Pion Bump: Astrophysical & Nuclear Signatures

• Pion bump is a term that denotes distinct spectral structures associated with pion production, decay, or propagation, with its definition varying by observational context.
• It serves as a diagnostic tool in gamma-ray astrophysics, enabling the identification of hadronic processes and candidate neutrino sources through its characteristic spectral turnover.
• In nuclear collisions and lepton–nucleus scattering, the bump highlights in-medium pion effects, resonance structures, and modifications in pion propagation.

Searching arXiv for papers on the pion bump across astrophysics, nuclear collisions, and lepton–nucleus scattering. “Pion bump” denotes several distinct spectral and cross-section structures associated with pion production, decay, or propagation, and its precise meaning is strongly context dependent. In high-energy astrophysics the term most commonly refers to the MeV–GeV gamma-ray turnover produced by $\pi^0 \to 2\gamma$, and by extension to the MeV spectral break used to select candidate hadronic neutrino sources. In nuclear and hadronic reaction studies, the same term is used for a peak in the kinetic-energy dependence of $\pi^-/\pi^+$, for a broad isoscalar enhancement near $\sqrt{s}\approx 2.31$ GeV in single- or double-pion production, and for the $\Delta(1232)$-region enhancement above the quasielastic peak in lepton–nucleus scattering (Yang et al., 2018, Granados et al., 30 Jul 2025, Feng, 2016, Clement et al., 2020, Isaacson et al., 26 Aug 2025).

1. Terminological scope and principal usages

The term is not attached to a single universal observable. Its meaning is fixed by the reaction channel, the plotted variable, and the underlying pion-production mechanism.

Context	Observable called “pion bump”	Characteristic scale
Hadronic gamma-ray astrophysics	$\pi^0$-decay turnover in $\gamma$-ray spectrum	$\sim 10$ MeV to a few GeV
Galactic neutrino source selection	Fermi-LAT spectral break interpreted as $\pi^0$-decay onset	$50$ MeV to $1$ GeV; often near $\pi^-/\pi^+$0 MeV
Heavy-ion collisions	Peak in $\pi^-/\pi^+$1	$\pi^-/\pi^+$2 MeV
Isoscalar single-/double-pion production	Broad enhancement in total cross section	$\pi^-/\pi^+$3–$\pi^-/\pi^+$4 GeV
Lepton–nucleus scattering	Inclusive enhancement above QE peak from single-$\pi^-/\pi^+$5 production	$\pi^-/\pi^+$6 GeV
Flaring blazars	Narrow $\pi^-/\pi^+$7-decay feature in VHE $\pi^-/\pi^+$8-ray spectrum	$\pi^-/\pi^+$9 TeV

This diversity is not merely terminological. In the astrophysical case the bump is a decay-kinematics signature of $\sqrt{s}\approx 2.31$0, whereas in heavy-ion and hadronic spectroscopy it is a resonance- and medium-structured feature in pion production or propagation. Precision usage therefore requires explicit specification of the observable.

2. The classical $\sqrt{s}\approx 2.31$1-decay bump in gamma-ray astrophysics

In hadronic cosmic-ray scenarios, high-energy protons interact with ambient gas through $\sqrt{s}\approx 2.31$2–$\sqrt{s}\approx 2.31$3 or $\sqrt{s}\approx 2.31$4–$\sqrt{s}\approx 2.31$5 collisions and produce pions. Neutral pions decay promptly via $\sqrt{s}\approx 2.31$6, while charged pions decay through $\sqrt{s}\approx 2.31$7 followed by $\sqrt{s}\approx 2.31$8 (Granados et al., 30 Jul 2025). The $\sqrt{s}\approx 2.31$9 mass, $\Delta(1232)$0, fixes the rest-frame photon energy to $\Delta(1232)$1. The threshold kinetic energy for $\Delta(1232)$2 is given as

$\Delta(1232)$3

while a complementary treatment describes the threshold as near $\Delta(1232)$4 (Granados et al., 30 Jul 2025, Yang et al., 2018).

Because astrophysical $\Delta(1232)$5 are produced with a distribution of energies and angles, the monochromatic $\Delta(1232)$6 MeV photons in the pion rest frame are Doppler broadened in the observer frame. Together with the production threshold, this generates a broad turnover or “bump” in the $\Delta(1232)$7-ray spectrum. One formulation describes the resulting feature in $\Delta(1232)$8 as a bell-type structure between about $\Delta(1232)$9 MeV and a few GeV; another emphasizes that in practice Fermi-LAT identifies it as a spectral break between $\pi^0$0 MeV and $\pi^0$1 GeV, with a characteristic turnover near $\pi^0$2 MeV (Yang et al., 2018, Granados et al., 30 Jul 2025). A standard hadronic emissivity representation is

$\pi^0$3

or equivalently through the intermediate pion distribution (Liu et al., 2024).

The detailed shape below the bump maximum is not unique to $\pi^0$4-decay alone. Secondary $\pi^0$5, produced through $\pi^0$6-meson decays, generate bremsstrahlung that can distort the spectrum below $\pi^0$7 MeV, and sub-relativistic heavy ions can contribute additional $\pi^0$8-ray flux in the same band (Yang et al., 2018). In dense, calorimetric environments, secondary bremsstrahlung can dominate below $\pi^0$9 MeV; at $\gamma$0 MeV it can exceed $\gamma$1-decay $\gamma$2-rays by an order of magnitude once $\gamma$3 (Yang et al., 2018).

The supernova remnant W44 provides a standard case study for the observational ambiguity. Using broken power laws in momentum, a hadronic fit employs proton indices $\gamma$4, $\gamma$5, and $\gamma$6, while a pure leptonic bremsstrahlung fit uses electron indices $\gamma$7, $\gamma$8, $\gamma$9, and a low-energy cutoff $\sim 10$0 (Liu et al., 2024). Both fit the GeV-band data, but they diverge strongly in the MeV band: the hadronic model predicts a sharp downturn across $\sim 10$1–$\sim 10$2 MeV, whereas the leptonic model remains bright and smooth. A MeGaT-like instrument with $\sim 10$3, $\sim 10$4 PSF, and a $\sim 10$5-month exposure is forecast to achieve $\sim 10$6–$\sim 10$7 bin-by-bin detections across $\sim 10$8–$\sim 10$9 MeV and to separate the models decisively (Liu et al., 2024).

3. The pion bump as a multimessenger hadronic tag

In Galactic-source neutrino searches, the pion bump is used not primarily as an end in itself but as a source-selection criterion. A recent IceCube analysis targets $\pi^0$0 Galactic Plane sources from the Fermi-LAT 4FGL catalog that exhibit the spectral break between $\pi^0$1 MeV and $\pi^0$2 GeV associated with the pion bump (Granados et al., 30 Jul 2025). These sources are treated as candidate hadronic emitters because the same hadronic interactions that generate $\pi^0$3-decay $\pi^0$4-rays also produce $\pi^0$5, and hence neutrinos.

The catalog contains $\pi^0$6 SNRs, $\pi^0$7 HMBs, $\pi^0$8 PWNe, $\pi^0$9 SFR, $50$0 SNR/PWN/composite sources, $50$1 binary, $50$2 unidentified sources, and $50$3 unknown sources; $50$4 of the $50$5 have a TeV counterpart within $50$6, including IC 443, W28, W49B, W51, MSH 15−52, HESS J1857+026, LSI+61 303, Eta Carinae, and the Cocoon (Granados et al., 30 Jul 2025). The IceCube search uses $50$7 years of data, combines track-like and cascade-like events while removing overlaps, and applies the standard unbinned point-source likelihood

$50$8

with test statistic

$50$9

Background is estimated by right-ascension scrambling with the Galactic Plane masked, using $1$0 scrambled pseudo-experiments per test (Granados et al., 30 Jul 2025).

The physical link from the $1$1-ray bump to neutrinos is standard but not one-to-one. For $1$2–$1$3 interactions at GeV–TeV energies, the approximate production ratio is $1$4, and after oscillations the flavor composition at Earth approaches $1$5 (Granados et al., 30 Jul 2025). For optically thin sources, the all-flavor neutrino flux and the $1$6-ray flux are approximately proportional at comparable energies, with a proportionality factor $1$7 of order unity to a few. However, a MeV bump only confirms hadronic $1$8-ray production at low energies. IceCube sensitivity is in the TeV–PeV range, and for $1$9–$\pi^-/\pi^+$00 interactions a neutrino typically carries $\pi^-/\pi^+$01–$\pi^-/\pi^+$02; TeV neutrinos therefore require proton acceleration to tens of TeV or higher (Granados et al., 30 Jul 2025).

The reported sensitivities are approximately two orders of magnitude below the diffuse Galactic Plane neutrino flux measured by IceCube in 2023, implying sensitivity to source populations contributing at the $\pi^-/\pi^+$03 level of the Galactic Plane emission. The analysis is explicitly framed as a search under active analysis: no detections, TS values, $\pi^-/\pi^+$04-values, or final upper limits are reported (Granados et al., 30 Jul 2025).

4. In-medium “pion bumps” in heavy-ion collisions

In heavy-ion transport theory, “pion bump” can refer to a different structure: a peak in the kinetic-energy dependence of the charged-pion ratio

$\pi^-/\pi^+$05

Within the Lanzhou Quantum Molecular Dynamics model, simulations of $\pi^-/\pi^+$06 at $\pi^-/\pi^+$07nucleon show a pronounced bump in $\pi^-/\pi^+$08 at $\pi^-/\pi^+$09, identified with the $\pi^-/\pi^+$10 resonance region (Feng, 2016).

Near threshold, pion production proceeds predominantly through $\pi^-/\pi^+$11, followed by $\pi^-/\pi^+$12, and is strongly coupled to $\pi^-/\pi^+$13 reabsorption cycles in dense matter (Feng, 2016). The key medium ingredient is an isospin-dependent pion–nucleon potential based on the $\pi^-/\pi^+$14-hole model. The in-medium pion energy is written as

$\pi^-/\pi^+$15

with $\pi^-/\pi^+$16 for $\pi^-/\pi^+$17, $\pi^-/\pi^+$18, and $\pi^-/\pi^+$19 (Feng, 2016). The $\pi^-/\pi^+$20-hole self-energy splits the pion mode into $\pi^-/\pi^+$21-like and $\pi^-/\pi^+$22-like branches, which cross near the $\pi^-/\pi^+$23 energy and generate a “pocket” in the optical potential $\pi^-/\pi^+$24. At $\pi^-/\pi^+$25 and $\pi^-/\pi^+$26, the quoted values are $\pi^-/\pi^+$27 for $\pi^-/\pi^+$28, $\pi^-/\pi^+$29 for $\pi^-/\pi^+$30, and $\pi^-/\pi^+$31 for $\pi^-/\pi^+$32; at $\pi^-/\pi^+$33 they become $\pi^-/\pi^+$34, $\pi^-/\pi^+$35, and $\pi^-/\pi^+$36, respectively (Feng, 2016).

This isospin splitting modifies pion propagation and reabsorption differently for $\pi^-/\pi^+$37 and $\pi^-/\pi^+$38, producing the local enhancement in $\pi^-/\pi^+$39 at $\pi^-/\pi^+$40 MeV. The feature appears only when the in-medium pion potential is included; without it, the bump is not emphasized in the displayed spectra. By contrast, the stiffness of the nuclear symmetry energy has negligible influence on $\pi^-/\pi^+$41 around the bump region, even though neutron/proton squeeze-out ratios remain sensitive to $\pi^-/\pi^+$42 at higher momenta (Feng, 2016).

5. Broad bumps in isoscalar single- and double-pion production

In hadronic spectroscopy, “pion bump” can denote a broad enhancement in the isoscalar part of single-pion production. The isoscalar cross section is extracted through

$\pi^-/\pi^+$43

After consolidation of available data and small $\pi^-/\pi^+$44–$\pi^-/\pi^+$45 renormalizations within quoted systematics, the energy dependence is described by a broad Lorentzian-like structure with peak position $\pi^-/\pi^+$46–$\pi^-/\pi^+$47 and width $\pi^-/\pi^+$48 (Clement et al., 2020). The bump peaks about $\pi^-/\pi^+$49 below the nominal $\pi^-/\pi^+$50 threshold at $\pi^-/\pi^+$51, and the isoscalar $\pi^-/\pi^+$52 invariant mass peaks at $\pi^-/\pi^+$53 with apparent width $\pi^-/\pi^+$54, a pattern interpreted as consistent with a bound or quasi-bound $\pi^-/\pi^+$55 configuration (Clement et al., 2020).

The same work argues that the observed shape is incompatible with a pure $\pi^-/\pi^+$56-channel opening of Roper production, which would rise with increasing phase space, and also with the narrow $\pi^-/\pi^+$57 Breit–Wigner proposed elsewhere (Clement et al., 2020). Instead, it proposes molecular-like $\pi^-/\pi^+$58–$\pi^-/\pi^+$59 dibaryon states with $\pi^-/\pi^+$60 and $\pi^-/\pi^+$61, overlapping near threshold. The final-state partial-wave content is central: the isoscalar $\pi^-/\pi^+$62 spectrum peaks at low $\pi^-/\pi^+$63, consistent with exit waves $\pi^-/\pi^+$64 and $\pi^-/\pi^+$65, not $\pi^-/\pi^+$66 (Clement et al., 2020).

A related but distinct usage appears in double-pionic fusion. A sequential single-pion production chain,

$\pi^-/\pi^+$67

was proposed as an explanation of the $\pi^-/\pi^+$68 peak. A corrected treatment instead yields a broad enhancement near $\pi^-/\pi^+$69–$\pi^-/\pi^+$70 with width $\pi^-/\pi^+$71, not a narrow $\pi^-/\pi^+$72-like structure (Bashkanov et al., 2023). In the $\pi^-/\pi^+$73 channel the residual bump after subtraction of the $\pi^-/\pi^+$74 contribution has peak cross section $\pi^-/\pi^+$75; the authors identify it with the sequential mechanism and/or broad isoscalar dibaryonic excitations rather than with the genuine narrow resonance at $\pi^-/\pi^+$76, $\pi^-/\pi^+$77, $\pi^-/\pi^+$78 (Bashkanov et al., 2023). The channel relation

$\pi^-/\pi^+$79

makes the neutral channel particularly clean for isolating the isoscalar bump (Bashkanov et al., 2023).

6. The $\pi^-/\pi^+$80-region pion bump in lepton–nucleus scattering

In inclusive lepton–nucleus scattering, the pion bump denotes the enhancement just above the quasielastic peak in distributions such as $\pi^-/\pi^+$81 or $\pi^-/\pi^+$82, centered near the $\pi^-/\pi^+$83 resonance. For a nucleon initially at rest,

$\pi^-/\pi^+$84

so the bump region corresponds to $\pi^-/\pi^+$85, with nuclear motion and removal energy smearing this mapping (Isaacson et al., 26 Aug 2025). In the Achilles event generator, the feature arises from single-$\pi^-/\pi^+$86 production via $\pi^-/\pi^+$87 excitation and nearby $\pi^-/\pi^+$88 resonances, followed by pion propagation and final-state interactions in the nucleus.

The electroweak vertex is modeled with the ANL–Osaka Dynamical Coupled-Channels framework, in which the hadronic current is the coherent sum of nonresonant background and resonant terms. The exclusive $\pi^-/\pi^+$89 hadron tensor is folded with realistic hole spectral functions $\pi^-/\pi^+$90, so shell structure, correlations, and removal energy broaden the bump already at the production level (Isaacson et al., 26 Aug 2025). Final-state interactions are treated by a semi-classical intranuclear cascade with DCC meson–baryon amplitudes; pion absorption is handled either by an Oset–Salcedo optical-potential mode,

$\pi^-/\pi^+$91

or by explicit propagation of intermediate resonances such as the $\pi^-/\pi^+$92 (Isaacson et al., 26 Aug 2025).

Initial-state smearing, elastic and inelastic rescattering, charge exchange, and absorption reshape the free-nucleon $\pi^-/\pi^+$93 peak into the nuclear pion bump. The net effect is to widen the $\pi^-/\pi^+$94-distribution, shift some strength to lower $\pi^-/\pi^+$95, and suppress low-$\pi^-/\pi^+$96 yields (Isaacson et al., 26 Aug 2025). Achilles reproduces the qualitative structure of inclusive electron scattering on $\pi^-/\pi^+$97 and $\pi^-/\pi^+$98, namely the QE peak followed by the $\pi^-/\pi^+$99-region component, and it gives favorable comparisons to T2K, MINER$\sqrt{s}\approx 2.31$00A, and MicroBooNE semi-inclusive data. The distinction between the Virtual Resonances and Propagating Resonances cascade modes brackets the amount of absorption and migration between $\sqrt{s}\approx 2.31$01 and $\sqrt{s}\approx 2.31$02 samples (Isaacson et al., 26 Aug 2025).

7. Transient TeV $\sqrt{s}\approx 2.31$03 bumps, ambiguities, and future measurements

A distinct high-energy usage of the term appears in flaring blazars. Here a $\sqrt{s}\approx 2.31$04 bump is a narrow, quasi-line-like excess in the very-high-energy $\sqrt{s}\approx 2.31$05-ray spectrum, produced when protons of tens of TeV interact with hard X-ray photons through the $\sqrt{s}\approx 2.31$06-resonance channel $\sqrt{s}\approx 2.31$07, followed by $\sqrt{s}\approx 2.31$08 (Petropoulou et al., 2023). For Mrk 501 during the extreme X-ray flare on MJD $\sqrt{s}\approx 2.31$09, MAGIC data showed a Gaussian-like excess at $\sqrt{s}\approx 2.31$10 at the $\sqrt{s}\approx 2.31$11–$\sqrt{s}\approx 2.31$12 level. The threshold relation

$\sqrt{s}\approx 2.31$13

shows why such features require hard X-ray target photons and are favored during synchrotron-dominated flares with peak energies $\sqrt{s}\approx 2.31$14 (Petropoulou et al., 2023). CTA simulations indicate that a $\sqrt{s}\approx 2.31$15 bump of this type could be detected at $\sqrt{s}\approx 2.31$16 with a $\sqrt{s}\approx 2.31$17-minute exposure, but the same modeling also requires an optically thin $\sqrt{s}\approx 2.31$18 region whose energy content is dominated by relativistic protons and whose jet power is highly super-Eddington (Petropoulou et al., 2023).

Across all domains, the literature does not treat the pion bump as a complete diagnostic on its own. In W44, current systematics below $\sqrt{s}\approx 2.31$19 MeV still allow electron bremsstrahlung to mimic the low-energy break, which is why future MeV detectors are treated as decisive (Liu et al., 2024). In Galactic neutrino searches, a MeV $\sqrt{s}\approx 2.31$20-decay signature establishes hadronic $\sqrt{s}\approx 2.31$21-ray production at low energies but does not itself guarantee IceCube-detectable neutrino emission, since the proton spectrum may be steep or cut off below the multi-TeV range (Granados et al., 30 Jul 2025). In heavy-ion collisions, the $\sqrt{s}\approx 2.31$22 bump diagnoses in-medium pion optics more directly than the stiffness of the symmetry energy (Feng, 2016).

Future work therefore separates into domain-specific programs. In MeV astrophysics, MeGaT-, COSI-, and AMEGO-class missions are motivated by the need to resolve the $\sqrt{s}\approx 2.31$23–$\sqrt{s}\approx 2.31$24 MeV turnover and the bremsstrahlung floor below it (Liu et al., 2024). In multimessenger Galactic studies, IceCube-Gen2 and KM3NeT are expected to improve sensitivity to pion-bump-selected source populations (Granados et al., 30 Jul 2025). In neutrino-event simulation, extensions such as coherent single-$\sqrt{s}\approx 2.31$25 production, $\sqrt{s}\approx 2.31$26–$\sqrt{s}\approx 2.31$27 mechanisms, and medium-modified resonances are identified as the next steps for standardizing the $\sqrt{s}\approx 2.31$28-region pion bump across targets and channels (Isaacson et al., 26 Aug 2025). The persistence of the term across these fields reflects a common link to pion physics, but the observable itself remains irreducibly context specific.