# Local Modular Hamiltonians from the Quantum Null Energy Condition

###### Abstract

The vacuum modular Hamiltonian K of the Rindler wedge in any relativistic quantum field theory is given by the boost generator. Here we investigate the modular Hamiltoninan for more general half-spaces which are bounded by an arbitrary smooth cut of a null plane. We derive a formula for the second derivative of the modular Hamiltonian with respect to the coordinates of the cut which schematically reads K^{\prime\prime}=T_{vv}. This formula can be integrated twice to obtain a simple expression for the modular Hamiltonian. The result naturally generalizes the standard expression for the Rindler modular Hamiltonian to this larger class of regions. Our primary assumptions are the quantum null energy condition — an inequality between the second derivative of the von Neumann entropy of a region and the stress tensor — and its saturation in the vacuum for these regions. We discuss the validity of these assumptions in free theories and holographic theories to all orders in 1/N.

## I Introduction and Summary

The reduced density operator \rho for a region in quantum field theory encodes all of the information about observables localized to that region. Given any \rho, one can define the modular Hamiltonian K by

\displaystyle\rho=e^{-K}. | (1) |

Knowledge of this operator is equivalent to knowledge of \rho, but the modular Hamiltonian frequently appears in calculations involving entanglement entropy. In general, i.e. for arbitrary states reduced to arbitrary regions, K is a complicated non-local operator. However, in certain cases it is known to simplify.

The most basic example where K simplifies is the vacuum state of a QFT in Rindler space, i.e. the half-space t=0,x\geq 0. The Bisognano–Wichmann theorem Bisognano:1976za () states that in this case the modular Hamiltonian is

\displaystyle\Delta K=\frac{2\pi}{\hbar}\int d^{d-2}y\int_{0}^{\infty}x\,T_{tt% }\,dx | (2) |

where \Delta K\equiv K-\braket{K}_{\rm vac} defines the vacuum-subtracted modular Hamiltonian, and y are d-2 coordinates parametrizing the transverse directions. The vacuum subtraction generally removes regulator-dependent UV-divergences in K. Other cases where the modular Hamiltonian is known to simplify to an integral of local operators are obtained via conformal transformation of Eq. (2), including spherical regions in CFTs Casini:2011kv (), regions in a thermal state of 1+1 CFTs Cardy:2016fqc (), and null slabs Bousso:2014sda (); Bousso:2014aa ().

Using conservation of the energy-momentum tensor, one can easily re-express the Rindler modular Hamiltonian in Eq. (2) as an integral over the future Rindler horizon u\equiv t-x=0 which bounds the future of the Rindler wedge:

\displaystyle\Delta K=\frac{2\pi}{\hbar}\int d^{d-2}y\int_{0}^{\infty}v\,T_{vv% }\,dv, | (3) |

where v\equiv t+x. It is important to note that standard derivations of (2) or (3), e.g. Bisognano:1976za (); Casini:2011kv (), do not apply when the entangling surface is defined by a non-constant cut of the Rindler horizon (see Fig. 1). One of the primary goals of this paper is to provide such a derivation.

For a large class of quantum field theories satisfying a precise condition specified momentarily, we will show that the vacuum modular Hamiltonian for the region \mathcal{R}[V(y)] above an arbitrary cut v=V(y) of a null plane is given by

\displaystyle\Delta K=\frac{2\pi}{\hbar}\int d^{d-2}y\int_{V(y)}^{\infty}(v-V(% y))\,T_{vv}\,dv | (4) |

This equation has been previously derived by Wall for free field theories Wall:2011hj () building on Burkardt:1995ct (); Sewell:1982zz (), and to linear order in the deformation away from V(y)=\rm const in general QFTs by Faulkner et al. Faulkner:2016mzt (). In CFTs, conformal transformations of Eq. (4) yield versions of the modular Hamiltonian for non-constant cuts of the causal diamond of a sphere.

The condition leading to Eq. (4) is that the theory should satisfy the quantum null energy condition (QNEC) Bousso:2015mna (); Bousso:2015wca (); Koeller:2015qmn (); Akers:2016ugt () — an inequality between the stress tensor and the von Neumann entropy of a region — and saturate the QNEC in the vacuum for regions defined by cuts of a null plane. We will review the statement of the QNEC in Sec. II.

The QNEC has been proven for free and superrenormalizable Bousso:2015wca (), as well as holographic Koeller:2015qmn (); Akers:2016ugt () quantum field theories. We take this as reasonable evidence that the QNEC is a true fact about relativistic quantum field theories in general, and for the purposes of this paper take it as an assumption. In Sec. II we will show how saturation of the QNEC in a given state leads to an operator equality relating certain derivatives of the modular Hamiltonian of that state to the energy-momentum tensor. Applied to the case outlined above, this operator equality will be integrated to give Eq. (4).

Given the argument in Sec. II, the only remaining question is whether the QNEC is in fact saturated in the vacuum state for entangling surfaces which are cuts of a null plane. This has been shown for free theories in Bousso:2015wca (). In Sec. III, we prove that this is the case for holographic theories to all orders in 1/N. We emphasize that Eq. (4) holds purely as a consequence of the validity of the QNEC and the saturation in the vacuum for \mathcal{R}, two facts which are potentially true in quantum field theories much more generally than free and holographic theories.

Finally, in Sec. IV we will conclude with a discussion of possible extensions to curved backgrounds and more general regions, connections between the relative entropy and the QNEC, and relations to other work.

## II Main Argument

### II.1 Review of QNEC

The von Neumann entropy of a region in quantum field theory can be regarded as a functional of the entangling surface. We will primarily be interested in regions to one side of a cut of a null plane in flat space, for which the entangling surface can be specified by a function V(y) which indicates the v-coordinate of the cut as a function of the transverse coordinates, collectively denoted y. See Fig. 1 for the basic setup. Each cut V(y) defines a half-space, namely the region to one side of the cut. We will pick the side towards the future of the null plane. For the purposes of this section we are free to consider the more general situation where the entangling surface is only locally given by a cut of a null plane. Thus the von Neumann entropy can be considered as a functional of a profile V(y) which defines the shape of the entangling surface, at least locally.

Suppose we define a one-parameter family of cuts V(y;\lambda)\equiv V(y;0)+\lambda\dot{V}(y), with \dot{V}(y)>0 to ensure that \mathcal{R}(\lambda_{1})\subset\mathcal{R}(\lambda_{2}) if \lambda_{1}>\lambda_{2}. If S(\lambda) is the entropy of region \mathcal{R}(\lambda), then the QNEC in integrated form states that

\int d^{d-2}y\,\braket{T_{vv}(y)}\dot{V}(y)^{2}\geq\frac{\hbar}{2\pi}\frac{d^{% 2}S}{d\lambda^{2}}. | (5) |

In general there would be a \sqrt{h} induced metric factor weighting the integral, but here and in the rest of the paper we will assume that the y coordinates have been chosen such that \sqrt{h}=1.

By taking advantage of the arbitrariness of \dot{V}(y) we can derive from this the local form of the QNEC. If we take a limit where \dot{V}(y^{\prime})^{2}\to\delta(y-y^{\prime}), then the l.h.s. reduces to \braket{T_{vv}}. We define S^{\prime\prime}(y) as the limit of d^{2}S/d\lambda^{2} in the same situation:

\frac{d^{2}S}{d\lambda^{2}}\to S^{\prime\prime}(y)~{}~{}~{}~{}\text{when}~{}~{% }~{}~{}\dot{V}(y^{\prime})^{2}\to\delta(y-y^{\prime}). | (6) |

Taking the limit of the nonlocal QNEC then gives the local one:

\braket{T_{vv}}\geq\frac{\hbar}{2\pi}S^{\prime\prime}. | (7) |

The local QNEC together with strong subadditivity can likewise be used to go backward and derive the nonlocal QNEC Bousso:2015mna (); Bousso:2015wca (); Koeller:2015qmn (). The details of that argument are not important here. In the next section we will discuss the consequences of the saturation of the QNEC, and will have to distinguish whether we mean saturation of the nonlocal inequality Eq. (5) or the local inequality Eq. (7), the latter condition being weaker.

### II.2 The QNEC under state perturbations

In this section we consider how the QNEC behaves under small deformations of the state. We begin with a reference state \sigma and consider the deformed state \rho=\sigma+\delta\rho, with \delta\rho traceless but otherwise arbitrary.

Consider a one-parameter family of regions \mathcal{R}(\lambda) as in the previous section. Define \overline{\mathcal{R}}(\lambda) to be the complement of \mathcal{R}(\lambda) within a Cauchy surface. The reduced density operator for any given region \mathcal{R}(\lambda) given by

\rho(\lambda)=\sigma(\lambda)+\delta\rho(\lambda)={\rm Tr}_{\overline{\mathcal% {R}}(\lambda)}\sigma+{\rm Tr}_{\overline{\mathcal{R}}(\lambda)}\delta\rho. | (8) |

By the First Law of entanglement entropy, the entropy of \rho(\lambda) is given by

S(\rho(\lambda))=S(\sigma(\lambda))-{\rm Tr}_{\mathcal{R}(\lambda)}\delta\rho(% \lambda)\log\sigma(\lambda)+o(\delta\rho^{2}). | (9) |

The second term can be written in a more useful way be defining the modular Hamiltonian K_{\sigma}(\lambda) as

K_{\sigma}(\lambda)\equiv-\mathds{1}_{\overline{\mathcal{R}}(\lambda)}\otimes% \log\sigma(\lambda). | (10) |

Defining K_{\sigma}(\lambda) this way makes it a global operator, which makes taking derivatives with respect to \lambda formally simpler. Using this definition, we can write Eq. (9) as

S(\rho(\lambda))=S(\sigma(\lambda))+{\rm Tr}\,\delta\rho K_{\sigma}(\lambda)+o% (\delta\rho^{2}). | (11) |

Now in the second term the trace is over the global Hilbert space, and the \lambda-dependence has been isolated to the operator K_{\sigma}(\lambda). Taking two derivatives, and simplifying the notation slightly, we find

\frac{d^{2}S}{d\lambda^{2}}(\rho)=\frac{d^{2}S}{d\lambda^{2}}(\sigma)+{\rm Tr}% \,\delta\rho\frac{d^{2}K_{\sigma}}{d\lambda^{2}}+o(\delta\rho^{2}). | (12) |

Suppose that the nonlocal QNEC, Eq. (5), is saturated in the state \sigma for all profiles \dot{V}(y). Then, using Eq. (12), the nonlocal QNEC for the state \rho can be written as

\displaystyle\int d^{d-2}y\,\left({\rm Tr}\,\delta\rho\,T_{vv}\right)\dot{V}^{% 2}\geq\frac{\hbar}{2\pi}{\rm Tr}\,\delta\rho\frac{d^{2}K_{\sigma}}{d\lambda^{2% }}+o(\delta\rho^{2}). | (13) |

The operator \delta\rho was arbitrary, and in particular could be replaced by -\delta\rho. Then the only way that Eq. 13 can hold is if we have the operator equality

\frac{d^{2}K_{\sigma}}{d\lambda^{2}}=C+\frac{2\pi}{\hbar}\int d^{d-2}y\,T_{vv}% \dot{V}^{2}. | (14) |

Here C is a number that we cannot fix using this method that is present because of the tracelessness of \delta\rho.

Eq. (14) can be integrated to derive the full modular Hamiltonian K_{\sigma} if we have appropriate boundary conditions. Up until now we have only made use of local properties of the entangling surface, but in order to provide boundary conditions for the integration of Eq. (14) we will assume that the entangling surface is globally given by a cut of a null plane, and that V(y;\lambda=0)=0. We will also make \sigma the vacuum state. In that situation it is known that the QNEC is saturated for free theories, and in the next section we will show that this is also true for holographic theories at all orders in the large-N expansion.

Our first boundary condition is at \lambda=\infty.^{1}^{1}1It is not always possible to consider the \lambda\to\infty limit of a null perturbation to an entangling surface because parts of the entangling surface may become timelike related to each other at some finite value of \lambda, at which point the surface is no longer the boundary of a region on a Cauchy surface. However, when the entangling surface is globally equal to a cut of a null plane this is not an issue. Since we expect that K_{\sigma}(\lambda) should have a finite expectation value in any state as \lambda\to\infty, it must be that dK_{\sigma}/d\lambda\to 0 as \lambda\to\infty. Then integrating Eq. (14) gives

\frac{dK_{\sigma}}{d\lambda}=-\frac{2\pi}{\hbar}\int d^{d-2}y\int_{V(y;\lambda% )}^{\infty}dv\,T_{vv}\dot{V}. | (15) |

Note that this equation implies that the vacuum expectation value \langle K_{\sigma}(\lambda)\rangle_{\rm vac} is actually \lambda-independent, which makes vacuum subtraction easy.

Our second boundary condition is Eq. (3), valid at \lambda=0 when V(y;\lambda)=0. Integrating once morenand making use of this boundary condition, we find

\Delta K_{\sigma}(\lambda)=\frac{2\pi}{\hbar}\int d^{d-2}y\int_{V(y;\lambda)}^% {\infty}(v-V(y;\lambda))\,T_{vv}\,dv | (16) |

which is Eq. (4). Note that the l.h.s. of this equation is now the vacuum-subtracted modular Hamiltonian.

Before moving on, we will briefly comment on the situation where the local QNEC, Eq. (7), is saturated but the nonlocal QNEC, Eq. (5), is not. Then, analogously to S^{\prime\prime} in Eq. (6), one may define a local second derivative of K_{\sigma}:

\frac{d^{2}K_{\sigma}}{d\lambda^{2}}\to K^{\prime\prime}_{\sigma}(y)~{}~{}~{}~% {}\text{when}~{}~{}~{}~{}\dot{V}(y^{\prime})^{2}\to\delta(y-y^{\prime}). | (17) |

Very similar manipulations then show that saturation of the local QNEC implies the equality

K_{\sigma}^{\prime\prime}=\frac{2\pi}{\hbar}T_{vv}. | (18) |

This equation is weaker than Eq. (14), which is meant to be true for arbitrary profiles of \dot{V}(y), but it may have a greater regime of validity. We will comment on this further in Sec. IV.

## III Holographic Calculation

In the previous section we argued that the form of the modular Hamiltonian could be deduced from saturation of the QNEC. In this section we will use the holographic entanglement entropy formula Ryu:2006bv (); Ryu:2006ef (); Hubeny:2007xt (); Faulkner:2013ana () to show that the QNEC is saturated in vacuum for entangling surfaces defined by arbitrary cuts v=V(y) of the null plane u=0 in holographic theories. Our argument applies to any holographic theory defined by a relevant deformation to a holographic CFT, and will be at all orders in the large-N expansion. To reach arbitrary order in 1/N we will assume that the all-orders prescription for von Neumann entropy is given by the quantum extremal surface proposal of Engelhardt and Wall Engelhardt:2014gca (). This is the same context in which the holographic proof of the QNEC was extended to all orders in 1/N Akers:2016ugt ().^{2}^{2}2It is crucial that we demonstrate saturation beyond leading order in large-N. The argument in the previous section used exact saturation, and an error that is naïvely subleading when evaluated in certain states may become very large in others.

As before, the entangling surface in the field theory is given by the set of points \partial\mathcal{R}=\{(u,v,y):v=V(y),u=0\} with null coordinates u=t-x and v=t+x, and the region \mathcal{R} is chosen to lie in the u<0 portion of spacetime. Here y represents d-2 transverse coordinates. The bulk quantum extremal surface anchored to this entangling surface is parameterized by the functions \bar{V}(y,z) and \bar{U}(y,z). It was shown in Koeller:2015qmn (); Akers:2016ugt () that if we let the profile V(y) depend on a deformation parameter \lambda, then the second derivative of the entropy is given by

\frac{d^{2}S}{d\lambda^{2}}=-\frac{d}{4G\hbar}\int d^{d-2}y\,\frac{d\bar{U}_{(% d)}}{d\lambda}~{}, | (19) |

to all orders in 1/N, where \bar{U}_{(d)}(y) is the coefficient of z^{d} in the small-z expansion of \bar{U}(z,y). We will show that \bar{U}=0 identically for any profile V(y), which then implies that d^{2}S/d\lambda^{2}=0, which is the statement of QNEC saturation in the vacuum.

One way to show that \bar{U} vanishes is to demonstrate that \bar{U}=0 solves the quantum extremal surface equations of motion in the bulk geometry dual to the vacuum state of the boundary theory. The quantum extremal surface is defined by having the sum of the area plus the bulk entropy on one side be stationary with respect to first-order variations of its position. One can show that \bar{U}=0 is a solution to the equations of motion if any only if

\frac{\delta S_{\rm bulk}}{\delta\bar{V}(y,z)}=0 | (20) |

in the vacuum everywhere along the extremal surface. This would follow from null quantization if the bulk fields were free Bousso:2015wca (), but that would only allow us to prove the result at order-one in the 1/N expansion.

For an all-orders argument, we opt for a more indirect approach using subregion duality, or entanglement wedge reonstruction Czech:2012bh (); Headrick:2014cta (); Dong:2016eik (); Harlow:2016vwg ().^{3}^{3}3The entanglement wedge of a boundary region is the set of bulk points which are spacelike- or null-related to that region’s quantum extremal surface on the same side of the quantum extremal surface as the boundary region itself. A version of this argument first appeared in Akers:2016ugt (), and we elaborate on it here.

Entanglement wedge reconstruction requires two important consistency conditions in the form of constraints on the bulk geometry which must hold at all orders in 1/N: The first constraint, entanglement wedge nesting (EWN), states that if one boundary region is contained inside the domain of dependence of another, then the quantum extremal surface associated to the first boundary region must be contained within the entanglement wedge of the second boundary region Czech:2012bh (); Wall:2012uf (). The second constraint, \mathcal{C}\subseteq\mathcal{E}, demands that the causal wedge of a boundary region be contained inside the entanglement wedge of that region Czech:2012bh (); Headrick:2014cta (); Wall:2012uf (); Engelhardt:2014gca (); Hubeny:2012wa (). Equivalently, it says that no part of the quantum extremal surface of a given boundary region can be timelike-related to the (boundary) domain of dependence of that boundary region. It was shown in Akers:2016ugt () that \mathcal{C}\subseteq\mathcal{E} follows from EWN, and EWN itself is simply the statement that a boundary region should contain all of the information about any of its subregions. We will now explain the consequences of these two constraints for \bar{U}(y,z).

Without loss of generality, suppose the region \mathcal{R} is defined by a coordinate profile which is positive, V(y)>0. Consider a second region \mathcal{R}_{0} which has an entangling surface at v=u=0 and whose domain of dependence (i.e., Rindler space) contains \mathcal{R}. The quantum extremal surface associated to \mathcal{R}_{0} is given by \bar{U}_{0}=\bar{V}_{0}=0. This essentially follows from symmetry.^{4}^{4}4One might worry that the quantum extremal surface equations display spontaneous symmetry breaking in the vacuum, but this can be ruled out using \mathcal{C}\subseteq\mathcal{E} with an argument similar to the one we present here. The entanglement wedge of \mathcal{R}_{0} is then a bulk extension of the boundary Rindler space, namely the set of bulk points satisfying u\leq 0 and v\geq 0. Then EWN implies that \bar{U}\leq 0 and \bar{V}\geq 0.

The only additional constraint we need from \mathcal{C}\subseteq\mathcal{E} is the requirement that the quantum extremal surface for \mathcal{R} not be in the past of the domain of dependence of \mathcal{R}. From the definition of \mathcal{R}, it is clear that a bulk point is in the past of the domain of dependence of \mathcal{R} if and only if it is in the past of the region u<0 on the boundary, which is the same as the region u<0 in the bulk. Therefore it must be that \bar{U}\geq 0. Combined with the constraint from EWN above, we then conclude that the only possibility is \bar{U}=0. This completes the proof that the QNEC is saturated to all orders in 1/N.

## IV Discussion

We conclude by discussing the generality of our analysis, some implications and future directions, and connections with previous work.

### IV.1 Generalizations and Future Directions

#### General Killing horizons

Though we restricted to cuts of Rindler horizons in flat space for simplicity, all of our results continue to hold for cuts of bifurcate Killing horizons for QFTs defined in arbitrary spacetimes, assuming the QNEC is true and saturated in the vacuum in this context. In particular, Eq. (4) holds with v a coordinate along the horizon. For holographic theories, entanglement wedge nesting (EWN) and the entanglement wedge being outside of the causal wedge (\mathcal{C}\subseteq\mathcal{E}) continue to prove saturation of the QNEC. To see this, note that a Killing horizon on the boundary implies a corresponding Killing horizon in the bulk. Now take the reference region \mathcal{R}_{0} satisfying V(y)=U(y)=0 to be the boundary bifurcation surface. By symmetry, the associated quantum extremal surface lies on the bifurcation surface of the bulk Killing horizon. Then the quantum extremal surface of the region \mathcal{R} defined by V(y)\geq 0 must lie in the entanglement wedge of \mathcal{R}_{0} — inside the bulk horizon — by entanglement wedge nesting, but must also lie on or outside of the bulk horizon by \mathcal{C}\subseteq\mathcal{E}. Thus it lies on the bulk horizon, \bar{U}=0, and the QNEC remains saturated by Eq. (19).

#### Future work

In this work, we have only established the form of K_{\mathcal{R}} for regions \mathcal{R} bounded by arbitrary cuts of a null plane. A natural next direction would be to understand if and how we can extend Eq. (18) to more general entangling surfaces. As discussed above, the QNEC was shown to hold for locally flat entangling surfaces in holographic, free and super-renormalizable field theories Bousso:2015wca (); Koeller:2015qmn (); Akers:2016ugt (). Thus, if we could prove saturation, i.e. that S_{\rm vac}^{\prime\prime}=0 at all orders in 1/N, then we would establish (18) for all regions with a locally flat boundary.

One technique to probe this question is to perturb the entangling surface away from a flat cut and compute the contributions to the QNEC order-by-order in a perturbation parameter \epsilon. Preliminary calculations futurework () have revealed that for holographic theories at leading order in large N, S_{\rm vac}^{\prime\prime}=0 at all orders in \epsilon.

Another interesting problem is to show that in a general QFT vacuum, null derivatives of entanglement entropy across arbitrary cuts of null planes vanish. That, along with a general proof of QNEC will establish (18) as a consequence. We will leave this to future work.

### IV.2 The QNEC as S(\rho\|\sigma)^{\prime\prime}\geq 0

There is a connection between the QNEC and relative entropy, first pointed out in Akers:2016ugt (), that we elaborate on here. The relative entropy S(\rho\|\sigma) between two states \rho and \sigma is defined as

\displaystyle S(\rho\|\sigma)=\operatorname{Tr}{\rho\log{\rho}}-\operatorname{% Tr}{\rho\log{\sigma}} | (21) |

and provides a measure of distinguishability between the two states Nielsen:2011:QCQ:1972505 (). Substituting the definition of K, Eq. (1), into Eq. (21) provides a useful alternate presentation:

\displaystyle S(\rho\|\sigma)=\langle K_{\sigma}\rangle_{\rho}-S(\rho). | (22) |

If Eq. (4) is valid, then taking two derivatives with respect to a deformation parameter, as in the main text, shows that the nonlocal QNEC, Eq. (5), is equivalent to

\displaystyle\partial_{\lambda}^{2}S(\rho(\lambda)\|\sigma(\lambda))\geq 0. | (23) |

For comparison, monotonicity of relative entropy for the types of regions and deformations we have been discussing can be written as

\displaystyle\partial_{\lambda}S(\rho(\lambda)\|\sigma(\lambda))\leq 0. | (24) |

Eq. (23) is a sort of “convexity” of relative entropy.^{5}^{5}5This is distinct from the well-known convexity of relative entropy, which says that S(t\rho_{1}+(1-t)\rho_{2}\|\sigma)\leq tS(\rho_{1}\|\sigma)+(1-t)S(\rho_{2}\|\sigma). Unlike monotonicity of relative entropy, which says that the first derivative is non-positive, there is no general information-theoretic reason for the second derivative to be non-negative. In the event that Eq. (18) holds but not Eq. (4), we would still have

S(\rho\|\sigma)^{\prime\prime}\geq 0. | (25) |

where the {}^{\prime\prime} notation denotes a local deformation as in Sec. II.

It would be extremely interesting to characterize what about quantum field theory and null planes makes (23) true. We can model the null deformation as a non-unitary time evolution in the space of states, with the vacuum state serving as an equilibrium state for this evolution. Then an arbitrary finite-energy state will relax toward the equilibrium state, with the relative entropy S(\rho\|\sigma) characterizing the free energy as a function of time. Monotonicity of relative entropy is then nothing more than the statement that free energy decreases, i.e. the second law of thermodynamics. The second derivative statement gives more information about the approach to equilibrium. If that approach is of the form of exponential decay, then all successive derivatives would alternate in sign. However, for null deformations in quantum field theory we do not expect to have a general bound on the behavior of derivatives of the energy-momentum tensor, meaning that the third derivative of the free energy should not have a definite sign.^{6}^{6}6We thank Aron Wall for a discussion of this point. Perhaps there is some way of characterizing the approach to equilibrium we have here, which is in some sense smoother than the most general possibility but not so constrained as to force exponential behavior.

### IV.3 Relation to previous work

Faulkner, Leigh, Parrikar and Wang Faulkner:2016mzt () have discussed results very similar to the ones presented here. They demonstrated that for first-order null deformations \delta V(y) to a flat cut of a null plane, the perturbation to the modular Hamiltonian takes the form

\braket{K_{\mathcal{R}}}_{\psi}-\braket{K_{\mathcal{R}_{0}}}_{\psi}=-\frac{2% \pi}{\hbar}\int d^{d-2}y\int_{V(y)}dv\,T_{vv}(y)\,\delta V(y) | (26) |

This is precisely the form expected from our equation (4). Faulkner et al. went on to suggest that the natural generalization of the modular Hamiltonian to finite deformations away from a flat cut takes the form of Eq. (4). In the context of holography they showed that this conclusion applied both on the boundary and in the bulk is consistent with JLMS Jafferis:2015del (). In the present paper, we have shown that Eq. (4) holds for theories which obey the QNEC, and for which the QNEC is saturated in the vacuum. A non-perturbative, field theoretic proof of these assumptions remains a primary goal of future work.

###### Acknowledgements.

We would like to thank Chris Akers, Raphael Bousso, Aron Wall and Zach Fisher for discussions and correspondence. This work was supported in part by the Berkeley Center for Theoretical Physics, by the National Science Foundation (award numbers 1521446, and 1316783), by FQXi, and by the US Department of Energy under contract DE-AC02-05CH11231.## References

- (1) J. J. Bisognano and E. H. Wichmann, “On the Duality Condition for Quantum Fields,” J. Math. Phys. 17 (1976) 303–321.
- (2) H. Casini, M. Huerta, and R. C. Myers, “Towards a Derivation of Holographic Entanglement Entropy,” JHEP 05 (2011) 036, arXiv:1102.0440 [hep-th].
- (3) J. Cardy and E. Tonni, “Entanglement hamiltonians in two-dimensional conformal field theory,” J. Stat. Mech. 1612 (2016) no. 12, 123103, arXiv:1608.01283 [cond-mat.stat-mech].
- (4) R. Bousso, H. Casini, Z. Fisher, and J. Maldacena, “Proof of a Quantum Bousso Bound,” Phys. Rev. D90 (2014) no. 4, 044002, arXiv:1404.5635 [hep-th].
- (5) R. Bousso, H. Casini, Z. Fisher, and J. Maldacena, “Entropy on a null surface for interacting quantum field theories and the Bousso bound,” 1406.4545. https://arxiv.org/abs/1406.4545.
- (6) A. C. Wall, “A proof of the generalized second law for rapidly changing fields and arbitrary horizon slices,” Phys. Rev. D85 (2012) 104049, arXiv:1105.3445 [gr-qc]. [Erratum: Phys. Rev.D87,no.6,069904(2013)].
- (7) M. Burkardt, “Light front quantization,” Adv. Nucl. Phys. 23 (1996) 1–74, arXiv:hep-ph/9505259 [hep-ph].
- (8) G. L. Sewell, “Quantum fields on manifolds: PCT and gravitationally induced thermal states,” Annals Phys. 141 (1982) 201–224.
- (9) T. Faulkner, R. G. Leigh, O. Parrikar, and H. Wang, “Modular Hamiltonians for Deformed Half-Spaces and the Averaged Null Energy Condition,” arXiv:1605.08072 [hep-th].
- (10) R. Bousso, Z. Fisher, S. Leichenauer, and A. C. Wall, “A Quantum Focussing Conjecture,” arXiv:1506.02669 [hep-th].
- (11) R. Bousso, Z. Fisher, J. Koeller, S. Leichenauer, and A. C. Wall, “Proof of the Quantum Null Energy Condition,” arXiv:1509.02542 [hep-th].
- (12) J. Koeller and S. Leichenauer, “Holographic Proof of the Quantum Null Energy Condition,” arXiv:1512.06109 [hep-th].
- (13) C. Akers, J. Koeller, S. Leichenauer, and A. Levine, “Geometric Constraints from Subregion Duality Beyond the Classical Regime,” arXiv:1610.08968 [hep-th].
- (14) S. Ryu and T. Takayanagi, “Holographic derivation of entanglement entropy from AdS/CFT,” Phys. Rev. Lett. 96 (2006) 181602, arXiv:hep-th/0603001 [hep-th].
- (15) S. Ryu and T. Takayanagi, “Aspects of Holographic Entanglement Entropy,” JHEP 08 (2006) 045, arXiv:hep-th/0605073 [hep-th].
- (16) V. E. Hubeny, M. Rangamani, and T. Takayanagi, “A Covariant holographic entanglement entropy proposal,” JHEP 07 (2007) 062, arXiv:0705.0016 [hep-th].
- (17) T. Faulkner, A. Lewkowycz, and J. Maldacena, “Quantum Corrections to Holographic Entanglement Entropy,” JHEP 11 (2013) 074, arXiv:1307.2892 [hep-th].
- (18) N. Engelhardt and A. C. Wall, “Quantum Extremal Surfaces: Holographic Entanglement Entropy Beyond the Classical Regime,” JHEP 01 (2015) 073, arXiv:1408.3203 [hep-th].
- (19) B. Czech, J. L. Karczmarek, F. Nogueira, and M. Van Raamsdonk, “The Gravity Dual of a Density Matrix,” Class. Quant. Grav. 29 (2012) 155009, arXiv:1204.1330 [hep-th].
- (20) M. Headrick, V. E. Hubeny, A. Lawrence, and M. Rangamani, “Causality & Holographic Entanglement Entropy,” JHEP 12 (2014) 162, arXiv:1408.6300 [hep-th].
- (21) X. Dong, D. Harlow, and A. C. Wall, “Reconstruction of Bulk Operators within the Entanglement Wedge in Gauge-Gravity Duality,” Phys. Rev. Lett. 117 (2016) no. 2, 021601, arXiv:1601.05416 [hep-th].
- (22) D. Harlow, “The Ryu-Takayanagi Formula from Quantum Error Correction,” arXiv:1607.03901 [hep-th].
- (23) A. C. Wall, “Maximin Surfaces, and the Strong Subadditivity of the Covariant Holographic Entanglement Entropy,” Class. Quant. Grav. 31 (2014) no. 22, 225007, arXiv:1211.3494 [hep-th].
- (24) V. E. Hubeny and M. Rangamani, “Causal Holographic Information,” JHEP 06 (2012) 114, arXiv:1204.1698 [hep-th].
- (25) J. Koeller, S. Leichenauer, A. Levine, and A. Shahbazi Moghaddam, “In Progress,”.
- (26) M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information: 10th Anniversary Edition. Cambridge University Press, New York, NY, USA, 10th ed., 2011.
- (27) D. L. Jafferis, A. Lewkowycz, J. Maldacena, and S. J. Suh, “Relative entropy equals bulk relative entropy,” JHEP 06 (2016) 004, arXiv:1512.06431 [hep-th].