Gibbons–Hawking–York boundary term explained

In general relativity, the Gibbons–Hawking–York boundary term is a term that needs to be added to the Einstein–Hilbert action when the underlying spacetime manifold has a boundary.

The Einstein–Hilbert action is the basis for the most elementary variational principle from which the field equations of general relativity can be defined. However, the use of the Einstein–Hilbert action is appropriate only when the underlying spacetime manifold

l{M}

is closed, i.e., a manifold which is both compact and without boundary. In the event that the manifold has a boundary

\partiall{M}

, the action should be supplemented by a boundary term so that the variational principle is well-defined.

The necessity of such a boundary term was first realised by James W. York and later refined in a minor way by Gary Gibbons and Stephen Hawking.

For a manifold that is not closed, the appropriate action is

l{S}EH+l{S}GHY=

1
16\pi
4
\int
l{M}d

x\sqrt{-g}R+

1
8\pi

\int\partial

} \mathrm^3 y \, \epsilon \sqrtK,

where

l{S}EH

is the Einstein–Hilbert action,

l{S}GHY

is the Gibbons–Hawking–York boundary term,

hab

is the induced metric (see section below on definitions) on the boundary,

h

its determinant,

K

is the trace of the second fundamental form,

\epsilon

is equal to

+1

where the normal to

\partiall{M}

is spacelike and

-1

where the normal to

\partiall{M}

is timelike, and

ya

are the coordinates on the boundary. Varying the action with respect to the metric

g\alpha\beta

, subject to the condition

\deltag\alpha|\partial

} = 0,

gives the Einstein equations; the addition of the boundary term means that in performing the variation, the geometry of the boundary encoded in the transverse metric

hab

is fixed (see section below). There remains ambiguity in the action up to an arbitrary functional of the induced metric

hab

.

That a boundary term is needed in the gravitational case is because

R

, the gravitational Lagrangian density, contains second derivatives of the metric tensor. This is a non-typical feature of field theories, which are usually formulated in terms of Lagrangians that involve first derivatives of fields to be varied over only.

The GHY term is desirable, as it possesses a number of other key features. When passing to the Hamiltonian formalism, it is necessary to include the GHY term in order to reproduce the correct Arnowitt–Deser–Misner energy (ADM energy). The term is required to ensure the path integral (a la Hawking) for quantum gravity has the correct composition properties. When calculating black hole entropy using the Euclidean semiclassical approach, the entire contribution comes from the GHY term. This term has had more recent applications in loop quantum gravity in calculating transition amplitudes and background-independent scattering amplitudes.

In order to determine a finite value for the action, one may have to subtract off a surface term for flat spacetime:

SEH+SGHY,0=

1
16\pi
4
\int
l{M}d

x\sqrt{-g}R+

1
8\pi

\int\partial

} \mathrm^3 y \, \epsilon \sqrt K - \int_ \mathrm^3 y \, \epsilon \sqrt K_0,

where

K0

is the extrinsic curvature of the boundary imbedded flat spacetime. As

\sqrt{h}

is invariant under variations of

g\alpha

, this addition term does not affect the field equations; as such, this is referred to as the non-dynamical term.

Introduction to hyper-surfaces

Defining hyper-surfaces

In a four-dimensional spacetime manifold, a hypersurface is a three-dimensional submanifold that can be either timelike, spacelike, or null.

A particular hyper-surface

\Sigma

can be selected either by imposing a constraint on the coordinates

f(x\alpha)=0,

or by giving parametric equations,

x\alpha=x\alpha(ya),

where

ya(a=1,2,3)

are coordinates intrinsic to the hyper-surface.

For example, a two-sphere in three-dimensional Euclidean space can be described either by

f(x\alpha)=x2+y2+z2-r2=0,

where

r

is the radius of the sphere, or by

x=r\sin\theta\cos\phi,y=r\sin\theta\sin\phi,z=r\cos\theta,

where

\theta

and

\phi

are intrinsic coordinates.

Hyper-surface orthogonal vector fields

We take the metric convention (-,+,...,+). We start with the family of hyper-surfaces given by

f(x\alpha)=C

where different members of the family correspond to different values of the constant

C

. Consider two neighbouring points

P

and

Q

with coordinates

x\alpha

and

x\alpha+dx\alpha

, respectively, lying in the same hyper-surface. We then have to first order

C=f(x\alpha+dx\alpha)=f(x\alpha)+{\partialf\over\partialx\alpha}dx\alpha.

Subtracting off

C=f(x\alpha)

from this equation gives

{\partialf\over\partialx\alpha}dx\alpha=0

at

P

. This implies that

f,

is normal to the hyper-surface. A unit normal

n\alpha

can be introduced in the case where the hyper-surface is not null. This is defined by

n\alphan\alpha\equiv\epsilon=\begin{cases}-1&if\Sigmaisspacelike\ +1&if\Sigmaistimelike\end{cases}

and we require that

n\alpha

point in the direction of increasing

f:n\alphaf,>0

. It can then easily be checked that

n\alpha

is given by

n\alpha={\epsilon{}f,\over(\epsilon{}g\alphaf,f,)1

}

if the hyper-surface either spacelike or timelike.

Induced and transverse metric

The three vectors

\alpha
e
a

=\left({\partialx\alpha\over\partialya}\right)\partial

} \quad a=1,2,3

are tangential to the hyper-surface.

The induced metric is the three-tensor

hab

defined by

hab=g\alpha

\alpha
e
a
\beta
e
b

.

This acts as a metric tensor on the hyper-surface in the

ya

coordinates. For displacements confined to the hyper-surface (so that

x\alpha=x\alpha(ya)

)

\begin{align} ds2&=g\alphadx\alphadx\beta\\ &=g\alpha\left(

\partialx\alpha
\partialya

dya\right)\left(

\partialx\beta
\partialyb

dyb\right)\\ &=\left(g\alpha

\alpha
e
a
\beta
e
b

\right)dyadyb\\ &=habdyadyb \end{align}

Because the three vectors

\alpha
e
1,
\alpha
e
2,
\alpha
e
3
are tangential to the hyper-surface,

n\alpha

\alpha
e
a

=0

where

n\alpha

is the unit vector (

n\alphan\alpha=\pm1

) normal to the hyper-surface.

We introduce what is called the transverse metric

h\alpha=g\alpha-\epsilonn\alphan\beta.

It isolates the part of the metric that is transverse to the normal

n\alpha

.

It is easily seen that this four-tensor

\alpha}
{h
\beta

=

\alpha}
{\delta
\beta

-\epsilonn\alphan\beta

projects out the part of a four-vector transverse to the normal

n\alpha

as
\alpha}
{h
\beta

n\beta=

\alpha}
({\delta
\beta

-\epsilonn\alphan\beta)n\beta=(n\alpha-\epsilon2n\alpha)=0andifw\alphan\alpha=0then

\alpha}
{h
\beta

w\beta=w\alpha.

We have

hab=h\alpha

\alpha
e
a
\beta
e
b.

If we define

hab

to be the inverse of

hab

, it is easy to check

h\alpha=hab

\alpha
e
a
\beta
e
b

where

h\alpha=g\alpha-\epsilonn\alphan\beta.

Note that variation subject to the condition

\deltag\alpha|\partial

} = 0,

implies that

hab=g\alpha

\alpha
e
a
\beta
e
b
, the induced metric on

\partiall{M}

, is held fixed during the variation. See also [1] for clarification on

\deltah\alpha

and

\deltan\alpha

etc.

On proving the main result

In the following subsections we will first compute the variation of the Einstein–Hilbert term and then the variation of the boundary term, and show that their sum results in

\deltaSTOTAL=\deltaSEH+\deltaSGHY=

1
16\pi

\intl{M}G\alpha\deltag\alpha\sqrt{-g}d4x

where

G\alpha=R\alpha-{1\over2}g\alphaR

is the Einstein tensor, which produces the correct left-hand side to the Einstein field equations, without the cosmological term, which however is trivial to include by replacing

SEH

with

{1\over16\pi}\intl{M}(R-2Λ)\sqrt{-g}d4x

where

Λ

is the cosmological constant.

In the third subsection we elaborate on the meaning of the non-dynamical term.

Variation of the Einstein–Hilbert term

We will use the identity

\delta\sqrt{-g}\equiv-{1\over2}\sqrt{-g}g\alpha\deltag\alpha,

and the Palatini identity:

\deltaR\alpha\equiv\nabla\mu(\delta

\mu
\Gamma
\alpha\beta

)-\nabla\beta(\delta

\mu
\Gamma
\alpha\mu

),

which are both obtained in the article Einstein–Hilbert action.

We consider the variation of the Einstein–Hilbert term:

\begin{align} (16\pi)\deltaSEH&=\intl{M}\delta\left(g\alphaR\alpha\sqrt{-g}\right)d4x\\ &=\intl{M}\left(R\alpha\sqrt{-g}\deltag\alpha+g\alphaR\alpha\delta\sqrt{-g}+\sqrt{-g}g\alpha\deltaR\alpha\right)d4x\\ &=\intl{M}\left(R\alpha-{1\over2}g\alphaR\right)\deltag\alpha\sqrt{-g}d4x+

\alpha\beta
\int
l{M}g

\deltaR\alpha\sqrt{-g}d4x. \end{align}

The first term gives us what we need for the left-hand side of the Einstein field equations. We must account for the second term.

By the Palatini identity

g\alpha\deltaR\alpha=\delta

\mu}
{V
;\mu

,    \deltaV\mu=g\alpha\delta

\mu
\Gamma
\alpha\beta

-g\alpha\delta

\beta
\Gamma
\alpha\beta

.

We will need Stokes theorem in the form:

\begin{align}

\mu}
\int
;\mu

\sqrt{-g}d4x&=\intl{M}(\sqrt{-g}

\mu)
A
,\mu

d4x\\ &=\oint\partial

} A^\mu d \Sigma_\mu \\& = \oint_ \epsilon A^\mu n_\mu \sqrt
d^3y\end

where

n\mu

is the unit normal to

\partiall{M}

and

\epsilon\equivn\mun\mu=\pm1

, and

ya

are coordinates on the boundary. And

d\Sigma\mu=\epsilonn\mud\Sigma

where

d\Sigma=|h|1d3y

where

h=\det[hab]

, is an invariant three-dimensional volume element on the hyper-surface. In our particular case we take

A\mu=\deltaV\mu

.

We now evaluate

\deltaV\mun\mu

on the boundary

\partiall{M}

, keeping in mind that on

\partiall{M},\deltag\alpha=0=\deltag\alpha

. Taking this into account we have

\delta

\mu
\Gamma
\alpha\beta

|\partial

} = \frac g^ (\delta g_ + \delta g_ - \delta g_). It is useful to note that

\begin{align} g\alpha\delta

\beta
\Gamma
\alpha\beta

|\partial

} & = g^ g^ (\delta g_ + \delta g_ - \delta g_) \\& = g^ g^ (\delta g_ + \delta g_ - \delta g_)\end

where in the second line we have swapped around

\alpha

and

\nu

and used that the metric is symmetric. It is then not difficult to work out

\deltaV\mu=g\mug\alpha(\deltag\nu-\deltag\alpha)

.

So now

\begin{align} \deltaV\mun\mu|\partial

} & = n^\mu g^ (\delta g_ - \delta g_) \\& = n^\mu (\epsilon n^\alpha n^\beta + h^) (\delta g_ - \delta g_) \\& = n^\mu h^ (\delta g_ - \delta g_)\end

where in the second line we used the identity

g\alpha=\epsilonn\alphan\beta+h\alpha

, and in the third line we have used the anti-symmetry in

\alpha

and

\mu

. As

\deltag\alpha

vanishes everywhere on the boundary

\partiall{M},

its tangential derivatives must also vanish:

\deltag\alpha

\gamma
e
c

=0

. It follows that

h\alpha\deltag\mu=hab

\alpha
e
a
\beta
e
b

\deltag\mu=0

. So finally we have

n\mu\deltaV\mu|\partial

} = - h^ \delta g_ n^\mu.

Gathering the results we obtain

(16\pi)\deltaSEH=\intl{M}G\alpha\deltag\alpha\sqrt{-g}d4x-\oint\partial

} \epsilon h^ \delta g_ n^\mu \sqrt d^3 y \quad Eq 1.

We next show that the above boundary term will be cancelled by the variation of

SGHY

.

Variation of the boundary term

We now turn to the variation of the

SGHY

term. Because the induced metric is fixed on

\partiall{M},

the only quantity to be varied is

K

is the trace of the extrinsic curvature.

We have

\begin{align} K&=

\alpha}
{n
;\alpha

\\ &=g\alphan\alpha\\ &=\left(\epsilonn\alphan\beta+h\alpha\right)n\alpha\\ &=h\alphan\alpha\\ &=h\alpha(n\alpha,-

\gamma
\Gamma
\alpha\beta

n\gamma) \end{align}

where we have used that

0=(n\alphan\alpha);

implies

n\alphan\alpha;=0.

So the variation of

K

is

\begin{align} \deltaK&=-h\alpha\delta

\gamma
\Gamma
\alpha\beta

n\gamma\\ &=-h\alphan\gamma

1
2

g\gamma\left(\deltag\sigma+\deltag\sigma-\deltag\alpha\right)\\ &=-{1\over2}h\alpha\left(\deltag\mu+\deltag\mu-\deltag\alpha\right)n\mu\\ &=

1
2

h\alpha\deltag\alphan\mu \end{align}

where we have used the fact that the tangential derivatives of

\deltag\alpha

vanish on

\partiall{M}.

We have obtained

(16\pi)\deltaSGHY=\oint\partial

} \epsilon h^ \delta g_ n^\mu \sqrt d^3 y

which cancels the second integral on the right-hand side of Eq. 1. The total variation of the gravitational action is:

\deltaSTOTAL={1\over16\pi}\intl{M}G\alpha\deltag\alpha\sqrt{-g}d4x.

This produces the correct left-hand side of the Einstein equations. This proves the main result.

This result was generalised to fourth-order theories of gravity on manifolds with boundaries in 1983[2] and published in 1985.[3]

The non-dynamical term

We elaborate on the role of

S0={1\over8\pi}\oint\partial

} \epsilon K_0 |h|^ d^3y

in the gravitational action. As already mentioned above, because this term only depends on

hab

, its variation with respect to

g\alpha

gives zero and so does not effect the field equations, its purpose is to change the numerical value of the action. As such we will refer to it as the non-dynamical term.

Let us assume that

g\alpha

is a solution of the vacuum field equations, in which case the Ricci scalar

R

vanishes. The numerical value of the gravitational action is then

S={1\over8\pi}\oint\partial

} \epsilon K |h|^ d^3y,

where we are ignoring the non-dynamical term for the moment. Let us evaluate this for flat spacetime. Choose the boundary

\partiall{M}

to consist of two hyper-surfaces of constant time value

t=t1,t2

and a large three-cylinder at

r=r0

(that is, the product of a finite interval and a three-sphere of radius

r0

). We have

K=0

on the hyper-surfaces of constant time. On the three cylinder, in coordinates intrinsic to the hyper-surface, the line element is

\begin{align} ds2&=-dt2+

2
r
0

d\Omega2\\ &=-dt2+

2
r
0

(d\theta2+\sin2\thetad\phi2) \end{align}

meaning the induced metric is

hab=\begin{bmatrix}-1&0&0\ 0&

2
r
0

&0\ 0&0&

2
r
0

\sin2\theta\end{bmatrix}.

so that

|h|1=

2
r
0

\sin\theta

. The unit normal is

n\alpha=\partial\alphar

, so

K=

\alpha}
{n
;\alpha

=2/r0

. Then

\oint\partial

} \epsilon K |h|^ d^3y = \int_^ dt \int_0^ d \varphi \int_0^\pi d \theta \left(\right) (r_0^2 \sin \theta) = 8 \pi r_0 (t_2 - t_1)

and diverges as

r0\toinfty

, that is, when the spatial boundary is pushed to infinity, even when the

l{M}

is bounded by two hyper-surfaces of constant time. One would expect the same problem for curved spacetimes that are asymptotically flat (there is no problem if the spacetime is compact). This problem is remedied by the non-dynamical term. The difference

SGHY-S0

will be well defined in the limit

r0\toinfty

.

Variation of modified gravity terms

See main article: Alternatives to general relativity. There are many theories which attempt to modify General Relativity in different ways, for example f(R) gravity replaces R, the Ricci scalar in the Einstein–Hilbert action with a function f(R). Guarnizo et al. found the boundary term for a general f(R) theory.[4] They found that the "modified action in the metric formalism of f(R) gravity plus a Gibbons–York–Hawking like boundary term must be written as:"

Smod=

1
2\kappa

\intVd4x\sqrt{-g}f(R)+2\int\partiald3y\epsilon|h|f'(R)K

where

f'(R)\equiv

df(R)
dR
.

By using the ADM decomposition and introducing extra auxiliary fields, in 2009 Deruelle et al. found a method to find the boundary term for "gravity theories whose Lagrangian is an arbitrary function of the Riemann tensor."[5] This method can be used to find the GHY boundary terms for Infinite derivative gravity.[6]

A path-integral approach to quantum gravity

As mentioned at the beginning, the GHY term is required to ensure the path integral (a la Hawking et al.) for quantum gravity has the correct composition properties.

This older approach to path-integral quantum gravity had a number of difficulties and unsolved problems. The starting point in this approach is Feynman's idea that one can represent the amplitude

\langleg2,\phi2,\Sigma2|g1,\phi1,\Sigma1\rangle

to go from the state with metric

g1

and matter fields

\phi1

on a surface

\Sigma1

to a state with metric

g2

and matter fields

\phi2

on a surface

\Sigma2

, as a sum over all field configurations

g

and

\phi

which take the boundary values of the fields on the surfaces

\Sigma1

and

\Sigma2

. We write

\langleg2,\phi2,\Sigma2|g1,\phi1,\Sigma1\rangle=\intl{D}[g,\phi]\exp(iS[g,\phi])

where

l{D}[g,\phi]

is a measure on the space of all field configurations

g

and

\phi

,

S[g,\phi]

is the action of the fields, and the integral is taken over all fields which have the given values on

\Sigma1

and

\Sigma2

.

It is argued that one need only specify the three-dimensional induced metric

h

on the boundary.

Now consider the situation where one makes the transition from metric

h1

, on a surface

\Sigma1

, to a metric

h2

, on a surface

\Sigma2

and then on to a metric

h3

on a later surface

\Sigma3

One would like to have the usual composition rule

\langleh3,\Sigma3|h1,\Sigma1\rangle=

\sum
h2

\langleh3,\Sigma3|h2,\Sigma2\rangle\langleh2,\Sigma2|h1,\Sigma1\rangle

expressing that the amplitude to go from the initial to final state to be obtained by summing over all states on the intermediate surface

\Sigma2

.

Let

g1

be the metric between

\Sigma1

and

\Sigma2

and

g2

be the metric between

\Sigma2

and

\Sigma3

. Although the induced metric of

g1

and

g2

will agree on

\Sigma2

, the normal derivative of

g1

at

\Sigma2

will not in general be equal to that of

g2

at

\Sigma2

. Taking the implications of this into account, it can then be shown that the composition rule will hold if and only if we include the GHY boundary term.[7]

In the next section it is demonstrated how this path integral approach to quantum gravity leads to the concept of black hole temperature and intrinsic quantum mechanical entropy.

Calculating black-hole entropy using the Euclidean semi-classical approach

See main article: Euclidean quantum gravity.

Application in loop quantum gravity

See main article: Loop quantum gravity.

Transition amplitudes and the Hamilton's principal function

In the quantum theory, the object that corresponds to the Hamilton's principal function is the transition amplitude. Consider gravity defined on a compact region of spacetime, with the topology of a four dimensional ball. The boundary of this region is a three-dimensional space with the topology of a three-sphere, which we call

\Sigma

. In pure gravity without cosmological constant, since the Ricci scalar vanishes on solutions of Einstein's equations, the bulk action vanishes and the Hamilton's principal function is given entirely in terms of the boundary term,

S[q]=\int\SigmaKab[q]qab\sqrt{q}d3\sigma

where

Kab

is the extrinsic curvature of the boundary,

qab

is the three-metric induced on the boundary, and

\sigma

are coordinates on the boundary.

The functional

S[q]

is a highly non-trivial functional to compute; this is because the extrinsic curvature

Kab[q]

is determined by the bulk solution singled out by the boundary intrinsic geometry. As such

Kab[q]

is non-local. Knowing the general dependence of

Kab

from

qab

is equivalent to knowing the general solution of the Einstein equations.

Background-independent scattering amplitudes

Loop quantum gravity is formulated in a background-independent language. No spacetime is assumed a priori, but rather it is built up by the states of theory themselves however scattering amplitudes are derived from

n

-point functions (Correlation function (quantum field theory)) and these, formulated in conventional quantum field theory, are functions of points of a background space-time. The relation between the background-independent formalism and the conventional formalism of quantum field theory on a given spacetime is far from obvious, and it is far from obvious how to recover low-energy quantities from the full background-independent theory. One would like to derive the

n

-point functions of the theory from the background-independent formalism, in order to compare them with the standard perturbative expansion of quantum general relativity and therefore check that loop quantum gravity yields the correct low-energy limit.

A strategy for addressing this problem has been suggested;[8] the idea is to study the boundary amplitude, or transition amplitude of a compact region of spacetime, namely a path integral over a finite space-time region, seen as a function of the boundary value of the field.[9] [10] In conventional quantum field theory, this boundary amplitude is well-defined[11] [12] and codes the physical information of the theory; it does so in quantum gravity as well, but in a fully background-independent manner.[13] A generally covariant definition of

n

-point functions can then be based on the idea that the distance between physical points arguments of the

n

-point function is determined by the state of the gravitational field on the boundary of the spacetime region considered.

The key observation is that in gravity the boundary data include the gravitational field, hence the geometry of the boundary, hence all relevant relative distances and time separations. In other words, the boundary formulation realizes very elegantly in the quantum context the complete identification between spacetime geometry and dynamical fields.

See also

References

External links

Notes and References

  1. Feng, J. C., Matzner R. A. The Weiss variation of the gravitational action. Theory Group, Department of Physics, University of Texas at Austin. arXiv:1708.04489v3 [gr-qc]. 24 July 2018 https://arxiv.org/pdf/1708.04489
  2. Web site: Second and fourth order gravitational actions on manifolds with boundaries. ResearchGate . en. 2017-05-08.
  3. Barth . N H . The fourth-order gravitational action for manifolds with boundaries . Classical and Quantum Gravity . IOP Publishing . 2 . 4 . 1985-07-01 . 0264-9381 . 10.1088/0264-9381/2/4/015 . 497–513. 1985CQGra...2..497B . 250893849 .
  4. 1002.0617. Boundary Term in Metric f(R) Gravity: Field Equations in the Metric Formalism. General Relativity and Gravitation . 42. 11. 2713–2728. Guarnizo . Alejandro . Castaneda. Leonardo. Tejeiro. Juan M.. 2010. 10.1007/s10714-010-1012-6. 2010GReGr..42.2713G. 119099298.
  5. 0908.0679. Hamiltonian formulation of f(Riemann) theories of gravity. Progress of Theoretical Physics. 123. 169–185. Deruelle. Nathalie. Nathalie Deruelle . Sasaki. Misao . Sendouda. Yuuiti. Yamauchi. Daisuke. 2010. 1 . 10.1143/PTP.123.169. 2010PThPh.123..169D. 118570242.
  6. 1606.01911 . Teimouri . Ali . Generalised Boundary Terms for Higher Derivative Theories of Gravity. Journal of High Energy Physics. 2016 . 8 . Talaganis. Spyridon. Edholm . James. Mazumdar . Anupam. 2016. 144 . 10.1007/JHEP08(2016)144. 2016JHEP...08..144T. 55220918 .
  7. For example see the book "Hawking on the big bang and black holes" by Stephen Hawking, chapter 15.
  8. Modesto . Leonardo . Rovelli . Carlo . Particle Scattering in Loop Quantum Gravity . Physical Review Letters . 95 . 19 . 2005-11-01 . 0031-9007 . 10.1103/physrevlett.95.191301 . 191301. 16383970 . gr-qc/0502036 . 2005PhRvL..95s1301M . 46705469 .
  9. Oeckl . Robert . A "general boundary" formulation for quantum mechanics and quantum gravity . Physics Letters B . Elsevier BV . 575 . 3–4 . 2003 . 0370-2693 . 10.1016/j.physletb.2003.08.043. hep-th/0306025 . 318–324. 2003PhLB..575..318O . free.
  10. Oeckl . Robert . Schrödinger's cat and the clock: lessons for quantum gravity . Classical and Quantum Gravity . 20 . 24 . 2003-11-03 . 0264-9381 . 10.1088/0264-9381/20/24/009 . 5371–5380. gr-qc/0306007 . 2003CQGra..20.5371O . 118978523 .
  11. Conrady . Florian . Rovelli . Carlo . Generalized Schrödinger equation in Euclidean field theory . International Journal of Modern Physics A . World Scientific Pub Co Pte Lt . 19 . 24 . 2004-09-30 . 0217-751X . 10.1142/s0217751x04019445 . 4037–4068. hep-th/0310246. 2004IJMPA..19.4037C . 18048123 .
  12. Doplicher . Luisa . Generalized Tomonaga-Schwinger equation from the Hadamard formula . Physical Review D . American Physical Society (APS) . 70 . 6 . 2004-09-24 . 1550-7998 . 10.1103/physrevd.70.064037 . 064037. gr-qc/0405006. 2004PhRvD..70f4037D . 14402915 .
  13. Conrady . Florian . Doplicher . Luisa . Oeckl . Robert . Rovelli . Carlo . Testa . Massimo . Minkowski vacuum in background independent quantum gravity . Physical Review D . American Physical Society (APS) . 69 . 6 . 2004-03-18 . 1550-7998 . 10.1103/physrevd.69.064019 . 064019. gr-qc/0307118. 2004PhRvD..69f4019C . 30190407 .