Evidence of cosmic-ray acceleration up to sub-PeV energies in the supernova remnant IC 443

Observations of the supernova remnant IC 443 by the Large High Altitude Air Shower Observatory (LHAASO) reveal high-energy gamma-ray emission extending beyond 30 TeV without a cutoff, providing compelling evidence that supernova remnant shocks can efficiently accelerate cosmic rays to sub-PeV energies.

The Gist

Imagine our galaxy, the Milky Way, as a vast, cosmic highway. On this highway, tiny particles called cosmic rays zoom around at incredible speeds, carrying more energy than anything we can create in our particle accelerators on Earth. For decades, scientists have been trying to figure out: Where do these super-fast particles get their boost?

The leading theory is that Supernova Remnants (SNRs)—the expanding, glowing debris fields left behind after massive stars explode—are the "cosmic accelerators" responsible. But there's a catch: while we know they can speed up particles, we haven't been able to prove they can push them all the way to the "knee" of the energy spectrum (a massive energy level called the PeV, or peta-electronvolt). Finding a source that does this is like finding the "Holy Grail" of astrophysics, often called a PeVatron.

This paper is about a major breakthrough in that search, using a giant cosmic-ray detector in China called LHAASO to study a specific cosmic graveyard known as IC 443.

The Detective Work: Two Sources in One

Think of IC 443 not as a single object, but as a messy crime scene where two different suspects might be hiding. When the LHAASO team looked at the data, they didn't just see one blob of energy; they resolved two distinct sources of high-energy gamma rays (the "smoke" left behind by the cosmic rays):

1. The Compact Source (C0): This is a tight, bright spot right where the supernova shockwave is crashing into a dense cloud of gas. It's like a high-pressure hose spraying water into a brick wall.
2. The Extended Source (C1): This is a huge, fuzzy halo of energy stretching out far beyond the main shockwave. It's like the mist and spray that drifts far away from the hose.

The Smoking Gun: Protons vs. Electrons

In the universe, high-energy light (gamma rays) can be made in two main ways:

• The "Leptonic" Way: Fast-moving electrons bounce off light waves (like a billiard ball hitting a cue ball).
• The "Hadronic" Way: Fast-moving protons (the heavy nuclei of atoms) crash into other gas particles, creating a particle called a "pion" that decays into gamma rays. This is the "smoking gun" for cosmic ray acceleration because it proves heavy particles are being accelerated, not just light electrons.

The team found that the Compact Source (C0) has a very specific signature (a "bump" in its energy spectrum) that only happens when protons crash into gas. This confirms that IC 443 is indeed a cosmic particle accelerator.

The Big Discovery: Breaking the Energy Barrier

Here is the most exciting part. Scientists had a theory that supernova shocks could only accelerate protons up to about 10–50 TeV (Tera-electronvolts). But when the LHAASO team looked at the energy of the protons in the Compact Source, they found no sign of a limit.

They calculated that the protons are being pushed to at least 300 TeV, and likely even higher (approaching the PeV range).

The Analogy:
Imagine a rollercoaster designed to go 100 mph. For years, physicists thought the track ended there. But this new data shows the rollercoaster is actually going 300 mph, and the track keeps going! This proves that supernova shocks are much more powerful than we thought, capable of acting as the "PeVatrons" we've been searching for.

What About the Fuzzy Halo?

The Extended Source (C1) is a bit of a mystery. It could be:

• Leaking Protons: Protons that escaped the main shockwave and are wandering through the galaxy, crashing into gas far away.
• A Different Accelerator: It might be related to a nearby pulsar (a spinning dead star) or even a different, older supernova remnant nearby.
• Electrons: It might just be fast electrons doing the "billiard ball" dance, rather than protons.

The data is consistent with both ideas, so the team needs more observations to solve this specific puzzle.

Why This Matters

This paper is a landmark because it provides the first solid evidence that supernova remnants can accelerate particles to "sub-PeV" energies (just under the PeV limit). It's like finally finding the engine that powers the cosmic highway.

In summary:

• The Problem: We didn't know if supernovae could create the most energetic particles in the galaxy.
• The Tool: LHAASO, a giant detector array in the mountains of China, acted like a high-resolution camera.
• The Result: They looked at IC 443 and found a "compact" accelerator pushing protons to at least 300 TeV.
• The Takeaway: Supernova explosions are the cosmic factories that build the universe's most energetic particles, and we have finally caught them in the act.

Technical Summary

Here is a detailed technical summary of the paper "Evidence of cosmic-ray acceleration up to sub-PeV energies in the supernova remnant IC 443" by the LHAASO Collaboration.

1. Problem Statement

Supernova remnants (SNRs) are the leading candidates for the origin of Galactic cosmic rays (CRs), capable of accelerating particles via diffusive shock acceleration. However, it remains uncertain whether SNRs can accelerate protons to the "knee" region of the CR spectrum (several PeV), a threshold required to explain the bulk of Galactic CRs. While interactions between SNRs and dense molecular clouds (MCs) are expected to produce characteristic $\pi^0$-decay gamma-ray signatures, previous observations of IC 443 (a middle-aged SNR interacting with MCs) by MAGIC, VERITAS, and HAWC had not definitively confirmed acceleration up to sub-PeV energies due to limited energy reach or spectral cutoffs. The specific maximum energy of accelerated protons in IC 443 and the nature of its extended gamma-ray emission remained open questions.

2. Methodology

The study utilizes data from the Large High Altitude Air Shower Observatory (LHAASO), a hybrid detector array in China comprising the Kilometer Square Array (KM2A) and the Water Cherenkov Detector Array (WCDA).

• Data Selection:
  • WCDA: March 5, 2021 – July 31, 2024 (~1136 days livetime).
  • KM2A: July 20, 2021 – December 31, 2024 (~1228 days livetime).
• Background Modeling:
  • A Region of Interest (ROI) was defined as a fan-shaped area centered on the Geminga pulsar (6° away from IC 443) to account for the extended Geminga halo.
  • Diffuse gamma-ray emission was modeled using gas templates traced by PLANCK dust opacity and gas surveys.
  • A 3D-likelihood method was employed to simultaneously fit the morphology and spectrum of the target sources, the Geminga halo, and the diffuse background.
• Analysis Strategy:
  • The analysis tested single-source vs. two-source hypotheses for the IC 443 region.
  • Spectral fitting was performed using Power-Law (PL) and Exponentially Cutoff Power-Law (ECPL) models.
  • Hadronic and Leptonic models were constructed to interpret the spectral energy distributions (SEDs) and constrain the maximum energy of accelerated protons.

3. Key Contributions and Results

A. Morphological Decomposition

The LHAASO data resolved the gamma-ray emission from the IC 443 region into two distinct components:

1. Source C0 (Compact):
   • Significance: $10.5\sigma$.
   • Morphology: Point-like (upper limit on extension $< 0.27^\circ$).
   • Association: Coincides with the Fermi-LAT compact source (4FGL J0617.2+2234e) and the densest part of the shocked molecular clouds.
   • Spectrum: Well-described by a Power-Law ($\Gamma \approx 2.95$) extending beyond 30 TeV with no apparent spectral cutoff.
2. Source C1 (Extended):
   • Significance: $13.1\sigma$.
   • Morphology: Extended with a radius $R_{39} = 0.67^\circ \pm 0.07^\circ$.
   • Association: Overlaps with a newly reported Fermi source (2FGES J0618.3+2227) and potentially the SNR G189.6+3.3 or the pulsar wind nebula (PWN) CXOU J061705.3+222127.
   • Spectrum: Best fit by an ECPL model with a cutoff energy $E_{cut} \approx 19.65 \pm 8.67$ TeV, though a broken power-law (BPL) model without a cutoff is also statistically viable.

B. Evidence for Sub-PeV Acceleration (Source C0)

Assuming a hadronic origin (proton-proton collisions producing $\pi^0$ decay) for Source C0:

• The lack of a spectral cutoff in the LHAASO data (up to ~30 TeV) allows for a robust constraint on the maximum proton energy.
• By fitting the proton spectrum (modeled as a broken power-law with an exponential cutoff) to the combined Fermi-LAT and LHAASO data, the authors derived a 95% lower limit for the proton cutoff momentum ($p_{cut}$) of $\sim 300$ TeV.
• This implies that the SNR shock in IC 443 is efficiently accelerating protons to sub-PeV energies, challenging standard test-particle acceleration models which typically predict limits in the tens of TeV range.

C. Interpretation of Source C1

The nature of Source C1 remains ambiguous but offers two plausible scenarios:

• Hadronic Scenario: If associated with SNR G189.6+3.3, the higher break energy ($\sim 22$ TeV) compared to C0 suggests a higher shock speed or different acceleration efficiency. If associated with IC 443, the spectral difference may result from particle propagation effects (diffusion of escaping protons).
• Leptonic Scenario: The emission could be produced by inverse Compton scattering of accelerated electrons. This model fits the data with an electron cutoff energy of $\sim 125$ TeV, consistent with cooling timescales in the ambient magnetic field.

4. Significance

• Confirmation of PeVatron Candidates: This study provides compelling evidence that SNR shocks can accelerate cosmic rays to sub-PeV energies, supporting the hypothesis that SNRs are the primary sources of Galactic CRs up to the "knee."
• Resolution of Morphology: The ability to separate the compact, shock-interacting component (C0) from the extended component (C1) clarifies previous conflicting results from other observatories (MAGIC, VERITAS, HAWC).
• Physical Constraints: The derived lower limit of 300 TeV for proton acceleration necessitates mechanisms beyond simple diffusive shock acceleration, such as magnetic field amplification, turbulence generation, or specific shock configurations.
• Multi-wavelength Synergy: The work successfully bridges the gap between GeV (Fermi-LAT) and multi-TeV (LHAASO) observations, demonstrating the power of wide-field, high-sensitivity detectors in studying complex SNR-MC interaction regions.

In conclusion, the LHAASO observations of IC 443 represent a milestone in high-energy astrophysics, offering the most stringent evidence to date that a specific SNR shock is capable of accelerating protons to energies approaching the PeV scale.