New holographic Weyl superconductors in Lifshitz gravity

Jun-Wang Lu,Ya-Bo Wu,Huai-Fan Li,Hao Liao,Yong Zheng and Bao-Ping Dong

1 School of Physics and Electronics,Qiannan Normal University for Nationalities,Duyun 558000,China

2 CAS Key Laboratory of Theoretical Physics,Institute of Theoretical Physics,Chinese Academy of Sciences,Beijing 100190,China

3 Department of Physics,Liaoning Normal University,Dalian 116029,China

4 Department of Physics,Shanxi Datong University,Datong 037009,China

Abstract We build holographic p-wave conductor(insulator)/superconductor models via the numerical method with a new form of Weyl coupling in five-dimensional Lifshitz gravity,and then investigate how the Weyl coupling parameter γ and the Lifshitz scaling parameter z affect the superconductor models.In the conductor/superconductor model,an increase in the Weyl correction(Lifshitz scaling)enhances(inhibits)the superconductor phase transition.Meanwhile,both the Weyl correction(when the Lifshitz parameter is large enough and fixed)and the Lifshitz scaling suppress the growth of the real part of the conductivity.The Weyl correction used here (CB2) shows weaker effects on the critical value than the previous Weyl correction(CF2).In the insulator/superconductor model,larger vaules of the Weyl parameter hinder the formation of condensate.However,in increase in the Lifshitz scaling enhances the appearance of condensate.In addition,the calculation suggests that a competitive relation may exist between the Weyl correction and the Lifshitz scaling.

Keywords: Weyl correction,Lifshitz gravity,Holographic superconductor

1.Introduction

Since the gauge/gravity duality was proposed [1,2],it has been widely applied in various fields of modern physics involving strong interactions[3-9],and,in particular,to hightemperature superconductors [10,11].

In 2008,a holographic s-wave conductor/superconductor model was first realized by coupling Maxwell’s complex scalar field to Einstein’s gravity in the probe limit[10].It followed that the model was able to satisfyingly mimic the scalar condensate and infinite DC conductivity as the temperature decreased below the critical value.Subsequently,the holographic superconductor model was further investigated from the perspectives of the magnetic phenomenon[12],departure from the probe limit[13],and the analytical method [14].The results of these works showed that the holographic tool is indeed an effective and powerful way to probe high-temperature superconductors.

Subsequently,various holographic superconductor models were constructed.For example,to construct superconductors that were closer to the real world and to display more and more of the universal features observed in condensed matter,p-and d-wave conductor(insulator)/superconductor models,superfluids,and multiple order parameters as well as lattice effects were successfully simulated from the holographic perspective[5,15-27].To model the anisotropy of spacetime that appears in many real physics systems,a (d+2)-dimensional nonrelativistic gravity with a Lifshitz fixed point was proposed[28,29],which was then generalized to the finite temperature case(2).Based on the Lifshitz background(2),s- and p-wave conductor(insulator)/superconductor models were constructed in [30-35].The results showed that for both superconductor models,stronger Lifshitz scaling hinders the superconductor phase transition and suppresses the growth of frequencydependent conductivity.In addition,via a Maxwell-complexvector (MCV) field,the authors of [36] realized a magneticfield-induced condensate,which was confirmed to be a massive generalized version of the SU(2) p-wave model [37-39].Following this,[22,24,40-48] described the construction of electric-field-induced p-wave superconductor models.It followed that the MCV superconductor model displays abundant phase structure in the presence of a back reaction.However,the stronger Lifshitz scaling still makes the superconductor phase transition more difficult.

On the other hand,to broaden the applied range of gauge/gravity duality,many works have focused on the effects of the(the constant λ denotes the ’t Hooft coupling) correction on holographic superconductor models.In particular,these investigations usually include nonlinear electrodynamic corrections [49,50] and the higher curvature corrections [38,51-54] which were found to inhibit the phase transition.In addition,there is another type of correction which involves the combination of the curvature tensor/scalar(or Weyl tensorCμνρσ)and the Maxwell field strength Fμν,for instance,the Weyl form CF2[55] and its generalized high-derivative version CnF2[56] as well as the RF2form [42].The related superconductor models implied that the MCV p-wave insulator/superconductor model(in particular,the critical chemical potential) is independent of the Weyl correction(CF2).However,the four derivative Weyl(CF2) as well as the six derivative Weyl(C2F2) corrections always promote both the swave and MCV p-wave conductor/superconductor models[57-59].In addition,as one comprehensively takes into account the Lifshitz scaling and the Weyl correction in the MCV superconductor model,both corrections exhibit interesting competitive effects on the critical value and the condensate [60].In fact,from the viewpoint of effective field theory,apart from the familiar Weyl correction CF2,there are also many other kinds of possible coupling form,such as the couplingwhere the antisymmetric tensor Bμνis made up of the MCV field ρμ;for details,see the next section.For simplicity,we call this new Weyl correction CB2,to distinguish it from the previous Weyl correction CF2.Bearing in mind the interesting CF2Weyl effects as well as its competing effects on the Lifshitz scaling of the MCV superconductor[60],in this paper,we will systematically investigate how the CB2Weyl correction influences the MCV p-wave conductor(insulator)/superconductor models in the Lifshitz black hole(soliton) background.It follows that for some parameter spaces,the influences of the Weyl parameter and the Lifshitz parameter on the conductor/superconductor model are in contrast to their influences on the insulator/superconductor phase transition.

The organization of the rest of this paper is as follows.In section II,we numerically construct a p-wave conductor/superconductor model.Based on a superconducting background,frequency-dependent conductivity is studied.Following a similar procedure,an insulator/superconductor model is built in section III.Conclusions and discussions are summarized in the final section.

2.Conductor/superconductor model

To study the Weyl and Lifshtiz effects on the p-wave conductor/superconductor phase transition,we first introduce the MCV field with the Weyl coupling term CB2in five-dimensional Lifshitz gravity.The corresponding full action,consisting of the gravitational(Lg) and matter (Lm) parts reads[28,29,36,40,41,56,61]

where φ,FμνinLgstand for the massless scalar field and an Abelian gauge field where the parameter ι is related to the Lifshitz dynamic exponent z.The gravitational partLgadmits the so-called five-dimensional Lifshitz black hole solution as

where r+is the location of the event horizon and the Lifshitz parameter z characterizes the degree of anisotropy of spacetime by the scaling symmetryMeanwhile,the Hawking temperature is given byThe solution of (2) returns standard anti-de Sitter (AdS)spacetime in the case of z=1.In addition,the matter partLmthat describes the vector condensate contains a Maxwell field strength Fμν=?μAν-?νAμand a vector field ρμcoupled to the Weyl tensorμνρσC.Specifically,the second-order tensor Bμν=Dμρν-Dνρμand Dμ=?μ-iqAμ.m and q are regarded as the mass and the charge of the vector field ρμ,respectively,and also an identity matrixWithout losing generality,we will take L=1.In addition,we construct our model in the probe limit which is still believed to maintain the main characters of superconductors and simplifies the calculation [10,11,15,36].

By a variation of action (1) with respect to the matter fields ρμand Aμ,respectively,we can write the equations of motion as follows:

The above equations have the same form as the equations in[36,40,41,44] when γ=0 and z=1 and the equations in[48] when γ=0 and z ≠1.

Figure 1.The condensate as a function of the temperature corresponding t o(black solid),0(red dashes),(blue dots and dashes)for z=1(left) and (right).

To build the p-wave superconductor,we turn on the components of the matter fields as described in [36,40,41],

substituting this into equations (3) and (4) yields the coupled differential equations for ψx(r) and φ(r).The familiar approach for obtaining the condensate is to solve the differential equations by the shooting method and the imposition of suitable boundary conditions [10,11,13,15,17,36].Specifically,at the horizon(r=r+),we require ψx(r+) to be regular and φ(r+)to be zero to obtain the finite norm of Aμ.At infinity (r →∞),the asymptotical expansions of ψx(r) and φ(r) read

Plots of the vector condensate curves forandare shown in the left(z=1)and rightpanels of figure 1,from which we observe that a critical temperature exists for all cases.Below the critical value,the vector condensate starts to grow.Fitting the curves near the critical point,we find there is always a square-root behavior,such thatwhich implies that the system may suffer from a second-order phase transition.Meanwhile,when the temperature is decreased,the condensate reaches a saturated value at z=1 but continuously increases forwhich seems to be universal for Lifshitz superconductors[33].What is more,at the low temperaturethe condensate decreases for z=1 but increases forwith larger values of γ.By probing the Lifshitz parameter carefully,the location of the demarcation for the qualitatively different behaviors of the condensate is at aboutwhere the condensate curves are confirmed to be almost the same for different Weyl parametersComparing the effects of the new Weyl coupling,CB2,with the one of the previous Weyl couplings,CF2,on the vector condensate,we find that the effect of the new coupling,CB2,is weaker than that of the old Weyl coupling CF2,[60]based on the left panel in figure 1.Furthermore,when γ=0,the present model is the same as the one in[60]with a vanishing CF2correction and the model in [42] without an RF2correction.

In the left panel of figure 2,we display the critical temperature with respect to the Weyl parameter γ for different Lifshitz parameters z,regarding which,we have the following remarks.First of all,when the Weyl parameter γ is fixed and the Lifshitz parameter z increases,the critical temperature clearly decreases,which implies that the improving anisotropy of spacetime hinders the appearance of the vector condensate.The suppressing effect of the Lifshitz scaling on the conductor/superconductor models was also observed in[33,35,42,48].Secondly,when the Lifshitz parameter z is fixed and the Weyl correction parameter γ gradually increases,the critical temperature increases monotonously,which suggests that the larger the Weyl correction,the more easily the condensate emerges in the conductor/superconductor model.In addition,it is worth mentioning that the critical value increases very slowly with the larger Weyl coupling parameter γ,which again agrees with the weaker influence of the Weyl correction on the condensate in the left panel of figure 1 compared with the left panel of figure 3 in [60].In addition,for γ=0,we read off the critical temperatures0.07958 (z=1) ,0.061330.04266 (z=2)and 0.02244which agree with the corresponding values in the left panel of figure 3 and Tab.1 in [60].Last but not least,rescaling the value of ψ2,μ as well as Tcby the charge density ρ,we can translate the system from a grand canonical ensemble with a fixed chemical potential to a canonical ensemble with a fixed charge density; correspondingly,the critical temperatureis rewritten asfor the case of z=1 and γ=0,which obviously restores to the value in [58] without a CF2correction.

Next,we check the thermodynamic stability as well as the order of the phase transition by calculating the grand potential density of the system,where Sosis the onshell action,and V3=∫d3x.The on-shell action,given by the Minkowski action(1),is

where we have taken ∫d3x=V3,into consideration and ignored the factorObviously,both the Weyl parameter γ and the Lifshitz parameter z contribute to the grand potential.

We representatively plot the grand potential versus temperature withz=andγ=in the right-hand panel of figure 2.It follows that at the critical point,the two curves merge smoothly together but not in the form of a swallowtail,which matches the condensed behaviors in figure 1 and,in particular,implies that the phase transition is indeed second order.In addition,the value of the black dashed curve is always larger than that of the red solid curve,which verifies that the hairy state is thermodynamically stable below the critical point.Besides,we have also calculated the behaviors of the grand potential for the other Weyl cases(such asas well as the other Lifshitz casesand obtained similar results to the one in figure 2.Therefore,we believe that our numerical results are quite reliable in the parameter space

In what follows,we take a look at the signals of superconductors,i.e.,their infinite DC conductivity,which requires us to calculate the perturbation of the Maxwell field in the bulk from the AdS/conformal field theory (CFT) correspondence.Taking the perturbation along the y direction with the ansatz δAy(t,r)dy=Ay(r)e-iωtdy,and substituting it into equation(4)yields the linearized differential equation of Ay(r)as well as its general falloff as

where A(i),and η are both constants.From equations (1)and(8),the retarded Green’s function reads

where the logarithmic divergent term for z=1 is canceled by the holographic renormalization [11,33,42,60].By means of the Kubo formula,the frequency-dependent conductivity is given byσ(ω) =from which we can see that the key point in the calculation of the conductivity is to extract the coefficients A(0)and A(z+1)in equation (8).The effective method is still to numerically solve the differential equation of Ay(r) with the ingoing boundary condition at the horizon asAy(r)=(r-r+)-iω(3+z)(1 +Ay1(r-r+) +Ay2(r-r+)2+Ay3(r-r+)3+…).

We show a typical plot of the AC conductivity atfor the fixed Weyl parameter(left) and the Lifshitz parameter(right),respectively,in figure 3.There are clear poles for the imaginary part of the conductivity Im[σ] at the lower frequencies in all cases.From the Kramers-Kronig relation,these poles correspond to a delta function of the real part of the conductivity Re[σ] at the vanishing frequency and thus indicate infinite DC conductivity.Moreover,Re[σ]increases with frequency and diverges at large frequencies,which is also its universal character in five-dimensional spacetime[11,33,42,60].Furthermore,in the case of z=1,there is a minimum for Im[σ]in the intermediate frequency region,which corresponds to the energy gap of the superconductor ωg.After the calculations,we find that the larger the Weyl correction,the smaller the value offor instance,Meanwhile,the variation ofwith respect to the Weyl parameter is very slow,which corresponds with the weak Weyl effects on the condensate.As we know,the location of the minimum of Im[σ]corresponds to the frequency region where Re[σ] increases quickly with frequency.As a result,for z=1,the phenomenon whereby the energy gap moves leftwards with increasing γ corresponding to Re[σ] increases more quickly with larger γ.Interestingly,for a big enough Lifshitz parameter z,such ason the right-hand panel of figure 3,we find that Re[σ] increases more and more slowly with increases in γ.These different effects of the Lifshitz parameter z on conductivity seem to be consistent with its effects on the condensate shown in figure 1.In addition,as displayed in the left-hand panel of figure 3,Re[σ] is suppressed more and more obviously with an increase in the Lifshitz scaling,which is also observed in[33,42,60].In addition,we also calculated the conductivity for other cases with fixed γ,for example,0and qualitatively obtained the same results as the ones shown in the left-hand panel of figure 3.

3.Insulator/superconductor model

First of all,by double Wick rotating the Lifshitz black hole(2),i.e.,t →iχ and w →it,we read off the five-dimensional Lifshitz soliton

where r0represents the tip f(r0)=0.In order to guarantee smoothing of the geometry,we have to impose Scherk-Schwarz periodicityon the spatial direction χ.As argued in[17],the Lifshitz soliton solution has a negative mass gap,which corresponds to a dual field theory with a mass gap.Due to a similar phenomenon whereby the mass gap also exists in the insulator,the current soliton solution is believed to provide a gravitational background to construct the holographic insulator/superconductor model [17,33,41,58].

Having the Lifshitz soliton background in hand,to build the MCV p-wave superconductor model with the CB2Weyl correction,we take the Lagrangian density of the matter field Lmto be equal to the one in equation (1).Considering the same ansatzs of the matter field(5),we can read off the coupled differential equations of motion of ψx(r) and φ(r) in the Lifshitz soliton geometry(10).To numerically solve the differential equations of motion,we add the Neumann-like boundary condition to ψx(r) and φ(r) at the tip; meanwhile the general falloff of ψx(r) is the same as the falloff in equation (6),and the asymptotic solution of φ(r) reads

Moreover,the physical translations of the coefficients(ψ1,ψ2,μ and ρ)are the same as the ones in the Lifshitz black hole.In addition,to conveniently compare the models with different parameters,we fix the periodicity Γ=π by rescalling the values of ψ2,μ,and ρ.

After a series of calculations,we display the condensate versus the chemical potential in figure 4 with z=1 (on the left),(on the right),respectively,and the left-hand panel in figure 5 withMeanwhile,we show the dependence of charge density on the chemical potential within the right-hand panel of figure 5.It follows that a critical chemical potential exists for all cases.Above the critical value,the condensate and the charge density appear and then grow monotonously with increasing chemical potential.In particular,near the critical value,all condensates behave according toThe critical exponentindicates that the system suffers from a second-order phase transition [17,41].Meanwhile,the charge density increases linearly,which seems to be a universal property at the probe approximation[17,33].What is more,as the Weyl parameter γ increases for z=1,the critical chemical potential moves towards the right,which implies that the improving Weyl correction inhibits the insulator/superconductor phase transition.On the other hand,the critical chemical value moves towards the left with increasing values of γ forwhich suggests the larger the Weyl parameter γ,the easier the phase transition of the system.Obviously,the effects of the Weyl correction on the insulator/superconductor phase transition are quite different from its effects on the conductor/superconductor phase transition.In addition,we observe that the location of the critical chemical potential moves towards the left with improving Lifshitz scaling forwhich suggests that stronger Lifshitz scaling allows the insulator/superconductor phase transition to occur more easily.

Motivated by the above abundant phenomena arising from the Lifshitz and Weyl parameters,we display the critical chemical potential as a function of the Weyl correction with different Lifshitz scaling in figure 6.It can be observed that larger values of γ inhibit the superconductor phase transition for z=1 andand promote the superconductor phase transition for z=2 andMeanwhile,the critical chemical potential varies almost linearly with the Weyl parameter.By calculation,we find that in thecase of the Lifshitz parameter,the critical value has nearly no dependence on the Weyl correction.In addition,since the Weyl correction is sufficiently large and fixed,such thatthe critical value decreases with an increase in the Lifshitz parameter,which means that larger Lifshitz scaling promotes the phase transition.In summary,the Weyl correction and the Lifshitz scaling exhibit competitive behavior in the superconductor model.

Similarly to the procedure in section 2,we calculate the grand potential densityfrom the defined ‘temperature’ of the soliton system and the on-shell action.The concrete on-shell action reads

which is quite different from the one given in(7)for the black hole.We representatively show the density of the grand potential forin figure 6.Obviously,as the chemical potential gradually increases from the critical point,all curves with hairy vectorsdevelop smoothly from the horizontal dashed line characterized by no hair.Furthermore,except for the critical point,all hairy curves are always lower than the horizontal line,but do not intersect it,since μ ＞μc.These characters are not only consistent with the condensed curves in figures 4 and 5,but also indicate that a second-order phase transition occurs at the critical point and the hairy state is thermodynamically stable,which are expected in the superconducting state.

To further check that the hairy state is a superconducting state,it is useful to calculate the conductivity from the investigation of gauge perturbation δAy(t,r)=Ay(r)e-iωtin the hairy bulk.Following the same procedure in the conductor/superconductor model,we read off the linearized differential equation for Ay(r) and further obtain its asymptotic solution as

the form of which is the same as that shown in [33,42].After analysis,we find the concrete expression of the retarded Green’s function related equation(13)is the same as the one for the black hole(9),so that the AC conductivity is again formally equal to the one in the conductor/superconductor model.Given the Neumann-like boundary conditionAy(r) = 1+Ay1(r-r0) +Ay2(r-r0)2+Ay3(r-r0)3+…,we can obtain the imaginal part of the conductivity Im[σ] by solving the differential equation for Ay(r).The frequency-dependent conductivity Im[σ]with z=1 andis plotted in figure 7,from which we observe the pole in the low-frequency region.As discussed in the conductor/superconductor model,this pole corresponds to the delta function for Re[σ] and thus indicates infinite DC conductivity,which not only closely matches the behaviors of the condensate and the grand potential,but also provides a signal for the superconductor [17,33].What is more,when z=1,we can see clearly that the curves for the conductivity nearly overlap with each other forThis means that when the chemical potential is sufficiently large,the effects of the CB2Weyl on the conductivity are rather weak,which is consistent with the phenomenon in figure 4,where the condensate curves for different values of γ tend to coincide at the large chemical potential.Moreover,for a sufficiently large Weyl parameter(such asnot displayed)or Lifshtz scaling,for example,shown on the right panel of figure 7,the location of the second pole always moves towards the right with the other parameter,z or γ,which seems to agree with the increasing differences of condensates with different parameters in figures 4 and 5.

4.Conclusions and discussions

In the probe limit,we have numerically built the MCV pwave conductor/superconductor and insulator/superconductor phase transitions with CB2Weyl coupling in fivedimensional Lifshitz gravity.We mainly studied the influences of the Weyl(γ) and Lifshitz(z) parameters on the superconductor models.The main results are summarized as follows.

For the conductor/superconductor model,as the Weyl parameter γ increases or the Lifshitz parameter z decreases,the critical temperature becomes larger and larger,which suggests that the phase transition becomes easier and easier with an increase in the Weyl correction or a weakening of the Lifshitz scaling.At low temperatures,as the Weyl parameter increases,the condensate value decreases(increases)below(above) the point of demarcationIn addition,we find that Re[σ] grows more and more quickly (slowly)with frequency for smaller(larger) values of z in the intermediate region,which is consistent with the behaviors of the condensate with the corresponding Lifshitz parameter z.Compared with the CF2Weyl correction [60],the CB2Weyl correction used here has weaker influences on the condensate and the critical temperature.

For the insulator/superconductor model,as the Weyl parameter increases,the phase transition becomes increasingly difficult(easy) below(above) the point of demarcationMeanwhile,when the Weyl correction is sufficiently strong,the superconductor phase transition becomes increasingly easy with an improvement of the Lifshitz parameter z.On the other hand,compared with the standard case,i.e.,z=1 and γ=0,we can understand the above results as showing that larger Weyl corrections inhibit the phase transition,while an increase in the Lifshitz scaling enhances the phase transition.Meanwhile,the abundant effects of the Weyl correction and the Lifshitz scaling on the superconductor result from the competition between the two corrections.

In fact,the phenomenon whereby the insulator/superconductor phase transition becomes easier with larger Lifshitz parameters is not strange,and has been observed in the swave and MCV p-wave insulator/superconductor models in[33,42].In particular,if we redefine the periodicity as Γ=π(with a scaling factor ofor even the Γ=2π(with a scaling factor ofused in [42],we can find the critical value without RF2correction is the same as the one in the model used here with γ=0.Meanwhile,for any Weyl and Lifshitz parameters in the rangeandin the conductor/superconductor model,Re[σ] is always positive though it is very small at low frequencies below the critical temperature,which means that the current superconductor model satisfies the unitary requirement for the present parameter space.Even so,it is meaningful to check whether phenomena such as superluminal velocities or the violation of causality will happen or not for the current Weyl parameter space.In addition,we also calculated the frequency-dependent conductivity near the critical temperature,but did not see the pronounced peak and/or the Drude-like peak that appeared in the case with C2F2correction [35,61,62],so would be interesting to investigate how the sixth-order derivative term C2B2influences the model.In the near future,we will focus on solving the above problems to improve the understanding of the impact of the Weyl correction on the superconductor model.

Acknowledgments

We would like to thank Prof.L.Li for his helpful discussion.This work is supported in part by NSFC (Grant Nos.11865012,12 075 109,12 075 143 and 11 747 615),the Foundation of Scientific Innovative Research Team of Education Department of Guizhou Province (201329).

ORCID iDs