Antiadiabatic Theory of Superconducting State Transition: Phonons and Strong Electron Correlations—The Old Physics and New Aspects

Complex electronic ground state of molecular and solid state system is analyzed on the ab initio level beyond the adiabatic BornOppenheimer approximation (BOA). The attention is focused on the band structure fluctuation (BSF) at Fermi level, which is induced by electron-phonon coupling in superconductors, and which is absent in the non-superconducting analogues. The BSF in superconductors results in breakdown of the adiabatic BOA. At these circumstances, chemical potential is substantially reduced and system is stabilized (effect of nuclear dynamics) in the antiadiabatic state at broken symmetry with a gap(s) in one-particle spectrum. Distorted nuclear structure has fluxional character and geometric degeneracy of the antiadiabatic ground state enables formation of mobile bipolarons in real space. It has been shown that an effective attractive e-e interaction (Cooper-pair formation) is in fact correction to electron correlation energy at transition from adiabatic into antiadiabatic ground electronic state. In this respect, Cooper-pair formation is not the primary reason for transition into superconducting state, but it is a consequence of antiadiabatic state formation. It has been shown that thermodynamic properties of system in antiadiabatic state correspond to thermodynamics of superconducting state. Illustrative application of the theory for different types of superconductors is presented.


Introduction
Superconductivity, an amazing physical phenomenon discovered nearly 100 years ago by Kamerlingh Onnes [1] and his assistant Gilles Holst whose name has been basically forgotten by history (see, e.g., [2]), has been one of the most important research fields of solid-state physics of last century and it remains until the present days.What is extremely irritating is the fact that microscopic mechanism of superconducting state transition, in spite of enormous attention which has been paid to this effect, remains still unclear and represents an open challenge for theory.
Until the discovery of high-temperature superconductivity of cuprates by Bednorz and Muller in 1986 [3] and synthesis of first 90 K superconductor [4] in 1987, understanding of microscopic mechanism of superconducting (SC) state transition formulated within the BCS theory in 1957 [5] was generally accepted as a firm theoretical basis behind the physics of this phenomenon.Here, the basic physics is the idea of Cooper-pairs formation, that is, formation of boson-like particles in momentum space, which are stable in thin layer above the Fermi level and drive system into more stable-superconducting state.Sufficient condition of pair formation is whatever weak, but attractive, interaction between electrons.Real possibility of effective attractive electron-electron interactions was well known.Some years, ago it was derived by Fröhlich [6,7] as a consequence of electron-phonon (e-p) interactions.
The range of validity of the BCS theory with respect to e-p interactions has been specified by Migdal [8] and Eliashberg [9,10].It can be interpreted as Migdal's theorem and Eliashberg's restriction (ME approximation).The first is related to validity of the condition ωλ/E F 1 and the second one restricts the validity only for λ ≤ 1, where λ is e-p coupling strength and ω and E F are characteristic phonon and electron energy scales, respectively.Expressed explicitly, BCS-like theories are valid only for adiabatic systems that obey the adiabatic Born-Oppenheimer approximation (BOA): ω/E F 1. While for conventional (low-temperature) superconductors, the BCS theory within the ME approximation (i.e., weak coupling regime) is an excellent extension of standard theory of metals, for high-temperature cuprates in order to interpret high critical temperature and ensure the pairs condensation, beside (or instead of) the e-p interactions the important role of other interaction mechanisms has been advocated (see e.g., [11][12][13]).Since the copper, a transition metal with incompletely filled d-shell when chemically bounded, is a central atom of hight-T c cuprates, it is quite natural that the attention has been focused on strong electron correlations (in a sense of standard Coulomb-repulsive e-e interactions), magnetic interactions, and/or spin fluctuation effects.At the present, the effect of Coulomb repulsion is usually incorporated via Coulomb pseudopotential μ * and critical temperature is calculated according to McMillan formula [14].
The e-p interactions, which have been accepted to be responsible for electron pairing that drive transition into superconducting state for classical low-T c superconductors, have became nearly abandoned and considered to be rather harmful for superconductivity in high-T c cuprates (see, e.g., [12]).Some aspects of d-wave superconductivity can be described within the models of strongly correlated electrons, for example, Hubbard-like or t-J models (e.g., [11,[15][16][17][18][19]), even without explicit account for e-p interactions.The underlying leitmotiv behind the electron correlation treatments has been to understand the phase diagram of high-T c cuprates, that is, the doping process.Introduction of charge carriers (holes or electrons) into the parent antiferomagnetic insulator that causes transition to superconductor (or metal) has been generally accepted to be a universal feature of high-T c cuprates and believed to be a matter intimately related to microscopic mechanism of superconductivity.
Bell-like-shaped dependence of T c on hole doping in the doping range 0.05 ≤ x < 0.27 for family of high-T c cuprates is well known (see, e.g., [20] and references therein).With the exception of YBCO (x = 0.05), the optimal hole doping with maximal T c is x ≈ 0.16.The electron doping is usually less favorable, but basically either hole or electron doping can induce superconductivity.It is, for instance, the case of infinite-layer compound CaCuO 2 (itself is an insulator) where by hole or electron doping in field-induced transistor configuration superconductivity can be induced with T c up to 89 K (or 34 K for electron doping) at about 0.15 charge carriers per CuO 2 .
Without any doubts, charge doping has no-negligible impact on e-e interactions and influences, to some extent, also more subtle spin interactions.Question remains if these are the key effects behind the physics which causes superconducting state transition upon doping?
In this respect, one has to realize that like T c , the lattice parameters and, what is very important, the lattice dynamics is strongly influenced by doping.For instance, dependence of T c on lattice parameter a in case of Hg-based cuprates follows similar dependence as it does T c on hole doping [21].In spite that isotope effect coefficient α for optimally doped cuprates is very small (α ≈ 0.05, exception is YBCO-α ≈ 0.8), doping in underdoped as well in overdoped region for Oisotope effect (all-over the cuprates family), results in a huge changes of isotope effect coefficient α (see, e.g., [22,23].It is important experimental evidence that doping induces very strong changes in lattice dynamics, in particular in dynamics of CuO 2 layers.From theoretical stand point, it means that mechanism of superconducting state transition is a complex matter of electron and nuclear degrees of freedom. The results of high-resolution ARPES study [24,25] of wide family of different high-T c cuprates have brought the other experimental evidence indicating that beside doping, an abrupt change (decrease) of electron velocity near Fermi level, at about 50-80 meV in nodal direction, is the other feature common to high-T c cuprates.The kink in nodal direction is temperature independent.More important, from the stand point of microscopic theory of superconductivity, seems to be formation of temperature-dependent kink on momentum distribution curve (dispersion renormalization) close to Fermi level (∼ 60 meV) in off-nodal direction at transition to superconducting state.It has been reported for Bi2212 [26][27][28][29].Recently, presence of the kink in off-nodal direction has been observed at VUV ARPES study with sub-meV resolution of optimally doped untwinned YBCO in superconducting state [30].
Formation of the off-nodal kink (dispersion renormalization) has been attributed by the authors [26] to the coupling of electrons to a bosonic excitations, preferably they consider a magnetic resonance mode such as observed in some inelastic neutron scattering experiments.The inconsistency in this interpretation has been pointed out by the authors of [27,29].The main points are described in [29] as follows: (a) magnetic resonance has not yet been observed by neutron scattering in such a heavily doped cuprates, and (b) magnetic resonance has little spectral weight and may be too weak to cause the effect seen by ARPES.They agree, however, with the opinion of the authors [26] that the renormalization effect seen by ARPES in cuprates may indeed be related directly to the microscopic mechanism of superconductivity.The authors of [27,29] instead of magnetic resonance mode, attribute the dispersion renormalization to coupling with phonon mode, in particular, with B 1g -buckling mode of CuO 2 plane.The temperature dependence of dispersion renormalization they attribute to the DOS enhancement due to SC-gap opening and to the thermal broadening of the phonon self-energy in normal state.
These results along with the results of neutron scattering [31,32] indicate that also for high-T c cuprates the ep coupling has to be considered as a crucial element of microscopic mechanism of SC state transition.Expressed explicitly, also in case of cuprates the role of phonons at superconducting state transition must not be overlooked.As soon as low-Fermi energy situation occurs (ω ≤ E F ), one can expect important contribution of nonadiabatic vertex corrections at SC state transition.It is beyond the standard ME approximation and this problem has been studied within the nonadiabatic theory of superconductivity [33][34][35].On the other hand, as the ARPES results indicate, electron kinetic energy is decreased and importance of proper treatment of electron-electron Coulomb interactions is increased.The competition between Coulomb and e-p interactions has been intensively studied within the Holstein-Hubbard models [36][37][38][39][40] with both interactions of shortrange character.The obtained results are not satisfactory since heavy-mass polarons are formed that yield low values of T c .It has been improved within the Frohlich-Coulomb model [41] that introduces long-range repulsion between charge-carriers and also long-range e-p interactions.The results show that there is a narrow window of parameters of Coulomb repulsion V c and e-p interactions E P (V c / E P ) resulting in the light-mass bipolarons formation.In this case, according to bipolaron theory of superconductivity [42][43][44], coherent motion of bipolarons represents the supercarrier motion and high T c can be reached.
The McMillan formula, which is very good approximation for T c of elementary metals and their alloys [45], is often used also for calculation of critical temperature of high-T c superconductors within the BCS-generic framework.It has been shown [46] that in strong-coupling regime λ 1, T c can be as large as k B T c = Dωλ 1/2 /2π.However, in reality there is problem with correct estimation of Coulomb pseudopotential μ * and with unrealistically large values of λ that would match high experimental T c of novel superconductors.It has to be stressed, however, that strong coupling regime λ > 1 violates adiabatic condition ω/E F 1 of the ME approximation, which is behind the derivation of the McMillan formula.
More over, new class of superconductors, for example, cuprates, fullerides, and MgB 2 are systems that are rather pseudoadiabatic with sizeable adiabatic ratio ω/E F < 1 [47], in contrast to elementary metals, where adiabatic condition ω/E F 1 is perfectly fulfilled.This situation indicates the importance of nonadiabatic contributions at calculation of ep interactions within the BOA, an effect which is beyond the standard ME approach.As mentioned above, formulation of the nonadiabatic theory of superconductivity by Pietronero et al. [33-35, 48, 49], which accounts for vertex corrections and cross phonon scattering (beyond ME approximation), has solved this nontrivial problem by generalization of Eliashberg equations.The theory, which is nonperturbative in λ and perturbative in λω D /E F , has been applied at simulation and interpretation of different aspects of high-T c superconductivity [50][51][52][53][54]. Basically, it can be concluded that accounts for nonadiabatic effects in quasiadiabatic state ω/E F ≤ 1 is able to simulate different properties of high-T c superconductors, including high-value of T c , already at relatively moderate value of e-p coupling, λ ≈ 1.Moreover, it has also been shown that increased electron correlation is important factor that makes corrections to vertex function positive, which is in this context crucial for increasing T c .
Nonetheless, sophisticated treatment of high-T c superconductivity within the nonadiabatic theory faces serious problem related to possibility of polaron collapse of the band and bipolaron formation.According to bipolaron theory of Alexandrov [42-44, 47, 55-58], polaron collapse occurs already at λ ≈ 0.5 for uncorrelated polarons and even at smaller value for a bare e-p coupling in strongly correlated systems.For ω/E F ≤ 1, or λ ≥ 1, and for ω/E F ≥ 1 at whatever small value λ 1, the nonadiabatic polaron theory has been shown to be basically exact [58].Bipolarons can be simultaneously small and light in suitable range of Coulomb repulsion and e-p interaction [59].These results have important physical consequences.There are serious arguments that effect of polaron collapse cannot be covered through calculation of vertex corrections due to translation symmetry breaking and mainly, polaron collapse changes possible mechanism of pair formation, that is, instead of BCS scenario with Cooper pair formation in momentum space, the BEC with mobile bipolarons (charged bosons) in real space becomes operative.
Discovery of superconductivity in a simple compound MgB 2 at 40 K [60] has been very surprising and has started a new revitalization of superconductivity research.Beside the many interesting aspects, discovery of the MgB 2 superconductivity is, in my opinion, crucial for general theoretical understanding of SC state transition on microscopic level.It is related to band structure (BS) fluctuation and dramatic changes of BS topology at e-p coupling.
The σ bands splitting at coupling to E 2g mode in MgB 2 has been reported [61] already in 2001 but, with exception of possible impact of anharmonicity [62], no special attention has been paid to this effect.Superconductivity in MgB 2 has been straightforwardly interpreted [63] shortly after the discovery as a standard BCS-like, even of intermediatestrong coupling character.For clumped nuclear equilibrium geometry, the BS is of adiabatic metal-like character.The E σ F of σ band electrons (chemical potential μ) is relatively small, ≈ 0.4 eV, but still great enough comparing to vibration energy of E 2g phonon mode, ω 2g ≈ 0.07 eV.Thought the adiabatic ratio ω/E F ≈ 0.15 is sizable, it is small enough in order to interpret superconductivity within the BCSgeneric framework.It is supposed that nonadiabatic effects, anharmonic contributions and/or Coulomb interactions within generalized Eliashberg approach should be important in this case.On the other hand, the value of e-p coupling, λ ≈ 0.7 indicates that polaron collapse can be expected and superconductivity should be of nonadiabatic bipolaron character rather than the BCS-like.
Nevertheless, the matter is even more complicated.It has been shown [64,65] that analytic critical point (ACP-maximum, minimum, or saddle point of dispersion, in case of MgB 2 it is maximum) of σ band at Γ point crosses Fermi level (FL) at vibration displacement ≈ 0.016 A • /B-atom, that is, with amplitude ≈ 0.032 A • , which is smaller than root-mean square (rms) displacement (≈0.036A • ) for zeropoint vibration energy in E 2g mode.It means, however, that in vibration when ACP approaches FL on the distance less than ±ω, the adiabatic Born-Oppenheimer approximation (BOA) is not valid.In this case, Fermi energy E σ F of σ band electrons (chemical potential μ) close to Γ point is smaller than E 2g mode vibration energy E σ F < ω 2g and at the moment when the ACP of the band touches Fermi level, the Fermi energy is reduced to zero, E F → 0.Moreover, shift of the ACP substantially increases density of states (DOS) at FL, Advances in Condensed Matter Physics EF , and induces corresponding decrease of effective electron velocity (∂ε 0 σ /∂k) EF of fluctuating band in this region of k-space.From the physical stand point, it represents transition of the system from adiabatic ω ≤ E F into true nonadiabatic ω > E F , or even into strong antiadiabatic state with ω E F .This effect has crucial theoretical impact.At these circumstances, not only ME approximation is not valid (including impossibility to calculate nonadiabatic vertex corrections which represent off-diagonal corrections to adiabatic ground state) but adiabatic BOA itself does not hold as well.
The BOA is crucial approximation of theoretical molecular as well as of solid-state physics.It enables to solve many-body problem via separation of electronic and nuclear motion and to study electronic problem in a field which is created by fixed nuclei.
On the level of the BOA, the motion of the electrons is a function of the instantaneous nuclear coordinates Q, but is not dependent on the instantaneous nuclear momenta P. Usually, and in solid state physics basically always, only parametric dependence is considered, that is, nuclear coordinates are only parameters in solution of electronic problem within the clamped nuclei Hamiltonian treatment.Nuclear coordinate-dependence, when explicitly treated, modifies nuclear potential energy by so called diagonal BO correction (DBOC) that reflects an influence of small nuclear displacements out of the equilibrium positions and corrects the electronic energy of clamped nuclear structure.The DBOC enters directly into the potential energy term of nuclear motion (but leaves unchanged the nuclear kinetic energy) and in this way modifies vibration frequencies.The off-diagonal terms of the nuclear part of system Hamiltonian that mix electronic and nuclear motion through the nuclear kinetic energy operator term are neglected and it enables independent diagonalization of electronic and nuclear motion (adiabatic approximation).Neglecting the off-diagonal terms is justified only if these are very small, that is, if the energy scales of electron and nuclear motion are very different and when adiabatic condition holds, that is, ω/E 1.If necessary, small contribution of the off-diagonal terms can be calculated by perturbation methods as so called nonadiabatic correction to the adiabatic ground state.
Situation for superconductors seems to be substantially different, at least in case of the MgB 2 .There is considerable reduction of electron kinetic energy, which for antiadiabatic state results even in dominance of nuclear dynamics (ω E F ) in some region of k-space.In this case, it is necessary to study electronic motion as explicitly dependent on the operators of instantaneous nuclear coordinates Q as well as on the operators of instantaneous nuclear momenta P. It is a new aspect for many-body theory.
The electronic theory of solids has been developed with the assumption of validity of the adiabatic BOA.In this respect, it is natural that different theoretical-microscopic treatments of superconductivity based on model Hamiltonians which stress importance of one or the other type of interaction mechanism, implicitly assume validity of the BOA, and it is very seldom that possibility of the BOA breakdown at transition to SC state is risen.The notion "nonadiabatic" effects in relation to electronic structure is commonly used for contributions of the off-diagonal matrix elements of interaction Hamiltonian (e.g., e-p coupling, ee correlations, etc.) to the adiabatic ground state electronic energy calculated in second and higher orders of perturbation theory and does not account for true nonadiabaticantiadiabatic situation, ω > E F .
In this connection, a lot of important questions arise, as the following examples.
(i) How to treat antiadiabatic state?
(ii) Can be system stable in antiadiabatic state?(iii) Are the physical properties of the system in antiadiabatic state different from the corresponding properties in adiabatic state?
(iv) What is the driving force for adiabatic ↔ antiadiabatic state transition, that is, which type of interaction mechanism and at which circumstances trigger this type of transition?Theoretical aspects related to the above problems have been elaborated and discussed in details within "Ab initio theory of complex electronic ground state of superconductors", which has been published in the papers, [66][67][68].The main theoretical point is generalization of the BOA by sequence of canonical-base functions transformations.The final electronic wave function is explicitly dependent on nuclear coordinates Q and nuclear momenta P, or alternatively, emerging new quasiparticles-nonadiabatic fermions are explicitly dependent on nuclear dynamics.As a result, the effect of nuclear dynamics can be calculated in a form of corrections to the clamped nuclei ground state electronic energy, to the one-particle spectrum and to the two-particle term, that is, to the electron correlation energy.

Note.
To avoid confusion, it should be stressed that the notion electron correlation energy as used in this paper stands for improvement of e-e interaction term contribution beyond the Hartree-Fock (HF) level, E corr = E exact − E HF (E exact < E HF ), as it can be calculated, for example, by (1/r)perturbation theory in the second and the higher orders, or by configurations interaction method.In condensed matter physics, electron correlation usually stands for an account for Coulomb e-e interaction at least on Hartree or HF level.On the HF level not only repulsive e-e term is present, like on Hartree level where spin is not considered at all, but also exchange term, fermion Coulomb-hole only for electrons with parallel spins.Correlation energy improves an account for unbalanced treatment of e-e interaction for electrons with parallel and antiparalel spins on HF level.
It has been shown that due to e-p interactions, which drive system from adiabatic into antiadiabatic state, adiabatic symmetry is broken and system is stabilized in the antiadiabatic state at distorted geometry with respect to the adiabatic equilibrium high symmetry structure.Stabilization effect is due to participation of nuclear kinetic energy term, that is, it is the effect of nuclear dynamics (dependence on P) which is absent in the adiabatic state on the level of the BOA.The antiadiabatic ground state at distorted geometry is geometrically degenerate with fluxional nuclear configuration in the phonon modes that drive system into this state.It has been shown that while system remains in antiadiabatic state, nonadiabatic polaron-renormalized phonon interactions are zero in well defined k-region of reciprocal lattice.Along with geometric degeneracy of the antiadiabatic state, it enables formation of mobile bipolarons (in a form of polarized intersite charge density distribution) that can move over lattice in external electric potential as supercarriers without dissipation.Moreover, it has been shown that due to e-p interactions at transition into antiadiabatic state, k-dependent gap in one-electron spectrum has been opened.Gap opening is related to the shift of the original adiabatic Hartree-Fock orbital energies and to the kdependent change of density of states of particular band(s) at the Fermi level.The shift of orbital energies determines in a unique way one-particle spectrum and thermodynamic properties of system.It has been shown that resulting oneparticle spectrum yields all thermodynamic properties that are characteristic for system in superconducting state, that is, temperature dependence of the gap, specific heat, entropy, free energy, and critical magnetic field.The k-dependent change of the density of states at the Fermi level in transition from adiabatic (nonsuperconducting) into antiadiabatic state (superconducting) can be experimentally verified by ARPES or tunneling spectroscopy as spectral weight transfer at cooling superconductor from temperatures above T c down to temperatures below T c .
Results of the ab initio theory of antiadiabatic state have shown that Fröhlich's effective attractive electron-electron interaction term represents correction to electron correlation energy in transition from adiabatic into antiadiabatic state due to e-p interactions.Analysis of this term has shown that increased electron correlation is a consequence of stabilization of the system in superconducting electronic ground state, but not the reason of its formation.
In the present article, the key points of the theory are recapitulated and the adiabatic ↔ antiadiabatic state transition is shown to be operative for different types of superconductors.

Electronic Ground State Beyond the Born-Oppenheimer Approximation
2.1.Preliminary Remarks.Development of the theory of molecules and solids has been enabled due to fundamental approximation, the Born-Oppenheimer approximation (BOA).With respect to electronic structure of superconductors and transition to superconducting state, some aspects of this approximation should be outlined at the beginning.Solution of the Schrödinger equation of many-body system composed of N e electrons and N n nuclei (total system) with the Hamiltonian and wave function Ψ(r, R), which is a general function of the sets of electron {r} and nuclear {R} Cartesian coordinates, is possible at the assumption of validity of the Born-Oppenheimer approximation (BOA).The BOA, originally formulated in the work in [69] by power expansion of potential surface for nuclear motion at equilibrium geometry with respect to displacement and electron/nuclear mass ratio (m e /M n ) 1/4 , has been reformulated later by Born [70,71] in a more practical form.According to it, the wave function of the total system (1) can be expressed in the factorized form as a linear combination of known adiabatic electronic wave functions {Φ m (r, R)} that are the eigenfunctions of clumped nuclei electronic Schrödinger equation with the electronic Hamiltonian: at fixed nuclear configuration [R].Expansion coefficients {χ m (R)} in (3), regarded as unknown, are nuclear wave functions for nuclear configuration R with the electronic subsystem in particular adiabatic electronic state Φ m (r, R).
The R-dependence in ( 5) is only parametric and in general it should be calculated over the full configuration space {R} in order to calculate the adiabatic potential hypersurface {E e m (R)} of nuclear motion.With respect to (3) and (4), the Schrödinger equation of the total system (1) for electronic state Φ n (r, R) can be written in the form with the term Λ nm (R), which couples electronic and nuclear motions, on the right-hand side (rhs): The term E te n in ( 6): is total adiabatic electronic energy, that is, adiabatic electronic energy plus nuclear Coulomb repulsion at nuclear configuration R.
Until the Born approach (3) is valid, ( 6) is exact and it still describes coupled motion of electrons and nuclei over the term (7), Λ nm (R), which represents possibility of transitions between different adiabatic electronic states, Φ n (r, R) ↔ Φ m (r, R), due to nuclear motion (Rdependence).If such transitions are forbidden from the symmetry reasons, or if there is physically reasonably justified assumption that contributions of such transitions are negligibly small, then one can omit the rhs term, and ( 6) can be written in the diagonal form: where B nn is the only a nonzero diagonal contribution of the Λ term The term (10) is the mean-value of the nuclear kinetic energy for adiabatic electronic state Φ n (r, R) at nuclear configuration {R} and represents so called adiabatic diagonal Born-Oppenheimer (DBOC) correction to the total adiabatic electronic energy E te n (R).However, ( 9) is then the equation of motion of nuclei and it has the form of Schrödinger equation with Hamiltonian: The effective-adiabatic potential for nuclear motionV eff NN (R): is represented by the total electronic energy (8), that is, Coulomb potential energy of the bare nuclei repulsion E NN (R) and adiabatic electronic energy E e n (R), which is corrected by mean-value of the nuclear kinetic energy (DBOC) for the particular adiabatic electronic state Φ n (r, R) (10).Contribution of the adiabatic DBOC is usually neglected as very small quantity, mean value of kinetic energy of slowly moving heavy nuclei in contrast to fast motion of light adiabatic electrons.
At these circumstances, the motion of electrons and nuclei is effectively decoupled, that is, it is possible to realize an independent diagonalization of the electronic Schrödinger equation ( 4) and nuclear Schrödinger equation (9), electrons and nuclei of the system behave like two statistically independent sets.Assumed small contribution of the rhs term of ( 6) can be calculated by some approximate way, usually by perturbation theory.Then, problem with the Hamiltonian: is studied, where is unperturbed part, and is a small perturbation.
In practice, physical and/or chemical properties of a many-body system in its ground electronic state Φ 0 (r, R) are of the prime interest.In this case, the Born approach (3) is usually restricted to the single term and the total wave function of system is a simple product of the adiabatic electronic wave function Φ 0 (r, R) and corresponding nuclear wave function: The Born approach in the form ( 16) is the adiabatic approximation.In a common sense, what is usually called the Born-Oppenheimer approximation (BOA, or clamped nuclei-crude adiabatic BOA) is the adiabatic approximation where the contribution of the DBOC (i.e., B 00 (R) term in (12)) is also neglected and wave function (16) has the form Contributions of the off-diagonal terms Λ 0m (R) in (15) which are calculated as a small perturbation to the Hamiltonian (14) represent a nonadiabatic correction to the unperturbed adiabatic ground state.The conditions at which the nonadiabatic correction can be expected to be small, and the BOA (16, and in general 3) is valid, can be estimated by analysis of the second-order contributions to the energy of the total system, E TS 0 (R), which are small providing that It can be derived (see, e.g., [72]), based on the expansion of the effective nuclear potential (V eff NN (R) = E te 0 (R)) at equilibrium nuclear geometry R eq (where the total electronic ground-state energy reaches its absolute minimum) at least up to the quadratic term in a displacement ΔR: The result is that (18) holds and the BOA ( 16), ( 3) is valid, if inequality is fulfilled for electronic and vibration (phonon) energy spectrum of a system.The meaning of ( 20) is clear, the electronic spectrum, that is, the differences between the total electronic energies of the excited electronic states and the ground-state energy has to be much greater than vibration (phonon) energy spectrum of system.The long-time experience of theoretical molecular and solid-state physics has shown that for a ground electronic state of vast majority of molecular systems and solids at equilibrium geometry R eq , the BOA is absolutely good approximation.In case of solids, it has enabled to derive (Bloch's assumption of small perturbation of the periodic lattice potential for small nuclear displacement out of equilibrium) the field theory Hamiltonian of total system (13) in a well-known form, H = H 0 + H ep , with The perturbation Hamiltonian H ep represents now, instead of nuclear kinetic energy Λ-perturbation term (15), an electron-phonon (e-p) interaction term: In this case, no special attention has been paid to the diagonal correction (10) to the total electronic ground-state energy, and at derivation of ( 21) and (22), this correction has been omitted as negligibly small quantity.Moreover, since perturbation H is now introduced as an e-p interaction H ep , it is immediately seen from the form of ( 16) and ( 22) that the first-order correction, that is, diagonal perturbation term, equals zero: Φ 0 | H ep |Φ 0 = 0.All interesting physics is then related to higher-order contributions with participation of excited electronic states, that is, the first possible nonzero contributions are in the second order of perturbation theory, that is, terms of the form Hamiltonian (21), (22) has been the starting point for theoretical study of the effects connected to e-p interactions in solids.By means of the famous unitary transformation of this Hamiltonian, Fröhlich has derived [6,7] an effective electron-electron interaction term which is a crucial element at formulation of the BCS [5] or Migdal-Eliashberg (ME) theory of superconductivity [8][9][10].The extension of the ME theory, that is, inclusion of higher-order contributions of e-p interactions by means of Feynman-Dyson perturbation expansion of H ep (vertex corrections) is the basis for so called "nonadiabatic" theory of superconductivity, as it has been developed by 48,49].It should be pointed-out that not only the above-mentioned standard and "nonadiabatic" theories of superconductivity but also other theoretical models of superconductivity which do not consider explicitly e-p interactions assume tacitly validity of the BOA as soon as the model Hamiltonian is written in the form no matter what perturbation H represents.The reason for short sketch of the BOA in the introductory part is to attract an attention toward some aspects, which at study of solids, in particular of superconductors, are tacitly assumed to hold implicitly and seemingly there are no indications raising doubts that this class of solids should be an exception.
The main aspect is the R-dependence of the BOA.Whatever trivial it seems to be, it has to be stressed that the study of many-body system by means of the clumped nuclei Hamiltonian ( 21)-( 23) requires validity of the Born ansatz ( 16), (3), that is, fulfillment of (20) to hold not only for equilibrium nuclear geometry R eq but also for displaced geometry (R eq + ΔR).Displacement ΔR has to be as large as to cover full configuration space {R} experienced by nuclei in the vibration (phonon) modes of system.At 0 K temperature, ΔR has to be greater than (or at least equal to) root-mean square displacement x 0(ν) of vibration (phonon) mode (ν) at zero-point energy: It means that for validity of the BOA, the inequality for the energy spectrum which, along with (20), has to be also fulfilled is In this context, the adiabatic electronic energies of the ground and excited electronic states {E te 0 , (E te n )}, at fixed nuclear geometry [R], need also some comment.Without the loss of generality, the Hartree-Fock (HF) calculation scheme in direct-real space orbital representation can be assumed.The adiabatic wave function of the electronic ground state is represented by the Slater determinant that satisfies requirement of antisymmetry: where for simplicity a closed-shell, N e = 2n e electron system, is considered.In (26), lowest laying n e spatial orbitals are occupied, each being occupied twice, once by an electron with α spin, ϕ n , and once with β spin, ϕ n .Adiabatic total electronic energy of the ground state is then the sum of adiabatic electronic energy and Coulomb repulsion E NN (R 0 ) of bare nuclei at fixed geometry [R] = R 0 : Summation indices {I, J} in ( 27) run over the occupied states {ϕ I , ϕ J }.
Set of {ε 0 I } in ( 27), is one-electron spectrum of the adiabatic electronic ground state (26), that is, orbital energies {ε 0 I } are eigenvalues: of the Fock operator and {ϕ Q } is the set of corresponding eigenfunctions (orthonormal set of optimized orbitals).In (28), h PP is diagonal element of one-electron core Hamiltonian h core (integrals of electronic kinetic energy term plus electronnuclear Coulomb attraction term), and J Q , K Q are 2-electron (1,2,3, or 4 nuclear center) Coulomb repulsion and exchange integrals.
In the next step, instead of independent calculation of electronic excited states energies {E e n( / = 0) (R)}, which would require new optimization of excited state wave functions Φ m (r, R) (it should be extremely complicated since excited state wave functions have to be orthogonal to the ground state wave function), in practice a fairly good approach is used which is based on the orbitals {ϕ Q } already optimized for the ground-state wave function Φ 0 (r, R 0 ).Optimization is always performed within a basis set with some finite number n of basis set functions |μ n , while n > N e .Diagonalization of (29) then yields {n} eigenfunctions {ϕ n (R)} and corresponding eigenenergies {ε 0 n (R)}, but only n e lowest-lying spatial orbitals {ϕ I (R)} (ne) are occupied and remaining (n − n e ) are virtual-unoccupied orbitals {ϕ A (R)} (n−ne) .By promotion of electron(s) from occupied orbital(s) {I, J, . ..} to virtual-unoccupied orbital(s) {A, B, . ..}, excited state(s) configuration(s) {Φ A } as a linear combination of corresponding Slater determinants {Φ I → A } can be constructed.It can be shown that, for example, single-electron excitations yield two excited state electronic configurations-lowest lying excited state, that is, singly excited triplet state 3 ΦA(I → A) , and singly excited singlet state 1 ΦA(I → A) .Differences in the electronic energies of these excited state configurations with respect to the electronic energy of the ground state are as follows: (i) for singly excited triplet state, (i) for singly excited singlet state, From ( 30), (31), it is clear that for approximations which do not consider explicitly for two-electron terms, the differences in energies of singly excited triplet and singlet states with respect to the ground-state energy are the same: Multiple electronic excitations can be calculated in a similar way, by generation of Slater determinants of p-particle, h(≡ p)-hole states in the language of particle-hole formalism.
The important point is that clumped nuclei electronic ground-state energy calculation provides approximate information about electronic excited states over the one-electron spectrum which corresponds to the electronic ground state.In particular, over the optimized set of occupied (ϕ I ) and unoccupied-virtual (ϕ A ) spinorbitals.
In the simplest form, with respect to (32), the inequalities (20), (25) which have to be valid for save application of the BOA, are in the terms of the optimized ground state orbital energies expressed as follows: (i) for system at equilibrium geometry R eq , has to hold for the couple of frontier orbitals, that is, for highest occupied ϕ I ≡ ϕ HOMO and lowest unoccupied ϕ A ≡ ϕ LUMO orbitals; (i) more over, has to hold for the same system at displaced nuclear configuration with respect to particular vibration(phonon) mode ν.
In case of solids with quasicontinuum of states, in momentum k-space representation, these inequalities can be rewritten in the form for equilibrium nuclear geometry, and Dω r (36) for the same system at displaced nuclear configuration in particular phonon mode r.These inequalities have to be valid in a multiband system for each band S and energies {ε 0 S (k c )} of analytic critical points (ACP-i.e., absolute and local minima and maxima or inflex points) located on the particular band-dispersion curve of the first BZ with respect to the energy of the Fermi level ε 0 F .An important aspect should be reminded at this place.Inequalities ( 33)-( 36) have been derived on the basis of the expansion ( 19) of the effective nuclear potential (i.e., total electronic energy E te 0 (R)) at R eq .With respect to the inner structure of the total electronic energy (27), it is clear that in principle, the expansion (19) can be performed (i.e., smallness in energy change ΔE te 0 (R eq +ΔR) as to ensure at least harmonic approximation) even if displacement ΔR yields significant changes in the one-electron spectrum (orbital energies) of the ground state provided that these changes are well balanced by changes Δh II and ΔE NN , eventually also by changes Δ{J IJ , K IJ } if two-electron terms are explicitly included.The main point is that for the electronic ground state, the one-electron spectrum at R eq , {ε 0 Q } Req , can be significantly different from the one-electron spectrum {ε 0 Q } Req+x0 at displaced nuclear configuration (R eq + x 0,ν ) for some phonon mode ν.
In respect of it, possibility that displacement x 0,(ν) yields for some band situation when ≤ Dω ν (38) even if for the equilibrium geometry R eq (35) holds, is not excluded at all.The reason of possibility for such substantial changes of the electronic (band) structure is hidden in the chemical composition and structure of particular system.With respect to ( 27)- (31), it is connected to the changes of one-electron core part (Δh IJ ) and/or two-electron Coulomb repulsion and exchange integrals (ΔJ IJ , ΔK IJ ).Which of these are dominant (hopping terms, onsite/intersite repulsion, electron correlations, dynamic screening, etc.) can hardly be determined on the basis of model Hamiltonians without invoking for complete chemical composition and structure of the particular system.
Occurrence of such a situation means that nuclear motion (nuclear vibration, in particular, phonon mode) has induced sudden decrease of effective electron velocity (∂|ε 0,a (k)|/∂k) kc , in particular band close to analytic critical point k c at the Fermi level ε 0,F .Or, what is equivalent, the effective mass of electrons have been increased in this region.In any case, however, it means that at displaced nuclear geometry (R eq + x 0,(ν) ), electrons, and nuclei of the system do not move as independent particles, that is, in this region, electrons are not able to follow nuclear motion adiabaticaly.The motion of electrons and nuclei is now strongly correlated.The new situation occurs, the adiabatic BOA is not valid for nuclear geometry (R eq + x 0,(ν) ), and wave function cannot be factorized in the form ( 16) and (3), that is, and standard adiabatic form of the Hamiltonian ( 21)-( 23) cannot be applied to study an electron-phonon problem.
What is crucial, however, is the fact that the nonadiabatic effect cannot be calculated as a perturbation, that is, as a nonadiabatic correction to the adiabatic ground state at R eq .In this case, system is in true nonadiabatic or antiadiabatic state.

Generalization of the BOA Beyond Adiabatic Regime.
In order to solve the problems which are sketched above, one needs to study electronic motion as explicitly dependent on nuclear dynamics, that is, dependent on instantaneous nuclear coordinates Q and nuclear momenta P. It means that electronic state has to be explicitly dependent on nuclear operators Q and P. Approximate solution of this problem can be found in [66].It is based on the sequence of basis functions transformations, starting with fixed basis set (Q, P-independent) over Q-dependent base up to QP-dependent basis set.A base function transformation is canonical transformation and it ensures that each transformation step preserves corresponding statistics.Moreover, since we start from fixed basis set (clamped nuclei situation) which enables factorization of total system wave function, this property is required to be preserved at each transformation step up to quadratic terms in Q and P. Consequently, it makes the solution of Q, P-dependence tractable.It should be stressed that base functions transformation is equivalent to the quasiparticle transformation which we have formulated in 1992 [67].
Main aspects and the results of this treatment are presented below.
General nonrelativistic form of system Hamiltonian (2) can be written in second quantization formalism as Nuclear potential energy E NN and one-electron core term h PQ (electron kinetic energy plus electron-nuclear coulomb attraction term) are functions of the nuclear coordinate Q r operators (normal modes {r} nuclear displacements out of fixed nuclear geometry R 0 ) and nuclear kinetic energy T N is a quadratic function of the corresponding Advances in Condensed Matter Physics nuclear momenta operators P r .For nuclear coordinate and momentum operators hold: The (a, a + ), (b, b + ) are standard fermion and boson annihilation and creation operators, respectively.The Q-dependence of terms E NN (Q) and h PQ (Q) in ( 40) can be expressed through the Taylor's expansion at fixed nuclear configuration R 0 : The term E 0 NN (R 0 ) represents potential energy of bare nuclei at fixed nuclear configuration R 0 , h 0 PQ (R 0 ) is one-electron core term (electron kinetic energy plus electron-nuclear coulomb attraction term) at fixed nuclear configuration R 0 and terms {u (i)  PQ (Q)} are related to matrix elements of electron-vibration (phonon) coupling (u (i) PQ ), that is, Two-electron terms v 0 PQRS (electron-electron coulomb repulsion and exchange integrals) do not depend explicitly on the nuclear operators.

Crude-Adiabatic Approximation: Clamped Nuclei Solution.
Since the crude-adiabatic approximation is the reference level for study the effect of nuclear dynamics on electronic structure, some details should be introduced.
With respect to sequence of transformations, which will be presented in the following parts, let as distinguish particular representation; that is, dependence of electronic states on nuclear operator Q only by double-bar over particular operator symbol (e.g., a,b), dependence on both Q and P will be without bar (a, b) and for fixed basis (independent on Q and QP, resp.) the operators will be with a single-bar (a, b) over the particular symbol.In this sense, the operators in crude-adiabatic approximation are written as single-bar operators, that is, (a, b, Q, P).
Provided that phonon and electronic energy spectrum are well separated and (20) holds for relevant configuration space R near to R eq , crude-adiabatic (clamped nuclei) treatment is justified.Electronic Hamiltonian is now only parametrically dependent on nuclear configuration, that is, nuclear geometry is fixed at nuclear configuration R 0 : Application of Wick's theorem to the product of creation and annihilation operators yields for particular terms the normal product form (N-product) with corresponding contractions: Application of Wick's theorem introduces renormalized Fermi vacuum, that is, the total set of orthonormal base orbitals (ϕ P , ϕ Q , . ..) is divided on two distinct groups; the set of occupied (ϕ I , ϕ J , . ..) spinorbitals which form a Fermi see (FS) and set of unoccupied-virtual (ϕ A , ϕ B , . ..) spinorbitals which are above the FS.
The electronic Hamiltonian (44) can now be written in a quasiparticle form as a sum of zero, one-and two-particle terms: (I) Electronic ground state energy.
The scalar quantity in this Hamiltonian, that is, zeroparticle term H (0) (R 0 ), is the result of the operators contractions and has the form This term represents ground-state electronic energy calculated by unrestricted Hartree-Fock SCF (HF-SCF) procedure at fixed nuclear configuration R 0 .It is obvious that addition of nuclear repulsion E NN (R 0 ) to this term yields total electronic energy E te 0 (R 0 ).Electronic ground state is represented by renormalized Fermi vacuum Φ 0 .It is an antisymmetric electronic wave function, that is, expressed in the form of single Slater determinant constituted by lowest laying occupied spinorbitals {ϕ I } of complete orthonormal base {ϕ P }: (II) Electronic energy spectrum.
The one-particle term H (1) (R 0 ) = H (cad) (R 0 ) of the electronic Hamiltonian (47) has the form, The elements F PQ are calculated as matrix elements of Hartree-Fock one-particle operator: Diagonalization of ( 50), F PQ = ε 0 P δ PQ , that is, solution of electronic Hartree-Fock equations by HF-SCF procedure: yields set of eigenvalues, that is, HF-orbital energies: and corresponding set of eigenfunctions {ϕ P }, the orthonormal set of optimized spinorbitals.It means that one-particle term (50) can always be written in diagonal form and represents a set of unperturbed eigenenergies of system at fixed nuclear geometry R 0 , that is, for excited states hold: In an approximate way, the one-particle Hamiltonian H (1) (R 0 ) represents complete electronic spectrum expressed over occupied and unoccupied-virtual spinorbitals calculated for electronic ground state Φ 0 (r, R 0 ) by HF-SCF procedure.In particular, n-electron excited state wave function Φ hn pn can be constructed, as it has already been mentioned, by promotion of n electrons from n occupied spinorbitals to n unoccupied spinorbitals (i.e., the same number n of holes (h n ) and particles (p n ) are created).Electronic energy of such excited state can be calculated through HF-eigenenergies that correspond to optimized spinorbitals of the ground state.It should be pointed that it is only an approximate way for description of true electronic excited states.
In terms of orbital energies (51c), for total electronic ground state energy, E te 0 (R 0 ) holds: For closed-shell electronic systems, within restricted HF approximation (doubly occupied spatial orbitals with α and β spins, (27) holds.Total electronic ground-state energy of the system reaches the minimum at some equilibrium nuclear configuration R 0 = R eq .Corresponding Slater determinant (49) represents wave function of the electronic ground state at equilibrium nuclear configuration.
(III) Electron correlation energy.(see Note in Section 1).The third term of the electronic Hamiltonian (47), that is, the two-particle term has the form Formally it looks like standard coulomb electron-electron interaction term in (44).With respect to the fact that after application of Wick's theorem (45) the renormalized Fermi vacuum has been introduced and zero-particle (scalar) quantity (48) represents electronic energy of the ground state that accounts also for coulomb electron-electron interactions (48) and one-particle term is diagonal (51d) and represents unperturbed HF-orbital energies {ε 0 P } of the system (oneelectron energy spectrum), then two-particle term (53) represents perturbation part of the electronic Hamiltonian (47).Since perturbation (53) contains only electron-electron interaction term, contributions of this term represent electron correlation energy of the system in its ground electronic state.In this respect, electron correlation energy is treated as a perturbation.Calculation of the electron correlation energy up to higher order of perturbation theory is usually done by diagrammatic many-body perturbation theory.For correlation energy in second order of perturbation theory, analytic expression for closed-shell system can be derived in a simple form Related to energy of the Fermi level, energies of unoccupied states are {ε 0 A } > 0 and energies of occupied states are {ε 0 I } < 0. It means that E corr is negative, that is, it decreases electronic energy of the ground state (this contribution corrects and stabilizes total electronic energy of the system ( 27), ( 52)).This holds for arbitrary nuclear geometry R 0 until system remains in a bound state.

Advances in Condensed Matter Physics
Solution of nuclear problem (9) on crude-adiabatic level is now straightforward.
Nuclear part of Hamiltonian is For different, but fixed nuclear configurations [R] at R 0 = R eq , E te 0 (R) is calculated and effective potential for nuclear motion is derived V eff NN (R).It has to be stressed that contributions of DBOC-B 00 (R), that is, (10), to V eff NN (R) are in this case neglected.
Energy and wave function of the total system in the ground electronic state are Nuclear wave function is

Base Transformation-Introduction of New Dynamical
Variables.Each eigenfunction ϕ J in Slater determinant (49) can be expressed as a linear combination of atomic orbital (AO) basis functions {μ}, ϕ P = μ c μP |μ , which are fixed at the positions of the particular nuclei in clamped nuclei configuration R 0 .It represents fixed basis set {μ(x, 0) ≡ μ(x, Q)}.
In second quantization, with single-bar a + being creation operator of the crude-adiabatic electron (clamped nuclei approximation, independent on Q and QP), base function μ(x, 0) can be written with respect to "absolute" Fermi vacuum |0 as and For solids, in electronic quasimomentum k-space, the basis functions are Bloch-periodic orbitals {μ(k, x, 0)}: In (60b), tR is translation vector and μ(x − tR) = μ(x).The set of {μ(x)} is fixed basis set {μ(x, 0)} at frozen nuclear configuration R 0 (Q = 0).An eigenfunction ϕ P is then crystal orbital CO-ϕ P (band ϕ P ), which is a linear combination of the Bloch-periodic basis functions {μ(k, x, 0)}.In a second quantization, it has the form In this case, occupancy of the band is not distinguished by the index (subscript) of the band, but it is determined by the value of k-vector of particular band dispersion ε 0 P (k) with respect to the energy of Fermi level ε 0 F .It is obvious that Hamiltonian of system ( 40), ( 44) is now Hamiltonian which corresponds to the composition of unit cell and total electronic energy calculation results in corresponding band structure.
For fermion and boson creation and annihilation operators, the standard anticommutation, and commutation relations hold: Simultaneous diagonalization of electronic and nuclear part of system Hamiltonian, that is, factorized form of system wave function (57) within adiabatic BOA implies also validity of the following commutation relations: In case of crude-adiabatic approximation, the electrons "see" the nuclei at theirs instantaneous positions at rest and nuclei do not "feel" internal dynamics of electrons.Within the spirit of the BOA, it would be correct if the electrons follow nuclear motion instantaneously, that is, electronic state has to dependent explicitly on instantaneous nuclear positions.In this case, the wave function of the system, instead of the form (57) with Q-independent electronic part should be replaced by Q-dependent form, that is, Advances in Condensed Matter Physics 13 Adiabatic, nuclear displacement Q-dependent electronic wave function Φ 0 (r, Q) in ( 62) assumes existence of complete orthonormal basis set {ϕ R (x, Q)}, that is, validity of the following relations: Now, electron creation and annihilation operators which correspond to the Q-dependent moving base are denoted as double-bar operators (a + , a).Also, the boson operators related to the Q-dependent moving base are written as double-bar operators: Then, Since adiabatic electrons remain fermions, the operators have to obey standard fermion anticommutation relations: (66).
Crude-adiabatic electronic wave function Φ 0 (r, 0) which does not depend on the nuclear displacements Q is expanded over fixed basis set with spinorbitals {ϕ R (x, 0)} that are eigenfunctions of clamped nuclear electronic Hartree-Fock equations (30).This is complete and orthonormal basis set Crude-adiabatic fermion creation and annihilation operators that correspond to the fixed basis set are written as single-bar operators (a + , a), that is, or in short-hand notation These operators pertain to ordinary (crude-adiabatic) electrons and the standard anticommutation relations (61a) hold.
Due to properties ( 63), (67), the two bases are interconnected by the base transformation of the following form: Now, for fermion operators in Q-dependent moving base one can write 69), (70) are Since then due to closure property and orthonormality ( 63), (67) of the bases, it can be derived that base transformation matrix C(Q) is an unitary matrix: It can be shown [66] that the base transformation is identical with canonical transformation of operators (see (70) and ( 74)) It ensures preservation of statistics, that is, validity of the anticommutation relations (66) for new, Q-dependent, adiabatic electrons.The exponential form of canonical transformation (74) legitimates Taylor's expansion of the matrix elements ( 71) and ( 72) of base transformation matrix, that is, The form of transformation relations for boson operators of system Hamiltonian has to respect the factorized form of the total system wave function (62).It implies possibility of simultaneous, independent diagonalization of electron and boson subsystems.It means that transformed fermion and transformed boson operators obey not only standard anticommutation and commutation relations within the individual subsystems: but also transformed operators of both subsystems have to commute mutually like the original operators, that is, also the following commutation relations are required to hold:

Advances in Condensed Matter Physics
With respect to the fermion transformation relations (70), the form of transformation relations for boson operators is expressed as Also for matrix elements d(Q) of transformation matrix D(Q), the Taylor' expansion is defined as In order to ensure possibility of practical solution, in what follows, important restriction is imposed.The commutation relation ( 77) is required to hold up to quadratic terms in Taylor's expansions.It enables to express transformation coefficients d(Q) RS through coefficients c(Q) RS (see [66]-Appendix A).It will be shown that c r PR covers the strength of electron-vibration (phonon) coupling up to the first order of Taylor's expansion and determine also adiabatic correction to the electronic energy of the ground state ΔE 0 (ad) .The adiabatic transformation preserves total number of electrons, and nuclear coordinate operator is invariant under the transformation, that is, N e = P a + P a P = P a + P a P = N e , Up to the first order of Taylor's expansion, the crudeadiabatic momentum operator P r = (b r − b + ȓ ) is transformed as 81) is nuclear momentum operator on adiabatic level.
For adiabatic Q-dependent spinorbitals |ϕ P (x, Q) , which are the basis functions of the adiabatic Q-dependent electronic wave function of the ground state Φ 0 (r, Q), expressed over crude-adiabatic orbitals can be derived: As it is seen from ( 82), adiabatic wave function is modulated by the instantaneous nuclear coordinates {Q r = Q} of particular vibration (phonon) modes {r} with the weight proportional to transformation coefficients c r PR (coefficients of transformation matrix in the first order of Taylor's expansion-c r PR = ∂c PR (Q)/∂Q r ).At solution of the problem on adiabatic (Q-dependent) level, we have restricted ourselves to study total system in its electronic ground state with (62) representing wave function Ψ 0 (r, R) = χ 0,ν (R)Φ 0 (r, R).The Schrödinger equation of the total system (6) for electronic ground state Φ 0 (r, R) is then of the diagonal form (9). Seemingly we have lost the effect of electron-nuclear coupling through nuclear kinetic energy operator which is covered by Λ terms on the rhs of equation ( 6).The off-diagonal Λ terms are absent and from the diagonal Λ terms, the only non-zero element is the mean value of nuclear kinetic energy in the electronic ground state which is expected to be negligibly small for systems in adiabatic state.Nonetheless, it will be shown that this term (B 00 ) covers the same effect of e-p coupling as it is routinely calculated in solid-state physics by perturbation theory with e-p coupling Hamiltonian (22) when electronic excited states are approximated by promotion of electrons to virtual orbitals of the electronic ground state.
As it has been mentioned in Section 1, study of band structure of superconductors indicates that e-p coupling induces fluctuation of some band through Fermi level.At the moment when ACP of such a band approaches Fermi level, there is considerable reduction of electron kinetic energy, which for antiadiabatic state results even for dominance of nuclear dynamics (ω E F ) in some region of kspace.Electrons at these circumstances are not able to follow nuclear motion adiabatically.It means that electronic wave function, in order to respect this fact, should be dependent not only on instantaneous nuclear coordinates Q but it should also be an explicit function of the instantaneous nuclear momenta P, that is, Φ 0 ≡ Φ 0 (r, Q, P).
Let us assume that wave function of total system can be found in the following factorized form: The form of the wave function ( 83) is basically P-dependent modification of the original Q-dependent BOA (3).Like in adiabatic case, solution of the problem will be restricted to electronic ground state, that is, for total system, we have It means that effect of nuclear momenta will be covered only in the form of QP-dependent diagonal correction Λ 00 (Q, P) = Φ 0 (r, Q, P)| T N |Φ 0 (r, Q, P) , that is, in a similar way as it has been covered the effect of instantaneous nuclear coordinates Q on the adiabatic level, that is, Advances in Condensed Matter Physics 15 Solution of this problem is similar to the transition from crude-adiabatic to adiabatic level as presented above.Now, however, the transition from adiabatic to antiadiabatic level is established.
Nonadiabatic, nuclear displacements, and momenta (QP)-dependent electronic wave function Φ 0 (r, Q, P) in (84) assume existence of complete orthonormal basis set {ϕ R (x, Q, P)}, that is, validity of the following relations: Electron creation and annihilation operators that correspond to the (QP)-dependent moving base are written as barless operators (a + , a).Boson operators related to the (QP)dependent moving base are denoted also as bar-less operators, Since nonadiabatic electrons remain fermions, the operators obey standard fermion anticommutation relations (it follows from canonical transformation, see (91)): In (87), shorthand notation is used, that is, a R ≡ a R (x, Q, P) and a + R ≡ a + R (x, Q, P).Since adiabatic Q-dependent moving base derived by adiabatic transformation is complete and orthonormal (63), then due to (85), the base transformation to nonadiabatic (QP)-dependent moving base can be established over the base transformation relation: For fermion operators of second quantization, it follows that Elements of the P-dependent transformation matrix C(P) are The P-dependent transformation matrix C(P) is also unitary matrix, that is, the relations hold: The form of transformation relations for boson operators of system Hamiltonian has to respect again the factorized form of the total system wave function (84).Also in this case, it expresses possibility of simultaneous, independent diagonalization of electron and boson subsystems.It means that transformed-nonadiabatic fermion and transformed nonadiabatic boson operators obey not only standard anticommutation and commutation relations within the individual subsystems, but also, like the original and adiabatic operators, transformed nonadiabatic operators of both subsystems have to commute mutually, that is, also commutation relations have to hold: With respect to the fermion transformation relations (89), the form of transformation relations for boson operators that fully respects conditions (92) is Again, in order to enable practical solution, the commutation relation ( 93) is required to hold up to quadratic terms in Taylor's expansions.This restriction enables to express transformation coefficients d RS trough coefficients c RS , see [66]-Appendix B.
It can be shown that also this transformation preserves the total number of particles, that is, (96) However, coordinate operator is transformed up to first order of Taylor's expansion, as Advances in Condensed Matter Physics that is, For nonadiabatic (QP)-dependent spinorbitals |ϕ P (x, Q, P) which are the basis functions of the nonadiabatic (QP)-dependent electronic wave function of the ground state Φ 0 (r, Q, P), expressed over crude-adiabatic orbitals can be derived: Nonadiabatic wave function (99) in contrast to adiabatic wave function ( 82) is modulated not only by the instantaneous nuclear coordinates {Q r }of particular vibration (phonon) modes {r} but modulation is also over corresponding instantaneous nuclear momenta {P ȓ }.The weight of momentum modulation is proportional to the P-dependent transformation coefficients c r PR .It represents first derivative of c PR matrix element with respect to nuclear momentum P r , c r PR = ∂ c PR (P)/∂P r , that is, coefficient of transformation matrix in the first order of Taylor's expansion.It will be shown that these coefficients reflect not only the strength of e-p coupling but mainly the magnitude of nonadiabaticity.For true nonadiabatic situation, that is, for antiadiabatic state |ε 0 I (R) − ε 0 A (R)| < Dω r , the weight of such P-modulated state can be significant.

Solution of Nonadiabatic Problem:
Corrections to Energy Terms (1) Transformations of System Hamiltonian Base functions transformations have incorporated dependence of electronic states on operators of nuclear motion and vice versa.It implies, before the system Hamiltonian transformations, necessity to rearrange starting crude-adiabatic Hamiltonian: into more convenient form.
Let us formally divide this Hamiltonian on two parts, H A and H B .
(A) The nuclear part H N , as we already know, is quantized on crude-adiabatic level as In general, this part can be written as the sum of nuclear kinetic and nuclear potential energy: whereas The standard, usually harmonic-quadratic part of the nuclear potential energy, E (2)  NN (Q) is corrected now by some, yet unknown, potential energy term V (2)  N (Q), that is, supposed to be also quadratic function of nuclear coordinate operators.Evidently, the V (2)  N (Q) term is absent on crude-adiabatic level.This correction originates from the interaction of vibrating nuclei with electrons on adiabatic Q-dependent level.In general, kinetic energy of vibration motion can also be corrected by some, yet unknown quadratic function of nuclear momenta operators-W (2)  N (P), that is, On the adiabatic Q-dependent level kinetic energy correction is negligibly small and it is neglected.It becomes important only when the BOA is broken, |E te 0 (R) − E te n (R eq )| ≤ Dω ν , that is, in the case when electrons due to increased effective mass are not able to follow nuclear motion adiabatically and electronic states are QP-dependent.It should be stressed that both corrections are absent on crude-adiabatic level, they have been introduced just with respect to subsequent adiabatic and nonadiabatic transformation of system Hamiltonian.
(B) The second part of divided Hamiltonian is It is evident that the division on the two parts is only formal The reason of it is mainly pragmatic with respect to transformations and final solution.In this way the total system Hamiltonian is divided on quasi-bosonic (nuclear part H B ) and quasifermionic (electronic H A ) Hamiltonians.
The system Hamiltonian (in the form (105)) is now prepared for canonical transformations.The frequently used form of canonical transformation, is applied in order to get diagonal, or "more diagonal", final form at least up to first order of commutation expansion, [H 0 , S] + H int = 0.This kind of transformation changes the form of the Hamiltonian, but leaves unchanged the system variables.
Base functions transformations, as presented in preceding parts, introduce new dynamical variables; starting from Unlike to (106), now the Hamiltonian is not transformed itself, it remains of the original form, but its variables (A ν ) are replaced by new variables (A ν ): The Hamiltonian written in new variables is Since at this transformation, there is no any requirement for fulfillment of condition like at (106), the transformation does not make Hamiltonian "more" diagonal.Usually it is more complex, but very often it discloses physical aspects of the problem that are not obvious from nontransformed form with the original variables of system Hamiltonian.
The system Hamiltonian (105) is now subjected to this type of transformation.
In the first step, transformation from crude-adiabatic to adiabatic quasiparticles is realized.
The adiabatic quasi-particle transformations, up to the second order of Taylor's expansion, generate terms H Ai ( j,k) and H Bi( j,k) whereas i, j, k in subscript i( j, k) stand now for the ith order of Taylor' series expansion in adiabatic representation, the jth order of the coordinate operator Q, and the kth order of momentum operator P within the given order i of Taylor's expansion; that is, i = ( j + k).
Up to the second-order expansion in Q, P and in secondorder expansion of the operators of quasiparticle transformation (c PQ , d rPQ ), the H A and H B parts are transformed in the following way: The particular form of the terms (110a)-(110e) and details of adiabatic transformation are in [66, Appendix A], [67].
In the next step, the adiabatic form (terms (110a)-(110e)) is transformed to nonadiabatic form, nonadiabatic quasiparticles are introduced (up to the 2nd order of Taylor' series expansion): The particular form of the terms (111a)-(111h) and details of nonadiabatic transformation are in [66, Appendix B], [67]. (

2) Nonadiabatic Solution-Corrections to the Crude-Adiabatic Energy Terms
Like in case of the Hamiltonian transformations, also details of nonadiabatic solution are published in [66,67].
The main results of the solution can be written in the form of corrections to particular crude-adiabatic terms: Advances in Condensed Matter Physics (A) Correction to the ground state electronic energy: correction to H (0) (R) of (47) it is correction to zero-particle term of system Hamiltonian, that is, correction to the ground state electronic energy (48) calculated in clamped-nuclei approximation at fixed nuclear configuration [R] due to the effect of nuclear dynamics (QP-dependence).For total electronic energy then holds: The correction ΔE 0 (R) is {I}-occupied orbitals, {A}-unoccupied orbitals ( 113) are coefficients of adiabatic transformation (70) and nonadiabatic transformation (89) in first order of Taylor's expansion in Q or P, respectively.Approximate solution (see [66, Appendix C], [67]) yields for these terms: The final form of the correction The Ω matrix is a symmetric matrix of the form From the final form of the correction (117), the two limiting situations can be extracted.
(i) True nonadiabatic or antiadiabatic state, (ε 0 In this case the correction is negative, ΔE 0 (nad) < 0, and represents stabilization contribution to the electronic ground state energy.It is evident from (117) that the groundstate energy stabilization is exclusively the effect of nuclear dynamics, which is expressed through the coefficients c r AI of P-dependent transformation.This contribution can be considerably large and reach the extreme negative value for left-hand side limit toward singular point in (117).Singular point itself is excluded [66, Appendix C].The correction is always negative for the extreme case of strong nonadiabatic limit, Dω r /|ε 0 However, the contribution in this case does not represent the largest possible negative value and it is equal to ΔE 0 2 /Dω r .(ii) The second limit corresponds to adiabatic level, when electronic state is only Q-dependent.In this case the coefficients c r AI of P-dependent transformation are absent in (118) and for the correction holds: In case of strong adiabaticity, (ε 0 A − ε 0 I ) Dω r , which is standard condition for the BOA we have In contrast to nonadiabatic level when the correction is negative, ΔE 0 (nad) < 0, on the adiabatic level the correction is whatever small, but always positive (119).
An important aspect has to be mentioned in relation to adiabatic level.Introduction of adiabatic quasiparticles, that is, Q-dependent transformation, in solution of total system problem has to reproduce the basic results of standard treatment of solid state physics (or molecular physics) at solution of electron-phonon interactions by means of perturbation theory (see, e.g., (21), (22)).It can be shown that the base transformation is in this respect equivalent to standard treatment.
The exact adiabatic correction to the electronic ground state is, according to (9), the term B 00 , that is, DBOC-the mean-value of nuclear kinetic energy in adiabatic electronic ground state: We have shown [73] that (121) is the same as the adiabatic correction (119) to the electronic energy of the ground state Φ 0 , calculated through the expansion coefficients c PQ (Q) of Q-dependent transformation, that is, Within the single Slater determinant representation of the ground electronic state, this relation is exact since the eigenfunctions-spinorbitals {ϕ Q (R)} of the Hartree-Fock equations ( 29) are orthonormal and form complete basis set, that is, closure property holds It means that both subsets, that is, occupied {ϕ I (r, R)} and unoccupied {ϕ A (r, R)} orbitals are included at calculation of (119), (121).As it is seen from ( 119), (122), electronic ground state energy correction is due to virtual transitions between occupied {ϕ I (r, R)} and unoccupied orbitals {ϕ A (r, R)} at nuclear vibration motion.It should be stressed that at these virtual transitions the system remains in its electronic ground state Φ 0 (r, R eq ).In this respect, even ΔE 0 (ad) represents exactly the DBOC, it covers basically "nonadiabatic" (off-diagonal) corrections in the sense as these are calculated in a second order of perturbation theory when excited electronic states are approximated through virtual orbitals optimized for electronic ground-state Φ 0 (r, R eq ) and perturbation is an electron-vibration coupling H ep , as it has been discussed at the beginning.It can be seen very clearly from the derived expression for correction to frequency of normal modes: 2  2 . (124) The final effect is softening of normal mode vibration frequency: The normal mode energy correction (124) is identical to the normal mode energy correction which can be derived as nonadiabatic correction in second order of perturbation theory assuming single boson excitation processes, that is, ( The term nonadiabatic is related to the fact that within the perturbation treatment contributions to the correction arise due to virtual transitions between occupied {I} and unoccupied {A} states, that is, correction is calculated through off-diagonal elements of the form For quasimomentum k, q-space representation of multiband solids, with respect to (k;−k) symmetry, the temperature-dependent form of the equation for correction to electronic ground-state energy (117) can be expressed in the following form: Summations in ( 127) are over all bands {ϕ R , ϕ S } and k points of the first BZ of multiband system, while In (128), the wave vector k corresponds to states fully occupied below Fermi level ( f k = 1), and wave vector k corresponds to empty-virtual states above the Fermilevel( f k = 0).It should be noticed that u k −k is matrix element of e-p coupling (43).(B) Correction to one-particle term: correction to H (1) (R) of ( 47).
The one-particle correction has been derived in the form It should be noticed that the second contribution in (129) originates from boson excitation dependent terms of nonadiabatic Hamiltonian and expression in the second line of (129) is for boson vacuum.The first term in (129) originates from those terms of nonadiabatic Hamiltonian which are not dependent on boson operators (pure fermionic terms).At derivation of (129), restriction has been to the first order Advances in Condensed Matter Physics of e-p coupling and the terms that are the products of e-p coupling and coulomb two-electron interactions (e.g., (ν 0 PRPA − ν 0 PRAP )(c r IR c r * IA + c r IR c r * IA )) have been expected to be negligibly small comparing to pure e-p coupling terms and were neglected.Details of derivation can be found in [66,67].
(i) Nonadiabatic polarons The diagonal form of the one-particle correction (129) is Substitution for transformation coefficients (115), (116) and rearrangement of the expression yields for e-p interaction part of the Hamiltonian: Transcription of (130b) to quasimomentum k,q-space of multiband solids has the form Expression (130c) represents total one-electron energy correction on the general Q, P−dependent level due to e-p interactions.The first term of (130c) is standard adiabatic (Q-dependent) polaron as it can be derived from Fröhlich Hamiltonian by the Lee-Low-Pines transformation [74].The second term of (130c) is the correction to polaron energy that arises due to dependence of electronic motion not only on nuclear coordinates but also on the nuclear momenta P (nonadiabatic modification of the BOA).This term can be interpreted as a correction to the energy of individual polaron by an effective field created by all other polarons of the system.(ii) Correction to orbital energies: gap opening in oneparticle spectrum of quasidegenerate states at Fermi level.
For quasicontinuum of states at Fermi level, which is characteristic for metal-like band structures, contribution of second term in (130a) can be expected to be negligibly small.Summation in this term runs over all states, occupied and unoccupied.Since the term is odd function of (ε 0 P − ε 0 R ), contributions from occupied and unoccupied states will tend to cancel mutually.On the other hand, character of one-particle spectrum can be changed at the e-p coupling through contribution of the first term in (130a) since summations in this term run separately over occupied and unoccupied states.Moreover, magnitude of the change in orbital positions is temperature dependent through the temperature dependence of the population of states.
At these circumstances, for investigation of possible changes in the character of one-electron spectrum of system due to e-p interactions on Q,P-dependent level, the first term of (130a) is crucial.For correction to orbital energy Δε P of general state ε 0 P , from the first term in (130a), it follows that Finally, corrected orbital energy is Let us consider only the couple of quasidegenerate states at Fermi level, occupied state ε 0 I , and unoccupied state ε 0 A ,a situation that can characterize couple of states in antiadiabatic state, 131) and from the structure of Ω matrix, (118) follows directly: It means that orbital energy of unoccupied state has been increased, ε A > ε 0 A , and orbital energy of the occupied state has been decreased, ε I < ε 0 I .It can be generalized for a set of quasidegenerate occupied {J} and unoccupied {B} states (quasicontinuum of states) at Fermi level.Then, with respect to the fact that for antiadiabatic state correction to the ground electronic state (117) is negative AI Ω AI < 0, the following relations can be derived: In particular, At finite temperature T, for a correction Δε P to an arbitrary state ε 0 P , from the set of quasidegenerate occupied and unoccupied states at Fermi level can be written: The occupation factor f Q obeys the Fermi-Dirac statistics: It is evident that for temperature 0 K, (135) is reduced to (131).From (135), temperature dependence of energy gap which is open in one-particle spectrum at Fermi level can be derived in the form The gap is defined as the energy difference of lowest lying corrected unoccupied state ε B(LUMO) and highest lying corrected occupied state ε J(HOMO), , that is, as an energy gap in one-particle spectrum.At temperature 0 K holds: Factor 4 in the denominator of the argument of tgh in (137) follows from the assumption that at Fermi level density of quasidegenerate occupied {ε 0 J } and unoccupied {ε 0 B } states of the band with gap opening is the same and consequently ε B(LUMO) = |ε J(HOMO) |.This factor can be larger or smaller than 4, depending on the actual difference in the density of states of occupied and unoccupied states of band at FL.
It is trivial to show that the corrections to orbital energies are negligibly small (basically zero) for a system in adiabatic state when (ε 0 A − ε 0 I ) ω. (C) Correction to two-particle term -correction to H (2) (R) of (47).
Correction to H (2) (R) represents correction to electron correlation energy ( see Note in Section 1) due to e-p interactions, in particular, it is dependence of electronic motion on nuclear vibration displacements Q and momenta P. It has been derived in the form [66,67] Substituting for transformation coefficients (115), ( 116) yields Transcription of (140) to quasimomentum k,q-space representation of multiband solids has the form

The Effect of e-p Coupling on Band Structure of Superconductors and Nonsuperconducting Analogues.
In Figures 1(a), 1(c), 1(e), 1(g), (the pictures on the left), there are band structures (BS) of four superconducting compounds of different types and crystal structures at optimized equilibrium geometries.The details of calculation are published elsewhere [64,75,76].

Advances in Condensed Matter Physics
As it can be seen from Figure 1, all the band structures are of adiabatic metal-like character with a relatively low density of states at the FL (indicated by a dashed line).Coupling to the respective phonon mode(s) in particular compounds seemingly does not change the metal-like character of BS.In all cases, however, e-p coupling induces BS fluctuation (see the pictures on the right), which is characteristic by fluctuation of the analytic critical point (ACP) of some band across the FL (cf.a-b, c-d, e-f, g-h).
In particular, for MgB 2 (hP3, P6mmm, #191; AlB 2 ωtype) coupling to the E 2g phonon mode (in-plane stretching vibration of B-B) results in splitting of σ bands (p x , p y electrons of B atoms in a-b plane) in Γ point of the first Brillouin zone (BZ)-Figure 1(b).Related to band topology, the analytic critical point (ACP, maximum) of σ bands is located at Γ point and, for the displacement ≈ 0.016 A • /Batom out of equilibrium position, the ACP crosses FL.This means periodic fluctuation of the BS between topologies 1a ↔ 1b in coupling to vibration in the E 2g mode.
The situation is similar for YB 6 (cP7, Pm3m, #221; CaB 6 type).In this case, BS fluctuation is related to the T 2g mode (valence vibration of B atoms in basal a-b plane of Boctahedron).At the displacement ≈0.017A • /B-atom out of equilibrium position, the ACP (saddle point) of the band with dominance of B-p and Y-d electrons crosses FL in the M point and the BS fluctuates between topologies 1c ↔1d.
As it can be seen from Figure 1 in all cases a particular ACP of the band crosses FL at a displacement which is smaller than root-mean square (rms) displacement for zero-point vibration energy in respective phonon mode.This means, however, that in vibration where ACP approaches FL at the distance less than ±ω, the Fermi energy E F (chemical potential μ) of the particular band-electrons close to the point where ACP crosses FL is smaller than the vibration energy of corresponding phonon mode, E F < ω 2g .At these circumstances the adiabatic BOA is not valid and standard adiabatic theories cannot be applied.Moreover, shift of the ACP increases substantially the density of states (DOS) at FL, n σ (E F ) = (∂ε 0 σ /∂k) −1 EF , and induces the corresponding decrease in the effective electron velocity (∂ε 0 σ /∂k) EF of fluctuating band in this region of k-space.Under these circumstances, the system is in the intrinsic nonadiabatic state, or even in the antiadiabatic state, E F ω 2g , and electronic motion depends not only on nuclear coordinates Q but is strongly influenced by nuclear dynamics-momenta P. The main part of the effect of nuclear kinetic energy on electronic motion can be derived as diagonal correction by sequential Q, P-dependent base transformations as it has been presented in preceding parts.
Instability of the electronic structure at e-p coupling is absent in respective non-superconducting analogues, such as XB 2 (X≡Al, Sc, Y, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, W, Mn, etc.), CaB 6 , YBa 2 Cu 3 O 6 , and Mg.Even in the case of the XB 2 , coupling to the E 2g mode induces splitting of σ bands in the Γ point, the systems remain stable in the adiabatic state.For these systems, the ACP of σ band does not fluctuate across FL.As an illustration, in Figure 2 are band structures of corresponding nonsuperconducting analogues at equilibrium high-symmetry structure (Figures 2(a), 2(c), 2(e), and 2(g)) and at distorted geometry (Figures 2(b), 2(d), 2(f), and 2(h)) with the same displacements in respective phonon modes as those in the case of corresponding superconductors at the transition to the antiadiabatic state.
In the case of AlB 2 , in spite of σ bands splitting and nearly the same value of the e-p interaction strength (the calculated mean value is u ≈ 1.01 eV/u.cell) as that of MgB 2 (u ≈ 0.98 eV/u.cell),AlB 2 remains in e-p coupling in the adiabatic state as a nonsuperconducting compound.In this case, BS fluctuation (bands splitting in EP coupling) does not decrease chemical potential, it remains in e-p coupling still larger than the vibration energy (μ ad > Dω) and, consequently, there is no driving force for transition to the antiadiabatic state.
The CaB 6 is an insulator, and coupling to the T 2g mode (valence vibration of B atoms in basal a-b plane of Boctahedron) does change the character of the BS (Figures 2(c) and 2(d)).
In the case of deoxygenated YBCO, YBa 2 Cu 3 O 6 , in contrast to the superconducting YBa 2 Cu 3 O 7 , a combination of electron coupling to A g , B 2g , and B 3g phonon modes leaves band structure without substantial change (Figures 2(e) and 2(f)).In the case of YBa 2 Cu 3 O 7 , the ACP (saddle point) at Y point fluctuates across FL (see Figures 1(e) and 1(f)), which yields substantial reduction of chemical potential → μ antiad < Dω.For YBa 2 Cu 3 O 6 the ACP does not fluctuate across FL and chemical potential remains larger than the phonon energy spectrum, μ ad > Dω, and the system remains in the adiabatic state.
Also the band structure of Mg remains without changes at e-p coupling (Figures 2(g) and 2(h)).

Nonadiabatic e-p Coupling Effects in
the Antiadiabatic State

Formation of Antiadiabatic Ground State and Gap
Opening.On crude-adiabatic level, total ground state electronic energy E te 0 (R eq ) has minimum at R eq , (dE te 0 /dR) Req = 0, and within the HF-SCF approximation it is equal to (27).The only correction to this energy is electron correlation energy (53), or (54) in second order of PT, which is negative and contributes to stabilization of the ground state at equilibrium geometry R eq : Related to any phonon mode, nuclear displacements at vibration motion increase total electronic energy (potential energy of nuclear motion in particular phonon mode).For displaced geometry R d on crude-adiabatic level holds: Since two-electron coulomb interactions {ν 0 PQRS } do not depend explicitly on nuclear variables it can be expected Advances in Condensed Matter Physics  that electron correlation energy at vibration displacement R d has not been changed significantly, that is, (E corr ) Req ≈ (E corr ) Rd .On crude-adiabatic level, for an increase of the total electronic energy ΔE d due to nuclear displacement R d then holds: In principle, two situations can occur; nuclear displacements related to some phonon mode(s) induce forma-tion of antiadiabatic state, or system remains in adiabatic state with respect to vibration motion in all phonon modes.The band structures in Figure 1 indicate possibility of antiadiabatic state formation for superconductors.The question if antiadiabatic state can be a stable state, that is, to be a ground electronic state of system at distorted nuclear configuration R d , depends on the value of the ground state energy correction ΔE 0 (R d ) (117).Since for antiadiabatic state this correction is negative, holds, the electronic state of the system is stabilized at distorted geometry R d .The reason of it is significant participation of the nuclear kinetic energy term expressed through contribution of the coefficients c r AI of P-dependent transformation (see structure of (113)), which stabilizes fermionic ground state energy in antiadiabatic state at distorted nuclear configuration R d .
Stabilization (condensation) energy at transition from adiabatic into antiadiabatic state is Correction (127) to the electronic ground-state energy in the k-space representation due to interaction of pair of states mediated by the phonon mode r can be rewritten as In general, all bands of 1st BZ of a multiband system are covered, including intraband terms, that is, stands for matrix element of e-p coupling and n εk , n εk are DOS of interacting bands at ε 0 k and ε 0 k .For adiabatic systems, such as metals, this correction is positive and negligibly small (DBOC).As it has been mentioned, only for systems in the antiadiabatic state the correction is negative and its absolute value depends on the magnitudes of u r k−k and n εk , n εk at displacement for FL crossing.At the moment when ACP approach FL, the system not only undergoes transition to the antiadiabatic state but also DOS, n σ (E F ) = (∂ε 0 σ /∂k) −1 EF , of the fluctuating band is considerably increased at FL.It is a situation invoking possibility of van Hove singularity formation.Important role of van Hove singularity formation for high-T c superconductors has been proposed by Bok and his collaborators [77,78].
For all the presented systems at 0 K, ΔE 0 (na) which covers the effect of nuclear momenta prevails in absolute value the electronic energy increase ΔE d (R dcr ), that is, (145) is fulfilled at nuclear displacements d cr when ACP crosses FL as calculated for clumped nuclear adiabatic structures.
The highest value of the "condensation energy" (146), E 0 cond ≈ −34 meV/unit cell, has been calculated [76] for YBa 2 Cu 3 O 7 .In the case of MgB 2 , E 0 cond ≈ −13 meV/unit cell has been obtained.For YB 6 , and Zn the values of E 0 cond are smaller but in both cases are also negative.
Under these circumstances, each of the studied systems is stabilized in the antiadiabatic electronic ground state at broken symmetry with respect to the adiabatic equilibrium high-symmetry structure.In my opinion, it can be identified experimentally by ARPES as a kink on momentum distribution curve at FL in form of the band curvature at the ACP when it approaches FL (see Figure 3(a), which presents calculated results for YBa 2 Cu 3 O 7 ).The crucial influence on the final value of the correction ΔE 0 (na) (R d ) is related to an increase of the DOS (n εk , n εk in (147)) connected to the ACP when it approaches (±ω) FL at transition into antiadiabatic state.Situation for YBa 2 Cu 3 O 7 is presented in Figure 3(b).In the ARPES spectra of YBa 2 Cu 3 O 7 a kink on momentum distribution curve in off-nodal direction (Γ-Y line, near to Y point) should appear in low-energy region at FL for temperatures below T c and it will be shifted away from FL, or will disappear above T c .It is related to fluctuation of the ACP (inflex point) of one of CuO 2 layer (d x 2 −y 2 − pσ) band when it approaches FL at coupling to mentioned phonon modes.In my opinion, it is an effect of dispersion renormalization as it has been discussed in Section 1.
Nontrivial and very important property of system which is stabilized in the antiadiabatic state should be stressed.Due to translation symmetry of the lattice, the created antiadiabatic electronic ground state is geometrically degenerate in distorted geometry (R d,cr ) with fluxional nuclear configuration in particular phonon mode(s) (see, e.g., [68,Figure 5] for case of MgB 2 ).It means that the ground-state energy is the same for different positions of the involved atoms (in phonon modes which drive the system into this state).Expressed in the other way, the involved atoms can circulate over perimeters of the circles with the radii equal to displacements d cr at FL crossing, without the energy dissipation.
In transition to the antiadiabatic state, k-dependent gap Δ k (T) in quasicontinuum of adiabatic one-electron spectrum is opened.The gap opening is related to shift Δε Pk of the original adiabatic orbital energies ε 0 Pk , ε Pk = ε 0 Pk + Δε Pk , and to the k-dependent change of DOS of particular band(s) at Fermi level.Shift of orbital energies in band ϕ P (k) has the form (134a), (134b).In quasimomentum k-space of multiband solids, these equations are for k > k F , and Replacement of discrete summation by integration, . It is of crucial importance in relation to fluctuating band (see Figure 3(b)).For corrected DOS n(ε k ), which is the consequence of shift Δε k of orbital energies, the following relation can be derived: The term n 0 (ε 0 k ) stands for uncorrected DOS of the original adiabatic states of particular band: Close to the k-point where the original band which interacts with fluctuating band intersects FL, the occupied states near FL are shifted downward below FL and unoccupied states are shifted upward, above FL.The gap is identified as the energy distance between created peaks in corrected DOS above FL (half-gap) and below FL.The formation of peaks is related to the spectral weight transfer, that is, observed by ARPES or tunneling spectroscopy in cooling below T c .For some of the studied compounds, the calculated corrected DOS (149) of particular band(s) with gap opening are in Figure 4.
In particular, YBa 2 Cu 3 O 7 exhibits an asymmetric gap in two directions: for Cu-O1 chain-derived (d-pσ) band the gap in one-particle spectrum is Δ b (0) ≈ 35.7 meV in the Γ-Y direction (4a) and Δ a (0) ≈ 24.2 meV in the Γ-X direction (4b), respectively.The calculated asymmetry, that is, the ratio (Δ a (0)/Δ b (0)) theor ≈ 0.68, or (Δ b (0)/Δ a (0)) theor ≈1.47 is very close to the experimental values, (Δ a /Δ b ≈ 0.66 and Δ b /Δ a ≈ 1.5) which have been recorded [79,80] for untwined YBa 2 Cu 3 O 7 .Inspection of the topology (Figures 1(e)-1(f)) of the Cu-O 1 chain-derived (d-pσ) band suggests that one can expect gap opening in one-particle spectrum of this band not only in Γ-Y(b) and Γ-X (a) direction, but it can be expected to be opened also in nodal, that is, Γ-S direction.Based on k-dependence of DOS of this band in particular directions at FS, it can be expected that Δ Γ−S ≤ Δ a(Γ−X) , that is, top of the dispersion should appear ≤ 15 meV below FS.One can expect opening of very small gap for this band also in Γ-T(U) direction [76].However, due to relatively large distance of CuO 2 layers (≈ 8.3 A • ), e-p coupling which mediates interactions between electrons of these layers is small, which suggests that gap in one-particle spectrum of CuO 2 layer-based (d x 2 −y 2 − pσ) bands is not expected to be opened.
The corrections to orbital energies (148a), (148b) and to the ground state energy (147) are temperature-dependent and decrease with increasing T. At a critical value T c , the gap in one-particle spectrum (137) which is formed at 0 K disappears, Δ(T c ) = 0 (continuum of states is established at FL), while |ΔE 0 (na) (T c )| ≤ ΔE d,cr and the system undergoes transition from the antiadiabatic into adiabatic state, which is stable for equilibrium high-symmetry structure above T c .With respect to Δ(0), from (137) a simple approximate relation follows for critical temperature T c of the adiabatic ↔ antiadiabatic state transition: The calculated values of T c for presented compounds are YBa 2 Cu 3 O 7 , T c ≈ 92.8 K, MgB 2 , T c ≈ 39.5 K, and YB 6, T c ≈ 11 K.It is in a good agreement with the experimental values of T c for superconducting state transition of these compounds [4,60,83].

Thermodynamic Properties of System in the Antiadiabatic
State.Presented results naturally provoke the question if the antiadiabatic state is related in some way to superconducting state.Crucial in this respect are thermodynamic properties of antiadiabatic state, that is, at temperatures T < T c .

(I) Character of condensation in antiadiabatic state
As it has already been mentioned, in the antiadiabatic state, ground-state total electronic energy of system is geometrically degenerate.Distorted nuclear structure, related to couple of nuclei in the phonon mode r that induces transition into antiadiabatic state, has fluxional character.There exist an infinite number of different, distorted configurations of this couple of nuclei in the phonon mode r and all these configurations, due to translation symmetry of the lattice, have the same ground state energy.Position of the involved displaced couple of nuclei is on the perimeter of circles with the centers at R eq (equilibrium on crude-adiabatic level) and with radii equal to ΔR = |R eq − R d,cr |.The R d,cr is distorted geometry at which ACP approaches FL and system undergoes transition from adiabatic into antiadiabatic state.
Due to the geometric degeneracy of the ground state energy, the involved atoms can circulate over perimeters of the circles without the energy dissipation.

Advances in Condensed Matter Physics
The dissipation-less motion of the couple of nuclei implies, however, that e-p coupling of involved phonon mode and electrons of corresponding band has to be zero.
Let as shown that in the antiadiabatic state this aspect is fulfilled.
The effective e-p interactions which cover the Q, Pdependence has the form that is for boson vacuum represented by the second term in (130a) is After substitution for transformation coefficients c r PR , c r PR (115), (116) and algebraic manipulation, the quasimomentum form of this term is It means that for electrons that satisfy condition of extreme nonadiabaticity with respect to interacting phonon mode r in particular direction of reciprocal lattice where the gap in one-electron spectrum has been opened, the electron (nonadiabatic polaron)-renormalized phonon interaction energy equals zero.Expressed explicitly, in the presence of external electric potential, dissipation-less motion of relevant valence band electrons (holes) on the lattice scale can be induced at the Fermi level (electric resistance ρ = 0), while motion of nuclei remains bound to fluxional revolution over distorted, energetically equivalent, configurations.The electrons move in a form of itinerant bipolarons, that is, as a polarized cloud of intersite charge density distribution (see Figure 5).Due to temperature increase, thermal excitations of valence band electrons to conduction band induce sudden transition from the antiadiabatic state into adiabatic state at T = T c , that is, (145) does hold.For temperatures ) holds and the system becomes stable at equilibrium R eq as it is characteristic for adiabatic structure.
In the adiabatic state, properties of the electrons are in sharp contrast with the properties of electrons in antiadiabatic state.The electrons in this case are in a valence band more or less, tightly bound to respective nuclei at adiabatic equilibrium positions and theirs motion in conducting band is restricted by scattering with interacting phonon modes.It corresponds to situation at T > T c .
For extreme adiabatic limit Dω q /|ε 0 k − ε 0 k−q | → 0, from (153a) for electron-phonon interaction energy in this case it follows that Expression (153c) represents basically energy of standard adiabatic polarons (small, self-trapped) that contributes to the total energy of system.II.Electronic specific heat and entropy in antiadiabatic state.
It has been shown that at formation of the antiadiabatic ground-state, electronic energy is decreased and for involved band(s) the gap in one-particle spectrum has been opened (shift of orbital energies).This fact has to be reflected by change of related thermodynamic properties.In particular, for electronic specific heat, can be derived: At derivation of (155), it has been assumed that for the band with gap opening, uncorrected density of quasicontinuum of occupied and unoccupied states at Fermi level is the same.Expression (155) can be derived either from the ground-state energy correction (127), (147), or from the basic statistical relation for entropy: In (156), f k is the Fermi-Dirac occupation factor of state ε k = ε 0 k + Δε k .Then for entropy related to formation of antiadiabatic state can be derived: For temperature derivative of entropy (157a), it follows that It yields directly to final expression (155) for C V ,el (T).
From the (155), in the limit T 0 K → , the exponential behavior characteristic for superconductors can be derived: The density of states at Fermi level n(ε kF ) in the above equations represents mean value of corrected density of states close to the k-point where the peak in density of states has been formed.From practical reasons, it can be approximated by the mean value of density of states of the fluctuating band in ACP at the moment when it approaches (±ω/2) Fermi level (antiadiabatic state).In the antiadiabatic state, density of states at Fermi level is considerably increased since density of states of fluctuating band at ACP is usually high (possibility of van Hove singularity formation at Fermi level).III.Magnetic properties of system in the antiadiabatic statecritical magnetic field System in superconducting state can exhibit absolute diamagnetism and Meissner effect only if inside the system B = 0.In this case, there has to exist some critical value of external magnetic field H c which destroys superconducting state and induces transition of the system into normal state (characteristic by finite-nonzero value of electric resistance ρ / = 0 at finite-nonzero density of electric current j), like it occurs in case of temperature increase above T c .It also means that at critical temperature and above it, T ≥ T c , critical magnetic field has to be zero, H c (T ≥ T c ) = 0.
It can be shown that antiadiabatic state exhibits this property.
From thermodynamics for critical magnetic field in this case, it follows that In (159), F (ad) and F (na) stand for free energies of the system in adiabatic and nonadiabatic (antiadiabatic) state.Then, From ( 154), (155), it follows that After substitution of (160b), (157a) into (160a) and algebraic rearrangements, for critical magnetic field at finite temperature T, it follows that At temperature 0 K, for critical magnetic field results, The relation between critical magnetic fields at finite and zero temperature follows from (160c), (160d): At critical temperature T c , when due to thermal excitations the antiadiabatic state is suddenly changed into adiabatic state and above this temperature, critical magnetic field has to be zero, H c (T ≥ T c ) = 0. Since there is no gap in one-electron spectrum in metal-like adiabatic state, Δ(T ≥ T c ) = 0, then zero value of critical magnetic field follows directly from (160c), (161).Derived equations show that system in the antiadiabatic state, beside zero value of electric resistance ρ = 0 (dissipation-less motion of bipolarons), has also specific property that is, necessary for occurrence of the Meissner effect.

Electron Correlation, Cooper's Pairs and Bipolarons
Correction ΔH ep to the two-particle term H (2) (R) represents correction to the electron correlation ( see Note in Section 1) energy due to e-p coupling on QP-dependent nonadiabatic level and has the form (141).The electron correlation energy on crude-adiabatic level H (2) (R) is a perturbation contribution to the electronic ground-state energy H (0) (R 0 ) = E 0 SCF (R 0 )-, (48).It is an improvement of e-e interaction contribution beyond the level of HF method which treats electrons with α and β spins in an unbalanced way.The electron correlation energy is negative (see Note in Section 1) and contributes to stabilization of ground state electronic energy as a perturbation contribution, which in the second order of perturbation theory has the form ( 54), E corr (R 0 ).
Let us see if the correction to correlation energy due to e-p coupling ΔH ep represents stabilization or destabilization contribution [66,67] 141) are positive and negative value of the matrix elements of this two-particle correction is reached for a reduced form if nominators are negative, that is, Since q / = 0, it can be reached if q = −k − k .At these circumstances for reduced form which maximizes e-e attraction, it can be derived that that is, for system which is in antiadiabatic state.
In case of strong antiadiabatic regime 0, for correction to electron correlation energy, it follows that The second extreme case is for strong adiabatic regime Dω k −k /|ε 0 k − ε 0 k | → 0. At these circumstances, correction to electron correlation energy approaches zero value: General adiabatic expression on Q-dependent adiabatic level that is, the adiabatic form of this correction can be derived directly from (139) in the limit c = 0, The reduced form is Like for nonadiabatic Q, P-dependent case also for adiabatic Q-dependent level, the correction to electron correlation energy is negative.By comparing (162) with (163a) it can be seen that nonadiabatic correction to electron correlation energy is in absolute value larger than corresponding adiabatic correction.Moreover, for antiadiabatic state, when the system is close to singular point in (162), the nonadiabatic correction can be significant.
Along with crude-adiabatic correlation energy (54), also the correction to electron correlation energy (162) contributes to stabilization of system in antiadiabatic state, which is the ground state of the system at distorted nuclear geometry R d,cr .It should be reminded, however, that regardless of the correction to electron correlation energy (which represents contributions of biexcited configurations in second and higher orders of perturbation theory), the system has already been stabilized in antiadiabatic state at distorted geometry due to correction to the ground state electronic energy (128) which represents zero-order correction in terms of perturbation theory.In this respect, increased electron correlation is the consequence of e-p interactions which have induced transition of the system from adiabatic into antiadiabatic state and stabilized it in this state no matter if correction to electron correlation energy is accounted for or not.
At finite temperature, the product of Fermi-Dirac occupation factors has to be introduced into derived equations: kk . . .With increasing temperature from 0 K, the value of the correction (162) which is characteristic for antiadiabatic state decreases and at T = T c when system undergoes sudden transition from antiadiabatic state at distorted geometry R d,cr into adiabatic state at undistorted geometry R 0 = R eq and above this temperature, for correction to electron correlation energy holds corresponding temperature-dependent adiabatic form (165). Above T c , this correction along with crude-adiabatic correlation energy (54) stabilizes adiabatic ground state of system at undistorted-adiabatic equilibrium geometry R 0 = R eq .
Let us turn attention to the Fröhlich effective Hamiltonian of e-e interactions [7]: The reduced for of this Hamiltonian is This interaction term is either attractive or repulsive depending on the sign of denominator.For antiadiabatic represents effective attractive electron-electron interactions.It should be reminded that effective attractive e-e interactions are the crucial condition of Cooper's pair formation and the basis of the BCS theory.
In the limit of extreme antiadiabaticity It is immediately seen that this interaction, in contrast to (165), is always repulsive and (k, k ) symmetry is disturbed.With respect to the fact that effective attractive ee interaction is a pure correction to electron correlation energy, the Cooper's pair idea within the BCS treatment is without any doubts formally correct, nevertheless, it is rather a model treatment how to solve the instability problem with severe restriction to fixed nuclear framework, that is, within the crude-adiabatic BOA which requires for total system validity of the adiabatic condition Discussion of these and related aspects, the BCS treatment of superconducting state transition and theory of antiadiabatic state formation have been analyzed in [68].
From the theory of the antiadiabatic ground state, instead of Cooper's pairs, formation of mobile bipolarons results in a natural way as a consequence of translation symmetry breakdown on antiadiabatic level.The bipolarons arise as polarized intersite charge density distribution that can move over lattice without dissipation due to geometric degeneracy (fluxional structure) of the antiadiabatic ground state at distorted nuclear configurations.Formation of polarized intersite charge density distribution at transition from adiabatic into antiadiabatic state is reflected by corresponding change of the wave function.For spinorbital (band) ϕ R(k) relation (99) holds, that is, Dω r , coefficients c r RS of Q-dependent transformation matrix become negligibly small and absolutely dominant for modulation of crude-adiabatic wave function are in this case coefficients c r RS of P-dependent transformation matrix.For simplicity, consider that transition into antiadiabatic state is driven by coupling to a phonon mode r with stretching vibration of two atoms (e.g., B-B in E 2g mode of MgB 2 , valence T 2g mode vibration of B-B atoms in basal a-b plane of B-octahedron in YB 6 , or vibration motion of O2, O3 in Cu-O planes, B 2g , B 3g modes of YBCO).Let m 1 and m 2 are equilibrium site positions of involved nuclei on crudeadiabatic level and d 1 and d 2 are nuclear displacements at which crossing into antiadiabatic state occurs.At these circumstances, the original crude-adiabatic wave function ϕ 0 k (x, 0, 0) is changed in the following way: In (171), site approximation for momentum has been used, that is, P q ∝ (sign •q) m P m e iq•m .In antiadiabatic state, for particular k and proper q values, nonadiabatic prefactors under summation symbol in (171) can be large.The prefactors, that is, coefficients of P-dependent transformation matrix, reflect influence of nuclear kinetic energy on electronic structure.At the dominance of these contributions (antiadiabatic state), strong increase in localization of charge density appears at distorted site-positions for x equal to (m 1 −d 1 ) and (m 2 + d 2 ).It induces (or increases) intersite polarization of charge density distribution.As an example, calculated iso-density line for highest electron density in MgB 2 is presented in Figure 5.For equilibrium high-symetry structure (R eq ), the highest electron density is localized at equilibrium position of B atoms (Figure 5(a)).For distorted nuclear geometry (R d,cr ) in the E 2g mode, electron density is polarized and the highest value is shifted into the intersite positions, bipolarons are formed (Figure 5(b)).This kind of intersites polarization persists, and on a lattice scale, it has itinerant character at nuclear revolution over perimeters of the fluxional circles with the radius R d,cr , until the system remains in the antiadiabatic state.

Discussion and Conclusion
Experimental results are always crucial for any theory which aims to formulate basic physics behind observed phenomenon or property.However, an experiment always cover much wider variety of different influences which have impact on results of experimental observation than any  theory can account for, mainly if theory is formulated on microscopic level and some unnecessary approximations and assumptions are usually incorporated.On the other hand, interpretation of many experimental results is based on particular theoretical model.This is also the case of ARPES experiments at reconstruction of Fermi surface for electronic structure determination of high-T c cuprates.Interpretation of experimental results is based on band structure calculated for particular compound.Methods of band structure calculations are always approximate, with different level of sophistication.Calculated band structure, mainly its topology at FL, is a kind of reference frame for assignment of particular dispersion of energy distribution curve (EDC) or momentum distribution curve (MDC) to particular band of studied compound at interpretation of ARPES.This is in direct relation with theoretical understanding of crucial aspects of SC-state transition and with mechanism of superconductivity in general.
In case of high-T c cuprates, the crucial is considered to be the order parameter, that is, symmetry of superconducting gap which should be of (d x 2 −y 2 ) symmetry.Within a simple orbital model consideration of CuO 2 plane in square lattice configuration, this type of symmetry is generally accepted to simulate a-b plane superconductivity of cuprates in a best way.
In this sense, unexpected results have been published recently [30].The high-resolved (sub-meV) ARPES with VUV laser as a light source for optimally doped untwined YBCO reveals that in superconducting state, gap is opened in one-particle spectrum.The momentum dependence is very interesting, however.Like for other cuprates, strong band renormalization (kink formation at ≈60 meV) in offnodal direction has been observed.The gap is opened in offnodal (Γ-X and Γ-Y) directions with top of the dispersion at ≈ 20 meV below FL.In contrast to [79,80], the results of [30] do not indicate presence of a/b asymmetry.What is surprising, however, is the finding that gap remains finite also in nodal Γ-S direction, with top of dispersion at ≈ 12 meV below FL.Presence of finite "nodeless" gap (i.e., in nodal direction) sharply contradicts to commonly accepted idea of cuprates as superconductors of (d x 2 −y 2 ) symmetry.In normal state (100 K), the gap is closed in both, off-nodal and nodal directions.
The order parameter of (d x 2 −y 2 ) symmetry assumes gap opening in a band with dominant contribution of Cubased (d x 2 −y 2 ) orbital.In layered cuprates, it corresponds to (d x 2 −y 2 − pσ) bands of CuO 2 planes.Respected band structure calculation for YBCO is the DFT-based LDA published by Andersen et al. [84] to which the authors [30] refer.Inspection of BS in [84] reveals, however, that Cu2-O2O3 planes bands do not intersect FL in off-nodal Γ-X nor in Γ-Y direction, but the bands intersect FL in nodal Γ-S direction and in S-Y and S-X directions.It is an indication that experimentally detected gap [30], which is opened in off-nodal (Γ-X and Γ-Y) directions and also in nodal Γ-S direction, does not correspond probably to bands of CuO 2 planes.
In spite that there is a general agreement among different DFT-based band structure calculations in overall character of band structure for particular cuprate, there are small, but important differences in details concerning topology of some bands at FL.In particular, for YBCO, DFT-based all-electron band structure calculation [85] using the FLAPW method yields important differences in topology of Cu-O chainderived (d − pσ) band.In contrast to band structure [84], where the band intersects FL in S-X, Γ-Y and Γ-S directions, the BS [85] yields for Cu-O chain band intersections with FL in Γ-X, Γ-Y and Γ-S directions, respectively.Calculated Cu-O chain band topology [85] is in an excellent agreement with experimental electron-positron momentum density in optimally doped untwined YBCO detected [86] by 2Dangular correlation of electron-positron annihilation radiation (ACAR) technique.It should be stressed that ACAR is particularly sensitive for study of Cu-O chain Fermi sheets in YBCO.The character of ACAR, as experimentally detected, was predicted theoretically [87], and beside intersection of FL in Γ-X, Γ-Y directions, prediction has also been for FL intersection in Γ-T, Γ-U directions by Cu-O chain band.
With respect to topology of Cu-O chain band, as discussed above, this band should be the best candidate to be considered for gap opening in both, off-nodal (Γ-X, Γ-Y) and nodal (Γ-S) directions as seen in the ARPES results for YBCO [30].One should stress, however, that also in case of gap opening in Cu-O chain band, expected (d x 2 −y 2 ) symmetry is lost.Moreover, the authors in [30] attribute the effects seen in the ARPES to bands of Cu-O 2 planes and declare that in the spectra dispersion of Cu-O chain band is not present.There are a lot of experimental peculiarities of ARPES experiments with YBCO (twinning/untwining, surface states, CuO-chain problems, spectral dependence on ligh-source energy, etc.), so it is difficult for me to make any comment about assignment of particular EDC/MDC to CuO 2 plane or Cu-O chain bands as it has been done by the authors [30].
Nonetheless, an interpretation of the effects seen in the discussed ARPES from the stand-point of antiadiabatic theory should be of interest.The basic aspects concerning YBCO have been predicted [76] and shortly mentioned also in the Section 3.2.1.One should start with topology of the band structure (Figures 1(e)-1(f)), in particular with dispersion of Cu-O chain band at FL. ( It should be stressed that HF-SCF method with semiempirical INDO Hamiltonian used at band structure calculation (see reference 12 in [76]), overestimates bonding character and consequently bandwidth, which means that high-energy effects can hardly be studied, but for low-energy physics (like gap opening, kink formation, etc.) the method is reliable enough at least in a qualitative way.)It can be seen that dispersion of this band at FL corresponds qualitatively to Cu-O chain band dispersion as it has been calculated by FLAPW method [85] and to its experimental character as seen at FL by ACAR method [86].
Coupling to A g , B 2g , B 3g phonon modes induces fluctuation (Figures 1(e)-1(f)) of the ACP (inflex) of Cu-O 2 plane (d x 2 −y 2 − pσ) band at Y point across the FL.In the moment when the ACP approaches FL from the bonding side, strong renormalization of dispersion of this band (kink formation) could be seen by ARPES in offnodal Y-Γ direction (Figure 3(a), this paper) at about k ≈ 0.41/3.82≈ 0.1 A −1 (cf.[28, Figure 2a(c)]).At this situation, there is dramatic decrease of effective electron velocity and chemical potential (μ ad μ antiad < Dω), while DOS is increased (Figure 3(b), this paper) at FL and system undergoes transition into antiadiabatic state.Influence of nuclear dynamics on electronic structure is now significant, antiadiabatic theory of e-p coupling as presented in this paper.System is stabilized in antiadiabatic state at distorted nuclear configuration R d,cr and gap is opened in one-particle spectrum.The gap (Section 3.2.1,Figures 4(a ) and also in nodal Γ-S direction.The half-gaps in the off-nodal directions are Δ ΓY /2, ≈ 22 meV, Δ ΓX /2, ≈ 15 meV and gap in the nodal direction is expected to be Δ ΓS /2 ≤ 15 meV (cf.[28,Figures 1,7,8]).In contrast to [30], antiadiabatic theory yields a/b-gap asymmetry, which is in good agreement with other experimental results [79,80].At temperatures > T c , the gap extinct.It should be manifested by disappearing of peaks on EDC (MDC) at FL in the ARPES spectra (see [30,Figure 8(e)]).In this situation, the adiabatic DOS, which is of constant value at FL, is established and system is in adiabatic-nonsuperconducting state.
It should be mentioned that considerably smaller gap (Δ ΓT/U ≈ 5 meV) has been predicted [76] to be opened also in Γ-T/U directions.The prediction is related to the topology of the Cu-O chain (d − pσ) band in these directions ( [76,Figures 10,11]).Theoretical calculation of the electronpositron momentum density [87] confirms this character of Cu-O chain topology.It is obvious that symmetry of the Cu-O chain band gap is not of (d x 2 −y 2 ) character.Symmetry of the gap is, in my opinion, the matter of band structure topology at FL, which is far more complicated than the one emerging from simple model of CuO 2 plane confined in a square lattice.It is a complex mater of crystal structure and chemical composition of particular cuprate.It should be reminded that YBCO is the only hight-T c cuprate with Cu-O chain in its structure.The gap opening in Cu-O chain band does not mean, however, that superconductivity in YBCO is realized in Cu-O chains.Superconductivity is realized by bipolaron mechanism in Cu-O 2 planes, no matter if gap symmetry is (d x 2 −y 2 ) or any other (antiadiabatic theory, see also [76]).
In conclusion, it can be summarized that based on the ab initio theory of complex electronic ground state of superconductors, it can be concluded that e-p coupling in superconductors induces the temperature-dependent electronic structure instability related to fluctuation of analytic critical point (ACP--maximum, minimum, or saddle point of dispersion) of some band across FL, which results in breakdown of the adiabatic BOA.When ACP approaches FL, chemical potential μ ad is substantially reduced to μ antiad (μ ad μ antiad < Dω).Under these circumstances, the system is stabilized due to the effect of nuclear dynamics, in the antiadiabatic state at broken symmetry with a gap in one-particle spectrum.Distorted nuclear structure, which is related to couple of nuclei in the phonon mode r that induces transition into antiadiabatic state, has fluxional character.It has been shown that while system remains in antiadiabatic state, nonadiabatic polaron-renormalized phonon interactions are zero in well-defined k-region of reciprocal lattice.Along with geometric degeneracy of the antiadiabatic ground state it enables formation of mobile bipolarons (in a form of polarized intersite charge density distribution in real space) that can moveover lattice in external electric potential as supercarriers without dissipation.With increasing T, the stabilization effect of nuclear kinetic energy to the electronic ground-state energy decreases and at critical temperature T c the gap(s) extinct and system is stabilized in the adiabatic metal-like state with a continuum of states at FL, which is characteristic by high-symmetry structure.
An analysis of e-p interaction Hamiltonian has shown that an effective attractive e-e interaction, which is, the basis of Cooper's pair formation, is in fact the correction to electron correlation energy at transition from adiabatic into antiadiabatic ground electronic state.In this respect, increased electron correlation is not the primary reason for transition into superconducting state, but it is a consequence of antiadiabatic state formation which is stabilized by nonadiabatic e-p interactions at broken translation symmetry.It has been shown that thermodynamic properties of system in the antiadiabatic state correspond to thermodynamics of superconducting state.
crude-adiabatic (a, a + , b, b + ), new adiabatic (a, a + , b, b + ), and nonadiabatic (a, a + , b, b + ) quasiparticles are introduced.At these circumstances, canonical transformation of Hamiltonian means introduction of new dynamical variables

2
Cu 3 O 7 (oP14, Pmmm, #47, with chain oxygen in b-direction and vacancy in a-direction), the BS fluctuation is associated with coupling to three modes, A g , B 2g , B 3g , with the apical O(4) and CuO-plane O(2), O(3) atom displacements.At displacements ≈ 0.031 A • of apical O(4) in the A g mode and ≈ 0.022 A • of O(2) and O(3) in the B 2g , B 3g modes, the ACP (saddle point) of one of the Cu-O plane (d-pσ) band in Y point crosses FL and undergoes periodic fluctuation between topologies 1e ↔1f.

Figure 1 :
Figure 1: Band structures of MgB 2 (a, b), YB 6 (c, d), YBa 2 Cu 3 O 7 (e, f) and Zn (g, h).Pictures on the left (a, c, e, g) correspond to equilibrium high-symmetry structures.On the right are the band structures (b, d, f, h) of distorted geometry with atom displacements in the respective phonon modes.

Figure 3 :
Figure 3: Calculated dispersion of CuO 2 layer (d x 2 −y 2 − pσ) band in the ΓY direction with the kink formation indicated by the arrow (a) and increase in the (d x 2 −y 2 − pσ) band DOS when ACP approaches the FL (b) for YBa 2 Cu 3 O 7 at e-p coupling (see text).In both figures, the energy of Fermi level is rescaled to 0 eV.

Figure 5 :
Figure 5: Iso-density line of highest electron density at equilibrium nuclear geometry R eq of MgB 2 (a) and at the distorted geometry of B atoms in E 2g phonon mode R d,cr when nonadiabatic e-p interactions trigger transition into antiadiabatic state and intersites polarization is induced (b).

) 30
Advances in Condensed Matter PhysicsIn this expression, summation over bands is not explicitly indicated, but it should be understand implicitly.This correction is due to pairs of electrons with opposite quasimomentum and antiparallel spins (k ↑, −k ↓).It should be noticed that it is the contribution of biexcited configurations{Φ (k ,−k ) → (k,−k) } (i.e., two-particle (k ↑, −k ↓), two-hole (k ↑, −k ↓)excited singlet states) to the electronic ground state, that is, represented by renormalized Fermi vacuum Φ 0 .Expressed explicitly, first nonzero contributions are from the matrix elements of the typeΦ 0 |ΔH ep |Φ (k ,−k ) → (k,−k)2, that is, contributions in the second order of perturbation theory.Now, {ε k } represent particle states that are occupied above Fermi level and {ε k } are, due to excitations, empty-hole states below Fermi level.From (162), it follows that due to e-p interactions on the Q, P-dependent nonadiabatic level, the electron correlation energy (54) of the electronic ground state is increased by contribution of this correction.Inspection of (162) indicates that the largest contribution to the correlation energy correction is for pair of electrons at Fermi level which satisfy the condition |ε 0 k ΔH ep Δε/Dω → 0 = lim H eff (Fr) Δε/Dω → 0 the form of the Fröhlich two-particle effective Hamiltonian