Complete patient exposure during paediatric brain cancer treatment for photon and proton therapy techniques including imaging procedures.

Background: In radiotherapy, especially when treating children, minimising exposure of healthy tissue can prevent the development of adverse outcomes, including second cancers. In this study we propose a validated Monte Carlo framework to evaluate the complete patient exposure during paediatric brain cancer treatment. Materials and methods: Organ doses were calculated for treatment of a diffuse midline glioma (50.4 Gy with 1.8 Gy per fraction) on a 5-year-old anthropomorphic phantom with 3D-conformal radiotherapy, intensity modulated radiotherapy (IMRT), volumetric modulated arc therapy (VMAT) and intensity modulated pencil beam scanning (PBS) proton therapy. Doses from computed tomography (CT) for planning and on-board imaging for positioning (kV-cone beam CT and X-ray imaging) accounted for the estimate of the exposure of the patient including imaging therapeutic dose. For dose calculations we used validated Monte Carlo-based tools (PRIMO, TOPAS, PENELOPE), while lifetime attributable risk (LAR) was estimated from dose-response relationships for cancer induction, proposed by Schneider et al. Results: Out-of-field organ dose equivalent data of proton therapy are lower, with doses between 0.6 mSv (testes) and 120 mSv (thyroid), when compared to photon therapy revealing the highest out-of-field doses for IMRT ranging between 43 mSv (testes) and 575 mSv (thyroid). Dose delivered by CT ranged between 0.01 mSv (testes) and 72 mSv (scapula) while a single imaging positioning ranged between 2 µSv (testes) and 1.3 mSv (thyroid) for CBCT and 0.03 µSv (testes) and 48 µSv (scapula) for X-ray. Adding imaging dose from CT and daily CBCT to the therapeutic demonstrated an important contribution of imaging to the overall radiation burden in the course of treatment, which is subsequently used to predict the LAR, for selected organs. Conclusion: The complete patient exposure during paediatric brain cancer treatment was estimated by combining the results from different Monte Carlo-based dosimetry tools, showing that proton therapy allows significant reduction of the out-of-field doses and secondary cancer risk in selected organs.

Experimental Validation of an Analytical Program and a Monte Carlo Simulation for the Computation of the Far Out-of-Field Dose in External Beam Photon Therapy Applied to Pediatric Patients.

Background: The out-of-the-field absorbed dose affects the probability of primary second radiation-induced cancers. This is particularly relevant in the case of pediatric treatments. There are currently no methods employed in the clinical routine for the computation of dose distributions from stray radiation in radiotherapy. To overcome this limitation in the framework of conventional teletherapy with photon beams, two computational tools have been developed-one based on an analytical approach and another depending on a fast Monte Carlo algorithm. The purpose of this work is to evaluate the accuracy of these approaches by comparison with experimental data obtained from anthropomorphic phantom irradiations. Materials and Methods: An anthropomorphic phantom representing a 5-year-old child (ATOM, CIRS) was irradiated considering a brain tumor using a Varian TrueBeam linac. Two treatments for the same planned target volume (PTV) were considered, namely, intensity-modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT). In all cases, the irradiation was conducted with a 6-MV energy beam using the flattening filter for a prescribed dose of 3.6 Gy to the PTV. The phantom had natLiF : Mg, Cu, P (MCP-N) thermoluminescent dosimeters (TLDs) in its 180 holes. The uncertainty of the experimental data was around 20%, which was mostly attributed to the MCP-N energy dependence. To calculate the out-of-field dose, an analytical algorithm was implemented to be run from a Varian Eclipse TPS. This algorithm considers that all anatomical structures are filled with water, with the exception of the lungs which are made of air. The fast Monte Carlo code dose planning method was also used for computing the out-of-field dose. It was executed from the dose verification system PRIMO using a phase-space file containing 3x109 histories, reaching an average standard statistical uncertainty of less than 0.2% (coverage factor k = 1 ) on all voxels scoring more than 50% of the maximum dose. The standard statistical uncertainty of out-of-field voxels in the Monte Carlo simulation did not exceed 5%. For the Monte Carlo simulation the actual chemical composition of the materials used in ATOM, as provided by the manufacturer, was employed. Results: In the out-of-the-field region, the absorbed dose was on average four orders of magnitude lower than the dose at the PTV. For the two modalities employed, the discrepancy between the central values of the TLDs located in the out-of-the-field region and the corresponding positions in the analytic model were in general less than 40%. The discrepancy in the lung doses was more pronounced for IMRT. The same comparison between the experimental and the Monte Carlo data yielded differences which are, in general, smaller than 20%. It was observed that the VMAT irradiation produces the smallest out-of-the-field dose when compared to IMRT. Conclusions: The proposed computational methods for the routine calculation of the out-of-the-field dose produce results that are similar, in most cases, with the experimental data. It has been experimentally found that the VMAT irradiation produces the smallest out-of-the-field dose when compared to IMRT for a given PTV.

Validation of a Monte Carlo Framework for Out-of-Field Dose Calculations in Proton Therapy.

Proton therapy enables to deliver highly conformed dose distributions owing to the characteristic Bragg peak and the finite range of protons. However, during proton therapy, secondary neutrons are created, which can travel long distances and deposit dose in out-of-field volumes. This out-of-field absorbed dose needs to be considered for radiation-induced secondary cancers, which are particularly relevant in the case of pediatric treatments. Unfortunately, no method exists in clinics for the computation of the out-of-field dose distributions in proton therapy. To help overcome this limitation, a computational tool has been developed based on the Monte Carlo code TOPAS. The purpose of this work is to evaluate the accuracy of this tool in comparison to experimental data obtained from an anthropomorphic phantom irradiation. An anthropomorphic phantom of a 5-year-old child (ATOM, CIRS) was irradiated for a brain tumor treatment in an IBA Proteus Plus facility using a pencil beam dedicated nozzle. The treatment consisted of three pencil beam scanning fields employing a lucite range shifter. Proton energies ranged from 100 to 165 MeV. A median dose of 50.4 Gy(RBE) with 1.8 Gy(RBE) per fraction was prescribed to the initial planning target volume (PTV), which was located in the cerebellum. Thermoluminescent detectors (TLDs), namely, Li-7-enriched LiF : Mg, Ti (MTS-7) type, were used to detect gamma radiation, which is produced by nuclear reactions, and secondary as well as recoil protons created out-of-field by secondary neutrons. Li-6-enriched LiF : Mg,Cu,P (MCP-6) was combined with Li-7-enriched MCP-7 to measure thermal neutrons. TLDs were calibrated in Co-60 and reported on absorbed dose in water per target dose (µGy/Gy) as well as thermal neutron dose equivalent per target dose (µSv/Gy). Additionally, bubble detectors for personal neutron dosimetry (BD-PND) were used for measuring neutrons (>50 keV), which were calibrated in a Cf-252 neutron beam to report on neutron dose equivalent dose data. The Monte Carlo code TOPAS (version 3.6) was run using a phase-space file containing 1010 histories reaching an average standard statistical uncertainty of less than 0.2% (coverage factor k = 1) on all voxels scoring more than 50% of the maximum dose. The primary beam was modeled following a Fermi-Eyges description of the spot envelope fitted to measurements. For the Monte Carlo simulation, the chemical composition of the tissues represented in ATOM was employed. The dose was tallied as dose-to-water, and data were normalized to the target dose (physical dose) to report on absorbed doses per target dose (mSv/Gy) or neutron dose equivalent per target dose (µSv/Gy), while also an estimate of the total organ dose was provided for a target dose of 50.4 Gy(RBE). Out-of-field doses showed absorbed doses that were 5 to 6 orders of magnitude lower than the target dose. The discrepancy between TLD data and the corresponding scored values in the Monte Carlo calculations involving proton and gamma contributions was on average 18%. The comparison between the neutron equivalent doses between the Monte Carlo simulation and the measured neutron doses was on average 8%. Organ dose calculations revealed the highest dose for the thyroid, which was 120 mSv, while other organ doses ranged from 18 mSv in the lungs to 0.6 mSv in the testes. The proposed computational method for routine calculation of the out-of-the-field dose in proton therapy produces results that are compatible with the experimental data and allow to calculate out-of-field organ doses during proton therapy.