Weinberg Gravitational-Wave Damping Formula

• Weinberg gravitational-wave damping formula is a framework describing how tensor perturbations are attenuated by shear viscosity and anisotropic stress from free-streaming neutrinos.
• Studies show that viscous media produce exponential damping via factors like 16πGη, while collisionless neutrino damping yields suppression factors near 0.645 and 0.91032.
• Extensions of the framework include modifications from spatial curvature, dark matter interactions, and altered gravity theories, offering insights into early-universe and induced gravitational wave phenomena.

Searching arXiv for papers directly relevant to Weinberg gravitational-wave damping, including viscous-fluid damping, free-streaming neutrino damping, and later generalizations. Weinberg gravitational-wave damping formula denotes a class of results in relativistic cosmology in which gravitational-wave tensor perturbations do not propagate as freely oscillating modes once the cosmic medium develops dissipative or anisotropic stress. In the classic viscous-fluid form, shear viscosity adds a friction term to the tensor propagation equation and exponentially attenuates the wave amplitude with damping rate proportional to $G\eta$, while bulk viscosity does not contribute at leading order to transverse-traceless modes (Lu et al., 2018). In the early-universe free-streaming form associated with Weinberg’s later treatment of collisionless neutrinos, the damping is nonlocal in time: anisotropic stress sourced by decoupled relativistic particles feeds back through an integro-differential kernel and suppresses tensor amplitudes after horizon entry (Stefanek et al., 2012). Subsequent work has extended these results to viscous cosmologies, self-interacting dark matter, generalized free-streaming species, modified gravity, closed spatial geometry, and dissipative corrections relevant to induced gravitational waves (Lu et al., 2018, Dent et al., 2013, Scomparin et al., 2019, Khodagholizadeh et al., 2014, Domènech et al., 17 Mar 2025).

1. Historical definition and scope

The expression “Weinberg gravitational-wave damping formula” refers to two closely related strands of analysis. One strand concerns tensor perturbations in a viscous medium, where shear viscosity enters the stress-energy tensor and produces friction-like attenuation of gravitational waves during propagation. In the formulation used for cosmological applications, the damping rate is written as $\beta = 16\pi G\eta$, so that the strain amplitude acquires an exponential suppression $e^{-\beta D/2}$ over propagation distance $D$ (Lu et al., 2018). A later paper discussing viscous damping in a local Minkowski treatment uses $\beta = 8\pi G\eta/c^2$ in the wave equation $$\partial_t^2 h_{ij} + 2\beta \partial_t h_{ij} - \partial_z^2 h_{ij}=0$$, with the abstract stating a “damping rate of $16\pi G\eta$,” corresponding to intensity rather than amplitude damping (Ning et al., 2020).

The second strand concerns free-streaming relativistic species, especially neutrinos. In that setting the effective dissipation is not ordinary hydrodynamic viscosity but anisotropic stress generated by the collisionless Boltzmann hierarchy. Weinberg’s tensor-mode equation becomes an integro-differential relation whose kernel is built from spherical Bessel functions, and the observable damping depends on a reduced wave number $Q$ measuring horizon entry relative to matter-radiation equality (Stefanek et al., 2012). In the short-wavelength regime, the power suppression approaches $0.645$, while in the long-wavelength regime relevant to the tensor contribution to the CMB the suppression of the derivative-based observable is $0.91032$ (Stefanek et al., 2012).

These two strands are conceptually unified by the same principle: tensor modes couple to non-perfect-fluid structure in the cosmic medium. In a hydrodynamic description, the relevant quantity is shear viscosity; in a collisionless description, it is anisotropic stress carried by free-streaming species. A plausible implication is that “Weinberg damping” is best understood as a general tensor attenuation framework rather than a single formula.

2. Viscous-medium formulation

In the viscous-fluid description, the background metric is taken to be FRW with transverse-traceless perturbations,

$\beta = 16\pi G\eta$0

and the imperfect-fluid stress-energy tensor includes shear viscosity $\beta = 16\pi G\eta$1 and bulk viscosity. The cosmological application in a viscous universe states that “when a nonzero shear viscosity $\beta = 16\pi G\eta$2 is introduced to the fluid energy-momentum tensor, GWs would be dissipated by matter with a damping rate $\beta = 16\pi G\eta$3” (Lu et al., 2018). The same work emphasizes that the damping depends on shear viscosity rather than bulk viscosity or GW frequency (Lu et al., 2018).

The resulting strain formula is written as

$\beta = 16\pi G\eta$4

so the physically relevant attenuation law is

$\beta = 16\pi G\eta$5

This is the viscous form of the Weinberg damping formula (Lu et al., 2018). The interpretation is direct: the usual $\beta = 16\pi G\eta$6 geometric dilution is supplemented by exponential suppression from shear viscosity accumulated along the propagation path.

A complementary derivation in Minkowski background plus a viscous medium writes the tensor equation as

$\beta = 16\pi G\eta$7

with $\beta = 16\pi G\eta$8 in SI units (Ning et al., 2020). Its approximate solution,

$\beta = 16\pi G\eta$9

shows explicitly that amplitude damping is exponential, while intensity damping carries the factor $e^{-\beta D/2}$0 (Ning et al., 2020). This difference of a factor of two is purely a matter of whether one quotes amplitude damping or energy-flux damping.

The same viscous framework underlies later kinetic-theory work on damping by matter more generally. A unified Boltzmann treatment distinguishes collisional damping and Landau damping, showing that collisions lead to a shear viscosity but also erase anisotropic stress, thereby suppressing damping in most astrophysical settings; damping by intergalactic or interstellar matter is negligible for all but primordial radiation (Baym et al., 2017). In that sense, the viscous Weinberg formula is quantitatively significant mainly for cosmological or extreme-medium applications.

3. Effective-distance representation and observational use

A notable reformulation of the viscous result introduces an effective distance

$e^{-\beta D/2}$1

which rewrites the damped strain in the same formal form as the perfect-fluid case (Lu et al., 2018). The significance of this substitution is operational: luminosity-distance posteriors inferred from standard LIGO/Virgo waveform analyses can be reinterpreted directly as constraints on $e^{-\beta D/2}$2, rather than on the true luminosity distance $e^{-\beta D/2}$3 (Lu et al., 2018).

Using LIGO/Virgo events together with GW170817, the viscous-universe analysis found no evidence for damping and obtained an upper bound

$e^{-\beta D/2}$4

at $e^{-\beta D/2}$5 confidence level (Lu et al., 2018). The same paper notes that a source at luminosity distance $e^{-\beta D/2}$6 would be needed to probe dark-matter self-interactions through this damping channel under the assumed hydrodynamic description (Lu et al., 2018).

A related 2020 study introduces the notion of a “viscosity redshift”

$e^{-\beta D/2}$7

so that the amplitude suppression factor becomes $e^{-\beta D/2}$8 and is degenerate with luminosity distance in the observed waveform amplitude (Ning et al., 2020). In that treatment, viscosity alters amplitude but not frequency, unlike cosmological, Doppler, and gravitational redshifts, which affect both (Ning et al., 2020). This distinction is useful conceptually because it isolates the phenomenological signature of viscous damping: sources appear farther away in gravitational waves without an accompanying frequency shift.

The same paper applies the formalism to a dark-matter mini-spike around an intermediate-mass black hole, with a density-dependent viscosity $e^{-\beta D/2}$9, and concludes that damping can be negligible for $D$0 or $D$1, but can become strong for more extreme parameters (Ning et al., 2020). This suggests that environmental rather than truly cosmological viscosity may dominate in some astrophysical scenarios.

4. Free-streaming neutrinos and the integro-differential Weinberg equation

The collisionless-neutrino version of Weinberg damping is structurally different. Instead of a local friction term proportional to $D$2, one obtains an integro-differential equation for the tensor transfer function $D$3. Stefanek and Repko give the general-wavelength form in terms of a reduced wave number $D$4, where $D$5 is the comoving wave number entering the horizon at equality (Stefanek et al., 2012). Their $D$6-space equation is

$D$7

with $D$8 in the radiation era (Stefanek et al., 2012).

The kernel may be written as

$D$9

which makes a spherical-Bessel-series solution natural (Stefanek et al., 2012). The analytic representation

$\beta = 8\pi G\eta/c^2$0

provides a uniform description across all wavelengths (Stefanek et al., 2012).

Two limiting damping results are especially important. For $\beta = 8\pi G\eta/c^2$1, the relevant derivative ratio at last scattering gives

$\beta = 8\pi G\eta/c^2$2

corresponding to about $\beta = 8\pi G\eta/c^2$3 suppression in power (Stefanek et al., 2012). For $\beta = 8\pi G\eta/c^2$4, the asymptotic power suppression is

$\beta = 8\pi G\eta/c^2$5

equivalent to roughly $\beta = 8\pi G\eta/c^2$6 suppression in amplitude (Stefanek et al., 2012). These results refine the original Weinberg discussion by giving an explicit analytic series valid at all $\beta = 8\pi G\eta/c^2$7.

Generalizations of this free-streaming framework replace standard neutrinos by arbitrary collisionless species. One study derives the anisotropic-stress kernel for general mass and distribution function, showing that nonzero neutrino masses reduce damping for modes with $\beta = 8\pi G\eta/c^2$8–200, and that a nonthermal axion background can produce a distinct phase-dependent damping signature (Dent et al., 2013). Another paper embeds the Einstein-Boltzmann structure into restricted Horndeski gravity, finding that modified-gravity functions $\beta = 8\pi G\eta/c^2$9 alter the effective damping and can maintain non-vanishing free-streaming damping even during matter domination, unlike general relativity (Scomparin et al., 2019).

5. Geometry, modified gravity, and later-era extensions

Several works extend the Weinberg damping framework beyond the standard flat radiation-dominated setting. In spatially closed cosmology, the Laplacian spectrum is discrete and curvature modifies both the tensor equation and the neutrino Boltzmann hierarchy. For modes entering the horizon during radiation domination, the closed-universe calculation finds

$$\partial_t^2 h_{ij} + 2\beta \partial_t h_{ij} - \partial_z^2 h_{ij}=0$$0

so the squared amplitude is reduced by about $$\partial_t^2 h_{ij} + 2\beta \partial_t h_{ij} - \partial_z^2 h_{ij}=0$$1, substantially stronger than the flat-space Weinberg result $$\partial_t^2 h_{ij} + 2\beta \partial_t h_{ij} - \partial_z^2 h_{ij}=0$$2 (Khodagholizadeh et al., 2014). A subsequent study extending the analysis to matter- and $$\partial_t^2 h_{ij} + 2\beta \partial_t h_{ij} - \partial_z^2 h_{ij}=0$$3-dominated eras reports much smaller suppression for flat matter-era modes, but strong damping in the closed case, again highlighting the importance of geometry (Khodagholizadeh et al., 2016).

A different extension considers expansion itself as an effective damping mechanism in closed de Sitter spacetime. There the amplitude and frequency are extracted from the geodesic deviation equation rather than directly from the tensor Einstein equation, and the standard $$\partial_t^2 h_{ij} + 2\beta \partial_t h_{ij} - \partial_z^2 h_{ij}=0$$4 scaling of subhorizon modes in flat FRW is generalized to curvature- and $$\partial_t^2 h_{ij} + 2\beta \partial_t h_{ij} - \partial_z^2 h_{ij}=0$$5-dependent functions $$\partial_t^2 h_{ij} + 2\beta \partial_t h_{ij} - \partial_z^2 h_{ij}=0$$6 and $$\partial_t^2 h_{ij} + 2\beta \partial_t h_{ij} - \partial_z^2 h_{ij}=0$$7 (Khodagholizadeh et al., 2021). The paper argues that the cosmological-constant contribution matters mainly in the early universe, while the dominant late-time effect remains expansion-induced damping (Khodagholizadeh et al., 2021). This suggests that “Weinberg damping” can also encompass geometric redshifting and dilution, not only material dissipation.

Modified-gravity generalizations include $$\partial_t^2 h_{ij} + 2\beta \partial_t h_{ij} - \partial_z^2 h_{ij}=0$$8 gravity, where the tensor mode acquires extra friction terms proportional to $$\partial_t^2 h_{ij} + 2\beta \partial_t h_{ij} - \partial_z^2 h_{ij}=0$$9 and a scalar polarization appears with its own damping behavior (Feng et al., 2023). In that framework, Landau damping is absent in flat spacetime while collision damping is present, and for $16\pi G\eta$0 the tensor mode with $16\pi G\eta$1 decays faster than in the massless case, whereas the scalar mode with $16\pi G\eta$2 suppresses decay (Feng et al., 2023). The general lesson is that the Weinberg idea survives intact, but the effective damping kernel and the propagating degrees of freedom depend strongly on the gravitational theory.

6. Current reinterpretations and induced-gravitational-wave applications

Recent work on induced gravitational waves uses the same conceptual machinery but shifts the damping from the tensor sector to the scalar source sector. In radiation domination, finite mean free paths of photons and neutrinos damp scalar perturbations through shear viscosity, with damping scale determined by

$16\pi G\eta$3

explicitly invoking Weinberg’s 1971 viscosity treatment (Domènech et al., 17 Mar 2025). The induced-GW kernel is then computed with an exponentially damped scalar transfer function rather than a damped tensor Green’s function (Domènech et al., 17 Mar 2025).

The main effects are regularization of the resonant frequency and replacement of the usual far-infrared logarithmic running by a low-frequency tail without logarithmic running (Domènech et al., 17 Mar 2025). The paper interprets the damping parameter as an effective finite source lifetime, directly analogous to the finite-width regularization of resonances (Domènech et al., 17 Mar 2025). This is not Weinberg’s 2004 neutrino damping formula in the strict sense, but it is a clear descendant of the same dissipative logic.

The paper further notes that direct free-streaming damping of gravitational waves, as in Weinberg’s neutrino calculation, is not included in its main computation (Domènech et al., 17 Mar 2025). That separation is useful: it isolates the effect of dissipation in the source sector from dissipation in the propagation sector. A plausible implication is that a fully consistent small-scale induced-GW calculation should incorporate both.

A different modern direction considers microphysical GW-to-EM conversion in a magnetized intergalactic plasma. That mechanism is not part of Weinberg’s original anisotropic-stress formalism, but it can be cast phenomenologically as an effective attenuation coefficient in the GW dispersion relation, with strong frequency dependence at long wavelengths (Lieu et al., 2022). This suggests that the broader Weinberg paradigm now includes both hydrodynamic and microscopic damping channels.

7. Conceptual interpretation and common misconceptions

A common misconception is that “gravitational-wave damping” always means ordinary Hubble dilution. In fact the literature distinguishes at least three mechanisms. First, expansion alone redshifts subhorizon amplitudes roughly as $16\pi G\eta$4, which is geometric rather than dissipative (Khodagholizadeh et al., 2021). Second, free-streaming species generate nonlocal anisotropic stress that damps tensor modes even in otherwise standard cosmology (Stefanek et al., 2012). Third, imperfect-fluid effects such as shear viscosity produce genuine exponential attenuation $16\pi G\eta$5 with $16\pi G\eta$6 (Lu et al., 2018). These mechanisms can coexist but are not interchangeable.

A second misconception is that any viscosity contributes equally to tensor damping. The viscous-universe treatments are explicit that bulk viscosity does not attenuate transverse-traceless gravitational waves at leading order; only shear viscosity enters the damping rate (Lu et al., 2018, Ning et al., 2020).

A third misconception is that damping by matter should be substantial for ordinary astrophysical gravitational waves propagating through galaxies or the intergalactic medium. Kinetic-theory analyses show that both collisional and Landau damping by matter are negligible for essentially all astrophysical sources except primordial radiation, because maximal collisional damping is significant only when the wavelength is comparable to the horizon scale (Baym et al., 2017). This sharply limits the phenomenological importance of the Weinberg damping formula outside early-universe or deliberately exotic settings.

Finally, the various quoted damping factors should not be conflated. The viscous-fluid exponential law refers to amplitude suppression during propagation through a medium (Lu et al., 2018). The $16\pi G\eta$7 and $16\pi G\eta$8 numbers refer to tensor-mode suppression by free-streaming neutrinos in specific $16\pi G\eta$9 limits (Stefanek et al., 2012). The $Q$0 result refers to a spatially closed cosmology, not to flat FRW (Khodagholizadeh et al., 2014). Each number belongs to a distinct physical and geometric regime.

In contemporary usage, the Weinberg gravitational-wave damping formula is therefore best regarded as a family of tensor attenuation laws arising from anisotropic stress or dissipation in cosmological media, with the precise kernel determined by whether the medium is viscous, collisionless, curved-background, or modified-gravity in nature.