No Arabic abstract
First-principles calculations of $e$-ph interactions are becoming a pillar of electronic structure theory. However, the current approach is incomplete. The piezoelectric (PE) $e$-ph interaction, a long-range scattering mechanism due to acoustic phonons in non-centrosymmetric polar materials, is not accurately described at present. Current calculations include short-range $e$-ph interactions (obtained by interpolation) and the dipole-like Frohlich long-range coupling in polar materials, but lack important quadrupole effects for acoustic modes and PE materials. Here we derive and compute the long-range $e$-ph interaction due to dynamical quadrupoles, and apply this framework to investigate $e$-ph interactions and the carrier mobility in the PE material wurtzite GaN. We show that the quadrupole contribution is essential to obtain accurate $e$-ph matrix elements for acoustic modes and to compute PE scattering. Our work resolves the outstanding problem of correctly computing $e$-ph interactions for acoustic modes from first principles, and enables studies of $e$-ph coupling and charge transport in PE materials.
Electron-phonon ($e$-ph) interactions are pervasive in condensed matter, governing phenomena such as transport, superconductivity, charge-density waves, polarons and metal-insulator transitions. First-principles approaches enable accurate calculations of $e$-ph interactions in a wide range of solids. However, they remain an open challenge in correlated electron systems (CES), where density functional theory often fails to describe the ground state. Therefore reliable $e$-ph calculations remain out of reach for many transition metal oxides, high-temperature superconductors, Mott insulators, planetary materials and multiferroics. Here we show first-principles calculations of $e$-ph interactions in CES, using the framework of Hubbard-corrected density functional theory (DFT+$U$ ) and its linear response extension (DFPT+$U$), which can describe the electronic structure and lattice dynamics of many CES. We showcase the accuracy of this approach for a prototypical Mott system, CoO, carrying out a detailed investigation of its $e$-ph interactions and electron spectral functions. While standard DFPT gives unphysically divergent and short-ranged $e$-ph interactions, DFPT+$U$ is shown to remove the divergences and properly account for the long-range Frohlich interaction, allowing us to model polaron effects in a Mott insulator. Our work establishes a broadly applicable and affordable approach for quantitative studies of e-ph interactions in CES, a novel theoretical tool to interpret experiments in this broad class of materials.
The interaction between electrons and lattice vibrations determines key physical properties of materials, including their electrical and heat transport, excited electron dynamics, phase transitions, and superconductivity. We present a new ab initio method that employs atomic orbital (AO) wavefunctions to compute the electron-phonon (e-ph) interactions in materials and interpolate the e-ph coupling matrix elements to fine Brillouin zone grids. We detail the numerical implementation of such AO-based e-ph calculations, and benchmark them against direct density functional theory calculations and Wannier function (WF) interpolation. The key advantages of AOs over WFs for e-ph calculations are outlined. Since AOs are fixed basis functions associated with the atoms, they circumvent the need to generate a material-specific localized basis set with a trial-and-error approach, as is needed in WFs. Therefore, AOs are ideal to compute e-ph interactions in chemically and structurally complex materials for which WFs are challenging to generate, and are also promising for high-throughput materials discovery. While our results focus on AOs, the formalism we present generalizes e-ph calculations to arbitrary localized basis sets, with WFs recovered as a special case.
We present an ab initio theory of the spin-wave stiffness tensor for ordered and disordered itinerant ferromagnets with pair exchange interactions derived from a method of infinitesimal spin rotations. The resulting formula bears an explicit form of a linear-response coefficient which involves one-particle Greens functions and effective velocity operators encountered in a recent theory of electron transport. Application of this approach to ideal metal crystals yields more reliable values of the spin stiffness than traditional ill-convergent real-space lattice summations. The formalism can also be combined with the coherent potential approximation for an effective-medium treatment of random alloys, which leads naturally to an inclusion of disorder-induced vertex corrections to the spin stiffness. The calculated concentration dependence of the spin-wave stiffness of random fcc Ni-Fe alloys can be ascribed to a variation of the reciprocal value of alloy magnetization. Calculations for random iron-rich bcc Fe-Al alloys reveal that their spin-wave stiffness is strongly reduced owing to the atomic ordering; this effect takes place due to weakly coupled local magnetic moments of Fe atoms surrounded by a reduced number of Fe nearest neighbors.
We present a comprehensive ab initio study of structural, electronic, lattice dynamical and electron-phonon coupling properties of the Bi(111) surface within density functional perturbation theory. Relativistic corrections due to spin-orbit coupling are consistently taken into account. As calculations are carried out in a periodic slab geometry, special attention is given to the convergence with respect to the slab thickness. Although the electronic structure of Bi(111) thin films varies significantly with thickness, we found that the lattice dynamics of Bi(111) is quite robust and appears converged already for slabs as thin as 6 bilayers. Changes of interatomic couplings are confined mostly to the first two bilayers, resulting in super-bulk modes with frequencies higher than the optic bulk spectrum, and in an enhanced density of states at lower frequencies for atoms in the first bilayer. Electronic states of the surface band related to the outer part of the hole Fermi surfaces exhibit a moderate electron-phonon coupling of about 0.45, which is larger than the coupling constant of bulk Bi. States at the inner part of the hole surface as well as those forming the electron pocket close to the zone center show much increased couplings due to transitions into bulk projected states near Gamma_bar. For these cases, the state dependent Eliashberg functions exhibit pronounced peaks at low energy and strongly deviate in shape from a Debye-like spectrum, indicating that an extraction of the coupling strength from measured electronic self-energies based on this simple model is likely to fail.
We report first principles calculations of the phonon dispersions of PbTe both for its observed structure and under compression. At the experimental lattice parameter we find a near instability of the optic branch at the zone center, in accord with experimental observations.This hardens quickly towards the zone boundary. There is also a very strong volume dependence of this mode, which is rapidly driven away from an instability by compression. These results are discussed inrelation to the thermal conductivity of the material.