MPEXS-DNA, a new GPU-based Monte Carlo simulator for track structures and radiation chemistry at subcellular scale

Track structure simulation codes can accurately reproduce the stochastic nature of particle-matter interactions in order to evaluate quantitatively radiation damage in biological cells such as DNA strand breaks and base damage. Such simulations handle large numbers of secondary charged particles and molecular species created in the irradiated medium. Every particle and molecular species are tracked step-by-step using a Monte Carlo method to calculate energy loss patterns and spatial distributions of molecular species inside a cell nucleus with high spatial accuracy. The Geant4-DNA extension of the Geant4 general-purpose Monte Carlo simulation toolkit allows for such track structure simulations and can be run on CPUs. However, long execution times have been observed for the simulation of DNA damage in cells. We present in this work an improvement of the computing performance of such simulations using ultraparallel processing on a graphical processing unit (GPU).

A new Monte Carlo simulator named MPEXS-DNA, allowing high computing performance by using a GPU, has been developed for track structure and radiolysis simulations at the subcellular scale. MPEXS-DNA physics and chemical processes are based on Geant4-DNA processes available in Geant4 version 10.02 p03. We have reimplemented the Geant4-DNA process codes of the physics stage (electromagnetic processes of charged particles) and the chemical stage (diffusion and chemical reactions for molecular species) for microdosimetry simulation by using the CUDA language. MPEXS-DNA can calculate a distribution of energy loss in the irradiated medium caused by charged particles and also simulate production, diffusion, and chemical interactions of molecular species from water radiolysis to quantitatively assess initial damage to DNA. The validation of MPEXS-DNA physics and chemical simulations was performed by comparing various types of distributions, namely the radial dose distributions for the physics stage, and the G-value profiles for each chemical product and their linear energy transfer dependency for the chemical stage, to existing experimental data and simulation results obtained by other simulation codes, including PARTRAC.

For physics validation, radial dose distributions calculated by MPEXS-DNA are consistent with experimental data and numerical simulations. For chemistry validation, MPEXS-DNA can also reproduce G-value profiles for each molecular species with the same tendency as existing experimental data. MPEXS-DNA also agrees with simulations by PARTRAC reasonably well. However, we have confirmed that there are slight discrepancies in G-value profiles calculated by MPEXS-DNA for molecular species such as H₂ and H₂ O₂ when compared to experimental data and PARTRAC simulations. The differences in G-value profiles between MPEXS-DNA and PARTRAC are caused by the different chemical reactions considered. MPEXS-DNA can drastically boost the computing performance of track structure and radiolysis simulations. By using NVIDIA's GPU devices adopting the Volta architecture, MPEXS-DNA has achieved speedup factors up to 2900 against Geant4-DNA simulations with a single CPU core.

The MPEXS-DNA Monte Carlo simulation achieves similar accuracy to Monte Carlo simulations performed using other codes such as Geant4-DNA and PARTRAC, and its predictions are consistent with experimental data. Notably, MPEXS-DNA allows calculations that are, at maximum, 2900 times faster than conventional simulations using a CPU.