EPW is an open-source software for ab initio calculations of electron–phonon interactions and related materials properties. The code combines density functional perturbation theory and maximally localized Wannier functions to efficiently compute electron–phonon coupling matrix elements, and to perform predictive calculations of temperature-dependent properties and phonon-assisted quantum processes in bulk solids and low-dimensional materials. Here, we report on significant developments in the code since 2016, namely: a transport module for the calculation of charge carrier mobility under electric and magnetic fields using the Boltzmann transport equation; a superconductivity module for calculations of phonon-mediated superconductors using the anisotropic multi-band Eliashberg theory; an optics module for calculations of phonon-assisted indirect transitions; a module for the calculation of small and large polarons without supercells; and a module for calculating band structure renormalization and temperature-dependent optical spectra using the special displacement method. For each capability, we outline the methodology and implementation and provide example calculations.

Ab initio calculations of the phonon-induced band structure renormalization are currently based on the perturbative Allen-Heine theory and its many-body generalizations. These approaches are unsuitable to describe materials where electrons form localized polarons. Here, we develop a self-consistent, many-body Green’s function theory of band structure renormalization that incorporates localization and self-trapping. We show that the present approach reduces to the Allen-Heine theory in the weak-coupling limit, and to total energy calculations of self-trapped polarons in the strong-coupling limit. To demonstrate this methodology, we reproduce the path-integral results of Feynman and diagrammatic Monte Carlo calculations for the Fröhlich model at all couplings, and we calculate the zero point renormalization of the band gap of an ionic insulator including polaronic effects.

We present a theoretical framework to describe polarons from first principles within a many-body Green's function formalism. Starting from a general electron-phonon Hamiltonian, we derive a self-consistent Dyson equation in which the phonon-mediated self-energy is composed by two distinct terms. One term is the Fan-Migdal self-energy and describes dynamic electron-phonon processes, the other term is a contribution to the self-energy originating from the static displacements of the atomic nuclei in the polaronic ground state. The lowest-order approximation to the present theory yields the standard many-body perturbation theory approach to electron-phonon interactions in the limit of large polarons, and the ab initio polaron equations introduced [Sio et al., Phys. Rev. B 99, 235139 (2019); Phys. Rev. Lett. 122, 246403 (2019)] in the limit of small polarons. A practical recipe to implement the present unifying formalism in first-principles calculations is outlined. We apply our method to the Fröhlich model, and obtain remarkably accurate polaron energies at all couplings, in line with Feynman's polaron theory and diagrammatic Monte Carlo calculations. We also recover the well-known results of Fröhlich and Pekar at weak and strong coupling, respectively. The present approach enables predictive many-body calculations of polarons in real materials at all couplings.

Excitons consist of electrons and holes held together by their attractive Coulomb interaction. Although excitons are neutral excitations, spatial fluctuations in their charge density couple with the ions of the crystal lattice. This coupling can lower the exciton energy and lead to the formation of a localized excitonic polaron or even a self-trapped exciton in the presence of strong exciton-phonon interactions. Here, we develop a theoretical and computational approach to compute excitonic polarons and self-trapped excitons from first principles. Our methodology combines the many-body Bethe-Salpeter approach with density-functional perturbation theory and does not require explicit supercell calculations. As a proof of concept, we demonstrate our method for a compound of the halide perovskite family.

Excitons are neutral excitations that are composed of electrons and holes bound together by their attractive Coulomb interaction. The electron and the hole forming the exciton also interact with the underlying atomic lattice, and this interaction can lead to a trapping potential that favors exciton localization. The quasiparticle thus formed by the exciton and the surrounding lattice distortion is called excitonic polaron. Excitonic polarons have long been thought to exist in a variety of materials, and are often invoked to explain the Stokes shift between the optical absorption edge and the photoluminescence peak. However, quantitative ab initio calculations of these effects are exceedingly rare. In this manuscript, we present a theory of excitonic polarons that is amenable to first-principles calculations. We first apply this theory to model Hamiltonians for Wannier excitons experiencing Fröhlich or Holstein electron-phonon couplings. We find that, in the case of Fröhlich interactions, excitonic polarons only form when there is a significant difference between electron and hole effective masses. Then, we apply this theory to calculating excitonic polarons in lithium fluoride ab initio. The key advantage of the present approach is that it does not require supercells, therefore it can be used to study a variety of materials hosting either small or large excitonic polarons. This work constitutes the first step toward a complete ab initio many-body theory of excitonic polarons in real materials.

We present a comprehensive first-principles analysis of the nonadiabatic effects due to the electron-phonon interaction on the vibrational spectrum of the electron-doped monolayer ${\mathrm{MoS}}_2$. Deep changes in the Fermi surface upon doping cause the linewidth broadening of the normal modes governing the spin-conserving intervalley electronic scattering, which become unstable with the population of all the spin-split conduction valleys. We find that the nonadiabatic spectral effects modify dramatically the adiabatic dispersion of the long-wavelength optical phonon modes, responsible for intravalley scattering, as soon as inequivalent valleys get populated. These results are illustrated by means of a simple analytical model. Finally, we explain the emergence of an intricate dynamical structure for the strongly interacting out-of-plane polarized $A_1^{\prime }$ optical vibrational mode spectrum by means of a multiple-phonon quasiparticle picture defined in the full complex frequency plane, showing that this intriguing spectral structure originates from the splitting of the original adiabatic branch induced by the electron-phonon coupling.