Fact-checked by Grok 2 weeks ago

Speed of gravity

The speed of gravity refers to the propagation velocity of gravitational influences across spacetime, which, according to Einstein's general theory of relativity, equals the speed of light in vacuum, c ≈ 299,792 km/s.^[1] This finite speed contrasts with Newtonian gravity's instantaneous action and implies that gravitational effects, such as those from accelerating masses, manifest as ripples in spacetime known as gravitational waves. The prediction has been rigorously tested and confirmed through direct observations of these waves, establishing general relativity's framework for gravity as a dynamic, wave-like phenomenon rather than a static force.^[2] In the early 20th century, Albert Einstein's development of general relativity resolved inconsistencies between Newtonian mechanics and special relativity by describing gravity as the curvature of spacetime induced by mass and energy, with perturbations propagating at c.^[1] Prior theoretical alternatives, like those proposing superluminal or subluminal speeds, were largely ruled out by subsequent analyses. Direct experimental verification emerged in 2015 with the Laser Interferometer Gravitational-Wave Observatory (LIGO) detection of GW150914, a binary black hole merger producing waves that arrived simultaneously across detectors separated by thousands of kilometers, consistent with travel at c.^[2] Further confirmation came from the 2017 multimessenger event GW170817, involving a binary neutron star merger observed via both gravitational waves and electromagnetic signals, including gamma rays. The near-simultaneous arrival of these signals—separated by just 1.7 seconds over 140 million light-years—constrained the difference between the gravitational-wave speed and c to less than 1 part in 10¹⁵, ruling out many modified gravity theories.^[3] Ongoing detections by LIGO, Virgo, and KAGRA, including those from the O4 observing run as of 2025, continue to refine these bounds and support general relativity while probing potential deviations in extreme regimes like cosmology or quantum gravity.^[4]^[5]

Historical Development

Newtonian Gravitation

In Newtonian mechanics, gravity is conceptualized as an attractive force acting instantaneously between masses, without the need for an intervening medium or propagation delay. This action-at-a-distance framework posits that the gravitational influence of a mass is felt immediately by all other masses, regardless of separation.^[6] Newton's law of universal gravitation quantifies this force as directly proportional to the product of the interacting masses and inversely proportional to the square of the distance between their centers. The key equation is

F = G \frac{m_1 m_2}{r^2},

where F is the magnitude of the gravitational force, m_1 and m_2 are the masses, r is the distance, and G is the gravitational constant. This formulation describes static gravitational fields, in which any change in a mass distribution—such as the motion of a celestial body—instantaneously alters the field everywhere, affecting distant objects without time lag.^[7]^[8] Isaac Newton introduced this theory in his seminal work Philosophiæ Naturalis Principia Mathematica, published in 1687, which unified the motions of terrestrial and celestial bodies under a single inverse-square law. The instantaneous nature of this action raised philosophical concerns regarding causality and the structure of space, as it implied forces traversing infinite distances in zero time, a concept Newton himself viewed with discomfort but accepted for its explanatory power in mechanics.^[7]^[6] This classical instantaneous model laid the groundwork for later efforts to reconcile gravity with the finite speeds observed in other physical phenomena.

Laplace's Finite-Speed Theory

In the mid-18th century, Pierre-Simon Laplace began exploring modifications to Newtonian gravity to reconcile its instantaneous action with the known finite speed of light, initially proposing in 1776 that gravitational effects might propagate via impulses from a fluid medium directed toward the attracting body. This idea addressed potential dynamical inconsistencies in celestial mechanics, particularly for systems like the Earth-Moon orbit influenced by solar motion. Laplace's early consideration appeared in a memoir examining lunar perturbations and ether influences, where he noted that a finite propagation speed could introduce subtle effects on orbital stability. Laplace expanded this concept in the 1805 volume of his Traité de mécanique céleste, performing a detailed calculation to determine the minimum propagation speed required for gravity to maintain observed planetary and satellite stabilities without detectable deviations from Newtonian predictions. Assuming gravity acts along straight lines at finite speed u, he analyzed the Earth-Moon system, where the Moon's gravitational attraction from Earth would be retarded by the time delay \tau = R/u, with R the Earth-Moon distance of approximately $3.84 \times 10^8 m. Because Earth orbits the Sun at velocity V_E \approx 3 \times 10^4 m/s, the retarded position of Earth lags behind its instantaneous position by a distance V_E \tau, introducing an aberration angle \theta \approx V_E / u. This misalignment produces a tangential component to the gravitational force on the Moon, F_\theta \approx -(GMm / R^2) (V_E / u), where G is the gravitational constant, M and m are the masses of Earth and Moon, directed opposite to the Moon's orbital velocity v_m \approx 10^3 m/s. The resulting power loss is dE/dt \approx -(GMm v_m V_E)/(R^2 u), leading to orbital energy dissipation and a gradual inward spiral of the Moon's radius according to R^{3/2} \approx R_0^{3/2} \left(1 - \frac{6\pi V_E t}{u T_0}\right), where R_0 is the initial radius and T_0 \approx 2.36 \times 10^6 s is the lunar orbital period.^[9]^[10] To ensure the Moon's orbit remains stable over the estimated age of the solar system (on the order of billions of years), without spiraling into Earth within observable timescales, Laplace derived that u must exceed roughly $7 \times 10^6 times the speed of light c \approx 3 \times 10^8 m/s, yielding a minimum u \gtrsim 2 \times 10^{15} m/s. This value, derived specifically from lunar dynamics, implied that any phase lags in gravitational propagation would be negligible for solar system scales, preserving the stability of planetary orbits against destabilizing torques. Laplace's analysis thus demonstrated that while a finite speed was conceptually appealing for consistency with light propagation, it required an extraordinarily high value to avoid contradicting astronomical observations of orbital regularity.^[9]^[11]

Electrodynamical Analogies

In the mid-18th century, Georges-Louis Le Sage proposed a kinetic theory of gravity known as push gravity, positing that gravitational attraction arises from streams of tiny, ultra-mundane corpuscles impinging on massive bodies from all directions, with the particles traveling at a finite velocity through space and creating a shadowing effect between bodies. This model treated gravity as mediated by a fluid-like medium of high-speed particles, inherently implying a finite propagation speed for gravitational influences, though Le Sage did not quantify it explicitly.^[12] Building on such mechanical ideas, Bernhard Riemann advanced a potential theory of gravity in 1853, envisioning the gravitational field as disturbances in an incompressible elastic aether where massive bodies act as sinks, leading to wave-like propagation at a finite speed determined by the medium's properties. Riemann's approach emphasized mathematical rigor in describing gravitational potentials, drawing parallels to fluid dynamics and foreshadowing field-theoretic treatments.^[13] The publication of James Clerk Maxwell's electrodynamics in the 1860s provided a powerful framework for analogies, as physicists began modeling gravity as a field governed by similar differential equations, incorporating retarded potentials to account for finite propagation delays in changing fields.^[14] These retarded effects, analogous to those in the Liénard-Wiechert potentials for accelerating charges, implied that gravitational influences from moving sources would arrive with a time lag, distorting the field's direction for accelerated masses—a phenomenon known as gravitational aberration.^[15] A seminal example is Oliver Heaviside's 1893 work, which explicitly drew an electromagnetic analogy for gravity, formulating equations where gravitational disturbances propagate as waves at the speed of light, with retarded potentials ensuring causality.^[16] Heaviside argued that "the velocity of propagation of gravitational effects is the same as that of light," predicting wave-like behavior and aberration for dynamic sources, thus linking gravity's speed directly to electromagnetic propagation.^[16] This built briefly on Pierre-Simon Laplace's earlier 1805 suggestion of finite-speed gravity to resolve orbital stability issues, though Laplace estimated a much faster speed of approximately $7 \times 10^6 times that of light.^[11]

Lorentz Covariant Models

In the early 1900s, Hendrik Lorentz sought to reconcile Newtonian gravity with his electron theory and the ether model, proposing that gravitational potentials behave analogously to electromagnetic four-vectors under Lorentz transformations. In his 1900 work, Lorentz formulated gravity as a field effect propagated through the ether at the speed of light, ensuring compatibility with the relativity principle while reducing to Newton's law for low velocities and stationary sources. This approach implied finite-speed propagation for gravitational influences, addressing inconsistencies in instantaneous action at a distance under special relativity.^[17] Building on Lorentz's ideas, Henri Poincaré extended the framework in 1905 by explicitly applying Lorentz transformations to gravitational forces in his paper "Sur la dynamique de l'électron." He argued that all forces, including gravity, must transform covariantly to preserve the relativity principle, leading to retarded potentials that propagate at the speed of light and eliminate absolute simultaneity issues in gravitational interactions. In 1908, Poincaré further refined this in his analysis of Lorentz's gravitational theory, confirming its Lorentz invariance and suggesting a unified four-dimensional treatment, though he retained an ether for interpretation while emphasizing mathematical covariance. These efforts marked the first attempts to create a fully relativistic, flat-spacetime theory of gravity.^[18]^[19] Central to these models is the covariant Poisson equation for the gravitational potential φ in weak fields:

\square \phi = \nabla_\mu \nabla^\mu \phi = -4\pi G \rho,

where □ is the d'Alembertian operator in Minkowski spacetime, G is the gravitational constant, and ρ is the mass-energy density. This equation describes a scalar gravitational field satisfying the wave equation, confirming propagation at speed c without introducing curvature. In tensor form for weak fields, it extends to a linearized metric perturbation, but remains linear overall. These theories predict no transverse gravitational waves, as the scalar nature lacks the tensor degrees of freedom present in general relativity; they agree with GR in weak fields but diverge in strong regimes due to neglected nonlinearities.^[19]

General Relativity

Theoretical Background

The theoretical foundation for the speed of gravity in general relativity originates from Albert Einstein's development of the theory between 1907 and 1916, beginning with a pivotal thought experiment on the equivalence principle and culminating in the formulation of field equations that describe gravitational propagation. In 1907, Einstein introduced the equivalence principle, stating that the effects of a uniform gravitational field are locally indistinguishable from those of a uniformly accelerated reference frame, which implies that gravitational influences cannot propagate instantaneously but must be finite to preserve causality in a relativistic framework.^[20] This principle, articulated in his review article, marked a departure from Newtonian instantaneous action at a distance and set the stage for a relativistic theory of gravity.^[21] By 1915, Einstein finalized the field equations of general relativity, which relate the curvature of spacetime to the distribution of mass and energy:

G_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}

where G_{\mu\nu} is the Einstein tensor encoding spacetime geometry, T_{\mu\nu} is the stress-energy tensor representing matter and energy, G is the gravitational constant, and c is the speed of light. These equations, presented in Einstein's November 25 paper, establish that gravity arises from the geometry of spacetime, curved by matter, and their hyperbolic structure ensures that gravitational effects propagate at a finite speed.^[22] In the weak-field limit, applicable to nearly flat spacetimes with small perturbations, general relativity linearizes to reveal wave-like behavior. Einstein's 1916 review derived this approximation, where the metric perturbation h_{\mu\nu} satisfies a wave equation propagating at the speed of light c:

\square \bar{h}_{\mu\nu} = -\frac{16\pi G}{c^4} T_{\mu\nu}

with \square the d'Alembertian operator and \bar{h}_{\mu\nu} the trace-reversed perturbation. This prediction of gravitational waves traveling at c directly follows from the relativistic invariance of the theory, confirming the finite propagation speed implied by the equivalence principle. Linearized models, such as those explored by Hendrik Lorentz in flat spacetime, served as precursors but lacked the full nonlinear curvature of general relativity.^[23]

Propagation Mechanism

In general relativity, gravitational disturbances manifest as gravitational waves, which are dynamic ripples in the curvature of spacetime generated by the asymmetric acceleration of massive objects, such as orbiting binary systems. These waves propagate outward from their sources, carrying energy and angular momentum away while altering the geodesic paths of test particles in a characteristic pattern. Unlike static gravitational fields, these perturbations are dynamic and reveal the wave-like nature of gravity.^[24] Gravitational waves are transverse, meaning their oscillations occur perpendicular to the direction of propagation, and they exhibit quadrupolar polarization modes—specifically, plus (+) and cross (×) polarizations—that stretch and squeeze spacetime in two independent directions. This quadrupolar character arises because the leading-order radiation from mass distributions vanishes for monopolar (breathing) and dipolar (oscillating) modes due to conservation laws, leaving higher-order quadrupole moments as the dominant source. In the far field, these waves behave as plane waves with wavelength much larger than the source size.^[25] The propagation speed of these waves is derived from the linearized form of the Einstein field equations, applicable in weak-field regimes where spacetime curvature is small. The metric is expressed as g_{\mu\nu} = \eta_{\mu\nu} + h_{\mu\nu}, with |h_{\mu\nu}| \ll 1, where \eta_{\mu\nu} is the Minkowski flat metric. Substituting into the full Einstein equations G_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu} and retaining first-order terms in h yields the linearized equations. Imposing the Lorenz gauge condition \partial^\mu \bar{h}_{\mu\nu} = 0, where \bar{h}_{\mu\nu} = h_{\mu\nu} - \frac{1}{2} \eta_{\mu\nu} h and h = h^\lambda_\lambda, simplifies the system to the wave equation

\square \bar{h}_{\mu\nu} = -\frac{16\pi G}{c^4} T_{\mu\nu},

with the d'Alembertian operator \square = \eta^{\alpha\beta} \partial_\alpha \partial_\beta = -\frac{1}{c^2} \partial_t^2 + \nabla^2. In vacuum (T_{\mu\nu} = 0), this reduces to the homogeneous wave equation \square \bar{h}_{\mu\nu} = 0. The general solution includes plane-wave forms \bar{h}_{\mu\nu} \propto \epsilon_{\mu\nu} e^{i(k^\alpha x_\alpha)}, where the dispersion relation \omega = c |\mathbf{k}| (with \omega the angular frequency and \mathbf{k} the wave vector) enforces propagation at exactly the speed of light c. The transverse-traceless (TT) gauge further constrains the tensor to two physical degrees of freedom, confirming the wave's tensorial structure.^[25]^[24] This tensor nature fundamentally distinguishes gravitational waves from electromagnetic waves, which are vector perturbations described by Maxwell's equations and primarily sourced by dipolar charge accelerations, resulting in two polarization states but stronger coupling via the electric charge rather than the gravitational constant G. Gravitational radiation is inherently quadrupolar, requires non-spherical mass distributions for efficient emission, and interacts far more weakly, with the characteristic power scaling as L_{GW} \sim \frac{G}{c^5} \left( \frac{d^3 Q}{dt^3} \right)^2 compared to electromagnetic L_{EM} \sim \frac{q^2 a^2}{c^3}.^[26] Albert Einstein first predicted in 1916 that gravitational disturbances propagate precisely at the speed of light, deriving this from approximate solutions to the field equations using retarded potentials. This prediction received theoretical confirmation in the 1930s, notably through the work of Einstein, Infeld, and Hoffmann, whose equations of motion for point masses emerge directly from the field equations and incorporate finite-speed propagation consistent with c.

Field Aberration and Conventions

In general relativity, the finite propagation speed of gravitational disturbances at the speed of light introduces an aberration effect in the direction of the gravitational field experienced by a moving observer relative to a source. For a weakly accelerated observer or a source with transverse velocity, the field direction lags behind the instantaneous position of the source due to the retardation in signal arrival, resulting in an angular displacement approximately given by \theta \approx v/c, where v is the relative velocity and c is the speed of light. This effect mirrors the aberration of light, where the apparent position shifts because the observer moves during the time the signal travels from the source.^[27] In the weak-field limit of general relativity, the effective aberration in the gravitational field direction can be more precisely expressed as \delta \phi = (v_\perp / c) \sin \phi, where v_\perp is the component of the source's velocity perpendicular to the line of sight from the observer to the source's retarded position, and \phi is the angle between the source's velocity vector and the line connecting the observer to the retarded position. This formula arises from expanding the retarded potentials in the post-Newtonian approximation, where the metric perturbation h_{00} and related components depend on the retarded time, leading to a directional shift proportional to the velocity over c. The derivation involves solving the linearized Einstein field equations with retarded Green's functions, incorporating the source's motion, and projecting the resulting acceleration vector onto the observer's frame; the lag \delta \phi represents the difference between the direction to the retarded and instantaneous positions, scaled by the geometric factor \sin \phi. However, this aberration is largely canceled by additional velocity-dependent terms, such as gravitomagnetic contributions from the off-diagonal metric components h_{0i}, ensuring the net field aligns closely with the instantaneous source position to leading order in v/c.^[27]^[28] To accurately model this propagation and aberration, specific conventions for gauges and coordinate systems are essential. The harmonic gauge, defined by the condition \partial^\mu \bar{h}_{\mu\nu} = 0 where \bar{h}_{\mu\nu} = h_{\mu\nu} - \frac{1}{2} \eta_{\mu\nu} h is the trace-reversed perturbation, simplifies the wave equation for metric perturbations to \square \bar{h}_{\mu\nu} = -16\pi [G](/page/G) T_{\mu\nu}/c^4, facilitating calculations of retarded solutions that propagate at speed c. For far-field gravitational waves, the transverse-traceless (TT) gauge is often adopted, where the perturbations satisfy \bar{h}_{0\mu} = 0, \partial^i \bar{h}_{ij} = 0, and \bar{h}^i_i = 0, rendering the waves transverse to the propagation direction and traceless, which isolates the physical tensor modes without spurious coordinate effects. These gauges are chosen within Lorentz-covariant coordinate systems, such as asymptotically flat Cartesian coordinates, to ensure consistency with the null propagation of gravitational signals.^[29] In stark contrast to general relativity, the Newtonian theory of gravity assumes instantaneous action at a distance, with the gravitational field determined solely by the instantaneous positions and masses of sources, resulting in no aberration or retardation effects regardless of relative motion. This instantaneous propagation eliminates any directional lag, as the field updates simultaneously across space, avoiding the velocity-induced shifts inherent in finite-speed theories.^[27]

Experimental and Observational Tests

Binary Pulsar Measurements

The Hulse-Taylor binary pulsar PSR B1913+16, discovered in 1974, consists of a 59-millisecond pulsar orbiting a compact companion neutron star in a highly eccentric 7.75-hour orbit, providing a laboratory for testing strong-field general relativity through high-precision radio timing measurements.^[30] Long-term observations have measured the relativistic periastron advance at a rate of 4.226 ± 0.001 degrees per year, precisely matching the general relativistic prediction based on the system's Keplerian parameters and total mass.^[31] Over subsequent decades, timing observations tracked the secular decay of the orbit due to energy dissipation via gravitational wave emission, with the measured orbital period derivative \dot{P}_b = -2.423 \times 10^{-12} agreeing with the general relativistic quadrupole formula prediction to within 0.2%.^[32] This decay manifests as a cumulative shift in periastron passage time, accumulating to over 40 seconds since discovery, directly evidencing the reality of gravitational radiation. The precision of these measurements, refined through observations from 1974 to 2010, confirmed the post-Keplerian parameters and earned Hulse and Taylor the 1993 Nobel Prize in Physics for the discovery and its implications for gravitational wave theory.^[33] In general relativity, the orbital frequency change rate \dot{f} due to gravitational wave energy loss for a circular orbit is given by

\dot{f} = \frac{96}{5} \pi^{8/3} \left( \frac{G \mathcal{M}}{c^3} \right)^{5/3} f^{11/3},

where f is the orbital frequency, \mathcal{M} is the chirp mass, G is the gravitational constant, and c is the speed of light; for the eccentric orbit of PSR B1913+16, additional factors account for eccentricity, yielding the observed \dot{P}_b.^[32] This formula assumes gravitational waves propagate at speed c; deviations in the propagation speed v_g would alter the radiated power by a factor of (c/v_g)^5, modifying the decay rate accordingly. The close agreement between observed and predicted decay rates thus constrains the speed of gravitational waves. Analysis of PSR B1913+16 timing data, incorporating effects like the Shklovskii acceleration, bounds the gravitational wave speed to c_T = c (1 \pm 10^{-2}) at the 1\sigma level in modified gravity models where screening mechanisms fail to fully suppress anomalies. Other binary pulsar systems provide even tighter constraints. For example, the double pulsar PSR J0737−3039, discovered in 2003, shows orbital decay agreeing with GR predictions to better than 0.1%, offering stringent tests of gravitational wave emission and propagation assumptions.^[34]

Astronomical Occultations

Astronomical occultations offer a unique opportunity to test the speed of gravity by examining the timing and effects of gravitational deflection on light or radio signals from a background source as a massive body passes in front. Early proposals in the 1960s and 1970s utilized radar ranging to measure the Shapiro time delay—the additional travel time for signals passing near the Sun to planets such as Venus—yielding lower bounds on the gravitational propagation speed of approximately v_g > 0.8c.^[35] These experiments provided indirect constraints by comparing observed delays with predictions assuming finite gravitational speeds, though the bounds were model-dependent and limited by measurement precision at the time. A landmark test occurred during the 2002–2003 Jovian occultation of the quasar J0842+1835, where Jupiter passed within 3.7 arcminutes of the quasar's line of sight on September 8, 2002. Astronomers employed very long baseline interferometry (VLBI) using the National Radio Astronomy Observatory's Very Long Baseline Array (VLBA) and additional global radio telescopes to monitor changes in the quasar's apparent radio position and scintillation patterns as Jupiter's gravitational field influenced the signal. The setup targeted potential timing discrepancies: if gravity propagated at a finite speed different from light's, the deflection would lag or lead relative to Jupiter's instantaneous position, altering the observed signal path. Observations spanned several days around the event to capture the deflection's onset and evolution, with data reduced to achieve sub-milliarcsecond positional accuracy despite challenges from Earth's rotation and instrumental effects.^[36] The analysis revealed a deflection consistent with general relativity's prediction of twice the Newtonian value, with no significant timing mismatch indicative of superluminal or subluminal gravity propagation. Interpreting the results under the assumption of finite gravitational speed yielded v_g = (0.95 \pm 0.25)c, implying the speed of gravity equals the speed of light to within 20–25% accuracy—a key empirical support for relativity if the gravitational interpretation holds. However, the measurement has faced substantial critique. Prominent relativists argued that the experiment effectively probes the propagation of light in a moving gravitational potential, not the independent speed of gravity, as the field's configuration is tied to light-speed dragging in standard general relativity frameworks; thus, deviations would reflect preferred-frame effects (parameterized post-Newtonian \alpha_1 \neq 0) rather than v_g \neq c. Solar system tests already constrain such effects tightly, rendering the result insensitive to v_g at leading order.^[37] Further contention arose from potential plasma effects in Jupiter's magnetosphere, which could induce refractive scintillation in the radio signal, mimicking or masking gravitational deflections through ionized particle scattering. While Fomalont and Kopeikin modeled and subtracted plasma contributions using multi-frequency data, skeptics questioned whether residual ionospheric and magnetospheric interferences fully accounted for observed position shifts, potentially biasing the timing analysis toward light-speed equivalence without confirming gravitational propagation directly. Despite these debates, the event remains a pioneering effort to parameterize gravitational speed via occultation dynamics, highlighting the interplay of aberration and field propagation in relativistic tests.^[38]

Gravitational Wave Events

The detection of gravitational waves from the binary neutron star merger GW170817 on August 17, 2017, by the LIGO and Virgo observatories marked the first multimessenger astronomical event, providing a direct probe of the speed of gravity through its association with the gamma-ray burst GRB 170817A.^[39] The gravitational wave signal arrived at Earth, followed 1.74 ± 0.05 seconds later by the gamma-ray emission from GRB 170817A, detected by the Fermi Gamma-ray Burst Monitor and the Integral SPI-ACS instrument.^[39] The source was localized to the galaxy NGC 4993 at a luminosity distance of approximately 40 Mpc.^[40] This short time delay between the gravitational and electromagnetic signals enabled a precise constraint on the propagation speed of gravitational waves, v_g, relative to the speed of light, c. Assuming the signals originated from the same astrophysical source, the time delay \Delta t arises from the difference in propagation speeds over the distance D, given by the approximation

\Delta t = D \left| \frac{1}{v_g} - \frac{1}{c} \right|

.^[39] Using the observed \Delta t \approx 1.7 s and D \approx 40 Mpc (corresponding to a light-travel time of about $1.3 \times 10^8 years), the analysis yields |v_g - c| / c < 5 \times 10^{-16}, specifically between -3 \times 10^{-15} and +7 \times 10^{-16}.^[39] This bound, consistent with general relativity's prediction that v_g = c, rules out numerous modified gravity theories that predict significant deviations in gravitational wave speeds, such as certain scalar-tensor models and massive graviton theories. Subsequent gravitational wave detections during the LIGO-Virgo-KAGRA observing runs have not yielded confirmed electromagnetic counterparts capable of similarly constraining the speed of gravity. As of November 2025, the fourth observing run (O4, spanning 2023–2025) has identified approximately 300 events, but none with verified multimessenger associations.^[41] For instance, GW230529, a candidate neutron star–black hole merger detected in May 2023, prompted extensive electromagnetic follow-up searches, yet no counterpart was confirmed, limiting further direct tests of propagation speed.^[42]

Implications and Future Prospects

Theoretical Consequences

If the speed of gravity, v_g, were greater than the speed of light c, it would permit superluminal signaling, enabling closed timelike curves and causality violations that contradict the principles of special relativity.^[43] Such a scenario would undermine the causal structure of spacetime, allowing effects to precede causes in certain reference frames, a fundamental inconsistency in Lorentz-invariant theories.^[44] Conversely, if v_g < c, as in theories of massive gravity, it introduces the van Dam-Veltman-Zakharov (vDVZ) discontinuity, where the weak-field limit fails to recover general relativity, leading to unphysical discrepancies in solar system dynamics such as incorrect planetary deflection angles. This discontinuity manifests as an extra scalar mode that alters gravitational interactions at low energies, causing instabilities like ghost fields or tachyonic modes in the theory's spectrum.^[45] In cosmology, deviations from v_g = c would profoundly impact large-scale structure formation and the evolution of the universe. Slower gravitational propagation would delay the clustering of matter perturbations, suppressing the growth of density fluctuations and altering the power spectrum observed in cosmic microwave background anisotropies or galaxy distributions.^[46] Faster propagation, meanwhile, could accelerate black hole mergers and gravitational wave emission rates, modifying the stochastic gravitational wave background and the merger rates inferred from observations.^[47] These effects would disrupt the standard \LambdaCDM model's predictions for cosmic expansion and structure, requiring significant revisions to inflationary scenarios and dark energy parameters. Modified gravity theories, such as scalar-tensor models like Brans-Dicke theory, generally predict v_g \approx c for tensor gravitational waves, with an additional scalar mode also propagating at c, but potential deviations arise from the scalar field's coupling strength parameterized by \omega_{BD}. Observations bound these deviations tightly; for instance, the Brans-Dicke parameter must satisfy \omega_{BD} > 40000 to align with solar system tests, ensuring negligible differences from general relativity.^[48] In these frameworks, any residual speed mismatch is constrained by multimessenger events, confirming consistency with v_g = c within $10^{-15}. The empirical verification that v_g = c philosophically reinforces general relativity's unification of gravity as spacetime curvature rather than an independent force, affirming Einstein's vision where gravitational influences are inseparable from the causal limits of relativity. This equivalence underscores the theory's elegance, integrating gravity into the relativistic framework without ad hoc adjustments for propagation speed.

Ongoing and Proposed Tests

The LIGO-Virgo-KAGRA (LVK) collaboration's fourth observing run (O4), which began in May 2023 and concluded in November 2025, detected over 200 gravitational wave events by March 2025 and approximately 300 by the run's end due to improved detector sensitivities.^[49]^[50] These observations include systematic searches for electromagnetic counterparts, which aim to refine bounds on the speed of gravity beyond the benchmark established by GW170817 in 2017, by comparing arrival times of gravitational waves and photons from the same cosmic events.^[51] Pulsar timing arrays, such as those operated by the NANOGrav collaboration since 2023, provide indirect constraints on the speed of gravity through analyses of the stochastic gravitational wave background. By examining timing residuals in millisecond pulsar signals influenced by nanohertz-frequency waves, recent studies using the 15-year dataset have derived stringent upper limits on the graviton mass, implying that the propagation speed of gravitational waves is very close to the speed of light.^[52] Ongoing data collection and international PTA collaborations continue to enhance these bounds as more years of observations accumulate. Proposed future experiments include the DECi-hertz Interferometer Gravitational-wave Observatory (DECIGO), a space-based detector targeting the 0.1–10 Hz band, which could enable direct, model-independent measurements of the gravitational wave speed by resolving arrival time differences in signals from binary neutron star mergers or other sources.^[53] Additionally, time-delay cosmography using the James Webb Space Telescope (JWST) offers prospects for joint estimation of the Hubble constant and gravitational wave speed through precise measurements of arrival time delays in strongly lensed supernovae and quasars, which can be cross-correlated with gravitational wave detections in multimessenger frameworks.^[54]^[55] Recent 2025 proposals emphasize model-independent tests of the gravitational wave speed via galaxy-scale strong lensing and supernova timing delays, leveraging datasets from surveys like DESI and JWST to compare electromagnetic time delays with predicted gravitational wave propagations, potentially achieving sub-percent precision in cosmological parameter constraints.^[56]^[57] These approaches build on multi-messenger lensing techniques to probe deviations from general relativity without assuming specific propagation models.^[58]

References

[1]
[PDF] Observation of Gravitational Waves from a Binary Black Hole Merger
Feb 11, 2016 · He found that the linearized weak-field equations had wave solutions: transverse waves of spatial strain that travel at the speed of light, ...
[2]
Observation of Gravitational Waves from a Binary Black Hole Merger
Feb 11, 2016 · He found that the linearized weak-field equations had wave solutions: transverse waves of spatial strain that travel at the speed of light, ...
[3]
Bounding the Speed of Gravity with Gravitational Wave Observations
Of order 20 detections by the two LIGO detectors will constrain the speed of gravity to within 20% of the speed of light, while just five detections by the LIGO ...
[4]
Tests of gravitational wave propagation with LIGO-Virgo catalog - arXiv
Apr 23, 2024 · Abstract:In the framework of general relativity (GR), gravitational waves (GWs) travel at the speed of light across all frequencies.
[5]
[PDF] Isaac Newton on the action at a distance in gravity - PhilArchive
Feb 16, 2019 · In Opticks Query 29, Newton states: "Pellucid. Substances act upon the Rays of Light at a distance in refracting, reflecting, and inflecting ...
[6]
Newton's Philosophiae Naturalis Principia Mathematica
Dec 20, 2007 · The historical context in which Newton wrote the Principia involved a set of issues that readers of the first edition saw it as addressing ...
[7]
Philosophiae naturalis principia mathematica : Newton, Isaac, 1642 ...
Jan 19, 2017 · This is the first edition of Newton's Principia, in which he elucidates the universal laws of gravitation and motion that underlay the phenomena described by ...
[8]
[PDF] Laplace and the Speed of Gravity 1 Problem 2 Solution
May 5, 2015 · We assume that the force of gravity acts along straight lines, but with finite speed u. Then, the force exerted at time t = 0 on the Moon by the ...
[9]
http://kirkmcd.princeton.edu/examples/laplace.pdf
[10]
Laplace on the Speed of Gravity - MathPages
Therefore, the finite speed of the posited gravitational particles in a shadow theory pulls the satellite forward, adding angular momentum, which (as discussed ...
[11]
Lesage's Shadows - MathPages
One of the most famous of the attempts to identify a mechanism for gravity was the "shadow theory" proposed by Georges-Louis Lesage (1724-1803), who is usually ...
[12]
Bernhard Riemann, a(rche)typical mathematical-physicist? - Frontiers
He chose for the first the subject of Fourier series, presenting an essay in 1853, in which he gave a criterion for a function to be (Riemann) integrable, and ...
[13]
[PDF] Gravitomagnetism and gravitational waves - arXiv
Feb 11, 2010 · Abstract. After extensively reviewing general relativistic gravitomagnetism, both historically and phenomenologically, we review in detail ...
[14]
[PDF] Aberration and the Speed of Gravity - arXiv
the speed of gravity? The answer, unfortunately, is that so far they say fairly little. In the absence of direct measurements of propagation speed ...
[15]
Electromagnetic theory : Heaviside, Oliver, 1850-1925
Mar 27, 2007 · Appendix B: A gravitational and electromagnetic analogy.--v.2. V ... FULL TEXT download · download 1 file · HOCR download · download 1 file.
[16]
Lorentz ether theory - chemeurope.com
In 1900 Lorentz tried to explain gravity on the basis of the Maxwell equations. He first considered a Le Sage type model and argued that there possibly exists ...<|control11|><|separator|>
[17]
https://www.chemeurope.com/en/encyclopedia/Lorentz_ether_theory.html
[18]
(PDF) Breaking in the 4-Vectors: The Four-Dimensional Movement ...
PDF | In July, 1905, Henri Poincaré (1854–1912) proposed two laws of gravitational attraction compatible with the principle of relativity and all.Missing: gravity | Show results with:gravity
[19]
[PDF] on the relativity principle and the conclusions drawn from it
In the section "Principle of relativity and gravitation", a reference system at rest situated in a temporally constant, homogeneous gravitational field is ...
[20]
[PDF] Einstein's Principle of Equivalence
gravity must be confined in a finite region; otherwise the light speed as the maximum velocity would be violated. In other words,. 22 uniform gravity like an ...
[21]
[PDF] Die Feldgleichungen der Gravitation.
799, 1915. Page 2. EINSTEIN: Die Feldgleichungen der Gravitation. 845. Die allgemein kovarianten zehn Gleichungen des Gravitationsfeldes in Räumen, in denen ...
[22]
[PDF] Die Grundlage der allgemeinen Relativitätstheorie
1916. ANNALEN DER PHYSIK. VIERTE FOLGE. BAND 49. 1. Die Grundlage der allgemeinen Relativitätstheorie; von A. Einstein. № 7. +. Die im nachfolgenden ...
[23]
[PDF] Gravitation - physicsgg
This is a textbook on gravitation physics (Einstein's "general relativity" or "geo- metrodynamics"). It supplies two tracks through the subject.
[24]
https://physicsgg.me/wp-content/uploads/2023/05/misner_thorne_wheeler_gravitation_freema.pdf
[25]
Comparison of electromagnetic and gravitational radiation: What we ...
Aug 1, 2013 · There are many similarities between electromagnetic (EM) radiation and gravitational radiation: both travel at the speed of light; both carry ...
[26]
[PDF] Comparison of electromagnetic and gravitational radiation: What we ...
This difference can be partially ascribed to the difference between the nature of eigenlines for a tensor field and “lines of force” for a vector field. But ...
[27]
[PDF] Näherungsweise Integration der Feldgleichungen der Gravitation
Albert Einstein (1916). Näherungsweise Integration der Feldgleichungen der Gravitation. Königlich Preußische Akademie der Wissenschaften (Berlin).
[28]
[PDF] Albert Einstein, Leopold Infeld and Banesh Hoffmann (1938)
Albert Einstein, Leopold Infeld and Banesh Hoffmann (1938). The Gravitational Equa- tions and the Problem of Motion. The Annals of Mathematics, 39: 65–100.
[29]
[gr-qc/9909087] Aberration and the Speed of Gravity - arXiv
Sep 29, 1999 · Title:Aberration and the Speed of Gravity. Authors:S. Carlip. View a PDF of the paper titled Aberration and the Speed of Gravity, by S. Carlip.
[30]
Gravimagnetism, causality, and aberration of gravity in the ...
Aug 11, 2007 · We prove that the extrinsic gravimagnetic field is generated by the relativistic effect of the aberration of the gravity force caused by the ...<|control11|><|separator|>
[31]
6. Weak Fields and Gravitational Radiation
... of the harmonic gauge known as the transverse traceless gauge (or sometimes "radiation gauge"). The name comes from the fact that the metric perturbation is ...
[32]
https://iopscience.iop.org/article/10.1088/0004-637X/722/2/1030
[33]
Measurements of general relativistic effects in the binary pulsar ...
Feb 8, 1979 · Measurements of second- and third-order relativistic effects in the orbit of binary pulsar PSR1913 + 16 have yielded self-consistent estimates of the masses of ...
[34]
TIMING MEASUREMENTS OF THE RELATIVISTIC BINARY ...
Sep 24, 2010 · 1. INTRODUCTION. Pulsar PSR B1913+16 (PSR J1915+1606) was the first binary pulsar to be discovered (Hulse & Taylor 1975).
[35]
Press release: The 1993 Nobel Prize in Physics - NobelPrize.org
Hulse's and Taylor's discovery in 1974 of the first binary pulsar, called PSR 1913 + 16 (PSR stands for pulsar, and 1913 + 16 specifies the pulsar's position in ...
[36]
The Confrontation between General Relativity and Experiment
3 Shapiro time delay and the speed of gravity. In 2001, Kopeikin [147] suggested that a measurement of the time delay of light from a quasar as the light ...
[37]
Aberration and the Fundamental Speed of Gravity in the Jovian ...
Abstract: We describe our explicit Lorentz-invariant solution of the Einstein and null geodesic equations for the deflection experiment of 2002 September 8 ...Missing: (v_perp / sin φ
[38]
Propagation Speed of Gravity and the Relativistic Time Delay - arXiv
Authors:Clifford M. Will. View a PDF of the paper titled Propagation Speed of Gravity and the Relativistic Time Delay, by Clifford M. Will. View PDF. Abstract ...
[39]
First speed of gravity measurement revealed | New Scientist
Jan 7, 2003 · Fomalont and Kopeikin combined observations from a series of radio telescopes across the Earth to measure the apparent change in the ...
[40]
Gravitational Waves and Gamma-Rays from a Binary Neutron Star ...
The association of GW170817 and GRB 170817A provides new insight into fundamental physics and the origin of short GRBs.
[41]
Early spectra of the gravitational wave source GW170817 - Science
Oct 16, 2017 · Two neutron stars merging together generate a gravitational wave signal and have also been predicted to emit electromagnetic radiation.
[42]
LIGO-Virgo-KAGRA Announce the 200th Gravitational Wave ...
Mar 20, 2025 · Since the “advanced” era of gravitational wave detection began in 2015, the LVK network has now (as of March 19, 2025) recorded 290 GW events, ...
[43]
LIGO-Virgo-KAGRA (LVK) Collaboration Detected a Remarkable ...
Apr 5, 2024 · What makes this signal, called GW230529, intriguing is the mass of the heavier object. It falls within a possible mass-gap between the heaviest ...
[44]
[hep-th/0011152] No Van Dam-Veltman-Zakharov Discontinuity in ...
Nov 16, 2000 · We prove that the van Dam-Veltman-Zakharov discontinuity arising in the massless limit of massive gravity theories is peculiar to Minkowski space.
[45]
"Speed" of Gravity and Speed of Light - Physics Stack Exchange
Mar 16, 2011 · The special theory of relativity is really enough to see that gravitational signals have to propagate by the speed c which we call "speed of ...
[46]
Van Dam-Veltman-Zakharov discontinuity in topologically new ...
Dec 31, 2012 · We study van Dam-Veltman-Zakharov discontinuity in the topologically new massive gravity (TNMG). The reduction from 2 degrees of freedom to one ...
[47]
Measuring the speed of cosmological gravitational waves - NASA ADS
We investigate the effects of a modified speed of gravity, cT2, on the B modes of the cosmic microwave background (CMB) anisotropy in polarization. We find that ...
[48]
Gravitational waves in Brans-Dicke Theory with a cosmological ...
Jan 10, 2021 · Weak field gravitational wave solutions are investigated in Brans-Dicke (BD) theory in the presence of a cosmological constant.
[49]
Gravitational waves in Brans–Dicke theory with a cosmological ...
Apr 17, 2021 · Weak field gravitational wave solutions are investigated in Brans–Dicke (BD) theory in the presence of a cosmological constant.
[50]
Record detection of 200 gravitational waves in the current ... - Virgo
Mar 20, 2025 · The observing period O4 has been running since May 2023, with a short break at the beginning of 2024 after which the Virgo detector joined the ...Missing: speed | Show results with:speed
[51]
Gravitational wave events hint at 'second-generation' black holes
Oct 28, 2025 · The current run started in late May 2023 and is expected to continue through mid-November of this year. To date, approximately 300 black hole ...
[52]
GWTC-4.0: Updated Gravitational-Wave Catalog Released | LIGO Lab
Aug 26, 2025 · August 26th 2025 marks an important date for the international LIGO, Virgo, and KAGRA (LVK) gravitational-wave network. The collaborations are ...
[53]
Unveiling the graviton mass bounds through the analysis of 2023 ...
Mar 18, 2024 · In this work, we investigate the graviton mass bounds by analyzing the observational data of the overlap reduction functions from the NANOGrav 15-year data ...
[54]
Direct measurements of the speed of gravitational waves using the ...
Oct 7, 2025 · In this Letter, we focus on a new approach to directly measuring vg with the DECi-hertz Interferometer Gravitational-wave Observatory (DECIGO), ...
[55]
(PDF) SN H0pe: The First Measurement of H 0 from a Multiply ...
Follow-up JWST observations of all three SN images enabled photometric and rare spectroscopic measurements of the two relative time delays. Following strict ...
[56]
Multi-messenger gravitational lensing - Journals
May 1, 2025 · Multi-messenger gravitational lensing advances in cosmology will be driven by the complementary arrival time and angular position accuracies of ...
[57]
https://ui.adsabs.harvard.edu/abs/2025SCPMA..6880407R/abstract
[58]
Forecast of gravitationally lensed Type Ia supernovae time delay ...
Strong lensing time delay measurement is a promising method to address the Hubble tension, offering a completely independent approach compared with both the ...
[59]
[PDF] Multi-messenger gravitational lensing | WRAP: Warwick
2017 Testing the speed of gravitational waves over cosmological dis- tances with strong gravitational lensing. Phys. Rev. Lett. 118, 091101. (doi:10.1103 ...<|control11|><|separator|>