A Monte Carlo study of different LET definitions and calculation parameters for proton beam therapy

Absolute energy-dependent scintillating screen calibration for real-time detection of laser-accelerated proton bunches

Laser-plasma accelerators (LPAs) can deliver pico- to nanosecond long proton bunches with ≳100 nC of charge dispersed over a broad energy spectrum. Increasing the repetition rates of today's LPAs is a necessity for their practical application. This, however, creates a need for real-time proton bunch diagnostics. Scintillating screens are one detector solution commonly applied in the field of electron LPAs for spatially resolved particle and radiation detection. Yet their establishment for LPA proton detection is only slowly taking off, also due to the lack of available calibrations. In this paper, we present an absolute proton number calibration for the scintillating screen type DRZ High (Mitsubishi Chemical Corporation, Düsseldorf, Germany), one of the most sensitive screens according to calibrations for relativistic electrons and x rays. The presented absolute light yield calibration shows an uncertainty of the proton number of 10% and can seamlessly be applied at other LPA facilities. For proton irradiation of the DRZ High screen, we find an increase in light yield of >60% compared to reference calibration data for relativistic electrons. Moreover, we investigate the scintillating screen light yield dependence on proton energy since many types of scintillators (e.g., plastic, liquid, and inorganic) show a reduced light yield for increased local energy deposition densities, an effect termed ionization quenching. The ionization quenching can reduce the light yield for low-energy protons by up to ∼20%. This work provides all necessary data for absolute spectral measurements of LPA protons with DRZ High scintillating screens, e.g., when used in the commonly applied Thomson parabola spectrometers.

A linear energy transfer distributions computation method for inhomogeneous medium by using the water equivalent ratio approximation

Dose-averaged linear energy transfer (LET), LETd is widely used in proton therapy. Compared with analytical models, Monte Carlo (MC) simulations are more accurate in obtaining LETd distributions, but they are time-consuming. This study used the 3D LETd distributions of proton beam spots in water by MC simulations as a benchmark data set. Subsequently, by combining the water equivalent ratio of various human tissues, the 3D LETd distributions of clinical cases could be quickly obtained. Our method was applied to a single spot of 160 MeV proton beam in a water-bone phantom and a pelvic case. We also computed the 3D LETd distributions for multiple proton beam spots in the pelvic case and a lung case. The results of our method were compared with the results of MC simulations, demonstrating that our method can rapidly provide 3D LETd distributions of clinical cases with acceptable differences from MC simulations.

Evaluation of monte carlo to support commissioning of the treatment planning system of new pencil beam scanning proton therapy facilities

Objective. To demonstrate the potential of Monte Carlo (MC) to support the resource-intensive measurements that comprise the commissioning of the treatment planning system (TPS) of new proton therapy facilities.Approach. Beam models of a pencil beam scanning system (Varian ProBeam) were developed in GATE (v8.2), Eclipse proton convolution superposition algorithm (v16.1, Varian Medical Systems) and RayStation MC (v12.0.100.0, RaySearch Laboratories), using the beam commissioning data. All models were first benchmarked against the same commissioning data and validated on seven spread-out Bragg peak (SOBP) plans. Then, we explored the use of MC to optimise dose calculation parameters, fully understand the performance and limitations of TPS in homogeneous fields and support the development of patient-specific quality assurance (PSQA) processes. We compared the dose calculations of the TPSs against measurements (DDTPSvs.Meas.) or GATE (DDTPSvs.GATE) for an extensive set of plans of varying complexity. This included homogeneous plans with varying field-size, range, width, and range-shifters (RSs) (n= 46) and PSQA plans for different anatomical sites (n= 11).Main results. The three beam models showed good agreement against the commissioning data, and dose differences of 3.5% and 5% were found for SOBP plans without and with RSs, respectively. DDTPSvs.Meas.and DDTPSvs.GATEwere correlated in most scenarios. In homogeneous fields the Pearson's correlation coefficient was 0.92 and 0.68 for Eclipse and RayStation, respectively. The standard deviation of the differences between GATE and measurements (±0.5% for homogeneous and ±0.8% for PSQA plans) was applied as tolerance when comparing TPSs with GATE. 72% and 60% of the plans were within the GATE predicted dose difference for both TPSs, for homogeneous and PSQA cases, respectively.Significance. Developing and validating a MC beam model early on into the commissioning of new proton therapy facilities can support the validation of the TPS and facilitate comprehensive investigation of its capabilities and limitations.

Clinical benefit of range uncertainty reduction in proton treatment planning based on dual-energy CT for neuro-oncological patients

OBJECTIVE

Several studies have shown that dual-energy CT (DECT) can lead to improved accuracy for proton range estimation. This study investigated the clinical benefit of reduced range uncertainty, enabled by DECT, in robust optimisation for neuro-oncological patients.

METHODS

DECT scans for 27 neuro-oncological patients were included. Commercial software was applied to create stopping-power ratio (SPR) maps based on the DECT scan. Two plans were robustly optimised on the SPR map, keeping the beam and plan settings identical to the clinical plan. One plan was robustly optimised and evaluated with a range uncertainty of 3% (as used clinically; denoted 3%-plan); the second plan applied a range uncertainty of 2% (2%-plan). Both plans were clinical acceptable and optimal. The dose-volume histogram parameters were compared between the two plans. Two experienced neuro-radiation oncologists determined the relevant dose difference for each organ-at-risk (OAR). Moreover, the OAR toxicity levels were assessed.

RESULTS

For 24 patients, a dose reduction >0.5/1 Gy (relevant dose difference depending on the OAR) was seen in one or more OARs for the 2%-plan; e.g. for brainstem D0.03cc in 10 patients, and hippocampus D40% in 6 patients. Furthermore, 12 patients had a reduction in toxicity level for one or two OARs, showing a clear benefit for the patient.

CONCLUSION

Robust optimisation with reduced range uncertainty allows for reduction of OAR toxicity, providing a rationale for clinical implementation. Based on these results, we have clinically introduced DECT-based proton treatment planning for neuro-oncological patients, accompanied with a reduced range uncertainty of 2%.

ADVANCES IN KNOWLEDGE

This study shows the clinical benefit of range uncertainty reduction from 3% to 2% in robustly optimised proton plans. A dose reduction to one or more OARs was seen for 89% of the patients, and 44% of the patients had an expected toxicity level decrease.

A Systematic Review of LET-Guided Treatment Plan Optimisation in Proton Therapy: Identifying the Current State and Future Needs

The well-known clinical benefits of proton therapy are achieved through higher target-conformality and normal tissue sparing than conventional radiotherapy. However, there is an increased sensitivity to uncertainties in patient motion/setup, proton range and radiobiological effect. Although recent efforts have mitigated some uncertainties, radiobiological effect remains unresolved due to a lack of clinical data for relevant endpoints. Therefore, RBE optimisations may be currently unsuitable for clinical treatment planning. LET optimisation is a novel method that substitutes RBE with LET, shifting LET hotspots outside critical structures. This review outlines the current status of LET optimisation in proton therapy, highlighting knowledge gaps and possible future research. Following the PRISMA 2020 guidelines, a search of the MEDLINE® and Scopus databases was performed in July 2023, identifying 70 relevant articles. Generally, LET optimisation methods achieved their treatment objectives; however, clinical benefit is patient-dependent. Inconsistencies in the reported data suggest further testing is required to identify therapeutically favourable methods. We discuss the methods which are suitable for near-future clinical deployment, with fast computation times and compatibility with existing treatment protocols. Although there is some clinical evidence of a correlation between high LET and adverse effects, further developments are needed to inform future patient selection protocols for widespread application of LET optimisation in proton therapy.

Proposing a Clinical Model for RBE Based on Proton Track-End Counts

PURPOSE

In proton therapy, the clinical application of linear energy transfer (LET) optimization remains contentious, in part because of challenges associated with the definition and calculation of LET and its exact relationship with relative biological effectiveness (RBE) because of large variation in experimental in vitro data. This has raised interest in other metrics with favorable properties for biological optimization, such as the number of proton track ends in a voxel. In this work, we propose a novel model for clinical calculations of RBE, based on proton track end counts.

METHODS AND MATERIALS

We developed an effective dose concept to translate between the total proton track-end count per unit mass in a voxel and a proton RBE value. Dose, track end, and dose-averaged LET (LETd) distributions were simulated using Monte Carlo models for a series of water phantoms, in vitro radiobiological studies, and patient treatment plans. We evaluated the correlation between track ends and regions of elevated biological effectiveness in comparison to LETd-based models of RBE.

RESULTS

Track ends were found to correlate with biological effects in in vitro experiments with an accuracy comparable to LETd. In patient simulations, our track end model identified the same biological hotspots as predicted by LETd-based radiobiological models of proton RBE.

CONCLUSIONS

These results suggest that, for clinical optimization and evaluation, an RBE model based on proton track end counts may match LETd-based models in terms of information provided while also offering superior statistical properties.

The OpenGATE ecosystem for Monte Carlo simulation in medical physics

This paper reviews the ecosystem of GATE, an open-source Monte Carlo toolkit for medical physics. Based on the shoulders of Geant4, the principal modules (geometry, physics, scorers) are described with brief descriptions of some key concepts (Volume, Actors, Digitizer). The main source code repositories are detailed together with the automated compilation and tests processes (Continuous Integration). We then described how the OpenGATE collaboration managed the collaborative development of about one hundred developers during almost 20 years. The impact of GATE on medical physics and cancer research is then summarized, and examples of a few key applications are given. Finally, future development perspectives are indicated.