⚛️ phenomenology

High Energy Neutrinos, Gravitational Waves, and Dark Matter from a Cosmological First Order Phase Transition

This paper proposes a multimessenger framework where superheavy dark matter produced during a cosmological first-order phase transition via a filtered mechanism simultaneously explains observed high-energy neutrino events and generates a stochastic gravitational wave background, with specific mass ranges linking the signals to current or future observatories like IceCube, Cosmic Explorer, and KM3NeT.

Original authors: James M. Cline, Savas Stoica, Yong Xu

Published 2026-10-05

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: James M. Cline, Savas Stoica, Yong Xu

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 universe is filled with invisible particles that rarely interact with anything, known as dark matter. For decades, scientists have tried to figure out what these particles are made of and how they came to exist. One leading idea is that dark matter was created in the very early universe during a dramatic event called a phase transition. This is similar to how water freezes into ice, but happening to the fundamental fields of the universe itself. When such a transition occurs, it creates bubbles of a new state of matter that expand and collide, releasing energy and creating ripples in space-time called gravitational waves. At the same time, other theories suggest that heavy, unstable dark matter particles could decay over billions of years, breaking apart into high-energy neutrinos. These neutrinos are ghostly particles that travel across the cosmos and can be detected by massive observatories on Earth. The big question is whether these two phenomena—the ripples in space-time and the high-energy particles—could be linked to the same cosmic event.

A team of researchers has explored this possibility by building a detailed model of how dark matter could be created during such a phase transition. They focused on a specific mechanism where dark matter particles are filtered out as the universe changes state. Imagine a wall moving through a hot soup of particles. As the wall passes, the particles on the other side suddenly become much heavier. Only the fastest, most energetic particles can push through this wall; the rest are left behind. This process naturally limits how much dark matter exists today, matching the amount astronomers observe. The researchers used this idea to calculate exactly how fast the wall of the phase transition must have moved to produce the right amount of dark matter. By doing this, they connected the speed of the wall to the amount of dark matter and the energy of the particles it might produce later.

The study found that if the dark matter particles are extremely heavy, with a mass around 440 million times that of a proton, they could explain a specific high-energy neutrino event recently detected by the KM3NeT observatory. This event, known as KM3-230213A, involved a neutrino with an energy far higher than what is typically seen. The researchers showed that dark matter with this specific mass, created during a phase transition, would decay and produce neutrinos with just the right energy to match this observation. However, this scenario also predicts a background of gravitational waves. For this heavy mass, the waves would oscillate at a frequency of about 10,000 hertz. This is far too high for current or planned gravitational wave detectors to hear, meaning that while the neutrino signal fits, the gravitational wave signal remains hidden for now.

The researchers then looked at what would happen if the dark matter were lighter. They found that as the mass of the dark matter decreases, the energy of the resulting neutrinos drops, moving into the range that the IceCube observatory in Antarctica is designed to detect. Crucially, a lighter dark matter mass also means the phase transition happened at a lower temperature in the early universe. This shift moves the frequency of the gravitational waves down to a range where future detectors, such as the Cosmic Explorer, the Big Bang Observer, and DECIGO, could potentially hear them. This creates a powerful link: if we see a specific pattern of neutrinos from IceCube and simultaneously detect a matching gravitational wave signal from these future experiments, it would strongly suggest that both signals came from the same phase transition that created the dark matter.

The paper emphasizes that this connection relies on precise calculations of how the bubble walls of the phase transition move. Previous studies often guessed this speed, but this team used advanced hydrodynamic models to determine it based on the physics of the transition itself. This makes their prediction more reliable. They also noted that for their model to work, the lightest known neutrino must be essentially massless, a condition that could be tested by other experiments searching for the absolute mass of neutrinos. While the specific scenario for the heavy dark matter event is currently untestable via gravitational waves, the framework offers a clear path forward. It suggests that by combining data from neutrino telescopes and gravitational wave detectors, scientists could eventually confirm the origin of dark matter and witness the echoes of the universe's earliest moments.

Technical Summary: High Energy Neutrinos, Gravitational Waves, and Dark Matter from a Cosmological First Order Phase Transition

Problem Statement
The paper investigates the potential for a cosmological first-order phase transition (FOPT) to simultaneously generate three distinct observables: the relic abundance of super-heavy dark matter (DM), a high-energy neutrino (HE $\nu$ ) flux from the decay of that DM, and a stochastic gravitational wave (GW) background. While previous studies have explored connections between GWs and HE $\nu$ s (e.g., via cosmic strings or primordial black holes), a direct link where the same phase transition event determines both the DM production mechanism and the GW generation has been less explored. Specifically, the authors address the challenge of producing super-heavy DM (masses $\gg$ the thermal scale) while reproducing the observed relic abundance, and the difficulty of predicting the bubble wall velocity ( $v_w$ ), a critical parameter that governs both the DM filtering efficiency and the GW spectrum.

Methodology
The authors adopt the "filtered DM" mechanism, where a right-handed neutrino (RHN) serves as the DM candidate. This particle is thermalized in the symmetric phase of a scalar field $S$ but acquires a large mass in the broken phase. As bubbles of the true vacuum expand, only a fraction of the DM particles possess sufficient energy to penetrate the bubble wall, effectively filtering the population and suppressing the relic density.

Key methodological advancements include:

Hydrodynamic Wall Velocity: Unlike previous filtered DM studies that treated $v_w$ as a free parameter, this work calculates $v_w$ self-consistently using recent hydrodynamic results (Ref. [22]). This involves solving for the wall velocity based on thermodynamic quantities: the phase transition strength ( $\alpha$ ), the inverse duration ( $\beta/H_n$ ), and the enthalpy ratio ( $\Psi_n$ ) across the wall.
Scale Invariance and Parameter Scan: The authors exploit the scale invariance of the bounce action to perform a numerical scan over dimensionless parameters ( $\lambda_s, D, \epsilon, y_s$ ). They fix the DM mass ( $m_{DM}$ ) to match specific neutrino event energies (e.g., KM3-230213A) and the observed relic abundance to determine the required phase transition parameters.
Multimessenger Calculation:
- Neutrinos: They compute the differential neutrino flux from DM decay ( $N \to \nu h$ ), accounting for both the Milky Way and extragalactic components, including secondary neutrinos from Higgs decay showers.
- Gravitational Waves: They calculate the GW spectrum arising from sound waves, magnetohydrodynamic turbulence, and bubble collisions, utilizing the hydrodynamic inputs derived from the phase transition dynamics.

Key Contributions and Results
The paper presents a coherent framework linking DM production, HE $\nu$ s, and GWs through a single FOPT. The primary results are:

Benchmark for KM3-230213A: For a DM mass of $m_{DM} \simeq 4.4 \times 10^8$ GeV (motivated by the KM3-230213A event), the authors identify parameter regions where the filtered mechanism reproduces the observed relic abundance. The resulting neutrino flux, assuming a DM lifetime $\tau_{DM} \gtrsim 5 \times 10^{29}$ s, is consistent with the KM3-230213A event and IceCube high-energy starting event (HESE) data.
GW Signal for Heavy DM: The corresponding phase transition for this heavy DM mass produces a GW spectrum peaking at frequencies around $10^4$ Hz. However, the amplitude ( $\Omega_{GW}h^2 \sim 10^{-14}$ ) is below the projected sensitivities of current and planned detectors (LISA, BBO, DECIGO, Cosmic Explorer).
Lower Mass Regime and Observability: The authors demonstrate that for lighter DM masses ( $m_{DM} \lesssim 10^6$ GeV), the neutrino spectrum shifts into the energy range relevant for IceCube observations. Crucially, the associated GW spectrum shifts to lower frequencies, entering the projected sensitivity bands of future detectors such as Cosmic Explorer (CE), Big Bang Observer (BBO), and DECIGO.
Self-Consistency: By determining $v_w$ from hydrodynamics rather than assuming it, the paper establishes a constrained relationship between the predicted DM abundance and the GW signal. The study shows that for the benchmark points, the sound wave contribution dominates the GW spectrum, while bubble collisions are negligible due to subluminal wall velocities.

Significance and Claims
The paper claims to provide a "multimessenger connection" between high-energy neutrinos and gravitational waves, where both signals originate from the production and subsequent decay of dark matter generated during a first-order phase transition.

The authors emphasize that their model offers a concrete test of the phase transition origin of super-heavy DM. They note that the model is directly related to the seesaw mechanism for light neutrino mass generation, as the DM candidate is a heavy sterile neutrino. The longevity of the DM candidate depends on the lightest neutrino mass being essentially zero; thus, future constraints on absolute neutrino masses from cosmology or neutrinoless double beta decay could potentially falsify the model.

The work is modest in its claims regarding immediate detection: for the specific benchmark matching the KM3-230213A event, the GW signal is currently unobservable. However, the framework suggests that a correlated observation of neutrinos and GWs in the lower mass regime ( $m_{DM} \lesssim 10^6$ GeV) would strongly support this specific production mechanism, distinguishing it from scenarios where GWs and DM production are decoupled.

More like this