🔭 astrophysics

Tomographic constraints on the high-energy cosmic neutrino emission rate

This study utilizes cross-correlations between the IceCube 10-year neutrino dataset and tomographic galaxy maps to place upper bounds on the bias-weighted high-energy neutrino emissivity from extragalactic sources up to redshift $z\sim3$ , finding no significant correlation but establishing a new method for constraining source models.

Original authors: Alberto Gálvez Ureña, Federico Urban, David Alonso

Published 2026-02-04

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Alberto Gálvez Ureña, Federico Urban, David Alonso

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: Hunting for the Source of Cosmic "Ghost Particles"

Imagine the universe is filled with a constant, invisible rain of tiny, ghostly particles called neutrinos. These particles are so light and weak that they pass through planets and stars without stopping. A giant detector buried deep in the ice at the South Pole, called IceCube, has been catching these particles for years.

Scientists know these high-energy neutrinos come from outside our galaxy (extragalactic), but they don't know exactly where. It's like hearing a constant hum of traffic from a distant city but not knowing which specific streets the cars are coming from.

The Detective Work: Looking for Patterns in the Stars

This paper asks a simple question: Do these neutrinos come from the same places where the galaxies are clustered?

Think of the universe like a giant, 3D web. Galaxies aren't scattered randomly; they hang out in clumps and filaments, like cities on a map. If the neutrinos are coming from powerful cosmic engines (like black holes or exploding stars) that live inside these galaxies, then the neutrinos should also be clustered in those same spots.

The authors tried to find this connection by doing a "cross-correlation."

The Analogy: Imagine you have a map of where the rain is falling (the neutrinos) and a separate map of where the cities are (the galaxies). If the rain is falling mostly over the cities, you have a match. If the rain is falling randomly everywhere, there is no connection.

How They Did It: The "Tomographic" Slice

To make this comparison, the scientists didn't just look at the whole sky at once. They used a technique called tomography.

The Analogy: Think of a loaf of bread. Instead of looking at the whole loaf, you slice it. You look at the bottom slice (nearby galaxies), the middle slice (medium-distance galaxies), and the top slice (very far-away galaxies).
They used four different catalogs of galaxies, each covering a different "slice" of the universe's history (from very close to us to very far away, up to 3 billion light-years back).
They compared the IceCube neutrino map against each of these galaxy slices to see if the neutrinos were "hanging out" with the galaxies in any specific layer.

The Models: Guessing the Recipe

Since they didn't know exactly how the neutrino sources behave, they tested two main "recipes" (models) for how the neutrino production might change over time:

The Power Law: A simple rule where the number of neutrinos changes steadily as you go back in time (like a volume knob turning up or down).
The Star Formation History: A model that assumes neutrinos are made in proportion to how many stars are being born. Since we know star formation peaked in the past and is slowing down now, this model predicts a specific curve of neutrino production.

They also tried a "non-parametric" approach, which is like asking, "How many neutrinos are coming from each specific slice of the bread?" without assuming a specific rule for the whole loaf.

The Results: The Silence of the Neutrinos

After crunching the numbers, the result was... nothing significant.

The Finding: They found no strong evidence that the neutrinos are clustering with the galaxies. The maps looked like random noise.
The "Signal": The strongest hint of a connection was a "1.9 sigma" deviation. In the world of science, this is like hearing a faint whisper in a noisy room. It's interesting, but it's not loud enough to say, "I definitely heard something!" (Scientists usually need a "5 sigma" result to claim a discovery).
The Conclusion: The neutrinos might be coming from places we haven't mapped yet, or they might be coming from everywhere so evenly that we can't see a pattern.

What They Did Learn: Setting the Limits

Even though they didn't find a "smoking gun," they learned something valuable: Upper Limits.

The Analogy: Imagine you are trying to guess how much sugar is in a mystery cake. You can't taste it, but you can say, "It definitely doesn't have more than 500 grams of sugar." You haven't found the exact amount, but you've ruled out the possibility that it's a giant sugar bomb.
The authors calculated the maximum possible amount of neutrinos that could be coming from these galaxy clusters without us noticing a pattern. They set strict "speed limits" on how bright these cosmic neutrino sources can be.

Summary

In short, the scientists took a 10-year dataset of cosmic neutrinos and compared it to a detailed map of galaxies across different eras of the universe. They used clever math to look for a pattern, like matching rain to cities. They didn't find a match. However, by proving that the neutrinos aren't clustering heavily with the galaxies, they have placed strict limits on how powerful these mysterious cosmic sources can be. This helps narrow down the list of suspects for what is creating these high-energy particles.

Technical Summary: Tomographic Constraints on the High-Energy Cosmic Neutrino Emission Rate

Problem Statement
Despite the detection of a diffuse flux of high-energy astrophysical neutrinos by the IceCube Observatory, the specific origins of the majority of these events remain unidentified. While specific sources have been associated with individual neutrino events, they account for only a minor fraction of the total flux. The dominant population appears isotropic, leaving the characteristics of the source populations undetermined. Identifying these sources is a primary challenge in contemporary astrophysics. A promising approach involves searching for statistical anisotropies in the arrival directions of neutrinos, hypothesizing that they trace the large-scale structure (LSS) of the Universe. However, previous attempts to correlate neutrino maps with galaxy surveys have yielded mixed or non-significant results, often lacking the ability to constrain the absolute physical parameters of the emission rate.

Methodology
This work employs a cross-correlation analysis between the IceCube 10-year muon-track dataset and tomographic maps of galaxy overdensity to constrain the bias-weighted high-energy neutrino emissivity. The methodology proceeds as follows:

Data Construction:
- Neutrino Data: Utilized 10 years of IceCube muon-track data (403,529 events) in the energy range $10^3$ to $10^6$ GeV. Maps were constructed for both number density ( $I_N$ ) and energy density ( $I_E$ ). To mitigate atmospheric neutrino contamination, events with declinations $\delta < -5^\circ$ were excluded, and monopole and dipole components were subtracted from the maps.
- Galaxy Data: Four distinct galaxy and quasar catalogues were used to cover a broad redshift range ( $z \lesssim 3$ ): 2MPZ ( $z \sim 0.06$ ), WI $\times$ SC ( $z \sim 0.23$ ), DECaLS ( $z \sim 0.50$ ), and Quaia ( $z \sim 1.72$ ). These were processed into HEALPix maps ( $N_{side}=256$ ) with masks applied to remove galactic plane contamination.
Theoretical Framework:
- The analysis defines neutrino intensity in terms of comoving neutrino emissivity ( $\dot{n}_\nu$ ) and comoving neutrino energy emissivity ( $\dot{\rho}_\nu$ ), weighted by the linear bias of the sources ( $b_\nu$ ).
- Angular cross-power spectra ( $C_\ell^{NG}$ ) between the neutrino maps and galaxy overdensities are calculated using the Limber approximation. The theoretical spectra depend on the product of the galaxy bias ( $b_g$ ), the neutrino bias ( $b_\nu$ ), the mean emissivity, and the matter power spectrum.
- Galaxy biases ( $b_g$ ) were determined independently via auto-correlation analysis and treated as fixed parameters in the cross-correlation likelihood.
Statistical Analysis:
- A Bayesian inference approach using nested sampling (UltraNest) was employed to constrain model parameters.
- Three models were tested for the redshift evolution of the bias-weighted emission rate ( $b_\nu \dot{n}_\nu$ $b_{ν} \overset{n}{˙}_{ν}$ and $b_\nu \dot{\rho}_\nu$ $b_{ν} \overset{ρ}{˙}_{ν}$ ):
  - Power Law: $\propto (1+z)^a$ .
  - Peak Model: Tracking the star formation history (Madau & Dickinson 2014).
  - Tomographic (Non-parametric): Constraining the emission rate independently in the mean redshift bins of the four galaxy catalogues.
- Signal-to-noise ratios (SNR) were calculated to assess detection significance.

Key Results

Correlation Significance: The analysis found no significant correlation between the IceCube neutrino maps and the large-scale structure tracers. The strongest deviation from a zero-signal model was a $1.9\sigma$ detection (specifically for the WI $\times$ SC sample in the energy-weighted map). All other cross-correlations were consistent with zero within statistical uncertainties.
Parameter Constraints:
- Parametric Models: For the power-law model, the amplitude parameter $b_\nu \dot{n}_\nu$ was constrained to $5.8 \pm 3.6 \times 10^{45} \text{ Mpc}^{-3} \text{ yr}^{-1}$ (compatible with zero at $1.6\sigma$ ). For the peak model, the value was $4.0 \pm 2.1 \times 10^{45} \text{ Mpc}^{-3} \text{ yr}^{-1}$ (compatible with zero at $1.9\sigma$ ). The power-law index $a$ remained unconstrained within the prior range.
- Tomographic Constraints: In a model-independent approach, upper bounds were placed on the emission rate in four redshift bins. While the central values were positive, they remained consistent with zero within $2\sigma$ across all bins (e.g., $4.7 \pm 3.5 \times 10^{45} \text{ Mpc}^{-3} \text{ yr}^{-1}$ at $z \approx 0.06$ and $310 \pm 250 \times 10^{45} \text{ Mpc}^{-3} \text{ yr}^{-1}$ at $z \approx 1.7$ ).
Robustness Tests: The results were verified against variations in energy bins, declination cuts ( $\delta \ge 0^\circ$ and $\delta \ge 5^\circ$ ), and the subtraction of an atmospheric neutrino template. These tests confirmed that the lack of significant detection is robust and not an artifact of atmospheric contamination or analysis choices.

Significance and Claims
The authors claim this work provides the first tomographic constraints on the emissivity of the diffuse astrophysical background of high-energy neutrinos. Unlike previous studies that focused on relative anisotropies or normalized maps, this analysis preserves the absolute flux scale, allowing for the direct derivation of upper bounds on the physical comoving neutrino emissivity ( $\dot{n}_\nu$ ) and energy emissivity ( $\dot{\rho}_\nu$ ), albeit degenerate with the unknown source bias $b_\nu$ .

The paper positions this analysis as a new probe for testing extragalactic neutrino source models. While the current data does not yield a definitive detection, the methodology establishes a framework to place quantitative upper limits on emission rates as a function of redshift. The authors note that future datasets from next-generation neutrino observatories (IceCube-Gen2, KM3NeT) and deep galaxy surveys (LSST, Euclid) will be required to increase the constraining power and potentially achieve a detection of this correlation.