Title:
Study of the secondary neutral radiation in proton therapy: toward an indirect in vivo dosimetry.

Abstract:
Secondary particles produced in the collision of protons with beam modifiers are of concern in proton therapy. Nevertheless, secondary radiation can provide information on the dosimetric parameters through its dependency on the modulating accessories (range shifter and range modulating wheel). Relatively little data have been reported in the literature for low-energy proton beams. The present study aims at characterizing the neutron and photon secondary radiation at the low-energy proton therapy facility of the Centre Antoine Lacassagne (CAL), and studying their correlation to the dosimetric parameters to explore possible practical uses of secondary radiation in the treatment quality for proton therapy.
The Monte Carlo code MCNPX was used to simulate the proton therapy facility at CAL. Neutron and photon fluence, Φ, and ambient dose equivalent per proton dose, H∗(10)∕D, were determined across the horizontal main plane spanning the whole treatment room. H∗(10)∕D was also calculated at two positions of the treatment room where dosimetric measurements were performed for validation of the Monte Carlo calculations. Calculations and measurements were extended to 100 clinical spread-out Bragg Peaks (SOBPs) covering the whole range of therapeutic dose rates (D∕MU) employed at CAL. In addition, the values of D and MU were also calculated for each SOBP and the results analyzed to study the relationship between secondary radiation and dosimetric parameters.
The largest production of the secondary particles takes place at the modulating devices and the brass collimators located along the optical bench. Along the beam line and off the beam axis to 2.5 m away, H∗(10)∕D values ranged from 5.4 μSv∕Gy to 5.3 mSv∕Gy for neutrons, and were 1 order of magnitude lower for photons. H∗(10)∕D varied greatly with the distance and angle to the beam axis. A variation of a factor of 5 was found for the different range of modulations (SOBPs). The ratios between calculations and measurements were 2.3 and 0.5 for neutrons and photons, respectively, and remained constant for all the range of SOBPs studied, which provided validation for the Monte Carlo calculations. H∗(10)∕D values were found to correlate to the proton dose rate D∕MU with a power fit, both for neutrons and photons. This result was exploited to implement a system to obtain D∕MU values from the measurement of the integrated photon ambient dose equivalent H∗(10) during treatment, which provides a method to control the dosimetric parameters D∕MU and D.
The treatment room at CAL is moderately polluted by secondary particles. The constant ratio between measurements and calculations for all SOBPs showed that simulations correctly predict the dosimetric parameters and the dependence of the production of secondary particles on the modulation. The correlation between H∗(10)∕D and D∕MU is a useful tool for quality control and is currently used at CAL. This system works as an indirect in vivo dosimetry method, which is so far not feasible in proton therapy. This tool requires very simple instrumentation and can be implemented from the measurement of either photons or neutrons.