Analytical and numerical study of backreacting one-dimensional holographic superconductors in the presence of Born–Infeld electrodynamics

We analytically as well as numerically study the effects of Born–Infeld nonlinear electrodynamics on the properties of (1+1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(1+1)$$\end{document}-dimensional s-wave holographic superconductors. We relax the probe limit and further assume the scalar and gauge fields to affect the background spacetime. We thus explore the effects of backreaction on the condensation of the scalar hair. For the analytical method, we employ the Sturm–Liouville eigenvalue problem, and for the numerical method, we employ the shooting method. We show that these methods are powerful enough to analyze the critical temperature and phase transition of the one-dimensional holographic superconductor. We find that increasing the backreaction as well as the nonlinearity makes the condensation harder to form. In addition, this one-dimensional holographic superconductor faces a second order phase transition and the critical exponent has the mean field value β=1/2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta ={1}/{2}$$\end{document}.


Introduction
The best-known theory for describing the mechanism behind superconductivity from a microscopic perspective is the BCS theory proposed by Bardeen, Cooper and Schrieffer. According to BCS theory, the condensation of Cooper pairs into a boson-like state, at low temperature, is responsible for an infinite conductivity in the solid state system [1]. However, when the temperature increases, the Cooper pair decouples and thus the BCS theory is unable to explain the mechanism of superconductivity for high temperature superconductors [1]. The correspondence between gravity in an anti-de Sitter (AdS) spacetime and Conformal Field Theory (CFT) living on the a e-mail: asheykhi@shirazu.ac.ir b e-mail: mkzangeneh@scu.ac.ir boundary of the spacetime provides a powerful tool for calculating correlation functions in a strongly interacting field theory using a dual classical gravity description [2]. According to the AdS/CFT duality proposal an n-dimensional conformal field theory on the boundary is equivalent to gravity theory in (n+1)-dimensional AdS bulk [2][3][4][5][6][7]. The dictionary of AdS/CFT duality implies that each quantity in the bulk has a dual on the boundary. For example, the energy-momentum tensor T μν on the boundary corresponds to the bulk metric g μν [3,4]. Based on this duality, Hartnoll et al. proposed a model for a holographic superconductor in 2008 [5]. Their motivation was to shed light on the problem of high temperature superconductors. According to the theory of holographic superconductors, we need a hairy black hole on the gravity side to describe a superconductor on its boundary. During the past decade, the investigation of the holographic superconductor has got a lot of attention (see e.g. [6][7][8][9][10][11][12][13][14][15][16][17][18][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36][37]).
On the other hand, BTZ (Bandos-Teitelboim-Zanelli) black holes, the well-known solutions of general relativity in (2+1)-dimensional spacetime, provide a simplified model to investigate some conceptual issues in black hole thermodynamics, quantum gravity, string theory, gauge theory and the AdS/CFT correspondence [38][39][40][41][42]. It has been shown that the quasinormal modes in this spacetime coincide with the poles of the correlation function in the dual CFT. This gives quantitative evidence for AdS/CFT [43]. In addition, BTZ black holes play a crucial role for improving our perception of gravitational interaction in low-dimensional spacetimes [44]. These kinds of solutions have been studied from different point of views [45][46][47][48].
Holographic superconductors have also been studied extensively in the presence of nonlinear electrodynamics (see e.g. [24,25,[28][29][30][31][32][33]36,37]). The most famous nonlinear electrodynamics is Born-Infeld electrodynamics. This model was presented for the first time to solve the problem of the divergence of the electrical field at the position of a point particle [60][61][62][63][64]. It was later showed that this model could be reproduced by string theory. In the present work, we would like to extend the investigation of the (1 + 1)-dimensional holographic superconductors by taking into account the nonlinear Born-Infeld (BI) electrodynamics, as our gauge field. Furthermore, we will study the effects of backreaction on our holographic superconductors. We perform our investigation both analytically and numerically and shall compare the result of two methods. Our analytical approach is based on the Sturm-Liouville variational method. In the latter study, we find the relation between the critical temperature and the chemical potential. Moreover, in order to study our holographic superconductors numerically, we use the shooting method. We show that the analytical results are in good agreement with numerical ones, which implies that the Sturm-Liouvile variation method is still powerful enough for studying the (1 + 1)-dimensional holographic superconductor.
The structure of our paper is as follows. In Sect. 2, the basic field equations of one-dimensional holographic superconductors with backreaction in the presence of BI nonlinear electrodynamics is introduced. In Sect. 3, we describe the procedure of an analytical study of the one-dimensional holographic superconductor based on Sturm-Liouvile method and obtain the relation between the critical temperature and the chemical potential. In Sect. 4, the numerical approach to the study of our holographic superconductors will be presented. Finally, we summarize our results in Sect. 5.

Basic field equations
The action of three-dimensional AdS gravity in the presence of a gauge and a scalar field is given by where m and q are the mass and the charge of the scalar field, κ 2 = 8π G 3 and G 3 is the three-dimensional Newton gravitation constant. Also, g, R and l are the metric determinant, Ricci scalar and AdS radius, respectively. In (1), is the electrodynamics field tensor and A μ is the vector potential. L(F) stands for the Lagrangian density of the BI nonlinear electrodynamics defined as where b is the nonlinear parameter. When b → 0, L reduces to −F μν F μν /4, which is the standard Maxwell Lagrangian [5]. Variation of the above action with respect to the scalar field ψ, the gauge field A μ and the metric g μν yields the following equations of motion: where L F = ∂L/∂F. Since we would like to consider the effect of the backreaction on the holographic superconductor, we take a metric ansatz as follows [53]: The Hawking temperature of the three-dimensional black hole on the outer horizon r + (where r + is the horizon radius obtained as the largest root of f (r + ) = 0) may be obtained through the use of the general definition of the surface gravity [65], where κ sg is the surface gravity and χ = ∂/∂t is the null Killing vector of the horizon. Taking χ ν = (−1, 0, 0), we have χ ν = ( f (r + )e −χ(r + ) , 0, 0) and hence on the horizon where f (r + ) = 0, we find (∇ μ χ ν )(∇ μ χ ν ) = − 1 2 f (r + ) 2 e −χ(r + ) . Thus, the temperature is obtained: We also choose the scalar and the gauge fields as [5] A Substituting (6) and (8) into the field Eqs. (3)- (5), we arrive at where the prime denotes the derivative with respect to r . Note that in the presence of the nonlinear BI electrodynamics equations (9) and (12) do not change compared to the linear Maxwell case. In the limiting case where b → 0 the equations of motion (10) and (11) turn into the corresponding equations of the one-dimensional holographic superconductor with a Maxwell field [53]. The field Eqs. (9)-(12) enjoy the symmetries Hereafter, we set q and l equal to unity by virtue of these symmetries. The behavior of our model functions for large r (near the boundary) read where μ and ρ are the chemical potential and the charge density of the field theory on the boundary and ± = 1 ± √ 1 + m 2 , which implies m 2 ≥ −1. Actually, χ could be a constant near the boundary but by using the symmetry of the field equation, e χ → a 2 e χ , φ → φ/a, it could be set to zero there. We consider ψ − and ψ + as the source and the dual of the order parameter, respectively. One may think that either ψ + or ψ − can be dual to the expectation value of the condensation operator (or order operator) O , while the other is dual to its source. However, in general this is not the case. While ψ + can always be dual to O , ψ − can be dual to it only if −1 ≤ m 2 < 0 where − is positive. We give ψ − the role of source and ψ + the role of the expectation value of the order parameter O + in this work. Since we seek for study the effects of b and κ on our holographic superconductors and different values of the scalar field mass do not influence this behavior qualitatively, we consider m 2 = 0 in this work. With this choice, we have + = 2, − = 0 and thus near the boundary. We set ψ − = 0 at the boundary and consider ψ + as the dual of order parameter O + . In Appendix A, we present the calculation for O + and show that it is proportional to ψ + and the other parameters do not appear in its formula explicitly. It is remarkable that the asymptotic solution for ψ given in Eq. (16) does not depend on the type of electrodynamics and thus for the Maxwell case in three dimensions the solution is the same as in Eq. (16). However, the solution depends on the spacetime dimensions. This is due to the fact that the equation for ψ given in (9) is independent of the type of electrodynamics but depends on the spacetime dimensions and the mass parameter m [34,53,66]. The next step is to solve the coupled nonlinear field equations (9)- (12) simultaneously and obtain the behavior of the model functions. Then we could figure out the behavior of different parameters of holographic superconductor, specially the order parameter O + and the critical temperature, by using these functions. In this work, we use both analytical and numerical methods for studying the holographic superconductor. For analytical study, we perform the Sturm-Liouville method, while for numerical study, we use the shooting method.

Sturm-Liouville method
In this section, we employ the Sturm-Liouville eigenvalue problem to investigate analytically the phase transition of the one-dimensional s-wave holographic superconductor in the presence of BI nonlinear electrodynamics. In addition, we calculate the relation between the critical temperature T c and the chemical potential μ near the horizon. Furthermore, we study the effect of backreaction and BI nonlinear electrodynamics on the critical temperature. For future convenience, we define a new variable z = r + /r (∈ [0, 1]). With this new coordinate, the field equations (9)-(12) could be rewritten as where ϒ = 1−bz 4 e χ φ 2 /r 2 + and now the prime indicates the derivative with respect to z. Since in the vicinity of critical temperature the order parameter is small, we can consider it as an expansion parameter where i = + or −. We focus on solutions for small values of condensation parameter , therefore we can expand the model functions as where 1 near the critical temperature. Moreover, by considering δμ 2 > 0, the chemical potential can be expressed as During a phase transition, μ c = μ 0 , thus the order parameter vanishes. Meanwhile, the critical exponent β = 1 2 is in good agreement with the mean field theory result.
At zeroth order of , the gauge field equation of motion (18) which could be rewritten as a first order Bernoulli differential equation by taking φ as a new function [67]. Therefore, one obtains for small values of b where we define λ = μ/r + and fix the integration constants by looking at the behavior of φ near the boundary given in (15). Integrating (22) and using the fact that φ(z = 1) = 0, 1 we can obtain When b = 0 the above equation reduces to the one of [50]. Note that at the zeroth order with respect to , ψ 0 = χ 0 = 0. Substituting (22) in (19), the equation for f at zeroth order with respect to has the following form: The asymptotic behavior of the scalar field ψ was given in (15). Near the boundary, we define the function F(z) so that Inserting Eq. (25) in Eq. (17) yields We can rewrite this equation in the Sturm-Liouville form as where the functions T , P, Q are defined as We can consider the trial function F(z) = 1 − αz 2 , which satisfies the required boundary conditions F(0) = 1 and F (0) = 0. Then the eigenvalue problem could be solved for (27) by minimizing the expression with respect to α. For the backreaction parameter, we could use the iteration method and define [68] κ n = n κ, n = 0, 1, 2, . . . , where κ = κ n+1 − κ n . Here, we take κ = 0.05. Since we are interested in finding the effects of nonlinearity on backreaction up to the order κ 2 , we have where we take κ −1 = 0 and λ 2 | κ −1 = 0. We shall also retain the linear terms with respect to the nonlinearity parameter b, and therefore Then the minimum eigenvalue of Eq. (31) can be obtained. At the critical point, the temperature is defined as (see Eq. (7) and note that at zeroth order with respect to , χ is zero.) Using Eqs. (11) and (23), we obtain and thus As an example, if b = κ 2 = 0 we have Inserting α = 0.759, λ 2 min = 13.76 and T c = 0.429μ. The latter result perfectly agrees with [50].
The values of T c /μ for different backreaction and nonlinearity parameters are listed in Table 1. As shown, the effect of increasing the backreaction parameter κ for a fixed value of nonlinearity parameter b follows the same trend as raising b for a fixed value of κ. In both cases, the critical temperature T c diminishes by the growth of the backreaction or nonlinearity parameters. It shows that the presence of backreaction and Born-Infeld nonlinear electrodynamics makes the scalar hair harder to form. In the next section, we will re-study the problem numerically using the shooting method.

Shooting method
In this section, we will study our holographic superconductor numerically. In order to do this, we use the shooting method [8]. In this method, the boundary values are found by setting appropriate initial conditions. So, for doing this, we need to know the behavior of the equations of motion both at horizon and at the boundary. Using a Taylor expansion at the horizon around z = 1, we get Note that φ = 0 at horizon, otherwise it will be ill-defined there. In our procedure, we find all coefficients in terms of φ 1 , ψ 0 and χ 0 by using the equations of motion. Varying them at the horizon, we try to get ψ − = χ = 0 at the boundary. So, the values of ψ + and μ are obtained. In addition, we will set r + = 1 by virtue of the equations of motion's symmetry, Performing a numerical solution, we can find the values of T c /μ for different backreaction and nonlinearity parameters. In order to compare the latter results with the analytical ones, we list both of them next to each other in Table 1. It is obvious that there is a reasonable agreement between the results of both methods. Moreover, in Table 1, the results of [53] for b = 0 have been recovered for different values of the backreaction parameter. As one could see in this table, increasing the backreaction parameter for a fixed value of b decreases the critical temperature. This means that larger values of the backreaction parameter make the condensation harder to form. Similarly, for a fixed value of κ, increasing the nonlinearity of the electrodynamic model makes scalar hair harder to form because it diminishes the critical temperature.   Figures 1 and 2 give information about the effect of the backreaction and nonlinear electrodynamics on condensation. All curves follow a same trend. As b → 0, we regain the results of the Maxwell case presented in [53]. As the figures show, the condensation gap increases by making the backreaction and nonlinearity parameters larger, while the other one is fixed. So, it can be understood that it is harder to form a superconductor. This is in agreement with the results obtained from the behavior of the critical temperature before.

Summary and discussion
In this work, by using the Sturm-Liouville eigenvalue problem, we analytically investigated the properties of the (1+1)dimensional holographic superconductor developed in a BTZ black hole background in the presence of BI nonlinear electrodynamics. We have relaxed the probe limit and further assumed that the gauge and scalar fields do backreact on the background metric. We determined the critical temperature for different values of backreaction and nonlinear parameters. We have continued our study by using the numerical shooting method and confirmed that the analytical results are in agreement with the numerical approach. We observed that the formation of scalar hair is harder in the presence of BI nonlinear electrodynamics and backreaction and it becomes harder and harder to form by increasing the strength of either the nonlinear or the backreaction parameters.
Finally, it would be of interest to extend this procedure for other nonlinear electrodynamics like the Power-Maxwell version and logarithmic cases and investigate the effects of nonlinear electrodynamics on the critical temperature and condensation operator of one-dimensional holographic superconductors. These issues are now under investigation and the results will appear elsewhere. raz University Research Council. The work of AS has been supported