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.

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.