## Antiadiabatic Phonons and Superconductivity in Eliashberg McMillan Theory

M.V. Sadovskii

Received: date / Accepted: date

Abstract The standard Eliashberg - McMillan theory of superconductivity is essentially based on the adiabatic approximation. Here we present some simple estimates of electron - phonon interaction within Eliashberg - McMillan approach in non - adiabatic and even antiadiabatic situation, when characteristic phonon frequency Ω 0 becomes large enough, i.e. comparable or exceeding the Fermi energy E F . We discuss the general definition of Eliashberg - McMillan (pairing) electron phonon coupling constant λ , taking into account the finite value of phonon frequencies. We show that the mass renormalization of electrons is in general determined by different coupling constant ˜ , λ which takes into account the finite width of conduction band, and describes the smooth transition from the adiabatic regime to the region of strong nonadiabaticity. In antiadiabatic limit, when Ω 0 /greatermuch E F , the new small parameter of perturbation theory is λ E F Ω 0 ∼ λ D Ω 0 /lessmuch 1 ( D is conduction band half - width), and corrections to electronic spectrum (mass renormalization) become irrelevant. However, the temperature of superconducting transition T c in antiadiabatic limit is still determined by Eliashberg - McMillan coupling constant λ . We consider in detail the model with discrete set of (optical) phonon frequencies. A general expression for superconducting transition temperature T c is derived, which is valid in situation, when one (or several) of such phonons becomes antiadiabatic. We also analyze the contribution of such phonons into the Coulomb pseudopotential µ /star and show, that antiadiabatic phonons do not contribute

M.V. Sadovskii

Institute for Electrophysics, Russian Academy of Sciences, Ural Branch, Amundsen str. 106, Ekaterinburg 620016, Russia to Tolmachev's logarithm and its value is determined by partial contributions from adiabatic phonons only.

Keywords Eliashberg - McMillan theory · Electron phonon interaction · Antiadiabatic phonons · Coulomb pseudopotential · Critical temperature

## 1 Introduction

Eliashberg - McMillan superconductivity theory is currently the basis for microscopic description of Cooper pairing and all general properties of conventional superconductors [1,2,3,4,5]. It is is essentially based on adiabatic approximation and Migdal's theorem [6], which allows to neglect the vertex corrections to electron -phonon coupling in typical metals. The actual small parameter of perturbation theory is λ Ω 0 E F /lessmuch 1, where λ is the dimensionless Eliashberg - McMillan electron phonon coupling constant, Ω 0 is characteristic phonon frequency and E F is Fermi energy of electrons. This leads to the widely accepted opinion, that vertex corrections can be neglected even for the case of λ &gt; 1, due to the fact, that in common metal Ω 0 E F /lessmuch 1 . The possible breaking of Migdal's theorem for the case of λ ∼ 1 due to polaronic effects was widely discussed in the literature [7,8]. In the following we consider only the case of λ &lt; 1 where we can safely neglect these effects [8].

Recently a number of superconductors was discovered, where the adiabatic approximation is not necessarily valid, and characteristic frequencies of phonons are of the order or even greater than Fermi energy. In this respect we can mention single - atomic layers of FeSe on the SrTiO 3 substrate (FeSe/STO) [9], as well as record breaking hydride based superconductors at high pressures [10]. This is also the case in the long -

standing puzzle of superconductivity in doped StTiO 3 [11]. The role of nonadiabatic phonons was recently analyzed in important papers by Gor'kov [12,13] within the standard BCS - like weak - coupling approach, and directly addressed to these new superconductors. Here we review some further estimates, derived by us in Refs. [14,15] in the framework of Eliashberg - McMillan theory.

## 2 Electron self - energy and electron - phonon coupling constant

Let us consider first a metal in normal (non superconducting) state, which is sufficient to introduce some basic notions of Eliashberg - McMillan theory [2,3]. The second - order (in electron - phonon coupling) diagram is shown in Fig. 1. Making all calculations in finite temperature technique, after the analytic continuation from Matsubara to real frequencies iω n → ± ε iδ and in the limit of T = 0 (i.e. E F /greatermuch T ), the contribution of diagram in Fig. 1 can be written [1,2] as:

$$\text{gram in fig. 1 can be written } [ 1, 2 ] \text{ as: } \text{ ada} \\ \Sigma ( \varepsilon, \mathbf p ) = \sum _ { \mathbf p ^ { \prime }, \alpha } | g _ { \mathbf p \mathbf p ^ { \prime } } ^ { \alpha } | ^ { 2 } \left \{ \frac { f _ { \mathbf p ^ { \prime } } } { \varepsilon - \varepsilon _ { \mathbf p ^ { \prime } } + \Omega _ { \mathbf p - \mathbf p ^ { \prime } } ^ { \alpha } - i \delta } \text{ \quad } \text{ theiv} \\ + \frac { 1 - f _ { \mathbf p ^ { \prime } } } { \varepsilon - \varepsilon _ { \mathbf p ^ { \prime } } - \Omega _ { \mathbf p - \mathbf p ^ { \prime } } ^ { \alpha } + i \delta } \right \} \text{ \quad } ( 1 ) \text{ initial} \\ \text{where in notations of Fig. 1 } \mathbf p ^ { \prime } = \mathbf p + \mathbf q. \text{ Here } g _ { \mathbf n. \mathbf n ^ { \prime } } ^ { \alpha } \text{ is } \text{ ores}$$

where in notations of Fig, 1 p ′ = p + q . Here g α p p , ′ is Fr¨hlich electron - phonon coupling constant, o ε p is electronic spectrum with energy zero taken at the Fermi level, Ω α q represents the phonon spectrum, and f p is Fermi distribution (step - function at T = 0). In these expressions index α enumerates the branches of phonon spectrum, which below is just dropped for brevity.

Now we can essentially follow the analysis, presented in Ref. [2,3]. Eq. (1) can be identically rewritten as:

$$& \Pi \intertext { n e. } [ 4, \vartheta ]. \, \varkappa \Delta \, \{ 1 \, \text{ can use memory new use} \, as. \\ & \Sigma ( \varepsilon, \mathbf p ) = \int d \omega \sum _ { \mathbf p ^ { \prime } } | g _ { \mathbf p \mathbf p ^ { \prime } } | ^ { 2 } \delta ( \omega - \Omega _ { \mathbf p - \mathbf p ^ { \prime } } ) \times \quad & \text{After} \\ & \times \left \{ \frac { f _ { \mathbf p ^ { \prime } } } { \varepsilon - \varepsilon _ { \mathbf p ^ { \prime } } + \omega - i \delta } + \frac { 1 - f _ { \mathbf p ^ { \prime } } } { \varepsilon - \varepsilon _ { \mathbf p ^ { \prime } } - \omega + i \delta } \right \} \quad & ( 2 ) \quad & \text{int fact} \\ & \text{average} \\ & \text{To simplify calculations we can get rid of explicit $\mathbf$mo$} \quad \text{which}$$

To simplify calulations we can get rid of explicit momentum dependencies here by averaging the matrix element of electron - phonon interaction over surfaces of constant energies, corresponding to initial and final momenta p and p ′ , which usually reduces to the averaging over corresponding Fermi surfaces, as phonon scattering takes place only within the narrow energy interval close to the Fermi level, with effective width of the order of double characteristic frequency of phonons 2 Ω 0 , and taking into account that in typical metals we always have Ω 0 /lessmuch E F .

Fig. 1 Second - order diagram for self - energy. Dashed line -phonon Green's function D (0) , continuous line - electron Green's function G in Matsubara representation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00718v1-figure-089b5245-8125-44ca-90c8-525c70e52e32.jpg)

This averaging can be achieved by the following replacement in Eq. (2):

$$\text{pacement in req. (4)} \\ \text{rcon-} \\ \text{e ba-} \\. The } ^ { 1 } _ { \text{P} } \sum _ { \text{P} } \frac { 1 } { N ( 0 ) } \sum _ { \text{P} ^ { \prime } } | g _ { \text{P} ^ { \prime } } | ^ { 2 } \delta ( \omega - \Omega _ { \text{P} - \text{P} ^ { \prime } } ) \delta ( \varepsilon _ { \text{P} } ) \delta ( \varepsilon _ { \text{P} ^ { \prime } } ) \\ \text{gram} \\ \text{item-} } ^ { 1 } \equiv \frac { 1 } { N ( 0 ) } \alpha ^ { 2 } ( \omega ) F ( \omega ) \ ( 3 ) \\ \text{from} \quad \text{where in the last expression we have introduced the}$$

where in the last expression we have introduced the definition of Eliashberg function α 2 ( ω ) and F ω ( ) = ∑ q δ ω ( -Ω q ) is the phonon density of states. In non adiabatic case, when phonon energy becomes comparable with or even exceeds the Fermi energy, electron scattering is effective not only in the narrow energy layer around the Fermi surface, but in much wider energy interval of the order of Ω 0 ∼ E F . Then, for the case of initial | p | ∼ p F the averaging over p ′ in expression like (3) should be done over the surface of constant energy, corresponding to E F + Ω p -p ′ , as is shown in Fig. (2). Now the Eq. (3) is directly generalized as:

$$s \text{ ele- \quad } & \text{ Now the Eq. (3) is directly generalized as: } \\ \ f _ { \text{p is } } & \text{$|g_{\text{p} ^ { \prime } |^{2} \delta(\omega - \Omega_{\text{p} - \text{p} ^ { \prime } } \Longrightarrow \\ \text{these} & \text{$\frac{1}{N(0)} \sum_{p} \frac{1}{N(0)} \sum_{p^{\prime}}|g_{\text{p} ^ { \prime } |^{2} \times \\ y$.} \\ \text{ented} \\ \text{as: } } & \quad \cong \frac { 1 } { N(0)} \alpha ^ { 2 } ( \omega ) F ( \omega ) \\ & \text{After the replacement like (3) or (4) the explicit mo-}$$

After the replacement like (3) or (4) the explicit momentum dependence of the self - energy disappears and in fact in the following we are dealing with Fermi surface average of self - energy Σ ε ( ) ≡ 1 N (0) ∑ p δ ε ( p ) Σ ε, ( p ), which is now written as:

$$\Phi \text{mov} \quad \text{ which is now wrttten as:} \\ \text{rix el-} \\ \text{rfaces} \quad \Sigma ( \varepsilon ) = \int d \varepsilon ^ { \prime } \int d \omega \alpha ^ { 2 } ( \omega ) F ( \omega ) & \left \{ \frac { f ( \varepsilon ^ { \prime } ) } { \varepsilon - \varepsilon ^ { \prime } + \omega - i \delta } \\ \text{e av-} \\ \text{horon} \\ \text{energy} \\ \text{dth of} \quad \text{ This expression forms the basis of Eliashberg - McMil-}$$

This expression forms the basis of Eliashberg - McMillan theory and determines the structure of Eliashberg equations for the description of superconductivity. Now the self - energy is dependent only on frequency (and

Fig. 2 (a) Elementary act of electron - phonon scattering in the vicinity of the Fermi surface. (b) Surfaces of constant energy for initial and final states of an electron scattered by an optical phonon with energy comparable to Fermi energy. Non - trivial contribution to the average of the matrix element in (11) or (13) comes here from the intersection of these surfaces.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00718v1-figure-1e9df6e1-b81b-4ddf-9069-bbccbf84c7cc.jpg)

not on momentum) and we can use the following simple expressions, relating mass renormalization of an electron to the residue a the pole of the Green's function [16]:

$$\stackrel { \bullet } { Z ^ { - 1 } } & = 1 - \frac { \partial \Sigma ( \varepsilon ) } { \partial \varepsilon } \Big | _ { \varepsilon = 0 } & ( 6 ) & \quad \text{where} \\ m ^ { * } & = \frac { m } { \tau } = m \, \Big ( 1 - \frac { \partial \Sigma ( \varepsilon ) } { \bigcap _ { \bigcap _ { \bigcap _ { \bigcap } } } } & \quad \Big ) & ( 7 ) & \quad \text{more}$$

∣ Then from Eq. (5) by direct calculations we obtain:

$$Z ^ { - 1 } = 1 - \frac { \dots \nwarrow } { \partial \varepsilon } \Big | _ { \varepsilon = 0 } & & ( 6 ) & & \text{where} \\ m ^ { * } = \frac { m } { Z } = m \left ( 1 - \frac { \partial \Sigma ( \varepsilon ) } { \partial \varepsilon } \Big | _ { \varepsilon = 0 } \right ) & & ( 7 ) & & \text{more} \\ \text{Then from Eq. (5) by direct calculations we obtain:} \\ \partial \varepsilon / \wedge \vdash 1 & & \text{$e\infty$ $A$} \,.$$

$$& \Gamma ( \varepsilon ) \\ & - \frac { \partial \Sigma ( \varepsilon ) } { \partial \varepsilon } \Big | _ { \varepsilon = 0 } = 2 \int _ { 0 } ^ { \infty } \frac { d \omega } { \omega } \alpha ^ { 2 } ( \omega ) F ( \omega ) & ( 8 ) \quad \text{$start$} \\ & \text{and introducing the dimensionless Eliasberg - McMil-} \\ & \text{an electron - nhonon coupling constant as:} & \quad \lambda = \frac { \pi } { \pi }$$

∣ and introducing the dimensionless Eliashberg - McMillan electron - phonon coupling constant as:

$$& \dots \dots \dots \quad r \dots \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ & \lambda = 2 \int _ { 0 } ^ { \infty } \frac { d \omega } { \omega } \alpha ^ { 2 } ( \omega ) F ( \omega ) & ( 9 ) \quad \text{As in} \\ & \text{we immediately obtain the standard expression for elec-} \quad \text{$[18,20}$}$$

we immediately obtain the standard expression for electron mass renormalization due to electron - phonon interaction:

$$m ^ { * } & = m ( 1 + \lambda ) & ( 1 0 ) \quad \text{ w}$$

The function α 2 ( ω F ω ) ( ) in the expression for Eliashberg - McMillan electron - phonon coupling constant

(9) should be calculated according to (3) or (4) depending on the relation between Fermi energy E F and characteristic phonon frequency Ω 0 As long as Ω 0 /lessmuch E F we can use the standard expression (3), while in case of Ω 0 ∼ E F we should use (4).

Using Eq. (4) we can rewrite (9) in the following form:

$$& \stackrel { \dots } { \lambda } = \frac { 2 } { N ( 0 ) } \int \frac { d \omega } { \omega } \sum _ { \mathbf p } \sum _ { \mathbf p ^ { \prime } } | g _ { \mathbf p \mathbf p ^ { \prime } } | ^ { 2 } \times \\ & \times \delta ( \omega - \Omega _ { \mathbf p - \mathbf p ^ { \prime } } ) \delta ( \varepsilon _ { \mathbf p } ) \delta ( \varepsilon _ { \mathbf p ^ { \prime } } - \Omega _ { \mathbf p - \mathbf p ^ { \prime } } ) \quad \text{(11)} \\ & \text{which gives the most general version to calculate the}$$

which gives the most general expression to calculate the electron phonon constant λ , determining pairing in Eliashberg - McMillan theory. Implicitly this result was contained already in Ref. [17]. Below we shall present some simple estimates, based on this general relation.

## 3 Estimates of electron - phonon coupling with non - adiabatic phonons

Let us consider the simplest possible model of electrons interacting with a single optical (Einstein - like) phonon mode with high - enough frequency Ω 0 . The general qualitative picture of such scattering is shown in Fig. 2. In this case in Eq. (11) the density of phonon states is simply F ω ( ) = δ ω ( -Ω 0 ). Just for orientation we may take the possible momentum dependence of interaction with this optical phonon in the form proposed in Refs. [18,19] to describe nearly 'forward' scattering by optical phonons at FeSe/STO interface, as a possible mechanism of strong T c enhancement in this system:

$$g ( \mathbf q ) = g _ { 0 } \exp ( - | \mathbf q | / q _ { 0 } ),$$

where the typical value of q 0 /lessmuch p F ( p F is the Fermi momentum) to ensure the nearly 'forward' nature of scattering. This model allows explicit estimates, which may illustrate the general situation.

Now we can write the dimensionless pairing constant of electron -phonon interaction in Eliashberg theory as:

$$\mathfrak { I c } \text{Mil-} \quad & \lambda = \frac { 2 } { N ( 0 ) \Omega _ { 0 } } \sum _ { \mathbf p } \sum _ { \mathbf q } | g _ { \mathbf q } | ^ { 2 } \delta ( \varepsilon _ { \mathbf p } ) \delta ( \varepsilon _ { \mathbf p + \mathbf q } - \Omega _ { 0 } ) \quad \quad ( 1 3 ) \\ ( 9 ) \quad & \mathfrak { I c } \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-} \text{-$$

As in FeSe/STO with rather shallow conduction band [18,20,21], where in fact we have Ω 0 &gt; E F , the finite value of Ω 0 in the second δ -function here should be definitely taken into account.

We can make our estimates assuming the simplest linearized form of electronic spectrum near the Fermi surface ( v F is Fermi velocity): ε p ≈ v F ( | p | -p F ), which allows us to perform all calculations analytically. Using

(12) in (13) and considering the two - dimensional case, after the calculation of all integrals we obtain [21]:

$$& \lambda = \frac { g _ { 0 } ^ { 2 } a ^ { 2 } } { \pi ^ { 2 } v _ { F } ^ { 2 } } K _ { 1 } \left ( \frac { 2 \Omega _ { 0 } } { v _ { F } q _ { 0 } } \right ), \quad \ \ \ ( 1 4 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad & \text{where $K$.} ( r ) $ is Ressel function of \maxinarray argument $ \quad. $ \quad \quad \quad$$

where K x 1 ( ) is Bessel function of imaginary argument (McDonald function). Using the asymptotic form of K x 1 ( ) and dropping a number of irrelevant constants of the order of unity, we get:

$$\lambda \sim \lambda _ { 0 } \frac { q _ { 0 } } { 4 \pi p _ { F } }, & & ( 1 5 ) & & 4 \, \text{Am} \\ \partial \dots \, a$$

for Ω 0 v F q 0 /lessmuch 1, and

$$& \lambda \sim \lambda _ { 0 } \frac { \Omega _ { 0 } } { \pi E _ { F } } \sqrt { \frac { v _ { F } q _ { 0 } } { \Omega _ { 0 } } } \exp \left ( - \frac { 2 \Omega _ { 0 } } { v _ { F } q _ { 0 } } \right ), \quad ( 1 6 ) \quad \text{$\text{$two$}$} \\ & \text{for $\underline{$\Omega_{0}$} \ll 1$                                                                                                                                                                                                        
                                                                                                                                                                                                        --band                                                                                                                                                                                                        -band$$

for Ω 0 v F q 0 /greatermuch 1. Here we introduced the standard dimensionless electron - phonon coupling constant:

$$\lambda _ { 0 } = \frac { 2 g _ { 0 } ^ { 2 } } { \Omega _ { 0 } } N ( 0 ), \quad \ \ & \quad ( 1 7 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where N (0) is the density of electronic states at the Fermi level per single spin projection.

The result (15) is by itself rather unfavorable for significant T c enhancement in model under discussion, where q 0 /lessmuch p F . Even worse is the situation if we take into account the large values of Ω 0 , as pairing constant becomes exponentially suppressed for Ω 0 v F q 0 &gt; 1, which is typical for FeSe/STO interface, where Ω 0 &gt; E F /greatermuch v F q 0 [9]. This makes the enhancement of T c due to interaction of FeSe electrons with optical phonons of STO rather improbable, as was stressed in Ref. [21].

However, this is not our main point here. Actually, using (12) we can also make estimates for generally more typical case, when the optical phonon scatters electrons not only in nearly 'forward' direction, but in a wider interval of transferred momenta. To do that we have simply to use in Eq. (12)the larger values of parameter q 0 . Choosing e.g. q 0 ∼ 4 πp F and using the low frequency limit of (15) we immediately obtain λ ≈ λ 0 , i.e. the standard result. Similarly, parameter q 0 can be taken of the order of inverse lattice vector 2 π/a (where a is the lattice constant). Then for q 0 ∼ 2 π/a from (15) we obtain:

$$\stackrel { \cdots \cdots \cdots \cdots \cdots } { \lambda \sim \lambda _ { 0 } \frac { 1 } { 2 p _ { F } a } } \sim \lambda _ { 0 } & & ( 1 8 ) & & \\ \stackrel { \cdots \cdots } { \nu \neq \nu }$$

for the typical case of p F ∼ 1 / 2 a . In general there always remains the dependence on the value of Fermi momentum and cutoff parameter (cf. similar analysis in Ref. [16]). These particular estimates are valid for the adiabatic case.

In antiadiabatic limit of (16), assuming q 0 ∼ p F we immediately obtain:

$$& \text{immediately obtain: } & \stackrel { \sim \dots } { \lambda } \sim \frac { \sqrt { 2 } } { \pi } \lambda _ { 0 } \sqrt { \frac { \Omega _ { 0 } } { E _ { F } } } \exp \left ( - \frac { \Omega _ { 0 } } { E _ { F } } \right ), & ( 1 9 ) & - \frac { \partial \Sigma } { \partial }$$

which simply signifies the effective interaction cutoff for Ω 0 &gt; E F in the antiadiabatic limit. This fact was already noted by Gor'kov in Refs. [12,13], where it was stressed that in antiadiabatic limit the cutoff in the Cooper channel is determined not by the average phonon frequency, but by Fermi energy.

## 4 Antiadiabatic limit and mass renormalization

Our discussion up to now implicitly assumed the conduction band of an infinite width. However, it is obvious that in case of large enough characteristic phonon frequency it may become comparable with conduction band width, which in typical metal case is of the order of Fermi energy E F . Now we will show that in the strongly nonadiabatic (antiadiabatic) limit, when Ω 0 /greatermuch E F ∼ D (here D is the conduction band halfwidth), we are in fact dealing with the situation, when there appears a new small parameter of perturbation theory λD/Ω 0 ∼ λE F /Ω 0 .

Consider the case of conduction band of the finite width 2 D with constant density of states (which formally corresponds to two - dimensional case). The Fermi level as always is considered as an origin of energy scale and for simplicity we assume the case of half - filled band. Then (5) reduces to:

$$E _ { F } \gg \text{ \emph{-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-S-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-b-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-n-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-a-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-ss-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-ns-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-r-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s- s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-t-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-g-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-l-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-sl-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-x-n-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s- s- s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s- s-ss-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-s-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-s-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-SS-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-ss-$$

For the model of a single optical phonon F ω ( ) = δ ω ( -Ω 0 ) and we immediately obtain:

$$\Gamma _ { \L } \nu \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ ( 1 8 ) \\ \Gamma _ { \L } \nu & \Sigma ( \varepsilon ) = \alpha ^ { 2 } ( \Omega _ { 0 } ) F ( \Omega _ { 0 } ) & \left \{ \ln \frac { \varepsilon + D + \Omega _ { 0 } - i \delta } { \varepsilon - D - \Omega _ { 0 } + i \delta } \\ \Gamma _ { \L } \nu & \quad - \ln \frac { \varepsilon + \Omega _ { 0 } - i \delta } { \varepsilon - \Omega _ { 0 } + i \delta } \right \} \\ \Gamma _ { \L } \nu & \quad \sim \quad.. \quad \nu \quad \widehat { \lambda } \quad \widehat { \lambda } \quad.$$

Correspondingly, from (20) we get:

$$( 1 9 ) \quad - \frac { \partial \Sigma ( \varepsilon ) } { \partial \varepsilon } \Big | _ { \varepsilon = 0 } = 2 \int _ { 0 } ^ { \infty } d \omega \alpha ^ { 2 } ( \omega ) F ( \omega ) \frac { D } { \omega ( \omega + D ) } \quad ( 2 2 )$$

and we can define the new generalized coupling constant as:

$$\tilde { \lambda } & = 2 \int _ { 0 } ^ { \infty } \frac { d \omega } { \omega } \alpha ^ { 2 } ( \omega ) F ( \omega ) \frac { D } { \omega + D } \quad \quad ( 2 3 ) \quad \text{All an} \\ & \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \$$

which for D → ∞ reduces to the usual Eliashberg McMillan constant (9), while for D → 0 ( D /lessmuch Ω 0 ) it gives the 'antiadiabatic' coupling constant:

$$\lambda _ { D } = 2 D \int \frac { d \omega } { \omega ^ { 2 } } \alpha ^ { 2 } ( \omega ) F ( \omega ) \quad \ \ \ ( 2 4 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
- \quad \quad \quad$$

Eq. (23) describes the smooth transition between the limits of wide and narrow conduction bands. Mass renormalization in general case is determined by ˜ : λ

$$m ^ { * } & = m ( 1 + \tilde { \lambda } ) & ( 2 5 ) & \quad \overset { \mu \sigma \dots } { \text{lating} }$$

For the model of a single optical phonon with frequency Ω 0 we have:

$$\tilde { \lambda } = \frac { 2 } { \Omega _ { 0 } } \alpha ^ { 2 } ( \Omega _ { 0 } ) \frac { D } { \Omega _ { 0 } + D } = \lambda \frac { D } { \Omega _ { 0 } + D } = \lambda _ { D } \frac { \Omega _ { 0 } } { \Omega _ { 0 } + D } \ ( 2 6 ) \quad G ( \varepsilon, \mathfrak { p }$$

where Eliashberg - McMillan constant is:

$$\lambda & = 2 \int _ { 0 } ^ { \infty } \frac { d \omega } { \omega } \alpha ^ { 2 } ( \omega ) F ( \omega ) = \alpha ^ { 2 } ( \Omega _ { 0 } ) \frac { 2 } { \Omega _ { 0 } } \quad \quad ( 2 7 ) \quad \text{where} \\ & \cdot \cdot \cdot \cdot$$

and λ D reduces to:

$$\lambda _ { D } = 2 \alpha ^ { 2 } ( \Omega _ { 0 } ) \frac { D } { \Omega _ { 0 } ^ { 2 } } = 2 \alpha ^ { 2 } ( \Omega _ { 0 } ) \frac { 1 } { \Omega _ { 0 } } \frac { D } { \Omega _ { 0 } } = \lambda \frac { D } { \Omega _ { 0 } } \quad \ \ ( 2 8 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where in the last expression we explicitly introduced the new small parameter D/Ω 0 /lessmuch 1, appearing in strong antiadiabatic limit. Correspondingly, in this limit we always have:

$$\lambda _ { D } = \lambda \frac { D } { \Omega _ { 0 } } \sim \lambda \frac { E _ { F } } { \Omega _ { 0 } } \ll \lambda$$

so that for reasonable values of λ (even up to a strong coupling region of λ ∼ 1) 'antiadiabatic' coupling constant remains small. Obviously, all vertex corrections here are also small, as was shown rather long ago by direct calculations in Ref. [22]. Thus we come to an unexpected conclusion - in the limit of strong nonadiabaticity the electron - phonon coupling becomes weak and we obtain a kind of 'anti - Migdal' theorem.

Physically, the weakness of electron - phonon coupling in strong nonadiabatic limit is more or less clear - when ions move much faster than electrons, these rapid oscillation are just averaged in time as electrons can not follow the very rapidly changing configuration of ions.

## 5 Eliashberg equations and the temperature of superconducting transition

All analysis above was performed for the normal state of a metal. Now let us turn to the superconducting phase. The problem arises, to what extent the results obtained can be generalized for the case of a metal in superconducting state? In particular, what coupling constant ( λ or ˜ ) λ determines the temperature of superconducting transition T c in antiadiabatic limit? Let us analyze this situation within appropriate generalization of Eliashberg equations.

Taking into account that in antiadiabatic approximation vertex corrections are are again irrelevant and neglecting the direct Coulomb repulsion, Eliashberg equations can be derived in the usual way by calculating the diagram of Fig. 1, where electronic Green's function in superconducting state is taken in Nambu's matrix representation. For real frequencies this Green's function is written in the following standard form [3,4]:

$$\mathfrak { i } \ ( 2 6 ) \quad G ( \varepsilon, \mathbf p ) = \frac { Z ( \varepsilon ) \varepsilon \tau _ { 0 } + \varepsilon _ { \mathbf p } \tau _ { 3 } + Z ( \varepsilon ) \Delta ( \varepsilon ) \tau _ { 1 } } { Z ^ { 2 } ( \varepsilon ) \varepsilon ^ { 2 } - Z ^ { 2 } ( \varepsilon ) \Delta ^ { 2 } ( \varepsilon ) - \varepsilon _ { \mathbf p } ^ { 2 } } \quad ( 3 0 ) \\ \quad \text{which a common code to the matrix of cell } \quad \text{anonym}.$$

which corresponds to the matrix of self - energy:

$$\Sigma ( \varepsilon ) = [ 1 - Z ( \varepsilon ) ] \varepsilon \tau _ { 0 } + Z ( \varepsilon ) \Delta ( \varepsilon ) \tau _ { 1 } \quad \quad ( 3 1 ) \\ ( 2 7 ) \quad. \quad. \quad.. \quad.. \quad.. \quad.. \quad.. \quad..$$

where τ i are standard Pauli matrices, while functions of mass renormalization Z ε ( ) and energy gap ∆ ε ( ) are determined from solution of integral Eliashberg equations [3,4]. For us now it is sufficient to consider only the linearized Eliashberg equations, determining superconducting transition temperature T c , which for the case of real frequencies are written as [3,4]:

$$\ e u n e \quad \infty \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \quad \ e u n e \\ \quad ( 2 9 ) \quad \ e ^ { D } \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{array}$$

$$& ( 2 9 ) & \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots.$$

In difference with the standard approach [4], we have introduced the finite integration limits, determined by the (half)bandwidth D . To simplify the analysis we again assume the half-filled band of degenerate electrons in two dimensions, so that D = E F /greatermuch T c , with constant density of states.

Situation is considerably simplified [14,15], if we consider these equations in the limit of ε → 0 and look

for the solutions 1 Z (0) = Z and ∆ (0) = ∆ . Then from (32) we obtain:

$$& [ 1 - Z ] \varepsilon = - 2 \varepsilon \int _ { 0 } ^ { \infty } d \omega \alpha ^ { 2 } ( \omega ) F ( \omega ) \frac { D } { \omega ( \omega + D ) } \quad ( 3 4 ) \quad 1 + \tilde { \lambda } \\ & \text{and we get the mass renormalization factor as} \colon$$

and we get the mass renormalization factor as:

$$Z = 1 + \tilde { \lambda }$$

where constant ˜ was defined above in Eq. (23), which λ for D → ∞ reduces to the usual Eliasberg - McMillan constant (9), while for D significantly smaller than characteristic phonon frequencies it gives the 'antiadiabatic' coupling constant (24). Mass renormalization is again determined by this generalized coupling constant ˜ λ as in Eq. (25). In particular, in the strong antiadiabatic limit this renormalization is quite small and determined by the limiting expression λ D given by Eq. (24).

Situation is quite different in Eq. (33). In the limit of ε → 0, using (35) we immediately obtain from (33) the following equation for T c :

$$1 + \tilde { \lambda } & = 2 \int _ { 0 } ^ { \mu } d \omega \alpha ^ { 2 } ( \omega ) F ( \omega ) \int _ { 0 } ^ { D } \frac { d \varepsilon ^ { \prime } } { \varepsilon ^ { \prime } ( \varepsilon ^ { \prime } + \omega ) } t h \frac { \varepsilon ^ { \prime } } { 2 T _ { c } } \ ( 3 6 ) \quad \text{Now} \ \, \\ & \text{where} \ \lambda \text{ is the standard Fliashher\sigma} - \text{McMillan counlin\sigma} } \quad \text{freque}$$

where λ is the standard Eliashberg - McMillan coupling constant as defined above in Eq. (9). Thus, in general case, different coupling constants determine mass renormalization and T c .

Let us consider rather general model with discrete set of dispersionless phonon modes (Einstein phonons). In this case the phonon density of states is written as:

$$F ( \omega ) & = \sum _ { i } \delta ( \omega - \Omega _ { i } ) & ( 3 7 ) & \quad \text{interp} \\ & \text{where } \Omega \colon \text{are discrete frequencies modeling the optical} & \quad \tau \ \dots$$

where Ω i are discrete frequencies modeling the optical branches of the phonon spectrum. Then from Eqs. (9) and (23) we get:

$$\lambda = 2 & \sum _ { i } \frac { \alpha ^ { 2 } ( \Omega _ { i } ) } { \Omega _ { i } } \equiv \sum _ { i } \lambda _ { i } \\ \cdots & \quad - \ \wedge ^ { 2 } ( O \cdot \, N \quad - \quad n \quad - \cdots \quad \ln ( \Omega )$$

$$& \tilde { \lambda } = 2 \sum _ { i } \frac { \alpha ^ { 2 } ( \Omega _ { i } ) D } { \Omega _ { i } ( \Omega _ { i } + D ) } = \sum _ { i } \lambda _ { i } \frac { D } { \Omega _ { i } + D } \equiv \sum _ { i } \tilde { \lambda } _ { i } \quad ( 3 9 ) \quad \text{$\ln \langle \Omega \rangle$} \\ & C o r r e s n o n d i n o v \text{ in this case} \cdot \quad \text{$\Gamma$} \text{man} \quad.$$

Correspondingly, in this case:

$$\text{expresspanny, in us case.} & & \text{quenc} \\ \alpha ^ { 2 } ( \omega ) F ( \omega ) = & \sum _ { i } \alpha ^ { 2 } ( \Omega _ { i } ) \delta ( \omega - \Omega _ { i } ) & & \ln ( \Omega ) \\ & = \sum _ { i } \frac { \lambda _ { i } } { 2 } \Omega _ { i } \delta ( \omega - \Omega _ { i } ) & & ( 4 0 ) & & \text{where} \\ \text{The standard Eliashberg equation (in adiabatic limit)} & & \text{mic f}$$

The standard Eliashberg equation (in adiabatic limit) for such model were consistently solved in Ref. [23]. For

1 To avoid confusion note, that according to standard notations of Eliashberg - McMillan theory the renormalization factor Z as defined here is just the inverse of a similar factor defined in Eq. (6) for the normal state

our purposes it is sufficient to analyze only Eq. (36), which takes now the following form:

$$( 3 4 ) \quad & 1 + \tilde { \lambda } = 2 \sum _ { i } \alpha ^ { 2 } ( \Omega _ { i } ) \int _ { 0 } ^ { D } \frac { d \varepsilon ^ { \prime } } { \varepsilon ^ { \prime } ( \varepsilon ^ { \prime } + \Omega _ { i } ) } t h \frac { \varepsilon ^ { \prime } } { 2 T _ { c } } \quad \quad ( 4 1 ) \\ & \text{This equation is easily solved to obtain:}$$

This equation is easily solved to obtain:

$$\text{which} _ { \text{IcMil-r than} } \quad & T _ { c } \sim \prod _ { i } \left ( \frac { D } { 1 + \frac { D } { \Omega _ { i } } } \right ) ^ { \frac { \lambda _ { i } } { \lambda } } \exp \left ( - \frac { 1 + \tilde { \lambda } } { \lambda } \right ) \quad & ( 4 2 ) \\ \text{$r$ than} \quad & \text{In the simple case of two optical phonons with frequen-}$$

In the simple case of two optical phonons with frequencies Ω 1 and Ω 2 we have:

$$\begin{smallmatrix} \text{tion is} \\ \text{istant} \\ \text{iadia} \\ \text{id de} \end{smallmatrix} \quad T _ { c } \sim \left ( \frac { D } { 1 + \frac { D } { \mathfrak { M } _ { 1 } } } \right ) ^ { \frac { \lambda _ { 1 } } { \lambda } } \left ( \frac { D } { 1 + \frac { D } { \mathfrak { M } _ { 2 } } } \right ) ^ { \frac { \lambda _ { 2 } } { \lambda } } \exp \left ( - \frac { 1 + \tilde { \lambda } } { \lambda } \right ) \quad ( 4 3 )$$

where ˜ λ = ˜ λ 1 + ˜ λ 2 and λ = λ 1 + λ 2 . For the case of Ω 1 /lessmuch D (adiabatic phonon), and Ω 2 /greatermuch D (antiadiabatic phonon) Eq. (43) is immediately reduced to:

$$\intertext { n ( 3 3 ) } _ { \intertext { \nu c } } T _ { c } \sim ( \Omega _ { 1 } ) ^ { \frac { \lambda _ { 1 } } { \lambda } } ( D ) ^ { \frac { \lambda _ { 2 } } { \lambda } } \exp \left ( - \frac { 1 + \tilde { \lambda } } { \lambda } \right )$$

Now we can see, that in the preexponential factor the frequency of antiadiabatic phonon is replaced by band half-width (Fermi energy), which plays a role of the cutoff for logarithmic divergence in Cooper channel in antiadiabatic limit [14,12,13].

Our general result (42) gives the general expression for T c for the model with discrete set of optical phonons, valid both in adiabatic and antiadiabatic regimes and interpolating between these limits in intermediate region. Actually Eq. (42) can be easily rewritten as:

$$\text{ptical} _ { \mathbb { S }. \ ( 9 ) } \quad T _ { c } \sim \langle \Omega \rangle \exp \left ( - \frac { 1 + \tilde { \lambda } } { \lambda } \right )$$

where we have introduced the average logarithmic frequency 〈 Ω 〉 as:

In the limit of continuous distribution of phonon frequencies this last expression reduces to:

$$\ln \langle \Omega \rangle = \ln \prod _ { i } \left ( \frac { D } { 1 + \frac { D } { \Omega _ { i } } } \right ) ^ { \frac { \lambda _ { i } } { \lambda } } = \sum _ { i } \frac { \lambda _ { i } } { \lambda } \ln \frac { D } { 1 + \frac { D } { \Omega _ { i } } } \quad ( 4 6 ) \\ \ln \text{ the limit of continuous distribution of phonon fre-}$$

$$& \ln \langle \Omega \rangle = \frac { 2 } { \lambda } \int \frac { d \omega } { \omega } \alpha ^ { 2 } ( \omega ) F ( \omega ) \ln \frac { D } { 1 + \frac { D } { \omega } } \quad \quad ( 4 7 ) \\ & \omega \in \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon \infty \infty$$

where λ is given by the usual expression (9). Eq. (47) generalizes the standard definition of average logarithmic frequency of Eliasberg - McMillan theory [4] for the case of finite bandwidth. Obviously, it reduces to the standard expression in adiabatic limit of phonon frequencies much lower than D , and gives 〈 Ω 〉 ∼ D in extreme antiadiabatic limit, when all phonon frequencies are much larger than D .

## 6 Coulomb pseudopotential

Up to now we have neglected the direct Coulomb repulsion of electrons, which in the standard approach [1,2,3, 4,5] is described by Coulomb pseudopotential µ /star , which is effectively suppressed by large Tolmachev's logarithm. As we noted in Ref. [14] antiadiabatic phonons actually suppress Tolmachev's logarithm, which can probably lead to rather strong suppression of the temperature of superconducting transition. To clarify this situation we consider the simplified version of integral equation for the gap (33), writing it in the standard form:

$$Z ( \varepsilon ) \Delta ( \varepsilon ) = \int _ { 0 } ^ { D } d \varepsilon ^ { \prime } K ( \varepsilon, \varepsilon ^ { \prime } ) \frac { 1 } { \varepsilon ^ { \prime } } t h \frac { \varepsilon ^ { \prime } } { 2 T _ { c } } \Delta ( \varepsilon ^ { \prime } ) \quad \quad ( 4 8 ) \quad \text{Now 1} \\ \text{where to internal normal is a combination of two action}$$

where the integral kernel is a combination of two step - functions:

$$K ( \varepsilon, \varepsilon ^ { \prime } ) & = \lambda \theta ( \langle \Omega \rangle - | \varepsilon | ) \theta ( \langle \Omega \rangle - | \varepsilon ^ { \prime } | ) & \quad \text{$stants$} \\ - \mu \theta ( D - | \varepsilon | ) \theta ( D - | \varepsilon ^ { \prime } | ) & \quad ( 4 9 ) & \quad \text{where} \\ \dots \dots \dots \vdots \vdots \vdots \dots \vdots \dots \vdots \varepsilon \dots \vdots \dots \vdots \dots \end{cases} \, \frac { \text{$stants$} } { \text{$\dots$} \dots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots \dots \vdots }$$

where µ is the dimensionless (repulsive) Coulomb potential, while the parameter 〈 Ω 〉 , determining the energy width of attraction region due to phonons is determined by preexponential factor (average logarithmic frequency) of Eqs. (42),(45).

$$\langle \Omega \rangle = \prod _ { i } \left ( \frac { D } { 1 + \frac { D } { \Omega _ { i } } } \right ) ^ { \frac { \lambda _ { i } } { \lambda _ { i } } } & & \mu ^ { * } = \cdot \\ \text{It is important that we always have } \langle \Omega \rangle < D \text{ } F o ( 4 8 ) & & \text{We ci}$$

It is important that we always have 〈 Ω &lt; D 〉 . Eq. (48) is now rewritten as:

$$\text{is now even as.} \\ Z ( \varepsilon ) \Delta ( \varepsilon ) = ( \lambda - \mu ) \int _ { 0 } ^ { \langle \Omega \rangle } \frac { d \varepsilon ^ { \prime } } { \varepsilon ^ { \prime } } t h \frac { \varepsilon ^ { \prime } } { 2 T _ { c } } \Delta ( \varepsilon ^ { \prime } ) & \text{the ra} \\ - \mu \int _ { \langle \Omega \rangle } ^ { D } \frac { d \varepsilon ^ { \prime } } { \varepsilon ^ { \prime } } \Delta ( \varepsilon ^ { \prime } ) & \text{the 51} \text{$persis} \\ & \text{the p\delta$} \\ \text{Writing the mass renormalization due to phonons as} \colon & \text{phone}$$

Writing the mass renormalization due to phonons as:

$$Z ( \varepsilon ) & = \begin{cases} 1 + \tilde { \lambda } \quad \text{for} \quad \varepsilon < \langle \Omega \rangle \\ 1 \quad \text{for} \quad \varepsilon > \langle \Omega \rangle \end{cases} \quad \begin{pmatrix} & \text{real} \ c \\ ( 5 2 ) & \text{cores} \\ & \text{by coin} \\ & \text{bastic} \end{pmatrix}$$

we look for the solution of Eq. (48) for ∆ ε ( ), as usual, in the following form [2,4,5]:

$$\Delta ( \varepsilon ) & = \begin{cases} \Delta _ { 1 } & \text{for} \quad \varepsilon < \langle \Omega \rangle \\ \Delta _ { 2 } & \text{for} \quad \varepsilon > \langle \Omega \rangle \end{cases} & \text{(53)} & \text{consid} \\ \tau \quad \psi \quad \psi \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

Then Eq. (51) is transformed into the system of two homogeneous linear equations for constants ∆ 1 and ∆ 2 :

$$\stackrel { \dots } { \Delta } _ { 1 } = ( \lambda - \mu ) \ln \frac { \langle \Omega \rangle } { T _ { c } } \Delta _ { 1 } - \mu \ln \frac { D } { \langle \Omega \rangle } \Delta _ { 2 } \quad \begin{pmatrix} \ln \, \text{pr} \\ \text{phone} \\ \Delta _ { 2 } = - \mu \ln \frac { \langle \Omega \rangle } { T _ { c } } \Delta _ { 1 } - \mu \ln \frac { D } { \langle \Omega \rangle } \Delta _ { 2 } \quad \begin{pmatrix} 5 4 \end{pmatrix} \text{tag} \\ \text{enoug}$$

The condition of the existence of nontrivial solution here is:

$$\text{repul-} \\ \text{which} \quad \text{$1+\tilde{\lambda}=\left(\lambda-\frac{\mu}{1+\mu\ln\frac{D}{\langle\Omega\rangle}} \right) \ln\frac{\langle\Omega\rangle} T_{c}$}$$

Then the transition temperature is given by:

$$\underset { \overset { \text{max} } { \text{the can} } } { \text{the team} } \quad T _ { c } = \langle \Omega \rangle \exp \left ( - \frac { 1 + \tilde { \lambda } } { \lambda - \mu ^ { * } } \right ) \quad \quad \quad \ \ ( 5 6 )$$

where the Coulomb pseudopotential is determined as:

$$\ n d a r d \begin{array} { c } \mu ^ { * } = \frac { \mu } { 1 + \mu \ln \frac { D } { \langle \Omega \rangle } } = \frac { \mu } { 1 + \mu \ln \prod _ { i } \left ( 1 + \frac { D } { \Omega _ { i } } \right ) ^ { \frac { \lambda _ { i } } { \lambda } } } \\ \quad \text{Now the phonon frequencies enter Tolmachev's logo-} \\ \quad \text{with} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$} \quad \text{$the$me$}$$

Now the phonon frequencies enter Tolmachev's logarithm as the product of partial contributions, with its values determined also by corresponding coupling constants. Similar structure of Tolmachev's logarithm was first obtained (in somehow different model) in Ref. [24], where the case of frequencies going outside the limits of adiabatic approximation was not considered. In this sense, Eq. (57) has a wider region of applicability. In particular, for the model of two optical phonons with frequencies Ω 1 /lessmuch D (adiabatic phonon) and Ω 2 /greatermuch D , from Eq. (57) we get:

$$& \dots \, - \, \psi \, \psi \, \cdots \, \circ \cdots \\ & \mu ^ { * } = \frac { \mu } { 1 + \mu \ln \left ( \frac { D } { \Omega _ { 1 } } \right ) ^ { \frac { \lambda _ { 1 } } { \lambda } } } = \frac { \mu } { 1 + \mu \frac { \lambda _ { 1 } } { \lambda } \ln \frac { D } { \Omega _ { 1 } } } \\ & \quad \text{We can see, that the contribution of antiadiabatic}$$

We can see, that the contribution of antiadiabatic phonon drops out of Tolmachev's logarithm, while the logarithm itself persists, with its value determined by the ratio of the band halfwidth (Fermi energy) to the frequency of adiabatic (low frequency) phonon. The general effect of suppression of Coulomb repulsion also persists, though it becomes somehow weaker due to the partial interaction of electrons with corresponding phonon. This situation is conserved also in the general case - the value of Tolmachev's logarithm and corresponding Coulomb pseudopotential is determined by contributions of adiabatic phonons, while antiadiabatic phonons drop out. Thus, in general case, situation becomes more favorable for superconductivity, as compared to the case of a single antiadiabatic phonon, considered in Ref. [14].

## 7 Conclusions

In present paper we have considered the electron -phonon coupling in Eliashberg -McMillan theory, taking into account antiadiabatic phonons with high enough frequency (comparable or exceeding the Fermi

energy E F ). The value of mass renormalization, in general case, was shown to be determined by the new coupling constant ˜ , λ while the value of the pairing interaction is always determined by the standard coupling constant λ of Eliashberg - McMillan theory, appropriately generalized by taking into account the finite value of phonon frequency [14]. Mass renormalization due to strongly antiadiabatic phonons is in general small and determined by the coupling constant λ D /lessmuch λ . In this sense, in the limit of strong antiadiabaticity, the coupling of such phonons with electrons becomes weak and corresponding vertex correction again become irrelevant [22,14], creating a kind of 'anti - Migdal' situation. This fact allows us to use Eliashberg - McMillan approach in the limit of strong antiadiabaticity. In the intermediate region all our expressions just produce a smooth interpolation between adiabatic and antiadiabatic limits.

The cutoff of pairing interaction in Cooper channel in antiadiabatic limit becomes effective at energies ∼ E F ∼ D , as was previously noted in Refs. [12,13, 14]), so that corresponding phonons do not contribute to Tolmachev's logarithm in Coulomb pseudopotential. However, the large enough values of this logarithm (and corresponding smallness of µ /star ) can be guaranteed due to contributions from adiabatic phonons [15].

Note that above we have used rather simplified analysis of Eliashberg equations. However, in our opinion, more elaborate approach, e.g. along the lines of Ref. [23], will not lead to qualitative change of our results. Some simple estimates for FeSe/STO system, based on these results, can be found in Refs. [14,15].

Acknowledgements It was a pleasure and honor to contribute this work to Ted Geballe's Festschrift, as his work in the field of superconductivity has influenced all of us for decades. This work was partially supported by RFBR grant No. 17-02-00015 and the program of fundamental research No. 12 of the RAS Presidium 'Fundamental problems of high temperature superconductivity'.

## References

- 1. J.R. Schrieffer. Theory of Superconductivity, WA Benjamin, NY, 1964
- 2. D.J. Scalapino. In 'Superconductivity', p. 449, Ed. by R.D. Parks, Marcel Dekker, NY, 1969
- 3. R.M. White, T.H. Geballe. Long Range Order in Solids. Academic Press, New York, San Francisco, London, 1979
- 4. S.V. Vonsovsky, Yu.A. Izyumov, E.Z. Kurmaev. Superconductivity of Transition metals, Their Alloys and Compounds, Springer, Berlin - Heidelberg, 1982
- 5. P.B. Allen, B. Mitrovi´. Solid State Physics, Vol. Vol. 37 c (Eds. F. Seitz, D. Turnbull, H. Ehrenreich), Academic Press, NY, 1982, p. 1
- 6. A.B. Migdal. Zh. Eksp. Teor. Fiz. 34 , 1438 (1958) [Sov. Phys. JETP 7 , 996 (1958)]
- 7. A.S. Alexandrov. A.B. Krebs. Usp. Fiz. Nauk 162 , 1 (1992) [Physics Uspekhi 35 , 345 (1992)]
- 8. I. Esterlis, B. Nosarzewski, E.W. Huang, D. Moritz, T.P. Devereux, D.J. Scalapino, S.A. Kivelson. Phys. Rev. B 97 , 140501(R) (2018)
- 9. M.V. Sadovskii. Usp. Fiz. Nauk 178 , 1243 (2008) [Physics Uspekhi 51 , 1243 (2008)]
- 10. L.P. Gor'kov, V.Z. Kresin. Rev. Mod. Phys. 90 , 011001 (2018)
- 11. L.P. Gor'kov. PNAS 113 , 4646 (2018)
- 12. L.P. Gor'kov. Phys. Rev. B 93 , 054517 (2016)
- 13. L.P. Gor'kov. Phys. Rev. B 93 , 060507 (2016)
- 14. M.V. Sadovskii.Zh. Eksp. Teor. Fiz 155 , 527 (2019) [JETP 128 , 455 (2019)]
- 15. M.V. Sadovskii. Pis'ma Zh. Eksp. Teor. Fiz. 109 , 165 (2010) [JETP Letters 109 , 166 (2019)]
- 16. M.V. Sadovskii. Diagrammatics. World Scientific, Singapore, 2006
- 17. P.B. Allen. Phys. Rev. B 6 , 2577 (1972)
- 18. J.J. Lee, F.T. Schmitt, R.G. Moore, S. Johnston, Y.T. Cui,W. Li, Z.K. Liu, M. Hashimoto, Y. Zhang, D.H. Lu, T.P. Devereaux, D.H. Lee, Z.X. Shen. Nature 515 , 245 (2014)
- 19. M.L. Kuli´, O.V. Dolgov. New J. Phys. c 19 , 013020 (2017)
- 20. I.A. Nekrasov, N.S. Pavlov, V.V. Sadovskii. Pis'ma Zh. Eksp. Teor. Fiz. 105 , 354 (2017) [JETP Letters 105 , 370 (2017)]
- 21. I.A. Nekrasov, N.S. Pavlov, M.V. Sadovskii. Zh. Eksp. Teor. Fiz. 153 , 590 (2018) [JETP 126 , 485 (2018)] Phys. Rev. B 93 , 134513 (2016)
- 22. M.A. Ikeda, A. Ogasawara, M. Sugihara. Phys. Lett. A 170 , 319 (1992)
- 23. A.E. Karakozov, E.G. Maksimov, S.A. Mashkov. Zh. Eksp. Teor. Fiz. 68 , 1937 (1975) [JETP 41 , 971 (1975)]
- 24. D.A. Kirzhnits, E.G. Maksimov, D.I. Khomskii. J. Low. Temp. Phys. 10 , 79 (1973)
M. V. Sadovskii
Antiadiabatic Phonons and Superconductivity in Eliashberg-McMillan Theory
The standard Eliashberg - McMillan theory of superconductivity is essentially
based on the adiabatic approximation. Here we present some simple estimates of
electron - phonon interaction within Eliashberg - McMillan approach in non -
adiabatic and even antiadiabatic situation, when characteristic phonon
frequency $\Omega_0$ becomes large enough, i.e. comparable or exceeding the
Fermi energy $E_F$. We discuss the general definition of Eliashberg - McMillan
(pairing) electron - phonon coupling constant $\lambda$, taking into account
the finite value of phonon frequencies. We show that the mass renormalization
of electrons is in general determined by different coupling constant
$\tilde\lambda$, which takes into account the finite width of conduction band,
and describes the smooth transition from the adiabatic regime to the region of
strong nonadiabaticity. In antiadiabatic limit, when $\Omega_0\gg E_F$, the new
small parameter of perturbation theory is
$\lambda\frac{E_F}{\Omega_0}\sim\lambda\frac{D}{\Omega_0}\ll 1$ ($D$ is
conduction band half -- width), and corrections to electronic spectrum (mass
renormalization) become irrelevant. However, the temperature of superconducting
transition $T_c$ in antiadiabatic limit is still determined by Eliashberg -
McMillan coupling constant $\lambda$. We consider in detail the model with
discrete set of (optical) phonon frequencies. A general expression for
superconducting transition temperature $T_c$ is derived, which is valid in
situation, when one (or several) of such phonons becomes antiadiabatic. We also
analyze the contribution of such phonons into the Coulomb pseudopotential
$\mu^{\star}$ and show, that antiadiabatic phonons do not contribute to
Tolmachev's logarithm and its value is determined by partial contributions from
adiabatic phonons only.