FLASH Radiotherapy: Advances & Challenges
- FLASH-RT is defined by delivering therapeutic radiation at ultra-high dose rates (≥40 Gy/s) in sub-second pulses, offering potential normal-tissue sparing while maintaining tumor control.
- Research spans radiobiology, accelerator physics, dosimetry, and treatment planning, using metrics like instantaneous, average, and dose-averaged dose rates to optimize the FLASH effect.
- Advanced methodologies, including Monte Carlo simulations, real-time beam monitoring, and nonconvex optimization algorithms, are central to refining clinical translation and quality assurance in FLASH-RT.
Searching arXiv for recent and foundational papers on FLASH radiotherapy to ground the article in current literature. arXiv query: FLASH radiotherapy proton review dose rate optimization commissioning dosimetry mechanisms FLASH radiotherapy (FLASH-RT) is the ultra-high dose-rate delivery of a therapeutic radiation dose in a fraction of a second, with a commonly used operational threshold of and, in several formulations, single large doses delivered in sub-second time windows (Boucher et al., 2022). In preclinical literature, this delivery regime is associated with the “FLASH effect,” namely reduced normal-tissue toxicity with apparently preserved tumor control in many models, although the exact physical and biological conditions that trigger the effect remain incompletely characterized and the underlying mechanism is still unclear (Ma et al., 2024). Contemporary FLASH-RT research therefore spans radiobiology, accelerator physics, dosimetry, treatment planning, quality assurance, and translational beam-delivery engineering.
1. Definition, dose-rate criteria, and scope
FLASH-RT is conventionally distinguished from conventional dose-rate radiotherapy by both average and instantaneous dose-rate scales. One formulation gives conventional RT as and FLASH-RT as , with instantaneous dose rates , total dose per fraction , and total delivery time (Boucher et al., 2022). In proton and scanned-particle work, a voxelwise FLASH criterion is often defined by a joint dose and dose-rate condition. For example, one proton Bragg peak formulation treats a voxel as FLASH-qualified only if both and (Zhang et al., 9 May 2025).
Several dose-rate metrics coexist. In pencil-beam-scanning proton therapy, the literature distinguishes instantaneous dose rate, average dose rate, dose-averaged dose rate (DADR), average dose rate (ADR), and dose-threshold dose rate (DTDR), and explicitly notes the need for standardization of PBS dose-rate definitions (Wei et al., 2022). This directly bears on the interpretation of “FLASH coverage,” because different metrics can report materially different irradiated volumes as FLASH-compliant.
A common misconception is that FLASH-RT is defined by a single universal threshold. The literature instead presents thresholds as modality- and assay-dependent operational criteria. In vitro studies summarized in proton SFRT-FLASH planning report tissue-sparing beginning at and 0 per fraction, whereas accelerator-oriented overviews describe current estimates in which large doses of 1 or more delivered in 2 or less produce normal-tissue sparing effects (Zhang et al., 9 May 2025). This suggests that “FLASH conditions” should be treated as an experimentally grounded multidimensional parameter space rather than a single scalar cutoff.
2. Physicochemical and biological hypotheses
Current mechanistic work organizes the FLASH effect around oxygen depletion, radical chemistry, DNA-damage signaling, mitochondrial pathways, and immune modulation. A minimal oxygen-depletion model writes
3
with empirical depletion constants reported as 4 in vitro and 5 in vivo in one review (Ma et al., 2024). However, the same review emphasizes that in many FLASH experiments the net 6 is only a few mmHg, too small to explain the full normal-tissue sparing by a simple oxygen-fixation argument alone (Ma et al., 2024).
Radical-recombination hypotheses address this limitation by making the relevant state variable the transient concentration of radiolysis products rather than only mean oxygen level. A representative model writes
7
with analogous kinetics for peroxyl radicals (Ma et al., 2024). A more explicit 2025 study adopts a two-compartment radiochemical model for peroxyl radicals and reports that under CONV-RT at 8 recombination is negligible (9) even at 0, whereas under FLASH-RT at 1 recombination efficiency reaches 2 at 3 delivered in 4; the recombination-versus-delivery-time curve yields a 5-maximal recombination at 6, corresponding to a mean dose-rate threshold of 7 (Zhang et al., 28 Aug 2025).
Temporal microstructure also matters. Monte Carlo IRT modeling and physicochemical measurements indicate that compared with multi-pulse irradiation, single-pulse mode with a width less than 8 of the radical lifetime, a repetition interval longer than the radical lifetime, and a dose exceeding 9 can lead to rapid radical consumption, reducing residual content; when the single-pulse dose exceeded 0, the overlap probability approached 1 (Sun et al., 28 Apr 2025). This suggests that average dose rate alone is insufficient to specify the relevant biological regime.
At the biological level, the literature describes reduced clustering of double-strand breaks, lower cytosolic dsDNA release, suppression of cGAS–STING signaling in normal intestine, altered mitochondrial signaling, modulation of immune response, and circulating lymphocyte sparing as candidate components of the downstream response cascade (Ma et al., 2024). In the intestinal-injury model that links radical kinetics to transcriptomics, FLASH-RT activated the NRF2 antioxidant pathway, suppressed ERK signaling, increased the GSH/GSSG ratio from 2 to 3, increased CAT activity from 4 to 5, reduced cleaved caspase-3 to 6 relative to CONV, and abolished the protection after NRF2 knockdown (Zhang et al., 28 Aug 2025).
Mechanistic controversy persists. In silico tumor-control modeling based on radiolytic oxygen depletion predicts systematically lower tumor control probability for FLASH-RT than for conventional RT when extrapolated to TCP curves, even though preclinical tumor-volume curves may appear iso-effective (González-Crespo et al., 2023). A plausible implication is that any clinically robust account of the FLASH effect likely requires multiple coupled mechanisms, not ROD alone.
3. Modalities, beam delivery, and platform-specific implementations
FLASH-RT is not a single beam modality. Electron systems remain the most mature preclinical and early translational platform, but proton, ion, synchrotron, and laser-driven implementations are under active development. Accelerator-oriented summaries describe electrons at 7, very-high-energy electrons at 8, photons requiring high beam power to reach 9 at isocenter, and proton shoot-through beams from 0 iso-cyclotrons with average beam current on the order of 1 for FLASH intensities (Boucher et al., 2022).
In PBS proton therapy, delivery constraints are dominated by beam current, minimum monitor units per spot, spot-switch overhead, and energy-layer switching. Reported room beam currents include values below 2 and up to approximately 3, minimum spot weights around 4 in Varian ProBeam FLASH modes, transverse switch times of about 5, and energy switching around 6 per layer for energy-degrading cyclotrons and greater than 7 for synchrotrons (Wei et al., 2022). These machine constraints explain the strong interest in single-energy transmission beams, pullback systems with universal range shifter and range compensator, and other strategies that avoid slow multi-layer delivery (Wei et al., 2022).
An important recent proton development is the explicit combination of FLASH-RT with spatially fractionated radiation therapy (SFRT). One treatment-planning study proposes “SFRT-FLASH” in two forms: pGRID-FLASH implemented as scissor-beam FLASH (“SB-FLASH”) and proton minibeam-FLASH (“MB-FLASH”). The rationale is depth complementarity: FLASH-RT achieves ultra-high dose-rate sparing mostly in deep tissue near the Bragg peak, whereas SFRT with protons achieves high peak-to-valley dose ratio sparing in shallow-to-intermediate entrance regions (Zhang et al., 9 May 2025). Across four anatomical sites, MB-FLASH and SB-FLASH achieved high FLASH effect coverage of approximately 8 in the CTV1cm ring while preserving PVDR values of about 9 at shallow-to-intermediate depths (Zhang et al., 9 May 2025).
Research beamlines have also been adapted specifically for mechanistic and preclinical FLASH studies. At the Bern Medical Cyclotron, an 18 MeV proton beamline was modified to support both conventional and FLASH regimes spanning dose rates from 0 to 1 in passive-scattered configuration, with minibeam collimators supporting patterns such as 2 slits with 3 pitch and PVDR for the 4 grid falling from approximately 5 at 6 separation to approximately 7 at 8 (Kasanda et al., 6 May 2026). This kind of platform is central for controlled exploration of dose-rate effects, LET interplay, and SFRT parameters.
At the extreme end of temporal compression, laser-driven proton accelerators deliver ultra-high instantaneous dose rates in ultrashort bunches. In the first in vivo investigation of normal tissue response to laser-driven proton irradiation, each proton bunch delivered 9 in an 0 full width at 1 amplitude, giving an instantaneous dose rate of approximately 2 when accounting for 3 of the pulse dose in that interval, while the mean dose rate remained approximately 4 because one bunch was delivered every 5 (Obst-Huebl et al., 24 Feb 2026). This sharp separation of instantaneous and mean dose-rate scales is relevant to ongoing debates over which temporal variable is biologically determinant.
4. Treatment planning and optimization
Because FLASH-RT introduces dose-rate objectives in addition to dose objectives, treatment planning departs from standard DVH-only optimization. A proton planning formulation for SFRT-FLASH defines fraction-wise spot weights 6, dose-influence matrices 7, fraction doses 8, and total dose 9, then solves
0
subject to minimum-monitor-unit constraints and fieldwise FLASH coverage in a CTV1cm ring (Zhang et al., 9 May 2025). The CTV1cm structure is a 1 isotropic expansion of the CTV excluding the CTV interior, used as a surrogate organ-at-risk to enforce FLASH-rate objectives (Zhang et al., 9 May 2025).
Clinical trade-offs are explicit. In the same study, conformity index to CTV decreased under MB-FLASH and SB-FLASH; an example given for head-and-neck is CI changing from 2 in conventional planning to 3 in MB-FLASH, reflecting dose-rate and MMU constraints (Zhang et al., 9 May 2025). However, the effective FLASH-modified dose using a dose-modifying factor for FLASH restored plan conformity, with 4 (Zhang et al., 9 May 2025). This is not a proof of biological equivalence, but it shows how biological-effect modeling can materially alter the interpretation of physical-dose conformity.
For scanned proton delivery, spot order is itself an optimization variable because voxelwise dose rate depends on the temporal ordering of spot delivery. A TSP-based heuristic represents each proton spot as a node in a complete graph and optimizes the delivery permutation to maximize FLASH-rate voxels in a region of interest. In a 26-patient prostate cohort, the approximate TSP heuristic required 5 per beam, compared with an average of 6 for the simulated-annealing global optimizer, while recovering essentially the same FLASH volume: 7 versus 8 for ROI = patient minus prostate (Wase et al., 2024). Since the method does not alter spot weights or the nominal dose distribution, it can be applied as a post-processing step to an existing PBS plan (Wase et al., 2024).
FLASH optimization has also motivated new nonconvex and nonsmooth algorithms. A stochastic three-operator splitting formulation rewrites FLASH planning as minimization of
9
where 0, 1 is the indicator of the nonconvex zero-or-2 minimum-MU set, and 3 is the indicator of the dose/dose-rate feasibility set (Bian et al., 2023). In reported lung and brain patient cases, STOS improved conformity index and dose-rate coverage relative to a convex-relaxation ADMM baseline; for the brain case, 4 increased from 5 to 6 (Bian et al., 2023).
A further planning complication is that the toxicity gain from FLASH may depend on organ architecture when dose distributions are heterogeneous. An in silico NTCP study using both a radiolytic oxygen depletion model and a phenomenological logistic model found that sparing increased with decreasing LKB volume-effect parameter 7, meaning more important sparing for serial organs. In one specific calculation with conventional NTCP fixed at 8, the corresponding FLASH NTCP ranged from 9 for 0 to 1 for 2 (Pardo-Montero, 28 Apr 2026). This suggests that local high-dose voxels may be disproportionately relevant for FLASH benefit in serial structures.
5. Dosimetry, commissioning, and quality assurance
FLASH dosimetry is difficult because detectors developed for conventional radiotherapy can saturate, become nonlinear, or lose temporal fidelity under high dose per pulse and high instantaneous dose rate. A review of FLASH dosimetry describes ion chambers as suffering recombination when 3, with charge-collection efficiency losses greater than 4 unless high bias and complex corrections are applied (Ashraf et al., 2020). The same review emphasizes the appeal of scintillation and Cherenkov methods because of nanosecond-scale response, sub-millimeter spatial resolution, and dose-rate independence up to at least 5 (Ashraf et al., 2020).
Commissioning frameworks for electron FLASH units therefore include non-standard beam parameters such as pulse width, pulse repetition frequency, dose per pulse, and instantaneous dose rate. A general protocol for acceptance testing, commissioning, and routine QA specifies, among other criteria, reproducibility of 6 consecutive irradiations at 7, dose-monitoring proportionality versus pulse number at 8, 9 independence versus PRF at 00, 01 linearity versus pulse width at 02, angular output stability at 03, and PDD parameters with 04 within 05 and 06 within 07 or 08 (Palmiero et al., 2024). The same framework advocates redundant dosimetry using ion chambers, beam current transformers, and passive dose-rate-independent dosimeters such as EBT3 film, TLD, and OSLD (Palmiero et al., 2024).
Monte Carlo beam models are emerging as the practical basis for electron FLASH treatment planning because commercial treatment-planning systems for electron FLASH radiotherapy are unavailable. For the Mobetron UHDR system at 09, pulse-width-specific phase-space modeling showed that mean energy decreased exponentially from 10 to 11 as pulse width increased from 12 to 13, while energy spread increased quadratically, both with 14 (Carballeira et al., 6 May 2026). A universal reference pulse width of 15 yielded 16 mean energy and reduced computational burden by 17 while maintaining clinical accuracy within AAPM TG-106 tolerances (Carballeira et al., 6 May 2026).
Real-time beam monitoring is another central requirement. A point scintillator detector coupled to a gated amplifier and FPGA controller was shown to be linear with mean dose rate from 18 to 19 and dose per pulse from 20 to 21 to within 22, enabling pulse-resolved dose integration and dose-based gating of a modified LINAC (Ashraf et al., 2021). The same system revealed a ramp-up of 23 pulses in vivo, during which average dose per pulse was only about 24 before stabilizing at 25 (Ashraf et al., 2021). This is important because early underdose during ramp-up may compromise the intended FLASH regime for low total doses.
Large-area two-dimensional scintillator monitors address complementary QA needs. A prototype FLASH Beam Scintillator Monitor using a proprietary inorganic hybrid scintillator achieved a small 26 signal decrease after a cumulative dose of 27, dose-per-pulse linearity with 28, spatial-resolution residual RMS of 29 in one configuration, and FPGA beam-parameter computation in less than 30 after the end of frame at 31 (Levin et al., 2023). Flexible 2D scintillating coatings further extend this approach to non-homogeneous preclinical electron beams and curved ex vivo surfaces, with linear response up to at least 32 pulses and setup robustness within 33, though signal becomes noisy for field sizes at or below about 34 (Vanreusel et al., 22 Apr 2025).
6. Experimental evidence, controversies, and translational outlook
The empirical basis for FLASH-RT remains strongest for normal-tissue sparing in preclinical systems. A white-paper synthesis cites reduced lung fibrosis, improved neurocognitive outcomes, and reduced skin toxicity in several electron-beam models, with normal tissues “universally spared” in preclinical studies while tumors are not (Boucher et al., 2022). More recent mechanistic work sharpens the dependence of this sparing on dose and delivery conditions. In a murine whole-abdomen model, a total-dose series at fixed mean dose rate of 35 found that only at 36 did FLASH-RT significantly reduce histological scores compared to CONV-RT, and a mean dose-rate series at fixed dose of 37 found the FLASH effect emerging at 38, deepening up to 39 and saturating thereafter (Zhang et al., 28 Aug 2025).
The same study also demonstrated that hybrid schedules can preserve or abolish sparing depending on the FLASH component: at total dose 40, FLASH component 41 retained sparing, whereas 42 FLASH plus 43 conventional did not (Zhang et al., 28 Aug 2025). Objective injury markers tracked these schedule dependencies. For example, tissue MDA at 44 after 45 was 46 for C15 and 47 for F15, a 48 reduction; TUNEL-positive crypt cells at 49 were reduced from 50 to 51 (Zhang et al., 28 Aug 2025).
Proton evidence is expanding but remains heterogeneous across delivery platforms. The Bern 18 MeV beamline supports systematic pre-clinical proton radiobiology studies over five decades of dose rate, while the first in vivo laser-driven proton study reported reduced tissue swelling compared with conventional-dose-rate X-rays and distinct RNA-sequencing signatures after ultra-high instantaneous dose-rate proton irradiation (Kasanda et al., 6 May 2026). In the laser-driven mouse ear model, peak swelling after 52 protons was 53 versus 54 for 55 X-rays, approximately 56 lower (Obst-Huebl et al., 24 Feb 2026). However, the authors also note that definitive attribution of a FLASH effect in that regime would require conventional-dose-rate protons with matching LET as the true dosimetric control (Obst-Huebl et al., 24 Feb 2026).
Several controversies follow from these data. First, the exact molecular targets and the weighting of oxygen depletion, radical recombination, immune modulation, and mitochondrial or nuclear signaling remain unresolved (Ma et al., 2024). Second, the clinically relevant dose-rate descriptor is unsettled, especially in pulsed and scanned beams where mean dose rate, instantaneous dose rate, dose per pulse, and temporal microstructure can dissociate strongly (Sun et al., 28 Apr 2025). Third, the assumption of tumor iso-effectiveness is not mechanistically secure: radiolytic oxygen depletion models predict lower TCP for FLASH-RT than for conventional RT unless compensating mechanisms not included in the model contribute in real tumors (González-Crespo et al., 2023).
The translational agenda implied by current work is therefore technically specific. It includes collimator-free minibeam generation to reduce monitor-unit burden in proton SFRT-FLASH, preclinical organ-specific validation of combined SFRT-FLASH biological effects, prospective clinical trials evaluating toxicity reduction and tumor control, verification of beam-model regressions on additional Mobetron installations, incorporation of universal and pulse-width-specific phase space files into Monte Carlo or GPU-accelerated dose engines, continued standardization of commissioning and QA for electron UHDR units, and end-to-end treatment-planning studies that jointly evaluate dose, dose rate, delivery accuracy, and normal-tissue sparing (Zhang et al., 9 May 2025).
FLASH-RT thus occupies a distinctive position in radiotherapy research. Its central claim is not merely geometric dose conformity, but biologically favorable temporal delivery. The field’s present state is defined by that promise, by substantial preclinical support for normal-tissue sparing, and by equally substantial uncertainty regarding the precise conditions under which the effect is robust, safe, and clinically generalizable.