No Arabic abstract
We investigate the applicability of Migdal--Eliashberg (ME) theory by revisiting Migdals analysis within the dynamical mean-field theory framework. First, we compute spectral functions, the quasi-particle weight, the self energy, renormalised phonon frequency and resistivity curves of the half-filled Holstein model. We demonstrate how ME theory has a phase-transition-like instability at intermediate coupling, and how the Engelsberg--Schrieffer (ES) picture is complicated by low-energy excitations from higher order diagrams (demonstrating that ES theory is a very weak coupling approach). Through consideration of the lowest-order vertex correction, we analyse the applicability of ME theory close to this transition. We find a breakdown of the theory in the intermediate coupling adiabatic limit due to a divergence in the vertex function. The region of applicability is mapped out, and it is found that ME theory is only reliable in the weak coupling adiabatic limit, raising questions about the accuracy of recent analyses of cuprate superconductors which do not include vertex corrections.
The superconducting (SC) and charge-density-wave (CDW) susceptibilities of the two dimensional Holstein model are computed using determinant quantum Monte Carlo (DQMC), and compared with results computed using the Migdal-Eliashberg (ME) approach. We access temperatures as low as 25 times less than the Fermi energy, $E_F$, which are still above the SC transition. We find that the SC susceptibility at low $T$ agrees quantitatively with the ME theory up to a dimensionless electron-phonon coupling $lambda_0 approx 0.4$ but deviates dramatically for larger $lambda_0$. We find that for large $lambda_0$ and small phonon frequency $omega_0 ll E_F$ CDW ordering is favored and the preferred CDW ordering vector is uncorrelated with any obvious feature of the Fermi surface.
We show that vertex corrections to Migdals theorem in general induce an odd-frequency spin-triplet superconducting order parameter, which coexists with its more commonly known even-frequency spin-singlet counterpart. Fully self-consistent vertex-corrected Eliashberg theory calculations for a two dimensional cuprate model, isotropically coupled to an Einstein phonon, confirm that both superconducting gaps are finite over a wide range of temperatures. The subordinate $d$-wave odd-frequency superconducting gap is found to be one order of magnitude smaller than the primary even-frequency $d$-wave gap. Our study provides a direct proof of concept for a previously unknown generation mechanism of odd-frequency superconductivity as well as for the generic coexistence of both superconducting states in bulk materials.
We formulate an efficient scheme to perform Migdal-Eliashberg calculation considering the retardation effect from first principles. While the conventional approach requires a huge number of Matsubara frequencies, we show that the intermediate representation of the Greens function [H. Shinaoka et al., Phys. Rev. B 96, 035147 (2017)] dramatically reduces the numerical cost to solve the linearized gap equation. Without introducing any empirical parameter, we demonstrate that we can successfully reproduce the experimental superconducting transition temperature of elemental Nb ($sim 10$ K) very accurately. The present result indicates that our approach has a superior performance for many superconductors for which $T_{rm c}$ is lower than ${mathcal O}(10)$ K
We simulate spectral functions for electron-phonon coupling in a filled band system - far from the asymptotic limit often assumed where the phonon energy is very small compared to the Fermi energy in a parabolic band and the Migdal theorem predicting 1+lambda quasiparticle renormalizations is valid. These spectral functions are examined over a wide range of parameter space through techniques often used in angle-resolved photoemission spectroscopy (ARPES). Analyzing over 1200 simulations we consider variations of the microscopic coupling strength, phonon energy and dimensionality for two models: a momentum-independent Holstein model, and momentum-dependent coupling to a breathing mode phonon. In this limit we find that any `effective coupling, lambda_eff, inferred from the quasiparticle renormalizations differs from the microscopic dimensionless coupling characterizing these Hamiltonians, lambda, and could drastically either over- or under-estimate it depending on the particular parameters and model. In contrast, we show that perturbation theory retains good predictive power for low coupling and small momenta, and that the momentum-dependence of the self-energy can be revealed via the relationship between velocity renormalization and quasiparticle strength. Additionally we find that (although not strictly valid) it is often possible to infer the self-energy and bare electronic structure through a self-consistent Kramers-Kronig bare-band fitting; and also that through lineshape alone, when Lorentzian, it is possible to reliably extract the shape of the imaginary part of a momentum-dependent self-energy without reference to the bare-band.
We perform single- and multi-band Migdal-Eliashberg (ME) calculations with parameters exctracted from density functional theory (DFT) simulations to study superconductivity in the electric-field-induced 2-dimensional hole gas at the hydrogenated (111) diamond surface. We show that according to the Eliashberg theory it is possible to induce a high-T$_{text{c}}$ superconducting phase when the system is field-effect doped to a surface hole concentration of $6times10^{14},$cm$^{-2}$, where the Fermi level crosses three valence bands. Starting from the band-resolved electron-phonon spectral functions $alpha^2F_{jj}(omega)$ computed ab initio, we iteratively solve the self-consistent isotropic Migdal-Eliashberg equations, in both the single-band and the multi-band formulations, in the approximation of a constant density of states at the Fermi level. In the single-band formulation, we find T$_{text{c}}approx40,$K, which is enhanced between $4%$ and $8%$ when the multi-band nature of the system is taken into account. We also compute the multi-band-sensistive quasiparticle density of states to act as a guideline for future experimental works.