Secondary neutron dose contribution from pencil beam scanning, scattered and spatially fractionated proton therapy

Secondary neutron dosimetry for conformal FLASH proton therapy

BACKGROUND

Cyclotron-based proton therapy systems utilize the highest proton energies to achieve an ultra-high dose rate (UHDR) for FLASH radiotherapy. The deep-penetrating range associated with this high energy can be modulated by inserting a uniform plate of proton-stopping material, known as a range shifter, in the beam path at the nozzle to bring the Bragg peak within the target while ensuring high proton transport efficiency for UHDR. Aluminum has been recently proposed as a range shifter material mainly due to its high compactness and its mechanical properties. A possible drawback lies in the fact that aluminum has a larger cross-section of producing secondary neutrons compared to conventional plastic range shifters. Accordingly, an increase in secondary neutron contamination was expected during the delivery of range-modulated FLASH proton therapy, potentially heightening neutron-induced carcinogenic risks to the patient.

PURPOSE

We conducted neutron dosimetry using simulations and measurements to evaluate excess dose due to neutron exposure during UHDR proton irradiation with aluminum range shifters compared to plastic range shifters.

METHODS

Monte Carlo simulations in TOPAS were performed to investigate the secondary neutron production characteristics with aluminum range shifter during 225 MeV single-spot proton irradiation. The computational results were validated against measurements with a pair of ionization chambers in an out-of-field region ( ≤ $\le$ 30 cm) and with a Proton Recoil Scintillator-Los Alamos rem meter in a far-out-of-field region (0.5-2.5 m). The assessments were repeated with solid water slabs as a surrogate for the conventional range shifter material to evaluate the impact of aluminum on neutron yield. The results were compared with the International Electrotechnical Commission (IEC) standards to evaluate the clinical acceptance of the secondary neutron yield.

RESULTS

For a range modulation up to 26 cm in water, the maximum simulated and measured values of out-of-field secondary neutron dose equivalent per therapeutic dose with aluminum range shifter were found to be( 0.57 ± 0.02 ) mSv/Gy $(0.57\pm 0.02)\ \text{mSv/Gy}$ and( 0.46 ± 0.04 ) mSv/Gy $(0.46\pm 0.04)\ \text{mSv/Gy}$ , respectively, overall higher than the solid water cases (simulation:( 0.332 ± 0.003 ) mSv/Gy $(0.332\pm 0.003)\ \text{mSv/Gy}$ ; measurement:( 0.33 ± 0.03 ) mSv/Gy $(0.33\pm 0.03)\ \text{mSv/Gy}$ ). The maximum far out-of-field secondary neutron dose equivalent was found to be (8.8 ± 0.5 $8.8 \pm 0.5$ ) μ Sv / Gy $\umu {\rm Sv/Gy}$ and (1.62 ± 0.02 $1.62 \pm 0.02$ ) μ Sv / Gy $\umu {\rm Sv/Gy}$ for the simulations and rem meter measurements, respectively, also higher than the solid water counterparts (simulation: (3.3 ± 0.3 $3.3 \pm 0.3$ ) μ Sv / Gy $\umu {\rm Sv/Gy}$ ; measurement: (0.63 ± 0.03 $0.63 \pm 0.03$ ) μ Sv / Gy $\umu {\rm Sv/Gy}$ ).

CONCLUSIONS

We conducted simulations and measurements of secondary neutron production under proton irradiation at FLASH energy with range shifters. We found that the secondary neutron yield increased when using aluminum range shifters compared to conventional materials while remaining well below the non-primary radiation limit constrained by the IEC regulations.

Secondary cancer risk assessments following the proton therapy of lung cancer as the functions of field characteristics and patient age

INTRODUCTION

Radiation-induced secondary cancers relevant to proton therapy are still a main concern among cancer survivors. This study aims to determine the effects of age at exposure and treatment field size on radiation-induced secondary tumors following the proton therapy of lung cancer within out of field organs through the Monte Carlo (MC) simulation approach.

MATERIAL AND METHODS

A full MC model of ICRP-110 male phantom was simulated to calculate the absorbed dose corresponding to secondary radiations within distant organs from the tumor volume. Then, the risks of secondary malignancies were estimated by employing the recommended risk model by the Committee of Biological Effects of Ionizing Radiation (BEIR) for different treatment field sizes and various patient ages at exposure.

RESULTS

The results revealed that by increasing the patient age from 25 to 45 years, lifetime attributable risk (LAR) values were decreased. Maximum and minimum mortality rates were obtained for the liver and thyroid at the fixed age of 25 years, respectively. Calculated risk values for most near organs to the tumor were higher than those for distant organs. Changing the aperture size from 5 × 5 cm2 to 8 × 10 cm2 resulted in LAR increments with maximum variations of 12.5% for the stomach and a rough variation of 1.12 times in LAR for all exposure ages.

CONCLUSION

Our work on whole-body phantom addresses the impact of age at exposure and aperture size on LAR during the proton therapy of lung cancer. To minimize secondary cancer risks relevant to proton therapy of lung cancer, extra attention should be considered.

Modelling the influence of radiosensitivity on development of second primary cancer in out-of-field organs following proton therapy for paediatric cranial cancer

OBJECTIVE

Radiobiological modelling the risks of second primary cancer (SPC) after proton therapy (PT) for childhood cranial cancer remains largely unknown. Organ-specific dose-response risk factors such as radiosensitivity require exploration. This study compared the influence of radiosensitivity data (slope of βEAR) on children's lifetime attributable risks (LAR) of SPC development in out-of-field organs following cranial scattering and scanning PT.

METHODS

Out-of-field radiosensitivity parameter estimates for organs (α/β and βEAR) were sourced from literature. Physical distances for 13 out-of-field organs were measured and input into Schneider's SPC model. Sensitivity analyses were performed as a function of radiosensitivity (α/β of 1-10 Gy) and initial slope (βEAR) from Japanese/UK data to estimate the influence on the risk of radiation-induced SPC following scattering and scanning PT.

RESULTS

Models showed similar LAR of SPC estimates for age and sex-matched paediatric phantoms, however, for breast there was a significant increase using Japanese βEAR data. For most organs, scattering PT demonstrated a larger risk of LAR for SPC which increased with α/β.

CONCLUSION

Breast tissue exhibited the highest susceptibility in calculated LAR risk, demonstrating the importance for accurate data input when estimating LAR of SPC.

ADVANCES IN KNOWLEDGE

The findings of this study demonstrated younger female patients undergoing cranial proton therapy have a higher risk of developing second primary cancer of the breast tissue. Long-term multicenter registries are important to improve predictive radiobiological modelling studies of side effects.

Improved lateral penumbra for proton ocular treatments on a general-purpose spot scanning beamline

PURPOSE

An ocular applicator that fits a commercial proton snout with an upstream range shifter to allow for treatments with sharp lateral penumbra is described.

MATERIALS AND METHODS

The validation of the ocular applicator consisted of a comparison of range, depth doses (Bragg peaks and spread out Bragg peaks), point doses, and 2-D lateral profiles. Measurements were made for three field sizes, 1.5, 2, and 3 cm, resulting in 15 beams. Distal and lateral penumbras were simulated in the treatment planning system for seven range-modulation combinations for beams typical of ocular treatments and a field size of 1.5 cm, and penumbra values were compared to published literature.

RESULTS

All the range errors were within 0.5 mm. The maximum averaged local dose differences for Bragg peaks and SOBPs were 2.6% and 1.1%, respectively. All the 30 measured point doses were within +/-3% of the calculated. The measured lateral profiles, analyzed through gamma index analysis and compared to the simulated, had pass rates greater than 96% for all the planes. The lateral penumbra increased linearly with depth, from 1.4 mm at 1 cm depth to 2.5 mm at 4 cm depth. The distal penumbra ranged from 3.6 to 4.4 mm and increased linearly with the range. The treatment time for a single 10 Gy (RBE) fractional dose ranged from 30 to 120 s, depending on the shape and size of the target.

CONCLUSIONS

The ocular applicator's modified design allows lateral penumbra similar to dedicated ocular beamlines while enabling planners to use modern treatment tools such as Monte Carlo and full CT-based planning with increased flexibility in beam placement.

Comparison of out-of-field normal tissue dose estimates for pencil beam scanning proton therapy: MCNP6, PHITS, and TOPAS

Monte Carlo (MC) methods are considered the gold-standard approach to dose estimation for normal tissues outside the treatment field (out-of-field) in proton therapy. However, the physics of secondary particle production from high-energy protons are uncertain, particularly for secondary neutrons, due to challenges in performing accurate measurements. Instead, various physics models have been developed over the years to reenact these high-energy interactions based on theory. It should thus be acknowledged that MC users must currently accept some unknown uncertainties in out-of-field dose estimates. In the present study, we compared three MC codes (MCNP6, PHITS, and TOPAS) and their available physics models to investigate the variation in out-of-field normal tissue dosimetry for pencil beam scanning proton therapy patients. Total yield and double-differential (energy and angle) production of two major secondary particles, neutrons and gammas, were determined through irradiation of a water phantom at six proton energies (80, 90, 100, 110, 150, and 200 MeV). Out-of-field normal tissue doses were estimated for intracranial irradiations of 1-, 5-, and 15-year-old patients using whole-body computational phantoms. Notably, the total dose estimates for each out-of-field organ varied by approximately 25% across the three codes, independent of its distance from the treatment volume. Dose discrepancies amongst the codes were linked to the utilized physics model, which impacts the characteristics of the secondary radiation field. Using developer-recommended physics, TOPAS produced both the highest neutron and gamma doses to all out-of-field organs from all examined conditions; this was linked to its highest yields of secondary particles and second hardest energy spectra. Subsequent results when using other physics models found reduced yields and energies, resulting in lower dose estimates. Neutron dose estimates were the most impacted by physics model choice, and thus the variation in out-of-field dose estimates may be even larger than 25% when considering biological effectiveness.

EURADOS REM-COUNTER INTERCOMPARISON AT MAASTRO PROTON THERAPY CENTRE: COMPARISON WITH LITERATURE DATA

The Maastro Proton Therapy Centre is the first European facility housing the Mevion S250i Hyperscan synchrocyclotron. The proximity of the accelerator to the patient, the presence of an active pencil beam delivery system downstream of a passive energy degrader and the pulsed structure of the beam make the Mevion stray neutron field unique amongst proton therapy facilities. This paper reviews the results of a rem-counter intercomparison experiment promoted by the European Radiation Dosimetry Group at Maastro and compares them with those at other proton therapy facilities. The Maastro neutron H*(10) in the room (100-200 μSv/Gy at about 2 m from the isocentre) is in line with accelerators using purely passive or wobbling beam delivery modalities, even though Maastro shows a dose gradient peaked near the accelerator. Unlike synchrotron- and cyclotron-based facilities, the pulsed beam at Maastro requires the employment of rem-counters specifically designed to withstand pulsed neutron fields.

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 10¹⁰ 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.

Technical aspects of proton minibeam radiation therapy: Minibeam generation and delivery

Proton minibeam radiation therapy (pMBRT) is a novel therapeutic strategy that combines the normal tissue sparing of sub-millimetric, spatially fractionated beams with the improved ballistics of protons. This may allow a safe dose escalation in the tumour and has already proven to provide a remarkable increase of the therapeutic index for high-grade gliomas in animal experiments. One of the main challenges in pMBRT concerns the generation of minibeams and the implementation in a clinical environment. This article reviews the different approaches for generating minibeams, using mechanical collimators and focussing magnets, and discusses the technical aspects of the implementation and delivery of pMBRT.