Fact-checked by Grok 2 weeks ago

Geodesics in general relativity

In general relativity, geodesics are the worldlines traced by test particles in free fall, generalizing the concept of straight lines to curved spacetime and serving as the paths of extremal proper time or length determined solely by the geometry of spacetime.^[1] These trajectories embody the principle that gravity is not a force but a manifestation of spacetime curvature, such that freely falling objects follow the "straightest" possible paths without acceleration relative to local inertial frames.^[2] The geodesic hypothesis, a foundational postulate of the theory, asserts that all free particles—regardless of their mass or composition—trace out these geodesics, unifying inertial and gravitational mass in a way that resolves long-standing puzzles in Newtonian gravity.^[2] The motion along geodesics is mathematically described by the geodesic equation, \frac{d^2 x^\mu}{d\lambda^2} + \Gamma^\mu_{\alpha\beta} \frac{dx^\alpha}{d\lambda} \frac{dx^\beta}{d\lambda} = 0, where x^\mu are spacetime coordinates, \lambda is an affine parameter (such as proper time for massive particles), and \Gamma^\mu_{\alpha\beta} are the Christoffel symbols encoding the metric's curvature.^[1] This equation arises from the requirement that the tangent vector to the curve undergoes parallel transport, preserving its direction relative to the local geometry, and can be derived variationally by extremizing the proper time functional \tau = \int \sqrt{-g_{\mu\nu} \frac{dx^\mu}{d\lambda} \frac{dx^\nu}{d\lambda}} \, d\lambda, where g_{\mu\nu} is the metric tensor.^[3] In coordinate systems adapted to local inertial frames, the equation simplifies to straight-line motion at constant velocity, but in curved coordinates reflecting gravitational fields, it captures deviations like orbital precession or light bending. Geodesics are classified by the sign of their tangent vector's norm: timelike geodesics (g_{\mu\nu} \frac{dx^\mu}{d\lambda} \frac{dx^\nu}{d\lambda} < 0) describe the paths of massive particles, maximizing proper time; null geodesics (= 0) govern photon trajectories, explaining phenomena like gravitational lensing; and spacelike geodesics (> 0) are less common but relevant for certain spacelike separations in spacetime.^[1] This classification underscores the theory's Lorentzian structure, where the causal structure of spacetime dictates allowable paths. The geodesic equation's universality stems from the equivalence principle, which equates gravitational fields to accelerated frames, ensuring that all observers agree on free-fall trajectories locally.^[2] Notable applications include the prediction of Mercury's perihelion advance, confirmed observationally as a 43 arcseconds per century deviation from Newtonian orbits due to spacetime curvature around the Sun.^[4] Null geodesics also account for the deflection of starlight during solar eclipses, verified in 1919 and foundational to general relativity's acceptance.^[5] In strong-field regimes, such as around black holes, geodesics delineate event horizons and photon spheres, influencing phenomena like gravitational wave signals from merging compact objects and the shadows observed by the Event Horizon Telescope.^[6] These paths not only test the theory but also form the basis for modeling cosmic structure formation and high-precision geodesy in relativistic contexts.^[7]

Introduction

Definition and Role in General Relativity

In general relativity, a geodesic is defined as an extremal curve in a pseudo-Riemannian manifold that generalizes the concept of a straight line from Euclidean space, representing the path of shortest or longest proper length between two points depending on the signature of the metric.^[8] These curves are solutions to the geodesic equation and embody the intrinsic geometry of spacetime.^[9] Geodesics play a central role in general relativity by describing the worldlines of test particles in free fall under the influence of gravity alone, with no non-gravitational forces acting on them.^[8] In this framework, the motion of such particles is determined solely by the curvature of spacetime, as dictated by the equivalence principle, where all freely falling bodies follow the same trajectory regardless of their mass or composition.^[10] The geometry of spacetime is encoded in the metric tensor g_{\mu\nu}, which defines distances and angles, allowing geodesics to predict the paths of particles without requiring solutions to the full Einstein field equations that couple matter to geometry.^[8] To ensure the geodesic equation is independent of the choice of parameter along the curve, an affine parameterization is employed, preserving the form of the equation under reparameterizations of the form \lambda \to a\lambda + b.^[8] In a coordinate basis, the geodesic equation takes the form

\frac{d^2 x^\mu}{d\tau^2} + \Gamma^\mu_{\alpha\beta} \frac{dx^\alpha}{d\tau} \frac{dx^\beta}{d\tau} = 0,

where \tau is an affine parameter (often proper time for timelike geodesics), and the Christoffel symbols \Gamma^\mu_{\alpha\beta} serve as the connection coefficients that quantify how the metric tensor varies across spacetime, capturing the effects of curvature on parallel transport.^[8]

Historical Development

The concept of geodesics originated in the mid-19th century within the framework of non-Euclidean geometry, where Bernhard Riemann introduced the idea of manifolds equipped with a metric structure in his 1854 habilitation lecture delivered at the University of Göttingen.^[11] In this foundational work, Riemann described continuous manifolds where metric relations determine distances and angles, laying the groundwork for curved spaces that would later underpin general relativity, though the full implications for geodesics as extremal paths emerged only with subsequent developments.^[12] The development of geodesics advanced significantly in the early 20th century through contributions in differential geometry, particularly Tullio Levi-Civita's work on parallel transport and affine connections between 1916 and 1917. In his 1916 paper, Levi-Civita explored intrinsic connections on hypersurfaces, providing a geometric framework for transporting vectors without torsion, which directly informed the notion of geodesics as auto-parallel curves.^[13] His 1917 publication further specified this parallelism in arbitrary manifolds, linking it to Riemann's curvature and enabling a precise definition of geodesic paths as those preserving direction under parallel transport, a concept essential for curved spacetime descriptions.^[14] Albert Einstein integrated geodesics into general relativity during 1915–1916, realizing that gravity manifests as the curvature of spacetime, with geodesics representing the inertial paths of freely falling objects analogous to straight lines in flat space. This insight, rooted briefly in the equivalence principle equating gravitational and inertial mass, allowed Einstein to predict the anomalous precession of Mercury's perihelion using approximate geodesic equations in a weak-field limit, achieving exact agreement with observations in his November 1915 paper. In his comprehensive 1916 exposition, Einstein formalized geodesics as solutions to the equations of motion derived from the spacetime metric, solidifying their role as the trajectories followed by matter and light in gravitational fields.^[15] A key milestone came in 1916 with Karl Schwarzschild's exact solution to Einstein's field equations for a spherically symmetric, non-rotating mass, which yielded the first precise geodesics describing orbits around such sources and later interpreted as paths around black holes. The predictive power of these null geodesics was empirically confirmed during the 1919 solar eclipse expedition led by Arthur Eddington, which measured the deflection of starlight grazing the Sun's edge, matching GR's forecast of 1.75 arcseconds to within experimental error and validating geodesics for light propagation.^[16] Post-Einstein refinements included John L. Synge's 1960 textbook, which clarified the geodesic deviation equation and its relation to tidal effects, providing a rigorous geometric interpretation of how nearby geodesics separate due to spacetime curvature, thus enhancing understanding of gravitational influences on extended bodies.^[17]

Mathematical Formulation

The Geodesic Equation

In general relativity, the geodesic equation governs the motion of test particles and light rays in curved spacetime, describing paths that extremalize the proper time or affine length. The standard form of this equation, in a coordinate basis with general coordinates x^\mu, is given by

\frac{d^2 x^\mu}{d\lambda^2} + \Gamma^\mu_{\alpha\beta} \frac{dx^\alpha}{d\lambda} \frac{dx^\beta}{d\lambda} = 0,

where \lambda is an affine parameter along the curve, and \Gamma^\mu_{\alpha\beta} are the Christoffel symbols of the second kind, which encode the geometry of spacetime via the metric tensor g_{\mu\nu}.^[18] This second-order differential equation arises from the requirement that the tangent vector to the curve is parallel transported along itself, ensuring the path is "straight" in the sense of the manifold's connection.^[19] The Christoffel symbols are uniquely determined by the metric and its partial derivatives through the formula

\Gamma^\mu_{\alpha\beta} = \frac{1}{2} g^{\mu\sigma} \left( \partial_\alpha g_{\beta\sigma} + \partial_\beta g_{\alpha\sigma} - \partial_\sigma g_{\alpha\beta} \right),

where g^{\mu\sigma} is the inverse metric.^[18] This expression follows from the Levi-Civita connection, which is the unique torsion-free connection compatible with the metric, satisfying \nabla_\rho g_{\mu\nu} = 0 to preserve distances under parallel transport.^[19] The symbols are symmetric in their lower indices, \Gamma^\mu_{\alpha\beta} = \Gamma^\mu_{\beta\alpha}, reflecting the absence of torsion.^[18] The parameter \lambda must be affine for the equation to hold in this form; under a reparameterization \lambda \to f(\lambda), the equation transforms to a non-affine version

\frac{d^2 x^\mu}{d\tilde{\lambda}^2} + \Gamma^\mu_{\alpha\beta} \frac{dx^\alpha}{d\tilde{\lambda}} \frac{dx^\beta}{d\tilde{\lambda}} = f'(\lambda) \frac{dx^\mu}{d\tilde{\lambda}},

where the extra term arises unless f(\lambda) is linear, \tilde{\lambda} = a\lambda + b.^[18] Affine parameterization scales invariantly, allowing \lambda to represent proper time for timelike geodesics or an arbitrary scaling for null ones, without altering the geometric path.^[19] In tensorial notation, the geodesic equation expresses the covariant derivative of the tangent vector along the curve vanishing:

\frac{D}{d\lambda} \frac{dx^\mu}{d\lambda} = 0,

or equivalently, \frac{dx^\sigma}{d\lambda} \nabla_\sigma \left( \frac{dx^\mu}{d\lambda} \right) = 0, where \nabla is the Levi-Civita covariant derivative.^[18] This form highlights the intrinsic, coordinate-independent nature of geodesics as autoparallels in the manifold.^[19] Given an initial position x^\mu(0) and initial velocity \frac{dx^\mu}{d\lambda}(0), standard ordinary differential equation theory guarantees the local existence and uniqueness of the solution to the geodesic equation in a sufficiently small neighborhood, assuming the metric is smooth.^[19] This theorem ensures that geodesics are well-defined locally as the unique curves satisfying the parallel transport condition for given initial data.^[18]

Parameter Choices and Equivalent Forms

In general relativity, the geodesic equation takes its simplest form when the curve is parameterized affinely, meaning the parameter \lambda is chosen such that the tangent vector satisfies g_{\mu\nu} \frac{dx^\mu}{d\lambda} \frac{dx^\nu}{d\lambda} = constant along the curve, ensuring no extraneous terms appear in the equation \frac{d^2 x^\mu}{d\lambda^2} + \Gamma^\mu_{\nu\rho} \frac{dx^\nu}{d\lambda} \frac{dx^\rho}{d\lambda} = 0.^[19] This affine parameterization preserves the parallel transport property of the tangent vector and allows linear reparameterizations \tilde{\lambda} = a\lambda + b (with a \neq 0) without altering the form of the equation.^[20] For timelike geodesics, which describe the worldlines of massive particles, the affine parameter is conventionally the proper time \tau, defined such that the normalization condition holds: g_{\mu\nu} \frac{dx^\mu}{d\tau} \frac{dx^\nu}{d\tau} = -1 in the mostly-plus signature (-,+,+,+).^[19] This choice yields the four-velocity u^\mu = \frac{dx^\mu}{d\tau}, satisfying u_\mu u^\mu = -1, which physically represents the invariant interval along the particle's path.^[20] Null geodesics, followed by massless particles such as photons, cannot use proper time since the interval vanishes, ds^2 = 0, so an arbitrary affine parameter \lambda is employed with the tangent vector k^\mu = \frac{dx^\mu}{d\lambda} obeying k^\mu k_\mu = 0.^[19] This parameterization maintains the geodesic equation in its standard form while ensuring the null condition is preserved under affine rescalings.^[21] In practical computations, especially for timelike geodesics in stationary spacetimes, coordinate time t (where x^0 = t) can parameterize the motion, leading to equations for spatial coordinates like \frac{d^2 x^i}{dt^2} + \Gamma^i_{\alpha\beta} \frac{dx^\alpha}{dt} \frac{dx^\beta}{dt} = 0, with \frac{dt}{dt} = 1 simplifying terms such as \Gamma^i_{00}.^[20] This form is particularly useful in weak-field approximations, where it recovers Newtonian equations of motion.^[19] The geodesic equation exhibits reparameterization freedom: under a non-affine change \lambda \to f(\lambda) with f'(\lambda) \neq 0, the equation modifies to \frac{d^2 x^\mu}{d\lambda^2} + \Gamma^\mu_{\nu\rho} \frac{dx^\nu}{d\lambda} \frac{dx^\rho}{d\lambda} = \phi(\lambda) \frac{dx^\mu}{d\lambda}, where \phi(\lambda) is a scalar function proportional to the derivative of the reparameterization.^[22] This extra term arises because non-affine parameters do not preserve the constant norm of the tangent vector, but the underlying curve remains a geodesic.^[21]

Derivations of the Geodesic Equation

From the Equivalence Principle

The equivalence principle, a cornerstone of general relativity, posits that in a sufficiently small region of spacetime, the effects of a uniform gravitational field are indistinguishable from those experienced in a uniformly accelerated reference frame, provided tidal forces—arising from spacetime curvature gradients—are negligible.^[23] In such a local domain, a freely falling observer's frame behaves as an inertial frame of special relativity, where test particles at rest relative to the frame remain at rest, or move with constant velocity along straight lines.^[23] This local flatness implies that the paths of freely falling particles are straight lines in these coordinates, establishing an intuitive physical basis for geodesic motion without invoking global variational principles.^[24] To formalize this, one can introduce local inertial coordinates at a chosen event P in spacetime, where the metric tensor g_{\mu\nu} takes the Minkowski form \eta_{\mu\nu} exactly at P, and the Christoffel symbols \Gamma^\lambda_{\mu\nu}, which measure deviations from flatness, vanish at P.^[24] Nearby, the metric deviates slightly, expressed as an expansion g_{\mu\nu} = \eta_{\mu\nu} + h_{\mu\nu}, where h_{\mu\nu} is small and captures the influence of curvature; the Christoffel symbols, proportional to the derivatives of h_{\mu\nu}, then quantify the tidal fields that cause paths to curve relative to local straight lines.^[24] In these coordinates, the equation of motion for a particle follows from the requirement that its four-acceleration vanishes in the inertial frame, but transforming back to general coordinates reveals the curvature-induced corrections. The derivation proceeds by considering the motion of a test particle in an accelerated frame equivalent to a gravitational field. In the local inertial frame, the particle satisfies the flat-space equation \frac{d^2 x'^\mu}{d\tau^2} = 0, indicating zero proper acceleration for free fall.^[24] Under a coordinate transformation to a non-inertial frame (mimicking gravity), the four-velocity v^\mu = dx^\mu / d\tau transforms, and the second derivative acquires additional terms from the acceleration of the frame. Specifically, the proper acceleration in curved coordinates becomes a^\mu = \frac{d^2 x^\mu}{d\tau^2} + \Gamma^\mu_{\alpha\beta} v^\alpha v^\beta, where the Christoffel symbols encode the "fictitious" forces due to the coordinate choice.^[24] For freely falling particles, the proper acceleration must vanish (a^\mu = 0), yielding the geodesic equation:

\frac{d^2 x^\mu}{d\tau^2} + \Gamma^\mu_{\alpha\beta} \frac{dx^\alpha}{d\tau} \frac{dx^\beta}{d\tau} = 0.

This equation describes how paths deviate from straight lines due to the varying gravitational field, with the Christoffel symbols serving as measures of the tidal distortions in the local frame.^[24] A classic illustration is Einstein's elevator thought experiment: imagine an elevator in free fall within a uniform gravitational field, where objects inside appear weightless and follow straight paths relative to the elevator.^[23] If the field varies spatially (as in a real gravitational source), tidal effects cause relative accelerations among particles, curving their worldlines in the elevator's coordinates and manifesting as geodesic motion; for instance, in Earth's field, a dropped ball follows a parabolic path in the lab frame but a straight line locally in the falling frame, with curvature arising from the field's inhomogeneity.^[23] This local equivalence thus bridges Newtonian intuition to relativistic geodesics, emphasizing that gravity dictates the "straightest" paths in curved spacetime.

Using the Variational Principle

The variational principle offers a systematic way to derive the geodesic equation by treating the worldline of a test particle as the path that extremizes the action functional in curved spacetime, rooted in the geometry of the metric tensor. For a massive test particle, the action is given by

S = -m \int \sqrt{-g_{\mu\nu} \, dx^\mu \, dx^\nu},

where m is the particle's rest mass, g_{\mu\nu} is the metric tensor, and the integral is taken along the path in spacetime with signature (-+++). This action corresponds to extremizing the proper time interval between events, as the proper time \tau satisfies d\tau^2 = -g_{\mu\nu} \, dx^\mu \, dx^\nu. For massless particles or light rays following null geodesics, the action simplifies to

S = \int g_{\mu\nu} \, dx^\mu \, dx^\nu,

since the proper time vanishes, and the path extremizes the null interval instead. These forms ensure reparameterization invariance, meaning the equations of motion are independent of the choice of affine parameter along the curve. To obtain the equations of motion, the action is varied with respect to the path x^\mu(\lambda), where \lambda is an arbitrary parameter. Equivalently, one employs the Lagrangian formulation with L = \frac{1}{2} g_{\mu\nu} \dot{x}^\mu \dot{x}^\nu, where the dot denotes differentiation with respect to the proper time \tau for timelike paths (or an affine parameter for null paths). The Euler-Lagrange equations then read

\frac{d}{d\tau} \left( \frac{\partial L}{\partial \dot{x}^\mu} \right) = \frac{\partial L}{\partial x^\mu}.

Computing the partial derivatives yields \frac{\partial L}{\partial \dot{x}^\mu} = g_{\mu\sigma} \dot{x}^\sigma and \frac{\partial L}{\partial x^\mu} = \frac{1}{2} \partial_\mu g_{\alpha\beta} \dot{x}^\alpha \dot{x}^\beta, so

\frac{d}{d\tau} (g_{\mu\sigma} \dot{x}^\sigma) = \frac{1}{2} \partial_\mu g_{\alpha\beta} \dot{x}^\alpha \dot{x}^\beta.

This equation can be rearranged using the compatibility of the Levi-Civita connection with the metric (\nabla_\rho g_{\mu\nu} = 0) and the definition of the Christoffel symbols \Gamma^\mu_{\alpha\beta} = \frac{1}{2} g^{\mu\rho} (\partial_\alpha g_{\rho\beta} + \partial_\beta g_{\rho\alpha} - \partial_\rho g_{\alpha\beta}). After lowering indices and contracting appropriately, followed by reparameterization to an affine parameter, the standard geodesic equation emerges:

\ddot{x}^\mu + \Gamma^\mu_{\alpha\beta} \dot{x}^\alpha \dot{x}^\beta = 0,

where the overdot now denotes differentiation with respect to the affine parameter. This test-particle action connects to the broader framework of general relativity through the Hilbert-Einstein action, which governs the gravitational field via S_\text{EH} = \frac{c^4}{16\pi G} \int \sqrt{-g} \, R \, d^4x, where R is the Ricci scalar. The full theory includes a matter action S_m, and for dilute test particles, the geodesic equation arises in the limit where the particle action S contributes negligibly to the field equations but determines the particle's motion independently. In this matterless limit for the background, free particles trace geodesics of the vacuum metric solutions. The variational approach has key advantages, including its ability to incorporate boundary conditions such as fixed endpoints for solving boundary-value problems in geodesic paths, and its extension to Hamiltonian mechanics in phase space, which facilitates treatments like quantization of geodesic motion or analysis in constrained systems.

Via Autoparallel Transport

In general relativity, the concept of parallel transport provides a geometric foundation for deriving the geodesic equation by considering how vectors are moved along curves without changing their direction relative to the spacetime manifold. Parallel transport of a vector V^\mu along a curve parameterized by \lambda is defined such that the covariant derivative along the tangent vector u = \frac{dx}{d\lambda} vanishes: \nabla_u V = 0, or in components, \frac{D V^\mu}{d\lambda} = u^\nu \nabla_\nu V^\mu = 0.^[25] This ensures that the vector is transported in a way that is consistent with the affine connection on the manifold, preserving its intrinsic properties under infinitesimal displacements along the curve.^[25] An autoparallel curve arises naturally from this transport when the curve's own tangent vector u^\mu = \frac{dx^\mu}{d\lambda} is parallel transported along itself, satisfying \nabla_u u = 0. This condition directly yields the geodesic equation in its covariant form: \frac{D u^\mu}{d\lambda} = 0.^[25] In the context of general relativity, geodesics are precisely these autoparallel curves with respect to the Levi-Civita connection, which is uniquely determined by being torsion-free and metric-compatible. The Levi-Civita connection \nabla is defined such that its Christoffel symbols \Gamma^\mu_{\nu\sigma} are symmetric in the lower indices (\Gamma^\mu_{\nu\sigma} = \Gamma^\mu_{\sigma\nu}), ensuring zero torsion, and it satisfies \nabla_\rho g_{\mu\nu} = 0 to preserve the metric g_{\mu\nu}. This torsion-free property guarantees that autoparallel curves coincide with metric geodesics, unlike in general affine connections where torsion can cause them to differ, leading to distinct notions of "straightest" paths. The metric compatibility ensures unique geodesics that extremize proper time or length, central to the description of free-falling particles in curved spacetime.^[25] To derive the coordinate form of the geodesic equation, expand the covariant derivative acting on the tangent vector: u^\nu \nabla_\nu u^\mu = u^\nu \left( \partial_\nu u^\mu + \Gamma^\mu_{\nu\sigma} u^\sigma \right) = 0. Assuming an affine parameterization where u^\nu u_\nu is constant, this simplifies to the second-order differential equation \frac{d^2 x^\mu}{d\lambda^2} + \Gamma^\mu_{\nu\sigma} \frac{dx^\nu}{d\lambda} \frac{dx^\sigma}{d\lambda} = 0, which governs the motion along geodesics.^[25] Geometrically, geodesics represent the "straightest" possible curves in curved spacetime, as they are the paths along which the direction of the tangent vector remains unchanged under parallel transport, generalizing the notion of straight lines in flat space while accounting for the manifold's intrinsic geometry.^[25]

Physical Interpretations and Properties

Geodesics as Paths of Stationary Proper Time

In general relativity, the proper time \tau experienced by an observer or test particle along a timelike worldline is given by the integral

\tau = \int \sqrt{ - g_{\mu\nu} \, dx^\mu \, dx^\nu },

where g_{\mu\nu} is the metric tensor describing the spacetime geometry, and the integral is taken along the path connecting two events, with the signature convention (-, +, +, +). Timelike geodesics represent the paths that extremize this proper time interval between fixed endpoints, specifically maximizing it in locally inertial frames, as derived from the variational principle applied to the worldline action. This extremization ensures that freely falling particles follow trajectories where no nearby path yields a longer proper time, reflecting the inertial nature of geodesic motion in curved spacetime. This property finds a direct analogy in the twin paradox, where one twin follows a geodesic path (e.g., inertial motion in flat spacetime or free fall in curved), accumulating the maximum proper time, while the other takes an accelerated, non-geodesic route, experiencing less proper time due to the shorter integrated interval along their worldline. In general relativity, this extends to scenarios like orbital motion around massive bodies, where the geodesic twin's clock advances more than that of a twin on a deviating path, underscoring how curvature alters the relative aging without violating causality.^[26] For null geodesics, which describe the propagation of light rays, the proper time vanishes, but the paths extremize the optical path length, analogous to Fermat's principle of least time in optics adapted to curved spacetime.^[27] Light follows null geodesics that make the arrival time stationary with respect to infinitesimal variations, leading to phenomena like gravitational lensing where rays bend to minimize effective travel time in the gravitational field. The stability of these extremal paths is analyzed through the second variation of the proper time functional, which determines whether a geodesic is a local maximum or minimum; the presence of conjugate points—where nearby geodesics intersect—signals instability and focusing due to spacetime curvature.^[28] In black hole spacetimes, such as Schwarzschild geometry, conjugate points along geodesics indicate regions of strong focusing, particularly near event horizons where null geodesics are trapped and timelike ones are drawn inexorably inward.^[29] Observationally, this maximization of proper time along geodesics is verified in systems like the Global Positioning System (GPS), where satellite clocks follow near-geodesic orbits and accumulate proper time at rates predicted by general relativity, ticking faster by about 38 microseconds per day relative to ground clocks due to weaker gravitational potential, despite special relativistic slowing from velocity.^[30] These corrections, accounting for the extremal proper time along satellite worldlines, ensure positional accuracy to within meters, demonstrating the practical impact of geodesic principles in navigation.^[30]

Relation to Field Equations in Vacuum

In vacuum regions of spacetime, where there is no matter or energy present, the Einstein field equations simplify to G_{\mu\nu} = 0, where G_{\mu\nu} = R_{\mu\nu} - \frac{1}{2} R g_{\mu\nu} is the Einstein tensor, R_{\mu\nu} is the Ricci tensor, R = g^{\mu\nu} R_{\mu\nu} is the Ricci scalar, and g_{\mu\nu} is the metric tensor. Taking the trace of this equation yields R = 0, so it is equivalent to the Ricci-flat condition R_{\mu\nu} = 0. The Ricci tensor itself is a contraction of the Riemann curvature tensor, given explicitly by

R_{\mu\nu} = \partial_{\lambda} \Gamma^{\lambda}_{\mu\nu} - \partial_{\nu} \Gamma^{\lambda}_{\mu\lambda} + \Gamma^{\lambda}_{\sigma\lambda} \Gamma^{\sigma}_{\mu\nu} - \Gamma^{\lambda}_{\sigma\nu} \Gamma^{\sigma}_{\mu\lambda},

where \Gamma^{\lambda}_{\mu\nu} are the Christoffel symbols determined by the metric.^[31] The second Bianchi identity, \nabla_{\lambda} R^{\rho}{}_{\sigma\mu\nu} + \nabla_{\mu} R^{\rho}{}_{\sigma\nu\lambda} + \nabla_{\nu} R^{\rho}{}_{\sigma\lambda\mu} = 0, when appropriately contracted (first on \rho and \nu, then on \mu and \lambda), implies the covariant divergence-free condition \nabla^{\mu} G_{\mu\nu} = 0. In the full theory with matter, the field equations G_{\mu\nu} = 8\pi T_{\mu\nu} (in units where G = c = 1) then enforce the covariant conservation law \nabla^{\mu} T_{\mu\nu} = 0 for the stress-energy tensor T_{\mu\nu}. In vacuum, where T_{\mu\nu} = 0, this conservation is automatically satisfied, but it plays a crucial role in the test particle limit.^[32] For a test particle of negligible mass (such that it does not source significant curvature), the stress-energy tensor is localized along its worldline and takes the form T^{\mu\nu} = m \int u^{\mu} u^{\nu} \frac{\delta^{4}(x - z(\tau))}{\sqrt{-g}} \, d\tau, where m is the rest mass, u^{\mu} = dz^{\mu}/d\tau is the four-velocity (with u_{\mu} u^{\mu} = -1), z^{\mu}(\tau) parameterizes the worldline by proper time \tau, and the integral is over the worldline. Outside the worldline, T_{\mu\nu} = 0, so the metric satisfies the vacuum equations R_{\mu\nu} = 0. The conservation law \nabla_{\mu} T^{\mu\nu} = 0 must hold globally; integrating it over a small tube surrounding the worldline and applying Stokes' theorem yields the equation of motion m \frac{D u^{\nu}}{d\tau} = 0, or equivalently \frac{D u^{\nu}}{d\tau} = u^{\mu} \nabla_{\mu} u^{\nu} = 0, which is the geodesic equation. This shows that the geodesic equation emerges directly as the compatibility condition for the particle's motion with the vacuum field equations.^[33] A related consequence is that, for any timelike vector u^{\mu} (such as the tangent to a geodesic), the vacuum condition implies R_{\mu\nu} u^{\mu} u^{\nu} = 0. Since R_{\mu\nu} = R^{\lambda}{}_{\mu\lambda\nu}, this contraction measures the Ricci curvature along the direction of u^{\mu}, and its vanishing ensures no net focusing or defocusing effect from the Ricci part of the curvature for a congruence of geodesics with that tangent; for a single geodesic, it is consistent with the autoparallel condition but follows from the overall vacuum geometry.^[32] Thus, in empty spacetime governed by the vacuum Einstein equations, the paths of test particles are precisely the geodesics of the Ricci-flat metric, requiring no separate equations of motion beyond the field equations themselves. This geometric interpretation underscores the equivalence principle, where gravity is encoded solely in the spacetime curvature.^[33]

Extensions and Applications

Geodesics for Charged Particles

In general relativity, the motion of uncharged test particles follows geodesics, but for charged particles, electromagnetic interactions modify this path, leading to a hybrid equation that incorporates both gravitational and Lorentz forces. The equation of motion for a charged particle of mass m and charge q is given by the covariant form

\frac{D u^\mu}{d\tau} = \frac{q}{m} F^\mu{}_\nu u^\nu,

where u^\mu = dx^\mu / d\tau is the four-velocity, \tau is the proper time, \frac{D}{d\tau} denotes the covariant derivative along the worldline, and F^\mu{}_\nu is the mixed electromagnetic field strength tensor. This equation generalizes the force-free geodesic equation by adding the electromagnetic contribution on the right-hand side, ensuring the path extremizes an action that includes both gravitational and electromagnetic terms. The derivation proceeds from the action principle, where the total action for the charged particle is

S = -m \int d\tau + q \int A_\mu \, dx^\mu.

Here, the first term -m \int d\tau = -m \int \sqrt{-g_{\alpha\beta} \, dx^\alpha dx^\beta} accounts for the gravitational interaction via the metric g_{\mu\nu}, while the second term represents the minimal coupling to the electromagnetic four-potential A_\mu. Varying this action with respect to the worldline x^\mu(\tau) and applying the Euler-Lagrange equations yields the modified equation of motion above, combining the geodesic deviation due to curvature with the Lorentz force. The electromagnetic field strength is defined by the Faraday tensor

F_{\mu\nu} = \partial_\mu A_\nu - \partial_\nu A_\mu,

which is antisymmetric and transforms as a tensor under general coordinate transformations; in curved spacetime, A_\mu couples minimally to the metric, ensuring gauge invariance under A_\mu \to A_\mu + \partial_\mu \Lambda. The raised index form F^\mu{}_\nu = g^{\mu\lambda} F_{\lambda\nu} appears in the force term, capturing both electric and magnetic effects in the observer's frame. In the weak-field limit, where the metric is nearly Minkowski (g_{\mu\nu} \approx \eta_{\mu\nu}) and curvature effects are negligible, this equation reduces to the flat-spacetime Lorentz force law \frac{d \mathbf{p}}{dt} = q (\mathbf{E} + \mathbf{v} \times \mathbf{B}), with \mathbf{p} = \gamma m \mathbf{v} the relativistic momentum. Conversely, in strong gravitational fields, such as near a charged black hole described by the Reissner-Nordström metric, the paths deviate significantly from pure geodesics; the electromagnetic term can cause spiraling orbits or altered accretion dynamics for charged particles, influenced by both the q^2/r^2 correction in the metric and the direct F^\mu{}_\nu u^\nu force. For consistency, the electromagnetic field must satisfy Maxwell's equations in curved spacetime, given covariantly by

\nabla_\mu F^{\mu\nu} = 4\pi J^\nu, \quad \nabla_{[\mu} F_{\nu\lambda]} = 0,

where J^\nu is the four-current, \nabla_\mu is the covariant derivative compatible with the metric, and the equations are in Gaussian units; the homogeneous equation follows from the definition of F_{\mu\nu}, while the inhomogeneous one sources the field via charges and currents. This framework ensures the overall dynamics respect both general covariance and the principles of electrodynamics in gravitational fields.

Examples in Curved Spacetime

In the Schwarzschild metric, which describes the spacetime geometry around a spherically symmetric, non-rotating mass M, the line element is given by

ds^2 = -\left(1 - \frac{2M}{r}\right) dt^2 + \left(1 - \frac{2M}{r}\right)^{-1} dr^2 + r^2 d\Omega^2,

where d\Omega^2 = d\theta^2 + \sin^2\theta \, d\phi^2 and units are such that G = c = 1^[34]. For timelike geodesics corresponding to orbital motion of massive particles, the equations of motion can be reduced using conserved quantities from the metric's symmetries, leading to an effective potential formulation for the radial motion: V_\text{eff} = -\frac{M}{r} + \frac{L^2}{2r^2} - \frac{M L^2}{r^3}, where L is the specific angular momentum. This potential allows analysis of bound orbits, such as nearly circular paths perturbed from Newtonian ellipses, with the innermost stable circular orbit at r = 6M marking the boundary beyond which orbits plunge into the central singularity. For null geodesics in the Schwarzschild metric, which trace the paths of photons, unstable circular orbits exist at the photon sphere radius r = 3M, where light can orbit the mass indefinitely before spiraling in or out^[35]. The weak-field deflection of light rays passing at large impact parameter b \gg M yields an angular deflection \delta \approx \frac{4M}{b}, twice the Newtonian prediction and a direct consequence of spacetime curvature. In the Friedmann-Lemaître-Robertson-Walker (FLRW) metric, modeling a homogeneous and isotropic expanding universe,

ds^2 = -dt^2 + a(t)^2 \left[ dr^2 + r^2 d\Omega^2 \right]

for flat spatial sections (with scale factor a(t)), radial timelike geodesics followed by comoving observers experience cosmological redshift z = \frac{a(t_0)}{a(t_e)} - 1, where t_0 and t_e are observation and emission times, respectively. Along these paths, the Hubble flow manifests as recession velocities v = H(t) d, with H(t) = \dot{a}/a the Hubble parameter and d the proper distance, illustrating how expansion stretches wavelengths without peculiar motion for distant galaxies^[36]. The geodesic deviation equation quantifies tidal effects in curved spacetime, describing the relative acceleration \xi^\mu between nearby geodesics with tangent vector u^\nu:

\frac{D^2 \xi^\mu}{d\tau^2} = - R^\mu{}_{\nu\rho\sigma} u^\nu \xi^\rho u^\sigma,

where R^\mu{}_{\nu\rho\sigma} is the Riemann curvature tensor and \tau the affine parameter; this leads to stretching along one direction and compression in others for freely falling observers, as in the tidal disruption near black holes or in cosmological voids. Numerical integration of geodesic equations is essential for simulating gravitational lensing, where ray-tracing algorithms solve the null geodesic ODEs in complex metrics to map source positions to observed images, accounting for multiple paths and magnification in strong-field regimes like black hole shadows^[37].

References

[1]
Lecture Notes on General Relativity - S. Carroll
The primary usefulness of geodesics in general relativity is that they are the paths followed by unaccelerated particles.
[2]
[PDF] 1. INTRODUCTION
1. INTRODUCTION. The geodesic hypothesis of general relativity asserts simply that A free particle follows a geodesic in spacetime. (where “a free particle” is ...
[3]
[PDF] General Relativity Fall 2019 Lecture 3: the geodesic equation
Sep 10, 2019 · This represents the equation of motion for a freely-moving particle in an arbitrary, non-inertial coordinate system. What it fundamentally means ...
[4]
[2503.09272] General Relativity and Geodesy - arXiv
Mar 12, 2025 · In this paper, we highlight the impact of General Relativity (GR) on geodesy: it provides corrections essential for the interpretation of high-precision ...
[5]
None
Below is a merged summary of the segments on "Geodesics in General Relativity" from the source https://www.damtp.cam.ac.uk/user/tong/gr/gr.pdf. To retain all information in a dense and organized manner, I will use a table in CSV format, followed by a concise narrative summary that integrates the key points. The table captures the details across all segments, including page references, definitions, roles, equations, and parameterization details.
[6]
[PDF] Unified Framework for Geodesic Dynamics with Conservative ... - arXiv
Sep 18, 2025 · In differential geometry and general relativity, a geodesic represents the extremal path between two space- time events, generalizing the notion ...
[7]
[PDF] Geodesics in Quantum Gravity - arXiv
A cornerstone of classical General Relativity (GR) is the geodesic equation, which governs the free-fall motion of test particles in curved spacetime.
[8]
1854: Riemann's classic lecture on curved space
Jun 1, 2013 · The geometric foundation for his work was laid some 60 years earlier, with the work of a German mathematician named Georg Friedrich Bernhard ...
[9]
On the Hypotheses Which Lie at the Bases of Geometry
In his inaugural lecture at Göttingen in 1854, Bernhard Riemann famously observed that "in a discrete manifold, the ground of its metric relations is given in ...
[10]
[PDF] On the history of Levi-Civita's parallel transport - arXiv
Aug 6, 2016 · This had never been considered before Levi-Civita's work, which therefore became rightly a milestone of differential geometry. Therefore, we ...
[11]
Schouten, Levi-Civita and the notion of parallelism in Riemannian ...
As is well known, Levi-Civita introduced his generalization of parallelism for contravariant vectors by considering the (local) embedding of a Riemannian ...
[12]
[PDF] EINSTEIN'S FIRST SYSTEMATIC EXPOSITION OF GENERAL ...
This was the first of many papers on general relativity in the Sit- zungsberichte, including such gems as (Einstein 1916c) on the relation between invariance ...
[13]
IX. A determination of the deflection of light by the sun's gravitational ...
The purpose of the expeditions was to determine what effect, if any, is produced by a gravitational field on the path of a ray of light traversing it.
[14]
Relativity: the general theory : Synge, J. L. (John Lighton), 1897
Jan 28, 2019 · Relativity: the general theory ; Publication date: 1964 ; Topics: Relativity (Physics) ; Publisher: Amsterdam, North-Holland Pub. Co. [1964].
[15]
[PDF] Lecture Notes on General Relativity - arXiv
them Christoffel symbols, but we won't use the funny notation. The study of ... This is the geodesic equation, another one which you should memorize.
[16]
[PDF] General Relativity - DAMTP - University of Cambridge
This is the relativistic version of the geodesic equation (1.6). We ... Proof: We start by showing uniqueness. Suppose that such a connection exists ...
[17]
[PDF] carroll-gr-textbook.pdf - Caltech Astronomy
Dec 3, 1997 · There are also null geodesics, which satisfy the same equation, except that the proper time cannot be used as a parameter (some set of allowed.
[18]
[PDF] Efficient computation of null affine parameters - arXiv
Dec 20, 2023 · It is only for the so-called affine parameterizations that the geodesic equation specifically reduces to Vkk = 0.Missing: reparameterization | Show results with:reparameterization
[19]
[PDF] Non-affinely parametrized geodesic equation - UMD Physics
The non-affinely parametrized geodesic equation is d/dλ (gαν ˙xν) − 1/gµν,α ˙xµ ˙xν = f(λ)gαν ˙xν, where f(λ) is a nonzero function.Missing: relativity | Show results with:relativity
[20]
Albert Einstein's 1916 Review Article on General Relativity - arXiv
May 13, 2004 · View PDF. Abstract: The first comprehensive overview of the final version of the general theory of relativity was published by Einstein in 1916 ...
[21]
[PDF] Geometrodynamics based on geodesic equation with Cartan ... - arXiv
Feb 9, 2020 · The most prominent application is in general relativity [4–7] that a geodesic generalizes the notion of a straight line to curved spacetime.
[22]
Notion of Parallelism on a Generic Manifold and Consequent ... - arXiv
Oct 20, 2022 · This paper written in 1916 by T. Levi-Civita introduces the notion of parallelism for a Riemannian manifold of arbitrary dimensions.
[23]
[PDF] The twin paradox and the principle of relativity - arXiv
Timelike geodesics have in general extremal proper time between two events, but the proper time along a geodesic curve can be either maximal or minimal. 5 ...
[24]
A variational principle for time of arrival of null geodesics - gr-qc - arXiv
Oct 20, 1998 · ... Fermat's Principle of Least Time which yields paths of extremal time. We here investigate a related but different problem. We consider an ...
[25]
[PDF] arXiv:2303.06788v1 [gr-qc] 13 Mar 2023
Mar 13, 2023 · The second variation of the length functional, used for timelike geodesics is called the index form, from which the method takes its name ...
[26]
[PDF] The singularity theorems of General Relativity and their low ... - arXiv
Nov 14, 2022 · Abstract. On the occasion of Sir Roger Penrose's 2020 Nobel Prize in Physics, we review the singularity theorems of General Relativity, ...
[27]
Relativity in the Global Positioning System - PMC - NIH
The Global Positioning System (GPS) uses accurate, stable atomic clocks in satellites and on the ground to provide world-wide position and time ...Missing: maximization | Show results with:maximization
[28]
4 The Einstein Equations‣ General Relativity by David Tong
As a quick sanity check, recall that the Ricci tensor takes the schematic form (3.131) R ∼ ∂ ⁡ Γ + Γ ⁢ Γ while the Levi-Civita connection itself is Γ ∼ ∂ ⁡ g . ...
[29]
[PDF] Lecture VII: Gravity gradients, the Ricci tensor, and the field equations
Oct 29, 2012 · This lecture will complete our description of curvature tensors. In particular, we will build the Einstein tensor Gµν.
[30]
https://pmc.ncbi.nlm.nih.gov/articles/PMC5253894/
[31]
Metamorphoses of a photon sphere | Phys. Rev. D
Oct 25, 2017 · There are circular planar null geodesics at r = 3 ⁢ M around a Schwarzschild black hole of mass M . These geodesics form a photon sphere.Missing: original | Show results with:original
[32]
Joining the Hubble flow: implications for expanding space
We conclude that particles in eternally expanding but otherwise arbitrary universes do not in general rejoin the Hubble flow.Joining The Hubble Flow... · 4 Joining The Hubble Flow · 4.2 Comparing The...<|separator|>
[33]
Numerical simulations of the light propagation in the gravitational ...
This paper describes numerical simulations of the light propagation in the gravitational field of moving gravitating bodies and summarizes the underlying the-.