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 cm3phantom 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.

A physics-based analytical model of absorbed dose from primary, leakage, and scattered photons from megavoltage radiotherapy with MLCs.

A burgeoning population of cancer survivors is at risk of late health effects following radiation therapy including second cancers, with the majority of these cancers occurring outside of the treatment volume of the primary cancer. Commercial radiotherapy treatment planning systems underestimate the out-of-field dose. Previous analytical models of out-of-field dose have assumed radial symmetry and/or approximated the dimensions of collimators as semi-infinite planes. The purpose of this work was to develop a physics-based analytical model of total absorbed dose from primary, scattered, and leakage radiation for square fields from a 6 MV beam at any arbitrary point in a phantom, including in-plane, cross-plane, and out-of-plane locations. The model includes the essential physics of radiation transport through beam-limiting-devices including rounded edges of MLC leaves. The model agreed well with measurements and Monte Carlo simulations of absorbed dose in a water-box phantom and was validated for field-sizes ranging from 2 [Formula: see text] 2 to 20 [Formula: see text] 20 cm2. The signed and unsigned average percent differences were [Formula: see text] and 15.9%, respectively, for all points and field-sizes considered. An extended gamma index analysis reveals a 92% pass rate with criteria of 3 mm distance-to-agreement, 3% relative dose difference in-field, and 3 mGy Gy-1 absolute dose difference out-of-field. The average wall-clock time to calculate dose to one million points was 3.3 min. These results suggest that it is feasible to calculate absorbed dose from both therapeutic and stray radiation using physics-based analytical models with good accuracy, thus overcoming a major obstacle to routinely assessing exposures. Additionally, this work demonstrates the importance of relaxing certain simplifications such as assuming a radially symmetric stray-dose distribution. Because the model is physics-based and may be reconfigured according to the dimensions of beam-limiting-devices, adapting it to other treatment units should be straight forward. Uses for such a model include clinical and research applications, such as clinical trials and epidemiological studies.

An objective method to evaluate radiation dose distributions varying by three orders of magnitude.

PURPOSE: Modern radiotherapy practices typically report the absorbed dose (D) within the 5% relative isodose volume (i.e., the therapeutic dose region) to an accuracy of 3%-5%. Gamma-index analysis, the most commonly used method to evaluate dosimetric accuracy, has low sensitivity to discrepancies that occur outside of this region. The objective of this study was to develop an evaluation method with high sensitivity across dose distributions spanning three orders of magnitude. METHODS: We generalized the gamma index to include an additional criterion, the absolute absorbed dose difference, specifically for the low-dose region (i.e., D ≤ 5%). We also proposed a method to objectively select the appropriate magnitudes for relative-dose-difference, absolute-dose-difference, and distance-to-agreement criteria. We demonstrated the generalized gamma-index method by first finding the appropriate generalized gamma-index agreement criteria at an interval of specified passing rates. Next, we used the generalized gamma index to evaluate one-, two-, and three-dimensional absorbed dose distributions in a water-box phantom and voxelized patient geometry. RESULTS: Generalized gamma-index passing rates for one-, two-, and three-dimensional dose distributions were 55.4%, 44.5%, and 8.9%, respectively. Traditional gamma-index passing rates were 100%, 97.8%, and 96.4%, respectively. These results reveal that the generalized method has adequate sensitivity in all regions (i.e., therapeutic and low dose). Additionally, the algorithmic determination of triplets of agreement criteria revealed that they are strong functions of the specified passing rate. CONCLUSIONS: The major finding of this work is that the proposed method provides an objective evaluation of the agreement of dose distributions spanning three orders of magnitude. In particular, this generalized method correctly characterized dosimetric agreement in the low-dose region, which was not possible by traditional methods. The proposed algorithmic selection of agreement criteria decreased subjectivity and requirements of user judgment and skill. This method could find utility in a variety of applications including dose-algorithm development and translation.

The Medical Physics Workforce.

The medical physics workforce comprises approximately 24,000 workers worldwide and approximately 8,200 in the United States. The occupation is a recognized, established, and mature profession that is undergoing considerable growth and change, with many of these changes being driven by scientific, technical, and medical advances. Presently, the medical physics workforce is adequate to meet societal needs. However, data are emerging that suggest potential risks of shortages and other problems that could develop within a few years. Some of the governing factors are well established, such as the increasing number of incident cancers thereby increasing workload, while others, such as the future use of radiation treatments and changes in healthcare economic policies, are uncertain and make the future status of the workforce difficult to forecast beyond the next several years. This review examines some of the major factors that govern supply and demand for medical physicists, discusses published projections and their uncertainties, and presents other information that may help to inform short- and long-term planning of various aspects of the future workforce. It includes a description of the general characteristics of the workforce, including information on its size, educational attainment, certification, age distribution, etc. Because the supply of new workers is governed by educational and training pathways, graduate education, post-doctoral training, and residency training are reviewed, along with trends in state and federal support for research and education. Selected professional aspects of the field also are considered, including professional certification and compensation. We speculate on the future outlook of the workforce and provide recommendations regarding future actions pertaining to the future medical physics workforce.

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.

ADVANTAGES OF MCNPX-BASED LATTICE TALLY OVER MESH TALLY IN HIGH-SPEED MONTE CARLO DOSE RECONSTRUCTION FOR PROTON RADIOTHERAPY.

Monte Carlo simulations are increasingly used to reconstruct dose distributions in radiotherapy research studies. Many studies have used the MCNPX Monte Carlo code with a mesh tally for dose reconstructions. However, when the number of voxels in the simulated patient anatomy is large, the computation time for a mesh tally can become prohibitively long. The purpose of this work was to test the feasibility of using lattice tally instead of mesh tally for whole-body dose reconstructions. We did this by comparing the dosimetric accuracy and computation time of lattice tallies with those of mesh tallies for craniospinal proton irradiation. The two tally methods generated nearly identical dosimetric results, within 1% in dose and within 1 mm distance-to-agreement for 99% of the voxels. For a typical craniospinal proton treatment field, simulation speed was 4 to 17 times faster using the lattice tally than using the mesh tally, depending on the numbers of proton histories and voxels. We conclude that the lattice tally is an acceptable substitute for the mesh tally in dose reconstruction, making it a suitable potential candidate for clinical treatment planning.

Impact of margin size on the predicted risk of radiogenic second cancers following proton arc therapy and volumetric modulated arc therapy for prostate cancer.

We previously determined that the predicted risk of radiogenic second cancer in the bladder and rectum after proton arc therapy (PAT) was less than or equal to that after volumetric modulated arc therapy (VMAT) with photons, but we did not consider the impact of margin size on that risk. The current study was thus conducted to evaluate margin size's effect on the predicted risks of second cancer for the two modalities and the relative risk between them. Seven treatment plans with margins ranging from 0 mm in all directions to 6 mm posteriorly and 8 mm in all other directions were considered for both modalities. We performed risk analyses using three risk models with varying amounts of cell sterilization and calculated ratios of risk for the corresponding PAT and VMAT plans. We found that the change in risk with margin size depended on the risk model but that the relative risk remained nearly constant with margin size, regardless of the amount of cell sterilization modeled. We conclude that while margin size influences the predicted risk of a second cancer for a given modality, it appears to affect both modalities in roughly equal proportions so that the relative risk between PAT and VMAT is approximately equivalent.

Risk of radiogenic second cancers following volumetric modulated arc therapy and proton arc therapy for prostate cancer.

Prostate cancer patients who undergo radiotherapy are at an increased risk to develop a radiogenic second cancer. Proton therapy has been shown to reduce the predicted risk of second cancer when compared to intensity modulated radiotherapy. However, it is unknown if this is also true for the rotational therapies proton arc therapy and volumetric modulated arc therapy (VMAT). The objective of this study was to compare the predicted risk of cancer following proton arc therapy and VMAT for prostate cancer. Proton arc therapy and VMAT plans were created for three patients. Various risk models were combined with the dosimetric data (therapeutic and stray dose) to predict the excess relative risk (ERR) of cancer in the bladder and rectum. Ratios of ERR values (RRR) from proton arc therapy and VMAT were calculated. RRR values ranged from 0.74 to 0.99, and all RRR values were shown to be statistically less than 1, except for the value calculated with the linear-non-threshold risk model. We conclude that the predicted risk of cancer in the bladder or rectum following proton arc therapy for prostate cancer is either less than or approximately equal to the risk following VMAT, depending on which risk model is applied.

The future workforce in cancer prevention: advancing discovery, research, and technology.

As part of a 2-day conference on October 15 and 16, 2009, a nine-member task force composed of scientists, clinicians, educators, administrators, and students from across the USA was formed to discuss research, discovery, and technology obstacles to progress in cancer prevention and control, specifically those related to the cancer prevention workforce. This article summarizes the task force's findings on the current state of the cancer prevention workforce in this area and its needs for the future. The task force identified two types of barriers impeding the current cancer prevention workforce in research, discovery, and technology from reaching its fullest potential: (1) limited cross-disciplinary research opportunities with underutilization of some disciplines is hampering discovery and research in cancer prevention, and (2) new research avenues are not being investigated because technology development and implementation are lagging. Examples of impediments and desired outcomes are provided in each of these areas. Recommended solutions to these problems are based on the goals of enhancing the current cancer prevention workforce and accelerating the pace of discovery and clinical translation.

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.

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.

Development and verification of an analytical algorithm to predict absorbed dose distributions in ocular proton therapy using Monte Carlo simulations.

Proton beam radiotherapy is an effective and non-invasive treatment for uveal melanoma. Recent research efforts have focused on improving the dosimetric accuracy of treatment planning and overcoming the present limitation of relative analytical dose calculations. Monte Carlo algorithms have been shown to accurately predict dose per monitor unit (D/MU) values, but this has yet to be shown for analytical algorithms dedicated to ocular proton therapy, which are typically less computationally expensive than Monte Carlo algorithms. The objective of this study was to determine if an analytical method could predict absolute dose distributions and D/MU values for a variety of treatment fields like those used in ocular proton therapy. To accomplish this objective, we used a previously validated Monte Carlo model of an ocular nozzle to develop an analytical algorithm to predict three-dimensional distributions of D/MU values from pristine Bragg peaks and therapeutically useful spread-out Bragg peaks (SOBPs). Results demonstrated generally good agreement between the analytical and Monte Carlo absolute dose calculations. While agreement in the proximal region decreased for beams with less penetrating Bragg peaks compared with the open-beam condition, the difference was shown to be largely attributable to edge-scattered protons. A method for including this effect in any future analytical algorithm was proposed. Comparisons of D/MU values showed typical agreement to within 0.5%. We conclude that analytical algorithms can be employed to accurately predict absolute proton dose distributions delivered by an ocular nozzle.

Risk of secondary malignant neoplasms from proton therapy and intensity-modulated x-ray therapy for early-stage prostate cancer.

PURPOSE: To assess the risk of a secondary malignant neoplasm (SMN) from proton therapy relative to intensity-modulated radiation therapy (IMRT) using X-rays, taking into account contributions from both primary and secondary sources of radiation, for prostate cancer. METHODS AND MATERIALS: A proton therapy plan and a 6-MV IMRT plan were constructed for 3 patients with early-stage adenocarcinoma of the prostate. Doses from the primary fields delivered to organs at risk of developing an SMN were determined from treatment plans. Secondary doses from the proton therapy and IMRT were determined from Monte Carlo simulations and available measured data, respectively. The risk of an SMN was estimated from primary and secondary doses on an organ-by-organ basis by use of risk models from the Committee on the Biological Effects of Ionizing Radiation. RESULTS: Proton therapy reduced the risk of an SMN by 26% to 39% compared with IMRT. The risk of an SMN for both modalities was greatest in the in-field organs. However, the risks from the in-field organs were considerably lower with the proton therapy plan than with the IMRT plan. This reduction was attributed to the substantial sparing of the rectum and bladder from exposure to the therapeutic beam by the proton therapy plan. CONCLUSIONS: When considering exposure to primary and secondary radiation, proton therapy can reduce the risk of an SMN in prostate patients compared with contemporary IMRT.

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.