Reversal of Band-Ordering Leads to High Hole Mobility in Strained p-type Scandium Nitride

Low hole mobility of nitride semiconductors is a significant impediment to realizing their high-efficiency device applications. Scandium nitride (ScN), an emerging rocksalt indirect band gap semiconductor, suffers from low hole mobility. Utilizing the ab initio Boltzmann transport formalism including spin-orbit coupling, here we show the dominating role of ionized impurity scattering in reducing the hole mobility in ScN thin films. We suggest a route to increase the hole mobility by reversing band ordering through strain engineering. Our calculation shows that the biaxial tensile strain in ScN lifts the split-off hole band above the heavy hole and light hole bands, leading to a lower hole-effective mass and increasing mobility. Along with the impurity scattering, the Fröhlich interaction also plays a vital role in the carrier scattering mechanism due to the polar nature of ScN. Increased hole mobility in ScN will lead to higher efficiencies in thermoelectric, plasmonics, and neuromorphic computing devices.

Electronic Structure of Few-Layer Black Phosphorus from μ-ARPES

Black phosphorus (BP) stands out among two-dimensional (2D) semiconductors because of its high mobility and thickness dependent direct band gap. However, the quasiparticle band structure of ultrathin BP has remained inaccessible to experiment thus far. Here we use a recently developed laser-based microfocus angle resolved photoemission (μ-ARPES) system to establish the electronic structure of 2-9 layer BP from experiment. Our measurements unveil ladders of anisotropic, quantized subbands at energies that deviate from the scaling observed in conventional semiconductor quantum wells. We quantify the anisotropy of the effective masses and determine universal tight-binding parameters, which provide an accurate description of the electronic structure for all thicknesses.

Accurate Prediction of Hall Mobilities in Two-Dimensional Materials through Gauge-Covariant Quadrupolar Contributions

Despite considerable efforts, accurate computations of electron-phonon and carrier transport properties of low-dimensional materials from first principles have remained elusive. By building on recent advances in the description of long-range electrostatics, we develop a general approach to the calculation of electron-phonon couplings in two-dimensional materials. We show that the nonanalytic behavior of the electron-phonon matrix elements depends on the Wannier gauge, but that a missing Berry connection restores invariance to quadrupolar order. We showcase these contributions in a MoS_{2} monolayer, calculating intrinsic drift and Hall mobilities with precise Wannier interpolations. We also find that the contributions of dynamical quadrupoles to the scattering potential are essential, and that their neglect leads to errors of 23% and 76% in the room-temperature electron and hole Hall mobilities, respectively.

Limits to Electrical Mobility in Lead-Halide Perovskite Semiconductors

Semiconducting polycrystalline thin films are cheap to produce and can be deposited on flexible substrates, yet high-performance electronic devices usually utilize single-crystal semiconductors, owing to their superior charge-carrier mobilities and longer diffusion lengths. Here we show that the electrical performance of polycrystalline films of metal-halide perovskites (MHPs) approaches that of single crystals at room temperature. Combining temperature-dependent terahertz conductivity measurements and ab initio calculations we uncover a complete picture of the origins of charge-carrier scattering in single crystals and polycrystalline films of CH₃NH₃PbI₃. We show that Fröhlich scattering of charge carriers with multiple phonon modes is the dominant mechanism limiting mobility, with grain-boundary scattering further reducing mobility in polycrystalline films. We reconcile the large discrepancy in charge-carrier diffusion lengths between single crystals and films by considering photon reabsorption. Thus, polycrystalline films of MHPs offer great promise for devices beyond solar cells, including light-emitting diodes and modulators.

Theory and Computation of Hall Scattering Factor in Graphene

The Hall scattering factor, r, is a key quantity for establishing carrier concentration and drift mobility from Hall measurements; in experiments, it is usually assumed to be 1. In this paper, we use a combination of analytical and ab initio modeling to determine r in graphene. Although at high carrier densities r ≈ 1 in a wide temperature range, at low doping the temperature dependence of r is very strong with values as high as 4 below 300 K. These high values are due to the linear bands around the Dirac cone and the carrier scattering rates due to acoustic phonons. At higher temperatures, r can instead become as low as 0.5 due to the contribution of both holes and electrons and the role of optical phonons. Finally, we provide a simple analytical model to compute accurately r in graphene in a wide range of temperatures and carrier densities.

Wannier90 as a community code: new features and applications

Wannier90 is an open-source computer program for calculating maximally-localised Wannier functions (MLWFs) from a set of Bloch states. It is interfaced to many widely used electronic-structure codes thanks to its independence from the basis sets representing these Bloch states. In the past few years the development of Wannier90 has transitioned to a community-driven model; this has resulted in a number of new developments that have been recently released in Wannier90 v3.0. In this article we describe these new functionalities, that include the implementation of new features for wannierisation and disentanglement (symmetry-adapted Wannier functions, selectively-localised Wannier functions, selected columns of the density matrix) and the ability to calculate new properties (shift currents and Berry-curvature dipole, and a new interface to many-body perturbation theory); performance improvements, including parallelisation of the core code; enhancements in functionality (support for spinor-valued Wannier functions, more accurate methods to interpolate quantities in the Brillouin zone); improved usability (improved plotting routines, integration with high-throughput automation frameworks), as well as the implementation of modern software engineering practices (unit testing, continuous integration, and automatic source-code documentation). These new features, capabilities, and code development model aim to further sustain and expand the community uptake and range of applicability, that nowadays spans complex and accurate dielectric, electronic, magnetic, optical, topological and transport properties of materials.

First-principles calculations of charge carrier mobility and conductivity in bulk semiconductors and two-dimensional materials

One of the fundamental properties of semiconductors is their ability to support highly tunable electric currents in the presence of electric fields or carrier concentration gradients. These properties are described by transport coefficients such as electron and hole mobilities. Over the last decades, our understanding of carrier mobilities has largely been shaped by experimental investigations and empirical models. Recently, advances in electronic structure methods for real materials have made it possible to study these properties with predictive accuracy and without resorting to empirical parameters. These new developments are unlocking exciting new opportunities, from exploring carrier transport in quantum matter to in silico designing new semiconductors with tailored transport properties. In this article, we review the most recent developments in the area of ab initio calculations of carrier mobilities of semiconductors. Our aim is threefold: to make this rapidly-growing research area accessible to a broad community of condensed-matter theorists and materials scientists; to identify key challenges that need to be addressed in order to increase the predictive power of these methods; and to identify new opportunities for increasing the impact of these computational methods on the science and technology of advanced materials. The review is organized in three parts. In the first part, we offer a brief historical overview of approaches to the calculation of carrier mobilities, and we establish the conceptual framework underlying modern ab initio approaches. We summarize the Boltzmann theory of carrier transport and we discuss its scope of applicability, merits, and limitations in the broader context of many-body Green's function approaches. We discuss recent implementations of the Boltzmann formalism within the context of density functional theory and many-body perturbation theory calculations, placing an emphasis on the key computational challenges and suggested solutions. In the second part of the article, we review applications of these methods to materials of current interest, from three-dimensional semiconductors to layered and two-dimensional materials. In particular, we discuss in detail recent investigations of classic materials such as silicon, diamond, gallium arsenide, gallium nitride, gallium oxide, and lead halide perovskites as well as low-dimensional semiconductors such as graphene, silicene, phosphorene, molybdenum disulfide, and indium selenide. We also review recent efforts toward high-throughput calculations of carrier transport. In the last part, we identify important classes of materials for which an ab initio study of carrier mobilities is warranted. We discuss the extension of the methodology to study topological quantum matter and materials for spintronics and we comment on the possibility of incorporating Berry-phase effects and many-body correlations beyond the standard Boltzmann formalism.

Route to High Hole Mobility in GaN via Reversal of Crystal-Field Splitting

A fundamental obstacle toward the realization of GaN p-channel transistors is its low hole mobility. Here we investigate the intrinsic phonon-limited mobility of electrons and holes in wurtzite GaN using the ab initio Boltzmann transport formalism, including all electron-phonon scattering processes and many-body quasiparticle band structures. We predict that the hole mobility can be increased by reversing the sign of the crystal-field splitting in such a way as to lift the split-off hole states above the light and heavy holes. We find that a 2% biaxial tensile strain can increase the hole mobility by 230%, up to a theoretical Hall mobility of 120  cm^{2}/V s at room temperature and 620  cm^{2}/V s at 100 K.

Polarons from First Principles, without Supercells

We develop a formalism and a computational method to study polarons in insulators and semiconductors 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 wave functions, formation energies, and spectral decomposition of polarons in LiF and Li_{2}O_{2}.

Dimensional Crossover in the Carrier Mobility of Two-Dimensional Semiconductors: The Case of InSe

Two-dimensional (2D) semiconductors are at the center of an intense research effort aimed at developing the next generation of flexible, transparent, and energy-efficient electronics. In these applications, the carrier mobility, that is the ability of electrons and holes to move rapidly in response to an external voltage, is a critical design parameter. Here, we show that the interlayer coupling between electronic wave functions in 2D semiconductors can be used to drastically alter carrier mobility and dynamics. We demonstrate this concept by performing state-of-the-art ab initio calculations for InSe, a prototypical 2D semiconductor that is attracting considerable attention, because of its exceptionally high electron mobility. We show that the electron mobility of InSe can be increased from 100 cm² V^-1 s^-1 to 1000 cm² V^-1 s^-1 by exploiting the dimensional crossover of the electronic density of states from two dimensions to three dimensions. By generalizing our results to the broader class of layered materials, we discover that dimensionality plays a universal role in the transport properties of 2D semiconductors.

Carrier Lifetimes and Polaronic Mass Enhancement in the Hybrid Halide Perovskite CH_{3}NH_{3}PbI_{3} from Multiphonon Fröhlich Coupling

We elucidate the nature of the electron-phonon interaction in the archetypal hybrid perovskite CH_{3}NH_{3}PbI_{3} using ab initio many-body calculations and an exactly solvable model. We demonstrate that electrons and holes near the band edges primarily interact with three distinct groups of longitudinal-optical vibrations, in order of importance: the stretching of the Pb-I bond, the bending of the Pb-I-Pb bonds, and the libration of the organic cations. These polar phonons induce ultrafast intraband carrier relaxation over timescales of 6-30 fs and yield polaron effective masses 28% heavier than the bare band masses. These findings allow us to rationalize previous experimental observations and provide a key to understanding carrier dynamics in halide perovskites.

Advanced capabilities for materials modelling with Quantum ESPRESSO

Quantum EXPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Quantum EXPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

Origin of Superconductivity and Latent Charge Density Wave in NbS_{2}

We elucidate the origin of the phonon-mediated superconductivity in 2H-NbS_{2} using the ab initio anisotropic Migdal-Eliashberg theory including Coulomb interactions. We demonstrate that superconductivity is associated with Fermi surface hot spots exhibiting an unusually strong electron-phonon interaction. The electron-lattice coupling is dominated by low-energy anharmonic phonons, which place the system on the verge of a charge density wave instability. We also provide definitive evidence for two-gap superconductivity in 2H-NbS_{2}, and show that the low- and high-energy peaks observed in tunneling spectra correspond to the Γ- and K-centered Fermi surface pockets, respectively. The present findings call for further efforts to determine whether our proposed mechanism underpins superconductivity in the whole family of metallic transition metal dichalcogenides.

Temperature Dependence of the Energy Levels of Methylammonium Lead Iodide Perovskite from First-Principles

Environmental effects and intrinsic energy-loss processes lead to fluctuations in the operational temperature of solar cells, which can profoundly influence their power conversion efficiency. Here we determine from first-principles the effects of temperature on the band gap and band edges of the hybrid pervoskite CH₃NH₃PbI₃ by accounting for electron-phonon coupling and thermal expansion. From 290 to 380 K, the computed band gap change of 40 meV coincides with the experimental change of 30-40 meV. The calculation of electron-phonon coupling in CH₃NH₃PbI₃ is particularly intricate as the commonly used Allen-Heine-Cardona theory overestimates the band gap change with temperature, and excellent agreement with experiment is only obtained when including high-order terms in the electron-phonon interaction. We also find that spin-orbit coupling enhances the electron-phonon coupling strength but that the inclusion of nonlocal correlations using hybrid functionals has little effect. We reach similar conclusions in the metal-halide perovskite CsPbI₃. Our results unambiguously confirm for the first time the importance of high-order terms in the electron-phonon coupling by direct comparison with experiment.

Temperature dependence of the electronic structure of semiconductors and insulators

The renormalization of electronic eigenenergies due to electron-phonon coupling (temperature dependence and zero-point motion effect) is sizable in many materials with light atoms. This effect, often neglected in ab initio calculations, can be computed using the perturbation-based Allen-Heine-Cardona theory in the adiabatic or non-adiabatic harmonic approximation. After a short description of the recent progresses in this field and a brief overview of the theory, we focus on the issue of phonon wavevector sampling convergence, until now poorly understood. Indeed, the renormalization is obtained numerically through a slowly converging q-point integration. For non-zero Born effective charges, we show that a divergence appears in the electron-phonon matrix elements at q → Γ, leading to a divergence of the adiabatic renormalization at band extrema. This problem is exacerbated by the slow convergence of Born effective charges with electronic wavevector sampling, which leaves residual Born effective charges in ab initio calculations on materials that are physically devoid of such charges. Here, we propose a solution that improves this convergence. However, for materials where Born effective charges are physically non-zero, the divergence of the renormalization indicates a breakdown of the adiabatic harmonic approximation, which we assess here by switching to the non-adiabatic harmonic approximation. Also, we study the convergence behavior of the renormalization and develop reliable extrapolation schemes to obtain the converged results. Finally, the adiabatic and non-adiabatic theories, with corrections for the slow Born effective charge convergence problem (and the associated divergence) are applied to the study of five semiconductors and insulators: α-AlN, β-AlN, BN, diamond, and silicon. For these five materials, we present the zero-point renormalization, temperature dependence, phonon-induced lifetime broadening, and the renormalized electronic band structure.