Prompt gamma-ray emission from biological tissues during proton irradiation: a preliminary study.

In this paper, we present the results of a preliminary study of secondary 'prompt' gamma-ray emission produced by proton-nuclear interactions within tissue during proton radiotherapy. Monte Carlo simulations were performed for mono-energetic proton beams, ranging from 2.5 MeV to 250 MeV, irradiating elemental and tissue targets. Calculations of the emission spectra from different biological tissues and their elemental components were made. Also, prompt gamma rays emitted during delivery of a clinical proton spread-out Bragg peak (SOBP) in a homogeneous water phantom and a water phantom containing heterogeneous tissue inserts were calculated to study the correlation between prompt gamma-ray production and proton dose delivery. The results show that the prompt gamma-ray spectra differ significantly for each type of tissue studied. The relative intensity of the characteristic gamma rays emitted from a given tissue was shown to be proportional to the concentration of each element in that tissue. A strong correlation was found between the delivered SOBP dose distribution and the characteristic prompt gamma-ray production. Based on these results, we discuss the potential use of prompt gamma-ray emission as a method to verify the accuracy and efficacy of doses delivered with proton radiotherapy.

Experimental validation of a Monte Carlo proton therapy nozzle model incorporating magnetically steered protons.

The purpose of this study is to validate the accuracy of a Monte Carlo calculation model of a proton magnetic beam scanning delivery nozzle developed using the Geant4 toolkit. The Monte Carlo model was used to produce depth dose and lateral profiles, which were compared to data measured in the clinical scanning treatment nozzle at several energies. Comparisons were also made between measured and simulated off-axis profiles to test the accuracy of the model's magnetic steering. Comparison of the 80% distal dose fall-off values for the measured and simulated depth dose profiles agreed to within 1 mm for the beam energies evaluated. Agreement of the full width at half maximum values for the measured and simulated lateral fluence profiles was within 1.3 mm for all energies. The position of measured and simulated spot positions for the magnetically steered beams agreed to within 0.7 mm of each other. Based on these results, we found that the Geant4 Monte Carlo model of the beam scanning nozzle has the ability to accurately predict depth dose profiles, lateral profiles perpendicular to the beam axis and magnetic steering of a proton beam during beam scanning proton therapy.

Commissioning dose computation models for spot scanning proton beams in water for a commercially available treatment planning system.

PURPOSE: To present our method and experience in commissioning dose models in water for spot scanning proton therapy in a commercial treatment planning system (TPS). METHODS: The input data required by the TPS included in-air transverse profiles and integral depth doses (IDDs). All input data were obtained from Monte Carlo (MC) simulations that had been validated by measurements. MC-generated IDDs were converted to units of Gy mm(2)/MU using the measured IDDs at a depth of 2 cm employing the largest commercially available parallel-plate ionization chamber. The sensitive area of the chamber was insufficient to fully encompass the entire lateral dose deposited at depth by a pencil beam (spot). To correct for the detector size, correction factors as a function of proton energy were defined and determined using MC. The fluence of individual spots was initially modeled as a single Gaussian (SG) function and later as a double Gaussian (DG) function. The DG fluence model was introduced to account for the spot fluence due to contributions of large angle scattering from the devices within the scanning nozzle, especially from the spot profile monitor. To validate the DG fluence model, we compared calculations and measurements, including doses at the center of spread out Bragg peaks (SOBPs) as a function of nominal field size, range, and SOBP width, lateral dose profiles, and depth doses for different widths of SOBP. Dose models were validated extensively with patient treatment field-specific measurements. RESULTS: We demonstrated that the DG fluence model is necessary for predicting the field size dependence of dose distributions. With this model, the calculated doses at the center of SOBPs as a function of nominal field size, range, and SOBP width, lateral dose profiles and depth doses for rectangular target volumes agreed well with respective measured values. With the DG fluence model for our scanning proton beam line, we successfully treated more than 500 patients from March 2010 through June 2012 with acceptable agreement between TPS calculated and measured dose distributions. However, the current dose model still has limitations in predicting field size dependence of doses at some intermediate depths of proton beams with high energies. CONCLUSIONS: We have commissioned a DG fluence model for clinical use. It is demonstrated that the DG fluence model is significantly more accurate than the SG fluence model. However, some deficiencies in modeling the low-dose envelope in the current dose algorithm still exist. Further improvements to the current dose algorithm are needed. The method presented here should be useful for commissioning pencil beam dose algorithms in new versions of TPS in the future.