Monte Carlo Range Research Articles

A systematic approach for calibrating a Monte Carlo code to a treatment planning system for obtaining dose, LET, variable proton RBE and out-of-field dose

Objective. While integration of variable relative biological effectiveness (RBE) has not reached full clinical implementation, the importance of having the ability to recalculate proton treatment plans in a flexible, dedicated Monte Carlo (MC) code cannot be understated . Here we provide a step-wise method for calibrating dose from a MC code to a treatment planning system (TPS), to obtain required parameters for calculating linear energy transfer (LET), variable RBE and in general enabling clinical realistic research studies beyond the capabilities of a TPS. Approach. Initially, Pristine Bragg peaks (PBP) were calculated in both the Eclipse TPS and the FLUKA MC code. A rearranged Bortfeld energy-range relation was applied to the initial energy of the beam to fine-tune the range of the MC code at 80% dose level distal to the PBP. The energy spread was adapted by dividing the TPS range by the MC range for dose level 80%–20% distal to the PBP. Density and relative proton stopping power were adjusted by comparing the TPS and MC for different Hounsfield units. To find the relationship of dose per primary particle from the MC to dose per monitor unit in the TPS, integration was applied to the area of the Bragg curve. The calibration was validated for spread-out Bragg peaks (SOBP) in water and patient treatment plans. Following the validation, variable RBE were calculated using established models. Main results. The PBPs ranges were within ± mm threshold, and a maximum of 5.5% difference for the SOBPs was observed. The patient validation showed excellent dose agreement between the TPS and MC, with the greatest differences for the lung tumor patient. Significance. A procedure for calibrating a MC code to a TPS was developed and validated. The procedure enables MC-based calculation of dose, LET, variable RBE, advanced (secondary) particle tracking and more from treatment plans.

Lacustrine leaf wax hydrogen isotopes indicate strong regional climate feedbacks in Beringia since the last ice age

The Late-Quaternary climate of Beringia remains unresolved despite the region's role in modulating glacial-interglacial climate and as the likely conduit for human dispersal into the Americas. Here, we investigate Beringian temperature change using an ∼32,000-year lacustrine record of leaf wax hydrogen isotope ratios (δ2Hwax) from Arctic Alaska. Based on Monte Carlo iterations accounting for multiple sources of uncertainty, the reconstructed summertime temperatures were ∼3 °C colder (range: −8 to +3 °C) during the Last Glacial Maximum (LGM; 21-25 ka) than the pre-industrial era (PI; 2–0.1 ka). This ice-age summer cooling is substantially smaller than in other parts of the Arctic, reflecting altered atmospheric circulation and increased continentality which weakened glacial cooling in the region. Deglacial warming was punctuated by abrupt events that are largely synchronous with events seen in Greenland ice cores that originate in the North Atlantic but which are also controlled locally, such as by the opening of the Bering Strait between 13.4 and 11 ka. Our reconstruction, together with climate modeling experiments, indicates that Beringia responds more strongly to North Atlantic freshwater forcing under modern-day, open-Bering Strait conditions than under glacial conditions. Furthermore, a 2 °C increase (Monte Carlo range: −1 to +5 °C) over the anthropogenic era reverses a 6 °C decline (Monte Carlo range: −10 to 0 °C) through the Holocene, indicating that recent warming in Arctic Alaska has not surpassed peak Holocene summer warmth.

Range versus density dependence of proton-therapy beams in human tissues

Abstract Range versus density dependence of proton-therapy beams in human tissues is discussed in the paper. The main attention is paid to energy-scaling of simple fitting functions that fit the range-to-density dependence for single proton-beam energies. These functions have been found by fitting ranges obtained with the aid of the Monte Carlo module of SRIM2013. The simulations have been run for 28 human tissues at 8 energies from 60 MeV to 220 MeV. It has been shown that a fitting function that depends solely on the target density can be found at each proton-beam energy. After the energy-scaling, the fitted ranges have been compared with the Monte Carlo ranges. We conclude that the final energy-scaled fitting function provides ranges within 0.4% of the Monte Carlo range (on average), which is less than 30% of the natural range-straggling (1.4% of the Monte Carlo range on average). The worst data-point has an off-set from the fitting function less than 1.7% of the Monte Carlo range.

SU‐E‐CAMPUS‐T‐02: Exploring Radiation Acoustics CT Dosimeter Design Aspects for Proton Therapy

Purpose:Investigate the design aspects and imaging dose capabilities of the Radiation Acoustics Computed Tomography (RA CT) dosimeter for Proton induced acoustics, with the objective to characterize a pulsed pencil proton beam. The focus includes scanner geometry, transducer array, and transducer bandwidth on image quality.Methods:The geometry of the dosimeter is a cylindrical water phantom (length 40cm, radius 15cm) with 71 ultrasound transducers placed along the length and end of the cylinder to achieve a weighted set of projections with spherical sampling. A 3D filtered backprojection algorithm was used to reconstruct the dosimetric images and compared to MC dose distribution. First, 3D Monte Carlo (MC) Dose distributions for proton beam energies (range of 12cm, 16cm, 20cm, and 27cm) were used to simulate the acoustic pressure signal within this scanner for a pulsed proton beam of 1.8×107 protons, with a pulse width of 1 microsecond and a rise time of 0.1 microseconds. Dose comparison within the Bragg peak and distal edge were compared to MC analysis, where the integrated Gaussian was used to locate the 50% dose of the distal edge. To evaluate spatial fidelity, a set of point sources within the scanner field of view (15×15×15cm3) were simulated implementing a low‐pass bandwidth response function (0 to 1MHz) equivalent to a multiple frequency transducer array, and the FWHM of the point‐spread‐function determined.Results:From the reconstructed images, RACT and MC range values are within 0.5mm, and the average variation of the dose within the Bragg peak are within 2%. The spatial resolution tracked with transducer bandwidth and projection angle sampling, and can be kept at 1.5mm.Conclusion:This design is ready for fabrication to start acquiring measurements. The 15 cm FOV is an optimum size for imaging dosimetry. Currently, simulations comparing transducer sensitivity, bandwidth, and proton beam parameters are being evaluated to assess signal‐to‐noise.

Variable-energy positron-beam studies of Ni implanted with He.

Variable-energy positron-beam studies have been made on well-annealed polycrystalline Ni samples implanted with 30-, 90-, and 180-keV $^{4}\mathrm{He}$ ions. The positron-annihilation characteristics were measured with a solid-state Ge detector at a number of different incident-positron energies and after isochronal annealing at various temperatures. The Doppler broadening of the annihilation photons was found to be strongly influenced by the $^{4}\mathrm{He}$ implantations. The data indicate that trapping of the positrons occurred predominantly at small He bubbles. The variation of the broadening with incident-positron energy was sensitive to the depth distribution of the traps. A diffusion model assuming a square concentration-defect profile was developed and analytically fitted to the parametrized momentum data. These fitted results were compared to Monte Carlo range calculations for $^{4}\mathrm{He}$ in Ni, and fairly good agreement was found. This investigation demonstrates the capabilities of positron annihilation for nondestructive depth profiling in ion-implanted systems. In addition, it establishes parallels between the trapping behavior of positrons and that reported elsewhere for hydrogen, thereby augmenting the present level of understanding of the technologically important trapping of hydrogen by the bubbles.

Near‐Surface Damage and Contamination after CF 4 / H 2 Reactive Ion Etching of Si

Silicon surfaces which had been exposed to reactive ion etching have been characterized by x‐ray photoelectron emission spectroscopy, He ion channeling, H profiling, and Raman scattering techniques. Plasma exposure of a clean Si surface leads to the deposition of a thin (∼30–50Å thick) C,F containing film. The near‐surface region (∼30–50Å) of the Si substrate is heavily disordered (“amorphized”), as found by ion channeling and Raman scattering. A modified, less damaged Si region extends from about 30–50Å from the surface to a depth in excess of 250Å. This layer contains a high concentration (∼5 atom percent) of H as shown by hydrogen profiling techniques. From the observation of Si‐H and vibrational modes by Raman scattering, it has been shown that the H is bonded to the Si lattice. In ion scattering, the extended modified Si layer appears to cause a strongly enhanced background in the energy spectra. Results of Monte Carlo range calculations are reported and compared with the damage depth found experimentally.