TOPAS simulation of photoneutrons in radiotherapy: accuracy and speed with variance reduction.

Objective. We provide optimal particle split numbers for speeding up TOPAS Monte Carlo simulations of linear accelerator (linac) treatment heads while maintaining accuracy. In addition, we provide a new TOPAS physics module for simulating photoneutron production and transport.Approach.TOPAS simulation of a Siemens Oncor linac was used to determine the optimal number of splits for directional bremsstrahlung splitting as a function of the field size for 6 MV and 18 MV x-ray beams. The linac simulation was validated against published data of lateral dose profiles and percentage depth-dose curves (PDD) for the largest square field (40 cm side). In separate simulations, neutron particle split and the custom TOPAS physics module was used to generate and transport photoneutrons, called 'TsPhotoNeutron'. Verification of accuracy was performed by comparing simulations with published measurements of: (1) neutron yields as a function of beam energy for thick targets of Al, Cu, Ta, W, Pb and concrete; and (2) photoneutron energy spectrum at 40 cm laterally from the isocenter of the Oncor linac from an 18 MV beam with closed jaws and MLC.Main results.The optimal number of splits obtained for directional bremsstrahlung splitting enhanced the computational efficiency by two orders of magnitude. The efficiency decreased with increasing beam energy and field size. Calculated lateral profiles in the central region agreed within 1 mm/2% from measured data, PDD curves within 1 mm/1%. For the TOPAS physics module, at a split number of 146, the efficiency of computing photoneutron yields was enhanced by a factor of 27.6, whereas it improved the accuracy over existing Geant4 physics modules.Significance.This work provides simulation parameters and a new TOPAS physics module to improve the efficiency and accuracy of TOPAS simulations that involve photonuclear processes occurring in high-Zmaterials found in linac components, patient devices, and treatment rooms, as well as to explore new therapeutic modalities such as very-high energy electron therapy.

Spatiotemporal optimisation of prostate intensity modulated proton therapy (IMPT) treatments.

Objective.In intensity modulated particle therapy (IMPT), the adoption of spatially and temporally heterogeneous dose distributions allows to decouple the fractionation scheme from the patient anatomy, so that an hypofractionated schedule can be selectively created inside the tumour, while simultaneously exploiting the fractionation effect in the healthy tissues. In this paper, the authors show the reproducibility of the method on a set of prostate patients, quantifying the dependencies of the achievable benefit with respect to conventional and hypofractionated schemes and the sensitivity of the method to setup errors and range uncertainty.Approach.On a cohort of 9 patients, non-uniform IMPT plans were optimised and compared to conventional and hypofractionated schedules. For each patient, the comparison of the three strategies has been based on the output of the cost function used to optimise the treatments. The analysis has been repeated considering differentα/ßratios for the tumour, namely 1.5, 3 and 4.5 Gy. For a single patient, setup errors and beam range uncertainty have been analysed: the plans, for each optimisation strategy, have been iteratively forward planned 500 times with randomly varying the patient position in each fraction, and 200 times for systematic range shift.Main results.An average 10% benefit has been shown for the lowestα/ßratio considered for the tumour, where the non-uniform schedule generally converges to hypofractionation; the benefit decreases to 5%-7% for higherα/ßratios, for which the non-uniform schedule always showed better outcomes with respect to the other fractionation schedules. An increased sensitivity to uncertainty, especially for setup errors, has been shown, which can be associated to the spatial non-uniformity of the dose distributions peculiar of the spatiotemporal plans.Significance.This work represents the first investigation of spatiotemporal fractionation for prostate cancer and the beginning of further investigations before clinical implementation can be considered.

TOPAS-nBio validation for simulating water radiolysis and DNA damage under low-LET irradiation.

The chemical stage of the Monte Carlo track-structure simulation code Geant4-DNA has been revised and validated. The root-mean-square (RMS) empirical parameter that dictates the displacement of water molecules after an ionization and excitation event in Geant4-DNA has been shortened to better fit experimental data. The pre-defined dissociation channels and branching ratios were not modified, but the reaction rate coefficients for simulating the chemical stage of water radiolysis were updated. The evaluation of Geant4-DNA was accomplished with TOPAS-nBio. For that, we compared predicted time-dependentGvalues in pure liquid water for·OH, e-aq, and H2with published experimental data. For H2O2and H·, simulation of added scavengers at different concentrations resulted in better agreement with measurements. In addition, DNA geometry information was integrated with chemistry simulation in TOPAS-nBio to realize reactions between radiolytic chemical species and DNA. This was used in the estimation of the yield of single-strand breaks (SSB) induced by137Csγ-ray radiolysis of supercoiled pUC18 plasmids dissolved in aerated solutions containing DMSO. The efficiency of SSB induction by reaction between radiolytic species and DNA used in the simulation was chosen to provide the best agreement with published measurements. An RMS displacement of 1.24 nm provided agreement with measured data within experimental uncertainties for time-dependentGvalues and under the presence of scavengers. SSB efficiencies of 24% and 0.5% for·OH and H·, respectively, led to an overall agreement of TOPAS-nBio results within experimental uncertainties. The efficiencies obtained agreed with values obtained with published non-homogeneous kinetic model and step-by-step Monte Carlo simulations but disagreed by 12% with published direct measurements. Improvement of the spatial resolution of the DNA damage model might mitigate such disagreement. In conclusion, with these improvements, Geant4-DNA/TOPAS-nBio provides a fast, accurate, and user-friendly tool for simulating DNA damage under low linear energy transfer irradiation.

TOPAS-nBio: An Extension to the TOPAS Simulation Toolkit for Cellular and Sub-cellular Radiobiology.

The TOPAS Monte Carlo (MC) system is used in radiation therapy and medical imaging research, having played a significant role in making Monte Carlo simulations widely available for proton therapy related research. While TOPAS provides detailed simulations of patient scale properties, the fundamental unit of the biological response to radiation is a cell. Thus, our goal was to develop TOPAS-nBio, an extension of TOPAS dedicated to advance understanding of radiobiological effects at the (sub-)cellular, (i.e., the cellular and sub-cellular) scale. TOPAS-nBio was designed as a set of open source classes that extends TOPAS to model radiobiological experiments. TOPAS-nBio is based on and extends Geant4-DNA, which extends the Geant4 toolkit, the basis of TOPAS, to include very low-energy interactions of particles down to vibrational energies, explicitly simulates every particle interaction (i.e., without using condensed histories) and propagates radiolysis products. To further facilitate the use of TOPAS-nBio, a graphical user interface was developed. TOPAS-nBio offers full track-structure Monte Carlo simulations, integration of chemical reactions within the first millisecond, an extensive catalogue of specialized cell geometries as well as sub-cellular structures such as DNA and mitochondria, and interfaces to mechanistic models of DNA repair kinetics. We compared TOPAS-nBio simulations to measured and published data of energy deposition patterns and chemical reaction rates (G values). Our simulations agreed well within the experimental uncertainties. Additionally, we expanded the chemical reactions and species provided in Geant4-DNA and developed a new method based on independent reaction times (IRT), including a total of 72 reactions classified into 6 types between neutral and charged species. Chemical stage simulations using IRT were a factor of 145 faster than with step-by-step tracking. Finally, we applied the geometric/chemical modeling to obtain initial yields of double-strand breaks (DSBs) in DNA fibers for proton irradiations of 3 and 50 MeV and compared the effect of including chemical reactions on the number and complexity of DSB induction. Over half of the DSBs were found to include chemical reactions with approximately 5% of DSBs caused only by chemical reactions. In conclusion, the TOPAS-nBio extension to the TOPAS MC application offers access to accurate and detailed multiscale simulations, from a macroscopic description of the radiation field to microscopic description of biological outcome for selected cells. TOPAS-nBio offers detailed physics and chemistry simulations of radiobiological experiments on cells simulating the initially induced damage and links to models of DNA repair kinetics.

A New Standard DNA Damage (SDD) Data Format.

Our understanding of radiation-induced cellular damage has greatly improved over the past few decades. Despite this progress, there are still many obstacles to fully understand how radiation interacts with biologically relevant cellular components, such as DNA, to cause observable end points such as cell killing. Damage in DNA is identified as a major route of cell killing. One hurdle when modeling biological effects is the difficulty in directly comparing results generated by members of different research groups. Multiple Monte Carlo codes have been developed to simulate damage induction at the DNA scale, while at the same time various groups have developed models that describe DNA repair processes with varying levels of detail. These repair models are intrinsically linked to the damage model employed in their development, making it difficult to disentangle systematic effects in either part of the modeling chain. These modeling chains typically consist of track-structure Monte Carlo simulations of the physical interactions creating direct damages to DNA, followed by simulations of the production and initial reactions of chemical species causing so-called "indirect" damages. After the induction of DNA damage, DNA repair models combine the simulated damage patterns with biological models to determine the biological consequences of the damage. To date, the effect of the environment, such as molecular oxygen (normoxic vs. hypoxic), has been poorly considered. We propose a new standard DNA damage (SDD) data format to unify the interface between the simulation of damage induction in DNA and the biological modeling of DNA repair processes, and introduce the effect of the environment (molecular oxygen or other compounds) as a flexible parameter. Such a standard greatly facilitates inter-model comparisons, providing an ideal environment to tease out model assumptions and identify persistent, underlying mechanisms. Through inter-model comparisons, this unified standard has the potential to greatly advance our understanding of the underlying mechanisms of radiation-induced DNA damage and the resulting observable biological effects when radiation parameters and/or environmental conditions change.

Online adaption approaches for intensity modulated proton therapy for head and neck patients based on cone beam CTs and Monte Carlo simulations.

To develop an online plan adaptation algorithm for intensity modulated proton therapy (IMPT) based on fast Monte Carlo dose calculation and cone beam CT (CBCT) imaging. A cohort of ten head and neck cancer patients with an average of six CBCT scans were studied. To adapt the treatment plan to the new patient geometry, contours were propagated to the CBCTs with a vector field (VF) calculated with deformable image registration between the CT and the CBCTs. Within the adaptive planning algorithm, beamlets were shifted following the VF at their distal falloff and raytraced in the CBCT to adjust their energies, creating a geometrically adapted plan. Four geometric adaptation modes were studied: unconstrained geometric shifts (Free), isocenter shift (Iso), a range shifter (RS), or isocenter shift and range shifter (Iso-RS). After evaluation of the geometrical adaptation, the weights of a selected subset of beamlets were automatically tuned using MC-generated influence matrices to fulfill the original plan requirements. All beamlet calculations were done with a fast Monte Carlo running on a GPU (graphics processing unit). Geometrical adaptation alone only worked with small anatomy changes. The weight-tuned adaptation worked for every fraction, with the Free and Iso modes performing similarly and being superior than the two range shifters modes. The mean V95 and V107 were 99.4 ± 0.9 and 6.4% ± 4.7% in the Free mode with weight tuning. The calculation time per fraction was ~5 min, but further task parallelization could reduce it to ~1-2 min for delivery adaptation right after patient setup. An online adaptation algorithm was developed that significantly improved the treatment quality for inter-fractional geometry changes. Clinical implementation of the algorithm would allow delivery adaptation right before treatment and thus allow planning margin reductions for IMPT.

Monte Carlo simulation of chemistry following radiolysis with TOPAS-nBio.

Simulation of water radiolysis and the subsequent chemistry provides important information on the effect of ionizing radiation on biological material. The Geant4 Monte Carlo toolkit has added chemical processes via the Geant4-DNA project. The TOPAS tool simplifies the modeling of complex radiotherapy applications with Geant4 without requiring advanced computational skills, extending the pool of users. Thus, a new extension to TOPAS, TOPAS-nBio, is under development to facilitate the configuration of track-structure simulations as well as water radiolysis simulations with Geant4-DNA for radiobiological studies. In this work, radiolysis simulations were implemented in TOPAS-nBio. Users may now easily add chemical species and their reactions, and set parameters including branching ratios, dissociation schemes, diffusion coefficients, and reaction rates. In addition, parameters for the chemical stage were re-evaluated and updated from those used by default in Geant4-DNA to improve the accuracy of chemical yields. Simulation results of time-dependent and LET-dependent primary yields Gx (chemical species per 100 eV deposited) produced at neutral pH and 25 °C by short track-segments of charged particles were compared to published measurements. The LET range was 0.05-230 keV µm-1. The calculated Gx values for electrons satisfied the material balance equation within 0.3%, similar for protons albeit with long calculation time. A smaller geometry was used to speed up proton and alpha simulations, with an acceptable difference in the balance equation of 1.3%. Available experimental data of time-dependent G-values for [Formula: see text] agreed with simulated results within 7% ± 8% over the entire time range; for [Formula: see text] over the full time range within 3% ± 4%; for H2O2 from 49% ± 7% at earliest stages and 3% ± 12% at saturation. For the LET-dependent Gx, the mean ratios to the experimental data were 1.11 ± 0.98, 1.21 ± 1.11, 1.05 ± 0.52, 1.23 ± 0.59 and 1.49 ± 0.63 (1 standard deviation) for [Formula: see text], [Formula: see text], H2, H2O2 and [Formula: see text], respectively. In conclusion, radiolysis and subsequent chemistry with Geant4-DNA has been successfully incorporated in TOPAS-nBio. Results are in reasonable agreement with published measured and simulated data.

Experimental assessment of proton dose calculation accuracy in inhomogeneous media.

PURPOSE: Proton therapy with Pencil Beam Scanning (PBS) has the potential to improve radiotherapy treatments. Unfortunately, its promises are jeopardized by the sensitivity of the dose distributions to uncertainties, including dose calculation accuracy in inhomogeneous media. Monte Carlo dose engines (MC) are expected to handle heterogeneities better than analytical algorithms like the pencil-beam convolution algorithm (PBA). In this study, an experimental phantom has been devised to maximize the effect of heterogeneities and to quantify the capability of several dose engines (MC and PBA) to handle these. METHODS: An inhomogeneous phantom made of water surrounding a long insert of bone tissue substitute (1×10×10 cm3) was irradiated with a mono-energetic PBS field (10×10 cm2). A 2D ion chamber array (MatriXX, IBA Dosimetry GmbH) lied right behind the bone. The beam energy was such that the expected range of the protons exceeded the detector position in water and did not attain it in bone. The measurement was compared to the following engines: Geant4.9.5, PENH, MCsquare, as well as the MC and PBA algorithms of RayStation (RaySearch Laboratories AB). RESULTS: For a γ-index criteria of 2%/2mm, the passing rates are 93.8% for Geant4.9.5, 97.4% for PENH, 93.4% for MCsquare, 95.9% for RayStation MC, and 44.7% for PBA. The differences in γ-index passing rates between MC and RayStation PBA calculations can exceed 50%. CONCLUSION: The performance of dose calculation algorithms in highly inhomogeneous media was evaluated in a dedicated experiment. MC dose engines performed overall satisfactorily while large deviations were observed with PBA as expected.

Comparing stochastic proton interactions simulated using TOPAS-nBio to experimental data from fluorescent nuclear track detectors.

Whilst Monte Carlo (MC) simulations of proton energy deposition have been well-validated at the macroscopic level, their microscopic validation remains lacking. Equally, no gold-standard yet exists for experimental metrology of individual proton tracks. In this work we compare the distributions of stochastic proton interactions simulated using the TOPAS-nBio MC platform against confocal microscope data for Al2O3:C,Mg fluorescent nuclear track detectors (FNTDs). We irradiated [Formula: see text] mm3 FNTD chips inside a water phantom, positioned at seven positions along a pristine proton Bragg peak with a range in water of 12 cm. MC simulations were implemented in two stages: (1) using TOPAS to model the beam properties within a water phantom and (2) using TOPAS-nBio with Geant4-DNA physics to score particle interactions through a water surrogate of Al2O3:C,Mg. The measured median track integrated brightness (IB) was observed to be strongly correlated to both (i) voxelized track-averaged linear energy transfer (LET) and (ii) frequency mean microdosimetric lineal energy, [Formula: see text], both simulated in pure water. Histograms of FNTD track IB were compared against TOPAS-nBio histograms of the number of terminal electrons per proton, scored in water with mass-density scaled to mimic Al2O3:C,Mg. Trends between exposure depths observed in TOPAS-nBio simulations were experimentally replicated in the study of FNTD track IB. Our results represent an important first step towards the experimental validation of MC simulations on the sub-cellular scale and suggest that FNTDs can enable experimental study of the microdosimetric properties of individual proton tracks.

Methodologies in the modeling of combined chemo-radiation treatments.

The variety of treatment options for cancer patients has increased significantly in recent years. Not only do we combine radiation with surgery and chemotherapy, new therapeutic approaches such as immunotherapy and targeted therapies are starting to play a bigger role. Physics has made significant contributions to radiation therapy treatment planning and delivery. In particular, treatment plan optimization using inverse planning techniques has improved dose conformity considerably. Furthermore, medical physics is often the driving force behind tumor control and normal tissue complication modeling. While treatment optimization and outcome modeling does focus mainly on the effects of radiation, treatment modalities such as chemotherapy are treated independently or are even neglected entirely. This review summarizes the published efforts to model combined modality treatments combining radiation and chemotherapy. These models will play an increasing role in optimizing cancer therapy not only from a radiation and drug dosage standpoint, but also in terms of spatial and temporal optimization of treatment schedules.

Dose enhancement effects to the nucleus and mitochondria from gold nanoparticles in the cytosol.

Gold nanoparticles (GNPs) have shown potential as dose enhancers for radiation therapy. Since damage to the genome affects the viability of a cell, it is generally assumed that GNPs have to localise within the cell nucleus. In practice, however, GNPs tend to localise in the cytoplasm yet still appear to have a dose enhancing effect on the cell. Whether this effect can be attributed to stress-induced biological mechanisms or to physical damage to extra-nuclear cellular targets is still unclear. There is however growing evidence to suggest that the cellular response to radiation can also be influenced by indirect processes induced when the nucleus is not directly targeted by radiation. The mitochondrion in particular may be an effective extra-nuclear radiation target given its many important functional roles in the cell. To more accurately predict the physical effect of radiation within different cell organelles, we measured the full chemical composition of a whole human lymphocytic JURKAT cell as well as two separate organelles; the cell nucleus and the mitochondrion. The experimental measurements found that all three biological materials had similar ionisation energies ∼70 eV, substantially lower than that of liquid water ∼78 eV. Monte Carlo simulations for 10-50 keV incident photons showed higher energy deposition and ionisation numbers in the cell and organelle materials compared to liquid water. Adding a 1% mass fraction of gold to each material increased the energy deposition by a factor of ∼1.8 when averaged over all incident photon energies. Simulations of a realistic compartmentalised cell show that the presence of gold in the cytosol increases the energy deposition in the mitochondrial volume more than within the nuclear volume. We find this is due to sub-micron delocalisation of energy by photoelectrons, making the mitochondria a potentially viable indirect radiation target for GNPs that localise to the cytosol.

Radiological protection in ion beam radiotherapy: practical guidance for clinical use of new technology.

Recently introduced technologies in radiotherapy have significantly improved the clinical outcome for patients. Ion beam radiotherapy, involving proton and carbon ion beams, provides excellent dose distributions in targeted tumours, with reduced doses to the surrounding normal tissues. However, careful treatment planning is required in order to maximise the treatment efficiency and minimise the dose to normal tissues. Radiation exposure from secondary neutrons and photons, particle fragments, and photons from activated materials should also be considered for radiological protection of the patient and medical staff. Appropriate maintenance is needed for the equipment and air in the treatment room, which may be activated by the particle beam and its secondary radiation. This new treatment requires complex procedures and careful adjustment of parameters for each patient. Therefore, education and training for the personnel involved in the procedure are essential for both effective treatment and patient protection. The International Commission on Radiological Protection (ICRP) has provided recommendations for radiological protection in ion beam radiotherapy in Publication 127 Medical staff should be aware of the possible risks resulting from inappropriate use and control of the equipment. They should also consider the necessary procedures for patient protection when new technologies are introduced into clinical practice.

Improved efficiency in Monte Carlo simulation for passive-scattering proton therapy.

The aim of this work was to improve the computational efficiency of Monte Carlo simulations when tracking protons through a proton therapy treatment head. Two proton therapy facilities were considered, the Francis H Burr Proton Therapy Center (FHBPTC) at the Massachusetts General Hospital and the Crocker Lab eye treatment facility used by University of California at San Francisco (UCSFETF). The computational efficiency was evaluated for phase space files scored at the exit of the treatment head to determine optimal parameters to improve efficiency while maintaining accuracy in the dose calculation. For FHBPTC, particles were split by a factor of 8 upstream of the second scatterer and upstream of the aperture. The radius of the region for Russian roulette was set to 2.5 or 1.5 times the radius of the aperture and a secondary particle production cut (PC) of 50 mm was applied. For UCSFETF, particles were split a factor of 16 upstream of a water absorber column and upstream of the aperture. Here, the radius of the region for Russian roulette was set to 4 times the radius of the aperture and a PC of 0.05 mm was applied. In both setups, the cylindrical symmetry of the proton beam was exploited to position the split particles randomly spaced around the beam axis. When simulating a phase space for subsequent water phantom simulations, efficiency gains between a factor of 19.9 ± 0.1 and 52.21 ± 0.04 for the FHTPC setups and 57.3 ± 0.5 for the UCSFETF setups were obtained. For a phase space used as input for simulations in a patient geometry, the gain was a factor of 78.6 ± 7.5. Lateral-dose curves in water were within the accepted clinical tolerance of 2%, with statistical uncertainties of 0.5% for the two facilities. For the patient geometry and by considering the 2% and 2mm criteria, 98.4% of the voxels showed a gamma index lower than unity. An analysis of the dose distribution resulted in systematic deviations below of 0.88% for 20% of the voxels with dose of 20% of the maximum or more.

A framework for implementation of organ effect models in TOPAS with benchmarks extended to proton therapy.

The aim of this work was to develop a framework for modeling organ effects within TOPAS (TOol for PArticle Simulation), a wrapper of the Geant4 Monte Carlo toolkit that facilitates particle therapy simulation. The DICOM interface for TOPAS was extended to permit contour input, used to assign voxels to organs. The following dose response models were implemented: The Lyman-Kutcher-Burman model, the critical element model, the population based critical volume model, the parallel-serial model, a sigmoid-based model of Niemierko for normal tissue complication probability and tumor control probability (TCP), and a Poisson-based model for TCP. The framework allows easy manipulation of the parameters of these models and the implementation of other models. As part of the verification, results for the parallel-serial and Poisson model for x-ray irradiation of a water phantom were compared to data from the AAPM Task Group 166. When using the task group dose-volume histograms (DVHs), results were found to be sensitive to the number of points in the DVH, with differences up to 2.4%, some of which are attributable to differences between the implemented models. New results are given with the point spacing specified. When using Monte Carlo calculations with TOPAS, despite the relatively good match to the published DVH's, differences up to 9% were found for the parallel-serial model (for a maximum DVH difference of 2%) and up to 0.5% for the Poisson model (for a maximum DVH difference of 0.5%). However, differences of 74.5% (in Rectangle1), 34.8% (in PTV) and 52.1% (in Triangle) for the critical element, critical volume and the sigmoid-based models were found respectively. We propose a new benchmark for verification of organ effect models in proton therapy. The benchmark consists of customized structures in the spread out Bragg peak plateau, normal tissue, tumor, penumbra and in the distal region. The DVH's, DVH point spacing, and results of the organ effect models are provided. The models were used to calculate dose response for a Head and Neck patient to demonstrate functionality of the new framework and indicate the degree of variability between the models in proton therapy.

Monitoring proton therapy with PET.

Protons are being used in radiation therapy because of typically better dose conformity and reduced total energy deposited in the patient as compared with photon techniques. Both aspects are related to the finite range of a proton beam. The finite range also allows advanced dose shaping. These benefits can only be fully utilized if the end of range can be predicted accurately in the patient. The prediction of the range in tissue is associated with considerable uncertainties owing to imaging, patient set-up, beam delivery, interfractional changes in patient anatomy and dose calculation. Consequently, a significant range (of the order of several millimetres) is added to the prescribed range in order to ensure tumour coverage. Thus, reducing range uncertainties would allow a reduction of the treatment volume and reduce dose to potential organs at risk.

Assessment of uncertainties in radiation-induced cancer risk predictions at clinically relevant doses.

PURPOSE: Theoretical dose-response models offer the possibility to assess second cancer induction risks after external beam therapy. The parameters used in these models are determined with limited data from epidemiological studies. Risk estimations are thus associated with considerable uncertainties. This study aims at illustrating uncertainties when predicting the risk for organ-specific second cancers in the primary radiation field illustrated by choosing selected treatment plans for brain cancer patients. METHODS: A widely used risk model was considered in this study. The uncertainties of the model parameters were estimated with reported data of second cancer incidences for various organs. Standard error propagation was then subsequently applied to assess the uncertainty in the risk model. Next, second cancer risks of five pediatric patients treated for cancer in the head and neck regions were calculated. For each case, treatment plans for proton and photon therapy were designed to estimate the uncertainties (a) in the lifetime attributable risk (LAR) for a given treatment modality and (b) when comparing risks of two different treatment modalities. RESULTS: Uncertainties in excess of 100% of the risk were found for almost all organs considered. When applied to treatment plans, the calculated LAR values have uncertainties of the same magnitude. A comparison between cancer risks of different treatment modalities, however, does allow statistically significant conclusions. In the studied cases, the patient averaged LAR ratio of proton and photon treatments was 0.35, 0.56, and 0.59 for brain carcinoma, brain sarcoma, and bone sarcoma, respectively. Their corresponding uncertainties were estimated to be potentially below 5%, depending on uncertainties in dosimetry. CONCLUSIONS: The uncertainty in the dose-response curve in cancer risk models makes it currently impractical to predict the risk for an individual external beam treatment. On the other hand, the ratio of absolute risks between two modalities is less sensitive to the uncertainties in the risk model and can provide statistically significant estimates.

Characterizing a proton beam scanning system for Monte Carlo dose calculation in patients.

The presented work has two goals. First, to demonstrate the feasibility of accurately characterizing a proton radiation field at treatment head exit for Monte Carlo dose calculation of active scanning patient treatments. Second, to show that this characterization can be done based on measured depth dose curves and spot size alone, without consideration of the exact treatment head delivery system. This is demonstrated through calibration of a Monte Carlo code to the specific beam lines of two institutions, Massachusetts General Hospital (MGH) and Paul Scherrer Institute (PSI). Comparison of simulations modeling the full treatment head at MGH to ones employing a parameterized phase space of protons at treatment head exit reveals the adequacy of the method for patient simulations. The secondary particle production in the treatment head is typically below 0.2% of primary fluence, except for low-energy electrons (<0.6 MeV for 230 MeV protons), whose contribution to skin dose is negligible. However, there is significant difference between the two methods in the low-dose penumbra, making full treatment head simulations necessary to study out-of-field effects such as secondary cancer induction. To calibrate the Monte Carlo code to measurements in a water phantom, we use an analytical Bragg peak model to extract the range-dependent energy spread at the two institutions, as this quantity is usually not available through measurements. Comparison of the measured with the simulated depth dose curves demonstrates agreement within 0.5 mm over the entire energy range. Subsequently, we simulate three patient treatments with varying anatomical complexity (liver, head and neck and lung) to give an example how this approach can be employed to investigate site-specific discrepancies between treatment planning system and Monte Carlo simulations.

Dosimetric feasibility of real-time MRI-guided proton therapy.

PURPOSE: Magnetic resonance imaging (MRI) is a prime candidate for image-guided radiotherapy. This study was designed to assess the feasibility of real-time MRI-guided proton therapy by quantifying the dosimetric effects induced by the magnetic field in patients' plans and identifying the associated clinical consequences. METHODS: Monte Carlo dose calculation was performed for nine patients of various treatment sites (lung, liver, prostate, brain, skull-base, and spine) and tissue homogeneities, in the presence of 0.5 and 1.5 T magnetic fields. Dose volume histogram (DVH) parameters such as D95, D5, and V20 as well as equivalent uniform dose were compared for the target and organs at risk, before and after applying the magnetic field. The authors further assessed whether the plans affected by clinically relevant dose distortions could be corrected independent of the planning system. RESULTS: By comparing the resulting dose distributions and analyzing the respective DVHs, it was determined that despite the observed lateral beam deflection, for magnetic fields of up to 0.5 T, neither was the target coverage jeopardized nor was the dose to the nearby organs increased in all cases except for prostate. However, for a 1.5 T magnetic field, the dose distortions were more pronounced and of clinical concern in all cases except for spine. In such circumstances, the target was severely underdosed, as indicated by a decrease in D95 of up to 41% of the prescribed dose compared to the nominal situation (no magnetic field). Sites such as liver and spine were less affected due to higher tissue homogeneity, typically smaller beam range, and the choice of beam directions. Simulations revealed that small modifications to certain plan parameters such as beam isocenter (up to 19 mm) and gantry angle (up to 10°) are sufficient to compensate for the magnetic field-induced dose disturbances. The authors' observations indicate that the degree of required corrections strongly depends on the beam range and direction relative to the magnetic field. This method was also applicable to more heterogeneous scenarios such as skull-base tumors. CONCLUSIONS: This study confirmed the dosimetric feasibility of real-time MRI-guided proton therapy and delivering a clinically acceptable dose to patients with various tumor locations within magnetic fields of up to 1.5 T. This work could serve as a guide and encouragement for further efforts toward clinical implementation of hybrid MRI-proton gantry systems.

Automation and uncertainty analysis of a method for in-vivo range verification in particle therapy.

We introduce the automation of the range difference calculation deduced from particle-irradiation induced ß(+)-activity distributions with the so-called most-likely-shift approach, and evaluate its reliability via the monitoring of algorithm- and patient-specific uncertainty factors. The calculation of the range deviation is based on the minimization of the absolute profile differences in the distal part of two activity depth profiles shifted against each other. Depending on the workflow of positron emission tomography (PET)-based range verification, the two profiles under evaluation can correspond to measured and simulated distributions, or only measured data from different treatment sessions. In comparison to previous work, the proposed approach includes an automated identification of the distal region of interest for each pair of PET depth profiles and under consideration of the planned dose distribution, resulting in the optimal shift distance. Moreover, it introduces an estimate of uncertainty associated to the identified shift, which is then used as weighting factor to 'red flag' problematic large range differences. Furthermore, additional patient-specific uncertainty factors are calculated using available computed tomography (CT) data to support the range analysis. The performance of the new method for in-vivo treatment verification in the clinical routine is investigated with in-room PET images for proton therapy as well as with offline PET images for proton and carbon ion therapy. The comparison between measured PET activity distributions and predictions obtained by Monte Carlo simulations or measurements from previous treatment fractions is performed. For this purpose, a total of 15 patient datasets were analyzed, which were acquired at Massachusetts General Hospital and Heidelberg Ion-Beam Therapy Center with in-room PET and offline PET/CT scanners, respectively. Calculated range differences between the compared activity distributions are reported in a 2D map in beam-eye-view. In comparison to previously proposed approaches, the new most-likely-shift method shows more robust results for assessing in-vivo the range from strongly varying PET distributions caused by differing patient geometry, ion beam species, beam delivery techniques, PET imaging concepts and counting statistics. The additional visualization of the uncertainties and the dedicated weighting strategy contribute to the understanding of the reliability of observed range differences and the complexity in the prediction of activity distributions. The proposed method promises to offer a feasible technique for clinical routine of PET-based range verification.

Site-specific range uncertainties caused by dose calculation algorithms for proton therapy.

The purpose of this study was to assess the possibility of introducing site-specific range margins to replace current generic margins in proton therapy. Further, the goal was to study the potential of reducing margins with current analytical dose calculations methods. For this purpose we investigate the impact of complex patient geometries on the capability of analytical dose calculation algorithms to accurately predict the range of proton fields. Dose distributions predicted by an analytical pencil-beam algorithm were compared with those obtained using Monte Carlo (MC) simulations (TOPAS). A total of 508 passively scattered treatment fields were analyzed for seven disease sites (liver, prostate, breast, medulloblastoma-spine, medulloblastoma-whole brain, lung and head and neck). Voxel-by-voxel comparisons were performed on two-dimensional distal dose surfaces calculated by pencil-beam and MC algorithms to obtain the average range differences and root mean square deviation for each field for the distal position of the 90% dose level (R90) and the 50% dose level (R50). The average dose degradation of the distal falloff region, defined as the distance between the distal position of the 80% and 20% dose levels (R80-R20), was also analyzed. All ranges were calculated in water-equivalent distances. Considering total range uncertainties and uncertainties from dose calculation alone, we were able to deduce site-specific estimations. For liver, prostate and whole brain fields our results demonstrate that a reduction of currently used uncertainty margins is feasible even without introducing MC dose calculations. We recommend range margins of 2.8% + 1.2 mm for liver and prostate treatments and 3.1% + 1.2 mm for whole brain treatments, respectively. On the other hand, current margins seem to be insufficient for some breast, lung and head and neck patients, at least if used generically. If no case specific adjustments are applied, a generic margin of 6.3% + 1.2 mm would be needed for breast, lung and head and neck treatments. We conclude that the currently used generic range uncertainty margins in proton therapy should be redefined site specific and that complex geometries may require a field specific adjustment. Routine verifications of treatment plans using MC simulations are recommended for patients with heterogeneous geometries.