First detection of ultra-high energy emission from gamma-ray binary LS I +61 303

Using data from the Large High Altitude Air Shower Observatory (LHAASO), researchers report the first detection of ultra-high-energy gamma-ray emission (exceeding 100 TeV) from the binary system LS I +61 303, providing compelling evidence for extreme particle acceleration and orbital modulation that suggests a composite leptonic and hadronic emission scenario.

The Gist

Here is an explanation of the paper, translated from complex astrophysics into a story you can understand over a cup of coffee.

The Big Picture: Catching the "Super-Bullets"

Imagine the universe as a giant, dark ocean. Most of the time, it's quiet. But every now and then, a cosmic "storm" happens, shooting out particles so fast and so energetic that they are like super-bullets traveling at nearly the speed of light.

Scientists have long known about a specific "storm system" in our galaxy called LS I +61°303. It's a cosmic dance between two partners:

1. A massive, young star (a Be-type star) that is spinning so fast it's flinging a giant ring of gas around itself (like a figure skater with a skirt of water).
2. A tiny, invisible partner (likely a neutron star) that is a dense, dead core of a star, orbiting the big one.

For years, telescopes could see the "rain" (lower energy light) from this storm. But this new paper, using a giant detector in China called LHAASO, has finally caught the "super-bullets" (Ultra-High Energy gamma rays). These are particles with energies over 100 trillion electron volts—energies so high they were previously thought impossible for this kind of system to produce.

The Detective Work: How They Found It

The LHAASO observatory is like a massive, high-tech net spread across a mountain in Sichuan, China. When a super-bullet from space hits the atmosphere, it creates a shower of secondary particles, like a splash in a pond. LHAASO catches these splashes.

The team looked at the data from two different parts of their net:

• The "Wide Net" (WCDA): Caught a huge number of medium-energy bullets.
• The "Deep Net" (KM2A): Caught the rare, super-powerful bullets.

The Result: They found a clear signal. It wasn't just background noise; it was a distinct "ping" coming from the direction of LS I +61°303. They even found one single particle with an energy of 267 TeV (that's 267 trillion electron volts!). To put that in perspective, if a single proton had that much energy, it would be like a baseball thrown at 100 mph, but compressed into a subatomic speck.

The Mystery: Why Does It Happen?

The most exciting part of the paper isn't just that they found these bullets, but when they found them.

Imagine the two stars dancing in an ellipse. Sometimes they are far apart (apastron), and sometimes they are very close (periastron).

• The Old Theory: Scientists thought the "super-bullets" would be shot out when the stars were close together, crashing into each other's gas.
• The New Clue: The data shows something weird. The lower-energy bullets come out when the stars are close. But the super-bullets (the ones over 100 TeV) seem to cluster when the stars are farther apart.

The Solution: A Two-Engine Car

The authors propose a clever explanation using a "Two-Engine" analogy.

Think of the binary system as a car with two engines:

1. The Electric Engine (Leptonic): This uses electrons. It's great for speed but runs out of steam quickly when things get too hot or crowded. This engine powers the lower-energy light we see when the stars are close together.
2. The Diesel Engine (Hadronic): This uses protons (heavy particles). It's slower to start but can handle the heavy lifting and go much further.

The Theory:

• When the stars are close, the "Electric Engine" revs up, creating a burst of light.
• But as the stars move apart, the environment changes. The "Electric Engine" gets choked by the intense radiation (like trying to run a fan in a hurricane). However, the "Diesel Engine" (protons) keeps going. Because the magnetic fields are weaker when the stars are apart, the protons don't lose energy as fast. They get accelerated to insane speeds and shoot out those Ultra-High Energy bullets.

Why Does This Matter?

This discovery is a game-changer for two reasons:

1. The "PeVatron" Candidate: It suggests that LS I +61°303 is a "PeVatron"—a natural particle accelerator in our galaxy capable of boosting particles to energies we can't even create in our biggest man-made machines (like the Large Hadron Collider). It proves that nature has found a way to build a super-rocket engine.
2. The Cosmic Mystery: It solves a puzzle about how these systems work. They aren't just simple collisions; they are complex machines where different types of particles take the lead at different times and distances.

The Bottom Line

The LHAASO team has opened a new window into the universe. They looked at a cosmic dance partner system and realized it's not just a simple waltz; it's a high-speed chase where the stars switch drivers depending on how far apart they are. They found the "smoking gun" of ultra-powerful particles, proving that our galaxy is full of hidden factories capable of creating the most extreme energy in the universe.

Technical Summary

Here is a detailed technical summary of the paper "First detection of ultra-high energy emission from gamma-ray binary LS I +61◦303" by the LHAASO Collaboration.

1. Problem and Scientific Context

Gamma-ray binaries are systems consisting of a compact object (neutron star or black hole) orbiting a massive O- or Be-type star. They are known for emitting variable gamma rays, modulated by their orbital periods, driven by interactions between the compact object's relativistic outflow and the companion star's wind or disk.

• The Gap: While these systems have been detected in the TeV range (up to ~10–20 TeV) by instruments like MAGIC and VERITAS, it remains unclear whether they can accelerate particles to Ultra-High Energies (UHE; >100 TeV) or even PeV energies.
• The Question: Can gamma-ray binaries act as Galactic "PeVatrons" (sources capable of accelerating cosmic rays to PeV energies)? Specifically, does LS I +61◦303, a prototypical gamma-ray binary, emit photons beyond 100 TeV, and what physical mechanisms (leptonic vs. hadronic) drive this emission?

2. Methodology

The study utilizes data from the Large High Altitude Air Shower Observatory (LHAASO), located in Sichuan, China. The analysis combines data from two sub-arrays:

• WCDA (Water Cherenkov Detector Array): Used for the energy range 1.4 – 30.5 TeV. Data collected from March 2021 to July 2024. Energy estimation relies on the number of triggered detectors ($N_{hit}$).
• KM2A (Kilometer Square Array): Used for the energy range 25 TeV to several hundred TeV. Data collected from December 2019 to July 2024. Energy reconstruction is based on particle density at 50m from the shower axis.
• Analysis Techniques:
  • 3D Likelihood Fitting: A joint fit was performed over a $4^\circ$ Region of Interest (ROI) to simultaneously model the source (LS I +61◦303), a nearby source (1LHAASO J0249+6022), and Galactic Diffuse Emission (GDE).
  • Background Estimation: The "direct integration method" was used to estimate cosmic-ray background.
  • Orbital Modulation: Events were binned into orbital phases ($\phi$) using the known period $P = 26.496$ days. A likelihood-ratio test compared constant flux vs. phase-dependent flux hypotheses.
  • High-Energy Event Selection: For the >100 TeV regime, specific photon-like events were identified against the background using the Li & Ma formula.

3. Key Contributions and Results

A. First UHE Detection

• Significance: The paper reports the first definitive detection of gamma-ray emission from LS I +61◦303 in the UHE regime.
• Statistical Significance:
  • WCDA (1.4–30.5 TeV): $9.2\sigma$ detection.
  • KM2A (25–267 TeV): $6.2\sigma$ detection.
  • >100 TeV Sub-sample: 16 photon-like events detected against an estimated 5.1 background events, yielding a $3.8\sigma$ significance.
• Highest Energy Event: A single photon-like event with an estimated energy of **$267 \pm 61$TeV** was detected at an angular distance of$0.09^\circ$ from the source, consistent with the detector's Point Spread Function (PSF).

B. Spectral Properties

• The combined spectrum from 1 TeV to hundreds of TeV is well-described by a single power law:
  $$dN/dE = N_0 (E/10 \text{ TeV})^{-\Gamma}$$
• Best-fit parameters:
  • Photon index: $\Gamma = 3.00 \pm 0.05$.
  • Flux normalization: $N_0 = (2.18 \pm 0.20) \times 10^{-15} \text{ TeV}^{-1} \text{ cm}^{-2} \text{ s}^{-1}$.
  • Integral flux ($>1$ TeV): $1.10 \times 10^{-12} \text{ cm}^{-2} \text{ s}^{-1}$ (approx. 3.8% of the Crab Nebula flux).
• No significant spectral cutoff was observed up to 267 TeV, challenging standard acceleration models.

C. Orbital Modulation and Energy Dependence

• Low Energy (1.4–30.5 TeV): Emission is detected across most phases, peaking near $\phi \approx 0.6$ (near apastron/inferior conjunction). Significance of modulation: $\sim 2.6\sigma$.
• Medium Energy (25–100 TeV): Emission is more confined to phases $\phi = 0.3–0.6$. Significance of modulation: $3.9\sigma$.
• Ultra-High Energy (>100 TeV): Due to low statistics, data was split into two bins. There is a hint ($2.4\sigma$) that UHE photons cluster in the **$\phi = 0.6–0.2$ range** (near apastron), which differs from the lower-energy distribution that peaks closer to periastron.

4. Physical Interpretation and Discussion

The authors propose a composite lepto-hadronic scenario to explain the observations:

1. Leptonic Component (Dominant at lower energies):

   • Electrons upscatter stellar photons via Inverse Compton (IC) scattering.
   • This explains the X-ray and lower-energy TeV emission.
   • However, pure leptonic models struggle to explain >100 TeV emission due to the Klein-Nishina suppression (reduced IC cross-section at high energies) and synchrotron cooling limits, which typically cap electron energies around 40–130 TeV.

2. Hadronic Component (Dominant at UHE):

   • Protons are accelerated and interact with dense matter (stellar wind or Be-disk) via $p-p$ collisions, producing neutral pions ($\pi^0$) that decay into gamma rays.
   • Disk Crossing: The compact object likely crosses the Be-star's circumstellar disk near periastron, creating a high-density environment ($n \sim 10^{11} \text{ cm}^{-3}$) that boosts hadronic interaction efficiency.
   • Energy Budget: A purely hadronic model for the entire spectrum would require a proton power exceeding the pulsar's spin-down luminosity. However, a mixed scenario works: electrons power the bulk of the emission (<25 TeV), while protons power the high-energy tail (>25 TeV). The required proton power ($\sim 2.5 \times 10^{33}$ erg s$^{-1}$) fits within the pulsar's spin-down budget ($\sim 6 \times 10^{33}$ erg s$^{-1}$).

3. Orbital Visibility:

   • The detection of UHE photons near apastron suggests these photons escape a region where $\gamma\gamma$ absorption (pair production) is minimized, whereas emission near superior conjunction might be absorbed by the dense stellar radiation field.

5. Significance

• New Window into Acceleration: This is the first time a gamma-ray binary has been confirmed to emit photons up to ~267 TeV, proving these systems are capable of accelerating particles to energies far beyond previous limits.
• PeVatron Candidate: The results strongly suggest LS I +61◦303 is a Galactic PeVatron candidate, capable of accelerating protons to PeV energies.
• Mechanism Validation: The energy-dependent orbital modulation provides crucial evidence for a mixed lepto-hadronic acceleration mechanism, where the dominant particle population shifts with energy.
• Future Outlook: The detection sets a precedent for future multi-wavelength campaigns to study the interplay between shock acceleration, magnetic fields, and dense matter in compact binary systems.

In conclusion, the LHAASO collaboration has fundamentally expanded the known energy range of gamma-ray binaries, providing compelling evidence for extreme particle acceleration and a complex, energy-dependent emission mechanism in LS I +61◦303.