Monte Carlo Modeling of a Protontherapy Beam Line Dedicated to Ophthalmologic Treatments

Configuration and validation of an analytical model predicting secondary neutron radiation in proton therapy using Monte Carlo simulations and experimental measurements

PURPOSE

This study focuses on the configuration and validation of an analytical model predicting leakage neutron doses in proton therapy.

METHODS

Using Monte Carlo (MC) calculations, a facility-specific analytical model was built to reproduce out-of-field neutron doses while separately accounting for the contribution of intra-nuclear cascade, evaporation, epithermal and thermal neutrons. This model was first trained to reproduce in-water neutron absorbed doses and in-air neutron ambient dose equivalents, H*(10), calculated using MCNPX. Its capacity in predicting out-of-field doses at any position not involved in the training phase was also checked. The model was next expanded to enable a full 3D mapping of H*(10) inside the treatment room, tested in a clinically relevant configuration and finally consolidated with experimental measurements.

RESULTS

Following the literature approach, the work first proved that it is possible to build a facility-specific analytical model that efficiently reproduces in-water neutron doses and in-air H*(10) values with a maximum difference less than 25%. In addition, the analytical model succeeded in predicting out-of-field neutron doses in the lateral and vertical direction. Testing the analytical model in clinical configurations proved the need to separate the contribution of internal and external neutrons. The impact of modulation width on stray neutrons was found to be easily adjustable while beam collimation remains a challenging issue. Finally, the model performance agreed with experimental measurements with satisfactory results considering measurement and simulation uncertainties.

CONCLUSION

Analytical models represent a promising solution that substitutes for time-consuming MC calculations when assessing doses to healthy organs.

Secondary neutron doses in proton therapy treatments of ocular melanoma and craniopharyngioma

Monte Carlo simulations were used to assess secondary neutron doses received by patients treated with proton therapy for ocular melanoma and craniopharyngioma. MCNPX calculations of out-of-field doses were done for ∼20 different organs considering realistic treatment plans and using computational phantoms representative of an adult male individual. Simulations showed higher secondary neutron doses for intracranial treatments, ∼14 mGy to the salivary glands, when compared with ocular treatments, ∼0.6 mGy to the non-treated eye. This secondary dose increase is mainly due to the higher proton beam energy (178 vs. 75 MeV) as well as to the impact of the different beam parameters (modulation, collimation, field size etc.). Moreover, when compared with published data, the assessed secondary neutron doses showed similar trends, but sometimes with sensitive differences. This confirms secondary neutrons to be directly dependent on beam energy, modulation technique, treatment configuration and methodology.

Secondary neutron doses received by paediatric patients during intracranial proton therapy treatments

This paper's goal is to assess secondary neutron doses received by paediatric patients treated for intracranial tumours using a 178 MeV proton beam. The MCNPX Monte Carlo model of the proton therapy facility, previously validated through experimental measurements for both proton and neutron dosimetry, was used. First, absorbed dose was calculated for organs located outside the clinical target volume using a series of hybrid computational phantoms for different ages and considering a realistic treatment plan. In general, secondary neutron dose was found to decrease as the distance to the treatment field increases and as the patient age increases. In addition, secondary neutron doses were studied as a function of the beam incidence. Next, neutron equivalent dose was assessed using organ-specific energy-dependent radiation weighting factors determined from Monte Carlo simulations of neutron spectra at each organ. The equivalent dose was found to reach a maximum value of ∼155 mSv at the level of the breasts for a delivery of 49 proton Gy to an intracranial tumour of a one-year-old female patient. Finally, a thorough comparison of the calculation results with published data demonstrated the dependence of neutron dose on the treatment configuration and proved the need for facility-specific and treatment-dependent neutron dose calculations.