No Arabic abstract
We apply the quasiparticle self-consistent GW method (QSGW) to slab models of ionic materials, LiF, KF, NaCl, MgO, and CaO, under electric field. Then we obtain the optical dielectric constants E(Slab) from the differences of the slopes of the electrostatic potential in the bulk and vacuum regions. Calculated E(Slab) show very good agreements with experiments. For example, we have E(Slab)=2.91 for MgO, in agreement with the experimental value E(Experiment)=2.96. This is in contrast to E(RPA)=2.37, which is calculated in the random-phase approximation for the bulk MgO in QSGW. After we explain the difference between the quasiparticle-based perturbation theory and the Greens function based perturbation theory, we interpret the large difference E(Slab)-E(RPA)=2.91-2.37 as the contribution from the vertex correction of the proper polarization which determines the screened Coulomb interaction W. Our result encourages the theoretical development of self-consistent G0W approximation along the line of QSGW self-consistency, as was performed by Shishkin, Marsman and Kresse [Phys. Rev. Lett. 99, 246403(2007)].
We introduce a first principles approach to determine the strength of the electronic correlations based on the fully self consistent GW approximation. The approach provides a seamless interface with dynamical mean field theory, and gives good results for well studied correlated materials such as NiO. Applied to the recently discovered iron arsenide materials, it accounts for the noticeable correlation features observed in optics and photoemission while explaining the absence of visible satellites in X-ray absorption experiments and other high energy spectroscopies.
We present quasiparticle (QP) energies from fully self-consistent $GW$ (sc$GW$) calculations for a set of prototypical semiconductors and insulators within the framework of the projector-augmented wave methodology. To obtain converged results, both finite basis-set corrections and $k$-point corrections are included, and a simple procedure is suggested to deal with the singularity of the Coulomb kernel in the long-wavelength limit, the so called head correction. It is shown that the inclusion of the head corrections in the sc$GW$ calculations is critical to obtain accurate QP energies with a reasonable $k$-point set. We first validate our implementation by presenting detailed results for the selected case of diamond, and then we discuss the converged QP energies, in particular the band gaps, for a set of gapped compounds and compare them to single-shot $G_0W_0$, QP self-consistent $GW$, and previously available sc$GW$ results as well as experimental results.
We describe a finite-field approach to compute density response functions, which allows for efficient $G_0W_0$ and $G_0W_0Gamma_0$ calculations beyond the random phase approximation. The method is easily applicable to density functional calculations performed with hybrid functionals. We present results for the electronic properties of molecules and solids and we discuss a general scheme to overcome slow convergence of quasiparticle energies obtained from $G_0W_0Gamma_0$ calculations, as a function of the basis set used to represent the dielectric matrix.
Finding an accurate ab initio approach for calculating the electronic properties of transition metal oxides has been a problem for several decades. In this paper, we investigate the electronic structure of the transition metal monoxides MnO, CoO, and NiO in their undistorted rock-salt structure within a fully iterated quasiparticle self-consistent GW (QPscGW) scheme. We study the convergence of the QPscGW method, i.e., how the quasiparticle energy eigenvalues and wavefunctions converge as a function of the QPscGW iterations, and we compare the converged outputs obtained from different starting wavefunctions. We find that the convergence is slow and that a one-shot G$_0$W$_0$ calculation does not significantly improve the initial eigenvalues and states. It is important to notice that in some cases the path to convergence may go through energy band reordering which cannot be captured by the simple initial unperturbed Hamiltonian. When we reach a fully iterated solution, the converged density of states, band-gaps and magnetic moments of these oxides are found to be only weakly dependent on the choice of the starting wavefunctions and in reasonably good agreement with the experiment. Finally, this approach provides a clear picture of the interplay between the various orbitals near the Fermi level of these simple transition metal monoxides. The results of these accurate {it ab initio} calculations can provide input for models aiming at describing the low energy physics in these materials.
The Breit correction, the finite-light-speed correction for the Coulomb interaction of the electron-electron interaction in $ O left( 1/ c^2 right) $, is introduced to density functional theory (DFT) based on the non-relativistic reduction with the local density approximation. Using this newly developed relativistic DFT, it is found that the possible outer-most electron of lawrencium atom is the $ p $ orbital instead of the $ d $ orbital, which is consistent with the previous calculations based on wave-function theory. A possible explanation of the anomalous behavior of its first ionization energy is also given. This DFT scheme provides a practical calculation method for the study of properties of super-heavy elements.