The indirect effect of radiation plays an important role in radio-induced biological damages. Monte Carlo codes have been widely used in recent years to study the chemical evolution of particle tracks. However, due to the large computational efforts required, their applicability is typically limited to simulations in pure water targets and to temporal scales up to the µs. In this work, a new extension of TRAX-CHEM is presented, namely TRAX-CHEMxt, able to predict the chemical yields at longer times, with the capability of exploring the homogeneous biochemical stage. Based on the species coordinates produced around one track, the set of reaction-diffusion equations is solved numerically with a computationally light approach based on concentration distributions. In the overlapping time scale (500 ns-1 µs), a very good agreement to standard TRAX-CHEM is found, with deviations below 6% for different beam qualities and oxygenations. Moreover, an improvement in the computational speed by more than three orders of magnitude is achieved. The results of this work are also compared with those from another Monte Carlo-based algorithm and a fully homogeneous code (Kinetiscope). TRAX-CHEMxt will allow for studying the variation in chemical endpoints at longer timescales with the introduction, as the next step, of biomolecules, for more realistic assessments of biological response under different radiation and environmental conditions.
Keywords: chemical track structure; homogeneous biochemical stage; radical/molecule yields; reaction–diffusion equations.