A Dynamical Equilibrium Linking Nanohertz Stochastic Gravitational Wave Background to Cosmic Structure Formation

Abstract: The stochastic gravitational wave background (SGWB) is conventionally treated as a passive relic of its astrophysical and cosmological sources, with negligible back-reaction on the matter content of the Universe. Here we show that this assumption needs to be modified once the SGWB and matter are treated as a dynamically coupled non-equilibrium system. Combining linearized general relativity with the fluctuation-dissipation theorem, we derive a generalized Langevin framework that drives the coupled system toward a dynamical equilibrium, which is characterized by a distinctive strain spectrum with a high-frequency cutoff $\mathcal{W}$, and a scale-dependent coupling parameter that screens gravity progressively for the most massive structures. Three findings support this framework. Fitting the equilibrium spectrum to the NANOGrav 15-year dataset yields a Bayes factor of $48\pm 3.8$ over the supermassive black hole binary baseline, achieved entirely within general relativity and the Standard Model. The PTA-calibrated screening mass scale $m_{c}\sim 10^{12}\text{--}10^{14}\,M_{\odot}$ overlaps, with no free cosmological parameter, the $Λ$CDM-derived linear-to-nonlinear transition mass $M_{\rm NL}$ of cosmic structure at $\sim 8\,h^{-1}\,\mathrm{Mpc}$. Most strikingly, promoting this concordance to a structural identification expresses $\mathcal{W}$ entirely in terms of $M_{\rm NL}$, and its inverse acquires a transparent physical reading as a coherence threshold for SGWB-matter coupling. $\mathcal{W}$ is thereby a derived quantity linking nanohertz gravitational-wave observables to the late-time cosmological sector. The framework makes distinctive scale-dependent predictions testable by forthcoming large-scale structure surveys and space-borne gravitational-wave observatories.

• The paper presents a novel DEGWB model that uses a generalized Langevin approach to balance SGWB fluctuations with gravitational radiation dissipation.
• It demonstrates, through Bayesian analysis of PTA data, that the equilibrium strain spectrum with a physical cutoff is strongly supported over standard models.
• The model predicts a mass screening (~10^12–10^14 M⊙) that influences cosmic structure formation, offering clear observational tests for future surveys.

Dynamical Equilibrium between Nanohertz Stochastic Gravitational Wave Background and Cosmic Structure Formation

Introduction and Motivation

This work presents a theoretical and observational re-examination of the role of the nanohertz SGWB in cosmic structure formation, departing from the traditional paradigm where GWs are treated as a passive relic with negligible back-reaction on the matter sector. The conventional $\Lambda$CDM model assumes that gravitational waves, due to their weak coupling to matter, exert no significant dynamical influence on the evolution of cosmic structures. However, recent pulsar timing array (PTA) results revealing nanohertz SGWB signals, especially from the NANOGrav, EPTA, PPTA, and CPTA collaborations, prompt reconsideration of this assumption.

The authors propose a novel framework in which the SGWB and cosmic matter fields constitute a dynamically coupled, non-equilibrium system. By invoking concepts from non-equilibrium statistical mechanics, particularly the generalized Langevin equation and the fluctuation-dissipation theorem, the paper constructs a model of scale-dependent energy exchange that drives the system towards a dynamical equilibrium—DEGWB—characterized by a non-trivial strain spectrum and a frequency cutoff $\mathcal{W}$. The interplay between stochastic driving by the SGWB and dissipation through gravitational radiation reaction introduces a feedback mechanism absent from the standard structure formation narrative.

Theoretical Framework: Generalized Langevin Description

The central theoretical construct is the generalized Langevin equation applied to a test mass in a stochastic, homogeneous, isotropic SGWB. The equation includes:

• A fluctuation term arising from the stochastic tidal forces due to the GW background
• A non-Markovian dissipation kernel reflecting radiative back-reaction, whose scale dependence induces "gravitational screening"

The dynamical equilibrium is achieved when the mean kinetic energy injected by the SGWB fluctuations is balanced by energy loss due to GW radiation, leading to a fluctuation-dissipation relation for the matter-GW coupled system. Importantly, the dissipative kernel shows Lorentzian frequency dependence with a high-frequency cutoff $\mathcal{W}$, encapsulating the loss of phase coherence for interactions at timescales shorter than the Schwarzschild light-crossing time of massive structures.

The spectrum of the equilibrium SGWB strain, $S_h(\omega)$, is predicted to manifest this cutoff and, for $\mathcal{W}\to\infty$, the model would be ultraviolet-divergent and unphysical. The cutoff emerges as a phenomenological consequence of the interaction timescales between the SGWB and the largest non-linear cosmic structures, not as a regulator imposed ad hoc. The effective gravitational constant becomes mass-dependent: $$G_{\mathrm{eff}}(m) = \frac{G}{1 + 4Gm\mathcal{W}/c^3}$$ resulting in progressive screening for $m \gg c^3/(4G\mathcal{W})$. This uniquely links GW observables with non-linear structure formation thresholds.

Bayesian Evidence from PTA Data

The equilibrium SGWB strain spectrum, when fitted to the NANOGrav 15-year free spectrum under Hellings-Downs correlations, is shown to yield high Bayesian evidence, with a Bayes factor $\mathcal{B} = 48 \pm 3.8$ over the canonical supermassive black hole binary background (SMBHB), on par with the best-performing SIGW models, without introducing any exotic physics beyond general relativity and the Standard Model. The amplitude and cutoff frequency parameters are tightly constrained by the data.

Figure 1: Bayes factors for model comparisons between DEGWB, DEGWB with varying cutoff prescriptions, and alternative SGWB scenarios, revealing strong statistical evidence for the proposed framework.

To explore parameter constraints, two-dimensional posterior distributions for the amplitude and cutoff frequency are constructed, exhibiting minimal degeneracy and well-localized posterior support.

Figure 2: Marginalized posterior for $\log_{10}A$ and $\log_{10}\mathcal{W}$, indicative of tight parameter constraints and minimal degeneracy in the Bayesian inference.

Fitting the equilibrium spectrum to the full free spectrum from NANOGrav, the DEGWB model shows excellent agreement across the entire PTA frequency domain.

Figure 3: Bayesian free spectrum for the NG15 dataset under Hellings-Downs correlations (grey violins) compared with model predictions at posterior mean parameters, demonstrating the empirical adequacy of the DEGWB strain spectrum.

The model demonstrates robustness with respect to details of the cutoff function (Lorentzian, exponential, power-law smoothing), with statistical evidence always favoring the existence of a physical cutoff and decisively rejecting both hard cutoffs and no-cutoff scenarios.

Physical Interpretation: Screening and the Mass Cutoff

A key result is the empirical determination of the mass scale at which SGWB-matter coupling becomes significant—the screening mass $\mathcal{W}$0–$\mathcal{W}$1—directly fixed by the PTA-inferred $\mathcal{W}$2. This mass scale overlaps with the $\mathcal{W}$3CDM-predicted nonlinear-to-linear transition mass $\mathcal{W}$4 at comoving scale $\mathcal{W}$5, required to match the present-day $\mathcal{W}$6. No additional cosmological parameters are tuned or introduced, making this identification a direct, falsifiable prediction.

The cutoff frequency $\mathcal{W}$7, therefore, is recast as a function of the growth scale for non-linear structure: $\mathcal{W}$8 Its inverse sets a phase-coherence threshold for SGWB-matter coupling, with interaction and dissipation sharply suppressed above this frequency.

The authors stress that their effective screening is not a strong, abrupt departure from standard gravitational dynamics but a mild, progressive effect, remaining within the perturbative regime compatible with the underlying Langevin framework.

Implications and Outlook

The model predicts a scale-dependent effective gravitational constant, with structure growth suppressed only for the most massive halos. This specificity—particularly the lack of uniform rescaling across all scales—distinguishes the DEGWB scenario from generic modified-gravity or interacting dark sector models, offering an observational discriminant in next-generation large-scale structure surveys.

The theoretical framework suggests SGWB plays an active, thermodynamic role in regulating the universe’s homogeneous and isotropic large-scale structure, rather than being a mere fossil record of early universe physics. Among immediate observational consequences is the potential for distinctive features in the high-mass end of the cluster halo mass function and the mass-dependent evolution of $\mathcal{W}$9 that can be scrutinized by surveys such as Euclid, LSST, and space-based GW observatories.

The overlap of the screening mass with long-standing structure formation scales also hints at a novel route to addressing, or at least interpreting, extant cosmological tensions, for instance in $\mathcal{W}$0, without resort to exotic new physics.

Conclusion

This work defines a consistent, GR-grounded framework—DEGWB—in which the SGWB establishes a dynamical equilibrium with the cosmic matter field, with testable, non-trivial consequences for both gravitational wave and cosmic structure observations. The key parameters of the model, notably the cutoff frequency $\mathcal{W}$1, emerge as derived, non-fundamental quantities linked to cosmological observables, providing a conceptual and phenomenological bridge between nanohertz GW experiments and the macrophysics of cosmic structure formation. Future multi-band GW detectors and percent-level large-scale structure surveys are expected to further test and constrain this equilibrium paradigm.