Exact results for corner contributions to the
[2.5mm] entanglement entropy and Rényi entropies
[2.5mm] of free bosons and fermions in 3d
Henriette Elvang and Marios Hadjiantonis
Randall Laboratory of Physics, Department of Physics,
University of Michigan, Ann Arbor, MI 48109, USA
email@example.com and firstname.lastname@example.org
In the presence of a sharp corner in the boundary of the entanglement region, the entanglement entropy (EE) and Rényi entropies for 3d CFTs have a logarithmic term whose coefficient, the corner function, is scheme-independent. In the limit where the corner becomes smooth, the corner function vanishes quadratically with coefficient for the EE and for the Rényi entropies. For a free real scalar and a free Dirac fermion, we evaluate analytically the integral expressions of Casini, Huerta, and Leitao to derive exact results for and for all . The results for agree with a recent universality conjecture of Bueno, Myers, and Witczak-Krempa that in all 3d CFTs, where is the central charge. For the Rényi entropies, the ratios do not indicate similar universality. However, in the limit , the asymptotic values satisfy a simple relationship and equal times the asymptotic values of the free energy of free scalars/fermions on the -covered 3-sphere.
1 Introduction and Results
For a 3d conformal field theory (CFT) in the ground state, the entanglement entropy for a region whose boundary has a sharp corner with angle can be written as
Here is a length scale associated with the size of the entangling region, is a short distance cutoff, and is a non-universal constant. The corner contribution to the entanglement entropy is the scheme-independent positive function of the opening angle [1, 2, 3]. Since the entanglement entropy of the region equals that of the complement region, the corner contribution satisfies . If the curve bounding the entangling region is smooth, the logarithmic term is absent, hence must vanish in the limit and it does so quadratically as
The value of the corner coefficient depends on the theory.
For the theory of a free real scalar or a Dirac fermion, Casini, Huerta, and Leitao [4, 5, 2] derived expressions that give implicitly in terms of some rather involved integrals. In the limit, one can extract double-integral expressions for the corner coefficient in (1.2). These integrals have been evaluated numerically [5, 6] and the results indicate that the exact values are 
for the free boson and free fermion, respectively.
The conjecture (1.4) has passed non-trivial holographic tests for gravity models with a family of higher derivative corrections [6, 7]. The central charge is defined as the coefficient of the vacuum 2-point function of the stress tensor (see eq. (3) in ). For free bosons and fermions, Osborn and Petkou  found that and in 3d. So with the values (1.3), the ratio is indeed for both free bosons and fermions.
In this paper, we evaluate analytically the integral expressions  of Casini, Huerta, and Leitao [4, 5, 2] for and and prove that their exact values are indeed those in (1.3). This verifies the universality conjecture (1.4) for the case of free bosons and fermions. One way of viewing the conjecture is simply as the statement that the corner coefficient in (1.2) does not contain independent information about the CFT, but is fixed in terms of the central charge .
Turning to the Rényi entropies , one can define a similar corner contribution which in the smooth limit goes to zero as for . (The limit of the Rényi entropy is the entanglement entropy.) It is not known if has any universal properties.
We calculate analytically for the free boson and free fermion using integral expressions for derived in [4, 5, 2].111We are grateful to Horacio Casini for sharing with us the integral expression for . For the free scalar we find
Note that when is even, the contribution from must be taken carefully using .
The result for the free fermion is
where sum is to be taken in integer steps from to .
In the case of the scalar, the exact results were guessed by the authors of  based on their high precision numerical evaluation of the integrals.
Since the ratios of the central charges of free fermions and bosons differ only by a factor of 2, universality of the ratio would require that obeys some simple, possibly -dependent, relation. Based on our results above, there is no hint of such a simple relationship. Of course to fully exclude this, one would need values of for other 3d CFTs.
As a function of , the Rényi corner coefficient decreases monotonically, as shown on the left in figure 1. When is large, asymptotes to a constant value, which we calculate analytically:
The appearance of the Riemann zeta-function is intriguing since also shows up in the free energies and Rényi entropies for free scalars/fermions on a 3-sphere, as shown by Klebanov, Pufu, Sachdev, and Safdi . Specifically, the free energy of a free real scalar or free fermion on an -covered 3-sphere behaves as222The authors of  work with a complex scalar, so the free energy there is twice that of a real scalar. for with
Thus, for both free scalars and fermions we have
For finite , there is no apparent relation between and , however there are some similarities in the subleading large- behaviors, as we discuss in section 4. The plot on the right in figure 1 shows the large- behaviors of the Rényi corner coefficients . A priori it is not clear if there is any relation at large between and , but it would be curious to test (1.9) in other examples.
The remainder of the paper details the derivations of the results summarized above. In section 2, we derive the results (1.3) for the entanglement entropy corner coefficient . We then evaluate the Rényi entropy corner coefficients in section 3. In section 4, we discuss the asymptotic behavior at large .
2 Evaluation of the EE integrals
In this section we describe the procedure for analytically evaluating the integrals for the coefficients and of the entanglement entropy. Our starting point is the integrals [4, 5, 2] presented in equations (B1)-(B3) of . After a change of integration variable from to , the integrals take the form
where for the scalar and for the fermion. The functions and are defined as
with333We simplified the expression for in  by writing it in terms of .
The functions , , , , and are defined as follows:
Here denotes the digamma function, .
Our first line of attack involves calculating the quantities and that appear in in (2.2). Beyond the immediate cancellations that occur in these ratios, one can perform further simplifications using identities involving gamma functions. Namely, one can use the recurrence relation
and the reflection relation
Surprisingly, all the gamma functions cancel after a series of such substitutions, giving
It is suggestive that the pre-factors and the form of these two results are the same. We then proceed by adding them together as in (2.2). The linear combination of digamma functions that appears in the result can be simplified using properties easily derived from (2.5) and (2.6). In the form that is useful for our purpose, these identities are
Then the combination of and that appears in simplifies to
Further simplifications of depend on the nature of variable , as we will see when we specialize to the cases of the free scalar and the free fermion.
Free scalar. To proceed with the evaluation of the integral , we set as prescribed for the free scalar. It is furthermore convenient to change integration variable . Using that both and are positive, the integrand of simplifies dramatically and becomes
Next, we integrate by parts. Writing
we see that the boundary term vanishes and we get
We have extended the limits of integration to facilitate the change of integration variables
This separates the two integrations and reduces the expression to
This completes the derivation of the result (1.3) for the free scalar.
Free fermion. With given in terms of as in (2.3), we have already done most of the leg-work needed to compute . For the free fermion, we have to take and it is again convenient to change integration variable . After putting everything together, we have
We can express the integrand as a total derivative plus remaining terms as
As before, the boundary term vanishes and we are left with the expression (after extending the limits of integration)
In the last step, we changed integration variables using (2.13). Since is odd, that part of the integral vanishes and the result is therefore simply
Thus we have derived the result (1.3) for the free fermion.
3 Rényi entropies
We now proceed to calculate the corner coefficients for the Rényi entropies.
Free scalar. For the scalar field, the Rényi corner coefficient is given by the integral (B7) in . We change of the integration variable to to write it as
As before, we write the integrand as a total derivative plus the remaining terms:
The boundary term vanishes and the expression simplifies to
The contribution of is easy to calculate and is equal to
For , there are contributions from two integrals:
Above, we manipulated the tri-logarithm using the polylog identity
which holds for .
Thus, having evaluated the integral in (3.4), we can write as the finite sum
Note that taking the limit as described below (1.5), the summand evaluates precisely to the special case (3.5). The expression (3.10) is the result for the Rényi corner coefficient presented in (1.5), so this completes our evaluation for the free scalar.
Free fermion. For the fermion field, the Rényi corner coefficient is given by the integral
where the sum is over from ( even) or ( odd) in integer steps to . Substituting the expressions for and obtained earlier gives
We then use integration by parts to simplify the integral
The boundary term integrates to zero and the expression simplifies to
The result of the integral again involves a difference of two tri-logarithms and it can be simplified using equation (3.8). The result is even in and we can write the final answer as
4 Asymptotic behavior of the Rényi’s
Let us now study the large behavior of the Rényi entropy corner coefficients . In particular, we evaluate analytically the value for the coefficients in the limit where . This is done by introducing a new variable and multiplying by . Then in the limit, the sum becomes an integral and we have
On the right in figure 1, we illustrated the asymptotic behavior of the corner coefficient which we find to be
Numerical fits show that , , and are for the free boson while is and in (4.2) for the free fermion. In fact, fitting up to , we find numerical evidence that for both the scalar and fermion. This indicates that a factor of can be factored out of the function in (4.2), so that
It is also interesting to study the ratios of the Rényi corner coefficients at large : based on numerical fits in the range to we find
The value of the -coefficient is guessed based on the numerics. Specifically, we fit to the function
and find that , , , , etc. The vanishing of the odd powers in (4.5) is consistent with (4.3). Note also that we can now identify the number from the fit (4.2) of the free fermion Rényi entropy corner coefficient at large as ; this is the value given in the caption of figure 1.
Again, the value of the -coefficient is guessed based on the numerics which give , , , , etc. The behaviors of individually is, however, very different that that of the Rényi corner coefficients. We find that while .
It is not clear whether the similarities observed at large between and have any significance or if it is a coincidence. Perhaps future investigations will clarify this.
We are grateful to Horacio Casini for sharing with us the integral expression for the corner coefficient of the free fermion. We would like to thank Pablo Bueno, Finn Larsen, Rob Myers, Silviu Pufu, and William Witczak-Krempa for useful discussions. HE is supported in part by NSF CAREER Grant PHY-0953232 and she is a Cottrell Scholar of the Research Corporation for Science Advancement. MH is supported by a Fulbright Fellowship and by the Department of Physics at the University of Michigan.
-  E. Fradkin and J. E. Moore, “Entanglement entropy of 2D conformal quantum critical points: hearing the shape of a quantum drum,” Phys. Rev. Lett. 97, 050404 (2006) [cond-mat/0605683 [cond-mat.str-el]].
-  H. Casini and M. Huerta, “Universal terms for the entanglement entropy in 2+1 dimensions,” Nucl. Phys. B 764, 183 (2007) [hep-th/0606256].
-  T. Hirata and T. Takayanagi, “AdS/CFT and strong subadditivity of entanglement entropy,” JHEP 0702, 042 (2007) [hep-th/0608213].
-  H. Casini, M. Huerta and L. Leitao, “Entanglement entropy for a Dirac fermion in three dimensions: Vertex contribution,” Nucl. Phys. B 814, 594 (2009) [arXiv:0811.1968 [hep-th]].
-  H. Casini and M. Huerta, “Entanglement entropy in free quantum field theory,” J. Phys. A 42, 504007 (2009) [arXiv:0905.2562 [hep-th]].
-  P. Bueno, R. C. Myers and W. Witczak-Krempa, “Universality of corner entanglement in conformal field theories,” arXiv:1505.04804 [hep-th].
-  P. Bueno and R. C. Myers, “Corner contributions to holographic entanglement entropy,” arXiv:1505.07842 [hep-th].
-  H. Osborn and A. C. Petkou, “Implications of conformal invariance in field theories for general dimensions,” Annals Phys. 231, 311 (1994) [hep-th/9307010].
-  I. R. Klebanov, S. S. Pufu, S. Sachdev and B. R. Safdi, “Renyi Entropies for Free Field Theories,” JHEP 1204, 074 (2012) [arXiv:1111.6290 [hep-th]].