Abstract

Recently it was revealed that the whole Fermi surface is fully gapped for several families of underdoped cuprates. The existence of the finite energy gap along the -wave nodal lines (nodal gap) contrasts the common understanding of the -wave pairing symmetry, which challenges the present theories for the high- superconductors. Here we propose that the incommensurate diagonal spin-density-wave order can account for the above experimental observation. The Fermi surface and the local density of states are also studied. Our results are in good agreement with many important experiments in high- superconductors.


The energy gap is one of the most important properties in the studies of the high- superconductors. Recently, measurements of the energy gaps by angle-resolved photoemission spectroscopy (ARPES) in lightly doped high- materials revealed the existence of the nonzero energy gap along the diagonal directions of the Brillouin zone (also referred to as the “nodal gap”) [16]. The earliest indication of fully gapped single-particle excitation was reported in [7]. The existence of the nodal gap seems to be generic. It has been observed in several families of cuprates [17]. Moreover, it has been reported that the nodal gap exists in the antiferromagnetic (AF) state [5, 6], the spin glass region [3], and the superconducting and normal states for deeply underdoped region [1, 2, 4, 7]. This result is surprising and contrasts the usual understanding of the -wave superconducting pairing or the conventional pseudogap behavior; both should generate energy nodes along diagonal lines of the Brillouin zone. The opening of the nodal gap in the AF state is also intriguing. It was reported that the commensurate AF order forms at K, well above the temperature that the nodal gap opens, which is only K [6].

The above experimental observations challenge the present theory for high- superconductors. As is known, even in the superconducting state the quasiparticle energy gap is not necessarily tied to the superconducting order parameter. It is rather important to explore the physics behind the gap-like feature. Previously, it was proposed that the spin-density-wave order may account for the nodal gap [8]. Recently many groups have attempted to give other possible theoretical scenarios for the nodal gap. It was proposed that the Coulomb disorder effects [9] or strong disordered magnetism [10] can account for the nodal gap while it was argued that the experimentally observed spectra are sufficiently sharp near the nodal momentum so that the nodal gap is unlikely disorder driven [2, 11]. Very recently, it was proposed that the underdoped high- superconductors are topological superconductors and the topological phase is induced by the AF order. The nodal gap is due to an additional -wave pairing term [11, 12] or -pairing term [13]. However, the existence of the AF order in the superconducting phase needs further experimental evidence. And so far there is no direct experimental observation to support the -wave or pairing component in high- superconductors. Moreover, note that the nodal gap is also observed in the nonsuperconducting state [16]. It is temperature independent when crossing the superconducting transition temperature [2]. This seems to suggest that the nodal gap is not superconductivity driven; thus the topological superconductivity scenario is challenged. Therefore, so far the origin for the nodal gap is still far from obvious and a widely accepted theoretical explanation is still awaited.

Actually, now our understanding of the underdoped high- superconductors is still rather incomplete. There may exist various competing orders that can induce the gap-like feature. One possible candidate order is the spin order. Experimentally the presence of the spin order can be examined by neutron scattering experiments. It was reported that the incommensurate diagonal spin-density-wave order (ID-SDW) with the wave vectors forms in the deeply underdoped region [1419]. As the doping density increases, the spin modulation along the parallel direction with the incommensurate wave vectors and was observed, accompanied by the appearance of the superconductivity. It was revealed that the diagonal spin order could persist in the superconducting state and coexist with the parallel spin order [18]. Very recently, based on the muon spin rotation measurement, the ID-SDW order was also observed in the AF state, with the temperature about 30 K. It is much lower than the AF Neel temperature, while it is close to that of the nodal gap being observed [6]. Therefore, the ID-SDW order and the nodal gap appear in the same region experimentally.

Motivated by the above experimental observations and theoretical attempts, and in order to give a more definitive explanation for the nodal gap, in the present work, we start from a phenomenological model in the presence of an ID-SDW order with the wave vectors to elaborate its effect on the spectral function. Our numerical results show that the nodal gap is reproduced in the presence of the above ID-SDW order. The numerical results are in good agreement with the experiments; thus we propose that the ID-SDW order is the main reason for the nodal gap. We also discuss the origin of the ID-SDW order and propose that it is understandable within the Fermi surface nesting picture. The renormalized Fermi surface in the normal state and the modulated real space local density of states are also calculated and the results are qualitatively consistent with previous experiments, which may be the additional evidences for the presence of the ID-SDW order. We also emphasize that our present theory is fundamentally different from previous one considering the existence of the static commensurate AF order [1113]. Namely, in the present work, the superconducting pairing has always pure symmetry. The nodal gap will appear for even very weak ID-SDW order. It is induced directly by the spin order and would persist in the normal state.

Our starting phenomenological model includes the superconducting term and an ID-SDW order term, which is expressed aswhere the superconducting Hamiltonian is expressed asIn the present work, we assume phenomenologically the SDW ordered periodically, expressed as

The above Hamiltonian can be transformed to the momentum space by taking into account the -wave superconducting pairing, the nearest-neighbor, and next-nearest-neighbor hopping, which is rewritten aswhere and . is expressed aswith being the ID-SDW order magnitude.

In the present work, we consider the four incommensurate scattering wave vectors , consistent with previous experiments [6, 1419]. Hereafter, if not specified otherwise, the parameters are used as , , (corresponding to the doping ), , and . The incommensurability is taken as , which is obtained from the Fermi surface nesting vector, as discussed below. We have checked numerically that our results are not sensitive to the reasonable changes of the chosen parameters.

Before presenting our results, let us elaborate the numerical technique with the incommensurate order. Firstly one may rewrite the Hamiltonian to the matrix form with the basis including the particle operator and (). For commensurate density wave (e.g., ), the wave vectors and have the same reduced wave vectors in the first Brillouin zone. As a result, the Hamiltonian is expressed as a finite order matrix and can be solved strictly. However, for incommensurate density-wave order with irrational real space modulation, the order of the Hamiltonian matrix is indeed infinite and can only be solved approximately. Since generally the incommensurability is usually considered approximately as a finite decimal, then the Hamiltonian can be solved numerically. For the present incommensurability under consideration , the whole Brillouin zone is divided into parts. The Hamiltonian with the superconducting pairing can be written as matrix. Then the retarded Green’s function can be obtained through diagonalizing the Hamiltonian. The quasiparticle spectral function is given from the retarded Green’s function with .

We show in Figures 1(a) and 1(b) the spectral function along the diagonal direction (along the cut indicated in Figure 1(c)) with and without the ID-SDW order, respectively. The dashed lines are the quasiparticle dispersions. The spectral functions as a function of the energy (EDCs) along the diagonal direction are plotted in Figure 1(d). As is seen, the quasiparticle energy decreases as the wave vector moves towards the Fermi surface. An obvious gap exists in the presence of the ID-SDW order, as shown in Figure 1(a). We can also see clearly that the gap closes and the quasiparticle dispersion crosses the Fermi momentum for , as is seen in Figure 1(b). We also checked numerically that the above results are in fact independent of the -wave pairing magnitude and the nodal gap exists when we set (not shown here). Therefore, the above nodal gap should exist both in the superconducting state and the normal state. Our results for the nodal gap are qualitatively consistent with the experiments [17].

The momentum dependence of the energy gap along the Fermi surface is studied in Figure 2. The EDCs with different Fermi surface angle (defined in Figure 2(b)) are plotted in Figure 2(a). We define the energy gaps as the peak positions of EDCs. Then the energy gap as a function of the Fermi angle is shown in Figure 2(c). As is seen, the energy gap is significantly anisotropic. It reaches the maxima value at the Brillouin boundary and decreases when the wave vector moves towards the diagonal direction. It reaches the minimum value at the diagonal direction. The -wave gap magnitude is also plotted in Figure 2(c) for comparison. The observed energy gap and the -wave gap are nearly the same near the antinodal direction. Near the diagonal direction the gap is different from the -wave one; namely, an obvious finite gap exists due to the presence of the ID-SDW order. The above results are qualitatively consistent with the experimental observations in the superconducting state [2]. Note that here the wave vector is important for the momentum dependence of the energy gap. We also note that very recent ARPES experiment on the insulating samples has also revealed that the energy gap is anisotropic; that is, it reaches the maxima value at the Brillouin zone boundary and is minimum at the diagonal direction [5]. This experimental result can also be explained qualitatively based on our model. Note that the origin of the energy gap in deeply underdoped high- materials is indeed complicated. There may exist several candidate competing orders. We propose that the energy gap near the -wave nodal points is still due to the ID-SDW order. The gap near the Brillouin boundary is generated by another order (e.g., the -density wave order [20] or the AF order (at the mean-field level, the AF order can be obtained simply by setting in (5))). The coexistence of two orders in the insulating sample is also supported by very recent experiments [6]. Theoretically, the energy gap produced by the -density-wave order or AF order should be maximum near the hot spots (the crossing points between the normal state Fermi surface and the magnetic Brillouin zone). For hole doped samples, the hot spots are close to the antinodal points. Thus the anisotropic gap in nonsuperconducting materials [5] is understandable. We have also checked numerically (not shown here) that similar anisotropic behavior can be reproduced with the model including both the ID-SDW order and the -density-wave order (or the AF order).

The explanation of the ID-SDW order can be given based on the Fermi surface nesting picture. The normal state Fermi surface is shown in Figure 3. As is seen, the tangent lines of the Fermi surface curve are parallel at the -wave nodal points, revealing the Fermi surface nesting feature. The corresponding nesting vector is marked in Figure 3, for which the incommensurability can be obtained from the Fermi momentum along the diagonal direction. The ID-SDW order comes mainly from such node-to-node excitations. As such, the existence of the nodal gap can be immediately understood. Namely, the vector of the ID-SDW order connects different nodal points of Fermi surface. The electron hopping between these points can occur due to the ID-SDW order. This destroys the state of the quasiparticle near nodal points and an energy gap opens.

We now discuss the differences between our model and previous one in [8]. A commensurate SDW order with the wave vector was taken into account in [8]. For the ID-SDW order considered in our present work, the SDW gap is opened first at the nodal direction, while, for the commensurate SDW order, the SDW gap is opened first near the antinodal direction and the nodal gap appears for a relatively strong SDW order. As a result, the momentum space of the total energy gap (shown in Figure 2) may be different as the SDW wave vector changes to the commensurate one.

The signatures of this ID-SDW order can be probed by studying the normal state Fermi surface. Previously many interesting results for the Fermi surface of underdoped cuprates have been obtained by ARPES experiments. One important observation is that the Fermi surface is gapped near the antinodal direction and leaves an ungapped Fermi arc [21]. This conventional pseudogap behavior is still unsolved and not concerned with in the present work. On the other hand, the electronic structure along diagonal direction is also nontrivial. It was revealed that the spectral weight is low and the quasiparticle peak is broad near the nodal direction for underdoped samples [22, 23]. Another interesting result is the observation of the Fermi pocket in the underdoped samples [24]. It was revealed that the Fermi pocket coexists with the Fermi arc and exists only in the underdoped samples. Interestingly, the experimentally observed Fermi pocket is not symmetrical with respect to the to line. Thus the -density-wave order or AF order may not account for the Fermi pocket. It was also proposed in [24] that the incommensurate diagonal density-wave may potentially explain their results.

The numerical results of the normal state zero energy spectral function () is presented in Figure 4. The normal state Fermi surface can be obtained through the peaks of the spectral function. As is seen in Figure 4(a), the spectral weight near the diagonal direction is quite low, consistent with the experimental observations [22, 23]. This is due to the node-to-node scattering caused by the ID-SDW. When the spectral function is plotted in a logarithmic scale, the weak features of the spectral function are revealed more clearly. As is seen, besides the main Fermi surface, contributed by the normal state energy band , another band with much lower spectral weight can be seen clearly. Then a Fermi pocket forms. Here the Fermi pocket is nonsymmetrical and not centered at (). The above results are qualitatively consistent with the recent experimental observation [24].

The relationship of the SDW and charge-density-wave (CDW) orders has been an important point and attracted intensive attention previously. Experimentally the charge order could be detected through the STM experiments. One prominent feature is the “checkerboard structure” from the energy-dependent local density of states (LDOS). It was first reported to exist in vortex cores of optimally doped materials, with a two-dimensional modulation along Cu–O bond directions [2527]. Later experiments also observed similar modulation in the superconducting samples without the magnetic field [28, 29]. In the meantime, the charge order can also be revealed in more detail through the Fourier transform of the LDOS (FT-LDOS). The nondispersive peaks would be observed in the presence of the charge modulation. The periodicity can be determined through the peak positions in the momentum space. An incommensurate charge modulation with the periodicity of about 4.5a4.7a was reported in the normal and superconducting states [3032].

We now study numerically the real space modulation induced by the ID-SDW order. Diagonalizing the Hamiltonian (1), we can obtain the LDOS numerically. To compare with the STM experiments, we also define the Fourier transformation of the LDOS (FT-LDOS), which is expressed as . The numerical results for the LDOS and FT-LDOS are presented in Figure 5. The LDOS in the normal state with the energy is plotted in Figure 5(a). The checkerboard pattern is revealed clearly. Figure 5(b) is the FT-LDOS spectra. There exist four peaks at the wave vector and (indicated with circles). The above results are robust and qualitatively the same for different energies. The LDOS and FT-LDOS in the superconducting state are plotted in Figures 5(c) and 5(d). As is seen, the intensities decrease due to the existence of the superconducting gap while the main results are qualitatively the same as those of the normal state. Interestingly, although here we consider the ID-SDW order in the starting model, the modulations of LDOS are along the Cu–O bond directions. This is consistent with the STM observations. It is also worthwhile to point out that the relation between the SDW order and charge order is still an open question. If both have the same origins, then a simple relation for the incommensurability of the SDW order and the CDW order should satisfy [33]. Experimentally this relation is consistent with the observations in [33, 34] and [35] samples. However, it was also revealed that this relation is not satisfied in the [36, 37] and samples [38].

In summary, based on a phenomenological model, we elaborate that an ID-SDW order can cause a finite gap along the -wave nodal line. This is in good agreement with recent experimental observations. The origin of the ID-SDW order can be explained through the Fermi surface nesting picture. The normal state Fermi surface and the local density of states are also studied. The results are qualitatively consistent with the experiments.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the NSFC (Grants nos. 11374005 and 11204138), the NCET (Grant no. NCET-12-0626), NSF of Jiangsu Province of China (Grant no. BK2012450), Jiangsu Qingnan Engineering project, and U.S. DOE Office of Basic Energy Sciences.