A Millimeter-scale Single Charged Particle Dosimeter for Cancer Radiotherapy.

This paper presents a millimeter-scale CMOS 64×64 single charged particle radiation detector system for external beam cancer radiotherapy. A 1×1 µm2 diode measures energy deposition by a single charged particle in the depletion region, and the array design provides a large detection area of 512×512 µm2. Instead of sensing the voltage drop caused by radiation, the proposed system measures the pulse width, i.e., the time it takes for the voltage to return to its baseline. This obviates the need for using power-hungry and large analog-to-digital converters. A prototype ASIC is fabricated in TSMC 65 nm LP CMOS process and consumes the average static power of 0.535 mW under 1.2 V analog and digital power supply. The functionality of the whole system is successfully verified in a clinical 67.5 MeV proton beam setting. To our' knowledge, this is the first work to demonstrate single charged particle detection for implantable in-vivo dosimetry.

Clinical trials involving carbon-ion radiation therapy and the path forward.

To describe the international landscape of clinical trials in carbon-ion radiotherapy (CIRT), the authors reviewed the current status of 63 ongoing clinical trials (median, 47 participants) involving CIRT identified from the US clinicaltrials.gov trial registry and the World Health Organization International Clinical Trials Platform Registry. The objectives were to evaluate the potential for these trials to define the role of this modality in the treatment of specific cancer types and identify the major challenges and opportunities to advance this technology. A significant body of literature suggested the potential for advantageous dose distributions and, in preclinical biologic studies, the enhanced effectiveness for CIRT compared with photons and protons. In addition, clinical evidence from phase I/II trials, although limited, indicated the potential for CIRT to improve cancer outcomes. However, current high-level phase III randomized clinical trial evidence does not exist. Although there has been an increase in the number of trials investigating CIRT since 2010, and the number of countries and sites offering CIRT is slowly growing, this progress has excluded other countries. Several recommendations are proposed to study this modality to accelerate progress in the field, including: 1) increasing the number of multinational randomized clinical trials, 2) leveraging the existing CIRT facilities to launch larger multinational trials directed at common cancers combined with high-level quality assurance; and 3) developing more compact and less expensive next-generation treatment systems integrated with radiobiologic research and preclinical testing.

Creating uniform cluster dose spread-out Bragg peaks for proton and carbon beams.

BACKGROUND: Preliminary data have shown a close association of the generalized ionization cluster size dose (in short, cluster dose) with cell survival, independent of particle type, and energy, when cluster dose is derived from an ionization detail parameter preferred for its association with cell survival. Such results suggest cluster dose has the potential to replace RBE-weighted dose in proton and ion beam radiotherapy treatment plan optimization, should a uniform cluster dose lead to comparable biological effects. However, further preclinical investigations are warranted to confirm this premise. PURPOSE: To present an analytical approach to create uniform cluster dose spread-out Bragg peaks (SOBP) for evaluation of the potential of cluster dose to result in uniform biological effect. METHODS: We modified the coefficients of the Bortfeld and Schlegel weight formula, an analytical method typically used for the creation of radiation dose SOBP in particle therapy, to produce uniform cluster dose SOBP of different widths (1-5 cm) at relevant clinical proton and carbon beam energies. Optimum parameters were found by minimization of the ratio between the maximum and minimum cluster dose in the SOBP region using the Nelder-Mead method. RESULTS: The coefficients of the Bortfeld and Schlegel weight formula leading to uniform cluster dose SOBPs were determined for each combination of beam energy and SOBP width studied. The uniformity of the resulting cluster dose SOBPs, calculated as the relative difference between the maximum and minimum cluster dose within the SOBP, was within 0.3%-3.5% for the evaluated proton beams and 1.3%-3.4% for the evaluated carbon beams. CONCLUSIONS: The modifications to the analytical approach to create radiation dose SOBPs resulted in uniform cluster dose proton and carbon SOBPs over a wide range of beam energies and SOBP widths. Such SOBPs should prove valuable in preclinical investigations for the selection of nanodosimetric quantities to be used in proton and ion therapy treatment planning.

A digital male pelvis phantom series showing anatomical variations over the course of fractionated radiotherapy treatment.

BACKGROUND: Daily IGRT images show day-to-day anatomical variations in patients undergoing fractionated prostate radiotherapy. This is of particular importance in particle beam treatments. PURPOSE: To develop a digital phantom series showing variation in pelvic anatomy for evaluating treatment planning and IGRT procedures in particle radiotherapy. METHODS: A pelvic phantom series was developed from the planning MRI and kVCT (planning CT) images along with six of the daily serial MVCT images taken of a single patient treated with a full bladder on a Tomotherapy unit. The selected patient had clearly visible yet unexceptional internal anatomy variation. Prostate, urethra, bladder, rectum, bowel, bowel gas, bone and soft tissue were contoured and a single Hounsfield Unit was assigned to each region. Treatment plans developed on the kVCT for photon, proton and carbon beams were recalculated on each phantom to demonstrate a clinical application of the series. Proton plans were developed with and without robust optimization. RESULTS: Limited to axial slices with prostate, the bladder volume varied from 6 to 46 cm3, the rectal volume (excluding gas) from 22 to 52 cm3, and rectal gas volume from zero to 18 cm3. The water equivalent path length to the prostate varied by up to 1.5 cm . The variations resulted in larger changes in the RBE-weighted Dose Volume Histograms of the non-robust proton plan and the carbon plan compared to the robust proton plan, the latter similar to the photon plan. The prostate coverage (V100%) decreased by an average of 18% in the carbon plan, 16% in the non-robust proton plan, 1.8% in the robust proton plan, and 4.4% in the photon plan. The volume of rectum receiving 75% of the prescription dose (V75%) increased by an average of 3.7 cm3, 4.7 cm3, 1.9 cm3, and 0.6 cm3 in those four plans, respectively. CONCLUSIONS: The digital pelvic phantom series provides for quantitative investigation of IGRT procedures and new methods for improving accuracy in particle therapy and may be used in cross-institutional comparisons for clinical trial quality assurance.

Lithium inelastic cross-sections and their impact on micro and nano dosimetry of boron neutron capture.

OBJECTIVE: To present a new set of lithium-ion cross-sections for (i) ionization and excitation processes down to 700 eV, and (ii) charge-exchange processes down to 1 keV/u. To evaluate the impact of the use of these cross-sections on micro a nano dosimetric quantities in the context of boron neutron capture (BNC) applications/techniques. Approach: The Classical Trajectory Monte Carlo (CTMC) method was used to calculate Li ion charge-exchange cross sections in the energy range of 1 keV/u to 10 MeV/u. Partial Li ion charge states ionization and excitation cross-sections were calculated using a detailed charge screening factor. The cross-sections were implemented in Geant4-DNA v10.07 and simulations and verified using TOPAS-nBio by calculating stopping power and CSDA range against data from ICRU and SRIM. Further microdosimetric and nanodosimetric calculations were performed to quantify differences against other simulation approaches for low energy Li ions. These calculations were: lineal energy spectra (yf(y) and yd(y)), frequency mean lineal energy (y_F ) ̅, dose mean lineal energy (y_D ) ̅ and ionization cluster size distribution analysis. Microdosimetric calculations were compared against a previous MC study that neglected charge-exchange and excitation processes. Nanodosimetric results were compared against pure ionization scaled cross-sections calculations. Main Results: Calculated stopping power differences between ICRU and Geant4-DNA decreased from 33.78% to 6.9%. The CSDA range difference decreased from 621% to 34% when compared against SRIM calculations. Geant4-DNA/TOPAS calculated dose mean lineal energy differed by 128% from the previous Monte Carlo. Ionization cluster size frequency distributions for Li ions differed by 76% to 344.11% for 21 keV and 2 MeV respectively. With a decrease in the N1 within 9% at 10 keV and agreeing after the 100 keV. With the new set of cross-sections being able to better simulate low energy behaviors of Li ions. Significance: This work shows an increase in detail gained from the use of a more complete set of low energy cross-sections which include charge exchange processes. Significant differences to previous simulation results were found at the microdosimetric and nanodosimetric scales that suggest that Li ions cause less ionizations per path length traveled but with more energy deposits. Microdosimetry results suggest that the BNC's contribution to cellular death may be mainly due to alpha particle production when boron-based drugs are distributed in the cellular membrane and beyond and by Li when it is at the cell cytoplasm regions.

An integrated Monte Carlo track-structure simulation framework for modeling inter and intra-track effects on homogenous chemistry.

Objective. The TOPAS-nBio Monte Carlo track structure simulation code, a wrapper of Geant4-DNA, was extended for its use in pulsed and longtime homogeneous chemistry simulations using the Gillespie algorithm.Approach. Three different tests were used to assess the reliability of the implementation and its ability to accurately reproduce published experimental results: (1) a simple model with a known analytical solution, (2) the temporal evolution of chemical yields during the homogeneous chemistry stage, and (3) radiolysis simulations conducted in pure water with dissolved oxygen at concentrations ranging from 10µM to 1 mM with [H2O2] yields calculated for 100 MeV protons at conventional and FLASH dose rates of 0.286 Gy s-1and 500 Gy s-1, respectively. Simulated chemical yield results were compared closely with data calculated using the Kinetiscope software which also employs the Gillespie algorithm.Main results. Validation results in the third test agreed with experimental data of similar dose rates and oxygen concentrations within one standard deviation, with a maximum of 1% difference for both conventional and FLASH dose rates. In conclusion, the new implementation of TOPAS-nBio for the homogeneous long time chemistry simulation was capable of recreating the chemical evolution of the reactive intermediates that follow water radiolysis.Significance. Thus, TOPAS-nBio provides a reliable all-in-one chemistry simulation of the physical, physico-chemical, non-homogeneous, and homogeneous chemistry and could be of use for the study of FLASH dose rate effects on radiation chemistry.

Ionization detail parameters and cluster dose: a mathematical model for selection of nanodosimetric quantities for use in treatment planning in charged particle radiotherapy.

Objective. To propose a mathematical model for applying ionization detail (ID), the detailed spatial distribution of ionization along a particle track, to proton and ion beam radiotherapy treatment planning (RTP).Approach. Our model provides for selection of preferred ID parameters (Ip) for RTP, that associate closest to biological effects. Cluster dose is proposed to bridge the large gap between nanoscopicIpand macroscopic RTP. Selection ofIpis demonstrated using published cell survival measurements for protons through argon, comparing results for nineteenIp:Nk,k= 2, 3, …, 10, the number of ionizations in clusters ofkor more per particle, andFk,k= 1, 2, …, 10, the number of clusters ofkor more per particle. We then describe application of the model to ID-based RTP and propose a path to clinical translation.Main results. The preferredIpwereN4andF5for aerobic cells,N5andF7for hypoxic cells. Significant differences were found in cell survival for beams having the same LET or the preferredNk. Conversely, there was no significant difference forF5for aerobic cells andF7for hypoxic cells, regardless of ion beam atomic number or energy. Further, cells irradiated with the same cluster dose for theseIphad the same cell survival. Based on these preliminary results and other compelling results in nanodosimetry, it is reasonable to assert thatIpexist that are more closely associated with biological effects than current LET-based approaches and microdosimetric RBE-based models used in particle RTP. However, more biological variables such as cell line and cycle phase, as well as ion beam pulse structure and rate still need investigation.Significance. Our model provides a practical means to select preferredIpfrom radiobiological data, and to convertIpto the macroscopic cluster dose for particle RTP.

TOPAS-imaging: extensions to the TOPAS simulation toolkit for medical imaging systems.

Objective. The TOol for PArticle Simulation (TOPAS) is a Geant4-based Monte Carlo software application that has been used for both research and clinical studies in medical physics. So far, most users of TOPAS have focused on radiotherapy-related studies, such as modeling radiation therapy delivery systems or patient dose calculation. Here, we present the first set of TOPAS extensions to make it easier for TOPAS users to model medical imaging systems.Approach. We used the extension system of TOPAS to implement pre-built, user-configurable geometry components such as detectors (e.g. flat-panel and multi-planar detectors) for various imaging modalities and pre-built, user-configurable scorers for medical imaging systems (e.g. digitizer chain).Main results. We developed a flexible set of extensions that can be adapted to solve research questions for a variety of imaging modalities. We then utilized these extensions to model specific examples of cone-beam CT (CBCT), positron emission tomography (PET), and prompt gamma (PG) systems. The first of these new geometry components, the FlatImager, was used to model example CBCT and PG systems. Detected signals were accumulated in each detector pixel to obtain the intensity of x-rays penetrating objects or prompt gammas from proton-nuclear interaction. The second of these new geometry components, the RingImager, was used to model an example PET system. Positron-electron annihilation signals were recorded in crystals of the RingImager and coincidences were detected. The simulated data were processed using corresponding post-processing algorithms for each modality and obtained results in good agreement with the expected true signals or experimental measurement.Significance. The newly developed extension is a first step to making it easier for TOPAS users to build and simulate medical imaging systems. Together with existing TOPAS tools, this extension can help integrate medical imaging systems with radiotherapy simulations for image-guided radiotherapy.

Sensitivity analysis of an asymmetric Monte Carlo beam model of a Siemens Primus accelerator.

The assumption of cylindrical symmetry in radiotherapy accelerator models can pose a challenge for precise Monte Carlo modeling. This assumption makes it difficult to account for measured asymmetries in clinical dose distributions. We have performed a sensitivity study examining the effect of varying symmetric and asymmetric beam and geometric parameters of a Monte Carlo model for a Siemens PRIMUS accelerator. The accelerator and dose output were simulated using modified versions of BEAMnrc and DOSXYZnrc that allow lateral offsets of accelerator components and lateral and angular offsets for the incident electron beam. Dose distributions were studied for 40 × 40 cm² fields. The resulting dose distributions were analyzed for changes in flatness, symmetry, and off-axis ratio (OAR). The electron beam parameters having the greatest effect on the resulting dose distributions were found to be electron energy and angle of incidence, as high as 5% for a 0.25° deflection. Electron spot size and lateral offset of the electron beam were found to have a smaller impact. Variations in photon target thickness were found to have a small effect. Small lateral offsets of the flattening filter caused significant variation to the OAR. In general, the greatest sensitivity to accelerator parameters could be observed for higher energies and off-axis ratios closer to the central axis. Lateral and angular offsets of beam and accelerator components have strong effects on dose distributions, and should be included in any high-accuracy beam model.

TOPAS-nBio simulation of temperature-dependent indirect DNA strand break yields.

Current Monte Carlo simulations of DNA damage have been reported only at ambient temperature. The aim of this work is to use TOPAS-nBio to simulate the yields of DNA single-strand breaks (SSBs) and double-strand breaks (DSBs) produced in plasmids under low-LET irradiation incorporating the effect of the temperature changes in the environment. A new feature was implemented in TOPAS-nBio to incorporate reaction rates used in the simulation of the chemical stage of water radiolysis as a function of temperature. The implemented feature was verified by simulating temperature-dependentG-values of chemical species in liquid water from 20 °C to 90 °C. For radiobiology applications, temperature dependent SSB and DSB yields were calculated from 0 °C to 42 °C, the range of available published measured data. For that, supercoiled DNA plasmids dissolved in aerated solutions containing EDTA irradiated by Cobalt-60 gamma-rays were simulated. TOPAS-nBio well reproduced published temperature-dependentG-values in liquid water and the yields of SSB and DSB for the temperature range considered. For strand break simulations, the model shows that the yield of SSB and DSB increased linearly with the temperature at a rate of (2.94 ± 0.17) × 10-10Gy-1Da-1°C-1(R2 = 0.99) and (0.13 ± 0.01) × 10-10Gy-1Da-1°C-1(R2 = 0.99), respectively. The extended capability of TOPAS-nBio is a complementary tool to simulate realistic conditions for a large range of environmental temperatures, allowing refined investigations of the biological effects of radiation.

Impact of DNA Geometry and Scoring on Monte Carlo Track-Structure Simulations of Initial Radiation-Induced Damage.

Track structure Monte Carlo simulations are a useful tool to investigate the damage induced to DNA by ionizing radiation. These simulations usually rely on simplified geometrical representations of the DNA subcomponents. DNA damage is determined by the physical and physicochemical processes occurring within these volumes. In particular, damage to the DNA backbone is generally assumed to result in strand breaks. DNA damage can be categorized as direct (ionization of an atom part of the DNA molecule) or indirect (damage from reactive chemical species following water radiolysis). We also consider quasi-direct effects, i.e., damage originated by charge transfers after ionization of the hydration shell surrounding the DNA. DNA geometries are needed to account for the damage induced by ionizing radiation, and different geometry models can be used for speed or accuracy reasons. In this work, we use the Monte Carlo track structure tool TOPAS-nBio, built on top of Geant4-DNA, for simulation at the nanometer scale to evaluate differences among three DNA geometrical models in an entire cell nucleus, including a sphere/spheroid model specifically designed for this work. In addition to strand breaks, we explicitly consider the direct, quasi-direct, and indirect damage induced to DNA base moieties. We use results from the literature to determine the best values for the relevant parameters. For example, the proportion of hydroxyl radical reactions between base moieties was 80%, and between backbone, moieties was 20%, the proportion of radical attacks leading to a strand break was 11%, and the expected ratio of base damages and strand breaks was 2.5-3. Our results show that failure to update parameters for new geometric models can lead to significant differences in predicted damage yields.

Accounting for the fringe magnetic field from the bending magnet in a Monte Carlo accelerator treatment head simulation.

PURPOSE: Monte Carlo (MC) simulation can be used for accurate electron beam treatment planning and modeling. Measurement of large electron fields, with the applicator removed and secondary collimator wide open, has been shown to provide accurate simulation parameters, including asymmetry in the measured dose, for the full range of clinical field sizes and patient positions. Recently, disassembly of the treatment head of a linear accelerator has been used to refine the simulation of the electron beam, setting tightly measured constraints on source and geometry parameters used in simulation. The simulation did not explicitly include the known deflection of the electron beam by a fringe magnetic field from the bending magnet, which extended into the treatment head. Instead, the secondary scattering foil and monitor chamber were unrealistically laterally offset to account for the beam deflection. This work is focused on accounting for this fringe magnetic field in treatment head simulation. METHODS: The magnetic field below the exit window of a Siemens Oncor linear accelerator was measured with a Tesla-meter from 0 to 12 cm from the exit window and 1-3 cm off-axis. Treatment head simulation was performed with the EGSnrc/BEAMnrc code, modified to incorporate the effect of the magnetic field on charged particle transport. Simulations were used to analyze the sensitivity of dose profiles to various sources of asymmetry in the treatment head. This included the lateral spot offset and beam angle at the exit window, the fringe magnetic field and independent lateral offsets of the secondary scattering foil and electron monitor chamber. Simulation parameters were selected within the limits imposed by measurement uncertainties. Calculated dose distributions were then compared with those measured in water. RESULTS: The magnetic field was a maximum at the exit window, increasing from 0.006 T at 6 MeV to 0.020 T at 21 MeV and dropping to approximately 5% of the maximum at the secondary scattering foil. It was up to three times higher in the bending plane, away from the electron gun, and symmetric within measurement uncertainty in the transverse plane. Simulations showed the magnetic field resulted in an offset of the electron beam of 0.80 cm (mean) at the machine isocenter for the exit window only configuration. The fringe field resulted in a 3.5%-7.6% symmetry and 0.25-0.35 cm offset of the clinical beam R(max) profiles. With the magnetic field included in simulations, a single (realistic) position of the secondary scattering foil and monitor chamber was selected. Measured and simulated dose profiles showed agreement to an average of 2.5%/0.16 cm (maximum: 3%/0.2 cm), which is a better match than previously achieved without incorporating the magnetic field in the simulation. The undulations from the 3 stepped layers of the secondary scattering foil, evident in the measured profiles of the higher energy beams, are now aligned with those in the simulated beam. The simulated fringe magnetic field had negligible effect on the central axis depth dose curves and cross-plane dose profiles. CONCLUSIONS: The fringe magnetic field is a significant contributor to the electron beam in-plane asymmetry. With the magnetic field included explicitly in the simulation, realistic monitor chamber and secondary scattering foil positions have been achieved, and the calculated fluence and dose distributions are more accurate.

Dispersion calibration for the National Ignition Facility electron-positron-proton spectrometers for intense laser matter interactions.

Electron-positron pairs, produced in intense laser-solid interactions, are diagnosed using magnetic spectrometers with image plates, such as the National Ignition Facility Electron-Positron-Proton Spectrometers (EPPSs). Although modeling can help infer the quantitative value, the accuracy of the models needs to be verified to ensure measurement quality. The dispersion of low-energy electrons and positrons may be affected by fringe magnetic fields near the entrance of the EPPS. We have calibrated the EPPS with six electron beams from a Siemens Oncor linear accelerator (linac) ranging in energy from 2.7 MeV to 15.2 MeV as they enter the spectrometer. A Geant4 traveling-wave optical parametric amplifier of superfluorescence Monte Carlo simulation was set up to match depth dose curves and lateral profiles measured in water at 100 cm source-surface distance. An accurate relationship was established between the bending magnet current setting and the energy of the electron beam at the exit window. The simulations and measurements were used to determine the energy distributions of the six electron beams at the EPPS slit. Analysis of the scanned image plates together with the determined energy distribution arriving in the spectrometer provides improved dispersion curves for the EPPS.

Assessment of DNA damage with an adapted independent reaction time approach implemented in Geant4-DNA for the simulation of diffusion-controlled reactions between radio-induced reactive species and a chromatin fiber.

PURPOSE: Simulation of indirect damage originating from the attack of free radical species produced by ionizing radiation on biological molecules based on the independent pair approximation is investigated in this work. In addition, a new approach, relying on the independent pair approximation that is at the origin of the independent reaction time (IRT) method, is proposed in the chemical stage of Geant4-DNA. METHODS: This new approach has been designed to respect the current Geant4-DNA chemistry framework while proposing a variant IRT method. Based on the synchronous algorithm, this implementation allows us to access the information concerning the position of radicals and may make it more convenient for biological damage simulations. Estimates of the evolution of free species as well as biological hits in a segment of DNA chromatin fiber in Geant4-DNA were compared for the dynamic time step approach of the step-by-step (SBS) method, currently used in Geant4-DNA, and this newly implemented IRT. RESULTS: Results show a gain in computation time of a factor of 30 for high LET particle tracks with a better than 10% agreement on the number of DNA hits between the value obtained with the IRT method as implemented in this work and the SBS method currently available in Geant4-DNA. CONCLUSION: Offering in Geant4-DNA more efficient methods for the chemical step based on the IRT method is a task in progress. For the calculation of biological damage, information on the position of chemical species is a crucial point. This can be achieved using the method presented in this paper.

Comparison of patient megavoltage cone beam CT images acquired with an unflattened beam from a carbon target and a flattened treatment beam.

PURPOSE: To use an imaging beam line (IBL) to obtain the first megavoltage cone-beam computed tomography (MV CBCT) images of patients with a low atomic number (Z) target, and to compare these images to those taken of the same patients with the 6 MV flattened beam from the treatment beam line (TBL). METHODS: The IBL, which produces a 4.2 MV unflattened beam from a carbon target, was installed on a linear accelerator in use for radiotherapy. Provision was made for switching between the IBL and TBL for imaging the same patient with beams from the low-Z and high-Z targets. Dose was quoted as monitor units times the dose per monitor unit for the standard calibration geometry. Images were acquired with institutional approval and patient consent with both the IBL and TBL on a series of 23 patients undergoing radiotherapy. Patients were imaged daily to weekly and aligned to the planning CT using the images. Doses were reduced over the course of treatment to determine the minimum doses required for alignment. Images were assessed offline. RESULTS: IBL MV CBCT images of prostate, head and neck, lung, and abdomen showed improvement in soft tissue contrast for the same dose as the TBL images. Bony anatomy, air cavities, and fiducial markers were sharper. CBCT with a dose of 1 cGy was sufficient for alignment of prostate and head and neck patients based on bony anatomy or implanted gold seeds, 2-4 cGy for lung, abdomen, and pelvis. Photon scatter in the patient had minimal effect on image quality. The metallic hip prosthesis in one patient showed reduced artifacts compared to diagnostic CT. CONCLUSIONS: The IBL has the advantage of improved image quality at the same dose, or reduced dose for the same image quality, over the TBL.

Practical quantitative measurement of graticule misalignment relative to collimator axis of rotation.

We design a practical procedure for measuring translational and rotational misalignment of graticule with collimator axis of rotation and collimator jaws, respectively. The procedure's quantitative results are accurate to less than 0.2 mm (at isocenter) and do not assume alignment of radiation focal spot with collimator axis of rotation. When provided with these quantitative results, the manufacturer can custom-adjust graticules to the purchaser's collimator head.

LET-Dependent Intertrack Yields in Proton Irradiation at Ultra-High Dose Rates Relevant for FLASH Therapy.

FLASH radiotherapy delivers a high dose (≥10 Gy) at a high rate (≥40 Gy/s). In this way, particles are delivered in pulses as short as a few nanoseconds. At that rate, intertrack reactions between chemical species produced within the same pulse may affect the heterogeneous chemistry stage of water radiolysis. This stochastic process suits the capabilities of the Monte Carlo method, which can model intertrack effects to aid in radiobiology research, including the design and interpretation of experiments. In this work, the TOPAS-nBio Monte Carlo track-structure code was expanded to allow simulations of intertrack effects in the chemical stage of water radiolysis. Simulation of the behavior of radiolytic yields over a long period of time (up to 50 s) was verified by simulating radiolysis in a Fricke dosimeter irradiated by 60Co γ rays. In addition, LET-dependent G values of protons delivered in single squared pulses of widths, 1 ns, 1 µs and 10 µs, were obtained and compared to simulations using no intertrack considerations. The Fricke simulation for the calculated G value of Fe3+ ion at 50 s was within 0.4% of the accepted value from ICRU Report 34. For LET-dependent G values at the end of the chemical stage, intertrack effects were significant at LET values below 2 keV/µm. Above 2 keV/µm the reaction kinetics remained limited locally within each track and thus, effects of intertrack reactions remained low. Therefore, when track structure simulations are used to investigate the biological damage of FLASH irradiation, these intertrack reactions should be considered. The TOPAS-nBio framework with the expansion to intertrack chemistry simulation provides a useful tool to assist in this task.

The TOPAS tool for particle simulation, a Monte Carlo simulation tool for physics, biology and clinical research.

PURPOSE: This paper covers recent developments and applications of the TOPAS TOol for PArticle Simulation and presents the approaches used to disseminate TOPAS. MATERIALS AND METHODS: Fundamental understanding of radiotherapy and imaging is greatly facilitated through accurate and detailed simulation of the passage of ionizing radiation through apparatus and into a patient using Monte Carlo (MC). TOPAS brings Geant4, a reliable, experimentally validated MC tool mainly developed for high energy physics, within easy reach of medical physicists, radiobiologists and clinicians. Requiring no programming knowledge, TOPAS provides all of the flexibility of Geant4. RESULTS: After 5 years of development followed by its initial release, TOPAS was subsequently expanded from its focus on proton therapy physics to incorporate radiobiology modeling. Next, in 2018, the developers expanded their user support and code maintenance as well as the scope of TOPAS towards supporting X-ray and electron therapy and medical imaging. Improvements have been achieved in user enhancement through software engineering and a graphical user interface, calculational efficiency, validation through experimental benchmarks and QA measurements, and either newly available or recently published applications. A large and rapidly increasing user base demonstrates success in our approach to dissemination of this uniquely accessible and flexible MC research tool. CONCLUSIONS: The TOPAS developers continue to make strides in addressing the needs of the medical community in applications of ionizing radiation to medicine, creating the only fully integrated platform for four-dimensional simulation of all forms of radiotherapy and imaging with ionizing radiation, with a design that promotes inter-institutional collaboration.

Independent reaction times method in Geant4-DNA: Implementation and performance.

PURPOSE: The simulation of individual particle tracks and the chemical stage following water radiolysis in biological tissue is an effective means of improving our knowledge of the physico-chemical contribution to the biological effect of ionizing radiation. However, the step-by-step simulation of the reaction kinetics of radiolytic species is the most time-consuming task in Monte Carlo track-structure simulations, with long simulation times that are an impediment to research. In this work, we present the implementation of the independent reaction times (IRT) method in Geant4-DNA Monte Carlo toolkit to improve the computational efficiency of calculating G-values, defined as the number of chemical species created or lost per 100 eV of deposited energy. METHODS: The computational efficiency of IRT, as implemented, is compared to that from available Geant4-DNA step-by-step simulations for electrons, protons and alpha particles covering a wide range of linear energy transfer (LET). The accuracy of both methods is verified using published measured data from fast electron irradiations for • OH and e aq - for time-dependent G-values. For IRT, simulations in the presence of scavengers irradiated by cobalt-60 γ-ray and 2 MeV protons are compared with measured data for different scavenging capacities. In addition, a qualitative assessment comparing measured LET-dependent G-values with Geant4-DNA calculations in pure liquid water is presented. RESULTS: The IRT improved the computational efficiency by three orders of magnitude relative to the step-by-step method while differences in G-values by 3.9% at 1 µs were found. At 7 ps, • OH and e aq - yields calculated with IRT differed from recent published measured data by 5% ± 4% and 2% ± 4%, respectively. At 1 µs, differences were 9% ± 5% and 6% ± 7% for • OH and e aq - , respectively. Uncertainties are one standard deviation. Finally, G-values at different scavenging capacities and LET-dependent G-values reproduced the behavior of measurements for all radiation qualities. CONCLUSION: The comprehensive validation of the Geant4-DNA capabilities to accurately simulate the chemistry following water radiolysis is an ongoing work. The implementation presented in this work is a necessary step to facilitate performing such a task.