Renormalized Stress-Energy Tensor (RSET)

• RSET is the finite expectation value of the stress-energy tensor after subtracting universal ultraviolet singularities, capturing quantum backreaction in curved spacetimes.
• It is computed using methods like adiabatic regularization, Hadamard renormalization, and point-splitting, ensuring covariance, conservation, and compliance with the Hadamard structure.
• RSET plays a critical role in phenomena such as vacuum polarization, particle creation, trace anomalies, and the semiclassical equilibrium of black holes and cosmological models.

Searching arXiv for the cited papers to ground the article in current literature. arXiv search: (Rio et al., 2024) "Adiabatic regularization for spin 1/2 fields and the renormalized stress-energy tensor" — found. arXiv search: (Numajiri et al., 2024) "Boulware Vacuum vs. Regularity: Thoughts on Anomaly-Induced Effective Action" — found. arXiv search: (Taylor et al., 2022) "A mode-sum prescription for the renormalized stress energy tensor on black hole spacetimes" — found. arXiv search: (Arrechea et al., 2024) "The renormalized stress-energy tensor for scalar fields in the Boulware state with applications to extremal black holes" — found. The renormalized stress-energy tensor (RSET), denoted $\langle T_{\mu\nu}\rangle_{\mathrm{ren}}$, is the finite expectation value of the stress-energy tensor operator of a quantum field on a classical curved spacetime after subtraction of state-independent ultraviolet singularities. It is the quantity that appears on the right-hand side of the semiclassical Einstein equations and therefore encodes quantum backreaction, vacuum polarization, particle creation, and the imprint of the chosen quantum state. Across quantum field theory in curved spacetime, the RSET is defined so as to be finite, covariantly conserved, compatible with the Hadamard short-distance structure, and subject to the standard local curvature ambiguities of renormalization (Rio et al., 2024).

1. Definition and semiclassical role

In semiclassical gravity the metric is classical while matter is quantized, and the basic dynamical equation takes the form

$$G_{\mu\nu} + \text{(curvature counterterms)} = 8\pi G\,\langle T_{\mu\nu}\rangle_{\mathrm{ren}},$$

or, in a more general renormalized form,

$$G_{\mu\nu}+\Lambda\,g_{\mu\nu}+\alpha_{1} H^{(1)}_{\mu\nu}+\alpha_{2}H^{(2)}_{\mu\nu}=8\pi\,\langle\hat T_{\mu\nu}\rangle_{\mathrm{ren}}.$$

These forms make explicit that renormalization of the matter sector is accompanied by renormalization of the cosmological constant, Newton’s constant, and higher-curvature couplings (Rio et al., 2024).

For concrete fields, the unrenormalized object is obtained from the classical stress tensor by quantization and state expectation value. For a Dirac field in curved spacetime,

$$T_{\mu\nu}=\frac{i}{2}\left[\bar{\psi}\,\gamma_{a}e^a{}_{(\mu}\nabla_{\nu)}\psi-\bar{\psi}\overleftarrow{\nabla}_{(\nu} e^a{}_{\mu)}\gamma_{a}\psi\right].$$

In spatially flat FLRW spacetime, homogeneity and isotropy reduce the vacuum expectation value to an energy density and isotropic pressure,

$$\langle T_{00}\rangle = \frac{1}{2\pi^2 a^3}\int_{0}^{\infty} dk\,k^2\,\rho_k(t),\qquad \langle T_{ii}\rangle = \frac{1}{2\pi^2 a}\int_{0}^{\infty} dk\,k^2\,p_k(t),$$

with mode functions carrying the information about particle creation and vacuum polarization (Rio et al., 2024).

For scalar fields, the same conceptual structure appears in static black holes, wormholes, cosmic strings, and collapsing geometries. In each case the RSET is the finite tensor obtained after removing the local Hadamard or DeWitt-Schwinger singular part of the two-point function, and it is the object used to test energy conditions, determine fluxes, and source backreaction (Taylor et al., 2022).

2. Divergences, renormalization conditions, and ambiguities

The need for renormalization is ultraviolet. In FLRW, the bare Dirac-mode integrals diverge because for large comoving momentum $k$ the modes approach their Minkowski form and the integrands behave as positive powers of $k$. In point-splitting language, the coincidence limit of the two-point function is singular. In Hadamard form, the ultraviolet divergences reside solely in the universal geometric part of the two-point function, while the state-dependent remainder is smooth (Rio et al., 2024).

The renormalized tensor is required to satisfy a standard set of physical and mathematical conditions. The data emphasize the following requirements: local covariance, conservation,

$$\nabla^\mu \langle T_{\mu\nu}\rangle_{\mathrm{ren}}=0,$$

the correct flat-space limit, and agreement with the short-distance singularity structure prescribed by Hadamard form or equivalent constructions. In the axiomatic formulation, differences between any two admissible RSET definitions are local curvature tensors (Rio et al., 2024).

In four dimensions, the finite ambiguity has the form

$$c_1\, {^{(1)}H_{\mu\nu}} + c_2\, {^{(2)}H_{\mu\nu}} + c_3\, m^2 G_{\mu\nu}+ c_4\, m^4 g_{\mu\nu},$$

with $^{(1)}H_{\mu\nu}$ and $$G_{\mu\nu} + \text{(curvature counterterms)} = 8\pi G\,\langle T_{\mu\nu}\rangle_{\mathrm{ren}},$$0 obtained by functional differentiation of $$G_{\mu\nu} + \text{(curvature counterterms)} = 8\pi G\,\langle T_{\mu\nu}\rangle_{\mathrm{ren}},$$1 and $$G_{\mu\nu} + \text{(curvature counterterms)} = 8\pi G\,\langle T_{\mu\nu}\rangle_{\mathrm{ren}},$$2. These terms correspond respectively to renormalizations of quadratic-curvature couplings, Newton’s constant, and the cosmological constant (Rio et al., 2024).

A central quantum effect built into the RSET is the trace anomaly. Even when the classical tensor is traceless, the renormalized expectation value acquires a nonzero trace fixed by local curvature invariants. For Dirac fields in curved spacetime, the adiabatic construction reproduces the standard trace anomaly in the massless limit, and for two-dimensional conformal matter one has

$$G_{\mu\nu} + \text{(curvature counterterms)} = 8\pi G\,\langle T_{\mu\nu}\rangle_{\mathrm{ren}},$$3

in the anomaly-induced formulation (Numajiri et al., 2024).

3. Principal computational frameworks

Several renormalization frameworks recur in the literature, and the modern theory treats them as complementary rather than competing.

For spin-$$G_{\mu\nu} + \text{(curvature counterterms)} = 8\pi G\,\langle T_{\mu\nu}\rangle_{\mathrm{ren}},$$4 fields in FLRW spacetime, adiabatic regularization is a mode-by-mode subtraction scheme based on a large-$$G_{\mu\nu} + \text{(curvature counterterms)} = 8\pi G\,\langle T_{\mu\nu}\rangle_{\mathrm{ren}},$$5 expansion of the exact solutions. One introduces an adiabatic ansatz for the Dirac modes, expands the frequency and amplitudes in derivatives of the scale factor, and subtracts the divergent adiabatic orders. In four spacetime dimensions, the ultraviolet divergences of the Dirac energy density and pressure are completely captured by adiabatic orders $$G_{\mu\nu} + \text{(curvature counterterms)} = 8\pi G\,\langle T_{\mu\nu}\rangle_{\mathrm{ren}},$$6, $$G_{\mu\nu} + \text{(curvature counterterms)} = 8\pi G\,\langle T_{\mu\nu}\rangle_{\mathrm{ren}},$$7, and $$G_{\mu\nu} + \text{(curvature counterterms)} = 8\pi G\,\langle T_{\mu\nu}\rangle_{\mathrm{ren}},$$8, so the renormalized quantities are defined by

$$G_{\mu\nu} + \text{(curvature counterterms)} = 8\pi G\,\langle T_{\mu\nu}\rangle_{\mathrm{ren}},$$9

$$G_{\mu\nu}+\Lambda\,g_{\mu\nu}+\alpha_{1} H^{(1)}_{\mu\nu}+\alpha_{2}H^{(2)}_{\mu\nu}=8\pi\,\langle\hat T_{\mu\nu}\rangle_{\mathrm{ren}}.$$0

The subtraction terms are purely geometric, the resulting integrals are ultraviolet convergent, and the tensor is conserved (Rio et al., 2024).

The adiabatic scheme is explicitly equivalent to DeWitt-Schwinger point-splitting for Dirac fields in FLRW. The proof exploits the relation

$$G_{\mu\nu}+\Lambda\,g_{\mu\nu}+\alpha_{1} H^{(1)}_{\mu\nu}+\alpha_{2}H^{(2)}_{\mu\nu}=8\pi\,\langle\hat T_{\mu\nu}\rangle_{\mathrm{ren}}.$$1

together with the coincidence of the fourth-order adiabatic two-point function and the DeWitt-Schwinger two-point function truncated to four derivatives of the metric. The conclusion is

$$G_{\mu\nu}+\Lambda\,g_{\mu\nu}+\alpha_{1} H^{(1)}_{\mu\nu}+\alpha_{2}H^{(2)}_{\mu\nu}=8\pi\,\langle\hat T_{\mu\nu}\rangle_{\mathrm{ren}}.$$2

This places adiabatic regularization within the same covariant uniqueness class as point-splitting (Rio et al., 2024).

For scalar fields on black-hole backgrounds, Hadamard renormalization is often implemented numerically through mode sums. One line of development rewrites the Hadamard parametrix itself as a compatible mode sum using extended coordinates, allowing direct subtraction mode by mode. This was developed for Hartle-Hawking states and later generalized to the Boulware state, including the extremal Reissner-Nordström case where no Hartle-Hawking state exists (Taylor et al., 2022). Another line is pragmatic mode-sum regularization, which performs point splitting along a Killing direction, expands the Christensen counterterm in the same basis functions as the field modes, subtracts the singular mode contributions, and then removes residual oscillations with self-cancellation. This framework has been used for Schwarzschild, Reissner-Nordström, and Kerr, and it works directly in Lorentzian signature (Levi, 2016).

A distinct but related formalism is the anomaly-induced effective action. In two dimensions it introduces an auxiliary scalar $$G_{\mu\nu}+\Lambda\,g_{\mu\nu}+\alpha_{1} H^{(1)}_{\mu\nu}+\alpha_{2}H^{(2)}_{\mu\nu}=8\pi\,\langle\hat T_{\mu\nu}\rangle_{\mathrm{ren}}.$$3 satisfying $$G_{\mu\nu}+\Lambda\,g_{\mu\nu}+\alpha_{1} H^{(1)}_{\mu\nu}+\alpha_{2}H^{(2)}_{\mu\nu}=8\pi\,\langle\hat T_{\mu\nu}\rangle_{\mathrm{ren}}.$$4, while in four dimensions it introduces two auxiliary scalars obeying fourth-order equations driven by the Euler density and Weyl-squared invariant. Functional differentiation of the local auxiliary-field action yields a conserved RSET with the correct anomalous trace, and the homogeneous parts of the auxiliary fields encode the quantum state (Numajiri et al., 2024).

In a separate exact comparison for a massless scalar field in $$G_{\mu\nu}+\Lambda\,g_{\mu\nu}+\alpha_{1} H^{(1)}_{\mu\nu}+\alpha_{2}H^{(2)}_{\mu\nu}=8\pi\,\langle\hat T_{\mu\nu}\rangle_{\mathrm{ren}}.$$5 dimensions, the RSET obtained from normal-mode quantization was shown to coincide with the RSET derived from the Polyakov effective action once the vacuum-state choices are related by the map $$G_{\mu\nu}+\Lambda\,g_{\mu\nu}+\alpha_{1} H^{(1)}_{\mu\nu}+\alpha_{2}H^{(2)}_{\mu\nu}=8\pi\,\langle\hat T_{\mu\nu}\rangle_{\mathrm{ren}}.$$6, where $$G_{\mu\nu}+\Lambda\,g_{\mu\nu}+\alpha_{1} H^{(1)}_{\mu\nu}+\alpha_{2}H^{(2)}_{\mu\nu}=8\pi\,\langle\hat T_{\mu\nu}\rangle_{\mathrm{ren}}.$$7 is the timelike conformal Killing vector defining positive frequency. This makes explicit how state information is encoded differently but equivalently in operator and effective-action approaches (Barceló et al., 2011).

4. State dependence, regularity, and horizon structure

The RSET is state-dependent in a mathematically precise and physically consequential sense. In black-hole spacetimes, the Boulware, Unruh, and Hartle-Hawking states produce different horizon behavior and different asymptotic fluxes. In Schwarzschild, the Unruh state describes an evaporating black hole and yields a nonzero flux

$$G_{\mu\nu}+\Lambda\,g_{\mu\nu}+\alpha_{1} H^{(1)}_{\mu\nu}+\alpha_{2}H^{(2)}_{\mu\nu}=8\pi\,\langle\hat T_{\mu\nu}\rangle_{\mathrm{ren}}.$$8

while the Boulware state is vacuum-like at infinity and divergent on the horizon (Levi et al., 2016).

The relation between regularity and vacuum choice is especially sharp in static horizonless spacetimes. In two dimensions, because every static metric is conformally flat, a Boulware-like state can be both regular at the center and asymptotically Minkowski. In four dimensions, by contrast, the anomaly-induced analysis of a horizonless Bardeen-type spacetime finds that the preferred regular vacuum is not the Boulware vacuum but a nontrivial state with a different RSET profile, including asymptotic terms decaying as $$G_{\mu\nu}+\Lambda\,g_{\mu\nu}+\alpha_{1} H^{(1)}_{\mu\nu}+\alpha_{2}H^{(2)}_{\mu\nu}=8\pi\,\langle\hat T_{\mu\nu}\rangle_{\mathrm{ren}}.$$9 or $$T_{\mu\nu}=\frac{i}{2}\left[\bar{\psi}\,\gamma_{a}e^a{}_{(\mu}\nabla_{\nu)}\psi-\bar{\psi}\overleftarrow{\nabla}_{(\nu} e^a{}_{\mu)}\gamma_{a}\psi\right].$$0 rather than the Schwarzschild Boulware $$T_{\mu\nu}=\frac{i}{2}\left[\bar{\psi}\,\gamma_{a}e^a{}_{(\mu}\nabla_{\nu)}\psi-\bar{\psi}\overleftarrow{\nabla}_{(\nu} e^a{}_{\mu)}\gamma_{a}\psi\right].$$1 behavior (Numajiri et al., 2024).

Extremal horizons provide a contrasting example. A direct mode-sum computation of the scalar-field Boulware-state RSET in extremal Reissner-Nordström finds numerical evidence for regularity of the full tensor at the extremal horizon regardless of the field mass and coupling (Arrechea et al., 2024). For the conformally coupled massless field at the horizon,

$$T_{\mu\nu}=\frac{i}{2}\left[\bar{\psi}\,\gamma_{a}e^a{}_{(\mu}\nabla_{\nu)}\psi-\bar{\psi}\overleftarrow{\nabla}_{(\nu} e^a{}_{\mu)}\gamma_{a}\psi\right].$$2

and for massive fields the analytic and numerical pieces separately diverge logarithmically but cancel in the full RSET, leaving finite energy density for freely falling observers (Arrechea et al., 2024).

These examples underscore a general point: regularity is not determined by geometry alone. It is a joint property of geometry, renormalization prescription, and quantum state.

5. Representative regimes and physical behavior

The literature records several characteristic regimes in which the RSET acquires distinctive physical content. The table summarizes examples explicitly discussed in the cited works.

Setting	Field/state	Salient RSET behavior
Spatially flat FLRW	Dirac vacuum	Adiabatic subtraction to fourth order yields finite, conserved $$T_{\mu\nu}=\frac{i}{2}\left[\bar{\psi}\,\gamma_{a}e^a{}_{(\mu}\nabla_{\nu)}\psi-\bar{\psi}\overleftarrow{\nabla}_{(\nu} e^a{}_{\mu)}\gamma_{a}\psi\right].$$3
Radiation-dominated FLRW	Dirac vacuum with Whittaker modes	Early-time $$T_{\mu\nu}=\frac{i}{2}\left[\bar{\psi}\,\gamma_{a}e^a{}_{(\mu}\nabla_{\nu)}\psi-\bar{\psi}\overleftarrow{\nabla}_{(\nu} e^a{}_{\mu)}\gamma_{a}\psi\right].$$4 radiation-like behavior; late-time $$T_{\mu\nu}=\frac{i}{2}\left[\bar{\psi}\,\gamma_{a}e^a{}_{(\mu}\nabla_{\nu)}\psi-\bar{\psi}\overleftarrow{\nabla}_{(\nu} e^a{}_{\mu)}\gamma_{a}\psi\right].$$5 matter-like behavior
Schwarzschild	Unruh state, minimally coupled massless scalar	WEC and NEC violated from $$T_{\mu\nu}=\frac{i}{2}\left[\bar{\psi}\,\gamma_{a}e^a{}_{(\mu}\nabla_{\nu)}\psi-\bar{\psi}\overleftarrow{\nabla}_{(\nu} e^a{}_{\mu)}\gamma_{a}\psi\right].$$6 to the horizon; ANEC violated on the circular null geodesic at $$T_{\mu\nu}=\frac{i}{2}\left[\bar{\psi}\,\gamma_{a}e^a{}_{(\mu}\nabla_{\nu)}\psi-\bar{\psi}\overleftarrow{\nabla}_{(\nu} e^a{}_{\mu)}\gamma_{a}\psi\right].$$7
Extremal Reissner-Nordström	Boulware state, scalar field	Numerical evidence for RSET regularity at the extremal horizon
Zero-tidal wormhole	Vacuum of non-minimally coupled massive scalar	Three disconnected supporting regions in $$T_{\mu\nu}=\frac{i}{2}\left[\bar{\psi}\,\gamma_{a}e^a{}_{(\mu}\nabla_{\nu)}\psi-\bar{\psi}\overleftarrow{\nabla}_{(\nu} e^a{}_{\mu)}\gamma_{a}\psi\right].$$8 and two mass exclusion intervals
Nonzero-radius cosmic string	Scalar field vacuum	Smooth core, asymptotic approach to point-string $$T_{\mu\nu}=\frac{i}{2}\left[\bar{\psi}\,\gamma_{a}e^a{}_{(\mu}\nabla_{\nu)}\psi-\bar{\psi}\overleftarrow{\nabla}_{(\nu} e^a{}_{\mu)}\gamma_{a}\psi\right].$$9 behavior, explicit conformal anomaly

In cosmology, the Dirac-field RSET in a radiation-dominated FLRW universe interpolates between

$$\langle T_{00}\rangle = \frac{1}{2\pi^2 a^3}\int_{0}^{\infty} dk\,k^2\,\rho_k(t),\qquad \langle T_{ii}\rangle = \frac{1}{2\pi^2 a}\int_{0}^{\infty} dk\,k^2\,p_k(t),$$0

at early times and

$$\langle T_{00}\rangle = \frac{1}{2\pi^2 a^3}\int_{0}^{\infty} dk\,k^2\,\rho_k(t),\qquad \langle T_{ii}\rangle = \frac{1}{2\pi^2 a}\int_{0}^{\infty} dk\,k^2\,p_k(t),$$1

at late times, for adiabatically renormalizable vacua satisfying

$$\langle T_{00}\rangle = \frac{1}{2\pi^2 a^3}\int_{0}^{\infty} dk\,k^2\,\rho_k(t),\qquad \langle T_{ii}\rangle = \frac{1}{2\pi^2 a}\int_{0}^{\infty} dk\,k^2\,p_k(t),$$2

The RSET therefore reproduces radiation-like behavior before mass effects dominate and matter-like behavior afterwards (Rio et al., 2024).

In black-hole physics, the RSET governs both local energy conditions and asymptotic fluxes. For a minimally coupled massless scalar in Schwarzschild Unruh state, the tensor is type I at weak field and becomes type IV at $$\langle T_{00}\rangle = \frac{1}{2\pi^2 a^3}\int_{0}^{\infty} dk\,k^2\,\rho_k(t),\qquad \langle T_{ii}\rangle = \frac{1}{2\pi^2 a}\int_{0}^{\infty} dk\,k^2\,p_k(t),$$3; both WEC and NEC are violated throughout $$\langle T_{00}\rangle = \frac{1}{2\pi^2 a^3}\int_{0}^{\infty} dk\,k^2\,\rho_k(t),\qquad \langle T_{ii}\rangle = \frac{1}{2\pi^2 a}\int_{0}^{\infty} dk\,k^2\,p_k(t),$$4, and the averaged null energy condition is violated on the circular null geodesic at $$\langle T_{00}\rangle = \frac{1}{2\pi^2 a^3}\int_{0}^{\infty} dk\,k^2\,\rho_k(t),\qquad \langle T_{ii}\rangle = \frac{1}{2\pi^2 a}\int_{0}^{\infty} dk\,k^2\,p_k(t),$$5 (Levi et al., 2016). By contrast, in a Hartle-Hawking computation for a scalar field in Schwarzschild, the null energy density on the photon sphere is positive for the parameter ranges examined in the paper (Taylor et al., 2022).

In wormhole spacetimes, the RSET is examined as a source of exotic matter. For a non-minimally coupled massive scalar field in a zero-tidal Morris-Thorne wormhole, the renormalized quantities at the throat are

$$\langle T_{00}\rangle = \frac{1}{2\pi^2 a^3}\int_{0}^{\infty} dk\,k^2\,\rho_k(t),\qquad \langle T_{ii}\rangle = \frac{1}{2\pi^2 a}\int_{0}^{\infty} dk\,k^2\,p_k(t),$$6

and the wormhole-supporting conditions are $$\langle T_{00}\rangle = \frac{1}{2\pi^2 a^3}\int_{0}^{\infty} dk\,k^2\,\rho_k(t),\qquad \langle T_{ii}\rangle = \frac{1}{2\pi^2 a}\int_{0}^{\infty} dk\,k^2\,p_k(t),$$7 and $$\langle T_{00}\rangle = \frac{1}{2\pi^2 a^3}\int_{0}^{\infty} dk\,k^2\,\rho_k(t),\qquad \langle T_{ii}\rangle = \frac{1}{2\pi^2 a}\int_{0}^{\infty} dk\,k^2\,p_k(t),$$8. The analysis finds three disconnected regions in $$\langle T_{00}\rangle = \frac{1}{2\pi^2 a^3}\int_{0}^{\infty} dk\,k^2\,\rho_k(t),\qquad \langle T_{ii}\rangle = \frac{1}{2\pi^2 a}\int_{0}^{\infty} dk\,k^2\,p_k(t),$$9 where the RSET supports the throat and two mass exclusion intervals in which no value of $k$0 suffices (Jiang et al., 1 Mar 2026).

In static compact objects, vacuum polarization in the Boulware state can produce negative-energy cores. Using an order-reduced approximation to the scalar-field RSET, ultracompact horizonless stars with compactness beyond the Buchdahl limit were found, with a negative Misner-Sharp mass interior generated by the negative quantum energy density near the center (Arrechea et al., 2023).

6. Theoretical status and semiclassical implications

The modern status of the RSET is that of a locally constrained but state-sensitive observable, computable in several equivalent schemes and indispensable for semiclassical backreaction. Equivalence results are particularly important: adiabatic and DeWitt-Schwinger renormalization coincide for Dirac fields in FLRW (Rio et al., 2024), and in $k$1 dimensions the normal-mode and Polyakov-effective-action constructions yield the same scalar-field RSET once their vacuum choices are mapped appropriately (Barceló et al., 2011). These results reinforce the view that admissible renormalization schemes differ only by the standard finite local curvature terms.

At the same time, the RSET can qualitatively alter spacetime structure. In extremal Reissner-Nordström, static semiclassical perturbations sourced by the Boulware-state scalar RSET either de-extremalize the black hole or convert it into a horizonless object (Arrechea et al., 2024). In horizonless ultracompact stars, negative vacuum energy in the Boulware state can invalidate the hypotheses of Buchdahl’s theorem and support equilibrium configurations that remain regular while exceeding the classical compactness bound (Arrechea et al., 2023).

The RSET is therefore not merely a renormalized local observable. It is the operational bridge between quantum matter and classical geometry: finite only after subtraction of universal short-distance structure, constrained by conservation and anomaly, sensitive to the global quantum state, and decisive in questions of evaporation, regularity, energy-condition violation, and semiclassical equilibrium.