Comparison of whole-body phantom designs to estimate organ equivalent neutron doses for secondary cancer risk assessment in proton therapy

Secondary neutron fluence created during proton therapy can be a significant source of radiation exposure in organs distant from the treatment site, especially in pediatric patients. Various published studies have used computational phantoms to estimate neutron equivalent doses in proton therapy. In these simulations, whole-body patient representations were applied considering either generic whole-body phantoms or generic age- and gender-dependent phantoms. No studies to date have reported using patient-specific geometry information. The purpose of this study was to estimate the effects of patient–phantom matching when using computational pediatric phantoms. To achieve this goal, three sets of phantoms, including different ages and genders, were compared to the patients' whole-body CT. These sets consisted of pediatric age-specific reference, age-adjusted reference and anatomically sculpted phantoms. The neutron equivalent dose for a subset of out-of-field organs was calculated using the GEANT4 Monte Carlo toolkit, where proton fields were used to irradiate the cranium and the spine of all phantoms and the CT-segmented patient models. The maximum neutron equivalent dose per treatment absorbed dose was calculated and found to be on the order of 0 to 5 mSv Gy−1. The relative dose difference between each phantom and their respective CT-segmented patient model for most organs showed a dependence on how close the phantom and patient heights were matched. The weight matching was found to have much smaller impact on the dose accuracy except for very heavy patients. Analysis of relative dose difference with respect to height difference suggested that phantom sculpting has a positive effect in terms of dose accuracy as long as the patient is close to the 50th percentile height and weight. Otherwise, the benefit of sculpting was masked by inherent uncertainties, i.e. variations in organ shapes, sizes and locations. Other sources of uncertainty included errors associated with beam positioning, neutron weighting factor definition and organ segmentation. This work demonstrated the importance of hybrid phantom height matching for more accurate organ dose calculation in proton therapy and the potential limitations of reference phantoms released by regulatory bodies for radiation therapy applications. (paper)