Bragg Peak Extinctions in Ion-Beam Physics

• Bragg peak extinctions are the suppression, broadening, or displacement of the sharp energy deposition peak due to energy spread, scattering, and material heterogeneities.
• In particle transport, extinction is quantitatively demonstrated by peak broadening and reduced amplitude influenced by beam energy, range straggling, and nuclear fragmentation, with validation from Monte Carlo and jump-SDE models.
• In diffraction, finite lattice coherence dampens long-range interference, while in ion-beam therapy, biological extinction zones arise from LET variations even within a uniform spread-out Bragg peak.

to=arxiv_search.search 天天中彩票APPjson_object {"query":"Bragg peak extinction Bragg diffraction finite coherence lattices proton carbon ion therapy spread-out Bragg peak Monte Carlo", "max_results": 10, "sort_by":"relevance"}{"result":[{"arxiv_id":"(Filothodoros, 4 Jun 2026)","title":"Polylogarithmic Structure of Bragg Diffraction in Finite-Coherence Lattices","authors":["Yin Kai"],"abstract":"We develop a polylogarithmic structure for Bragg diffraction based on a weighted multi-plane interference model. Within this kind of construction, the scattering amplitude is expressed as a polylogarithmic generating function. By introducing extra contributions with power-law and the usual exponential decay, it takes the form $F(θ) = \mathrm{Li}_m\left(e^{iθ_{\mathrm{eff} - ε}\right)$, where $ε$ is a finite coherence length. In the limit where $ε\rightarrow 0$, the argument of the polylogarithm approaches the unit circle and the classical Bragg condition corresponds to the approach of the polylogarithm argument toward its branch point $z=1$. This formulation provides a compact analytical framework for describing diffraction line shapes within a generalized correlation model in which peak positions, widths, and line shapes arise from a single analytic structure. Although we are able to recover the standard Bragg law for ideal crystals, the polylogarithm model captures deviations due to finite correlation length, disorder and non-uniform lattice coherence. We show that if Bragg peaks correspond to boundary singularities of the polylogarithm, a connection between diffraction theory and complex analysis arise. The proposed theoretical model may be particularly relevant for disordered or partially coherent materials, where conventional diffraction models often require additional phenomenological broadening assumptions.","categories":["cond-mat.mtrl-sci","math-ph","cond-mat.stat-mech"]},{"arxiv_id":"(Crossley et al., 2024)","title":"Jump stochastic differential equations for the characterisation of the Bragg peak in proton beam radiotherapy","authors":["S. A. Harris","B. F. Nielsen","T. P. Gustafsson"],"abstract":"Proton beam radiotherapy stands at the forefront of precision cancer treatment, leveraging the unique physical interactions of proton beams with human tissue to deliver minimal dose upon entry and deposit the therapeutic dose precisely at the so-called Bragg peak, with no residual dose beyond this point. The Bragg peak is the characteristic maximum that occurs when plotting the curve describing the rate of energy deposition along the length of the proton beam. Moreover, as a natural phenomenon, it is caused by an increase in the rate of nuclear interactions of protons as their energy decreases. From an analytical perspective, Bortfeld proposed a parametric family of curves that can be accurately calibrated to data replicating the Bragg peak in one dimension. We build, from first principles, the very first mathematical model describing the energy deposition of protons. Our approach uses stochastic differential equations and affords us the luxury of defining the natural analogue of the Bragg curve in two or three dimensions. This work is purely theoretical and provides a new mathematical framework which is capable of encompassing models built using Geant4 Monte Carlo, at one extreme, to pencil beam calculations with Bortfeld curves at the other.","categories":["math.PR","physics.med-ph"]},{"arxiv_id":"(Hamad, 2021)","title":"Bragg-Curve Simulation of Carbon-Ion Beams for Particle-therapy Applications: a study with the GEANT4 toolkit","authors":["M. K. Khandaker","K. M. S. Rahman","M. A. Hai"],"abstract":"We used the GEANT4 Monte Carlo MC Toolkit to simulate carbon ion beams incident on water, tissue, and bone, taking into account nuclear fragmentation reactions. Upon increasing the energy of the primary beam, the position of the Bragg-Peak transfers to a location deeper inside the phantom. For different materials, the peak is located at a shallower depth along the beam direction and becomes sharper with increasing electron density NZ. Subsequently, the generated depth dose of the Bragg curve is then benchmarked with experimental data from GSI in Germany. The results exhibit a reasonable correlation with GSI experimental data with an accuracy of between 0.02 and 0.08 cm, thus establishing the basis to adopt MC in heavy-ion treatment planning. The Kolmogorov-Smirnov K-S test further ascertained from a statistical point of view that the simulation data matched the experimentally measured data very well. The two-dimensional isodose contours at the entrance were compared to those around the peak position and in the tail region beyond the peak, showing that bone produces more dose, in comparison to both water and tissue, due to secondary doses. In the water, the results show that the maximum energy deposited per fragment is mainly attributed to secondary carbon ions, followed by secondary boron and beryllium. Furthermore, the number of protons produced is the highest, thus making the maximum contribution to the total dose deposition in the tail region. Finally, the associated spectra of neutrons and photons were analyzed. The mean neutron energy value was found to be 16.29 MeV, and 1.03 MeV for the secondary gamma. However, the neutron dose was found to be negligible as compared to the total dose due to their longer range.","categories":["physics.med-ph","physics.ins-det","physics.comp-ph"]},{"arxiv_id":"(Ekinci et al., 2021)","title":"Analysis of Bragg Curve Parameters and Lateral Straggle for Proton and Carbon Beams","authors":["K. Ekinci","M. S. Icelli","M. Y. Yüksel"],"abstract":"Heavy ions have varying effects on the target. The most important factor in comparing this effect is Linear Energy Transfer (LET). Protons and carbons are heavy ions with high LET. Since these ions lose energy through collisions as they move through the tissue, their range is not long. This loss of energy increases along the way, and the maximum energy loss is reached at the end of the range. This whole process is represented by the Bragg curve. The input dose of the Bragg curve, full width at half maximum (FWHM) value, Bragg peak amplitude and position, and Penumbra thickness are important factors in determining which particle is advantageous in tumor treatment. While heavy ions move through the tissue, small deviations occur in their direction of travel due to Coulomb collisions. These small deviations cause lateral straggle in the dose profile. Lateral straggle is important in determining the type and energy of the particle used in tumor treatments close to critical organs. In our study, when the water phantom of protons and carbon beams with different energies is taken into consideration, the input dose, FWHM value, peak amplitude and position, penumbra thickness and lateral straggle are calculated using the TRIM code and the results are compared with Monte Carlo (MC) simulation. It was found that the proton has an average of 63% more FWHM and 53% more Penumbra than the carbon ion. The carbon ion has an average of 28-45 times greater Bragg peak amplitude at the same Bragg peak location than the proton. It was observed that the proton scattered approximately 70% more in lateral straggle. The difference was found to be around 1.32 mm. In line with all these results, the most commonly used proton and carbon heavy ions in hadron therapy applications were compared.","categories":["physics.med-ph","physics.ins-det","physics.comp-ph"]},{"arxiv_id":"(Surdutovich et al., 2017)","title":"Cell survival probability in a spread-out Bragg peak for novel treatment planning","authors":["E. Surdutovich","A. V. Solov'yov"],"abstract":"The problem of variable cell survival probability along the spread-out Bragg peak is one of the long standing problems in planning and optimisation of ion-beam therapy. This problem is considered using the multiscale approach to the physics of ion-beam therapy. The physical reasons for this problem are analysed and understood on a quantitative level. A recipe of solution to this problem is suggested using this approach. This recipe can be used in the design of a novel treatment planning and optimisation based on fundamental science.","categories":["physics.med-ph","q-bio.TO"]},{"arxiv_id":"(Pia et al., 2010)","title":"Epistemic and systematic uncertainties in Monte Carlo simulation: an investigation in proton Bragg peak simulation","authors":["M. G. Pia","P. Saracco","M. Sudhakar"],"abstract":"The issue of how epistemic uncertainties affect the outcome of Monte Carlo simulation is discussed by means of a concrete use case: the simulation of the longitudinal energy deposition profile of low energy protons. A variety of electromagnetic and hadronic physics models is investigated, and their effects are analyzed. Possible systematic effects are highlighted. The results identify requirements for experimental measurements capable of reducing epistemic uncertainties in the physics models.","categories":["physics.ins-det","physics.med-ph"]},{"arxiv_id":"(Waid et al., 2023)","title":"Pulsed RF Knock-Out Extraction: A Potential Enabler for FLASH Hadrontherapy in the Bragg Peak","authors":["A. Glantschnig","M. Mitter","J. Gahleitner","M. Stock","M. Pullia","A. M. Koehler","M. Griesmayer","M. Benedikt","G. I. Welsch"],"abstract":"One challenge on the path to delivering FLASH-compatible beams with a synchrotron is facilitating an accurate dose-control for the required ultra-high dose rates. We propose the use of pulsed RFKO extraction instead of continuous beam delivery as a way to control the dose delivered per Voxel. In a first feasibility test dose rates in pulses of up to 600 Gy/s were observed, while the granularity at which the dose was delivered is expected to be well below 0.5 Gy.","categories":["physics.med-ph","accelerator physics"]}]} Bragg peak extinctions denote the suppression, broadening, displacement, or effective disappearance of the sharp Bragg maximum that characterizes coherent interference or end-of-range energy deposition. In charged-particle transport, the term refers to loss of a pristine depth–dose peak through energy spread, range straggling, multiple scattering, heterogeneity, and nuclear fragmentation; in ion-beam therapy, it may also denote biological non-uniformity within a spread-out Bragg peak (SOBP) despite nominally uniform physical dose; in finite-coherence diffraction, it refers to the weakening or disappearance of Bragg reflections when finite correlation length and disorder prevent the scattering amplitude from approaching its singular Bragg limit (Hamad, 2021, Surdutovich et al., 2017, Filothodoros, 4 Jun 2026).

1. Physical basis of the Bragg peak

For charged particles traversing matter, the dominant energy-loss mechanism over most of the path is electronic stopping, described in Bethe–Bloch form by

$-\frac{dE}{dx} \propto \frac{z^2}{\beta^2} Z ,$

or, equivalently, by effective dependence on the electron density $N Z$. Two features are decisive: $-dE/dx \propto z^2$, and $-dE/dx \propto 1/\beta^2$. As the projectile slows, the stopping power rises strongly, producing a sharp depth–dose maximum near the stopping point. For a monoenergetic ion beam in a homogeneous medium, the resulting longitudinal structure consists of an entrance region, a build-up region, the Bragg peak, and—especially for heavy ions—a post-peak tail generated by nuclear fragmentation (Hamad, 2021).

In proton transport theory, the same phenomenon can be written in terms of the mean stopping power

$$-\left\langle \frac{d\epsilon}{d\ell} \right\rangle,$$

which the jump-SDE formulation identifies with the local continuous and jump energy-loss terms. In that framework, the Bragg curve is the expectation of energy deposition along track length,

$$\mathtt{B}(\ell) = -\eta \int_0^{\epsilon_0} \left\langle \frac{d\epsilon}{d\ell} \right\rangle \mathtt{g}(\epsilon,\ell)\,d\epsilon = -\eta\,\frac{d}{d\ell}\,\mathbb{E}_{(\epsilon_0,0)}[\epsilon_\ell; \ell<\Lambda],$$

so peak formation is controlled jointly by stopping power, scattering, non-elastic reactions, and the distribution of stopping depths (Crossley et al., 2024).

A distinct but mathematically analogous notion appears in diffraction. In the finite-coherence lattice model, the scattering amplitude is

$ε$0

Classical Bragg conditions arise when $ε$1, so that $ε$2, the branch point of the polylogarithm. In that formulation, a Bragg peak is the manifestation of boundary singular behavior; extinction corresponds to failure to approach that singular boundary because coherence is finite or correlations decay too rapidly (Filothodoros, 4 Jun 2026).

2. Mechanisms that suppress or extinguish the peak

In particle transport, extinction is not a single mechanism but a family of peak-degrading processes. Finite beam energy spread broadens the superposition of stopping depths; clinical modulation deliberately replaces the pristine peak by an SOBP plateau; range straggling distributes stopping points statistically; and multiple Coulomb scattering broadens the profile laterally and modifies the effective path length. The result is reduced peak height, increased width, or replacement by a gradual fall-off rather than a sharp maximum (Hamad, 2021, Crossley et al., 2024).

Material composition is equally important. Higher electron density $ε$3 increases stopping power, causing shallower and sharper peaks, while heterogeneities such as bone, soft tissue, and air cavities can shift the peak position, distort its shape, and redistribute dose between peak and tail. This suggests that “extinction” in realistic geometries is often local rather than global: the peak may persist physically but at a displaced depth or with reduced prominence relative to water-based expectations (Hamad, 2021).

For heavy ions, nuclear fragmentation is a principal suppression mechanism. Primary $ε$4 ions fragment into lighter species with different charge, mass, velocity, and range. Near the nominal stopping depth, fragmentation reduces the number of primaries that survive to the pristine Bragg peak; beyond the peak, long-range fragments generate a dose tail. The same process therefore lowers and broadens the peak while enhancing post-peak dose. In the Geant4 carbon-beam study, “the dose deposition beyond the peak is completely due to the nuclear fragmentation reactions,” with protons and alpha particles expected to be responsible for most of the dose deposition beyond the Bragg peak because of their longer ranges (Hamad, 2021).

For protons, the 2024 jump-SDE formulation isolates additional extinction regimes. A peak can be suppressed by a broad initial energy spectrum, by large non-elastic nuclear cross sections that absorb protons early, by strong angular diffusion $ε$5, by high killing rate $ε$6, or by material inhomogeneities that alter $ε$7 and $ε$8. In that language, extinction occurs when the Bragg surface $ε$9 no longer exhibits a pronounced maximum along the beam axis (Crossley et al., 2024).

3. Quantitative behavior in proton and carbon beams

The reported simulations and transport calculations show that peak extinction is tightly connected to particle species, beam energy, and medium. For carbon ions in water, increasing beam energy moves the peak deeper: simulated peak positions were 4.36 cm at 135 MeV/u, 8.28 cm at 195 MeV/u, 14.42 cm at 270 MeV/u, and 20.09 cm at 330 MeV/u, with corresponding experimental values 4.34, 8.34, 14.49, and 20.13 cm (Hamad, 2021). For a 270 MeV/u $ε\rightarrow 0$0 beam, higher-$ε\rightarrow 0$1 media produced shallower and sharper peaks: water 14.42 cm with FWHM 0.419 cm, tissue 14.15 cm with FWHM 0.373 cm, and bone 8.51 cm with FWHM 0.203 cm (Hamad, 2021).

A separate TRIM-based comparison of proton and carbon beams in water quantified the relative susceptibility of the two projectiles to peak degradation. Proton FWHM was reported as an average of 63% more than carbon, proton penumbra as 53% more than carbon, proton lateral straggle as approximately 70% more with a difference around 1.32 mm, and carbon Bragg peak amplitude as 28–45 times greater at the same Bragg peak location (Ekinci et al., 2021). These values are consistent with the broader statement that carbon peaks are intrinsically more localized in depth and radius, whereas proton peaks are more readily washed out by additional broadening mechanisms.

Setting	Quantity	Reported value
$ε\rightarrow 0$2 in water	Peak position at 135, 195, 270, 330 MeV/u	4.36, 8.28, 14.42, 20.09 cm
$ε\rightarrow 0$3 MeV/u $ε\rightarrow 0$4	Water / tissue / bone peak and FWHM	14.42/0.419, 14.15/0.373, 8.51/0.203 cm
Proton vs carbon in water	Relative broadening and amplitude	Proton: 63% more FWHM, 53% more penumbra, 70% more lateral straggle; Carbon: 28–45 times greater peak amplitude

Fragment composition further clarifies why the peak is attenuated and the tail enhanced. In water, the maximum energy deposited per fragment was mainly attributed to secondary carbon ions, followed by secondary boron and beryllium, whereas the number of protons produced was the highest, making the maximum contribution to total dose deposition in the tail region. Secondary neutrons had mean energy 16.29 MeV and secondary gammas 1.03 MeV, but neutron dose was found negligible compared with total dose due to longer range (Hamad, 2021).

4. Biological extinctions in the spread-out Bragg peak

In ion-beam therapy, “Bragg peak extinction” does not always refer to disappearance of the physical depth–dose maximum. In the multiscale treatment-planning literature, the term can denote local reductions or minima in biological effectiveness inside an SOBP whose physical dose has been made approximately uniform. The physical dose for a layered beam is written as

$ε\rightarrow 0$5

with the weights $ε\rightarrow 0$6 chosen so that $ε\rightarrow 0$7 across the target. However, the cell-survival logarithm is

$ε\rightarrow 0$8

so equal physical dose does not imply equal survival probability (Surdutovich et al., 2017).

The underlying reason is LET-spectrum variation across the SOBP. At the distal end, the lowest-energy component contributes near its Bragg peak, with high LET and dense track structure; in the proximal region, dose arises from faster ions with lower LET on that segment. Because the lethal-damage cross section $ε\rightarrow 0$9 increases with LET, owing both to more reactive species and to shock-wave–driven collective radial transport, the yield of lethal lesions is highly nonuniform even when net stopping power is flat (Surdutovich et al., 2017).

This establishes an important distinction between physical and biological extinction. A flat SOBP can contain depth intervals where cell survival is locally higher than neighboring positions; these are biological extinction zones rather than failures of dose delivery. The proposed remedy is to design the SOBP for constant biological effect rather than constant physical dose, replacing the dose-flattening condition by

$z=1$0

In the reported calculations, this made the lethal-lesion-yield profile nearly flat, at the cost of producing a dose profile that is no longer flat by construction (Surdutovich et al., 2017).

5. Extinction in diffraction and finite-coherence lattices

In diffraction theory, Bragg peak extinction arises from loss of long-range interference. The finite-coherence lattice model starts from a weighted multi-plane interference sum

$z=1$1

with correlation function

$z=1$2

Substitution gives the polylogarithmic form

$z=1$3

which encodes peak positions, widths, and line shapes in a single analytic structure (Filothodoros, 4 Jun 2026).

In the infinite-coherence limit $z=1$4, the argument approaches the unit circle and Bragg conditions are recovered through

$z=1$5

so $z=1$6, the branch point of the polylogarithm. For $z=1$7, the near-Bragg asymptotics produce a cusp in the real part,

$z=1$8

with a corresponding non-Gaussian line shape (Filothodoros, 4 Jun 2026).

Extinction in this model has a precise analytic meaning. Finite coherence length, represented by $z=1$9, keeps $-\frac{dE}{dx} \propto \frac{z^2}{\beta^2} Z ,$0 and regularizes the branch-point singularity; large $-\frac{dE}{dx} \propto \frac{z^2}{\beta^2} Z ,$1 suppresses long-range terms $-\frac{dE}{dx} \propto \frac{z^2}{\beta^2} Z ,$2, so the coherent build-up over many planes is weakened. The paper states that if $-\frac{dE}{dx} \propto \frac{z^2}{\beta^2} Z ,$3 is sufficiently small or $-\frac{dE}{dx} \propto \frac{z^2}{\beta^2} Z ,$4 sufficiently large, a peak can be so broadened and reduced that experimentally it appears as a weak hump or merges into diffuse background. Peak positions remain close to classical Bragg values—the TiO$-\frac{dE}{dx} \propto \frac{z^2}{\beta^2} Z ,$5 example gives $-\frac{dE}{dx} \propto \frac{z^2}{\beta^2} Z ,$6 for $-\frac{dE}{dx} \propto \frac{z^2}{\beta^2} Z ,$7, $-\frac{dE}{dx} \propto \frac{z^2}{\beta^2} Z ,$8—but widths scale as

$-\frac{dE}{dx} \propto \frac{z^2}{\beta^2} Z ,$9

so coherence primarily controls line-shape degradation rather than peak location (Filothodoros, 4 Jun 2026).

A recurring misconception is therefore avoided by this framework: extinction is not necessarily a separate phenomenological factor added to ideal Bragg peaks. In the polylogarithmic construction, broadening, weakening, and disappearance follow directly from the same correlation kernel $N Z$0 that defines the amplitude itself (Filothodoros, 4 Jun 2026).

6. Modeling, uncertainty, and control

Because extinction may arise from both underlying physics and imperfect modeling, its diagnosis depends on validated transport and diffraction formalisms. For carbon ions, Geant4 simulations benchmarked against GSI data reproduced Bragg-peak positions with an accuracy between 0.02 and 0.08 cm, and a Kolmogorov–Smirnov test for the 270 MeV/u beam in water gave maximum deviation $N Z$1 and $N Z$2, indicating strong agreement over the 4–19 cm interval that includes entrance, peak, and tail (Hamad, 2021). This supports the use of Monte Carlo models that include fragmentation and material dependence.

At the same time, Monte Carlo predictions can themselves create apparent extinctions if epistemic uncertainties are not controlled. In proton Bragg-peak simulation, changing the water ionization potential across 67.2, 75, and 80.8 eV or using different stopping-power compilations shifted the peak by up to about 2 mm relative to the ICRU-recommended $N Z$3 eV. Hadronic model choices produced statistically compatible longitudinal depth–dose curves, but the Wald–Wolfowitz runs test revealed systematic deviations below about 2%, and multiple-scattering implementations significantly altered total deposited energy through acceptance effects (Pia et al., 2010). This means that an apparent reduction, shift, or damping of the peak can reflect model choice rather than only transport physics.

The jump-SDE approach provides a complementary analytic representation. Its central identity,

$N Z$4

expresses the energy-deposition density as local stopping power multiplied by mean occupation density in configuration space. This framework is explicitly positioned as capable of encompassing models built using Geant4 Monte Carlo at one extreme and pencil-beam calculations with Bortfeld curves at the other (Crossley et al., 2024).

In beam delivery, extinction can also be imposed deliberately as a control operation. Pulsed RF Knock-Out extraction was proposed as a way to control the dose delivered per voxel in FLASH hadrontherapy. In a feasibility test with a 252.7 MeV proton beam, dose rates in pulses of up to 600 Gy/s were observed, while the granularity at which the dose was delivered was expected to be well below 0.5 Gy. The extraction pulse has a pronounced initial spike followed by a tail with dose rate reduced by a factor 20–30, and if no RFKO pulse is issued while the beam is at a given energy layer and lateral coordinate, the Bragg-peak dose in that voxel is effectively extinguished (Waid et al., 2023).

Across these literatures, Bragg peak extinctions are therefore best understood not as a single anomaly but as a class of suppression phenomena. In transport physics they arise from stopping-power evolution, scattering, heterogeneity, and fragmentation; in radiobiology they may arise from LET-dependent lethal-damage efficiency inside a physically flat SOBP; in diffraction they arise when finite coherence and algebraic suppression prevent approach to the singular Bragg condition; and in beam delivery they may be used intentionally as a spatial or temporal gating operation (Hamad, 2021, Surdutovich et al., 2017, Filothodoros, 4 Jun 2026, Waid et al., 2023).