Ultra-high energy event KM3-230213A as a cosmogenic neutrino in light of minimal UHECR flux models

This study proposes that the recently detected 220 PeV neutrino (KM3-230213A) is a cosmogenic neutrino consistent with minimal ultra-high energy cosmic ray flux models derived from Telescope Array data, while also satisfying constraints from gamma-ray observations.

The Gist

Imagine the universe as a giant, chaotic ocean. Most of the time, the waves are small and predictable. But occasionally, a "rogue wave" appears—a particle of energy so massive it defies our understanding.

In February 2023, a deep-sea telescope called KM3NeT (located in the Mediterranean Sea) spotted one of these cosmic rogue waves. It was a neutrino (a ghostly, nearly massless particle) with an energy of 220 PeV. To put that in perspective, if this particle were a baseball, it would be moving with the energy of a professional pitcher throwing a fastball, but compressed into a single subatomic speck. It is the most energetic neutrino ever detected.

The Mystery:
Here's the problem: This neutrino showed up alone. Other giant detectors around the world (like IceCube in Antarctica and the Pierre Auger Observatory in Argentina), which have been watching the sky for much longer, haven't seen any neutrinos this powerful. It's like a lighthouse in a storm seeing one massive wave, while every other lighthouse on the coast says, "The sea is calm."

This created a tension. Was this a fluke? Was it a signal from a nearby monster? Or was it something else entirely?

The Paper's Solution: The "Cosmic Cannonball" Theory
The authors of this paper propose a solution: This neutrino wasn't a fluke; it was a cosmogenic neutrino.

Think of the universe as a giant highway.

1. The Drivers (UHECRs): Ultra-High-Energy Cosmic Rays (UHECRs) are like super-fast cars (protons and atomic nuclei) zooming through space.
2. The Traffic (Background Light): The universe isn't empty; it's filled with a faint, ancient fog of light (photons) left over from the Big Bang.
3. The Crash: When these super-fast "cars" crash into the "fog," they don't just stop; they explode into debris. This debris includes pions, which decay into neutrinos and gamma rays.

The paper asks: If we assume the "cars" (cosmic rays) are mostly made of light stuff (like protons and helium), does the "debris" (neutrinos) match what KM3NeT saw?

The Twist: Two Different Maps
For years, scientists have argued about what these cosmic "cars" are made of.

• Map A (Pierre Auger): Suggests the cars are heavy trucks (heavy nuclei like iron).
• Map B (Telescope Array): Suggests the cars are mostly light motorcycles (protons and helium).

Previous studies tried to explain the KM3NeT event using Map A (heavy trucks). But to make the math work, they had to assume the "drivers" were coming from incredibly far away in the early universe, which felt like a stretch.

This paper says: "Let's try Map B." They used the Telescope Array's data, which suggests the cosmic rays are light. They ran a simulation of these light particles crashing into the cosmic fog.

The Results: A Perfect Match
Surprisingly, the math worked out beautifully without needing any "magic tricks."

• The Prediction: Their model predicted that we should see about one or two of these massive neutrinos in the KM3NeT detector.
• The Reality: KM3NeT saw exactly one.
• The Silence: The model also predicted that the other detectors (IceCube, Auger) should see zero or very few of these events in the same timeframe because they weren't looking at the right time or with the right sensitivity. This matches reality perfectly.

The "Ghost" Check (Gamma Rays)
If these cosmic cars are crashing, they should also produce a different kind of debris: high-energy gamma rays (light).

• The authors checked if their model produced too much "light pollution" (gamma rays) that would have been seen by the Fermi-LAT satellite.
• Result: No problem. The predicted gamma rays are faint enough to be hidden, fitting perfectly within the limits of what we can currently see.

The Big Picture
This paper is like solving a puzzle where everyone was trying to force a square peg into a round hole. By switching to a model where cosmic rays are light (protons) rather than heavy (iron), the authors showed that:

1. The single, massive neutrino seen by KM3NeT is likely a natural byproduct of cosmic rays traveling through the universe.
2. We don't need to invent new physics or assume extreme, exotic sources to explain it.
3. The "silence" from other detectors is actually expected, not a contradiction.

In Simple Terms:
The universe is a busy highway. We saw one massive crash (the neutrino). Other observers didn't see crashes because they weren't looking at the right spot or time. This paper proves that if the cars on the highway are mostly light (protons), the crash debris perfectly explains what we saw, without needing to break the laws of physics. It's a "minimalist" solution that fits the data like a glove.

Technical Summary

1. Problem Statement

On February 13, 2023, the KM3NeT experiment detected an exceptionally high-energy neutrino event, KM3-230213A, with an estimated energy of 220 PeV. This is the highest-energy neutrino ever recorded. However, this detection presents a significant tension with the null results from other major high-energy neutrino observatories (IceCube, Pierre Auger Observatory, and Baikal-GVD), which have much larger exposures but have detected no neutrinos in the same energy range (72 PeV – 2.6 EeV).

The origin of KM3-230213A is debated. Previous explanations include galactic sources, blazars, or cosmogenic origins (neutrinos produced by Ultra-High Energy Cosmic Rays, UHECRs, interacting with background photons). However, standard cosmogenic models based on Pierre Auger Observatory data (which suggest a heavy UHECR composition) struggle to explain the event without invoking extreme source evolution or additional proton components that conflict with other gamma-ray constraints.

The Core Question: Can the KM3-230213A event be explained as a cosmogenic neutrino using minimal UHECR flux models that assume a light mass composition (dominated by protons and light nuclei), as inferred by the Telescope Array (TA) experiment, without requiring exotic physics or extreme source evolution?

2. Methodology

The authors employed a multi-step simulation and statistical analysis approach:

• UHECR Flux Models: The study utilizes two specific UHECR source models derived from a combined fit of Telescope Array (TA) surface detector (SD) spectrum and fluorescence detector (FD) stereo composition data (Ref. [43]).

  • "Best-fit" Model: Represents the global minimum of the fit, dominated by Helium (99.2%) with a small Iron fraction.
  • "Local min." Model: Represents a local minimum of the fit, dominated by Protons (97.1%) with a small Iron fraction.
  • Assumptions: Both models assume a single source population with standard cosmological evolution ($SE(z) \propto (1+z)^3$ for $z < 1.5$) and a rigidity-dependent cutoff. They focus strictly on the high-energy component ($E > 10^{18.7}$ eV), ignoring the sub-ankle component.

• Propagation Simulation:

  • The CRPropa 3.2.1 framework was used to simulate the propagation of UHECRs through the intergalactic medium.
  • Interactions considered: Photopion production ($p\gamma$), photo-nuclear disintegration, and pair production on the Cosmic Microwave Background (CMB) and Extragalactic Background Light (EBL).
  • The simulation generated the resulting fluxes of cosmogenic neutrinos, UHECRs, and cascade gamma-rays.

• Data Comparison:

  • Neutrino Data: The predicted event rates were compared against observations from KM3NeT, IceCube (9-year and 12-year datasets), Pierre Auger, and Baikal-GVD.
  • Gamma-Ray Data: The predicted cosmogenic cascade gamma-ray flux (MeV–GeV range) was compared with Fermi-LAT Isotropic Diffuse Gamma-Ray Background (IGRB) data. The predicted Ultra-High Energy (UHE) gamma-ray flux was compared with upper limits from Auger and TA.

• Statistical Analysis:

  • The authors calculated the expected number of events ($n_{pred}$) for each experiment by integrating the predicted flux over the detector's effective area and exposure time.
  • A combined $p$-value and statistical significance were calculated using a likelihood ratio test (Poisson statistics) and Wilks' theorem approximation to assess the compatibility of the models with the "global neutrino observatory" (all experiments combined).

3. Key Contributions

• Re-evaluation of KM3-230213A Origin: The paper provides a rigorous test of the cosmogenic hypothesis using TA-derived models (light composition), contrasting with previous studies that relied on Auger-derived models (heavy composition).
• Minimal Model Validation: It demonstrates that a "minimal" UHECR model (single population, standard evolution, no extra proton component) is sufficient to explain the KM3-230213A event while remaining consistent with null results from other detectors.
• Multi-Messenger Consistency Check: The study simultaneously validates the models against neutrino data, Fermi-LAT IGRB constraints, and UHE gamma-ray limits, ensuring the proposed scenario is physically viable across different messengers.

4. Results

Neutrino Flux Compatibility

• KM3NeT-Only: Both the "Best-fit" and "Local min." models predict a neutrino flux that is lower than the flux required to produce exactly one event in KM3NeT alone, but the discrepancy is within the $2\sigma$ confidence level.
• Global Neutrino Observatory: When combining data from KM3NeT (1 event), IceCube (0 events), Auger (0 events), and Baikal-GVD (0 events), the models remain consistent with the observations at approximately the $2\sigma$ level.
  • Best-fit model: Combined $p$-value $\approx 0.019$ ($2.09\sigma$).
  • Local min. model: Combined $p$-value $\approx 0.035$ ($1.81\sigma$).
• Conclusion: The tension between KM3NeT's single detection and the null results of other experiments is significantly reduced when assuming a light-composition UHECR model.

Gamma-Ray Constraints

• Fermi-LAT IGRB: The predicted cosmogenic cascade gamma-ray fluxes for both models lie well below the Fermi-LAT IGRB measurements. This confirms that the models do not overproduce gamma rays in the GeV range, a common failure point for proton-dominated models with extreme evolution.
• UHE Gamma-Ray Limits: The predicted integral flux of UHE gamma-rays ($E > 10^{17}$ eV) is below current upper limits from Auger and TA. However, the "Best-fit" model prediction approaches the current limits at the highest energies ($\sim 10^{10}$ GeV), suggesting that future improvements in gamma-ray search sensitivity could probe these models.

5. Significance and Implications

• Resolution of Tension: The study suggests that the apparent conflict between KM3NeT's high-energy detection and the lack of detections by other experiments can be resolved without invoking Beyond the Standard Model (BSM) physics or exotic source populations. The discrepancy is likely due to statistical fluctuations and the specific composition of UHECRs (light vs. heavy).
• Validation of Telescope Array Models: The results lend strong support to the Telescope Array's inference of a light UHECR composition at the highest energies, which has historically been in tension with the Pierre Auger Observatory's heavy composition results.
• Future Prospects:
  • The "Best-fit" model is on the verge of being testable by current UHE gamma-ray experiments (Auger/TA) if their analysis efficiency improves by a factor of 3–5.
  • The study highlights the importance of joint analyses of neutrino and gamma-ray data to constrain UHECR source properties (composition, evolution, and maximum energy).
• Minimalism: The work emphasizes that complex adjustments (e.g., adding a separate proton population or extreme source evolution) are not necessary to explain the data, favoring Occam's razor in astrophysical modeling.

In summary, the paper successfully argues that the KM3-230213A event is a plausible cosmogenic neutrino arising from a standard population of UHECR sources with a light mass composition, consistent with all current multi-messenger constraints.