Holographic entanglement entropy of a 1 + 1 dimensional p-wave superconductor

We examine the behavior of entanglement entropy SAEE of a subsystem A in a fully backreacted holographic model of a 1 + 1 dimensional p wave superconductor across the phase transition. For a given temperature, the system goes to a superconducting phase beyond a critical value of the charge density. The entanglement entropy, considered as a function of the charge density at a given temperature, has a cusp at the critical point. In addition, we find that there are three different behaviors in the condensed phase, depending on the subsystem size. For a subsystem size l smaller than a critical size lc1, SAEE continues to increase as a function of the charge density as we cross the phase transition. When l lies between lc1 and another critical size lc2 the entanglement entropy displays a non-monotonic behavior, while for l > lc2 it decreases monotonically. At large charge densities SAEE appears to saturate. The non-monotonic behavior leads to a novel phase diagram for this system.


Introduction
The entanglement entropy S EE A of a subsystem A of a system is a useful non-local quantity in quantum field theories [1][2][3][4][5][6]. This quantity often provides interesting probes of the physics of phase transitions, most notably for quantum critical transitions where it typically diverges [7]. Even for thermal transitions, the entanglement entropy provides additional diagnostics of the nature of the phases, since it is sensitive to the number of available degrees of freedom. In recent years, the Ryu-Takayanagi prescription has become a key tool to compute entanglement entropy of strongly coupled systems which have gravity duals [8][9][10]. Indeed, starting with [11], there have been several studies of the change of behavior of the entanglement entropy in holographic models of critical phase transitions. [11] examined a model of s-wave holographic superconducting transition, [12] considered a model of an insulator-superfluid transition while [13] studies a model of holographic p-wave superconductor phase transition based on [14,15]. Recently, [16,17] extended the study of the holographic entanglement entropy near critical phase transitions. In all these models a suitable order parameter condenses in the superfluid phase, reducing the number of degrees of freedom. Indeed these studies found that the magnitude of the entropy in the condensed phase is always lower than what it would have been in the absence of condensation. In all these studies, the entanglement entropy as a function of the temperature or the chemical potential displayed a cusp at the location of the transition where the relevant derivative becomes discontinuous. As one goes into the condensed phase the entanglement entropy sometimes shows a non-monotonic behavior close to the transition, while far from the transition it keeps decreasing.
In this paper, we perform a calculation of the holographic entanglement entropy (HEE) in a toy model of a 1 + 1 dimensional p-wave superconductor at non-zero charge density and temperature. The dual theory is a 2 + 1 theory of gravity and SU(2) Yang-Mills field with a negative cosmological constant. The action is given by

JHEP09(2017)016
where A µ is a SU (2) gauge field with the field strengthF µν = ∂ µÃν − ∂ νÃµ − ig YM [Ã µ ,Ã ν ]. In this convention,Ã µ is dimensionless, while g YM and 1/κ 2 have the dimension 1. This model has been introduced in [18] following earlier work in higher dimensions [14]. An important well known aspect of this 1 + 1 dimensional model is that the gauge field has to be treated in alternative quantization [19]: the non-normalizable part of the field can be interpreted as the expectation value of a dynamical gauge field living on the boundary, which makes it morally closer to true superconductors, as opposed to superfluids. In the normal phase the dual geometry is a AdS 3 charged BTZ black hole, with the role of the usual Maxwell field played by the time component of the diagonal gauge field A (3) t . The boundary theory then has a nonzero temperature and charge density. As we increase the charge density ρ for a given temperature, there is a critical value ρ c beyond which the black hole acquires a vector hair, e.g. a spatial component A (1) x becomes nonzero: 1 its boundary value then becomes the expectation value of a vector order parameter in the dual field theory. In the probe approximation the corresponding bulk solutions have been obtained in [18].
The dual field theory lives on a circle with circumference 2πL (whereL is the AdS scale as in equation (1.1), and we calculate the holographic entanglement entropy of an interval of size l in this circle. In the normal phase the exact gravity solution is of course known. In the condensed phase the solution is known only in the probe approximation. To obtain the holographic entanglement entropy we first obtain the fully backreacted solution for the hairy charged BTZ black hole by numerically integrating the bulk equations of motion. Using this backreacted solution we compute the Bekenstein-Hawking entropy S BH = A 2κ 2 (where A is the horizon area and κ is the gravitational constant) and the entanglement entropy S EE using the Ryu-Takayanagi prescription where γ A is the length of the one-dimensional bulk geodesic whose end points coincide with the two end points of the interval on the boundary. We focus on the behavior of these quantities at some fixed temperature as a function of the charge density for various values of l. We find that the Bekenstein-hawking entropy S BH has a cusp at the critical charge density q = q c where ∂S BH ∂q is discontinuous. For q > q c this monotonically decreases, approaching a constant value. This decrease is essentially due to the reduction of degrees of freedom due to condensation.
For a given temperature and q < q c the entanglement entropy monotonically increases as a function of q and its derivative has a discontinuity at q = q c . This cusp-like behavior has been observed earlier in higher dimensional models. However the behavior for q > q c appears quite different from that in the higher dimensional examples. This depends on the subsystem size l and there are two critical sizes l c1 , l c2 where the behavior changes. For small enough size, l < l c1 the entanglement entropy keeps increasing monotonically. When l c2 < l < l c1 this quantity displays a non-monotonic behavior, decreasing at first to JHEP09(2017)016 a local minimum, and then increasing to a maximum value and finally decreasing again. The local minimum moves closer to ρ = ρ c as l approaches l c1 . At l = l c2 the local minima and maxima merge to a point of inflection, so that for l > l c2 the entanglement entropy decreases monotonically with increasing charge density.
While we do not have a complete understanding, it appears that this complex behavior arises from a competition between charge density and condensation. If there was no condensation, the entanglement entropy would have kept increasing with charge density for any subsystem size. Condensation reduces the number of effective degrees of freedom. When the subsystem size is small enough, the effects of condensation are small, and the tendency of the entanglement entropy to increase with charge density wins over. When the subsystem size is large enough, the effect of condensation wins over. The non-monotonic behavior seen for intermediate sizes reflects a competition between these two effects. This results in an interesting "phase diagram" where the sign of (∂S EE /∂q) l distinguishes different "phases".
In section (2) we describe the system and the associated critical transition, and calculate the backreaction by numerically solving the equations of motion. The backreacted solution is then used to calculate the Bekenstein-Hawking entropy. In section (3) we use the solution above to calculate the minimal geodesics and hence the holographic entanglement entropy. (4) contains a brief discussion.

The backreacted solution
The action of Einstein-Yang Mills system we consider in this paper is where the field strength of the SU(2) Yang-Mills term is given byF One of the spatial directions (called y below) is compact, y ∼ y + 2πL. In (2.1),Ã µ is dimensionless, while g YM and 1/κ 2 have the mass dimension 1. A standard field redefinitionÃ µ → A µ /g YM the field strength becomes while the Yang-Mills kinetic term becomes proportional to 1/(κ 2 g 2 YM ). In this section, we analyze the backreaction (whose strength is controlled by 1/g 2 YM ) of the SU(2) Yang-Mills term into the metric of the 2+1 dimensional gravity to compute the holographic entanglement entropy.
The equations of motion of the gauge field and the metric derived from the above action become where the energy momentum tensor is

JHEP09(2017)016
If we set the Yang-Mills field to zero, the solutions to Einstein equations are either pure AdS 3 or an uncharged BTZ black hole.
In the presence of a nontrivial gauge field which is assumed to depend only on the radial AdS direction, the Einstein equations (2.2) and (2.3) can be solved with an ansatz: where σ a (a = 1, 2, 3) are the Pauli matrices. The function f 2 (z) vanishes at z = z h , the location of the black hole horizon. The Einstein equations can be re-written as the following three equations: The last equation is a constraint. Using the ansatz the equations of motion for the Yang-Mills fields become Note that the equations of motion depend only on the product of AdS radiusL and g YM . This can be used to set the scaleL = 1.
In the uncondensed phase the y component of the gauge field w(z) vanishes. In this case, nonlinear terms are absent in the field strength and the system reduces to an Einstein-Maxwell system. The solution is then well known: it is a charged AdS 3 black hole considered in [19,21] and described in (2.6) with The black hole horizon is located at z = z 0 . The Hawking temperature is given by (2.9) Furthermore, we have chosen φ(z) to vanish at the horizon, as required by regularity of the solution in the euclidean domain. The regularized mass of the black hole, M 0 = (L/z 0 ) 2 , is inversely proportional to the squared horizon position and satisfies the BPS-like bound JHEP09(2017)016 [22] which is saturated at T H = 0. At the equality, the horizon position z h is equal to the value of z where f (z) has its extremal value.
As explained in [19] the presence of a logarithmic term implies that the gauge field has to be treated in alternative quantization, which requires us to fix boundary conditions which specify the charge density q. The chemical potential is then given by the constant part µ = −q log z 0 . The mass dimension of q and µ are both 1.
When the charge density is large enough there is another solution which has w(z) = 0, with a free energy lower than a charged AdS 3 black hole.
The behavior of this solution near the AdS boundary z = 0 is where once again q is the charge density of the boundary field theory while the chemical potential is given by µ ≡ −q log(z p ). 2 The parameter w c is now interpreted as the expectation value of a vector order parameter in the boundary theory while J w is the source for this order parameter. Both these parameters have mass dimension 1. h 0 and r 0 are constants. The black hole horizon is now located at z = z h as specified by the condition f 2 (z h ) = 0. Regularity at the horizon requires as usual φ(z h ) = 0. The fields can be now expanded near the horizon as follows: The coefficients (a 1 , b 1 , b 2 , c 1 , c 2 , d 2 ) are constants. In terms of these parameters the Hawking temperature is Substitution of (2.11) into equations of motion (2.6) and (2.7) yields four independent parameters (a 1 , b 1 , c 1 , z h ). Other parameters like b 2 , c 2 and d 2 are then represented by these 4 parameters and g YM : We now obtain a backreacted solution by numerically solving the equations of motion. Recall that the last equation of (2.6) is a constraint equation. This can be solved by choosing regularity conditions at the horizon. We then solve the first two Einstein equations JHEP09(2017)016 of (2.6) and two equations (2.7) numerically, starting from the horizon and proceeding to the AdS boundary.
To look for a superconducting phase we need to find solutions with a vanishing source in the boundary theory. This implies that at the AdS boundary, J w has to vanish. A nontrivial solution for w(z) then signifies spontaneous symmetry breaking in the boundary theory. This also breaks residual bulk U(1) gauge symmetry generated by A 3 µ spontaneously. The parameters w c , J w (= 0), h 0 are then specified by using 4 parameters (a 1 , b 1 , c 1 , z h ) of the horizon expansion.
To go further, consider the scaling symmetry in the equations of motion as follows: 14) The first symmetry can be used to fix z h = 1. Using second symmetry, the leading coefficient of h 2 (z) in (2.10) can be fixed to set (h 0 = 1): this yields the standard AdS 3 metric near the boundary. Under scaling symmetry (2.13), the parameters of the solution transform as follows: where we have used J w = 0 becomes zero in the last equality. Note that the scaling transformation of log(z) in (2.10) produces the shift log λ s in the chemical potential µ.
The scaling transformation is important to fix parameters like the temperature or the charge. We can not fix these parameters without changing z h . When we fix the inverse temperature to be β, we need the scaling transformation λ s = β| z h =1 /β and then the horizon position z h is changed into z h = β/(β| z h =1 ) by using the above transformation. To consider the system of the finite charge q 0 , instead, we need to perform the scaling transformation λ s = q 0 /q| z h =1 . The horizon position is then changed into z h = q| z h =1 /q 0 .
We find that with a vanishing J w a solution with non-zero w c exists only when the charge density q exceeds a critical value q c . This solution is a hairy black hole. As discussed above w c is the expectation value of a vector order parameter in the boundary field theory: the hairy black hole is then the gravity representation of a superconducting phase.
In figure 1 we plot the behavior of w c /q as a function of q c /q at a fixed temperature T H = 0.15 with varying g 2 YM . The critical charge q c becomes 21.7T H , 33.5T H , and 45.1T H for g 2 YM = 5 × 10 5 , 10, and 6, respectively. The condensate vanishes for q < q c . The point q = q c is a critical phase transition. Near the critical charge q c , the condensate behaves as w c ∼ 1.18q 1 − q c /q so that we get a mean field critical exponent.
The probe approximation corresponds to large g 2 YM : this is the situation when the gravity backreaction can be ignored. Our results clearly shows that backreaction decreases w c /q at large charge q, while near the critical point the results approach those of the probe approximation. This is expected since near the critical point, w c is small so that the backreaction is small as well even for finite g 2 YM .
JHEP09 (2017) In figure 2, we plot the Bekenstein-Hawking entropy of the hairy black hole S = (2π) 2 /(κ 2 z h ) in units of 2π/κ 2 (= c 2d /6), where c 2d = 12π/κ 2 in units ofL = 1 is the central charge of the CFT. In the uncondensed phase, the entropy grows like q 2 when q/(g YM T H ) 1 (the left figure) and grows like a linear function of T H when q/(g YM T H ) 1 (the right figure). We find a cusp at the critical phase transition where q = q c (= 33.5T H ).

Holographic entanglement entropy and the phase transition
In this section we will calculate the entanglement entropy of an interval on the boundary in the y direction of length l, using the Ryu-Takayanagi formula [8][9][10]. The arc on the boundary is chosen to be −l/2 ≤ y ≤ l/2. We then need to calculate the length of a geodesic with minimum length in the bulk metric on a constant time slice which joins the JHEP09(2017)016 two endpoints of the interval. 3 Here z(y) denotes a curve in the bulk 2+1 dimensional space-time on a constant time slice.
Since we are dealing with a static situation, the action is independent of the time slice. It is useful to consider the action I EE as the integral of a Lagrangian with y being considered as the time. Since translations of y is a symmetry, the corresponding Hamiltonian is conserved. Thus, the equation of motion which follows from this action becomes a first order differential equation for z(y): where z = z * is the turning point, i.e. the point where z (y) vanishes. The curve z(y) is assumed to be smooth everywhere. The length l of the interval can be now easily calculated We have assumed that the curve is symmetric under reflections about the turning pointthis gives the factor of 2 in front of l curve A solution to the equation (3.2) with specified end points is a minimal length geodesic joining those points. According to the Ryu-Takayanagi formula the holographic entanglement entropy is the on-shell action of I EE divided by the gravitational constant One can check that l → λ −1 s l under the rescaling (z, z h , z * , ) → λ −1 s (z, z h , z * , ). It is then convenient to introduce the scale invariant quantities lq or lT . Note that the holographic entanglement entropy is invariant under the scaling transformation varying (z, z h , z * , ). In equation (3.4) the boundary z = 0 has been replaced by a cutoff boundary z = (which is essentially a UV cutoff of the boundary field theory). This regulates the UV divergence of the entanglement entropy, as discussed below.
We are interested in calculating S EE A as a function of the charge q for a given temperature. As discussed above, for q < q c the relevant bulk geometry is a charged black hole, while for q > q c the background is given by a hairy black hole with a nonzero condensate w c , with the metric (2.5) in section 2. In both cases, the UV divergent (area law) term of I EE behaves like I EE ∼ −2 log( ). This motivates the definition of the finite part of the minimal length, I m,f in is independent of the cutoff. Under the scaling transformation discussed above, I m,f in changes only by the additive constant −2 log(λ s ), which will not affect our analysis qualitatively.
We now use the metric obtained in the previous section for both the condensed and uncondensed phase to calculate I m,f in . This involves evaluation of the integral (3.4), which we perform numerically.
In the uncondensed phase T > T c or q < q c , we find that the holographic entanglement entropy is proportional to the subsystem size l (a volume law) when the length l is large (lT H 1) exactly as in higher dimensions [25,26]. In the opposite limit T H l 1 and ql 1, I m,f in approaches the value in pure AdS 3 , namely, S EE A ∼ 4π/κ 2 ·log(l/a) = c/3·log(l/a). These behaviors in the extreme limits are of course what is expected. Interesting non-trivial behavior is expected for T or q in the intermediate range [27], particularly when they are close to their values at the critical superconducting transition. Consequently, there is another critical interval size l c2 . For l > l c2 the quantity I m,f in decreases monotonically for q > q c . For g 2 YM = 10 we get l c2 T H ∼ 0.22. At this value of l the minimum of I m,f in is pushed beyond the largest value of q for which we could do the numerics, while the plateau value becomes very close to the minimum value.
This complex behavior is possibly a result of two competing trends. First, an increasing charge density tends to increase the entanglement entropy. This is clear in the uncondensed phase for all values of l. However condensation beyond q = q c results in a depletion of the number of degrees of freedom, and therefore tend to decrease the entanglement entropy. For small enough interval sizes, the effect of condensation is not very pronounced, so that the first trend wins leading to a monotonically increasing S EE A . In the dual gravity picture the condensate is pronounced near the horizon -the minimal geodesic stays far from the horizon for small intervals. This is consistent with the fact that at length scales which are significantly smaller than the scale of symmetry breaking, effects of symmetry breaking are invisible. However as the interval size increases, the minimal geodesic goes deep into the bulk and the effect of the condensate becomes pronounced. The intermediate region l c1 ≤ l ≤ l c2 is possibly characterized by the regime where these two effects are of the same order and therefore compete. Finally, for l > l c2 the effect of condensation overcomes the effect of increasing charge density. The degrees of freedom keep decreasing as we increase the charge density resulting in a monotonic decrease of the S EE A . While we do not have a quantitative calculation to back the above scenario we can check its consistency by calculating I m,f in for different values of g 2 YM and therefore changing the strength of backreaction. As we will see, this turns out to give further insights into JHEP09(2017)016 the nature of the critical values of the critical sizes of l c1 and l c2 and the mechanism of the competition. We present the g 2 YM = 10 case first, followed by g 2 YM = 6 and g 2 YM = 50. The smaller the value, the stronger the back-reaction and the effects of the competition mentioned above become more pronounced.
The entanglement entropy S EE A , normalized as I m,f in (I, q)/I m,f in (I, q c ), for g 2 YM = 10 with 6 different sub-system sizes is presented in figure 5, so that all the S EE A coincide at the critical point q = q c . In the left figure of figure 5, we show the behavior of entanglement entropy as we increase the sub-system size from the top plot to the bottom plot. The right figure is the slope of the entanglement entropy with respect to the charge, dI m,f in (I, q)/dq. In the range q > q c , the slope is positive for lT H = 0.14 and thus entanglement entropy increases monotonically, while the slope is negative for lT H = 0.23 with monotonically decreasing entanglement entropy. The entanglement entropy for 0.21 lT H 0.22 have non-monotonic behaviors revealing two zeros of its derivative. This is shown in the right panel of figure 6. As we increase the charge, entanglement entropy decreases, hits a minimum (first zero of the derivative), increases, arrives at a maximum (second zero of the derivative), and eventually reaches a plateau. Typically, the second zero corresponding to a local maximum happens at a very large value of q/q c . Finally, we clearly see that the minimum of the entanglement entropy moves to q = q c for lT H = 0. 15 As we increase the strength of the back reaction, decreasing g 2 YM , both critical lengths l c1 T H and l c2 T H decrease while the invariant distance between them, ∆lT H = l c2 T H − l c1 T H , increases.

A phase diagram
As discussed above, the non-monotonic behavior of the entanglement entropy is a result of a competition between two effects: the effect of charge density and the effect of condensation. It is useful to re-cast this in terms of two length scales which emerge in the problem. Let us recapitulate the salient features of our results in Fig 3 -Fig 8. These are all plots of the entanglement entropy as a function of the charge density, the different curves being results for different values of (lT H ). For a given value of g YM we have the following: • For l > l c2 the entanglement entropy decreases monotonically as a function of q/q c • At l = l c2 there is a point of inflextion at q = q • For l c2 > l > l c1 the entanglement entropy has a minimum at some value q min (l) and a maximum at a larger value q max (l).
• The minimum q min (l) increases as a function of l. At l = l c1 one has q min (l c1 ) = q c . On the other hand, q max (l) decreases with increasing l, reaching its minimum value q at l = l c2 . Thus at l = l c2 the maximum and the minimum merge into a point of inflection.
• At l = l c1 the minimum is at q = q c • For l < l c1 the entanglement entropy increases as a function of q at least upto q ∼ 5q c which is the maximum value of q used in our calculation.
It is then clear that the minimum of the entanglement entropy as a function of the charge can occur only for q < q , while the maximum of the entanglement entropy as a function of the charge can occur only for q > q .
To rephrase the behavior in terms of length scales it is useful to plot the quantity (∂I m,f in /∂q) l as a function of lT H . This is shown in figure (9).
We now define two length scales L c , L q , which depend on the charge density, as follows.
(∂I m,f in /∂q) l=Lc = 0, (∂ 2 I m,f in /∂q 2 ) l=Lc > 0 q < q (∂I m,f in /∂q) l=Lq = 0, (∂ 2 I m,f in /∂q 2 ) l=Lq < 0 q > q  Figure 9. ∂I m,f in /∂q as a function of the sub-system size l for g YM = 10 and for various values of q/q c . The left panel shows the results for q < q , where q is in between q/q c = 2 and q/q c = 3. The zero of ∂I m,f in /∂q in this regime of charge is a minimum (as a function of q), and the length scale L c is defined as the value of l at which ∂I m,f in /∂q = 0. Clearly, L c increases as the charge density increases. The right panel shows the results for q > q . Now the zero of ∂I m,f in /∂q is a maximum (as a function of q). The length scale L q is defined as the value of l at which ∂I m,f in /∂q = 0. Clearly L q decreases with increasing q.  Figure 10. Phase space (lT H , q/q c ). The two curves L c (q/q c ) and L q (q/q c ) meet at q = q . In the tent like region below these curves the derivative (∂I m,f in /∂q) l > 0, outside the tent this is negative.
From figure (9) it is clear that L c increases with the charge density while L q decreases with the charge density. At q = q one has L c (q ) = L q (q ) = l c2 , which is the maximum possible value of either of L c , L q .
These two length scales and their dependence on the charge density can be used to chart out a "phase diagram" for the entanglement entropy. Figure (10) shows the general nature of the phase diagram, parameterized by the length of a subsystem (vertical axis) and the charge density q/q c (horizontal axis). The solid lines denote the functions L c (q/q c ) and L q (q/q c ) -as discussed above the former quantity exists for q < q while the latter JHEP09(2017)016 of the quench rate [33]. The setup used in this paper should be useful in extending the discussion to entanglement entropy.
Finally, it will be interesting to explore in detail the stringy embedding of the 3d toy model (2.1). In the probe limit of flavor branes, the toy model (2.1) corresponds to a D3 − D3 system of type IIB string theory [18] without backreaction. In this limit the dilaton is a constant, which is why it is omitted in (2.1). In the presence of backreaction, the dilaton runs, making the problem much more difficult. It would be interesting to see if the strategy of [34,35] which expresses the leading order correction to the entanglement entropy in terms of the energy momentum tensors of the probe brane and the minimal surface and the bulk graviton propagator can be used to calculate this. Such a calculation would also need to include the backreaction of the four form gauge field sourced by the D3 branes. We hope to pursue this in the near future.