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.

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

The purpose of this study was to independently apply an analytical model for equivalent dose from neutrons produced in a passive-scattering proton therapy treatment unit, H. To accomplish this objective, we applied the previously-published model to treatment plans of two pediatric patients. Their model accounted for neutrons generated by mono-energetic proton beams stopping in a closed aperture. To implement their model to a clinical setting, we adjusted it to account for the area of a collimating aperture, energy modulation, air gap between the treatment unit and patient, and radiation weighting factor. We used the adjusted model to estimate H per prescribed proton absorbed dose, D Rx , for the passive-scattering proton therapy beams of two children, a 9-year-old girl and 10-year-old boy, who each received intracranial boost fields as part of their treatment. In organs and tissues at risk for radiation-induced subsequent malignant neoplasms, T, we calculated the mass-averaged H, H T , per D Rx . Finally, we compared H T /D Rx values to those of previously-published Monte Carlo (MC) simulations of these patients' fields. H T /D Rx values of the adjusted model deviated from the MC result for each organ on average by 20.8 ± 10.0% and 44.2 ± 17.6% for the girl and boy, respectively. The adjusted model underestimated the MC result in all T of each patient, with the exception of the girl's bladder, for which the adjusted model overestimated H T /D Rx by 3.1%. The adjusted model provided a better estimate of H T /D Rx than the unadjusted model. That is, between the two models, the adjusted model reduced the deviation from the MC result by approximately 37.0% and 46.7% for the girl and boy, respectively. We found that the previously-published analytical model, combined with adjustment factors to enhance its clinical applicability, predicted H T /D Rx in out-of-field organs and tissues at risk for subsequent malignant neoplasms with acceptable accuracy. This independent application demonstrated that the analytical model may be useful broadly for clinicians and researchers to calculate equivalent dose from neutrons produced externally to the patient in passive-scattering proton therapy.

Supplemental computational phantoms to estimate out-of-field absorbed dose in photon radiotherapy.

The purpose of this study was to develop a straightforward method of supplementing patient anatomy and estimating out-of-field absorbed dose for a cohort of pediatric radiotherapy patients with limited recorded anatomy. A cohort of nine children, aged 2-14 years, who received 3D conformal radiotherapy for low-grade localized brain tumors (LBTs), were randomly selected for this study. The extent of these patients' computed tomography simulation image sets were cranial only. To approximate their missing anatomy, we supplemented the LBT patients' image sets with computed tomography images of patients in a previous study with larger extents of matched sex, height, and mass and for whom contours of organs at risk for radiogenic cancer had already been delineated. Rigid fusion was performed between the LBT patients' data and that of the supplemental computational phantoms using commercial software and in-house codes. In-field dose was calculated with a clinically commissioned treatment planning system, and out-of-field dose was estimated with a previously developed analytical model that was re-fit with parameters based on new measurements for intracranial radiotherapy. Mean doses greater than 1 Gy were found in the red bone marrow, remainder, thyroid, and skin of the patients in this study. Mean organ doses between 150 mGy and 1 Gy were observed in the breast tissue of the girls and lungs of all patients. Distant organs, i.e. prostate, bladder, uterus, and colon, received mean organ doses less than 150 mGy. The mean organ doses of the younger, smaller LBT patients (0-4 years old) were a factor of 2.4 greater than those of the older, larger patients (8-12 years old). Our findings demonstrated the feasibility of a straightforward method of applying supplemental computational phantoms and dose-calculation models to estimate absorbed dose for a set of children of various ages who received radiotherapy and for whom anatomies were largely missing in their original computed tomography simulations.

Reducing the cost of proton radiation therapy: the feasibility of a streamlined treatment technique for prostate cancer.

Proton radiation therapy is an effective modality for cancer treatments, but the cost of proton therapy is much higher compared to conventional radiotherapy and this presents a formidable barrier to most clinical practices that wish to offer proton therapy. Little attention in literature has been paid to the costs associated with collimators, range compensators and hypofractionation. The objective of this study was to evaluate the feasibility of cost-saving modifications to the present standard of care for proton treatments for prostate cancer. In particular, we quantified the dosimetric impact of a treatment technique in which custom fabricated collimators were replaced with a multileaf collimator (MLC) and the custom range compensators (RC) were eliminated. The dosimetric impacts of these modifications were assessed for 10 patients with a commercial treatment planning system (TPS) and confirmed with corresponding Monte Carlo simulations. We assessed the impact on lifetime risks of radiogenic second cancers using detailed dose reconstructions and predictive dose-risk models based on epidemiologic data. We also performed illustrative calculations, using an isoeffect model, to examine the potential for hypofractionation. Specifically, we bracketed plausible intervals of proton fraction size and total treatment dose that were equivalent to a conventional photon treatment of 79.2 Gy in 44 fractions. Our results revealed that eliminating the RC and using an MLC had negligible effect on predicted dose distributions and second cancer risks. Even modest hypofractionation strategies can yield substantial cost savings. Together, our results suggest that it is feasible to modify the standard of care to increase treatment efficiency, reduce treatment costs to patients and insurers, while preserving high treatment quality.

A comparative study on the risks of radiogenic second cancers and cardiac mortality in a set of pediatric medulloblastoma patients treated with photon or proton craniospinal irradiation.

PURPOSE: To compare the risks of radiogenic second cancers and cardiac mortality in 17 pediatric medulloblastoma patients treated with passively scattered proton or field-in-field photon craniospinal irradiation (CSI). MATERIAL/METHODS: Standard of care photon or proton CSI treatment plans were created for all 17 patients in a commercial treatment planning system (TPS) (Eclipse version 8.9; Varian Medical Systems, Palo Alto, CA) and prescription dose was 23.4 or 23.4 Gy (RBE) to the age specific target volume at 1.8 Gy/fraction. The therapeutic doses from proton and photon CSI plans were estimated from TPS. Stray radiation doses were determined from Monte Carlo simulations for proton CSI and from measurements and TPS for photon CSI. The Biological Effects of Ionization Radiation VII report and a linear model based on childhood cancer survivor data were used for risk predictions of second cancer and cardiac mortality, respectively. RESULTS: The ratios of lifetime attributable risk (RLARs) (proton/photon) ranged from 0.10 to 0.22 for second cancer incidence and ranged from 0.20 to 0.53 for second cancer mortality, respectively. The ratio of relative risk (RRR) (proton/photon) of cardiac mortality ranged from 0.12 to 0.24. The RLARs of both cancer incidence and mortality decreased with patient's age at exposure (e), while the RRRs of cardiac mortality increased with e. Girls had a significantly higher RLAR of cancer mortality than boys. CONCLUSION: Passively scattered proton CSI provides superior predicted outcomes by conferring lower predicted risks of second cancer and cardiac mortality than field-in-field photon CSI for all medulloblastoma patients in a large clinically representative sample in the United States, but the magnitude of superiority depends strongly on the patients' anatomical development status.

Comparison of risk of radiogenic second cancer following photon and proton craniospinal irradiation for a pediatric medulloblastoma patient.

Pediatric patients who received radiation therapy are at risk of developing side effects such as radiogenic second cancer. We compared proton and photon therapies in terms of the predicted risk of second cancers for a 4 year old medulloblastoma patient receiving craniospinal irradiation (CSI). Two CSI treatment plans with 23.4 Gy or Gy (RBE) prescribed dose were computed: a three-field 6 MV photon therapy plan and a four-field proton therapy plan. The primary doses for both plans were determined using a commercial treatment planning system. Stray radiation doses for proton therapy were determined from Monte Carlo simulations, and stray radiation doses for photon therapy were determined from measured data. Dose-risk models based on the Biological Effects of Ionization Radiation VII report were used to estimate the risk of second cancer in eight tissues/organs. Baseline predictions of the relative risk for each organ were always less for proton CSI than for photon CSI at all attained ages. The total lifetime attributable risk of the incidence of second cancer considered after proton CSI was much lower than that after photon CSI, and the ratio of lifetime risk was 0.18. Uncertainty analysis revealed that the qualitative findings of this study were insensitive to any plausible changes of dose-risk models and mean radiation weighting factor for neutrons. Proton therapy confers lower predicted risk of second cancer than photon therapy for the pediatric medulloblastoma patient.

External-beam accelerated partial breast irradiation using multiple proton beam configurations.

PURPOSE: To explore multiple proton beam configurations for optimizing dosimetry and minimizing uncertainties for accelerated partial breast irradiation (APBI) and to compare the dosimetry of proton with that of photon radiotherapy for treatment of the same clinical volumes. METHODS AND MATERIALS: Proton treatment plans were created for 11 sequential patients treated with three-dimensional radiotherapy (3DCRT) photon APBI using passive scattering proton beams (PSPB) and were compared with clinically treated 3DCRT photon plans. Monte Carlo calculations were used to verify the accuracy of the proton dose calculation from the treatment planning system. The impact of range, motion, and setup uncertainty was evaluated with tangential vs. en face beams. RESULTS: Compared with 3DCRT photons, the absolute reduction of the mean of V100 (the volume receiving 100% of prescription dose), V90, V75, V50, and V20 for normal breast using protons are 3.4%, 8.6%, 11.8%, 17.9%, and 23.6%, respectively. For breast skin, with the similar V90 as 3DCRT photons, the proton plan significantly reduced V75, V50, V30, and V10. The proton plan also significantly reduced the dose to the lung and heart. Dose distributions from Monte Carlo simulations demonstrated minimal deviation from the treatment planning system. The tangential beam configuration showed significantly less dose fluctuation in the chest wall region but was more vulnerable to respiratory motion than that for the en face beams. Worst-case analysis demonstrated the robustness of designed proton beams with range and patient setup uncertainties. CONCLUSIONS: APBI using multiple proton beams spares significantly more normal tissue, including nontarget breast and breast skin, than 3DCRT using photons. It is robust, considering the range and patient setup uncertainties.

Predicted risks of second malignant neoplasm incidence and mortality due to secondary neutrons in a girl and boy receiving proton craniospinal irradiation.

The purpose of this study was to compare the predicted risks of second malignant neoplasm (SMN) incidence and mortality from secondary neutrons for a 9-year-old girl and a 10-year-old boy who received proton craniospinal irradiation (CSI). SMN incidence and mortality from neutrons were predicted from equivalent doses to radiosensitive organs for cranial, spinal and intracranial boost fields. Therapeutic proton absorbed dose and equivalent dose from neutrons were calculated using Monte Carlo simulations. Risks of SMN incidence and mortality in most organs and tissues were predicted by applying risks models from the National Research Council of the National Academies to the equivalent dose from neutrons; for non-melanoma skin cancer, risk models from the International Commission on Radiological Protection were applied. The lifetime absolute risks of SMN incidence due to neutrons were 14.8% and 8.5%, for the girl and boy, respectively. The risks of a fatal SMN were 5.3% and 3.4% for the girl and boy, respectively. The girl had a greater risk for any SMN except colon and liver cancers, indicating that the girl's higher risks were not attributable solely to greater susceptibility to breast cancer. Lung cancer predominated the risk of SMN mortality for both patients. This study suggests that the risks of SMN incidence and mortality from neutrons may be greater for girls than for boys treated with proton CSI.

A GPU implementation of a track-repeating algorithm for proton radiotherapy dose calculations.

An essential component in proton radiotherapy is the algorithm to calculate the radiation dose to be delivered to the patient. The most common dose algorithms are fast but they are approximate analytical approaches. However their level of accuracy is not always satisfactory, especially for heterogeneous anatomical areas, like the thorax. Monte Carlo techniques provide superior accuracy; however, they often require large computation resources, which render them impractical for routine clinical use. Track-repeating algorithms, for example the fast dose calculator, have shown promise for achieving the accuracy of Monte Carlo simulations for proton radiotherapy dose calculations in a fraction of the computation time. We report on the implementation of the fast dose calculator for proton radiotherapy on a card equipped with graphics processor units (GPUs) rather than on a central processing unit architecture. This implementation reproduces the full Monte Carlo and CPU-based track-repeating dose calculations within 2%, while achieving a statistical uncertainty of 2% in less than 1 min utilizing one single GPU card, which should allow real-time accurate dose calculations.

Adjustment of the lateral and longitudinal size of scanned proton beam spots using a pre-absorber to optimize penumbrae and delivery efficiency.

In scanned-beam proton therapy, the beam spot properties, such as the lateral and longitudinal size and the minimum achievable range, are influenced by beam optics, scattering media and drift spaces in the treatment unit. Currently available spot scanning systems offer few options for adjusting these properties. We investigated a method for adjusting the lateral and longitudinal spot size that utilizes downstream plastic pre-absorbers located near a water phantom. The spot size adjustment was characterized using Monte Carlo simulations of a modified commercial scanned-beam treatment head. Our results revealed that the pre-absorbers can be used to reduce the lateral full width at half maximum (FWHM) of dose spots in water by up to 14 mm, and to increase the longitudinal extent from about 1 mm to 5 mm at residual ranges of 4 cm and less. A large factor in manipulating the lateral spot sizes is the drift space between the pre-absorber and the water phantom. Increasing the drift space from 0 cm to 15 cm leads to an increase in the lateral FWHM from 2.15 cm to 2.87 cm, at a water-equivalent depth of 1 cm. These findings suggest that this spot adjustment method may improve the quality of spot-scanned proton treatments.

Assessment of targeting accuracy of a low-energy stereotactic radiosurgery treatment for age-related macular degeneration.

Age-related macular degeneration (AMD), a leading cause of blindness in the United States, is a neovascular disease that may be controlled with radiation therapy. Early patient outcomes of external beam radiotherapy, however, have been mixed. Recently, a novel multimodality treatment was developed, comprising external beam radiotherapy and concomitant treatment with a vascular endothelial growth factor inhibitor. The radiotherapy arm is performed by stereotactic radiosurgery, delivering a 16 Gy dose in the macula (clinical target volume, CTV) using three external low-energy x-ray fields while adequately sparing normal tissues. The purpose of our study was to test the sensitivity of the delivery of the prescribed dose in the CTV using this technique and of the adequate sparing of normal tissues to all plausible variations in the position and gaze angle of the eye. Using Monte Carlo simulations of a 16 Gy treatment, we varied the gaze angle by ±5° in the polar and azimuthal directions, the linear displacement of the eye ±1 mm in all orthogonal directions, and observed the union of the three fields on the posterior wall of spheres concentric with the eye that had diameters between 20 and 28 mm. In all cases, the dose in the CTV fluctuated <6%, the maximum dose in the sclera was <20 Gy, the dose in the optic disc, optic nerve, lens and cornea were <0.7 Gy and the three-field junction was adequately preserved. The results of this study provide strong evidence that for plausible variations in the position of the eye during treatment, either by the setup error or intrafraction motion, the prescribed dose will be delivered to the CTV and the dose in structures at risk will be kept far below tolerance doses.

An analytic model of neutron ambient dose equivalent and equivalent dose for proton radiotherapy.

Stray neutrons generated in passively scattered proton therapy are of concern because they increase the risk that a patient will develop a second cancer. Several investigations characterized stray neutrons in proton therapy using experimental measurements and Monte Carlo simulations, but capabilities of analytical methods to predict neutron exposures are less well developed. The goal of this study was to develop a new analytical model to calculate neutron ambient dose equivalent in air and equivalent dose in phantom based on Monte Carlo modeling of a passively scattered proton therapy unit. The accuracy of the new analytical model is superior to a previous analytical model and comparable to the accuracy of typical Monte Carlo simulations and measurements. Predictions from the new analytical model agreed reasonably well with corresponding values predicted by a Monte Carlo code using an anthropomorphic phantom.

The risk of developing a second cancer after receiving craniospinal proton irradiation.

The purpose of this work was to compare the risk of developing a second cancer after craniospinal irradiation using photon versus proton radiotherapy by means of simulation studies designed to account for the effects of neutron exposures. Craniospinal irradiation of a male phantom was calculated for passively-scattered and scanned-beam proton treatment units. Organ doses were estimated from treatment plans; for the proton treatments, the amount of stray radiation was calculated separately using the Monte Carlo method. The organ doses were converted to risk of cancer incidence using a standard formalism developed for radiation protection purposes. The total lifetime risk of second cancer due exclusively to stray radiation was 1.5% for the passively scattered treatment versus 0.8% for the scanned proton beam treatment. Taking into account the therapeutic and stray radiation fields, the risk of second cancer from intensity-modulated radiation therapy and conventional radiotherapy photon treatments were 7 and 12 times higher than the risk associated with scanned-beam proton therapy, respectively, and 6 and 11 times higher than with passively scattered proton therapy, respectively. Simulations revealed that both passively scattered and scanned-beam proton therapies confer significantly lower risks of second cancers than 6 MV conventional and intensity-modulated photon therapies.

Stray radiation dose and second cancer risk for a pediatric patient receiving craniospinal irradiation with proton beams.

Proton beam radiotherapy unavoidably exposes healthy tissue to stray radiation emanating from the treatment unit and secondary radiation produced within the patient. These exposures provide no known benefit and may increase a patient's risk of developing a radiogenic cancer. The aims of this study were to calculate doses to major organs and tissues and to estimate second cancer risk from stray radiation following craniospinal irradiation (CSI) with proton therapy. This was accomplished using detailed Monte Carlo simulations of a passive-scattering proton treatment unit and a voxelized phantom to represent the patient. Equivalent doses, effective dose and corresponding risk for developing a fatal second cancer were calculated for a 10-year-old boy who received proton therapy. The proton treatment comprised CSI at 30.6 Gy plus a boost of 23.4 Gy to the clinical target volume. The predicted effective dose from stray radiation was 418 mSv, of which 344 mSv was from neutrons originating outside the patient; the remaining 74 mSv was caused by neutrons originating within the patient. This effective dose corresponds to an attributable lifetime risk of a fatal second cancer of 3.4%. The equivalent doses that predominated the effective dose from stray radiation were in the lungs, stomach and colon. These results establish a baseline estimate of the stray radiation dose and corresponding risk for a pediatric patient undergoing proton CSI and support the suitability of passively-scattered proton beams for the treatment of central nervous system tumors in pediatric patients.

Monte Carlo fast dose calculator for proton radiotherapy: application to a voxelized geometry representing a patient with prostate cancer.

The Monte Carlo method is used to provide accurate dose estimates in proton radiation therapy research. While it is more accurate than commonly used analytical dose calculations, it is computationally intense. The aim of this work was to characterize for a clinical setup the fast dose calculator (FDC), a Monte Carlo track-repeating algorithm based on GEANT4. FDC was developed to increase computation speed without diminishing dosimetric accuracy. The algorithm used a database of proton trajectories in water to calculate the dose of protons in heterogeneous media. The extrapolation from water to 41 materials was achieved by scaling the proton range and the scattering angles. The scaling parameters were obtained by comparing GEANT4 dose distributions with those calculated with FDC for homogeneous phantoms. The FDC algorithm was tested by comparing dose distributions in a voxelized prostate cancer patient as calculated with well-known Monte Carlo codes (GEANT4 and MCNPX). The track-repeating approach reduced the CPU time required for a complete dose calculation in a voxelized patient anatomy by more than two orders of magnitude, while on average reproducing the results from the Monte Carlo predictions within 2% in terms of dose and within 1 mm in terms of distance.

A TRACK-REPEATING ALGORITHM FOR FAST MONTE CARLO DOSE CALCULATIONS OF PROTON RADIOTHERAPY.

Monte Carlo codes are utilized for accurate dose calculations in proton radiation therapy research. While they are superior in accuracy to commonly used analytical dose calculations, they require significantly longer computation times. The aim of this work is to characterize a Monte Carlo track-repeating algorithm to increase computation speed without compromising dosimetric accuracy. The track-repeating approach reduced the CPU time required for a complete dose calculation in voxelized patient anatomy by more than two orders of magnitude, while on average reproducing the results from the traditional Monte Carlo approach within 4% dose difference and within 1-mm distance to agreement.

Reducing stray radiation dose to patients receiving passively scattered proton radiotherapy for prostate cancer.

Proton beam radiotherapy exposes healthy tissue to stray radiation emanating from the treatment unit and secondary radiation produced within the patient. These exposures provide no known benefit and may increase a patient's risk of developing a radiogenic second cancer. The aim of this study was to explore strategies to reduce stray radiation dose to a patient receiving a 76 Gy proton beam treatment for cancer of the prostate. The whole-body effective dose from stray radiation, E, was estimated using detailed Monte Carlo simulations of a passively scattered proton treatment unit and an anthropomorphic phantom. The predicted value of E was 567 mSv, of which 320 mSv was attributed to leakage from the treatment unit; the remainder arose from scattered radiation that originated within the patient. Modest modifications of the treatment unit reduced E by 212 mSv. Surprisingly, E from a modified passive-scattering device was only slightly higher (109 mSv) than from a nozzle with no leakage, e.g., that which may be approached with a spot-scanning technique. These results add to the body of evidence supporting the suitability of passively scattered proton beams for the treatment of prostate cancer, confirm that the effective dose from stray radiation was not excessive, and, importantly, show that it can be substantially reduced by modest enhancements to the treatment unit.

A comparison of the measured responses of a tissue-equivalent proportional counter to high energy heavy (HZE) particles and those simulated using the Geant4 Monte Carlo code.

Monte Carlo simulations of heavy ion interactions using the Geant4 toolkit were compared with measurements of energy deposition in a spherical tissue-equivalent proportional counter (TEPC). A spherical cavity with a physical diameter of 12.7 mm was filled with propane-based tissue-equivalent gas surrounded by a wall of A-150 tissue-equivalent plastic that was 2.54 mm to thick. Measurements and Monte Carlo simulations were used to record the energy deposition and the trajectory of the incident particle on an event-by-event basis for ions ranging in atomic number from 2 ((4)He) to 26 ((56)Fe) and in energy from 200 MeV/nucleon to 1000 MeV/nucleon. In the simulations, tracking of secondary electrons was terminated when the range of an electron was below a specified threshold. The effects of range cuts for electrons at 0.5 µm, 1 µm, 10 µm, and 100 µm were evaluated. To simulate an energy deposition influenced by large numbers of low energy electrons with large transverse momentum, it was necessary to track electrons down to range cuts of 10 µm or less. The Geant4 simulated data closely matched the measured data acquired using a TEPC for incident particles traversing the center of the detector as well as near the gas-wall interface. Values of frequency mean lineal energy and dose mean lineal energy were within 8% of the measured data. The production of secondary particles in the aluminum vacuum chamber had no effect on the response of the TEPC for (56)Fe at 1000 MeV/nucleon. The results of this study confirm that Geant4 can simulate patterns of energy deposition for existing microdosimeters and is valuable for improving the design of a new generation of detectors used for space dosimetry and for characterizing particle beams used in hadron radiotherapy.