Einstein Field Equations

The Einstein field equations (EFE) are a set of 16 coupled, hyperbolic-elliptic, nonlinear partial differential equations at the heart of general relativity. They describe how mass and energy curve the fabric of spacetime, and that curvature is what you experience as gravity.

At a high level, the EFE relate the geometry of spacetime (the left side of the equation) to the distribution of matter and energy within it (the right side). Because the Einstein tensor G_μν and the stress-energy tensor T_μν are both symmetric, the 16 equations reduce to 10 independent ones. Apply the four Bianchi identities and you’re left with six truly independent equations, giving the metric four gauge-fixing degrees of freedom.

Albert Einstein first published these equations in 1915 as a tensor equation connecting local spacetime curvature (the Einstein tensor) with the stress, momentum, and energy present in that region (the stress-energy tensor). In practical terms, they determine the metric tensor of spacetime for a given arrangement of stress-energy-momentum. Think of it as the gravitational analog of Maxwell’s equations, which relate electromagnetic fields to the distribution of charges and currents.

Because the Einstein tensor depends on the metric tensor, you can rewrite the EFE as a set of nonlinear partial differential equations whose solutions are the components of the metric tensor itself. Once you have the metric, you can trace out the paths that particles and light follow through curved spacetime using the geodesic equation.

Mathematical form of the Einstein field equations

The Einstein field equations are written as:

G_μν + Λg_μν = κT_μν

Here, G_μν is the Einstein tensor, Λ is the cosmological constant, g_μν is the metric tensor, κ is the Einstein gravitational constant, and T_μν is the stress-energy tensor. Each symbol carries precise geometric or physical meaning, and understanding what they represent is the key to reading these equations.

You can expand the Einstein tensor as:

G_μν = R_μν – ½Rg_μν

Here, R_μν is the Ricci curvature tensor and R is the scalar curvature. The Einstein tensor is a symmetric second-rank tensor that depends only on the metric tensor and its first and second derivatives. That constraint is what makes it the natural geometric object to appear on the left side of the field equations.

The Einstein gravitational constant is defined as:

κ = 8πG/c⁴ ≈ 2.077 × 10^-43 N^-1

Here, G is the Newtonian constant of gravitation and c is the speed of light in vacuum. Notice how tiny κ is. That’s why you need enormous amounts of mass-energy to produce any noticeable curvature of spacetime.

Substituting the definition of G_μν, the full expanded form reads:

R_μν – ½Rg_μν + Λg_μν = κT_μν

Every term on the left side has units of 1/length². The left side encodes the curvature of spacetime as determined by the metric, while the right side represents the stress-energy-momentum content. You can read the equation as a precise statement: “matter tells spacetime how to curve.”

Together with the geodesic equation, the EFE form the mathematical backbone of general relativity. Each side of the equation is a symmetric 4×4 tensor with 10 independent components. After applying the four Bianchi identities, six independent equations remain, leaving the metric with four gauge-fixing degrees of freedom.

Although Einstein originally wrote these equations in four dimensions, theorists have explored their consequences in arbitrary n dimensions. When you write the equations out in full component form, you see that they are a system of 10 coupled, nonlinear, hyperbolic-elliptic partial differential equations. That complexity is exactly why finding exact solutions is so challenging.

Equivalent formulations

If you take the trace of both sides of the EFE with respect to the metric, you get:

R – R + DΛ = κT

Here, D is the spacetime dimension. Solving for R and substituting back into the original EFE gives you the “trace-reversed” form:

R_μν – Λg_μν = κ(T_μν – T_μν)

In D = 4 dimensions, this simplifies to:

R_μν – Λg_μν = κ(T_μν – ½Tg_μν)

You can reverse the trace again to recover the original EFE. In certain problems, the trace-reversed form turns out to be far more convenient to work with, particularly when you want to isolate the Ricci tensor on one side.

Linearized form

The nonlinearity of the EFE makes finding exact solutions difficult, but you can simplify things with an approximation. Far from a source of gravitating matter, the gravitational field becomes extremely weak and spacetime closely resembles flat Minkowski space.

In this regime, you write the metric as the sum of the Minkowski metric and a small perturbation term, then drop higher-order powers of the perturbation. This process of linearization is exactly how physicists study gravitational radiation, the ripples in spacetime that LIGO and Virgo have now famously detected.

Polynomial form

The standard EFE contain the inverse of the metric tensor, but you can rearrange them into a purely polynomial form. To do this, first express the determinant of the metric in four dimensions using the Levi-Civita symbol:

det(g) = ε^αβγδε^κλμνg_ακg_βλg_γμg_δν

Then the inverse metric can be written as:

g^ακ = (ε^αβγδε^κλμνg_ακg_βλg_γμg_δν) / det(g)

Substitute this into the field equations and multiply both sides by a suitable power of det(g) to clear the denominator. What you get is a set of polynomial equations in the metric tensor and its first and second derivatives. You can also rewrite the Einstein-Hilbert action in polynomial form using appropriate field redefinitions.

The cosmological constant

The term Λg_μν in the EFE wasn’t part of Einstein’s original 1915 publication. He added the cosmological constant Λ later to allow for a static universe, one that neither expands nor contracts. That attempt failed for two reasons:

• The static solution it produces is unstable. Any small perturbation causes the universe to expand or collapse.
• Edwin Hubble’s observations in the late 1920s showed that the universe is actually expanding.

Einstein reportedly called the cosmological constant the “biggest blunder” of his life and dropped it from his equations. However, including it doesn’t introduce any mathematical inconsistency, and for decades most physicists simply assumed Λ = 0.

That changed in 1998 when observations of distant Type Ia supernovae revealed that the expansion of the universe is accelerating. Explaining this acceleration requires a small, positive value of Λ. Its effect is negligible at the scale of galaxies or smaller, which is why you don’t notice it in everyday physics, but it dominates the large-scale dynamics of the cosmos.

Solutions of the Einstein field equations

Solutions to the EFE are metrics of spacetime. They describe the geometry of spacetime and, through the geodesic equation, the inertial motion of objects within it. Because the equations are nonlinear, finding complete solutions is rarely straightforward. For instance, no exact solution exists for a spacetime containing two massive bodies, such as a binary star system.

In practice, physicists make post-Newtonian approximations to handle such cases. Still, many important exact solutions do exist, and studying them has been one of the most productive activities in theoretical physics. Exact solutions have led to the prediction of black holes, the development of cosmological models for the evolution of the universe, and the theoretical framework for gravitational waves.

The method of orthonormal frames, introduced by Ellis and MacCallum, offers another route to new solutions. This approach reduces the EFE to a set of coupled, nonlinear ordinary differential equations. Hsu and Wainwright showed that self-similar solutions correspond to fixed points of the resulting dynamical system, and researchers like LeBlanc, Kohli, and Haslam have used these techniques to discover entirely new solutions.

Frequently asked questions