Fact-checked by Grok 2 weeks ago

Einstein field equations

The Einstein field equations (EFE), also known as the Einstein equations, are a set of ten nonlinear partial differential equations in the theory of that relate the local to the local , , and stress within that . First derived and published by on November 25, 1915, in the Sitzungsberichte der Preußischen Akademie der Wissenschaften, they represent the gravitational analog of in Newtonian gravity and form the cornerstone of by expressing how mass and dictate the geometry of . In their standard form, the EFE are written as
G_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu},
where G_{\mu\nu} is the (derived from the tensor and scalar), g_{\mu\nu} is the describing , T_{\mu\nu} is the stress-energy tensor encoding the of and , G is the , c is the , and \Lambda is the (introduced by Einstein in to permit static cosmological models but later revised). This formulation generalizes the , treating gravity not as a force but as the curvature of four-dimensional caused by energy-momentum.
The EFE were the culmination of Einstein's eight-year quest for a relativistic theory of , building on earlier work with and overcoming challenges like ensuring . Their successful prediction of the anomalous of Mercury's perihelion, announced by Einstein in 1915, provided immediate empirical validation, while the 1919 observations confirming light deflection by the Sun's cemented their acceptance. Solutions to the EFE underpin key phenomena in , including black holes (e.g., the Schwarzschild solution), (detected in 2015), the in cosmology, and the large-scale structure of the cosmos. Despite their complexity—solving them exactly remains challenging and often requires approximations or numerical methods—the EFE remain the most accurate description of gravitational interactions on cosmic scales.

Historical Development

Origins in General Relativity

The development of the Einstein field equations is rooted in Albert Einstein's quest to extend to include , beginning with the formulation of the in 1907. While reviewing the status of in a survey article for the Jahrbuch der Radioaktivität und Elektronik, Einstein realized that a person in experiences no gravitational force, equating the effects of with in a local frame. This insight, later formalized as the , provided the cornerstone for generalizing to accelerated frames and curved , resolving the incompatibility between Newtonian and special relativity's constant . By 1911, Einstein had predicted the deflection of light by , setting the stage for empirical tests, though the full theoretical framework remained elusive. A pivotal milestone occurred in 1912 upon Einstein's return to Zurich as a professor at the ETH, where he began exploring curved spacetime in private notes known as the Zurich Notebook. These calculations marked the first use of the metric tensor to describe variable geometry, introducing concepts like geodesics for particle motion and initial attempts at gravitational field equations, while grappling with the Newtonian limit. This work was deeply influenced by Riemannian geometry, which Marcel Grossmann introduced to Einstein, enabling the representation of gravity through spacetime curvature rather than forces. From 1913 to 1915, Einstein collaborated closely with Grossmann, his former classmate and mathematician, on the "Entwurf" (sketch) theory, published in 1913, which employed tensor calculus to formulate provisional field equations covariant under restricted coordinate transformations. Grossmann's guidance on differential geometry, including the absolute differential calculus developed by Gregorio Ricci-Curbastro and Tullio Levi-Civita, was crucial; Levi-Civita's concept of parallel transport and the torsion-free connection later helped Einstein refine the mathematical structure for general covariance. Einstein's core motivation was to geometrize , transforming Newtonian action-at-a-distance into the intrinsic of caused by and , thereby eliminating absolute space and aligning with of relational motion. This approach sought to unify inertia and gravitation dynamically, making geometry contingent on matter distribution. In parallel during , independently pursued a leading to similar field equations, with correspondence between the two scholars in November highlighting the competitive yet convergent efforts that spurred Einstein to achieve full . By late , these foundations culminated in the completion of , with the field equations expressing the precise relation between and energy-momentum.

Publication and Initial Reception

On November 25, 1915, presented his definitive form of the field equations to the in , in a paper titled "The Field Equations of Gravitation," marking the culmination of his work on . This presentation followed a series of four papers delivered to the academy that month, where Einstein had iteratively refined his gravitational theory, resolving earlier inconsistencies with the . The equations linked curvature to the distribution of mass and energy, providing a covariant framework that Einstein had sought since 1912. Nearly simultaneously, mathematician submitted a paper on November 20, 1915, to the Göttingen Academy, deriving mathematically equivalent field equations through a variational approach based on an action principle. This led to ongoing historical debates about , with analyses of Hilbert's proofs and Einstein's indicating that both arrived independently, though Einstein's physical interpretation preceded Hilbert's formal derivation. Hilbert's work, published in March 1916, acknowledged Einstein's contributions and focused on axiomatic foundations rather than empirical predictions. The equations received immediate empirical support through their successful prediction of the anomalous of Mercury's perihelion, which Einstein calculated as an advance of arcseconds per century, matching observations unexplained by Newtonian . This verification, detailed in Einstein's November 18, 1915, , generated excitement among astronomers and physicists, providing the first quantitative of the . However, early reactions were mixed: critics, including some traditional physicists, highlighted the equations' mathematical complexity and the apparent lack of uniqueness in selecting the specific tensor form from generally covariant possibilities, arguing it relied too heavily on aesthetic criteria over empirical necessity. Prominent endorsements helped counter these concerns; Hendrik Lorentz, a leading figure in relativity's development, praised the theory's elegance in private correspondence and reports, noting its consistency with while awaiting further tests. Figures like Max Abraham and Gustav Mie offered constructive critiques on uniqueness but acknowledged the framework's potential, fostering gradual acceptance within the . In 1916, Einstein solidified the equations' status with his comprehensive review article, "The Foundation of the General Theory of Relativity," published in , which systematically outlined the theory's principles, derivations, and implications, becoming a foundational text for subsequent .

Mathematical Formulation

The Core Equations

The Einstein field equations form the foundational set of ten coupled, nonlinear partial differential equations in that relate the geometry of four-dimensional to the distribution of mass and energy within it. These equations, first presented by in their final form on November 25, 1915, are expressed in tensor notation as G_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}, where Greek indices \mu, \nu run from 0 to 3, corresponding to the time and three spatial coordinates. The left-hand side involves the Einstein tensor G_{\mu\nu}, defined as G_{\mu\nu} = R_{\mu\nu} - \frac{1}{2} R g_{\mu\nu}, which captures the curvature of spacetime. Here, R_{\mu\nu} is the Ricci curvature tensor, obtained by contracting the Riemann curvature tensor R^\rho_{\ \sigma\mu\nu} via R_{\mu\nu} = R^\rho_{\ \mu\rho\nu}; R = g^{\mu\nu} R_{\mu\nu} is the Ricci scalar curvature, the trace of the Ricci tensor; and g_{\mu\nu} is the metric tensor that defines distances and angles in the spacetime manifold. The right-hand side features the stress-energy tensor T_{\mu\nu}, a that encodes the , flux, and stresses associated with , , and other forms in . The coupling constant \frac{8\pi [G](/page/Gravitational_constant)}{c^4} incorporates the Newtonian gravitational constant [G](/page/Gravitational_constant) \approx 6.67430 \times 10^{-11} \, \mathrm{m^3 \, kg^{-1} \, s^{-2}} and the c \approx 2.99792 \times 10^8 \, \mathrm{m \, s^{-1}}, ensuring dimensional consistency between geometric and physical quantities. In theoretical work and numerical simulations, it is common to adopt geometrized units where G = c = 1, reducing the equations to the simplified form G_{\mu\nu} = 8\pi T_{\mu\nu} and measuring lengths in units of time (or ). This streamlines calculations while preserving the physical content. Conceptually, the equations embody of equivalence between inertial and gravitational mass by equating curvature (left side) to the local presence of energy-momentum (right side), such that massive bodies "tell how to curve" and curved "tells matter how to move" along geodesics.

Sign Conventions

The Einstein field equations are formulated within the framework of a Lorentzian metric on spacetime, where the signature convention determines the relative signs of the temporal and spatial components in the line element. The two predominant conventions are the mostly-minus signature (+ − − −), in which the Minkowski line element takes the form ds^2 = c^2 \, dt^2 - dx^2 - dy^2 - dz^2, and the mostly-plus signature (− + + +), given by ds^2 = -c^2 \, dt^2 + dx^2 + dy^2 + dz^2. These differ in the sign of the time component, with the former making timelike intervals positive and the latter making spacelike intervals positive. This choice impacts the signs appearing in the , particularly through adjustments in the computation of the , which is contracted from the . The Riemann tensor's definition often includes a sign flip between conventions (e.g., the term involving ∂_λ Γ^σ_μν may change sign), necessitating corresponding changes to preserve the equation's form, such as altering the overall sign on one side of G_{\mu\nu} = 8\pi G T_{\mu\nu}. Despite these adjustments, the physical implications remain invariant, as the conventions are equivalent up to an overall rescaling. Historically, the (+ − − −) signature has been more common in classical general relativity literature, aligning with Einstein's original presentations where positive proper time for timelike paths was emphasized, while the (− + + +) signature is preferred in particle physics contexts for compatibility with positive-definite kinetic terms in quantum field theories. This divergence arose from differing emphases: relativity on geometric positivity for spacelike separations versus particle physics on positive energies for massive particles. Consistency across the formulation is essential; switching signatures requires flipping signs in the stress-energy tensor T_{\mu\nu} (e.g., from positive in (+ − − −) to negative in (− + + +) without adjustment) to avoid altering physical predictions like attractive or the positive mass-energy condition. The does not affect observable outcomes, as it is a notational freedom, but uniform application within a given work prevents errors in derivations.

Alternative Expressions

The trace-reversed form of the Einstein field equations expresses the Ricci curvature tensor directly in terms of the stress-energy tensor, providing an alternative to the standard Einstein tensor formulation. This form is derived by taking the trace of the original equations and rearranging, yielding R_{\mu\nu} = \frac{8\pi G}{c^4} \left( T_{\mu\nu} - \frac{1}{2} T g_{\mu\nu} \right), where T = g^{\alpha\beta} T_{\alpha\beta} is the trace of the stress-energy tensor T_{\mu\nu}. This representation is particularly advantageous in contexts such as , where it simplifies the treatment of perturbations around flat , and in numerical simulations of gravitational systems. Contracting the trace-reversed equations with the g^{\mu\nu} further reduces them to a scalar relating the Ricci scalar R to the trace of the stress-energy tensor: R = -\frac{8\pi G}{c^4} T. This scalar form highlights the direct proportionality between curvature and the total energy-momentum trace, serving as a useful condition in derivations and approximations of solutions. In the Palatini formulation, the Einstein field equations emerge from varying an with respect to both the and an independent , rather than assuming the connection is the derived from the alone. The resulting equations enforce that the connection is metric-compatible and torsion-free, reproducing the standard field equations while allowing for extensions in modified gravity theories. This approach underscores the geometric structure of by treating the as a variable. Index-free notations offer compact ways to express the Einstein field equations without explicit coordinate indices, emphasizing their tensorial nature. In abstract index notation, the equations appear as G_{ab} = \frac{8\pi G}{c^4} T_{ab}, where a and b are abstract placeholders indicating the tensor type rather than specific components, facilitating manipulations in differential geometry./06%3A_Waves/6.07%3A_Abstract_Index_Notation) Alternatively, using differential forms, the equations can be formulated in terms of the curvature 2-form \Omega and the volume form, such as \epsilon(\Omega) = \frac{8\pi G}{c^4} *T, where * denotes the Hodge dual, providing a coordinate-independent perspective suited to exterior calculus methods.

Physical Interpretation

Curvature and Energy-Momentum

The Einstein tensor G_{\mu\nu}, which forms the geometric side of the field equations, encapsulates the of through its construction from the R^\rho_{\sigma\mu\nu}. The Riemann tensor quantifies the extent to which nearby geodesics deviate from each other, a directly linked to tidal forces acting on extended bodies in gravitational fields. Specifically, the geodesic deviation equation, \frac{D^2 \xi^\mu}{d\tau^2} = -R^\mu_{\nu\rho\sigma} u^\nu \xi^\rho u^\sigma, describes how the separation vector \xi^\mu between two infinitesimally close geodesics with tangent u^\mu evolves, with the Riemann tensor term representing the relative acceleration due to spacetime curvature. This tidal effect generalizes Newtonian gravitational forces, illustrating how curvature influences the motion of matter along free-fall paths. On the physical side, the equations relate this to the stress-energy tensor T_{\mu\nu}, which describes the distribution of energy, , and stress in the . For a —characterized by isotropic and no —the tensor takes the form T_{\mu\nu} = (\rho + p) u_\mu u_\nu + p g_{\mu\nu}, where \rho is the (with T^{00} \approx \rho c^2 in the , capturing mass-energy content), p is the , u^\mu is the , and g_{\mu\nu} is the . The spatial components T^{ij} encode flux and stresses, such as forces perpendicular to the flow. In this framework, the core relation equates G_{\mu\nu} to a constant multiple of T_{\mu\nu}, balancing geometric structure with material sources. This interplay embodies the foundational insight of : matter and energy dictate the curvature of , while the resulting curved geometry governs the trajectories of matter. As physicist succinctly phrased it, "matter tells how to curve; curved tells matter how to move." This bidirectional relationship unifies as an intrinsic feature of geometry rather than a force. The uniqueness of these equations among possible gravitational theories arises from their derivation via the applied to the Hilbert-Einstein , S = \frac{c^4}{16\pi G} \int R \sqrt{-g} \, d^4x + S_m, where R is the Ricci scalar (contracted from the Riemann tensor), g is the determinant, and S_m is the matter . Varying this with respect to the g^{\mu\nu} yields the field equations after and using the Bianchi identities to ensure . Specifically, the variation \delta S = 0 leads to \delta(\sqrt{-g} R) / \sqrt{-g} = - \frac{8\pi G}{c^4} T_{\mu\nu} \delta g^{\mu\nu}, producing G_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}. This approach, introduced by in 1915, confirms the equations as the unique diffeomorphism-invariant theory linear in the second derivatives of the and quadratic in the first derivatives, consistent with .

Role of the Cosmological Constant

The Einstein field equations can be modified to include a term \Lambda, resulting in the form G_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}, where \Lambda is a universal constant with dimensions of inverse length squared, added to the left-hand side as a geometric contribution proportional to the g_{\mu\nu}. This term was introduced by in 1917 to enable a static, closed model consistent with , where \Lambda > 0 balances gravitational attraction with a repulsive effect to maintain spatial stability against collapse. In this model, the positive \Lambda provides the necessary condition for a finite, eternal without or , reflecting Einstein's initial cosmological vision before observational of emerged. Einstein retracted the cosmological constant in 1931, deeming it superfluous after Alexander Friedmann's 1922 derivation of expanding universe solutions from the original field equations without , and Georges Lemaître's 1927 independent work confirming dynamic models supported by Edwin Hubble's observations of galactic recession. Friedmann demonstrated that homogeneous, isotropic spacetimes could evolve over time, challenging the need for a static and prompting Einstein to acknowledge the viability of expansion. Lemaître further developed this into a "primeval atom" hypothesis, aligning theory with emerging data on cosmic . The term remained marginalized for decades until 1998, when observations of distant Type Ia supernovae by the Supernova Cosmology Project (led by ) and the High-Z Supernova Search Team (led by and ) revealed an , best explained by reintroducing a positive as the dominant component in late-time cosmology. Physically, the cosmological constant represents a uniform vacuum energy density \rho_\Lambda = \frac{\Lambda c^2}{8\pi G}, contributing to the total budget without dilution during expansion, unlike or . This density exerts an effective p_\Lambda = -\rho_\Lambda c^2, driving repulsive that accelerates cosmic expansion, as captured by the equation of state parameter w = -1. In the standard field equations without \Lambda, relates solely to and via the stress-energy tensor, but the constant term introduces a baseline repulsion independent of local sources, influencing large-scale and the universe's ultimate fate. However, the presents a profound theoretical challenge known as the "." predicts that , arising from quantum fluctuations, should contribute a term to \rho_\Lambda approximately 120 orders of magnitude larger than the observed value inferred from \Lambda. This enormous discrepancy remains one of the most significant unsolved problems in physics as of 2025, prompting ongoing research into possible resolutions such as principles, modified theories, or dynamical models.

Stress-Energy Tensor

The stress-energy tensor T_{\mu\nu}, also known as the energy-momentum tensor, encodes the distribution of , , and associated with and fields within , acting as the primary source term on the right-hand side of the Einstein field equations. For classical fields, T_{\mu\nu} is symmetric and defined such that its components T_{00} correspond to , T_{0i} to momentum density, T_{i0} to momentum flux, and T_{ij} to . This tensor arises naturally from the variation of the matter S_m with respect to the via the Hilbert prescription: T_{\mu\nu} = -\frac{2}{\sqrt{-g}} \frac{\delta S_m}{\delta g^{\mu\nu}}, where g is the of the g_{\mu\nu}. Alternatively, it can be derived from applied to the invariance of the under , yielding the canonical form through spacetime translations, though the Hilbert form is preferred in for its symmetry and compatibility with the metric. In both approaches, the tensor satisfies the local conservation law \nabla^\mu T_{\mu\nu} = 0, reflecting the absence of external sources and the invariance of the full theory. A common realization of T_{\mu\nu} for macroscopic matter is the form, which assumes isotropic p and \rho in the fluid's , with no or conduction: T_{\mu\nu} = (\rho + p) u_\mu u_\nu + p g_{\mu\nu}, where u^\mu is the normalized such that u^\mu u_\mu = -1 (in the mostly-plus ). This expression is obtained by varying the appropriate fluid or averaging microscopic contributions. Special cases include , modeling pressureless matter like non-relativistic particles or collisionless gases, where p = 0, simplifying to T_{\mu\nu} = \rho u_\mu u_\nu. Another example is or relativistic particles, characterized by the equation of state p = \rho/3, as in gas or massless fields, where the isotropic equals one-third of the due to the traceless nature of the underlying stress contributions. For electromagnetic fields, the stress-energy tensor is derived from the Maxwell Lagrangian \mathcal{L} = -\frac{1}{4} F_{\mu\nu} F^{\mu\nu}, where F_{\mu\nu} is the Faraday tensor. Applying to Lorentz transformations or the Hilbert variation yields the canonical symmetric form: T_{\mu\nu} = F_{\mu\alpha} F_\nu{}^\alpha - \frac{1}{4} g_{\mu\nu} F_{\alpha\beta} F^{\alpha\beta}. This tensor is traceless (T^\mu{}_\mu = 0) and satisfies \nabla^\mu T_{\mu\nu} = -j^\mu F_{\mu\nu}, where j^\mu is the four-current, reducing to the in source-free regions. The derivation traces back to early applications of Noether's ideas, as explicitly computed for electrodynamics. In relativistic contexts, T_{\mu\nu} is normalized such that its components carry units of . In natural units where c = 1 and \hbar = 1, all components have dimensions of energy per volume (or inverse length to the fourth power), with \rho representing rest mass energy density divided by c^2 in the non-relativistic limit. In SI units, the energy-momentum components scale as joules per cubic meter (J/m³), while momentum fluxes involve kg/(m s²), ensuring consistency with the coupling \frac{8\pi [G](/page/G)}{c^4} T_{\mu\nu} in the field equations, where G and c restore dimensional balance. This normalization underscores the equivalence of inertial and gravitational mass, central to .

Key Properties

Energy-Momentum Conservation

The Einstein field equations ensure a conservation principle for the distribution of and in through the inherent symmetries of the . Specifically, the twice-contracted second Bianchi , a fundamental property of the , implies that the covariant divergence of the vanishes identically: \nabla^\lambda G_{\lambda\mu} = 0. This holds regardless of the presence of matter, as it arises purely from the metric compatibility of the and the antisymmetry of the Riemann tensor. Given the Einstein field equations G_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu} (in the absence of a ), contracting with the and taking the covariant divergence yields the corresponding for the stress-energy tensor: \nabla^\mu T_{\mu\nu} = 0. This equation is derived directly by substituting the field equations into the Bianchi , ensuring that the matter content satisfies the automatically without additional constraints. The stress-energy tensor T_{\mu\nu}, which encodes the , flux, and stresses of matter and fields, thus obeys this relation as a consequence of the invariance of the theory. The form \nabla^\mu T_{\mu\nu} = 0 represents local of energy-momentum: at every point in , the of the energy-momentum vanishes, meaning no net creation or destruction occurs locally, with flows balanced along timelike and spacelike geodesics. To obtain global statements, one integrates this over a spacelike or volume, applying the to relate the volume integral of the to the surface of T_{\mu\nu} through the boundary. However, such global integrals yield conserved quantities only in spacetimes with suitable symmetries, like a timelike , which allows defining a . A key limitation arises in curved spacetimes without such symmetries, particularly those undergoing dynamical , such as cosmological models. While the local law \nabla^\mu T_{\mu\nu} = 0 always holds, global energy is not absolutely because the expanding alters the interpretation of integrated quantities— for instance, redshifts without violating the local equation, as the "lost" energy is accounted for by the work done against gravitational . This reflects that in is tied to the local rather than a fixed background, precluding a universal global energy in time-dependent or non-asymptotically flat spacetimes.

Nonlinearity

The nonlinearity of the Einstein field equations arises fundamentally from the structure of the Ricci tensor, which is computed from the metric tensor and its derivatives; specifically, the Christoffel symbols involve inverse metric components multiplied by metric derivatives, and the Ricci tensor includes quadratic terms in these symbols, such as \Gamma^\sigma_{\nu\mu} \Gamma^\rho_{\rho\sigma}. This results in the left-hand side of the equations being a nonlinear functional of the metric, reflecting how spacetime curvature depends quadratically on the gravitational field itself. A key physical consequence of this nonlinearity is the failure of the , which holds in linear theories like for but not in . In the latter, gravitational fields do not simply add; instead, the field produced by one source curves in a way that influences the propagation of other fields, embodying the principle that "gravity gravitates." This self-sourcing effect manifests dramatically in the interactions of , where waves can scatter off each other due to their own energy-momentum, producing higher-order effects observable in strong-field events like mergers. For instance, during such mergers, nonlinearities generate the characteristic ringdown phase of the waveform, as confirmed by simulations and gravitational-wave detections. The nonlinear nature also poses significant challenges in solving the equations analytically; unlike linear systems, there exists no general closed-form solution for arbitrary matter distributions, necessitating case-by-case approaches via symmetries or approximations. In cosmology, this nonlinearity drives the evolution of structure from primordial quantum fluctuations: initial small density perturbations, linear at early times, enter a nonlinear regime where gravitational instability amplifies them, leading to the collapse and hierarchical formation of galaxies, filaments, and voids observed in the universe today. Full general relativistic simulations reveal how these nonlinear dynamics deviate from Newtonian predictions, particularly on sub-horizon scales.

Newtonian Correspondence

In the weak field, limit, the Einstein field equations recover the equations of Newtonian gravity, providing a smooth correspondence between and . This limit applies when velocities are much smaller than the (v \ll c) and gravitational potentials are weak (|\Phi| \ll c^2), conditions typical of solar system dynamics. Under these assumptions, the simplifies significantly: the time-time component approximates g_{00} \approx -(1 + 2\Phi/c^2), where \Phi is the Newtonian , while spatial components remain nearly Minkowski flat to leading order. Substituting this metric form into the 00-component of the Einstein field equations yields the Poisson equation \nabla^2 \Phi = 4\pi G \rho, where \rho is the mass-energy density, directly matching Newton's law of gravity. To probe beyond this leading-order correspondence, the systematically includes higher-order corrections parameterized by powers of v/c and \Phi/c^2. In , these expansions are governed by specific values of the parametrized post-Newtonian (PPN) parameters, such as \gamma = 1 and \beta = 1, which ensure consistency with the theory's nonlinear structure while recovering Newtonian gravity at zeroth order. Solar system tests, including perihelion precession of Mercury and Shapiro time delays, validate this correspondence by measuring deviations from Newtonian predictions that align with post-Newtonian terms at the level of $10^{-5} to $10^{-6}. This Newtonian limit demonstrates how the nonlinearity of the Einstein field equations, which complicates exact solutions, is effectively linearized under the specified conditions, bridging relativistic and classical regimes. Historically, the equivalence was underscored in the context of the Eddington expedition, where observations during a confirmed general relativity's predictions for light deflection—twice the Newtonian value—while affirming the underlying Newtonian behavior for massive bodies.

Special Cases

Vacuum Equations

In the absence of and , the stress-energy tensor vanishes, T_{\mu\nu} = 0, simplifying the Einstein field equations to their form. Without a , contracting the equations yields the R = 0, which in turn implies that the Ricci tensor must satisfy R_{\mu\nu} = 0. This condition defines Ricci-flat spacetimes, where the Ricci tensor measures the local volume distortion due to , and its vanishing indicates no net or expansion effects from gravitational sources. Ricci-flat spacetimes encompass both the flat Minkowski , representing empty, special-relativistic space, and nontrivial curved geometries sustained by the global distribution of rather than local . In these regions, the absence of energy-momentum sources ensures that the motion of test particles is purely geometric, with freely falling bodies tracing geodesics determined by the alone. This reflects the principle that gravity in arises from , allowing propagation of gravitational influences without localized energy densities. When a cosmological constant \Lambda is included to account for vacuum energy, the equations modify to R_{\mu\nu} - \frac{1}{2} R g_{\mu\nu} + \Lambda g_{\mu\nu} = 0, with the trace giving R = 4\Lambda and thus R_{\mu\nu} = \Lambda g_{\mu\nu}. This form describes spacetimes of constant , where the geometry is homogeneous and isotropic on large scales, independent of specific matter distributions. Such solutions model idealized cosmological backgrounds, highlighting the role of \Lambda in balancing expansive or contractive tendencies in .

Einstein-Maxwell Equations

The Einstein-Maxwell equations form a coupled system that incorporates the as a source for curvature in , extending the pure gravitational description to include . This system comprises the Einstein field equations, where the stress-energy tensor T_{\mu\nu} is provided by the electromagnetic field, alongside the Maxwell equations formulated covariantly in curved . In the absence of charges and currents, the inhomogeneous Maxwell equation is \nabla_\alpha F^{\mu\alpha} = 0, while the homogeneous equation follows from the Bianchi identity \nabla_{[\lambda} F_{\mu\nu]} = 0, ensuring the field strength tensor F_{\mu\nu} satisfies the appropriate differential constraints. The electromagnetic stress-energy tensor, which sources the Einstein field equations G_{\mu\nu} = 8\pi T_{\mu\nu} (in units where G = c = 1), is given by T_{\mu\nu} = -F_{\lambda\mu} F_{\nu}{}^{\lambda} + \frac{1}{4} g_{\mu\nu} F_{\alpha\beta} F^{\alpha\beta}, where F_{\mu\nu} = \partial_\mu A_\nu - \partial_\nu A_\mu is the strength derived from the four-potential A_\mu, and g_{\mu\nu} is the metric tensor. This tensor captures the energy density, momentum flux, and stresses associated with electromagnetic fields, such as those from electric and magnetic contributions. A key application of the Einstein-Maxwell equations arises in the study of charged black holes, where the Reissner-Nordström metric provides an exact, static, spherically symmetric solution describing the spacetime around a non-rotating, charged , generalizing the Schwarzschild solution to include electromagnetic effects. When electromagnetic fields vanish, the Einstein-Maxwell system reduces to the vacuum Einstein equations. The coupled equations are consistent, as both the gravitational and electromagnetic dynamics derive from a unified : the total is the sum of the Einstein-Hilbert for and the Maxwell for the , S = \int d^4x \sqrt{-g} \left( \frac{R}{16\pi} - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} \right), whose variation yields the full set of field equations.

Solutions and Approximations

Exact Solutions

The represents the simplest and most influential exact solution to the vacuum Einstein field equations, describing the geometry surrounding a spherically symmetric, non-rotating mass distribution. This metric, derived by in early 1916 just weeks after Einstein's final formulation of , yields a static, asymptotically flat that features an at the r_s = 2GM/c^2, beyond which light cannot escape, characterizing an eternal . The in standard is ds^2 = \left(1 - \frac{r_s}{r}\right) c^2 dt^2 - \left(1 - \frac{r_s}{r}\right)^{-1} dr^2 - r^2 (d\theta^2 + \sin^2\theta \, d\phi^2), where G is the gravitational constant, M the mass, and c the speed of light. Birkhoff's theorem establishes that this solution is unique among all spherically symmetric vacuum solutions, implying that any such configuration evolves without radiating gravitational waves and remains static if initially so. Physically, the Schwarzschild metric underpins predictions of phenomena like the bending of light by the Sun, the precession of Mercury's orbit, and the existence of black holes, confirmed observationally through gravitational wave detections and imaging of supermassive black holes. The Friedmann–Lemaître–Robertson–Walker (FLRW) metric furnishes an exact class of solutions to the Einstein field equations for homogeneous and isotropic spacetimes filled with , , or a \Lambda, forming the backbone of modern . First proposed by in 1922, it demonstrates that the need not be static but can expand or contract dynamically, with subsequent refinements by in 1927 interpreting expansion as a origin, and kinematic completions by Howard Robertson and Arthur Walker in the 1930s. The metric takes the form ds^2 = -c^2 dt^2 + a(t)^2 \left[ \frac{dr^2}{1 - k r^2} + r^2 (d\theta^2 + \sin^2\theta \, d\phi^2) \right], where a(t) is the time-dependent scale factor dictating expansion, and k parameterizes spatial curvature (k = 0 for flat, k > 0 for closed, k < 0 for open universes). Substituting into the Einstein equations yields the Friedmann equations, which relate \dot{a}^2/a^2 to energy density, pressure, and \Lambda, enabling models of cosmic evolution that align with observations such as the cosmic microwave background uniformity and accelerated expansion driven by dark energy. These solutions assume perfect fluid sources, briefly referencing the matter-filled cases from the stress-energy tensor formulations. The extends the Schwarzschild solution to rotating, asymptotically flat s in vacuum, capturing the effects of on geometry. Discovered by in , it describes a stationary, axisymmetric with an and an where objects must co-rotate with the due to . In Boyer-Lindquist coordinates, the metric involves mass M and angular momentum per unit mass a = J/(Mc), reducing to the Schwarzschild case when a = 0. This solution is essential for modeling astrophysical phenomena like accreting s in X-ray binaries and the supermassive s at galactic centers, with its predictions verified through signals from binary mergers exhibiting spin effects. A systematic classification of exact solutions relies on the Petrov scheme, which algebraically categorizes the —the traceless, conformally invariant part of the representing the free decoupled from local sources. Introduced by Alexei Petrov in , this classification identifies six types (I, II, III, N, O, D) based on the multiplicity of principal null directions aligned with the Weyl tensor's eigenvectors, with type D denoting algebraically special metrics having two repeated directions, as in the Schwarzschild and Kerr solutions. The Weyl tensor's role is pivotal, as its algebraic speciality simplifies the Einstein equations, enabling the generation of new exact solutions through techniques like the Newman-Penrose formalism, and it distinguishes (type N) from more general distortions (type I). This framework has guided the cataloging of thousands of known solutions, emphasizing their geometric properties over source terms.

Linearized Form

In the weak-field limit, where gravitational effects are small compared to the squared, the Einstein field equations can be linearized by the around the flat Minkowski background: g_{\mu\nu} = \eta_{\mu\nu} + h_{\mu\nu}, with |h_{\mu\nu}| \ll 1. This approximation neglects higher-order terms in h_{\mu\nu}, simplifying the nonlinear Ricci tensor and to first order in the perturbation. To further simplify the equations, a choice is imposed, known as the or de Donder gauge, where \partial^\mu \bar{h}_{\mu\nu} = 0 and \bar{h}_{\mu\nu} = h_{\mu\nu} - \frac{1}{2} \eta_{\mu\nu} h, with h = h^\lambda_\lambda. Under this condition, the linearized Einstein field equations reduce to a sourced : \square \bar{h}_{\mu\nu} = -\frac{16\pi G}{c^4} T_{\mu\nu}, where \square = \partial^\alpha \partial_\alpha is the d'Alembertian operator in . This form reveals that gravitational perturbations propagate as waves at the , analogous to electromagnetic waves, and in the quantum context, these waves are mediated by massless spin-2 gravitons. In vacuum, where T_{\mu\nu} = 0, the equation becomes the homogeneous \square \bar{h}_{\mu\nu} = 0, admitting plane-wave solutions. For propagating in the z-direction, the transverse-traceless (TT) gauge further constrains the perturbations such that h_{0\mu} = 0, \partial^i h_{ij} = 0, h^i_i = 0, and h_{ij} has only two independent polarization modes: the plus (+) and (\times) modes. These solutions describe quadrupole emitted by accelerating masses, such as merging holes or . The linearized equations also capture effects in weak fields, where rotating masses induce a gravitomagnetic field analogous to in , causing inertial frames to be dragged along with the . This Lense-Thirring arises from the off-diagonal components of h_{\mu\nu} sourced by the stress-energy tensor's density. Direct evidence for these predictions came from the and collaborations, which detected from mergers in 2015, confirming the transverse-traceless modes with strains matching the linearized theory's plus and cross polarizations. The first such observation, GW150914, exhibited a signal consistent with general relativity's weak-field wave propagation.

Numerical and Polynomial Approaches

Numerical relativity employs computational methods to solve the Einstein field equations in regimes where exact solutions are unavailable, particularly for highly dynamical involving strong . A foundational approach is the 3+1 decomposition, which foliates into a family of spatial evolving along a time direction, transforming the four-dimensional equations into a of three-dimensional spatial equations coupled with evolution equations. This decomposition underpins the Arnowitt-Deser-Misner () formalism, which recasts the Einstein equations in form, yielding explicit evolution equations for the spatial and extrinsic on each . The formalism facilitates by treating the gravitational field as a amenable to finite-difference or spectral methods, though early implementations suffered from instabilities due to the hyperbolic-parabolic nature of the equations. To address these stability issues, the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation refines the system by introducing auxiliary variables that enforce constraint propagation and improve hyperbolicity. Originally developed through modifications to the ADM variables, BSSN decomposes the spatial metric into a conformal factor and a transverse-traceless part, allowing for better control of coordinate choices and damping of constraint violations during long-term evolutions. This formulation has become the standard for stable numerical simulations of mergers, enabling accurate modeling of nonlinear gravitational dynamics. A pivotal advancement occurred in 2005, when achieved stable, long-term evolutions of spacetimes, overcoming decades of challenges with instabilities and coordinate pathologies. This , demonstrated through moving-puncture techniques combined with BSSN, allowed simulations of the inspiral, merger, and ringdown phases, providing waveforms essential for gravitational-wave detection. Applications extend to astrophysical scenarios such as inspirals, where numerical models predict orbital decay rates and energy losses matching observations from /, and core-collapse supernovae, where general relativistic simulations reveal the role of rotation and in explosion mechanisms. Beyond full numerical evolution, polynomial approaches provide approximate solutions via series expansions in powers of the strength or , valid in weakly nonlinear regimes. These expansions, such as the post-Newtonian series, iteratively solve the field equations order by order in small parameters like v/c or Gm/rc^2, capturing corrections to Newtonian gravity for systems like binary pulsars. By truncating at desired orders, these methods yield analytic expressions for metrics and trajectories, bridging exact solutions and full numerics while informing modeling for inspiraling compact objects.

References

  1. [1]
    [PDF] Die Feldgleichungen der Gravitation.
    EINSTEIN: Die Feldgleichungen der Gravitation. 845. Die allgemein kovarianten zehn Gleichungen des Gravitationsfeldes in Räumen, in denen » Materie « fehlt ...
  2. [2]
    [PDF] The Meaning of Einstein's Equation - UCR Math Department
    Jan 4, 2006 · We have already seen that it explains how gravity attracts. We sketch a few other consequences below.
  3. [3]
    None
    ### Summary of Einstein Field Equations
  4. [4]
    [PDF] JOHN NORTON - How Einstein found his field equations: 1912-1915
    In the. Entwurf paper, Einstein and Grossmann had come within a hair's breadth of the generally covariant field equations of the final theory. *University of ...Missing: primary | Show results with:primary
  5. [5]
    One Hundred Years of Einstein's Field Equations, and of the ...
    100 years ago Albert Einstein submitted a paper to the Royal Prussian Academy of the Natural Sciences which contained what we today call the Einstein field ...
  6. [6]
    Einstein's Pathway to General Relativity - University of Pittsburgh
    The immediate outcome of this reflection was Einstein's "principle of equivalence." It formed the basis of the concluding Part V of his 1907 Jahrbuch article.
  7. [7]
    [PDF] Einstein's Pathway to the Equivalence Principle 1905-1907 - arXiv
    Between 1905 and 1907, Einstein first tried to extend the special theory of relativity in such a way so as to explain gravitational phenomena.
  8. [8]
    [PDF] einstein's quest for general relativity, 1907–1920 - Physics - UMD
    He thereupon modified his field equations (without compromising their general covariance) by adding a term with the so-called cosmological constant (Einstein ...<|control11|><|separator|>
  9. [9]
    Einstein's Zurich Notebook - University of Pittsburgh
    Its spacetime trajectory is a geodesic of the spacetime. That is, it traces a curve of extremal spacetime interval in spacetime.
  10. [10]
    Topics
    1913: The Entwurf theory. In collaboration with his mathematician friend, Marcel Grossmann, the two wrote a sketch of a generalized theory of relativity and ...
  11. [11]
    [PDF] Correspondence between Einstein and Levi- Civita Galina Weinstein
    Jan 31, 2012 · Before developing his 1915 General Theory of Relativity, Einstein held the "Entwurf" theory. Tullio Levi-Civita from Padua, one of the ...Missing: connection | Show results with:connection
  12. [12]
  13. [13]
    [physics/9811050] The Relativity of Discovery: Hilbert's First Note on ...
    Nov 26, 1998 · An intense correspondence between Hilbert and Einstein in the crucial month of November 1915, furthermore, confronts the historian with a case ...
  14. [14]
    The Field Equations of Gravitation - Wikisource, the free online library
    Aug 9, 2025 · In German: Die Feldgleichungen der Gravitation, Preussische Akademie der Wissenschaften, Sitzungsberichte, 1915 (part 2), 844–847, Internet ...
  15. [15]
  16. [16]
    Belated Decision in the Hilbert-Einstein Priority Dispute - jstor
    Nov 11, 2024 · November 1915 but published only on 31 March 1916, presents a generally covariant theory of gravitation, including field equations essentially ...
  17. [17]
    General relativity - MacTutor History of Mathematics
    In December 1915 Ehrenfest wrote to Lorentz referring to the theory of November 25, 1915. ... A J Kox, Einstein, Lorentz, Leiden and general relativity, Les ...
  18. [18]
    Albert Einstein's 1916 Review Article on General Relativity - arXiv
    May 13, 2004 · The first comprehensive overview of the final version of the general theory of relativity was published by Einstein in 1916.
  19. [19]
    [PDF] EINSTEIN'S EQUATIONS, COSMOLOGY AND ASTROPHYSICS
    This section is devoted to the genesis and properties of the field equations. The notation is standard, so x0,123 are the coordinates of time and ordinary space ...
  20. [20]
    [PDF] General Relativity - DAMTP - University of Cambridge
    A really excellent book that will satisfy your geometrical and topological needs for this course and much beyond. It is particularly useful for Sections 2 ...
  21. [21]
    [PDF] Lecture Notes on General Relativity - arXiv
    in a Euclidean signature metric they are orthogonal transformations, while in a. Lorentzian signature metric they are Lorentz transformations. We ...
  22. [22]
    [PDF] Sign Conventions in General Relativity
    Dec 11, 2023 · This article highlights some sign conventions that appear to be standard and describes some of the effects of changing the signs of the metric ...
  23. [23]
    Particle physics signature versus relativity signature of metric tensor
    Mar 26, 2017 · The physics doesn't change because both conventions give the same definition for all physical quantities.Advantages of using different metric signatures in relativity and QFTWhy are there two metric signature conventions? [duplicate]More results from physics.stackexchange.com
  24. [24]
    6. Weak Fields and Gravitational Radiation
    (The Einstein tensor is simply the trace-reversed Ricci tensor.) In terms of $ \bar{h}_{\mu\nu}^{} the harmonic gauge condition becomes. Equation 6.21, (6.21) ...
  25. [25]
    [PDF] arXiv:2405.06035v1 [gr-qc] 9 May 2024
    May 9, 2024 · We begin with the trace-reversed Einstein equations and, as above, define auxiliary constraints, this time called Za. We then replace the ...
  26. [26]
    [PDF] arXiv:1703.06331v1 [gr-qc] 18 Mar 2017
    Mar 18, 2017 · ... trace-reversed form, we have after taking the trace. Rg = 2 (3Λ′ − κtr(Tµν). (7). Thus the trace-reversed Einstein's equations are. ¯Rµν = 2Λ ...
  27. [27]
    [1112.1572] Palatini Approach Beyond Einstein's Gravity - arXiv
    Dec 7, 2011 · I review recent results obtained for extensions of general relativity formulated within the Palatini formalism, an approach in which metric and ...
  28. [28]
    [PDF] Palatini's formulation of General Relativity - arXiv
    Nov 28, 2023 · This paper aims to construct a pre-symplectic structure within General Relativity, specifically in Palatini's approach, which can be formulated ...
  29. [29]
    [PDF] Geodesic Deviation and Riemann Curvature Tensor
    Clearly, the vaccum Einstein equations can also be written as,. Since. Gur. == Gμr =o <=> Rar =0 -. Alternative form of vacuum. Einstein field equations.
  30. [30]
    [PDF] Perfect Fluids in General Relativity - MPG.PuRe
    Dec 15, 1970 · Stress-Energy Tensor and Equations of Motion. The relativistic one-component perfect fluid is defined by its equation of state, Eq. (2.6) ...
  31. [31]
    Stress-energy tensor - Cosmological Dynamics - E. Bertschinger
    For a perfect fluid the stress-energy tensor takes the well-known form ... We may already anticipate that these sources are responsible in the Einstein equations ...
  32. [32]
    GP-B — Einstein's Spacetime - Gravity Probe B
    In 1918, Einstein described Mach's principle as a philosophical pillar of general relativity, along with the physical principle of equivalence and the ...
  33. [33]
    Einstein's 1917 Static Model of the Universe: A Centennial Review
    Jan 25, 2017 · We present a historical review of Einstein's 1917 paper 'Cosmological Considerations in the General Theory of Relativity' to mark the centenary of a key work.
  34. [34]
    Einstein's conversion from his static to an expanding universe - arXiv
    Nov 12, 2013 · In 1922 Alexander Friedman showed that Einstein's fundamental equation also allowed dynamical worlds, and in 1927 Geroges Lemaitre, backed by ...Missing: Friedmann | Show results with:Friedmann
  35. [35]
    [2403.01555] Prediction of the Expansion of the Universe made by ...
    Mar 3, 2024 · This review devoted to the centenary of Alexander Friedmann's prediction of the Universe expansion presents the results obtained by him in 1922 and 1924.
  36. [36]
    [astro-ph/9805201] Observational Evidence from Supernovae for an ...
    May 15, 1998 · Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant. Authors:Adam G. Riess, Alexei V. Filippenko, ...
  37. [37]
    Cosmological constant and vacuum energy: old and new ideas - arXiv
    Jun 6, 2013 · The cosmological constant (CC) term in Einstein's equations, Lambda, was first associated to the idea of vacuum energy density.
  38. [38]
    [astro-ph/0203330] On the Cosmological Constant Problems ... - arXiv
    Mar 20, 2002 · After a short history of the cosmological term it is explained why the (effective) cosmological constant is expected to obtain contributions ...
  39. [39]
    [PDF] Noether's Theorems and Energy in General Relativity - arXiv
    Mar 31, 2021 · This paper has three main aims: first, to give a pedagogical introduction to Noether's two theorems and their implications for energy ...
  40. [40]
    [PDF] Noether's first theorem and the energy-momentum tensor ... - arXiv
    Jul 21, 2021 · The derivation of the canonical stress-energy tensor fails to use the full power of the mathematical machinery that Noether has given us by ...
  41. [41]
    [PDF] On the conservation laws of electrodynamics - Neo-classical physics
    By. Erich Bessel-Hagen in Göttingen. Translated by D. H. ... and on the basis of the definition of the electromagnetic energy-impulse tensor Sik that.
  42. [42]
    [PDF] From Noether's Theorem to Bremsstrahlung - arXiv
    Dec 10, 2021 · To do this, we will use the stress energy tensor of the EM field. Tµν = −F α. µ Fνα + 1. 4. gµνFαβFαβ. (3.33). We will use that the energy ...
  43. [43]
    [PDF] Energy conditions in general relativity and quantum field theory - arXiv
    Jun 5, 2020 · Abstract. This review summarizes the current status of the energy conditions in general relativity and quantum field theory.
  44. [44]
    [PDF] Revisiting the Bianchi identity 12.2 Making a theory of gravity - MIT
    The way we have derived the Einstein field equation here parallels how Einstein originally derived it: deduce that the source should be Tµν, infer that the ...
  45. [45]
    [PDF] Physics 225A: General Relativity
    Dec 18, 2024 · Stress-energy tensor in curved space. Earlier, we applied the Noether method to translation symmetry to find. T. µ flatν = −. ∂L. ∂(∂µφ).
  46. [46]
    [PDF] 4. The Einstein Equations - DAMTP
    The force of gravity is mediated by a gravitational field. The glory of general relativity is that this field is identified with a metric gµν(x) on a 4d ...
  47. [47]
    No-Go Theorem for Energy-Momentum Conservation in Curved ...
    Feb 10, 2007 · The law of conservation of energy-momentum holds valid only approximately when spacetime is not seriously curved.
  48. [48]
    [PDF] Newtonian approximation in gravitational physics
    In general relativity, the gravitational field itself is a source of gravity, a reflection of the nonlinearity of Einstein's equations, and in contrast to the.
  49. [49]
    [PDF] Exact Solutions of Einstein's Field Equations
    A revised edition of the now classic text, Exact Solutions of Einstein's Field Equations gives a unique survey of the known solutions of Einstein's field ...
  50. [50]
    [PDF] THE FOUNDATION OF THE GENERAL THEORY OF RELATIVITY
    he “HE special theory of relativity is based on the following postulate, which is also satisfied by the mechanics of Galileo and Newton.
  51. [51]
    On the unreasonable effectiveness of the post-Newtonian ... - PNAS
    The post-Newtonian approximation is a method for solving Einstein's field equations for physical systems in which motions are slow compared to the speed of ...
  52. [52]
    [PDF] GRAVITATION F10 Lecture 13 1. The Vacuum Einstein's Equation ...
    1.1. The metric tensor of space-time satisfies a nonlinear Partial Differ- ential Equation which determines it given initial conditions. The metric tensor ...
  53. [53]
    Introduction to Einstein-Maxwell equations and the Rainich conditions
    Jun 28, 2016 · Abstract:The first results of Einstein-Maxwell equations established by Raincih in 1925 are therefore called the Raincih conditions.
  54. [54]
    [PDF] Stress-Energy tensor for Maxwell Theory - Rutgers Physics
    Here the field strength tensor is defined by. Fµν = ∂µAν − ∂νAµ, where, as we have seen in Homework 1, Problem 3, consistency requires. ∂νjν = 0. Consider first ...
  55. [55]
    Über die Eigengravitation des elektrischen Feldes nach der ...
    Über die Eigengravitation des elektrischen Feldes nach der Einsteinschen Theorie. H. Reissner, ... First published: 1916. https://doi.org/10.1002/andp.19163550905.Missing: original | Show results with:original
  56. [56]
  57. [57]
    [physics/9905030] On the gravitational field of a mass point ... - arXiv
    May 12, 1999 · Plain TeX, 7 pages, English translation of the original paper by K. Schwarzschild. Subjects: History and Philosophy of Physics (physics.hist ...
  58. [58]
    Exact Solutions of Einstein's Field Equations
    Cambridge Core - Theoretical Physics and Mathematical Physics - Exact Solutions of Einstein's Field Equations.
  59. [59]
    [PDF] 15.1 Solving the linearized Einstein equation for dynamical sources
    15.1 Solving the linearized Einstein equation for dynamical sources. When sources are not dynamical, then the linearized Einstein field equation. D¯hαβ ...
  60. [60]
    [PDF] Introduction. Weak Field General Relativity.
    Far from the source the Einstein field equations reduce to linearized gravity. These limitations of the linearized gravity approach show up in Eq. 6 [1,2,3].
  61. [61]
    Observation of Gravitational Waves from a Binary Black Hole Merger
    Feb 11, 2016 · This is the first direct detection of gravitational waves and the first observation of a binary black hole merger. Comments: 16 pages, 4 figures.
  62. [62]
    Republication of: The dynamics of general relativity
    Aug 8, 2008 · This article—summarizing the authors' then novel formulation of General Relativity—appeared as Chap. 7, pp. 227–264, in Gravitation.
  63. [63]
    Evolution of Binary Black-Hole Spacetimes | Phys. Rev. Lett.
    The study uses a numerical code to simulate binary black hole evolution, including orbit, merger, and ringdown, and the resulting black hole's properties.
  64. [64]
    Core-collapse supernovae as supercomputing science: A status ...
    This is a status report on our endeavor to reveal the mechanism of core-collapse supernovae (CCSNe) by large-scale numerical simulations.