In-silico calculations of DNA damage induced by α-particles in the 224Ra DaRT decay chain for a better understanding of the radiobiological effectiveness of this treatment.

Diffusing alpha-emitters radiation Therapy (DaRT) is an interstitial brachytherapy technique using 224Ra seeds. For accurate treatment planning a good understanding of the early DNA damage due to α-particles is required. Geant4-DNA was used to calculate the initial DNA damage and radiobiological effectiveness due to α-particles with linear energy transfer (LET) values in the range 57.5-225.9 keV/µm from the 224Ra decay chain. The impact of DNA base pair density on DNA damage has been modelled, as this parameter varies between human cell lines. Results show that the quantity and complexity of DNA damage changes with LET as expected. Indirect damage, due to water radical reactions with the DNA, decreases and becomes less significant at higher LET values as shown in previous studies. As expected, the yield of complex double strand breaks (DSBs), which are harder for a cell to repair, increases approximately linearly with LET. The level of complexity of DSBs and radiobiological effectiveness have been found to increase with LET as expected. The quantity of DNA damage has been shown to increase for increased DNA density in the expected base pair density range of human cells. The change in damage yield as a function of base pair density is largest for higher LET α-particles, an increase of over 50% for individual strand breaks between 62.7 and 127.4 keV/µm. This change in yield shows that the DNA base pair density is an important parameter for modelling DNA damage particularly at higher LET where the DNA damage is greatest and most complex.

Geant4-DNA simulation of human cancer cells irradiation with helium ion beams.

PURPOSE: This study aimed to develop a computational environment for the accurate simulation of human cancer cell irradiation using Geant4-DNA. New cell geometrical models were developed and irradiated by alpha particle beams to induce DNA damage. The proposed approach may help further investigation of the benefits of external alpha irradiation therapy. METHODS: The Geant4-DNA Monte Carlo (MC) toolkit allows the simulation of cancer cell geometries that can be combined with accurate modelling of physical, physicochemical and chemical stages of liquid water irradiation, including radiolytic processes. Geant4-DNA is used to calculate direct and non-direct DNA damage yields, such as single and double strand breaks, produced by the deposition of energy or by the interaction of DNA with free radicals. RESULTS: In this study, the "molecularDNA" example application of Geant4-DNA was used to quantify early DNA damage in human cancer cells upon irradiation with alpha particle beams, as a function of linear energy transfer (LET). The MC simulation results are compared to experimental data, as well as previously published simulation data. The simulation results agree well with the experimental data on DSB yields in the lower LET range, while the experimental data on DSB yields are lower than the results obtained with the "molecularDNA" example in the higher LET range. CONCLUSION: This study explored and demonstrated the possibilities of the Geant4-DNA toolkit together with the "molecularDNA" example to simulate the helium beam irradiation of cancer cell lines, to quantify the early DNA damage, or even the following DNA damage response.

The general-purpose Geant4 Monte Carlo toolkit and its Geant4-DNA extension to investigate mechanisms underlying the FLASH effect in radiotherapy: Current status and challenges.

FLASH radiotherapy is a promising approach to cancer treatment that offers several advantages over conventional radiotherapy. With this novel technique, high doses of radiation are delivered in a short period of time, inducing the so-called FLASH effect - a phenomenon characterized by healthy tissue sparing without alteration of tumor control. The mechanisms behind the FLASH effect remain unknown. One way to approach this problem is to gain insight into the initial parameters that can distinguish FLASH from conventional irradiation by simulating particle transport in aqueous media using the general-purpose Geant4 Monte Carlo toolkit and its Geant4-DNA extension. This review article discusses the current status of Geant4 and Geant4-DNA simulations to investigate mechanisms underlying the FLASH effect, as well as the challenges faced in this research field. One of the primary challenges is to accurately simulate the experimental irradiation parameters. Another challenge is the temporal extension of the simulations. This review also focuses on two hypotheses to explain the FLASH effect - namely the oxygen depletion hypothesis and the inter-track interactions hypothesis - and discusses how the Geant4 toolkit can be used to investigate them. The aim of this review is to provide an overview of Geant4 and Geant4-DNA simulations for FLASH radiotherapy and to highlight the challenges that need to be overcome in order to better study the FLASH effect.

Modeling of scavenging systems in water radiolysis with Geant4-DNA.

PURPOSE: This paper presents the capabilities of the Geant4-DNA Monte Carlo toolkit to simulate water radiolysis with scavengers using the step-by-step (SBS) or the independent reaction times (IRT) methods. It features two examples of application areas: (1) computing the escape yield of H2O2 following a 60Co γ-irradiation and (2) computing the oxygen depletion in water irradiated with 1 MeV electrons. METHODS: To ease the implementation of the chemical stage in Geant4-DNA, we developed a user interface that helps define the chemical reactions and set the concentration of scavengers. The first application area example required two computational steps to perform water radiolysis using NO2- and NO3- as scavengers and a 60Co irradiation. The oxygen depletion computation technique for the second application area example consisted of simulating track segments of 1 MeV electrons and determining the radio-induced loss and gain of oxygen molecules. RESULTS: The production of H2O2 under variable scavenging levels is consistent with the literature; the mean relative difference between the SBS and IRT methods is 7.2 % ± 0.5 %. For the oxygen depletion 1 µs post-irradiation, the mean relative difference between both methods is equal to 9.8 % ± 0.3 %. The results in the microsecond scale depend on the initial partial pressure of oxygen in water. In addition, the computed oxygen depletions agree well with the literature. CONCLUSIONS: The Geant4-DNA toolkit makes it possible to simulate water radiolysis in the presence of scavengers. This feature offers perspectives in radiobiology, with the possibility of simulating cell-relevant scavenging mechanisms.

Sub-keV corrections to binary encounter cross section models for electron ionization of liquid water with application to the Geant4-DNA Monte Carlo code.

INTRODUCTION: The electron ionization cross section of water is one of the most important input in Monte Carlo studies of cellular radiobiological effects. Analytical cross section models of the binary-encounter type have the potential of reducing simulation time and facilitate application to a variety of biological materials (other than water). The Binary-Encounter-Bethe (BEB) and Binary-Encounter-Dipole (BED) models of NIST are perhaps the most popular of such models giving reliable results for atoms and molecules in the gas-phase over a wide energy range. However, the use of such models to sub-keV electron energies in liquid water raises concerns due to the neglect of condensed phase effects that leads to a significant overestimation when compared to medium-specific dielectric models. PURPOSE: To modify the BEB and BED models towards better agreement with the recommended low-energy dielectric model of Geant4-DNA (Option 4). To implement the new modifications to the existing BEB model of the Option 6 physics constructor of Geant4-DNA and re-evaluate fundamental transport quantities for sub-keV electrons. METHODS: In analogy to a Yukawa potential a simple, yet physically-motivated, modification of the Burgess correction term is proposed to account for the reduction of the Coulomb interaction due to the polarizability of the target. The magnitude of the correction is guided by the dielectric-based ionization cross section implemented in Option 4. RESULTS: Differential, total and stopping ionization cross sections for low-energy electrons in liquid water are presented. When combined with the Vriens correction (which is not included in Option 6), the proposed modification to the BEB and BED models brings the ionization and stopping cross sections in much better agreement against those used in the Option 4 dielectric model of Geant4-DNA, with up to 30% and 10% deviation, respectively. Implementation of the new correction to the Option 6 constructor of Geant4-DNA and re-evaluation of fundamental transport quantities, such as electron penetration ranges and dose-point-kernels, reduced the discrepancies from Option 4 at sub-keV energies from 20 to 100% (or more) to well below 10% in most cases. CONCLUSIONS: A simple modification to the BEB and BED analytic models was found to improve their performance for sub-keV electrons in liquid water medium. Implementation of the new modification to the Option 6 constructor of Geant4-DNA significantly improved the agreement with the recommended low-energy Option 4 constructor for a variety of fundamental quantities related to electron transport.

Prediction of DNA rejoining kinetics and cell survival after proton irradiation for V79 cells using Geant4-DNA.

PURPOSE: Track structure Monte Carlo (MC) codes have achieved successful outcomes in the quantitative investigation of radiation-induced initial DNA damage. The aim of the present study is to extend a Geant4-DNA radiobiological application by incorporating a feature allowing for the prediction of DNA rejoining kinetics and corresponding cell surviving fraction along time after irradiation, for a Chinese hamster V79 cell line, which is one of the most popular and widely investigated cell lines in radiobiology. METHODS: We implemented the Two-Lesion Kinetics (TLK) model, originally proposed by Stewart, which allows for simulations to calculate residual DNA damage and surviving fraction along time via the number of initial DNA damage and its complexity as inputs. RESULTS: By optimizing the model parameters of the TLK model in accordance to the experimental data on V79, we were able to predict both DNA rejoining kinetics at low linear energy transfers (LET) and cell surviving fraction. CONCLUSION: This is the first study to demonstrate the implementation of both the cell surviving fraction and the DNA rejoining kinetics with the estimated initial DNA damage, in a realistic cell geometrical model simulated by full track structure MC simulations at DNA level and for various LET. These simulation and model make the link between mechanistic physical/chemical damage processes and these two specific biological endpoints.

Intercomparison of nanodosimetric distributions in nitrogen simulated with Geant4 and PTra track structure codes.

To facilitate the use of Geant4-DNA for radiation transport simulations in micro- and nanodosimeters, which are physically operated with tissue-equivalent gases such as nitrogen (and propane), this work aims to extend the cross section data available in Geant4-DNA to include those of nitrogen for electron energies ranging from 1 MeV down to the ionisation threshold. To achieve this, interaction cross section data for nitrogen that have been used with the in-house PTB PTra track structure code have been implemented in the current state-of-the-art Geant4-DNA simulation toolkit. An intercomparison has been performed between the two codes to validate this implementation. To quantify the agreement between the cross section models for nitrogen adopted in PTra and those implemented in Geant4-DNA, the simulation results of both codes were analysed using three physical parameters describing the ionisation cluster size distribution (ICSD): mean ionisation cluster size, variance of the cluster size and the probability to obtain a single ionisation within the target. Statistical analysis of the results indicates that the interaction cross section models for nitrogen used in PTra (elastic scattering, impact ionisations and electronic excitations) have been successfully implemented in Geant4-DNA. In addition, simulated ICSDs were compared to those measured with the Jet Counter nanodosimeter for energies between 100 and 2000 eV. For greater energies, the ICRP data for LET and particle range were used as a reference. The modified Geant4-DNA code and data successfully passed all these benchmarks fulfilling the requirement for their public release in the next version of the Geant4 toolkit.

Nanodosimetric Calculations of Radiation-Induced DNA Damage in a New Nucleus Geometrical Model Based on the Isochore Theory.

Double-strand breaks (DSBs) in nuclear DNA represents radiation-induced damage that has been identified as particularly deleterious. Calculating this damage using Monte Carlo track structure modeling could be a suitable indicator to better assess and anticipate the side-effects of radiation therapy. However, as already demonstrated in previous work, the geometrical description of the nucleus and the DNA content used in the simulation significantly influence damage calculations. Therefore, in order to obtain accurate results, this geometry must be as realistic as possible. In this study, a new geometrical model of an endothelial cell nucleus and DNA distribution according to the isochore theory are presented and used in a Monte Carlo simulation chain based on the Geant4-DNA toolkit. In this theory, heterochromatin and euchromatin compaction are distributed along the genome according to five different families (L1, L2, H1, H2, and H3). Each of these families is associated with a different hetero/euchromatin rate related to its compaction level. In order to compare the results with those obtained using a previous nuclear geometry, simulations were performed for protons with linear energy transfers (LETs) of 4.29 keV/µm, 19.51 keV/µm, and 43.25 keV/µm. The organization of the chromatin fibers at different compaction levels linked to isochore families increased the DSB yield by 6-10%, and it allowed the most affected part of the genome to be identified. These new results indicate that the genome core is more radiosensitive than the genome desert, with a 3-8% increase in damage depending on the LET. This work highlights the importance of using realistic distributions of chromatin compaction levels to calculate radio-induced damage using Monte Carlo simulation methods.

Implementation of the EPICS2017 database for photons in Geant4.

This paper describes in detail the implementation of Geant4 Livermore electromagnetic physics models based on the EPICS2017 database for the low energy transport of photons. These models describe four photon processes: gamma conversion, Compton scattering, photoelectric effect and Rayleigh scattering. New parameterizations based on EPICS2017 were performed for scattering functions of Compton effect, subshell cross-sections of the photoelectric effect and form factors of Rayleigh scattering, in order to improve the precision of fitted values compared to tabulated values. Comparisons between new and old parameterizations were also carried out to evaluate the precision of the new parameterizations. The models were tested through a comparative study, in which the mass attenuation coefficient was calculated for both total photon interaction and each process using Geant4 simulations based on EPICS2017 and EPDL97 respectively. The results obtained from the simulations were found in good agreement with the XCOM reference data.

Geant4-DNA Modeling of Water Radiolysis beyond the Microsecond: An On-Lattice Stochastic Approach.

In this work, we use the next sub-volume method (NSM) to investigate the possibility of using the compartment-based ("on-lattice") model to simulate water radiolysis. We, first, start with a brief description of the reaction-diffusion master equation (RDME) in a spatially discretized simulation volume ("mesh"), which is divided into sub-volumes (or "voxels"). We then discuss the choice of voxel size and merging technique of a given mesh, along with the evolution of the system using the hierarchical algorithm for the RDME ("hRDME"). Since the compartment-based model cannot describe high concentration species of early radiation-induced spurs, we propose a combination of the particle-based step-by-step ("SBS") Brownian dynamics model and the compartment-based model ("SBS-RDME model") for the simulation. We, finally, use the particle-based SBS Brownian dynamics model of Geant4-DNA as a reference to test the model implementation through several benchmarks. We find that the compartment-based model can efficiently simulate the system with a large number of species and for longer timescales, beyond the microsecond, with a reasonable computing time. Our aim in developing this model is to study the production and evolution of reactive oxygen species generated under irradiation with different dose rate conditions, such as in FLASH and conventional radiotherapy.

Geant4-DNA simulation of the pre-chemical stage of water radiolysis and its impact on initial radiochemical yields.

This paper demonstrates the impact of the pre-chemical stage, especially the dissociation scheme and the associated probabilities, on water radiolysis simulation using the Geant4-DNA Monte Carlo track structure simulation toolkit. The models and parameters provided by TRACs have been collected and implemented into Geant4-DNA. In order to evaluate their influence on water radiolysis simulation, the radiochemical yields (G-values) are evaluated as a function of time and LET using the "chem6" Geant4-DNA example, and they are compared with published experimental and calculated data. The new pre-chemical models lead to a better agreement with literature data than the default pre-chemical models of Geant4-DNA, especially for OH radicals and H2O2. The revised chemistry constructor "G4EmDNAChemistry_option3" is available in Geant4-DNA version 10.7.

A benchmarking study of Geant4 for Auger electrons emitted by medical radioisotopes.

Auger emitting radioisotopes are of great interest in targeted radiotherapy because, once internalised in the tumour cells, they can deliver dose locally to the radiation sensitive targets, while not affecting surrounding cells. Geant4 is a Monte Carlo code widely used to characterise the physics mechanism at the basis of targeted radiotherapy. In this work, we benchmarked the modelling of the emission of Auger electrons in Geant4 deriving from the decay of 123I, 124I, 125I radionuclides against existing theoretical approaches. We also compared Geant4 against reference data in the case of 131Cs, which is of interest for brachytherapy. In the case of 125I and 131Cs, the simulation results are compared to experimental measurements as well. Good agreement was found between Geant4 and the reference data. As far as we know, this is the first study aimed to benchmark against experimental measurements the emission of Auger electrons in Geant4 for radiotherapy applications.

Monte Carlo Simulation of SARS-CoV-2 Radiation-Induced Inactivation for Vaccine Development.

Immunization with an inactivated virus is one of the strategies currently being tested towards developing a SARS-CoV-2 vaccine. One of the methods used to inactivate viruses is exposure to high doses of ionizing radiation to damage their nucleic acids. While gamma (γ) rays effectively induce lesions in the RNA, envelope proteins are also highly damaged in the process. This in turn may alter their antigenic properties, affecting their capacity to induce an adaptive immune response able to confer effective protection. Here, we modeled the effect of sparsely and densely ionizing radiation on SARS-CoV-2 using the Monte Carlo toolkit Geant4-DNA. With a realistic 3D target virus model, we calculated the expected number of lesions in the spike and membrane proteins, as well as in the viral RNA. Our findings showed that γ rays produced significant spike protein damage, but densely ionizing charged particles induced less membrane damage for the same level of RNA lesions, because a single ion traversal through the nuclear envelope was sufficient to inactivate the virus. We propose that accelerated charged particles produce inactivated viruses with little structural damage to envelope proteins, thereby representing a new and effective tool for developing vaccines against SARS-CoV-2 and other enveloped viruses.

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.

Fully integrated Monte Carlo simulation for evaluating radiation induced DNA damage and subsequent repair using Geant4-DNA.

Ionising radiation induced DNA damage and subsequent biological responses to it depend on the radiation's track-structure and its energy loss distribution pattern. To investigate the underlying biological mechanisms involved in such complex system, there is need of predicting biological response by integrated Monte Carlo (MC) simulations across physics, chemistry and biology. Hence, in this work, we have developed an application using the open source Geant4-DNA toolkit to propose a realistic "fully integrated" MC simulation to calculate both early DNA damage and subsequent biological responses with time. We had previously developed an application allowing simulations of radiation induced early DNA damage on a naked cell nucleus model. In the new version presented in this work, we have developed three additional important features: (1) modeling of a realistic cell geometry, (2) inclusion of a biological repair model, (3) refinement of DNA damage parameters for direct damage and indirect damage scoring. The simulation results are validated with experimental data in terms of Single Strand Break (SSB) yields for plasmid and Double Strand Break (DSB) yields for plasmid/human cell. In addition, the yields of indirect DSBs are compatible with the experimental scavengeable damage fraction. The simulation application also demonstrates agreement with experimental data of [Formula: see text]-H2AX yields for gamma ray irradiation. Using this application, it is now possible to predict biological response along time through track-structure MC simulations.

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.

Advances in modelling gold nanoparticle radiosensitization using new Geant4-DNA physics models.

Gold nanoparticles have demonstrated significant radiosensitization of cancer treatment with x-ray radiotherapy. To understand the mechanisms at the basis of nanoparticle radiosensitization, Monte Carlo simulations are used to investigate the dose enhancement, given a certain nanoparticle concentration and distribution in the biological medium. Earlier studies have ordinarily used condensed history physics models to predict nanoscale dose enhancement with nanoparticles. This study uses Geant4-DNA complemented with novel track structure physics models to accurately describe electron interactions in gold and to calculate the dose surrounding gold nanoparticle structures at nanoscale level. The computed dose in silico due to a clinical kilovoltage beam and the presence of gold nanoparticles was related to in vitro brain cancer cell survival using the local effect model. The comparison of the simulation results with radiobiological experimental measurements shows that Geant4-DNA and local effect model can be used to predict cell survival in silico in the case of x-ray kilovoltage beams.

Monte Carlo dosimetry of a realistic multicellular model of follicular lymphoma in a context of radioimmunotherapy.

PURPOSE: Small-scale dosimetry studies generally consider an artificial environment where the tumors are spherical and the radionuclides are homogeneously biodistributed. However, tumor shapes are irregular and radiopharmaceutical biodistributions are heterogeneous, impacting the energy deposition in targeted radionuclide therapy. To bring realism, we developed a dosimetric methodology based on a three-dimensional in vitro model of follicular lymphoma incubated with rituximab, an anti-CD20 monoclonal antibody used in the treatment of non-Hodgkin lymphomas, which might be combined with a radionuclide. The effects of the realistic geometry and biodistribution on the absorbed dose were highlighted by comparison with literature data. Additionally, to illustrate the possibilities of this methodology, the effect of different radionuclides on the absorbed dose distribution delivered to the in vitro tumor were compared. METHODS: The starting point was a model named multicellular aggregates of lymphoma cells (MALC). Three MALCs of different dimensions and their rituximab biodistribution were considered. Geometry, antibody location and concentration were extracted from selective plane illumination microscopy. Assuming antibody radiolabeling with Auger electron (125 I and 111 In) and ß- particle emitters (177 Lu, 131 I and 90 Y), we simulated energy deposition in MALCs using two Monte Carlo codes: Geant4-DNA with "CPA100" physics models for Auger electron emitters and Geant4 with "Livermore" physics models for ß- particle emitters. RESULTS: MALCs had ellipsoid-like shapes with major radii, r, of ~0.25, ~0.5 and ~1.3 mm. Rituximab was concentrated in the periphery of the MALCs. The absorbed doses delivered by 177 Lu, 131 I and 90 Y in MALCs were compared with literature data for spheres with two types of homogeneous biodistributions (on the surface or throughout the volume). Compared to the MALCs, the mean absorbed doses delivered in spheres with surface biodistributions were between 18% and 38% lower, while with volume biodistribution they were between 15% and 29% higher. Regarding the radionuclides comparison, the relationship between MALC dimensions, rituximab biodistribution and energy released per decay impacted the absorbed doses. Despite releasing less energy, 125 I delivered a greater absorbed dose per decay than 111 In in the r ~ 0.25 mm MALC (6.78·10-2 vs 6.26·10-2  µGy·Bq-1 ·s-1 ). Similarly, the absorbed doses per decay in the r ~ 0.5 mm MALC for 177 Lu (2.41·10-2  µGy·Bq-1 ·s-1 ) and 131 I (2.46·10-2  µGy·Bq-1 ·s-1 ) are higher than for 90 Y (1.98·10-2  µGy·Bq-1 ·s-1 ). Furthermore, radionuclides releasing more energy per decay delivered absorbed dose more uniformly through the MALCs. Finally, when considering the radiopharmaceutical effective half-life, due to the biological half-life of rituximab being best matched by the physical half-life of 177 Lu and 131 I compared to 90 Y, the first two radionuclides delivered higher absorbed doses. CONCLUSION: In the simulated configurations, ß- emitters delivered higher and more uniform absorbed dose than Auger electron emitters. When considering radiopharmaceutical half-lives, 177 Lu and 131 I delivered absorbed doses higher than 90 Y. In view of real irradiation of MALCs, such a work may be useful to select suited radionuclides and to help explain the biological effects.

Development of integrated prompt gamma imaging and positron emission tomography system for in vivo 3-D dose verification: a Monte Carlo study.

An accurate knowledge of in vivo proton dose distribution is key to fully utilizing the potential advantages of proton therapy. Two representative indirect methods for in vivo range verification, namely, prompt gamma (PG) imaging and positron emission tomography (PET), are available. This study proposes a PG-PET system that combines the advantages of these two methods and presents detector geometry and background reduction techniques optimized for the PG-PET system. The characteristics of the secondary radiations emitted by a water phantom by interaction with a 150 MeV proton beam were analysed using Geant4.10.00, and the 2-D PG distributions were obtained and assessed for different detector geometries. In addition, the energy window (EW), depth-of-interaction (DOI), and time-of-flight (TOF) techniques are proposed as the background reduction techniques. To evaluate the performance of the PG-PET system, the 3-D dose distribution in the water phantom caused by two proton beams of energies 80 MeV and 100 MeV was verified using 16 optimal detectors. The thickness of the parallel-hole tungsten collimator of pitch 8 mm and width 7 mm was determined as 200 mm, and that of the GAGG scintillator was determined as 30 mm, by an optimization study. Further, 3-7 MeV and 2-7 MeV were obtained as the optimal EWs when the DOI and both the DOI and TOF techniques were applied for data processing, respectively; the detector performances were improved by about 38% and 167%, respectively, compared with that when applying only the 3-5 MeV EW. In this study, we confirmed that the PG distribution can be obtained by simply combining the 2-D parallel hole collimator and the PET detector module. In the future, we will develop an accurate 3-D dose evaluation technique using deep learning algorithms based on the image sets of dose, PG, and PET distributions for various proton energies.