Excited states of holographic superconductors with backreaction

In this paper we investigate the model of the anti-de Sitter gravity coupled to a Maxwell field and a free, complex scalar field, and construct a fully back-reacted holographic model of superconductor with excited states. With the fixed charge $q$, there exist a series of excited states of holographic superconductor with the corresponding critical chemical potentials. The condensates as functions of the temperature for the two operators $\mathcal{O}_1$ and $\mathcal{O}_2$ of excited states are also studied. For the optical conductivity in the excited states, we find that there exist the additional peaks in the imaginary and real parts of the conductivity. Moreover, the number of peaks corresponding to $n$-th excited state is equal to $n$.


Introduction
In the past ten years, the anti-de Sitter/conformal field theory (AdS/CFT) correspondence [1,2,3] has been used to study the strongly correlated systems in condensed matter physics, and received a great deal of attention. In the remarkable papers [4,5,6], with the U (1) symmetry breaking in a four-dimensional Schwarzschild-AdS black hole background, the condensate of a scalar field could be interpreted as the holographic realization of superconductor condensate. When replacing the scalar field with other matter fields, one could obtain the holographic models for various kinds of superconductors. For example, by introducing a U (1) gauge field coupled to a symmetric, traceless second-rank tensor field in the bulk spacetime, the holographic realization of d-wave supeconductor was obtained in [7,8,9]. With the breaking of SU (2) gauge symmetry, the holographic model of p-wave superconductor was also studied in [10]. In addition, two alternative holographic models of p-wave superconductor could been realized from the condensates of a two-form field [11] and a complex, massive vector field with U (1) charge [12,13], respectively. Considering the weak-link barrier between two superconductor materials [14], the holographic model of Josephson junction was studied in [15,16,17,18,19,20,21,22,23,24,25,26,27]. A topdown construction of holographic superconductor from superstring theory was discussed in [28], and a similar construction using an M-theory truncation was proposed in [29,30]. For reviews of holographic superconductors, see [31,32,33,34].
Recently, the authors in [35] constructed the novel solutions of holographic superconductor, in which excited states of a scalar field were explored in the probe limit. In contrast to the ground state of holographic superconductors, the excited states could have some nodes along the radial coordinator, in which the value of the scalar field could change its sign. Moreover, there exist a series of excited states of holographic superconductor with the corresponding critical chemical potentials. It is interesting to mention that the conductivity in the excited states of holographic superconductors has an additional pole in imaginary part and a delta function in real part arising at the low temperature inside the gap. Similar behaviours for conductivities were found in the holographic model of s-wave superconductor with mass m 2 = −9/4 [40] and the type II Goldstone bosons [37]. These behaviors indicate the emergence of new bound states, which embodies the existence of excited states. Away from the large charge limit, the backreaction of the matter fields on the bulk metric needs to be considered [6]. In the present paper, we study a holographic s-wave superconductor model with full back-reaction, focusing on the excited states of scalar field.
After numerically solving the coupled equations, condensate of scalar field and the optical conductivity will be studied.
The paper is organized as follows. In Sec. 2, we review the gravity dual model of a holographic superconductor with backreaction. We explore the relations between the critical chemical potential and the charge, and study the characteristics of the condensate and conductivity of the excited states with backreaction in Sec. 3. The conclusions and discussions are given in the last section.

Review of holographic superconductors
In this section, we review the model of a Maxwell field and a charged complex scalar field coupled to Einstein gravity with a negative cosmological constant Λ in the four-dimensional spacetime [6]. The bulk action reads on the gravity will decrease when q increases, and the large q limit (q → ∞) corresponds to the probe limit (non-backreaction) of the matter fields. In this paper beyond the probe limit, we will numerically solve the full set of equations with finite charge q. The equations of the scalar and Maxwell fields can be derived from the action (2.1): and variation of the action with respect to the metric leads to the Einstein equations with the stress-energy tensor of the scalar and Maxwell fields given by For simplicity, we shall work with units in which ℓ = 1. When the charged scalar field ψ = 0, the solution of Einstein equations (2.4) is the well-known Reissner-Nordström-AdS (RN-AdS) black hole. This solution with a planar symmetric horizon can be written as follows (2.6) and the gauge field is given by 4r 2 , and µ 0 is the chemical potential for U (1) charge. Thus, the event horizon of the black hole is located at r = 1 and the boundary of the asymptotical AdS spacetime is at r → 0. The Hawking temperature is given by which can be regarded as the temperature of the holographic superconductors. Thus, the solution of a RN-AdS black hole could exist only for µ 0 ≤ µ ext = √ 12, which describes the extremal RN-AdS black hole with zero temperature.
It is well known that there exists a critical temperature T c , below which the black hole has a scalar hair that breaks the U (1) gauge symmetry spontaneously. As for T > T c , the black hole with a scalar hair degrades into RN-AdS black hole. In order to build a holographic model of holographic superconductors, one could introduce the following ansatz of metric: and matter fields with Subscript n is the principal quantum number of the scalar field, and n = 0 is regarded as the ground state and n ≥ 1 as the excited states.
With the above ansatzs (2.9) and (2.10), the independent equations of motion can be deduced as follows To obtain the numerical solutions for the four functions ψ n , φ, g and χ, one must set suitable boundary conditions imposed at the event horizon r = 1: In addition, at the asymptotic boundary r → ∞, the asymptotic behaviors of the four functions take the following forms where ψ (i) (i = 1, 2) are the corresponding expectation values of the dual scalar operators O i . According to AdS/CFT dictionary, the parameters ρ and µ are the charge density and chemical potential in the dual conformal field theory, respectively. In the four-dimensional spacetime, the values of mass m need to satisfy the Breitenlohner-Freedman (BF) bound (m 2 ≥ −9/4) [38]. To simplify, we will set m 2 = −2 in this paper.

Numerical results
In order to numerically solve the above coupled equations (2.11)-(2.14), we introduce a new radial variable z = 1/r which maps the semi-infinite region [1, ∞) to the finite region [1,0].
Thus the event horizon and the asymptotic boundary could be fixed at z = 1 and z = 0, respectively. We focus on pseudospectral methods based on the Chebyshev expansion, and the number of grid points ranges between 60 to 300 in the integration region 0 ≤ z ≤ 1.
The iterative process has been performed by using the Newton-Raphson method, and the relative error for the numerical solutions in our paper is estimated to be below 10 −5 .
With the fixed charge q, there exists a critical value of the chemical potential µ, below which the RN-AdS black hole solution is stable and the scalar field vanishes. Above this value, the charge static solution becomes unstable, the scalar condensate could form and lead to the spontaneously breaking of U (1) gauge symmetry. Before solving the full dynamic equations of motion including Einstein equations numerically, we would adapt the same method as in [39], and find the critical value of µ with the fixed charge q. It is convenient to   Table 1. For each curve, the charge q decreases with the increasing chemical potential µ. We see that for the fixed amplitude µ, the charge q of excited state is larger than the ground state.
In Table 1, we present the results of the critical chemical potential µ c for the operator O 2 from the ground to fifth excited states with the charges q = 6, 10, 14, respectively. It is obvious that an excited state has a higher critical chemical potential than the ground state, which means the excited state has a higher critical charge density ρ c . Due that the critical

Condensates
According to AdS/CFT duality, the expectation value of the condensate operator O i could be connected to the scalar field ψ (i) : The values of condensates as functions of the temperature for the two operators O 1 and O 2 of excited states with the charge q = 6 are shown in the left and right panels of Fig. 2, respectively. In both plots, the black, red and blue lines correspond to the ground state, first and second states, respectively. In the left panel, the condensate of the ground state for the operator O 1 starts to curve upwards as ones approach low temperatures and appears similar to that of the probe limit. However, the condensate of excited states appears to converge as T → 0, and there exists the constant condensate near zero temperature. Moreover, in contrast to the behavior of condensate for the operator O 1 in probe limit [35], we find the condensate of the first-excited state is also smaller than that of the ground state in the left panel, however, it is interesting that the second-excited and higher excited states have larger condensate than the first excited state, which is different from the probe limit case.
To clarify this issue, we plot an extra orange curve corresponding to the third excited state.
In the right panel, the condensate of the excited states for the operator O 2 appears to converge as the temperature T → 0, and the condensate of each excited state is larger than the ground state, which appears similar to that of the probe limit.

Conductivity
In this subsection we study the optical conductivity of the excited states of holographic superconductors with backreaction. According to the holographic duality, we can calculate electromagnetic fluctuations in the bulk geometry. Considering a time dependence perturbation of δA = A x (r)e −iωt dx, the linearized equation of the Maxwell equation is given as When implementing ingoing wave boundary conditions at the horizon, we could obtain the asymptotic behaviour of the gauge field A x at the boundary: With the AdS/CFT correspondence and Ohm's law, the conductivity can be computed by the following formula In Eq. (3.20), the scalar field ψ 0 corresponds to the ground state, and the conductivity of the ground state has been studied in [6].
In Fig. 3, we present the AC conductivity as a function of the frequency ω at the low temperature T /T c ≃ 0.2 for the operator O 2 from the ground state to second excited state. The real and imaginary parts of conductivity are plotted in the left and right panels, respectively. In both plots, the black, red and blue lines correspond to the ground state, first and second states, respectively. From the imaginary part in the right panel, we could see the optical conductivity of excited states also develops a gap at certain frequency which can be identified as the minimum of the imaginary part of the optical conductivity [40].
Moreover, the gaps of excited states are located inside the gap of the ground state. It is interesting to note that there exists an additional peak in the imaginary part of the conductivity for the first excited state, and two additional peaks in Im[σ] of the second excited state. Similar behaviours also appear in the real part of the conductivity. In the left panel of Fig. 3, there also exist additional peaks in the excited states. Moreover, the number of peaks corresponding to n-th excited state is equal to n. Recent work in [35] shows that in the probe limite case, there exists an additional pole in Im[σ] and a delta function in Re[σ] arising at low temperature inside the gap. In contrast to the probe limit case, we find that due to the introduction of backreaction, the pole and delta function can be broaden into the peaks with finite width.

Discussion
In this paper, we numerically solved the full dynamic equations of motion including Einstein equations, and constructed the holographic model of excited states of holographic superconductors with backreaction. When the temperature drops below the critical temperature, the condensate of holographic superconductor begins to appear, which could be regarded as the ground state solution. Further decreasing the temperature to another critical temperature, a new branch of solution with a node in radial direction begins to develop, which is the first excited state solution. As the temperature continues to decrease to lower values, we could obtain a series of excited states of holographic superconductor with corresponding critical temperatures, which is similar to the probe limit case. We also studied the optical conductivity in the excited states of holographic superconductors with backreaction. In contrast to the conductivity of the excited states in the probe limit case, it is very interesting to note that due to the introduction of backreaction, the pole and delta function in optical conductivity at the probe limit case can be transformed into the peaks with finite values at the backreaction case. Moreover, the number of peaks corresponding to n-th excited state is equal to n.
There are some interesting extensions of our work which we plan to investigate in future projects. First, in order to break the translational symmetry which means the momentum of charge carriers could be dissipated, the authors in [39] introduced a gravitational background lattice to the holographic model. It would be very interesting to study the optical conductivity of the excited states of holographic superconductor with holographic lattices.
The second extension of our study is to explore the properties of holographic superconductors with backreaction using the semi-analytical techniques [41]. Finally, we are planning to investigate the excited states of the p-wave holographic superconductor with backreaction and study the excited vector condensates and optical conductivity in future work.