Holographic Subregion Complexity for Singular Surfaces

Recently holographic prescriptions are proposed to compute quantum complexity of a given state in the boundary theory. A specific proposal known as `holographic subregion complexity' is supposed to calculate the the complexity of a reduced density matrix corresponding to a static subregion. We study different families of singular subregions in the dual field theory and find the divergence structure and universal terms of holographic subregion complexity for these singular surfaces. We find that there are new universal terms, logarithmic in the UV cutoff, due to the singularities of a family of surfaces including a kink in (2+1)-dimension and cones in even dimensional field theories. We find examples of new divergent terms such as square logarithm and negative powers times the logarithm of the UV cut-off parameter.


Introduction
Quantum entanglement has been widely studied in the context of holographic field theories after the pioneering Ryu-Takayanagi (RT) proposal [1,2]. Quantum complexity is another notion in quantum information theory which has been recently included in the context of holographic field theories. Roughly speaking, quantum complexity of a state is the minimum number of information gates needed to prepare a state from a given reference state. There exist some efforts to develop a holographic dual for quantities related to this notion in the context of AdS/CFT correspondence [3][4][5][6][7][8][9][10][11].
From a more geometrical point of view, it is well-established that the von Neumann entropy of a subregion in a given state corresponds to the area of a co-dimension two surface in the gravity solution dual of the state. People have also tried to find geometrical duals for other quantities in the context of information theory; such as Renyi entropies [12,13], information metric (fidelity susceptibility) [14,16,17] 1 , fisher information [18], etc.. Some of these geometrical objects are still co-dimension two objects in the dual theory but some are not.
There are two distinct proposals to compute complexity of a state in the dual gravity theory. The first one, which is sometimes called the 'complexity=volume' proposal, states that the complexity of a given state at a given time in the boundary theory is given by the volume of an extremal codimension one surface in the bulk which meets the corresponding time slice. To be more concrete, one can state this proposal as where the maximum is chosen among those co-dimension one surfaces which end on the corresponding time slice on the conformal boundary. In this proposal ℓ is some length scale which should be identified case by case, e.g. the radius of the asymptotically AdS solution or the radius of the horizon in case of AdS black-hole geometries. This non recognized length scale seems to be a disadvantage of this proposal.
The other proposal, which is sometimes called 'complexity=action', states that the complexity of a given state at a given time is equal to the on-shell action of the dual (Einstein) gravity theory computed in the domain of dependence of any Cauchy surface in the bulk which ends on the given time slice at the conformal boundary. 2 This region is known as the Wheeler-DeWitt patch corresponding to the given boundary time slice. Although this proposal (in contrast with the previous one) does not need any length scale by definition, it has its own challenges due to surface terms and corner contributions of the Wheeler-DeWitt patch (see [11,20]). We will come back to this point in the next section.
A natural generalization of the 'complexity=volume' proposal concerns with generic mixed states. A specific way of constructing a mixed state out of the entire state of a system is to trace out a part of the space-like manifold of the dual field theory. The mixed state constructed in this way is described by what is well-known as the reduced density matrix. Then the complexity of such a (static) state is proposed to be given by the volume enclosed by the Ryu-Takayanagi surface and the corresponding subregion in the boundary theory. 3 To be more concrete the subregion complexity is defined as [16] C subregion = max V (γ) 8πℓG N , (1.2) where γ is the RT surface of the corresponding subregion and ℓ is a length scale of the dual geometry.
The maximization is among volumes enclosed by surfaces ending on the same subregion. This proposal (up to a numerical factor) reduces to 'complexity=volume' given in (1.1) if the subregion is chosen to be the whole time slice of the dual theory.
Different proposals for complexity all lead to UV divergent results since they all contain a volume of a surface which reaches the conformal boundary of an asymptotically AdS geometry. This is the same as what happened in the case of holographic entanglement entropy. Natural questions about such quantities are: "What is the divergent structure of this quantity?", "How it can be regularized?", "What kind of universal information can be extracted from it?", and "Is it possible to find any monotonic function out of this quantity under the RG flow of the dual theory?". Specifically for the case of subregion complexity one may also ask about the (subregion) shape dependence of the divergence structure.
Some of the above questions has been recently addressed for different proposals of complexity and even for complexity of reduced states due to smooth subregions [11]. The goal of this paper is to investigate the divergence structure of subregion complexity when the subregion is a singular surface. Similar to the case of entanglement entropy we expect new divergent (sometimes new universal) terms due to singularities in the subregion. There has been done a considerable amount of efforts to investigate the role of singularities of entangling regions in the context of (mostly holographic) entanglement entropy [21][22][23][24][25][26][27][28][29][30][31][32][33][34]. We will consider the simplest case of a singular surface in a (2 + 1)-dimensional field theory and its generalizations to enough symmetric singular surfaces in higher dimensions (see [24,25] for a similar analysis for entanglement entropy) and study the divergent structure due to subregion complexity proposal [16].
The rest of this paper is organized as follows: in section 2 we define different families of singular surfaces which we study. If the reader is just interested in the final results, we have summarized our subsequent results in this section. In the following sections we study complexity of different subregions and we finalize in the last section with addressing interesting directions for future studies.

Singular Subregions and Summary of Results
We are interested in asymptotically AdS solutions of Einstein gravity with a negative cosmological constant in d + 1 dimensions. The simplest case which we study in this section is the pure AdS d+1 solution in the Poincare patch with the following coordinates where z is the radial coordinate and L is the AdS radius. Here dθ 2 + sin 2 θdΩ 2 n is the metric on a unit sphere S n+1 and the term m i=1 (dx i ) 2 indicates a flat R m space in Cartesian coordinates. The conformal boundary of this solution is achieved in the z → 0 limit. Hence, the boundary metric reads For the whole manifold of the bulk, as well as the boundary, the range of the parameter θ is (−π, π) for n = 0 and (0, π) for n > 0. However, throughout this paper we consider different kinds of singular subregions, i.e. the conic singular subregions, in which −Ω < θ < Ω for n = 0 and 0 < θ < Ω for n > 0. The simplest conical geometry is a kink (k) in d = 3 where n = m = 0, as the following subregion of the boundary We also study singular surfaces in asymptotically AdS d+1 geometries given by In these cases f 1 and f 2 are functions which are determined by the gravity equations of motion. We study different cones and creases in these asymptotically AdS geometries.
In Ref. [24] the holographic entanglement entropy for the above singular surfaces is calculated in Einstein gravity and also some specific higher derivative gravity theories. In this paper we calculate the holographic complexity in each case by using the proposal of Ref. [16]. As we have mentioned in the previous section, according to this proposal the volume of a co-dimension one surface enclosed by the subregion in the boundary theory and the RT co-dimension two surface in the bulk is proportional to the complexity of the (mixed) state corresponding to the subregion. To do so, one should find the RT surface corresponding to subregion A which we denote by γ A and calculate the volume V (γ A ) enclosed by γ A . The holographic complexity is proposed to be given by Eq. (1.2) [16]. We choose ℓ in the asymptotically AdS gravity solutions to be identified with the AdS radius. In what follows we will study this quantity in different singular subregions.

Summary of Results
Since the detailed calculations presented in next sections may be involved, here we briefly summarize our results. We study the divergent structure of holographic subregion complexity and find new divergences due to singular subregions which in some cases lead to new universal terms.
In the case of a crease entangling region in a (2+1)-dimensional boundary theory (see the left panel of Fig. 2) we find that there is a new divergent term of the form log δ which is a universal term. The entanglement entropy for the same subregion also leads to a logarithmic universal term.
For the case of a crease entangling region with a flat locus, which we denote by k × R m (see the right panel of Fig. 2) there is no universal term due to the singularity and even no actual new divergent term, although the subleading divergent term gets corrections from the singularity. This resembles to the entanglement entropy in having no new universal term. Even for the case of k × S 1 , which again the locus of the singularity is flat, there is no new universal term and no new divergent contribution from the singularity.
In the case of creases with a curved locus we again find that there is no new divergent term.
This is in contrast with what happens for entanglement entropy of these surfaces. We study the case of k × S 2 and k × S 3 and also k × R × S 2 and in all of them although there is a log δ term but it is suppressed with a positive power of δ resulting in no new divergent term.
The most interesting behavior happens for conical subregions which we show by c n (see the middle panel of Fig. 2). For these subregions we find that there is new universal log δ term for odd n and log 2 δ for even n's. We have worked out a few examples of this for n = 1, 2, 3, 4. In comparison with entanglement entropy of these surfaces we find a shift from odd to even n's where log 2 δ and log δ appear respectively. It would be very interesting to find out whether these universal terms are related to some characteristic feature of the dual field theory.
The other family of singular surfaces which we have studied are conical creases of the form c n × R m and c n × S m . Among these surfaces the only case which we find that a universal log δ term appears is c 1 × S 2 . In other cases new divergent terms appear due to the singularity which have the form of 1/δ log δ or 1/δ 2 log δ. These are very similar to what has been recently found from the 'complexity=action' proposal [11]. This similarity may be due to the singularities within the Wheeler-DeWitt patch. We have summarized our results in the above table.
3 Flat Locus Singular Surfaces

Kink k
The simplest case is a kink in a 2+1 dimensional boundary theory. The bulk metric dual to the vacuum state is given by IR cut-off. The corresponding Ryu-Takayanagi surface can be described by z = z(ρ, θ), hence the entanglement entropy is given by Since there is no length scale except ρ, the radial coordinate z depends on ρ linearly [23], i.e.
and h(θ) should be found such that it minimizes the entropy (area) functional and is anchored to the kink in the asymptotic boundary. Applying this into Eq. (3.2) gives does not depend on θ explicitly, we have the following conserved quantity along θ translation To find the holographic subregion complexity we should write the volume V (γ) of the subregion of the bulk where ǫ is a short distance cut-off in the boundary corresponding to δ in the bulk. To clarify the singular terms of Eq. (3.6) we convert θ integration to an integral over h as follows One can easily find the following expression from Eq. (3.5) Using the coordinate transformation y = 1 , where y → ∞ as we approach the boundary via θ → Ω, we have . (3.9) In the limit δ → 0 and hence y → ∞ the integrand is finite. So we can find it just for y → ∞. We have finally where α(h 0 ) is the cut-off independent term given by , (3.11) which vanishes in the smooth region limit (i.e. Ω → π). Thus the divergent structure of holographic complexity of kink is given by

Cone c n
As indicated in the previous section, to consider a conical subregion c n with n = d − 3, we use the following form of the bulk metric where dΩ n is the metric of a unit sphere S n . The subregion in the boundary is defined by ρ ∈ [0, H] and θ ∈ [0, Ω]. The extension of this region in the bulk is denoted by the function z(ρ, θ). One should find the profile of this extension via minimizing the following area functional where Ω n is the volume of the unit n-sphere andż = ∂ θ z, z ′ = ∂ ρ z.
As in the previous case, z can depend on ρ only linearly, i.e. z(ρ, θ) = ρ h(θ). Using this assumption, and change of variable y = sin θ = y(h) which giveṡ the equation of motion for the case d = 4 read as follows where y ′ = dy dh and y ′′ = d 2 y dh 2 . Since we are interested in the singular behavior of the complexity near the boundary, where h → 0, let us concentrate on this limit (still for d = 4). For this reason we consider a power law expansion for y(h) in terms of h and put it in Eq. (3.15). Then using the boundary condition y(0) = sin Ω we find the following result The expansion forḣ follows consequently fromḣ = The corresponding volume is given by Using asymptotic expansions (3.16) and (3.17) the integrand of (3.18) has the following behavior near the boundary Let us divide singular parts of V (γ) into I 1 and L 2 where the latter contains the singularities due to the integrand while the former shows the contribution of the limits of the integrations, i.e. So the singular part of the complexity is given by In the limit h → δ/ρ there is no singular term from integration over h (neither from the integrand nor from the integration limits); we have just a logarithmic singularity from the lower limit of the integration over ρ as follows  The singular terms in I 2 can be calculated directly. Hence we have (3.25) One can perform similar computations for cones in higher dimensions. We have done this for c 2 and c 3 in CFT 5 and CFT 6 respectively. The method is similar to what we have presented in d = 4, so we will skip the details and report the results in these cases.
In the case of c 2 one finds two family of divergent terms proportional to log δ and log 2 δ which are given by One should note that the C log 2 5,c 2 is not a universal term. For the case c 3 we find

Crease k × R m
Consider the following metric for a AdS d+1 space-time in the bulk (3.28) The area functional to be minimized is given by Again one can use the scaling property z = ρh(θ), to find the equation of motion as Eq. (3.30) can be integrated over to find the following constant along the θ variation Noticing that h is a decreasing function near the boundary, we have from Eq. (3.31) One can find the volume as In the limit h → 0 the integrand in the last term behaves as So we can write

(3.34)
We can separate the divergent term as follows Let us denote it is clear from Eq. (3.32) that J(h) ∼ (h 2 ) as h → 0. We can find the integral (I 1 ) by parts (3.37) Now for finding the divergences of I 2 , we make a change of variable from ρ to q = δ ρ and then Taylor expand the terms around δ = 0 (3.38) From Eq. (3.32) q d−3 J(q) ∼ q d−1 for small q, hence in the above expression the integral over q is finite. We have also (3.39) So the singular terms of the volume is as follows The complexity is finally given by where g [ab] (S n ) is the metric of the sphere (S n ). The surface function to be extrimized is the following S d,cn×R m =L d−1H m Ω n dρdθ ρ n sin n (θ) ż 2 + ρ 2 (1 + z ′2 ) z d−1 . (3.43) The equation of motion for z(ρ, θ) after imposing the scaling relation z = ρh(θ) reads h(1 + h 2 )ḧ + n cot(θ)hḣ 3 + (d + nh 2 −1)ḣ 2 + n cot(θ)h(1 + h 2 )ḣ  Now we can use it to make the h integral in holographic complexity finite, i.e.  Let us indicate the integrand in I 1 by J 5 (h) and integrate it by parts  For c 1 × R 2 the result is as follows For c 2 × R 1 similar steps leads to (3.57)

Curved Locus Singular Surfaces
In this section, we consider several singular embeddings which have curved locus such as k × Σ and c n × Σ , where locus Σ will take the form S m or S m−p × R p .

Crease k × Σ
Consider the geometries k × S 2 , k × S 3 and k × R × S 2 . We will see that singularities with even dimensional locus will contribute through a logarithmic term. To begin with, let us consider d = 5 CFT on background R 3 × S 2 . The action for six-dimensional dual Einstein gravity reads We consider the following ansatz for the solution, where dΩ 2 2 = dξ 2 0 + sin 2 (ξ 0 )dξ 2 1 represents a two-sphere metric and f 1 and f 2 are functions of the radial coordinate. The boundary of this solution is R 3 × S 2 with R 1 the radius of S 2 ; so we can recover the flat boundary results in the limit R 1 → ∞. Using the Fefferman-Graham expansion near the boundary to find f 1 and f 2 leads to The subregion of interest here is ρ ∈ [0, H] and θ ∈ [−Ω, Ω] where H is again a IR cut-off. The coordinates are (z, θ, ξ 0 , ξ 1 ) on the minimal surface and ρ = ρ(z, θ) on the sphere. In the limit R 1 → ∞ one may expect from the case of entanglement entropy that leading order correction to the holographic subregion complexity would be O(1/R 2 1 ), however, we show that in this case there is no new divergent term up to O(1/R 4 1 ). We first work out the solution ρ(z, θ) in this approximation with the following ansatz Using the ansatz (4.5) in the equation of motion of ρ(z, θ) leads to vanishing of even terms g 2n . In order to separate the logarithmic divergence, we impose ρ = ρ 0 (z, θ) + ρ 1 (z, θ)/R 2 1 + ρ 2 (z, θ)/R 4 1 , where ρ 0 = z/h(θ) and ρ 1 = z 3 g 3 (θ), ρ 2 = z 5 g 5 (θ) are higher corrections in the large R 1 regime. Now we come back to the metric (4.2) and find the volume holographic complexity as Now, we can insert the ansatz ρ = ρ 0 + ρ 1 /R 2 1 + ρ 2 /R 4 1 that ρ 0 = z/h(θ), ρ 1 = z 3 g 3 (θ) and ρ 2 = z 5 g 5 (θ) and use (4.4) in the integrand to simplify the results as and where δ is the UV cut-of. We have also changed the integration limits from (−Ω, Ω) to (0, Ω) and then changed the integration variable in V 1 to h(θ). It is instructive to use the following constant of motion which is related to h(0) at the turning point. To find the logarithmic divergent parts it is enough to find the asymptotic behavior of h and g 3 . Solving g 3 in terms of h in the limit of small h leads to 10) where b 3 can be fixed by demanding g 3 to have an extremum at θ = 0. We will need to find the series expansion of h 1c in terms of δ as follows where h 1c = h(Ω − ǫ). The result is obtained for the leading corrections in R 1 at any order of δ [?].
Now we look at (4.7) to analyze the divergent terms in the asymptotic limit and (4.14) Now we differentiate each of them with respect to the UV cut-off and look for 1/δ divergent terms.
One can easily find So from (4.17)-(4.27) we can find the logarithmic divergences in the holographic complexity for k × S 2 geometry as follows Note that in this case no new divergent term appears due to the singular surface. All new log δ terms are suppressed with a factor of δ α where α ≥ 1.

Subregion k × S 3
Now we want to find the holographic subregion complexity for k × S 3 geometry in a CFT on R 3 ×S 3 . We will show that in this case the singularity gives no logarithmic contribution to subregion complexity. Consider the following metric where dΩ 2 3 = dξ 2 0 + sin 2 (ξ 0 )dξ 2 1 + sin 2 (ξ 0 ) sin 2 (ξ 1 )dξ 2 2 is the unit S 3 and we find f 1 and f 2 as Similar to the previous case the induced coordinates on the RT surface are (z, θ, ξ 0 , ξ 1 , ξ 2 ) and ρ = ρ(z, θ). Using the equation of motion for h we can find the following constant of motion which can be fixed in terms of the boundary data. Using the metric (4.35) we find the holographic complexity as Inserting the ansatz ρ = ρ 0 + ρ 1 /R 2 1 , ρ 0 = z/h(θ)and ρ 1 = z 3 g 3 (θ) and using the expansions (4.36) in the integrand, simplifies the result as and where δ is the UV cut-of, such that ρ(z, Ω − ǫ) = H and z m is defined such that ρ(z m , 0) = H. We have also changed the integration limits from (−Ω, Ω) to (0, Ω) and then changed the integration variable in V 1 to h(θ).
Similar to what we have done in the previous sections in details, one can work out the logarithmic divergence in this case. Here we step the details and report to the final result (4.41) In this case no new divergent term appears due to the singular surface and all new log δ terms are suppressed with a factor of δ α where α ≥ 1.
In the following we give another example showing that odd dimensional locus does not contribute to logarithmic singularities, althogh it has non-zero curvature. We consider a CFT defined on R 4 × S 2 .
The bulk metric is given by where dΩ 2 is the line element over S 2 and f 1 and f 2 have the following expansions Returning to the metric (4.42) we find the holographic complexity as We then insert the ansatz ρ = ρ 0 + ρ 1 /R 2 1 , ρ 0 = z/h(θ) and ρ 1 = z 3 g 3 (θ) and use (4.43) in the integrand to simplify the expressions as follows and Again we step the details of the rest of this calculation we find (4.48) As the case of k × S 2 and k × S 3 new logarithmic divergent term in this case are also suppressed with a factor of δ α with a positive power.

Conical Crease c n × Σ
In this section, we will calculate holographic complexity for subregions with conical singularities of the form c n × S m .
Subregion c 1 × S 1 To begin with, we concider the simplest case with m = 1. In this case, the background geometry for CFT is R 4 × S 1 . The dual bulk geometry is then given by ds 2 = L 2 z 2 dz 2 + f 1 (z)(dt 2 + dρ 2 + ρ 2 dθ 2 + ρ 2 sin 2 (θ)dφ 2 ) + f 2 (z)R 2 1 dξ 2 0 ,  c n in even dimensional field theories a new universal log δ terms appears. In odd dimensional field theories the singularity of a cone c n gives rise to a log 2 δ divergent term. We also showed that surprisingly crease singularities of any type do not give rise to any universal term or even any new divergent term. For generalized conical singularities the situation is completely different. There are examples which new power law divergences appear but there is no new universal term due to the singularity. We found also an example, i.e. c 1 × S 2 , with a curved locus that has a new universal term. Another type of conical singularity has 1 δ log δ and 1 δ 2 log δ divergent terms for even and odd dual field theories respectively. The latter family is very similar to what has been recently found using 'complexity=action' proposal on the Wheeler-DeWitt patch which also posses corners. We have summarized all of these results in a table in section 2.
There are several directions to follow in future works. Regarding the divergence structure of subregion complexity, the most important question is whether one can define any monotonic function from the universal terms which leads to a kind of 'c-function' in higher odd-dimensional dual field

theories?
Another interesting open question is how to generalize complexity proposals beyond Einstein gravity. Recently there have been some proposals trying to address this question (see e.g. [35]).
A natural question about this work is how to study the role of singularities of subregions in the 'complexity=action' proposal. Recently some progress have been made in [11] for spherical subregions. The authors have proposed the intersection between the "entanglement wedge" and the corresponding WDW patch for 'complexity=action' for mixed states constructed from subregions.
It would be instructive to understand this proposal by considering more complicated examples.