Independent application of an analytical model for secondary neutron equivalent dose produced in a passive-scattering proton therapy treatment unit

Stray neutron radiation exposures from proton therapy: physics-based analytical models of neutron spectral fluence, kerma and absorbed dose

Objective. Patients who receive proton beam therapy are exposed to unwanted stray neutrons. Stray radiations increase the risk of late effects in normal tissues, such as second cancers and cataracts, and may cause implanted devices such as pacemakers to malfunction. Compared to therapeutic beams, little attention has been paid to modeling stray neutron exposures. In the past decade, substantial progress was made to develop semiempirical models of stray neutron dose equivalent, but models to routinely calculate neutron absorbed dose and kerma are still lacking. The objective of this work was to develop a new physics based analytical model to calculate neutron spectral fluence, kerma, and absorbed dose in a water phantom. Approach. We developed the model using dosimetric data from Monte Carlo simulations and neutron kerma coefficients from the literature. The model explicitly considers the production, divergence, scattering, and attenuation of neutrons. Neutron production was modeled for 120–250 MeV proton beams impinging on a variety of materials. Fluence, kerma and dose calculations were performed in a 30 × 180 × 44 cm3 phantom at points up to 43 cm in depth and 80 cm laterally. Main Results. Predictions of the analytical model agreed reasonably with corresponding values from Monte Carlo simulations, with a mean difference in average energy deposited of 20%, average kerma coefficient of 21%, and absorbed dose to water of 49%. Significance. The analytical model is simple to implement and use, requires less configuration data that previously reported models, and is computationally fast. This model appears potentially suitable for integration in treatment planning system, which would enable risk calculations in prospective and retrospective cases, providing a powerful tool for epidemiological studies and clinical trials.

Determining Out-of-Field Doses and Second Cancer Risk From Proton Therapy in Young Patients—An Overview

Proton therapy has the potential to provide survival and tumor control outcomes comparable and frequently superior to photon therapy. This has led to a significant concern in the medical physics community on the risk for the induction of second cancers in all patients and especially in younger patients, as they are considered more radiosensitive than adults and have an even longer expected lifetime after treatment. Thus, our purpose is to present an overview of the research carried out on the evaluation of out-of-field doses linked to second cancer induction and the prediction of this risk. Most investigations consisted of Monte Carlo simulations in passive beam facilities for clinical scenarios. These works established that equivalent doses in organs could be up to 200 mSv or 900 mSv for a brain or a craniospinal treatment, respectively. The major contribution to this dose comes from the secondary neutrons produced in the beam line elements. Few works focused on scanned-beam facilities, but available data show that, for these facilities, equivalent doses could be between 2 and 50 times lower. Patient age is a relevant factor in the dose level, especially for younger patients (by means of the size of the body) and, in addition, in the predicted risk by models (due to the age dependence of the radiosensitivity). For risks, the sex of the patient also plays an important role, as female patients show higher sensitivity to radiation. Thus, predicted risks of craniospinal irradiation can range from 8% for a 15-year-old male patient to 58% for a 2-year-old female patient, using a risk model from a radiological protection field. These values must be taken with caution due to uncertainties in risk models, and then dosimetric evaluation of stray radiation becomes mandatory in order to complement epidemiological studies and be able to model appropriate dose–response functions for this dose range. In this sense, analytical models represent a useful tool and some models have been implemented to be used for young patients. Research carried out so far confirmed that proton beam therapy reduces the out-of-field doses and second cancer risk. However, further investigations may be required in scanned-beam delivery systems.

Proton Radiotherapy Could Reduce the Risk of Fatal Second Cancers for Children with Intracranial Tumors in Low- and Middle-Income Countries

PURPOSE

To test our hypothesis that, for young children with intracranial tumors, proton radiotherapy in a high-income country does not reduce the risk of a fatal subsequent malignant neoplasm (SMN) compared with photon radiotherapy in low- and middle-income countries.

MATERIALS AND METHODS

We retrospectively selected 9 pediatric patients with low-grade brain tumors who were treated with 3-dimensional conformal radiation therapy in low- and middle-income countries. Images and contours were deidentified and transferred to a high-income country proton therapy center. Clinically commissioned treatment planning systems of each academic hospital were used to calculate absorbed dose from the therapeutic fields. After fusing supplemental computational phantoms to the patients' anatomies, models from the literature were applied to calculate stray radiation doses. Equivalent doses were determined in organs and tissues at risk of SMNs, and the lifetime attributable risk of SMN mortality (LAR) was predicted using a dose-effect model. Our hypothesis test was based on the average of the ratios of LARs from proton therapy to that of photon therapy ()(H0: = 1; H A : < 1).

RESULTS

Proton therapy reduced the equivalent dose in organs at risk for SMNs and LARs compared with photon therapy for which the for the cohort was 0.69 ± 0.10, resulting in the rejection of H0 (P < .001, α = 0.05). We observed that the younger children in the cohort (2-4 years old) were at a factor of approximately 2.5 higher LAR compared with the older children (8-12 years old).

CONCLUSION

Our findings suggest that proton radiotherapy has the strong potential of reducing the risk of fatal SMNs in pediatric patients with intracranial tumors if it were made available globally.

A Monte Carlo based analytic model of the in-room neutron ambient dose equivalent for a Mevion gantry-mounted passively scattered proton system

The goal of this study was to develop a Monte Carlo (MC)-based analytical model that can predict the in-room ambient dose equivalent from a Mevion gantry-mounted passively scattered proton system. The Mevion S250 and treatment vault were simulated using the MCNPX MC code. The results of the in-room neutron dose measurements, using an FHT 762 WENDI-II detector, were employed to benchmark the MC-derived values. After tuning the MCNPX MC code, for the same beam delivery parameters, the code was used to calculate the neutron spectra and ambient dose equivalent in the vault and at varying angles from the isocenter. Then, based on the calculations, an analytical model was reconstructed and data were fitted to derive the model parameters at 95% confidence intervals (CI). The MCNPX codes were tuned to within about 19% of the measured values for most of the measurements in the vault. For the maze, up to 0.08 mSv Gy-1 discrepancies were found between the experimental measurements and MCNPX calculated results. The analytical model showed up to 18% discrepancy for distances between 100 and 600 cm from the isocenter compared to the MC calculations. The model may underestimate the neutron ambient dose equivalent up to 21% for distances less than 100 cm from the isocenter. The proposed analytical model can be used to estimate the contribution of the secondary neutron dose from the Mevion S250 for the design of local shielding inside the proton therapy treatment vault.

A Monte Carlo-based analytic model of neutron dose equivalent for a mevion gantry-mounted passively scattered proton system for craniospinal irradiation

PURPOSE

To calculate in- and out-of-field neutron spectra and dose equivalent, using Monte Carlo (MC) simulation, for a Mevion gantry-mounted passively scattered proton system in craniospinal irradiation. An analytical model based on the MC calculations that estimates in- and out-of-field neutron dose equivalent from proton Craniospinal irradiation (CSI) was also developed.

METHODS

The MCNPX MC code was used to simulate a Mevion S250 proton therapy system. The simulated proton depth doses and profiles for pristine and spread-out Bragg peaks were benchmarked against the measured data. Previous measurements using extended-range Bonner spheres were used to verify the calculated neutron spectra and dose equivalent. Using the benchmarked results as a reference condition, a correction-based analytical model was reconstructed by fitting the data to derive model parameters at 95% confidence interval. Sensitivity analysis of brass aperture opening, thickness of the Lucite (PMMA) range compensator, and modulation width was performed to obtain correction parameters for nonreference conditions.

RESULTS

For the neutron dose equivalent per therapeutic proton dose, the MCNPX calculated dose equivalent matched the measured values to within 8%. The benchmarked neutron dose equivalent at the isocenter was 41.2 and 20.8 mSv/Gy, for cranial and spinal fields, respectively. For in- and out-of-field neutron dose calculations, the correction-based analytical model showed up to 17% discrepancy compared to the MC calculations. The correction factors may provide a conservative estimation of neutron dose, especially for depth ≤ 5 cm and regions underneath the brass aperture.

CONCLUSION

The proposed analytical model can be used to estimate the contribution of the neutron dose to the overall CSI treatment dose. Moreover, the model can be employed to estimate the neutron dose to the implantable cardiac electronic devices.

Neutron dose and its measurement in proton therapy-current State of Knowledge

Proton therapy has shown dosimetric advantages over conventional radiation therapy using photons. Although the integral dose for patients treated with proton therapy is low, concerns were raised about late effects like secondary cancer caused by dose depositions far away from the treated area. This is especially true for neutrons and therefore the stray dose contribution from neutrons in proton therapy is still being investigated. The higher biological effectiveness of neutrons compared to photons is the main cause of these concerns. The gold-standard in neutron dosimetry is measurements, but performing neutron measurements is challenging. Different approaches have been taken to overcome these difficulties, for instance with newly developed neutron detectors. Monte Carlo simulations is another common technique to assess the dose from secondary neutrons. Measurements and simulations are used to develop analytical models for fast neutron dose estimations. This article tries to summarize the developments in the different aspects of neutron dose in proton therapy since 2017. In general, low neutron doses have been reported, especially in active proton therapy. Although the published biological effectiveness of neutrons relative to photons regarding cancer induction is higher, it is unlikely that the neutron dose has a large impact on the second cancer risk of proton therapy patients.

ANALYTICAL MODEL TO ESTIMATE EQUIVALENT DOSE FROM INTERNAL NEUTRONS IN PROTON THERAPY OF CHILDREN WITH INTRACRANIAL TUMORS

This study developed a computationally efficient and easy-to-implement analytical model to estimate the equivalent dose from secondary neutrons originating in the bodies ('internal neutrons') of children receiving intracranial proton radiotherapy. A two-term double-Gaussian mathematical model was fit to previously published internal neutron equivalent dose per therapeutic absorbed dose versus distance from the field edge calculated using Monte Carlo simulations. The model was trained using three intracranial proton fields of a 9-year-old girl. The resulting model was tested against two intracranial fields of a 10-year-old boy by comparing the mean doses in organs at risk of a radiogenic cancer estimated by the model versus those previously calculated by Monte Carlo. On average, the model reproduced the internal neutron organ doses in the 10-year-old boy within 13.5% of the Monte Carlo at 3-10 cm from the field edge and within a factor of 2 of the Monte Carlo at 10-20 cm from the field edge. Beyond 20 cm, the model poorly estimated H/DRx, however, the values were very small, at <0.03 mSv Gy-1.