Story
Fast Real-Axis Eliashberg Calculations: Full-bandwidth solutions beyond the constant density of states approximation
Key takeaway
Researchers developed a new algorithm to model superconductivity more accurately, which could lead to improved understanding and prediction of real-world superconducting materials.
Quick Explainer
The core idea of this work is to efficiently solve the finite-temperature Migdal-Eliashberg equations directly on the real-frequency axis, without relying on the common constant density of states approximation. This approach enables the computation of real-frequency observables like spectral functions and transport properties, while accurately capturing the energy-dependent electronic structure - a crucial factor for materials like H₃S that exhibit strong particle-hole asymmetry near the Fermi level. The key innovations are numerical techniques that reduce the computational complexity of the real-axis integrals, and a semi-analytic quadrature method to accurately evaluate integrals of the electronic spectral function. This direct real-axis solution method provides significant efficiency gains over the conventional approach of analytic continuation from the imaginary axis.
Deep Dive
Technical Deep Dive: Fast Real-Axis Eliashberg Calculations
Overview
This work presents an efficient numerical method for solving the finite-temperature Migdal-Eliashberg equations directly on the real-frequency axis. This approach allows for the computation of real-frequency observables like spectral functions, optical response, and transport properties without the need for ill-conditioned analytic continuation from the imaginary axis.
The key innovations are:
- Formulation of the real-axis Eliashberg equations that accounts for an energy-dependent electronic density of states, avoiding the common constant density of states approximation.
- Numerical techniques that reduce the computational complexity of the real-axis integrals from quadratic to linear in the number of sampling points.
- Demonstration of the approach on the superconductor H₃S, which exhibits a prominent van Hove singularity near the Fermi level.
Methodology
- Derived the finite-temperature Migdal-Eliashberg equations in a form suitable for direct real-frequency calculations, retaining the full energy dependence of the electronic density of states and screened Coulomb interaction.
- Introduced an efficient numerical algorithm to evaluate the integral kernels, reducing the computational complexity from quadratic to linear in the number of sampling points.
- Employed a semi-analytic quadrature method to accurately evaluate integrals of the electronic spectral function, which has sharp peaks near the poles of the Green's function.
- Implemented fixed-point iterations to solve the coupled real-axis Eliashberg equations, enforcing the expected symmetry relations to maintain causality of the Green's function.
Results
- Applied the method to the superconductor H₃S, which features a prominent van Hove singularity near the Fermi level.
- The full-bandwidth solution for the superconducting gap at low temperatures is ~60 meV, in good agreement with experimental tunneling measurements, compared to ~75 meV predicted by the constant density of states approximation.
- Computed the temperature-dependent spectral function, quasiparticle density of states, and occupancies, clearly capturing the strong particle-hole asymmetry induced by the van Hove singularity.
- Compared the direct real-axis solutions to those obtained via analytic continuation from the imaginary axis, showing that the ill-conditioned nature of the continuation obscures fine spectral details.
Interpretation
- The efficiency of the real-axis approach, with typical runtimes of milliseconds to minutes, enables the computation of transport coefficients, optical conductivities, and time-dependent response functions.
- This opens the door to more systematic investigations of transport properties and nonequilibrium dynamics in conventional superconductors within the Migdal-Eliashberg framework.
- The treatment of the full electronic structure, rather than the constant density of states approximation, is crucial for accurately capturing the superconducting properties of materials like H₃S that exhibit strong particle-hole asymmetry near the Fermi level.
Limitations & Uncertainties
- The method relies on the validity of the Migdal-Eliashberg theory, which assumes a weak electron-phonon coupling and neglects vertex corrections.
- While the real-axis solutions agree qualitatively with analytic continuation, the latter suffers from numerical instabilities that can obscure fine spectral details.
- The computations were performed at zero pressure; the effects of pressure on the electronic structure and electron-phonon coupling in H₃S were not explored.
What Comes Next
- Extend the real-axis Eliashberg approach to study the nonequilibrium dynamics and time-dependent response of superconducting materials, leveraging the efficiency of the direct real-frequency solutions.
- Apply the method to other superconducting systems with notable particle-hole asymmetry near the Fermi level to further assess the importance of treating the full electronic structure.
- Investigate ways to improve the numerical stability and efficiency of the real-axis Eliashberg solver, potentially exploring more advanced quadrature techniques or alternative formulations.