arXiv:cond-mat/0307660 v2 25 Nov 2003
Collective modes and quasiparticle interference on the loc
al density of states of
cuprate superconductors
C.-T. Chen
1
and N.-C. Yeh
1
1
Department of Physics, California Institute of Technology
, Pasadena, CA 91125, USA
(Dated: May 30, 2005)
The energy, momentum and temperature dependence of the quas
iparticle local density of states
(LDOS) of a two-dimensional
d
x
2
−
y
2
-wave superconductor with random disorder is investigated
using the first-order T-matrix approximation. The results s
uggest that collective modes such as
spin/charge density waves are relevant low-energy excitat
ions of the cuprates that contribute to the
observed LDOS modulations in recent scanning tunneling mic
roscopy studies of Bi
2
Sr
2
CaCu
2
O
x
.
PACS numbers: 74.50.+r, 74.62.Dh, 74.72.-h
One of the most widely debated issues in cuprate su-
perconductivity is the possibility of preformed Cooper
pairs
1,2,3,4
and the origin of the pseudogap phe-
nomenon.
5,6
Recent experiments have demonstrated
that the pseudogap phenomenon is unique to the
hole-doped (p-type) cuprates and is absent above
T
c
in electron-doped (n-type) cuprates.
7,8,9
Furthermore,
in the quasiparticle tunneling spectra of the double-
layer Bi
2
Sr
2
CaCu
2
O
8+
δ
(Bi-2212)
10
and the one-layer
Bi
2
(Sr
2
−
x
La
x
)Cu
2
O
6+
δ
(Bi-2201)
11
systems, it is shown
that the pseudogap can be distinguished from the super-
conducting gap: the former evolves smoothly with in-
creasing temperature whereas the latter vanishes at
T
c
.
These phenomena suggest that the pseudogap may be
associated with a competing order
6,12
that coexists with
the superconducting phase for
T < T
c
and persists above
T
c
until a pseudogap temperature
T
∗
. The competing
quantum ordered phase
13,14
can be manifested in the
form of collective modes such as charge- and spin-density
waves (CDW and SDW) in the superconducting state,
as inferred from neutron scattering experiments in a va-
riety of p-type cuprates.
15,16,17,18,19
However, whether
these collective modes are closely correlated with super-
conductivity remain controversial. Recent scanning tun-
neling spectroscopic studies of the Fourier transformed
(FT) quasiparticle local density of states (LDOS) of Bi-
2212
20,21,22
have stimulated further discussions on the
relevance of collective modes.
23,24,25,26,27,28
While Bogoli-
ubov quasiparticle interference apparently plays an im-
portant role in the observed FT-LDOS in the supercon-
ducting state, certain spectral details of the LDOS can-
not be accounted for unless collective modes are consid-
ered.
23,24,25
In particular, the findings of 4 high-intensity
Bragg peaks remaining above
T
c
in the FT-LDOS map of
Bi-2212
29
cannot be reconciled with quasiparticles being
the sole low-energy excitations. These new developments
motivate us to reexamine the role of collective modes in
cuprate superconductors by considering the energy (
E
),
momentum transfer (
q
) and temperature (
T
) dependence
of the resulting FT-LDOS modulations.
We begin our model construction by noting that sub-
stantial quasiparticle gap inhomogeneities are observed
in the low-temperature tunneling spectroscopy of under-
and optimally doped Bi-2212 single crystals,
30,31
suggest-
ing at least two types of spatially separated regions, one
with sharp quasiparticle coherence peaks at smaller en-
ergies ∆
d
and the other with rounded hump-like features
at larger energies ∆
∗
. On the other hand, low-energy
LDOS (for
E <
0
.
5∆
d
) of Bi-2212 exhibit long-range
spectral homogeneity. We therefore conjecture that dy-
namic SDW or CDW coexist with cuprate superconduc-
tivity and that they are only manifested in the quasipar-
ticle LDOS when pinned by disorder. Thus, regions with
rounded hump features in the quasiparticle spectra are
manifestation of localized charge modulations due to pin-
ning of collective modes by disorder, and the wavevector
of the charge modulation is twice of that for the collinear
SDW order, as proposed in Refs. 23,24,32. In con-
trast, regions with sharp quasiparticle spectral peaks are
representative of generic Bogoliubov quasiparticle spec-
tra with a well-defined
d
-wave pairing order parameter
∆
k
≈
∆
d
cos 2
θ
k
, where ∆
d
is the maximum gap value
and
θ
k
is the angle between the quasiparticle wavevector
k
and the antinode direction. Our model therefore as-
sumes ‘puddles’ of spatially confined ‘pseudogap regions’
with a quasiparticle scattering potential modulated at a
periodicity of 4 lattice constants along the Cu-O bonding
directions, and the spatial modulations can be of either
the ‘checkerboard’ pattern
23,24
or ‘charge nematic’ with
short-range stripes.
12,14
In the limit of weak perturba-
tions, we employ the first-order T-matrix approximation
and consider a (400
×
400) sample area with either 24
randomly distributed point impurities or 24 randomly
distributed puddles of charge modulations that cover ap-
proximately 6% of the sample area. For simplicity, we do
not consider the effect of disorder on either suppressing
the local pairing potential ∆
d
(
r
) or altering the nearest-
neighbor hopping coefficient (
t
) in the band structure of
Bi-2212, although such effects reflect the internal struc-
tures of charge modulations.
25,26
Specifically, the Hamiltonian of the two-dimensional
superconductor is given by
H
=
H
BCS
+
H
imp
, where
H
BCS
denotes the unperturbed BCS Hamiltonian of the
d
-wave superconductor,
H
BCS
=
∑
k
σ
(
ǫ
k
−
)
c
†
k
σ
c
k
σ
+
∑
k
∆
k
[
c
†
k
↑
c
†
−
k
↓
+
c
−
k
↓
c
k
↑
], and
H
imp
is the perturbation
Hamiltonian associated with impurity-induced quasipar-
2
ticle scattering potential.
24,26,27,33
Using the T-matrix
method, the Green’s function
G
associated with
H
is
given by
G
=
G
0
+
G
0
T
G
0
, where
G
0
is the Green’s func-
tion of
H
BCS
and
T
=
H
imp
/
(1
− G
0
H
imp
). The Hartree
perturbation potential for single scattering events in the
diagonal part of
H
and for non-interacting identical
point impurities at locations
r
i
is
V
α
(
q
) =
∑
i
V
s,m
e
i
q
·
r
i
for non-magnetic (
V
s
) and magnetic (
V
m
) impurities
27
,
whereas that for puddles with short stripe-like modula-
tions centering at
r
j
is
33
V
β
(
q
) =
∑
j
V
0
e
i
q
·
r
j
2 sin(
q
y,x
R
j
) sin(
q
x,y
R
j
)
q
y,x
sin(2
q
x,y
)
,
(1)
and that for checkerboard modulations is
V
γ
(
q
) =
∑
j
V
0
e
i
q
·
r
j
[
2 sin(
q
y
R
j
) sin(
q
x
R
j
)
q
y
sin(2
q
x
)
+ (
q
x
↔
q
y
)
]
(2)
Here all lengths are expressed in units of the lattice con-
stant
a
,
R
j
is the averaged radius of the
j
-th puddle,
and
V
0
denotes the magnitude of the scattering poten-
tial by pinned collective modes. For simplicity, we ne-
glect the energy dependence of
V
α,β,γ
and assume that
V
s
,
V
m
and
V
0
are sufficiently small so that no resonance
occurs in the FT-LDOS.
26
For sufficiently large scatter-
ing potentials, full T-matrix calculations become neces-
sary as in Ref. 27. However, large
V
s,m
would result in
strong spectral asymmetry between positive and negative
bias voltages,
27
which differs from experimental observa-
tion.
20,21,22
We also note that the energy dependence of
V
β,γ
reflects the spectral characteristics of the collective
modes and their interaction with quasiparticles and im-
purities. For instance, we expect
V
γ
∼
ζγ
2
for pinned
SDW, where
ζ
is the impurity pinning strength and
γ
is the coupling amplitude of quasiparticles with SDW
fluctuations.
23,24
Empirically for nearly optimally doped
Bi-2212,
R
j
ranges from 5
∼
10.
31
Here we take different
values for
R
j
with a mean value
h
R
j
i
= 10.
Given the Hamiltonian and the scattering potentials
V
α,β,γ
(
q
), we find that for infinite quasiparticle lifetime
and in the first-order T-matrix approximation, the FT
of the LDOS
ρ
(
r
, E
) that involves elastic scattering of
quasiparticles from momentum
k
to
k
+
q
is:
ρ
q
(
ω
) =
−
1
πN
2
lim
δ
→
0
∑
k
V
α,β,γ
(
q
)
×
{
u
k
+
q
u
k
(
u
k
+
q
u
k
∓
v
k
+
q
v
k
)
ℑ
[
1
(
ω
−
E
k
+
iδ
)(
ω
−
E
k
+
q
+
iδ
)
]
+
u
k
+
q
v
k
(
u
k
+
q
v
k
±
v
k
+
q
u
k
)
ℑ
[
1
(
ω
+
E
k
+
iδ
)(
ω
−
E
k
+
q
+
iδ
)
]
+
v
k
+
q
u
k
(
u
k
+
q
v
k
±
v
k
+
q
u
k
)
ℑ
[
1
(
ω
−
E
k
+
iδ
)(
ω
+
E
k
+
q
+
iδ
)
]
−
v
k
+
q
v
k
(
u
k
+
q
u
k
∓
v
k
+
q
v
k
)
ℑ
[
1
(
ω
+
E
k
+
iδ
)(
ω
+
E
k
+
q
+
iδ
)
]
}
.
(3)
Here
N
is the total number of unit cells in the sam-
ple, and
ℑ
[
. . .
] denotes the imaginary part of the quan-
tity within the brackets, which is related to the equal-
energy quasiparticle joint density of states. The up-
per (lower) sign in the coherence factor applies to
spin-independent (spin-dependent) interactions,
u
k
and
v
k
are the Bogoliubov quasiparticle coefficients,
u
2
k
+
v
2
k
= 1,
u
2
k
= [1 + (
ξ
k
/E
k
)]
/
2,
ξ
k
≡
ǫ
k
−
,
ǫ
k
is the tight-binding energy of the normal state of Bi-
2212 according to Norman
et al.
34
,
ǫ
k
=
t
1
(cos
k
x
+
cos
k
y
)
/
2 +
t
2
cos
k
x
cos
k
y
+
t
3
(cos 2
k
x
+ cos 2
k
y
)
/
2 +
t
4
(cos 2
k
x
cos
k
y
+ cos
k
x
cos 2
k
y
)
/
2 +
t
5
cos 2
k
x
cos 2
k
y
,
t
1
−
5
=
−
0
.
5951
,
0
.
1636
,
−
0
.
0519
,
−
0
.
1117
,
0
.
0510 eV,
is the chemical potential, and
E
k
=
√
ξ
2
k
+ ∆
2
k
.
FIG. 1: Calculated energy-dependent Fourier transform (FT
)
maps of quasiparticle LDOS in the first Brillouin zone
with randomly distributed non-magnetic point defects us-
ing Eq. (3) and
V
α
: (a) ∆
d
= 40 meV and (
ω/
∆
d
) =
±
0
.
15
,
±
0
.
45
,
±
0
.
75 (up and down from left to right); (b)
∆
d
= 20 meV and (
ω/
∆
d
) = 0
.
15
,
0
.
45
,
0
.
75 (left to right).
(c) Schematic illustration of the equal-energy contours an
d
representative modulation wavevectors
q
A
,
q
B
and
q
C
, which
correspond to
q
1
,
q
7
and
q
2
in Refs. 20,21.
Using Eq. (3) and
V
α,β,γ
(
q
), we obtain the energy-
dependent FT-LDOS maps in the first Brillouin zone for
non-magnetic point impurities in Fig. 1 with two differ-
ent ∆
d
values and for pinned SDW (with spin-dependent
coherence factor) in Fig. 2, whereas the corresponding
LDOS modulations due to
V
α,β,γ
(
q
) in real space is
shown in Figs. 3(a)-(c). For non-magnetic point impurity
scattering at
T
≪
T
c
, the intensities associated with
q
B
and
q
C
are much stronger than those of
q
A
, as shown in
Fig. 1 and also in Fig. 4(a). The results in Fig. 1 differ
from the STM observation
20,21
that reveals comparable
3
FIG. 2: Energy-dependent FT-LDOS maps with randomly
distributed pinned SDW using Eq. (3) and
V
γ
: (a) ∆
d
= 40
meV and (
ω/
∆
d
) =
±
0
.
15
,
±
0
.
45
,
±
75, up and down from left
to right; (b) ∆
d
= 20 meV and (
ω/
∆
d
) = 0
.
15
,
0
.
45
,
0
.
75, from
left to right. The FT-LDOS does not exhibit discernible dif-
ferences in the spectral characteristics except the total i
nten-
sities if we simply replace
V
γ
by
V
β
and assume non-magnetic
coherence factors in Eq. (3).
intensities associated with
q
A
and
q
B
, and weaker inten-
sities with
q
C
. Interestingly, the intensities of
q
A
and
q
B,C
become reversed if one assumes magnetic point im-
purity scattering, as illustrated in Fig. 4(b). However,
there is no evidence of magnetic scattering in the samples
used in Refs. 20,21. In contrast, the presence of pinned
collective modes, regardless of CDW or SDW, gives rise
to much stronger intensities for
q
A
(by about two orders
of magnitude), as shown in Fig. 2. Thus, the empirical
FT-LDOS maps
20,21
cannot be solely attributed to quasi-
particle scattering by non-magnetic point impurities.
FIG. 3: Real space quasiparticle LDOS for a (400
×
400)
area at
T
= 0 due to scattering by (a)non-magnetic point
impurities, (b) pinned CDW and (c) pinned SDW, for ∆
d
=
40 meV and
ω
= 30 meV.
The relevance of collective modes become indisputable
when we consider the temperature dependence of the FT-
LDOS. As shown in Fig. 5(a), in the limit of
T
→
T
−
c
,
the
q
-values contribute to the FT-LDOS map become sig-
FIG. 4: Evolution of the relative intensities of FT-LDOS wit
h
energy (
ω
) for
q
A
,
q
B
and
q
C
as defined in Fig.1(c) and
V
s
,
V
m
and
V
0
all taken to be unity: quasiparticle scattering by (a)
single non-magnetic point impurity, and (b) single magneti
c
point impurity.
FIG. 5: The FT-LDOS maps at
T
= 0, 0
.
75
T
c
and
T
c
(from
left to right) for (a) point impurities
V
α
(
q
) and (b) pinned
SDW
V
γ
(
q
). We assume ∆
d
(
T
) = ∆
d
(0)[1
−
(
T/T
c
)]
1
/
2
,
∆
d
(0) = 40 meV, tunneling biased voltage = 18 mV,
and
T
c
= 80 K. Besides temperature dependent coher-
ence factors, the thermal smearing of quasiparticle tunnel
-
ing conductance (
dI/dV
) is obtained by using (
dI/dV
)
∝
|
∫
ρ
q
(
E
)(
df/dE
)
|
(
E
−
eV
)
dE
|
, where
f
(
E
) denotes the Fermi
function. (c)
|
q
A
|
-vs.-
V
(biased voltage) dispersion relation
for pinned SDW at
T
= 0 and
T
c
.
nificantly extended and smeared for point-impurity scat-
tering. In contrast, pinned SDW yields strong intensi-
ties in the FT-LDOS map only at
q
A
for
T >
∼
T
c
, as
shown in Fig. 5(b). The overall energy dispersion due to
pinned SDW is weaker than that due to point impuri-
ties, as shown in Fig. 5(c) for
|
q
A
|
-vs.-
V
(biased voltage)
at both
T
= 0 and
T
=
T
c
. In particular, we note that
the dispersion is further reduced at
T
c
. These findings
are consistent with recent experimental observation by
4
Yazdani
et al
.
29
The energy, momentum and temperature dependence
of our calculated FT-LDOS in Figs. 1-5 is supportive
of spatially modulated collective modes being relevant
low energy excitations in cuprates besides quasiparticles
.
In particular, only pinned collective modes can account
for the observation in the FT-LDOS map above
T
c
. Al-
though our simplified model cannot exclude CDW, we
note that pinned CDW would have coupled directly to
the quasiparticle spectra and resulted in stripe-like peri
-
odic local conductance modulation, which has not been
observed in STM studies. On the other hand, various
puzzling phenomena seem reconcilable with the SDW
scenario. For instance, the nano-scale gap variations ob-
served in Bi-2212
31
may be understood by noting that
the LDOS in regions with disorder-pinned SDW contains
information of disorder potential coupled with quasipar-
ticles and SDW, so that the hump-like spectral features
at
±
∆
∗
represent neither the SDW gap nor the super-
conducting gap ∆
d
, and the values of ∆
∗
vary in accor-
dance with the disorder potential. The long-range spatial
homogeneity of quasiparticle spectra in YBa
2
Cu
3
O
7
−
δ
(YBCO)
35
as opposed to the strong spatial inhomogene-
ity in Bi-2212 can also be reconciled in a similar context.
That is, SDW can be much better pinned in extreme two-
dimensional cuprates like Bi-2212 than in more three-
dimensional cuprates such as YBCO. Furthermore, SDW
can be stabilized by magnetic fields,
23,24,36
which natu-
rally account for the checkerboard-like spectral modula-
tions around the vortex cores of Bi-2212.
19,37
Finally, the
smooth evolution of the pseudogap phase with tempera-
ture through
T
c
may contribute to the anomalously large
Nernst effect under a c-axis magnetic field above
T
c
,
38
with spin fluctuations responsible for the excess entropy.
In summary, we employ first-order T-matrix approx-
imation to study modulations in the quasiparticle FT-
LDOS of cuprates as a function of energy, momentum and
temperature. Our results suggest that a full account for
all aspects of experimental observation below
T
c
must in-
clude collective modes as relevant low energy excitations
besides quasiparticles, and that only collective modes can
account for the observed FT-LDOS above
T
c
.
Acknowledgments
We thank Professors Subir Sachdev, Doug Scalapino
and C. S. Ting and Mr. Yuan-Yu Jau for useful dis-
cussions. This research was supported by NSF Grant
#DMR-0103045.
1
P. W. Anderson, Science
235
, 1196 (1987).
2
E. Dagotto, J. Riera, and D. Scalapino, Phys. Rev. B
45
,
5744 (1992).
3
V. J. Emery, S. A. Kivelson, and O. Zachar, Phys. Rev. B
56
, 6120 (1997).
4
Q. J. Chen, I. Kosztin, B. Janko, and K. Levin, Phys. Rev.
B
59
, 7083 (1999).
5
T. Timusk and B. Statt, Rep. Prog. Phys.
62
, 61 (1999).
6
N.-C. Yeh, Bulletin of Assoc. Asia Pacific Phys. Soc.
12
, 2
(2002), cond-mat/0210656.
7
C.-T. Chen, P. Seneor, N.-C. Yeh, R. Vasquez, L. Bell,
C. Jung, J. Kim, M. Park, H. Kim, and S. Lee, Phys. Rev.
Lett.
88
, 227002 (2002).
8
S. Kleefisch, B. Welter, A. Marx, L. Alff, R. Gross, and
M. Naito, Phys. Rev. B
6310
, 100507 (2001).
9
N.-C. Yeh, C.-T. Chen, R. P. Vasquez, C. U. Jung, S. I.
Lee, K. Yoshida, and S. Tajima, J. Low Temp. Phys.
131
,
435 (2003).
10
V. M. Krasnov, A. Yurgens, D. Winkler, P. Delsing, and
T. Claeson, Phys. Rev. Lett.
84
, 5860 (2000).
11
A. Yurgens, D. Winkler, T. Claeson, S. Ono, and Y. Ando,
Phys. Rev. Lett.
90
, 147005 (2003).
12
N.-C. Yeh and C.-T. Chen, Int. J. Mod. Phys. B (2003),
cond-mat/0302217.
13
S. Sachdev, Rev. Mod. Phys.
75
, 913 (2003).
14
S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan,
J. M. Tranquada, A. Kapitulnik, and C. Howald, Rev.
Mod. Phys.
75
(2003), cond-mat/0210683; and references
therein.
15
M. A. Kastner, R. J. Birgeneau, G. Shirane, and Y. Endoh,
Rev. Mod. Phys.
70
, 897 (1998), and references therein.
16
H. F. Fong, P. Bourges, Y. Sidis, L. P. Regnault, A. Ivanov,
G. Gu, N. Koshizuka, and B. Keimer, Nature
398
, 588
(1999).
17
H. A. Mook, P. C. Dai, and F. Dogan, Phys. Rev. Lett.
88
, 097004 (2002).
18
M. Fujita, H. Goka, K. Yamada, and M. Matsuda, Phys.
Rev. Lett.
88
, 167008 (2002).
19
B. Lake, G. Aeppli, K. N. Clausen, D. McMorrow, K. Lef-
mann, N. Hussey, N. Mangkorntong, M. Nohara, H. Tak-
agi, T. Mason, et al., Science
291
, 1759 (2001).
20
J. E. Hoffman, K. McElroy, D. H. Lee, K. M. Lang,
H. Eisaki, S. Uchida, and J. Davis, Science
297
, 1148
(2002).
21
K. McElroy, R. Simmonds, J. Hoffman, D. Lee, J. Oren-
stein, H. Eisaki, S. Uchida, and J. Davis, Nature
422
, 592
(2003).
22
C. Howald, H. Eisaki, N. Kaneko, M. Greven, and A. Ka-
pitulnik, Phys. Rev. B
67
, 014533 (2003).
23
A. Polkovnikov, M. Vojta, and S. Sachdev, Phys. Rev. B
65
, 220509 (2002).
24
A. Polkovnikov, S. Sachdev, and M. Vojta, Physica C
388
,
19 (2003).
25
D. Podolsky, E. Demler, K. Damle, and B. I. Halperin,
Phys. Rev. B
67
, 094514 (2003).
26
D. G. Zhang and C. S. Ting, Phys. Rev. B
67
, 100506
(2003).
27
Q. H. Wang and D. H. Lee, Phys. Rev. B
67
, 020511
(2003).
28
L. Capriotti, D. J. Sacalapino, and R. D. Sedgewick, cond-
mat/0302563.
29
A. Yazdani, private communications (2003).
5
30
C. Howald, R. Fournier, and A. Kapitulnik, Phys. Rev. B
64
, 100504 (2001).
31
K. M. Lang, V. Madhavan, J. Hoffman, E. Hudson,
H. Eisaki, S. Uchida, and J. Davis, Nature
415
, 412 (2002).
32
O. Zachar, S. A. Kivelson, and V. J. Emery, Phys. Rev. B
57
, 1422 (1998).
33
For coexisting superconductivity and CDW,
V
β
(
q
) repre-
sents the second-order effect of quasiparticle interferenc
e
with pinned CDW. The first-order effect of CDW has been
discussed in Ref.25.
34
M. R. Norman, M. Randeria, H. Ding, and J. C. Cam-
puzano, Phys. Rev. B
52
, 615 (1995).
35
N.-C. Yeh, C.-T. Chen, G. Hammerl, J. Mannhart,
A. Schmehl, C. Schneider, R. Schulz, S. Tajima,
K. Yoshida, D. Garrigus, et al., Phys. Rev. Lett.
87
,
087003 (2001).
36
E. Demler, S. Sachdev, and Y. Zhang, Phys. Rev. Lett.
8706
, 067202 (2001).
37
J. E. Hoffman, E. Hudson, K. Lang, V. Madhavan,
H. Eisaki, S. Uchida, and J. Davis, Science
295
, 466 (2002).
38
Y. Y. Wang, N. P. Ong, Z. Xu, T. Kakeshita, S. Uchida,
D. Bonn, R. Liang, and W. Hardy, Phys. Rev. Lett.
88
,
257003 (2002).