Microquasar Remnants as Pevatrons Illuminating the Galactic Cosmic Ray Knee

arXiv:2602.08940 · astro-ph.HE · Submitted 2026-02-09 · Read on arXiv

Listen

Transcript

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

Vera: Next we'll be talking about the paper "Microquasar Remnants as Pevatrons Illuminating the Galactic Cosmic Ray Knee".

Jocelyn: The paper was written by Bing Theodore Zhang and Shiqi Yu from Key Laboratory of Particle Astrophysics and Experimental Physics Division and Computing Center, Institute of High Energy Physics, Chinese Academy of Sciences and TIANFU Cosmic Ray Research Center and University of Utah.

Vera: Stay tuned as we take you through the paper and discuss its implications.

Paper discussion segment 1 — Title and Authors: Vera: We've seen how microquasars are candidates for these high-energy PeVatrons, but the paper by Zhang and Yu provides a much more detailed physical picture of how this actually works within the Galactic Cosmic Ray Knee. They’re not just saying they *can* produce the flux; they are describing *how* it is flowing.

Jocelyn: The authors are showing us that cosmic ray transport isn't random at all, which is a huge observational confirmation of something we've long suspected. They are modeling anisotropic diffusion through the Galactic magnetic field, and that’s a key physical insight.

Subrahmanyian: The magnetic field acts like a highway system for these charged particles. This paper demonstrates that the geometry of the GMF isn't just an accessory; it determines how much flux we actually see at Earth.

Vera: It does, and this is where they make some very specific claims about connectivity. They are comparing sources that are physically close but magnetically disconnected with sources that might be slightly further away but perfectly aligned with the local field lines.

Jocelyn: That’s a powerful distinction, because the paper is showing us that distance isn' not the only thing governing our local flux—the magnetic architecture is just as important. It’s a really important lesson for how we interpret data from surveys.

Subrahmanyian: The implication here is that our observational campaigns need to be aware of this complex connectivity, recognizing that spatial proximity alone won't always tell the whole story of cosmic ray origins.

Vera: That certainly changes the way we approach source identification, knowing that the magnetic path dictates what we observe. Let’s see how they quantify this in their summary and what does it mean for our observations next segment.

Paper discussion segment 2 — Summary: Jocelyn: The authors provide a very clear summary of their findings, which really confirms that the Galactic magnetic field is the primary driver of the local flux we measure. They used simulations to show how this anisotropy affects particle propagation dramatically.

Subrahmanyian: It’s a stark contrast between what happens with V616 Mon and what happens with Cygnus X-one for example, which really drives home the point about magnetic channeling being crucial. The physics of the transport is dominating the simple geometry.

Vera: And they use these two specific sources to illustrate that distance alone doesn't matter. V616 Mon is closer than Cygnus X-one but its lack of connectivity means its flux is heavily suppressed, which would be a huge shock to any models based on proximity.

Jocelyn: The way the GMF channels the rays toward us, it creates this specific filamentary morphology in the CR clouds that we can now use to interpret our own observational data. It gives us a physical map of how these particles arrive at Earth.

Subrahmanyian: The key insight is that they show how this anisotropic transport fundamentally dictates the observed flux, making it a highly efficient mechanism for local particle delivery.

Vera: It's not just one source contributing either, but the cumulative effect of how they all align with our specific patch of the magnetic field. It sets up a great example of how complex these local environments are.

Jocelyn: This finding is really going to change how we interpret the "local" nature of the cosmic ray signal, moving beyond simple source strength assumptions. Let's look at how they handle the time evolution in this paper next segment.

Paper discussion segment 3 — Improvements: Subrahmanyian: The authors suggest a significant improvement to our understanding by focusing on the cumulative contribution of microquasar remnants, not just the active ones, which is a huge shift in perspective for cosmic ray modeling.

Vera: This is where it really gets interesting. They are showing us that while active sources might be intense locally, their continuous presence doesn' not match what we see at the PeV knee. The history of the remnants must play a major role.

Jocelyn: It turns out that the expanded clouds of CRs from these remnants have had sufficient time to diffuse through the GMF and fill that observed "sea" of diffuse PeV CR flux over a much larger area. That's a massive change in our interpretation of the sky.

Subrahmanyian: The simulation shows that for this cumulative contribution to work, we need a specific source age of about one million years, favoring the remnant phase over younger, active systems. This is a strong constraint on our models.

Vera: And they also found that using a hard injection spectrum—Model A—is much more efficient at explaining the observed flux than standard power-law assumptions, which is a huge technical victory for our modeling efforts.

Jocelyn: It’s not just about the age of these sources; it' about their energy injection profile, showing that we need to look for specific spectral shapes to match the PeV bump.

Subrahmanyian: The implications are vast—we're seeing a signature of the integrated history of the galaxy, where those ancient jets are still contributing long after their central engines have shut off.

Vera: It's definitely a nuanced picture that we aren't just seeing in one spot, but across the entire galactic population. Let’s wrap up this discussion and talk about what it all means for our future observations.

Conclusion: Jocelyn: We’ve covered so much ground today, from the magnetic highways of our galaxy to the crucial role of remnants in explaining the cosmic ray knee, thanks to this paper. It's really cemented microquasars as PeVatrons, but not in a simplistic way.

Vera: The authors have shown that this whole picture is driven by a combination of factors—the injection spectrum, the age, and the magnetic field—and their findings are extremely robust. We've learned that relying on simple proximity will not work for our observations.

Subrahmanyian: I think the most exciting thing is how it frames the local environment as a dynamic snapshot in time, where these remnants provide a lasting signature of galactic history. It really moves us beyond steady-state assumptions about the cosmos.

Jocelyn: It’s a huge shift in perspective, realizing that for future studies, we need to fully account for both the spatial distribution and the particle escape processes of these microquasars.

Vera: We're so excited about this work by Zhang and Yu, and we hope that our discussion has helped listeners understand the profound implications of Microquasar Remnants as Pevatrons Illuminating the Galactic Cosmic Ray Knee.

Subrahmanyian: I agree, it’s a clear path forward for understanding that complex cosmic ray landscape.

Jocelyn: It's certainly a story of accumulated history, Vera.

Vera: Absolutely. And that’ all for us today, everyone; we're looking forward to the next paper we get to discuss with you all soon.

astro-ph.HE

Submitted: 2026-02-09

Updated: 2026-09-01

Comments: 4 + 4 figures and 2 table

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 82/100

The gist: This paper investigates the role of microquasars as primary candidates for Galactic PeVatrons and their collective contribution to the cosmic ray (CR) "knee." By simulating anisotropic diffusive

Key concepts

Galactic Magnetic Field (GMF): The GMF acts like a highway system for charged particles. This field dictates how cosmic ray flux reaches Earth, meaning the magnetic architecture is as important as spatial proximity when interpreting observations.
Anisotropic Transport: This describes particle movement that is not random. The GMF channels rays in specific directions, leading to filamentary patterns in cosmic ray clouds and fundamentally dictating the observed flux at Earth.
Microquasar Remnants: These are the remnants of microquasars, not just their active phases. Simulations suggest these remnants have diffused through the GMF over time (around one million years), contributing to the overall diffuse PeV flux.

Terminology

Anisotropy — 6× in this paper · explained in 3 episodes
Alignment — 3× in this paper · explained in 1 episode
Anisotropic Transport — 1× in this paper · explained in 1 episode

Summary

This paper investigates the role of microquasars as primary candidates for Galactic PeVatrons and their collective contribution to the cosmic ray (CR) knee. By simulating anisotropic diffusive propagation through the Galactic magnetic field (GMF), the authors explore how the interplay between source age, magnetic connectivity, and injection spectra shapes the observed spectral bump in the proton spectrum.

Methodology and Modeling

The researchers simulate CR propagation using the DiffusiveSDE module within the CRPropa 3.2 framework, representing the GMF via the JF12 model. The expected CR flux at Earth is estimated by recording particles entering a spherical volume of radius 200 pc centered on the solar neighborhood. To model the injection, the authors adopt a physically motivated Model A, which is characterized by a hard spectral shape with a steep cutoff, and compare it against a generic power-law template with a spectral index of 2.

The simulation incorporates several critical parameters:

A time-dependent injection rate Q(E) integrated over the source's history.
Anisotropic diffusion where the perpendicular component is D = epsilon D.
A cylindrical Galactic halo with R G = 20 kpc and H G = 2 kpc as boundary conditions.
The assumption that CRs are injected into the interstellar medium immediately after escaping from the sources.

The Impact of Magnetic Connectivity

A central finding is that the GMF channels CRs preferentially along field lines, creating a transport regime where magnetic connectivity between sources and the solar neighborhood determines the local flux. This implies that two equidistant sources can yield vastly different observed fluxes depending on their alignment with the GMF. For instance, active sources aligned with local GMF lines, such as Cygnus X-1, exhibit significant flux enhancements, whereas magnetically disconnected sources like V616 Mon are strongly suppressed due to inefficient cross-field transport.

The study demonstrates that:

The impact of the anisotropy diffusion parameter epsilon is more significant than the physical distance of the source.
Sources in the Cygnus region are magnetically connected to the solar neighborhood, funneling CRs toward Earth.
The spatial distribution of 10 PeV CRs exhibits a filamentary morphology driven by the magnetic field.

Microquasar Remnants and the CR Knee

The paper argues that the observed proton bump at the CR knee is best reproduced by the cumulative contribution of microquasar remnants rather than active sources alone. These are defined as fossil systems where relativistic jet activity has ceased, leaving behind expanding clouds of CRs. While active sources could potentially create localized spectral features, the absence of such features above 10 PeV suggests that the current local CR environment is shaped by these remnants, which represent an integrated history of Galactic microquasar activity.

The contribution of these remnants is characterized by:

The ability of remnant clouds to soften and fill the observed diffuse PeV CR 'sea'.
A preference for a harder injection spectrum (Model A), which requires a CR loading efficiency of only eta cr about 0.5% compared to 3% for standard power-law models.
The realization that the local PeV CR environment is a unique snapshot of the stochastic history and proximity of nearby microquasar activities.

Improvements for AI systems

1. Physics-Informed Neural Networks (PINNs) for Anisotropic Diffusion Modeling

Improvement: Integrate the anisotropic diffusion tensor (D and D) and the JF12 Galactic Magnetic Field (GMF) model directly into the loss function of a PINN. This moves beyond isotropic assumptions by encoding the relationship where transverse transport is driven by field-line random walks.
Capability: The AI can rapidly simulate cosmic ray (CR) density distributions in complex, non-uniform magnetic topologies, bypassing the massive computational cost of traditional Stochastic Differential Equation (SDE) solvers like CRPropa.

2. Graph Neural Networks (GNNs) for Magnetic Connectivity Mapping

Improvement: Represent the Galaxy as a dynamic graph where nodes are astrophysical sources (microquasars, remnants) and observers (Solar neighborhood), and edges represent magnetic field lines weighted by the anisotropy parameter epsilon.
Capability: The system can perform real-time connectivity analysis to predict the flux contribution of any given source. It can instantly determine if a source like Cygnus X-1 will provide a flux enhancement or if a source like V616 Mon will be suppressed due to magnetic disconnection, without running full propagation simulations.

3. Stochastic Generative Models (Diffusion/GANs) for Snapshot Temporal Reconstruction

Improvement: Train generative models on the stochastic birth rates (R MQ) and active lifetimes (tau dur) of microquasars to model the integrated history of the Galaxy.
Capability: Instead of predicting a steady-state average, the AI can generate high-fidelity probabilistic snapshots of the local CR environment. This allows the system to quantify the uncertainty and variance inherent in a dynamic background, simulating how a few recent or nearby events (remnants) dominate the observed knee spectral feature.

4. Multi-Scale Transformer Architectures for Inverse Spectral Problem-Solving

Improvement: Develop a Transformer-based architecture capable of ingesting multi-modal inputs: 3D spatial coordinates (Galactic l, b), energy spectra (e.g., LHAASO/IceTop data), and temporal injection histories.
Capability: The system can perform sophisticated inverse modeling. Given an observed spectral bump or knee, the AI can backtrack to identify the most likely spatial distribution, age, and injection spectrum (e.g., distinguishing between Model A hard spectra and standard E-2 power laws) of the contributing microquasar remnants.