Self-consistent theory of phonon renormalization and electron-phonon ...

1 downloads 0 Views 240KB Size Report
two-dimensional metals with a cylindrical Fermi surface cannot be applied near the wave vector corresponding to the Kohn singularity. Instead, the Dyson ...
PHYSICAL REVIEW B 77, 014517 共2008兲

Self-consistent theory of phonon renormalization and electron-phonon coupling near a two-dimensional Kohn singularity O. V. Dolgov,1 O. K. Andersen,1 and I. I. Mazin2 1Max-Planck-Institut

für Festkörperforschung, Heisenbergstrasse 1, 70569 Stuttgart, Germany Center for Computational Materials Science, Naval Research Laboratory, Washington, DC 20375, USA 共Received 3 October 2007; revised manuscript received 12 December 2007; published 23 January 2008兲

2

We show that the usual expression for evaluating electron-phonon coupling and the phonon linewidth in two-dimensional metals with a cylindrical Fermi surface cannot be applied near the wave vector corresponding to the Kohn singularity. Instead, the Dyson equation for phonons has to be solved self-consistently. If a self-consistent procedure is properly followed, there is no divergency in either the coupling constant or the phonon linewidth near the offending wave vectors, in contrast to the standard expression. DOI: 10.1103/PhysRevB.77.014517

PACS number共s兲: 63.20.K⫺, 74.70.Ad

INTRODUCTION

First principles calculations of the phonon spectra and electron-phonon coupling in MgB2 共see Ref. 1 for a review兲 have determined that the interaction mainly responsible for superconductivity in this material is coupling of small-q, high-energy optical phonons of a particular symmetry with approximately parabolic nearly two-dimensional hole bands forming practically perfect circular cylinders occupying only a small fraction of the Brillouin zone. It was realized early enough2 that the case of an ideal two-dimensional 共2D兲 cylinder leads to a divergency in the calculated phonon linewidth at the 2D Kohn singularity, q = 2kF, which presents serious difficulties in calculating the electron-phonon coupling function. One option that was exploited was to use an analytical integration for the wave vectors comparable with or smaller than 2kF,2,3 and a numerical one for the larger vectors. It was also pointed out1,4 that this singularity gets stronger when the Fermi surface gets smaller, while the integrated electron-phonon coupling 共for 2D parabolic bands兲 does not change. While a perfectly cylindrical Fermi surface is an idealized construction, deviations may be quite small, it seems, on the first glance, unphysical that all phonons with 兩q兩 = 2kF have infinite linewidth. Note that the problem is not specific for MgB2: it occurs in doped graphenes such as CaC6 and, YbC6,5 and, in fact, in any 2D material sporting Kohn singularities. In particular, the hypothetical hexagonal LiB, a subject of substantial recent interest, has ␴ bands that are even more 2D that those in MgB2 共Ref. 6兲 and the described problem is even more pronounced. This intuition is correct. In this paper, we show that close to a Kohn anomaly, standard formulas for calculating electron-phonon interaction 共EPI兲 become incorrect and new, self-consistent expressions replace them. These expressions have no singularities, and exhibit a much more natural, reasonably smooth, q dependence of the phonon self-energy. To start with, we shall remind the readers of the standard formalism. We first define the retarded phonon Green function ␤ 共0兲兴典, D␣␤共q,t兲 ⬅ − i␪共t兲具关uq␣共t兲u−q

where 关…兴 is a commutator and 具…典 denotes statistical averaging. The displacement operator uq␣ in the ␣ direction can 1098-0121/2008/77共1兲/014517共7兲

be expressed via the phonon eigenvectors eq␯ and frequencies squared ␻q2 ␯: uq = 兺 ␯





1 † eq␯共aq␯ + a−q ␯兲. 2M ␻q␯

共1兲

For simplicity, a primitive lattice with a single kind of ions with a single mass M will be considered below. Also, atomic 共hartree兲 units will be used throughout the paper. In this case, the “bare” phonon Green function has a form D0共q, ␻兲 =





1 1 1 − , 2M ␻q0 ␻ − ␻q0 + i␦ ␻ + ␻q0 + i␦

共2兲

where ␻q0 is the bare phonon frequency, before accounting for electron-phonon coupling 共screening by electrons兲. Without losing generality, it can be assumed to be q independent, ␻q0 = ␻0. Correspondingly, the full Green function is D共q, ␻兲 =





1 1 1 − , 2M ␻q ␻ − ␻q + i⌫q ␻ + ␻q + i⌫q

共3兲

where ␻q is the renormalized 共observable兲 frequency, and ⌫q is damping 共phonon linewidth兲 due to EPI.7 The Dyson equation reads D−1共q, ␻兲 = D−1 0 共q, ␻兲 − ⌸共q, ␻兲,

共4兲

where the polarization operator in the lowest approximation 共as usual, the Migdal theorem8 allows neglecting the vertex corrections兲 along the real frequency axis at T = 0 has a form 共a is the lattice constant in the plane兲 ⌸共q, ␻兲 = − 2i





冉 冊

␻ q 0 兩gk,k+q 兩2G0 k + ,␧ + 2 2

␻ a2d2k d␧ q ⫻G0 k − ,␧ − . 2 2 共2␲兲2 2␲

冊 共5兲

0 Here, gk,k+q is the bare electron-ion scattering matrix element 共the commonly used EPI matrix element differs in that the potential gradient is replaced by the derivatives with respect to the normal phonon coordinates兲

014517-1

©2008 The American Physical Society

PHYSICAL REVIEW B 77, 014517 共2008兲

DOLGOV, ANDERSEN, AND MAZIN

0 gk,k+q =



1 a2

a2

* ␺k+q 共r兲 ⵜ V共r兲␺k共r兲d2r,

共6兲

where G0共k,␧兲 =

1

共7兲

␧ − ␧0共k兲 + i␦ sign共k − kF兲

is the bare electron Green function, and ␧0共k兲 = 共k2x + k2y − kF2 兲 / 2m. The renormalized phonon frequency and the phonon linewidth ⌫q are determined by the pole of the phonon Green function D共q , ␻兲 or D−1 0 共q, ␻q

D 共q, ␻q + i⌫q兲 ⬅ −1

+ i⌫q兲 − ⌸共q, ␻q + i⌫q兲 = 0

I. COMPLEX POLARIZATION OPERATOR

Let us consider a model with a cylindrical Fermi surface of radius kF, whose electrons interact with an optical phonon with a bare frequency ␻0, and a momentum-independent matrix element g0 共we also neglect possible warping of the Fermi-surface cylinder, cf. Ref. 11兲. In this case, the imaginary part of Eq. 共10兲 reads Im ⌸共q, ␻ + i␦兲 = − 2␲兩g0兩2 兺 兵␪关␧0共k + q兲兴 − ␪关␧0共k兲兴其 k

⫻␦关␧0共k + q兲 − ␧0共k兲 − ␻兴 or

关note that in the experiment, the phonon linewidth is usually defined as the half-width of the peak in Im D共q , ␻兲, which differs from our definition by terms of the order 共⌫ / ␻兲4兴. This leads to

␻q2

=

1 Re ⌸共q, ␻q + i⌫q兲 + ⌫q2 + M

␻20

Im ⌸共q, ␻兲 = −

共8兲

=− 1 Im ⌸共q, ␻q + i⌫q兲. 2M ␻q

0 兩2 ⌸共q, ␻ + i␦兲 = − 2 兺 兩gk,k+q k

⌸共q,Z兲 =

共10兲

2 ␻appq

and ⌫qapp = −

冋 冉 冊册

␻2 1 = ⌸共q,0兲 1 + O 20 M ␧F



1 d Im ⌸共q, ␻兲 2M d␻





冑1 − 共q/2k + m␻/kq兲2

m兩g0兩2 关Re 冑1 − 共Q − ⍀兲2 ␲ a 2Q 共13兲

1 P ␲





−⬁

dE Im ⌸共q,E兲. E−Z

The result is ⌸共q,Z兲 = −

共11兲

m兩g0兩2 ␲ a 2Q



dxdy

C

x , − 共Z/4␧F兲2/Q2 + x2

共14兲

where ␧F = kF2 / 2m, x = kx / kF, y = ky / kF, and N共0兲 = m / 2␲a2 is the density of states at the Fermi level, per spin. The integration is performed over the range 共x − Q兲2 + y 2 ⬍ 1. The substitution x − Q = r cos ␸ and y = r sin ␸ leads to

␻=0

␲ 兺 兩g0 兩2␦关␧0共k + q兲兴␦关␧0共k兲兴, M k k,k+q

max兵kF,m␻/q+q/2其

where Q = q / 2kF and ⍀ = ␻ / qvF 共vF is the Fermi velocity兲. To find the polarization operator for a complex frequency Z = ␻ + i␥, we use the Hilbert transformation

nk − nk+q ␧0共k + q兲 − ␧0共k兲 − ␻ − i␦

to first order in frequency just above the real axis of the complex frequency ␻. In this case,

2 kF +2m␻

− Re 冑1 − 共Q + ⍀兲2兴,

共9兲

The next standard step, following Ref. 9, is to expand the polarization operator

冕冑

dk

⫻Re

and ⌫q = −

m兩g0兩2 2 ␲ a 2q



共12兲

where the factor 共nk − nk+q兲␦关␧0共k + q兲 − ␧0共k兲 − ␻兴 has been changed to −␻␦关␧0共k + q兲兴␦关␧0共k兲兴. According to Ref. 10, “Except for extremely pathological energy bands, it is an excellent approximation.” Unfortunately, MgB2 and some other recently discovered superconductors are examples where the energy bands are, in some aspects, pathological. Formally, expression 共12兲 for a 2D system is divergent near a Kohn anomaly q → 2kF, and we have ⌫q 艌 ␻q, i.e., phonons are not well defined quasiparticles. To describe electronphonon interaction in these systems, we have to calculate the polarization operator ⌸共q , Z兲 for a complex frequency Z and solve Eqs. 共8兲 and 共9兲.

C

⌸共q,Z兲 = −

冕 冕 1

dxdy ⇒

冋 冉 册

0

m兩g0兩2 Z 2 2− 1− 2␲a 4␧FQ2

− 兵Z → − Z其 .

2␲

rdr

d␸ ,

0

冊冑 冉

1− Q−

Z 4␧FQ



−2

共15兲

The branches of the square roots are chosen so as to get the correct behavior ⌸ ⬀ q2VF2 / ␻2 at large frequencies ␻ 共␻ = Re Z兲. For a 2D system for ⌸共q , ␻兲 on the real frequency axis, one can write:12–15

014517-2

PHYSICAL REVIEW B 77, 014517 共2008兲

SELF-CONSISTENT THEORY OF PHONON…

Re ⌸共q, ␻ + i␦兲 = −

m兩g0兩2 关2Q − 共Q − ⍀兲Re 冑1 − 共Q − ⍀兲−2 2 ␲ a 2Q

− 共Q + ⍀兲Re 冑1 − 共Q + ⍀兲−2兴, Im ⌸共q, ␻ + i␦兲 = −

m兩g0兩2 关Re 冑1 − 共Q − ⍀兲2 2 ␲ a 2Q

− Re 冑1 − 共Q + ⍀兲2兴.

共16兲

A three-dimensional 共3D兲 plot of the functions 共16兲 can be found, for instance, in Ref. 16. Note that while the expression for real frequencies 关Eq. 共16兲兴 has been published multiple times, we have not found in the literature the general expression for an arbitrary complex frequency argument 关Eq. 共15兲兴. First, we see that the imaginary part is finite for all values of the wave vector q and vanishes inside the Landaudamping cone q ⬍ ␻ / vF 关more exactly, at Q ⬍ 共␻ / 4␧F兲 − 共␻ / 4␧F兲2兴. It has two maxima: one is rather small, Im ⌸共␻兲 ⯝ −m兩g0兩2冑␻ / 4␧F / ␲a2, at Q ⯝ 1 − 共␻ / 4␧F兲, while the other has an antiadiabatical behavior, Im ⌸共␻兲 ⯝ −m兩g0兩2冑2␧F / ␻ / ␲a2, and occurs at a very low frequency Q ⯝ 共␻ / 4␧F兲 + 共␻ / 4␧F兲2. However, if we expand the polarization operator at small frequencies, we recover a standard result 共see, e.g., Refs. 2 and 4兲 ⌸app共q, ␻ + i␦兲 ⯝ − −



FIG. 1. 共Color online兲 The imaginary and real parts of the normalized polarization operator 2␲兩⌸共q , ␻ = ␻0兲兩a2 / m兩g0兩2 as a function of the reduced wave vector q / kBZ, for different fillings kF / kBZ. Solid lines represent the exact results, and the dashed lines the approximate solution 关Eq. 共17兲兴. The four different sets correspond, from left to right, to four increasing values of kF / kBZ.

Re ⌸共q,0兲 ⯝ −

m兩g0兩2 i␻/4␧F 1+ ␪共1 − Q兲 2 ␲a Q冑1 − Q2



1−



1 ␪共Q − 1兲 , Q2

=−



1 . ␧0共k + q兲 − ␧0共k兲

In this case,

⌸ M 共q,i␻n兲 = −

m兩g0兩2 关1 − ␪共Q − 1兲冑1 − Q−2兴. ␲a2

Re ⌸共q → 0, ␻兲 ⬇ −

共18兲

In the opposite limit, Im ⌸共q , ␻兲 ⬅ 0 for q 艋 m␻ / 冑kF2 − m␻. For the real part, one can set ␻ = 0 in Eq. 共17兲: k

2kF q

1−

2

共19兲

In the opposite limit,

02

Re ⌸共q,0兲 = − 2兩g0兩2 兺 ␪共兩k兩 − kF兲␪共kF − 兩k + q兩兲

冑 冉 冊册

共17兲

where the imaginary part diverges at q → 0 and q → 2kF. The real part of Eq. 共16兲 practically coincides with the real part of Eq. 共17兲 except inside the Landau-damping region. At ␻ → 0 and finite q, we get

␻ m兩g 兩 ␪共1 − Q兲 . Im ⌸app共q, ␻兲 ⯝ − 4␧F ␲a2 Q冑1 − Q2



m兩g0兩2 1 − ␪共q − 2kF兲 ␲a2

冉 冊

m兩g0兩2 qkF 2␲a2 m␻

2

.

The momentum dependence of the absolute values of the imaginary 共the upper panel兲 and real parts 共the lower panel兲 of Eqs. 共16兲 共solid lines兲 and 共17兲 共short-dash lines兲 at ␻ = ␻0 = 90 meV as the functions of the reduced wave vector q / kBZ 关kBZ = ␲ / a is the Brillouin zone 共BZ兲 vector, or the radius of the Wigner-Seitz cylinder兴 is shown in Fig. 1. Fermi vectors kF / kBZ = 0.075, kF / kBZ = 0.1, kF / kBZ = 0.15, and kF / kBZ = 0.2 correspond to ␧F = 0.15, 0.27, 0.60, and 1.07 eV, respectively. Along the imaginary 共Matsubara兲 axis, the polarization operator has the following form:





冑Q4 − Q2 − 共␻n/4Q␧F兲 + 冑关Q4 − Q2 − 共␻n/4Q␧F兲兴2 + 共Q␻n/4␧F兲2 m兩g0兩2 1+ , 2 冑2Q2 ␲a 014517-3

共20兲

PHYSICAL REVIEW B 77, 014517 共2008兲

DOLGOV, ANDERSEN, AND MAZIN

where ␻n = 2␲nT . T is temperature and n = 0 , ⫾ 1, ⫾2 , . . ., ⫾⬁. ⌸ M 共q , i␻n → 0兲 coincides with Eq. 共19兲. II. PHONON RENORMALIZATION IN TWO-DIMENSIONAL SYSTEMS

First, let us consider the approximate polarization operator from Eq. 共17兲. Then

␻qapp = ␻0冑1 − 2␨关1 − ␪共q − 2kF兲冑1 − 共2kF/q兲2兴

共21兲

and, according to Eqs. 共12兲 and 共17兲, ⌫qapp =

m兩g0兩2 ␨␻20 ␪共1 − Q兲 ␪共2kF − q兲 = , 2 2␲a MVFq 冑1 − 共q/2kF兲2 4␧F Q冑1 − Q2 共22兲

where we have introduced, following Ref. 17, an auxiliary coupling constant ␨ = N共0兲兩g0兩2 / M ␻20 共some authors use another dimensionless constant ␭0 = 2␨兲. ⍀q = ␻qapp + i⌫qapp gives the renormalized frequency and the damping 关see Eqs. 共11兲 and 共12兲兴. For q 艋 2kF, ⌫qapp

=

␻qapp

␵␻0

1

4␧F冑1 − 2␨ Q冑1 − Q2

共23兲

.

This expression diverges in the limits Q → 0 and Q → 1. Turning now to the exact Eq. 共15兲, we observe that in the lowest order in ␻ / ␧F ⌸共q → 2kF, ␻兲 = −



m兩g0兩2 1 − 共1 + i兲 ␲a2

冑 册

␻ . 8␧F

This leads to

␻2kF ⯝ ␻0冑1 − 2␨ , ⌫2kF =

␨␻20 ␻2kF

and ⌫2kF

␻2kF

=

␨ 1 − 2␨





␻2kF 8␧F

,

␻0冑1 − 2␨ . 8␧F

共24兲

This ratio remains finite in the limit Q → 1, although the approximate expression of Eq. 共23兲 diverges for any system with a cylindrical Fermi surface. Of course, the singularity at Q → 0 is also unphysical and, in principle, can be treated in a similar way. The common point is that, in both cases, the well-known popular formula 兩g20兩 ⌫q 兩g20兩 = ␦共␧k兲␦共␧k+q兲 = 兺 ␻q M k M共2␲a兲2



dk 兩v共k兲 ⫻ v共k + q兲兩 共25兲

is not valid near the singularity. We do not consider the Q → 0 case in this paper, but it is worth noting that there are other issues relevant in that limit, but not near the Kohn singularity, such as the Landau threshold and applicability of the Migdal theorem.

FIG. 2. 共Color online兲 The linewidth ⌫q and the renormalized phonon frequency ␻q obtained by using the approximation of Eq. 共17兲 and the exact expression Eq. 共15兲, for the following parameters: the bare phonon frequency ␻0 = 90 meV and the bare constant of EPI ␨ = 1 / 4. The filling corresponds to kF / kBZ = 0.17 and ␧F = 0.2 eV.

The results for the linewidth ⌫q and the renormalized phonon frequency ␻q obtained by using the approximate polarization operator as functions of the reduced wave vector Q = q / 2kF are shown in Fig. 2 by red lines. Parameters are the following: the bare phonon frequency ␻0 = 90 meV and the bare constant of EPI ␨ = 1 / 4. The ratio kF / kBZ is equal to 0.17. It corresponds to ␧F = 0.2 eV. The approximate result agrees with that of Ref. 4 共their Fig. 4兲. The exact result 共black lines兲 has been obtained by a numerical solution of Eqs. 共8兲 and 共9兲 using the polarization operator from Eq. 共15兲. The latter, in contrast to the approximate expression, shows two shoulders in the wave vector dependence of the renormalized frequency ␻q 共see also 兩Re ⌸共q , ␻0兲兩 in the bottom panel of Fig. 2兲. The first one corresponds to the maximum of 兩Im ⌸共q , ␻兲兩 共or ⌫q兲 and the second one to the vanishing of these values. III. ELECTRON SELF-ENERGY

The electron self-energy is expressed via the electron and phonon Green functions: ⌺共p兲 = ig0



G共p − k兲D共k兲⌫共p,p − k,k兲

dD+1k , 共2␲兲D+1

where p = 兵p , ␧其. It was shown by Migdal8 that one can neglect the vertex corrections ⌫共p , p − k , k兲 ⯝ g0共1 + 冑m / M兲 and that the function G共␧ , k兲 differs from the bare electron Green function 关Eq. 共7兲兴 only in the narrow interval of momenta 兩k − kF兩 ⱗ ␻ ph / VF and frequencies 兩␻兩 ⱗ ␻ ph. Thus, the full electron Green function G共p兲 can be substituted by the corresponding function for noninteracting electrons 关Eq. 共7兲兴. Using Eq. 共3兲, the electron self-energy for T = 0 becomes 共see, e.g., Refs. 18 and 19兲

014517-4

PHYSICAL REVIEW B 77, 014517 共2008兲

SELF-CONSISTENT THEORY OF PHONON…

⌺共k, ␻兲 = 兺 ␦共␧k+q兲 q

+

兩g0兩2 2M ␻q

冕 冋 d␰



␪共␰兲 ␻ − ␰ − ␻q + i⌫q

˜ 共0兲 = N q

␪共− ␰兲 . ␻ − ␰ + ␻q − i⌫q

1 兺 ␦共␧k兲⌺共k, ␻兲 N共0兲 k

=−



兩g0兩2 ⌫2 + 共␻q − ␻兲2 1 ␦共␧k兲␦共␧k+q兲 ln 2q 兺 兺 N共0兲M k q 4␻q ⌫q + 共␻q + ␻兲2

冋 冉 冊 冉 冊册冎

+ 2i tan−1

␻q − ␻ ␻q + ␻ − tan−1 ⌫q ⌫q

˜ 共0兲 = N q

q

兩g0兩2 1 , M ⌫q2 + ␻q2

02

共26兲

=





␻q2 ,

⌫qapp 1 ␦共␻ − ␻q兲, = 兺 2␲N共0兲 q ␻q





共30兲

the same function Im ⌸⬘共q , 0兲 ⬅ d Im ⌸共q , ␻兲 / 兩d␻兩␻=0. In a general case, these functions can be different. This result, without using the pole approximation for the phonon Green function, can be trivially obtained in the Matsubara formalism. In this case, one does not need to solve the Dyson equation. In the lowest order in coupling for T = 0, the self-energy has a form

⌺共i␧,k兲 = − i

共31兲

where in the last equality we have used the approximate Eq. 共12兲. Note that Eq. 共31兲 is a consequence of the fact that the damping ⌫qapp, according to Eq. 共12兲, can be expressed via the nesting function 关Eq. 共27兲兴 and both are determined by

共29兲



˜ 共0兲兩g0兩2 1 N 1 ⌫q ⌫q q − . 兺 2␻q 2 共␻ − ␻q兲2 + ⌫q2 共␻ + ␻q兲2 + ⌫q2 2␲ M q

␦共␧k兲兩g0兩2␦共␧k+q兲 1 ␦ 共 ␻ − ␻ q兲 兺 2␻q N共0兲M k,q

共28兲

but for ⌫q Ⰷ ␻q, the contribution of strongly damped phonons to total ␭q is suppressed. The result 共26兲 we can get also if we introduce, according to Eq. 共3兲, a generalized Eliashberg function

␦共␧k兲兩g0兩2␦共␧k+q兲 1 1 ⌫q ⌫q 兺 2 − 2 2 共␻ − ␻q兲 + ⌫q 共␻ + ␻q兲2 + ⌫q2 2␻q 2␲ MN共0兲 k,q

The second term in this expression cancels out the nonphysical behavior at low and high frequencies. Otherwise, ␭⌫ would have been divergent. Equations 共26兲 and 共30兲 are general and valid not only for the 2D systems, where phase space factor 关Eq. 共27兲兴 is divergent. For ⌫q Ⰶ ␻q, we have 2 ␣app 共␻兲F共␻兲 =

,

In the weak-damping approximation, we recover the standard formula ˜ 共0兲 兩g 兩 ␭qapp = N q M

where we introduced the phase space function 共sometimes called “nesting function”兲

␣⌫2 共␻兲F共␻兲 =

4␲␧FQ冑1 − Q2

02 1 ˜ 共0兲 兩g 兩 . ␭q = N q 2 M ⌫q + ␻q2

1 1 兺 兺 ␦共␧k兲兩g0兩2␦共␧k+q兲 ⌫2 + ␻2 MN共0兲 k q q q

˜ 共0兲 =兺N q

␪共1 − Q兲

which diverges at Q → 0 and Q → 1. Following Ref. 9, we can introduce the “mode ␭” via the expression ␭ = 兺q␭q. Then

.

The limit ␭ = −lim␻→0 Re ⌺共␻兲 / ␻ is nothing but the standard electron-phonon coupling constant ␭⌫ =

共27兲

For a 2D cylindrical Fermi surface, we have

This is a trivial generalization of the standard expressions on the finite phonon linewidth case. Let us average the selfenergy over the Fermi surface ⌺共␻兲 =

1 兺 ␦共␧k兲␦共␧k+q兲. N共0兲 k



where

014517-5



d共i␻兲 兺 兩g0兩2DM␯共i␻ + i␧,k,k⬘兲 2 ␲ k ,␯ ⬘

1 , i␻ − ␧k⬘ − ⌺共i␻,k⬘兲

PHYSICAL REVIEW B 77, 014517 共2008兲

DOLGOV, ANDERSEN, AND MAZIN

D M ␯共i␻,k,k⬘兲 = 共1/␲兲





d⍀ Im D共⍀ + i␦,k,k⬘兲关共i␻ − ⍀兲−1

0

− 共i␻ + ⍀兲−1兴.

共32兲

We can also average the self-energy over the Fermi surface ⌺共i␧兲 = 兺k␦共␧k兲⌺共i␧ , k兲 / N共0兲: ⌺共i␧兲 =

1 兺 兺 ␦共␧k兲 N共0兲 k q ⫻兩g0兩2␦共␧k+q兲





−⬁

d␻ D M 共q,i␻兲 2␲





−⬁

d␰ . i共␧ − ␻兲 − ␰

pling constants ␭ and ␨ are that they are measures of the renormalization of the phonon frequency from ␻0 to ␻q 共cf. a discussion for 3D systems in Ref. 20兲. Turning back to the 2D case, according to Eq. 共26兲, we can neglect all divergent contributions near Q = 0 and Q = 1. Elsewhere, we can use Eq. 共29兲. One should keep in mind that the conventional coupling constant is ␭ = ␨ / 共1 − 2␨兲. This parameter determines electronic properties 共Fermi velocities, Tc, etc.兲. The other parameter, ␨ = ␭ / 共1 + 2␭兲 ⬍ 1 / 2, defines the observable phonon frequency, ␻q = ␻0冑1 − 2␨. CONCLUSIONS

d␰ ⬁ 兰−⬁ i共␧−␻兲−␰ = −i2␲

The integral sign共␧ − ␻兲 allows us to calculate the physical coupling constant ␭ ␭=−



⳵⌺el共i␧兲 ⳵i␧



= − lim ␧→0

⫻兩g0兩2␦共␧k+q兲



␧→0

␣2共␻兲F共␻兲 =

1 兺 兺 ␦共␧k兲 N共0兲 k q



d␻ D共q,i␻兲2␲␦共␧ − ␻兲 2␲

1 =− 兺 兺 ␦共␧k兲兩g0兩2␦共␧k+q兲DM共q,i0兲 N共0兲 k q =−

First of all, the standard well-known expression

1 1 . ␦共␧k兲兩g0兩2␦共␧k+q兲 −1 兺 兺 N共0兲 k q D M0共i0兲 − ⌸M 共q,i0兲

On the Matsubara axes, for a 2D system, according to Eq. 共20兲, the phase space factor vanishes for q ⬎ 2kF and ⌸ M 共q 艋 2kF , i0兲 = −2␨␻20 is a constant 关see Eqs. 共20兲 and 共19兲兴. Us2 ing Eqs. 共2兲 and 共32兲, we get D−1 M0共i0兲 = −M ␻0 and ␭ = ␨/共1 − 2␨兲.

1 ⌫q ␦ 共 ␻ − ␻ q兲 兺 2␲N共0兲 q ␻q

␦共␧k兲兩g0兩2␦共␧k+q兲 1 ␦ 共 ␻ − ␻ q兲 兺 兺 2␻q N共0兲M k q

is valid only in the lowest order in the bare phonon linewidth, ⌫q, which is not an acceptable approximation in the case of strong Kohn singularities, and particularly for a cylindrical Fermi surface. ⌫q in this approximation is not the actual phonon linewidth; as opposed to ⌫q, determined by the oversimplified Eq. 共12兲, the real phonon linewidth does not diverge even for an ideally cylindrical Fermi surface. Second, following Eq. 共33兲 derived above, the more accurate treatment presented here still respects the old sum rules that the total ␭ for a cylindrical Fermi surface and parabolic bands do not depend on filling.

共33兲

ACKNOWLEDGMENT

In a 3D case, ⌸ M 共q , i0兲 is a rather complicated function of q and Eq. 共33兲 is only an approximation 共as was probably first mentioned by Fröhlich17兲. The physical meanings of the cou-

I.I.M. would like to thank the Max Planck Society for hospitality during his visit to the Max Planck Institute for Solid State Research, where part of this work was done.

I. I. Mazin and V. P. Antropov, Physica C 385, 49 共2003兲. Kong, O. V. Dolgov, O. Jepsen, and O. K. Andersen, Phys. Rev. B 64, 020501共R兲 共2001兲. 3 A. Y. Liu, I. I. Mazin, and J. Kortus, Phys. Rev. Lett. 87, 087005 共2001兲. 4 J. M. An, S. Y. Savrasov, H. Rosner, and W. E. Pickett, Phys. Rev. B 66, 220502共R兲 共2002兲; W. E. Pickett, J. M. An, H. Rosner, and S. Y. Savrasov, Physica C 387, 117 共2003兲. 5 For a review, see I. I. Mazin, L. Boeri, O. V. Dolgov, A. A. Golubov, G. B. Bachelet, M. Giantomassi, and O. K. Andersen, Physica C 460, 116 共2007兲. 6 A. Y. Liu and I. I. Mazin, Phys. Rev. B 75, 064510 共2007兲; M. Calandra, A. N. Kolmogorov, and S. Curtarolo, ibid. 75, 144506 共2007兲. 7 In publications, one can sometimes find the phonon Green functions which differ from ours in some factors. This ambiguity can be traced down to the definitions of the phonon field operators

关Eq. 共1兲兴. In this case, the corresponding factors appear in the electron-ion matrix element 关Eq. 共6兲兴. However, all physical quantities such as EPI coupling constant, Eliashberg functions, etc., are determined by the unique combination 兩g0兩2D, which is independent of definitions. 8 A. B. Migdal, Sov. Phys. JETP 7, 996 共1958兲. 9 P. B. Allen, in Dynamical Properties of Solids, edited by G. K. Horton and A. A. Maradudin, Vol. 3 共North-Holland, Amsterdam, 1980兲, p. 157. 10 P. B. Allen, Phys. Rev. B 6, 2577 共1972兲. 11 M. Calandra and F. Mauri, Phys. Rev. B 71, 064501 共2005兲. 12 F. Stern, Phys. Rev. Lett. 18, 546 共1967兲. 13 T. Ando, A. B. Fowler, and F. Stern, Rev. Mod. Phys. 54, 445 共1982兲. 14 Y. Zhang, V. M. Yakovenko, and S. Das Sarma, Phys. Rev. B 71, 115105 共2005兲. 15 A. Isihara and T. Toyoda, Z. Phys. B 23, 389 共1976兲; A. Isihara,

1

2 Y.

014517-6

PHYSICAL REVIEW B 77, 014517 共2008兲

SELF-CONSISTENT THEORY OF PHONON… in Solid State Physics, edited by H. Erenreich, F. Zeitz, and D. Turnbull 共Academic, New York, 1989兲, Vol. 42, p. 271. 16 M. Dressel and G. Grüner, Electrodynamics of Solids: Optical Properties of Electrons in Matter 共Cambridge University Press, Cambridge, England, 2002兲. 17 H. Fröhlich, Phys. Rev. 79, 845 共1950兲. 18 G. D. Mahan, Many-Particle Physics 共Plenum, New York, 1981兲.

19 P.

B. Allen and B. Mitrović, in Solid State Physics, edited by H. Erenreich, F. Zeitz, and D. Turnbull 共Academic, New York, 1982兲, Vol. 37, p. 1. 20 E. G. Maksimov and D. I. Khomskii, in High-Temperature Superconductivity, edited by V. L. Ginzburg and D. A. Kirzhnits 共Consultants Bureau, New York, 1982兲, Chap. 3.

014517-7