No Arabic abstract
The electronic transport behaviour of materials determines their suitability for technological applications. We develop an efficient method for calculating carrier scattering rates of solid-state semiconductors and insulators from first principles inputs. The present method extends existing polar and non-polar electron-phonon coupling, ionized impurity, and piezoelectric scattering mechanisms formulated for isotropic band structures to support highly anisotropic materials. We test the formalism by calculating the electronic transport properties of 16 semiconductors and comparing the results against experimental measurements. The present work is amenable for use in high-throughput computational workflows and enables accurate screening of carrier mobilities, lifetimes, and thermoelectric power.
The bulk photovoltaic effect (BPVE) has attracted an increasing interest due to its potential to overcome the efficiency limit of traditional photovoltaics, and much effort has been devoted to understanding its underlying physics. However, previous work has shown that theoretical models of the shift current and the phonon-assisted ballistic current in real materials do not fully account for the experimental BPVE photocurrent, and so other mechanisms should be investigated in order to obtain a complete picture of BPVE. In this Letter, we demonstrate two approaches that enable the ab initio calculation of the ballistic current originating from the electron-hole interaction in semiconductors. Using BaTiO$_3$ and MoS$_2$ as two examples, we show clearly that for them the asymmetric scattering from electron-hole interaction is less appreciable than that from electron-phonon interaction, indicating more scattering processes need to be included to further improve the BPVE theory. Moreover, our approaches build up a venue for predicting and designing materials with larger ballistic current due to electron-hole interactions.
Efficient ab initio computational methods for the calculation of thermoelectric transport properties of materials are of great avail for energy harvesting technologies. The BoltzTraP code has been largely used to efficiently calculate thermoelectric coefficients. However, its current version that is publicly available is based only on the constant relaxation time (RT) approximation, which usually does not hold for real materials. Here, we extended the implementation of the BoltzTraP code by incorporating realistic k-dependent RT models of the temperature dependence of the main scattering processes, namely, screened polar and nonpolar scattering by optical phonons, scattering by acoustic phonons, and scattering by ionized impurities with screening. Our RT models are based on a smooth Fourier interpolation of Kohn-Sham eigenvalues and its derivatives, taking into account non-parabolicity (beyond the parabolic or Kane models), degeneracy and multiplicity of the energy bands on the same footing, within very low computational cost. In order to test our methodology, we calculated the anisotropic thermoelectric transport properties of low temperature phase (Pnma) of intrinsic p-type and hole-doped tin selenide (SnSe). Our results are in quantitative agreement with experimental data, regarding the evolution of the anisotropic thermoelectric coefficients with both temperature and chemical potential. Hence, from this picture, we also obtained the evolution and understanding of the main scattering processes of the overall thermoelectric transport in p-type SnSe.
Identification and design of defects in two-dimensional (2D) materials as promising single photon emitters (SPE) requires a deep understanding of underlying carrier recombination mechanisms. Yet, the dominant mechanism of carrier recombination at defects in 2D materials has not been well understood, and some outstanding questions remain: How do recombination processes at defects differ between 2D and 3D systems? What factors determine defects in 2D materials as excellent SPE at room temperature? In order to address these questions, we developed first-principles methods to accurately calculate the radiative and non-radiative recombination rates at defects in 2D materials, using h-BN as a prototypical example. We reveal the carrier recombination mechanism at defects in 2D materials being mostly dominated by defect-defect state recombination in contrast to defect-bulk state recombination in most 3D semiconductors. In particular, we disentangle the non-radiative recombination mechanism into key physical quantities: zero-phonon line (ZPL) and Huang-Rhys factor. At the end, we identified strain can effectively tune the electron-phonon coupling at defect centers and drastically change non-radiative recombination rates. Our theoretical development serves as a general platform for understanding carrier recombination at defects in 2D materials, while providing pathways for engineering of quantum efficiency of SPE.
Light nuclei at room temperature and below exhibit a kinetic energy which significantly deviates from the predictions of classical statistical mechanics. This quantum kinetic energy is responsible for a wide variety of isotope effects of interest in fields ranging from chemistry to climatology. It also furnishes the second moment of the nuclear momentum distribution, which contains subtle information about the chemical environment and has recently become accessible to deep inelastic neutron scattering experiments. Here we show how, by combining imaginary time path integral dynamics with a carefully designed generalized Langevin equation, it is possible to dramatically reduce the expense of computing the quantum kinetic energy. We also introduce a transient anisotropic Gaussian approximation to the nuclear momentum distribution which can be calculated with negligible additional effort. As an example, we evaluate the structural properties, the quantum kinetic energy, and the nuclear momentum distribution for a first-principles simulation of liquid water.
Inelastic scattering experiments are key methods for mapping the full dispersion of fundamental excitations of solids in the ground as well as non-equilibrium states. A quantitative analysis of inelastic scattering in terms of phonon excitations requires identifying the role of multi-phonon processes. Here, we develop an efficient first-principles methodology for calculating the {it all-phonon} quantum mechanical structure factor of solids. We demonstrate our method by obtaining unprecedented agreement between measurements and calculations of the diffuse diffraction patterns of black phosphorus, showing that multi-phonon scattering plays a substantial role. The present approach constitutes a pivotal advancement in the interpretation of static and time-resolved electron, X-ray, and neutron inelastic scattering data.