Objective: In FLASH radiotherapy (dose rates ≥ 40 Gy/s), a reduced normal tissue toxicity has been observed, while maintaining the same tumor control compared to conventional radiotherapy (dose rates ≤ 0.03 Gy/s). This protecting effect could not be fully explained yet. One assumption is that interactions between the chemicals of different primary ionizing particles, so-called inter-track interactions, trigger this outcome. In this work, we included inter-track interactions in Monte Carlo track structure simulations and investigated the yield of chemicals (G-value) produced by ionizing particles. 
Approach: For the simulations, we used the Monte Carlo toolkit TOPAS, in which inter-track interactions cannot be implemented without further effort. Thus, we developed a method enabling the simultaneous simulation of N original histories in one event allowing chemical species to interact with each other. To investigate the effect of inter-track interactions we analyzed the G-value of different chemicals using various radiation sources. We used electrons with an energy of 60 eV in different spatial arrangements as well as a 10 MeV and 100 MeV proton source. For electrons we set N between 1 and 60, for protons between 1 and 100.
Main results: In all simulations, the total G-value decreases with increasing N. In detail, the G-value for ●OH, H3O and eaq decreases with increasing N, whereas the G-value of OH- , H2O2 and H2 increases slightly. The reason is that with increasing N, the concentration of chemical radicals increases allowing for more chemical reactions between the radicals resulting in a change of the dynamics of the chemical stage.
Signifcance: Inter-track interactions resulting in a variation of the yield of chemical species, may be a factor explaining the FLASH effect. To verify this hypothesis, further simulations are necessary in order to evaluate the impact of varying G-values on the yield of DNA damages.