No Arabic abstract
We present first-principles many-body perturbation theory calculations of the quasiparticle electronic structure and of the optical response of HfO$_2$ polymorphs. We use the $GW$ approximation including core electrons by the projector augmented wave (PAW) method and performing a quasiparticle self-consistency also on wavefunctions (QS$GW$). In addition, we solve the Bethe-Salpeter equation on top of $GW$ to calculate optical properties including excitonic effects. For monoclinic HfO$_2$ we find a fundamental band gap of $E_g = 6.33$ eV (with the direct band gap at $E_g^d = 6.41$ eV), and an exciton binding energy of 0.57 eV, which situates the optical gap at $E^o_g = 5.85$ eV. The latter is in the range of spectroscopic ellipsometry (SE) experimental estimates (5.5-6 eV), whereas our electronic band gap is well beyond experimental photoemission (PE) estimates ($< 6$ eV) and previous $GW$ works. Our calculated density of states and optical absorption spectra compare well to raw PE and SE spectra. This suggests that our predictions of both optical and electronic gaps are close to, or at least lower bounds of, the real values.
In the development of highly efficient photovoltaic cells, solid perovskite systems have demonstrated unprecedented promise, with the figure of merit exceeding nineteen percent of efficiency. In this paper, we investigate the optical and vibrational properties of organometallic cubic perovskite CH3NH3PbI3 using first-principles calculations. For accurate theoretical description, we go beyond conventional density functional theory (DFT), and calculated optical conductivity using relativist quasi-particle (GW) correction. Incorporating these many-body effects, we further solve Bethe-Salpeter equations (BSE) for excitons, and found enhanced optical conductivity near the gap edge. Due to the presence of organic methylammonium cations near the center of the perovskite cell, the system is sensitive to low energy vibrational modes. We estimate the phonon modes of CH3NH3PbI3 using small displacement approach, and further calculate the infrared absorption (IR) spectra. Qualitatively, our calculations of low-energy phonon frequencies are in good agreement with our terahertz measurements. Therefore, for both energy scales (around 2 eV and 0-20 meV), our calculations reveal the importance of many-body effects and their contributions to the desirable optical properties in the cubic organometallic perovskites system.
The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as Thiels set. Literature reference data based on Gaussian-type orbitals are reproduced to about one meV precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted tier2+aug2 is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.
We present a systematic investigation of the role and importance of excitonic effects on the optical properties of transitions metal oxide perovskites. A representative set of fourteen compounds has been selected, including 3$d$ (SrTiO$_3$, LaScO$_3$, LaTiO$_3$, LaVO$_3$, LaCrO$_3$, LaMnO$_3$, LaFeO$_3$ and SrMnO$_3$), 4$d$ (SrZrO$_3$, SrTcO$_3$ and Ca$_2$RuO$_4$) and 5$d$ (SrHfO$_3$, KTaO$_3$ and NaOsO$_3$) perovskites, covering a band gap ranging from 0.1 eV to 6.1 eV and exhibiting different electronic, structural and magnetic properties. Optical conductivities and optical transitions including electron-hole interactions are calculated through the solution of the Bethe-Salpeter equation (BSE) with quasi-particle energies evaluated by single-shot $G_0W_0$ approximation. The exciton binding energies are computed by means of a model-BSE (mBSE), carefully benchmarked against the full BSE method, in order to obtain well-converged results in terms of k-point sampling. The predicted results are compared with available measured data, with an overall satisfactory agreement between theory and experiment.
We present a hybrid approach for GW/Bethe-Salpeter Equation (BSE) calculations of core excitation spectra, including x-ray absorption (XAS), electron energy loss spectra (EELS), and non-resonant inelastic x-ray scattering (NRIXS). The method is based on {it ab initio} wavefunctions from the plane-wave pseudopotential code ABINIT; atomic core-level states and projector augmented wave (PAW) transition matrix elements; the NIST core-level BSE solver; and a many-pole GW self-energy model to account for final-state broadening and self-energy shifts. Multiplet effects are also accounted for. The approach is implemented using an interface dubbed OCEAN (Obtaining Core Excitations using ABINIT and NBSE). To demonstrate the utility of the code we present results for the K-edges in LiF as probed by XAS and NRIXS, the K-edges of KCl as probed by XAS, the Ti L_2,3-edge in SrTiO_3 as probed by XAS, and the Mg L_2,3-edge in MgO as probed by XAS. We compare the results to experiments and results obtained using other theoretical approaches.
We present a method to compute optical spectra and exciton binding energies of molecules and solids based on the solution of the Bethe-Salpeter equation (BSE) and the calculation of the screened Coulomb interaction in finite field. The method does not require the explicit evaluation of dielectric matrices nor of virtual electronic states, and can be easily applied without resorting to the random phase approximation. In addition it utilizes localized orbitals obtained from Bloch states using bisection techniques, thus greatly reducing the complexity of the calculation and enabling the efficient use of hybrid functionals to obtain single particle wavefunctions. We report exciton binding energies of several molecules and absorption spectra of condensed systems of unprecedented size, including water and ice samples with hundreds of atoms.