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

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>.