---
title: 'KM3-230213A: Ultra-High-Energy Neutrino Event'
url: https://www.emergentmind.com/topics/km3-230213a
type: topic
---

# KM3-230213A: Ultra-High-Energy Neutrino Event

Searching arXiv for recent papers on KM3-230213A to ground the article in the literature.
KM3-230213A is an ultra-high-energy neutrino event reported by the KM3NeT Collaboration and widely treated as the most energetic neutrino observed to date, with an inferred neutrino energy of approximately \(220\) PeV and a reconstructed muon energy of \(120^{+110}_{-60}\,\mathrm{PeV}\) [2502.09548]. It has become a focal object in multimessenger astrophysics because a single event at this energy probes the extreme tail of neutrino production models, the proton content and evolution of ultra-high-energy cosmic-ray source populations, and several beyond-standard-model scenarios, while also exposing a quantitative tension with null results from IceCube and the Pierre Auger Observatory [2605.07641]. The event has therefore been used both as a phenomenological source probe and as a constraint on Lorentz-invariance violation, in-vacuo dispersion, dark-matter decay, and related frameworks [2502.18256].

## 1. Discovery, reconstruction, and basic observational status

KM3-230213A was recorded by the 21-line configuration of KM3NeT/ARCA on 2023-02-13 and published in February 2025 [2605.07641]. It was detected as an almost horizontal or near-horizontal track event, with best-fit reconstructed direction \(\mathrm{RA}=94.3^\circ\), \(\mathrm{Dec}=-7.8^\circ\) and a reported \(68\%\) containment radius \(R(68\%)=1.5^\circ\); a \(90\%\) containment radius of \(2.2^\circ\) is quoted in VERITAS follow-up work [2502.13093]. The event is associated with a reconstructed muon energy \(E_\mu = 120^{+110}_{-60}\,\mathrm{PeV}\), while the inferred neutrino energy is typically given as \(220^{+570}_{-110}\,\mathrm{PeV}\), with a 68% confidence interval of \([110\,{\rm PeV},\,790\,{\rm PeV}]\) and a 90% confidence interval of \([72\,{\rm PeV},\,2.6\,{\rm EeV}]\) [2502.18256].

The event is treated as astrophysical rather than atmospheric because the amount of material traversed along the reconstructed trajectory strongly disfavors an atmospheric-muon origin. One review states that the path implies roughly \(300\ \mathrm{km\ water\ equivalent}\), making an atmospheric muon explanation “highly unlikely” [2507.16392]. A dedicated PBR study similarly argues that atmospheric muons are not viable once the required water-equivalent depth and survival probability are taken into account [2507.02246].

Its exceptional status derives not only from absolute energy, but also from comparison with previous neutrino detections. One Lorentz-invariance-violation analysis emphasizes that the event energy is “almost two orders of magnitude above the previous most energetic neutrino seen by IceCube” [2502.09548]. Another review describes it as the first detected astrophysical neutrino above \(100\) PeV and the first direct detection of a neutrino in the ultra-high-energy domain, defined there as \(E_\nu \gtrsim 50\ \mathrm{PeV}\) [2605.07641].

## 2. Flux interpretation and tension with other observatories

KM3-230213A quickly became central to the global ultra-high-energy neutrino landscape because a naive flux inference from a single event appears difficult to reconcile with non-observations elsewhere. One synthesis reports a single-flavor flux normalization
\[
E^2\Phi = 7.5^{+13.1}_{-4.7}\times 10^{-10}\ \mathrm{GeV\ cm^{-2}\ s^{-1}\ sr^{-1}}
\]
from a joint fit involving KM3NeT, IceCube, and Auger, with a tension of about \(2.5\sigma\) relative to the null detections by IceCube and Auger [2605.07641]. Another review summarizes a higher event-implied flux level,
\[
5.8^{+10.1}_{-3.7}\times 10^{-8}\ \mathrm{GeV\ cm^{-2}\ s^{-1}},
\]
and then notes that the discrepancy softens once a joint isotropic fit is used [2507.16392].

Several papers frame this mismatch in terms of source class. Under a diffuse isotropic interpretation, the tension with IceCube is quoted as \(3.5\sigma\), while transient-source interpretations reduce it to \(2.6\sigma\) [2507.02246]. A review gives a similar hierarchy: \(3.1\)–\(3.6\sigma\) for diffuse or cosmogenic scenarios, \(2.9\sigma\) for stationary point sources, and \(2.0\sigma\) for transient point sources [2507.16392]. This has made transience, anisotropy, or detector-geometry dependence recurrent themes in the literature.

That tension is methodologically important because it forces model comparisons to use not just the observed event, but also the absence of comparable events in IceCube and Auger. This is explicit in cosmogenic analyses, GRB-blastwave fits, and several beyond-standard-model papers, all of which treat KM3-230213A not as an isolated datum but as a constraint coupled to wider exposure and null results [2603.13476].

## 3. Source classes and astrophysical interpretations

Astrophysical interpretations of KM3-230213A divide broadly into Galactic, extragalactic transient, cosmogenic, and population-level diffuse scenarios. A review states that a Galactic origin is effectively ruled out because counterpart searches found no powerful Galactic emitters and no relevant gamma-ray sources in the error region seen by HAWC or LHAASO [2605.07641]. Another review similarly describes a Galactic origin as highly unlikely [2507.16392].

Extragalactic transient scenarios remain prominent because they reduce the KM3NeT–IceCube tension. In hadronic models, however, a neutrino source should also produce gamma rays, leading to dedicated follow-up studies. A Fermi-LAT cascade search found no convincing GeV counterpart and therefore constrained combinations of source redshift and intergalactic magnetic field, excluding roughly \(B \lesssim 10^{-14}\,\mathrm{G}\) for \(z<1\) in its benchmark transient hadronic picture [2503.16606]. VERITAS follow-up observations likewise found no statistically significant very-high-energy gamma-ray excess and derived a \(99\%\) CL upper limit at the event position of
\[
F(>550\ \mathrm{GeV}) < 9.49 \times 10^{-13}\ \mathrm{cm}^{-2}\,\mathrm{s}^{-1},
\]
for an assumed spectrum with \(\Gamma=-2.5\) [2509.25422].

Within the localization region, several papers discuss blazars as candidate counterparts. A review notes that 17 blazars were identified within the uncertainty region, with three showing flaring behavior near the event time, but no conclusive association was established [2605.07641]. A source-model study argues that PMN J0606-0724 is the most plausible counterpart because of a radio flare peaking 5 days after the neutrino, while interpreting the event through a jet–red-giant interaction and shock acceleration scenario [2511.01729]. In that framework, the baryon injection is attributed to the red giant, but the dominant photon field setting the neutrino energy scale is argued to be colder infrared emission from the dusty torus rather than the red-giant photosphere, since the latter would more naturally yield \(\sim 1\)–10 PeV neutrinos instead of \(\sim 220\) PeV [2511.01729].

Cosmogenic interpretations have also been developed in detail. One study fits a two-population UHECR model and finds that, if the event is interpreted using only KM3NeT exposure, the proton-source evolution must be strongly positive, with \(5.2 \le m^p \le 6.6\) at 68% CL, while inclusion of Auger and IceCube null detections relaxes this to \(-1.2 \le m^p \le 4.1\) [2603.13476]. Another multimessenger transport calculation argues that a steep sub-ankle proton component is disfavored because it saturates the Fermi-LAT isotropic gamma-ray background, whereas a hard proton spectrum extending beyond \(10^{20}\,\mathrm{eV}\) with evolution \(\propto (1+z)^3\) can reproduce KM3-230213A without violating gamma-ray or neutrino limits, provided the proton fraction is \(\lesssim 10\%\) at \(3\times 10^{19}\,\mathrm{eV}\) [2507.11993]. A different cosmogenic-point-source study proposes that a nearby transient cosmogenic-neutrino source with weak intergalactic magnetic field could be more favorable than diffuse cosmogenic emission, particularly because it reduces the contrast with IceCube non-detection [2504.10378].

Long-duration GRB blastwaves have likewise been tested as a population explanation. A Poisson-likelihood analysis using the joint exposure of KM3NeT/ARCA, IceCube, and Auger constrains the baryon loading to
\[
f_b \le \{392,131,39,13\}
\]
at \(90\%\) confidence for ISM densities \(n_0=\{1,3,10,30\}\,\mathrm{cm^{-3}}\), while for wind-like environments it finds
\[
f_b \le \{20,50,100\}
\]
for \(A_*=\{0.05,0.06,0.07\}\) [2509.14895]. This suggests that GRB blastwaves remain viable contributors to the diffuse UHE neutrino background, but only in restricted parameter regions [2509.14895].

## 4. Fundamental-physics tests enabled by the event

Because many Lorentz-violating effects scale steeply with energy, KM3-230213A immediately became a powerful probe of Lorentz invariance in the neutrino sector. One analysis considers superluminal neutrino Lorentz-invariance violation with modified dispersion relation
\[
E^2=m^2+k^2\left(1+s_n\left(\frac{k}{\Lambda}\right)^n\right),
\]
for \(n=1,2\), and uses the absence of neutrino splitting \(\nu\to \nu\nu\bar\nu\) to derive lower limits on the LIV scale [2502.09548]. Under an extragalactic source-distance assumption \(L=10\) Mpc, it finds
\[
\Lambda = \begin{cases}
5.4\times10^{30}\,\mathrm{GeV}, & n=1,\\
3.5\times10^{19}\,\mathrm{GeV}, & n=2,
\end{cases}
\]
and identifies neutrino splitting, rather than \(\nu\to \nu e^+e^-\), as the dominant constraint [2502.09548].

A second LIV study adopts an isotropic modified-dispersion framework and propagates neutrino spectra including the vacuum decay channels \(\nu\to \nu e^+e^-\) and \(\nu\to \nu\nu\bar\nu\). It concludes that the most conservative acceptable quadratic limit is
\[
\Lambda_2 > 5.0\times10^{19}\ {\rm GeV}\qquad (90\%\ {\rm C.L.}),
\]
while stating explicitly that the existing data do not provide a high-confidence constraint on the linear \(n=1\) scale in that setup [2502.18256]. A review further summarizes related velocity-based constraints as \(\delta=c_\nu^2-1<4.2\times10^{-22}\) [2507.16392].

KM3-230213A has also been used as a case study for in-vacuo dispersion. One paper investigates whether the neutrino could be associated with GRB090401B, observed about 14 years earlier, under a linear delay law
\[
\Delta t = D(z)\frac{E}{M_{QG}},
\]
and reports directional overlap, energy consistency, and a Monte Carlo p-value of \(0.0092\) for the GRB090401B–KM3-230213A pairing [2502.13093]. The authors explicitly state that this is “intriguing but far from conclusive,” and frame the exercise as proof of principle for archival transient searches with decade-scale delays rather than as evidence for quantum gravity [2502.13093].

## 5. Dark-matter and other beyond-standard-model interpretations

KM3-230213A has generated an extensive literature on superheavy or very heavy dark matter. The simplest class of models uses two-body dark-matter decay to neutrinos, exploiting the kinematic relation \(E_\nu \simeq M_{\rm DM}/2\). A paper based on a singlet-fermion dark matter extension of type-I seesaw interprets the event through \(\chi\to \nu h\), deriving an allowed dark-matter mass range
\[
1.7\times 10^8\ {\rm GeV} \lesssim M_{\rm DM} \lesssim 5.5\times10^9\ {\rm GeV},
\]
with corresponding lifetimes
\[
6.3\times10^{29}\ {\rm s}\ \text{to}\ 3.6\times10^{29}\ {\rm s},
\]
while also predicting stochastic gravitational waves from the disappearance of domain walls formed when a \(Z_2\) symmetry is spontaneously broken [2503.04464].

Another paper proposes a vector dark matter model based on a new \(U(1)_X\) gauge symmetry, with benchmark dark matter mass \(m_X=440\) PeV and viable lifetime range
\[
7.3\times 10^{28}\ \mathrm{s} \lesssim \tau_{\rm DM} \lesssim 2.9\times 10^{30}\ \mathrm{s}.
\]
In that model, the same \(U(1)_X\) breaking produces cosmic strings with predicted tension
\[
4.5 \times 10^{-11} \lesssim G\mu \lesssim 1.2 \times 10^{-10},
\]
which the authors describe as consistent with recent pulsar timing array observations [2507.21534].

More cosmological dark-sector scenarios have also been explored. One study links KM3-230213A and high-energy IceCube events to metastable superheavy dark matter produced through primordial black hole evaporation, finding viability for
\[
m_{\rm DM}\sim 10^8-10^{10}\,\mathrm{GeV},
\quad
M_{\rm BH,0}\sim 10-5\times10^8\,\mathrm{g},
\quad
\beta\sim 10^{-30}-10^{-14},
\]
with a benchmark \(m_{\rm DM}=5000\) PeV and \(a_d=10^{-2}\) [2510.26126]. These models do not require a local counterpart and instead interpret the event as part of a diffuse cosmological neutrino spectrum from late decays [2510.26126].

Not all dark-matter interpretations remain favorable. A detector-level analysis of heavy dark matter decay as the origin of the event finds preferred masses larger than about \(100\) PeV and best-fit lifetimes in the range \(10^{26}\)–\(10^{27}\) s, but concludes that these regions are in tension with existing gamma-ray and neutrino bounds [2606.09986]. It identifies best-fit points such as \((m_{\rm DM},\tau_{\rm DM})=(0.2~{\rm EeV},\,4\times10^{26}\,{\rm s})\) for \({\rm DM}\to \nu\bar\nu\), \((1~{\rm EeV},\,7\times10^{26}\,{\rm s})\) for \({\rm DM}\to \tau\bar\tau\), and \((100~{\rm EeV},\,4\times10^{26}\,{\rm s})\) for \({\rm DM}\to b\bar b\), but states that the preferred regions in all three channels are excluded by global gamma-ray lower bounds on the lifetime at \(95\%\) CL [2606.09986].

The event has additionally motivated proposals involving sterile-neutrino reconversion over detector-specific path lengths, speculative electroweak-vacuum-turbulence scenarios in merging black-hole binaries, and other geometry-dependent or transient new-physics mechanisms [2507.02246]. A PBR study is explicit that if KM3-230213A was “not beginner’s luck,” detector-dependent beyond-standard-model scenarios become compelling targets for future Earth-skimming neutrino follow-up [2507.02246].

## 6. Multimessenger follow-up, caveats, and continuing significance

The observational response to KM3-230213A has been overtly multimessenger. Fermi-LAT searches found no convincing cascade counterpart in GeV gamma rays [2503.16606]. VERITAS found no significant TeV emission during its 2025 follow-up campaign [2509.25422]. HAWC is cited in follow-up work as having reported no significant emission in multiple windows up to one year after the event [2509.25422]. Counterpart searches in radio, X-ray, and blazar catalogs have produced plausible associations, but none is secure [2605.07641].

Several caveats recur across the literature. The strongest interpretations are often model dependent in source distance, source duration, proton composition, or hidden-photon-field assumptions. Many analyses rely on the median event energy \(E=220\) PeV despite the broad uncertainty interval [2502.09548]. Source-specific counterpart studies are limited by the event’s angular uncertainty and by incomplete contemporaneous X-ray or gamma-ray coverage [2507.16392]. Population-level fits are intrinsically constrained by single-event statistics [2606.09986]. Even the most straightforward source comparisons are conditioned on whether the event is treated as diffuse, transient, point-like, cosmogenic, or exotic [2605.07641].

Despite these uncertainties, KM3-230213A has already altered the structure of the field. It has sharpened multimessenger tests of hadronic source models, constrained the proton fraction and source evolution of UHECR populations, strengthened neutrino-sector Lorentz-invariance tests, and motivated detector-complementarity studies for future instruments such as IceCube-Gen2, GRAND, POEMMA-Balloon with Radio, and the full KM3NeT/ARCA array [2507.11993]. A plausible implication is that the event’s long-term importance will depend less on whether any single current interpretation survives unchanged, and more on whether subsequent UHE neutrino detections reveal a population, a transient class, a cosmogenic component, or a pattern of anomalies requiring new physics.

Source: https://www.emergentmind.com/topics/km3-230213a