Holographic conductivity of holographic superconductors with higher-order corrections

22 pages, 38 figures. New references added. Text was revised

section_content: The correspondence between the gravity in a d-dimensional anti-de Sitter (AdS) spacetime and the conformal field theory (CFT) residing on the (d -1)-dimensional boundary of this spacetime, known as the AdS/CFT correspondence, provides an established method for calculating correlation functions in a e-mail: asheykhi@shirazu. ac. ir a strongly interacting field theory using a dual classical gravity description [1] [2] [3]. It has been confirmed that this duality can be applied to solving the problem of high temperature superconductors in condensed matter physics [4]. This is due to the fact the high temperature superconductors are basically in a strong coupling regime, and thus one expects that the holographic method could give some insights into the pairing mechanism in these systems. Understanding the mechanism of high temperature superconductors has long been a mysteries problem in modern condensed matter physics. Recently, it was suggested that it is logical to understand the properties of high temperature superconductors on the boundary of spacetime by considering classical general relativity in one higher dimension. This idea is called holographic superconductors (HSCs) [5] [6] [7] and has got a lot of attention in the past decade. According to the HSC proposal, on the gravity side, a Maxwell field and a charged scalar field are introduced to describe the U (1) symmetry and the scalar operator in the dual field theory, respectively. This holographic model undergoes a phase transition from a black hole with no hair (normal phase/conductor phase) to the case with scalar hair at low temperatures (superconducting phase) [8]. 

Nowadays, the investigations on the HSC have attracted considerable attention and have become an active field of research. Let us review some work in this direction. In the background of Schwarzschild AdS black holes in Einstein gravity, the properties of HSCs have been explored in [9] [10] [11] [12] [13] [14] [15] [16]. The studies were also generalized to higher-order gravity theories such as Gauss-Bonnet gravity [17] [18] [19] [20] [21] [22] [23]. It was argued that the critical temperature of the HSC decreases with increasing the backreaction, although the effect of the Gauss-Bonnet coupling is more subtle: the critical temperature first decreases then increases as the coupling tends towards the Chern-Simons value in a backreaction dependent fashion [20]. It was confirmed that the critical exponent of the condensation in Gauss-Bonnet HSCs still obeys the mean field theory and has the value 1/2 [19]. Other studies 123 on the holographic superconductor have been carried out (see for example and the references therein). 

It is also interesting to investigate the electrical conductivity of HSC in the dual CFT as a function of frequency. In the AdS/CFT correspondence, the electrical conductivity can be computed by looking at the linear response of the system to fluctuations of the fields A x and g t x in the bulk. These fluctuations are dual to the electric current J x and energy current T t x operators in the CFT. In the context of linear Maxwell field, the conductivity of HSC was computed in [4, 6, 7]. In the presence of nonlinear electrodynamics, the conductivity of HSC has been investigated in [50, 51]. Also, in the context of Born-Infeld nonlinear electrodynamics, the optical properties of Lifshitz HSC have been explored in [52]. It was demonstrated that this superconductor exhibits a metamaterial property for low frequency of the external electric field for a certain region of the nonlinear parameter. The effects of the Weyl coupling parameter and the Lifshitz dynamic exponent on the conductivity of HSC have been explored in [53]. In Ref. [54], rotating BTZ black holes were considered as the gravity dual to a (1 + 1) dimensional superconductor. In this case, the dependence on the angular momentum on the conductivity has been investigated. Recently, the authors of [55] have analytically computed the holographic conductivity of HSCs in the presence of Born-Infeld nonlinear electrodynamics by considering the backreaction of the matter field on the bulk metric. Further investigations on the holographic conductivity of HSCs have been performed in [56, 57]. 

In this work, we will address the effects of the higherorder corrections on the holographic conductivity of the swave HSC. On the gravity side, we will consider the Gauss-Bonnet curvature correction terms which is the most general action in the 5D spacetime and on the gauge field side we add the quadratic nonlinear gauge term. We shall investigate the effects of these correction terms on the imaginary and real parts of the electrical conductivity of the system. With these correction terms, especially including a Gauss-Bonnet correction to the 5D action, we have the most general action with second-order field equations in 5D [58], which provides the most general model for a s-wave HSC. Furthermore, in an effective action approach to string theory, the Gauss-Bonnet term corresponds to the leading-order quantum corrections to gravity, and its presence guarantees a ghost-free action [59]. The purpose of this work is to analytically and numerically explore the effects of these correction terms on the properties of the s-wave HSC. 

The plan of the work is as follows. In Sect. 2, we set up our model of the HSC in Gauss-Bonnet gravity with nonlinear electrodynamics in the probe limit and derive the equations of motion. In Sect. 3, we analytically and numerically compute the relationship between the critical temperature and the charge density of Gauss-Bonnet HSC. In Sect. 4, we study condensation operator near the critical temperature using analytical and numerical method. In Sect. 5 we investigate the electrical conductivity of the HSC in Gauss-Bonnet gravity with a nonlinear correction term to the Maxwell field. In particular, we shall find the ratio of the gap frequency in conductivity to the critical temperature. Finally, we summarize and discuss our results in Sect. 6.
section_name: Introduction

section_content: We consider the 5D Einstein-Gauss-Bonnet gravity in the background of AdS spaces which is described by the action [60]

where = -6/l 2 is the cosmological constant of 5dimensional AdS spacetime with radius l, α is the Gauss-Bonnet coefficient with dimension (length) 2, R μνρσ, R μν and R are the Riemann curvature tensor, Ricci tensor, and the Ricci scalar, respectively. For convenience, hereafter we set the AdS radius l = 1. We consider the Lagrangian density of the matter field, L M, as

where ψ is a scalar field, q and m are, respectively, the charge and the mass of the scalar field, and the Lagrangian density of the nonlinear electrodynamics is given by [61, 62 ]

where F = -1 4 F μν F μν is the Maxwell Lagrangian and b is a parameter. The term bF 2 is the first-order leading nonlinear correction term to the Maxwell field. There are several motivations for choosing the nonlinear Lagrangian in the form of (3). First, the series expansion of the three wellknown Lagrangians of nonlinear electrodynamics such as Born-Infeld, logarithmic and exponential nonlinear electrodynamics have the form of (3) [63]. Second, calculating a one-loop approximation of QED, it was shown [64] that the effective Lagrangian is given by (3). Besides, if one neglects all other gauge fields, one may arrive at the effective quadratic order of U (1) as F 2 [65, 66]. Furthermore, considering the next-order correction terms in the heterotic string effective action one can obtain the F 2 term as a correction to the bosonic sector of supergravity, which has the same order as the Gauss-Bonnet term [65] [66] [67] [68],

The field equations can be obtained by varying the action (1) with respect to the metric g μν, the scalar field ψ, and the gauge field A μ. We find

where T μν is the matter-stress tensor

The metric of a planar Schwarzschild-AdS black hole in 5D is [69]

with

The Hawking temperature at the horizon can be written in the form

It is worthwhile to note that in the limit r → ∞, we can obtain

so we can introduce the effective AdS radius by

We choose the gauge and the scalar fields in the form [4] A μ = (φ(r ), 0, 0, 0, 0),

Inserting the metric ( 9 ) and the gauge and scalar fields (14) in the field equations ( 6 ) and ( 7 ), we arrive at

The horizon radius is defined as the root of f (r + ) = 0. The regularity condition for the gauge field A t on the horizon r + implies the boundary condition φ(r + ) = 0, which substituting in Eq. ( 15 ) yields

Near the AdS boundary (r → ∞) the asymptotic behaviors of the solutions are given by

where

The value of ± depend on m2 = m 2 L 2 eff. For example, setting m2 = -3, we have + = 3 and -= 1. The coefficients ψ ± correspond to the vacuum expectation values of the condensate operator, namely ψ ± =< O ± >, where O ± is the dual operator to the scalar field with the conformal dimension ±. Following [4], we can impose the boundary condition in which either ψ -or ψ + vanishes, so that the theory is stable in the asymptotic AdS region. In what follows, we set ψ -= 0 and take ψ + = O + non-zero. The interpretation of the parameters μ and ρ, also comes from the gauge/gravity dictionary and they are, respectively, interpreted as the chemical potential and charge density of the conformal field theory on the boundary.
section_name: HSC in Gauss-Bonnet gravity with nonlinear electrodynamics

section_content: In this section, we are going to study the critical temperature of HSC, when the higher-order corrections to the gravity side and the gauge field are taken into account. We shall continue our study both analytically and numerically and compare the results of the two methods.
section_name: Relation between critical temperature and charge density

section_content: Here, we analytically obtain the relation between the critical temperature and the charge density of Gauss-Bonnet HSC. 

To do this, we first transform the coordinate r to z, such that z = r + /r. Using these new coordinates, the equations of motion ( 15 ) and ( 16 ) can be rewritten as

Near the critical temperature (T = T c ) we have ψ = 0, and thus Eq. ( 20 ) reduces to

Solving the above equation for a small value of the nonlinear parameter b, we find

where

Next, we consider the boundary conditions for ψ near the critical point (T → T c ). We assume ψ has the following form [70] :

where F(z) is the trial function near the boundary z = 0, which satisfies the boundary conditions F(0) = 1 and F (0) = 0. Substituting Eqs. ( 23 ) and (25) in Eq. ( 21 ) one arrives at

where the prime now indicates the derivative with respect to z, and p(z), q(z) and w(z) read

It is a matter of calculations to convert Eq. ( 26 ) to the standard form of the Sturm-Liouville equation,

where

In the above equations we have only kept the terms up to order α 3/2. Next, we perform a perturbative expansion bλ 2 and retain only the terms that are linear in b such that

where λ 2 | b=0 is the value of λ 2 for b = 0. Thus we can rewrite Eq. ( 33 ) as

Employing the Sturm-Liouville eigenvalues problem, the eigenvalues of Eq. ( 30 ) can be obtained by varying the following function:

where we also choose F(z) = 1az 2 and m 2 = -3/L 2 eff to assess this expression. Finally, using Eqs. ( 11 ) and ( 24 ), for T ∼ T c, one can obtain

where

and λ min is the minimum eigenvalue which can be obtained by variation of Eq. (36). Our strategy, in the analytical method, for calculating the critical temperature for condensation is to minimize the function (36) with respect to the coefficient a by fixing other parameters of the model such as b and α. Then we obtain λ min and hence the maximum value of T c /ρ 1/3 can be deduced through Eq. (37). As an example, we present the details of our calculation for α = 0. 01 and b = 0. 01. In this case Eq. ( 36 ) reduces to whose minimum is λ min = 25. 6427 at a = 0. 747087. And thus according to Eq. ( 37 ), the critical temperature becomes T c = 0. 185363ρ 1/3. In Tables 1, 2 and 3, we summarize our results for λ min and ζ for different values of the parameters α, a and b. This table shows that, for a small and fixed value of α, with increasing the nonlinear parameter b, the value of ζ = T c /ρ 1/3 decreases as well. As we shall see in the next section this result is in a very good agreement with the numerical results.
section_name: Analytical method

section_content: Now, we numerically investigate the critical behavior of the HSC in Gauss-Bonnet gravity with a quadratic correction term to the gauge field. For the numerical study we employ the shooting method [71]. For simplicity we assume r + = 1, and thus Eqs. ( 20 ) and (21) for φ and ψ reduce to

Near the horizon (z = 1), we can expand φ and ψ as

while near the AdS boundary (z = 0), they behave like

We calculate φ (1), ψ (1) and ψ (1) in the term of ψ(1) and φ (1) by using the equations of motion for φ and ψ, namely Eqs. ( 39 ) and (40), respectively. Since near the critical point ψ is very small, thus we choose ψ(1) = 0. 0001. Our strategy for using the shooting method is as follows. For a specific value of the reduced scalar field mass m2, we can perform a numerical calculation near the horizon boundary with one shooting parameter φ (1) to get proper solutions at the infinite boundary. For specific values of φ (1), we impose the boundary condition ψ -= 0. We also calculate the analytical values and numerical values of ζ for different b. We compare our numerical results with analytical results in Tables 1, 2 and 3. 

In Fig. 1, we plot ψ versus z for three first boundary condition φ ( 1 ), m2 = -3 and different values of the Gauss-Bonnet coefficient α and nonlinear parameter b. In the absence of a quadratic correction term (b = 0), our results exactly coincide with those presented in [19, 24]. The acceptable diagram for us is the red one in each plot, since there is nothing in the bulk to have effect on the speed of the wave, so the diagram of ψ will be stable. From Tables 1, 2 and 3, it is evident that when b becomes larger the condensation gets harder. A similar behavior can be seen for the fixed value of b and different values of α, namely the critical temperature reduces and condensation becomes harder when the Gauss-Bonnet coupling parameter α gets larger.
section_name: Numerical method

section_content: In this section, our aim is to calculate the critical exponent of HSC with first-order correction terms in gravity and gauge field. Again, we continue our studying both analytically and numerically.
section_name: Critical exponent and condensation values

section_content: We would like to obtain the critical exponent and the condensation values of the condensation operator near the critical temperature using the analytical method. Inserting Eq. ( 25 ) into Eq. ( 20

and

Near the critical temperature, O + 2 r 4

+ is a very small and thus we can expand φ(z) as φ(z)

where χ satisfies the following boundary condition:

With the help of Eq. ( 47 ), and comparing the coefficient of

+ on both sides, Eq. ( 45 ) leads to

From Eq. ( 49 ) we figure out, in the limit z → 0, the following equation:

It is a matter of calculations to show that Eq. ( 49 ) can be written

where (z) = 1+z 2 +z 4 +25z 6. Integrating both sides of the above equation in the interval [0, 1] and using the boundary condition (48), we arrive at

where

Combining Eqs. ( 18 ) and ( 47 ), we obtain

where in the last step we have expanded χ(z) around z = 0. Equating the coefficients of z 2 on both sides of Eq. ( 53 ), we find ρ

where

From Eq. ( 55 ) we observe that the critical exponent has the mean field value 1/2, which is independent of the nonlinear parameter b and Gauss-Bonnet parameter α. It is worth noting that O is zero at T = T c and condensation occurs for T < T c. We shall come back to calculation of the condensation value γ in the next subsection.
section_name: Analytical method

section_content: We use the numerical method to explore the behavior of the condensate operator O + in terms of temperature for different values of α and b (see Fig. 2 ). These curves are obtained by the shooting method which we described in the previous section. As one can see from this figure there is a critical temperature T c below which the condensate appears, then it increases quickly as the system is cooled and finally goes to a constant for sufficiently low temperatures. This behavior is qualitatively similar to that obtained in BCS theory and observed in many materials. Now we are going to study the condensation operator O + in a close neighborhood of the superconductor's critical temperature to compute the critical exponents and the condensation value γ of the Gauss-Bonnet HSC with quadratic nonlinear electromagnetic. For this purpose, we first take the logarithm of Eq. (55). We arrive at log

We have plotted the behavior of the above function in Fig. 3. From this figure, we observe that the numerical results are fitted to the above analytic form in the vicinity of the critical temperature. We summarize our results in Fig. 3 and also Tables 4, 5 and 6 for different values of b and α. We see that, for a fixed value of α, the condensation operator γ increases with increasing b, while for a fixed value of b, it decreases with increasing α.
section_name: Numerical method

section_content: In this section, we study the energy gap in the holographic superconductor phase which is constructed on the boundary of the background spacetime. In particular, we investigate the influence of the Gauss-Bonnet and nonlinear parameters on the superconducting gap. In order to do this, we must compute the electrical conductivity of holographic superconductor by turning on a small perturbation δ A x = A x (r ) exp(-iωt) to the gauge field in the bulk where ω is the frequency. At linearized order in the perturbation (δ A x ), the equation of motion for the gauge field A x (r ), which obeys Eq. ( 7 ), is

In the absence of the nonlinear correction (b = 0), this differential equation reduces to the Maxwell case as presented in [4, 5]. To determine the conductivity, we need the asymptotic (r → ∞) form of the second-order differential equation (58),

which admits the following solution near the boundary:

where A (0)

x, A (1) x are two constants and K is also a constant parameter with length dimension which is considered so as to have a dimensionless logarithmic argument. 

According to the AdS/CFT correspondence, the two point correlation function of the current operators in a system is given by its on-shell action where the action is evaluated on the equations of motion. Here, the on-shell action is

which, in the quadratic approximation for the gauge field perturbation, becomes

After performing an integration by parts and using Eq. ( 58 ), we get

Substituting Eqs. ( 12 ), ( 18 ) and (60) in the above expression, one arrives at

thus we can obtain S o. s. as follows:

in which logarithmic divergences appears. In order to cancel out this divergence, we obtain the boundary counterterm as described in the appendix by using Skenderis' method of holographic renormalization [72]. Therefore, the finite onshell action may be written as

where the gauge invariant counterterm S c. t. is given by Eq. ( 79 ). Now, we can obtain the current operator in the boundary field theory [4, 5] as

A (1) x -

According to Ohm's law, the electrical conductivity can be expressed as

where E x = -∂ t δ A x. Hence, using the current (67), the holographic conductivity is given by

Consequently, the holography conductivity is calculated by solving numerically a differential equation ( 58 ) such that the in-falling boundary condition is imposed at the event horizon

in which T is the Hawking temperature and

where a 1, a 2,. . . are calculated by Taylor series expansion of Eq. ( 58 ) around the horizon. The numerical results for holographic conductivity associated with the condensation operator < O + > are plotted in Figs. 4, 5, 6 and 7. The real and imaginary parts of the electrical conductivity versus frequency are illustrated at different temperatures below T c in Figs. 4 and 5, respectively. As one can see from Fig. 4, the superconducting gap has opened below the critical temperature, which became deeper with decreasing the temperature. Besides, for various values of α and b, the real part of the conductivity is proportional to the frequency per temperature, ω/T, as the frequency is large enough. According to Fig. 5 the divergence in the imaginary part, at ω = 0, points out a delta function in the real part, Re[σ ], at ω = 0, which is not displayed in Fig. 4. As one can see from Fig. 5, the holographic conductivity of the HSC has a minimum in the imaginary part. Thus with decreasing the temperature the minimum in the imaginary part shifts toward greater frequency for various values of the Gauss-Bonnet and nonlinear parameters. 

To study the formation of the superconducting gap with changing α and b at low temperature, e. g., T ≈ 0. 1T c, the real and the imaginary parts of holographic conductivity as a function of ω/< O + > 1/3 are plotted in Figs. 6 and 7, respectively. For a fixed value of the Gauss-Bonnet coefficient α, the energy gap (ω/< O + > 1/3 ) enlarges with increasing nonlinear parameter b. It is evident from Figs. 6 and 7 that the energy gap of HSC for various α exhibits different behaviors based on the nonlinear correction b. In the case of the Maxwell field (b = 0), the superconducting energy gap decreases as α increases at low temperature (see Fig. 6d ). When we take into account the nonlinear correction b, the energy gap of the HSC increases with increasing α (see Fig. 6d, g ). From Fig. 7, we see that, for a fixed value of α, the minimum of the imaginary part of the conductivity goes to the larger value of ω/ < O + > 1/3 when b increases. In the absence of a correction (b = 0), it decreases with enhancing Gauss-Bonnet coefficient (Fig. 6d ). Besides, for a fixed
section_name: Holographic conductivity

section_content: In this paper, we continue the studies on the s-wave holographic superconductors (HSC) by taking into account the higher correction terms both on the gravity side and on the gauge field side of the system. We considered the Gauss-Bonnet HSC when the Maxwell Lagrangian has a nonlinear correction term and is written in the form L = F + bF 2, where F is the Maxwell Lagrangian. We have provided several motivations for choosing this kind of Lagrangian for the gauge field. For example, all well-known nonlinear Lagrangians have a series expansion of which the first two terms are exactly in the above form. First, we have analytically and numerically investigated the relation between the critical temperature of the phase transition and charge density which depends on both the Gauss-Bonnet parameter α and the nonlinear parameter b. For this purpose, we employed the analytical Sturm-Liouville eigenvalue problem and the numerical shooting method. We find that, for a fixed value of α, with increasing nonlinear parameter b, the value of T c /ρ 1/3 decreases. This implies that when b becomes larger, the condensation gets harder. A similar behavior can be seen for the fixed value of b and different values of α, namely the critical temperature decreases and the condensation becomes harder when the Gauss-Bonnet coupling parameter α gets larger. We confirmed that these results are in a very good agreement with our numerical results. Then we obtained the critical exponent of the Gauss-Bonnet HSC with a nonlinear gauge field. We observed that the critical exponent has the mean field value 1/2, which is independent of the nonlinear parameter b and the Gauss-Bonnet parameter α. 

Then we explored, numerically, the holographic conductivity of the system. For this purpose, we plotted the real and imaginary parts of the electrical conductivity versus ω/T and ω/ < O + > 1/3 for T < T c. We observed that the superconducting gap has opened below the critical temperature, which became deeper with decreasing the temperature. Interestingly enough, we found that, for different values of α and b, and for large frequency, the real part of the conductivity is proportional to ω/T. We observed that the holographic conductivity of HSC has a minimum in the imaginary part. Besides, with decreasing the temperature the minimum in the imaginary part shifts toward greater frequency for various values of the Gauss-Bonnet parameter α and the nonlinear gauge field parameter b. Furthermore, for a fixed value of b, the minimum of Im[σ ] increases with increasing α, while for a fixed value of α, it increases with increasing b. 

Hence, the asymptotically metric becomes ds 2 = L 2 eff d 2 4 X ( ) + h μν dx μ dx ν, (73) where

According to the electromagnetic contribution, one can evaluate the on-shell action as

where is a small constant parameter. On the new coordinates (72), the gauge field equation of the bulk motion near the boundary is given by

where ∂ 2 0 refers to the wave operator of the boundary metric γ 0 μν and the general solution of this equation is

where ψ = -1/4∂ 2 0 A 0 i. With A i at hand, the on-shell electromagnetic action can be written as

which is logarithmically divergent. Now, according to Ref. [72], in order to determine the counterterm action, we first invert the solution (77) to find A 0 i = A i + O( ). Thus, the counterterm action is obtained:

It is notable that, since we consider A i = A i ( ) exp(-iωL eff τ ), one can calculate
section_name: Conclusions

section_content: Acknowledgements We thank the Research Council of Shiraz University. The work of AS has been supported financially by Research Institute for Astronomy and Astrophysics of Maragha (RIAAM), Iran. Open Access This article is distributed under the terms of the Creative Commons Attribution 4. 0 International License (http://creativecomm ons. org/licenses/by/4. 0/),which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3.

section_name: Appendix: Holographic renormalization

section_content: To construct the boundary counterterm action, we utilize the holographic renormalization method of Skenderis, which was presented in [72] [73] [74]. To apply this method, the spacetime metric takes the form

which relates to the metric (9) via τ = t/L eff, = L 2 eff /r and the asymptotic ( → 0) metric function is X ( ) →0 = 1.
section_name: Appendix: Holographic renormalization