On the shape dependence of Entanglement Entropy

We study the shape dependence of entanglement entropy (EE) by deforming symmetric entangling surfaces. We show that entangling surfaces with a rotational or translational symmetry extremize (locally) the EE with respect to shape deformations that break some of the symmetry (i.e. the 1st order correction vanishes). This result applies to EE and Renyi entropy for any QFT in any dimension. Using Solodukhin’s formula in 4d and holography in any d, we calculate the 2nd order correction to the universal EE for CFTs and simple symmetric entangling surfaces. For several entangling surfaces we find that the 2nd order correction is positive for any perturbation, and thus the corresponding symmetric entangling surface is a local minimum. Some of the results are extended to free massive fields and to 4d Renyi entropy.


Introduction
Entanglement entropy is a measure of the quantum correlations of a system. It has a very wide range of applications from condensed matter physics to quantum field theory  and AdS/CFT . The EE in QFT is generally hard to calculate, and most computations have been done in simple setups, such as: free fields, CFTs, and symmetric entangling surfaces (e.g. spheres and planes). Thus there is a need to obtain analytical results for interacting theories, non-CFTs, and less symmetrical entangling surfaces. This work aims to make a step in this direction by studying the shape dependence of EE. Previous works on the shape dependence of EE include [63][64][65][66][67][68][69][70][71][72][73][74][75][76][77][78][79][80][81].
The divergent structure of entanglement entropy for a CFT in d-dimensions is: where R is the scale of the entangling region, and δ is the UV cutoff. The leading divergence is the area law, and all of the power law divergences are non-universal. We will be interested in the universal term S (univ) , which in d = even is the coefficient of the log divergence and in d = odd it is the finite term.
JHEP12(2015)043 Figure 1. Illustration of a perturbed circle r(φ) = 1 + n a n cos(nφ). Left: a perturbation without a zero mode a 0 = 0. Right: a perturbation with a zero mode a 0 = 0.3. Consider a QFT parametrized by coordinates (t, y i ,r), where i = 1 . . . , d − 2. For instance take r to be the radial coordinate in spherical coordinates, and y i to be angles parameterizing the entangling surface. We will always work in a constant time slice t = 0. Consider a codimension-2 entangling surface defined by: where r 0 (y i ) is some given function of y i . Thus we choose r to be the dependent coordinate, and y i as the independent coordinates. We denote the entanglement entropy corresponding to the entangling surface r 0 (y i ) as S 0 . Now we slightly perturb the entangling surface: where is a small parameter, and f (y i ) is some arbitrary perturbation function. 1 The entanglement entropy will change as a result of the perturbation of the entangling surface, and it can generally be written as an expansion in : The above procedure was carried out in [63,64] for the case of a perturbed sphere for a CFT in d dimensions. They start with a sphere entangling surface in a flat space-time background, and perturb the sphere as follows: where R is the sphere radius, the a's are constants, and the Y 's are (real) hyper-spherical harmonics. Then they calculate the resulting change in the universal term of the EE. They find that the the change in the universal EE vanishes at 1 st order, namely S (univ) 1 = 0 (1.6) 1 As an example, a perturbed circle in d = 3 is shown in figure 1. In this case, (1.3) is given by: n an cos(nφ)], where we Fourier expanded the perturbation f . R is the radius of the circle, φ is the angle in polar coordinates, and an are the Fourier coefficients.

JHEP12(2015)043
Thus, for a CFT the sphere is a local 2 extremum with respect to perturbations of the entangling surface. Additionally, for holographic CFTs [64] uses the Ryu-Takayanagi formula [35,36] (more specifically, its generalization to higher derivative gravity [83][84][85][86]) to calculate the 2 nd order correction S (univ) 2 where C T is the (positive) central charge appearing in the 2-point function: could have depended on the three parameters C t , t 2 , and t 4 of the 3-point function T T T , but it turns out that it depends just on C T . S (univ) 2 is clearly positive, implying that the sphere is a local minimum for holographic CFTs. (1.7) was compared to S (univ) 2 obtained from Solodukhin's formula in 4d CFTs, and precise agreement was found.
In this work we will generalize the above results of [64] to less symmetric entangling surfaces, and (in some cases) to non-CFTs. We will show that entangling surfaces with a rotational or translational symmetry in some direction (see figure 2), extremize the universal EE with respect to shape deformations that break some of the symmetry, 3 i.e.: The proof of this result will be purely geometrical, and hence will apply to any QFT (and also to Renyi entropy). A simple corollary is that for d = even, (1.8) is true also for multiply connected entangling surfaces. 4 Additionally, we will calculate the 2 nd order correction S (univ) 2 (whose sign determines if the local extremum is a minimum or a maximum) for some simple entangling surfaces. We perform this calculation using holography (the Ryu-Takayanagi formula), and also by using Solodukhin's formula (for 4d CFTs). For several examples of symmetric entangling surfaces, we find that the 2 nd order correction is positive, and this corresponds to a local minimum. We conjecture this to hold more generally to symmetric entangling surfaces. We also comment on results for free massive fields, and 4d Renyi entropy. 2 If the topology of the entangling surface is allowed to change, then the EE can become unbounded from below, as shown in 4d in [82]. 3 It should be emphasized that the family of shape perturbations that extremize the EE, are only those which break some of the symmetry of the original entangling surface. For example, for a 3d circle the perturbation in figure 1-left is allowed, whereas the perturbation in figure 1-right is not allowed because it has a component in the radial direction (a zero mode of the Fourier expansion). Such perturbations cannot deform a sphere into a larger/smaller sphere. More generally, such perturbations cannot deform a surface of revolution into a different surface of revolution (containing the same number of symmetries). 4 For d = even the universal term of the EE is a log divergence which is determined locally by the shape of the entangling surface. Thus for a multiply connected entangling surface the log term is a superposition of the contribution from each separate piece of the entangling surface. r(φ, y i ) = r 0 (y i ) + A n a n (y i ) cos(nφ) , n = 0 (2.4)

JHEP12(2015)043
where r 0 (y i ) doesn't depend on the symmetry direction φ, the a n (y i ) are functions of y i , and A n are constants. The resulting EE can be expanded as: Now lets consider the same perturbation but with negative sign, i.e. → − : The resulting EE will be (just flipping the sign of in (2.5): 5 A general perturbation (without a zero mode) can be written as: Then the EE can be expanded: At linear order in the modes don't mix, and their contributions to S1 add up linearly: where U1(n) and U2(n) are some functions of n. Because the modes don't mix at linear order, we can compute the contribution of a single mode and then sum over all modes.

JHEP12(2015)043
But the two perturbations (2.6) and (2.4) describe precisely the same entangling surface, only rotated. This can be seen by performing φ → φ + π n on (2.6), which gives (2.4). Since the two entangling surfaces are the same, they have same EE and therefore from (2.5) and (2.7) we have:S and we proved what we wanted. An important point to make is that the proof above used only the rotation symmetry, and did not use any specific property of entanglement entropy or of the QFT. Therefore stationarity will hold for any quantity which is a function of the surface (e.g. Renyi entropy), and for any QFT.

The second order correction in holography
Consider a boundary QFT parametrized by coordinates (t, y i , r), i = 1 . . . , d − 2. z is the holographic coordinate, and we set t = 0. The holographic EE according to the Ryu-Takayanagi formula is: The corresponding equation of motion for the bulk surface is: where there is a summation convention on y i . Consider an entangling surface defined by: Now perturb the entangling surface: where is a small parameter. Then the bulk surface will also get perturbed: 6 For the rest of this section we assume d = even dimensions, and we will want to compute the universal log term. We can then write the resulting perturbed L and S as:

JHEP12(2015)043
We can derive a formula for the 2 nd order correction S 2 . Assuming that the entangling surface has a symmetry in all directions (e.g. the cylinder surface in the next section), terms containing r 2 will fall and we get:

Cylinder entangling surface
Now let us assume that the entangling surface is a cylinder S p × R d−2−p in flat spacetime. We denote y j as cartesian directions along R d−2−p (where −L ≤ y j ≤ L) and Ω p as directions along the sphere S p (of radius R). The metric in cylindrical coordinates is: where α(z) and β(z) are some functions of z. The holographic EE written in these cylindrical coordinates is (see also (A.20)): where: For an AdS metric we have: F 2 (z) = 1, and where Y lp (Ω p ) are (real) hyperspherical harmonics, the n j are integers, the a {n j ,lp} are the coefficients of the Fourier expansion, and for conciseness we defined {n j } ≡ n 1 , . . . , n d−2−p and {l p } ≡ l, m 1 . . . , m p−1 . The bulk surface perturbation r 1 of (3.5) can also be Fourier expanded (see (A.25)): where r (1) {n j ,lp} (z) are functions of z, which must obey the boundary condition: r

JHEP12(2015)043
Now we can derive an expression for the boundary term (the last term on the r.h.s. of (3.7)): where L is the length of the R d−2−p directions, and we take the limit L → ∞. This formula is a boundary term which is evaluated at the boundary z = δ. We will want to extract the universal log divergence from this formula in d = even dimensions.
In principle, r 0 (z) and r The terms in (3.13) can then be written as: where we extracted the log terms. Using (3.13), and also the fact that the background is where we used u 0 = 1 and q 0 = R, and multiplied by (−1) d−2 2 , see (1.1). It can be shown that for the two special cases p = 0, 1, the boundary term above is the only contribution to S (univ) 2 . Thus the p = 1 case gives: We see from (3.18) that the sign of S

Plane entangling surface
The plane entangling surface R d−2 is a special case of the cylinder with p = 0. Therefore (3.17) becomes: 9 S (univ) 2 So the sign of S (univ) 2 depends solely on the coefficientũ d , which we shall now determine. For a plane entangling surface and Einstein gravity in the bulk, we will find an explicit solution 10 (A.19): where N is a normalization constant, andñ 2 ≡ i n 2 i . The small z expansion of this function gives (see (A. 16 Since the above expression is positive, we have proved that a plane is a local minimum for Einstein gravity in the bulk. It is now natural to conjecture that (3.21) holds for any CFT. (3.21) was derived for d = even, but it would be simple task to obtain the analogous d = odd result. The d = odd result will differ from (3.21) only by its d dependence, and for d = 3 it was obtained in [77].

(3.22)
This precisely matches (4.13) which we will obtain in the next section for 4d CFTs via Solodukhin's formula.

4d CFT: Solodukhin's formula
For a d = even CFT there is a universal log term (1.1): and as in (1.4), we can expand in small : 9 Note that (3.17), (3.18), and (3.19) apply to any asymptotically AdS background. 10 This result was also derived in [52] in the context of "entanglement density".

JHEP12(2015)043
Solodukhin's formula [65] (which applies to a CFT in 4d) is: 11 where E 2 is the Euler density, and the integrals are over the entangling surface Σ, and The second fundamental form and extrinsic curvature are defined as: where γ µν is the induced metric on the entangling surface. We consider shape perturbations that leave the topology of the entangling surface fixed, therefore the change in the Euler term above is zero, and S (univ) 2 of (4.2) will be given by the integral of I 2 : In (4.4) I 2 is always positive [78,82] and is zero for the sphere. Therefore the sphere locally minimizes I 2 and S (univ) . This argument also works for flat entangling surfaces (plane, strips) since these have I 2 = 0.
In the following, we calculate S (univ) 2 for several examples. See also [63,64].
Example 1: plane entangling surface. The metric in cartesian coordinates is: Consider a plane at x = 0, where y 2 , y 3 are coordinates along the surface. Now perturb it as follows: x = 0 + f (y 2 , y 3 ) (4.8) The vector normal to the surface is: The second fundamental form at order O( 2 ) is: k a µν starts at order since the unperturbed plane is flat. Plugging this in (4.4) gives: 11 a4 and c4 are the a and c anomalies in 4d. Both are normalized such that for a real scalar field their value is 1. We will sometimes use CT = c 4 3π 4 instead of c4.

JHEP12(2015)043
Note that I 2 starts at order 2 , therefore the order correction vanishes even before integration over the surface. Now Fourier expand the perturbation: f (y 2 , y 3 ) = ∞ n 2 ,n 3 =1 a n 2 n 3 cos(n 2 y 2 ) cos(n 3 y 3 ) (4.12) and plug (4.11), (4.12) in (4.6): where we used c 4 = 3π 4 C T , and the integral L −L dy cos 2 (ny) = L, where the width of the plane is very large L → ∞. So we got a positive result, and therefore the universal term of a plane in 4d is a local minimum. We see that (4.13) precisely matches the result (3.22) obtained in holography.
Example 3: cylinder entangling surface. The metric in cylindrical coordinates is: Consider a perturbed cylinder entangling surface with radius R: We get:
Example 5: 4d waveguide surface. Consider a general waveguide entangling surface (e.g. figure 2-right) given by: r(y, φ) = r 0 (φ) (4.25) where we use cylindrical coordinates (r, y, φ). This surface has translational symmetry in the y direction. Now we perturb the surface as follows: where the G i (r 0 , f 2 ) are some functions of r 0 , f 2 and their derivatives. Explicit calculation gives the m 4 coefficient: may still turn out to be positive).

Renyi entropy
There is a generalization of (4.3) to Renyi entropy in a 4d a CFT [70,79]: where S (univ) q is the universal log term for the q-th Renyi entropy. Q 1,2 (q) are functions of q such that Q 1 (1) = Q 2 (1) = 1, so it matches (4.3) in the EE limit q → 1. Up to these functions, (4.29) and (4.3) are the same, hence the results in the previous sections can be used. In particular, if Q 2 (q) is positive then for any entangling surface the Renyi entropy will have the same sign as the EE. Note that Q 2 (q) is positive for free fields, and there is strong evidence that it is positive for holographic CFTs [70,82,89]. If this turns out to be correct, then the Renyi entropy will be a local minimum whenever the EE is.

Free massive fields
Certain universal EE terms for free massive field theories have been found [16][17][18][19]. The so called "universal area law" for free scalars or fermions with mass m has the following form: where γ d is a positive constant that depends only on the dimension d. All of the shape dependence is contained in A Σ , the area of the entangling surface. Table 1 lists the sign factors (−1) # in (4.30) for a Dirac fermion, a conformally coupled scalar, and a minimally coupled scalar.
As an example, let us now compute the "universal area law" term for a deformed sphere entangling surface: The area of the deformed sphere is: where we plugged (4.31) and performed the integrals. The first term on the r.h.s. is the area of the undeformed sphere, the O( ) correction vanishes as expected, and the O( 2 ) correction is positive. Therefore the O( 2 ) correction in S has the same sign as the zeroth order, which can be read from table 1. Generalizing to non-spheres, it is easy to see that the area of a perturbed surface of revolution is larger than that of the unperturbed surface of revolution. It would be interesting to perform a similar analysis to higher curvature terms for free massive fields (i.e to curvature corrections to the "universal area law").
A similar analysis applies to "universal area law" terms in interacting theories [19,21,69], with a shape dependence that comes only from the area A Σ . The EE for a CFT perturbed by a relevent operator of dimension ∆ = d+2 2 , contains the following term: where λ is the coupling constant. Such log terms occur both in odd and even dimensions.

Discussion
In this work we studied the shape dependence of entanglement entropy by deforming symmetric entangling surfaces. We showed that entangling surfaces with a rotational or translational symmetry locally extremize the EE with respect to shape deformations that break some of the symmetry. This result applies to EE and Renyi entropy for any QFT in any dimension. Using Solodukhin's formula and holography, we calculated the 2nd order correction to the EE for CFTs and simple symmetric entangling surfaces. In several cases we found that the 2nd order correction is positive, and thus the corresponding symmetric entangling surface is a local minimum. Perhaps this result can be shown to hold more generally for other symmetric entangling surfaces.
• The calculation in section 3 considered only Einstein gravity in the bulk, and it would be interesting to consider also higher derivative gravity. For spheres, [64] found that S (univ) 2 depends just on C T and not on t 2 or t 4 (these are the three parameters in the 3-point function of stress tensors). It would be nice to check if this continues to hold for other entangling surfaces.
• Computing the FG coefficientsũ d andq d in (3.17) by solving the bulk EOMs for the cylinder (A.23) and (A.26). Maybe it is also possible to show, using a more general principle, thatq d andũ d in (3.18) must have a definite sign. It would also be interesting to consider (in holography) more general entangling surfaces with a symmetry.
• The work of [69] attempted to compute S univ 2 for a plane entangling surface in a d = 4 CFT using the perturbative formalism of EE [66]. They were not able to obtain the I 2 term of Solodukhin's formula (4.3), and the current situation is somewhat puzzling. Our result (3.21) as well as (1.7) (obtained in [64]) might help in resolving this puzzle.
• It would be interesting to repeat the analysis of section 4.1 for 6d CFTs using the results of [80,81].
• Other possible extensions are to compute higher orders S (univ) j for j > 2, and also to perform computations in a curved space-time background.