Super PeVatrons: Extreme Galactic Accelerators
- Super PeVatrons are extreme Galactic accelerators that exceed the typical ∼1 PeV energy threshold, often evidenced by gamma-ray emissions above 100 TeV.
- They are identified through hard, nearly curvature-free gamma-ray spectra and de-attenuation analyses that indicate hadronic acceleration near or beyond the cosmic-ray knee.
- Observations from sources like the Cygnus X bubble and advances with instruments such as CTAO and LHAASO are driving the frontier in distinguishing these ultra-extreme accelerators.
“Super PeVatrons” are best understood as the most extreme subset of Galactic PeVatrons: sources whose observed or inferred particle acceleration extends beyond the ordinary PeV scale and enters a regime relevant to the knee and, in some usages, well beyond it. The term is not fully standardized across the literature. In one usage, it refers to Galactic systems whose gamma-ray emission implies parent protons reaching at least PeV, as argued for the Cygnus X bubble; in another, it denotes sources extending above $1$ PV and up to PV in rigidity; and in a propagation-oriented usage it marks sources whose intrinsic spectra and luminosities may be substantially more extreme than their attenuated observed spectra suggest (Collaboration, 2023, Wang et al., 28 Jul 2025, Zhang et al., 2024). Across these usages, the common idea is an accelerator that is not merely capable of reaching the PeV domain, but that probes the upper end of Galactic particle acceleration, where source physics, transport, and radiative diagnostics all become nontrivial (Wilhelmi et al., 2024).
1. Terminology and conceptual scope
Cristofari’s review emphasizes that a Galactic PeVatron is, in the strict sense, an astrophysical source capable of accelerating particles to energies of order , but that the field often conflates “finding PeVatrons” with “solving the origin of Galactic cosmic rays” (Cristofari, 2021). That distinction remains central for super-PeVatrons. A source can be a PeVatron without being the dominant source class of Galactic cosmic rays; the Crab Nebula is the canonical example, being a convincing electron PeVatron but not an accepted source of the bulk of Galactic CR protons (Cristofari, 2021).
In the strongest hadronic usage, super-PeVatrons are sources whose gamma-ray output implies parent hadrons extending well beyond the ordinary PeVatron threshold. The clearest explicit example is the LHAASO study of the Cygnus X bubble, which interprets photons up to a few PeV, including events above $1$ PeV, as evidence for “Super PeVatron(s)” accelerating protons to at least PeV (Collaboration, 2023). A broader rigidity-based definition appears in the proposal that Galactic super-accreting X-ray binaries may supply sources extending from above $1$ PV to PV, motivated by the possibility that the Galactic CR component extends well above the knee (Wang et al., 28 Jul 2025).
A more cautious, source-classification-oriented usage appears in attenuation studies of LHAASO sources. There the phrase “super PeVatrons” refers to sources whose intrinsic spectra, after correcting for propagation losses on the ISRF and CMB, may imply harder spectra, higher luminosities, and higher parent-particle energies than the observed spectra alone suggest (Zhang et al., 2024). This usage does not require an already established PeV hadronic accelerator; it identifies a bias in the observational ranking of extreme Galactic sources.
2. Physical criteria and observational diagnostics
The standard observational logic begins from the approximate hadronic relation 0, which makes gamma rays near 1 TeV natural tracers of parent protons in the PeV domain (Cristofari, 2021). The literature repeatedly uses this as a practical heuristic, but also warns that it is not an absolute discriminator because inverse-Compton emission from electrons is suppressed above 2 TeV by the Klein–Nishina effect without being eliminated (Cristofari, 2021). For that reason, neutrinos remain the cleanest discriminator of hadronic acceleration.
The theoretical confinement and acceleration requirements are commonly summarized by a Hillas-type estimate. The review by Gabici et al. writes
3
showing that PeV acceleration requires a sufficiently large product of size, flow speed, and magnetic field (Wilhelmi et al., 2024). The same review gives the maximal diffusive-shock-acceleration rate
4
and stresses that the actual maximum energy is set by balancing acceleration against losses or finite lifetime,
5
rather than by confinement alone (Wilhelmi et al., 2024).
For electrons, the synchrotron-loss limit is especially severe. Gabici et al. give
6
which is why electron PeV acceleration is disfavored in ordinary SNR shocks but remains plausible in relativistic pulsar environments or low-field cluster systems (Wilhelmi et al., 2024).
Operationally, PeVatron and super-PeVatron searches rely on a hard gamma-ray spectrum extending to the highest observed energies with little or no evidence for curvature. One formalization is the PeVatron Test Statistic,
7
with 8 PeV as a threshold for a proton PeVatron test; 9 would then constitute a robust hadronic PeVatron detection (Angüner, 2023). In practice, however, most current claims rest on a combination of UHE extension, absence of an observed cutoff, environmental plausibility, and exclusion of simple leptonic interpretations rather than on a single decisive statistic.
3. Source classes and acceleration environments
The search began with supernova remnants. SNRs remain attractive because of energetics, diffusive shock acceleration, and evidence for magnetic-field amplification, but Cristofari’s review is explicit that no clearly detected shell-type SNR has yet been shown to be an active proton PeVatron (Cristofari, 2021). A population-synthesis study of SNR PeVatrons argues that if SNRs produce Galactic CRs up to the knee, their active PeVatron phase should be short, of order a few $1$0 years, and detectable sources should be mostly young core-collapse remnants with compact angular sizes and hard spectra above $1$1 TeV (Cristofari et al., 2018). This suggests that SNR super-PeVatrons, if they exist, are likely rare and short-lived rather than the dominant observed UHE population.
A second class comprises young massive stellar clusters and superbubbles. Vink shows analytically that rapid second-order Fermi acceleration in superbubbles can accelerate particles to $1$2 eV within $1$3–$1$4 Myr, provided that $1$5, $1$6 at $1$7 TeV, and $1$8 pc, with diffusive escape limiting the maximum energy to a few PeV for typical superbubble sizes (Vink, 2024). A later time-dependent composition study goes further and argues that collective winds of massive star clusters, rather than individual stellar winds, naturally reproduce the common rigidity break of protons and helium near $1$9 PV reported by LHAASO, making collective cluster winds plausible dominant knee-region PeVatrons (Qiu et al., 29 May 2026). Morphological simulations with CTAO and the ASTRI Mini-Array predict that YMSCs should differ from TeV halos by an off-center radial emission peak and characteristic radial-profile parameters, providing a practical route to source identification in the LHAASO era (Bonollo et al., 24 Sep 2025).
Pulsars and PWNe define the clearest electron-PeVatron channel. The Crab Nebula is presented in multiple reviews as a convincing PeVatron, specifically an electron PeVatron, and many LHAASO UHE sources are associated with high-spin-down pulsars (Cristofari, 2021). The same literature stresses that this does not by itself solve the Galactic CR proton problem. For super-PeVatrons, PWNe are therefore a dual-use category: they prove that Galactic systems can reach the PeV domain, but they do not automatically imply hadronic acceleration to the knee.
Jet- and wind-powered binaries have emerged as additional candidates. The LHAASO attenuation study points out that microquasars such as Cyg X-3 and GRS 1915+105 can lie in the same attenuation regime as distant LHAASO PeVatron candidates, with strong deabsorption corrections required at the highest energies (Zhang et al., 2024). More aggressively, a ULX-wind-bubble model argues that super-Eddington ULX winds can accelerate protons to several PeV and, in the most powerful cases, to several tens of PeV, with SS 433 treated as a Galactic example whose 0 TeV morphology may be explained by hadronic emission from a wind bubble (Peretti et al., 2024). A related proposal extends the idea to super-accreting X-ray binaries more generally, defining super-PeVatrons as Galactic sources accelerating particles from above 1 PV to 2 PV and arguing that baryon-loaded trans-relativistic jets or winds with kinetic luminosity 3 can satisfy that requirement (Wang et al., 28 Jul 2025).
A speculative extension invokes exotic compact objects and black-hole–dark-sector systems. In that literature, ultra-spinning black-hole vortex-string systems, boson stars, axion stars, and Q-balls are proposed as “exotic PeVatrons,” with power controlled by quantized magnetic flux 4 or analogous relations and spin-down-like scaling 5, allowing PeV to beyond-PeV gamma rays in favorable parameter space (Addazi et al., 30 Sep 2025). The photon-production calculations there remain qualitative, but a plausible implication is that the phrase super-PeVatron is also beginning to denote a frontier where standard accelerator phenomenology and beyond-standard physics meet.
4. Representative Galactic systems
The present observational landscape is defined not by a single super-PeVatron, but by a small set of benchmark systems whose evidentiary status differs substantially. The table summarizes representative cases.
| System | Key observational fact | Interpretive status |
|---|---|---|
| Crab Nebula | PeV photons detected; robust PeVatron | Strong electron PeVatron; hadronic component uncertain |
| Galactic center ridge | Parent particle spectrum extends at least to 6 PeV at 95% CL in the classic H.E.S.S. analysis | Strong hadronic PeVatron candidate |
| HESS J1702-420A | 7, no curvature, 8–9 TeV bin at 0 | One of the strongest H.E.S.S. PeVatron candidates |
| G106.3+2.7 | Gamma rays extend to and above 1 TeV; morphology aligned with a molecular cloud | Strong SNR-related hadronic candidate, but not definitive |
| Cygnus X bubble | Eight events above 2 PeV within 3; implied protons at least 4 PeV | Explicit super-PeVatron claim |
| 1LHAASO J1857+0203u | 5 detection above 6 TeV | Good UHE accelerator candidate, but not especially extreme on current hadronic modeling |
The Crab Nebula is the cleanest proof that the Galaxy contains PeV particle accelerators. Reviews of Galactic PeVatrons treat it as a robust PeVatron, but specifically a leptonic one, while noting that a hadronic interpretation of its highest-energy hardening would imply proton cutoff energies of 7–8 PeV (Angüner, 2023). The Galactic center remains the strongest traditional hadronic case. In the H.E.S.S. interpretation summarized by the reviews, diffuse gamma rays from the central molecular zone imply a parent particle spectrum extending at least to 9 PeV at 95% CL, with absorption corrections potentially pushing the lower limit higher, although the identity of the central accelerator remains unresolved (Angüner, 2023).
HESS J1702-420A is the most compelling H.E.S.S. example of a compact hard component hidden inside a softer extended system. Deep H.E.S.S. observations resolved HESS J1702-420 into components A and B, with the physically decisive result being that A has $1$0, no detected curvature, and a $1$1 detection in the $1$2–$1$3 TeV interval; under hadronic priors, the 95% CL lower limit on the proton cutoff remains above $1$4 PeV (Giunti et al., 2021). The source therefore qualifies as a top-tier PeVatron candidate even though a leptonic origin cannot yet be excluded.
G106.3+2.7 has become the prototype of an SNR-related PeVatron candidate. Tibet AS$1$5MD reported gamma-ray emission extending to and above $1$6 TeV, with the $1$7 TeV centroid offset from PSR J2229+6114 and coincident with a molecular cloud, favoring hadronic emission (Amenomori et al., 2021). A later 12-year Fermi-LAT analysis strengthened this interpretation by showing a hard GeV–TeV spectrum with nondetection below 10 GeV, favoring a hybrid model in which electrons explain the radio-to-X-ray synchrotron emission while protons dominate the GeV–TeV gamma rays; the best-fit proton cutoff is $1$8 eV, with a 1$1$9 interval that straddles 0 PeV (Fang et al., 2022). The source is thus one of the best links between SNR environments and the knee-energy problem, though the papers stop short of claiming definitive proof of ongoing PeV acceleration at the current shock.
The most explicit super-PeVatron claim comes from the Cygnus X bubble. LHAASO reported a giant UHE gamma-ray bubble spanning at least 1, with 66 photon-like events above 2 TeV within 3, eight above 4 PeV, and two above 5 PeV in the central 6 core. The measured spectrum is well described by a log-parabola with
7
and the paper interprets the morphology, gas correlation, and PeV photons as evidence for “Super PeVatron(s)” accelerating protons to at least 8 PeV (Collaboration, 2023). In the present literature, this is the clearest published observational case for a Galactic source class entering a regime beyond the ordinary PeVatron threshold.
Not every 9 TeV source is equally extreme. The detailed LHAASO analysis of 1LHAASO J1857+0203u found $1$0 significance above $1$1 TeV, but the preferred hadronic fit for the HII-region interpretation gives $1$2 and $1$3 TeV, making it a credible UHE Galactic accelerator but not a compelling super-PeVatron in the strongest hadronic sense (Collaboration, 2024). This contrast is instructive: $1$4 TeV detection is necessary for a serious super-PeVatron discussion, but not sufficient.
5. Propagation, attenuation, and interpretive caveats
A central development in the field is the realization that Galactic propagation can no longer be neglected once one enters the LHAASO band. For source photons above tens of TeV, the relation between observed and intrinsic flux is
$1$5
with $1$6 receiving both ISRF and CMB contributions (Zhang et al., 2024). Using GALPROP v57 radiation fields, the LHAASO attenuation study finds that attenuation begins around $1$7 TeV, becomes material by $1$8 TeV for distant sources, and can reach values reported as 30% at $1$9 TeV and 80% at 3 PeV in the abstract, while the detailed text gives about 20% and 70% for the representative source 1LHAASO J1959+1129u. The paper explicitly notes this mismatch in presentation, but the qualitative conclusion is robust: observed curvature or softening in a distant inner-Galaxy source is not automatically intrinsic (Zhang et al., 2024).
This matters directly for super-PeVatron identification. Deabsorbed spectra of SS 433 and 1LHAASO J1959+1129u become harder and higher above 0 TeV than the measured spectra. The authors state that “the correction for the pair absorption can produces harder intrinsic spectra higher than 10 TeV than observed” and that “the hardening of the intrinsic spectrum around PeV may implicate different mechanism for the emissions” (Zhang et al., 2024). A plausible implication is that some of the most extreme Galactic accelerators are currently under-ranked because their 1 TeV flux is suppressed in transit.
A second caveat concerns air-shower interaction physics. The paper on the hadronic nature of high-energy photons argues that LHAASO’s nearly background-free 2 TeV photon sample opens a new regime for constraining the poorly known photoproduction cross section, with model differences exceeding 50% by 3 eV and already becoming relevant from 4 TeV (Sciascio, 2024). Because gamma/hadron separation in shower arrays depends strongly on muon content, uncertainties in photonuclear interactions propagate into flux reconstruction, cutoff determination, and source classification above 5 TeV (Sciascio, 2024).
These issues compound more familiar uncertainties: source confusion in crowded regions, uncertain distances, uncertain gas targets, and the persistent hadronic-versus-leptonic ambiguity emphasized in the reviews (Cristofari, 2021). For binaries and microquasars, internal absorption may be substantial but is not included in purely interstellar attenuation treatments (Zhang et al., 2024). For many LHAASO sources, angular resolution remains insufficient to decide whether the UHE emission is from the accelerator itself, from a nearby illuminated cloud, or from multiple overlapping emitters. Super-PeVatrons are therefore not defined by photon energy alone; they are defined by a chain of inference in which propagation, shower physics, and environment all matter.
6. Present frontier and future directions
The reviews now agree that gamma-ray astronomy has entered the PeVatron era and perhaps the super-PeVatron era, but without yet resolving the identity of the dominant Galactic sources above the knee (Cristofari, 2021). The field’s corrective lesson is that finding extreme accelerators is not identical to explaining the bulk CR spectrum. Even if some systems accelerate particles well beyond 6 PeV, it remains open whether those systems are numerous, efficient, and compositionally suitable enough to dominate the Galactic CR budget (Cristofari, 2021).
The most immediate experimental needs are improved angular resolution in crowded fields, improved spectral precision above 7–8 TeV, continued PeV-photon statistics, and multimessenger confirmation. Reviews and instrumentation studies converge on CTAO, SWGO, the ASTRI Mini-Array, continued LHAASO and HAWC exposure, and neutrino facilities such as IceCube-Gen2 and KM3NeT as the decisive tools for the next phase (Wilhelmi et al., 2024, Angüner, 2023, Bonollo et al., 24 Sep 2025). The YMSC morphology study makes this especially concrete: many LHAASO-like extended sources should become classifiable as cluster-wind systems or TeV halos only after tens to hundreds of hours of high-resolution IACT exposure (Bonollo et al., 24 Sep 2025).
The source landscape is also broadening. Massive star clusters and superbubbles are now serious candidates for long-lived hadronic super-PeVatrons; microquasar and ULX wind bubbles have become credible PeV-to-tens-of-PeV accelerators in dedicated models; super-accreting X-ray binaries have been proposed as a population extending to 9 PV rigidity; and exotic compact-object scenarios have been advanced for photons beyond the reach of conventional Galactic accelerators (Vink, 2024, Peretti et al., 2024, Wang et al., 28 Jul 2025, Addazi et al., 30 Sep 2025). This suggests that “super PeVatrons” are less a single source class than a frontier category in which several distinct acceleration environments compete.
At present, the most conservative statement is that the Galaxy contains multiple UHE gamma-ray sources, at least one robust electron PeVatron, several compelling hadronic PeVatron candidates, and at least one explicit observational super-PeVatron claim in the Cygnus X bubble (Collaboration, 2023). The strongest open question is no longer whether such extreme accelerators exist, but which of them are truly hadronic, which of them reach far beyond the ordinary PeV scale, and which—if any—dominate the Galactic cosmic-ray population around and above the knee.