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}
\partiall{M}
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
where
l{S}EH
l{S}GHY
hab
h
K
\epsilon
+1
\partiall{M}
-1
\partiall{M}
ya
g\alpha\beta
\deltag\alpha|\partial
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
hab
That a boundary term is needed in the gravitational case is because
R
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
where
K0
\sqrt{h}
g\alpha
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
f(x\alpha)=0,
or by giving parametric equations,
x\alpha=x\alpha(ya),
where
ya(a=1,2,3)
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
x=r\sin\theta\cos\phi, y=r\sin\theta\sin\phi, z=r\cos\theta,
where
\theta
\phi
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
P
Q
x\alpha
x\alpha+dx\alpha
C=f(x\alpha+dx\alpha)=f(x\alpha)+{\partialf\over\partialx\alpha}dx\alpha.
Subtracting off
C=f(x\alpha)
{\partialf\over\partialx\alpha}dx\alpha=0
at
P
f,
n\alpha
n\alphan\alpha\equiv\epsilon=\begin{cases}-1&if\Sigmaisspacelike\ +1&if\Sigmaistimelike\end{cases}
and we require that
n\alpha
f:n\alphaf,>0
n\alpha
n\alpha={\epsilon{}f,\over(\epsilon{}g\alphaf,f,)1
if the hyper-surface either spacelike or timelike.
The three vectors
\alpha | |
e | |
a |
=\left({\partialx\alpha\over\partialya}\right)\partial
are tangential to the hyper-surface.
The induced metric is the three-tensor
hab
hab=g\alpha
\alpha | |
e | |
a |
\beta | |
e | |
b |
.
This acts as a metric tensor on the hyper-surface in the
ya
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 |
n\alpha
\alpha | |
e | |
a |
=0
where
n\alpha
n\alphan\alpha=\pm1
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
\alpha} | |
{h | |
\beta |
n\beta=
\alpha} | |
({\delta | |
\beta |
-\epsilonn\alphan\beta)n\beta=(n\alpha-\epsilon2n\alpha)=0 and if w\alphan\alpha=0 then
\alpha} | |
{h | |
\beta |
w\beta=w\alpha.
We have
hab=h\alpha
\alpha | |
e | |
a |
\beta | |
e | |
b. |
If we define
hab
hab
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
implies that
hab=g\alpha
\alpha | |
e | |
a |
\beta | |
e | |
b |
\partiall{M}
\deltah\alpha
\deltan\alpha
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
SEH
{1\over16\pi}\intl{M}(R-2Λ)\sqrt{-g}d4x
where
Λ
In the third subsection we elaborate on the meaning of the non-dynamical 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
where
n\mu
\partiall{M}
\epsilon\equivn\mun\mu=\pm1
ya
d\Sigma\mu=\epsilonn\mud\Sigma
d\Sigma=|h|1d3y
h=\det[hab]
A\mu=\deltaV\mu
We now evaluate
\deltaV\mun\mu
\partiall{M}
\partiall{M},\deltag\alpha=0=\deltag\alpha
\delta
\mu | |
\Gamma | |
\alpha\beta |
|\partial
\begin{align} g\alpha\delta
\beta | |
\Gamma | |
\alpha\beta |
|\partial
where in the second line we have swapped around
\alpha
\nu
\deltaV\mu=g\mug\alpha(\deltag\nu-\deltag\alpha)
So now
\begin{align} \deltaV\mun\mu|\partial
where in the second line we used the identity
g\alpha=\epsilonn\alphan\beta+h\alpha
\alpha
\mu
\deltag\alpha
\partiall{M},
\deltag\alpha
\gamma | |
e | |
c |
=0
h\alpha\deltag\mu=hab
\alpha | |
e | |
a |
\beta | |
e | |
b |
\deltag\mu=0
n\mu\deltaV\mu|\partial
Gathering the results we obtain
(16\pi)\deltaSEH=\intl{M}G\alpha\deltag\alpha\sqrt{-g}d4x-\oint\partial
We next show that the above boundary term will be cancelled by the variation of
SGHY
We now turn to the variation of the
SGHY
\partiall{M},
K
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);
n\alphan\alpha;=0.
K
\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
\partiall{M}.
(16\pi)\deltaSGHY=\oint\partial
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]
We elaborate on the role of
S0={1\over8\pi}\oint\partial
in the gravitational action. As already mentioned above, because this term only depends on
hab
g\alpha
Let us assume that
g\alpha
R
S={1\over8\pi}\oint\partial
where we are ignoring the non-dynamical term for the moment. Let us evaluate this for flat spacetime. Choose the boundary
\partiall{M}
t=t1,t2
r=r0
r0
K=0
\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
n\alpha=\partial\alphar
K=
\alpha} | |
{n | |
;\alpha |
=2/r0
\oint\partial
and diverges as
r0\toinfty
l{M}
SGHY-S0
r0\toinfty
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]
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
\phi1
\Sigma1
g2
\phi2
\Sigma2
g
\phi
\Sigma1
\Sigma2
\langleg2,\phi2,\Sigma2|g1,\phi1,\Sigma1\rangle=\intl{D}[g,\phi]\exp(iS[g,\phi])
where
l{D}[g,\phi]
g
\phi
S[g,\phi]
\Sigma1
\Sigma2
It is argued that one need only specify the three-dimensional induced metric
h
Now consider the situation where one makes the transition from metric
h1
\Sigma1
h2
\Sigma2
h3
\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
\Sigma1
\Sigma2
g2
\Sigma2
\Sigma3
g1
g2
\Sigma2
g1
\Sigma2
g2
\Sigma2
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.
See main article: Euclidean quantum gravity.
See main article: Loop quantum gravity.
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
S[q]=\int\SigmaKab[q]qab\sqrt{q} d3\sigma
where
Kab
qab
\sigma
The functional
S[q]
Kab[q]
Kab[q]
Kab
qab
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
n
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
n
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.