FLASH radiotherapy is a method of irradiating tumors with ultra-short pulses at a high dose rate, which spares normal tissues while maintaining tumor control. The mechanisms behind this effect are actively debated: hypotheses include transient oxygen depletion, reduced yield of reactive oxygen species, radical recombination, and others. However, all biological responses begin with the absorption of radiation energy, and the spatial distribution of that energy is the first step toward understanding the phenomenon. This work aligns with studies modeling the influence of oxygen heterogeneity and oxygen depletion on the radiobiological response to FLASH irradiation but differs by addressing a purely physical aspect—nonlinear energy absorption in the material—without invoking specific biochemical pathways.
A two-dimensional model of a homogeneous medium was built on a 250×250 pixel grid with a constant energy diffusion coefficient. The diffusion equation with absorption was solved using an explicit finite-difference method with a time step of 0.05 (in arbitrary units) over 8000 steps until reaching a steady state. The energy source was placed in a central circular region of radius 70 pixels, and the measurement zone was an annulus from radius 70 to 100 pixels. Four absorption coefficient models were considered: constant, linearly increasing with energy, saturating (Hill function with exponent 0.5), and threshold (sigmoid). Two delivery modes of the same total energy were compared: continuous low-intensity (CONV-like) and short high-intensity pulse (FLASH-like). The effectiveness was quantified by the relative reduction in absorbed energy in the annular region (RER).
With a constant absorption coefficient, both modes produced identical energy distributions, confirming the correctness of the numerical scheme. For all nonlinear models, the FLASH-like mode reduced the total absorbed energy in the surrounding region compared to the continuous mode: absorption in the central region increased by about 1.5%, while in the annulus it decreased, with the effect being stronger for the linear and sigmoid models than for the sublinear Hill model. Analysis of the parameter space showed that the dependence of RER on pulse intensity is non-monotonic: at a low energy threshold (Eth), increasing intensity can even increase absorption in the surrounding region, while at a high threshold the benefit plateaus. This qualitatively matches experimental data showing that the FLASH effect weakens or disappears above a critical dose.
The work demonstrates that nonlinear energy absorption alone, without any biological mechanisms, can reproduce the qualitative picture of the tissue-sparing effect of FLASH therapy. Unlike neighboring studies that focus on oxygen depletion and radiation-chemical processes, this work proposes a minimal physical mechanism based on material absorption properties. This does not rule out other hypotheses but provides a simple, testable foundation for interpreting the heterogeneous response of different tissues and tumors.
The model can be extended by incorporating tissue heterogeneity, spatial variations in absorption and diffusion parameters, and more realistic nonlinear response functions. Comparing model predictions with experimental dose distribution and radiobiological outcome data for different tissue types is promising. A transition to three-dimensional geometry and inclusion of real parameters from measurements is possible.
The results may help optimize FLASH therapy protocols: tailoring pulse intensity and duration to the absorbing medium's characteristics (e.g., tumor type) could enhance treatment selectivity. The model also applies to other systems with nonlinear energy absorption, such as laser material processing.
Plans include adding spatial parameter variations (diffusion and absorption threshold) to simulate heterogeneous biological tissue and testing whether the qualitative conclusions hold. Also, comparing predictions with published experimental data on FLASH irradiation is intended.
In radiobiology, the physical mechanism of the FLASH effect remains an open question: why ultra-fast irradiation spares normal tissues but not tumors. This model offers a possible explanation through nonlinear energy absorption, which may differ between tissue types due to their composition and microstructure. This links radiation transport physics to observed biological selectivity.
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
This study presents a minimal model of energy transfer in a homogeneous medium with nonlinear absorption. Numerical simulations show that under pulsed FLASH-like irradiation, energy is redistributed: the central zone absorbs more energy, while the surrounding region absorbs less than under continuous irradiation with the same total dose. The effect depends on the type of nonlinearity and pulse parameters, providing a physical explanation for the tissue-sparing effect of FLASH therapy.
Where its strength is
The strength of this work lies in demonstrating that nonlinear energy absorption alone, without any biological mechanisms, reproduces the qualitative picture of the FLASH effect. This provides a minimal physical foundation onto which more complex models can be built.
Where it could go next
Worth trying
Why this may matter
If nonlinear absorption is indeed a fundamental physical mechanism behind the FLASH effect, it could simplify our understanding of the phenomenon and aid in optimizing irradiation parameters without resorting to complex biological models. Moreover, it opens the door to engineering solutions aimed at enhancing spatial energy redistribution.
Nearby works in our archive
The closest work in the archive sits at 0.70 by meaning — the usual level for our corpus.
The model used a 250×250 pixel grid — imagine a battlefield the size of a football field divided into 10-centimeter squares, with each square recording how much energy landed there.