Stress-energy tensor as functional derivative of the matter action

Intuition Beginner

Einstein's equation has two sides. The left side measures how spacetime is curved. The right side measures what is doing the curving: the energy and momentum carried by matter and fields. For the equation to make sense, you need a single object that packages all of that source information -- energy density, momentum, pressure, internal stresses -- into one tidy bundle. That object is the stress-energy tensor.

Where does it come from? Every kind of matter is described by its own action, a single number assigned to each possible history of the matter. A remarkable fact is that you do not have to invent the stress-energy tensor separately. It is already hidden inside the matter action. You read it off by asking a precise question: if I gently change the shape of spacetime, by how much does the matter action change? The answer, suitably arranged, is exactly the stress-energy tensor.

This is a clean and economical idea. The same metric that appears in the curvature side also appears inside the matter action, because matter feels the geometry it lives in. Probing the matter action with a small change of metric extracts the source term, and it does so automatically with the right symmetries built in. David Hilbert used this route in 1915 to write gravity and matter as one variational system.

Visual Beginner

The picture shows the extraction. The matter action depends on the spacetime shape. Nudging the shape and watching the response returns a four-by-four table whose entries are energy density, momentum, pressure, and shear. That table is the stress-energy tensor.

Worked example Beginner

Take the simplest field: a single number $\varphi$ defined at every point of spacetime, a scalar field. Think of it as a temperature-like quantity stored everywhere. Its energy comes from two sources: how fast it changes from place to place, and how large it is (a mass-like cost for the field to be nonzero at all).

When you carry out the extraction -- nudge the spacetime shape and read off the response -- the energy density that comes out is a sum of these contributions. In a region where the field is changing rapidly in time, the energy density is large. In a region where the field is uniform and small, the energy density is small. The momentum entries record how energy flows from one place to another, which happens whenever the field changes in both time and space at once.

Here is the takeaway worth holding onto. You did not have to guess the energy density of the scalar field or postulate a formula for its pressure. Both fell out of one computation applied to the field's action. The stress-energy tensor is not an extra ingredient bolted onto the theory. It is the matter action, viewed through the lens of "how does this respond to a change of geometry." The same procedure works for the electromagnetic field, for a fluid, for any matter whose action you can write down.

Check your understanding Beginner

Formal definition Intermediate+

Fix an oriented four-dimensional manifold $M$ carrying a Lorentzian metric $g_{\mu \nu }$ of signature $(-,+,+,+)$, with Levi-Civita connection $\nabla$, and write $g=\det (g_{\mu \nu })$, so that $\sqrt{-g}{d}^{4}x$ is the invariant volume element. Matter is described by fields $\psi$ (collectively denoting scalars, vectors, tensors, or spinors) and a matter action

$$S_m[g,\psi ]={\int }_M{\mathcal{L}}_m(g,\psi ,\nabla \psi )\sqrt{-g}{d}^{4}x,$$

a diffeomorphism-invariant scalar functional of the metric and the matter fields. The symbol $\nabla$ denotes covariant differentiation; $\delta$ denotes a first-order field variation; $\int$ is integration against $d^{4}x$; $\partial$ denotes partial differentiation.

The Hilbert (metric) stress-energy tensor is defined by the metric variation of the matter action:

$$T_{\mu \nu }:=-\frac{2}{\sqrt{-g}}\frac{\delta S_m}{\delta g^{\mu \nu }}=-\frac{2}{\sqrt{-g}}\frac{\delta (\sqrt{-g}{\mathcal{L}}_m)}{\delta g^{\mu \nu }}.$$

Carrying the variation of the volume density through with $\delta \sqrt{-g}=-\frac{1}{2}\sqrt{-g}{g}_{\mu \nu }\delta g^{\mu \nu }$ from 13.04.02 gives the equivalent working form

$$T_{\mu \nu }=-2\frac{\partial {\mathcal{L}}_m}{\partial g^{\mu \nu }}+g_{\mu \nu }{\mathcal{L}}_m,$$

valid when ${\mathcal{L}}_m$ depends on the metric only algebraically (no derivatives of $g$); the functional-derivative form is the general definition. The factor $-2$ and the sign are fixed so that $T_{\mu \nu }$ appears as the source on the right of the Einstein equation $G_{\mu \nu }=\kappa T_{\mu \nu }$ with $\kappa =8\pi G/c^4$, and so that the energy density $T_{00}$ is positive for ordinary matter.

Sign and convention notes

The signature is $(-,+,+,+)$. With this signature and the definition above, a scalar field with positive kinetic and potential energy has $T_{00}>0$. Texts using the opposite signature $(+,-,-,-)$ absorb the change into the overall sign of ${\mathcal{L}}_m$ and of the definition, so the physical tensor is unchanged. The functional derivative $\delta S_m/\delta g^{\mu \nu }$ is the coefficient of $\delta g^{\mu \nu }$ in $\delta S_m$ after every total-divergence term has been integrated away, on the matter fields' equations of motion.

Counterexamples to common slips

• The Hilbert tensor is defined by varying with respect to the inverse metric $g^{\mu \nu }$, not $g_{\mu \nu }$; the two differ by a sign through $\delta g_{\mu \nu }=-g_{\mu \alpha }{g}_{\nu \beta }\delta g^{\mu \nu }$, and confusing them flips the sign of $T_{\mu \nu }$.
• The Hilbert tensor is not in general equal to the canonical (Noether) stress tensor $T_{\mu \nu }^{\mathrm{can}}=\frac{\partial {\mathcal{L}}_m}{\partial ({\partial }^{\mu }\psi )}{\partial }_{\nu }\psi -g_{\mu \nu }{\mathcal{L}}_m$. They agree only up to a Belinfante-Rosenfeld improvement term that is itself identically conserved. The Hilbert tensor is symmetric by construction; the canonical tensor need not be.
• Tracelessness $T=g^{\mu \nu }{T}_{\mu \nu }=0$ is a property of conformally invariant matter (such as the electromagnetic field in four dimensions), not a general feature. A massive scalar field has $T\ne 0$.

Key theorem with proof Intermediate+

Theorem (symmetry and on-shell conservation of the Hilbert tensor). Let $S_m[g,\psi ]$ be a diffeomorphism-invariant matter action on a four-dimensional Lorentzian manifold, and let $T_{\mu \nu }=-\frac{2}{\sqrt{-g}}\frac{\delta S_m}{\delta g^{\mu \nu }}$. Then (i) $T_{\mu \nu }$ is symmetric, $T_{\mu \nu }=T_{\nu \mu }$, and (ii) on solutions of the matter field equations $\delta S_m/\delta \psi =0$, it is covariantly conserved, ${\nabla }^{\mu }{T}_{\mu \nu }=0$. ^{[Wald Appendix E.1]}

Proof. Symmetry. The inverse metric $g^{\mu \nu }$ is a symmetric tensor: $g^{\mu \nu }=g^{\nu \mu }$. The functional derivative of a scalar with respect to a symmetric tensor is itself symmetric, because only the symmetric part of $\delta g^{\mu \nu }$ contributes to $\delta S_m$. Concretely, write $\delta S_m=\int A_{\mu \nu }\delta g^{\mu \nu }{d}^{4}x$; since $\delta g^{\mu \nu }$ is symmetric, replacing $A_{\mu \nu }$ by its symmetric part $A_{(\mu \nu )}=\frac{1}{2}(A_{\mu \nu }+A_{\nu \mu })$ leaves $\delta S_m$ unchanged, so the functional derivative is well-defined only as a symmetric tensor. Hence $T_{\mu \nu }$, being $-\frac{2}{\sqrt{-g}}{A}_{\mu \nu }$ read as a symmetric object, satisfies $T_{\mu \nu }=T_{\nu \mu }$. The symmetry is built into the definition; it requires no improvement term.

Conservation. Consider an infinitesimal diffeomorphism generated by a compactly supported vector field ${\xi }^{\mu }$. Under its flow the dynamical fields change by their Lie derivatives: the metric by $\delta g_{\mu \nu }=({\mathcal{L}}_{\xi }g{)}_{\mu \nu }={\nabla }_{\mu }{\xi }_{\nu }+{\nabla }_{\nu }{\xi }_{\mu }$, and the matter fields by $\delta \psi ={\mathcal{L}}_{\xi }\psi$. Because $S_m$ is a diffeomorphism-invariant scalar, its value is unchanged by this active transformation of all fields: ${\delta }_{\xi }{S}_m=0$ identically. Splitting the variation into its metric and matter parts,

$$0={\delta }_{\xi }{S}_m={\int }_M\frac{\delta S_m}{\delta g^{\mu \nu }}\delta g^{\mu \nu }{d}^{4}x+{\int }_M\frac{\delta S_m}{\delta \psi }\delta \psi d^{4}x.$$

On solutions of the matter field equations, $\delta S_m/\delta \psi =0$, so the second integral vanishes and the first must vanish on its own. Insert the definition $\frac{\delta S_m}{\delta g^{\mu \nu }}=-\frac{1}{2}\sqrt{-g}{T}_{\mu \nu }$ and the variation $\delta g^{\mu \nu }=-({\nabla }^{\mu }{\xi }^{\nu }+{\nabla }^{\nu }{\xi }^{\mu })$ (obtained by raising both indices on $\delta g_{\mu \nu }$ with a sign, since $\delta g^{\mu \nu }=-g^{\mu \alpha }{g}^{\nu \beta }\delta g_{\alpha \beta }$):

$$0=-\frac{1}{2}{\int }_M{T}_{\mu \nu }\delta g^{\mu \nu }\sqrt{-g}{d}^{4}x={\int }_M{T}_{\mu \nu }{\nabla }^{\mu }{\xi }^{\nu }\sqrt{-g}{d}^{4}x,$$

where the factor of two from the symmetric combination ${\nabla }^{\mu }{\xi }^{\nu }+{\nabla }^{\nu }{\xi }^{\mu }$ is absorbed using the symmetry $T_{\mu \nu }=T_{\nu \mu }$ established in part (i). Integrate by parts with the covariant divergence theorem; the boundary term vanishes because ${\xi }^{\nu }$ has compact support:

$$0=-{\int }_M({\nabla }^{\mu }{T}_{\mu \nu }){\xi }^{\nu }\sqrt{-g}{d}^{4}x.$$

Since ${\xi }^{\nu }$ is an arbitrary compactly supported vector field, the fundamental lemma of the calculus of variations forces ${\nabla }^{\mu }{T}_{\mu \nu }=0$ pointwise. $■$

Bridge. This result builds toward 13.04.01, where covariant conservation ${\nabla }^{\mu }{T}_{\mu \nu }=0$ is the consistency condition that lets the Einstein equation $G_{\mu \nu }=\kappa T_{\mu \nu }$ be solved, since the left side is divergence-free by the contracted Bianchi identity. It appears again in 13.08.01, where the conservation law applied to the perfect-fluid tensor yields the continuity equation governing the energy density of the expanding universe, and in 13.05.01, where the vanishing of $T_{\mu \nu }$ outside a star defines the vacuum region whose Schwarzschild geometry is then determined. The diffeomorphism-invariance argument used here is the matter-side counterpart of the contracted Bianchi identity of 05.05.08: the same Noether-identity mechanism that makes the geometric side automatically conserved makes the matter source conserved on-shell, and the two halves must match for the field equation to be consistent.

Exercises Intermediate+

Lean formalization Intermediate+

The definition and properties of the Hilbert stress-energy tensor rest on a tower of structures Mathlib does not provide, so this unit ships with lean_status: none. From the bottom:

1. Lorentzian metrics and the volume density. Mathlib has positive-definite Riemannian metrics but no indefinite Lorentzian metric as a first-class object. The density $\sqrt{-g}$ and its variation $\delta \sqrt{-g}=-\frac{1}{2}\sqrt{-g}{g}_{\mu \nu }\delta g^{\mu \nu }$, which carry the $g_{\mu \nu }{\mathcal{L}}_m$ term of the tensor, are unavailable.

2. The functional derivative on the space of metrics. The definition $T_{\mu \nu }=-\frac{2}{\sqrt{-g}}\delta S_m/\delta g^{\mu \nu }$ is a Gateaux derivative of a functional on the infinite-dimensional space of metrics. Mathlib's variational-calculus and variational-bicomplex coverage for fields on manifolds is essentially absent.

3. Covariant divergence and the Lie derivative of the metric. The conservation theorem uses $\delta g_{\mu \nu }={\nabla }_{\mu }{\xi }_{\nu }+{\nabla }_{\nu }{\xi }_{\mu }$ and an integration by parts with the covariant divergence. Neither the Levi-Civita connection of a Lorentzian metric nor the Lie derivative of the metric tensor is present as a named object.

A realistic first target is purely tensor-algebraic and frame-independent: given ${\mathcal{L}}_m$, $g^{\mu \nu }$, and the matter fields as data, define $T_{\mu \nu }$ by the working form $-2\partial {\mathcal{L}}_m/\partial g^{\mu \nu }+g_{\mu \nu }{\mathcal{L}}_m$ and prove its symmetry from the symmetry of $g^{\mu \nu }$, then verify tracelessness for the Maxwell Lagrangian as a finite linear-algebra identity in four dimensions. The conservation statement waits on items 1-3.

Advanced results Master

The canonical tensor and the Belinfante-Rosenfeld improvement. Noether's first theorem applied to spacetime translation symmetry in flat space produces the canonical stress-energy tensor

$$T^{\mathrm{can}}{}^{\mu }{}_{\nu }=\frac{\partial {\mathcal{L}}_m}{\partial ({\partial }_{\mu }\psi )}{\partial }_{\nu }\psi -{\delta }_{\nu }^{\mu }{\mathcal{L}}_m,$$

derived in the Noether framework of 05.05.08. This tensor is conserved, ${\partial }_{\mu }{T}^{\mathrm{can}}{}^{\mu }{}_{\nu }=0$, but for fields carrying spin it is generally not symmetric: $T_{\mu \nu }^{\mathrm{can}}\ne T_{\nu \mu }^{\mathrm{can}}$. The electromagnetic canonical tensor $\frac{\partial {\mathcal{L}}_m}{\partial ({\partial }_{\mu }{A}_{\lambda })}{\partial }_{\nu }{A}_{\lambda }-{\delta }_{\nu }^{\mu }{\mathcal{L}}_m=-F^{\mu \lambda }{\partial }_{\nu }{A}_{\lambda }-{\delta }_{\nu }^{\mu }{\mathcal{L}}_m$ is neither symmetric nor gauge-invariant. The Belinfante-Rosenfeld construction repairs both defects by adding an improvement term built from the spin current $S^{\lambda \mu \nu }$:

$$T_{\mu \nu }^{\mathrm{BR}}=T_{\mu \nu }^{\mathrm{can}}+{\partial }_{\lambda }{B}^{\lambda }{}_{\mu \nu },B^{\lambda }{}_{\mu \nu }=\frac{1}{2}(S^{\lambda }{}_{\mu \nu }+S_{\mu \nu }{}^{\lambda }+S_{\nu \mu }{}^{\lambda }),$$

with $B^{\lambda }{}_{\mu \nu }$ antisymmetric in $(\lambda ,\mu )$, so that ${\partial }_{\mu }{\partial }_{\lambda }{B}^{\lambda \mu }{}_{\nu }=0$ and conservation is preserved by the improvement. The resulting $T_{\mu \nu }^{\mathrm{BR}}$ is symmetric and gauge-invariant, and for the electromagnetic field it reproduces $F_{\mu \beta }{F}_{\nu }{}^{\beta }-\frac{1}{4}{g}_{\mu \nu }{F}^2$ exactly ^{[Belinfante 1940] [Rosenfeld 1940]}.

The Hilbert tensor equals the improved tensor. The deep fact is that the metric variation performs the Belinfante-Rosenfeld improvement automatically. Coupling the matter to a general metric and varying it tests the response of the action to the symmetric deformation $\delta g_{\mu \nu }$, which is exactly the combination that the symmetrized improvement extracts. One proves that $T_{\mu \nu }^{\mathrm{Hilb}}=T_{\mu \nu }^{\mathrm{BR}}$ for any matter whose action is diffeomorphism-invariant and minimally coupled: the curved-space covariantization supplies the spin-connection terms whose variation generates precisely the improvement ${\partial }_{\lambda }{B}^{\lambda }{}_{\mu \nu }$. The Hilbert prescription thus delivers the symmetric, gauge-invariant, conserved tensor in one step, with no need to identify and add the improvement term by hand. This is the technical content of the statement that the metric stress-energy tensor is "the right one."

The trace and conformal invariance. Under a Weyl rescaling $g_{\mu \nu }\mapsto e^{2\omega }{g}_{\mu \nu }$, the matter action changes by $\delta S_m=\int \frac{\delta S_m}{\delta g_{\mu \nu }}\delta g_{\mu \nu }{d}^{4}x$, and with $\delta g_{\mu \nu }=2\omega g_{\mu \nu }$ at first order this is $\delta S_m=\int (-\frac{1}{2}\sqrt{-g}{T}^{\mu \nu })(2\omega g_{\mu \nu })d^{4}x=-\int \omega T^{\mu }{}_{\mu }\sqrt{-g}{d}^{4}x$. Hence the action is invariant under all Weyl rescalings if and only if the stress-energy tensor is traceless, $T=T^{\mu }{}_{\mu }=0$. The massless, conformally coupled scalar and the four-dimensional electromagnetic field are conformally invariant and have traceless $T_{\mu \nu }$; a mass scale or a dimensionful coupling breaks the symmetry and gives $T\ne 0$. In the quantum theory this classical statement acquires the conformal (trace) anomaly: even a classically traceless tensor develops $⟨T^{\mu }{}_{\mu }⟩\propto (c{C}^2-aE)$ built from the Weyl tensor squared and the Euler density, a one-loop effect with no classical analogue.

Synthesis. The Hilbert prescription resolves several classical difficulties at once and reorganizes the source side of gravity around a single definition. First, the stress-energy tensor is no longer postulated case by case; it is the metric functional derivative of whatever matter action one writes, and the same formula serves scalars, gauge fields, fluids, and spinors. Second, symmetry is automatic, not engineered: the tensor is built from the symmetric pair $(g,g^{\mu \nu })$, so it inherits symmetry without a Belinfante-Rosenfeld step, and in fact equals the improved canonical tensor for minimally coupled diffeomorphism-invariant matter. Third, conservation is automatic on-shell, following from diffeomorphism invariance of the matter action by the same Noether-identity mechanism of 05.05.08 that makes the geometric side divergence-free, so the two sides of the Einstein equation are conserved in lockstep. Fourth, the trace encodes conformal structure: classical Weyl invariance is exactly tracelessness, and its quantum breaking is the conformal anomaly. These four facts -- universality of the definition, automatic symmetry, automatic conservation, and the conformal trace -- are the structural content of treating the matter source as a functional derivative rather than a separate input.

Full proof set Master

Proposition (symmetry of the Hilbert tensor). For any scalar matter action $S_m[g,\psi ]$, the tensor $T_{\mu \nu }=-\frac{2}{\sqrt{-g}}\frac{\delta S_m}{\delta g^{\mu \nu }}$ is symmetric in $\mu ,\nu$.

Proof. The metric variation $\delta g^{\mu \nu }$ is a symmetric $(2,0)$-tensor, since $g^{\mu \nu }=g^{\nu \mu }$ and the space of metrics is the space of symmetric nondegenerate tensors. Write $\delta S_m={\int }_M{A}_{\mu \nu }\delta g^{\mu \nu }\sqrt{-g}{d}^{4}x$ for some coefficient $A_{\mu \nu }$. Decompose $A_{\mu \nu }=A_{(\mu \nu )}+A_{[\mu \nu ]}$ into symmetric and antisymmetric parts. The contraction of an antisymmetric tensor with the symmetric $\delta g^{\mu \nu }$ vanishes identically: $A_{[\mu \nu ]}\delta g^{\mu \nu }=0$. Hence $\delta S_m={\int }_M{A}_{(\mu \nu )}\delta g^{\mu \nu }\sqrt{-g}{d}^{4}x$, and only the symmetric part $A_{(\mu \nu )}$ is determined by the functional derivative. The functional derivative is therefore well-defined only as a symmetric tensor, and $T_{\mu \nu }=-2A_{(\mu \nu )}$ satisfies $T_{\mu \nu }=T_{\nu \mu }$. No improvement term, gauge choice, or field-dependent correction is needed; symmetry is a property of varying with respect to a symmetric object. $■$

Proposition (tracelessness of Maxwell theory in four dimensions). The Maxwell stress-energy tensor $T_{\mu \nu }=F_{\mu \beta }{F}_{\nu }{}^{\beta }-\frac{1}{4}{g}_{\mu \nu }{F}_{\alpha \beta }{F}^{\alpha \beta }$ is traceless in spacetime dimension $n=4$, and only then.

Proof. Contract with $g^{\mu \nu }$. The first term gives $g^{\mu \nu }{F}_{\mu \beta }{F}_{\nu }{}^{\beta }=F^{\nu }{}_{\beta }{F}_{\nu }{}^{\beta }=F_{\alpha \beta }{F}^{\alpha \beta }$, writing $F^2:=F_{\alpha \beta }{F}^{\alpha \beta }$. The second term gives $-\frac{1}{4}(g^{\mu \nu }{g}_{\mu \nu })F^2=-\frac{n}{4}{F}^2$, using $g^{\mu \nu }{g}_{\mu \nu }=n$ in dimension $n$. Therefore

$$T=g^{\mu \nu }{T}_{\mu \nu }=F^2-\frac{n}{4}{F}^2=(1-\frac{n}{4})F^2.$$

For generic field configurations $F^2\ne 0$, so $T=0$ requires $1-n/4=0$, i.e. $n=4$. In any other dimension the trace is nonzero. The coefficient $1-n/4$ is the obstruction to Weyl invariance of Maxwell theory away from four dimensions, matching the fact that the gauge field $A_{\mu }$ has conformal weight zero precisely in $n=4$. $■$

Proposition (conservation implies the geometric consistency of the Einstein equation). If ${\nabla }^{\mu }{T}_{\mu \nu }=0$ holds, then it is consistent to set $G_{\mu \nu }=\kappa T_{\mu \nu }$, since the left side is identically divergence-free.

Proof. The contracted second Bianchi identity gives ${\nabla }^{\mu }{G}_{\mu \nu }=0$ for every metric, as an identity of the Levi-Civita connection (the geometric Noether identity of 05.05.08). Taking the covariant divergence of the proposed field equation $G_{\mu \nu }=\kappa T_{\mu \nu }$ yields $0={\nabla }^{\mu }{G}_{\mu \nu }=\kappa {\nabla }^{\mu }{T}_{\mu \nu }$. This is satisfied automatically because ${\nabla }^{\mu }{T}_{\mu \nu }=0$ was established from diffeomorphism invariance of the matter action. The two divergence laws are independent in origin -- one geometric, one from matter symmetry -- yet they must agree for the field equation to admit solutions, and the variational framework guarantees that they do: both are instances of the same diffeomorphism-Noether identity, one for the gravitational action and one for the matter action. $■$

Connections Master

The stress-energy tensor developed here is the source on the right of the Einstein field equations 13.04.01; covariant conservation ${\nabla }^{\mu }{T}_{\mu \nu }=0$ is the consistency condition matching the contracted Bianchi identity of the geometric side, without which the equation $G_{\mu \nu }=\kappa T_{\mu \nu }$ could not be solved.

The definition $T_{\mu \nu }=-\frac{2}{\sqrt{-g}}\delta S_m/\delta g^{\mu \nu }$ is exactly the source term that emerged from the matter variation in the variational derivation of 13.04.02; this unit develops the properties of that term, where the prior unit established only its definition in passing.

Covariant conservation of the Hilbert tensor is the matter-side instance of Noether's second theorem for the diffeomorphism gauge group 05.05.08; the same local-symmetry mechanism that forces the geometric Euler-Lagrange tensor to be divergence-free forces the matter source to be conserved on-shell.

The perfect-fluid form $T_{\mu \nu }=(\rho +p)u_{\mu }{u}_{\nu }+p{g}_{\mu \nu }$, together with its continuity equation, supplies the matter content of the Friedmann equations of 13.08.01, where the equation of state $p=w\rho$ controls the expansion history of the universe.

The vanishing of $T_{\mu \nu }$ outside a static star defines the vacuum region whose spherically symmetric solution is the Schwarzschild geometry of 13.05.01; the interior solution requires the fluid tensor matched across the stellar surface.

Historical & philosophical context Master

David Hilbert, in his November 1915 communication to the Göttingen Academy, obtained the gravitational field equations from a variational principle and identified the matter source through the metric variation of the matter action, the construction now bearing his name ^{[Hilbert 1915]}. The prescription gives a symmetric, gauge-invariant tensor directly, in contrast to the canonical tensor that Noether's first theorem produces from translation invariance, which for fields with spin is asymmetric and, for gauge fields, gauge-dependent. The problem of symmetrizing the canonical tensor was resolved independently in 1940 by Frederik Belinfante and by Léon Rosenfeld, who constructed the improvement term built from the spin current that converts the canonical tensor into a symmetric one while preserving its conservation ^{[Belinfante 1940] [Rosenfeld 1940]}. Rosenfeld's analysis established that the symmetrized canonical tensor coincides with the Hilbert tensor obtained by coupling the field to a metric and varying, which is why the curved-space construction is regarded as the physically primary definition.

Sternberg's treatment situates the energy-momentum tensor within the geometry of invariant variational problems, deriving its conservation as a Noether identity of the diffeomorphism symmetry ^{[Sternberg Ch. 20]}. The perfect-fluid form as a phenomenological stress-energy tensor predates general relativity, entering the relativistic literature through the early work on relativistic hydrodynamics and codified in the textbook treatments of Misner, Thorne, and Wheeler ^{[MTW §21.3]}. The conformal trace condition, classical at the level treated here, becomes the conformal anomaly in the quantum theory, a result of the 1970s that remains central to the study of conformal field theory and the renormalization of fields in curved spacetime.

Bibliography Master

@article{Hilbert1915,
  author  = {Hilbert, David},
  title   = {Die Grundlagen der Physik (Erste Mitteilung)},
  journal = {Nachrichten von der K\"oniglichen Gesellschaft der Wissenschaften zu G\"ottingen, Mathematisch-Physikalische Klasse},
  year    = {1915},
  pages   = {395--407}
}

@article{Rosenfeld1940,
  author  = {Rosenfeld, L\'eon},
  title   = {Sur le tenseur d'impulsion-\'energie},
  journal = {M\'emoires de l'Acad\'emie Royale de Belgique, Classe des Sciences},
  volume  = {18},
  number  = {6},
  year    = {1940},
  pages   = {1--30}
}

@article{Belinfante1940,
  author  = {Belinfante, Frederik J.},
  title   = {On the current and the density of the electric charge, the energy, the linear momentum and the angular momentum of arbitrary fields},
  journal = {Physica},
  volume  = {7},
  number  = {5},
  year    = {1940},
  pages   = {449--474}
}

@book{Wald1984,
  author    = {Wald, Robert M.},
  title     = {General Relativity},
  publisher = {University of Chicago Press},
  year      = {1984}
}

@book{MTW1973,
  author    = {Misner, Charles W. and Thorne, Kip S. and Wheeler, John A.},
  title     = {Gravitation},
  publisher = {W. H. Freeman},
  year      = {1973}
}

@book{Sternberg2012,
  author    = {Sternberg, Shlomo},
  title     = {Curvature in Mathematics and Physics},
  publisher = {Dover Publications},
  year      = {2012}
}

@book{Carroll2004,
  author    = {Carroll, Sean M.},
  title     = {Spacetime and Geometry: An Introduction to General Relativity},
  publisher = {Addison-Wesley},
  year      = {2004}
}