Microdistribution and local dosimetry of 226Ra in trabecular bone of the beagle

Sections of lumbar vertebral bodies of young adult beagle dogs have been analyzed autoradiographically to characterize and quantify the local distribution of 226Ra by means of a scanning microscope photometer. The animals received a single injection of 355 kBq/kg body weight and were serially sacrificed at 5 to 1381 days postinjection. Hotspot concentrations decreased from about 51 kBq/g bone at 5 days to 20 kBq/g at 1381 days postinjection. The diffuse concentration changed from 8.3 to 1.9 kBq/g. The mean 226Ra concentration in the trabecular areas scanned was initially higher and at the end of the observation period lower than the average calculated for the whole lumbar vertebral column. Density and area of, and fraction of bone activity in, hotspots virtually remained constant. With time hotspots tended to become translocated into bone volume. Mean dose rates to lining cells from both hotspots and diffuse labels decreased from about 210 mGy/d at early postinjection times to 105 mGy/d. This corresponds to 2.5 to 1.1 times the average skeletal dose rate. A discussion of the level of irradiation in terms of hit frequencies shows that osteoblasts in the initial phase of hotspot formation receive about 60 hits to their nucleus for the duration of bone formation. After about 6 months, however, the 226Ra concentration in new bone and the corresponding hit frequency appears to be low enough that interference with bone formation is unlikely. Morphometric measurements showed that abnormal bone accretion and thickening of trabeculae occurred. This was interpreted as an imbalance between bone formation and resorption. Both formation and resorption seem to be substantially lowered compared to control animals.     

Radial secondary electron dose profiles and biological effects in light-ion beams based on analytical and Monte Carlo calculations using distorted wave cross sections

To speed up dose calculation, an analytical pencil-beam method has been developed to calculate the mean radial dose distributions due to secondary electrons that are set in motion by light ions in water. For comparison, radial dose profiles calculated using a Monte Carlo technique have also been determined. An accurate comparison of the resulting radial dose profiles of the Bragg peak for (1)H(+), (4)He(2+) and (6)Li(3+) ions has been performed. The double differential cross sections for secondary electron production were calculated using the continuous distorted wave-eikonal initial state method (CDW-EIS). For the secondary electrons that are generated, the radial dose distribution for the analytical case is based on the generalized Gaussian pencil-beam method and the central axis depth-dose distributions are calculated using the Monte Carlo code PENELOPE. In the Monte Carlo case, the PENELOPE code was used to calculate the whole radial dose profile based on CDW data. The present pencil-beam and Monte Carlo calculations agree well at all radii. A radial dose profile that is shallower at small radii and steeper at large radii than the conventional 1/r(2) is clearly seen with both the Monte Carlo and pencil-beam methods. As expected, since the projectile velocities are the same, the dose profiles of Bragg-peak ions of 0.5 MeV (1)H(+), 2 MeV (4)He(2+) and 3 MeV (6)Li(3+) are almost the same, with about 30% more delta electrons in the sub keV range from (4)He(2+)and (6)Li(3+) compared to (1)H(+). A similar behavior is also seen for 1 MeV (1)H(+), 4 MeV (4)He(2+) and 6 MeV (6)Li(3+), all classically expected to have the same secondary electron cross sections. The results are promising and indicate a fast and accurate way of calculating the mean radial dose profile.     

The effect of 238Pu alpha-particles on the mouse fibroblast cell line C3H 10T1/2: characterization of source and RBE for cell survival

Considerable interest has been aroused in recent years by reports that the transforming and carcinogenic effectiveness of low doses of high LET radiations can be increased by reducing the dose rate, especially for transformation of 10T1/2 cells in vitro by fission-spectrum neutrons. We report on conditions which have been established for irradiation of 10T1/2 cells with high LET monoenergetic alpha-particles (energy of 3.2 MeV, LET of 124 keV microns-1) from 238Pu. The alpha-particle irradiator allows convenient irradiation of multiple dishes of cells at selectable high or low dose rates and temperatures. The survival curves of irradiated cells showed that the mean lethal dose of alpha-particles was 0.6 Gy and corresponded to an RBE, at high dose rates, of 7.9 at 80 per cent survival and 4.6 at 5 per cent survival, relative to 60Co gamma-rays. The mean areas of the 10T1/2 nuclei, perpendicular to the incident alpha-particles, was measured as 201 microns2, from which it follows that, on average, only one in six of the alpha-particle traversals through a cell nucleus is lethal. Under the well-characterized conditions of these experiments the event frequency of alpha-particle traversals through cell nuclei is 9.8 Gy-1.     

Monte Carlo modeling of the effects of injected short-, medium- and longer-range alpha-particle emitters on human marrow at various ages

The effects of injected short-, medium- and longer-range alpha-particle emitters ((149)Tb, (211)At/(211)Po and (213)Bi/(213)Po, respectively) on the total hemopoietic stem cell population of active normal bone marrow in humans of various ages has been estimated using Monte Carlo modeling. The fraction of the normal hemopoietic stem cells that are hit and survive has been calculated as a first step toward estimating the risk of development of therapy-induced leukemia. The fraction was lowest for the shorter-range alpha-particle emitter ((149)Tb) and highest for the longer-range alpha-particle emitter ((213)Bi/(213)Po), with the value for the medium-range alpha-particle emitter (211)At/(211)Po being intermediate between these. There was little variation in the data with the age of the subject within each alpha-particle emitter. This lack of age dependence provides reassurance that the fraction of cells hit in any subject of any age with normal marrow can be estimated by modeling newborn marrow (which requires little computing time) despite age-related differences in microarchitecture.     

Impact of photon cross section uncertainties on Monte Carlo-determined depth-dose distributions

This work studies the impact of systematic uncertainties associated to interaction cross sections on depth dose curves determined by Monte Carlo simulations. The corresponding sensitivity factors are quantified by changing cross sections by a given amount and determining the variation in the dose. The influence of total and partial photon cross sections is addressed. Partial cross sections for Compton and Rayleigh scattering, photo-electric effect, and pair production have been accounted for. The PENELOPE code was used in all simulations. It was found that photon cross section sensitivity factors depend on depth. In addition, they are positive and negative for depths below and above an equilibrium depth, respectively. At this depth, sensitivity factors are null. The equilibrium depths found in this work agree very well with the mean free path of the corresponding incident photon energy. Using the sensitivity factors reported here, it is possible to estimate the impact of photon cross section uncertainties on the uncertainty of Monte Carlo-determined depth dose curves.     

Microdosimetry of electron microbeams

Track structures of 25, 50 and 80 keV primary electrons, simulated by the detailed-history Monte Carlo method, were analyzed for the frequency distributions of energy deposited in spheres with a diameter of 1 microm, placed in a cylindrically symmetrical array around the projected initial direction of the primary electron. The frequency mean of specific energy, the dose mean of lineal energy, and the parameters of lognormal functions fit to the dose distributions were calculated as a function of beam penetration and radial distance from the projected beam axis. Given these data, the stochastics of dose and radiation quality for micrometer-scale sites targeted by a medium-energy electron microbeam can be predicted as a function of the site's location relative to the beam entry point.     

Dosimetric characterization of an 192Ir brachytherapy source with the Monte Carlo code PENELOPE

Monte Carlo calculations are highly spread and settled practice to calculate brachytherapy sources dosimetric parameters. In this study, recommendations of the AAPM TG-43U1 report have been followed to characterize the Varisource VS2000 ¹⁹²Ir high dose rate source, provided by Varian Oncology Systems.In order to obtain dosimetric parameters for this source, Monte Carlo calculations with PENELOPE code have been carried out. TG-43 formalism parameters have been presented, i.e., air kerma strength, dose rate constant, radial dose function and anisotropy function. Besides, a 2D Cartesian coordinates dose rate in water table has been calculated. These quantities are compared to this source reference data, finding results in good agreement with them.The data in the present study complement published data in the next aspects: (i) TG-43U1 recommendations are followed regarding to phantom ambient conditions and to uncertainty analysis, including statistical (type A) and systematic (type B) contributions; (ii) PENELOPE code is benchmarked for this source; (iii) Monte Carlo calculation methodology differs from that usually published in the way to estimate absorbed dose, leaving out the track-length estimator; (iv) the results of the present work comply with the most recent AAPM and ESTRO physics committee recommendations about Monte Carlo techniques, in regards to dose rate uncertainty values and established differences between our results and reference data.The results stated in this paper provide a complete parameter collection, which can be used for dosimetric calculations as well as a means of comparison with other datasets from this source.     

A Study on the Dose Distributions in Various Materials from an Ir-192 HDR Brachytherapy Source

Dose distributions of (192)Ir HDR brachytherapy in phantoms simulating water, bone, lung tissue, water-lung and bone-lung interfaces using the Monte Carlo codes EGS4, FLUKA and MCNP4C are reported. Experiments were designed to gather point dose measurements to verify the Monte Carlo results using Gafchromic film, radiophotoluminescent glass dosimeter, solid water, bone, and lung phantom. The results for radial dose functions and anisotropy functions in solid water phantom were consistent with previously reported data (Williamson and Li). The radial dose functions in bone were affected more by depth than those in water. Dose differences between homogeneous solid water phantoms and solid water-lung interfaces ranged from 0.6% to 14.4%. The range between homogeneous bone phantoms and bone-lung interfaces was 4.1% to 15.7%. These results support the understanding in dose distribution differences in water, bone, lung, and their interfaces. Our conclusion is that clinical parameters did not provide dose calculation accuracy for different materials, thus suggesting that dose calculation of HDR treatment planning systems should take into account material density to improve overall treatment quality.     

Distributions of specific energy in sensitive layers of the human respiratory tract.

The purpose of the present work was to calculate the specific energy distribution in sensitive cells of the human lung. Specific energy distributions were calculated by applying Monte Carlo methods and the ICRP 66 model of the human respiratory tract. Specific energies were calculated at various depths in the epithelium and combined according to the relative cell abundance. Distributions are given for various combinations of sources, alpha-particle energies, targets and regions of the lung. The chord length does not follow the triangular distribution when the particle range is comparable to the diameter of the target. The notion of "effective volume" is introduced and defined, which is needed for estimation of hit frequency in some particular targets. It has been shown that basal cells are subjected to a larger proportion of alpha-particle hits with small energy transfer than secretory cells. Small energy transfer events that lead to minor damage of the DNA are more efficient in cancer induction than are hits with large energy deposition that lead to cell killing.     

Dosimetry of 252Cf sources for neutron radiotherapy with and without augmentation by boron neutron capture therapy.

Interstitial and intracavity 252Cf sources have been used to treat a number of tumor types with encouraging results. In particular these tumors include a variety of cervical, head-and-neck, and oral-cavity cancers and possible malignant gliomas. As a neutron source, 252Cf offers certain theoretical advantages over photon therapy (i.e., in treating tumors with significant hypoxic or necrotic components). With the recent availability of 10B-labeled tumor-seeking compounds, the usefulness of 252Cf may be further improved by augmenting the 252Cf dose to the tumor with an additional dose due to the fission (following thermal neutron capture) of 10B located in the tumor itself. While the high mean neutron energy permits 252Cf to deliver a high-LET, low-OER dose to the tumor on a macroscopic scale, thermalization of neutrons followed by 10B capture may augment this dose at the cellular level if adequate loading of tumor cells with 10B is possible. This paper presents results of a Monte Carlo simulation study investigating the dosimetric characteristics of linear 252Cf sources both with and without the quantitative increase in tumor dose possible with the addition of 10B. Results are displayed in the form of "along and away" tables and dose profiles in a water phantom. Comparisons of Monte Carlo results with experimental and analytical dosimetry data available in the literature are also presented.     

An EGS4 based Monte Carlo code for the calculation of organ equivalent dose to a modified Yale voxel phantom.

Organ or tissue equivalent dose, the most important quantity in radiation protection, cannot be measured directly. Therefore it became common practice to calculate the quantity of interest with Monte Carlo methods applied to so-called human phantoms, which are virtual representations of the human body. The Monte Carlo computer code determines conversion coefficients, which are ratios between organ or tissue equivalent dose and measurable quantities. Conversion coefficients have been published by the ICRP (Report No. 74) for various types of radiation, energies and fields, which have been calculated, among others, with the mathematical phantoms ADAM and EVA. Since then progress of image processing, and of clock speed and memory capacity of computers made it possible to create so-called voxel phantoms, which are a far more realistic representation of the human body. Voxel (Volume pixel) phantoms are built from segmented CT and/or MRI images of real persons. A complete set of such images can be joined to a 3-dimensional representation of the human body, which can be linked to a Monte Carlo code allowing for particle transport calculations. A modified version of the VOX_TISS8 human voxel phantom (Yale University) has been connected to the EGS4 Monte Carlo code. The paper explains the modifications, which have been made, the method of coupling the voxel phantom with the code, and presents results as conversion coefficients between organ equivalent dose and kerma in air for external photon radiation. A comparison of the results with published data shows good agreement.     

Lek mating system and its impact on male annihilation technique

Simulation models were presented which describe the matings inside a lek. The size of the lek was assumed to follow a poisson distribution, the dominance factor of each of the males was drawn from a gamma distribution, and the mating probability for each male was determined by its mating effectiveness raised to female mating factor divided by the sum of all the mating effectivenesses raised to this power. The mating probatility for each male and the actual matings were determined by Monte Carlo simulation. Based on the simulation data, the effects of female mating factor and size of the lek on the frequency of actual matings accomplished by the males inside the lek were discussed. In the case of male annihilation, the mean number of effective matings per male in the lek, calculated from the frequency of actual matings weighted by the effective mating ratio, was used to evaluate the efficiency of male annihilation method.     

A quantitative morphological analysis of macrophage stimulation

Summary The principles of stereology have been applied to a stereometric analysis of the ultrastructural composition of normal rat peritoneal macrophages and of cells stimulated five days previously by a single intraperitoneal dose of Freund's adjuvants. Material obtained by a systematic random sampling regime was analysed and the data groups for various morphological parameters were compared.Estimates were made of cell numbers. From electron micrographs the volume proportions of nuclei, mitochondria, heterogeneous granules and rough endoplasmic reticulum were determined. The average numbers and dimensions of mitochondria, granules and free ribosomes were also evaluated. The volume-to-surface ratios of cells and their nuclei were computed. By the independent measurement of cellular volumes it was possible to obtain estimates of nuclear and cytoplasmic volumes, to evaluate the membrane surface areas of the average cells and their nuclei, and to calculate numbers of organelles per average cell.The structural alterations induced by the complete adjuvant were similar to those evoked by the incomplete adjuvant. Stimulation was characterized by a cellular hypertrophy accompanied by an increase in mean granule size and in the number of mitochondria per average cell. At the same time, there was a substantial decrease in the number of granules and a considerable depletion of both nuclear and plasma membrane surface areas.     

Microbeams of heavy charged particles.

We have established a single cell irradiation system, which allows selected cells to be individually hit with defined number of heavy charged particles, using a collimated heavy-ion microbeam apparatus at JAERI-Takasaki. This system has been developed to study radiobiological processes in hit cells and bystander cells exposed to low dose and low dose-rate high-LET radiations, in ways that cannot be achieved using conventional broad-field exposures. Individual cultured cells grown in special dishes were irradiated in the atmosphere with a single or defined numbers of 18.3 MeV/amu 12C, 13.0 MeV/amu 20Ne, and 11.5 MeV/amu 40Ar ions. Targeting and irradiation of the cells were performed automatically at the on-line microscope of the microbeam apparatus according to the positional data of the target cells obtained at the off-line microscope before irradiation. The actual number of particle tracks that pass through cell nuclei was detected with prompt etching of the bottom of the cell dish made of ion track detector TNF-1 (modified CR-39), with alkaline-ethanol solution at 37 degrees C for 15-30 minutes. Using this system, separately inoculated Chinese hamster ovary cells, confluent normal human fibroblasts, and single plant cells (tobacco protoplasts) have been irradiated. These are the first studies in which single-ion direct hit effect and the bystander effect have been investigated using a high-LET heavy particle microbeam.     

Two-step validation of a Monte Carlo dosimetry framework for general radiology

The Monte Carlo technique is considered gold standard when it comes to patient-specific dosimetry. Any newly developed Monte Carlo simulation framework, however, has to be carefully calibrated and validated prior to its use. For many researchers this is a tedious work. We propose a two-step validation procedure for our newly built Monte Carlo framework and provide all input data to make it feasible for future related application by the wider community. The validation was at first performed by benchmarking against simulation data available in literature. The American Association of Physicists in Medicine (AAPM) report of task group 195 (case 2) was considered most appropriate for our application. Secondly, the framework was calibrated and validated against experimental measurements for trunk X-ray imaging protocols using a water phantom. The dose results obtained from all simulations and measurements were compared. Our Monte Carlo framework proved to agree with literature data, by showing a maximal difference below 4% to the AAPM report. The mean difference with the water phantom measurements was around 7%. The statistical uncertainty for clinical applications of the dosimetry model is expected to be within 10%. This makes it reliable for clinical dose calculations in general radiology. Input data and the described procedure allow for the validation of other Monte Carlo frameworks.     

Test of the Inverse Monte Carlo Method for the Calculation of Interatomic Potential Energies in Atomic Liquids

Abstract

The principle purpose of this paper is to demonstrate the use of the Inverse Monte Carlo technique for calculating pair interaction energies in monoatomic liquids from a given equilibrium property. This method is based on the mathematical relation between transition probability and pair potential given by the fundamental equation of the “importance sampling” Monte Carlo method. In order to have well defined conditions for the test of the Inverse Monte Carlo method a Metropolis Monte Carlo simulation of a Lennard Jones liquid is carried out to give the equilibrium pair correlation function determined by the assumed potential. Because an equilibrium configuration is prerequisite for an Inverse Monte Carlo simulation a model system is generated reproducing the pair correlation function, which has been calculated by the Metropolis Monte Carlo simulation and therefore representing the system in thermal equilibrium. This configuration is used to simulate virtual atom displacements. The resulting changes in atom distribution for each single simulation step are inserted in a set of non-linear equations defining the transition probability for the virtual change of configuration. The solution of the set of equations for pair interaction energies yields the Lennard Jones potential by which the equilibrium configuration has been determined.     

Radiation dosimetry by automatic image analysis of dicentric chromosomes.

A system for scoring dicentric chromosomes by image analysis comprised fully automatic location of mitotic cells, automatic retrieval, focus and digitization at high resolution, automatic rejection of nuclei and debris and detection and segmentation of chromosome clusters, automatic centromere location, and subsequent rapid interactive visual review of potential dicentric chromosomes to confirm positives and reject false positives. A calibration set of about 15,000 cells was used to establish the quadratic dose response for 60Co gamma-irradiation. The dose-response function parameters were established by a maximum likelihood technique, and confidence limits in the dose response and in the corresponding inverse curve, of estimated dose for observed dicentric frequency, were established by Monte Carlo techniques. The system was validated in a blind trial by analysing a test set comprising a total of about 8000 cells irradiated to 1 of 10 dose levels, and estimating the doses from the observed dicentric frequency. There was a close correspondence between the estimated and true doses. The overall sensitivity of the system in terms of the proportion of the total population of dicentrics present in the cells analysed that were detected by the system was measured to be about 40%. This implies that about 2.5 times more cells must be analysed by machine than by visual analysis. Taking this factor into account, the measured review time and false positive rates imply that analysis by the system of sufficient cells to provide the equivalent of a visual analysis of 500 cells would require about 1 h for operator review.     

Apoptotic bone cells may be engulfed by osteoclasts during alveolar bone resorption in young rats

The alveolar bone is a suitable in vivo physiological model for the study of apoptosis and interactions of bone cells because it undergoes continuous, rapid and intense resorption/remodelling, during a long period of time, to accommodate the growing tooth germs. The intensity of alveolar bone resorption greatly enhances the chances of observing images of the extremely rapid events of apoptosis of bone cells and also of images of interactions between osteoclasts and osteocytes/osteoblasts/bone lining cells. To find such images, we have therefore examined the alveolar bone of young rats using light microscopy, the TUNEL method for apoptosis, and electron microscopy. Fragments of alveolar bone from young rats were fixed in Bouin and formaldehyde for morphology and for the TUNEL method. Glutaraldehyde-formaldehyde fixed specimens were processed for transmission electron microscopy. Results showed TUNEL positive round/ovoid structures on the bone surface and inside osteocytic lacunae. These structures--also stained by hematoxylin--were therefore interpreted, respectively, as osteoblasts/lining cells and osteocytes undergoing apoptosis. Osteoclasts also exhibited TUNEL positive apoptotic bodies inside large vacuoles; the nuclei of osteoclasts, however, were always TUNEL negative. Ultrathin sections revealed typical apoptotic images--round/ ovoid bodies with dense crescent-like chromatin--on the bone surface, corresponding therefore to apoptotic osteoblasts/lining cells. Osteocytes also showed images compatible with apoptosis. Large osteoclast vacuoles often contained fragmented cellular material. Our results provide further support for the idea that osteoclasts internalize dying bone cells; we were however, unable to find images of osteoclasts in apoptosis.     

Proton-HZE-particle sequential dual-beam exposures increase anchorage-independent growth frequencies in primary human fibroblasts

The radiation field in deep space contains high levels of high-energy protons and substantially lower levels of high-atomic-number, high-energy (HZE) particles. Calculations indicate that cellular nuclei of human space travelers will be hit during a 3-year Mars mission by approximately 400 protons and approximately 0.4 HZE particles. Thus most cells in astronauts will be hit by a proton(s) before being hit by an HZE particle. To investigate effects of dual ion irradiations on human cells, we irradiated primary human neonatal fibroblasts with protons (1 GeV/nucleon, 20 cGy) followed from 2.5 min to 48 h later by iron or titanium ions (1 GeV/nucleon, 20 cGy) and then measured clonogenic survival and frequency of anchorage-independent growth. This frequency depends on the interval between hydrogen- and iron-ion irradiation, with a critical window between 2.5 min and 1 h producing about three times more anchorage-independent colonies per survivor than expected from simple addition of the two ions separately. The hydrogen-titanium-ion dual-beam irradiation produced similar increases that persisted to approximately 6 h. At longer intervals, anchorage-independent growth frequencies were similar to those expected for additivity. However, irradiation of cells with either an iron or a titanium particle first followed by protons produced only additive levels.