Ten Decades of Neutrino Searches, One Combined Constraint

High-energy neutrinos are the ghosts of the cosmos. They cross galaxies, magnetic fields, and planets without stopping. That is what makes them valuable: a neutrino arriving from deep space points back to some of the most violent processes in the Universe. It is also what makes them hard to catch.

Above PeV energies—one PeV is 10^15 electron volts—the neutrino sky becomes especially important. These particles are expected to accompany ultra-high-energy cosmic rays, the most energetic protons and nuclei known. Some may be born inside or near cosmic accelerators. Others, called cosmogenic neutrinos, may be produced when cosmic rays interact with the faint radiation filling intergalactic space during their long journeys.

A patchwork of limits, not a stack of curves

The problem is that no single detector can see this whole energy range. Cherenkov neutrino telescopes such as IceCube and Baikal-GVD watch for faint light in Antarctic ice or deep water. The Pierre Auger Observatory looks for very inclined air showers made by neutrinos hitting the atmosphere. Radio experiments in and above polar ice listen for coherent radio pulses from particle cascades. Lunar searches use the Moon itself as a target, hunting for radio flashes from neutrino interactions in the lunar regolith.

Each technique has produced upper limits. Plotted together, they look like a patchwork of curves stretching from 10^16 to 10^26 eV. But those curves cannot simply be stacked. Some quote limits for one neutrino flavor; others for all flavors. Some use different energy bins. Some use different statistical conventions to turn event counts into limits. A limit is not the same as a detection; it is a statement about what a detector could have seen if a certain flux were present.

A Poisson combination from 10^16 to 10^26 eV

That is the gap a new preprint by S. V. Troitsky addresses. Rather than combining the published limit curves, the paper goes back to what each experiment could actually see. For each search, Troitsky reconstructs the energy-dependent exposure: a measure of how much detector area, observing time, and solid angle were available to a neutrino flux at each energy. The exposures are converted to a common convention for the total all-flavor neutrino plus antineutrino flux, assuming an equal mix of electron, muon, and tau neutrinos at Earth.

Then the published event counts and expected backgrounds are combined channel by channel in a Poisson likelihood. This keeps track of which experiment saw each surviving event, instead of collapsing everything into a single number. The result is a set of 90% confidence upper limits for one-decade-wide test spectra, spanning 10^16 eV to 10^26 eV.

The strongest bound appears in the lowest energy bin, near 3 x 10^16 eV. There, the 90% upper limit on the flux, written in the common scaled form E^2 dΦ/dE, is 6.1 x 10^-9 GeV cm^-2 s^-1 sr^-1. At higher energies the bound is generally weaker; the paper's table shows the one-decade limits increasing overall from the lowest bin.

The combined exposure also shows how the search methods fit together. Optical Cherenkov telescopes dominate the low-energy end. Auger and radio detectors become important around EeV energies, or 10^18 eV. ANITA controls part of the intermediate radio range. Lunar observations extend the reach to the highest energies considered.

Which models sit above the combined ceiling

The paper then uses the same reconstructed exposures to test theoretical spectra. Instead of comparing a model curve by eye with a plot of limits, it folds each model through the detector exposures and asks how much the model's normalization would have to be scaled down to satisfy the combined data. A scale factor below 1 means the model's nominal flux is already above the 90% upper endpoint of this count-level test.

Several cosmogenic-neutrino benchmarks remain below the combined ceiling, including the broad envelope associated with Pierre Auger fits and a joint-fit two-population model. But two benchmarks sit above it: the Telescope Array best-fit model and a KM3NeT-only fit to the ultra-high-energy event KM3-230213A. The paper treats the Alhebsi et al. benchmarks—including that KM3NeT-only fit—as compatibility diagnostics because those spectra were inferred using overlapping neutrino data. Three cosmic-string models also remain well below the present sensitivity.

Not every exotic scenario can be tested this way. Superheavy dark matter would produce an anisotropic neutrino signal, because the Milky Way dark-matter halo is not uniform, rather than the isotropic glow assumed here. A meaningful search would need direction-dependent exposures from each experiment, which are not uniformly public. The paper therefore does not quote a numerical constraint on those models.

One observational wrinkle is the set of four near-horizon air-shower events seen by the ANITA-IV balloon experiment. They are omitted from the main combination because the ANITA-IV analysis did not define an internal diffuse-background expectation for them. If all four were treated as possible diffuse neutrino candidates, the combined limits weaken by as much as a factor of 2.5. If they are treated as background, the change in the primary one-decade bounds is at most about 11%. That spread is a reminder that a handful of ambiguous events can shape a limit in a statistically limited regime.

The analysis is deliberately conservative and not a replacement for the dedicated likelihoods built by each collaboration. It uses public event counts and backgrounds, not reconstructed energies, directions, or detector-specific nuisance parameters. Some exposures had to be reconstructed from published flux limits, and the results are quoted to two significant figures because experimental systematics dominate. The paper is a preprint and has not yet completed peer review.

Still, its value is clear. It turns a heterogeneous collection of searches into a reproducible benchmark across ten decades of neutrino energy. Projected sensitivities—including IceCube-Gen2's radio array, the full Askaryan Radio Array, PUEO, GRAND, Trinity, and lunar radio concepts—are expected to push deeper; the paper notes that ARA and PUEO have already collected relevant data but have not yet published final diffuse-flux limits. When those limits come, they will have a common scale against which to measure their progress. The ultra-high-energy neutrino sky is still mostly quiet, but it is becoming a better-bounded, more comparable kind of quiet.