No Arabic abstract
We construct a reference database of materials properties calculated using density-functional theory in the local or generalized-gradient approximation, and an all-electron or a projector augmented-wave (PAW) formulation, for verification and validation of first-principles simulations. All-electron calculations use the full-potential linearised augmented-plane wave method, as implemented in the texttt{Elk} open-source code, while PAW calculations use the datasets developed by some of us in the open-source texttt{PSlibrary} repository and the texttt{Quantum ESPRESSO} distribution. We first calculate lattice parameters, bulk moduli, and energy differences for alkaline metals, alkaline earths, and $3d$ and $4d$ transition metals in three ideal, reference phases (simple cubic, fcc, and bcc), representing a standardized crystalline monoatomic solid-state test. Then, as suggested by K. Lejaeghere {it et al.}, [Critical Reviews in Solid State and Material Sciences 39, p 1 (2014)], we compare the equations of state for all elements, except lanthanides and actinides, in their experimental phase (or occasionally a simpler, closely related one). PAW and all-electron energy differences and structural parameters agree in most cases within a few meV/atom and a fraction of a percent, respectively. This level of agreement, comparable with the previous study, includes also other PAW and all-electron data from the electronic-structure codes texttt{VASP} and texttt{WIEN2K}, and underscores the overall reliability of current, state-of-the-art electronic-structure calculations. At the same time, discrepancies that arise even within the same formulation for simple, fundamental structural properties point to the urgent need of establishing standards for verification and validation, reference data sets, and careful refinements of the computational approaches used.
We present valence electron Compton profiles calculated within the density-functional theory using the all-electron full-potential projector augmented-wave method (PAW). Our results for covalent (Si), metallic (Li, Al) and hydrogen-bonded ((H_2O)_2) systems agree well with experiments and computational results obtained with other band-structure and basis set schemes. The PAW basis set describes the high-momentum Fourier components of the valence wave functions accurately when compared with other basis set schemes and previous all-electron calculations.
We have implemented the so called GW approximation (GWA) based on an all-electron full-potential Projector Augmented Wave (PAW) method. For the screening of the Coulomb interaction W we tested three different plasmon-pole dielectric function models, and showed that the accuracy of the quasiparticle energies is not sensitive to the the details of these models. We have then applied this new method to compute the quasiparticle band structure of some small, medium and large-band-gap semiconductors: Si, GaAs, AlAs, InP, SiMg$_2$, C and (insulator) LiCl. A special attention was devoted to the convergence of the self-energy with respect to both the {bf k}-points in the Brillouin zone and to the number of reciprocal space $bf G$-vectors. The most important result is that although the all-electron GWA improves considerably the quasiparticle band structure of semiconductors, it does not always provide the correct energy band gaps as originally claimed by GWA pseudopotential type of calculations. We argue that the decoupling between the valence and core electrons is a problem, and is some what hidden in a pseudopotential type of approach.
We present a Projector Augmented-Wave~(PAW) method based on a wavelet basis set. We implemented our wavelet-PAW method as a PAW library in the ABINIT package [http://www.abinit.org] and into BigDFT [http://www.bigdft.org]. We test our implementation in prototypical systems to illustrate the potential usage of our code. By using the wavelet-PAW method, we can simulate charged and special boundary condition systems with frozen-core all-electron precision. Furthermore, our work paves the way to large-scale and potentially order-N simulations within a PAW method.
Theoretical calculations of core electron binding energies are important for aiding the interpretation of experimental core level photoelectron spectra. In previous work, the $Delta$-Self-Consistent-Field ($Delta$-SCF) method based on density functional theory has been shown to yield highly accurate 1s and 2p binding energies in free molecules. However, most experimental work is concerned with solids, not gases. In this study, we demonstrate the application of the all-electron $Delta$-SCF method to periodic systems. A consideration of the experimentally accessible points of reference leads to the definition of a core electron binding energy in a solid as the difference between the total energies of two $N-1$ electron systems: one with an explicit, localized core hole, and one with an electron removed from the highest occupied state. The calculation of each of these quantities is addressed. In addition, the analogy between a localized core hole and a charged defect in a solid is highlighted, and the extrapolation of calculated core electron binding energies to the infinite supercell limit is discussed. It is found that the extrapolated values of the core electron binding energies from periodic $Delta$-SCF calculations agree well with experimental results for both metallic and insulating systems, with a mean absolute error of 0.24 eV for the 15 core levels considered in this study.
The self-consistent evaluation of Hubbard parameters using linear-response theory is crucial for quantitatively predictive calculations based on Hubbard-corrected density-functional theory. Here, we extend a recently-introduced approach based on density-functional perturbation theory (DFPT) for the calculation of the on-site Hubbard $U$ to also compute the inter-site Hubbard $V$. DFPT allows to reduce significantly computational costs, improve numerical accuracy, and fully automate the calculation of the Hubbard parameters by recasting the linear response of a localized perturbation into an array of monochromatic perturbations that can be calculated in the primitive cell. In addition, here we generalize the entire formalism from norm-conserving to ultrasoft and projector-augmented wave formulations, and to metallic ground states. After benchmarking DFPT against the conventional real-space Hubbard linear response in a supercell, we demonstrate the effectiveness of the present extended Hubbard formulation in determining the equilibrium crystal structure of Li$_x$MnPO$_4$ (x=0,1) and the subtle energetics of Li intercalation.