Self-confinement of ultra-high-energy nuclei in cosmic filaments: implications for the UHECR spectrum and composition

arXiv:2605.14869 · astro-ph.HE · Submitted 2026-05-14 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: We're looking at the paper: "Self-confinement of ultra-high-energy nuclei in cosmic filaments".

Jocelyn: Here's the gist of the paper. The gist: Self-confinement scenarios, driven by self-generated magnetic turbulence in cosmic filaments, provide a physically motivated mechanism to transform standard acceleration spectra into the hard effective spectra inferred from ultra-high-energy cosmic ray observations.

Vera: Model Setup and Instability Growth The scenario extends the self-confinement framework to a mixed nuclear composition and considers a steady UHECR source embedded in the intergalactic medium (IGM) The electric current carried by escaping UHECRs amplifies magnetic turbulence in the surrounding intergalactic medium, reduces the diffusion coefficient, and delays particle escape in a rigidity-dependent way The growth of the instability is controlled by this current The transverse size of the escaping flux is taken to be the larger of the source size RS and the particle Larmor radius in the background field RL(R) At saturation, the fastest-growing mode reaches a scale comparable to the Larmor radius of the particles carrying the current in the amplified field: kmax(R) = 1/RL(R) = eδB R.

Jocelyn: Confinement and Transport Physics The transport within the confinement region is modeled using a Leaky Box model, where the equilibrium density of each nuclear species is determined by the balance between injection, escape, and interactions The escape time includes both diffusive transport through the self-generated turbulence and advection with the background plasma The equilibrium density of confined particles is given by n(A,Z)I(E,z) = q(A,Z)acc(E) τ eff(E,z), where the effective confinement time is τ eff(E,z) = τ esc(E,z) τ int(E,z).

Vera: The low-rigidity suppression arises because particles cannot travel across one magnetic-field coherence length within the available source lifetime.

Paper summary: Vera: So let's talk a bit more about how they set up this self-confinement model, looking at the physics behind that initial instability growth. They start by saying that when UHECRs escape, they generate magnetic turbulence in the surrounding intergalactic medium. This turbulence actually reduces the diffusion coefficient around the source and delays particle escape in a rigidity-dependent way.

Jocelyn: That delay is governed by an electric current carried by those escaping UHECRs, which amplifies the magnetic field fluctuations in a specific way. They take the transverse size of that escaping flux to be the larger of the source size RS and the particle Larmor radius RL(R) in the background field.

Subrahmanyan: The instability itself is controlled by this current, and they determine its growth rate by looking at k max equal one over RL(R), which simplifies to e delta B R. This tells us exactly what scale the turbulence reaches when it saturates, based on the local magnetic field strength.

Vera: And once that saturation point is found, they use a leaky box model to describe the transport within that confinement region. They determine the equilibrium density of each nuclear species by balancing injection from acceleration against escape and interactions with ambient photon fields.

Jocelyn: The effective confinement time is defined as tau eff equal tau esc times tau int, where the escape time includes both transport through the self-generated turbulence and advection with the background plasma. This is how they figure out how long a particle actually stays trapped in that environment.

Subrahmanyan: They then show that the low-rigidity suppression happens because particles simply cannot travel across one magnetic-field coherence length within the available lifetime of the source. This is captured by an exponential attenuation factor inspired by the Green function of a one-dimensional advection diffusion equation.

Vera: That exponential factor is what translates those physical transport rules into a spectral shape, showing that this trapping suppresses the escaping flux below a characteristic rigidity of about one EV <ref:2605.14869#pg1>. The paper notes that this is important because it allows them to match the observed suppression without needing an extremely hard injection spectrum from the source itself.

Jocelyn: It’s really about how the magnetic field structure—the coherence length lambda B—plays a critical role in setting that cutoff rigidity Rc, and they show different coherence lengths like fifteen Mpc, thirty Mpc, and fifty Mpc lead to different spectral shifts.

Subrahmanyan: The math shows that for lambda B equal to thirty or fifty Mpc, the model achieves a reasonable agreement with observed data while having a negligible contribution from primary protons in the escaping flux. This suggests that in those filament environments, this confinement is quite effective at shaping the observed spectrum.

Vera: But they also mention that if you look at different coherence lengths, like fifteen Mpc or thirty Mpc, the escaping emissivity shifts toward higher energies accessible at the source and can show an apparent hardening just below that cutoff energy. It’s a bit of a trade-off there.

Jocelyn: So to summarize this part, they've shown how the instability grows based on current currents, how transport is modeled using leaky boxes with effective confinement times, and what the spectral consequences are for different magnetic field scales.

Conclusion: Vera: Now we move into what these changes actually mean for the spectrum and composition of those cosmic rays. They look at how the cutoff rigidity Rc depends on all those model parameters, which highlights some degeneracies in this approach, but they find that for different coherence lengths, like fifteen Mpc or thirty Mpc or fifty Mpc, the escaping emissivity is shifted toward higher energies accessible at the source.

Jocelyn: This shifting means that for a given source luminosity and magnetic field setup, you can get an apparent hardening just below the cutoff energy without needing an exceptionally hard injection spectrum from the accelerator. It’s a way to explain what we observe without assuming something weird about the source physics.

Subrahmanyan: From a cosmological perspective, this is significant because it links local magnetic field properties and source environments—like cosmic filaments—to the observed high-energy phenomena, suggesting that structure formation plays a role in shaping these spectra.

Vera: The model explores how heavy nuclei interact with the ambient radiation fields during this confinement phase. They find that heavy nuclei efficiently photodisintegrate, which produces secondary protons that contribute to the observed flux below the ankle energy.

Jocelyn: That’s interesting because it directly addresses the composition puzzle. It suggests that even if we see a mix of nuclei at high energies, there’s a secondary component of protons coming from these interactions during confinement that helps account for the composition we measure at lower energies.

Subrahmanyan: This secondary proton contribution is important because it modifies the overall mass composition observed and helps bridge the gap between what we expect from pure acceleration and what we actually detect.

Vera: They also compute the associated cosmogenic neutrino and gamma-ray emission resulting from these confined particles interacting with the background radiation fields. The predicted neutrino flux remains compatible with current limits, which is a good sign for this scenario.

Jocelyn: However, they also note that the diffuse gamma-ray background provides a potentially strong constraint on the most extreme configurations of this confinement model. If you push the parameters to their limit, those gamma rays get pretty strong and put pressure on the theory.

Subrahmanyan: The paper’s main conclusion here is that self-generated turbulence can suppress that escaping flux below the EV scale for source luminosities and magnetic-field coherence lengths compatible with UHECR sources hosted in galaxy clusters and propagating through cosmic filaments.

Vera: So, essentially, this mechanism works when you have these specific conditions—the right source power and the right magnetic field coherence length—to naturally suppress the flux below one EV <ref:2605.14869#pg1>.

Gran Sasso Science Institute (GSSI) · INFN-Laboratori Nazionali del Gran Sasso (LNGS)

astro-ph.HE

Submitted: 2026-05-14

Updated: 2026-05-14

Comments: 13 pages (11 main text + 2 appendix), 7 figures

DOI: 10.1051/0004-6361/202660944

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 72/100

The gist: The gist: Self-confinement scenarios, driven by self-generated magnetic turbulence in cosmic filaments, provide a physically motivated mechanism to transform standard acceleration spectra into the

Key concepts

Self-confinement: This is a scenario where the UHECRs themselves generate magnetic turbulence in their surrounding environment. This turbulence acts like a cage, trapping the particles near their sources and delaying their escape. This process modifies the standard acceleration spectrum, making it appear harder at lower energies.
Leaky Box Model: This is a transport model used to describe how particles move within the confinement region. It balances particle injection from the source against particle escape through turbulence and interactions like photodisintegration. The equilibrium density of each nuclear species is found by balancing these rates.
Photodisintegration: When heavy nuclei are confined, they interact with ambient radiation fields, such as photons in the cosmic background. This interaction causes the nuclei to break apart into lighter particles, specifically secondary protons. These secondary protons help explain the observed composition of UHECRs below the ankle.
Magnetic Field Coherence Length ($\lambda_B$): This refers to the typical distance over which magnetic field lines are correlated or coherent in the intergalactic medium. The paper finds that confinement is most effective when this length is in the range of 5 to 50 Mpc, suggesting that filaments with these specific magnetic properties are ideal for this mechanism.

Terminology

Coherence — 1× in this paper · explained in 2 episodes
Coherence Length — 1× in this paper · explained in 2 episodes
Streaming Instability — 1× in this paper · explained in 3 episodes

Summary

The gist: Self-confinement scenarios, driven by self-generated magnetic turbulence in cosmic filaments, provide a physically motivated mechanism to transform standard acceleration spectra into the hard effective spectra inferred from ultra-high-energy cosmic ray observations.

Model Setup and Instability Growth

The scenario extends the self-confinement framework to a mixed nuclear composition and considers a steady UHECR source embedded in the intergalactic medium (IGM) The electric current carried by escaping UHECRs amplifies magnetic turbulence in the surrounding intergalactic medium, reduces the diffusion coefficient, and delays particle escape in a rigidity-dependent way The growth of the instability is controlled by this current The transverse size of the escaping flux is taken to be the larger of the source size RS and the particle Larmor radius in the background field RL(R) At saturation, the fastest-growing mode reaches a scale comparable to the Larmor radius of the particles carrying the current in the amplified field: kmax(R) = 1/RL(R) = eδB R.

Confinement and Transport Physics

The transport within the confinement region is modeled using a Leaky Box model, where the equilibrium density of each nuclear species is determined by the balance between injection, escape, and interactions The escape time includes both diffusive transport through the self-generated turbulence and advection with the background plasma The equilibrium density of confined particles is given by n(A,Z)I(E,z) = q(A,Z)acc(E) τ eff(E,z), where the effective confinement time is τ eff(E,z) = τ esc(E,z) τ int(E,z). The low-rigidity suppression arises because particles cannot travel across one magnetic-field coherence length within the available source lifetime. This is captured by an exponential attenuation factor inspired by the Green function of the one-dimensional advection–diffusion equation.

Spectral Consequences and Composition

The model explores how the cutoff rigidity Rc depends on model parameters, highlighting degeneracies in this approach. For different coherence lengths, such as λB = 15 Mpc, 30 Mpc, and 50 Mpc, the escaping emissivity is shifted toward the highest energies accessible at the source and can exhibit an apparent hardening just below the cutoff energy. A suppression of R ∼ 1 EV can reproduce combined fits to UHECR spectrum and composition without requiring extremely hard injection spectra. For λB = 30–50 Mpc, the model achieves a reasonable agreement with observed data, with a negligible contribution from primary protons.

Multi-messenger Signatures

Particles confined within magnetized regions can continue to interact with the ambient background radiation fields, producing cosmogenic neutrinos and gamma-rays. The production of secondary protons from photodisintegration during confinement is expected to enhance the composition observed below the ankle. The predicted neutrino flux remains compatible with current limits, while the diffuse gamma-ray background provides a potentially strong constraint on the most extreme confinement configurations.

Astrophysical Viability and Constraints

The scenario is viable for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The required source density is of the order 10−6–10−4 Mpc−3, which is consistent with arrival-direction analyses. The main limitation set by the pre-existing intergalactic magnetic field is that particles remain sufficiently tied to magnetic-field lines over the full energy range of interest. The mechanism appears most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The resulting allowed parameter space is shown in the left panel of Fig. 7.

Future Directions

Future measurements of the mass composition at the highest energies, in particular with the AugerPrime upgrade of the Pierre Auger Observatory, together with improved sensitivity to diffuse UHE neutrinos expected from next-generation observatories such as GRAND, GCOS, and IceCube-Gen2 will be crucial to further probe the viable parameter space of this self-confinement scenario. The mechanism appears to be most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The mechanism appears to be most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The mechanism appears to be most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The mechanism appears to be most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The mechanism appears to be most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The mechanism appears to be most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The mechanism appears to be most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The mechanism appears to be most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The mechanism appears to be most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The mechanism appears to be most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc. The mechanism appears to be most efficient for sources with a total cosmic-ray power of the order of 1043–1046 erg s−1, averaged over the age of the Universe, and for magnetic-field coherence lengths in the range λB ∼ 5–50 Mpc.

Improvements for AI systems

Precision modeling of UHECR source injection spectra can be improved by incorporating self-confinement effects, which transform standard acceleration spectra into the hard effective spectra inferred from observations. This allows AI systems to distinguish between intrinsic accelerator properties and observable transport imprints on the spectrum.
Composition inference algorithms can be enhanced to account for secondary particle production within the confinement region, as heavy nuclei efficiently photodisintegrate, producing secondary protons that naturally contribute below the ankle and help account for the observed composition evolution. This enables more accurate mass composition reconstruction from shower maximum depth measurements (Xmax).
Multi-messenger prediction capabilities can be expanded to include cosmogenic neutrinos and gamma rays generated by confined particles, as the predicted neutrino flux remains compatible with current limits, while the diffuse gamma-ray background provides a potentially strong constraint on the most extreme confinement configurations. This allows AI to test physical scenarios using neutrino and gamma-ray data alongside cosmic ray observations.
Parameter space exploration tools can be refined to map the relationship between source luminosity and magnetic-field coherence length by incorporating constraints derived from transport regimes, such as the blue curves show the coherence length λB(L) required, for a given luminosity, to suppress the escaping flux at 500 PV (dashed), 1 EV (solid), and 2 EV (dotted). This helps identify astrophysically viable regions for self-confinement.
Source identification algorithms can be improved to prioritize candidate sources based on their inferred cosmic ray power, focusing on those that satisfy the source luminosity must exceed L min = cΛR2sB204 to trigger the non-resonant streaming instability (NRSI). This steers searches toward environments like galaxy clusters and filaments.
Spectral shape analysis can be refined to account for the spectral filtering effect of confinement, where the confinement effect acts as a spectral filter, reshaping the spectrum emerging from the source environment even when the injected spectrum itself remains unchanged. This prevents misinterpretation of source acceleration mechanisms based solely on observed escape spectra.

Abstract

The spectrum and composition of ultra-high-energy cosmic rays (UHECRs) suggest that the population dominating above the ankle releases particles with an unusual hard spectrum at low rigidity, below the EV scale. In self-confinement scenarios, such an apparent hardening arises from transport: escaping UHECRs generate magnetic turbulence that delays their own release from the magnetized environments surrounding their sources. We extend the self-confinement scenario based on the non-resonant streaming instability to a mixed nuclear composition. We describe the confinement region with an effective leaky-box model including escape, photodisintegration, and secondary production. We then compare the resulting spectrum and composition with Auger measurements and compute the associated cosmogenic neutrino and gamma-ray emission. We find that self-generated turbulence can suppress the escaping flux below the EV scale for source luminosities and magnetic-field coherence lengths compatible with UHECR sources hosted in galaxy clusters and propagating through cosmic filaments. During confinement, heavy nuclei efficiently photodisintegrate, producing secondary protons that contribute below the ankle and help account for the observed composition. The predicted neutrino flux remains compatible with current limits, while the diffuse gamma-ray background provides a potentially strong constraint on the most extreme configurations.