More on phase transition and Rényi entropy

In this paper, we study the scalar field condensation around the hyperbolic black hole solutions in the Einstein and Gauss–Bonnet gravities. We investigate the Rényi entropy and inequalities governing on it under this phase transition. Our numerical computations show that for the positive values of the Gauss–Bonnet coupling and below a critical temperature one of these inequalities is violated. This puts more restrictions on the allowed values of the Gauss–Bonnet coupling.

Generally speaking, when a quantum system is in a pure state, EE is a good measure for the degree of entanglement between any two subsystems. If the composite system is in a mixed state, other metrics such as the mutual information have to be used. To define EE more carefully, consider a system in a pure quantum state, composed of two subsystems A and its complement A c , whose state can be described by a density matrix ρ. The entangling surface ∂ A is defined as a boundary surface of the spatial region A in the quantum field theory (QFT) under consideration.
By tracing over the degrees of freedom in A c , one can construct a reduced density matrix ρ A = tr A c ρ. The EE between these two subsystems is measured by the Von Neumann entropy of the reduced density matrix i.e. S E E = −trρ A log ρ A . The computation of log ρ A is a hard task even in the simplest cases of two dimensional QFTs. Instead, one a e-mail: s.qolibikloo@mail.um.ac.ir b e-mail: a-ghodsi@ferdowsi.um.ac.ir can compute trρ n A and use it for another measure, known as Entanglement Rényi Entropy (ERE) [16,17] (1.1) Taking the limit n → 1 one yields the entanglement entropy, i.e. S E E = lim n→1 S n . ERE contains more useful information about the spectrum of ρ A than EE. Indeed, one can recover the whole spectrum of the reduced density matrix ρ A only by knowing the Rényi entropies for all integers n > 0 [18]. For example, S ∞ = − log λ max , where λ max is the largest eigenvalue of ρ A or S 0 = log D, in which, D is the number of non-vanishing eigenvalues of ρ A . In particular, S 2 = − log trρ 2 A = − ln P, where P is the probability of finding two systems in the same state, after the measurement in a diagonalized basis. These special cases are called the min-entropy, the Hartley or the max-entropy and the collision entropy respectively.
As we mentioned, the computation of trρ n A is much easier than the computation of tr log ρ A . In the context of Conformal Field Theories (CFT)s, a standard way to calculate trρ n A (and consequently S E E ) is the "Replica Trick". In this method, after a Wick rotation from the flat Minkowski metric where the CFT lives to a metric with Euclidean signature, the desired trρ n A operator is given in terms of the path integral on an nsheeted Riemann surface where Z is the partition function of the original space-time and Z n is the partition function on a singular space, which is constructed by gluing n copies of the original space along the boundary ∂ A. At the end, trρ n A essentially becomes product of the two-point correlation functions of the twisted vertex operators [2,19] which is very hard to calculate, except in very specific cases, for example see [20][21][22] or [8] and references therein. In any case, after the computation of trρ n A , we can substitute it in the following equation to compute S E E , . (1.3) In the context of holography, EE first appeared in [6] where the authors, Ryu and Takayanagi (RT), gave a simple yet ingenious way for the computation of the S E E . In their method, the EE between the two spatial regions A and its complement in the d-dimensional boundary CFT, is found to be proportional to the extremum of the surface area of the bulk hypersurfaces , which are homologous to the region A, i.e. ∂ A = ∂ .
The resemblance between the formula above and the Bekenstein-Hawking formula for black hole (BH) thermal entropy S B H = A 4G , is extraordinary. However, the surface does not need to coincide with any BH event horizon in general. This striking similarity prompted physicists to look for a derivation of the RT conjecture, for example, see [23].
In paper [24], authors give a novel method for the computation of the entanglement entropy of a Spherical Entangling Surfaces (SES) in a CFT, on a flat Minkowski space-time. In what follows we present a quick review of their approach.
First, let us consider a quantum system in a d-dimensional flat Minkowski space-time R 1,d−1 , which is divided into two subsystems, A and its complement A c . The subsystem A consists of a spatial region inside a ball of radius R, whose boundary (at time t = 0) is denoted by ∂ A. The subsystem A can be described by a reduced density matrix ρ A = tr A c ρ, where ρ is the density matrix of the vacuum state of the whole system. The Cauchy development D(A), by definition is the set of all space-time events p ∈ R 1,d−1 , through which, every non-space-like curve intersects with the region A at least once. After a Wick rotation the metric is where t E is the Euclidean time, r is the radial coordinate on the constant time slice (t E = 0) and d 2 d−2 is the metric on a unit (d − 2)-sphere. Notice that, SES in this metric is the region (t E , r ) = (0, R). After a series of conformal transformations (first z = r + it E , second exp(−w) = R−z R+z and finally w = u + iτ E R ) the metric can be written as where d 2 d−1 is the metric on a unit (d − 1)-dimensional hyperbolic plane and is a conformal factor = 2R 2 |R 2 −z 2 | = |1 + cosh w|.
In short, D(A) is conformally mapped to S 1 × H d−1 . In other words after the inverse Wick rotation, we will see that the causal development of the ball enclosed by the entangling surface (SES) is conformally mapped into a hyperbolic cylinder R × H d−1 , or equivalently the minimal surface of the RT conjecture gets mapped to the horizon of the topological BH. This mapping also translates the vacuum of the CFT in the original Minkowski space-time into a thermal bath with the temperature T 0 = 1 2π R in the hyperbolic cylinder. Note that, the curvature scale on the hyperbolic spatial slice, is equal to the radius R of the original SES. Just like any other operator in a CFT, the density matrix in the new spacetime R 1 × H d−1 , can be achieved by a unitarity transformation of the density matrix in the original geometry R 1,d−1 i.e. ρ therm = Uρ A U −1 . More explicitly, we may write the density matrix on D(A) as By considering the n copies of ρ A and after taking its trace Now by using the definition of free energy, F(T ) = −T log(Z (T )), we can write the Rényi entropy as follow Finally by utilizing the thermodynamical identity S = −∂ F/∂ T , one can rewrite the Rényi entropy of the above equation to (1.10) Here, S n is the desired entanglement Rényi entropy between the two subsystems A and its complement in the vacuum of the CFT, and S therm (T ) denotes the thermal entropy of the CFT. Upon taking the n → 1 limit of the above formula, S E E = lim n→1 S n = S therm (T 0 ), which means that the conformal transformations and the corresponding unitarity transformation of the density matrix discussed above, relate the EE of the SES to the thermal entropy of the same CFT at a temperature T 0 in R 1 × H d−1 space. This insight might not be particularly useful in the computation of the EE for a generic CFT unless using the AdS/CFT correspondence. We can relate the thermal bath in the boundary to a topological black hole in the AdS bulk space whose event horizon has a hyperbolic cross-section [25][26][27][28][29][30]. According to the AdS/CFT dictionary The right hand side of this equation is easy to compute by using the Wald's formula for entropy [31,32] in any general gravitational theory.
The authors of [33], use this method to find the Rényi entropies of the dual CFTs by finding the thermal entropy of black holes in the Einstein, Gauss-Bonnet (GB) and quasitopological gravitational theories. Moreover, the authors of [34,35] expand the aforementioned method to compute the ERE for the grand canonical ensembles. For these cases, one requires to consider a U (1) charged hyperbolic black hole. The computation of the charged Rényi entropies in [36,37] shows that the Rényi entropy inequalities put a new restriction on the allowed values of the coupling constants of the GB gravitational theory. The Rényi entropies obey the following inequalities [38] As discussed in [33], while the second and the third inequalities are yielded from the positivity of the black hole thermal entropy the first and the last inequalities are correct as long as the black hole has a positive specific heat. For a recent study of these inequalities see [39].
In [40] it was shown that adding a scalar field to the Einstein gravitational action creates a "hairy" black hole solution below a critical temperature and this encourages one to study the Rényi entropies of the dual CFT at this phase transition. It was also shown that at this critical temperature where the scalar field condenses, the second derivative of the Rényi entropy becomes discontinuous. Similar phase transitions occur in the study of the holographic superconductors [41], where it might cause a phase transition in the dual boundary CFT. For a review on holographic superconductors see [42].
Inspired by the works of [37,40] we are going to add a scalar field to the GB gravity to study the phase transition in the presence of the new gravitational coupling. We also find new restrictions on the allowed region of the GB coupling constant by studying the Rényi entropy inequalities.
The organization of the paper is as follow: In Sect. 2 we introduce the Einstein and the hairy black holes and compute their physical quantities such as the temperature, energy, and thermal entropy and also discuss the condensation of the scalar field. We review the related holographic Rényi entropy and its behavior under the phase transition. In Sect. 3 we add the Gauss-Bonnet terms and investigate the effect of these higher derivative terms to the results of Sect. 2. In the last section, we discuss our results and the inequalities of the HRE.

Phase transition and Rényi entropy
In this section we review the phase transition between Einstein and hairy black holes in five dimensions which have been already studied in [40,43]. We also look at the behavior of the Rényi entropies under the scalar field condensation.

The Einstein black hole
We begin this section by introducing the Einstein-Hilbert action together with a cosmological constant term in a five dimensional space-time To study the entanglement entropy of a spherical region with radius R in a quantum field theory side, which is supposed to exist as a gauge theory on R 1 × H d−1 geometry, we need to consider a black hole in the AdS space-time with characteristic scale of L and a hyperbolic spatial boundary. Therefore we start from the following metric where d 2 3 is the metric on a unit hyperboloid in three dimensions By inserting (2.2) into the equations of motion, we find the following solutions for the unknown functions of the metric [43] where r H is the radius of horizon. In a series of straightforward steps, we will compute the temperature, energy and thermal entropy of this black hole. By knowing the metric we may find the temperature of the black hole by using the following definition To compute the black hole energy (mass) in an asymptotically AdS space-time, one could use the Astekhar-Das formalism [44]. A quick but equivalent way is to expand the g 00 in the metric [43] g 00 (r ) = −N 2 r 2 For simplicity we assume that R = L (we can get rid of this coefficient by a rescaling of time). Therefore the black-hole energy would be where V is the regularized volume of the hyperbolic space. Finally to compute the thermal entropy we use the Wald's formula for entropy [31,32] In this formula h is the determinant of the induced metric on the horizon andε ab is a binormal, composed from the Killing vectors ξ b and the normal vectors to the horizon η a , i.e.ε ab = η a ξ b − η b ξ a . At the level of the Einstein-Hilbert action, the Wald's formula coincides with the Bekenstein-Hawking area formula For future purposes it would be useful to introduce a dimensionless parameter x ≡ r H L , so that the temperature, energy and thermal entropy of the black hole can be written as where we have assumed R = L = V L 2 G N = 1 for simplicity. It is easy to show that the first law of thermodynamics, i.e. dẼ =T dS, is satisfied by the above values.
Now we can go one step further and find the Rényi entropies which we have introduced in Eq. (1.10) by integrating the thermal entropy where T (x) x=x n = T 0 n and T 0 is related to the length scale of the hyperbolic space, here T 0 = 1 2π . In order to draw the behavior of the S n versus n, we should substitute the value of x n in S n . It is given by the real positive root of the following quadratic equation (2.12) Substitution gives (2.13) In the limit n → 1,S 1 = 1 4 gives the entanglement entropy and when n → ∞ the Rényi entropy goes to theS ∞ → 5 32 (see Fig. 1). Moreover one can simply check that all the inequalities in (1.12a)-(1.12d) hold for the Rényi entropy in (2.13).

The hairy black hole
In this section we review the scalar field condensation in the background of an AdS black hole with a hyperbolic spatial slicing. We also study the effect of condensation (phase transition) on the behavior of the Rényi entropies. To begin, we add a real scalar field ψ which is coupled to gravity as (2.14) For simplicity we just consider a mass and a kinetic term. To find the desired black hole solution, we choose the following ansatz for the metric and scalar field The back reaction of the scalar filed on the metric will appear through F and χ functions. The equations of motion are where the prime stands for differentiation with respect to the radial coordinate r . To solve this system of coupled nonlinear differential equations numerically, we need to impose boundary conditions on the metric and scalar field at the horizon (r = r H ) and at the space-like boundary (r → ∞) as follows: • F(r ): The definition for the location of the horizon demands that F(r H ) = 0. For asymptotically AdS black holes and at large distances near the infinite boundary, In order to have a hairy black hole, we need a regular scalar field at the horizon of the black hole, so we suppose that ψ(r H ) = O(1). As is well known in the dictionary of AdS/CFT, our gravitational bulk theory must be dual to a conformal field theory that lives at the boundary of the AdS space. This CFT contains a scalar operator with the conformal dimension so that μ 2 L 2 = ( − 4). In order to have an asymptotically AdS black hole, unstable against the scalar field condensation, has to be sufficiently small. We use the benefits of this instability to study the phase transition due to the scalar field condensation. The instability condition requires that the scalar mass takes the values between the two Breitenlohner-Freedman (BF) bounds of AdS 5 and AdS 2 [40,[45][46][47]] To specify the behavior of the scalar field at the AdS boundary, it will be sufficient to study its equation of motion around the AdS-Schwarzschild background, which assumes the form of a Klein-Gordon equation. The solution to this equation behaves asymptotically as where ± = 2 ± 4 + μ 2 L 2 , and A (±) are the expectation values of the conformal operators with conformal dimensions ± , i.e. O ± ≡ A (±) . If we assume the Dirichlet boundary condition then we can keep the fastest falling off mode near the boundary and ignore the other one, hence we choose, • χ(r ): By using the third equation of (2.16) and applying the boundary conditions ψ(r ) ≈ , it is easy to obtain the boundary conditions for χ(r ) In this paper, we are going to find the numerical solutions for equations of motion by using the shooting method. To do this it would be easier to change the radial variable from r to z via z = r H r . This substitution maps the region r H ≤ r < ∞ to the region 1 ≥ z > 0. Here we also introduce two dimensionless parameters, m = −μ 2 L 2 and z 0 = r H L . By these changes the equations of motion assume the form and boundary conditions yield In the shooting method, we need to calculate the series expansions of the scalar and metric functions around the horizon and boundary. By inserting these series into the equations of motion we can find the unknown coefficients of the expansions. After that, we match the series smoothly at some small distance, say , from one end of the interval 0 < z < 1. Note that for a given ψ(z) one can find χ(z) from the last equation of motion in (2.16), simply by an integration.
Since the equations of motion are regular at the horizon, the expansion of functions are given by a Taylor series near the horizon The boundary conditions (2.21) require that the F (0) = 0. By substituting the near horizon expansions of (2.22) into the equations of motion (2.20) and expanding again near the z = 1, we can find the values of F (n) and ψ (n) as functions of free parameters z 0 , m and ψ (0) ≡ ψ 1 . In the precision region of our computations it would be enough to keep expansions up to n = 6. For example the first few terms are 1 The expansion near the boundary at z = 0 is a little bit trickier, because the equations of motion are irregular at this point. Here we have a power series expansion and the power of the leading term depends on the value of m. We denote this, by a mass dependent parameter δ Here we have expanded the functions until the resulting algebraic system of equations gives a non-trivial solution. In Table 1 we have presented the values of m, + , δ and i min which we have used to solve the equations of motion numerically By substituting the expansions in (2.24) for every mass parameter m from Table 1 into the equations of motion, we can find the expansion coefficients, a i and b j in terms of three free parameters C + , C m and z 0 . Here C + is the coefficient Eq. (2.24) (note that C m z 2 always exists in the near boundary expansion of F(z)).
Let us explain the shooting method a little bit. At first step, we read the initial values of the fields and their derivatives near the horizon at z = 1 − from the series expansion in (2.23), remembering that the coefficients are functions of C + , C m and z 0 . Then we solve the first two differential equations in (2.20) numerically to find the values of the fields near the boundary at z = (note that we need the very small parameter as a regulator). Now we can compare these values of the fields (two equations) and their first derivatives (two equations) with those which are computed from the near boundary expansion in (2.24). At the end of the day, we will have four equations for four unknown parameters C + , C m , ψ 1 and z 0 to solve. This set of equations has innumerable numerical solutions which can be found and put into the power series expansions of the fields both at the horizon (2.23) and boundary (2.24) to compute various physical and thermodynamical quantities.
The first quantity that we consider, is again the temperature of the hairy black hole is computed from Eq. (2.23) and can be computed numerically by imposing the boundary conditions of (2.21). The next quantity is the thermal entropy of the hairy black hole where we have calculated the entropy again in units where To compute the energy of the hairy black hole we can use two different ways again. Either we can expand the g 00 and pick up the coefficient of z 2 z 2 0 , therefore the energy yields As we discussed previously, in the AdS/CFT dictionary the coefficient of r − + gives the expectation value of the conformal operator O with a conformal dimension + . This coefficient is called the condensate. In the z coordinates at leading order of ψ(z) ∼ C + z + therefore the value of the condensate is O = C + z + 0 . Figure 2 shows the behavior of the O 1 + in terms of the temperature, for various possible values of the + in the interval 2 ≤ + ≤ 2 + √ 3. The condensation happens for each value of the + at different critical temperaturesT c , which decreases when + increases.
When the AdS 5 BF bound saturates, i.e. + = 2 or m = 4, the critical temperature T c reaches to its maximum value at T 0 = 1 2π ≈ 0.159. This is the Rindler temperature (the temperature of massless black hole).
In the following figures, we have sketched the behavior of various physical quantities as a function of the energy or the temperature.
In each diagram in Fig. 3 the dashed curve represents the Einstein black hole. TheS −T diagram shows how the condensation happens. When temperature decreases there is a critical temperature for each value of the + (the point where solid curves meet the dashed one) in which, the Einstein black hole is replaced by a hairy black hole. As we see the critical temperature decreases when + increases. One can convert the temperature to energy using theT −Ẽ diagram. TheS −Ẽ diagram shows that the entropy of a hairy black hole at a specific value of the energy is larger than that of the Einstein black hole, therefore, the hairy black holes are more favorable (stable) after the condensation. The O −Ẽ diagram displays where the condensa-tion starts and reaches to its maximum value in terms of the energy.
To compute the Rényi entropies from (1.10) according to theS −T diagram in the Fig. 3, we should take into account the phase transition from the Einstein black hole (EBH) to the hairy black hole (HBH) [40] (2.30) The second term above, differs from (2.11) only on its lower limit. Here we can define T (x) x=x c = T c in accordance with the definition of the x n . So theS n yields where x c is the real positive root of the 1 By a numerical computation, we can draw theS n in terms of n, see Fig. 4.
As we observe in Fig. 4, one can define a critical number n c =T T c , where the curve of HBH approaches that of EBH. At this point the second derivative ofS n suddenly changes. This discontinuity of the second derivative ofS n with respect to n confirms that the condensation of the scalar field is a second order phase transition. Moreover, the value of n c increases when + increases as we depicted in Fig. 5. Here all four inequalities (1.12a)-(1.12d) hold again.

A note on the normalizability of the scalar modes
The scalar modes that we have considered so far, are constant on the hyperboloid i.e. ψ = ψ(r ). These modes are called the non-normalizable modes in the context of AdS/CFT. Since these modes preserve the symmetries of the hyperboloid, the construction of fully backreacting hairy black holes in the gravitational background is relatively easy; and for this rea- However, the more physical cases would be those which are not constant on the hyperboloid ψ = ψ(t, r, u, θ, φ), i.e. the normalizable modes. In these cases, the construction of the fully backreacting hairy black holes is considerably more difficult (if not impossible). As an example, the authors of [40] perform a linearized analysis of the Klein-Gordon equation for these modes and manage to show the existence of a phase transition. The main result of this linearized analysis is a phase transition, that occurs if the scalar mass was in the range −4 ≤ μ 2 L 2 ≤ −3 which is more restricting than −4 ≤ μ 2 L 2 ≤ −1.
Conversely, one could say that the linearized analysis restricts the admissible region of conformal dimension from 2 ≤ ≤ 2 + √ 3 to 1 ≤ ≤ 3. Indeed, if one compares the Fig. 5 of our paper to the Fig. 6 in [40], one would see that both curves generally behave similarly, especially in the lower admissible region of the conformal dimension. The reason for the similarity of these two curves is hidden in the fact that if the scalar mass is just above the AdS 5 BF bound and below the unitarity bound (a bound imposed by unitarity conditions) i.e. − 4 L 2 < μ 2 < − 3 L 2 , then the constant modes on the hyperboloid are normalizable, because it is possible to quotient the hyperboloid to form a compact space, see for example [40,43,48,49].
Therefore, despite the difference in the length of the allowed region and the difference in the behavior of the aforementioned curves in the upper part of the allowed region, at least in the lower parts, the constant modes of the scalar field could be regarded as a good approximation to the complicated problem of finding the fully backreacting normalizable modes.
All in all, since our main objective in this paper is to study the effect of the higher derivative terms on the entanglement Rényi entropy in the process of the phase transition, we will consider the constant modes as a good toy model and focus on them in the remainder of this paper. We will leave the study of normalizable modes to future works.

The ERE from the Gauss-Bonnet gravity
In this section, we are going to explore the effect of higher curvature bulk theories of gravity on the condensation and the holographic Rényi entropy. We consider the Gauss-Bonnet (GB) gravity in five dimensions where the equations of motion are of the second order in derivatives and the numerical method of the previous section works here as well without imposing any further boundary condition. The GB gravity has been studied extensively in the literature [50][51][52][53][54], which we will review in the following subsection. Afterwards, in Sect. 3.2, we will present our own original work on the study of the hairy black holes for this gravitational theory.

The modified Einstein Gauss-Bonnet black hole
Similarly to the previous section, we start from a pure gravitational action in five dimensions where λ is a free dimensionless coupling. To find the asymptotically AdS black hole solutions with a hyperbolic horizon, we use again the ansatz in Eq. (2.2) and insert it into the equations of motion. We will find with the following Einstein Gauss-Bonnet (EGB) black hole solution where the constant of integration is fixed in terms of the horizon radius by demanding that F(r ) vanishes at the horizon, i.e. C = r 4 H −r 2 H L 2 +λL 4 . By comparing (3.3) with Eq. (2.4) we can fix the constant value of N (r ) = N . Since asymptot- where L eff is an effective asymptotic AdS scale. Note that asymptotically F(r ) → r 2 In order to avoid an imaginary AdS scale or a naked singularity, the Gauss-Bonnet coupling λ has to be limited to the region λ ≤ 1/4. On the other hand the unitarity of the boundary theory dual to the Gauss-Bonnet gravity in this background, demands that − 7 36 ≤ λ ≤ 9 100 , for example see [55,56], or [57]. Therefore we will restrict our numerical computations to this interval from now on. Now we follow analogous steps, similar to the Sect. 2, to calculate the thermodynamical quantities and then compute the holographic Rényi entropy. The black hole temperature is given by By introducing the X = r H L eff and T 0 = 1 2π R , the temperature can be expressed as [33] To compute the black hole energy, we use again the procedure discussed in the Sect. 2. By expanding the metric near the boundary and assuming that R = L eff we have (3.7) Therefore the black hole energy assumes the form whereẼ is the dimensionless energy of the EGB black hole.
As it was mentioned in Sect. 2, for the higher derivative gravitational theories we must use the Wald's formula (2.8) to compute the thermal entropy. After some algebra we find from which, we can express the thermal entropy as By knowing the thermal entropy we can compute the holographic Rényi entropy. Using the Eq. (2.11) and defining the T (X ) X =X n = T 0 n , we find [33] S n = n T 0 (n − 1) Note that, X n stands for the real positive root of the following third order equation In [36,37] the authors show that in the Gauss-Bonnet gravity there is a new bound on the coupling λ when we demand the positivity of the thermal entropy. In five dimensions, this is − 7 36 ≤ λ ≤ 1 12 . We will discuss this in the next section.

The modified hairy black hole
In Sect. 2, we reviewed the condensation of the scalar field around a five-dimensional asymptotically AdS black hole with the hyperbolic horizon in the Einstein gravity. In this section, we are going to study the effect of higher-derivative terms on the scalar condensation, by considering the Gauss-Bonnet gravity. We first prove the existence of a non-trivial scalar solution below the critical temperature and then compare its thermodynamical properties and holographic Rényi entropy with those of the hairy black hole of the Einstein gravity, in various values of the scalar mass μ and coupling λ. Let us start from the following action where a real scalar field is coupled to the gravity The metric and the scalar field ansatz for the Gauss-Bonnet gravity are the same as those in Eq. (2.15). The equations of motion are where we have defined, V (λ, r ) = 2λL 2 r 2 1 + F(r ) . Since in the Gauss-Bonnet gravity the equations of motion remain the second order differential equations, the number of boundary conditions are the same as in Einstein gravity. Similarly to the previous section, we need to know the behavior of every unknown function near the horizon r = r H and near the boundary when r → ∞.
As always, the definition of the location of horizon F(r = r H ) = 0 gives the first boundary condition on F(r ). On the other hand, as we mentioned previously, the asymptotic behavior of this function is given by We also demand that the scalar field is a regular function at the horizon so the ψ(r H ) = O(1). If we consider the scalar field as a quantum field living in the curved spacetime outside the black hole horizon, even in the presence of the Gauss-Bonnet terms, it should satisfy the Klein-Gordon equation. The solution of this equation decays at infinity as ψ(r ) ∼ are the expectation values of the conformal operators with ± conformal dimensions.
To find a relation for , one needs to substitute F(r ) = r 2 L 2 eff − 1 and ψ(r ) = A r − into the equation of motion for the scalar field (the first equation in (3.14)) in the limit of r → ∞. In this regard, we will find the following values for the conformal dimension In what follows, we again assume the Dirichlet boundary conditions, A (−) = 0, therefore we expect ψ(r → ∞) ∼ A r + asymptotically. It is also easy to derive the boundary conditions for χ(r ) just by looking at the last equation of motion in (3.14) We can use the advantage of working with the dimensionless parameters which is desirable in the numerical methods. To do this we change the radial coordinate from r to z via z = r H r and introduce a dimensionless mass parameter m = −μ 2 L 2 .
We also use the scaling symmetry of the equations of motion to define, z 0 = r H L eff . Then the equations of motion will be the following The boundary condition also can be written as (3.20) These boundary conditions are identical to those of (2.21). As we discussed in Sect. 2, in order to find a numerical solution for equations of motion, we need to compute the expansions of the metric and the scalar functions near the horizon and boundary. Once again, we substitute the near horizon expansion (2.22) into the first two equations of motion in (3.18) and expand them near the z = 1. We find Similarly, we use the near boundary expansion at z = 0 in the Eq. (2.24), but here we note that the leading order is given by the value of + in the Eq. (3.16) or Here we have a two-parameter family of the solutions which are controlled by the values of m and f ∞ or equivalently by μ and λ. In our numerical method, it would be simpler to fix the value of + firstly and then choose some appropriate values for the dimensionless mass and coupling.
Once again, in shooting method we first read the initial values of the fields and their derivatives near the horizon and we find The next physical quantity will be the thermal entropy, which we have calculated for the EGB black hole from the Wald's entropy formula in (3.10). For the modified hairy black hole we can use again this equation just by replacing r H → z 0 L eff . Hence (3.23) The computation of the energy is straightforward and we can use the same formula that we found in Eqs. (2.28) and (2.29) again Since we are interested in the thermal phase transition, we need to find the condensate or the O = C + ( z 0 √ f ∞ ) + . By plotting the condensate versus the temperature (see Fig. 6), we observe that the scalar field condensates at every value of the λ within the allowed interval − 7 36 ≤ λ ≤ 9 100 . As we see in Fig. 6, for the fixed value of + = 2, the critical temperature decreases when the Gauss-Bonnet coupling increases.
To see the relationship between the various thermodynamical quantities and moreover to follow how the condensation happens in this theory, we have depicted their graphs in Fig. 7. The behavior of various quantities is the same as in the case of Einstein gravity. We can see this by comparing the curves of λ = 0 with the curve λ = 0 (black curve) in each graph. For each value of the coupling, there is a critical temperature where a phase transition from the EGB black hole to the modified hairy black hole happens.
But there is an important point which can be seen easily in theS −Ẽ or theS −T diagrams. For the positive values of the coupling λ, there is a temperature or energy below the critical point, where the value of the entropy becomes negative. The value of these temperatures, T v , is depicted in Fig. 8. The same behavior is happening for the EGB thermal entropy in a narrow region of the couplings 1 12 < λ < 9 100 , this has been already reported in [36,37].
Finally the Rényi entropies should be computed using the EGB temperature in the Eq. (3.6). By considering the dimensionless temperature,X c will be the real positive root of the following equatioñ (3.25) Then according to (3.11), the dimensionless entanglement Rényi entropies yield whereT 0 = 1 2π . Once again the first term of this expression should be computed numerically. SinceT c ≥T 0 n we can define a critical value, n c =T 0 T c , where the Rényi entropy corresponding to the modified hairy black hole approaches that of the EGB black hole, see Fig. 9. These diagrams show that in the presence of the Gauss-Bonnet terms the condensation of the scalar field is again a second order phase transition.
To show the behavior of n c under the change of the coupling λ, we have sketched the Fig. 10 (left). As we see, this critical value increases by moving from the negative values of the coupling to the positive values.
Since we have the numerical behavior of the Rényi entropies, it is worth checking the inequalities of (1.12a)-(1.12d). We have sketched the corresponding diagrams in terms of the various values of n, see Fig. 11.
As we observe, the inequalities (1.12a) (up-left) and (1.12d) (down-right) hold here but the second inequality (1.12b) (up-right) is violated for the positive values of the coupling λ (there is an interval for n, where the curves get the negative values). This is completely consistent with the fact that for the positive values of the coupling, there is a temperature where the thermal entropy becomes negative. By a   numerical study, one can find a violation number, n v , where the second inequality governing the EREs is violated, see the right diagram in Fig. 10. The third inequality (1.12c) (downleft) also holds here, and as we see, there is a discontinuity in the slope of the curves at n = 1. This shows a second order phase transition. We will discuss these results in the next section.

Results for + > 2
To complete our results, we have drawn the Fig. 12 for other values of + > 2. In this figure we have presented the critical temperatureT c , the critical Rényi parameter n c and the violation number n v for each value of the Gauss-Bonnet coupling λ. The behavior of these parameters are similar to those of + = 2 but is more pronounced as + increases.

Discussion
In Sect. 2, we have reviewed the phase transition from the Einstein black hole to a hairy black hole due to the condensation of a scalar field. Both solutions exist in a fivedimensional space-time and have a hyperbolic spatial boundary. The numerical results of the condensate in terms of the temperature are sketched in Fig. 2 (the left diagram). We have also found the critical temperatures for this phase transition within the unitarity bound 2 ≤ + ≤ 2 + √ 3, see Fig. 2 (the right diagram). These diagrams show that, by increasing the value of the conformal dimension of the scalar field, the critical temperature of the phase transition decreases. Moreover, the value of the condensate decreases too.
One can explain this condensation in terms of the energy, by finding the relation between the energy of the black holes and their temperatures, via theT −Ẽ diagram in Fig. 3. Thẽ S −Ẽ diagram shows that the entropy of a hairy black hole at a specific value of the energy is larger than the entropy of the Einstein black hole, so the hairy black hole is preferred to the Einstein black hole after the condensation. The O −Ẽ curves display a phase transition when the energy decreases.
As we reviewed in the introduction, by knowing the thermal entropy, one can compute the entanglement Rényi entropies S n , for a spherical region in the dual gauge theory. The results for these entropies are given in theS n − n diagrams of Fig. 4. They show the existence of a critical value, n c , where a discontinuity is happening for the second derivative of the S n . This suggests that the condensation is a second order phase transition.
In Sect. 3 we have extended the study of the phase transition in Sect. 2, by adding the Gauss-Bonnet terms to the action. We have found the modified black hole solutions before and after the condensation, which now depend on the GB coupling λ in the allowed region − 7 36 ≤ λ ≤ 9 100 . The effect of adding this coupling on the condensation is presented in Figs. 6 and 7. Here we have two important parameters, the conformal dimension of the scalar field, + and the coupling λ. We have presented the results for + = 2 completely and for + > 2 briefly in this paper.
We show that by going from the positive values of the coupling to the negative ones the value of the condensate and the critical temperature are increasing, see Fig. 6. By looking at the diagramsS −Ẽ orS −T in Fig. 7, we observe a new behavior for λ > 0. For the positive values of the coupling λ, there is a temperature or energy, below the critical point, where the value of the entropy becomes negative. The value of these temperatures, T v , is depicted in Fig. 8.
As it was mentioned, the violation of the second inequality of the Rényi entropy (1.12b), is related to the negativity of the black hole entropy. In references [36,37], the authors have studied the same violation in the Gauss-Bonnet gravitational theory in the presence of a gauge field and found a new bound on the coupling of the theory i.e. − 7 36 ≤ λ ≤ 1 12 in 5 dimensions, instead of the interval demanded by the unitarity of the theory i.e. − 7 36 ≤ λ ≤ 9 100 . In this paper, we show that for the Gauss-Bonnet gravity and in the presence of a scalar field, if we demand the positivity of the black hole thermal entropy, then there will be a novel bound on the coupling of the theory i.e. − 7 36 ≤ λ ≤ 0. The entanglement Rényi entropy is calculated in the Eq. (3.26) and for the various values of the coupling are depicted in Fig. 9. As we see, a second order phase transition is happening again. The critical points of n c , where the ERE of the hairy black holes meet the ERE of the EGB black holes, have been drawn in the left diagram of the Fig. 10.
One of the important subjects in the study of the entanglement Rényi entropy is the existence of the inequalities. We have checked some of these inequalities, (1.12a)-(1.12d), and the reader can see the results of our numerical computations in Fig. 11.
Let us explain our results according to the Eq. where ξ is a positive number independent of n. As it was discussed, the inequalities (1.12a) and (1.12d) hold in our case. To explain the first inequality we can write (4. 2) The first term in the above equation is negative, because the S HBH is a monotonically increasing function of the temperature, see theS −T diagram in the Fig. 7. Since the ξ is a positive constant, the total value of the right hand side of the above expression will be negative. The fourth inequality can be written in terms of the specific heat, to see this we notice that TheT −Ẽ diagram in Fig. 7 easily shows that the right hand side of the above equation is always negative. Now let us look at the second inequality (1.12b), which is violated for the positive values of the coupling λ, as depicted in the Fig. 11. This is completely consistent with the fact that for the positive values of the coupling there is a temperature where the thermal entropy becomes negative. In fact, we have TheS −T diagram in the Fig. 7 proves that whenever theT 0 n is smaller than a specific temperature then the value of the thermal entropy becomes negative and therefore the second inequality is violated. We have sketched the violation number n v , in the right diagram of the Fig. 12 for the positive values of the coupling. In fact, this violation never occurs for λ ≤ 0, because as we see from Eq. (3.23), for the negative values of the coupling, the thermal entropy is always positive. The third inequality (1.12c) has a more complicated form (4.5) The first and the second terms above, may have positive or negative values, depending on the value of theT 0 n . The last term is always positive, therefore it would be difficult to find the total sign of the above expression, otherwise we use the numerical data. Our numerical analysis shows that, this inequality holds for the most areas of the parameters of the theory as depicted in the Fig. 11 (down-left).
By more numerical analysis for the other values of the + > 2, we conclude that all the above arguments and discussions are + independent, see Fig. 12.