We have developed the recent investigations on the second-order phase transition in the holographic superconductor using the probe limit for a nonlinear Maxwell field strength coupled to a massless scalar field. By analytical methods, based on the variational Sturm-Liouville minimization technique, we study the effects of the spacetime dimension and the nonlinearity parameter on the critical temperature and the scalar condensation of the dual operators on the boundary. Further, as a motivated result, we analytically deduce the DC conductivity in the low and zero temperatures regime. Especially in the zero temperature limit and in two dimensional toy model, we thoroughly compute the conductivity analytically. Our work clarifies more features of the holographic superconductors both in different space dimensions and on the effect of the nonlinearity in Maxwell's strength field.
1. Introduction
In the recent years, using the holographic picture of the world, the AdS/CFT (anti de Sitter/conformal field theory) correspondence [1–3] has been applied to study some strongly correlated systems in condensed matter physics, especially for strongly coupled systems with the scale-invariance. Particularly, people studied the low temperature, quantum critical systems near critical point (see, e.g., [4, 5] and references therein). The critical phenomena, which happen here, is a second-order phase transition from normal phase to the superconducting phase, in which below a specific temperature Tc, the DC conductivity becomes infinite. Such second-order phase transitions happen in the high-temperature superconductors and can be described very well by the AdS/CFT dictionary [6, 7]. From the classical and phenomenological point of view, superconductivity, in the high-temperature type II superconductors, modeled by a phenomenological Landau-Ginzburg Lagrangian. This Lagrangian contains a general complex value scalar field Ψ, plays the role of a condensate in a superconductive phase. Basically, to have a scalar condensation in the boundary quantum field theory using CFT on the boundary of the bulk, Hartnoll et al. [8] introduced a U(1) abelian gauge field Aμ and a typical conformal coupled charged complex scalar field in the bulk black hole background. The conformal mass is above the Breitenlohner-Freedman (BF) bound [9]. To solve the negative mass problem, finally Gubser [10] showed that the vector potential modifies the mass term of scalar field and we have a possibility to have hairy black holes in some parts of the parameter space.
The full description of the superconductivity in the probe limit or away this limit needs to provide the numerical solutions of a couple of nonlinear differential equations. By simplicity, they can be solved using the shooting approach by expanding in series the functions and matching these by varying the free parameters of the series in a typical point between the horizon and the spatial infinity. Parallel to the numerical studies, recently some analytical approaches have been proposed to find the universal properties of second-order phase transitions in holographic superconductors [11–18]. In particular, the authors in [18] used the variational functional method on the real valued functions for the Sturm-Liouville eigenvalue problem to compute analytically the basic properties of the holographic superconductors in a three-dimensional boundary using the solutions of the bulk gauge fields in the four dimensional bulk theory. In [18, 19] the authors have shown that also it is possible to calculate critical exponent and critical temperature by this approach. The eigenvalue of this variational problem is a function of the critical chemical potential μc, and consequently, it is related to the Tc. Different modes of super criticality s-wave, p-wave and d-wave have been studied [20]. Also, one can apply this method to superconductors with external magnetic fields, which gives the analytical results in good agreement with numerical results produced previously [21, 22].
Furthermore, a number of aspects of external and bulk magnetic fields in holographic superconductors have been investigated [23–26]. The phase transition can be interpreted in terms of the string interactions [27–31]. The effects of the nonlinear electrodynamics in the holographic superconductors have been investigated recently [32–36].
There are many interests in the modified gravity theories. For example, on Gauss-Bonnet and Weyl corrected superconductors, on which, we are working with a higher derivative corrected bulk black hole, like Weyl corrections [37] numerically. Furthermore, we have studied the Weyl corrections to the superconductors analytically [38, 39]. Moreover, we showed that there exists a family of p-wave holographic with Weyl corrections [40].
In the present paper, we would like to study the ((d-2)+1)-dimensional holographic superconductor in the probe limit for a power law Maxwell field strength (FμνFμν)δ coupled to a scalar field. We focus just on the s-wave cases. We must clarify the motivation of the s-wave approximation in holographic models of superconductors. In the relativistic models of the gravity, it is highly known that s-wave approximation is not a good approximation, for example, in the cosmological models and black holes [41]. The meaning of the s-wave here does not back to the reduction of the action from four to two dimensional like the dilatonic action from the four dimensional spinor (Majorana) action. We mean by s-wave, in the context of holographic superconductors, a scalar order parameter, whose expectation value breaks the U(1) but not rotational symmetry. Moreover, we can have the Yang-Mills fields with SU(2) symmetry which additionally they can generate another symmetry breaking of an axial vector type. The last case resembles the p-wave models. We mention here that the three dimensional non-linear model of the superconductors, which we used in this paper, is a realistic model and it will be more interesting that we can find a direct relation between this nonlinear model and the results of a higher dimensional model, by a principle like the detailed balance.
Another additional point is to restrict ourselves just to the case of a single horizon. The problem of the multihorizon cases needs more investigation, for example, the case of the Nariai black holes. The theory here will be so different. This later appeared in the lower dimensional models. For example, the case of the quantum corrected BTZ like black hole is a good example [42]. In the holographic set up for superconductors one must identify a temperature in his gravitational bulk model to the CFT temperature on the boundary. If the black hole has only one horizon, in this case, we can use the Hawking-Bekenstein (horizon) or Kodama- Hayward temperature [43] as a reasonable candidate. But if our asymptotically AdS bulk has more than one horizon, for example, in the case of the charged BTZ like black holes, then we take the temperature of the real physical horizon (the temperature which is obtained by calculation the surface gravity of the biggest null hypersurface orthogonal surface) as the candidate for temperature of the CFT. In fact the effects of the quantum corrections and charged Maxwell field on the background of the bulk are very interesting problems and can be investigated in more details. Also it is possible to relate the instability of such charged dilaton configurations in the AdS spacetime [44] to the symmetry breaking mechanism of the superconductors. The idea has motivation enough as a new work. In this paper we investigate analytically the effect of the spacetime dimension d and the power δ on the critical temperature Tc. Although our problem is the especial massless case of the model which has been investigated recently [45], we study these corrections to the superconductors analytically. Additionally, we want to compute the DC conductivity for this kind of the superconductor using the perturbation method. In this approach, we apply an external linear electromagnetic field. This field is periodic in time. By calculating, the response in the first order linear approximation, we compute the conductivity for the low temperature case, especially for zero temperature configuration.
Our plan in this paper is as the following. In Section 2, we clarify our motivation for considering the nonlinear Maxwell action instead of the linear theory. In Section 3, we introduce our model for holographic superconductors. In Section 4, we apply the variational method to obtain the critical temperature of the system. In Section 5, we calculate the critical exponent for the condensation operator. In Sections 6 and 7 we compute the conductivity for low and zero temperature cases. We summarize and conclude in the final section.
2. Motivation for Nonlinear Maxwell Effects in Holographic Superconductors
The linear approximations in the mathematical physics, as we know, have limitations, both in predictions of the model and especially on matching with the full description of the model using the numerical results. The linear Maxwell theory fails in some domains, and it is needed to consider the general form of the Lagrangian instead of the linear one. The Lagrangian of the Maxwell model is
(1)ℒM=-14F2,
where F2=FμνFμν. A natural extension of (1) is obtained by replacing a general function of F in the form Φ(F). Recently, it has been shown that such nonlinear general forms have a rich family of black holes in f(R) gravity [46]. There are different reasons for investigating these forms. The oldest one may be the Born-Infeld (BI) alternative for linear Maxwell's theory. In string theory language, this BI action can be replaced by the tachyonic action. The BI Lagrangian reads [47]
(2)ℒBI=η2(1-1+F22η2).
Here η is the string's tension parameter. This form reduces to (1) in the limit of η→∞. Indeed, we can expand in series (2) in the following form:
(3)ℒBI=∑n=1∞cnη2(1-n)F2n.
The leading order term n=1 is just in the form of (1). Even, if we do not work with BI theory, this nonlinearity meets us from a geometrical point of view. Suppose that we want to write a conformal invariance (CI) Lagrangian, constructed from the U(1) gauge fields Aμ. Such CI is the invariance of the whole theory under geometrical transformation gμν→e2σgμν in a d≥4 Riemannian manifold without torsion or nonmetricity fields. However, previously, a version of such Lagrangian has been found in Weitznbock spacetime with torsion and nonmetricity fields [48]. It is easy to show that the Lagrangian ℒ∝Fd/2 is invariant under CI transformations. When d=4, the proper Lagrangian is ℒ∝F2 but in d=5 the suitable form is ℒ∝F5/2. The last form belongs to the nonlinear, noninteger Maxwell family. Thus from geometrical view, the nonlinearity is welcome in our Lagrangian dynamical theory. Further, even if we do not know any on BI or CI, when we are working with vacuum effects, there is a simple generalization of Maxwell Lagrangian in a logarithmic form
(4)ℒlog=-η2log(1+F24η2).
As BI case, in the limit of |F|≪η, by expanding this equation the leading order term n=1, is the linear theory (1). In brief, according to the above discussion, it seems that consideration of the nonlinear effects of the Maxwell field as Fδ is important, and the significant differences between the usual linear theory F2 and nonlinear theory Fδ can be shown. Holographic superconductors provide a rich background for testing such type of new physics. In this paper, we will describe a (d-1) dimensional holographic superconductor (HSC) via d-dimensional gravity dual, described by a d-dimensional AdS black hole on the static patch. We set the Stuckelberg field to zero and work with massless scalar fields.
3. Field Equations
We write the following action [32–36] for a power Maxwell field Fμν, which it is coupled minimally to a massless scaler field ψ in a d-dimensional asymptotic AdSd spacetime [45]
(5)S=∫ddx-g[R+(d-1)(d-2)-ξ(FμνFμν)δ-|Dμψ|2).
Here, we set the radius of the AdS, L=1, Dμ=∂μ-ieAμ, Fμν≡2∂[μAν]=∂μAν-∂νAμ, Aμ is U(1) vector potential. The exponent (power) δ can be noninteger as we explained it in the previous section. Further, when we work with Maxwell linear electrodynamics δ=1, d=4, then ξ=1/4. However, in our case with δ≠1, in general, ξ∈ℛ remains as a free parameter in our model. We set the electric charge e=1. In normal phase and in the absence of the scalar field we put ψ=0. This action has been used before in literature for full description of the nonlinear Maxwell field in the gravitational action. In the probe limit when the matter action and the gravitational part decouple, the system of the field equations has an exact solution which will be discussed here. The solution of the generalized Maxwell-Einstein equations is the simple d-dimensional static Reissner-Nordstrom-Anti-de Sitter (RNAdS) black hole with the following metric form:
(6)gμν=diag(-f(r),1f(r),r2Σd-2),
where Σd-2 is the metric on a (d-2) dimensional sphere and the metric function is
(7)f(r)=r2(1-(r+r)d-1).
Here r+ is the black hole horizon which in general is the largest root of the algebraic equation f(r+)=0. If we set δ=1 in (5), we recover the usual holographic superconductors. In the probe limit by ignoring the backreaction effects of the matter fields in the background metric gμν, and with spherically symmetric static metric and by choosing a suitable gauge fixing for the gauge field Aμ we can take the functions (ψ,ϕ)∈ℛ and one variable functions. We assume that the functions ψ and ϕ have finite numbers of poles on the real axis. It means the analytical solutions can be written in the forms of the algebraic series expressions. To obtain the field equations, we assume that Aμ=ϕ(r)δμt, so the nonzero components of the Fμν read as
(8)Frt=-Ftr=ϕ′.
Hence, we have
(9)FμνFμν=-2ϕ′2.
Also, the Ricci R of a d-dimensional spacetime with spherical symmetry reads
(10)R=-(f′′+2(d-2)f′r-(d-2)(d-3)fr2).
Further, we compute
(11)|Dμψ|2=fψ′2-ϕ′2f.
Finally, by plugging the above expressions into the action (5), making a partial integration, we get the following effective Lagrangian:
(12)ℒ=rd-2[-(f′′+2(d-2)f′r-(d-2)(d-3)fr2)+(d-1)(d-2)-ξ(-2ϕ′2)δ-(fψ′2-ϕ′2f)].
As a first step, it is necessary to eliminate the f′′ term, by integration part by part. After it, to write the field equation, we use the Euler-Lagrange equation as the following:
(13)ddr(∂ℒ∂q,r)=∂ℒ∂q,q={ϕ,ψ}.
The field equations, derived from (12), read
(14)ψ′′+(f′f+d-2r)ψ′+ϕ2f2ψ=0,(15)ϕ′′+d-22δ-1ϕ′r-Cδψ2ξfϕ(ϕ′)2(1-δ)=0,
here Cδ=((-2)2-δδ(2δ-1))-1. To avoid the pure complex numbers in our field equations, we assume that δ≠2-(1/2N). The case with δ=1, D=4 is the usual four-dimensional HSC which describes the three dimensional superconductors.
These field equations are the special massless case of the model which has been investigated recently [45, 49]. To avoid the divergence near the singularity f(r+)=0 we write the boundary conditions for (14), (15) by
(16)ϕ(r+)=ψ′(r+)=0.
The asymptotic solutions for system (14), (15) on the AdS boundary r→∞, are
(17)ψ=D-+D+rd-1,(18)ϕ=A-Brη,
in (18), D±=〈𝒪±〉, η=(d-2δ-1)/(2δ-1), where 〈𝒪±〉 denotes the vacuum expectation value of the dual operator 𝒪± on the boundary, and the chemical potential and charge density of the dual theory are A=μ, B=ρ(2δ-1)-1, respectively.
We must clarify the reason for the modification of the form of the electric potential ϕ modified. In the asymptotic regime, we know that the metric function f behaves like f~r2. Also, the scalar field has the following asymptotic form ψ~0, so (15) gives us
(19)ϕ∞′′~-d-22δ-1ϕ∞′r.
The solution (19) reads
(20)ϕ∞(r)=c1+c0(2δ-1)2δ-d+11r(d-2δ-1)/(2δ-1).
This solution coincides completely with the solution presented in (18). We mention here that the above function ϕ∞(r) in the limit of the linear electrodynamic theory δ=1 has the true asymptotic form of ϕ∞(r)~r3-d, d≠2. For d=2 the expression of ϕ∞(r) is in the form of a diverging log term ϕ∞(r)~log(r) and the application of the AdS/CFT fails. At least, we do not know the unique and true dictionary of the AdS/CFT in this lower dimensional bulk theory.
The only thing left is to identify the parameters with the physical quantities in the dual theory, that is, the chemical potential μ and the charge density ρ. We did it before, so the asymptotic behaviors have the same forms.
The asymptotic solutions for ϕ, ψ are the same as the previous expressions which have been presented in [45]. For normalization purposes, we set D-=0.
4. Variational Method
To solve the solutions of the field equations given by (14) and (15), the well-known technique is numerical algorithms. However, from these numerical solutions, it is not so easy and straightforward to read 〈𝒪〉. Another method is using the matching method. It is a potentially powerful method. Even so, the results must be interpreted very carefully near the boundaries. The first step for solving system (14), (15) using variational approach [18], is rewriting the equations in a new dimensionless coordinate z=r+/r in the following forms:
(21)ψ′′+(f′f-d-4z)ψ′+(r+z2)2ϕ2f2ψ=0,(22)ϕ′′-η-1zϕ′-Cδr+2δξz4δψ2fϕ(ϕ′)2(1-δ)=0.
Now, the prime is for differentiation with respect to the dimensionless radial coordinate z. Near the critical point T=Tc, the following solution is valid:
(23)ϕ≈ϕB(1-zη).
Here ϕB=μc=ρ(2δ-1)-1/(r+)η+1 is the value of the ϕ at the horizon r=r+,δ≠(d-1)/2 and the critical point corresponds by the critical chemical potential μ=μc. Further, we can write
(24)ψ≈〈O+〉(zr+)d-1F(z),
near the AdS boundary z→0 with F(0)=1, F′(0)=0. The trial variational function satisfies the following second-order Sturm-Liouville self sdjoint differential equation:
(25)[μ(z)F′(z)]′-Q(z)F(z)+(ϕBr+)2P(z)F(z)=0,
where
(26)μ(z)=z2d-2(zd-1-1),P(z)=μ(z)(1-zη1-zd-1)2,Q(z)=-μ(z)[(2+(d-3)zd-1z(1-zd-1)+d-4z)(d-1)(d-2)z2-d-1z-μ(z)×(2+(d-3)zd-1z(1-zd-1)+d-4z)],
our strategy is to obtain the minimum value of ϕB/r+ from the minimization of the following functional:
(27)Θ(δ,d)≡(ϕBr+)Min2=∫01(μ(z)F′(z)2+Q(z)F2(z))dz∫01P(z)F2(z)dz.
The minimum of the critical temperature Tc is obtained from the T=(d-1)r+/4π. It reads
(28)Tc=γρη+1,
and γ=((d-1)/4π)((ϕB/r+)Min)-(η+1). Here we take trial function F(z)=1-αz2 in (27). It is useful to set ρ=1. The values of the critical temperature Tc for d=4,5, by minimizing the functional (27) with trial function F(z) and using (28) (in case ρ=1) are given by the following.
For d=4, δ=1: Tc=0.0844, for δ=3/4, Tc=0.1692.
For d=5, δ=1: Tc=0.01676, for δ=3/4, Tc=0.2503, for δ=5/4, Tc=0.0954.
These values are in good agreements with the numerical values [50] and also coincide to the analytical values given in [45].
5. Calculating the Critical Exponent
We begin from (22) by writing it near the critical point (CP) and in limit T→Tc. The first step is rewriting the solution in a perturbative scheme with respect to the perturbation parameter ϵ=〈O+〉2, as it was described by Kanno [51]. Near the CP, the solution of the field ϕ is written as
(29)ϕ(z)=ϕ0+ϵχ(z),
where ϕ0=κTc(1-zη). Using (24) and (29) in first order with respect to the O(ϵ) we obtain the following ordinary differential equation:
(30)χ′′(z)-η-1zχ′(z)=E(z),
where
(31)E(z)=Cδr+2z-δ-2F2(z)ξf(z)ϕ02(ϕ0)′2(1-δ),F(z)|z→0≈1.
Since ϕ(z)=(A(Tc)η+1/Tη)(1-zη),A=(d-1)/4π writing the solution for ϕ in z=0, we have
(32)ϕ(0)=ϕ0(0)+ϵχ(0).
The general solution for (30) is given by
(33)χ(z)=z(c1J1(x)+c2Y1(x))+zπη2-2δTc4-2δκ4-2δf(x),
with
(34)f(x)=αJ1(x)∫Y1(x)(zη-1)2z2η-1/2+δ-2δηzd-zdz+βY1(x)∫J1(x)(zη-1)2z2η-1/2+δ-2δηzd-zdz.
Here, x=2(1-η)z, {Jn(x),Yn(x)} are Bessel functions of first and second kinds.
Finding the value of χ(0) from (33), and solving it for 〈O+〉=ϵ, we obtain (we take Tc=1)
(35)〈O+〉∝Td-δ-η/2[1-Tη]1/2.
When the power δ decreases, the value of the 〈O+〉 increases. Thus we conclude that the effect of the power δ in d=4 model is in the direction of the increase of 〈O+〉.
However, in the five dimensions, the analysis is a little bit different. As we observe, in d=5, when the power δ increases, the value of the 〈O+〉 increases. Thus, we can say that the effect of the power δ in d=5 model is in the direction of the increase of 〈O+〉.
6. Calculating the Low Temperature DC Conductivity
In this section, we compute the low-temperature DC conductivity. We concentrate on the general space time dimension d, and we will try to calculate the conductivity σ as a function of the rescaled frequency as follows:
(36)ω^=ω〈O+〉1/Δ,Δ=d-1.
In this limit, the asymptotic form of the scalar field is
(37)ψ(z)=bd-12zd-1F(z),b≡〈O+〉.
We follow the method in [18]. Assuming that there exists an external magnetic field A(r,t)=A(r)e-iωt. Note that here, the applied Maxwell field is linear. It is not related to the nonlinear Maxwell's field in the bulk action. In fact, the nonlinearity of the Maxwell field now is stored in the background metric. As we know, the field equations have some terms which involve the exponent δ. This parameter denotes the nonlinearity, which is hidden in the structure of the background metric and through it. Moreover, it diffuses to the dynamics of the scalar field and the Abelian gauge field U(1), so to compute the conductivity the applied external magnetic field is linear and satisfies the usual linear Maxwell field Fμν;μ=0, which is nothing but the linear wave equation.
The linear wave equation for this field reduces to
(38)-d2Adr-2+V(r)A=ω2A.
Equation (38) is written in the Schrodinger's form and in terms of the new tortoise coordinate r-=-(1/r+)Σn=0∞(zn(d-1)+1/(n(d-1)+1)). The horizon r=r+ is located at z=1 or r-→-∞. Here the potential is V=2f(r)ψ2(r). The ingoing waves in horizon behave as a boundary condition (BC) for solving this wave equation (38) read as
(39)A~e-iωr-~(1-z)-iω/(d-1)r+.
The electromagnetic wave equation (38), in the coordinate z, reads
(40)Azz+(2z+f′2f)Az+(-2r+2ψ(z)2z4+ω2r+2z4f)A=0.
By replacing (37), we obtain
(41)Azz+p1(z)Az+Q1(z)A=0,p1(z)=(2z+f′2f),Q1(z)=-r+2b2(d-1)z2(d-1)F2(z)z4+ω2r+2z4f.
To keep the boundary conditions, we put
(42)A=(1-z)-iω/(d-1)r+e-iωz/(d-1)r+Θ(z).
Substituting this ansatz in (41), we obtain
(43)Θ′′+P(z)Θ′+Q(z)Θ=0,(44)P(z)=r+(1-z)(d-1)p1(z)+2iωzr+(1-z)(d-1),(45)Q(z)=((1-z)2iz(1-z)ωr+(d-1)p1(z)+r+2(d-1)2(1-z)2Q1(z)+ω(i(d-1)r+-ωz2))×(r+2(d-1)2(1-z)2)-1.
By imposing the regularity condition on wave function Θ at the black hole horizon z=1, we obtain the following auxiliary boundary condition:
(46)0=Θ′(1)limz→1(12r+(1-z)(d-1)p1(z)+iωz)+Θ(1)limz→1((d-1)2(i(d-1)r+-ωz2)iz(1-z)ωr+(d-1)p1(z)Θ(1)limz→1+r+2(d-1)2(1-z)2Q1(z)Θ(1)limz→1+ω(i(d-1)r+-ωz2)).
Explicitly, by computing the limits, we have
(47)i(-r++iω+r+d)ωΘ(1)=0.
One possibility is Θ(1)=0. Another ω=ir+(d-1),but the last case from (39) leads to the
(48)A~e-iωr-~,(1-z)
which has no meaning as the ingoing wave toward the horizon, so we impose Θ(1)=0. Since we are working in the low temperature limit, we take the limit b→∞, so we rescale the coordinate z by z→z/b. Also we must put one suitable trial form for F(z), for the case of Δ=d-1>3/2,d>2 we put F(z/b)→F(0)=1. We rescale (43), so we obtain
(49)Θ′′+bP(zb)Θ′+b2Q(zb)Θ=0,′=d/d(z/b).
Finally, we obtain
(50)Θ′′+zr+2(1+2ir+ω^d-1)Θ′+3iω^r+(d-1)Θ=0.
The general solution for (50) reads
(51)Θ(z)=ze-z2(1/2+iω^r+)/r+2×(c+M(μ,ν,12(1+2iω^r+)z2r+2)+c-U(μ,ν,(1+2iω^r+)z2r+2))(52)μ=122d-2+4iω^r+d-7iω^r+(1+2iω^r+)(d-1),ν=32.
The DC conductivity in the low temperature limit is defined by
(53)σ(ω^)=iω^Θ′(-ω^2)Θ(-ω^2).
We need the asymptotic limit of (51). Indeed, we guess that limz→∞Θ(z)≈Θ(-ω^2) where here the ω^ is the quasinormal modes and locates on the real axis, so we guess ℜ𝔢[σ(ω^)]=0 except at the poles of ℑ𝔪[σ(ω^)] where as ℜ𝔢[σ(ω^)]≈δ(ω^). First, by imposing Θ(1)=0 we have
(54)Θ(z)=c-ze-z2(1/2+iω^r+)/r+2U(μ,ν,(1+2iω^r+)/r+2)×[(μ,ν,(1+2iω^r+)z2r+2)M(μ,ν,(1+2iω^r+)z2r+2)U(μ,ν,1+2iω^r+r+2)-U(μ,ν,(1+2iω^r+)z2r+2)×M(μ,ν,1+2iω^r+r+2)(μ,ν,(1+2iω^r+)z2r+2)].
Here M(μ,ν,x), U(μ,ν,x) denote the Kummer functions [52].
The quasinormal modes are the solutions of the following equation:
(55)U(μ,ν,1+2iω^r+r+2)=M(μ,ν,1+2iω^r+r+2),
which has no closed, analytical solution and can be solved just numerically.
7. Calculating the Zero Temperature Conductivity
In this section, we compute analytically the conductivity σ(ω^) in the zero temperature T=0 limit. This case corresponds to the limiting case b→∞. We begin by writing (38) in the coordinate z′=z/b(56)Az′z′+12z′[2-(d+1)(bz′)d-11-(bz′)d-1]×Az′+r+2(bz′)4[b2z′2ω^21-(bz′)d-1+r+2(bz′)4-(b2z′)2(d-1)F(bz′)2b2z′2ω^21-(bz′)d-1]A=0.
By taking the limit b→∞ we have
(57)Az′z′+d+12z′Az′+r+2[-ω^2(bz′)d+1-b4(d-2)z′2(d-3)]A=0.
In (57) we set limb→∞F(bz′)=F(0)=1. The closed form of the exact solution for (57) depends on the value of the d. Below we list the solutions for special dimensions d=2, d=3 as follows:
(58)A(z′)=Cz′-1/4N-1+16r+2/2(2iω^r+b3/2z′),d=2(59)A(z′)=H(0,β,δ,γ,(z′2+1)/(z′2-1))z′×(c1+c2∫dz′zH(0,β,δ,γ,(z′2+1)/(z′2-1))),d=3.
Here
(60)β=-144b8r+2+b4+4r+2ω^2b4,δ=-2r+2(b8-ω^2)b4,γ=-144r+2ω^2-b4+4b8r+2b4.
Here Nν(x) is the second kind of the Bessel function and H(0,β,δ,γ,x) denotes the Heun double confluent function [53], by expending (58), (59) in series in the form
(61)A(z′)=A0+A2(bz′)2-⋯
The conductivity can be computed via the following simple formula:
(62)σ(ω^)=2iω^A2A0+iω^2.
So by computing the asymptotic series, the conductivity will be determined by the analytical expression.
It is appropriate here to present the explicit form of the conductivity for one case. We choose the case given by (58). In the zero temperature limit, we treat the horizon size very tiny, so in application we take 1+16r+2≈1. Also, it is better we define x=2ω^r+/b3/2z′, which by the definition of the z′ reads as x=2ωr+/bz. When b→∞, then x≪1, so, we can expand the Bessel function in terms of the infinitesimal argument x, and hence rewrite (58) as the following form:
(63)A(z)=Cbz4N-ν(ix),ν=12,x=2ωr+bz.
Remembering the following identities
(64)N-ν(ix)=i-ν-1I-ν(x)cosh(νx)-iν-1Iνsinh(νx).
Also, we know that
(65)I1/2(x)=2πxcosh(x),I-(1/2)(x)=2πxsinh(x),
so we have
(66)N-(1/2)(ix)=2iπx[icosh(x)cosh(x/2)-sinh(x)sinh(x/2)].
Also, in limit x≪1,we have
(67)icosh(x)cosh(x/2)-sinh(x)sinh(x/2)=2ix-2+76ix-14x2+59360ix3+O(x4).
So for the electromagnetic field we have
(68)A(z)≈Cb4iπr+ω^[izω^r+-2+73iω^r+z-ω^2r+2z].
Finally, we obtain the conductivity using (62) by the following simple formula:
(69)σ(ω^)=12iω^+iω^2,
where as we guess ℜ𝔢[σ(ω^)]=0. Figure 1 shows ℑ𝔪[σ(ω^)] as a function of the ω^ for zero temperature case and in d=2. Equation (69) gives the expression for DC conductivity in the zero temperature limit in AdS2/CFT1 model. Indeed, it is comparable with the numerical results of the previous papers about one dimensional holographic superconductors [54].
Variation of the log(ℑ𝔪[σ(ω^)]) for zero temperature limit in d=2.
8. Conclusion
In this paper, we investigated the analytical properties of a holographic superconductor with power Maxwell's field. We studied the problem in the probe limit. We observed that it is possible to find the critical temperature Tc and the condensation 〈O+〉 and the conductivity σ(ω) via Sturm-Liouville variational approach. We concluded that in d=4, when the power δ decreases, the value of the 〈O+〉 increases. Thus we can say that the effect of the power δ in d=4 model is in the direction of the increase of 〈O+〉. In d=5 when the power δ increases, the value of the 〈O+〉 increases. Thus, we can say that the effect of the power δ in d=5 model is in the direction of the increase of 〈O+〉. Further, we analytically deduced the low temperature and the zero temperature DC conductivity σ as a function of the ω^=ω/〈O+〉1/Δ. Our work helps to give a better understanding of some unfamiliar effects of the holographic superconductors both in different space dimensions and nonlinearity in Maxwell's strength field.
MaldacenaJ. M.The large N limit of superconformal field theories and supergravity19982231252GubserS. S.KlebanovI. R.PolyakovA. M.Gauge theory correlators from non-critical string theory19984281-21051142-s2.0-0347811692WittenE.Anti-de sitter space and holography19982253291HerzogC. P.Lectures on holographic superfluidity and superconductivity2009423434300110.1088/1751-8113/42/34/343001HartnollS. A.Lectures on holographic methods for condensed matter physics2009262222400210.1088/0264-9381/26/22/224002PolicastroG.SonD. T.StarinetsA. O.Shear viscosity of strongly coupled N=4 supersymmetric yang-mills plasma2001878408160110.1103/PhysRevLett.87.081601KovtunP.SonD. T.StarinetsA. O.Holography and hydrodynamics: diffusion on stretched horizons2003710167517012-s2.0-23044505243HartnollS. A.HerzogC. P.HorowitzG. T.Building a holographic superconductor20081013403160110.1103/PhysRevLett.101.031601BreitenlohnerP.FreedmanD. Z.Stability in gauged extended supergravity198214422492812-s2.0-33750157541GubserS. S.Breaking an abelian gauge symmetry near a black hole horizon2008781506503410.1103/PhysRevD.78.065034GeX. H.WangB.WuS. F.YangG. H.Analytical study on holographic superconductors in external magnetic field2010100810810.1007/JHEP08(2010)108LiH. F.CaiR. G.ZhangH. Q.Analytical studies on holographic superconductors in Gauss-Bonnet gravity2011110402810.1007/JHEP04(2011)028CaiR. G.LiH. F.ZhangH. Q.Analytical studies on holographic insulator/superconductor phase transitions20118312712600710.1103/PhysRevD.83.126007ChenC. M.WuM. F.An analytic analysis of phase transitions in holographic superconductors201112638739510.1143/PTP.126.387GeX. H.LengH. Q.Analytical calculation on critical magnetic field in holographic superconductors with backreaction20121281211122810.1143/PTP.128.1211SetareM. R.MomeniD.MyrzakulovR.RazaM.Holographic superconductors in the AdS black hole with a magnetic charge20128604500510.1088/0031-8949/86/04/045005MomeniD.SetareM. R.MajdN.Holographic superconductors in a model of non-relativistic gravity20110511810.1007/JHEP05(2011)118SiopsisG.TherrienJ.MusiriS.Holographic superconductors near the Breitenlohner-Freedman bound29808500710.1088/0264-9381/29/8/085007ZengH. B.GaoX.JiangY.ZongH. S.Analytical computation of critical exponents in several holographic superconductors201110.1007/JHEP05(2011)002GangopadhyayS.RoychowdhuryD.Analytic study of properties of holographic p-wave superconductors20120810410.1007/JHEP08(2012)104NakanoE.WenW. Y.Critical magnetic field in a holographic superconductor2008784304600410.1103/PhysRevD.78.046004AlbashT.JohnsonC. V.A holographic superconductor in an external magnetic field2008200812110.1088/1126-6708/2008/09/121DomenechO.MontullM.PomarolA.SalvioA.SilvaP. J.Emergent gauge fields in holographic superconductors2010201003310.1007/JHEP08(2010)033MontullM.PujolasO.SalvioA.SilvaP. J.Flux periodicities and quantum hair on holographic superconductors201110718181601MontullM.PujolasO.SalvioA.SilvaP. J.Magnetic response in the holographic insulator/superconductor transition20121204135Salvio AHolographic superfluids and superconductors in dilaton-gravity20121209134BobevN.KunduA.PilchK.WarnerN. P.Minimal holographic superconductors from maximal supergravity20112510.1007/JHEP03(2012)064BobevN.HalmagyiN.PilchK.WarnerN. P.Supergravity instabilities of non-supersymmetric quantum critical points20104610.1088/0264-9381/27/23/235013BobevN.KunduA.PilchK.WarnerN. P.Minimal holographic superconductors from maximal supergravity20122510.1007/JHEP03(2012)064FischbacherT.PilchK.WarnerN. P.New supersymmetric and stable, non-supersymmetric phases in supergravity and holographic field theoryhttp://arxiv.org/abs/1010.4910BobevN.KunduA.PilchK.WarnerN. P.Minimal holographic superconductors from maximal supergravity201110.1007/JHEP03(2012)064RoychowdhuryD.AdS/CFT superconductors with power maxwell electrodynamics: reminiscent of the meissner effect201371831089109410.1016/j.physletb.2012.11.019RoychowdhuryD.Effect of external magnetic field on holographic superconductors in presence of nonlinear corrections201286101210600910.1103/PhysRevD.86.106009BanerjeeR.GangopadhyayS.RoychowdhuryD.LalaA.Holographic s-wave condensate with non-linear electrodynamics: a nontrivial boundary value problemhttp://arxiv.org/abs/1208.5902GangopadhyayS.RoychowdhuryD.Analytic study of gauss-bonnet holographic superconductors in born-Infeld electrodynamics20120515610.1007/JHEP05(2012)156GangopadhyayS.RoychowdhuryD.Analytic study of properties of holographic superconductors in born-infeld electrodynamics20120500210.1007/JHEP05(2012)002WuJ. P.CaoY.KuangX. M.LiW. J.The 3+1 holographic superconductor with Weyl corrections2011697215315810.1016/j.physletb.2011.01.045MomeniD.SetareM. R.A note on holographic superconductors with Weyl corrections20112638288910.1142/S0217732311037169MomeniD.SetareM. R.MyrzakulovR.Condensation of the scalar field with stuckelberg and weyl corrections in the background of a planar AdS-schwarzschild black hole201227125012810.1142/S0217751X1250128XMomeniD.MajdN.MyrzakulovR.p-wave holographic superconductors with Weyl corrections20129766100110.1209/0295-5075/97/61001NojiriS.OdintsovS. D.Anomaly-induced effective actions in even dimensions and reliability of s-wave approximation19994631576299041410.1016/S0370-2693(99)00946-6NojiriS.OdintsovS. D.Can quantum-corrected btz black hole anti-evaporate?19981333269598060310.1142/S0217732398002862HaywardS. A.Unified first law of black-hole dynamics and relativistic thermodynamics19981510314710.1088/0264-9381/15/10/017NojiriS.OdintsovS. D.Quantum (in)stability of 2D charged dilaton black holes and 3D rotating black holes19995941104400310.1103/PhysRevD.59.044003JingJ.PanQ.ChenS.Holographic superconductors with Power-Maxwell field2011201104510.1007/JHEP11(2011)045MazharimousaviS. H.HalilsoyM.Black hole solutions in f(R) gravity coupled with non-linear Yang-Mills field20118406403210.1103/PhysRevD.84.064032BornM.InfeldL.Foundations of the new field theory193414485242545110.1098/rspa.1934.0059ObukhovY. N.Conformal invariance and space-time torsion1982901-2131610.1016/0375-9601(82)90037-8PanQ.JingJ.WangB.Holographic superconductor models with the Maxwell field strength corrections20118410.1103/PhysRevD.84.126020HorowitzG. T.RobertsM. M.Holographic superconductors with various condensates20087812812600810.1103/PhysRevD.78.126008KannoS.A note on gauss-bonnet holographic superconductors2011281212700110.1088/0264-9381/28/12/127001AbramowitzM.StegunI.1972New York, NY, USADoverRonveauxA.1995Oxford University PressRenJ.One-dimensional holographic superconductor from AdS3/CFT2 correspondence2010101105510.1007/JHEP11(2010)055