TOPAS Monte Carlo model of MD anderson scanning proton beam for simulation studies in proton therapy

Fetal dose from proton pencil beam scanning craniospinal irradiation during pregnancy: a Monte Carlo study

Objective. We conducted a Monte Carlo study to comprehensively investigate the fetal dose resulting from proton pencil beam scanning (PBS) craniospinal irradiation (CSI) during pregnancy.Approach. The gestational-age dependent pregnant phantom series developed at the University of Florida (UF) were converted into DICOM-RT format (CT images and structures) and imported into a treatment planning system (TPS) (Eclipse v15.6) commissioned to a IBA PBS nozzle. A proton PBS CSI plan (prescribed dose: 36 Gy) was created on the phantoms. The TOPAS MC code was used to simulate the proton PBS CSI on the phantoms, for which MC beam properties at the nozzle exit (spot size, spot divergence, mean energy, and energy spread) were matched to IBA PBS nozzle beam measurement data. We calculated mean absorbed doses for 28 organs and tissues and whole body of the fetus at eight gestational ages (8, 10, 15, 20, 25, 30, 35, and 38 weeks). For contextual purposes, the fetal organ/tissue doses from the treatment planning CT scan of the mother's head and torso were estimated using the National Cancer Institute dosimetry system for CT (NCICT, Version 3) considering a low-dose CT protocol (CTDIvol: 8.97 mGy).Main results. The majority of the fetal organ/tissue doses from the proton PBS CSI treatment fell within a range of 3-6 mGy. The fetal organ/tissue doses for the 38 week phantom showed the largest variation with the doses ranging from 2.9 mGy (adrenals) to 8.2 mGy (eye lenses) while the smallest variation ranging from 3.2 mGy (oesophagus) to 4.4 mGy (brain) was observed for the doses for the 20 week phantom. The fetal whole-body dose ranged from 3.7 mGy (25 weeks) to 5.8 mGy (8 weeks). Most of the fetal doses from the planning CT scan fell within a range of 7-13 mGy, approximately 2-to-9 times lower than the fetal dose equivalents of the proton PBS CSI treatment (assuming a quality factor of 7).Significance. The fetal organ/tissue doses observed in the present work will be useful for one of the first clinically informative predictions on the magnitude of fetal dose during proton PBS CSI during pregnancy.

Monte Carlo methods for device simulations in radiation therapy

Monte Carlo (MC) simulations play an important role in radiotherapy, especially as a method to evaluate physical properties that are either impossible or difficult to measure. For example, MC simulations (MCSs) are used to aid in the design of radiotherapy devices or to understand their properties. The aim of this article is to review the MC method for device simulations in radiation therapy. After a brief history of the MC method and popular codes in medical physics, we review applications of the MC method to model treatment heads for neutral and charged particle radiation therapy as well as specific in-room devices for imaging and therapy purposes. We conclude by discussing the impact that MCSs had in this field and the role of MC in future device design.

Monte Carlo framework for commissioning a synchrotron-based discrete spot scanning proton beam system and treatment plan verification

This study aimed to develop a Monte Carlo (MC) framework for commissioning the narrow proton beams (spot size sigma, 5.2 mm 2 mm at isocenter for 69.4 MeV-221.3 MeV for the main beam option and 4.1 mm 1.3 mm for the minibeam option respectively) of a synchrotron-based proton therapy system and design an independent absolute dose calculation engine for intensity-modulated proton treatments. A proton therapy system (Hitachi PROBEAT-V) was simulated using divergent and convergent beam models at the nozzle entrance. The innovative source weighting scheme for the MC simulation with TOPAS (TOol for PArticle Simulations) was implemented using dose output data for the absolute dose calculations. The results of the MC simulation were compared to the experimental data, analyzed and used to commission the treatment planning system. Two MC models, divergent and convergent beams were implemented. The convergent beam model produced a high level of agreement when MC and measurements were analyzed. The beam ellipticity did not result in significant differences between MC simulated and treatment planning system calculated doses. A model of a synchrotron-based spot scanning proton therapy system has been developed and implemented in the TOPAS MC transport code framework. The dose computation engine is useful for treatment plan verification with primary and minibeam beam option.

A dose voxel kernel method for rapid reconstruction of out-of-field neutron dose of patients in pencil beam scanning (PBS) proton therapy

Monte Carlo (MC) radiation transport methods are used for dose calculation as 'gold standard.' However, the method is computationally time-consuming and thus impractical for normal tissue dose reconstructions for the large number of proton therapy patients required for epidemiologic investigations of late health effects. In the present study, we developed a new dose calculation method for the rapid reconstruction of out-of-field neutron dose to patients undergoing pencil beam scanning (PBS) proton therapy. The new dose calculation method is based on neutron dose voxel kernels (DVKs) generated by MC simulations of a proton pencil beam irradiating a water phantom (60 × 60 × 300 cm³), which was conducted using a MC proton therapy simulation code, TOPAS. The DVKs were generated for 19 beam energies (from 70 to 250 MeV with the 10 MeV interval) and three range shifter thicknesses (1, 3, and 5 cm). An in-house program was written in C++ to superimpose the DVKs onto a patient CT images according to proton beam characteristics (energy, position, and direction) available in treatment plans. The DVK dose calculation method was tested by calculating organ/tissue-specific neutron doses of 1- and 5-year-old whole-body computational phantoms where intracranial and craniospinal irradiations were simulated. The DVK-based doses generally showed reasonable agreement with those calculated by direct MC simulations with a detailed PBS model that were previously published, with differences mostly less than 30% and 10% for the intracranial and craniospinal irradiations, respectively. The computation time of the DVK method for one patient ranged from 1 to 30 min on a single CPU core of a personal computer, demonstrating significant improvement over the direct MC dose calculation requiring several days on high-performance computing servers. Our DVK-based dose calculation method will be useful when dosimetry is needed for the large number of patients such as for epidemiologic or clinical research.

A Monte Carlo model for organ dose reconstruction of patients in pencil beam scanning (PBS) proton therapy for epidemiologic studies of late effects

Significant efforts such as the Pediatric Proton/Photon Consortium Registry (PPCR) involving multiple proton therapy centers have been made to conduct collaborative studies evaluating outcomes following proton therapy. As a groundwork dosimetry effort for the late effect investigation, we developed a Monte Carlo (MC) model of proton pencil beam scanning (PBS) to estimate organ/tissue doses of pediatric patients at the Maryland Proton Treatment Center (MPTC), one of the proton centers involved in the PPCR. The MC beam modeling was performed using the TOPAS (TOol for PArticle Simulation) MC code and commissioned to match measurement data within 1% for range, and 0.3 mm for spot sizes. The established MC model was then tested by calculating organ/tissue doses for sample intracranial and craniospinal irradiations on whole-body pediatric computational human phantoms. The simulated dose distributions were compared with the treatment planning system dose distributions, showing the 3 mm/3% gamma index passing rates of 94%-99%, validating our simulations with the MC model. The calculated organ/tissue doses per prescribed doses for the craniospinal irradiations (1 mGy Gy^-1 to 1 Gy Gy^-1) were generally much higher than those for the intracranial irradiations (2.1 μGy Gy^-1 to 0.1 Gy Gy^-1), which is due to the larger field coverage of the craniospinal irradiations. The largest difference was observed at the adrenal dose, i.e. ∼3000 times. In addition, the calculated organ/tissue doses were compared with those calculated with a simplified MC model, showing that the beam properties (i.e. spot size, spot divergence, mean energy, and energy spread) do not significantly influence dose calculations despite the limited irradiation cases. This implies that the use of the MC model commissioned to the MPTC measurement data might be dosimetrically acceptable for patient dose reconstructions at other proton centers particularly when their measurement data are unavailable. The developed MC model will be used to reconstruct organ/tissue doses for MPTC pediatric patients collected in the PPCR.