IceCube Search for MeV Neutrinos from Mergers using Gravitational Wave Catalogs

This paper reports a search by the IceCube Neutrino Observatory for MeV neutrinos associated with compact binary mergers containing at least one neutron star during the LIGO-Virgo-KAGRA O1–O3 runs, finding no significant excess and establishing stringent upper limits on neutrino emission, including for the confirmed GW170817 event.

The Gist

Imagine the universe is a giant, dark ocean. For a long time, we've only been able to see the "waves" on the surface (light) or feel the "tremors" in the water (gravitational waves). But there's another kind of signal hiding in the deep: neutrinos. These are ghostly, tiny particles that zip through everything without stopping.

This paper is about a team of scientists using a giant detector called IceCube (buried deep in the ice at the South Pole) to try and catch a specific type of these ghost particles: MeV neutrinos. These are "low-energy" neutrinos, which are like the gentle ripples compared to the massive tsunamis of high-energy neutrinos IceCube usually hunts.

Here is the story of their search, broken down simply:

1. The Big Idea: Listening for the "Pop"

When two heavy cosmic objects crash into each other—specifically when a Neutron Star (an ultra-dense dead star) hits another Neutron Star or a Black Hole—it creates a massive explosion.

• The Analogy: Think of it like two cars crashing at high speed. The crash creates a huge cloud of smoke and heat. In space, this "smoke" is a burst of hot, dense energy that spits out a flood of these ghostly neutrinos.
• The Goal: The scientists wanted to see if, right at the moment the LIGO/Virgo detectors heard the "crash" (gravitational waves), IceCube would also see a sudden "pop" of these neutrinos.

2. The Challenge: Finding a Needle in a Haystack (That's on Fire)

Detecting these specific neutrinos is incredibly hard.

• The Problem: IceCube is made of thousands of light sensors (like lightbulbs) frozen in the ice. Usually, these sensors are quiet. But they are also noisy. Radioactive glass in the sensors and cosmic rays from space make them "twitch" randomly, like a room full of people constantly dropping their phones.
• The Solution: You can't look at one sensor to find a neutrino; it's just too much noise. Instead, the scientists looked at the whole room at once.
• The Analogy: Imagine you are in a crowded stadium where everyone is whispering. If one person shouts, you might not hear them. But if everyone suddenly starts clapping at the exact same millisecond, you know something big happened. IceCube looks for that "sudden collective clap" across all its sensors.

3. The Hunt: The "On" and "Off" Switch

The scientists used a list of 83 cosmic crashes (gravitational wave events) detected between 2015 and 2020.

• The Strategy: For every crash, they looked at the IceCube data for a few seconds before and after the event.
  • The "On" Time: The exact moment of the crash.
  • The "Off" Time: The hours before and after, when no crash was happening.
• The Test: They compared the noise level during the crash to the noise level when nothing was happening. They asked: "Did the sensors get significantly louder exactly when the crash happened?"

They tried looking at four different time windows (0.5 seconds, 1.5 seconds, 4 seconds, and 10 seconds) because they weren't sure how long the "neutrino burst" would last. It's like trying to catch a firework: you don't know if it explodes instantly or burns for a few seconds, so you watch for a while.

4. The Results: The Silence

After checking all 83 events, the answer was: Nothing.

• The Finding: There was no "collective clap." The sensors didn't get any louder during the crashes than they were at random times.
• The Population Check: They also looked at the group of crashes involving Neutron Stars (the ones most likely to make neutrinos) versus Black Hole crashes (which shouldn't make them). Even when grouping them together, there was no signal.

5. Why This Matters (Even if they found nothing)

You might think, "If they found nothing, why write a paper?"

• Setting the Rules: In science, finding nothing is still a discovery. It tells us what doesn't happen.
• The New Limit: Before this, we didn't know how bright these neutrino bursts could be. Now, we know they can't be brighter than a certain limit.
• The Analogy: Imagine you are looking for a specific type of firefly in a forest. You look everywhere and don't see one. You can't say fireflies don't exist, but you can say, "If they are there, they are much dimmer than we thought, or they are much rarer than we hoped."

The Bottom Line:
The IceCube team acted like cosmic detectives, listening for a specific "ghostly whisper" from the most violent crashes in the universe. They listened to 83 crashes and heard silence. While they didn't find the neutrinos this time, they have now drawn a very tight boundary around where those neutrinos could be, helping future scientists know exactly what to look for next time.

Technical Summary

1. Problem Statement

The paper addresses the challenge of detecting low-energy (MeV-scale) neutrino bursts associated with compact binary mergers, specifically those involving neutron stars (Binary Neutron Star [BNS] and Neutron Star–Black Hole [NSBH] systems).

• Physical Context: Mergers involving neutron stars are predicted to generate hot, dense conditions that produce intense bursts of quasi-thermal neutrinos (average energy $\langle E_\nu \rangle \approx 10–30$ MeV) via electron antineutrino emission. Theoretical models predict luminosities on the order of $10^{51}–10^{53}$ erg/s.
• Detection Challenge: Unlike high-energy neutrinos, MeV neutrinos interact primarily via Inverse Beta Decay (IBD: $\bar{\nu}_e + p \to e^+ + n$) in ice. The resulting positrons produce short tracks and Cherenkov photons that often do not reach the Digital Optical Modules (DOMs) individually. Consequently, the signal manifests not as distinct tracks but as a collective, transient increase in the single-photoelectron hit rate across the entire detector array, which must be distinguished from a high background of noise (radioactive decays, atmospheric muons).
• Motivation: While Core-Collapse Supernovae (CCSNe) also emit thermal neutrinos, no nearby CCSN has coincided with gravitational wave (GW) observations in the current era. Therefore, GW-triggered mergers represent the most viable targets for multi-messenger searches for MeV neutrinos.

2. Methodology

The analysis utilizes data from the IceCube Neutrino Observatory during the LIGO-Virgo-KAGRA (LVK) O1, O2, and O3 observing runs.

• Data Source:

  • GW Catalogs: 83 GW alerts were analyzed, comprising 77 Binary Black Hole (BBH) mergers and 6 mergers containing at least one neutron star (BNS/NSBH).
  • Neutrino Data: The Supernova Data Acquisition System (SNDAQ) provided a dedicated dataset ("SN Data") containing 0.5-second binned detector rates from all 5,160 DOMs.

• Signal Detection Algorithm:

  • Test Statistic ($\xi$): The search looks for a statistically significant deviation ($\Delta\mu$) in the collective hit rate relative to the expected background mean ($\mu_i$) and standard deviation ($\sigma_i$).
  • Correction for Noise: To mitigate broadening of the test statistic distribution caused by correlated noise (e.g., atmospheric muons), a deadtime of $\approx 250$ $\mu$s is applied during acquisition. A linear correction is subsequently applied to the test statistic: $\xi_{corr} = \xi - b \cdot R_\mu - a$, where $R_\mu$ is the muon-induced hit rate.
  • Time Windows: Four distinct time windows centered on the GW trigger time were analyzed to account for uncertainties in emission duration: 0.5 s, 1.5 s, 4.0 s, and 10 s.

• Statistical Framework:

  • Individual Event Analysis: A frequentist "on-time" vs. "off-time" approach was used. The "on-time" is the window containing the GW trigger; the "off-time" spans $\pm 24$ hours (excluding the on-time) to establish a background distribution.
  • Population-Level Analysis: A binomial test was performed on the set of p-values from all events. The GW catalog was split into two populations: (i) BNS/NSBH (expected to emit neutrinos) and (ii) BBH (expected not to emit). This tests if the number of significant events exceeds background fluctuations.

3. Key Contributions

• First Constraints: This work presents the first IceCube constraints on MeV neutrino emission temporally associated with gravitational wave sources.
• Model-Independent Search: The analysis employs a model-independent approach, scanning multiple time windows to capture various transient models (from short-lived merger emissions to time-of-flight delays).
• Population-Level Sensitivity: By applying a binomial test to the entire catalog, the study assesses the potential for a cumulative signal even if individual events are below detection thresholds.
• Specific Limits on GW170817: The paper provides specific upper limits for GW170817, the first confirmed BNS merger, establishing one of the strongest limits to date for such a source.

4. Results

• Individual Events: No significant excess of neutrinos was found for any of the 83 GW events. After correcting for the trials factor (selecting the most significant window and the number of sources), all p-values exceeded 0.65.
• Population Analysis:
  • BNS/NSBH Population: The corrected p-value was 0.12.
  • BBH Population: The corrected p-value was 0.81.
  • Conclusion: There is no statistically significant correlation between MeV neutrino emission and GW events in either population.
• Upper Limits: Since no signal was detected, the collaboration set upper limits on the isotropic equivalent energy of electron antineutrinos at the source.
  • For the BNS/NSBH population, limits were derived assuming a monochromatic spectrum between 3 and 100 MeV.
  • These limits constrain the neutrino luminosity of merger events, effectively ruling out the most optimistic theoretical emission models for the observed sample.

5. Significance

• Multi-Messenger Astrophysics: This study demonstrates IceCube's capability to act as a real-time monitor for low-energy neutrino bursts, complementing high-energy neutrino searches.
• Physics Constraints: The non-detection places stringent constraints on the physics of neutron star mergers. It suggests that either the neutrino emission is less luminous than some theoretical predictions, the emission is beamed away from Earth, or the duration/intensity is insufficient to trigger a collective rate increase in IceCube given current background levels.
• Future Prospects: The methodology established here provides a robust framework for future searches as the LVK detector sensitivity improves and the number of detected mergers increases, potentially allowing for the detection of the first MeV neutrino burst from a compact binary merger.