Irreducible Gravitational Wave Background as a Particle Detector

Irreducible Gravitational Wave Background as a Particle Detector

Motivation and Conceptual Framework

The stochastic gravitational wave background (GWB) offers a unique, direct probe of the early Universe, accessing pre-BBN cosmological epochs that are not observable through electromagnetic or other standard messengers. The sensitivity of upcoming gravitational wave observatories will allow scrutiny of GW spectra across more than twelve decades in frequency. Among potential GW sources, transient processes—including primordial tensor perturbations from inflation, first-order phase transitions, and decaying domain walls—are particularly salient. These sources generate a GWB that can propagate freely once production ceases.

A central finding of this work is that a range of beyond-the-Standard-Model (BSM) frameworks generically include long-lived heavy particles (LLPs), which, if sufficiently long-lived, induce a finite epoch of early matter domination (EMD) in the Universe’s thermal history. This period, even if short, alters the cosmic equation of state and leaves indelible signatures in the gravitational wave spectrum, manifesting as two parametric spectral features that correspond to the onset and termination of EMD. These spectral breaks serve as direct encoders of the mass, decay rate, and initial abundance of the LLP responsible for the nonstandard expansion history.

The key logical structure of the inference from GW observables to microscopic BSM parameters is illustrated as follows: measurements of these two characteristic frequencies fully determine, up to model uncertainties in initial abundance, both the mass and decay width of the LLP, and hence the Lagrangian coupling governing the new particle’s decay.

(Figure 1)

Figure 1: Schematic of the inference chain from gravitational wave observables to Lagrangian parameters; spectral features $f_1$ and $f_2$ encode the EMD period's end and onset respectively, and thus the underlying microphysics.

Early Matter Domination and Gravitational Wave Propagation

Considering a generic fermionic or bosonic LLP—decoupled from the thermal bath at high temperature with initial yield $Y_i$, mass $M$, and decay width $\Gamma$—one finds that the onset ($f_2$) and termination ($f_1$) of EMD imprint corresponding breaks in the GW spectrum. The temperature at which EMD begins is proportional to $Y_i M$, whereas its end is fixed by $\Gamma$ via the equality $\Gamma = H$ (Hubble parameter). Accurate mapping of these relations requires solving the Boltzmann system for LLP and radiation densities and tracking $f_2$0, the evolving equation of state.

The GW modes themselves are sensitive to changes in the background only after they enter the horizon. This leads to a suppression of the GW energy density at frequencies that re-entered the horizon during EMD, resulting in a broken power-law feature: modes above $f_2$1 are suppressed, while those below $f_2$2 remain unaffected, and the intermediate region marks the transition period.

Analytical and Numerical Mapping—Spectral Features and Microphysical Parameters

The analysis establishes that the transfer function encoding EMD imprint on the GWB, $f_2$3, can be fitted accurately by the functional form

$f_2$4

where $f_2$5 and $f_2$6 are, respectively, the characteristic frequencies associated with the end and onset of EMD. Solving for a wide range of $f_2$7, the following best-fit expressions are obtained:

• $f_2$8
• $f_2$9

  Figure 2: Benchmark GW backgrounds from inflation (red), first-order phase transitions (orange), and domain wall collapse (purple), showing the impact (solid) and absence (dashed) of an EMD epoch. The frequency positions of spectral breaks map directly to particle decay lengths accessible to colliders.

  

  Figure 3: Numerical scan showing the dependence of the feature frequencies $Y_i$0 (left: decay rate dependence only) and $Y_i$1 (right: strong initial yield and mass dependence, mild decay rate dependence); model benchmarks indicated.

These fits are robust for regimes with pronounced EMD, breaking down smoothly as the LLP decay rate approaches the threshold where EMD is minimized. Notably, $Y_i$2 is a clean probe of the decay rate, with $Y_i$3 encoding both the initial abundance and mass.

By recasting $Y_i$4 in terms of laboratory-observable decay lengths $Y_i$5, one finds that frequencies probed by current and next-generation pulsar timing arrays, such as NANOGrav, SKA, and Theia, are directly sensitive to the same decay lengths targeted by laboratory LLP searches (FASER, DUNE, SHiP).

Gravitational Wave Mode Evolution and Validity of the Horizon-Entry Approximation

An explicit check—comparing the instantaneous horizon-entry approximation to the full tensor evolution in both constant and smoothly changing equation-of-state backgrounds—demonstrates that this approximation tightly reproduces the physics away from the immediate neighborhood of horizon crossing. It accurately captures the suppression of GW modes during EMD, ensuring the robustness of the mapping between the observed spectral breaks and the underlying cosmological transitions.

Figure 4: Comparison of tensor-mode evolution in radiation and matter domination; the approximation reproduces both super-horizon freezing and sub-horizon amplitude decay.

Figure 5: Full tensor-mode solution and instantaneous entry approximation for a smooth matter-to-radiation transition. The equation of state $Y_i$6 and crossover features are tracked.

Detector Prospects and Complementarity with Laboratory Searches

The presence of EMD-induced features in the stochastic GWB places strong constraints on or allows the direct inference of BSM parameter space. The frequency range of $Y_i$7 aligns with the decay lengths of LLPs that current (LHC-based) and upcoming dedicated experiments (such as FASER and MATHUSLA) are seeking. Thus, non-observation, or characterization, of these GW spectral features in PTA data directly constrains models of LLPs, often probing mass and coupling regimes that are far beyond the reach of terrestrial experiments.

Figure 6: Detector reach and frequency overlap of upcoming GW observatories for $Y_i$8 backgrounds and corresponding decay length mapping; direct test of parameter space inaccessible at colliders.

Additionally, scans over the allowed decay rates confirm the precise regime of validity for the analytic fits and highlight the necessity of a prolonged EMD epoch (two orders of magnitude below the onset threshold) for clean spectral feature extraction.

Figure 7: Sensitivity of the lower spectral feature ($Y_i$9) to proximity to the maximal decay rate allowing EMD; only for $M$0 is the analytic scaling precisely applicable.

Implications and Outlook

The paper demonstrates that GW observations—specifically the detection of broken power-law features in the stochastic background—function as cosmic detectors with direct access to BSM properties traditionally considered to be intractable in collider or direct/indirect detection searches. The method is agnostic to the details of the initial transient GW source, relying only on the subsequent universal evolution during a departure from standard radiation domination.

Such features, if observed, would enable the reconstruction of fundamental Lagrangian parameters from cosmological data. This provides an unprecedented complementarity between collider-based LLP programs and GW cosmology, with GW data accessing much higher mass scales and weaker couplings.

Conclusion

This study establishes that the cosmic gravitational wave background, when analyzed for specific spectral features reflecting early matter domination induced by long-lived heavy particles, provides an irreducible and robust probe of BSM physics. The two-parameter mapping between GWB observables and particle properties is numerically validated and applicable across all transient GW production mechanisms. These results motivate interdisciplinary efforts connecting gravitational wave astrophysics, particle collider phenomenology, and early Universe cosmology, confirming the role of gravitational wave observatories as de facto particle physics detectors. Future observations have the potential to exclude or decisively measure large portions of parameter space for non-thermal cosmological histories and exotic particle sectors.

Reference: "Irreducible Gravitational Wave Background as a Particle Detector" (2604.20792)