No Arabic abstract
We extend the two leading methods for the emph{ab initio} computational descrip tion of phonon-mediated superconductors, namely Eliashberg theory and density fu nctional theory for superconductors (SCDFT), to include plasmonic effects. Furth ermore, we introduce a hybrid formalism in which the Eliashberg approximation fo r the electron-phonon coupling is combined with the SCDFT treatment of the dynam ically screened Coulomb interaction. The methods have been tested on a set of we ll-known conventional superconductors by studying how the plasmon contribution a ffects the phononic mechanism in determining the critical temperature (tc). Our simulations show that plasmonic SCDFT leads to a good agreement between predict ed and measured tcs, whereas Eliashberg theory considerably overestimates the plasmon-mediated pairing and, therefore, tc. The hybrid approach, on the other hand, gives results close to SCDFT and overall in excellent agreement with exper iments.
We present a first-principles approach to describe magnetic and superconducting systems and the phenomena of competition between these electronic effects. We develop a density functional theory: SpinSCDFT, by extending the Hohenberg-Kohn theorem and constructing the non-interacting Kohn- Sham system. An exchange-correlation functional for SpinSCDFT is derived from the Sham Schluter connection between the SpinSCDFT Kohn-Sham and a self-energy in Eliashberg approximation. The reference Eliashberg equations for superconductors in the presence of magnetism are also derived and discussed.
We address an issue of how to accurately include the self energy effect of the screened electron-electron Coulomb interaction in the phonon-mediated superconductors from first principles. In the Eliashberg theory for superconductors, self energy is usually decomposed using the $2times 2$ Pauli matrices in the electron-hole space. We examine how the diagonal ($sigma_{0}$ and $sigma_{3}$) components, which results in the quasiparticle correction to the normal state, behave in the homogeneous electron gas in order to establish a norm of treating those components in real metallic systems. Within the $G_{0}W_{0}$ approximation, we point out that these components are non-analytic near the Fermi surface but their directional derivatives and resulting corrections to the quasiparticle velocity are nevertheless well defined. In the low-energy spectrum, we observe large cancellation between effects of these components and, without the numerically more tedious $sigma_{3}$ component, the effective mass is incorrectly increased. Feasible paths to manage this cancellation in the ab initio Eliashberg calculations are discussed.
The newly discovered iron pnictide superconductors apparently present an unusual case of interband-channel pairing superconductivity. Here we show that, in the limit where the pairing occurs within the interband channel, several surprising effects occur quite naturally and generally: different density-of-states on the two bands lead to several unusual properties, including a gap ratio which behaves inversely to the ratio of density-of-states; the weak-coupling limit of the Eliashberg and the BCS theory, commonly taken as equivalent, in fact predict qualitatively different dependence of the $Delta_{1}/Delta_{2}$ and $Delta/T_{c}$ ratios on coupling constants. We show analytically that these effects follow directly from the interband character of superconductivity. Our results show that in the interband-only pairing model the maximal gap ratio is $sqrt{N_{2}/N_{1}}$ as strong-coupling effects act only to reduce this ratio. This suggests that if the large experimentally reported gap ratios (up to a factor 2) are correct, the pairing mechanism must include more intraband interaction than is usually assumed.
The standard Eliashberg - McMillan theory of superconductivity is essentially based on the adiabatic approximation. Here we present some simple estimates of electron - phonon interaction within Eliashberg - McMillan approach in non - adiabatic and even antiadiabatic situation, when characteristic phonon frequency $Omega_0$ becomes large enough, i.e. comparable or exceeding the Fermi energy $E_F$. We discuss the general definition of Eliashberg - McMillan (pairing) electron - phonon coupling constant $lambda$, taking into account the finite value of phonon frequencies. We show that the mass renormalization of electrons is in general determined by different coupling constant $tildelambda$, which takes into account the finite width of conduction band, and describes the smooth transition from the adiabatic regime to the region of strong nonadiabaticity. In antiadiabatic limit, when $Omega_0gg E_F$, the new small parameter of perturbation theory is $lambdafrac{E_F}{Omega_0}simlambdafrac{D}{Omega_0}ll 1$ ($D$ is conduction band half -- width), and corrections to electronic spectrum (mass renormalization) become irrelevant. However, the temperature of superconducting transition $T_c$ in antiadiabatic limit is still determined by Eliashberg - McMillan coupling constant $lambda$. We consider in detail the model with discrete set of (optical) phonon frequencies. A general expression for superconducting transition temperature $T_c$ is derived, which is valid in situation, when one (or several) of such phonons becomes antiadiabatic. We also analyze the contribution of such phonons into the Coulomb pseudopotential $mu^{star}$ and show, that antiadiabatic phonons do not contribute to Tolmachevs logarithm and its value is determined by partial contributions from adiabatic phonons only.
We numerically investigate the Spin Density Functional theory for superconductors (SpinSCDFT) and the approximated exchange-correlation functional, derived and presented in the preceding paper I. As a test system we employ a free electron gas featuring an exchange-splitting, a phononic pairing field and a Coulomb repulsion. SpinSCDFT results are compared with Sarma, the Bardeen Cooper and Schrieffer theory and with an Eliashberg type of approach. We find that the spectrum of the superconducting Kohn-Sham SpinSCDFT system is not in agreement with the true quasi particle structure. Therefore, starting from the Dyson equation, we derive a scheme that allows to compute the many body excitations of the superconductor and represents the extension to superconductivity of the G0W0 method in band structure theory. This superconducting G0 W0 method vastly improves the predicted spectra.