Analysis of the track- and dose-averaged LET and LET spectra in proton therapy using the geant4 Monte Carlo code

Technical note: Determination of proton linear energy transfer from the integral depth dose

BACKGROUND

Proton linear energy transfer (LET) is associated with the relative biological effectiveness of radiation on tissues. Monte Carlo (MC) simulations have been known to be the preferred method to calculate LET. Detectors have also been built to measure LET, but they need to be calibrated with MC simulations.

PURPOSE

To propose and test a MC-free method for determining LET from the measured integral depth dose (LFI) of the protons of interest.

METHOD AND MATERIALS

LFI consists of three steps: (1) IDD measurements, (2) extraction of energy spectrum (ES) from the IDD, and (3) LET determination from the extracted ES and the stopping power of each energy. To validate the accuracy of the extraction of ES, we use Gaussian ES to synthesize IDD, extract ES from the synthesized IDD, and then compare the original (ground truth) and extracted ES. LETs calculated from the original and extracted ES are also compared. To obtain the LET of protons of interest, we measure IDDs by a large-area plane-parallel ionization chamber in water. Finally, TOPAS MC is employed to simulate IDDs, ES, and LETs. From the simulated IDD, the extracted ES and LET are compared with the simulations from TOPAS MC.

RESULTS

From the synthesized IDDs, the LETs agreed excellently when the peak energies ≥10 and 1.25 MeV with depth resolutions 0.1 and 0.01 mm, respectively. For energy <1.25 MeV, even higher depth resolution than 0.01 mm is required. From the MC simulated IDDs, our track-averaged LET excellently agreed with MC simulation, but not the LETd. Our LETd was smaller than MC simulated LETd in the shallow region but larger in the distal Bragg peak region.

CONCLUSION

LET can be accurately determined from the IDD. This method can be used in the clinic to commission or validate LETs from other measurement methods or a treatment planning system.

Comparison of linear energy transfer measurement for therapeutic carbon beam using CR-39 and TLD

The measurement of linear energy transfer (LET) is crucial for the evaluation of the radiation effect in heavy ion therapy. As two detectors which are convenient to implant into the phantom, the performance of CR-39 and thermoluminescence detector (TLD) for LET measurement was compared by experiment and simulation in this study. The results confirmed the applicability of both detectors for LET measurements, but also revealed that the CR-39 detector would lead to potential overestimation of dose-averaged LET compared with the simulation by PHITS, while the TLD would have a large uncertainty measuring ions with LET larger than 20 keVμm-1. The results of this study were expected to improve the detection method of LET for therapeutic carbon beam and would finally be benefit to the quality assurance of heavy ion radiotherapy.

Biological optimization for hybrid proton-photon radiotherapy

Objective.Hybrid proton-photon radiotherapy (RT) is a cancer treatment option to broaden access to proton RT. Additionally, with a refined treatment planning method, hybrid RT has the potential to offer superior plan quality compared to proton-only or photon-only RT, particularly in terms of target coverage and sparing organs-at-risk (OARs), when considering robustness to setup and range uncertainties. However, there is a concern regarding the underestimation of the biological effect of protons on OARs, especially those in close proximity to targets. This study seeks to develop a hybrid treatment planning method with biological dose optimization, suitable for clinical implementation on existing proton and photon machines, with each photon or proton treatment fraction delivering a uniform target dose.Approach.The proposed hybrid biological dose optimization method optimized proton and photon plan variables, along with the number of fractions for each modality, minimizing biological dose to the OARs and surrounding normal tissues. To mitigate underestimation of hot biological dose spots, proton biological dose was minimized within a ring structure surrounding the target. Hybrid plans were designed to be deliverable separately and robustly on existing proton and photon machines, with enforced uniform target dose constraints for the proton and photon fraction doses. A probabilistic formulation was utilized for robust optimization of setup and range uncertainties for protons and photons. The nonconvex optimization problem, arising from minimum monitor unit constraint and dose-volume histogram constraints, was solved using an iterative convex relaxation method.Main results.Hybrid planning with biological dose optimization effectively eliminated hot spots of biological dose, particularly in normal tissues surrounding the target, outperforming proton-only planning. It also provided superior overall plan quality and OAR sparing compared to proton-only or photon-only planning strategies.Significance.This study presents a novel hybrid biological treatment planning method capable of generating plans with reduced biological hot spots, superior plan quality to proton-only or photon-only plans, and clinical deliverability on existing proton and photon machines, separately and robustly.

Methodology for small animals targeted irradiations at conventional and ultra-high dose rates 65 MeV proton beam

As part of translational research projects, mice may be irradiated on radiobiology platforms such as the one at the ARRONAX cyclotron. Generally, these platforms do not feature an integrated imaging system. Moreover, in the context of ultra-high dose-rate radiotherapy (FLASH-RT), treatment planning should consider potential changes in the beam characteristics and internal movements in the animal. A patient-like set-up and methodology has been implemented to ensure target coverage during conformal irradiations of the brain, lungs and intestines. In addition, respiratory cycle amplitudes were quantified by fluoroscopic acquisitions on a mouse, to ensure organ coverage and to assess the impact of respiration during FLASH-RT using the 4D digital phantom MOBY. Furthermore, beam incidence direction was studied from mice µCBCT and Monte Carlo simulations. Finally,in vivodosimetry with dose-rate independent radiochromic films (OC-1) and their LET dependency were investigated. The immobilization system ensures that the animal is held in a safe and suitable position. The geometrical evaluation of organ coverage, after the addition of the margins around the organs, was satisfactory. Moreover, no measured differences were found between CONV and FLASH beams enabling a single model of the beamline for all planning studies. Finally, the LET-dependency of the OC-1 film was determined and experimentally verified with phantoms, as well as the feasibility of using these filmsin vivoto validate the targeting. The methodology developed ensures accurate and reproducible preclinical irradiations in CONV and FLASH-RT without in-room image guidance in terms of positioning, dose calculation andin vivodosimetry.

The Effects of Particle LET and Fluence on the Complexity and Frequency of Clustered DNA Damage

Motivation

Clustered DNA-lesions are predominantly induced by ionizing radiation, particularly by high-LET particles, and considered as lethal damage. Quantification of this specific type of damage as a function of radiation parameters such as LET, dose rate, dose, and particle type can be informative for the prediction of biological outcome in radiobiological studies. This study investigated the induction and complexity of clustered DNA damage for three different types of particles at an LET range of 0.5-250 keV/μm.

Methods

Nanometric volumes (36.0 nm3) of 15 base-pair DNA with its hydration shell was modeled. Electron, proton, and alpha particles at various energies were simulated to irradiate the nanometric volumes. The number of ionization events, low-energy electron spectra, and chemical yields for the formation of °OH, H°, e aq - , and H2O2 were calculated for each particle as a function of LET. Single- and double-strand breaks (SSB and DSB), base release, and clustered DNA-lesions were computed from the Monte-Carlo based quantification of the reactive species and measured yields of the species responsible for the DNA lesion formation.

Results

The total amount of DNA damage depends on particle type and LET. The number of ionization events underestimates the quantity of DNA damage at LETs higher than 10 keV/μm. Minimum LETs of 9.4 and 11.5 keV/μm are required to induce clustered damage by a single track of proton and alpha particles, respectively. For a given radiation dose, an increase in LET reduces the number of particle tracks, leading to more complex clustered DNA damage, but a smaller number of separated clustered damage sites.

Conclusions

The dependency of the number and the complexity of clustered DNA damage on LET and fluence suggests that the quantification of this damage can be a useful method for the estimation of the biological effectiveness of radiation. These results also suggest that medium-LET particles are more appropriate for the treatment of bulk targets, whereas high-LET particles can be more effective for small targets.

Deep learning-based synthetic dose-weighted LET map generation for intensity modulated proton therapy

The advantage of proton therapy as compared to photon therapy stems from the Bragg peak effect, which allows protons to deposit most of their energy directly at the tumor while sparing healthy tissue. However, even with such benefits, proton therapy does present certain challenges. The biological effectiveness differences between protons and photons are not fully incorporated into clinical treatment planning processes. In current clinical practice, the relative biological effectiveness (RBE) between protons and photons is set as constant 1.1. Numerous studies have suggested that the RBE of protons can exhibit significant variability. Given these findings, there is a substantial interest in refining proton therapy treatment planning to better account for the variable RBE. Dose-average linear energy transfer (LETd) is a key physical parameter for evaluating the RBE of proton therapy and aids in optimizing proton treatment plans. Calculating precise LETddistributions necessitates the use of intricate physical models and the execution of specialized Monte-Carlo simulation software, which is a computationally intensive and time-consuming progress. In response to these challenges, we propose a deep learning based framework designed to predict the LETddistribution map using the dose distribution map. This approach aims to simplify the process and increase the speed of LETdmap generation in clinical settings. The proposed CycleGAN model has demonstrated superior performance over other GAN-based models. The mean absolute error (MAE), peak signal-to-noise ratio and normalized cross correlation of the LETdmaps generated by the proposed method are 0.096 ± 0.019 keVμm-1, 24.203 ± 2.683 dB, and 0.997 ± 0.002, respectively. The MAE of the proposed method in the clinical target volume, bladder, and rectum are 0.193 ± 0.103, 0.277 ± 0.112, and 0.211 ± 0.086 keVμm-1, respectively. The proposed framework has demonstrated the feasibility of generating synthetic LETdmaps from dose maps and has the potential to improve proton therapy planning by providing accurate LETdinformation.

Incorporating variable RBE in IMPT optimization for ependymoma

PURPOSE

To study the dosimetric impact of incorporating variable relative biological effectiveness (RBE) of protons in optimizing intensity-modulated proton therapy (IMPT) treatment plans and to compare it with conventional constant RBE optimization and linear energy transfer (LET)-based optimization.

METHODS

This study included 10 pediatric ependymoma patients with challenging anatomical features for treatment planning. Four plans were generated for each patient according to different optimization strategies: (1) constant RBE optimization (ConstRBEopt) considering standard-of-care dose requirements; (2) LET optimization (LETopt) using a composite cost function simultaneously optimizing dose-averaged LET (LETd ) and dose; (3) variable RBE optimization (VarRBEopt) using a recent phenomenological RBE model developed by McNamara et al.; and (4) hybrid RBE optimization (hRBEopt) assuming constant RBE for the target and variable RBE for organs at risk. By normalizing each plan to obtain the same target coverage in either constant or variable RBE, we compared dose, LETd , LET-weighted dose, and equivalent uniform dose between the different optimization approaches.

RESULTS

We found that the LETopt plans consistently achieved increased LET in tumor targets and similar or decreased LET in critical organs compared to other plans. On average, the VarRBEopt plans achieved lower mean and maximum doses with both constant and variable RBE in the brainstem and spinal cord for all 10 patients. To compensate for the underdosing of targets with 1.1 RBE for the VarRBEopt plans, the hRBEopt plans achieved higher physical dose in targets and reduced mean and especially maximum variable RBE doses compared to the ConstRBEopt and LETopt plans.

CONCLUSION

We demonstrated the feasibility of directly incorporating variable RBE models in IMPT optimization. A hybrid RBE optimization strategy showed potential for clinical implementation by maintaining all current dose limits and reducing the incidence of high RBE in critical normal tissues in ependymoma patients.

The dirty and clean dose concept: Towards creating proton therapy treatment plans with a photon-like dose response

BACKGROUND

Applying tolerance doses for organs at risk (OAR) from photon therapy introduces uncertainties in proton therapy when assuming a constant relative biological effectiveness (RBE) of 1.1.

PURPOSE

This work introduces the novel dirty and clean dose concept, which allows for creating treatment plans with a more photon-like dose response for OAR and, thus, less uncertainties when applying photon-based tolerance doses.

METHODS

The concept divides the 1.1-weighted dose distribution into two parts: the clean and the dirty dose. The clean and dirty dose are deposited by protons with a linear energy transfer (LET) below and above a set LET threshold, respectively. For the former, a photon-like dose response is assumed, while for the latter, the RBE might exceed 1.1. To reduce the dirty dose in OAR, a MaxDirtyDose objective was added in treatment plan optimization. It requires setting two parameters: LET threshold and max dirty dose level. A simple geometry consisting of one target volume and one OAR in water was used to study the reduction in dirty dose in the OAR depending on the choice of the two MaxDirtyDose objective parameters during plan optimization. The best performing parameter combinations were used to create multiple dirty dose optimized (DDopt) treatment plans for two cranial patient cases. For each DDopt plan, 1.1-weighted dose, variable RBE-weighted dose using the Wedenberg RBE model and dose-average LETd distributions as well as resulting normal tissue complication probability (NTCP) values were calculated and compared to the reference plan (RefPlan) without MaxDirtyDose objectives.

RESULTS

In the water phantom studies, LET thresholds between 1.5 and 2.5 keV/µm yielded the best plans and were subsequently used. For the patient cases, nearly all DDopt plans led to a reduced Wedenberg dose in critical OAR. This reduction resulted from an LET reduction and translated into an NTCP reduction of up to 19 percentage points compared to the RefPlan. The 1.1-weighted dose in the OARs was slightly increased (patient 1: 0.45 Gy(RBE), patient 2: 0.08 Gy(RBE)), but never exceeded clinical tolerance doses. Additionally, slightly increased 1.1-weighted dose in healthy brain tissue was observed (patient 1: 0.81 Gy(RBE), patient 2: 0.53 Gy(RBE)). The variation of NTCP values due to variation of α/β from 2 to 3 Gy was much smaller for DDopt (2 percentage points (pp)) than for RefPlans (5 pp).

CONCLUSIONS

The novel dirty and clean dose concept allows for creating biologically more robust proton treatment plans with a more photon-like dose response. The reduced uncertainties in RBE can, therefore, mitigate uncertainties introduced by using photon-based tolerance doses for OAR.

Multiple Mesh-type Real Human Cell Models for Dosimetric Application Coupled with Monte Carlo Simulations

The mesh-type models are superior to voxel models in cellular dose assessment coupled with Monte Carlo codes. The aim of this study was to expand the micron-scale mesh-type models based on the fluorescence tomography of real human cells, and to investigate the feasibility of these models in the application of various irradiation scenarios and Monte Carlo codes. Six different human cell lines, including pulmonary epithelial BEAS-2B, embryonic kidney 293T, hepatocyte L-02, B-lymphoblastoid HMy2.CIR, Gastric mucosal GES-1, and intestine epithelial FHs74Int, were adopted for single mesh-type models reconstruction and optimization based on laser confocal tomography images. Mesh-type models were transformed into the format of polygon mesh and tetrahedral mesh for the GATE and PHITS Monte Carlo codes, respectively. The effect of model reduction was analyzed by dose assessment and geometry consideration. The cytoplasm and nucleus doses were obtained by designating monoenergetic electrons and protons as external irradiation, and S values with different "target-source" combinations were calculated by assigning radioisotopes as internal exposure. Four kinds of Monte Carlo codes were employed, i.e., GATE with "Livermore," "Standard" and "Standard and Geant4-DNA mixed" models for electrons and protons, as well as PHITS with "EGS" mode for electrons and radioisotopes. Multiple mesh-type real human cellular models can be applied to Monte Carlo codes directly without voxelization when combined with certain necessary surface reduction. Relative deviations between different cell types were observed among various irradiation scenarios. The relative deviation of nucleus S value reaches up to 85.65% between L-02 and GES-1 cells by 3H for the "nucleus-nucleus" combination, while that of 293T and FHs74Int nucleus dose for external beams at a 5.12 cm depth of water is 106.99%. Nucleus with smaller volume is far more affected by physical codes. There is a considerable deviation for dose within BEAS-2B at the nanoscale. The multiple mesh-type real cell models were more versatile than voxel models and mathematical models. The present study provided several models which can easily be extended to other cell types and irradiation scenarios for RBE estimations and biological effect predictions, including radiation biological experiments, radiotherapy and radiation protection.

A simulation study of in-beam visualization system for proton therapy by monitoring scattered protons

Recently, in-beam positron emission tomography (PET) has been actively researched for reducing biological washout effects and dose monitoring during irradiation. However, the positron distribution does not precisely reflect the dose distribution since positron production and ionization are completely different physical processes. Thus, a novel in-beam system was proposed to determine proton dose range by measuring scattered protons with dozens of scintillation detectors surrounding the body surface. While previous studies conducted a preliminary experiment with a simple phantom, we simulated more complex situations in this paper. Especially, we conducted three stepwise simulation studies to demonstrate the feasibility of the proposed method. First, a simple rectangular phantom was reproduced on simulation and irradiated with protons for obtaining current values and Monte Carlo (MC) dose. Next, we trained a deep learning model to estimate 2-dimensional-dose range (2D-DL dose) from measured current values for simulation (A). We simulated plastic scintillators as detectors to measure the scattered protons. Second, a rectangular phantom with an air layer was used, and 3D-DL dose was estimated in simulation (B). Finally, a cylindrical phantom that mimics the human body was used for confirming the estimation quality of the simulation (C). Consequently, the position of the Bragg peak was estimated with an error of 1.0 mm in simulation (A). In addition, the position of the air layer, as well as the verifying peak position with an error of 2.1 mm, was successfully estimated in simulation (B). Although the estimation error of the peak position was 12.6 mm in simulation (C), the quality was successfully further improved to 9.3 mm by incorporating the mass density distribution obtained from the computed tomography (CT). These simulation results demonstrated the potential of the as-proposed verification system. Additionally, the effectiveness of CT utilization for estimating the DL dose was also indicated.

Biophysical characterization of collimated and uncollimated fields in pencil beam scanning proton therapy

Objective. The lateral dose fall-off in proton pencil beam scanning (PBS) technique remains the preferred choice for sparing adjacent organs at risk as opposed to the distal edge due to the proton range uncertainties and potentially high relative biological effectiveness. However, because of the substantial spot size along with the scattering in the air and in the patient, the lateral penumbra in PBS can be degraded. Combining PBS with an aperture can result in a sharper dose fall-off, particularly for shallow targets.Approach. The aim of this work was to characterize the radiation fields produced by collimated and uncollimated 100 and 140 MeV proton beams, using Monte Carlo simulations and measurements with a MiniPIX-Timepix detector. The dose and the linear energy transfer (LET) were then coupled with publishedin silicobiophysical models to elucidate the potential biological effects of collimated and uncollimated fields.Main results. Combining an aperture with PBS reduced the absorbed dose in the lateral fall-off and out-of-field by 60%. However, the results also showed that the absolute frequency-averaged LET (LET_F) values increased by a maximum of 3.5 keVμm^-1in collimated relative to uncollimated fields, while the dose-averaged LET (LET_D) increased by a maximum of 7 keVμm^-1. Despite the higher LET values produced by collimated fields, the predicted DNA damage yields remained lower, owing to the large dose reduction.Significance. This work demonstrated the dosimetric advantages of combining an aperture with PBS coupled with lower DNA damage induction. A methodology for calculating dose in water derived from measurements with a silicon-based detector was also presented. This work is the first to demonstrate experimentally the increase in LET caused by combining PBS with aperture, and to assess the potential DNA damage which is the initial step in the cascade of events leading to the majority of radiation-induced biological effects.

Dosimetric response of Gafchromic™ EBT-XD film to therapeutic protons

EBT-XD model of Gafchromic™ films has a broader optimal dynamic dose range, up to 40 Gy, compared to its predecessor models. This characteristic has made EBT-XD films suitable for high-dose applications such as stereotactic body radiotherapy and stereotactic radiosurgery, as well as ultra-high dose rate FLASH radiotherapy. The purpose of the current study was to characterize the dependence of EBT-XD film response on linear energy transfer (LET) and dose rate of therapeutic protons from a synchrotron. A clinical spot-scanning proton beam was used to study LET dependence at three dose-averaged LET (LETd) values of 1.0 keV/µm, 3.6 keV/µm, and 7.6 keV/µm. A research proton beamline was used to study dose rate dependence at 150 Gy/second in the FLASH mode and 0.3 Gy/second in the non-FLASH mode. Film response data from LETd values of 0.9 keV/µm and 9.0 keV/µm of the proton FLASH beam were also compared. Film response data from a clinical 6 MV photon beam were used as a reference. Both gray value method and optical density (OD) method were used in film calibration. Calibration results using a specific OD calculation method and a generic OD calculation method were compared. The four-parameter NIH Rodbard function and three-parameter rational function were compared in fitting the calibration curves. Experimental results showed that the response of EBT-XD film is proton LET dependent but independent of dose rate. Goodness-of-fit analysis showed that using the NIH Rodbard function is superior for both protons and photons. Using the "specific OD + NIH Rodbard function" method for EBT-XD film calibration is recommended.

Linear Energy Transfer and Relative Biological Effectiveness Investigation of Various Structures for a Cohort of Proton Patients With Brain Tumors

Purpose

The current knowledge on biological effects associated with proton therapy is limited. Therefore, we investigated the distributions of dose, dose-averaged linear energy transfer (LETd), and the product between dose and LETd (DLETd) for a patient cohort treated with proton therapy. Different treatment planning system features and visualization tools were explored.

Methods and Materials

For a cohort of 24 patients with brain tumors, the LETd, DLETd, and dose was calculated for a fixed relative biological effectiveness value and 2 variable models: plan-based and phenomenological. Dose threshold levels of 0, 5, and 20 Gy were imposed for LETd visualization. The relationship between physical dose and LETd and the frequency of LETd hotspots were investigated.

Results

The phenomenological relative biological effectiveness model presented consistently higher dose values. For lower dose thresholds, the LETd distribution was steered toward higher values related to low treatment doses. Differences up to 26.0% were found according to the threshold. Maximum LETd values were identified in the brain, periventricular space, and ventricles. An inverse relationship between LETd and dose was observed. Frequency information to the domain of dose and LETd allowed for the identification of clusters, which steer the mean LETd values, and the identification of higher, but sparse, LETd values.

Conclusions

Identifying, quantifying, and recording LET distributions in a standardized fashion is necessary, because concern exists over a link between toxicity and LET hotspots. Visualizing DLETd or dose × LETd during treatment planning could allow for clinicians to make informed decisions.

Technical note: In silico benchmarking of the linear energy transfer-based functionalities for carbon ion beams in a commercial treatment planning system

BACKGROUND

The increasing number of studies dealing with linear energy transfer (LET)-based evaluation and optimization in the field of carbon ion radiotherapy (CIRT) indicates the rising demand for LET implementation in commercial treatment planning systems (TPS). Benchmarking studies could play a key role in detecting (and thus preventing) computation errors prior implementing such functionalities in a TPS.

PURPOSE

This in silico study was conducted to benchmark the following two LET-related functionalities in a commercial TPS against Monte Carlo simulations: (1) dose averaged LET (LET_d ) scoring and (2) physical dose filtration based on LET for future LET-based treatment plan evaluation and optimization studies.

METHODS

The LET_d scoring and LET-based dose filtering (in which the deposited dose can be separated into the dose below and above the user specified LET threshold) functionalities for carbon ions in the research version RayStation (RS) 9A-IonPG TPS (RaySearch Laboratories, Sweden) were benchmarked against GATE/Geant4 simulations. Pristine Bragg peaks (BPs) and cuboid targets, positioned at different depths in a homogeneous water phantom and a setup with heterogeneity were used for this study.

RESULTS

For all setups (homogeneous and heterogeneous), the mean absolute (and relative) LET_d difference was less than 1 keV/ μ $\umu$ m (3.5%) in the plateau and target and less than 2 keV/ μ $\umu$ m (8.3%) in the fragmentation tail. The maximum local differences were 4 and 6 keV/ μ $\umu$ m, respectively. The mean absolute (and relative) physical dose differences for both low- and high-LET doses were less than 1 cGy (1.5%) in the plateau, target and tail with a maximum absolute difference of 2 cGy.

CONCLUSIONS

No computation error was found in the tested functionalities except for LET_d in lateral direction outside the target, showing the limitation of the implemented monochrome model in the tested TPS version.

Short-lived radioactive8Li and8He ions for hadrontherapy: a simulation study

Purpose.Although charged particle therapy (CPT) for cancer treatment has grown these past years, the use of protons and carbon ions for therapy remains debated compared to x-ray therapy. While a biological advantage of protons is not clearly demonstrated, therapy using carbon ions is often pointed out for its high cost. Furthermore, the nuclear interactions undergone by carbons inside the patient are responsible for an additional dose delivered after the Bragg peak, which deteriorates the ballistic advantage of CPT. Therefore, a renewed interest for lighter ions with higher biological efficiency than protons was recently observed. In this context, helium and lithium ions represent a good compromise between protons and carbons, as they exhibit a higher linear energy transfer (LET) than protons in the Bragg peak and can be accelerated by cyclotrons. The possibility of accelerating radioactive⁸Li, decaying in 2α-particles, and⁸He, decaying in⁸Li byβ^-decay, is particularly interesting.Methods. This work aims to assess the interest of the use of⁸Li and⁸He ions for therapy by Monte Carlo simulations carried out withGeant4.Results. It was calculated that the⁸Li and⁸He decay results in an increase of the LET of almost a factor 2 in the Bragg peak compared to stable⁷Li and⁴He. This results also in a higher dose deposited in the Bragg peak without an increase of the dose in the plateau region. It was also shown that both⁸He and⁸Li can have a potential interest for prompt-gamma monitoring techniques. Finally, the feasibility of accelerating facilities delivering⁸Li and⁸He was also discussed.Conclusion. In this study, we demonstrate that both⁸Li and⁸He have interesting properties for therapy. Indeed, simulations predict that⁸Li and⁸He are a good compromise between proton and¹²C, both in terms of LET and dose.

An iterative convex relaxation method for proton LET optimization

Objective:A constant relative biological effectiveness of 1.1 in current clinical practice of proton radiotherapy (RT) is a crude approximation and may severely underestimate the biological dose from proton RT to normal tissues, especially near the treatment target at the end of Bragg peaks that exhibits high linear energy transfer (LET). LET optimization can account for biological effectiveness of protons during treatment planning, for minimizing biological proton dose and hot spots to normal tissues. However, the LET optimization is usually nonlinear and nonconvex to solve, for which this work will develop an effective optimization method based on iterative convex relaxation (ICR).Approach: In contrast to the generic nonlinear optimization method, such as Quasi-Newton (QN) method, that does not account for specific characteristics of LET optimization, ICR is tailored to LET modeling and optimization in order to effectively and efficiently solve the LET problem. Specifically, nonlinear dose-averaged LET term is iteratively linearized and becomes convex during ICR, while nonconvex dose-volume constraint and minimum-monitor-unit constraint are also handled by ICR, so that the solution for LET optimization is obtained by solving a sequence of convex and linearized convex subproblems. Since the high LET mostly occurs near the target, a 1 cm normal-tissue expansion of clinical target volume (CTV) (excluding CTV), i.e. CTV1cm, is defined to as an auxiliary structure during treatment planning, where LET is minimized.Main results: ICR was validated in comparison with QN for abdomen, lung, and head-and-neck cases. ICR was effective for LET optimization, as ICR substantially reduced the LET and biological dose in CTV1cm the ring, with preserved dose conformality to CTV. Compared to QN, ICR had smaller LET, physical and biological dose in CTV1cm, and higher conformity index values; ICR was also computationally more efficient, which was about 3 times faster than QN.Significance: A LET-specific optimization method based on ICR has been developed for solving proton LET optimization, which has been shown to be more computationally efficient than generic nonlinear optimizer via QN, with better plan quality in terms of LET, biological and physical dose conformality.

Characterisation and Quenching Correction for an Al2O3:C Optical Fibre Real Time System in Therapeutic Proton, Helium, and Carbon-Charged Beams

Real time radioluminescence fibre-based detectors were investigated for application in proton, helium, and carbon therapy dosimetry. The Al₂O₃:C probes are made of one single crystal (1 mm) and two droplets of micro powder in two sizes (38 μm and 4 μm) mixed with a water-equivalent binder. The fibres were irradiated behind different thicknesses of solid slabs, and the Bragg curves presented a quenching effect attributed to the nonlinear response of the radioluminescence (RL) signal as a function of linear energy transfer (LET). Experimental data and Monte Carlo simulations were utilised to acquire a quenching correction method, adapted from Birks' formulation, to restore the linear dose-response for particle therapy beams. The method for quenching correction was applied and yielded the best results for the '4 μm' optical fibre probe, with an agreement at the Bragg peak of 1.4% (160 MeV), and 1.5% (230 MeV) for proton-charged particles; 2.4% (150 MeV/u) for helium-charged particles and of 4.8% (290 MeV/u) and 2.9% (400 MeV/u) for the carbon-charged particles. The most substantial deviations for the '4 μm' optical fibre probe were found at the falloff regions, with ~3% (protons), ~5% (helium) and 6% (carbon).

Characterisation and Quenching Correction for an Al2O3:C Optical Fibre Real Time System in Therapeutic Proton, Helium, and Carbon-Charged Beams

Real time radioluminescence fibre-based detectors were investigated for application in proton, helium, and carbon therapy dosimetry. The Al2O3:C probes are made of one single crystal (1 mm) and two droplets of micro powder in two sizes (38 μm and 4 μm) mixed with a water-equivalent binder. The fibres were irradiated behind different thicknesses of solid slabs, and the Bragg curves presented a quenching effect attributed to the nonlinear response of the radioluminescence (RL) signal as a function of linear energy transfer (LET). Experimental data and Monte Carlo simulations were utilised to acquire a quenching correction method, adapted from Birks' formulation, to restore the linear dose-response for particle therapy beams. The method for quenching correction was applied and yielded the best results for the '4 μm' optical fibre probe, with an agreement at the Bragg peak of 1.4% (160 MeV), and 1.5% (230 MeV) for proton-charged particles; 2.4% (150 MeV/u) for helium-charged particles and of 4.8% (290 MeV/u) and 2.9% (400 MeV/u) for the carbon-charged particles. The most substantial deviations for the '4 μm' optical fibre probe were found at the falloff regions, with ~3% (protons), ~5% (helium) and 6% (carbon).

Targeting hippocampal neurogenesis to protect astronauts' cognition and mood from decline due to space radiation effects

Neurogenesis is an essential, lifelong process during which neural stem cells generate new neurons within the hippocampus, a center for learning, memory, and mood control. Neural stem cells are vulnerable to environmental insults spanning from chronic stress to radiation. These insults reduce their numbers and diminish neurogenesis, leading to memory decline, anxiety, and depression. Preserving neural stem cells could thus help prevent these neurogenesis-associated pathologies, an outcome particularly important for long-term space missions where environmental exposure to radiation is significantly higher than on Earth. Multiple developments, from mechanistic discoveries of radiation injury on hippocampal neurogenesis to new platforms for the development of selective, specific, effective, and safe small molecules as neurogenesis-protective agents hold great promise to minimize radiation damage on neurogenesis. In this review, we summarize the effects of space-like radiation on hippocampal neurogenesis. We then focus on current advances in drug discovery and development and discuss the nuclear receptor TLX/NR2E1 (oleic acid receptor) as an example of a neurogenic target that might rescue neurogenesis following radiation.

Adaptation and dosimetric commissioning of a synchrotron-based proton beamline for FLASH experiments

Objective. Irradiation with ultra-high dose rates (>40 Gy s−1), also known as FLASH irradiation, has the potential to shift the paradigm of radiation therapy because of its reduced toxicity to normal tissues compared to that of conventional irradiations. The goal of this study was to (1) achieve FLASH irradiation conditions suitable for pre-clinical i n vitro and in vivo biology experiments using our synchrotron-based proton beamline and (2) commission the FLASH irradiation conditions achieved. Approach. To achieve these suitable FLASH conditions, we made a series of adaptations to our proton beamline, including modifying the spill length and size of accelerating cycles, repurposing the reference monitor for dose control, and expanding the field size with a custom double-scattering system. We performed the dosimetric commissioning with measurements using an Advanced Markus chamber and EBT-XD films as well as with Monte Carlo simulations. Main results. Through adaptations, we have successfully achieved FLASH irradiation conditions, with an average dose rate of up to 375 Gy s−1. The Advanced Markus chamber was shown to be appropriate for absolute dose calibration under our FLASH conditions with a recombination factor ranging from 1.002 to 1.006 because of the continuous nature of our synchrotron-based proton delivery within a spill. Additionally, the absolute dose measured using the Advanced Markus chamber and EBT-XD films agreed well, with average and maximum differences of 0.32% and 1.63%, respectively. We also performed a comprehensive temporal analysis for FLASH spills produced by our system, which helped us identify a unique relationship between the average dose rate and the dose in our FLASH irradiation. Significance. We have established a synchrotron-based proton FLASH irradiation platform with accurate and precise dosimetry that is suitable for pre-clinical biology experiments. The unique time structure of the FLASH irradiation produced by our synchrotron-based system may shed new light onto the mechanism behind the FLASH effect.

Quantification of biological range uncertainties in patients treated at the Krakow proton therapy centre

BACKGROUND

Variable relative biological effectiveness (vRBE) in proton therapy might significantly modify the prediction of RBE-weighted dose delivered to a patient during proton therapy. In this study we will present a method to quantify the biological range extension of the proton beam, which results from the application of vRBE approach in RBE-weighted dose calculation.

METHODS AND MATERIALS

The treatment plans of 95 patients (brain and skull base patients) were used for RBE-weighted dose calculation with constant and the McNamara RBE model. For this purpose the Monte Carlo tool FRED was used. The RBE-weighted dose distributions were analysed using indices from dose-volume histograms. We used the volumes receiving at least 95% of the prescribed dose (V95) to estimate the biological range extension resulting from vRBE approach.

RESULTS

The vRBE model shows higher median value of relative deposited dose and D95 in the planning target volume by around 1% for brain patients and 4% for skull base patients. The maximum doses in organs at risk calculated with vRBE was up to 14 Gy above dose limit. The mean biological range extension was greater than 0.4 cm.

DISCUSSION

Our method of estimation of biological range extension is insensitive for dose inhomogeneities and can be easily used for different proton plans with intensity-modulated proton therapy (IMPT) optimization. Using volumes instead of dose profiles, which is the common method, is more universal. However it was tested only for IMPT plans on fields arranged around the tumor area.

CONCLUSIONS

Adopting a vRBE model results in an increase in dose and an extension of the beam range, which is especially disadvantageous in cancers close to organs at risk. Our results support the need to re-optimization of proton treatment plans when considering vRBE.

Quantifying Systematic RBE-Weighted Dose Uncertainty Arising from Multiple Variable RBE Models in Organ at Risk

PURPOSE

Relative biological effectiveness (RBE) uncertainties have been a concern for treatment planning in proton therapy, particularly for treatment sites that are near organs at risk (OARs). In such a clinical situation, the utilization of variable RBE models is preferred over constant RBE model of 1.1. The problem, however, lies in the exact choice of RBE model, especially when current RBE models are plagued with a host of uncertainties. This paper aims to determine the influence of RBE models on treatment planning, specifically to improve the understanding of the influence of the RBE models with regard to the passing and failing of treatment plans. This can be achieved by studying the RBE-weighted dose uncertainties across RBE models for OARs in cases where the target volume overlaps the OARs. Multi-field optimization (MFO) and single-field optimization (SFO) plans were compared in order to recommend which technique was more effective in eliminating the variations between RBE models.

METHODS

Fifteen brain tumor patients were selected based on their profile where their target volume overlaps with both the brain stem and the optic chiasm. In this study, 6 RBE models were analyzed to determine the RBE-weighted dose uncertainties. Both MFO and SFO planning techniques were adopted for the treatment planning of each patient. RBE-weighted dose uncertainties in the OARs are calculated assuming( α β ) x of 3 Gy and 8 Gy. Statistical analysis was used to ascertain the differences in RBE-weighted dose uncertainties between MFO and SFO planning. Additionally, further investigation of the linear energy transfer (LET) distribution was conducted to determine the relationship between LET distribution and RBE-weighted dose uncertainties.

RESULTS

The results showed no strong indication on which planning technique would be the best for achieving treatment planning constraints. MFO and SFO showed significant differences (P <.05) in the RBE-weighted dose uncertainties in the OAR. In both clinical target volume (CTV)-brain stem and CTV-chiasm overlap region, 10 of 15 patients showed a lower median RBE-weighted dose uncertainty in MFO planning compared with SFO planning. In the LET analysis, 8 patients (optic chiasm) and 13 patients (brain stem) showed a lower mean LET in MFO planning compared with SFO planning. It was also observed that lesser RBE-weighted dose uncertainties were present with MFO planning compared with SFO planning technique.

CONCLUSIONS

Calculations of the RBE-weighted dose uncertainties based on 6 RBE models and 2 different( α β ) x revealed that MFO planning is a better option as opposed to SFO planning for cases of overlapping brain tumor with OARs in eliminating RBE-weighted dose uncertainties. Incorporation of RBE models failed to dictate the passing or failing of a treatment plan. To eliminate RBE-weighted dose uncertainties in OARs, the MFO planning technique is recommended for brain tumor when CTV and OARs overlap.

Towards harmonizing clinical linear energy transfer (LET) reporting in proton radiotherapy: a European multi-centric study

BACKGROUND

Clinical data suggest that the relative biological effectiveness (RBE) in proton therapy (PT) varies with linear energy transfer (LET). However, LET calculations are neither standardized nor available in clinical routine. Here, the status of LET calculations among European PT institutions and their comparability are assessed.

MATERIALS AND METHODS

Eight European PT institutions used suitable treatment planning systems with their center-specific beam model to create treatment plans in a water phantom covering different field arrangements and fulfilling commonly agreed dose objectives. They employed their locally established LET simulation environments and procedures to determine the corresponding LET distributions. Dose distributions D1.1 and DRBE assuming constant and variable RBE, respectively, and LET were compared among the institutions. Inter-center variability was assessed based on dose- and LET-volume-histogram parameters.

RESULTS

Treatment plans from six institutions fulfilled all clinical goals and were eligible for common analysis. D1.1 distributions in the target volume were comparable among PT institutions. However, corresponding LET values varied substantially between institutions for all field arrangements, primarily due to differences in LET averaging technique and considered secondary particle spectra. Consequently, DRBE using non-harmonized LET calculations increased inter-center dose variations substantially compared to D1.1 and significantly in mean dose to the target volume of perpendicular and opposing field arrangements (p < 0.05). Harmonizing LET reporting (dose-averaging, all protons, LET to water or to unit density tissue) reduced the inter-center variability in LET to the order of 10-15% within and outside the target volume for all beam arrangements. Consequentially, inter-institutional variability in DRBE decreased to that observed for D1.1.

CONCLUSION

Harmonizing the reported LET among PT centers is feasible and allows for consistent multi-centric analysis and reporting of tumor control and toxicity in view of a variable RBE. It may serve as basis for harmonized variable RBE dose prescription in PT.

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

The strongin vitroevidence that proton Relative Biological Effectiveness (RBE) varies with Linear Energy Transfer (LET) has led to an interest in applying LET within treatment planning. However, there is a lack of consensus on LET definition, Monte Carlo (MC) parameters or clinical methodology. This work aims to investigate how common variations of LET definition may affect potential clinical applications. MC simulations (GATE/GEANT4) were used to calculate absorbed dose and different types of LET for a simple Spread Out Bragg Peak (SOBP) and for four clinical PBT plans covering a range of tumour sites. Variations in the following LET calculation methods were considered: (i) averaging (dose-averaged LET (LET_d) & track-averaged LET); (ii) scoring (LET_dto water, to medium and to mass density); (iii) particle inclusion (LET_dto all protons, to primary protons and to particles); (iv) MC settings (hit type and Maximum Step Size (MSS)). LET distributions were compared using: qualitative comparison, LET Volume Histograms (LVHs), single value criteria (maximum and mean values) and optimised LET-weighted dose models. Substantial differences were found between LET values in averaging, scoring and particle type. These differences depended on the methodology, but for one patient a difference of ∼100% was observed between the maximum LET_dfor all particles and maximum LET_dfor all protons within the brainstem in the high isodose region (4 keVμm^-1and 8 keVμm^-1respectively). An RBE model using LET_dincluding heavier ions was found to predict substantially different LET-weighted dose compared to those using other LET definitions. In conclusion, the selection of LET definition may affect the results of clinical metrics considered in treatment planning and the results of an RBE model. The authors' advocate for the scoring of dose-averaged LET to water for primary and secondary protons using a random hit type and automated MSS.

Investigating beam range uncertainty in proton prostate treatment using pelvic-like biological phantoms

This study aims to develop a method for verifying site-specific and/or beam path specific proton beam range, which could reduce range uncertainty margins and the associated treatment complications. It investigates the range uncertainties from both CT HU to relative stopping power conversion and patient positioning errors for prostate treatment using pelvic-like biological phantoms. Three 25 × 14 × 12 cm³phantoms, made of fresh animal tissues mimicking the pelvic anatomies of prostate patients, were scanned with a general electric CT simulator. A 22 cm circular passive scattering beam with 29 cm range and 8 cm modulation width was used to measure the water equivalent path lengths (WEPL) through the phantoms at multiple points using the dose extinction method with a MatriXXPT detector. The measured WEPLs were compared to those predicted by TOPAS simulations and ray-tracing WEPL calculations. For the three phantoms, the WEPL differences between measured and theoretical prediction (WDMT) are below 1.8% for TOPAS, and 2.5% for ray-tracing. WDMT varies with phantom anatomies by about 0.5% for both TOPAS and ray-tracing. WDMT also correlates with the tissue types of a specific treated region. For the regions where the proton beam path is parallel to sharp bone edges, the WDMTs of TOPAS and ray-tracing respectively reach up to 1.8% and 2.5%. For the region where proton beams pass through just soft tissues, the WDMT is mostly less than 1% for both TOPAS and ray-tracing. For prostate treatments, range uncertainty depends on the tissue types within a specific treated region, patient anatomies and the range calculation methods in the planning algorithms. Our study indicates range uncertainty is less than 2.5% for the whole treated region with both ray-tracing and TOPAS, which suggests the potential to reduce the current 3.5% range uncertainty margin used in the clinics by at least 1% even for single-energy CT data.

Uncertainty in tissue equivalent proportional counter assessments of microdosimetry and RBE estimates in carbon radiotherapy

Microdosimetry is an important tool for assessing energy deposition distributions from ionizing radiation at cellular and cellular nucleus scales. It has served as an input parameter for multiple common mathematical models, including evaluation of relative biological effectiveness (RBE) of carbon ion therapy. The most common detector used for microdosimetry is the tissue-equivalent proportional counter (TEPC). Although it is widely applied, TEPC has various inherent uncertainties. Therefore, this work quantified the magnitude of TEPC measurement uncertainties and their impact on RBE estimates for therapeutic carbon beams. Microdosimetric spectra and frequency-, dose-, and saturation-corrected dose-mean lineal energy (****) were calculated using the Monte Carlo toolkit Geant4 for five monoenergetic and three spread-out Bragg peak carbon beams in water at every millimeter along the central beam axis. We simulated the following influences on these spectra from eight sources of uncertainty: wall effects, pulse pile-up, electronics, gas pressure, W-value, gain instability, low energy cut-off, and counting statistics. Statistic uncertainty was quantified as the standard deviation of perturbed values for each source. Bias was quantified as the difference between default lineal energy values and the mean of perturbed values for each systematic source. Uncertainties were propagated to RBE using the modified microdosimetric kinetic model (MKM). Variance introduced by statistic sources iny¯Fandy¯Daveraged 3.8% and 3.4%, respectively, and 1.5% iny*across beam depths and energies. Bias averaged 6.2% and 7.3% iny¯Fandy¯D,and 4.8% iny*.These uncertainties corresponded to 1.2 ± 0.9% on average in RBE_MKM. The largest contributors to variance and bias were pulse pile-up and wall effects. This study established an error budget for microdosimetric carbon measurements by quantifying uncertainty inherent to TEPC measurements. It is necessary to understand how robust the measurement of RBE model input parameters are against this uncertainty in order to verify clinical model implementation.

Radiosensitization Effect of Gold Nanoparticles in Proton Therapy

The number of proton therapy facilities and the clinical usage of high energy proton beams for cancer treatment has substantially increased over the last decade. This is mainly due to the superior dose distribution of proton beams resulting in a reduction of side effects and a lower integral dose compared to conventional X-ray radiotherapy. More recently, the usage of metallic nanoparticles as radiosensitizers to enhance radiotherapy is receiving growing attention. While this strategy was originally intended for X-ray radiotherapy, there is currently a small number of experimental studies indicating promising results for proton therapy. However, most of these studies used low proton energies, which are less applicable to clinical practice; and very small gold nanoparticles (AuNPs). Therefore, this proof of principle study evaluates the radiosensitization effect of larger AuNPs in combination with a 200 MeV proton beam. CHO-K1 cells were exposed to a concentration of 10 μg/ml of 50 nm AuNPs for 4 hours before irradiation with a clinical proton beam at NRF iThemba LABS. AuNP internalization was confirmed by inductively coupled mass spectrometry and transmission electron microscopy, showing a random distribution of AuNPs throughout the cytoplasm of the cells and even some close localization to the nuclear membrane. The combined exposure to AuNPs and protons resulted in an increase in cell killing, which was 27.1% at 2 Gy and 43.8% at 6 Gy, compared to proton irradiation alone, illustrating the radiosensitizing potential of AuNPs. Additionally, cells were irradiated at different positions along the proton depth-dose curve to investigate the LET-dependence of AuNP radiosensitization. An increase in cytogenetic damage was observed at all depths for the combined treatment compared to protons alone, but no incremental increase with LET could be determined. In conclusion, this study confirms the potential of 50 nm AuNPs to increase the therapeutic efficacy of proton therapy.

Exploratory Investigation of Dose-Linear Energy Transfer (LET) Volume Histogram (DLVH) for Adverse Events Study in Intensity Modulated Proton Therapy (IMPT)

PURPOSE

We proposed a novel tool-a dose linear energy transfer (LET)-volume histogram (DLVH)-and performed an exploratory study to investigate rectal bleeding in prostate cancer treated with intensity modulated proton therapy.

METHODS AND MATERIALS

The DLVH was constructed with dose and LET as 2 axes, and the normalized volume of the structure was contoured in the dose-LET plane as isovolume lines. We defined the DLVH index, DLv%(d,l) (ie, v% of the structure) to have a dose of ≥d Gy and an LET of ≥l keV/μm, similar to the dose-volume histogram index Dv%. Nine patients with prostate cancer with rectal bleeding (Common Terminology Criteria for Adverse Events grade ≥2) were included as the adverse event group, and 48 patients with no complications were considered the control group. A P value map was constructed by comparison of the DLVH indices of all patients between the 2 groups using the Mann-Whitney U test. Dose-LET volume constraints (DLVCs) were derived based on the P value map with a manual selection procedure facilitated by Spearman's correlation tests. The obtained DLVCs were further cross-validated using a multivariate support vector machine (SVM)-based normal tissue complication probability (NTCP) model with an independent testing data set composed of 8 adverse event and 13 control patients.

RESULTS

We extracted 2 DLVC constraints. One DLVC was obtained, Vdose/LETboundary:2.5keVμmat 75 Gy to 3.2keVμmat8.65Gy <1.27% (DLVC1), revealing a high LET volume effect. The second DLVC, V(72.2Gy,0keVμm) < 2.23% (DVLC2), revealed a high dose volume effect. The SVM-based NTCP model with 2 DLVCs provided slightly superior performance than using dose only, with an area under the curve of 0.798 versus 0.779 for the testing data set.

CONCLUSIONS

Our results demonstrated the importance of rectal "hot spots" in both high LET (DLVC1) and high dose (DLVC2) in inducing rectal bleeding. The SVM-based NTCP model confirmed the derived DLVCs as good predictors for rectal bleeding when intensity modulated proton therapy is used to treat prostate cancer.

A Critical Review of LET-Based Intensity-Modulated Proton Therapy Plan Evaluation and Optimization for Head and Neck Cancer Management

In this review article, we review the 3 important aspects of linear-energy-transfer (LET) in intensity-modulated proton therapy (IMPT) for head and neck (H&N) cancer management. Accurate LET calculation methods are essential for LET-guided plan evaluation and optimization, which can be calculated either by analytical methods or by Monte Carlo (MC) simulations. Recently, some new 3D analytical approaches to calculate LET accurately and efficiently have been proposed. On the other hand, several fast MC codes have also been developed to speed up the MC simulation by simplifying nonessential physics models and/or using the graphics processor unit (GPU)–acceleration approach. Some concepts related to LET are also briefly summarized including (1) dose-weighted versus fluence-weighted LET; (2) restricted versus unrestricted LET; and (3) microdosimetry versus macrodosimetry. LET-guided plan evaluation has been clinically done in some proton centers. Recently, more and more studies using patient outcomes as the biological endpoint have shown a positive correlation between high LET and adverse events sites, indicating the importance of LET-guided plan evaluation in proton clinics. Various LET-guided plan optimization methods have been proposed to generate proton plans to achieve biologically optimized IMPT plans. Different optimization frameworks were used, including 2-step optimization, 1-step optimization, and worst-case robust optimization. They either indirectly or directly optimize the LET distribution in patients while trying to maintain the same dose distribution and plan robustness. It is important to consider the impact of uncertainties in LET-guided optimization (ie, LET-guided robust optimization) in IMPT, since IMPT is sensitive to uncertainties including both the dose and LET distributions. We believe that the advancement of the LET-guided plan evaluation and optimization will help us exploit the unique biological characteristics of proton beams to improve the therapeutic ratio of IMPT to treat H&N and other cancers.

Design and characterization of a prototype tertiary device for proton beam stereotactic radiosurgery

Though potentially beneficial, proton beam stereotactic radiosurgery has not been adopted widely secondary to the technical challenge of safely delivering multiple focused beams of proton radiation. In this study, we describe the design and characterization of a proton beam stereotactic radiosurgery system that can be adopted by existing passive scattering systems. This system utilizes a helmet-like device in which patient-specific brass apertures required for final beam collimation are positioned on a scaffold that is separate from the treatment gantry. The proton snout is then fitted with a generic aperture to focus the primary proton beam onto the patient specific apertures that are in the helmet-like device. The patient-specific apertures can all be placed at the start of the treatment, thus treatment with multiple beams can be accomplished without the delay of switching the apertures. In this report we describe a prototype design of this collimation system and dosimetric testing to verify efficacy. Subsequently, we describe a custom 3D printing of a prototype device and report on overall localization accuracy using Winston-Lutz tests. Our results show that it is possible to develop an add-on device for proton beam radiosurgery that is safe and efficient and capable of wide adoption on existing proton delivery systems.

Measurement of therapeutic12C beam in a water phantom using CR-39

The motivation for this study was to explore a new method to test the particle spatial distribution for a therapeutic carbon beam. CR-39 plastic nuclear track detectors were irradiated to a 276.5 MeV u^-1mono-energy carbon beam at the heavy ion facility in the Shanghai Proton and Heavy Ion Center. The spatial distribution of the primary carbon beam and secondary fragments in a water phantom were systematically analyzed both in the transverse direction (perpendicular to the projection direction of the primary beam) and at different depths in the longitudinal direction (along the projection direction of the primary beam) with measured tracks on the CR-39 detectors. Meanwhile, the theoretically spatial distribution and linear energy transfer (LET) spectra of the primary beam and secondary fragments were calculated using the Monte Carlo (MC) toolkit Geant4. The results showed that the CR-39 detectors are capable of providing high lateral resolution of carbon ion at different depths. In the range of the primary carbon beam, the beam width simulated with MC is in good agreement with that of experimental measurement. The track size registered in the CR-39 has a good correlation with the particle LET. These findings indicate that the CR-39 can be used for measuring both the particle flux and its spatial distribution of carbon ions.

First theoretical determination of relative biological effectiveness of very high energy electrons

Very high energy electrons (VHEEs, E > 70 MeV) present promising clinical advantages over conventional beams due to their increased range, improved penumbra and relative insensitivity to tissue heterogeneities. They have recently garnered additional interest in their application to spatially fractionated radiotherapy or ultra-high dose rate (FLASH) therapy. However, the lack of radiobiological data limits their rapid development. This study aims to provide numerical biologically-relevant information by characterizing VHEE beams (100 and 300 MeV) against better-known beams (clinical energy electrons, photons, protons, carbon and neon ions). Their macro- and microdosimetric properties were compared, using the dose-averaged linear energy transfer ([Formula: see text]) as the macroscopic metric, and the dose-mean lineal energy [Formula: see text] and the dose-weighted lineal energy distribution, yd(y), as microscopic metrics. Finally, the modified microdosimetric kinetic model was used to calculate the respective cell survival curves and the theoretical RBE. From the macrodosimetric point of view, VHEEs presented a potential improved biological efficacy over clinical photon/electron beams due to their increased [Formula: see text]. The microdosimetric data, however, suggests no increased biological efficacy of VHEEs over clinical electron beams, resulting in RBE values of approximately 1, giving confidence to their clinical implementation. This study represents a first step to complement further radiobiological experiments.

A systematic review on the usage of averaged LET in radiation biology for particle therapy

Linear Energy Transfer (LET) is widely used to express the radiation quality of ion beams, when characterizing the biological effectiveness. However, averaged LET may be defined in multiple ways, and the chosen definition may impact the resulting reported value. We review averaged LET definitions found in the literature, and quantify which impact using these various definitions have for different reference setups. We recorded the averaged LET definitions used in 354 publications quantifying the relative biological effectiveness (RBE) of hadronic beams, and investigated how these various definitions impact the reported averaged LET using a Monte Carlo particle transport code. We find that the kind of averaged LET being applied is, generally, poorly defined. Some definitions of averaged LET may influence the reported averaged LET values up to an order of magnitude. For publications involving protons, most applied dose averaged LET when reporting RBE. The absence of what target medium is used and what secondary particles are included further contributes to an ill-defined averaged LET. We also found evidence of inconsistent usage of averaged LET definitions when deriving LET-based RBE models. To conclude, due to commonly ill-defined averaged LET and to the inherent problems of LET-based RBE models, averaged LET may only be used as a coarse indicator of radiation quality. We propose a more rigorous way of reporting LET values, and suggest that ideally the entire particle fluence spectra should be recorded and provided for future RBE studies, from which any type of averaged LET (or other quantities) may be inferred.

On calculation of the average linear energy transfer for radiobiological modelling

Applying the concept of linear energy transfer (LET) to modeling of biological effects of charged particles usually involves calculation of the average LET. To calculate this, the energy distribution of particles is characterized by either the source spectrum or fluence spectrum. Also, the average can be frequency-or dose-weighted. This makes four methods of calculating the average LET, each producing a different number. The purpose of this note is to describe which of these four methods is best suited for radiobiological modelling. We focused on data for photons (x-rays and gamma radiation) because in this case differences in the four averaging methods are most pronounced. However, our conclusions are equally applicable to photons and hadrons. We based our arguments on recently emerged Monte Carlo data that fully account for transport of electrons down to very low energies comparable to the ionization potential of water. We concluded that the frequency average LET calculated using the fluence spectrum has better predictive power than does that calculated using any of the other three options. This optimal method is not new but is different from those currently dominating research in this area.

Biological effects of passive scattering and spot scanning proton beams at the distal end of the spread-out Bragg peak in single cells and multicell spheroids

PURPOSE

The present study investigated the biological effects of spot scanning and passive scattering proton therapies at the distal end region of the spread-out Bragg peak (SOBP) using single cell and multicell spheroids.

MATERIALS AND METHODS

The Geant4 Monte Carlo simulation was used to calculate linear energy transfer (LET) values in passive scattering and spot scanning beams. The biological doses of the two beam options at various points of the distal end region of SOBP were investigated using EMT6 single cells and 0.6-mm V79 spheroids irradiated with 6 and 15 Gy, respectively, by inserting the fractions surviving these doses onto dose-survival curves and reading the corresponding dose.

RESULTS

LET values in the entrance region of SOBP were similar between the two beam options and increased at the distal end region of SOBP, where the LET value of spot scanning beams was higher than that of passive scattering beams. Increases in biological effects at the distal end region were similarly observed in single cells and spheroids; biological doses at 2-10 mm behind the distal end were 4.5-57% and 5.7-86% higher than physical doses in passive scattering and spot scanning beams, respectively, with the biological doses of spot scanning beams being higher than those of passive scattering beams (p < .05).

CONCLUSIONS

In single cells and spheroids, the effects of proton irradiation were stronger than expected from measured physical doses at the distal end of SOBP and were correlated with LET increases.

The Emerging Potential of Multi-Ion Radiotherapy

Research into high linear energy transfer (LET) radiotherapy now spans over half a century, beginning with helium and deuteron treatment in 1952 and today ranging from fast neutrons to carbon-ions. Owing to pioneering work initially in the United States and thereafter in Germany and Japan, increasing focus is on the carbon-ion beam: 12 centers are in operation, with five under construction and three in planning. While the carbon-ion beam has demonstrated unique and promising suitability in laboratory and clinical trials toward the hypofractionated treatment of hypoxic and/or radioresistant cancer, substantial developmental potential remains. Perhaps most notable is the ability to paint LET in a tumor, theoretically better focusing damage delivery within the most resistant areas. However, the technique may be limited in practice by the physical properties of the beams themselves. A heavy-ion synchrotron may provide irradiation with multiple heavy-ions: carbon, helium, and oxygen are prime candidates. Each ion varies in LET distribution, and so a methodology combining the use of multiple ions into a uniform LET distribution within a tumor may allow for even greater treatment potential in radioresistant cancer.

Simultaneous Measurements of Dose and Microdosimetric Spectra in a Clinical Proton Beam Using a scCVD Diamond Membrane Microdosimeter

A single crystal chemical vapor deposition (scCVD) diamond membrane-based microdosimetric system was used to perform simultaneous measurements of dose profile and microdosimetric spectra with the Y1 proton passive scattering beamline of the Center of Proton Therapy, Institute Curie in Orsay, France. To qualify the performance of the set-up in clinical conditions of hadrontherapy, the dose, dose rate and energy loss pulse-height spectra in a diamond microdosimeter were recorded at multiple points along depth of a water-equivalent plastic phantom. The dose-mean lineal energy (y¯D) values were computed from experimental data and compared to silicon on insulator (SOI) microdosimeter literature results. In addition, the measured dose profile, pulse height spectra and y¯D values were benchmarked with a numerical simulation using TOPAS and Geant4 toolkits. These first clinical tests of a novel system confirm that diamond is a promising candidate for a tissue equivalent, radiation hard, high spatial resolution microdosimeter in beam quality assurance of proton therapy.

Characterization of a new scintillation imaging system for proton pencil beam dose rate measurements

The goal of this work was to create a technique that could measure all possible spatial and temporal delivery rates used in pencil-beam scanning (PBS) proton therapy. The proposed system used a fast scintillation screen for full-field imaging to resolve temporal and spatial patterns as it was delivered. A fast intensified CMOS camera used continuous mode with 10 ms temporal frame rate and 1 × 1 mm2 spatial resolution, imaging a scintillation screen during clinical proton PBS delivery. PBS plans with varying dose, dose rate, energy, field size, and spot-spacing were generated, delivered and imaged. The captured images were post processed to provide dose and dose rate values after background subtraction, perspective transformation, uniformity correction for the camera and the scintillation screen, and calibration into dose. The linearity in scintillation response with respect to varying dose rate, dose, and field size was within 2%. The quenching corrected response with varying energy was also within 2%. Large spatio-temporal variations in dose rate were observed, even for plans delivered with similar dose distributions. Dose and dose rate histograms and maximum dose rate maps were generated for quantitative evaluations. With the fastest PBS delivery on a clinical system, dose rates up to 26.0 Gy s-1 were resolved. The scintillation imaging technique was able to quantify proton PBS dose rate profiles with spot weight as low as 2 MU, with spot-spacing of 2.5 mm, having a 1 × 1 mm2 spatial resolution. These dose rate temporal profiles, spatial maps, and cumulative dose rate histograms provide useful metrics for the potential evaluation and optimization of dose rate in treatment plans.

Standardizing Monte Carlo simulation parameters for a reproducible dose-averaged linear energy transfer

OBJECTIVE

Dose-averaged linear energy transfer (LETD) is one of the factors which determines relative biological effectiveness (RBE) for treatment planning in proton therapy. It is usually determined from Monte Carlo (MC) simulation. However, no standard simulation protocols were established for sampling of LETD. Simulation parameters like maximum step length and range cut will affect secondary electrons production and have an impact on the accuracy of dose distribution and LETD. We aim to show how different combinations of step length and range cut in GEANT4 will affect the result in sampling of LETD using different MC scoring methods.

METHODS

In this work, different step length and range cut value in a clinically relevant voxel geometry were used for comparison. Different LETD scoring methods were established and the concept of covariance between energy deposition per step and step length is used to explain the differences between them.

RESULTS

We recommend a maximum step length of 0.05 mm and a range cut of 0.01 mm in MC simulation as this yields the most consistent LETD value across different scoring methods. Different LETD scoring methods are also compared and variation up to 200% can be observed at the plateau of 80 MeV proton beam. Scoring Method one has one of the lowest percentage differences compared across all simulation parameters.

CONCLUSION

We have determined a set of maximum step length and range cut parameters to be used for LETD scoring in a 1 mm voxelized geometry. LETD scoring method should also be clearly defined and standardized to facilitate cross-institutional studies.

ADVANCES IN KNOWLEDGE

Establishing a standard simulation protocol for sampling LETD would reduce the discrepancy when comparing data across different centres, and this can improve the calculation for RBE.

A DNA damage multiscale model for NTCP in proton and hadron therapy

PURPOSE

To develop a first principle and multiscale model for normal tissue complication probability (NTCP) as a function of dose and LET for proton and in general for particle therapy with a goal of incorporating nanoscale radio-chemical to macroscale cell biological pathways, spanning from initial DNA damage to tissue late effects.

METHODS

The method is a combination of analytical and multiscale computational steps including (a) derivation of functional dependencies of NTCP on DNA-driven cell lethality in nanometer and mapping to dose and LET in millimeter, and (b) three-dimensional-surface fitting to Monte Carlo data set generated based on postradiation image change and gathered for a cohort of 14 pediatric patients treated by scanning beam of protons for ependymoma. We categorize voxel-based dose and LET associated with development of necrosis in NTCP.

RESULT

Our model fits well the clinical data, generated for postradiation tissue toxicity and necrosis. The fitting procedure results in extraction of in vivo radio-biological α-β indices and their numerical values.

DISCUSSION AND CONCLUSION

The NTCP model, explored in this work, allows to correlate the tissue toxicities to DNA initial damage, cell lethality and the properties and qualities of radiation, dose, and LET.

Dose- rather than fluence-averaged LET should be used as a single-parameter descriptor of proton beam quality for radiochromic film dosimetry

PURPOSE

The dose response of Gafchromic EBT3 films exposed to proton beams depends on the dose, and additionally on the beam quality, which is often quantified with the linear energy transfer (LET) and, hence, also referred to as LET quenching. Fundamentally different methods to determine correction factors for this LET quenching effect have been reported in literature and a new method using the local proton fluence distribution differential in LET is presented. This method was exploited to investigate whether a more practical correction based on the dose- or fluence-averaged LET is feasible in a variety of clinically possible beam arrangements.

METHODS

The relative effectiveness (RE) was characterized within a high LET spread-out Bragg peak (SOBP) in water made up by the six lowest available energies (62.4-67.5 MeV, configuration " b 1 ") resulting in one of the highest clinically feasible dose-averaged LET distributions. Additionally, two beams were measured where a low LET proton beam (252.7 MeV) was superimposed on " b 1 ", which contributed either 50% of the initial particle fluence or 50% of the dose in the SOBP, referred to as configuration " b 2 " and " b 3 ," respectively. The proton LET spectrum was simulated with GATE/Geant4 at all measurement positions. The net optical density change differential in LET was integrated over the local proton spectrum to calculate the net optical density and therefrom the beam quality correction factor. The LET dependence of the film response was accounted for by an LET dependence of one of the three parameters in the calibration function and was determined from inverse optimization using measurement " b 1 ." This method was then validated on the measurements of " b 2 " and " b 3 " and subsequently used to calculate the RE at 900 positions in nine clinically relevant beams. The extrapolated RE set was used to derive a simple linear correction function based on dose-averaged LET ( L d ) and verify the validity in all points of the comprehensive RE set.

RESULTS

The uncorrected film dose deviated up to 26% from the reference dose, whereas the corrected film dose agreed within 3% in all three beams in water (" b 1 ", " b 2 " and " b 3 "). The LET dependence of the calibration function started to strongly increase around 5 keV/μm and flatten out around 30 keV/μm. All REs calculated from the proton fluence in the nine simulated beams could be approximated with a linear function of dose-averaged LET (RE = 1.0258-0.0211 μm/keV L d ). However, no functional relationship of RE- and fluence-averaged LET could be found encompassing all beam energies and modulations.

CONCLUSIONS

The film quenching was found to be nonlinear as a function of proton LET as well as of the dose-averaged LET. However, the linear relation of RE on dose-averaged LET was a good approximation in all cases. In contrast to dose-averaged LET, fluence-averaged LET could not describe the RE when multiple beams were applied.

Dose-average linear energy transfer of electrons released in liquid water and LiF:Mg,Ti by low-energy x-rays, 137Cs and 60Co gamma

For many years, track-average linear energy transfer (LET), [Formula: see text] has been used to quantify the radiation-induced phenomena in biological and physical systems. However, due to the need for including into the radiotherapy treatment planning system, parameters that are clinical and biologically relevant, a precise knowledge of the dose-average LET, [Formula: see text] becomes essential. Besides, several dosimetric studies have revealed that [Formula: see text] is fundamental to describe the dosimeter's response induced by photons. The most important data sets publicly available for [Formula: see text] of electron generated by photons are those reported for measurements performed in methane-based tissue-equivalent gas. However, comparing to liquid water, the electron spectra generated by low photon energy might not be similar due to the photoelectric effect. Thus, this work aimed at investigating the [Formula: see text] of electron spectra generated in liquid water and LiF:Mg,Ti by ten x-ray beams from 20 kV to 300 kV, ¹³⁷Cs and ⁶⁰Co gamma. The results suggest that [Formula: see text] is more sensitive to the surrounding environment than [Formula: see text] and consequently, it might be a more appropriate parameter to quantify the radiation effect and damage in matter induced by photons. Besides, good agreement (6% to 12% differences versus 10% to 15% uncertainties in the experiments) was observed between the data obtained in this work for liquid water and the experimental values published for methane-based tissue-equivalent gas at energies above 60 keV. Whereas at lowest energies, the minimum difference is around 18% which can be associated to the difference between the two media.

Physical advantages of particles: protons and light ions

Proton and ion beam therapy has been introduced in the Lawrence Berkeley National Laboratory in the mid-1950s, when protons and helium ions have been used for the first time to treat patients. Starting in 1972, the scientists at Berkeley also were the first to use heavier ions (carbon, oxygen, neon, silicon and argon ions). The first clinical ion beam facility opened in 1994 in Japan and since then, the interest in radiotherapy with light ion beams has been increasing slowly but steadily, with 13 centers in clinical operation in 2019. All these centers are using carbon ions for clinical application.The article outlines the differences in physical properties of various light ions as compared to protons in view of the application in radiotherapy. These include the energy loss and depth dose properties, multiple scattering, range straggling and nuclear fragmentation. In addition, the paper discusses differences arising from energy loss and linear energy transfer with respect to their biological effects.Moreover, the paper reviews briefly the existing clinical data comparing protons and ions and outlines the future perspectives for the clinical use of ions like oxygen and helium.

Exploring the advantages of intensity-modulated proton therapy: experimental validation of biological effects using two different beam intensity-modulation patterns

In current treatment plans of intensity-modulated proton therapy, high-energy beams are usually assigned larger weights than low-energy beams. Using this form of beam delivery strategy cannot effectively use the biological advantages of low-energy and high-linear energy transfer (LET) protons present within the Bragg peak. However, the planning optimizer can be adjusted to alter the intensity of each beamlet, thus maintaining an identical target dose while increasing the weights of low-energy beams to elevate the LET therein. The objective of this study was to experimentally validate the enhanced biological effects using a novel beam delivery strategy with elevated LET. We used Monte Carlo and optimization algorithms to generate two different intensity-modulation patterns, namely to form a downslope and a flat dose field in the target. We spatially mapped the biological effects using high-content automated assays by employing an upgraded biophysical system with improved accuracy and precision of collected data. In vitro results in cancer cells show that using two opposed downslope fields results in a more biologically effective dose, which may have the clinical potential to increase the therapeutic index of proton therapy.