Hadamard Tail in Curved Spacetime

• Hadamard tail is the part of wave propagation inside the light cone in curved spacetime, capturing curvature-induced scattering effects and deviations from the strong Huygens principle.
• It is mathematically formulated as the smooth biscalar V(x,x') in the Hadamard representation and appears as a logarithmic term in quantum two-point functions.
• The concept underpins practical insights in self-force computations and quantum detector communications, highlighting singular behavior near caustic regions.

Searching arXiv for papers directly using the term “Hadamard tail” and closely related curved-spacetime/QFT usage. Hadamard tail is the part of wave propagation on curved spacetime that lies inside the light cone rather than exactly on it. In the local Hadamard representation of the retarded Green function, it is the smooth biscalar $V(x,x')$; in a closely related short-distance quantum-field-theoretic usage, it is the curvature-dependent non-analytic term $v(x_1,x_2)\log\!\big(s^2+i\epsilon\big)$ in the two-point function of a scalar field. In both formulations, the term records the failure of sharp null propagation in curved geometry and the consequent appearance of propagation or singular structure at timelike separation (Aruquipa et al., 2022, Chavda et al., 29 Oct 2025).

1. Definition through the Hadamard form

In curved spacetime, the retarded Green function is not supported only on the null cone. In flat $3+1$ dimensions, massless propagation obeys the strong Huygens principle, so disturbances propagate sharply along null geodesics. In curved spacetime, curvature scatters the field, and part of the signal propagates effectively at sub-light speeds and contributes inside the light cone. The local Hadamard representation makes this precise as (Aruquipa et al., 2022)

$$G_{\mathrm{ret}}(x,x')= \big[U(x,x')\delta(\sigma(x,x'))+V(x,x')\theta(-\sigma(x,x'))\big]\theta_+(x,x').$$

Here $\sigma(x,x')$ is Synge’s world function, defined as one-half the squared geodesic distance along the unique geodesic joining $x$ and $x'$ in a normal neighbourhood. The term $U(x,x')\delta(\sigma)$ has support on the light cone and is the direct or null part. The term $V(x,x')\theta(-\sigma)$ has support inside the light cone and is the tail. In this sense, the Hadamard tail is the smooth biscalar $V(x,x')$ appearing in the local Hadamard form.

This definition is local: it requires a normal neighbourhood, so that the relevant geodesic between $v(x_1,x_2)\log\!\big(s^2+i\epsilon\big)$0 and $v(x_1,x_2)\log\!\big(s^2+i\epsilon\big)$1 is unique. Physically, the tail encodes the violation of the strong Huygens principle due to curvature backscatter. Even for a massless field, curvature produces support for the Green function at timelike separation.

2. Governing equations and characteristic data

For a scalar field with mass $v(x_1,x_2)\log\!\big(s^2+i\epsilon\big)$2 and curvature coupling $v(x_1,x_2)\log\!\big(s^2+i\epsilon\big)$3, the retarded Green function satisfies

$v(x_1,x_2)\log\!\big(s^2+i\epsilon\big)$4

The tail biscalar obeys the homogeneous equation

$v(x_1,x_2)\log\!\big(s^2+i\epsilon\big)$5

Its value on the null cone, $v(x_1,x_2)\log\!\big(s^2+i\epsilon\big)$6, is fixed by a transport equation along null geodesics,

$v(x_1,x_2)\log\!\big(s^2+i\epsilon\big)$7

supplemented by the coincidence condition

$v(x_1,x_2)\log\!\big(s^2+i\epsilon\big)$8

A standard local expansion is

$v(x_1,x_2)\log\!\big(s^2+i\epsilon\big)$9

On the light cone, $3+1$0, only $3+1$1 survives, so

$3+1$2

For higher-order characteristic schemes, derivatives of $3+1$3 on the cone involve $3+1$4 as well. This makes the tail a characteristic initial value problem: the null cone is the characteristic hypersurface, the transport equation supplies data on it, and the homogeneous wave equation propagates that data into the cone interior (Aruquipa et al., 2022).

3. Plebański–Hacyan spacetime and null-geodesic crossing

A central modern implementation of this program uses Plebański–Hacyan spacetime,

$3+1$5

with

$3+1$6

This background is a black-hole toy model: null geodesics from a point can cross and form caustics, as in Schwarzschild, but the direct-product structure reduces the scalar wave equation from a generic $3+1$7-dimensional PDE to a $3+1$8-dimensional PDE (Aruquipa et al., 2022).

The world function splits additively,

$3+1$9

with

$$G_{\mathrm{ret}}(x,x')= \big[U(x,x')\delta(\sigma(x,x'))+V(x,x')\theta(-\sigma(x,x'))\big]\theta_+(x,x').$$0

where

$$G_{\mathrm{ret}}(x,x')= \big[U(x,x')\delta(\sigma(x,x'))+V(x,x')\theta(-\sigma(x,x'))\big]\theta_+(x,x').$$1

Hence

$$G_{\mathrm{ret}}(x,x')= \big[U(x,x')\delta(\sigma(x,x'))+V(x,x')\theta(-\sigma(x,x'))\big]\theta_+(x,x').$$2

Before the first caustic, null geodesics satisfy $$G_{\mathrm{ret}}(x,x')= \big[U(x,x')\delta(\sigma(x,x'))+V(x,x')\theta(-\sigma(x,x'))\big]\theta_+(x,x').$$3, so $$G_{\mathrm{ret}}(x,x')= \big[U(x,x')\delta(\sigma(x,x'))+V(x,x')\theta(-\sigma(x,x'))\big]\theta_+(x,x').$$4. In this spacetime, null geodesics focus first at

$$G_{\mathrm{ret}}(x,x')= \big[U(x,x')\delta(\sigma(x,x'))+V(x,x')\theta(-\sigma(x,x'))\big]\theta_+(x,x').$$5

The future boundary of the maximal normal neighbourhood is

$$G_{\mathrm{ret}}(x,x')= \big[U(x,x')\delta(\sigma(x,x'))+V(x,x')\theta(-\sigma(x,x'))\big]\theta_+(x,x').$$6

while the direct null boundary is

$$G_{\mathrm{ret}}(x,x')= \big[U(x,x')\delta(\sigma(x,x'))+V(x,x')\theta(-\sigma(x,x'))\big]\theta_+(x,x').$$7

The tail is defined only between these boundaries.

4. Explicit characteristic formulation and computation in $$G_{\mathrm{ret}}(x,x')= \big[U(x,x')\delta(\sigma(x,x'))+V(x,x')\theta(-\sigma(x,x'))\big]\theta_+(x,x').$$8

The paper introduces

$$G_{\mathrm{ret}}(x,x')= \big[U(x,x')\delta(\sigma(x,x'))+V(x,x')\theta(-\sigma(x,x'))\big]\theta_+(x,x').$$9

In Plebański–Hacyan, $\sigma(x,x')$0, so

$\sigma(x,x')$1

Most of the analysis is massless, hence effectively $\sigma(x,x')$2. The direct coefficient on the light cone is

$\sigma(x,x')$3

and the tail data on the light cone are

$\sigma(x,x')$4

For the specific case $\sigma(x,x')$5,

$\sigma(x,x')$6

Using the symmetry reduction, the PDE for $\sigma(x,x')$7 becomes

$\sigma(x,x')$8

In null coordinates

$\sigma(x,x')$9

the operator takes the form

$x$0

with

$x$1

so the tail equation becomes

$x$2

The characteristic data on the null cone are

$x$3

For a fourth-order scheme one also uses

$x$4

and

$x$5

This implementation computes $x$6 everywhere it is defined, namely throughout the maximal normal neighbourhood of an arbitrary point. For a massless field, the study reports results for

$x$7

For the static-path example with $x$8 and $x$9, the magnitude of $x'$0 decreases as $x'$1 increases; for larger $x'$2, especially $x'$3 in the notation of the plots, the shape of $x'$4 changes more markedly near the caustic region $x'$5. For $x'$6 with $x'$7 on static paths, the characteristic-initial-data evolution agrees with earlier calculations based on infinite sums and integrals (Aruquipa et al., 2022).

5. Physical roles and short-distance quantum-field-theoretic usage

The tail is relevant to classical self-force calculations and to communication between quantum particle detectors. In the self-force context, the retarded Green function near a worldline must be separated into singular/direct and regular/tail pieces. In the detector context, the inside-the-light-cone support permits timelike signaling through quantum fields in curved spacetime (Aruquipa et al., 2022).

A second usage of the same term appears in the short-distance singularity structure of the Feynman two-point function of a scalar field in curved spacetime. In four dimensions, for a Hadamard state,

$x'$8

Here the term

$x'$9

is called the Hadamard tail. For a minimally coupled massless scalar,

$U(x,x')\delta(\sigma)$0

so the local curvature correction is

$U(x,x')\delta(\sigma)$1

With non-minimal coupling,

$U(x,x')\delta(\sigma)$2

For the retarded propagator,

$U(x,x')\delta(\sigma)$3

the Hadamard form implies

$U(x,x')\delta(\sigma)$4

Thus $U(x,x')\delta(\sigma)$5 gives support for $U(x,x')\delta(\sigma)$6, that is, inside the light cone. In this usage the Ricci-scalar-controlled logarithmic term is the short-distance manifestation of the same inside-the-light-cone phenomenon (Chavda et al., 29 Oct 2025).

6. Domain of validity, singular structure, and scope

The local Hadamard form requires a unique geodesic between the points, so the smooth biscalar $U(x,x')\delta(\sigma)$7 is defined only within a normal neighbourhood. In Plebański–Hacyan spacetime, the light-cone value $U(x,x')\delta(\sigma)$8 is regular for $U(x,x')\delta(\sigma)$9 but diverges at $V(x,x')\theta(-\sigma)$0 like

$V(x,x')\theta(-\sigma)$1

Since $V(x,x')\theta(-\sigma)$2 lies outside the maximal normal neighbourhood and corresponds to null-geodesic focusing, this divergence signals the onset of the global singular structure of the retarded Green function.

As one approaches the end of the maximal normal neighbourhood away from caustics, the tail is argued to diverge like

$V(x,x')\theta(-\sigma)$3

This connects the local Hadamard picture with the known global singularity sequence under repeated caustic crossings,

$V(x,x')\theta(-\sigma)$4

Beyond caustics, multiple geodesics connect the points, the standard local Hadamard form breaks down, and one must use other tools such as Kirchhoff representations, mode sums, or global Green-function methods. The characteristic-initial-data implementation on $V(x,x')\theta(-\sigma)$5 therefore serves as a proof of concept specifically inside the maximal normal neighbourhood, but in a spacetime where null geodesics from a point do cross. This suggests a route toward more realistic black-hole spacetimes, while also making clear that the numerical and geometric difficulties increase sharply in higher-dimensional settings such as Schwarzschild and Kerr (Aruquipa et al., 2022).