Generalized superconductors from the coupling of a scalar field to the Einstein tensor and their refractive index in massive gravity

We construct the generalized superconductors from the coupling of a scalar field to the Einstein tensor in the massive gravity and investigate their negative refraction in the probe limit. We observe that the larger graviton mass and Einstein tensor coupling parameters both hinder the formation of the condensation, but the larger graviton mass or smaller coupling parameter makes it easier for the emergence of the Cave of Winds. Furthermore, we see that the larger graviton mass but smaller coupling parameter make the range of frequencies or the range of temperatures larger for which a negative Depine-Lakhtakia index occurs, which indicates that the graviton mass and Einstein tensor have completely different effects on the negative refraction. In addition, we find that the larger graviton mass and coupling parameters both can reduce the dissipation and improve the propagation in the holographic setup. PACS numbers: 11.25.Tq, 04.70.Bw, 74.20.-z ∗ panqiyuan@126.com † jljing@hunnu.edu.cn


I. INTRODUCTION
It is well known that Maldacena conjectured a remarkable connection between a strongly coupled field theories to a gravitational theory in higher spacetime dimensions in 1998 [1], which states that a d-dimensional weakly coupled dual gravitational description in the bulk is equivalent to a (d − 1)-dimensional strongly coupled conformal field theory on the boundary. This is the so-called anti-de Sitter/conformal field theories (AdS/CFT) correspondence [2,3], which provides a novel approach to studying the strongly coupled systems in condensed matter physics in the past ten years, especially the high temperature superconductors [4]. With this holographic duality, it is shown that the holographic superconductor models exhibit many characteristic properties shared by real superconductors [5]; for reviews, see Refs. [6][7][8][9] and references therein. However, in most cases, the studies on the holographic superconductors focus on the superconductors without an impurity.
Turning on a coupling between the gauge field and a new massive gauge field, the authors of Ref. [10] analyzed the impurity effect in a holographic superconductor and observed that the mass gap in the optical conductivity disappears when the coupling is sufficiently large. More recently, Hu et al. constructed the holographic Josephson junction from the massive gravity and found that the increasing graviton mass parameter will change the superconductor phase into a normal metal phase for a fixed temperature (chemical potential), which indicates that the graviton mass parameter and the doping have a similar effect as from superconductivity to a normal metal [11]. Considering a holographic superconductor with a scalar field coupled kinematically to the Einstein tensor, Kuang and Papantonopoulos observed that, as the strength of the coupling increases, the critical temperature below which the scalar field condenses is lowering, the condensation gap decreases faster than the temperature, the width of the condensation gap is not proportional to the size of the condensate and at low temperatures the condensation gap tends to zero for the strong coupling, which suggests that the derivative coupling in the gravity bulk can have a dual interpretation on the boundary corresponding to impurities concentrations in a real material [12]. Other investigations based on the impurity effect on the holographic dual models can be found, for example, in Refs. [13,14].
On the other hand, since Veselago proposed in theory that the refractive index might be negative in some special material [15], there have been many investigations concerning this exotic electromagnetic phenomenonnegative refraction, which indicates that electromagnetic waves propagate in a direction opposite to that of the flow of energy in the so-called "metamaterials" [16][17][18]. Interestingly, it is believed that the superconductors are the ideal candidates for metamaterials due to their low loss, compact structure, extraordinary degree of nonlinearity and tunability, magnetic flux quantization and the Josephson effect, quantum effects in which photons interact with quantized energy levels in the meta-atom, and strong diamagnetism [19,20]. With the help of the AdS/CFT correspondence, the refractive index in the superconductors was investigated by holography in the probe limit [21] and away from the probe limit [22], which shows that in the superconducting phase there is negative refraction at low enough frequencies only in the backreacted case. It is interesting to note that Mahapatra et al. [23] constructed the generalized holographic superconductors via the Stückelberg mechanism [24][25][26] in the four dimensional R-charged black hole and observed that the superconducting phase can exhibit a negative Depine-Lakhtakia (DL) index [27] at low frequencies and below a cut-off value of the charge parameter even in the probe limit. Extending the investigation to the generalized superconductors with Born-Infeld electrodynamics, Cheng et al. found that the system has a negative DL index in the superconducting phase at small frequencies and the greater the Born-Infeld corrections the larger the range of frequencies or the range of temperatures for which the negative refraction occurs [28]. The refractive index in the holographic dual models can also be found, for example, in Refs. [29][30][31][32][33][34][35][36][37][38][39].
In this work, motivated by the recent studies in Refs. [11,12], we will use the AdS/CFT correspondence to investigate systematically the effect of the doping or impurity on the superconductors in strongly coupled condensed matter systems and their negative refraction in the probe limit. In order to obtain the rich phase structure, we construct the generalized superconductors from the coupling of a scalar field to the Einstein tensor in the massive gravity via the Stückelberg mechanism just as in [23]. We will observe that the larger graviton mass and Einstein tensor coupling parameters both hinder the formation of the condensation, but the larger graviton mass or smaller coupling parameter makes it easier for the emergence of the socalled Cave of Winds [40], i.e., the second-order transition occurs before the first-order transition to a new superconducting phase when the temperature decreases [41][42][43]. Furthermore, we assume that the boundary theory is weakly coupled to a dynamical electromagnetic field and one can calculate the refractive index of the system perturbatively. We will find that the larger graviton mass parameter or smaller coupling parameter makes the range of frequencies or the range of temperatures larger for which a negative DL index occurs, but the larger graviton mass and coupling parameters both can reduce the dissipation and improve the propagation in the holographic setup.
The organization of this work is as follows. In Sec. II we will use the AdS/CFT correspondence to construct the generalized superconductors from the coupling of a scalar field to the Einstein tensor in the massive gravity.
In Sec. III we will discuss the effects of the graviton mass and Einstein tensor coupling parameters on the scalar condensation and phase transition in the generalized holographic models. In Sec. IV we will analyze the effects of the graviton mass and Einstein tensor coupling parameters on the negative refraction in the generalized holographic models. We will conclude in the last section of our main results.

II. GENERALIZED SUPERCONDUCTORS FROM THE COUPLING OF A SCALAR FIELD TO THE EINSTEIN TENSOR IN MASSIVE GRAVITY
Since we are interested in the background of a massive gravity theory, we start with the following action for a four-dimensional ghost-free dRGT massive gravity [44,45] where f is the reference metric, m is the graviton mass parameter, c i are constants and U i are symmetric polynomials of the eigenvalues of the matrix K µ ν ≡ √ g µα f αν where the square root in K means ( The equation of motion of this action will be with We consider the following reference metric [44] f µν = diag(0, 0, with c 0 being a positive constant, and a general black hole solution is given by [45] with where m 0 is related to the mass of the black hole, and h ij dx i dx j is the line element for the two-dimensional spherical, flat or hyperbolic space with κ = −1, 0 or 1, respectively. In this work we are only interested in the plane-symmetric black hole, so we will set κ = 0 and the corresponding Hawking temperature is determined by where r h is the radius of the event horizon. For this considered solution, the nonzero components of the Einstein tensor G µν are where the prime denotes a derivative with respect to r. Without loss of generality, we will scale the AdS radius L = 1 in the following calculation.
In the probe limit, we consider a Maxwell field and a charged complex scalar field coupled to the Einstein tensor G µν with the generalized action where the charged scalar field Ψ with the mass m Ψ and the phase α both are real, and the local U (1) gauge symmetry is given by A µ → A µ + ∂ µ λ together with α → α + λ. When the coupling parameter η → 0 and the graviton mass m → 0, our model reduces to the generalized holographic superconductors of the Einstein's general relativity (GR) investigated in [23].
Using the gauge symmetry to fix the phase α = 0, from the action (10) we can obtain equations of motion for the matter fields Taking the ansatz of the matter fields as Ψ = Ψ(r) and A = Φ(r)dt, under the solution (6), the above equations of motion can be rewritten into Using the shooting method [46], we can solve Eqs. (13) and (14) numerically by doing integration from the horizon out to the infinity with the appropriate boundary conditions for Ψ and Φ. At the event horizon r = r h , the regularity leads to the boundary conditions At the infinity r → ∞, the solutions behave like with the characteristic exponent where µ and ρ are interpreted as the chemical potential and charge density in the dual field theory, respectively.
According to the AdS/CFT correspondence, provided ∆ − is larger than the unitarity bound, both Ψ − and Ψ + can be normalizable and they can be used to define operators in the dual field theory Ψ − = O − and Ψ + = O + , respectively. Considering that we focus on the condensate for the operator O + , we will impose boundary condition Ψ − = 0 in this work.

III. SCALAR CONDENSATION AND PHASE TRANSITION
In this section, we will solve the system of coupled differential equations (13) and (14) numerically and then discuss the scalar condensation and phase transition of these generalized superconductors. For concreteness, we will set m 2 Ψ = −2 since the choices of the scalar field mass will not qualitatively modify our results.
Following [23,28], we consider a particular form of G(Ψ) for simplicity, i.e., where ξ is the model parameter. In order to compare with the results given in Refs. [23,28], we will concentrate on the cases for the model parameters ξ = 0.2 and ξ = 0.5 in our discussion. For the convenience of numerics, we will set c 1 = r h and c 2 = −r 2 h /2, and choose the range of graviton mass parameters 0 ≤ m ≤ 1.2, just as in [11]. It is interesting to note that, from the equations of motion (13) and (14), there exist the useful scaling symmetries and the transformation of the relevant quantities where λ is a real positive number. We will use the scaling symmetries to set r h = 1 when performing numerical calculations.
On the other hand, we have to calculate the grand potential Ω = −T S os of the bound state in order to obtain the important information about the phase transition of the system, where S os is the Euclidean on-shell action. Working in the grand canonical ensemble, we have where we have used the integration dtdxdy = V T . Thus, we can express the grand potential in the superconducting phase as In the normal phase, i.e., Ψ = 0, we can easily get the grand potential A. Effect of the graviton mass on the scalar condensation and phase transition In Fig. 1, we plot the condensate and corresponding grand potential as a function of the chemical potential with the fixed coefficient ξ for different values of the graviton mass m in the case of the coupling parameter η = 0.0 by solving the equations of motion (13) and (14) numerically. Obviously, we can see clearly from left two panels that, for a fixed model parameter ξ, the critical chemical potential µ c increases with the increase of the graviton mass parameter, i.e., the critical temperature T c decreases as m increases, which indicates that the larger graviton mass parameters will hinder the formation of the condensation. This agrees well with the findings in the first holographic Josephson junction from the massive gravity introduced in Ref. [11], i.e., the graviton mass parameter and the doping have a similar effect as the phase transition from superconductivity to a normal metal.
What is more interesting is the rich phase structure of this system. In the case of ξ = 0 or small ξ, for example ξ = 0.2, the transition is of the second order and the condensate approaches zero as for all values of m considered here, which is supported by the behavior of the grand potential just as shown in the top-right panel for m = 1.2. Note that this conclusion still holds in the cases of m ≤ 1.0 with ξ = 0.5.
However, the story is obviously different if we consider the case of m = 1.2 in the bottom two panels of Fig.   1. We find that, as the chemical potential is increased, i.e., the temperature is lowered, the condensate O + becomes multivalued, and there is at first a second order phase transition from the normal to a superconducting phase and then a first order phase transition between two superconducting phases, which is the Cave of Winds Before proceeding, we first briefly review the electric permittivity ǫ(ω), effective magnetic permeability µ(ω) and their dependence on the retarded correlators in the linear response theory investigated by Amariti et al. [29]. Expanding the transverse part of the dielectric tensor for the small wave vector k as one can obtain the refractive index for the isotropic media [29] n 2 (ω) = k 2 ω 2 = ǫ(ω)µ(ω).
With the form of ǫ T (ω, k) specified by the transverse part of the retarded current-current correlator G T (ω, k) [47,48], i.e., and the expanded form of the retarded correlator in k the electric permittivity and effective magnetic permeability can be expressed as where C em is the electromagnetic coupling constant which will be set to unity when performing numerical calculations, and G 0 T (ω) and G 2 T (ω) are the expansion coefficients of the retarded correlators [49], respectively.
It should be noted that, when taking the dissipation of the system into account, the electric permittivity ǫ, magnetic permeability µ and refractive index n are generally complex functions of the frequency ω. In Ref.
[50], McCall et al. derived the condition for the negative refraction in the dissipative medium and pointed out that it is not necessary for real parts Re(ǫ) < 0 and Re(µ) < 0 simultaneously. Later, Depine and Lakhtakia found that, for the dissipative medium, the sign of the refractive index coincide with the sign of the DL index [27] n DL = Re[ǫ(ω)]|µ(ω)| + Re[µ(ω)]|ǫ(ω)|, where n DL < 0 indicates that the phase velocity in the medium is opposite to the direction of the energy flow, i.e., the system has a negative index of refraction. As a simple, convenient, and widely adopted criterion, we will use the DL index n DL to check if the medium has positive or negative refraction in this work.
In order to compute the refractive index of our holographic systems, we will use the AdS/CFT correspondence to derive the expression of G 0 T (ω) and G 2 T (ω). Assuming that the perturbed Maxwell field has a form

we get the equation of motion
Taking the same series expansion of A x as in Eq. (26) for G T (ω, k), i.e., from Eq. (29) we arrive at the equations of motion for A x0 and A x2 which can be numerically solved by the ingoing wave boundary conditions near the horizon and the behaviors near the asymptotic AdS boundary Thus, according to the AdS/CFT dictionary, we obtain the retarded correlators which can be used to calculate ǫ(ω) and µ(ω) in Eq. (27) and n DL in Eq. (28).
A. Effect of the graviton mass on the negative refractive index Just as pointed out in Ref. [29], the refraction may be negative when Re(ǫ) and Re(µ) are not simultaneously negative because of the presence of the imaginary parts. It should be noted that Im(µ) < 0 could imply some problem in the ǫ − µ approach [22,23,35], although Markel argued that Im(µ) can in fact be negative for diamagnetic materials [51]. We will leave this point for further study since this delicate issue goes beyond the scope of our paper.
In order to obtain the effect of the graviton mass m on the negative refraction, we plot the behavior of the obvious for the larger graviton mass in both cases, i.e., T = 0.81T c and T = 0.85T c , which indicates that the system with the larger graviton mass has a negative DL index in the superconducting phase even T < 0.81T c or T > 0.85T c . This means that the larger graviton mass also makes the range of temperatures larger for which negative refraction appears. Thus, it is interesting to note that we can use the graviton mass m to accommodate the range of frequencies and the range of temperatures for which negative refraction occurs. Note that |k 2 G 2 T (ω)/G 0 T (ω)| ≪ 1 is the validity condition of Eq. (26). With the help of Eq. (24), we present in Fig. 6 the curves of | T (ω) n 2 |ω 2 as a function of ω/T with the fixed temperatures T = 0.81T c and T = 0.85T c for different values of the graviton mass parameter m. Due to the appearance of a negative imaginary part of the magnetic permeability, the validity condition is not very strictly satisfied, just as shown in Fig. 6, which may be an unfortunate feature in the probe limit [23,35]. However, we find that within the negative refraction frequency range | T (ω) n 2 |ω 2 is always less than one for the fixed temperatures T and all values of m considered here, which means that the validity condition gets marginally satisfied. Thus, our results are reliable. It is also important to calculate the ratio Re(n)/Im(n) and analyze the propagation Re(n) and the dissipation Im(n) of electromagnetic waves in systems. In Fig. 7 . It is shown that, for the permittivity, Re(ǫ) is negative at low frequencies and Im(ǫ) is always positive, but for the permeability, Re(µ) is always positive and Im(µ) is negative, which implies that the refraction may be negative since the real parts of the permittivity Re(ǫ) and of the permeability Re(µ) are not simultaneously negative. It should be noted that we have defined an effective magnetic permeability that is not an observable, although there is a caveat, i.e., the appearance of a negative imaginary part of the magnetic permeability. We will proceed with this point in mind. that the larger Einstein tensor coupling parameter also reduces the range of temperatures for which negative refraction occurs. The behavior of the Einstein tensor is reminiscent of that seen for the graviton mass, so we conclude that the graviton mass and Einstein tensor have completely different effects on the negative refractive index of our holographic systems.
In order to verify the validity of the expansion used in Eq. (26), we plot in Fig. 11 the behavior of   increasing values of η for the fixed ω/T , which indicates that the larger Einstein tensor coupling parameter enhances the magnitude of Re(n)/Im(n) and can be used to improve the propagation in the holographic setup.
Interestingly, the effect of the Einstein tensor on the ratio Re(n)/Im(n) is similar to that of the graviton mass.

V. CONCLUSIONS
Recent studies showed that the graviton mass parameter has the similar effect on the phase transition from superconductivity to a normal metal as that of the doping [11] and the coupling of a scalar field to the Einstein tensor in the gravity bulk can have a dual interpretation on the boundary corresponding to impurities concentrations in a real material [12]. In order to investigate systematically the influence of the doping or impurity on the superconductors in strongly coupled condensed matter systems by holography, we constructed the generalized superconductors from the coupling of a scalar field to the Einstein tensor in the massive gravity and analyzed the effects of the graviton mass and Einstein tensor on the holographic superconductor models and their negative refraction. In the probe limit, we found that the critical chemical potential increases with the increase of the graviton mass or Einstein tensor coupling parameter, which supports the findings in Ref. [11,12] and indicates that the larger graviton mass or Einstein tensor parameters hinder the formation of the condensation. Interestingly, in our model we observed that the second-order transition occurs before the first-order transition to a new superconducting phase when the chemical potential increases (namely the temperature decreases), i.e, the so-called Cave of Winds, and the larger graviton mass but smaller coupling parameter make it easier for the emergence of the Cave of Winds. On the other hand, we investigated the optical properties of our generalized holographic superconductors and found the existence of the negative Depine-Lakhtakia index in the superconducting phase at small frequencies. We noted that the larger graviton mass but smaller Einstein tensor coupling parameter make the range of frequencies or the range of temperatures larger for which a negative DL index occurs, which means that the graviton mass and Einstein tensor have completely different effects on the negative refraction. Considering the diverse effects of the doping or impurity on the refractive index of the real material, for example LiNbO 3 [52,53], our result suggests that the doping or impurity will have various effects on the negative refraction in the high-temperature superconductor systems.
Furthermore, we found that the larger graviton mass and coupling parameters both can reduce the dissipation and improve the propagation in the holographic setup. Our analysis shows that the graviton mass and the coupling of the scalar field to the Einstein tensor both can play important roles in determining the optical properties of the boundary theory.