Modeling ultra-high dose rate electron and proton FLASH effect with the physicochemical approach

Abstract

Objective. A physicochemical model built on the radiochemical kinetic theory was recently proposed in (Labarbe et al ) to explain the FLASH effect. We performed extensive simulations to scrutinize its applicability for oxygen depletion studies and FLASH-related experiments involving both proton and electron beams. Approach. Using the dose and beam delivery parameters for each FLASH experiment, we numerically solved the radiochemical rate equations comprised of a set of coupled nonlinear ordinary differential equations to obtain the area under the curve (AUC) of radical concentrations. Main results. The modeled differences in AUC induced by ultra-high dose rates appeared to correlate well with the FLASH effect. (i) For the whole brain irradiation of mice performed in (Montay-Gruel et al ), the threshold dose rate values for memory preservation coincided with those at which AUC started to decrease much less rapidly. (ii) For the proton pencil beam scanning FLASH of (Cunningham et al ), we found linear correlations between radicals’ AUC and the biological endpoints: TGF-β1, leg contracture and plasma level of cytokine IL-6. (iii) Compatible with the findings of the proton FLASH experiment in (Kim et al ), we found that radicals’ AUC at the entrance and mid-Spread-Out Bragg peak regions were highly similar. In addition, our model also predicted ratios of oxygen depletion G-values between normal and UHDR irradiation similar to those observed in (Cao et al ) and (El Khatib et al ). Significance. Collectively, our results suggest that the normal tissue sparing conferred by UHDR irradiation may be due to the lower degree of exposure to peroxyl and superoxide radicals. We also found that the differential effect of dose rate on the radicals’ AUC was less pronounced at lower initial oxygen levels, a trait that appears to align with the FLASH differential effect on normal versus tumor tissues.