Monte Carlo simulation of a dynamic MLC based on a multiple source model

Detailed knowledge of the characteristics of the radiation field shaped by a multileaf collimator (MLC) is essential in intensity modulated radiotherapy (IMRT). A previously developed multiple source model (MSM) for a 6 MV beam was extended to a 15 MV beam and supplemented with an accurate model of an 80-leaf dynamic MLC. Using the supplemented MSM and the MC code GEANT, lateral dose distributions were calculated in a water phantom and a portal water phantom. A field which is normally used for the validation of the step and shoot technique and a field from a realistic IMRT treatment plan delivered with dynamic MLC are investigated. To assess possible spectral changes caused by the modulation of beam intensity by an MLC, the energy spectra in five portal planes were calculated for moving slits of different widths. The extension of the MSM to 15 MV was validated by analysing energy fluences, depth doses and dose profiles. In addition, the MC-calculated primary energy spectrum was verified with an energy spectrum which was reconstructed from transmission measurements. MC-calculated dose profiles using the MSM for the step and shoot case and for the dynamic MLC case are in very good agreement with the measured data from film dosimetry. The investigation of a 13 cm wide field shows an increase in mean photon energy of up to 16% for the 0.25 cm slit compared to the open beam for 6 MV and of up to 6% for 15 MV, respectively. In conclusion, the MSM supplemented with the dynamic MLC has proven to be a powerful tool for investigational and benchmarking purposes or even for dose calculations in IMRT.

VMC++ versus BEAMnrc: a comparison of simulated linear accelerator heads for photon beams

BEAMnrc, a code for simulating medical linear accelerators based on EGSnrc, has been bench-marked and used extensively in the scientific literature and is therefore often considered to be the gold standard for Monte Carlo simulations for radiotherapy applications. However, its long computation times make it too slow for the clinical routine and often even for research purposes without a large investment in computing resources. VMC++ is a much faster code thanks to the intensive use of variance reduction techniques and a much faster implementation of the condensed history technique for charged particle transport. A research version of this code is also capable of simulating the full head of linear accelerators operated in photon mode (excluding multileaf collimators, hard and dynamic wedges). In this work, a validation of the full head simulation at 6 and 18 MV is performed, simulating with VMC++ and BEAMnrc the addition of one head component at a time and comparing the resulting phase space files. For the comparison, photon and electron fluence, photon energy fluence, mean energy, and photon spectra are considered. The largest absolute differences are found in the energy fluences. For all the simulations of the different head components, a very good agreement (differences in energy fluences between VMC++ and BEAMnrc <1%) is obtained. Only a particular case at 6 MV shows a somewhat larger energy fluence difference of 1.4%. Dosimetrically, these phase space differences imply an agreement between both codes at the <1% level, making VMC++ head module suitable for full head simulations with considerable gain in efficiency and without loss of accuracy.

The Genesis Solar Wind Concentrator: Flight and Post-Flight Conditions and Modeling of Instrumental Fractionation

The Genesis mission Solar Wind Concentrator was built to enhance fluences of solar wind by an average of 20x over the 2.3 years that the mission exposed substrates to the solar wind. The Concentrator targets survived the hard landing upon return to Earth and were used to determine the isotopic composition of solar-wind—and hence solar—oxygen and nitrogen. Here we report on the flight operation of the instrument and on simulations of its performance. Concentration and fractionation patterns obtained from simulations are given for He, Li, N, O, Ne, Mg, Si, S, and Ar in SiC targets, and are compared with measured concentrations and isotope ratios for the noble gases. Carbon is also modeled for a Si target. Predicted differences in instrumental fractionation between elements are discussed. Additionally, as the Concentrator was designed only for ions ≤22 AMU, implications of analyzing elements as heavy as argon are discussed. Post-flight simulations of instrumental fractionation as a function of radial position on the targets incorporate solar-wind velocity and angular distributions measured in flight, and predict fractionation patterns for various elements and isotopes of interest. A tighter angular distribution, mostly due to better spacecraft spin stability than assumed in pre-flight modeling, results in a steeper isotopic fractionation gradient between the center and the perimeter of the targets. Using the distribution of solar-wind velocities encountered during flight, which are higher than those used in pre-flight modeling, results in elemental abundance patterns slightly less peaked at the center. Mean fractionations trend with atomic mass, with differences relative to the measured isotopes of neon of +4.1±0.9 ‰/amu for Li, between -0.4 and +2.8 ‰/amu for C, +1.9±0.7‰/amu for N, +1.3±0.4 ‰/amu for O, -7.5±0.4 ‰/amu for Mg, -8.9±0.6 ‰/amu for Si, and -22.0±0.7 ‰/amu for S (uncertainties reflect Monte Carlo statistics). The slopes of the fractionation trends depend to first order only on the relative differential mass ratio, Δ m/ m. This article and a companion paper (Reisenfeld et al. 2012, this issue) provide post-flight information necessary for the analysis of the Genesis solar wind samples, and thus serve to complement the Space Science Review volume, The Genesis Mission (v. 105, 2003).

The abundance and isotopic composition of Cd in iron meteorites

Cadmium is a highly volatile element and its abundance in meteorites may help better understand volatility-controlled processes in the solar nebula and on meteorite parent bodies. The large thermal neutron capture cross section of 113Cd suggests that Cd isotopes might be well suited to quantify neutron fluences in extraterrestrial materials. The aims of this study were (1) to evaluate the range and magnitude of Cd concentrations in magmatic iron meteorites, and (2) to assess the potential of Cd isotopes as a neutron dosimeter for iron meteorites. Our new Cd concentration data determined by isotope dilution demonstrate that Cd concentrations in iron meteorites are significantly lower than in some previous studies. In contrast to large systematic variations in the concentration of moderately volatile elements like Ga and Ge, there is neither systematic variation in Cd concentration amongst troilites, nor amongst metal phases of different iron meteorite groups. Instead, Cd is strongly depleted in all iron meteorite groups, implying that the parent bodies accreted well above the condensation temperature of Cd (i.e., ≈650 K) and thus incorporated only minimal amounts of highly volatile elements. No Cd isotope anomalies were found, whereas Pt and W isotope anomalies for the same iron meteorite samples indicate a significant fluence of epithermal and higher energetic neutrons. This observation demonstrates that owing to the high Fe concentrations in iron meteorites, neutron capture mainly occurs at epithermal and higher energies. The combined Cd-Pt-W isotope results from this study thus demonstrate that the relative magnitude of neutron capture-induced isotope anomalies is strongly affected by the chemical composition of the irradiated material. The resulting low fluence of thermal neutrons in iron meteorites and their very low Cd concentrations make Cd isotopes unsuitable as a neutron dosimeter for iron meteorites.

MeV ion beam lithography of biocompatible halogenated Parylenes using aperture masks

Parylenes are poly(p-xylylene) polymers that are widely used as moisture barriers and in biomedicine because of their good biocompatibility. We have investigated MeV ion beam lithography using 16O+ ions for writing defined patterns in Parylene-C, which is evaluated as a coating material for the Cochlear Implant (CI) electrode array, a neuroprosthesis to treat some forms of deafness. Parylene-C and -F on silicon and glass substrates as well as 50 μm thick PTFE were irradiated to different fluences (1×1013-1×10161×1013-1×1016 1 MeV 16O+ ions cm−2) through aperture masks under high vacuum and a low pressure (<10−3 mbar) oxygen atmosphere. Biocompatibility of the irradiated and unirradiated surfaces was tested by cell-counting to determine the proliferation of murine spiral ganglion cells. The results reveal that an oxygen ion beam can be used to pattern Parylene-C and -F without using a liquid solvent developer in a similar manner to PTFE but with a ∼25× smaller removal rate. Biocompatibility tests showed no difference in cell adhesion between irradiated and unirradiated areas or ion fluence dependence. Coating the Parylene surface with an adhesion-promoting protein mixture had a much greater effect on cell proliferation.