Gravity duals of boundary cones

The replica trick defines Renyi entropies as partition functions on conically singular geometries. We discuss their gravity duals: regular bulk solutions to the Einstein equations inducing conically singular metrics at the boundary. When the conical singularity is supported on a flat or spherical surface, these solutions are rewritings of the hyperbolic black hole. For more general shapes, these solutions are new. We construct them perturbatively in a double expansion in the distance and strength of the conical singularity, and extract the vacuum polarisation due to the cone. Recent results about the structure of logarithmic divergences of Renyi entropies are reproduced ---in particular, $f_b\neq f_c$. We discuss in detail the dynamical resolution of the singularity in the bulk. This resolution is in agreement with a previous proposal, and indicates a non-minimal settling to the `splitting problem': an apparent ambiguity in the holographic entropy formula of certain theories with higher derivatives.


Introduction and summary
The Ryu-Takayanagi area formula [1,2] seeds many recent insights from quantum information into quantum gravity. This formula calculates the entanglement entropy of a field theory region as the area of a minimal surface in the gravity dual, generalising the Bekenstein-Hawking entropy.
Building on work by Casini, Huerta and Myers [3], Lewkowycz and Maldacena have constructed a derivation of the Ryu-Takayanagi formula. This derivation, called 'generalized entropy' [4], uses the replica trick in field theory and its gravity dual, which is in the Euclidean quantum gravity regime.

JHEP09(2016)139
In this context, the replica trick defines entanglement entropy S as the n → 1 limit of Rényi entropies S n . Using Euclidean techniques, the S n can be related to field theory partition functions on conically singular manifolds [5]. The conical singularity has support on the entangling surface that bounds the partition whose entanglement is being calculated. The period of the cycle contracting on the entangling surface is 2πn, and the conical singularity is absent when n = 1.
In a CFT, conformal symmetry can be used to send this entangling surface to infinity. Doing so leaves behind an Euclidean field theory geometry with a non-contractible cycle. This makes the setup analogous to that of thermal field theory. When these thermal states have black holes as gravity duals, arguments from Euclidean quantum gravity apply. An area law emerges as a Gibbons-Hawking derivation of black hole entropy [6] -in which regularity of the bulk Euclidean geometry for all values of the temperature (the period of the thermal cycle) plays a central role. [3] implemented these ideas for a spherical entangling surface in the vacuum of a CFT in flat space. This conformally maps to the Euclidean thermal hyperboloid, whose gravity dual is well known -it is the hyperbolic black hole [7]. Entanglement entropy in the flat space picture is related to thermal entropy in the hyperboloid, in which case it is dual to black hole entropy. Using Bekenstein-Hawking for the thermal case, a Ryu-Takayanagi area follows for the entangling one. [4] extended this picture to more general states and entangling regions. This involves field theory geometries for which the entangling cycle no longer generates a symmetry, and thus the connection to thermal physics weakens. The subtleties of the replica trick become more prominent, and the analytic continuation of the dual geometries to non-integer values of n, relevant for the n → 1 limit, is less direct and needs to be discussed in detail.
Both [3] and [4] work in 'hyperbolic frames', in which the boundary conical singularity has been mapped to infinity. The primary goal of this paper is to carry out explicitly the construction in the 'entangling frame', in which case the conical singularity of the boundary at n = 1 remains within sight, but nevertheless has a regular gravity dual. We will exhibit in detail these geometries, and discuss how exactly gravity in the bulk dynamically regulates boundary conical singularities. 1 The new geometries are regular Euclidean solutions of the Einstein equations with a negative cosmological constant, subject to the boundary condition that the geometry induced at the conformal boundary has a conical singularity on a specified surface. To make analytic progress, we will find these geometries perturbatively in a double expansion in the distance and strength of the singularity, (n − 1). The distance is measured in units of the smallest lengthscale characterising the geometry of the background and surface supporting the singularity. We focus on the case of five bulk dimensions for convenience. An extension to general dimensions would be very interesting.
These geometries encode the vacuum polarisation of holographic conformal field theory due to conical singularities on general surfaces. Although strictly speaking there is no Fefferman-Graham expansion because the boundary metric is singular, we will extract an JHEP09(2016)139 expectation value for the stress tensor. In our approximations, and ignoring contact terms, the stress tensor is traceless.
The Euclidean action of these bulk geometries relates to Rényi entropy S n , that we extract to first order in (n − 1). We find that, contrary to the entanglement entropy term S 1 , the first order term does not have an area law. Also, in precise agreement with [10], the logarithmic divergence of the first order term turns out to have a different structure from the entanglement entropy one.
The specific way in which gravity regulates these singularities impacts the generalisation of the Ryu-Takayanagi formula to theories with higher derivatives. This formula extends Wald's black hole entropy to setups without U(1) symmetry. For theories without explicit derivatives of the Riemann tensor in their lagrangian, 2 it takes the schematic form [12,13]: where L is the gravity lagrangian, K is the extrinsic curvature, and q α are coefficients characterising how exactly the conical singularity is regulated in the bulk (this is reviewed in section 7). We will see below that they differ from a minimal prescription. The rest of this paper is organised as follows. Section 2 reviews the Casini-Huerta-Myers construction and exhibits that Euclidean hyperbolic black holes can be written as smooth gravity duals to straight conical singularities. In the rest of the paper we will deform this cone away from straightness and explore the consequences of its gravity dual. Section 3 reviews the construction of Fermi-like coordinates adapted to codimension two surfaces, and the natural implementation of the replica trick in these coordinates. Section 4 constructs explicitly the gravity duals, and is the core of the paper. There are many ways to deform a surface away from straightness (equivalently, many ways to squash a cone [14]), and, after a general overview, we proceed in a casuistic way. To the order of Riemann curvature, this results in thirteen subsections analysing different such deformations. 3 Section 5 summarises the results regarding the vacuum polarisation induced by these singularities, and section 6 explains how to reproduce the results of [10] regarding logarithmic divergences of Rényi entropies in CFT. In section 7 we discuss the consequences of section 4 for the entropy formula of theories of gravity with higher-derivative interactions. We conclude in section 8.

Hyperbolic black holes and boundary cones
This section reviews some aspects of the Casini-Huerta-Myers construction [3] and hyperbolic black holes. It also serves to set notation and discuss coordinates that we will be using throughout.

JHEP09(2016)139
Consider the Rényi entropy of a 4D conformal field theory in the vacuum across a straight plane. 4 The replica trick maps this quantity to the Euclidean partition function on the conically singular geometry: where we took complex coordinates ζ = σ 1 + iσ 2 on the entangling plane, at r = 0, and the geometry is singular for n = 1 because of the period of τ . A convenient way to write cone metrics uses complex coordinates also in the plane of the cone. Define z = r 1/n e iτ /n . This is a good complex coordinate when τ has the period in (2.1), and the cone becomes which can be obtained from Euclidean space by the 'quotient' z → z n . This conical geometry (2.1) is conformal to the Euclidean thermal hyperboloid S 1 ×H 3 : which is regular for all periods of τ .
If the CFT has a GR gravity dual, the geometry dual to the partition function on (2.3) is the hyperbolic black hole [7]: 5 The metric at the boundary, at ρ → ∞, is (2.3). τ closes smoothly at ρ = ρ h when τ ∼ τ + 2πn. Hence this is a good holographic dual to (2.3). The geometry becomes AdS when n = 1, when the temperature is 1/2π in units of the radius of the hyperboloid. Boundary conformal transformations are implemented by large diffeomorphisms in the bulk, so there is a change of coordinates that writes the geometry (2.4) as the gravity dual of the conical singularity (2.2), [9]. To leading order in (n − 1), one such diffeomorphism is: Now (2.4) reads: In these coordinates the boundary is at ρ → 0. Its geometry is indeed conically singular, because for ρ → 0: and we recover (2.2) as the boundary metric. The axis of the cone extends into the bulk in a regular manner: at fixed ρ, as zz → 0, a constant, and thus regular. The factors of (n−1) 1+x are bounded corrections that remain small everywhere for small (n − 1). At fixed ρ, eq. (2.8) is a regularised version of the boundary conical singularity, to which it tends far from the axis, x ≫ 1. Therefore, eq. (2.8) is a boundary conical singularity that gravity dynamically regularises in the bulk.

Fermi coordinates for field theory replicas
Following the program of generalized entropy [4], we will find the gravity duals to cones on entangling surfaces and backgrounds more general than the plane in flat space.
One can take coordinates adapted to a generic surface by shooting geodesics orthogonal to it. In an expansion near such surface the metric can be written as: where, as above, we parametrise the transverse geodesics in complex coordinates: z,z. The objects γ ij , K ijz , Q ijzz , etc., characterise the embedding and background geometry on the surface and may depend on its coordinates σ i . c.c. stands for complex conjugation within the square brackets. κ is a small book-keeping parameter counting powers of the distance to the surface in units of the characteristic lengthscale of the geometry. As in the previous section, it is convenient to take complex coordinates also in the surface and analogously expand around their origin. By appropriately constructing these coordinates, one can eliminate A izz and its symmetrised first derivatives at the origin, as well as the

JHEP09(2016)139
Christoffel symbols of the induced metric γ ij : where κ now also keeps track of the distance to the origin in the surface, ζζ = 0. R µνρσ denotes components of the background Riemann tensor, and R ijkl refers to the Riemann tensor of γ ij -the intrinsic curvature of the zz = 0 surface. A convenient feature of complex coordinates is that trace and traceless elements of symmetric tensors are readily distinguished. For example, K ζζz belongs to the trace of K ijz , whereas K ζζz belongs to its traceless part.
The use of these coordinates on entanglement entropy calculations in field theory was pioneered in [15,16], and [17] for Rényi entropy. We expand to O(κ 2 ) because we are interested in effects due to Riemann curvature. The Riemann tensor of (3.1)-(3.2) at ζζ = zz = 0 is fully captured at this order, and reads (setting κ = 1):

3)
R izzz and R zzzz coming directly from (3.1). Only some combinations of the above objects transform covariantly under conformal transformations. When working with CFTs, one may take advantage of such symmetries to eliminate non-covariant elements. For example, one can choose to eliminate the trace of the extrinsic curvature K ζζz , as well as the intrinsic Riemann R ζζζζ , and the traces Q ζζzz and Q ζζzz . For reasons that will be clear, we will drop K ζζz but will keep the R and Qs to keep track of the topology of the surface (via Gauss-Bonnet) and as a device to check the conformal covariance of our results. One check of such covariance will be the appearance of the trace of the projection of the bulk Weyl on the entangling surface: Now, as earlier, replicating around the zz = 0 surface is implemented by z → z n . We will then be after the gravity duals of partition functions of CFTs on with γ ij and A izz as in (3.2).

JHEP09(2016)139
z n is multivalued for generic n, and so (3.5) is a geometry only for n a positive integer (in particular, only then ds 2 is continuous for z in the complex plane). We will nevertheless treat n as a real number, expand ds 2 in powers of (n − 1), and speak about a CFT in (3.5) for real n. This is tantamount to speaking about a CFT on an 'analytic continuation of a geometry', and it is an abuse of language. Quantities of interest for n ∈ R should be thought as being analytically continued from n ∈ Z, as Rényi entropies in the replica trick. The usefulness of this picture is that it picks the right analytic continuation for the quantities of interest [4].

Gravity duals of squashed conical singularities
This section is the core of the paper, where the regular gravity duals to (3.5) are spelled out, to first order in (n − 1) and to second order in κ. We start by giving an overall picture of these geometries and their properties.
The expansion in κ is a derivative expansion, and for this reason our results may be reminiscent of other such expansions in gravity, as fluids/gravity [18,19] or blackfolds [20,21]. The strategy is to deform the boundary metric in (2.8) in the expansion of (3.5). This generically does not solve the Einstein equations by itself, and one needs to add a small correction to that end. We then solve for this correction subject to the boundary conditions of normalisability and regularity.
The notion of bulk regularity we are alluding to at non-integer n is not a standard one if there is squashing. The metrics we will call 'regular' have 1/r divergent curvature invariants at the bulk axis, 6 but (i) these divergences drop from the field equations, and (ii) they are altogether absent at integer n. Condition (ii) is the usual AdS/CFT bulk regularity at integer n, while (i) picks the right analytic continuation of the bulk metric to non-integer n [4]. In a sense, this is the closest we can get to the usual notion of regularity for n ∈ R.
We will assume that all the non-trivial dependence of the corrections is in the dimensionless x of eq. (2.7). The geometry then does not have any more dependence on the angle in which we approach the singularity other than the one following from the index structures in (3.5). This reduces the equations for the corrections to ODEs. To leading order in (n − 1), our notion of regularity at the bulk axis, x = 0, boils down to the expandability of the metric in non-negative powers of x 1 n and x n−1 n [22]; 7 and that for g zz such expansion has a constant term, (2.10). This gives smoothness at positive n ∈ Z. Assuming dependence just in x implies replica symmetry -a discrete rotational symmetry in the plane of the cone, z → z e ik/n for k ∈ Z n . 8

JHEP09(2016)139
This dependence on x implies covariance under diffeomorphisms in the surface, that we will maintain explicitly. This is a useful principle when writing ansatze for the bulk corrections, because it forbids appearances of the surface coordinates ζ other than the ones in (3.2).
The geometries we will find are exact in x ≡ (zz) n /ρ 2 . Given this and the expansion in powers of κ, the range of ρ -the holographic radial coordinate -needs also be small in units of the boundary curvature. This means that the expansion in κ is also an expansion around the conformal boundary. The solutions we will write down are then analogs to Fefferman-Graham expansions, but for conically singular boundary metrics. For this reason, the field theory properties we will extract from them are approximate and belong in a UV expansion.
We will also further fix the gauge by requiring that the metric has no derivative corrections at the axis in the legs involving dz. That is, e.g., that there are no κ corrections to eq. (2.10). As a gauge choice, this does not affect the geometry. Its advantage is that it is straightforward to read the geometric properties of the bulk axis.
To illustrate this language, we now describe how the solution (2.8) would look like if we were finding it in this way. Starting from AdS in Poincare coordinates, and introducing a cone in the boundary by z → z n , we would write with ds 2 1 a small correction. Taking the ansatz Einstein equations plus boundary conditions give: which is indeed the O(n − 1) expansion of (2.8). While the 1 1+x terms are clearly regular at x → 0 and normalisable at x → ∞, the log term may look problematic at the axis. However, this log is exactly what is needed to balance the zero of (zz) n−1 in (4.1), so that we are left with the finite result of eq. (2.10). The key point is that the (zz) n−1 and the log can be grouped into (1 + x) n−1 n of eq. (2.8), exhibiting the desired contrasting behaviours at small and large x -the axis is regular in the bulk but not in the boundary.
We will find similar log divergences below, and we will have to interpret whether they reflect singular behaviour or not. They may just indicate that the ansatz following from replicating at the boundary does not capture a specific regular behaviour near the bulk axis. Sometimes this can be anticipated, like for g zz above. Consider, e.g., the case of F ζζzz z dz.

JHEP09(2016)139
After replicating, this term is accompanied by a factor of (zz) n−1 at the boundary. In the bulk, a minimal replica symmetric ansatz near the axis does not have such factor, and we may expect that this is the behaviour chosen by dynamics [12,13]. For these cases, appropriate factors of (1 + x) n−1 n in the bulk ansatz can offset the appearance of logarithms in ds 2 1 . We will anticipate so in a few cases below, except for the R zzzz and Q ζζzz cases; exhibiting their logarithms will be useful for understanding regularity in the more delicate case of K ζζz Kζζz.
Each term in the expansion in κ of (3.5) needs its own derivative corrections. For readability, we will present all these contributions separately. It is straightforward to put them together.
We will start in 4.1 with the O(κ) corrections, due to traceless extrinsic curvaturewe remind the reader that we exploit conformal symmetry in the boundary to drop the trace of the extrinsic curvature. We will then move on to O(κ 2 ) terms. Some of these are seeded by squares of extrinsic curvatures, and others are sourced linearly by κ 2 terms in (3.5) -including those implicit in the derivatives of the extrinsic curvature.
We will present the O(κ 2 ) contributions in an order that groups them by their tensorial character in the parallel and transverse coordinates ζ and z. In 4.2 we present the correction due to K ζζz,ζ , which might be called the 3|1 contribution, because it has three holomorphic indices in ζ and one in z. 9 In 4.3 we present the correction due to F ζζzz , that may be called axial-axial because it is antisymmetric in both pairs of indices. In 4.4 we move on to the 2|2 term, Q ζζzz ; followed in 4.5 and 4.6 by the 0|2 ones, Q ζζzz and K ζζz Kζζ z ; 2|0 in 4.7, Q ζζzz ; 4|0 in 4.8, K ζζz K ζζz ; and 1|1 in 4.9 and 4.10, K ζζzζ and R ζzzz . The last four contributions 4.11-4.14 are the 0|0 ones: R ζζζζ , R zzzz , Q ζζzz and K ζζz Kζζz. The last one has the subtlest log structure, and its regularity has consequences for the splitting problem.
Sections 6 and 7 only use the results of 4.1 and 4.11-4.14, and some readers may want to focus on these.
Had we not used conformal symmetry to drop the trace of the extrinsic curvature, there would be five more cases at second order: K ζζz,ζ , K ζζz K ζζz , K ζζz K ζζz , K ζζz K ζζz , and K ζζz K ζζz .
Reality implies that any geometry containing, e.g., the 0|2 correction K ζζz Kζζ z must also contain the corrections leading to the 0|0 correction K ζζz Kζζz. We will leave the latter implicit when displaying the results for the former, which means that formally we will be writing down complex metrics. This is just an artefact of the presentation in terms of ζ and z tensor behaviour.
For all cases we will present the expectation value of the stress tensor due to the particular squashing of the cone. This is obtained by conventional holographic methods (see section 5 for more details, discussion, and a comprehensive expression for T ).

K ζζz
Before replicating, the boundary metric for this term is: We replicate by z → z n , leading to the bulk ansatz: In (4.5) we choose K ζζz to multiply the same 1 1+x factor as g ζζ for later convenience, but that is not significant, as these factors are precisely what f K2|1 ζζ (x) is designed to discover. Notice that, since ds 2 1 is linear in κ, it has to be proportional to K ζζz ; and, since this is traceless, it can only seed a dζ 2 leg if we want to avoid explicit appearances of ζwhich we do because of covariance in the entangling surface. The important feature of this ansatz (4.5) is that the boundary metric, at ρ → 0, is different from the one in (2.8) by the factor of K ζζz z n of (3.5).
Einstein's equations lead an ODE for f K2|1 ζζ (x), whose normalisable solution is: with C K2|1 an integration constant. 1/x is singular at the axis, and regularity sets C K2|1 = 0. Restoring the AdS radius ℓ, the stress tensor reads 10 This term follows directly the one we just analysed, from covariance in the ζ coordinate. Its boundary metric is, before replicating

F ζζzz
The boundary metric for this term is, before replicating, After z → z n , a natural ansatz for the bulk is: with Notice the factor of (1 + x) n−1 n , anticipating different behaviours of the F ζζzz term in the boundary and bulk axes. In the boundary, we require the (zz) n−1 behaviour from the replica trick (3.5); in the bulk, a minimal guess suggests (zz) 0 .
The normalisable solution we find is: and it should be set to zero because it is not regular at the bulk axis x → 0.

JHEP09(2016)139
The O(κ 2 ) contribution to the stress tensor coming from this term is:

Q ζζzz
This one comes from Then, with The normalisable solution to Einstein equations is: which regularity sets to zero. The contribution to the stress tensor vanishes:

Q ζζzz
This one comes from Then,

JHEP09(2016)139
with (4.25) This case has new features compared to the previous two. First, (4.24) has a κ 2 term in the dρ dz leg that remains finite as n → 1. This term does not change the boundary metric and therefore is a bulk response to Q ζζzz , even before introducing any conical singularity. Its origin is, in fact, well known. In the Fefferman-Graham expansion it is the Shouten term, 11 which is indeed non-zero for (4.23). As explained in the opening of this section, our gauge demands that g zz vanishes on the bulk axis. This places this Shouten term in the dρ dz leg (in Fefferman-Graham coordinates this would have had a dz 2 leg).
Second, the scalar character of Q ζζzz in the ζ directions allows for many more terms in ds 2 1 compared to earlier. In fact, covariance would also allow for a dz 2 leg that we have not written down. As it turns out, a gauge transformation can move this correction from the dz 2 leg to dz dρ, and, as explained, our gauge fixing places it in the latter.
There are two normalizable zero modes to this ansatz: 26) and they are banished by regularity at the axis. The regular solution we find is: The Fefferman-Graham expansion of AAdS spacetimes reads, for small ρ [24]: (0) gµν is the Shouten tensor of the boundary metric (for boundary dimension d > 2). Notice that the Fefferman-Graham expansion is a derivative expansion, so the expansion in κ has a FG character.

JHEP09(2016)139
The contribution of Q ζζzz to the stress tensor is: This one comes from Then, with Again, the third line in (4.30) comes from the Shouten tensor, and there is no dz 2 leg because of gauge fixing. Discarding the same zero mode as earlier, (4.26), the regular solution we find is: The contribution to the stress tensor: (4.33)

Q ζζzz
This one comes from Then, (n−1) 1+x with Again, notice a Shouten correction in the brackets. Notice also that we have anticipated a change of behaviour of Q ζζzz from the boundary, (zz) n , to the bulk axis, zz. f Q2|0 ζζ (x) = 0 is the regular solution we seek. There is also one singular normalisable zero mode: Stress tensor contribution: (4.38)

K ζζz K ζζz
This one comes from Then, but ds 2 1 = 0 because the four legs in ζ of K ζζz K ζζz would force at least two contractions with ζ, which clashes with covariance in ζ. This case, as the one in 4.2, follows from the one in 4.1.
Therefore, T = 0 . Then, Notice a Shouten correction in the fourth line. There are two normalisable zero-modes: x . (4.45) We set C 1|1 1 = 0 for regularity at the axis, and C 1|1 2 = 0 with the gauge condition g zµ = 0 at the axis.
The regular particular solution is: Stress tensor contribution: The first of these terms follows from covariance in (4.8). The rest are new.

R ζzzz
This one comes from

JHEP09(2016)139
Then, Notice a Shouten correction in the fourth line. There are two normalisable zero-modes: x . (4.52) We set C with the gauge condition g zµ = 0. The regular particular solution is: Stress tensor contribution: (4.55)

R ζζζζ
This one comes from

JHEP09(2016)139
Then, There are two normalizable zero modes to this ansatz, that we will omit henceforthalso for the three other 0|0 cases that follow. They look (4.59) C 0|0 2 = 0 with the gauge condition of no corrections in κ to g zz at the axis, and C 0|0 1 = 0 with regularity. Notice that C 0|0 1 = 0 would result in a logarithmic divergence at the axis in, e.g., the dz dz leg. This leg has no κ 2 term surviving the n → 1 limit, so no candidate to absorb the logarithmic divergence in a change of behaviour of the exponent, as is the case for the logarithmic divergences that we have been allowing.
The regular solution is: (4.60) Stress tensor contribution:

R zzzz
This one comes from Then, after replicating, we write the following ansatz As earlier, Shouten corrections appear in brackets. Notice that we have not written down (1 + x) n−1 n factor in the R zzzz factor second line, even though we anticipate it. The reason for not writing it will be clear in section 4.14.
The regular solution we find is: Notice the appearance of the logarithm, that could be absorbed in the ansatz by writing, instead of R zzzz (zz) n in the second line of (4.63), R zzzz zz (1 + x) n−1 n ρ 2 n−1 n . Thus, this logarithm does not reflect singular behaviour.
Stress tensor contribution:

Q ζζzz
This one comes from

JHEP09(2016)139
Then, after replicating, Notice again the Shouten correction inside the brackets, and the absence of a likely (1 + x) n−1 n factor in the third line of (4.68). After discarding the singular zero mode of (4.59), we find the regular solution: (4.70) The apparently singular log in f Q0|0 ζζ (x) can again be absorbed in the ansatz (4.68) by replacing Q ζζzz (zz) n in the third line by Q ζζzz zz (1 + x) n−1 n ρ 2 n−1 n . The contribution to the stress tensor: (zz) n−1 z 2 dz 2 +z 2 dz 2 zz . (4.71)

K ζζz Kζζz
This one comes from Then, n−1 n ρ 2/n dz dz The κ 2 ρ 2 terms correspond to the Shouten correction to AdS if we were using FG coordinates. This just depends on the boundary geometry and appears at n = 1.
Discarding the zero mode of (4.59), we find: As advertised at the beginning of the section, regularity at the bulk axis now appears to be subtler, as there are no obvious terms with which to absorb these logarithms into (1 + x) n−1 n factors.
It helps to notice, however, that the logarithms appear in the same legs as for the two previous cases (4.65), (4.70). They can then be cancelled by adding them to the current case, fine-tuned as leading to a regular solution. This solution implies, however, that whenever K ζζz Kζζz is not zero there will be a part of R zzzz and Q ζζzz that is not accompanied by the (1 + x) n−1 n factors, unlike in sections 4.12 and 4.13. This part, unlike the rest, will have a factor of (zz) n−1 near the bulk axis.
In practice, this amounts to absorbing the logs in (4.75) by adding to the ansatz (4.73) the following: (4.77) Thus, the logs in (4.75) are not signalling singular behaviour. We will elaborate on this in section 7.
The contribution to the stress tensor reads: (4.78)

Summary
In this section we have written down explicitly smooth gravity duals to all possible squashings of a boundary cone, to the order of Riemann curvature -with the exception of those involving the trace of the extrinsic curvature, that we set to zero without loss of generality using conformal symmetry. It is a rather lengthy section because of its exhaustiveness. The cases we have not written down explicitly follow straightforwardly from the ones presented by exchanging ζ ↔ζ and/or z ↔z.
In the remainder of the paper we will discuss the consequences of these bulk geometries and benchmark them against known features they should reproduce.

Vacuum polarisation
In the preceding section we presented the vacuum polarisation T induced by a squashed conical singularity in 4D holographic conformal field theory. This expectation value of the stress tensor was obtained by conventional holographic techniques [24], which for us simplify to reverting to Fefferman-Graham-like coordinates and selecting the coefficient of the ρ 4 term (see footnote 11).
However, strictly speaking the bulk geometries of this paper do not have well defined Fefferman-Graham expansions, because the boundary metric is singular. It is then necessary to explain in which sense the method used calculate this vacuum polarisation is conventional and valid.
Let us first recall one aspect of the boundary metrics (3.5). These are obtained by 'quotiening' the regular ones in (3.1), by z → z n . The quotient is, away from zz = 0, locally a change of coordinates. As such, it does not alter local properties of the geometry, as the curvature. Hence, the Riemann tensor is regular away from zz = 0. However, the quotient does introduce singular, delta-like, contributions to the curvature at the originthe tip of the cone.
Up to these contributions and the multivaluedness discussed at the end of section 3, the boundary geometry behaves regularly, and one can formally develop the Fefferman-Graham expansion and extract a stress tensor. But one needs to bear in mind that the stress tensor calculated in this way ignores contact terms.
As emphasised in the beginning of section 4, the T we have extracted should be thought of as belonging to a UV expansion. Indeed, notice that it diverges at the origin as where we momentarily reverted to polar coordinates r 2 = zz.
The T k contributions we have presented are all local functions of the geometry. Such local probes are characteristic of UV expansions. 12 The r → 0 divergences conform one such expansion. Generically, one also expects there to be finite, non-local dependence on

JHEP09(2016)139
The stress tensor induced by the introduction of the cone on (5.3) (by z → z n ) can be read from (5.2). Its dependence on κ also displays conformal flatness: providing a check of the good conformal covariance properties of (5.2). Similarly, the other two locally conformally flat cases, R ζζζζ = −R zzzz and Q ζζzz = 0, can also be seen to follow from the one without bending, κ = 0.

Logarithmic divergences of holographic Rényi entropy
The metrics dual to squashed cones reproduce the results of [10] regarding logarithmic divergences of Rényi entropy for holographic CFTs. Entanglement entropy is known to be UV divergent in field theory. This divergence is due to correlations across the entangling surface between infinitely many short distance degrees of freedom. It is therefore localised around the entangling surface. Taming it with a short distance cutoff, it reads, for the vacuum of a 4D CFT [25]: where the area and the integrals are on the entangling surface. R is the Ricci scalar of the induced metric γ ij ; W ij ij is the contraction of the projection of the Weyl tensor on the surface; and K {ij}a K {ij}a is the contraction of the square of the traceless part of the extrinsic curvature. The last two are conformal invariant. While the coefficient of the area term is sensitive to the choice of cutoff, the logarithmic divergence is not. It therefore has a physical character. a and c are the central chargesthe logarithmic divergence can be derived from the conformal anomaly, when the latter is written as the logarithmic divergence of the regulated effective action [25].
Rényi entropies are conjectured to have a similar UV behavior [26]: There are known relations between f a (n) and f c (n) [17,27], but less is known about f b (n), apart from f b (1) = c. Free field theory results prompted the conjecture that f b (n) = f c (n) [28]. However, this relation fails for holographic theories [10]. The method of [29] applied to our metrics for duals to squashed cones reproduces this failing. This method builds on [9] and [4] to argue that a certain derivative of Rényi entropy with respect to the index n is given, for theories holographically dual to General Relativity, by the area of the bulk axis:

JHEP09(2016)139 7 Splitting problem and singularity resolution in the bulk
This section discusses the impact of the analysis of section 4.14 on the holographic entanglement entropy formula of higher-derivative theories of gravity. We will conclude, in precise agreement with [32], that the 'splitting problem' has a non-minimal solution in a class of theories, resulting in a slightly different entropy formula from previously anticipated in [12,13]. This difference is visible only beyond curvature squared interactions, and does not impact the entropy of Lovelock nor f (R) interactions. We include a brief but self-contained description of the entropy formula for higher-derivative gravity and its splitting problem. The application of generalized entropy to higher-derivative theories of gravity results in a new holographic entanglement entropy formula. For the class of theories with a lagrangian depending on the Riemann tensor but not on its derivatives, 15 √ g d 5 x + Boundary terms, (7.1) this formula is Here R zzzz , R zizj , K ijz and their complex conjugates refer to an expansion of the type (3.1) around the bulk entangling surface, at zz → 0. This surface is what we called the 'bulk axis' in previous sections. In contrast to the use of the expansion (3.1) in section 3, i and j now run over three values, that in the coordinates of that section would be ζ,ζ and ρ.
The first term in this holographic entropy formula is Wald entropy [33], and the second one can be thought of as a correction to it. Wald entropy was constructed on bifurcation surfaces of event horizons, which necessarily have vanishing extrinsic curvature. The entropy of (7.2) applies also to situations in which the extrinsic curvature may be non-zero.
To explain the meaning of the sum in α in (7.2) we need to discuss some details of the application of generalized entropy to the class of theories (7.1). This application involves evaluating actions of bulk geometries 16 that regulate conical singularities. The prescription of [12,13] for the sum in α assumes a 'minimal' regulation, of the type discussed below eq. (4.3). Here minimal means that, in an expansion around the axis (at zz = 0), the metric of the regulated cone takes the form: This follows from taking (3.1) and promoting any holomorphic factors of z and dz that are not paired with antiholomorphic ones to z n and d(z n ) = nz n−1 dz, respectively. This

JHEP09(2016)139
achieves a replica symmetric metric (7.3) that is regular at the axis for integer n, and for which the exponents in z differ minimally from the ones before replicating (3.1).
The entropy following from the evaluation of the action of such cone-regulating geometries involves integrals of the type lim n→1 ∂ n ∞ 0 (n − 1) 2 r n r 2 2 r 2qα(n−1) e −r 2 r dr = 1 2 where we used polar coordinates r = √ zz. The role of the exponential function is to localise around the axis, and this function could be replaced without change in the r.h.s. by any other regulating function, interpolating smoothly between 1 at the origin and 0 at infinity -e.g., 1 1+r 2 . The integral in eq. (7.4) is dominated by a logarithmic divergence at the lower end as n → 1. That explains the independence from the regulating function. The outcome is sensitive to the details of the expansion of the geometry around the axis, that are encoded in the q α in (7.4). q α parametrises n−dependence in the power of r in the integrand, reflecting n−dependent exponents of z andz in the geometry (7.3).
Terms in the expansion (7.3) that are accompanied by different powers of z andz contribute differently to q α . That is what the sum in α in (7.2) captures. In this formula we need to decompose the second derivative of the lagrangian in monomials of the curvature, that α labels. These monomials are, however, not in the background Riemann tensor, as may appear natural for L(Riem). Rather, its constituents are the quantities appearing (7.3), in terms of which one can write the Riemann tensor (3.3): R ijkl , K ijz , Q ijzz , Q ijzz , R izzz , F ijzz and R zzzz (and complex conjugates). Each monomial α is then assigned a value of q α , and the sum is performed with the 1 qα+1 weight. Constituents of α contribute additively to q α with a weight that depends on the exponent of z andz that they are accompanied by around the axis of the regulated cone. For the regulation of (7.3), q α is contributed 1/2 for each K ijz and R izzz , 1 for Q ijzz , and 0 otherwise.
Note, however, that exchanging, e.g., the zz factor accompanying R zzzz in (7.3) for (zz) n would also achieve a regular replica symmetric metric, although with a different weight of this term in (7.2). Now, R zzzz would contribute 1 to q α , instead of 0. The obvious such ambiguities are in terms with a zz pair in their indices: Q ijzz , F ijzz , R izzz and R zzzz ; although there may be more [22]. These ambiguities have been called 'the splitting problem' [32].
A lesson that follows from the analysis in section 4.14 is that, in General Relativity, the expansion around the bulk axis does not take the form of eq. (7.3). Rather, the Q ζζzz and R zzzz terms look: 17 Notice that upon taking n → 1 the extrinsic curvature contributions cancel and we recover (3.1). Let us for a moment discuss what is the σ i −covariant version of eqs. (7.6) and (7.7). Recall that, in this section, σ i encompasses ζ,ζ and ρ. This is irrelevant for the second term, that can be written because for the configurations we studied in section 4, K ρiz = 0. Notice also that in this section K ijz can not have a trace -as the Ryu-Takayanagi surface is a minimal surface. For Q ′ a similar argument implies that the covariantisation should read: Equivalent expressions for the analogs of Q ′ and R ′ were found in [32] by solving the Einstein equations around the bulk axis. 18 Since the factors of (zz) n in (7.5) are different from those in (7.3), we conclude that the splitting problem has a non-minimal solution in General Relativity. This translates into the α sum of (7.2) meaning something different than it would in the minimal case of (7.3). Now, α labels monomials in terms of R ijkl , K ijz , Q ijzz , Q ′ ijzz , R izzz , F ijzz and R ′ zzzz ; instead of Q ijzz and R zzzz . Explicit factors of K ijz and R izzz still contribute 1/2 to q α ; Q ijzz contributes 1; and the rest, including Q ′ ζζzz and R ′ zzzz , contribute 0. This may appear irrelevant, since the lagrangian of GR has a vanishing second derivative in the Riemann and therefore no splitting problem; its entropy is just the area. However, this splitting does have consequences for the entropy formula of theories that contain perturbative higher-derivative corrections to General Relativity. For these corrections there is a splitting problem, and the splitting is fixed by the leading result -the GR one we just discussed.
This affects Riem k interactions for k ≥ 3, but does not have consequences for Lovelock interactions, because in those ∂ 2 L ∂R zizj ∂Rz kzl does not depend on Q ζζzz nor R zzzz ; 19 or for f (R) interactions, for which the second derivative vanishes identically.

JHEP09(2016)139
where the 3 is a symmetry factor and there is a factor of 2 3 from three g zz . From here we get ∂ 2 L ∂R zizj ∂Rz kzl K ij z K klz = −3λ R zζzζ K ζζz Kζζz + ζ ↔ζ , (7.12) and with the splitting of (7.5): in agreement with [32]. This could be simplified further using the 0th order background eoms (Einstein's). 20 For the second line we have used eq. (3.3), and for the third, (7.6). In the fourth line we summed over α with the splitting we have discussed, and in the fifth and sixth we have used eqs. (7.6) and (3.3) again. For the 'minimal splitting', we would perform the sum in α directly from the second line, getting in the end −3λ R zζzζ − K ζζz Kζζ z K ζζz Kζζz + ζ ↔ζ, for which the K 4 term is a different tensor structure altogether.

Outlook
This paper has described the regular GR duals to CFTs on squashed cones. These metrics show how bulk gravity regulates a conical singularity in the boundary. We have worked to first order in the strength of the cone (n − 1), and to second order in κ, parametrising an ultralocal expansion around the cone. Going to second order allows sensitivity to Riemann curvature -although only within a UV expansion. A quantity that follows from these geometries is (the UV expansion of) the vacuum polarisation in the presence of these cones (5.2), up to contact terms. This stress tensor bears some resemblance to those of fluids/gravity, and one can interpret its many coefficients as response coefficients to the squashing of the cone. The number of such second-order coefficients is large compared to [19] because the entangling surface breaks O(4) symmetry

JHEP09(2016)139
to two planes (parallel and transverse). Since we did not take full advantage of conformal symmetry, three of the contributions to (5.2) can be generated via covariance under conformal transformations, as we saw in section 5. For closure, it would be interesting to write down the missing contact terms, that may be interpreted as defect degrees of freedom.
Our setup should not be confused with that of holographic entanglement entropy across surfaces with singular shapes [34]. In that setup, the background boundary metric is regular.
The metric of a conical singularity at the boundary is simple both in complex coordinates (2.2) and in hyperbolic ones (2.3). (2.3) has a simple gravity dual for all values of n, eq. (2.4), while for (2.2) we have worked only to leading order in (n − 1), (2.8). It is natural to suspect that there should be a simple gravity dual to (2.2) for all values of n, and a correspondingly simple generalisation of our results non-linearly in n. Such generalisation would be applicable, e.g., to negativity as the n → 1/2 limit [35]. We plan to investigate this elsewhere. It should also be possible to generalise the results of this paper to other bulk dimensions, and to other theories of gravity.
The detailed mechanism by which the bulk regulates the boundary cone is in agreement with [32]. There, this structure was derived by solving the finite part of Einstein's equations around the Ryu-Takayanagi surface (the infinite part gives that the surface is minimal [4]). That suggests that the addition of matter may change the detailed regulation, and therefore the solution of the splitting problem we presented in section 7. It may be interesting to explore this possibility, and whether it impacts the log divergence of Rényi entropy we discussed in section 6 -perhaps there is after all a gravity dual for which f b (n) = f c (n) is realised.