Holographic s-wave superconductors with Horndeski correction

Via both numerical and analytical methods, we build the holographic s-wave insulator/superconductor model in the five-dimensional AdS soliton with the Horndeski correction in the probe limit and study the effects of Horndeski parameter k on the superconductor model. For the fixed mass squared of the scalar field (m2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m^2$$\end{document}), the critical chemical potential μc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _c$$\end{document} increases with the larger Horndeski parameter k, which means that the increasing Horndeski correction hinders the superconductor phase transition. Meanwhile, above the critical chemical potential, the obvious pole arises in the low frequency of the imaginal part of conductivity, which signs the appearance of superconducting state. What is more, the energy of quasiparticle excitation decreases with the larger Horndeski correction. Furthermore, the critical exponent of the condensate (charge density) is 12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{1}{2}$$\end{document} (1), which is independent of the Horndeski correction. In addition, the analytical results agree well with the numerical results. Subsequently, the conductor/superconductor model with Horndeski correction is analytically realized in the four- and five-dimensional AdS black holes. It is observed that the increasing Horndeski correction decreases the critical temperature and thus hinders the superconductor phase transition, which agrees with the numerical result in the previous works.


Introduction
Due to the clear strong/weak dual dictionary, the AdS/CFT correspondence provides a powerful tool to study (d + 1) dimensional strongly coupled systems by means of its ddimensional weak AdS gravity [1,2]. In the recent years, the AdS/CFT correspondence and its generalized version (gauge/gravity duality) have been widely applied in many strong correlated systems [3][4][5][6], especially the holographic a e-mail: lujunwang.2008@163.com (corresponding author) b e-mail: ybwu61@163.com construction of the high critical temperature superconductor [7,8].
In contrast with the conventional superconductor, the high critical temperature superconductor is believed to involve strong interaction and thus its pair mechanism is one of the most challenges in condensed matter physics. The AdS/CFT correspondence obviously opens a new window to understand the properties of superconductors. In 2008, by using the Einstein-Abelian-Higgs system, the authors in Ref. [7] built numerically a holographic description of the s-wave conductor/superconductor model in the probe limit, and produced the main characters of the superconductor, such as the appearance of the scalar condensate accompanied by spontaneous breaking of the U (1) symmetry below the critical temperature as well as the infinite DC conductivity on top of the hairy gravitational background. Soon after, the Meissner effect and the vortex lattice solution about the superconductor were managed to realize via the holographic duality [9], following which the numerical results about holographic superconductor model were upheld by the analytical Sturm-Liouville (S-L) eigenvalue method [10].
Of course, there are also many works which took into account simultaneously the effects of both directions of development on the holographic superconductors, for example, the superconductor model with backreaction in Gauss-Bonnet gravity [38], the superconductor with momentum relaxation and Weyl correction [51] and the MCV p-wave model with Weyl correction [47] or R F 2 correction [45] and the MCV p-wave superfluid in AdS soliton [44] as well as superfluid with vortex lattice [52]. In addition, there is a more natural gravity theory worthy of attention for the construction of holographic superconductor, i.e., the Horndeski gravity [53]. The Horndeski gravity has a series of advantages. First of all, although its lagrangian contain high order derivative about the metric and scalar field, the equations of motion and the energy-momentum tensor involve less than second order derivative of the fields which is similar to the Lovelock gravity. Secondly, as one constructs the s-wave superconductor model by means of the scalar field of the present theory, the coupled term between the Einsten tensor and the scalar provide a new kind of interaction compared with the above corrections. Furthermore, the s-wave conductor/superconductor phase transition in the Horndeski gravity showed that the increasing Horndeski coupling makes the condensate more difficult [54][55][56]. Meanwhile, the Horndeski coupling was argued to play the role of a high concentration of impurities in a material [54,57]. In terms of other investigations of the Horndeski coupling in holography and cosmology, see, for example, Refs. [58][59][60][61][62][63]. Based on the above works, it is interesting to ask how the Horndeski coupling parameter affects the insulator/superconductor phase transition. Therefore, we will generalized the insulator/superconductor model from the minimal coupling between the gravity and scalar field to the nonminimal case by introducing the Horndeski correction and study the effects of the Horndeski coupling parameter on the critical value as well as the conductivity. It follows that the increasing Horndeski coupling inhibits the scalar condensate to appear and suppresses the growth of the scalar condensate.
This paper is planed as follows. In Sect. 2, we will numerically build the s-wave insulator/superconductor phase transition and thus calculate the conductivity, following which we restudy the critical value and the critical behavior by the analytical S-L method. To testify numerical results of the s-wave conductor/superconductor model with the Horndeski correction, we study the corresponding superconductor model by the S-L method. The final section is devoted to the conclusions and discussions.

Insulator/superconductor phase transition
In this section, we firstly construct the holographic swave insulator/superconductor phase transition in the fivedimensional AdS soliton with the Maxwell complex scalar field coupled to the Einstein tensor (i.e., the Horndeski correction) via the numerical method. To verify that above the critical chemical potential the state with scalar hair is indeed stable, we define the grand potential of the system density and then compare the grand potential of the hairy state with the one of no hairy state, following which the frequency dependent conductivity is studied. In order to check the reasonability of the model, we reconstruct the holographic superconductor model by the analytical S-L method to back up the above numerical results.
As a preparation, we start with the review of the fivedimensional AdS soliton, the action of which consists of the Ricci scalar and a cosmological constant as [64,65] where = − 6 l 2 with l the radius of the AdS spacetime. For simplicity, we will set l = 1 throughout the paper. Varying the action (1) yields the Einstein equation of motion, by solving which we can obtain the five-dimensional AdS soliton solution as [29,[64][65][66][67][68][69] Some remarks about the soliton solution (2) are as follows: firstly, the parameter r 0 denotes the tip of the soliton geometry. To avoid the conical singularity at r = r 0 , we impose a periodicity for the Scherk-Schwarz circle on the spatial direction η, i.e., η ∼ η + π r 0 , which thus suggests that the present five-dimensional solition corresponds to a three-dimensional CFT. Secondly, different from the black hole, the soliton solution has no event horizon, and thus no the Hawking temperature defined from the horizon, of course the associated entropy vanishes. Thirdly, for the boundary topology (R × S 1 × R 2 ) to the soliton solution, there exist other two kinds of bulk solutions as for the action (1), i.e., the Ricci flat AdS black hole and the AdS space [64,65]. In fact, the soliton solution can be obtained from the double Wick rotation to the Ricci flat AdS black hole solution. Meanwhile, Horowitz and Myers have conjectured that the soliton solution is not only a negative energy solution but also the minimum energy solution compared with the AdS vacuum and the Ricci flat AdS black hole with the same boundary topology [67]. Choosing the AdS soliton as a reference background, the authors in Ref. [66] found that there exists a Hawking-Page phase transition between the Ricci flat AdS black hole and the AdS soliton if at least one of black hole horizon coordinate is compact, which corresponds to the confinement/deconfinement phase transition in the dual CFT [70]. Fourthly, due to the soliton's lower energy than the one of the AdS vacuum, the soliton gravitational background corresponds to the field theory with a "mass gap", which is similar to the insulator in the condensed physics and is thus believed to provide a gravitational background for modeling the holographic insulator/superconductor phase transition [29,35,68,69].
Subsequently, we take the Lagrangian density consisting of a Maxwell field and a complex scalar field coupled to the Einstein tensor as [54,55] where and m (q) is the mass (charge) of the scalar field ψ. Obviously, the parameter k characterizes the strength of the Horndeski correction. In particular, as k = 0, the above Lagrangian density restores to the case in Refs. [7,8,16,68,69]. Besides, to simplify the calculation, we regarded the matter field as the probe to the soliton, where the equations of motion for both the scalar and the gauge field decouple from the Einstein field equation and the main physics of the system is believed to be grasped at the same time. From Eq. (3), we can read off the equations of motion of scalar field and gauge field as Throughout the paper, we will set L = 1 and q = 1 without loss of generality.
Following the works in Refs. [7,8,29,54,55,68,69], we take the complex scalar field to be real and only turn on the time component of the Maxwell field, which are Combining the ansatz (6) with the soliton background (2), the concrete equations of motion in term of ψ(r ) and φ read where , the prime stands for the derivative with respect to r . Obviously, as k = 0, Eqs. (7) and (8) restore to the standard case in Refs. [29,35,68,69].

Numerical part
To solve the above coupled differential equations by the shooting method, we have to specify the boundary conditions for ψ(r ) and φ(r ). It should be noted that the constant μ is a trial solution of Eq. (8) which is obviously different from the AdS black holes requiring the gauge field vanishing [7,8,54]. Here, we only impose the Neumann-like boundary condition [29] to remove the logarithm term in order for both ψ(r 0 ) and φ(r 0 ) to be finite at the tip r = r 0 . At the boundary r → ∞, ψ(r ) and φ(r ) behave as where ± = 2 ± 4 + m 2 1+6k , ψ 1 , ψ 2 , μ and ρ are all constants. According to the gauge/gravity duality, ψ 1 (ψ 2 ) can be regarded as the source (the vacuum expectation value) of the dual operator O, and μ(ρ) is chemical potential (charge density) of the dual field theory. Since the U(1) symmetry is expected to be broken spontaneously, we impose the sourcefree condition ψ 1 = 0. Hereafter, we denote = + for simplicity. The mass squared m 2 of the scalar field has a lower bound as m 2 = −4 − 24k. In the present paper, we choose m 2 = − 63 16 and mainly consider the effects of the Horndeski correction with the range 0 ≤ k ≤ 1.  Table 1 The critical chemical potential and the fitted coefficient of the condensate (C 4 ) and the charge density (C 5 ) for various We still take r 0 = 1 in the numerical calculation. Thus the period of the spatial coordinate χ is = π . In order for the insulator/superconductor model with different Horndeski correction to be comparable, we fix the period of the spatial coordinate χ , which is another difference from the case of the black holes requiring either the charge density or the chemical potential to be fixed. 1 Fortunately, different from the Lifshitz case in Ref. [35], the period is independent of the Horndeski correction, which implies that the present work does not need rescale the numerical result if we fix = π . After series of calculation, we plot the scalar condensate O and the charge density ρ versus chemical potential with = 9 4 (black solid), 5 2 (red dashed), 3 (blue dotted) and 7 2 (purple dotdashed) in Fig. 1, from which we have the following remarks. As for the curves of scalar condensate, first of all, for all cases there always exists a critical chemical potential, above which the scalar hair starts to condense. Meanwhile, the critical chemical potential increases with the larger Horndeski parameter k, which means that the larger Horndeski correction hinders the insulator/superconductor phase transition. It should be mentioned that we have also listed critical chemical potential in Table 1 and plotted the value of μ c versus the Horndeski parameter k in the form the black solid points in the left panel of Fig. 2. It follows that when = 5 2 , 3, the results return to the ones in Refs. [29,35,68,69]. What is more, near the critical point, we have 1 In addition, as mentioned above, the soliton solution has no temperature, so the insulator/superconductor phase transition in the soliton background is usually restricted to occur at the zero temperature and the related condensate is always displayed as a function of the chemical potential but not the Hawking temperature [29,35,68,69].
√ μ − μ c by fitting the numerical curves. The critical exponent of the condensate 1 2 indicates that the system undergoes a second-order phase transition at the critical point. Meanwhile, we read off the coefficient C 4 in Table 1 and find it decreases with the increasing Horndeski parameter, which suggests that the larger correction suppresses the growth of the condensate and is consistent with the phenomenon that the critical chemical potential increases with the larger Horndeski correction. In term of the curves for the charge density, we observe that above the critical point, the charge density appears and increases with the chemical potential. By fitting the numerical results, we find the charge density has the linear dependent on the chemical potential as ρ ∼ C 5 (μ−μ c ), which agrees with the meaning field theory. Besides, the coefficient C 5 listed in Table 1 decreases with the improving Horndeski parameter, which agrees with the fact from the condensate curve that the increasing Horndeski correction makes the phase transition more difficult.
To check that above the critical chemical potential, the superconducting state with scalar hairy is indeed thermodynamically favored in contrast with the normal state, we define the 'temperature' 2 of the soliton as dt = 1 T , upon which the grand potential of the system is defined by the Euclidean onshell action S E timing the 'temperature' of the soliton, i.e., = T S E . Integrating the Minkowski action (3) by parts yields the on-shell part of action as where we have taken into account d 3 x = V 3 , and also Eqs. (4) and (5). Keep in mind that S E = −S os , the density of the grand potential reads Due to the existence of the coefficient C 3 , it is clear that the Horndeski parameter k will affect the grand potential. Especially, in the case of k = 0, Eq. (11) returns to the standard AdS case [29]. Next, we typically display the grand potential with respect to the chemical potential for the case of = 9 4 , 5 2 , 3, 7 2 in the right panel of Fig. 2, from which we can observe that near the critical point, the superconducting curve stretches out from the insulator horizontal line smoothly with the increasing chemical potential, which means that the system indeed suffers from a second-order phase transition at the critical chemical potential, and thus agrees with the behavior of the condensate in Fig. 1. What is more, the value of the grand potential of the superconducting state is always lower than the one of the normal state, which means that above the critical point, the superconducting state is indeed thermodynamically stable. In addition, the behaviors of the other values of the Horndeski parameter in 0 ≤ k ≤ 1 are similar to the case in Fig. 2. The above results uphold clearly the numerical results about the condensate and charge density.
As we all know, the infinite DC conductivity is a typical character of the superconductor. So it is meaningful to study the conductivity to check whether the hairy state is indeed superconducting state. According to the gauge/gravity duality, we should calculate an electromagnetic perturbation in the hairy soliton to obtain the retarded current-current two-point function. Concretely, we turn on the perturbation δ A = A x (r )e −iωt , and thus obtain the linear equation as In order for A x to be finite at the tip, we take the ansatz of A x near the tip where A x1 , A x2 and A x3 are all constants and the leading term is taken to be unity due to the linearity of the equation for A x . At the boundary r → ∞, the asymptotical expansion of A x is of the form where A (i) , and ξ are all constants. From the gauge field perturbation we can then obtain the Green function as It should be noted that we have removed the logarithmic divergence by the holographic renormalization [8]. Hence, the AC conductivity reads In Fig. 3, we typically show the imaginal part of conductivity as a function of the frequency for different values of chemical potential and , from which some remarks are in order. In term of the left panel, the imaginal part of the conductivity vanishes at the critical chemical potential, which corresponds to the finite conductivity. However, when the chemical potential increases away from the critical point, such as μ μ c ≈ 3 2 , 2, a clear pole appears in the low frequency region, which suggests the infinite DC conductivity as expected from the superconducting state. Meanwhile, the value of the location for the second pole of conductivity increases with the larger chemical potential, which implies that the larger chemical potential increases the energy of quasiparticle excitation. As for the right panel, due to the fact that all curves are from superconducting state, it is reasonable that there always exists a pole in the low frequency. Furthermore, for the fixed ratio of the chemical to the critical value, the location of the second pole moves toward left when one increases the value of . In addition, for all cases, we find the behavior of the conductivity is independent of the Horndeski correction.

Analytical part
To back up the above numerical results, especially, the effects of the Horndeski correction on the critical chemical potential and the condensate, we construct the s-wave insulator/superconductor model by the analytical S-L method. Concretely, we need resolve analytically the coupled differential equations (7) and (8) with the same boundary conditions as the ones in Sect. 2.1. For simplicity, by introducing a new variable z = r 0 r , Eqs. (7) and (8) can be rewritten as where the prime represents the derivative with respect to the new variable z.
In the normal phase, ψ(r ) = 0, the general solution φ(r ) from Eq. (18) is of the form where C 6 and C 7 are two constants. As we have analyzed in Sect. 2.1, we take C 7 = 0 in order for the gauge field to be finite at the tip via the Neumann-like boundary condition [29,68,69]. Hence, the constant C 6 = μ is regarded as the chemical potential in the dual field theory. Obviously, the charge density vanishes in the normal phase, which agrees well with the numerical results in Fig. 1.
When the chemical potential goes slightly beyond the critical point, the scalar condensate begins to condense and can be expressed as where F(z) is a function to be determined with the boundary condition F(0) = 1. Plugging Eq. (20) in to Eq. (17) yields the equation of F(z) as Multiplying the factor T = z 2 −3 (1 − z 4 ) to the above equation reads the S-L eigenvalue equation as where P and Q are given by The minimal eigenvalue μ 2 c is obtained by minimizing the following expression with respect to the coefficient α with the boundary condition F (0) = 0 [68,69]. According to the boundary conditions for F(z), i.e., F(0) = 1 and F (0) = 0, we introduce a trial function by substituting which into Eq. (24) we can read off the critical chemical potential. In order to compare the analytical results with the numerical ones, we plot the critical chemical potential μ c versus the Horndeski parameter k in the form of curve in the left plot of Fig. 2 and list the results in Table 1. It is observed that the analytical results agree well with the numerical ones, which means that S-L method is still powerful in the superconductor with Horndeski correction. When the chemical potential is slightly above the critical point, the condensate O is small, so we can expand the gauge field φ(z) in the form of the small parameter O as Combining with Eq. (20), we have the equation of χ(z) as Defining T (z) = 1−z 4 z , Eq. (27) can be rewritten as At the boundary z → 0, the function χ(z) can be further expanded as Substituting Eq. (29) into Eq. (26) and thus comparing it with Eq. (10), we can get Obviously, the following main task is to calculate the values of χ (0) and χ(0). By using the boundary condition χ (0) = 0 [68,69], integrating Eq. (28) reads Further integrating the above equation, we can obtain the value of χ(0) as Combining Eq. (33) with Eq. (30), we can obtain the scalar condensate Considering Eq. (27) up to the zero order of O , we can obtain Thus we can rewrite the charge density as  Table 1 to compare with the numerical results.
It is obvious that the analytical results agree with the numerical ones at the same order, especially, the trend of the effect for Horndeski parameter on the coefficient is consistent with each other. What is more, the results in the case of m 2 = − 15 4 recover the ones in Refs. [29,35,68,69].

Conductor/superconductor phase transition
The authors in Ref. [54] studied the s-wave conductor/superconductor phase transition model with Horndeski correction in the probe limit by the numerical method and found that the increasing Horndeski parameter decreases the critical temperature and the phase transition is always second order. To check the reasonability of the numerical results, we will restudy the s-wave model by the analytical S-L method in this section. Firstly, we give the four-dimensional Schwarzshild AdS planer black hole as where r + denotes the location of event horizon corresponding to the Hawking temperature T = 3r + 4π . Considering the Lagrangian density (3) and the ansatz in terms of the scalar field and the gauge field (6), we can derive the equations of motion as where , a prime stands for the derivative with respect to r . Obviously, Eqs. (38) and (39) are the same to Eqs. (13) and (14) with κ = 0, c 1 = 0, c 2 = 0, ξ = 0 in Ref. [56]. It should be noted that by redefining g(r ) = r 2 f (r ), Eqs. (38) and (39) restore to Eqs. (23) and (24) in Ref. [54] and Eqs. (14) and (15) with τ = 0 in Ref. [57]. At the boundary r → ∞, ψ(r ) and φ(r ) behave as where ± = 3 2 ± 1 2 9 + 4m 2 1+3k consistent with Eq. (17) in Refs. [56,57], ψ 1 , ψ 2 , μ and ρ are all constants. According to the gauge/gravity duality, the explanations of the constants (ψ 1 , ψ 2 , μ and ρ) are the same with the ones for the insulator/superconductor model in section 2. In order to compare the current model with the ones in Ref. [54], we choose m 2 = −2 and restrict the Horndeski parameter in the range −0.01 ≤ k ≤ 1.
By means of the new variable z = r + r , Eqs. (38) and (39) can be rewritten as where the prime represents the derivative with respect to the variable z.
At the critical point, the scalar condensate vanishes, we can thus read off the solution of the gauge field as where r +c represents the location of the horizon at the critical temperature. Near the critical point, similar to the soliton case (20), we can express the scalar field ψ(z) as Combining Eqs. (43) and (44), Eq. (41) can be written as the S-L eigenvalue equation where the coefficients T , P and Q are as follows Thus λ 2 is obtained by minimizing the following expression as Given the boundary conditions F(0) = 1, F (0) = 0 [10,35,49], we take also the form of trial function (25) and obtain the value of λ. Thus the critical temperature can be written as The concrete results of the critical temperature are listed in Table 2, from which we can see clearly that the critical temperature decreases with the increasing Horndeski parameter k, which means that the larger Horndeski correction makes the conductor/superconductor phase transition more difficult to occur. Meanwhile, the analytical results agree well with the numerical ones in Ref. [54]. In particular, as k = 0, the result restores to the standard case in Refs. [7,8,10,35,49] Below (but close to) the critical temperature, the condensate O r + is very small, by using which we expand the gauge field (43) in the form of the small parameter as with the boundary condition χ(1) = 0 = χ (1). Integrating Eq. (50), we can obtain By further expanding Eq. (49) at the infinite boundary (z → 0) and comparing the linear order of z with Eq. (40), we can Take notice of the value of χ (0) in Eq. (52), the condensate can be expressed as The coefficient C 10 are listed in Table 2, from which we can find that the coefficient agrees with the numerical results in Ref. [54] at the same order, especially, the monotone decreasing trend as a function of the Horndeski parameter k.
For completeness, we also calculate the critical temperature of the conductor/superconductor phase transition in the  five-dimensional AdS black hole with the Horndeski correction by the S-L method in the probe limit. Concretely, in order for our results to be comparable with the ones in the previous literature, we take the mass squared of the scalar field as m 2 = −3, and thus the scaling dimension of the scalar operator reads ± = 2 ± 4 + m 2 1+6k = 2 ± 4 − 3 1+6k . The critical temperature for various value of the Horndeski parameter k is listed in Table 3.
It is shown that the critical temperature decreases with the increasing Horndeski correction, which is similar to the four-dimensional case. In particular, the critical value for the vanishing Horndeski parameter k restores to the one in Refs. [8,35]. Meanwhile, we also consider the critical behavior of the scalar condensate. It is found that the critical exponent of the condensate is still 1 2 , which means the second order phase transition at the critical point. What is more, we derive the coefficient (C 11 ) of the condensate near the critical tempera- and also list the value of C 11 in Table 3. It is observed that the coefficient increases with the Horndeski correction on the whole, which is different from the case in the four-dimensional black hole. Even though, the value with k = 0 agrees well with the analytical one in Ref. [35] and at the same order with the numerical one in Ref. [8].

Conclusions and discussions
In this paper, at the probe approximation, we have constructed the holographic s-wave insulator/superconductor model in the five-dimensional AdS soliton with the Horndeski correction via both numerical and analytical methods. Then, we realized analytically the corresponding conductor/superconductor model in the four-and five-dimensional AdS black hole. The effects of the Horndeski parameter k on the superconductor models were studied and the main conclusions of which are as follows.
In term of the insulator/superconducotor model, when the mass squared is fixed, the critical chemical potential increases with the larger Horndeski parameter k, which means that the increasing Horndeski correction hinders the superconductor phase transition. Meanwhile, the critical exponent of the condensate is 1 2 , which suggests that the system suffers from a second-order phase transition at the critical point, and thus be upheld by the behavior of the grand potential. What is more, near the critical point, the charge density increases linearly with the chemical potential, which is the university of holographic insulator/superconductor model. Furthermore, beyond the critical point, the imaginal part of conductivity displays an obvious pole at the low frequency, which implies that the system is indeed at superconducting state above the critical chemical potential. In addition, the stability of the hairy state is checked by the comparison of the grand potential with the no hair state. The analytical results such as the value of the critical chemical potential, the critical exponents of condensate O and charge density ρ agree with the numerical results, and the coefficients of O and ρ are qualitatively the same with the numerical ones near the critical point. As for the conductor/superconductor model, our analytical results showed that the increasing Horndeski parameter k hinders the superconductor phase transition. Especially, the analytical results in the four-dimensional case are consistent with the previous numerical results in Ref. [54]. Even though the present calculation is restricted to some special cases of the parameter space of mass squared and Horndeski correction, we can obtain the qualitatively same results for other values of Horndeski parameter k.
Comprehensively speaking, the increasing Horndeski correction hinders both superconducting phase transition models for fixed mass squared of scalar field, which is reasonable, because in this case the larger k will increases the effective mass of the scalar field and the scaling dimension of the operator O and thus prevent the instability occurring. On the contrary, if we do not fix the mass squared m 2 but the scaling dimension during the calculation (i.e., the m 2 1+6k is constant), by rescaling (1 + 6k)ψ 2 →ψ 2 , we can see clearly that the equations of motion for our insulator/superconductor model do not depend on the Horndeski correction, and thus the critical value and grand potential as well as the conductivity are independent of Horndeski correction. However, it should be noted that the phenonemon with fixed scaling dimension might only occur in the probe limit. In the near future, we will try to generalize the insulator/superconductor model from the probe limit to the backreaction case, and further study the effect of the Horndeski correction on the phase transition for fixed mass squared and scaling dimension , respectively, which will shed light on the understanding of the effect of Horndeski correction on our superconductor model.