# Stress–energy tensor

__: Tensor describing the density and flux of energy in spacetime__

**Short description**General relativity |
---|

The **stress–energy tensor**, sometimes called the **stress–energy–momentum tensor** or the **energy–momentum tensor**, is a tensor quantity in physics that describes the density and flux of energy and momentum in spacetime, generalizing the stress tensor of Newtonian physics. It is an attribute of matter, radiation, and non-gravitational force fields. This density and flux of energy and momentum are the sources of the gravitational field in the Einstein field equations of general relativity, just as mass density is the source of such a field in Newtonian gravity.

## Definition

The stress–energy tensor involves the use of superscripted variables (*not* exponents; see tensor index notation and Einstein summation notation). If Cartesian coordinates in SI units are used, then the components of the position four-vector are given by: *x*^{0} = *t*, *x*^{1} = *x*, *x*^{2} = *y*, and *x*^{3} = *z*, where *t* is time in seconds, and *x*, *y*, and *z* are distances in meters.

The stress–energy tensor is defined as the tensor *T*^{αβ} of order two that gives the flux of the *α*th component of the momentum vector across a surface with constant *x*^{β} coordinate. In the theory of relativity, this momentum vector is taken as the four-momentum. In general relativity, the stress–energy tensor is symmetric,^{[1]}

- [math]\displaystyle{ T^{\alpha \beta} = T^{\beta \alpha}. }[/math]

In some alternative theories like Einstein–Cartan theory, the stress–energy tensor may not be perfectly symmetric because of a nonzero spin tensor, which geometrically corresponds to a nonzero torsion tensor.

## Identifying the components of the tensor

Because the stress–energy tensor is of order two, its components can be displayed in 4 × 4 matrix form:

- [math]\displaystyle{ (T^{\mu\nu})_{\mu,\nu=0,1,2,3} = \begin{pmatrix} T^{00} & T^{01} & T^{02} & T^{03} \\ T^{10} & T^{11} & T^{12} & T^{13} \\ T^{20} & T^{21} & T^{22} & T^{23} \\ T^{30} & T^{31} & T^{32} & T^{33} \end{pmatrix}. }[/math]

In the following, *i* and *k* range from 1 through 3.

The time–time component is the density of relativistic mass, i.e., the energy density divided by the speed of light squared.^{[2]} Its components have a direct physical interpretation. In the case of a perfect fluid this component is

- [math]\displaystyle{ T^{00} = \rho, }[/math]

where [math]\displaystyle{ \rho }[/math] is the relativistic mass per unit volume, and for an electromagnetic field in otherwise empty space this component is

- [math]\displaystyle{ T^{00} = {1 \over c^2}\left(\frac{1}{2} \epsilon_0 E^2 + \frac{1}{2\mu_0} B^2 \right), }[/math]

where *E* and *B* are the electric and magnetic fields, respectively.^{[3]}

The flux of relativistic mass across the *x ^{i}* surface is equivalent to the density of the

*i*th component of linear momentum,

- [math]\displaystyle{ T^{0i} = T^{i0}. }[/math]

The components

- [math]\displaystyle{ T^{ik} }[/math]

represent flux of *i*th component of linear momentum across the *x ^{k}* surface. In particular,

- [math]\displaystyle{ T^{ii} }[/math]

(not summed) represents normal stress, which is called pressure when it is independent of direction. The remaining components

- [math]\displaystyle{ T^{ik} \quad i \ne k }[/math]

represent shear stress (compare with the stress tensor).

In solid state physics and fluid mechanics, the stress tensor is defined to be the spatial components of the stress–energy tensor in the proper frame of reference. In other words, the stress energy tensor in engineering *differs* from the relativistic stress–energy tensor by a momentum convective term.

### Covariant and mixed forms

Most of this article works with the contravariant form, *T*^{μν} of the stress–energy tensor. However, it is often necessary to work with the covariant form,

- [math]\displaystyle{ T_{\mu \nu} = T^{\alpha \beta} g_{\alpha \mu} g_{\beta \nu}, }[/math]

or the mixed form,

- [math]\displaystyle{ T^\mu{}_\nu = T^{\mu \alpha} g_{\alpha \nu}, }[/math]

or as a mixed tensor density

- [math]\displaystyle{ \mathfrak{T}^\mu{}_\nu = T^\mu{}_\nu \sqrt{-g} \,. }[/math]

This article uses the spacelike sign convention (−+++) for the metric signature.

## Conservation law

### In special relativity

The stress–energy tensor is the conserved Noether current associated with spacetime translations.

The divergence of the non-gravitational stress–energy is zero. In other words, non-gravitational energy and momentum are conserved,

- [math]\displaystyle{ 0 = T^{\mu \nu}{}_{;\nu} = \nabla_\nu T^{\mu \nu}{}. \! }[/math]

When gravity is negligible and using a Cartesian coordinate system for spacetime, this may be expressed in terms of partial derivatives as

- [math]\displaystyle{ 0 = T^{\mu \nu}{}_{,\nu} = \partial_{\nu} T^{\mu \nu}. \! }[/math]

The integral form of this is

- [math]\displaystyle{ 0 = \int_{\partial N} T^{\mu \nu} \mathrm{d}^3 s_{\nu} \! }[/math]

where *N* is any compact four-dimensional region of spacetime; [math]\displaystyle{ \partial N }[/math] is its boundary, a three-dimensional hypersurface; and [math]\displaystyle{ \mathrm{d}^3 s_{\nu} }[/math] is an element of the boundary regarded as the outward pointing normal.

In flat spacetime and using Cartesian coordinates, if one combines this with the symmetry of the stress–energy tensor, one can show that angular momentum is also conserved:

- [math]\displaystyle{ 0 = (x^{\alpha} T^{\mu \nu} - x^{\mu} T^{\alpha \nu})_{,\nu} . \! }[/math]

### In general relativity

When gravity is non-negligible or when using arbitrary coordinate systems, the divergence of the stress–energy still vanishes. But in this case, a coordinate-free definition of the divergence is used which incorporates the covariant derivative

- [math]\displaystyle{ 0 = \operatorname{div} T = T^{\mu \nu}{}_{;\nu} = \nabla_{\nu} T^{\mu \nu} = T^{\mu \nu}{}_{,\nu} + \Gamma^{\mu}{}_{\sigma \nu}T^{\sigma \nu} + \Gamma^{\nu}{}_{\sigma \nu} T^{\mu \sigma} }[/math]

where [math]\displaystyle{ \Gamma^{\mu}{}_{\sigma \nu} }[/math] is the Christoffel symbol which is the gravitational force field.

Consequently, if [math]\displaystyle{ \xi^{\mu} }[/math] is any Killing vector field, then the conservation law associated with the symmetry generated by the Killing vector field may be expressed as

- [math]\displaystyle{ 0 = \nabla_\nu (\xi^{\mu} T_{\mu}^{\nu}) = \frac{1}{\sqrt{-g}} \partial_\nu( \sqrt{-g} \ \xi^{\mu} T_{\mu}^{\nu}) }[/math]

The integral form of this is

- [math]\displaystyle{ 0 = \int_{\partial N} \sqrt{-g} \ \xi^{\mu} T_{\mu}^{\nu} \ \mathrm{d}^3 s_{\nu} = \int_{\partial N} \xi^{\mu} \mathfrak{T}_{\mu}^{\nu} \ \mathrm{d}^3 s_{\nu} }[/math]

## In special relativity

In special relativity, the stress–energy tensor contains information about the energy and momentum densities of a given system, in addition to the momentum and energy flux densities.^{[4]}

Given a Lagrangian Density [math]\displaystyle{ \mathcal{L} }[/math] that is a function of a set of fields [math]\displaystyle{ \phi_{\alpha} }[/math] and their derivatives, but explicitly not of any of the spacetime coordinates, we can construct the tensor by looking at the total derivative with respect to one of the generalized coordinates of the system. So, with our condition

- [math]\displaystyle{ \frac{\partial \mathcal{L}}{\partial x_{\nu}} = 0 }[/math]

By using the chain rule, we then have

- [math]\displaystyle{ \frac{d \mathcal{L}}{dx_{\nu}} = \partial^{\nu}\mathcal{L} = \frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})}\frac{\partial(\partial_{\mu}\phi_{\alpha})}{\partial x_{\nu}} + \frac{\partial \mathcal{L}}{\partial \phi_{\alpha}}\frac{\partial \phi_{\alpha}}{\partial x_{\nu}} }[/math]

Written in useful shorthand,

- [math]\displaystyle{ \partial^{\nu}\mathcal{L} = \frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})}\partial^{\nu}\partial_{\mu}\phi_{\alpha} + \frac{\partial \mathcal{L}}{\partial \phi_{\alpha}}\partial^{\nu}\phi_{\alpha} }[/math]

Then, we can use the Euler–Lagrange Equation:

- [math]\displaystyle{ \partial_{\mu}(\frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})}) = \frac{\partial\mathcal{L}}{\partial \phi_{\alpha}} }[/math]

And then use the fact that partial derivatives commute so that we now have

- [math]\displaystyle{ \partial^{\nu}\mathcal{L} = \frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})}\partial_{\mu}\partial^{\nu}\phi_{\alpha} + \partial_{\mu}(\frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})})\partial^{\nu}\phi_{\alpha} }[/math]

We can recognize the right hand side as a product rule. Writing it as the derivative of a product of functions tells us that

- [math]\displaystyle{ \partial^{\nu}\mathcal{L} = \partial_{\mu}[\frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})}\partial^{\nu}\phi_{\alpha}] }[/math]

Now, in flat space, one can write [math]\displaystyle{ \partial^{\nu}\mathcal{L} = \partial_{\mu}g^{\mu\nu}\mathcal{L} }[/math]. Doing this and moving it to the other side of the equation tells us that

- [math]\displaystyle{ \partial_{\mu}[\frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})}\partial^{\nu}\phi_{\alpha}] - \partial_{\mu}(g^{\mu\nu}\mathcal{L}) = 0 }[/math]

And upon regrouping terms,

- [math]\displaystyle{ \partial_{\mu}[\frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})}\partial^{\nu}\phi_{\alpha}-g^{\mu\nu}\mathcal{L}] = 0 }[/math]

This is to say that the divergence of the tensor in the brackets is 0. Indeed, with this, we define the stress–energy tensor:

- [math]\displaystyle{ T^{\mu\nu} \equiv \frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})}\partial^{\nu}\phi_{\alpha}-g^{\mu\nu}\mathcal{L} }[/math]

By construction it has the property that

- [math]\displaystyle{ \partial_{\mu}T^{\mu\nu} = 0 }[/math]

Note that this divergenceless property of this tensor is equivalent to four continuity equations. That is, fields have at least four sets of quantities that obey the continuity equation. As an example, it can be seen that [math]\displaystyle{ T^{0}_0 }[/math] is the energy density of the system and that it is thus possible to obtain the Hamiltonian density from the stress–energy tensor.

Indeed, since this is the case, observing that [math]\displaystyle{ \partial_{\mu}T^{\mu 0} = 0 }[/math], we then have

- [math]\displaystyle{ \frac{\partial \mathcal{H}}{\partial t} + \nabla\cdot(\frac{\partial \mathcal{L}}{\partial\nabla \phi_{\alpha}}\dot{\phi}_{\alpha}) = 0 }[/math]

We can then conclude that the terms of [math]\displaystyle{ \frac{\partial \mathcal{L}}{\partial\nabla \phi_{\alpha}}\dot{\phi}_{\alpha} }[/math] represent the energy flux density of the system.

### The trace

Note that the trace is defined to be [math]\displaystyle{ T^{\mu}_{\mu} }[/math]. Note that

- [math]\displaystyle{ T^{\mu}_{\mu} = T^{\mu\nu}g_{\nu\mu} . }[/math]

When we use our formula for the stress–energy tensor found above,

- [math]\displaystyle{ T^{\mu}_{\mu} = \frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})}g_{\mu\nu}\partial^{\nu}\phi_{\alpha}-g_{\mu\nu}g^{\mu\nu}\mathcal{L} . }[/math]

Using the raising and lowering properties of the metric, and the fact that [math]\displaystyle{ g^{\mu\nu}g_{\mu\alpha} = \delta^{\nu}_{\alpha} }[/math],

- [math]\displaystyle{ T^{\mu}_{\mu} = \frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})}\partial_{\mu}\phi_{\alpha}-\delta^{\mu}_{\mu}\mathcal{L} . }[/math]

Since [math]\displaystyle{ \delta^{\mu}_{\mu} = 4 }[/math], so that we can conclude that

- [math]\displaystyle{ T^{\mu}_{\mu} = \frac{\partial \mathcal{L}}{\partial(\partial_{\mu}\phi_{\alpha})}\partial_{\mu}\phi_{\alpha}-4\mathcal{L} . }[/math]

## In general relativity

In general relativity, the symmetric stress–energy tensor acts as the source of spacetime curvature, and is the current density associated with gauge transformations of gravity which are general curvilinear coordinate transformations. (If there is torsion, then the tensor is no longer symmetric. This corresponds to the case with a nonzero spin tensor in Einstein–Cartan gravity theory.)

In general relativity, the partial derivatives used in special relativity are replaced by covariant derivatives. What this means is that the continuity equation no longer implies that the non-gravitational energy and momentum expressed by the tensor are absolutely conserved, i.e. the gravitational field can do work on matter and vice versa. In the classical limit of Newtonian gravity, this has a simple interpretation: kinetic energy is being exchanged with gravitational potential energy, which is not included in the tensor, and momentum is being transferred through the field to other bodies. In general relativity the Landau–Lifshitz pseudotensor is a unique way to define the *gravitational* field energy and momentum densities. Any such stress–energy pseudotensor can be made to vanish locally by a coordinate transformation.

In curved spacetime, the spacelike integral now depends on the spacelike slice, in general. There is in fact no way to define a global energy–momentum vector in a general curved spacetime.

### The Einstein field equations

In general relativity, the stress tensor is studied in the context of the Einstein field equations which are often written as

- [math]\displaystyle{ R_{\mu \nu} - \tfrac{1}{2} R\,g_{\mu \nu} + \Lambda g_{\mu \nu} = {8 \pi G \over c^4} T_{\mu \nu}, }[/math]

where [math]\displaystyle{ R_{\mu \nu} }[/math] is the Ricci tensor, [math]\displaystyle{ R }[/math] is the Ricci scalar (the tensor contraction of the Ricci tensor), [math]\displaystyle{ g_{\mu \nu}\, }[/math] is the metric tensor, Λ is the cosmological constant (negligible at the scale of a galaxy or smaller), and [math]\displaystyle{ G }[/math] is the universal gravitational constant.

## Stress–energy in special situations

### Isolated particle

In special relativity, the stress–energy of a non-interacting particle with rest mass *m* and trajectory [math]\displaystyle{ \mathbf{x}_\text{p}(t) }[/math] is:

- [math]\displaystyle{ T^{\alpha \beta}(\mathbf{x},t) = \frac{m \, v^{\alpha}(t) v^{\beta}(t)}{\sqrt{1 - (v/c)^2}}\;\, \delta(\mathbf{x} - \mathbf{x}_\text{p}(t)) = \frac{E}{c^2}\; v^{\alpha}(t) v^{\beta}(t)\;\, \delta(\mathbf{x} - \mathbf{x}_\text{p}(t)) }[/math]

where [math]\displaystyle{ (v^{\alpha})_{\alpha=0,1,2,3} \! }[/math] is the velocity vector (which should not be confused with four-velocity, since it is missing a [math]\displaystyle{ \gamma }[/math])

- [math]\displaystyle{ (v^{\alpha})_{\alpha=0,1,2,3} = \left(1, \frac{d \mathbf{x}_\text{p}}{dt}(t) \right) \,, }[/math]

δ is the Dirac delta function and [math]\displaystyle{ E = \sqrt{p^2 c^2 + m^2 c^4} }[/math] is the energy of the particle.

### Stress–energy of a fluid in equilibrium

For a perfect fluid in thermodynamic equilibrium, the stress–energy tensor takes on a particularly simple form

- [math]\displaystyle{ T^{\alpha \beta} \, = \left(\rho + {p \over c^2}\right)u^{\alpha}u^{\beta} + p g^{\alpha \beta} }[/math]

where [math]\displaystyle{ \rho }[/math] is the mass–energy density (kilograms per cubic meter), [math]\displaystyle{ p }[/math] is the hydrostatic pressure (pascals), [math]\displaystyle{ u^{\alpha} }[/math] is the fluid's four velocity, and [math]\displaystyle{ g^{\alpha \beta} }[/math] is the reciprocal of the metric tensor. Therefore, the trace is given by

- [math]\displaystyle{ T = 3p - \rho c^2 \,. }[/math]

The four velocity satisfies

- [math]\displaystyle{ u^{\alpha} u^{\beta} g_{\alpha \beta} = - c^2 \,. }[/math]

In an inertial frame of reference comoving with the fluid, better known as the fluid's proper frame of reference, the four velocity is

- [math]\displaystyle{ (u^{\alpha})_{\alpha=0,1,2,3} = (1, 0, 0, 0) \,, }[/math]

the reciprocal of the metric tensor is simply

- [math]\displaystyle{ (g^{\alpha \beta})_{\alpha,\beta=0,1,2,3} \, = \left( \begin{matrix} - c^{-2} & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{matrix} \right) \, }[/math]

and the stress–energy tensor is a diagonal matrix

- [math]\displaystyle{ (T^{\alpha \beta})_{\alpha,\beta=0,1,2,3} = \left( \begin{matrix} \rho & 0 & 0 & 0 \\ 0 & p & 0 & 0 \\ 0 & 0 & p & 0 \\ 0 & 0 & 0 & p \end{matrix} \right). }[/math]

### Electromagnetic stress–energy tensor

The Hilbert stress–energy tensor of a source-free electromagnetic field is

- [math]\displaystyle{ T^{\mu \nu} = \frac{1}{\mu_0} \left( F^{\mu \alpha} g_{\alpha \beta} F^{\nu \beta} - \frac{1}{4} g^{\mu \nu} F_{\delta \gamma} F^{\delta \gamma} \right) }[/math]

where [math]\displaystyle{ F_{\mu \nu} }[/math] is the electromagnetic field tensor.

### Scalar field

The stress–energy tensor for a complex scalar field [math]\displaystyle{ \phi }[/math] which satisfies the Klein–Gordon equation is

- [math]\displaystyle{ T^{\mu\nu} = \frac{\hbar^2}{m} (g^{\mu \alpha} g^{\nu \beta} + g^{\mu \beta} g^{\nu \alpha} - g^{\mu\nu} g^{\alpha \beta}) \partial_{\alpha}\bar\phi \partial_{\beta}\phi - g^{\mu\nu} m c^2 \bar\phi \phi , }[/math]

and when the metric is flat (Minkowski) its components work out to be:

- [math]\displaystyle{ \begin{align} T^{00} & = \frac{\hbar^2}{m c^4} \left(\partial_0 \bar{\phi} \partial_0 \phi + c^2 \partial_k \bar{\phi} \partial_k \phi \right) + m \bar{\phi} \phi, \\ T^{0i} = T^{i0} & = - \frac{\hbar^2}{m c^2} \left(\partial_0 \bar{\phi} \partial_i \phi + \partial_i \bar{\phi} \partial_0 \phi \right),\ \mathrm{and} \\ T^{ij} & = \frac{\hbar^2}{m} \left(\partial_i \bar{\phi} \partial_j \phi + \partial_j \bar{\phi} \partial_i \phi \right) - \delta_{ij} \left(\frac{\hbar^2}{m} \eta^{\alpha\beta} \partial_\alpha \bar{\phi} \partial_\beta \phi + m c^2 \bar{\phi} \phi\right). \end{align} }[/math]

## Variant definitions of stress–energy

There are a number of inequivalent definitions of non-gravitational stress–energy:

### Hilbert stress–energy tensor

The Hilbert stress–energy tensor is defined as the functional derivative

- [math]\displaystyle{ T_{\mu\nu} = \frac{-2}{\sqrt{-g}}\frac{\delta S_{\mathrm{matter}} }{\delta g^{\mu\nu}} =\frac{-2}{\sqrt{-g}}\frac{\partial(\sqrt{-g}\mathcal{L}_{\mathrm{matter}})}{\partial g^{\mu\nu}}= - 2 \frac{\partial \mathcal{L}_\mathrm{matter}}{\partial g^{\mu\nu}} + g_{\mu\nu} \mathcal{L}_\mathrm{matter} , }[/math]

where [math]\displaystyle{ S_{\mathrm{matter}} }[/math] is the nongravitational part of the action, [math]\displaystyle{ \mathcal{L}_{\mathrm{matter}} }[/math] is the nongravitational part of the Lagrangian density, and the Euler-Lagrange equation has been used. This is symmetric and gauge-invariant. See Einstein–Hilbert action for more information.

### Canonical stress–energy tensor

Noether's theorem implies that there is a conserved current associated with translations through space and time. This is called the canonical stress–energy tensor. Generally, this is not symmetric and if we have some gauge theory, it may not be gauge invariant because space-dependent gauge transformations do not commute with spatial translations.

In general relativity, the translations are with respect to the coordinate system and as such, do not transform covariantly. See the section below on the gravitational stress–energy pseudo-tensor.

### Belinfante–Rosenfeld stress–energy tensor

In the presence of spin or other intrinsic angular momentum, the canonical Noether stress energy tensor fails to be symmetric. The Belinfante–Rosenfeld stress energy tensor is constructed from the canonical stress–energy tensor and the spin current in such a way as to be symmetric and still conserved. In general relativity, this modified tensor agrees with the Hilbert stress–energy tensor.

## Gravitational stress–energy

By the equivalence principle gravitational stress–energy will always vanish locally at any chosen point in some chosen frame, therefore gravitational stress–energy cannot be expressed as a non-zero tensor; instead we have to use a pseudotensor.

In general relativity, there are many possible distinct definitions of the gravitational stress–energy–momentum pseudotensor. These include the Einstein pseudotensor and the Landau–Lifshitz pseudotensor. The Landau–Lifshitz pseudotensor can be reduced to zero at any event in spacetime by choosing an appropriate coordinate system.

## See also

- Cooperstock's energy-localization hypothesis
- Electromagnetic stress–energy tensor
- Energy condition
- Energy density of electric and magnetic fields
- Maxwell stress tensor
- Poynting vector
- Ricci calculus
- Segre classification

## Notes and references

- ↑ On pp. 141–142 of Misner, Thorne, and Wheeler, section 5.7 "Symmetry of the Stress–Energy Tensor" begins with "All the stress–energy tensors explored above were symmetric. That they could not have been otherwise one sees as follows."
- ↑ Charles W., Misner, Thorne, Kip S., Wheeler, John A., (1973). Gravitation. San Frandisco: W. H. Freeman and Company. ISBN:0-7167-0334-3.
- ↑ d'Inverno, R.A, (1992). Introducing Einstein's Relativity. New York: Oxford University Press. ISBN:978-0-19-859686-8.
- ↑ Landau, L.D.; Lifshitz, E.M. (2010).
*The Classical Theory of Fields*(4th ed.). Butterworth-Heinemann. pp. 84-85. ISBN 978-0-7506-2768-9.

- W. Wyss (2005). "The Energy–Momentum Tensor in Classical Field Theory". Colorado, USA. http://merlin.fic.uni.lodz.pl/concepts/2005_3_4/2005_3_4_295.pdf.

## External links

- Lecture, Stephan Waner
- Caltech Tutorial on Relativity — A simple discussion of the relation between the Stress–Energy tensor of General Relativity and the metric