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.

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.