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

The European Physical Journal C
August 2018, 78:654 | Cite as
Analytical and numerical study of backreacting one-dimensional holographic superconductors in the presence of Born–Infeld electrodynamics
AuthorsAuthors and affiliations
Mahya MohammadiAhmad SheykhiMahdi Kord Zangeneh
Open Access
Regular Article - Theoretical Physics
First Online: 17 August 2018
Received: 07 May 2018
Accepted: 04 August 2018
49 Downloads
Abstract
We analytically as well as numerically study the effects of Born–Infeld nonlinear electrodynamics on the properties of \((1+1)\)-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 \(\beta ={1}/{2}\).
1 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 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_{\mu \nu }\) on the boundary corresponds to the bulk metric \(g_{\mu \nu }\) [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 dual to asymptotic BTZ black holes have been explored widely (see e.g. [19, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59]). In order to construct the \((1+1)\)-dimensional holographic superconductors one should employ the AdS\(_{3}\)/CFT\(_{2}\) correspondence. In [49], the \((1+1)\)-dimensional holographic superconductors were explored in the probe limit and its distinctive features in both normal and superconducting phases were investigated. Employing the variational method of the Sturm–Liouville eigenvalue problem, the one-dimensional holographic superconductors have been analytically studied in [50, 51, 52]. It is also interesting to study the \((1+1)\)-dimensional holographic superconductor away from the probe limit by considering the backreaction. In [53], the effects of the backreaction have been studied for s-wave linearly charged one-dimensional holographic superconductors.
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.
2 Basic field equations
The action of three-dimensional AdS gravity in the presence of a gauge and a scalar field is given by
$$\begin{aligned} S= & {} \frac{1}{2\kappa ^{2}}\int \mathrm{d}^{3}x\sqrt{-g}\left( R+\frac{2}{l^{2}} \right) \nonumber \\&+\int \mathrm{d}^{3}x\sqrt{-g} [ \mathcal {L}(\mathcal {F})-|\nabla \psi -iqA\psi |^{2}-m^{2}|\psi |^{2}],\nonumber \\ \end{aligned}$$
(1)
where m and q are the mass and the charge of the scalar field, \(\kappa ^{2}={8\pi 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), \(\mathcal {F}=F_{\mu \nu }F^{\mu \nu }\) where \(F_{\mu \nu }=\nabla _{[\mu }A_{\nu ]}\) is the electrodynamics field tensor and \(A_{\mu }\) is the vector potential. \( \mathcal {L}(\mathcal {F})\) stands for the Lagrangian density of the BI nonlinear electrodynamics defined as
$$\begin{aligned} \mathcal {L}(\mathcal {F})=\frac{1}{b}\left( 1-\sqrt{1+\frac{b\mathcal {F}}{2}} \right) , \end{aligned}$$
(2)
where b is the nonlinear parameter. When \(b\rightarrow 0\), \(\mathcal {L}\) reduces to \(-F_{\mu \nu }F^{\mu \nu }/{4}\), which is the standard Maxwell Lagrangian [5]. Variation of the above action with respect to the scalar field \(\psi \), the gauge field \(A_{\mu }\) and the metric \(g_{\mu \nu }\) yields the following equations of motion:
$$\begin{aligned} 0= & {} \left( \nabla _{\mu }-i{q}A_{\mu }\right) \left( \nabla ^{\mu }-i{q} A^{\mu }\right) \psi -m^{2}\psi , \end{aligned}$$
(3)
$$\begin{aligned} 0= & {} \nabla ^{\mu }\left( 4\mathcal {L}_{\mathcal {F}}F_{\mu \nu }\right) \nonumber \\&-i{q}\left[ -\psi ^{*}(\nabla _{\nu }-i{q}A_{\nu })\psi +\psi (\nabla _{\nu }+i{q}A_{\nu })\psi ^{*}\right] \,, \end{aligned}$$
(4)
$$\begin{aligned} 0= & {} \frac{1}{2\kappa ^{2}}\left[ R_{\mu \nu }-g_{\mu \nu }\left( \frac{R}{2} +\frac{1}{l^{2}}\right) \right] +2F_{ac}F_{b}{}^{c}\mathcal {L}_{\mathcal {F}} \nonumber \\&-\frac{g_{\mu \nu }}{2}\left[ \mathcal {L}(\mathcal {F})-m^{2}|\psi |^{2}-{ |\nabla \psi -i{q}A\psi |^{2}}\right] \nonumber \\&-\frac{1}{2}\left[ (\nabla _{\mu }\psi -i{q}A_{\mu }\psi )(\nabla _{\nu }\psi ^{*}+i{q}A_{\nu }\psi ^{*})+\mu \leftrightarrow \nu \right] ,\nonumber \\ \end{aligned}$$
(5)
where \(\mathcal {L}_{\mathcal {F}}=\partial \mathcal {L}/\partial {\mathcal {F}}\).
Since we would like to consider the effect of the backreaction on the holographic superconductor, we take a metric ansatz as follows [53]:
$$\begin{aligned} {\mathrm{d}s}^{2}=-f(r)\mathrm{e}^{-\chi (r)}{\mathrm{d}t}^{2}+\frac{{\mathrm{d}r}^{2}}{f(r)}+\frac{r^{2}}{l^{2}} {\mathrm{d}x}^{2}. \end{aligned}$$
(6)
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],
$$\begin{aligned} T=\frac{\kappa _\mathrm{{sg}}}{2\pi }=\frac{1}{2\pi }\sqrt{-\frac{1}{2}(\nabla _{\mu }\chi _{\nu })(\nabla ^{\mu }\chi ^{\nu })}, \end{aligned}$$
where \(\kappa _\mathrm{{sg}}\) is the surface gravity and \(\chi =\partial /\partial t\) is the null Killing vector of the horizon. Taking \(\chi ^{\nu }=(-1,0,0)\), we have \(\chi _{\nu }=(f(r_{+})\mathrm{e}^{-\chi (r_{+})},0,0)\) and hence on the horizon where \(f(r_{+})=0\), we find \((\nabla _{\mu }\chi _{\nu })(\nabla ^{\mu }\chi ^{\nu })=-\frac{1}{2}\left[ f^{\prime }(r_{+})\right] ^{2}\mathrm{e}^{-\chi (r_{+})}\). Thus, the temperature is obtained:
$$\begin{aligned} T=\frac{\mathrm{e}^{-\chi (r_{+})/2}f^{^{\prime }}(r_{+})}{4\pi }. \end{aligned}$$
(7)
We also choose the scalar and the gauge fields as [5]
$$\begin{aligned} A_{\mu }=(\phi (r),0,0), \,\quad \psi =\psi (r). \end{aligned}$$
(8)
Substituting (6) and (8) into the field Eqs. (3)–(5), we arrive at
$$\begin{aligned} 0= & {} \psi ^{\prime \prime }+\psi ^{\prime }\left[ \frac{1}{r}+\frac{ f^{\prime }}{f}-\frac{\chi ^{\prime }}{2}\right] +\psi \left[ \frac{ q^{2}\phi ^{2}{\mathrm {e}}^{\chi }}{f^{2}}-\frac{m^{2}}{f}\right] , \end{aligned}$$
(9)
$$\begin{aligned} 0= & {} \phi ^{\prime \prime }+\phi ^{\prime }\left[ -\frac{b\mathrm{e}^{\chi }\phi ^{\prime 2}}{r}+\frac{\chi ^{\prime }}{2}+\frac{1}{r}\right] \nonumber \\&- \frac{2q^{2}\psi ^{2}\phi }{f}\left[ 1-b\mathrm{e}^{\chi }\phi ^{\prime 2}\right] ^{3/2}, \end{aligned}$$
(10)
$$\begin{aligned} 0= & {} f^{\prime }-\frac{2r}{l^{2}} +2\kappa ^{2}r \bigg [ \frac{q^{2}\mathrm{e}^{\chi }\psi ^{2}\phi ^{2}}{f}+f\psi ^{\prime 2}+m^{2}\psi ^{2} \nonumber \\&-\frac{1}{b}+\frac{1}{b\sqrt{1-b\mathrm{e}^{\chi }\phi ^{\prime 2}}} \bigg ] , \end{aligned}$$
(11)
$$\begin{aligned} 0= & {} \chi ^{\prime }+4\kappa ^{2}r\left[ \frac{q^{2}\phi ^{2}\psi ^{2}{ \mathrm {e}}^{\chi }}{f^{2}}+\psi ^{\prime 2}\right] , \end{aligned}$$
(12)
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\rightarrow 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
$$\begin{aligned} q\rightarrow q/a,\phi \rightarrow a\phi , \psi \rightarrow a\psi , \nonumber \\ \kappa \rightarrow \kappa /a, b\rightarrow b/a^{2}, \end{aligned}$$
(13)
$$\begin{aligned} l\rightarrow al, r\rightarrow ar, q\rightarrow q/a, \nonumber \\ m\rightarrow m/a, b\rightarrow a^{2}b. \end{aligned}$$
(14)
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
$$\begin{aligned} \chi \rightarrow 0,\ \ \ f(r)\sim r^{2}, \nonumber \\ \phi \sim \rho +\mu \ln (r),\ \ \ \psi \sim \psi _{-}r^{-\Delta _{-}}+\psi _{+}r^{-\Delta _{+}}, \end{aligned}$$
(15)
where \(\mu \) and \(\rho \) are the chemical potential and the charge density of the field theory on the boundary and \(\Delta _{\pm }=1\pm \sqrt{1+m^{2}}\), which implies \(m^{2}\ge -1\). Actually, \(\chi \) could be a constant near the boundary but by using the symmetry of the field equation, \(\mathrm {e}^{\chi }\rightarrow a^{2}\mathrm {e}^{\chi },\) \(\phi \rightarrow \phi /a\), it could be set to zero there. We consider \(\psi _{-}\) and \(\psi _{+}\) as the source and the dual of the order parameter, respectively. One may think that either \(\psi _{+}\) or \(\psi _{-}\) can be dual to the expectation value of the condensation operator (or order operator) \(\left\langle O\right\rangle \), while the other is dual to its source. However, in general this is not the case. While \(\psi _{+} \) can always be dual to \(\left\langle O\right\rangle \), \(\psi _{-}\) can be dual to it only if \(-1\le m^{2}<0\) where \(\Delta _{-}\) is positive. We give \(\psi _{-}\) the role of source and \(\psi _{+}\) the role of the expectation value of the order parameter \( \left\langle O_{+}\right\rangle \) in this work. Since we seek for study the effects of b and \(\kappa \) 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 \(\Delta _{+}=2\), \(\Delta _{-}=0\) and thus
$$\begin{aligned} \psi \sim \psi _{-}+\frac{\psi _{+}}{r^{2}}, \end{aligned}$$
(16)
near the boundary. We set \(\psi _{-}=0\) at the boundary and consider \(\psi _{+}\) as the dual of order parameter \(\left\langle O_{+}\right\rangle \). In Appendix A, we present the calculation for \( \left\langle O_{+}\right\rangle \) and show that it is proportional to \(\psi _{+}\) and the other parameters do not appear in its formula explicitly. It is remarkable that the asymptotic solution for \(\psi \) 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 \( \psi \) 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 \(\langle O_{+}\rangle \) 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.
3 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_\mathrm{{c}}\) and the chemical potential \(\mu \) 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\) \((\in \left[ 0,1\right] )\). With this new coordinate, the field equations (9)–(12) could be rewritten as
$$\begin{aligned} 0= & {} \psi ^{\prime \prime }+\psi ^{\prime }\left[ \frac{f^{\prime }}{f}- \frac{\chi ^{\prime }}{2}+\frac{1}{z}\right] +\psi \left[ \frac{ r_{+}^{2}\mathrm{e}^{\chi }\phi ^{2}}{z^{4}f^{2}}-\frac{m^{2}r_{+}^{2}}{z^{4}f}\right] , \nonumber \\\end{aligned}$$
(17)
$$\begin{aligned} 0= & {} \phi ^{\prime \prime }+\phi ^{\prime }\left[ \frac{bz^{3}\mathrm{e}^{\chi }\phi ^{\prime 2}}{r_{+}^{2}}+\frac{\chi ^{\prime }}{2}+\frac{1}{z}\right] -\frac{ 2r_{+}^{2}\psi ^{2}\phi }{z^{4}f}\Upsilon ^{\frac{3}{2}}, \end{aligned}$$
(18)
$$\begin{aligned} 0= & {} f^{\prime }+\frac{2r_{+}^{2}}{z^{3}}+\frac{2r_{+}^{2}\kappa ^{2}}{z^{3}} \nonumber \\&\times \left[ \frac{1}{b} (1-\Upsilon ^{-\frac{1}{2}}) -\frac{ z^{4}f\psi ^{\prime 2}}{r_{+}^{2}}-\frac{\mathrm{e}^{\chi }\psi ^{2}\phi ^{2}}{f} -m^{2}\psi ^{2}\right] , \nonumber \\ \end{aligned}$$
(19)
$$\begin{aligned} 0= & {} \chi ^{\prime }-4\kappa ^{2}\left[ \frac{r_{+}^{2}\mathrm{e}^{\chi }\psi ^{2}\phi ^{2}}{z^{3}f^{2}}+z\psi ^{\prime 2}\right] , \end{aligned}$$
(20)
where \(\Upsilon =1-bz^{4}\mathrm{e}^{\chi }\phi ^{\prime 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
$$\begin{aligned} \epsilon \equiv \langle O_{i} \rangle , \end{aligned}$$
where \(i=+\) or −. We focus on solutions for small values of condensation parameter \(\epsilon \), therefore we can expand the model functions as
$$\begin{aligned} \psi\approx & {} \epsilon \psi _{1}+\epsilon ^{3}\psi _{3}+\epsilon ^{5}\psi _{5}+\cdots , \\ \phi\approx & {} \phi _{0}+\epsilon ^{2}\phi _{2}+\epsilon ^{4}\phi _{4}+\cdots , \\ f\approx & {} f_{0}+\epsilon ^{2}f_{2}+\epsilon ^{4}f_{4}+\cdots , \\ \chi\approx & {} \epsilon ^{2}\chi _{2}+\epsilon ^{4}\chi _{4}+\cdots , \end{aligned}$$
where \(\epsilon \ll 1\) near the critical temperature. Moreover, by considering \(\delta \mu _{2}>0\), the chemical potential can be expressed as
$$\begin{aligned} \mu =\mu _{0}+\epsilon ^{2}\delta \mu _{2}+ \cdots \rightarrow \epsilon \thickapprox \Bigg (\frac{\mu -\mu _{0}}{\delta \mu _{2}}\Bigg )^{1/2}. \end{aligned}$$
During a phase transition, \(\mu _\mathrm{{c}}=\mu _{0}\), thus the order parameter vanishes. Meanwhile, the critical exponent \(\beta =\frac{1}{2}\) is in good agreement with the mean field theory result.
At zeroth order of \(\epsilon \), the gauge field equation of motion (18) reduces to
$$\begin{aligned} \phi ^{\prime \prime }+\frac{\phi ^{\prime }}{z}+\frac{bz^{3}\phi ^{\prime 3} }{r_{+}^{2}}=0, \end{aligned}$$
(21)
which could be rewritten as a first order Bernoulli differential equation by taking \(\phi ^{\prime }\) as a new function [67]. Therefore, one obtains
$$\begin{aligned} \phi ^{\prime }=\frac{\lambda r_{+}}{z\sqrt{b\lambda ^{2}z^{2}+1}} \end{aligned}$$
(22)
for small values of b where we define \(\lambda =\mu /r_{+}\) and fix the integration constants by looking at the behavior of \(\phi \) near the boundary given in (15). Integrating (22) and using the fact that \(\phi (z=1)=0\),1 we can obtain
$$\begin{aligned} \phi _{0}(z)= & {} \int _{1}^{z}\frac{\lambda r_{+}}{z}\left( 1-\frac{1}{2} b\lambda ^{2}z^{2}\right) \,\mathrm{d}z \nonumber \\= & {} \lambda r_{+}\log (z)-\frac{1}{4}b\lambda ^{3}r_{+} ( z^{2}-1). \end{aligned}$$
(23)
When \(b=0\) the above equation reduces to the one of [50]. Note that at the zeroth order with respect to \(\epsilon \), \(\psi _{0}=\chi _{0}=0\). Substituting (22) in (19), the equation for f at zeroth order with respect to \(\epsilon \) has the following form:
$$\begin{aligned} f_{0}(z)= & {} r_{+}^{2}g(z), \nonumber \\ g(z)= & {} \frac{1}{z^{2}}-1+\frac{1}{8}b\kappa ^{2}\lambda ^{4}\left( 1-z^{2}\right) +\kappa ^{2}\lambda ^{2}\log (z). \end{aligned}$$
(24)
The asymptotic behavior of the scalar field \(\psi \) was given in (15). Near the boundary, we define the function F(z) so that
$$\begin{aligned} \psi (z)=\frac{\left\langle O_{i}\right\rangle z^{\triangle _{i}}}{\sqrt{2} r_{+}^{\triangle _{i}}}F(z). \end{aligned}$$
(25)
Inserting Eq. (25) in Eq. (17) yields
$$\begin{aligned}&F^{\prime \prime }(z)+F^{\prime }(z)\left[ \frac{g^{\prime }(z)}{g(z)}+ \frac{2\Delta }{z}+\frac{1}{z}\right] \nonumber \\&\quad +F(z)\left[ \frac{\Delta g^{\prime }(z)}{zg(z)}-\frac{m^{2}}{z^{4}g(z)}+ \frac{\Delta ^{2}}{z^{2}}\right] \nonumber \\&\quad -\frac{F(z)}{2z^{4}g(z)^{2}}[ \lambda ^{2}\log (z) ( b\lambda ^{2}r_{+} ( z^{2}-1) -2\log (z)) ] =0.\nonumber \\ \end{aligned}$$
(26)
We can rewrite this equation in the Sturm–Liouville form as
$$\begin{aligned} \left[ T(z)F^{\prime }(z)\right] ^{\prime }+P(z)T(z)F(z)+\lambda ^{2}Q(z)T(z)F(z)=0, \end{aligned}$$
(27)
where the functions T, P, Q are defined as
$$\begin{aligned} T(z)=z^{2\Delta +1}\left[ \frac{1}{z^{2}}-1+\frac{1}{8}b\kappa ^{2}\lambda ^{4}\left( 1-z^{2}\right) +\kappa ^{2}\lambda ^{2}\log (z)\right] , \end{aligned}$$
(28)
$$\begin{aligned} P(z)=\frac{\Delta }{z}\left( \frac{g^{\prime }(z)}{g(z)}+\frac{\Delta }{z} \right) -\frac{m^{2}}{z^{4}g(z)}, \end{aligned}$$
(29)
$$\begin{aligned} Q(z)=-\frac{\log (z)\left( b\lambda ^{2}r_{+}\left( z^{2}-1\right) -2\log (z)\right) }{2z^{4}g(z)^{2}}. \end{aligned}$$
(30)
We can consider the trial function \(F(z)=1-\alpha z^{2}\), which satisfies the required boundary conditions \(F(0)=1\) and \(F^{^{\prime }}(0)=0\). Then the eigenvalue problem could be solved for (27) by minimizing the expression
$$\begin{aligned} \lambda ^{2}=\frac{\int _{0}^{1}T ( F^{\prime 2}-PF^{2}) \mathrm{d}z}{ \int _{0}^{1}TQF^{2}\mathrm{d}z}, \end{aligned}$$
(31)
with respect to \(\alpha \). For the backreaction parameter, we could use the iteration method and define [68]
$$\begin{aligned} \kappa _{n}=n\Delta \kappa ,\ \ \ n=0,1,2,\ldots , \end{aligned}$$
(32)
where \(\Delta \kappa =\kappa _{n+1}-\kappa _{n}\). Here, we take \(\Delta \kappa =0.05\). Since we are interested in finding the effects of nonlinearity on backreaction up to the order \(\kappa ^{2}\), we have
$$\begin{aligned} \kappa ^{2}\lambda ^{2}={\kappa _{n}}^{2}\lambda ^{2}={\kappa _{n}} ^{2}(\lambda ^{2}|_{\kappa _{n-1}})+O[(\Delta \kappa )^{4}], \end{aligned}$$
(33)
where we take \(\kappa _{-1}=0\) and \(\lambda ^{2}|_{\kappa _{-1}}=0\). We shall also retain the linear terms with respect to the nonlinearity parameter b, and therefore
$$\begin{aligned} b\lambda ^{2}=b\left( \lambda ^{2}|_{b=0}\right) +\mathcal {O}(b^{2}). \end{aligned}$$
(34)
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 \(\epsilon \), \(\chi \) is zero.)
$$\begin{aligned} T_\mathrm{{c}}=\frac{f^{\prime }\left( r_{+c}\right) }{4\pi }. \end{aligned}$$
(35)
Using Eqs. (11) and (23), we obtain
$$\begin{aligned} f^{\prime }\left( r_{+c}\right) =2r_{+c}+\frac{2\kappa ^{2}r_{+c}}{b}\left[ 1-\frac{1}{\sqrt{1-b\phi ^{\prime }\left( r_{+c}\right) {}^{2}}}\right] , \end{aligned}$$
(36)
and thus
Table 1
Analytical results of \({T_\mathrm{{c}}}/{\mu }\) for different values of backreaction and nonlinearity parameters.
\(b=0\)
\(b=0.04\)
\(b=0.08\)
Analytical
Numerical
Analytical
Numerical
Analytical
Numerical
\(\kappa ^{2}=0\)
0.0429
0.0460
0.0360
0.0410
0.0275
0.0362
\(\kappa ^{2}=0.05\)
0.0399
0.0369
0.0337
0.0326
0.0260
0.0286
\(\kappa ^{2}=0.1\)
0.0381
0.0295
0.0311
0.0260
0.0218
0.0227
\(\kappa ^{2}=0.15\)
0.0352
0.0236
0.0280
0.0207
0.0174
0.0180
\(\kappa ^{2}=0.2\)
0.0313
0.0189
0.0242
0.0165
0.0136
0.0143
\(\kappa ^{2}=0.25\)
0.0264
0.0151
0.0195
0.0131
0.0089
0.0114
Open image in new window
Fig. 1
The behavior of the condensation parameter as a function of the temperature for different values of backreaction
Open image in new window
Fig. 2
The behavior of the condensation parameter as a function of the temperature for different values of the nonlinearity parameter b
$$\begin{aligned} T_\mathrm{{c}}= & {} \frac{1}{4\pi }(\frac{\mu }{\lambda })[2-\kappa _{n}^{2}(\lambda ^{2}|_{\kappa _{n-1}}) \nonumber \\&-\frac{3}{4}b\kappa _{n}^{2}(\lambda ^{4}|_{\kappa _{n-1},b=0})+b\kappa _{n}^{2}(\lambda ^{4}|_{\kappa _{n-1},b=0})]. \end{aligned}$$
(37)
As an example, if \(b=\kappa ^{2}=0\) we have
$$\begin{aligned} \lambda ^{2}=\dfrac{\frac{2}{3}\alpha ^{2}-\frac{4}{3}\alpha +1}{-\frac{ 251\alpha ^{2}}{864}+\frac{9\alpha }{16}+\frac{\alpha ^{2}\zeta (3)}{4}- \frac{\alpha \zeta (3)}{2}+\frac{\zeta (3)}{4}-\frac{1}{4}}. \end{aligned}$$
Inserting \(\alpha =0.759\), \(\lambda _\mathrm{{min}}^{2}=13.76\) and \(T_\mathrm{{c}}=0.429\mu \). The latter result perfectly agrees with [50].
The values of \(T_\mathrm{{c}}/\mu \) for different backreaction and nonlinearity parameters are listed in Table 1. As shown, the effect of increasing the backreaction parameter \(\kappa \) for a fixed value of nonlinearity parameter b follows the same trend as raising b for a fixed value of \( \kappa \). In both cases, the critical temperature \(T_\mathrm{{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.
4 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
$$\begin{aligned} f(z)= & {} f_{1}\left( 1-z\right) +f_{2}\left( 1-z\right) {}^{2}+\cdots , \end{aligned}$$
(38)
$$\begin{aligned} \phi (z)= & {} \phi _{1}\left( 1-z\right) +\phi _{2}\left( 1-z\right) {}^{2}+\cdots , \end{aligned}$$
(39)
$$\begin{aligned} \psi (z)= & {} \psi _{0}+\psi _{1}\left( 1-z\right) +\psi _{2}\left( 1-z\right) {}^{2}+\cdots , \end{aligned}$$
(40)
$$\begin{aligned} \chi (z)= & {} \chi _{0}+\chi _{1}\left( 1-z\right) +\chi _{2}\left( 1-z\right) {}^{2}+\cdots . \end{aligned}$$
(41)
Note that \(\phi =0\) at horizon, otherwise it will be ill-defined there. In our procedure, we find all coefficients in terms of \(\phi _{1}\), \(\psi _{0}\) and \(\chi _{0}\) by using the equations of motion. Varying them at the horizon, we try to get \(\psi _{-}=\chi =0\) at the boundary. So, the values of \(\psi _{+}\) and \(\mu \) are obtained. In addition, we will set \(r_{+}=1\) by virtue of the equations of motion’s symmetry,
$$\begin{aligned} r\rightarrow ar,f\rightarrow a^{2}f,\phi \rightarrow a\phi . \end{aligned}$$
Performing a numerical solution, we can find the values of \(T_\mathrm{{c}}/\mu \) 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 \( \kappa \), 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\rightarrow 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.
5 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.
Footnotes
1.
It is necessary so that the norm of the gauge potential is finite at the horizon.
Notes
Acknowledgements
We thank the referee for constructive comments which helped us to improve our paper significantly. We also thank Shiraz University Research Council. The work of AS has been supported financially by Research Institute for Astronomy and Astrophysics of Maragha (RIAAM), Iran. MKZ would like to thank Shahid Chamran University of Ahvaz, Iran.
Appendix A: Holographic calculation of \(\left\langle O_{+}\right\rangle \)
Here, we present the calculation for the expectation value of the order parameter \(\left\langle O_{+} \right\rangle \). In this work, we have considered \(m^{2}=0\) and thus the behavior of the scalar field \(\psi \) near the boundary is
$$\begin{aligned} \psi \sim \psi _{-}+\frac{\psi _{+}}{r^{2}}. \end{aligned}$$
(A1)
\(\psi _{-}\) plays the role of source. From gauge/gravity duality, we know that the expectation value of the order parameter \(\left\langle O_{+} \right\rangle \) for this case is given by
$$\begin{aligned} \left\langle O_{+} \right\rangle =\left. \frac{\partial \mathfrak {L} }{\partial \left( \partial _{r}\psi \right) }\right| _{r\rightarrow \infty }, \end{aligned}$$
(A2)
where \(\mathfrak {L}\) is the Lagrangian density which could be obtained from Eq. (1) via the formula \(S=\int d^{3}x \mathfrak {L}\). Employing Eq. (A2) and the method of holographic renormalization [69, 70], one can obtain \(\left\langle O_{+} \right\rangle \).
The term \(\partial _{r}\psi \) appears in \(\mathfrak {L}\) as
$$\begin{aligned} \mathfrak {L}=-\mathrm{e}^{-\chi /2} r f \left( \partial _{r}\psi \right) ^2+\cdots . \end{aligned}$$
(A3)
Note that we have set l to unity by virtue of the symmetry (13). From (A2) and (A3), we have
$$\begin{aligned} \left\langle O_{+} \right\rangle =-2 \lim _{r \rightarrow \infty } \mathrm{e}^{-\chi /2} r f \partial _{r}\psi . \end{aligned}$$
(A4)
Using (15) and (16), one could easily show that \(\left\langle O_{+} \right\rangle =4\psi _{+}\). This result shows that the expectation value of the order parameter is proportional to \(\psi _{+}\) and the other parameters do not appear in its formula explicitly. It is worth to mention that the divergence terms which may arise when one calculates holographic quantities (via Eq. (A2) or other similar ones) do not affect the value of it and can easily be eliminated by using holographic renormalization approach. In the latter method, the divergence is canceled by adding appropriate counterterms to the action [69, 70].
References
1.
J. Bardeen, L.N. Cooper, J.R. Schrieer, Phys. Rev. 108, 1175 (1957)ADSMathSciNetCrossRefGoogle Scholar
2.
J.M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998). arXiv:hep-th/9711200v3 ADSMathSciNetCrossRefGoogle Scholar
3.
S.S. Gubser, I.R. Klebanov, A.M. Polyakov, Phys. Lett. B 428, 105 (1998). arXiv:hep-th/9802109 ADSMathSciNetCrossRefGoogle Scholar
4.
E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998). arXiv:hep-th/9802150 ADSMathSciNetCrossRefGoogle Scholar
5.
S.A. Hartnoll, C.P. Herzog, G.T. Horowitz, Phys. Rev. Lett. 101, 031601 (2008). arXiv:0803.3295 ADSCrossRefGoogle Scholar
6.
G.T. Horowitz, M.M. Roberts, Phys. Rev. D 78, 126008 (2008)ADSCrossRefGoogle Scholar
7.
J. Ren, JHEP. 1011, 055 (2010). arXiv:1008.3904 ADSCrossRefGoogle Scholar
8.
S.A. Hartnoll, Class. Quantum Gravity 26, 224002 (2009). arXiv:0903.3246 ADSCrossRefGoogle Scholar
9.
C.P. Herzog, J. Phys. A 42, 343001 (2009). arXiv:0904.1975 MathSciNetCrossRefGoogle Scholar
10.
G.T. Horowitz, Lect. Notes Phys. 828, 313 (2011). arXiv:1002.1722 ADSCrossRefGoogle Scholar
11.
S.S. Gubser, C.P. Herzog, S.S. Pufu, T. Tesileanu, Phys. Rev. Lett. 103, 141601 (2009). arXiv:0907.3510 ADSMathSciNetCrossRefGoogle Scholar
12.
S.A. Hartnoll, C.P. Herzog, G.T. Horowitz, JHEP 0812, 015 (2008). arXiv:0810.1563 ADSCrossRefGoogle Scholar
13.
J. Jing, S. Chen, Phys. Lett. B 686, 68 (2010). arXiv:1001.4227 ADSCrossRefGoogle Scholar
14.
R.G. Cai, L. Li, L.-F. Li, R.-Q. Yang, Sci China Phys. Mech. Astron. 58, 060401 (2015). arXiv:1502.00437 Google Scholar
15.
X.H. Ge, B. Wang, S.F. Wu, G.H. Yang, JHEP 1008, 108 (2010). arXiv:1002.4901 ADSCrossRefGoogle Scholar
16.
X.H. Ge, S.F. Tu, B. Wang, JHEP 09, 088 (2012). arXiv:1209.4272 ADSCrossRefGoogle Scholar
17.
X.M. Kuang, E. Papantonopoulos, G. Siopsis, B. Wang, Phys. Rev. D 88, 086008 (2013). arXiv:1303.2575 ADSCrossRefGoogle Scholar
18.
Q. Pan, J. Jing, B. Wang, JHEP 11, 088 (2011). arXiv:1105.6153 ADSCrossRefGoogle Scholar
19.
M. Kord Zangeneh, Y.C. Ong, B. Wang, Phys. Lett. B 771, 235 (2017). arXiv:1704.00557 ADSCrossRefGoogle Scholar
20.
R.G. Cai, H.F. Li, H.Q. Zhang, Phys. Rev. D 83, 126007 (2011)ADSCrossRefGoogle Scholar
21.
R.G. Cai, Z.Y. Nie, H.Q. Zhang, Phys. Rev. D 82, 066007 (2010)ADSCrossRefGoogle Scholar
22.
R.G. Cai, L. Li, L.F. Li, JHEP 1401, 032 (2014). arXiv:1309.4877 ADSCrossRefGoogle Scholar
23.
W. Yao, J. Jing, JHEP 1305, 101 (2013). arXiv:1306.0064 ADSCrossRefGoogle Scholar
24.
Z. Zhao, Q. Pan, S. Chen, J. Jing, Nucl. Phys. B 871, 98 (2013). arXiv:1212.6693 ADSCrossRefGoogle Scholar
25.
Y. Liu, Y. Gong, B. Wang, JHEP 1602, 116 (2016). arXiv:1505.03603 ADSCrossRefGoogle Scholar
26.
S. Gangopadhyay, D. Roychowdhury, JHEP 05, 002 (2012). arXiv:1201.6520 ADSCrossRefGoogle Scholar
27.
S. Gangopadhyay, D. Roychowdhury, JHEP 05, 156 (2012). arXiv:1204.0673 ADSCrossRefGoogle Scholar
28.
A. Sheykhi, H.R. Salahi, A. Montakhab, JHEP 1604, 058 (2016). arXiv:1603.00075 ADSGoogle Scholar
29.
H.R. Salahi, A. Sheykhi, A. Montakhab, Eur. Phys. J. C 76, 575 (2016). arXiv:1608.05025 ADSCrossRefGoogle Scholar
30.
A. Sheykhi, F. Shaker, Int. J. Mod. Phys. D 26, 1750050 (2017). arXiv:1606.04364 ADSCrossRefGoogle Scholar
31.
A. Sheykhi, A. Ghazanfari, A. Dehyadegari, Eur. Phys. J. C 78, 159 (2018)ADSCrossRefGoogle Scholar
32.
A. Sheykhi, F. Shaker, Can. J. Phys. 94, 1372 (2016). arXiv:1601.05817 ADSCrossRefGoogle Scholar
33.
A. Sheykhi, F. Shaker, Phys Lett. B 754, 281 (2016). arXiv:1601.04035 ADSCrossRefGoogle Scholar
34.
A. Sheykhi, D. Hashemi Asl, A. Dehyadegari, Phys. Lett. B 781, 139 (2018). arXiv:1803.05724 ADSMathSciNetCrossRefGoogle Scholar
35.
A. Sheykhi, A. Ghazanfari, A. Dehyadegari, Eur. Phys. J. C 78, 159 (2018). arXiv:1712.04331 ADSCrossRefGoogle Scholar
36.
M. Kord Zangeneh, S.S. Hashemi, A. Dehyadegari, A. Sheykhi, B. Wang. arXiv:1710.10162
37.
S. I. Kruglov. arXiv:1801.06905
38.
S. Carlip, Class. Quantum Gravity 12, 2853 (1995). arXiv:gr-qc/9506079 ADSMathSciNetCrossRefGoogle Scholar
39.
A. Ashtekar, J. Wisniewski, O. Dreyer, Adv. Theor. Math. Phys. 6, 507 (2002). arXiv:gr-qc/0206024 MathSciNetCrossRefGoogle Scholar
40.
T. Sarkar, G. Sengupta, B. NathTiwari, JHEP 0611, 015 (2006). arXiv:hep-th/0606084 ADSCrossRefGoogle Scholar
41.
E. Witten, Adv. Theor. Math. Phys. 2, 505 (1998). arXiv:hep-th/9803131 MathSciNetCrossRefGoogle Scholar
42.
S. Carlip, Class. Quantum Gravity 22, R85 (2005). arXiv:gr-qc/0503022 ADSCrossRefGoogle Scholar
43.
D. Birmingham, I. Sachs, S.N. Solodukhin, Phys. Rev. Lett. 88, 151301 (2002). arXiv:hep-th/0112055 ADSMathSciNetCrossRefGoogle Scholar
44.
E. Witten. arXiv:0706.3359
45.
G. Panotopoulos, A. Rincon, Phys. Rev. D 97, 085014 (2018). arXiv:1804.04684 ADSCrossRefGoogle Scholar
46.
A. Rincon, G. Panotopoulos, Phys. Rev. D 97, 024027 (2018). arXiv:1801.03248 ADSCrossRefGoogle Scholar
47.
G. Panotopoulos, A. Rincon, Int. J. Mod. Phys. D 27, 1850034 (2018). arXiv:1711.04146 ADSCrossRefGoogle Scholar
48.
A. Rincon, E. Contreras, P. Bargueño, B. Koch, G. Panotopoulos, A. Hernández-Arboleda, Eur. Phys. J. C 77, 494 (2017). arXiv:1704.04845 ADSCrossRefGoogle Scholar
49.
P. Chaturvedi, G. Sengupta, Phys. Rev. D 90, 046002 (2014). arXiv:1310.5128 ADSCrossRefGoogle Scholar
50.
R. Li, Mod. Phys. Lett. A. 27, 1250001 (2012)ADSCrossRefGoogle Scholar
51.
D. Momeni, M. Raza, M.R. Setare, R. Myrzakulov, Int. J. Theor. Phys. 52, 2773 (2013). arXiv:1305.5163 CrossRefGoogle Scholar
52.
Y. Peng, G. Liu, Int. J. Mod. Phys. A 32(26), 1750160 (2017)ADSCrossRefGoogle Scholar
53.
Y. Liu, Q. Pan, B. Wang, Phys. Lett. B 702, 94 (2011). arXiv:1106.4353 ADSCrossRefGoogle Scholar
54.
N. Lashkari, JHEP 1111, 104 (2011). arXiv:1011.3520 ADSMathSciNetCrossRefGoogle Scholar
55.
H.B. Zeng. arXiv:1204.5325
56.
Y. Bu, Phys. Rev. D 86, 106005 (2012). arXiv:1205.1614 ADSCrossRefGoogle Scholar
57.
Y. Peng. arXiv:1604.06990
58.
G. Alkac, S. Chakrabortty, P. Chaturvedi, Phys. Rev. D 96, 086001 (2017). arXiv:1610.08757 ADSCrossRefGoogle Scholar
59.
B. Binaei Ghotbabadi, M. Kord Zangeneh, A. Sheykhi, Eur. Phys. J. C 78, 381 (2018). arXiv:1804.05442 ADSCrossRefGoogle Scholar
60.
M. Born, L. Infeld, Proc. R. Soc. A 144, 425 (1934)ADSCrossRefGoogle Scholar
61.
B. Hoffmann, Phys. Rev. 47, 877 (1935)ADSCrossRefGoogle Scholar
62.
W. Heisenberg, H. Euler, Physics 98, 714 (1936). arXiv:physics/0605038 [physics.hist-ph]CrossRefGoogle Scholar
63.
H.P. de Oliveira, Class. Quantum Gravity 11, 1469 (1994)ADSCrossRefGoogle Scholar
64.
G.W. Gibbons, D.A. Rasheed, Nucl. Phys. B 454, 185 (1995). arXiv:hep-th/9506035 ADSCrossRefGoogle Scholar
65.
A. Sheykhi, A. Kazemi, Phys. Rev. D 90, 044028 (2014)ADSCrossRefGoogle Scholar
66.
D. Ghorai, S. Gangopadhyay, Eur. Phys. J. C 76, 146 (2016). arXiv:1511.02444 ADSCrossRefGoogle Scholar
67.
T.L. Chow, Mathematical Methods for Physicists—A Concise Introduction (Cambridge University Press, Cambridge, 2000)CrossRefzbMATHGoogle Scholar
68.
C. Lai, Q. Pan, J. Jing, Y. Wang, Phys. Lett. B 749, 437 (2015). arXiv:1508.05926 ADSCrossRefGoogle Scholar
69.
M. Bianchi, D.Z. Freedman, K. Skenderis, Nucl. Phys. B 631, 159 (2002). arXiv:hep-th/0112119 ADSCrossRefGoogle Scholar
70.
K. Skenderis, Class. Quantum Gravity 19, 5849 (2002). arXiv:hep-th/0209067 CrossRefGoogle Scholar
Copyright information
© The Author(s) 2018
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Funded by SCOAP3
Authors and Affiliations
Mahya Mohammadi1Ahmad Sheykhi12Email authorMahdi Kord Zangeneh341.Physics Department and Biruni ObservatoryShiraz UniversityShirazIran2.Research Institute for Astronomy and Astrophysics of Maragha (RIAAM)MaraghaIran3.Physics Department, Faculty of ScienceShahid Chamran University of AhvazAhvazIran4.Center for Research on Laser and PlasmaShahid Chamran University of AhvazAhvazIran