Analytical study of holographic p-wave superfluid models in Gauss–Bonnet gravity

In Gauss–Bonnet gravity, we analytically investigate the p-wave superfluid models in five dimensional AdS soliton and AdS black hole in order to explore the influences of the higher curvature correction on the holographic superfluid phase transition. We observe that the analytical findings are in good agreement with the numerical computations. Our results show that the critical chemical potential of the system increases with the increase of the Gauss–Bonnet parameter in AdS soliton background, while the critical temperature decreases as the Gauss–Bonnet factor grows if the phase transition of the system is of the second order in AdS black hole background, both of which indicate that the higher curvature correction hinders the formation of the condensation of the vector operator. Moreover, the critical exponent of the system takes the mean-field value 1/2, which is independent of the Gauss–Bonnet parameter and the spatial component of the gauge field.


Introduction
The emergence of the anti-de Sitter/conformal field theory (AdS/CFT) correspondence [1][2][3][4], which connects the gravity on asymptotically anti-de Sitter spacetime to the conformal field theory on the (d − 1)-dimensional boundary of this spacetime, provides a dual description for a strongly interacting system from the perspective of classical gravity. Through this correspondence, there is a great progress in comprehending the dynamics of strongly coupled gauge theories. One of the remarkable applications of the AdS/CFT duality is the study of high-temperature superconductor in a e-mail: laichuyu@gzhu.edu.cn (corresponding author) b e-mail: panqiyuan@126.com c e-mail: jljing@hunnu.edu.cn condensed matter physics [5]. The core mechanism that is responsible for the holographic superconductor is the spontaneous breaking of U (1) symmetry at the horizon of bulk black hole, which implies the superconducting phase transition in the dual CFTs. In the past decade, a variety of holographic superconductors have attracted considerable attention and been explored. The holographic superconductor in the background of the bulk AdS black hole can be used to model phase transition from normal state to superconducting state, the critical temperature and critical phenomena were studied in Refs. [6][7][8][9][10][11]. By considering the background of AdS soliton, the holographic models describing the insulator/superconductor transition at zero temperature have been constructed, and the investigations on the properties of these superconductors were carried out in Refs. [12][13][14][15][16][17]. The studies were also generalized to holographic models of superfluid, in which the spatial component of the U (1) gauge field on the boundary is turned on and it corresponds to the current along the same spatial direction, for reviews, see Refs. [18][19][20][21][22][23].
Most of the holographic models are built and investigated based on the Einstein gravity. In order to gain an insight into the influences of the higher curvature correction on the superconductor or superfluid phase transition, we will extend the investigation to the higher dimensional spacetime in Gauss-Bonnet (GB) gravity [24] where α is the GB coupling constant with dimension (length) 2 . The motivation for this is due to the contradiction between the Mermin-Wagner theorem that forbids spontaneously symmetry breaking at finite temperature in spatial dimensions d ≤ 2 on the gauge theory side and the indeed observation that the superconducting phase transition exists in the four-dimensional AdS black hole backgrounds [5,25]. It is possible that higher curvature correction should suppress condensation based on the suppression of the large fluctuations in the large N limit. According to the AdS/CFT dictionary, the higher derivative curvature terms correspond to the corrections of large N expansion of boundary CFTs in the strong coupling limit. In this direction, lots of work studying various holographic superconductors in Einstein-Gauss-Bonnet gravity have been performed and the effects of the GB correction on the phase transition have been discovered [26][27][28][29][30][31][32][33][34][35]. Recently, Nam studied d-dimensional holographic superconductors in the probe limit in the framework of Einstein-Gauss-Bonnet gravity as well as exponential nonlinear electrodynamics, and found that the GB correction makes the formation of condensation harder, and the superconducting energy gap becomes larger when increasing GB parameter [36]. Ref. [37] numerically investigated the holographic p-wave superconductors with GB curvature correction and nonlinear electrodynamics, and further disclosed the effect of the GB parameter on the behavior of conductivity. In the literature [38], we studied the holographic p-wave superfluid model in the background of soliton solution to Einstein-Maxwell theory. The chemical potential, as a parameter of the soliton solution, has a critical value where a second-order transition will be triggered and the AdS soliton will reach the superconductor (or superfluid) phase. Also from the analysis of the critical phenomenon, it was found that the spatial component of the U (1) gauge field will not bring up the first-order transition. It is of interest to further consider the holographic p-wave superfluid model in AdS soliton including higher order curvature correction. The variation method for the Sturm-Liouville (S-L) eigenvalue problem, which was firstly proposed by Siopsis et al. [39,40] to analytically study the critical phenomenon in holographic superconductor model, is generalized to the studies of various holographic models and proved to be an effective approach, see Refs. [10,14,15,29,33,[41][42][43][44]. The analytical results obtained from the S-L method are in great agreement with the numerical findings. We will employ the analytical S-L method to reveal some details of the holographic p-wave superfluid model in GB AdS soliton and disclose several general properties. On the other hand, a holographic superfluid solution was constructed based on the AdS black hole background by Herzog et al. [19], and they found that the second-order superfluid phase transition can change to the first order when the velocity of the superfluid component increases relative to the normal component. By coupling a Maxwell-complex vector field into the AdS black hole, Wu et al. investigated the holographic p-wave superfluid and detected that the lager mass of the vector field leads to the larger translating superfluid velocity from second order to first order, and the Cave of Winds only exists in the case of five-dimensioned spacetime [23]. The holographic p-wave superfluid in AdS black hole background in GB gravity has been studied numerically by Liu et al. recently [45]. It revealed that the GB parameter will make it easier for the appearance of translating point from the second-order transition to the first-order one or for the emergence of the Cave of Winds. In particular, for the sufficiently large mass of the vector field, the phase transition of the system is always of the second order. In order to back up this significant feature by analytical calculation, we use the S-L method to analytically reproduce the properties of the holographic p-wave superfluid model in GB AdS black hole with the large mass of the vector field fixed where the superfluid phase transition is of the second order. In our paper, we will work in the probe limit where the gravitational backreaction of the vector fields on the background geometry is neglected.
This paper is organized as follows. In the next section, we introduce holographic model of superfluid in GB AdS soliton background. We focus on the influences of the higher curvature correction on the critical chemical potential and the critical phenomena of the system via the S-L method. We also study the behavior of the spatial component of the gauge field near the critical point. In Sect. 3, we extend the analytical study to the p-wave superfluid model in GB AdS black hole with the large mass of the vector field. We calculate the critical temperature as well as the condensation of the vector operator to back up the numerical findings and analyze the behavior of the spatial component of the gauge field. In Sect. 4, we summarise the general properties of the holographic p-wave superfluid model in GB gravity with our main results.

Holographic p-wave superfluid model in Gauss-Bonnet AdS soliton background
Since we will discuss the holographic superfluid dual to the five-dimensional GB AdS soliton configuration [46], in the probe limit, we begin with the soliton solution to the action (1) for d = 5 in the form where the metric function with the tip r s is a conical singularity in this solution, α is the Gauss-Bonnet coupling constant and L is the AdS radius. In the asymptotic region (r → ∞), the function f (r ) behaves as so we can define the effective asymptotic AdS scale by where the upper bound of the GB parameter α = L 2 /4 is the so-called Chern-Simons limit [47]. Given the constraints of the causality via the holographic correspondence [48,49], we will take the range −7L 2 /36 ≤ α ≤ 9L 2 /100 for the GB coupling constant in the following calculation.
It should be noted that the solution (2) goes back to the Schwarzschild AdS soliton in the Einstein limit, i.e., α → 0. In order to remove the singularity at the tip, we impose a period β = π L 2 /r s for the coordinate ϕ [46], which is the same as that of Schwarzschild AdS soliton [50]. This period remains unchanged as the Gauss-Bonnet parameter varies, which means that we can fix the period of the spatial coordinate ϕ to discuss the influence of the high order correction on the holographic model. For simplicity, we will scale L = 1 in the following.
In the GB AdS soliton background, considering a Maxwell field and a charged complex vector field coupled, we build the holographic p-wave model of superfluidity via the action where the tensor ρ μν = D μ ρ ν − D ν ρ μ with the covariant derivative D μ = ∇ μ − iq A μ , q and m are the charge and mass of the vector field ρ μ , respectively. In this work the parameter γ , which is the magnetic moment of the vector field, will not be considered. In order to take the possibility of DC supercurrent into account, we adopt the ansatz of the guage field A μ as where both a time component A t and a spatial component A ϕ of the vector potential have been introduced. For the vector field ρ μ , we assume it to be real and take the following ansatz Therefore we can get a set of equations of motion for the holographic p-wave superfluid model where the prime denotes the derivative with respect to r . The effective mass of the vector field obtained from the equation of motion for ρ x reads Evidently, the expression of the effective mass implies that the p-wave superfluid phase transition depends on m 2 , A t and A ϕ . The specific results can be deduced from the following computation.
The boundary conditions at the tip of AdS soliton and the asymptotic AdS boundary have to be imposed to solve the nonlinear equations (9). At the tip r → r s , the solutions behave as At the asymptotic AdS boundary r → ∞, we have the constraints where the characteristic exponent ± = 1 ± 1 + m 2 L 2 e f f , μ and S ϕ are interpreted as the chemical potential and the superfluid velocity, while ρ and J ϕ stand for the charge density and the current in the dual field theory, respectively. The quantities ρ x− and ρ x+ , according to the AdS/CFT dictionary, correspond to the source and the vacuum expectation value of the vector operator O x , respectively. Since we expect that the condensate appears spontaneously, we concentrate on the case of ρ x− = 0 and O x = ρ x+ in our discussion.
In this system, there are several useful scaling symmetries which lead to the transformation of the relevant quantities where λ is a real positive number. Using the scaling symmetries (13) we can assume q = 1 and r s = 1 without losing the generality in the following numerical calculation.

Critical chemical potentials
It was found that [27,51,52], for the AdS soliton, the chemical potential μ of the system has a critical value μ c , above which the solution becomes unstable to develop a hair. It implies the appearance of the condensation of the vector field ρ x in gravity side and the dual field theory reaches a superconductor (or superfluid) phase. For lower chemical potential μ < μ c , the vector condensate vanishes, and the dual theory is in an insulator phase because in this model the normal phase is described by an AdS soliton where the system exhibits a mass gap. The critical chemical potential is a turning point of the insulator/superconductor phase transitions.
To proceed further, we would like to rewrite the equations of motion (9) in the new coordinates z = r s /r , so they yield where the prime denotes the derivative with respect to z, and f (z) is denoted by At the critical chemical potential μ c , the vector field ρ x is zero. So the equation of motion for the gauge field component A t turns out to be The above equation can not be solved directly, but considering the GB factor α is small, we can still get the effective general solution In order to keep A t finite, we have to fix c 1 = 0 in view of the Neumann-like boundary condition for A t at the tip of the soliton. Thus, it leads to that the physical solution of A t is a constant μ at the phase transition point.
Similarly, the equation of the spatial component A ϕ in Eq. (15) simplifies to and the solution is which fulfills the condition A ϕ (1) = 0 given in Eq. (11). Now we consider the case that the chemical potential is away from (but very close to) the critical value μ c , the field equation for the vector field ρ x approaches to From the boundary condition given in (12), we suppose that ρ x takes the form where we have defined a trial function F(z) with the boundary conditions F(0) = 1 and F (0) = 0, and (21) can be rewritten as Here we introduce a new parameter proposed in Ref. [53] which can help to prevent the problem of the divergent behavior for the larger values of the parameter k = S ϕ μ c that we had introduced in our previous work [38].
Following the Sturm-Liouville eigenvalue problem, the equation of motion for F(z) can be converted into a standard Sturm-Liouville form with  Therefore, we can obtain the minimum eigenvalue of = qμ/r s which minimizes the following functional To further estimate it, we assume the trial function to be F(z) = 1 − az 2 with a constant a.
In Tables 1, 2  19, which gives the critical chemical potential c = 1.62945. We can observe that the analytical results derived from the S-L method are greatly consistent with the numerical computations. This shows that the S-L method is an effective analyti-cal approach to investigate the holographic p-wave superfluid model in Gauss-Bonnet gravity.
From Tables 1, 2 and 3, we see that the critical chemical potential in our model, as a transition point, depends on the GB coupling constant α and the superfluid parameter K for different mass of the vector field. These attractive features can also be shown in Figs. 1 and 2, in which we plot the condensation of the vector operator O x = ρ x and the charge density ρ as a function of the chemical potential μ in the holographic p-wave superfluid model, respectively. The increase of the GB coupling constant α generally results in the rise of the critical chemical potential μ c for all values of K with the fixed mass, which means the higher curvature correction makes the phase transition more difficult to take place. Also the superfluid velocity is responsible for the condensation formation of the vector field in GB gravity. The spatial component of the gauge field modeling the superfluid hinders the phase transition, which is shown up phenomenally by the fact that the larger superfluid parameter K leads

Critical phenomena
In this subsection we move on to investigate the condensation of the vector operator and the relation between the charge density and the chemical potential, which can help us further understand the properties of the p-wave superfluid phase transition in GB AdS soliton background.
Near the critical point, i.e., μ → μ c , the vector field ρ x approximates to (22), so the equation of motion for A t presents It should be noted that the condensation of the vector operator O x is very small when μ → μ c , so we are able to expand A t (z) in O x as follow We have to impose χ(1) = 0 and χ (1) = constant to recover the boundary condition at the tip. For simplicity, the following function is introduced then using Eqs. (28) and (29), we can easily arrive at the equation of motion for ξ(z) with Considering the asymptotic behavior of A t in Eq. (12) and expanding ξ(z) near the boundary z = 0, we get with where the integration constants c 2 and c 3 can be determined by the boundary condition of χ(z). Concretely, we take the case of K = 0.25 and m 2 L 2 e f f = 5/4 as an example, Eq. (34) gives O x ≈ 1.98909(μ − μ c ) 1/2 for α = −0.19, which is consistent with the numerical result shown in Fig. 1. We can find that near the critical point, for various values of the GB coupling constant α, the order parameter O x always yields which reveals that in the GB AdS soliton background, the phase transition of the holographic p-wave superfluid model represents the second-order phase transition, with the critical exponent of system taking the mean-field value 1/2. Moreover, both the GB coupling constant α and the spatial component of the gauge field do not have any effect on the second-order phase transition. Figure 1 illustrates that there is a second-order phase transition when the chemical potential μ arrives at the critical value μ c , and the AdS soliton reaches the superconductor (or superfluid) phase for larger μ, which confirms the above analytical conclusion. Also we note that for the vanishing superfluid velocity, the situation is the same as the one discussed in Ref. [29]. On the other hand, from the coefficients of the z 2 term in Eq. (33), we know that the dependence of the charge density ρ on the chemical potential μ reads From the integration of both sides of Eq. (31) and the fact that ξ (0) = 0 revealed by the comparison of the coefficients of the z 1 term in Eq. (33), we have Then, inserting the above equation into Eq. (37), we attain ρ where (α, K , m) is given by The charge density ρ for the holographic p-wave superfluid depends on the GB factor α, the parameter K and the mass of the vector field m 2 via function , and it is proportional to (μ − μ c ).

A ϕ in holographic p-wave superfluid
Let us remind that the vector field approximates to ρ x ∼ O x r + s z + F(z) near the critical point, so the equation of motion for the spatial component A ϕ of the vector potential can be rewritten as where F(z) obeys F(0) = 1 and F (0) = 1. Similarly to the procedure in the preceding subsection, as μ → μ c , A ϕ can be expanded in small O x as From Eqs. (41) and (42), we easily find that the introduced function w(z) satisfies the following equation: We can expand w(z) close to the boundary z = 0 as Substituting the above expansion into Eq. (42), and comparing with the boundary condition (12) for A ϕ , from the z 2 term, we can deduce that the superfluid current reads This expression shows that, in soliton background, GB coupling constant α does not directly affect the current but through μ c (α) due to the relation O x ∼ (μ − μ c ) 1/2 .
Eventually the behavior of A ϕ near the critical point can be described as (46) which matches the behavior of A ϕ in Eq. (20) at the critical chemical potential μ c .

Holographic p-wave superfluid model in Gauss-Bonnet AdS black hole
The holographic p-wave model of superfluidity in AdS soliton background is used to describe a system at zero temperature. In this section we employ the S-L method to analytically investigate the holographic p-wave superfluid model with a certain temperature. For this purpose, our attention is concentrated on the Maxwell complex vector field model with the the same action of matter as Eq. (6), and taking into account the background of five-dimensional AdS black hole in the GB gravity where , in units in which the AdS radius is unity, i.e., L = 1. The Hawking temperature of the black hole is related to the horizon radius r + as T B H = r + π . We suppose that the gauge field has the form A = A t (r )dt + A y (r )dy and the vector field presents ρ μ dx μ = ρ x (r )dx. Then we have the field equations where the prime denotes the derivative with respect to r . Imposing the regularity conditions at the horizon r = r + , we have It should be noted that A t is vanish at r = r + , which is in strong contrast to the case in AdS soliton where time compo-nent A t (r s ) of the the vector potential is a nonzero constant at The parameter λ min = qμ r +c gives the value of horizon radius r + . Since the analysis is valid close to the transition point T c ∼ T B H , using the definition of the Hawking temperature T B H of the black hole, we can compute the critical temperature of the superconductor by the following relation  Table 4, we exhibit the typical values of the critical temperature T c which depends on the GB coupling constant α and the superfluid parameter K for the case of the fixed mass of the vector field m 2 L 2 e f f = 3. We can clearly see that the critical temperature T c diminishes as the GB parameter α increases for various of K . That is to say, the higher curvature correction will make the condensation of the vector operator harder to be formed. In addition, for the same strength of the curvature correction, the critical temperature T c decreases when the superfluid velocity K becomes bigger, which backs up the numerical findings in Ref. [45]. Furthermore, the comparison with the numerical findings indicates that the consistency of our analytical results derived from the S-L method with the ones from numerical calculation is impressive.

Condensation values
In this subsection, we aim to investigate the critical phenomena of the p-wave superfluid model in AdS black hole and figure out the correlation between the condensation operator and the GB coupling constant. For T below but not far from the critical one T c , we can expand the time component A t of the gauge field in small vector operator O x with the boundary conditions κ(1) = 0 resulting from A t (1) = 0 at the horizon. We redefine κ(z) by a new function ς(z) as so near the critical temperature, we easily get the equation of motion for ς(z) Near the boundary z = 0, from the asymptotic behavior of A t and the expansion ς(z) = ς(0) + ς (0)z + 1 2 ς (0)z 2 + · · · , Eq. (62) can be rewritten as Comparing the coefficients multiplying z 0 , z 1 and z 2 in both sides of the above formula, respectively, we have Obviously, there is a relation ς(0) = −ς (0)/2 to make the above equations ture. Therefore, for the further calculation we need the value of ς (0). On the other hand, we observe that the third part of Eq. (66) is consistent with the following relation by making integration of both sides of Eq. (64) In virtue of Eqs. (61), (66) and (67), the correlation between the condensation value of the vector operator O x and the temperature is given by increasing superfluid velocity, i.e., K = q S y /r + = 0 (black), 0.25 (blue), 0.50 (red), 0.75 (orange) and 1.00 (green), respectively. We scale q = 1 and r + = 1 in the numerical calculation where we have set To concretely detect the condensation operator of p-wave superfluid model with respect to T in GB AdS black hole with large mass of the vector field by numerical shooting method, we focus on the case of the mass fixed by m 2 L 2 e f f = 3. In Fig. 3, we plot the condensate of the vector operator as a function of temperature for several values of the superfluid velocity K = q S y /r + with the fixed GB factor α = −0.19, 0.0001, 0.09, respectively. We remark that in each panel, the curves for different K are similar. In other words, although the coefficent ϒ is different, near the critical point T c , the condensation operator O x always takes the form which implies that the phase transition of the holographic p-wave superfluid model belongs to the second order and the critical exponent of the system takes the mean-field value 1/2 for the case of the fixed mass m 2 L 2 e f f = 3 in AdS black hole background. This result is independent of the GB coupling constant α. Moreover, we see that for the fixed value of the superfluid velocity K , the critical temperature decreases when the GB factor α grows. Again, it states that the higher curvature correction makes the condensation of the vector operator harder to occur.

A y in holographic p-wave superfluid in AdS black hole
Having in mind that the vector field ρ x (z) when T tends to T c . The spatial component A y (z) satisfies the equation of motion in the form Near the critical point, A y takes the approximation Inserting Eq. (72) into Eq. (71), the equation of motion for τ (z) is given by where Q(z) has been introduced in Eq. (32). As we processed in the preceding subsection, the expansion of τ (z) near z = 0 presents Substituting the expansion into Eq. (72) and comparing the coefficient of the z 1 term with the one in the asymptotic behavior of A y on the boundary, we have τ (0) → 0. By virtue of this condition and integrating both sides of Eq. (73), the following relation is satisfied Consequently, we get where the term multiplying z 2 stands for the superfluid current. In view of the dependence of O x on the temperature, from Eq. (70), we can easily find that the current J y depends on the GB coupling constant α indirectly via T c (α), and there is a linear relation between the current and temperature, i.e., J y ∝ (1 − T T c ).

Conclusions
In this paper, we have analytically investigated the general properties of holographic p-wave model of superfluid in Gauss-Bonnet gravity by employing the Sturm-Liouville eigenvalue problem. We performed our analysis in the background of AdS soliton as well as AdS black hole in the probe limit and deduced the dependence of the critical chemical potential and the critical temperature on the GB coupling constant α, which are verified by numerical calculation and in good agreement with the numerical findings. Our results have shown that: for the case in AdS soliton with different mass of the vector field, the critical chemical potential increases with the increase of the GB coupling constant α, while for the case in AdS black hole with the larger mass of vector field, the critical temperature decreases for the larger GB factor. Both of the observations reveal the fact that the higher curvature correction hinders the formation of the condensation of the vector operator in holographic p-wave superfluid models. On the other hand, based on the S-L method, we attain the condensation of the vector operator and the charge density with respect to the chemical potential in soliton case, the relation between the condensation value and temperature in black hole case and the behavior of the spatial component of the gauge field near the phase transition point for the p-wave superfluid model. In the AdS soliton background, the superfluid phase transition is always of the second order and the critical exponent of the condensation operator takes the mean-field value 1/2, which cannot be affected by the GB coupling constant α. Especially, we analytically demonstrate that, for the p-wave superfluid model in AdS black hole, when the mass of the vector field is sufficiently large, the phase transition always belongs to the second order for different values of superfluid velocity, that is to say, the spatial component of the gauge field modeling the superfluid cannot bring up the first-order phase transition. And this conclusion is independent of the GB parameter. We also carry out a numerical study of the condensation of the vector operator by using the shooting method. All the analytical results obtained from the S-L method are perfectly consistent with the numerical computation in holographic pwave model of superfluid in GB gravity. Since in this paper we have only considered the probe limit, where the backreaction of matter fields on the metric background is neglected, it is worthy to extend the investigation to the Gauss-Bonnet holographic superfluid model away from the probe limit and take the backreaction into account. It would also be of interest to study other characteristics of these systems like the behavior of conductivity or optical features.