-
Polarons in two-dimensional atomic crystals
Authors:
Weng Hong Sio,
Feliciano Giustino
Abstract:
The polaron is the archetypal example of a quasiparticle emerging from the interaction between fermionic and bosonic fields in quantum field theory. In crystalline solids, polarons are formed when electrons and holes become dressed by the quanta of lattice vibrations. While experimental signatures of polarons in bulk three-dimensional materials abound, only rarely have polarons been observed in tw…
▽ More
The polaron is the archetypal example of a quasiparticle emerging from the interaction between fermionic and bosonic fields in quantum field theory. In crystalline solids, polarons are formed when electrons and holes become dressed by the quanta of lattice vibrations. While experimental signatures of polarons in bulk three-dimensional materials abound, only rarely have polarons been observed in two-dimensional atomic crystals. Here, we shed light on this asymmetry by developing a quantitative ab initio theory of polarons in atomically-thin crystals. Using this conceptual framework, we unravel the real-space structure of the recently-observed hole polaron in hexagonal boron nitride, we discover an unexpected critical condition for the existence of polarons in two-dimensional crystals, and we establish the key materials descriptors and the universal laws that underpin polaron physics in two dimensions.
△ Less
Submitted 16 February, 2023;
originally announced February 2023.
-
Electron-phonon physics from first principles using the EPW code
Authors:
Hyungjun Lee,
Samuel Poncé,
Kyle Bushick,
Samad Hajinazar,
Jon Lafuente-Bartolome,
Joshua Leveillee,
Chao Lian,
Francesco Macheda,
Hari Paudyal,
Weng Hong Sio,
Marios Zacharias,
Xiao Zhang,
Nicola Bonini,
Emmanouil Kioupakis,
Elena R. Margine,
Feliciano Giustino
Abstract:
EPW is an open-source software for $\textit{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 on ultra-fine Brillouin zone grids. This data is employed for predictive calculations of temperature-d…
▽ More
EPW is an open-source software for $\textit{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 on ultra-fine Brillouin zone grids. This data is employed for predictive calculations of temperature-dependent properties and phonon-assisted quantum processes in bulk solids and low-dimensional materials. Here, we report on significant new developments in the code that occurred during the period 2016-2022, namely: a transport module for the calculation of charge carrier mobility and conductivity under electric and magnetic fields within the $\textit{ab initio}$ Boltzmann transport equation; a superconductivity module for the calculation of critical temperature and gap structure in phonon-mediated superconductors within the $\textit{ab initio}$ 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 using the $\textit{ab initio}$ polaron equations; and a module for calculating electron-phonon couplings, 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. We describe recent code refactoring to prepare EPW for exascale architectures, we discuss efficient parallelization strategies, and report on extreme parallel scaling tests.
△ Less
Submitted 15 February, 2023;
originally announced February 2023.
-
Ab initio self-consistent many-body theory of polarons at all couplings
Authors:
Jon Lafuente-Bartolome,
Chao Lian,
Weng Hong Sio,
Idoia G. Gurtubay,
Asier Eiguren,
Feliciano Giustino
Abstract:
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…
▽ More
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 new 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 in [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.
△ Less
Submitted 12 August, 2022;
originally announced August 2022.
-
Unified approach to polarons and phonon-induced band structure renormalization
Authors:
Jon Lafuente-Bartolome,
Chao Lian,
Weng Hong Sio,
Idoia G. Gurtubay,
Asier Eiguren,
Feliciano Giustino
Abstract:
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…
▽ More
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.
△ Less
Submitted 12 August, 2022;
originally announced August 2022.
-
Unified ab initio description of Fröhlich electron-phonon interactions in two-dimensional and three-dimensional materials
Authors:
Weng Hong Sio,
Feliciano Giustino
Abstract:
\textit{Ab initio} calculations of electron-phonon interactions including the polar Fröhlich coupling have advanced considerably in recent years. The Fröhlich electron-phonon matrix element is by now well understood in the case of bulk three-dimensional (3D) materials. In the case of two-dimensional (2D) materials, the standard procedure to include Fröhlich coupling is to employ Coulomb truncation…
▽ More
\textit{Ab initio} calculations of electron-phonon interactions including the polar Fröhlich coupling have advanced considerably in recent years. The Fröhlich electron-phonon matrix element is by now well understood in the case of bulk three-dimensional (3D) materials. In the case of two-dimensional (2D) materials, the standard procedure to include Fröhlich coupling is to employ Coulomb truncation, so as to eliminate artificial interactions between periodic images of the 2D layer. While these techniques are well established, the transition of the Fröhlich coupling from three to two dimensions has not been investigated. Furthermore, it remains unclear what error one makes when describing 2D systems using the standard bulk formalism in a periodic supercell geometry. Here, we generalize previous work on the \textit{ab initio} Fröhlich electron-phonon matrix element in bulk materials by investigating the electrostatic potential of atomic dipoles in a periodic supercell consisting of a 2D material and a continuum dielectric slab. We obtain a unified expression for the matrix element, which reduces to the existing formulas for 3D and 2D systems when the interlayer separation tends to zero or infinity, respectively. This new expression enables an accurate description of the Fröhlich matrix element in 2D systems without resorting to Coulomb truncation. We validate our approach by direct \textit{ab initio} density-functional perturbation theory calculations for monolayer BN and MoS$_2$, and we provide a simple expression for the 2D Fröhlich matrix element that can be used in model Hamiltonian approaches. The formalism outlined in this work may find applications in calculations of polarons, quasiparticle renormalization, transport coefficients, and superconductivity, in 2D and quasi-2D materials.
△ Less
Submitted 15 March, 2022; v1 submitted 1 March, 2022;
originally announced March 2022.
-
Ab initio theory of polarons: formalism and applications
Authors:
Weng Hong Sio,
Carla Verdi,
Samuel Ponce,
Feliciano Giustino
Abstract:
We develop a theoretical and computational framework to study polarons in semiconductors and insulators from first principles. Our approach provides the formation energy, excitation energy, and wavefunction of both electron and hole polarons, and takes into account the coupling of the electron or hole to all phonons. An important feature of the present method is that it does not require supercell…
▽ More
We develop a theoretical and computational framework to study polarons in semiconductors and insulators from first principles. Our approach provides the formation energy, excitation energy, and wavefunction of both electron and hole polarons, and takes into account the coupling of the electron or hole to all phonons. An important feature of the present method is that it does not require supercell calculations, and relies exclusively on electron band structures, phonon dispersions, and electron-phonon matrix elements obtained from calculations in the crystal unit cell. Starting from the Kohn-Sham (KS) equations of density-functional theory, we formulate the polaron problem as a variational minimization, and we obtain a nonlinear eigenvalue problem in the basis of KS states and phonon eigenmodes. In our formalism the electronic component of the polaron is expressed as a coherent superposition of KS states, in close analogy with the solution of the Bethe-Salpeter equation for the calculation of excitons. We demonstrate the power of the methodology by studying polarons in LiF and Li2O2. We show that our method describes both small and large polarons, and seamlessly captures Frohlich-type polar electron-phonon coupling and non-Frohlich coupling to acoustic and optical phonons. To analyze in quantitative terms the electron-phonon coupling mechanisms leading to the formation of polarons, we introduce spectral decompositions similar to the Eliashberg spectral function. We validate our theory using both analytical results and direct calculations on large supercells. This study constitutes a first step toward complete ab initio many-body calculations of polarons in real materials.
△ Less
Submitted 19 June, 2019;
originally announced June 2019.
-
Polarons from first principles, without supercells
Authors:
Weng Hong Sio,
Carla Verdi,
Samuel Ponce,
Feliciano Giustino
Abstract:
We develop a formalism and a computational method to study polarons in insulators and semi-conductors from first principles. Unlike in standard calculations requiring large supercells, we solve a secular equation involving phonons and electron-phonon matrix elements from density-functional perturbation theory, in a spirit similar to the Bethe-Salpeter equation for excitons. We show that our approa…
▽ More
We develop a formalism and a computational method to study polarons in insulators and semi-conductors from first principles. Unlike in standard calculations requiring large supercells, we solve a secular equation involving phonons and electron-phonon matrix elements from density-functional perturbation theory, in a spirit similar to the Bethe-Salpeter equation for excitons. We show that our approach describes seamlessly large and small polarons, and we illustrate its capability by calculating wavefunctions, formation energies, and spectral decomposition of polarons in LiF and Li2O2.
△ Less
Submitted 19 June, 2019;
originally announced June 2019.
-
Cs$_2$InAgCl$_6$: A new lead-free halide double perovskite with direct band gap
Authors:
George Volonakis,
Amir A. Haghighirad,
Rebecca L. Milot,
Weng H. Sio,
Marina R. Filip,
Bernard Wenger,
Michael B. Johnston,
Laura M. Herz,
Henry J. Snaith,
Feliciano Giustino
Abstract:
A$_2$BB$^\prime$X$_6$ halide double perovskites based on bismuth and silver have recently been proposed as potential environmentally-friendly alternatives to lead-based hybrid halide perovskites. In particular, Cs$_2$BiAgX$_6$ (X = Cl, Br) have been synthesized and found to exhibit band gaps in the visible range. However, the band gaps of these compounds are indirect, which is not ideal for applic…
▽ More
A$_2$BB$^\prime$X$_6$ halide double perovskites based on bismuth and silver have recently been proposed as potential environmentally-friendly alternatives to lead-based hybrid halide perovskites. In particular, Cs$_2$BiAgX$_6$ (X = Cl, Br) have been synthesized and found to exhibit band gaps in the visible range. However, the band gaps of these compounds are indirect, which is not ideal for applications in thin film photovoltaics. Here, we propose a new class of halide double perovskites, where the B$^{3+}$ and B$^{+}$ cations are In$^{3+}$ and Ag$^{+}$, respectively. Our first-principles calculations indicate that the hypothetical compounds Cs$_2$InAgX$_6$ (X = Cl, Br, I) should exhibit direct band gaps between the visible (I) and the ultraviolet (Cl). Based on these predictions, we attempt to synthesize Cs$_2$InAgCl$_6$ and Cs$_2$InAgBr$_6$, and we succeed to form the hitherto unknown double perovskite Cs$_2$InAgCl$_6$. X-ray diffraction yields a double perovskite structure with space group $Fm\overline{3}m$. The measured band gap is 3.3 eV, and the compound is found to be photosensitive and turns reversibly from white to orange under ultraviolet illumination. We also perform an empirical analysis of the stability of Cs$_2$InAgX$_6$ and their mixed halides based on Goldschmidt's rules, and we find that it should also be possible to form Cs$_2$InAg(Cl$_{1-x}$Br$_{x}$)$_6$ for $x<1$. The synthesis of mixed halides will open the way to the development of lead-free double perovskites with direct and tunable band gaps.
△ Less
Submitted 25 January, 2017; v1 submitted 16 November, 2016;
originally announced November 2016.