This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

KUNS-2995
Coarse-graining black holes out of equilibrium
with boundary observables on time slice

Daichi Takeda Department of Physics, Kyoto University, Kyoto 606-8502, Japan
Abstract

In black hole thermodynamics, defining coarse-grained entropy for dynamical black holes has long been a challenge, and various proposals, such as generalized entropy, have been explored. Guided by the AdS/CFT, we introduce a new definition of coarse-grained entropy for a dynamical black hole in Lorentzian Einstein gravity. On each time slice, this entropy is defined as the horizon area of an auxiliary Euclidean black hole that shares the same mass, (angular) momenta, and asymptotic normalizable matter modes with the original Lorentzian solution. The entropy is shown to satisfy a generalized first law and, through holography, the second law as well. Furthermore, by applying this thermodynamics to several Vaidya models in AdS and flat spacetime, we discover a connection between the second law and the null energy condition.

000 [email protected]

1 Introduction and summary

Gravity is considered thermodynamic [1]. The zeroth law asserts the existence of intensive variables, which in black hole thermodynamics correspond to quantities constant over the event horizon, such as surface gravity and angular velocity. The first law is one of the key concept, stating the energy conservation among neighboring stationary black holes [2, 3] when the horizon area is viewed as the entropy [4]. Applied to the local Rindler patch, the first law reproduces the Einstein equation, allowing us to regard it as an equation of state [5]. Above all, not only does black hole thermodynamics apply as a formal analogy, but it also has physical reality in the Hawking radiation [6].

Thermodynamics helps us understand spacetime physics, because it gives a macroscopic constraints that the statistical mechanics of a quantum theory of gravity must adhere to. First, the microscopic counting of states in string theory (beginning with [7]) agrees with the horizon area, confirming the importance of exploring fundamental theories through macroscopic observations. Second, given the entropy as a function of sufficient extensive variables, thermodynamics dictates among what states transitions are allowed to happen. This is the role of the second law. Since general relativity itself does not teach what “physical” time evolutions are, energy conditions are necessary for describing “physical” systems. Here, the thermodynamics of spacetime is expected to provide an answer if established.

However, the second law remains under debate. While the Hawking area theorem [8] guarantees that the horizon area does not decrease toward the future, the question as to whether the horizon area can be regarded as the entropy even out of equilibrium is nontrivial. Besides, it is incompatible with the Hawking radiation occurring in the semiclassical regime. An alternative proposal that has been extensively studied is the generalized entropy [9, 10], which is the sum of the entropy of matter outside the horizon on a time slice and the area of the cross section between the horizon and the time slice. Despite numerous attempts, however, a proof of the generalized second law valid in any situation, or an agreement on a suitable definition of the entropy for dynamical black holes has not been achieved (see [11, 12, 13, 14] for review). The second law is one of the primary topics addressed in this paper.111 The third law seems not to hold in general for black hole physics [15, 16], but it is not a drawback since the third law is just a phenomenological observation in laboratories, and in principle, it is not necessary in the axiomatic construction of thermodynamics.

Gravity, on the other hand, is considered holographic. In the classical Einstein gravity, the bulk contribution of the on-shell action vanishes and meaningful values arise from the boundary term, by which we define the energy, (angular) momenta, and other charges. At the same time, the entropy in stationary cases is defined on the causal boundary, the horizon, and interestingly, the first law holds among those boundary values [3].

The holographic nature is a key concept also at the quantum level [17, 18, 19, 20], having provided hints towards quantum gravity. Among those, the AdS/CFT correspondence [21, 22, 23] has attracted significant attention so far. According to the correspondence, the gravitational degrees of freedom are mapped onto the quantum field theory defined on the AdS boundary. Notably, the entanglement entropy is expected to reveal the nature of the dual spacetime [24, 25, 26, 27]. The entanglement entropy of a boundary region is equivalent to the area of an extremal surface in the bulk, and in the semiclassical regime the von Neumann entropy of matter is added to it as quantum corrections [28, 29, 30]; the total entropy takes a form similar to the thermodynamic generalized entropy explained above. This direction recently developed into the Island formula [31, 32, 33, 34, 35] to study the Page curve [36, 37], with applications extending beyond AdS gravity [38].

When applying holography to black hole thermodynamics, particularly concerning the second law, it seems better to introduce coarse-grained entropies. This is because, for instance, the unitarity of the CFT ensures that the fine-grained entropy of a total system state is time-independent, which makes the HRT surface [26] insensitive to any probe sent from the boundary [39]. Within the framework of the AdS/CFT, a coarse-grained entropy associated with the second law in the bulk was first explored in [40], where the one-point entropy was introduced. The relationship with the first law was also investigated in [41]. In [40], the bulk dual of the one-point entropy was first conjectured to be the causal holographic information [42, 43], but unfortunately, counterexamples were identified in [44], leading the authors to conclude that the bulk counterpart seems more intricate than initially proposed. Other approaches, as seen in [45, 46, 47, 48], begin with defining their coarse-graining procedures in the bulk. While [45, 47] and [48] employ different ways of coarse-graining, both approaches are associated with the apparent horizon, presenting intriguing implications to gravity and its thermodynamics. However, as noted by the authors, none of those approaches has yet acquired a complete description within the boundary language.

In the same spirit, we propose via the AdS/CFT a new coarse-grained entropy, which is defined both on the boundary and in the bulk. It is valuable to study variety of ways of coarse-graining, as the coarse-grained entropy differs depending on what information one aims to preserve, and different coarse-graining methods will offer different perspectives on gravity. The entropy we propose satisfies a generalized version of the first law, and by the AdS/CFT the second law as well. In addition, it coincides with the Bekenstein-Hawking entropy for stationary solutions.

We begin with the positivity of the relative entropy in the boundary quantum theory.222 The usage of relative entropy in the AdS/CFT to survey the bulk is also seen in [49, 50, 51, 52]. Let ρ1\rho_{1} and ρ2\rho_{2} be normalized density operators, then the relative entropy is defined as

S(ρ1||ρ2):=Tr(ρ1lnρ1ρ1lnρ2),\displaystyle S(\rho_{1}||\rho_{2}):=\mathrm{Tr}(\rho_{1}\ln\rho_{1}-\rho_{1}\ln\rho_{2}), (1.1)

which is always non-negative. Suppose the initial state is a steady state ρ0\rho_{0}, which evolves to ρt\rho_{t} with a time-dependent Hamiltonian. The time-dependence is triggered by sources coupling to composite operators. On the other hand, we prepare a reference state ρref,t\rho_{\mathrm{ref},t} sharing with ρt\rho_{t} the same expectation values of composite operators {OI}\{O_{I}\} of our interest: Tr(ρref,tOI)=Tr(ρtOI)\mathrm{Tr}(\rho_{\mathrm{ref},t}O_{I})=\mathrm{Tr}(\rho_{t}O_{I}). At t=0t=0 we set the reference state as ρref,0=ρ0\rho_{\mathrm{ref},0}=\rho_{0}. Then, we define our coarse-grained entropy by St:=Tr(ρref,tlnρref,t)S_{t}:=-\mathrm{Tr}(\rho_{\mathrm{ref},t}\ln\rho_{\mathrm{ref},t}), and will show that StS0S_{t}\geq S_{0} holds for any tt, based on the positivity S(ρt||ρref,t)0S(\rho_{t}||\rho_{\mathrm{ref},t})\geq 0.

This coarse-grained entropy is in the bulk equivalent to the horizon area of an auxiliary Euclidean spacetime333 As is also the case in [47], an appropriate auxiliary spacetime different from the original is necessary, since the coarse-graining process changes the state on the boundary. realized as the dominant saddle point of the gravitational path integral dual to ρref,t\rho_{\mathrm{ref},t}. In the bulk language, Tr(ρref,tOI)=Tr(ρtOI)\mathrm{Tr}(\rho_{\mathrm{ref},t}O_{I})=\mathrm{Tr}(\rho_{t}O_{I}) means that the Euclidean solution shares the same asymptotic normalizable modes with the original Lorentzian time slice. Besides, StS0S_{t}\geq S_{0} means that the horizon area of the Euclidean black hole we refer to at tt is always greater than or equal to the horizon area of the initial black hole. This is the second law (for adiabatic process) in this paper. As explained below, this is not an assertion of monotonically increasing area as in Hawking’s area law. After obtaining the gravitational description, we will also derive the first law, which, in addition to the usual terms, contains contributions from the asymptotic matter modes.

We will also check the two laws explicitly in several null-ray collapse models and find that our second law implies the null energy condition — derived by the AdS/CFT, this can be viewed as a consequence of quantum gravity. The gravitational description we finally obtain can formally be applied to non-AdS spacetimes. To test its applicability, an asymptotically flat collapse is chosen as one of the examples, and the two laws actually hold in this example.

It is worth mentioning that the second law in this paper, though not monotonic, is compatible with the expressions of the second law in the context of non-equilibrium thermodynamics. The second law in thermodynamics serves as the criterion for possible transitions between steady states. In non-equilibrium thermodynamics, the second law is generalized to the statement that the entropy production444 The Shannon entropy change plus the heat absorption divided by temperature. never becomes negative. When the total system is deterministic, the second law follows from the result of [53, 54], where the initial state is prepared as a product state of the system and bath, with the bath in equilibrium. Since the initial state is prepared specially, the monotonicity is not guaranteed. The second law remains monotonic, for example, when the system is Markovian [55, 56], as a Markovian system does not care about its history, making the initial state not special anymore. It may not be obvious whether the monotonicity holds for dynamical gravitational systems as well. The quantum versions of the second law without the monotonicity were also shown via the relative entropy in [57, 58] (see [59] for review). The organization of this paper is as follows. In section 2, the second law is derived in the QFT. In section 3, we rewrite it in the bulk Einstein gravity, and also derive the (generalized) first law. In section 4, the derived thermodynamic laws are applied to Vaidya-type models, indicating that thermodynamics constraints gravity.

Summary of the results

Here, the main results in the gravity side is roughly overviewed. Let MM be a (d+1)(d+1)-dimensional Lorentzian manifold with timelike boundary M\partial M. The theory is supposed to be the Einstein theory with matter fields, whose action is denoted as IgravI_{\mathrm{grav}} (see (3.1)). Gauge fields like Maxwell field can be included.

We suppose that a stationary configuration is realized on the initial time slice. This is prepared by analytically continuing a Euclidean solution (figure 3). To define a coarse-grained entropy at time t>0t>0, we need to specify a set of fields to be respected, whose values on M\partial M we write as {wtI(θ)}\{w^{I}_{t}(\theta)\}. Here, II is the label of fields, and the boundary coordinate is written as z=(t,θA)z=(t,\theta^{A}) with A=1,,d1A=1,\cdots,d-1. The induced metric on M\partial M is assumed static, while each wtI(θ)w^{I}_{t}(\theta) can depend on time. The normalizable mode conjugate to wtI(θ)w^{I}_{t}(\theta) is defined to be

πI,t(θ):=δδwtI(θ)(Igrav|onshell[w]).\displaystyle\pi_{I,t}(\theta):=\frac{\delta}{\delta w^{I}_{t}(\theta)}\left(\left.I_{\mathrm{grav}}\right|_{\mathrm{on-shell}}[w]\right). (1.2)

We write the mass and the (angular) momenta in θA\theta^{A}-direction as hth_{t} and pA,tp_{A,t}, both of which are defined by the Brown-York tensor (3.8). Note that πI,t(θ)\pi_{I,t}(\theta) is a local quantity, dependent on the spatial coordinates θA\theta^{A}, while hth_{t} and pA,tp_{A,t} are independent of θA\theta^{A}.

Our coarse-grained entropy StS_{t} is defined by the following procedure. On each time slice, we find the auxiliary Euclidean solution which dominates the gravitational path integral, while respecting hth_{t}, pA,tp_{A,t} and πI,t\pi_{I,t} as its mass, (angular) momenta, and asymptotic modes. We assume that this Euclidean solution is realized stationary, independent of the imaginary time. Then, StS_{t} is defined to be the horizon area, the area of the cigar tip in figure 4. Since the initial configuration of MM is stationary, S0S_{0} is exactly the Bekenstein-Hawking entropy of the initial black hole.

The main claims of this paper are as follows. First, S˙t\dot{S}_{t} (the dot means the time-derivative) is given by the first law, with additional contributions consisting of {wI}\{w^{I}\} and {πI}\{\pi_{I}\}. Second, StS_{t} satisfies the second law, StS0S_{t}\geq S_{0}. While the first law is derived within the Einstein gravity, the second law is derived via the holographic dictionary. Nevertheless, we conjecture that it must hold also in generic cases under the setup mentioned in the beginning (see section 3.5.1 and 4.3).

2 Coarse-grained state in QFT

In this section, we first introduce the reference state ρref,t\rho_{\mathrm{ref},t}, a coarse-grained state to be compared with the original state at time tt through the relative entropy. The second law is derived from the positivity of the relative entropy. After that, in preparation for section 3, we move on to the path integral representation.

2.1 Coarse-grained state and second law

2.1.1 The simplest case

Let us grasp the essence with a warm-up example. We consider a generic quantum theory. The initial state at t=0t=0 is supposed to be

ρ0=eβH(0)Z0(β),Z0(β):=TreβH(0),\displaystyle\rho_{0}=\frac{e^{-\beta H(0)}}{Z_{0}(\beta)},\qquad Z_{0}(\beta):=\mathrm{Tr}e^{-\beta H(0)}, (2.1)

where H(0)H(0) is the Hamiltonian at t=0t=0 and β0\beta_{0} is an arbitrary inverse temperature. The system evolves with a time-dependent Hamiltonian H(t)H(t) as

ρt=U(t)ρ0U(t)1,U(t):=Texp(i0tdsH(s)).\displaystyle\rho_{t}=U(t)\rho_{0}U(t)^{-1},\qquad U(t):=\mathrm{T}\exp\left(-i\int_{0}^{t}\mathrm{d}s\,H(s)\right). (2.2)

We choose the following as the reference state:

ρref,t=eβH(t)Zt(β),Zt(β):=TreβH(t),\displaystyle\rho_{\mathrm{ref},t}=\frac{e^{-\beta^{\prime}H(t)}}{Z_{t}(\beta^{\prime})},\qquad Z_{t}(\beta^{\prime}):=\mathrm{Tr}e^{-\beta^{\prime}H(t)}, (2.3)

with β\beta^{\prime} being any positive parameter.

In this setup, we compute the relative entropy S(ρt||ρref,t)S(\rho_{t}||\rho_{\mathrm{ref},t}) defined in (1.1). First, since the evolution is unitary, we have

Tr(ρtlnρt)=Tr(ρ0lnρ0)=lnZ0(β)βTr(ρ0H(0)).\displaystyle\mathrm{Tr}(\rho_{t}\ln\rho_{t})=\mathrm{Tr}(\rho_{0}\ln\rho_{0})=-\ln Z_{0}(\beta)-\beta\mathrm{Tr}\left(\rho_{0}H(0)\right). (2.4)

The remaining piece in (1.1) is readily

Tr(ρtlnρref,t)=lnZt(β)βTr(ρtH(t)).\displaystyle\mathrm{Tr}(\rho_{t}\ln\rho_{\mathrm{ref},t})=-\ln Z_{t}(\beta^{\prime})-\beta^{\prime}\mathrm{Tr}\left(\rho_{t}H(t)\right). (2.5)

The inequality S(ρt||ρref,t)0S(\rho_{t}||\rho_{\mathrm{ref},t})\geq 0 reads

Tr(ρtlnρref,t)=βTr(ρtH(t))+lnZt(β)βTr(ρ0H(0))+lnZ0(β).\displaystyle-\mathrm{Tr}(\rho_{t}\ln\rho_{\mathrm{ref},t})=\beta^{\prime}\mathrm{Tr}\left(\rho_{t}H(t)\right)+\ln Z_{t}(\beta^{\prime})\geq\beta\mathrm{Tr}\left(\rho_{0}H(0)\right)+\ln Z_{0}(\beta). (2.6)

When β\beta^{\prime} is set to be β\beta, the inequality gives [57]

Tr(ρtH(t))Tr(ρ0H(0))[β1lnZt(β)][β1lnZ0(β)].\displaystyle\mathrm{Tr}\left(\rho_{t}H(t)\right)-\mathrm{Tr}\left(\rho_{0}H(0)\right)\geq[-\beta^{-1}\ln Z_{t}(\beta)]-[-\beta^{-1}\ln Z_{0}(\beta)]. (2.7)

The l.h.s. is regarded as the work by the time-dependent part of the Hamiltonian, not as the heat because the system is unitary. The r.h.s. can be seen as the free energy difference. This coincides with the maximum work principle.

In this paper, we would rather like to ask what β\beta^{\prime} makes (2.6) the tightest bound. The condition that extremizes the l.h.s. of (2.6) is found to be

Tr(ρtH(t))=Tr(ρref,tH(t)),\displaystyle\mathrm{Tr}\left(\rho_{t}H(t)\right)=\mathrm{Tr}\left(\rho_{\mathrm{ref},t}H(t)\right), (2.8)

where β\beta^{\prime}-dependence resides in ρref,t\rho_{\mathrm{ref},t}. With β\beta^{\prime} satisfying this, (2.6) is reduced to

S(ρref,t)S(ρ0),S(ρ):=Tr(ρlnρ).\displaystyle S(\rho_{\mathrm{ref},t})\geq S(\rho_{0}),\qquad S(\rho):=-\mathrm{Tr}\left(\rho\ln\rho\right). (2.9)

Here are some comments. First, the extremal condition (2.8) actually gives the minimum, since the l.h.s. of (2.6) goes to positive infinity as β\beta^{\prime}\to\infty. Second, our ansatz (2.3) is justified by considering coarse-graining ρt\rho_{t} as follows. When coarse-graining a state ρt\rho_{t} with a given operator set {OI}\{O_{I}\} respected, we find the maximum of S(ρ)S(\rho^{\prime}) under the following condition:

Tr(ρOI)=Tr(ρtOI)for all OI in {OI}.\displaystyle\mathrm{Tr}\left(\rho^{\prime}O_{I}\right)=\mathrm{Tr}\left(\rho_{t}O_{I}\right)\qquad\mbox{for all $O_{I}$ in $\{O_{I}\}$}. (2.10)

This problem is solved by the optimizing

S~(ρ)=S(ρ)λI[Tr(ρOI)Tr(ρOI)]+λ~(Trρ1),\displaystyle\tilde{S}(\rho^{\prime})=S(\rho^{\prime})-\lambda^{I}\left[\mathrm{Tr}\left(\rho^{\prime}O_{I}\right)-\mathrm{Tr}\left(\rho O_{I}\right)\right]+\tilde{\lambda}(\mathrm{Tr}\rho-1), (2.11)

where λI\lambda^{I} and λ~\tilde{\lambda} are the Lagrange multipliers, and summation is taken for the label II. Then the solution is found to be ρ=ρcg,texp(λIOI)\rho^{\prime}=\rho_{\mathrm{cg},t}\propto\exp(-\lambda^{I}O_{I}). When {OI}={H(t)}\{O_{I}\}=\{H(t)\}, the coarse-grained state ρcg,t\rho_{\mathrm{cg},t} is exactly (2.3), so S(ρref,t)S(\rho_{\mathrm{ref},t}) is viewed as the coarse-grained entropy. Since the coarse-graining is a process of treating the system thermodynamically, the ansatz (2.3) is appropriate.

2.1.2 Generalization

We consider more general situations in a relativistic quantum field theory on dd-dimensional spacetime, whose metric is supposed to be static. The coordinate is written as za=(t,θA)z^{a}=(t,\theta^{A}) with tt being the time coordinate, and the metric as

ds2=σabdzadzb=dt2+σAB(θ)dθAdθB(σ:=detσab=detσAB).\displaystyle\mathrm{d}s^{2}=\sigma_{ab}\mathrm{d}z^{a}\mathrm{d}z^{b}=-\mathrm{d}t^{2}+\sigma_{AB}(\theta)\mathrm{d}\theta^{A}\mathrm{d}\theta^{B}\qquad(\sigma:=\mathrm{det}\,\sigma_{ab}=-\mathrm{det}\,\sigma_{AB}). (2.12)

The action is supposed to take the form

IQFT[w]=ddzσ(+wI(z)OI(z)),\displaystyle I_{\mathrm{QFT}}[w]=\int\mathrm{d}^{d}z\sqrt{-\sigma}\left(\mathcal{L}+w^{I}(z)O_{I}(z)\right), (2.13)

where \mathcal{L} is the Lagrangian density without explicit time-dependence, each OI(z)O_{I}(z) is a composite operator (without time-derivative), and wI(z)w^{I}(z) is a source coupling to it. The Hamiltonian operator is then given by

H(t)=Hdd1θσwI(t,θ)OI(θ),\displaystyle H(t)=H-\int\mathrm{d}^{d-1}\theta\sqrt{-\sigma}\,w^{I}(t,\theta)O_{I}(\theta), (2.14)

where HH is the Hamiltonian when w=0w=0. In this expression and below, we adopt the Schrödinger picture. Let Tab(θ)T_{ab}(\theta) be the stress tensor operator when w=0w=0. Then, the Hamiltonian HH and momentum PAP_{A} without sources are given as

H=dd1θσuaTbatb,PA=dd1θσuaTbaeAb,\displaystyle H=\int\mathrm{d}^{d-1}\theta\sqrt{-\sigma}\,u_{a}T^{a}_{b}t^{b},\qquad P_{A}=\int\mathrm{d}^{d-1}\theta\sqrt{-\sigma}\,u_{a}T^{a}_{b}e^{b}_{A}, (2.15)

where uau_{a} is the unit normal to time slices, taa=tt^{a}\partial_{a}=\partial_{t}, and eAbb=Ae^{b}_{A}\partial_{b}=\partial_{A}. For the metric (2.12), we have ua=tau^{a}=t^{a}.

At t=0t=0, we prepare the initial state as

ρ0=1Z0exp[β(HωAPAdd1θσjI(θ)OI(θ))],\displaystyle\rho_{0}=\frac{1}{Z_{0}}\exp\left[-\beta\left(H-\omega^{A}P_{A}-\int\mathrm{d}^{d-1}\theta\sqrt{-\sigma}\,j^{I}(\theta)O_{I}(\theta)\right)\right], (2.16)

with Z0Z_{0} determined through Trρ0=1\mathrm{Tr}\rho_{0}=1 and dependent on β0,ω0A\beta_{0},\omega^{A}_{0} and j0Ij^{I}_{0}. This state evolves according to (2.14):

ρt=U(t)ρ0U(t)1,U(t):=Texp(i0tdsH(s)).\displaystyle\rho_{t}=U(t)\rho_{0}U(t)^{-1},\qquad U(t):=\mathrm{T}\exp\left(-i\int_{0}^{t}\mathrm{d}s\,H(s)\right). (2.17)

We choose the reference state as

ρref,t=1Ztexp[βt(HωtAPAdd1θσjtI(θ)OI(θ))].\displaystyle\rho_{\mathrm{ref},t}=\frac{1}{Z_{t}}\exp\left[-\beta_{t}\left(H-\omega^{A}_{t}P_{A}-\int\mathrm{d}^{d-1}\theta\,\sqrt{-\sigma}j^{I}_{t}(\theta)O_{I}(\theta)\right)\right]. (2.18)

The parameters βt\beta_{t}, ωtA\omega_{t}^{A}, and jtI(θ)j_{t}^{I}(\theta) are the Lagrange multipliers to be optimized below, and as the initial conditions, we require

βt=0=β,ωt=0A=ωA,jt=0I=jI.\displaystyle\beta_{t=0}=\beta,\qquad\omega^{A}_{t=0}=\omega^{A},\qquad j^{I}_{t=0}=j^{I}. (2.19)

The relative entropy between ρt\rho_{t} and ρref,t\rho_{\mathrm{ref},t} is calculated as

S(ρt||ρref,t)=S(ρt)+Tr(ρtlnρref,t)=S(ρ0)+Tr(ρtlnρref,t)0.\displaystyle S(\rho_{t}||\rho_{\mathrm{ref},t})=-S(\rho_{t})+\mathrm{Tr}\left(\rho_{t}\ln\rho_{\mathrm{ref},t}\right)=-S(\rho_{0})+\mathrm{Tr}\left(\rho_{t}\ln\rho_{\mathrm{ref},t}\right)\geq 0. (2.20)

The tightest bound is achieved by optimizing this with respect to (βt,ωt,jt)(\beta_{t},\omega_{t},j_{t}), and the conditions are found to be

Tr(ρtH)=Tr(ρref,tH),Tr(ρtPA)=Tr(ρref,tPA),Tr(ρtOI(θ))=Tr(ρref,tOI(θ)).\displaystyle\mathrm{Tr}\left(\rho_{t}H\right)=\mathrm{Tr}\left(\rho_{\mathrm{ref},t}H\right),~{}~{}\mathrm{Tr}\left(\rho_{t}P_{A}\right)=\mathrm{Tr}\left(\rho_{\mathrm{ref},t}P_{A}\right),~{}~{}\mathrm{Tr}\left(\rho_{t}O_{I}(\theta)\right)=\mathrm{Tr}\left(\rho_{\mathrm{ref},t}O_{I}(\theta)\right). (2.21)

Thus, we have again obtained the coarse-graining conditions. As explained before, conversely, maximizing the entropy with those conditions reveals that ρref,t\rho_{\mathrm{ref},t} must take the form of (2.18). Under these conditions, (2.20) is reduced to

S(ρref,t)S(ρ0).\displaystyle S(\rho_{\mathrm{ref},t})\geq S(\rho_{0}). (2.22)

Thus, our coarse-grained entropy S(ρref,t)S(\rho_{\mathrm{ref},t}) never gets smaller than the initial value. Note that the initial state is not arbitrary.

Although not written down, the set of differential equations for βt\beta_{t}, ωtA\omega^{A}_{t}, and jtIj^{I}_{t}, i.e, the differential equations to determine ρref,t\rho_{\mathrm{ref},t}, can be derived from the tt-derivative of (2.21). The solution is unique under the initial condition (2.19). However, we do not argue this point anymore, as the our target is the application to gravity.

2.1.3 Purification

For later convenience, we consider purifying ρt\rho_{t}. To purify ρ0\rho_{0}, we bring a copy of the QFT, and name the original one QFTR\mathrm{QFT}_{\mathrm{R}} and the copied one QFTL\mathrm{QFT}_{\mathrm{L}}. Let |ψ0\ket{\psi_{0}} be the pure state in QFTLQFTR\mathrm{QFT}_{\mathrm{L}}\otimes\mathrm{QFT}_{\mathrm{R}} given by

|ψ0=1Z0n|nLeβH~R/2|nR,\displaystyle\ket{\psi_{0}}=\frac{1}{\sqrt{Z_{0}}}\sum_{n}\ket{n}_{\mathrm{L}}\otimes e^{-\beta\tilde{H}_{\mathrm{R}}/2}\ket{n}_{\mathrm{R}}, (2.23)

where {|n}\{\ket{n}\} is any orthonormal basis, and H~\tilde{H} is defined to be,

H~:=HωAPAdd1θσjI(θ)OI(θ),\displaystyle\tilde{H}:=H-\omega^{A}P_{A}-\int\mathrm{d}^{d-1}\theta\sqrt{-\sigma}\,j^{I}(\theta)O_{I}(\theta), (2.24)

which appeared in the exponent of (2.16). We see ρ0=Tr|ψ0ψ0|\rho_{0}=\mathrm{Tr}\ket{\psi_{0}}\bra{\psi_{0}}, and hence |ψ0\ket{\psi_{0}} is the purified state of ρ0\rho_{0}. Here, note that H~\tilde{H} is Hermitian due to ρ0=ρ0\rho_{0}^{\dagger}=\rho_{0}.

The time evolution of the total system is defined by

ULR(t)=VL(t)UR(t),\displaystyle U_{\mathrm{LR}}(t)=V_{L}(t)\otimes U_{R}(t), (2.25)

where UR(t)U_{\mathrm{R}}(t) is the one in (2.17), and VL(t)V_{\mathrm{L}}(t) is any unitary operator. With this, |ψ0\ket{\psi_{0}} evolves to

|ψt:=ULR(t)|ψ0=1Z0nVL(t)|nLUR(t)eβH~R/2|nR.\displaystyle\ket{\psi_{t}}:=U_{\mathrm{LR}}(t)\ket{\psi_{0}}=\frac{1}{\sqrt{Z_{0}}}\sum_{n}V_{\mathrm{L}}(t)\ket{n}_{\mathrm{L}}\otimes U_{\mathrm{R}}(t)e^{-\beta\tilde{H}_{\mathrm{R}}/2}\ket{n}_{\mathrm{R}}. (2.26)

For any VLV_{\mathrm{L}} and any operator OO, we can show

TrL|ψtψt|=ρt,\displaystyle\mathrm{Tr}_{\mathrm{L}}\ket{\psi_{t}}\bra{\psi_{t}}=\rho_{t}, (2.27)
ψt|OR|ψt=Tr(ρtO).\displaystyle\braket{\psi_{t}}{O_{\mathrm{R}}}{\psi_{t}}=\mathrm{Tr}\left(\rho_{t}O\right). (2.28)

Since our focus is QFTR\mathrm{QFT}_{\mathrm{R}} and VLV_{\mathrm{L}} is any, we hereafter set

V(t)=eiHt,\displaystyle V(t)=e^{-iHt}, (2.29)

while continuing using the notation VV even below.

2.2 Path integral representation

Here, we write down the generating functionals for ρref,t\rho_{\mathrm{ref},t} and (2.28) to use the GKPW formula in section 3.

2.2.1 Coarse-grained state

Let us start with

ρ[Γ]:=1Z[Γ]exp[B(HΩAPAdd1θσJI(θ)OI(θ))],\displaystyle\rho[\Gamma]:=\frac{1}{Z[\Gamma]}\exp\left[-B\left(H-\Omega^{A}P_{A}-\int\mathrm{d}^{d-1}\theta\,\sqrt{-\sigma}J^{I}(\theta)O_{I}(\theta)\right)\right], (2.30)

where Γ:=(B,Ω,J)\Gamma:=(B,\Omega,J). First, notice that

HΩAPA=dd1θσuaTbaξb,ξa:=uaΩAeAa.\displaystyle H-\Omega^{A}P_{A}=\int\mathrm{d}^{d-1}\theta\sqrt{-\sigma}\,u_{a}T^{a}_{b}\xi^{b},\qquad\xi^{a}:=u^{a}-\Omega^{A}e^{a}_{A}. (2.31)

From the viewpoint of the ADM formalism, this is the Hamiltonian when the time direction is chosen as ξa\xi^{a}. In other words, it is the Hamiltonian on the background given by

ds2=dt2+σAB(dθAΩAdt)(dθBΩBdt).\displaystyle\mathrm{d}s^{2}=-\mathrm{d}t^{2}+\sigma_{AB}(\mathrm{d}\theta^{A}-\Omega^{A}\mathrm{d}t)(\mathrm{d}\theta^{B}-\Omega^{B}\mathrm{d}t). (2.32)

Then, the partition function in (2.30) is expressed as

Z[Γ]=𝒟φexp[0Bdτdd1θσ~(EJI(θ)OI(τ,θ))],\displaystyle Z[\Gamma]=\oint\mathcal{D}\varphi\,\exp\left[-\int_{0}^{B}\mathrm{d}\tau\int d^{d-1}\theta\sqrt{\tilde{\sigma}}\left(\mathcal{L}_{\mathrm{E}}-J^{I}(\theta)O_{I}(\tau,\theta)\right)\right], (2.33)

where E\mathcal{L}_{\mathrm{E}} is the Euclid Lagrangian on

ds2=σ~abdzadzb=dτ2+σAB(dθA+iΩAdτ)(dθB+iΩBdτ).\displaystyle\mathrm{d}s^{2}=\tilde{\sigma}_{ab}\mathrm{d}z^{a}\mathrm{d}z^{b}=\mathrm{d}\tau^{2}+\sigma_{AB}(\mathrm{d}\theta^{A}+i\Omega^{A}\mathrm{d}\tau)(\mathrm{d}\theta^{B}+i\Omega^{B}\mathrm{d}\tau). (2.34)

In the path integral, φ\varphi means the collection of the elemental fields, and the (anti-)periodic boundary condition is imposed on φ\varphi as τ+Bτ\tau+B\sim\tau. The normalization factor ZtZ_{t} in (2.18) is equal to Z[Γt]Z[\Gamma_{t}] with Γt=(βt,ωt,jt)\Gamma_{t}=(\beta_{t},\omega_{t},j_{t}). The expectation values are generated as

Tr(ρ[Γ]OI(θ))=1Bσ~δδJI(θ)lnZ[Γ].\displaystyle\mathrm{Tr}\left(\rho[\Gamma]O_{I}(\theta)\right)=\frac{1}{B\sqrt{\tilde{\sigma}}}\frac{\delta}{\delta J^{I}(\theta)}\ln Z[\Gamma]. (2.35)

By a coordinate transformation from θA\theta^{A} to ϑA=θA+iΩAτ\vartheta^{A}=\theta^{A}+i\Omega^{A}\tau, the metric is changed to

ds2=dτ2+σABdϑAdϑB.\displaystyle\mathrm{d}s^{2}=\mathrm{d}\tau^{2}+\sigma_{AB}\mathrm{d}\vartheta^{A}\mathrm{d}\vartheta^{B}. (2.36)

In this coordinate, the periodicity condition for τ\tau is modified to (τ+B,ϑA+iΩAB)(τ,ϑA)(\tau+B,\vartheta^{A}+i\Omega^{A}B)\sim(\tau,\vartheta^{A}), and additionally, the source term could depend on tt. In the following, we will rather use the coordinate (2.34). The choice of (2.36) is discussed in section 3.5.2.

Refer to caption

==          Refer to caption

Figure 1: (Left) The evolution of (2.23) to (2.26). (Right) The contour CC in (2.38) to give the wave function of (2.26). The two vertical lines are the Lorentzian time-evolution over [0,t][0,t] by VV and UU, while the half curve is the Euclidean evolution over [0,β/2][0,\beta/2]. In the figure, the Lorentzian time direction of VV is drawn opposite to UU, to match with figure 3.

2.2.2 Original state

Next, we find the path integral representation of the functional to generate (2.28). As will be shown in appendix A, the following holds for |ψt\ket{\psi_{t}}555 Equation (2.37) does not hold as it is, but as explained in appendix A, there is no problem in proceeding with (2.37). :

ϕ|Lϕ|RVL(t)UR(t)(n|nLeβH~R/2|nR)=ϕ|U(t)eβH~/2V(t)|ϕ.\displaystyle{}_{\mathrm{L}}\bra{\phi}\otimes{}_{\mathrm{R}}\bra{\phi^{\prime}}V_{\mathrm{L}}(t)\otimes U_{\mathrm{R}}(t)\left(\sum_{n}\ket{n}_{\mathrm{L}}\otimes e^{-\beta\tilde{H}_{\mathrm{R}}/2}\ket{n}_{\mathrm{R}}\right)=\braket{\phi^{\prime}}{U(t)e^{-\beta\tilde{H}/2}V(t)}{\phi}. (2.37)

Note that the l.h.s. is proportional to the wave function of the thermofield double state (2.26), and that the r.h.s. is written in terms of a single QFT. The r.h.s. can be expressed as

ϕ|U(t)eβH~/2V(t)|ϕ=ϕϕ𝒟φeiIQFT[φ;C],\displaystyle\braket{\phi^{\prime}}{U(t)e^{-\beta\tilde{H}/2}V(t)}{\phi}=\int_{\phi}^{\phi^{\prime}}\mathcal{D}\varphi\,e^{iI_{\mathrm{QFT}}[\varphi;C]}, (2.38)

where CC is the contour depicted in figure 1; first, +t+t in the real direction, then iβ-i\beta in the pure imaginary direction, and finally +t+t again in the real direction. The metric of the Lorentzian parts is (2.12), and the one for the Euclidean part is (2.34) with the replacement ΩAωA\Omega^{A}\to\omega^{A}. Regarding the continuity, the induced metric on constant time surface is the same between both metrics, but the extrinsic curvature is not analytically continuous. This is not a problem however, because there is no dynamical gravity on the QFT.

Refer to caption
Figure 2: The contour CC^{\prime} for the generating functional (2.39). The source ww^{\prime} in the inverse contour is different from ww, and UU^{\prime} is the time-evolution operator with source ww^{\prime}. Note again that the Lorentzian time direction of VV and V1V^{-1} is taken opposite to UU and U1U^{\prime-1}.

The generating functional to compute (2.28) is the path integral over the closed contour CC^{\prime} in figure 2:

ZC[w,w]=𝒟φeiIQFT[φ;C].\displaystyle Z_{C^{\prime}}[w,w^{\prime}]=\oint\mathcal{D}\varphi\,e^{iI_{\mathrm{QFT}}[\varphi;C^{\prime}]}. (2.39)

In the arguments, jj and ω\omega have been omitted, as they will be fixed hereafter. Since CC^{\prime} is closed, the path integral is subject to the (anti-)periodic boundary condition. This path integral corresponds to

Tr[(U(t)eβH~/2V(t))U(t)eβH~/2V(t)],\displaystyle\mathrm{Tr}\left[\left(U^{\prime}(t)e^{-\beta\tilde{H}/2}V(t)\right)^{\dagger}U(t)e^{-\beta\tilde{H}/2}V(t)\right], (2.40)

and thus when w=ww^{\prime}=w, it is equal to Z0Z_{0}, i.e, ZC[w,w]=Z0Z_{C^{\prime}}[w,w]=Z_{0}. The relation we need is

Tr(ρtOI(θ))=1iσδδwI(t,θ)lnZC[w,w]|ww.\displaystyle\mathrm{Tr}\left(\rho_{t}O_{I}(\theta)\right)=\left.\frac{1}{i\sqrt{-\sigma}}\frac{\delta}{\delta w^{I}(t,\theta)}\ln Z_{C^{\prime}}[w,w^{\prime}]\right|_{w^{\prime}\to w}. (2.41)

To compute the expectation values of HH and PAP_{A}, we need to include the stress tensor Tab(z)T_{ab}(z) in the source term of the Lagrangian, coupled to the source δσab\delta\sigma^{ab} as (1/2)δσabTab(1/2)\delta\sigma^{ab}T_{ab}. This corresponds to that the Lagrangian (2.13), which defines the evolution U(t)U(t), is replaced with

IQFT[w,δσ]\displaystyle I_{\mathrm{QFT}}[w,\delta\sigma] =IQFT[w]12ddzσδσab(z)Tab(z)\displaystyle=I_{\mathrm{QFT}}[w]-\frac{1}{2}\int\mathrm{d}^{d}z\sqrt{-\sigma}\,\delta\sigma^{ab}(z)T_{ab}(z)
=ddzdet(σab+δσab)[σab+δσab]+ddzσwI(z)OI(z),\displaystyle=\int\mathrm{d}^{d}z\,\sqrt{-\mathrm{det}\left(\sigma_{ab}+\delta\sigma_{ab}\right)}\mathcal{L}[\sigma_{ab}+\delta\sigma_{ab}]+\int\mathrm{d}^{d}z\sqrt{-\sigma}\,w^{I}(z)O_{I}(z), (2.42)

where δσab:=σacσbdδσcd\delta\sigma_{ab}:=-\sigma_{ac}\sigma_{bd}\delta\sigma^{cd}. Note that the metric in the source term is not perturbed. In the following, we include δσab/2-\delta\sigma^{ab}/2 in {wI}\{w^{I}\}, so that we can deal with them collectively. In the expression (2.41), we take δσ0\delta\sigma\to 0 together with www^{\prime}\to w.

3 Description in gravity

Since the path integral representations (2.33) and (2.39) have been obtained, the GKPW formula enables us to write them in terms of gravity. Especially, (2.33) corresponds to a classical Euclidean gravity solution, and we will derive its thermodynamic laws. Our coarse-grained entropy is, in the classical limit, equal to the horizon area of the Euclidean cigar geometry (figure 4). The application to non-AdS cases is discussed at the end of this section.

3.1 Setup

Here, a dd-dimensional holographic CFT is considered. The contents in section 2 apply as they are, and we use the holographic dictionary to rewrite them in the gravitational language. Particularly, we find the dual descriptions of the conditions (2.21) and the coarse-grained entropy S(ρref,t)S(\rho_{\mathrm{ref},t}).

The bulk configuration dual to (2.39) is constructed by the common procedure [60, 61]. As figure 12 of [61], the resulting spacetime fills the bulk of figure 2 as shown in figure 3. (In this figure, only the forward half is drawn, and it is up to the time slice drawn red that corresponds to the forward half of figure 2.) Our target is the unshaded single-sided black hole spacetime.

Refer to caption
Figure 3: The bulk spacetime. The half disc depicts the Euclidean time evolution over [0,β/2][0,\beta/2] to prepare the initial configuration at t=0t=0, which evolves to the Lorentzian part drawn vertically. The dotted diagonal lines are the event horizons. This classical solution is constructed by properly joining the saddle point of each part [61].

Our setup and assumptions in the dual gravity are as follows:

  • We consider that the bulk MM is governed by the Einstein gravity,

    Igrav=116πGMdd+1xg(R2Λ)+18πGMddzγK+Imatter+Ict,\displaystyle I_{\mathrm{grav}}=\frac{1}{16\pi G}\int_{M}\mathrm{d}^{d+1}x\sqrt{-g}\,(R-2\Lambda)+\frac{1}{8\pi G}\int_{\partial M}\mathrm{d}^{d}z\sqrt{-\gamma}\,K+I_{\mathrm{matter}}+I_{\mathrm{ct}}, (3.1)

    where ImatterI_{\mathrm{matter}} is the matter Lagrangian, IctI_{\mathrm{ct}} is the counterterm both for gravity and matter, and KabK_{ab} is the extrinsic curvature of the boundary and KK is its trace, K:=γabKabK:=\gamma^{ab}K_{ab}. We assume that ImatterI_{\mathrm{matter}} does not contain any boundary term.

  • We write the bulk metric as

    ds2=gμνdxμdxν=N2dt2+Nr2dr2+cAB(dθA+NAdt)(dθB+NBdt),\displaystyle\mathrm{d}s^{2}=g_{\mu\nu}\mathrm{d}x^{\mu}\mathrm{d}x^{\nu}=-N^{2}\mathrm{d}t^{2}+N_{r}^{2}\mathrm{d}r^{2}+c_{AB}\left(\mathrm{d}\theta^{A}+N^{A}\mathrm{d}t\right)\left(\mathrm{d}\theta^{B}+N^{B}\mathrm{d}t\right), (3.2)

    with the gauge-fixing condition grt=grA=0g_{rt}=g_{rA}=0. Throughout this paper, μ,ν,\mu,\nu,\cdots are used for the bulk spacetime coordinates, A,B,A,B,\cdots are for θ\theta-coordinates, and a,b,a,b,\cdots are for the combination of (t,θA)(t,\theta^{A}).

  • With the unit normal 𝚞μμt\mathtt{u}_{\mu}\propto\partial_{\mu}t and the pull-back eAμe_{A}^{\mu}, the vector of time direction is expressed as

    𝚝μμ:=t=(N𝚞μ+NAeAμ)μ=N𝚞aa+NAA.\displaystyle\mathtt{t}^{\mu}\partial_{\mu}:=\partial_{t}=(N\mathtt{u}^{\mu}+N^{A}e^{\mu}_{A})\partial_{\mu}=N\mathtt{u}^{a}\partial_{a}+N^{A}\partial_{A}. (3.3)
  • The boundary M\partial M is located at rr\to\infty, where the induced metric of the Lorentzian part becomes (2.12) after the conformal factor is removed:

    ds2r2L2(dt2+σABdθAdθB)+Nr2dr2.\displaystyle\mathrm{d}s^{2}\sim\frac{r^{2}}{L^{2}}\left(-\mathrm{d}t^{2}+\sigma_{AB}\mathrm{d}\theta^{A}\mathrm{d}\theta^{B}\right)+N_{r}^{2}\mathrm{d}r^{2}. (3.4)

    The induced metric on r=constr=\mathrm{const} is γab\gamma_{ab}, which asymptotes to the first term above.

  • The on-shell value of the action is expressed as Igrav[w,w]I_{\mathrm{grav}}[w,w^{\prime}]. The argument ww is the one appeared in (2.39), and related to the asymptotic value of the bulk fields ΦI\Phi^{I} (including gμνg_{\mu\nu}) as

    ΦI(x)(RL)ΔIdwI(z)=:WI(z),\displaystyle\Phi^{I}(x)\sim\left(\frac{R}{L}\right)^{\Delta_{I}-d}w^{I}(z)=:W^{I}(z), (3.5)

    where RR is the IR cut-off to be sent RR\to\infty, and ΔI\Delta_{I} is the conformal dimension of OIO_{I}.666 In the bulk definition, ΔI\Delta_{I} is determined by the asymptotic rr-dependence of ΦI\Phi^{I}. The argument ww^{\prime} is related to WW^{\prime} in the same say in the backward evolution of figure 2.

  • The counterterm IctI_{\mathrm{ct}} is usually determined as [62, 63].

3.2 Coarse-graining conditions in the bulk

According to the the GKPW formula, (2.41) is rewritten as

Tr(ρtOI(θ))=limR(RL)ΔI1γδIgrav[w,w]δWI(t,θ)|ww.\displaystyle\mathrm{Tr}\left(\rho_{t}O_{I}(\theta)\right)=\lim_{R\to\infty}\left(\frac{R}{L}\right)^{\Delta_{I}}\left.\frac{1}{\sqrt{-\gamma}}\frac{\delta I_{\mathrm{grav}}[w,w^{\prime}]}{\delta W^{I}(t,\theta)}\right|_{w^{\prime}\to w}. (3.6)

For the stress tensor, in particular, we have

Tr(ρtTab(θ))=limR(RL)d2Yab(z),\displaystyle\mathrm{Tr}\left(\rho_{t}T_{ab}(\theta)\right)=\lim_{R\to\infty}\left(\frac{R}{L}\right)^{d-2}Y_{ab}(z), (3.7)

where YabY_{ab} is the Brown-York tensor defined as [62]

Yab:=2γδIgravδγab=18πG(KγabKab)2γδIctδγab.\displaystyle Y_{ab}:=-\frac{2}{\sqrt{-\gamma}}\frac{\delta I_{\mathrm{grav}}}{\delta\gamma^{ab}}=\frac{1}{8\pi G}\left(K\gamma_{ab}-K_{ab}\right)-\frac{2}{\sqrt{-\gamma}}\frac{\delta I_{\mathrm{ct}}}{\delta\gamma^{ab}}. (3.8)

Therefore, all we need is

δwIgrav[w]:=\displaystyle\delta_{w}I_{\mathrm{grav}}[w]:= Igrav[w+δw,w]Igrav[w,w]\displaystyle I_{\mathrm{grav}}[w+\delta w,w]-I_{\mathrm{grav}}[w,w]
=\displaystyle= dd+1x(EOM)IδwΦI+ddzγπI(z)δWI(z).\displaystyle\int\mathrm{d}^{d+1}x\,(\mathrm{EOM})_{I}\delta_{w}\Phi^{I}+\int\mathrm{d}^{d}z\sqrt{-\gamma}\,\pi_{I}(z)\delta W^{I}(z). (3.9)

Since we are interested in the on-shell variation, only the boundary terms remain. By definition, the Brown-York tensor is considered as

πI=YabforδWI=δγab/2.\displaystyle\pi_{I}=Y_{ab}\qquad\mathrm{for}\qquad\delta W^{I}=-\delta\gamma^{ab}/2. (3.10)

When ImatterI_{\mathrm{matter}} takes the form dd+1xgmatter(ΦI,μΦI;g)\int\mathrm{d}^{d+1}x\sqrt{-g}\,\mathcal{L}_{\mathrm{matter}}(\Phi^{I^{\prime}},\nabla_{\mu}\Phi^{I^{\prime}};g) (II^{\prime} does not contain gμνg_{\mu\nu}), πI\pi_{I^{\prime}} for the matter fields are written as

πI=rμmatter(μΦI)+1γδIctδWIforδWI=δΦI|r=R,\displaystyle\pi_{I^{\prime}}=r_{\mu}\frac{\partial\mathcal{L}_{\mathrm{matter}}}{\partial\left(\nabla_{\mu}\Phi^{I^{\prime}}\right)}+\frac{1}{\sqrt{-\gamma}}\frac{\delta I_{\mathrm{ct}}}{\delta W^{I^{\prime}}}\qquad\mathrm{for}\qquad\delta W^{I^{\prime}}=\delta\Phi^{I^{\prime}}|_{r=R}, (3.11)

where rμ=Nr1δμrr_{\mu}=N_{r}^{-1}\delta_{\mu}^{r} is the unit normal of r=constr=\mathrm{const}. As seen above, πI\pi_{I} is like the “canonical conjugate momentum” (up to counterterm contributions) if rr were viewed as the time coordinate. This viewpoint agrees with the fact that differentiating the on-shell action with the final position generates the canonical conjugate momentum.

Therefore, (3.6) is reduced to

Tr(ρtOI(θ))=limR(RL)ΔIπI(t,θ)=:oI,t(θ),\displaystyle\mathrm{Tr}\left(\rho_{t}O_{I^{\prime}}(\theta)\right)=\lim_{R\to\infty}\left(\frac{R}{L}\right)^{\Delta_{I^{\prime}}}\pi_{I^{\prime}}(t,\theta)=:o_{I^{\prime},t}(\theta), (3.12)

and (3.9) reproduces the on-shell variation formula:

δwIgrav[w]=ddzσoI,t(θ)δwI(t,θ)12ddzγYab(t,θ)δγab(t,θ).\displaystyle\delta_{w}I_{\mathrm{grav}}[w]=\int\mathrm{d}^{d}z\sqrt{-\sigma}\,o_{I^{\prime},t}(\theta)\delta w^{I^{\prime}}(t,\theta)-\frac{1}{2}\int\mathrm{d}^{d}z\sqrt{-\gamma}\,Y_{ab}(t,\theta)\delta\gamma^{ab}(t,\theta). (3.13)

For HH and PAP_{A}, the following holds from (2.15) and (3.7):

Tr(ρtH)=limR(RL)d2dd1θσuaYabtb=limRdd1θcN𝚞aYab𝚞b=:ht,\displaystyle\mathrm{Tr}\left(\rho_{t}H\right)=\lim_{R\to\infty}\left(\frac{R}{L}\right)^{d-2}\int\mathrm{d}^{d-1}\theta\sqrt{-\sigma}\,u^{a}Y_{ab}t^{b}=\lim_{R\to\infty}\int\mathrm{d}^{d-1}\theta\sqrt{c}\,N\mathtt{u}^{a}Y_{ab}\mathtt{u}^{b}=:h_{t}, (3.14)
Tr(ρtPA)=limR(RL)d2dd1θσuaYabeAb=limRdd1θc𝚞aYabeAb=:pA,t.\displaystyle\mathrm{Tr}\left(\rho_{t}P_{A}\right)=\lim_{R\to\infty}\left(\frac{R}{L}\right)^{d-2}\int\mathrm{d}^{d-1}\theta\sqrt{-\sigma}\,u^{a}Y_{ab}e^{b}_{A}=\lim_{R\to\infty}\int\mathrm{d}^{d-1}\theta\sqrt{c}\,\mathtt{u}^{a}Y_{ab}e^{b}_{A}=:p_{A,t}. (3.15)

Here, we have used ua=taN𝚞au^{a}=t^{a}\sim N\mathtt{u}^{a} as rr\to\infty, which follows from (3.4).

It should be noted that oIo_{I^{\prime}} and YabY_{ab} are evaluated by Igrav[w,w]I_{\mathrm{\mathrm{grav}}}[w,w] (ww^{\prime} is set to be ww). According to [61], the bulk in such a case will be simply constructed as follows. First, we cut into half the dominant Euclidean black hole solution dual to Z0Z_{0} (see also next paragraph), which provides the Euclidean “half disc” in figure 3. Next, we solve the initial value problem from t=0t=0, with the initial configuration given by analytically continuing the Euclidean solution. The initial values obtained in this way are compatible with constrains (e.g. Hamiltonian constraint). In solving this problem toward time tt, the fields are subject to the Neumann boundary conditions specified by wIw^{I} on M\partial M. The backward evolution part is simply obtained by time-reversing it.

Next, we move on to the dual gravitational description for ρref,t\rho_{\mathrm{ref},t} in (2.18). We assume that the bulk solution that dominates Z[Γ]Z[\Gamma] does not depend on the imaginary time (i.e, stationary), and let Igrav(E)[Γ]I_{\mathrm{grav}}^{(\mathrm{E})}[\Gamma] be the bulk Euclidean action evaluated by the dominant saddle point for given Γ\Gamma.777 Although phase transitions can happen depending on Γ\Gamma, such as the Hawking-Page transition [64], we formally perform the following calculations. Similarly to the Lorentzian case, the same formulae hold for the Euclidean solution:

Tr(ρ[Γ]OI(θ))=limR(RL)ΔIπI(θ)=:oI[Γ](θ),\displaystyle\mathrm{Tr}\left(\rho[\Gamma]O_{I^{\prime}}(\theta)\right)=\lim_{R\to\infty}\left(\frac{R}{L}\right)^{\Delta_{I^{\prime}}}\pi_{I^{\prime}}(\theta)=:o_{I^{\prime}}[\Gamma](\theta), (3.16)
Tr(ρ[Γ]H)=limRdd1θc~N𝚞aYab𝚞b=:h[Γ],\displaystyle\mathrm{Tr}\left(\rho[\Gamma]H\right)=\lim_{R\to\infty}\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{c}}\,N\mathtt{u}^{a}Y_{ab}\mathtt{u}^{b}=:h[\Gamma], (3.17)
Tr(ρ[Γ]PA)=limRdd1θc~𝚞aYabeAb=:pA[Γ],\displaystyle\mathrm{Tr}\left(\rho[\Gamma]P_{A}\right)=\lim_{R\to\infty}\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{c}}\,\mathtt{u}^{a}Y_{ab}e^{b}_{A}=:p_{A}[\Gamma], (3.18)

where, πI\pi_{I^{\prime}} here is defined as

πI(θ):=1Bγ~δIgrav(E)[Γ]δ𝒥I(θ),𝒥I(θ):=(RL)ΔIdJI(θ).\displaystyle\pi_{I^{\prime}}(\theta):=-\frac{1}{B\sqrt{\tilde{\gamma}}}\frac{\delta I_{\mathrm{grav}}^{\left(\mathrm{E}\right)}[\Gamma]}{\delta\mathcal{J}^{I}(\theta)},\qquad\mathcal{J}^{I}(\theta):=\left(\frac{R}{L}\right)^{\Delta_{I^{\prime}}-d}J^{I}(\theta). (3.19)

The tildes on the metrics indicates the Euclidean signature. As opposed to the Lorentzian case, the time vector is twisted as τa=ua+iΩAeAa\tau^{a}=u^{a}+i\Omega^{A}e_{A}^{a} due to Ω\Omega in (2.34). However, in order to measure HH rather than the Hamiltonian that generates τ\tau-translation of (2.34), we had to choose uaTabubu^{a}T_{ab}u^{b}, not uaTabτbu^{a}T_{ab}\tau^{b}. This point has already been taken into account in (3.17).

From the above, (2.21) is equivalent in the bulk to finding the set Γ=(B,Ω,J)\Gamma=(B,\Omega,J) such that

oI[Γ]=oI,th[Γ]=ht,pA[Γ]=pA,t.\displaystyle o_{I^{\prime}}[\Gamma]=o_{I^{\prime},t}\qquad h[\Gamma]=h_{t},\qquad p_{A}[\Gamma]=p_{A,t}. (3.20)

We write the solution as Γt=(βt,ωt,jt)\Gamma_{t}=(\beta_{t},\omega_{t},j_{t}). The first condition is equivalent to equating the leading behaviors of πI\pi_{I}.

3.3 Coarse-grained entropy

After Γt\Gamma_{t} is found for each time tt, our coarse-grained entropy can be computed in the gravity. We write the bulk metric as

ds2=g~μνdxμdxν=N~2dτ2+N~r2dr2+c~AB(dθA+iN~Adτ)(dθB+iN~Bdτ),\displaystyle\mathrm{d}s^{2}=\tilde{g}_{\mu\nu}\mathrm{d}x^{\mu}\mathrm{d}x^{\nu}=\tilde{N}^{2}\mathrm{d}\tau^{2}+\tilde{N}_{r}^{2}\mathrm{d}r^{2}+\tilde{c}_{AB}\left(\mathrm{d}\theta^{A}+i\tilde{N}^{A}\mathrm{d}\tau\right)\left(\mathrm{d}\theta^{B}+i\tilde{N}^{B}\mathrm{d}\tau\right), (3.21)

similar to (3.2). We write the metric on r=constr=\mathrm{const} as γ~ab\tilde{\gamma}_{ab}. This metric asymptotes to (2.34) as

N~rL,N~AΩA,c~ABr2L2σAB,(r).\displaystyle\tilde{N}\sim\frac{r}{L},\qquad\tilde{N}^{A}\sim\Omega^{A},\qquad\tilde{c}_{AB}\sim\frac{r^{2}}{L^{2}}\sigma_{AB},\qquad(r\to\infty). (3.22)

We suppose N~r|r=rh=0\tilde{N}_{r}|_{r=r_{h}}=0 and rN~r|r=rh0\partial_{r}\tilde{N}_{r}|_{r=r_{h}}\neq 0 for some rh>0r_{h}>0, as figure 4.888 If N~0\tilde{N}\neq 0 anywhere, the entropy just vanishes and such a case is not our interest. In this case by a coordinate transformation dρ=Nrdr\mathrm{d}\rho=N_{r}\mathrm{d}r (ρ(rh)=0\rho(r_{h})=0), the metric near r=rhr=r_{h} must behave as

ds2dρ2+(2πBρ)2dτ2+c~ABdθAdθB,\displaystyle\mathrm{d}s^{2}\sim\mathrm{d}\rho^{2}+\left(\frac{2\pi}{B}\rho\right)^{2}\mathrm{d}\tau^{2}+\tilde{c}_{AB}\mathrm{d}\theta^{A}\mathrm{d}\theta^{B}, (3.23)

so as to avoid the conical singularity. Note that the shift vectors must also vanish (see below (2.34) and also section 3.5.2).

Refer to caption
Figure 4: The Euclidean cigar geometry and a time slice τ=const\tau=\mathrm{const}.

First, the entropy of ρ[Γ]\rho[\Gamma] in (2.30) is

S(ρ[Γ])\displaystyle S(\rho[\Gamma]) =Tr(ρ[Γ]lnρ[Γ])\displaystyle=-\mathrm{Tr}\left(\rho[\Gamma]\ln\rho[\Gamma]\right)
=lnZ[Γ]+BTr(ρ[Γ]{HΩAPAdd1θσJI(θ)OI(θ)})\displaystyle=\ln Z[\Gamma]+B\,\mathrm{Tr}\left(\rho[\Gamma]\left\{H-\Omega^{A}P_{A}-\int\mathrm{d}^{d-1}\theta\,\sqrt{-\sigma}J^{I}(\theta)O_{I}(\theta)\right\}\right) (3.24)
=B2B(B1lnZ[Γ]).\displaystyle=B^{2}\frac{\partial}{\partial B}\left(-B^{-1}\ln Z[\Gamma]\right). (3.25)

Using the holographic dictionary, lnZ[Γ]=Igrav(E)[Γ]-\ln Z[\Gamma]=I_{\mathrm{grav}}^{(\mathrm{E})}[\Gamma], we obtain

S(ρ[Γ])=Igrav(E)[Γ]BIgrav(E)[Γ]B=Bd^BIgrav(E)[Γ],\displaystyle S(\rho[\Gamma])=I_{\mathrm{grav}}^{(\mathrm{E})}[\Gamma]-B\frac{\partial I_{\mathrm{grav}}^{(\mathrm{E})}[\Gamma]}{\partial B}=B\hat{\mathrm{d}}_{B}I_{\mathrm{grav}}^{(\mathrm{E})}[\Gamma], (3.26)

where d^B\hat{\mathrm{d}}_{B} is BB-derivative that acts only on the fields (with the integration range of τ\tau fixed) and the last equation follows from the τ\tau-independence of the on-shell fields. Thus, d^BIgrav(E)[Γ]\hat{\mathrm{d}}_{B}I_{\mathrm{grav}}^{(\mathrm{E})}[\Gamma] is just the field variation triggered by BB, and the Euclidean version of (3.9) can be used with δwd^B\delta_{w}\to\hat{\mathrm{d}}_{B}. Since the bulk fields are subject to boundary conditions JI(θ)J^{I}(\theta) independent of BB, (3.9) seems to vanish. But, the variation perturbs BB in (3.23), generating a non-trivial conical singularity. Due to this, we have [65]

116πGdd+1xd^B(g~R)=14GBrhdd1θc~+(regular terms),\displaystyle-\frac{1}{16\pi G}\int\mathrm{d}^{d+1}x\,\hat{\mathrm{d}}_{B}(\sqrt{\tilde{g}}\,R)=\frac{1}{4GB}\int_{r_{h}}\mathrm{d}^{d-1}\theta\sqrt{\tilde{c}}+\left(\mbox{regular terms}\right), (3.27)

by which we obtain

S(ρ[Γ])=14Grhdd1θc~.\displaystyle S\left(\rho[\Gamma]\right)=\frac{1}{4G}\int_{r_{h}}\mathrm{d}^{d-1}\theta\sqrt{\tilde{c}}. (3.28)

The fact that the entropy is given by the area of the tip agrees with [66], where the replica trick is used.

If necessary, quantum corrections can be computed as in [28, 29, 30, 35, 67]. Also, since (3.28) is derived from the gravitational path integral Z[Γ]Z[\Gamma], the non-perturbative effects such as phase transitions are automatically taken into account when the saddle point approximation is taken for Z[Γ]Z[\Gamma]. If the dominant saddle point changes, the entropy is then given by the corresponding Euclidean solution. In any case, the second law introduced below in (3.29) must be satisfied.

3.4 Thermodynamic laws

The second law is guaranteed via the AdS/CFT correspondence by (2.22):

S(ρ[Γt])S(ρ[Γ0]).\displaystyle S(\rho[\Gamma_{t}])\geq S(\rho[\Gamma_{0}]). (3.29)

Here, S(ρ[Γ0])S(\rho[\Gamma_{0}]) coincides with the horizon area of the initial black hole, as we have set ρ0=ρref,0\rho_{0}=\rho_{\mathrm{ref},0}. Therefore, the horizon area of the auxiliary Euclidean black hole never gets smaller than the initial value. This second law is a requirement from the AdS/CFT and is expected to provide a thermodynamic constraint on quantum gravity. As we will see in section 4, it seems to be related to some energy condition.

Next, let us derive the first law. To find it, we perform a scaling transformation τs=B1τ\tau_{s}=B^{-1}\tau. Values in the scaled coordinate is put subscript ss below. The periodicity is modified to τs+1τs\tau_{s}+1\sim\tau_{s}, and (3.22) and (3.23) are changed to

N~sBrL,N~sABΩA,(r),\displaystyle\tilde{N}_{s}\sim B\frac{r}{L},\qquad\tilde{N}^{A}_{s}\sim B\Omega^{A},\qquad(r\to\infty), (3.30)
ds2dρ2+(2πρ)2dτs2+c~ABdθAdθB(r=rh).\displaystyle\mathrm{d}s^{2}\sim\mathrm{d}\rho^{2}+(2\pi\rho)^{2}\mathrm{d}\tau_{s}^{2}+\tilde{c}_{AB}\mathrm{d}\theta^{A}\mathrm{d}\theta^{B}\qquad(r=r_{h}). (3.31)

In this coordinate, the derivative d^B\hat{\mathrm{d}}_{B} does not give any conical singularity, but this time it act on the boundary metric. In addition, the volume factors transformed as γ~s=Bγ~\sqrt{\tilde{\gamma}_{s}}=B\sqrt{\tilde{\gamma}} and σ~s=Bσ~\sqrt{\tilde{\sigma}_{s}}=B\sqrt{\tilde{\sigma}} while YabuaubY_{ab}u^{a}u^{b} and YabuaeAbY_{ab}u^{a}e^{b}_{A} are invariant. Noting those points, we take a variation δΓ\delta\Gamma for Igrav(E)[Γ]I_{\mathrm{grav}}^{(\mathrm{E})}[\Gamma] (see (2.35), (3.13) and (3.19)):

δΓIgrav(E)[Γ]=\displaystyle\delta_{\Gamma}I_{\mathrm{grav}}^{(\mathrm{E})}[\Gamma]= dd1θσ~soI[Γ](θ)δJI(θ)12dd1θγ~sYab(θ)δΓγab(θ),\displaystyle-\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{\sigma}_{s}}\,o_{I^{\prime}}[\Gamma](\theta)\delta J^{I^{\prime}}(\theta)-\frac{1}{2}\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{\gamma}_{s}}\,Y_{ab}(\theta)\delta_{\Gamma}\gamma^{ab}(\theta), (3.32)
δΓγab=\displaystyle\delta_{\Gamma}\gamma_{ab}= 2B(δB𝚞a𝚞bδ(BΩA)𝚞aeAb).\displaystyle\frac{2}{B}\left(\delta B\mathtt{u}_{a}\mathtt{u}_{b}-\delta(B\Omega^{A})\mathtt{u}_{a}e_{Ab}\right). (3.33)

Therefore, going back to the original coordinate, we obtain

δΓIgrav(E)[Γ]=h[Γ]δBpA[Γ]δ(BΩA)dd1θσ~o^I[Γ](θ)δJI(θ),\displaystyle\delta_{\Gamma}I_{\mathrm{grav}}^{(\mathrm{E})}[\Gamma]=h[\Gamma]\delta B-p_{A}[\Gamma]\delta(B\Omega^{A})-\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{\sigma}}\,\hat{o}_{I^{\prime}}[\Gamma](\theta)\delta J^{I^{\prime}}(\theta), (3.34)

where we have defined

o^I[Γ]:=BoI[Γ].\displaystyle\hat{o}_{I^{\prime}}[\Gamma]:=Bo_{I^{\prime}}[\Gamma]. (3.35)

In (3.34), we have supposed that {ΦI}\{\Phi^{I^{\prime}}\} includes only fields that are invariant under the scaling. Soon later, we will see what happens if there are scaling fields such as the Maxwell field.

The entropy function is obtained by the Legendre transformation from Γ\Gamma to Ξ:=(h,p,o^)\Xi:=(h,p,\hat{o}):

S[Ξ]:=\displaystyle S[\Xi]:= B[Ξ](hpAΩA[Ξ])dd1θσ~o^I(θ)JI[Ξ](θ)Igrav(E)[Γ[Ξ]]\displaystyle B[\Xi]\left(h-p_{A}\Omega^{A}[\Xi]\right)-\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{\sigma}}\,\hat{o}_{I^{\prime}}(\theta)J^{I^{\prime}}[\Xi](\theta)-I_{\mathrm{grav}}^{(\mathrm{E})}[\Gamma[\Xi]] (3.36)
δΞS[Ξ]=\displaystyle\delta_{\Xi}S[\Xi]= B[Ξ](δhΩA[Ξ]δpA)dd1θσ~δo^I(θ)JI[Ξ](θ).\displaystyle B[\Xi]\left(\delta h-\Omega^{A}[\Xi]\delta p_{A}\right)-\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{\sigma}}\,\delta\hat{o}_{I^{\prime}}(\theta)J^{I^{\prime}}[\Xi](\theta). (3.37)

Now, the entropy has a new expression (3.36), which is seemingly different from (3.28). But, they must be equivalent:

S(ρ[Γt])=S[Ξt],Ξt:=(ht,pt,o^t).\displaystyle S(\rho[\Gamma_{t}])=S[\Xi_{t}],\qquad\Xi_{t}:=(h_{t},p_{t},\hat{o}_{t}). (3.38)

As a matter of fact, (3.36) is exactly the same as the corresponding CFT expression (3.24). By construction of Γt\Gamma_{t}, the argument Ξt\Xi_{t} is the one given in (3.12), (3.14), and (3.15) up to the hat symbol. The first law follows from (3.37):

S˙[Ξt]=B[Ξt](h˙tΩA[Ξt]p˙A,t)dd1θσ~o^˙I,t(θ)JI[Ξ](θ).\displaystyle\dot{S}[\Xi_{t}]=B[\Xi_{t}]\left(\dot{h}_{t}-\Omega^{A}[\Xi_{t}]\dot{p}_{A,t}\right)-\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{\sigma}}\,\dot{\hat{o}}_{I^{\prime},t}(\theta)J^{I^{\prime}}[\Xi](\theta). (3.39)

This equation contains the local terms from matter fields, in addition to the usual first law.

Finally, let us see how (3.34) and below are modified when there are scaling fields. To be concrete, we demonstrate with the time component of U(1)\mathrm{U}(1) gauge field A0A_{0}, whose JIJ^{I} and oIo_{I} we write as aa and qq, respectively. In this case, (3.34) is modified to

δΓIgrav(E)[Γ]=(3.34)dd1θσ~q[Γ](θ)δ(Ba(θ)).\displaystyle\delta_{\Gamma}I_{\mathrm{grav}}^{(\mathrm{E})}[\Gamma]=\eqref{eq: unscaled variation}-\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{\sigma}}\,q[\Gamma](\theta)\delta(Ba(\theta)). (3.40)

Accordingly, qq is added to Ξ\Xi, and (3.36) and (3.39) acquire additional terms:

S[Ξ]:=\displaystyle S[\Xi]:= (3.36)B[Ξ]dd1θσ~q(θ)a[Ξ](θ),\displaystyle~{}\eqref{eq: entropy function}-B[\Xi]\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{\sigma}}\,q(\theta)a[\Xi](\theta), (3.41)
S˙[Ξt]=\displaystyle\dot{S}[\Xi_{t}]= (3.39)B[Ξt]dd1θσ~q˙t(θ)a[Ξt](θ).\displaystyle~{}\eqref{eq: first law}-B[\Xi_{t}]\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{\sigma}}\,\dot{q}_{t}(\theta)a[\Xi_{t}](\theta). (3.42)

The additional term is the usual term consisting of the U(1)\mathrm{U}(1) electric charge and the gauge potential, except that they are locally treated in the current case. If there are more fields scaled by τs=B1τ\tau_{s}=B^{-1}\tau, similar terms are added depending on their transformation rules.

3.5 Comments on generalizations

3.5.1 Non-AdS cases

We have obtained the definition of our coarse-grained entropy purely in the gravitational language. If gravity is in general holographic in the sense that the quantum theory behind it is defined by specifying the boundary as assumed in [66], then we would reach the same results in the classical limit from the saddle point approximation applied to the fomal path integral of the quantum gravity. Thus, it is expected that the first and second law also holds for non-AdS cases.

Let us consider the Einstein gravity (3.1) on a Lorentzian manifold MM with a timelike boundary M\partial M. The induced metric on M\partial M is assumed static, but the boundary conditions for the other bulk fields can depend on the time. Then, (3.20) determines the reference Euclidean solution on each time, whose horizon area is the entropy we want. In solving (3.20), the Brown-York tensor is defined in the same way, and the matching condition of oIo_{I} is replaced with the matching of the leading term of πI\pi_{I}.

3.5.2 Choice of shift vectors

So far, we have adopted (2.34) as our gauge choice. As has been explained, it is possible to choose (2.36), where ϑ\vartheta is shifted in the pure imaginary direction along with the τ\tau-periodicity. If we had taken this gauge, (3.23) would be changed to

ds2dρ2+(2πBρ)2dτ2+c~AB(dϑAiΩAdτ)(dϑBiΩBdτ),\displaystyle\mathrm{d}s^{2}\sim\mathrm{d}\rho^{2}+\left(\frac{2\pi}{B}\rho\right)^{2}\mathrm{d}\tau^{2}+\tilde{c}_{AB}(\mathrm{d}\vartheta^{A}-i\Omega^{A}\mathrm{d}\tau)(\mathrm{d}\vartheta^{B}-i\Omega^{B}\mathrm{d}\tau), (3.43)

from the requirement of removing the conical singularity. Seemingly different, but physically same results must be reproduced also in this gauge. A difference, for example, is that the (angular) momenta are measured not on the boundary, but at r=rhr=r_{h}. In fact, Ω\Omega does not appear in (3.33), but in (3.23) as above, and hence BB enters to (3.31) accompanying Ω\Omega when τ\tau is scaled.

It was addressed in [68] for Kerr-AdS5 spacetime, how the choice of the shift vectors affect on the first law. Since we allow matter fields to exist, the story seems to be not so simple as [68]. In vacuum solutions, for example, the angular momentum (Komar integral) measured at infinity is equal to the one at the horizon by the Stokes’ theorem, but different by the volume integral when matter exists. Nevertheless, the first law will be properly derived in a similar way.

4 Demonstration in Vaidya-type spacetimes

In this section, we check thermodynamic laws derived in the previous section for null-ray collapse models. Since the second law was predicted from the CFT, it provides an energy constraint for those examples.

4.1 Schwarzschild-AdSd+1

Refer to caption
Figure 5: Vaidya-Schwarzschild-AdSd+1 spacetime. At t=0t=0, the Schwarzschild-AdS spacetime is prepared, and null rays exist in the yellow region.

Let us start with the simplest case where a null-ray is shot from the boundary to the Schwarzschild-AdSd+1 spacetime. The initial metric is given as

ds2=f(r)dt2+dr2f(r)+r2dΘ2,f(r):=1+r2L22μ(0)rd2,\displaystyle\mathrm{d}s^{2}=-f(r)\mathrm{d}t^{2}+\frac{\mathrm{d}r^{2}}{f(r)}+r^{2}\mathrm{d}\Theta^{2},\qquad f(r):=1+\frac{r^{2}}{L^{2}}-\frac{2\mu(0)}{r^{d-2}}, (4.1)

where dΘ2\mathrm{d}\Theta^{2} is the metric of the unit (d1)(d-1)-dimensional sphere. The boundary metric with the conformal factor removed is

ds2|bdy=dt2+r2dΘ2.\displaystyle\mathrm{d}s^{2}|_{\mathrm{bdy}}=\mathrm{d}t^{2}+r^{2}\mathrm{d}\Theta^{2}. (4.2)

The metric describing the null ray is

ds2=f(v,r)dv2+2dvdr+r2dΘ2,f(v,r)=1+r2L22μ(v)rd2.\displaystyle\mathrm{d}s^{2}=-f(v,r)\mathrm{d}v^{2}+2\mathrm{d}v\mathrm{d}r+r^{2}\mathrm{d}\Theta^{2},\qquad f(v,r)=1+\frac{r^{2}}{L^{2}}-\frac{2\mu(v)}{r^{d-2}}. (4.3)

At the boundary, both time coordinates coincide, namely, t=vt=v. In the bulk, the time coordinate tt is properly extended towards the future, but we do not concretely specify the foliation, since only the boundary observables at time tt are required in coarse-graining. Hence, we use vv for the time coordinate hereafter.

We take the mass hvh_{v} as the only respected observable on the boundary.999 In principle, we can add boundary values of the matter fields making up the null ray, but it is in general difficult to explicitly construct null ray with some field [69]. The mass hvh_{v} depends on the counterterm [62, 63], but most of time [62], it is computed as

hv=d18πGVol(𝕊d1)μ(v)+(independent of v).\displaystyle h_{v}=\frac{d-1}{8\pi G}\mathrm{Vol}\left(\mathbb{S}^{d-1}\right)\mu(v)+\left(\mbox{independent of $v$}\right). (4.4)

The Euclidean black hole that minimizes the Euclidean action, having the same mass, will be

ds2=f(v,r)dτ2+dr2f(v,r)+r2dΘ2.\displaystyle\mathrm{d}s^{2}=f(v,r)\mathrm{d}\tau^{2}+\frac{\mathrm{d}r^{2}}{f(v,r)}+r^{2}\mathrm{d}\Theta^{2}. (4.5)

Note that τ\tau is not related to the Lorentzian time vv, but just the imaginary time for the reference Euclidean solution that we refer to at time vv. It is explicitly confirmed that this geometry has the mass hvh_{v}. The Lagrange multiplier βv\beta_{v} is found to be

βv=4πrf(t,rh(v)),\displaystyle\beta_{v}=\frac{4\pi}{\partial_{r}f(t,r_{h}(v))}, (4.6)

where r=rh(v)r=r_{h}(v) is the largest root of f(v,r)f(v,r), related to μ(v)\mu(v) as

μ(v)=12rh(v)d2(1+rh(v)2L2).\displaystyle\mu(v)=\frac{1}{2}r_{h}(v)^{d-2}\left(1+\frac{r_{h}(v)^{2}}{L^{2}}\right). (4.7)

The coarse-grained entropy is the area of the cigar tip,

Sv=Vol(𝕊d1)4Grh(v)d1.\displaystyle S_{v}=\frac{\mathrm{Vol}\left(\mathbb{S}^{d-1}\right)}{4G}r_{h}(v)^{d-1}. (4.8)

We see that the first law is satisfied:

S˙v=βvh˙v.\displaystyle\dot{S}_{v}=\beta_{v}\dot{h}_{v}. (4.9)

On the other hand, the second law SvS0S_{v}\geq S_{0} is equivalent to rh(v)rh(0)r_{h}(v)\geq r_{h}(0), which is not satisfied for generic μ(v)\mu(v). We rather see this fact as a constraint from quantum gravity, since this must be satisfied in the dual description.

Refer to caption
Figure 6: The behavior of f(v,r)f(v,r) under the null energy condition.

Actually, the second law is satisfied when the null energy condition holds. The energy flow along a null vector =v+(f(v,r)/2)r\ell=\partial_{v}+(f(v,r)/2)\partial_{r} is

T:=Tμνμνvf(v,r).\displaystyle T_{\ell\ell}:=T_{\mu\nu}\ell^{\mu}\ell^{\nu}\propto-\partial_{v}f(v,r). (4.10)

If T0T_{\ell\ell}\geq 0, then f(v,r)f(v,r) decreases in vv, which implies that rh(v)r_{h}(v) increases (see figure 6). However, this monotonicity is stronger than required by our second law, and the following integrated version is sufficient for only saying rh(v)rh(0)r_{h}(v)\geq r_{h}(0):

v>0,0vdvT0.\displaystyle\forall v>0,\qquad\int_{0}^{v}\mathrm{d}v\,T_{\ell\ell}\geq 0. (4.11)

4.2 Rotating BTZ

Next, we consider the uncharged BTZ spacetime and add the angular momentum pϕp_{\phi} to the respected observable set. The metric with null ray is given as

ds2\displaystyle\mathrm{d}s^{2} =N(v,r)2dv2+2dvdr+r2(dϕ+Nϕ(v,r)dv)2,\displaystyle=-N(v,r)^{2}\mathrm{d}v^{2}+2\mathrm{d}v\mathrm{d}r+r^{2}(\mathrm{d}\phi+N^{\phi}(v,r)\mathrm{d}v)^{2}, (4.12)
N(v,r)2\displaystyle N(v,r)^{2} =8Gm(v)+r2L2+16G2j(v)2r2,Nϕ(v,r)=4Gj(v)r2.\displaystyle=-8Gm(v)+\frac{r^{2}}{L^{2}}+\frac{16G^{2}j(v)^{2}}{r^{2}},\qquad N^{\phi}(v,r)=\frac{4Gj(v)}{r^{2}}. (4.13)

The boundary metric is again (4.2) with dΘ2=dϕ2\mathrm{d}\Theta^{2}=\mathrm{d}\phi^{2}. We assume that N(v,r)2N(v,r)^{2} has two positive roots, i.e, L2m(v)2>j(v)2L^{2}m(v)^{2}>j(v)^{2}, and the larger one is

rh(v)=[4L2Gm(v)(1+1j(v)2L2m(v)2)]1/2.\displaystyle r_{h}(v)=\left[4L^{2}Gm(v)\left(1+\sqrt{1-\frac{j(v)^{2}}{L^{2}m(v)^{2}}}\right)\right]^{1/2}. (4.14)

The counterterm for gravity in d=2d=2 [62, 63] is known as

Ict=18πGLd2zγ,\displaystyle I_{\mathrm{ct}}=-\frac{1}{8\pi GL}\int\mathrm{d}^{2}z\sqrt{-\gamma}, (4.15)

with which the mass and angular momentum are calculated as

hv=m(v),pϕ,v=j(v).\displaystyle h_{v}=m(v),\qquad p_{\phi,v}=j(v). (4.16)

In considering the reference Euclidean geometry at vv, we have two parameters βv\beta_{v} and ωvϕ\omega^{\phi}_{v}. The metric that has hvh_{v} and pϕ,vp_{\phi,v}, and is compatible with the boundary metric (2.34), is found to be

ds2=\displaystyle\mathrm{d}s^{2}= N(v,r)2dτ2+dr2N(v,r)2+r2(dϕiN^ϕ(v,r)dτ)2,\displaystyle N(v,r)^{2}\mathrm{d}\tau^{2}+\frac{\mathrm{d}r^{2}}{N(v,r)^{2}}+r^{2}\left(\mathrm{d}\phi-i\hat{N}^{\phi}(v,r)\mathrm{d}\tau\right)^{2}, (4.17)
N^ϕ(v,r):=\displaystyle\hat{N}^{\phi}(v,r):= Nϕ(v,r)Nϕ(v,rh(v)).\displaystyle N^{\phi}(v,r)-N^{\phi}(v,r_{h}(v)). (4.18)

Thus, the parameters are determined as

βv=4πrN(v,rh(v))2,ωvϕ=Nϕ(v,rh(v)).\displaystyle\beta_{v}=\frac{4\pi}{\partial_{r}N(v,r_{h}(v))^{2}},\qquad\omega^{\phi}_{v}=N^{\phi}(v,r_{h}(v)). (4.19)

The coarse-grained entropy is

Sv=π2Grh(v),\displaystyle S_{v}=\frac{\pi}{2G}r_{h}(v), (4.20)

and the first law holds as

S˙v=βv(h˙vωvϕp˙ϕ,v).\displaystyle\dot{S}_{v}=\beta_{v}(\dot{h}_{v}-\omega^{\phi}_{v}\dot{p}_{\phi,v}). (4.21)

The second law means rh(v)rh(0)r_{h}(v)\geq r_{h}(0). Again, this is guaranteed under the null energy condition (4.11) for a null vector =v+(N2/2)rNϕϕ\ell=\partial_{v}+(N^{2}/2)\partial_{r}-N^{\phi}\partial_{\phi}, by the same logic.

4.3 Reissner-Nordström (4d flat)

Refer to caption
Figure 7: Vaidya-Bonner spacetime. The boundary is sent to infinity.

The final example is a 4-dimensional charged black hole. The solution we consider is the Vaidya-Bonner spacetime,

ds2=\displaystyle\mathrm{d}s^{2}= N(v,r)2dv2+2dvdr+r2dΘ2,N(v,r)2=12Gm(v)r+Gq(v)24πr2,\displaystyle-N(v,r)^{2}\mathrm{d}v^{2}+2\mathrm{d}v\mathrm{d}r+r^{2}\mathrm{d}\Theta^{2},\qquad N(v,r)^{2}=1-\frac{2Gm(v)}{r}+\frac{Gq(v)^{2}}{4\pi r^{2}}, (4.22)
A=\displaystyle A= q(v)2πr(dvdrN(v,r)2),\displaystyle-\frac{q(v)}{2\pi r}\left(\mathrm{d}v-\frac{\mathrm{d}r}{N(v,r)^{2}}\right), (4.23)

where dΘ2=dθ2+sin2θdϕ2\mathrm{d}\Theta^{2}=\mathrm{d}\theta^{2}+\sin^{2}\theta\mathrm{d}\phi^{2}, and rh(v)r_{h}(v) is the largest root of N(v,r)2N(v,r)^{2} given as

rh(v)=Gm(v)+(Gm(v))2Gq(v)24π,\displaystyle r_{h}(v)=Gm(v)+\sqrt{(Gm(v))^{2}-\frac{Gq(v)^{2}}{4\pi}}, (4.24)

with the quantity in the square root assumed positive. We choose r=Rr=R as the boundary, whose induced metric is

ds2|bdy=N(v,R)2dv2+R2dΘ2.\displaystyle\mathrm{d}s^{2}|_{\mathrm{bdy}}=-N(v,R)^{2}\mathrm{d}v^{2}+R^{2}\mathrm{d}\Theta^{2}. (4.25)

Since this induced metric is time-dependent, we send RR to infinity, in order for it to be time-independent. This infinity corresponds to the past null infinity \mathscr{I}^{-}. This time, the respected quantities are the mass and the electric charge.101010 We can also consider (3.11) for all components of AμA_{\mu}, which will make the coarse-graining conditions stronger. However, as the spacetime is spherically symmetric now, the entropy and its thermodynamic laws will be the same as what we are going to derive below. Observed at \mathscr{I}^{-}, the mass and charge are

hv=m(v),qv=q(v),\displaystyle h_{v}=m(v),\qquad q_{v}=q(v), (4.26)

where the counterterm is, as usual, taken to be the contribution from the Minkowski spacetime. The mass is the Bondi-Sachs type mass.

The reference Euclidean spacetime is

ds2=N(v,r)2dτ2+dr2N(v,r)2+r2dΘ2,A=q(v)2π(1r1rh(v))dτ,\displaystyle\mathrm{d}s^{2}=N(v,r)^{2}\mathrm{d}\tau^{2}+\frac{\mathrm{d}r^{2}}{N(v,r)^{2}}+r^{2}\mathrm{d}\Theta^{2},\qquad A=-\frac{q(v)}{2\pi}\left(\frac{1}{r}-\frac{1}{r_{h}(v)}\right)\mathrm{d}\tau, (4.27)

with the Lagrange multipliers given by

βt=4πrN(v,rh(v))2,λt=12πAτ|r,\displaystyle\beta_{t}=\frac{4\pi}{\partial_{r}N(v,r_{h}(v))^{2}},\qquad\lambda_{t}=\frac{1}{2\pi}A_{\tau}|_{r\to\infty}, (4.28)

where λt\lambda_{t} is the conjugate of qvq_{v}.111111 Here we take the gauge Aτ|r=rh=0A_{\tau}|_{r=r_{h}}=0 so that λt\lambda_{t} becomes a pure boundary value. Without this gauge, (3.42) does not hold as it is, but the term related to Aτ|r=rhA_{\tau}|_{r=r_{h}} will appear.

The entropy is, hence,

Sv=πGrh(v)2.\displaystyle S_{v}=\frac{\pi}{G}r_{h}(v)^{2}. (4.29)

The first law is satisfied in the form (3.42):

S˙v=βv(h˙vλvq˙v).\displaystyle\dot{S}_{v}=\beta_{v}(\dot{h}_{v}-\lambda_{v}\dot{q}_{v}). (4.30)

The second law again implies (4.11), with =v+(N2/2)r\ell=\partial_{v}+(N^{2}/2)\partial_{r}.

If we decide not to respect the charge, the resulting entropy must become larger. This choice is possible when q(0)=0q(0)=0, necessary for S0S_{0} to coincide with the initial entropy. In this case, the reference Euclidean geometry would be

ds2=f(v,r)dτ2+dr2f(v,r)+r2dΘ2,f(v,r)=12Gm(v)r.\displaystyle\mathrm{d}s^{2}=f(v,r)\mathrm{d}\tau^{2}+\frac{\mathrm{d}r^{2}}{f(v,r)}+r^{2}\mathrm{d}\Theta^{2},\qquad f(v,r)=1-\frac{2Gm(v)}{r}. (4.31)

Then, the horizon radius is 2Gm(v)2Gm(v), which is larger than rh(v)r_{h}(v), and hence the entropy becomes larger as well. This reflects the fact that coarse-grained entropy increases when fewer constraints are imposed.

5 Discussions

We have introduced a coarse-grained entropy that respects asymptotic values on each time slice, as a new measure of black hole thermodynamics. The entropy is the horizon area of a specific reference Euclidean black hole, faithful to the first and second laws. The second law, proven from the CFT, is expected to be a reflection of the thermal (or statistical) nature of quantum gravity. As a matter of fact, the second law is not always satisfied, providing a constraint on gravity. Through several examples, we have discovered that the constraint seems to be related to the integrated null energy condition.

The coarse-grained entropy introduced in this paper has its definition both on the boundary and in the bulk. As a newly established dictionary, this will enable us to study the non-equilibrium process from CFT to gravity, and vice versa. Although we have focused on the thermodynamics of gravity in this paper, it is of course possible to survey the evolution of the coarse-grained entropy in the CFT from the bulk analysis.

While the first law was derived within the Einstein gravity, we asked for the help of the holographic dictionary in showing the second law. It is valuable to further survey what the second law means in more complex situations, since the second law seems to imply something that the Einstein theory itself does not tell. Particularly, we did not specify the origin of the null rays in Vaidya models, but if it were possible, we could add the matter fields to the set of respected quantities. For this purpose, the solution in [70], which offers a null ray collapse model within given fields, will be helpful.

Also, studying the relation between the coarse-grained entropy and the area of the apparent horizon will be helpful to understand the origin of the second law in the bulk. In the Schwarzschild-AdS case in section 4, for example, our coarse-grained entropy coincides with the apparent horizon area when the time foliation is specified by the coordinate vv. However, the definition of our coarse-grained entropy is essentially different from the apparent horizon; while the apparent horizon is foliation-dependent, our entropy depends only on the boundary time foliation. To connect them, a proper rule to extend the boundary foliation to the bulk must be established. A definition of the entropy for dynamical black holes was recently proposed also in [71], where the apparent horizon appears as well. One possible future direction is to find a connection between those notions.

It is also interesting to consider other ways of coarse-graining to provide some thermodynamic constraints. In our coarse-graining process, especially, only the boundary observables on the time slice are referred to at each time tt, and hence it does not tell any difference that happens deep inside the bulk. To be more finely grained, the coarse-grained state must also respect observables off the time slice. By doing so, the information of some codimension-0 domain on the boundary is kept, preserving some bulk region causally related to it. This idea has been addressed in the literature [40, 45, 47].

Finally, one of the aims of non-equilibrium thermodynamics is to describe the relaxation process. If there is non-equilibrium thermodynamics in gravity, does it also have knowledge about the relaxation of gravitational systems? In some sense, we can say that the no-hair theorem teaches what states are possible as the final states, after matter has fallen behind the horizon. However, it is still a mystery how the system evolves to the final state, or whether the fate of the system must be thermodynamically controlled. We expect that studying non-equilibrium black hole thermodynamics will reveal those problems.

5.1 Erratum added

After the publication, I found some ambiguity regarding the dictionary (2.42) and (3.7), which leads some inconsistencies as explained below.121212I thank Sotaro Sugishita for discussion.

From (3.34), we obtain

δBIgrav(E)[Γ]=(h[Γ]ΩApA[Γ])δB,\displaystyle\delta_{B}I_{\mathrm{grav}}^{(E)}[\Gamma]=\left(h[\Gamma]-\Omega^{A}p_{A}[\Gamma]\right)\delta B, (5.1)

where Γ\Gamma was defined to be the collection of the Lagrange multipliers. On the other hand, from the holographic dictionary, Igrav(E)[Γ]=lnZ[Γ]I_{\mathrm{grav}}^{(E)}[\Gamma]=-\ln Z[\Gamma], with

Z[Γ]=Trexp[B(HΩAPAdd1θσJI(θ)OI(θ))],\displaystyle Z[\Gamma]=\mathrm{Tr}\exp\left[-B\left(H-\Omega^{A}P_{A}-\int\mathrm{d}^{d-1}\theta\sqrt{-\sigma}J^{I}(\theta)O_{I}(\theta)\right)\right], (5.2)

we also have

δB(lnZ[Γ])=(Tr(ρ[Γ]H)ΩATr(ρ[Γ]PA)dd1θσJITr(ρt[Γ]OI))δB.\displaystyle\delta_{B}\left(-\ln Z[\Gamma]\right)=\left(\mathrm{Tr}\left(\rho[\Gamma]H\right)-\Omega^{A}\mathrm{Tr}\left(\rho[\Gamma]P_{A}\right)-\int\mathrm{d}^{d-1}\theta\sqrt{-\sigma}J^{I}\mathrm{Tr}\left(\rho_{t}[\Gamma]O_{I}\right)\right)\delta B. (5.3)

Hence, (5.1) and (5.3) are contradictory when they are compared via (3.17) and (3.18). Here, (3.17) and (3.18) come from (3.7).

Therefore, although source insertion in (2.42) seems widely accepted as far as I recognize, we need to include the metric variation to the other source terms if they are not turned off. That means, (3.7) is modified as

Tr(ρtTab)Tr(ρtOI(z))wI(z)σab=limR(RL)d2Yab(z),\displaystyle\mathrm{Tr}\left(\rho_{t}T_{ab}\right)-\mathrm{Tr}\left(\rho_{t}O_{I}(z)\right)w^{I}(z)\sigma_{ab}=\lim_{R\to\infty}\left(\frac{R}{L}\right)^{d-2}Y_{ab}(z), (5.4)

which resolves the above contradiction. In fact, (5.4) can also be confirmed by comparing the computations in the Einstein-scalar theory and in the CFT. One may refer to [72] where such an example is demonstrated.

Unfortunately, however, the story was not so simple. For the U(1)\mathrm{U}(1) gauge field in the bulk, instead of (5.1), we have from (3.40)

δBIgrav(E)=(h[Γ]ΩApA[Γ]dd1θσ~q[Γ](θ)a(θ))δB,\displaystyle\delta_{B}I_{\mathrm{grav}}^{(E)}=\left(h[\Gamma]-\Omega^{A}p_{A}[\Gamma]-\int\mathrm{d}^{d-1}\theta\sqrt{\tilde{\sigma}}q[\Gamma](\theta)a(\theta)\right)\delta B, (5.5)

where recall σ~=σ\sqrt{\tilde{\sigma}}=\sqrt{-\sigma}. This was, as a matter of fact, derived also in [73] under more general situations. Now, (5.5) is consistent with the form in (5.3) without modifying (3.7).

Therefore, I think that the dictionary between the stress tensor and the Brown-York tensor must be studied more carefully.

Finally, let me explain the interpretation of the examples in section 4, in the light of the above discussion. Since we do not know the matter theory that makes up the null ray in the Vaidya spacetime, it is not obvious whether the Brown-York tensor follows (3.7) or (5.4). Thus, it is possible that hth_{t} is not dual to Tr(ρtH)\mathrm{Tr}(\rho_{t}H), while the masses of the Euclidean black holes in section 4 purely measure Tr(ρref,tH)\mathrm{Tr}(\rho_{\mathrm{ref},t}H) because we chose not to respect the matter field, meaning that the coarse-graining is not properly done. However, it becomes correct only while the source is turned off, in other words, during the periods when μ(v)\mu(v) and the other parameters are constant. Thus, by considering turning on and off the source, our coarse-graining is correct and its associated second law must be satisfied while the source is off.

Acknowledgement

I thank Masafumi Fukuma, Osamu Fukushima, Yuji Hirono, Takanori Ishii, Tomohiro Shigemura, Keito Shimizu, Sotaro Sugishita, Ryota Watanabe, and Takuya Yoda for discussion. I am also grateful for daily advice by Koji Hashimoto and Shigeki Sugimoto. My work was supported by Grant-in-Aid for JSPS Fellows No. 22KJ1944.

Appendix A Derivation of (2.37)

We derive (2.37) in this appendix.131313 Based on discussion with Osamu Fukushima, Takanori Ishii, Tomohiro Shigemura, Keito Shimizu, Sotaro Sugishita, and Ryota Watanabe. We denote the r.h.s. as TT, and will show that TT is deformed to the l.h.s. In the following, the CPT operator Θ\Theta plays a key role. The CPT conjugate state Θψ|\bra{\Theta\psi} of a state |ψ\ket{\psi} is defined as

|ψ,Θψ|ψ=(ψ|Θ|ψ),\displaystyle\forall\ket{\psi^{\prime}},\qquad\braket{\Theta\psi}{\psi^{\prime}}=(\braket{\psi^{\prime}}{\Theta}{\psi})^{*}, (A.1)

where the asterisk means the complex conjugation. Note that Θ\Theta is antiunitary.

With Θ\Theta and any orthonormal basis {|n}\{\ket{n}\}, TT is deformed as follows:

T=\displaystyle T= ϕ|U(t)eβH~/2ΘΘV(t)|ϕ=ϕ|U(t)eβH~/2Θ(n|nn|ΘV(t)|ϕ)\displaystyle\braket{\phi^{\prime}}{U(t)e^{-\beta\tilde{H}/2}\Theta^{\dagger}\Theta V(t)}{\phi}=\bra{\phi^{\prime}}U(t)e^{-\beta\tilde{H}/2}\Theta^{\dagger}\left(\sum_{n}\ket{n}\bra{n}\Theta V(t)\ket{\phi}\right)
=\displaystyle= nn|ΘV(t)|ϕϕ|U(t)eβH~/2|nΘ.\displaystyle\sum_{n}\braket{n}{\Theta V(t)}{\phi}^{*}\braket{\phi^{\prime}}{U(t)e^{-\beta\tilde{H}/2}}{n_{\Theta}}. (A.2)

Here, we have set |nΘ:=Θ|n\ket{n_{\Theta}}:=\Theta^{\dagger}\ket{n}. Assuming that the theory without sources is CPT-invariant, i.e, HH in (2.24) is invariant, we obtain

n|ΘV(t)|ϕ=n|V(t)Θ|ϕ=Θϕ|V(t)|n.\displaystyle\braket{n}{\Theta V(t)}{\phi}^{*}=\braket{n}{V(t)^{\dagger}\Theta}{\phi}^{*}=\braket{\Theta\phi}{V(t)}{n}. (A.3)

In addition, being a c-number, this can be seen calculated in a copied Hilbert space, which we call as QFTL. The remaining part is kept in the original Hilbert space, QFTR. Therefore we conclude

T=Θϕ|Lϕ|RVL(t)UR(t)(n|nLeβH~R/2|nΘR).\displaystyle T={}_{\mathrm{L}}\bra{\Theta\phi}\otimes{}_{\mathrm{R}}\bra{\phi^{\prime}}V_{\mathrm{L}}(t)\otimes U_{\mathrm{R}}(t)\left(\sum_{n}\ket{n}_{\mathrm{L}}\otimes e^{-\beta\tilde{H}_{\mathrm{R}}/2}\ket{n_{\Theta}}_{\mathrm{R}}\right). (A.4)

In this expression, Θϕ|\bra{\Theta\phi} and |nΘ\ket{n_{\Theta}} appear, instead of ϕ|\bra{\phi} and |n\ket{n}. However, those differences do not matter. First, ϕ\phi is just a dummy variable in the generating functional (2.39). Second, either |n\ket{n} or |nΘ\ket{n_{\Theta}} can guarantee (2.28), meaning that we could have started from (2.23) with |nR|nΘR\ket{n}_{\mathrm{R}}\to\ket{n_{\Theta}}_{\mathrm{R}} and reach the same results in the main text. Thus, for simplicity, we have ignored those subtleties in (2.37).

References