Irreversibility in quantum field theories with boundaries

# Irreversibility in quantum field theories with boundaries

Horacio Casini, Ignacio Salazar Landea, Gonzalo Torroba

Centro Atómico Bariloche and CONICET
S.C. de Bariloche, Río Negro, R8402AGP, Argentina
###### Abstract

We study conformal field theories with boundaries, and their boundary renormalization group (RG) flows, using methods from quantum information theory. Positivity of the relative entropy, together with unitarity and Lorentz invariance, give rise to bounds that characterize the irreversibility of such flows. This generalizes the recently proved entropic -theorem to higher dimensions. In dimensions with a boundary, we prove the entropic -theorem – the decrease of the two-dimensional Weyl anomaly under boundary RG flows. In higher dimensions, the bound implies that the leading area coefficient of the entanglement entropy induced by the defect decreases along the flow. Our proof unifies these properties, and provides an information-theoretic interpretation in terms of the distinguishability between the short distance and long distance states. Finally, we establish a sum rule for the change in the area term in theories with boundaries, which could have implications for models with localized gravity.

\@maketitle

## 1 Introduction

Boundaries are ubiquitous in broad areas of physics, including high energy, condensed matter physics and, of course, real experiments. In particular, systems that preserve scale invariance play a central role in the dynamics and phases of quantum field theories (QFTs) with boundaries. This is analogous to how fixed points of the renormalization group (RG) organize the dynamics of QFTs without boundaries [1]. In this work we consider fixed points of relativistic QFTs with boundaries, and how they are connected by RG flows. Our goal is to understand if these flows are irreversible. We will approach the problem using quantum information theory.

The irreversibility of relativistic QFTs without boundaries has been established in two [2, 3], three [4], and four [5, 6] spacetime dimensions. Methods from quantum information theory, based on strong subadditivity of the entanglement entropy (EE), have provided a unifying proof for all these irreversibility theorems [6]. A natural question is whether some of these properties survive in the presence of boundaries.

The earliest result in this direction is the -theorem: this establishes the irreversibility of boundary RG flows111A boundary RG flow is an RG flow triggered by a relevant deformation that is turned on only at the boundary. The bulk, as in all cases in this work, is conformal. in 2d CFTs with boundaries [7, 8]. The entropic version of this result was recently obtained in [9]. More generally, one can consider higher-dimensional quantum field theories with boundaries, that respect conformal invariance, and their boundary RG flows.222See [10] for a review with references. In this paper we study generalizations of the approach of [9] to higher-dimensional boundary conformal field theories (BCFTs). We will use methods from quantum information theory to characterize boundary RG flows.

The introduction of boundaries poses new challenges to this approach. One is related to the (partial) breaking of the Poincaré group, which makes relativistic constraints less powerful. Another new issue is that, while the RG flow occurs at the boundary, the dynamics of the boundary degrees of freedom by themselves is nonlocal. The nonlocality is induced by interactions with the bulk, and it could prevent the existence of irreversibility theorems for such flows. A surprising outcome of our results will be that the nonlocality cancels out from the measures of quantum information theory that we analyze (the entanglement entropy and the relative entropy). While locality, causality and conformal invariance in the bulk are ultimately responsible for this, this is a point that remains to be fully understood and exploited.

Before proceeding, let us mention some encouraging results in higher dimensions. One is the -theorem of [11], regarding the irreversibility of boundary RG flows in 3d. This was shown using dilaton methods [5]. There is also evidence for a boundary -theorem in 4d systems with boundaries [12]. On the other hand, holographic models also exhibit irreversibility of boundary RG flows in general [13]. Other examples, and a suggestion for an irreversible quantity, are explored in the recent work [14].

This paper is organized as follows. In Sec. 2 we discuss relevant aspects of boundary RG flows, and how they can be measured in terms of the entanglement entropy. In Sec.3 we study the relative entropy and its connection with the entanglement entropy. We argue that, by taking the limit where the Cauchy surface becomes null, the modular hamiltonian contribution vanishes (in a certain range for the relevant deformation dimension). This is one of our main results, which allows us to use positivity of the relative entropy to bound the change in the EE between the UV and the IR. It extends [9] to higher dimensions; its consequences are discussed in Sec. 4. We prove the entropic version of the -theorem in 2+1 dimensions, and “area theorems” for higher dimensional QFTs with boundaries. In Sec. 5 we derive a sum rule that provides an explicit expression for the change in the EE in terms of the two-point function of the trace of the boundary stress tensor. Finally, Sec. 6 contains our conclusions and future directions.

## 2 Boundary RG flows and entanglement entropy

We consider a quantum field theory (QFT) on a -dimensional spacetime with a boundary at . The first coordinates , while lives on the half-line . The ‘bulk’ fields have some prescribed boundary conditions at , and there may also be degrees of freedom localized at the boundary.

We will be interested in boundary conformal field theories (BCFTs) and the RG flows that connect them.333See e.g. [15] for a discussion of defects in CFTs. The bulk is always taken to be a CFT, and the boundary breaks the bulk conformal group to

 SO(2,d)→SO(2,d−1). (2.1)

In other words, a conformal boundary preserves the dilatation, Poincaré and special conformal transformations that keep the boundary fixed. We start from a UV boundary CFT , and trigger a boundary RG flow by turning on a relevant deformation on the boundary,

 ST1=ST0+∫∂Mdd−1xgO (2.2)

where the dimension of at the UV fixed point is . The theory with relevant deformation is denoted by , and we assume that the flow ends at an IR fixed point – a different BCFT.

We wish to understand if these flows are irreversible in some sense. In order to address this, we will employ methods from quantum information theory, generalizing the recent entropic proof of the -theorem given in [9]. The basic idea is to compute the entanglement entropy for a (half) sphere of radius , centered at the boundary, and then use the monotonicity of the relative entropy in order to characterize the irreversibility of the RG. Before describing this in more detail, let us review the simpler case of the -theorem.

For a 2d theory with boundary at , the contribution to the thermal entropy due to the boundary can only decrease along boundary RG flows [7, 8]. This contribution, denoted in general as , can also be captured by the entanglement entropy (EE) on an interval ,

 S(ρR)=−tr(ρRlogρR), (2.3)

where is the vacuum density matrix reduced to the interval. At a fixed point, this is of the form [16]

 S(ρR)=c6logRϵ+c0+logg, (2.4)

where is the bulk CFT central charge, is a short distance cutoff, and is some regulator-dependent constant. While the regulator dependence does not allow to obtain directly from the EE, we can measure changes in impurity entropy by subtracting the EE between the UV and IR fixed points, where the bulk contributions cancel out:

 S(mR≪1)−S(mR≫1)=loggUVgIR. (2.5)

Here is some characteristic mass scale for the RG flow, for instance, in (2.2). Ref. [9] showed that positivity of the relative entropy implies . Furthermore, monotonicity of this quantity also allows to define a monotonic interpolating -function. Our goal is to extend this entropic -theorem to higher dimensions.

Motivated by these results, we will consider the vacuum EE for (half) spheres

 (x1)2+…+(xd−2)2+y2=R2,y≥0, (2.6)

in boundary field theories that undergo RG flows (2.2). This is shown in Fig. 1, together with the associated causal domain. Unitarity and causality dictate that the EE depends on geometric properties of the boundary of the entangling region, and is the same for any Cauchy surface inside the causal diamond. As in [9, 17], our strategy will be to deform the Cauchy surface towards the light-cone in order to equate the EE with the relative entropy.

In more detail, let be the vacuum reduced density matrix for the UV fixed point theory ,

 σR=tr¯V(|0⟩⟨0|), (2.7)

where is the ground-state of , and is the region outside the half-sphere (2.6). And let be the corresponding reduced density matrix for the theory along the flow (2.2). The two theories are defined microscopically by the same operator content, but evolve with different hamiltonians. The relative entropy

 S(ρR|σR)=trV(ρR(logρR−logσR)) (2.8)

measures the distinguishability between the two states. It is positive and monotonic, namely it increases with the size of the region. Introducing the modular Hamiltonian of the UV BCFT,

 σR=e−HσtrHσ, (2.9)

the relative entropy can be rewritten as

 S(ρR|σR) = trV(ρHσ)−trV(σHσ)+trV(ρRlogρR−σRlogσR) (2.10) = Δ⟨Hσ⟩−ΔS

with and .

Following [9], the strategy will be to consider a Cauchy surface that approaches the light-cone of the causal diamond associated to the entangling region. We will show that, under certain circumstances, in this limit. Then , and so the change in the EE inherits the positivity and monotonicity properties of the relative entropy. Before turning to this, though, let us discuss the general structure of the EE for BCFTs.

### 2.1 Structure of the entanglement entropy

In order to understand the structure of the EE for a half-sphere of radius in a BCFT, let us first recall that the sphere entropy in a CFT without boundary is of the form

 S(R)=μd−2Rd−2+μd−4Rd−4+…+⎧⎨⎩(−)d−224Alog(R/ϵ)deven.(−)d−12Fdodd. (2.11)

For a CFT the coefficients are proportional to inverse powers of the cutoff. The geometric origin of these terms, for spheres or their null deformations, was recently explained in [18]. The area term proportional to comes from a volume integral over the entangling surface, the term proportional to arises from a similar integral containing the intrinsic curvature of the surface, and so on. The nonlocal contribution, the logarithmic term in even , comes from a Wess-Zumino action on the surface.

Moving on to the BCFT case, the entangling surface intersects the boundary on a sphere of radius . As in [18], this will give rise to new divergent terms in the EE, which can be written as integrals of local geometric quantities on this that respect Lorentz invariance. Therefore, for even (i.e. an odd - dimensional boundary), we have

 S(R)=μd−2Rd−2+~μd−3Rd−3+μd−4Rd−4+…+(−1)d−224Alog(R/ϵ)+(−1)d−22~F, (2.12)

where quantities that come from the boundary are denoted with a tilde. The last term behaves like a boundary F-function. Although the logarithmic term implies that will change when choosing different regulators, differences in for boundary RG flows will be physical. The simplest example is the case of the -theorem discussed before, where in (2.4) plays the role of . Examples of this for were studied in [12]. The different coefficients are powers of the cutoff for a BCFT but can contain other divergent or finite contributions depending on the dimensionful scales of the theory as we approach a fixed point.

For odd, the arguments in [18] show that it is possible to have a Wess-Zumino action, now localized on the , consistent with Lorentz invariance. Hence, we expect an entropy of the form

 S(R)=μd−2Rd−2+~μd−3Rd−3+μd−4Rd−4+…+(−1)d−324~Alog(R/ϵ)+(−1)d−12F. (2.13)

The logarithmic term, absent for CFTs in odd dimensions, comes from the Wess-Zumino action on . This should correspond to a Weyl anomaly in the BCFT. To see this, one can follow the arguments of [19] that related the logarithmic term in (2.11) to the Weyl anomaly of the CFT. This involves a conformal mapping between the causal diamond of the spherical region and the static patch of de Sitter. The EE entropy on the sphere then becomes the thermal entropy of de Sitter. Rotating to euclidean time, this entropy is obtained from the logarithm of the partition function on . The logarithmic term in the entropy is then seen to be determined directly by the Weyl anomaly

 ⟨Tμμ⟩=2(−1)d/2AEd+…, (2.14)

where is the Euler density,444The integral of on a -dimensional sphere gives . and we are not showing terms that vanish for a sphere. Applying this to BCFTs with (2.13) maps the causal diamond of (2.6) to half of a causal patch in de Sitter. In the euclidean version, this is a hemisphere of , with boundary at the equator. The logarithmic term in the thermal entropy should then come from a Weyl anomaly localized at the equator,

 ⟨Tμμ⟩=2(−1)d−12~AEd−1δ(θ−π2)+…, (2.15)

where is the azimuthal angle of the hemisphere. See also [20, 14, 21].

The case of a 3d theory with boundary will play an important role below. There is a “Graham-Witten” anomaly from the 2d boundary [22]

 ⟨Tμμ⟩=−b24πR(2)δ(θ−π2)+… (2.16)

Then we have

 ~A=b12. (2.17)

The -anomaly decreases along boundary RG flows [11]. Our goal will be to establish this, together with natural extensions to higher dimensions (area-theorems) using entropic methods.

## 3 Analysis of the relative entropy

In this section we analyze the relative entropy (2.10) between theories and . After determining the structure of the modular Hamiltonian and properties of the stress tensor, we take the limit where the Cauchy surface approaches the light-cone of the causal domain.

### 3.1 Modular Hamiltonian and stress-tensor

For a CFT and a spherical region, the modular Hamiltonian is local, and given by

 H=∫ΣdσημξνTμν, (3.1)

where is a unit vector normal to the Cauchy surface , and is the Killing vector for a conformal transformation that keeps the sphere fixed,

 ξν=πR(R2−(x0)2−(xi)2,−2x0xi). (3.2)

Recall that this can be proved by a conformal map between the causal domain of the sphere and the Rindler wedge [19, 23]; here the density matrix is thermal and the modular flow corresponds to boosts.

This approach can be applied to BCFTs as well to establish that (3.1) is valid for a half-sphere centered at the boundary. The conformal transformation of [19] in signature reads

 xμ=2Xμ+(X⋅X)Cμ1+2(X⋅C)+(X⋅X)(C⋅C)−R2Cμ. (3.3)

We can perform it along the direction transverse to the boundary, . The boundary is then uniformly accelerated inside the Rindler wedge, with acceleration . This respects the boost symmetry, and so the density matrix is again thermal with respect to the boost Hamiltonian. Therefore, the modular Hamiltonian is still given by the Rindler hamiltonian, and transforming back to the sphere gives (3.1).555See also [24] for this result; this work also discusses a map with tangential to the boundary. We stress that here is the stress-tensor operator of the BCFT , because we are using the conformal map between the sphere and the Rindler wedge.

Having established the form of the modular hamiltonian for the BCFT theory , let us analyze the contribution to the relative entropy, for a given Cauchy surface . Since is a conserved current, is a conserved charge on . This means that, as an operator, is independent of the choice of Cauchy surface inside the given domain of dependence. All the dependence on has to come from the state. The key point, discussed in [9, 17], is that the state brings in explicit dependence on . The reason is that, in order to compare both states and , we have to map the algebra of operators of to that of , , at a given Cauchy surface . In the Heisenberg representation, these algebras evolve with different Hamiltonians, which means that the map depends on ,

 UΣ~ϕλ(x)U†Σ=ϕλ(x),x∈Σ. (3.4)

The corresponding map for the state, , will then also depend on the choice of the Cauchy surface.

To make this more explicit, let us consider the change in the stress-tensor expectation value between the two states,

 Δ⟨Tμν(x)⟩=tr(ρRTμν(x))−tr(σRTμν(x)). (3.5)

Let us split the energy-momentum tensor into a bulk part, that is continuous at the boundary, and localized terms,

 Tμν(x)=Tbulkμν(x)+δ(y)tμν(x), (3.6)

with associated to the presence of the boundary. First, away from the boundary and evolve with the same Hamiltonian –they only differ by the relevant perturbation (2.2) localized at . This means that is independent on the choice of Cauchy surface. This is a local quantity, which can only depend on (there is translation invariance on the other coordinates), and on the metric and the normal . We parametrize

 Δ⟨Tbulkμν(x)⟩=α1(y)(gμν−NμNν)+α2(y)NμNν. (3.7)

Imposing the conservation condition away from the defect sets , so that is constant. But since the expectation value has to decay away from the defect, we have . Also, is the BCFT stress tensor and so it is traceless, ; this sets . Therefore,

 Δ⟨Tbulkμν(x)⟩=0, (3.8)

and we are only left with a possible contribution proportional to . This is a key property of codimension one defects; as discussed in the Appendix, it is not valid for other codimensions.

On the boundary , the state depends on the choice of Cauchy surface . We can write the stress-tensor expectation value in terms of local geometric objects , and curvatures:

 Δ⟨Tμν(x)⟩=δ(y)(α1(gμν−NμNν)+α2ημην+α3NμNν+α4(Nμην+Nνημ))+O(ϵ2K2,…), (3.9)

where we indicated that contributions from curvatures are suppressed by the short-distance cutoff [17]. In conformal perturbation theory around the UV fixed point we expect, on dimensional grounds,

 αi∼g2ϵ(d−1)−2ΔO (3.10)

for . For smaller operator dimension, a nonperturbative finite term of order can appear.

Plugging (3.9) into (3.1), and taking into account that at , we find that the integral is restricted to the component of the entangling surface at the boundary, and that all surviving terms are proportional to the contraction :

 Δ⟨Hσ⟩=∫Σ∩∂Mdσαηαξα+O(ϵ2K2). (3.11)

Here denotes a linear combination of the in (3.9); in what follows we will only need their order of magnitude (3.10). Note that becomes an integral of a local quantity on the boundary.

We conclude that the dependence of on is given simply by the flux of the conformal Killing vector through . This is analogous to what happened for bulk RG flows in [17], although here the bulk part is conformal, and all the contributions are restricted to the boundary.

### 3.2 The light-cone limit

The last step is to take a Cauchy surface that approaches the null boundary of the causal domain of dependence. The modular flow keeps the boundary of the causal domain fixed, so becomes null there. As a result,

 ξαηα∣∣Σnull→0. (3.12)

When this is valid, , and hence the change in the EE equals the relative entropy between the two states,

 S(ρR|σR)=−ΔS. (3.13)

However, the null limit has to be taken in a controlled way, because (3.12) multiplies a possibly divergent contribution (3.10). Following [17], we can approach the null limit with a Cauchy surface forming a hyperboloid of radius , obtaining

 Δ⟨Hσ⟩∼(g2ϵd−1−2ΔO)(a2Rd−3). (3.14)

The first factor here comes from the expansion (3.10), while is the contribution of for a Cauchy surface that is a hyperboloid of curvature scale of order . This scaling with should be contrasted with the dependence on a constant time slice, which always dominates over the area term associated to . The null limit reduces the scaling with by two powers, replacing this by .

Given a short distance cutoff , we can have at most , so (3.14) becomes

 Δ⟨Hσ⟩∼g2Rd−3ϵd+1−2ΔO. (3.15)

In this way, in the null limit for relevant perturbations in the range

 ΔO

For a -dimensional boundary, the perturbation is relevant if . For , (3.16) then imposes no restrictions (it covers the whole range of relevant deformations), while for higher it does. In this regime, (3.13) applies. For , the modular Hamiltonian contribution to the relative entropy will compete with the area term from the EE, and we cannot fix the sign of .

## 4 Consequences in different dimensions

We now restrict to (3.16), and analyze the consequences of (3.13).

First, in the EEs of the two states can be written as (see (2.12))

 S(σR) = c6logRϵ+c′0+loggUV S(ρR) = c6logRϵ+c′0+logg(R), (4.1)

where we have defined a running -function . In the IR, , the impurity entropy for the IR BCFT. The relative entropy then reads

 S(ρR|σR)=−logg(R)gUV. (4.2)

Monotonicity of the relative entropy requires . This is the entropic -theorem proved in [9].

In , corresponding to a 2d boundary, (2.12) gives the fixed point expression

 S(σR)=μ1Rϵ−F+bUV3logRϵ (4.3)

where we used , the boundary central charge in (2.16). The entropy has a more complicated radial dependence along the RG flow, but near the IR fixed point,

 S(ρR)=μ1Rϵ−F+bIR3log(mR)−bUV3log(mϵ). (4.4)

Here is a typical mass scale for the RG flow. Then the relative entropy for becomes

 S(ρR|σR)≈13(bUV−bIR)log(mR). (4.5)

Positivity of the relative entropy thus implies

 bUV≥bIR, (4.6)

so that the boundary -anomaly decreases along boundary RG flow. The result (4.6) was proved using dilaton methods in [11].

Lastly, in higher dimensions and for long distances , we have

 ΔS=(~μIRd−3−~μUVd−3)Rd−3+… (4.7)

In the range (3.16), Eq. (3.13) together with positivity of the relative entropy imply that

 ~μUVd−3≥~μIRd−3. (4.8)

This is a new result on the irreversibility of boundary RG flows, established with quantum information theory methods. It implies that the leading area term associated to the boundary can only decrease along boundary RG flows. From [25], the flow in the area term in relativistic QFTs is related to the renormalization of Newton’s constant. Thus, we expect that (4.8) may be relevant in theories of localized gravity [26].

In summary, we have found that, in QFTs with boundaries, positivity of the relative entropy on a null Cauchy surface leads to an inequality that unifies the -theorem, the -theorem and an area theorem in higher dimensions. The change in these quantities thus acquires a precise information-theoretic meaning as a meassure of distinguishability between the reduced states of the UV and IR fixed points in boundary RG flows.

## 5 A sum rule in boundary QFTs

For RG flows in relativistic QFTs without boundary, the change in the area term

 Δμd−2=μIRd−2−μUVd−2 (5.1)

between the UV and IR fixed points obeys a sum rule (see [27, 28, 25, 29]),

 Δμd−2=−πd(d−1)(d−2)∫ddxx2⟨Θ(x)Θ(0)⟩, (5.2)

with . Given that for boundary RG flows we just established that

 Δ~μd−3=~μIRd−3−~μUVd−3≤0, (5.3)

it is natural to ask whether a similar sum rule exists for this quantity.

In order to address this, let us follow [25], and consider a half-sphere of radius . The entangling region approaches a Rindler wedge along (say) . The boundary preserves boosts along , so the modular Hamiltonian is given by

 H=−2π∫w1≥0,wd−1≥0dd−1→ww1T00(w), (5.4)

for any boundary QFT (not necessarily at a fixed point). Note that here is the Rindler edge, while is the position of the boundary.

Taking , is dominated by the IR fixed point, and hence

 ΔS=S(ρR)−S(σR)≈(~μIRd−3−~μUVd−3)Rd−3; (5.5)

all the other terms are subleading. The change can be obtained by performing a small variation of ,

 RdΔSdR=(d−3)ΔS. (5.6)

Recall that under a small change of state , the first law allows to relate the variation in the entropy to the change in the modular hamiltonian, . We can view (5.6) as a small change in the state due to a dilatation. By applying the first law to the variation of entropy in each of the states we get

 ΔS = 1d−3∫ddx⟨Θρ(x)Hρ⟩−1d−3∫ddx⟨Θσ(x)Hσ⟩ (5.7) = 1d−3∫ddxδ(xd−1)⟨θρ(x)Hρ⟩,

where we used that the trace of the stress-tensor vanishes in the BCFT state , and introduced

 Θρ(x)=δ(xd−1)θρ(x),θρ=tr(tρ). (5.8)

Also, and are the modular Hamiltonians of the theories and respectively. Intuitively, is implementing a global scale transformation, and its only nonvanishing contribution comes from the scale dependence at the boundary.

Replacing (5.4) into (5.7),

 ΔS=−2πd−3∫ddxδ(xd−1)∫dd−1ww1⟨θ(x)T00(w)⟩, (5.9)

where the stress-tensors are evaluated in the theory , but we have eliminated the subindices to streamline the formulas.

As we just reviewed, the factor of is implementing a global scale transformation. The quantity is traceless and conserved in for bulk points, and therefore we can use the same arguments as in Sec. 3.1 to conclude it has to vanish for in the bulk. In fact this is the change of the bulk stress tensor expectation value under dilatations, and the bulk stress tensor expectation value vanishes identically. Only the boundary component of the stress-tensor will contribute to (5.9), so here we can replace . In this way, we end up with a purely boundary expression

 ΔS=−2πd−3∫ddxδ(xd−1)∫w1≥0dd−1wδ(wd−1)w1⟨θ(x)t00(w)⟩. (5.10)

Translation invariance along the boundary implies that the integrand is independent of the spatial coordinates transverse to the Rindler edge . These -integrals give simply a factor of the boundary area , and we arrive to

 Δ~μd−3=−2πd−3∫∂Mdd−1x∫w1>0dw1w1⟨θ(x)t00(w)⟩. (5.11)

Finally, we need to relate this to an integral of the two-point function of . Performing a diffeomorphism transformation tangent to the boundary (i.e. , with ) implies the conservation equation [8]

 ∂αtαβ=Tβ,d−1bulk. (5.12)

We argued in Sec. 3.1 that in vacuum expectation value. Since in (5.11) we have a global scale transformation on the vacuum expectation value , we can here use inside the correlation function. In more detail,

 ∂∂wα∫∂Mdd−1x⟨θ(x)tαβ(w)⟩=R∂∂R⟨∂∂wαtαβ(w)⟩=R∂∂R⟨Tβ,d−1bulk(w)⟩=0. (5.13)

This means that

 ∫∂Mdd−1x⟨θ(x)tαβ(w)⟩=∫∂Mdd−1x(gαβ∇2−∂α∂β)F(s), (5.14)

with , and we used that this correlator is symmetric in . In particular, taking the trace over the coordinates of the boundary gives

 ⟨θ(x)θ(w)⟩=(d−2)∇2F(s). (5.15)

Let us now perform the integral,

 ∫∂Mdd−1x∫w1>0dw1w1⟨θ(x)t00(w)⟩ = ∫∂Mdd−1x∫w1>0dw1w1→∇2F(s) (5.16) = ∫∂Mdd−1xF(x),

where we integrated by parts twice. On the other hand,

 ∫∂Mdd−1xx2⟨θ(x)θ(0)⟩ = (d−2)∫∂Mdd−1xx2∇2F(x) (5.17) = (d−1)(d−2)∫∂Mdd−1xF(x),

after integration by parts.

Comparing both expressions, we derive the sum rule

 Δ~μd−3=−2π(d−1)(d−2)(d−3)∫∂Mdd−1xx2⟨θ(x)θ(0)⟩. (5.18)

This sum rule agrees with the sign deduced from positivity of the relative entropy.

The case can be obtained by taking into account that in this case the behavior is logarithmic instead of a power law, or alternatively from the equation above by adimensionalizing the integral using a mass scale

 ΔS=−(mR)d−32π(d−1)(d−2)(d−3)∫∂Mdd−1xm−(d−3)x2⟨θ(x)θ(0)⟩, (5.19)

and expanding for we have

 ΔS=−πlog(mR)∫∂Md2xx2⟨θ(x)θ(0)⟩. (5.20)

Comparing with (2.17), (4.3), and (4.4), this gives

 bUV−bIR=3π∫∂Md2xx2⟨θ(x)θ(0)⟩. (5.21)

This coincides with the sum rule obtained in [11] by a different method.

Note that this applies to , so that it does not provide a sum rule for the change in in . The reason is that for , the contribution to does not scale like a power of , and so it cannot be deduced from the Rindler limit. We can obtain the change in by integrating , which makes it clear that the modular Hamiltonian over all scales would be needed.

## 6 Conclusions and future directions

In this paper we proved a series of theorems about the monotonicity of the first subleading term in the entanglement entropy of CFTs with planar boundaries. We accomplished this by considering a null Cauchy surface, and equating the change in the EE with the relative entropy between the UV and IR states. This provides a unified picture for the -theorem, -theorem, and boundary area theorems in higher dimensions. Several questions arise on possible generalizations of our results. We will summarize our main future directions in what follows.

One direction is to consider defects with higher codimension. Our approach here does not apply to these cases because, as we review in the Appendix, the bulk contribution (3.7) to the modular Hamiltonian does not vanish. It is not clear yet how to resolve this issue. In this context, the recent work [14] conjectured that the defect free energy may provide a monotonic quantity for defect RG flows. For codimension one defects this agrees with the change in the EE entropy, but for higher codimensions both quantities are different. It would be interesting to find an information-theoretic version of this conjecture, and analyze its validity.

We saw that for , the first subleading term in the EE of a BCFT corresponds to a logarithmic divergence. This is a universal term, in the sense that it is independent on the way we choose to regularize the theory, and has a nice interpretation in terms of the boundary Weyl anomaly . For the universal terms are even more subleading and we do not have access to their monotonicity properties using only the relative entropy. For CFTs without boundaries, one can access the universal terms in using the strong subaditivity of the EE and the Markovian property of the vacuum on the null cone. This gives the -theorem and the entropic -theorem [4, 6]. It would be interesting to explore if some of these results can be extended to BCFTs.

Finally, in Sec. 5 we have shown that the change of the subleading terms of the EE under an RG flow is related to a sum rule for the two point function of the trace of the boundary stress-energy tensor. In QFTs with full Poincare invariance, Adler and Zee [30, 31] showed that this sum rule for the two point function is related to the renormalization of the Newton’s constant. It would be interesting to understand if there are similar implications for the graviton effective action in field theories with boundaries.

## Acknowledgments

We thank Y. Sato for discussions. This work was partially supported by CONICET (PIP grant 11220150100299), CNEA, and Universidad Nacional de Cuyo, Argentina. H.C. acknowledges an “It From Qubit” grant of the Simons Foundation. G.T. is also supported by ANPCYT PICT grant 2015-1224.

## Appendix A Comments on higher-codimension defects

In this Appendix we briefly discuss some properties of the stress-tensor in theories with higher codimension defects. Most of this is standard and well-known; we include it to highlight how the approach in the main text fails when the codimension is not one.

We have a bulk CFT in -dimensional flat spacetime, coupled to a defect of codimension . We turn on a relevant deformation localized at the boundary; as before, denotes the density matrix of the UV fixed point, while is the density matrix for the theory along the flow.

We split the coordinates as

 xμ=(xα,ya),μ=0,…,d−1;α=0,…,d−n−1;a=d−n,…,d−1. (A.1)

The defect is placed at

 ya=0. (A.2)

In the main part of this work we focused on . Our goal is to evaluate , namely the change in expectation value of the defect CFT stress tensor between the states and . This will generalize the discussion in Sec. 3.1 to other codimensions.

As discussed before, is independent of the choice of Cauchy surface. Operatorially, is conserved and traceless away from the defect. Furthermore, by rotational invariance it can only depend on . The conservation condition requires

 Δ⟨Taα⟩=0. (A.3)

The remaining nonzero components can then be parametrized as

 Δ⟨Tαβ⟩ = h(y)ηαβ (A.4) Δ⟨Tab⟩ = f1(y)δab+f2(y)(yayby2−δabn).

Here we used Poincaré invariance along the defect to constrain , and rotational invariance in the transverse directions to fix .

Requiring vanishing trace relates

 f1(y)=−d−nnh(y). (A.5)

The conservation condition gives, on the other hand,

 n−1n(f′2(y)+nyf2(y))+f′1(y)=0. (A.6)

So in general we have an arbitrary free function, which we may take as .

The form of is further determined if the vacuum state is conformally invariant: since there is no dimensionful coupling, and the stress tensor has dimension , . As a result,

 Δ⟨Tαβ⟩ = h0ydηαβ (A.7) Δ⟨Tab⟩ = −d−nnh0ydδab+dn−1h0yd(yayby2−δabn) Δ⟨Taα⟩ = 0.

Only the constant is arbitrary, and this depends on the type of conformal defect.

For , the contribution in (A.4) vanishes identically, and then setting in (A.6) requires to be a constant. This constant has to vanish because the one-point function decays away from the defect. Then has to vanish as well, in order to satisfy (A.5). This recovers of the main text. On the other hand, we see that for higher codimensions is non-vanishing. This means that the modular Hamiltonian contribution to the relative entropy will be generically nonzero and larger than . The approach of Secs. 3, 4 does not apply here.

## References

You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters