The formal relationship between ultrasoft (US) Vanderbilt-type pseudopotentials and Bl\"ochl's projector augmented wave (PAW) method is derived. It is shown that the total energy functional for US pseudopotentials can be obtained by linearization of two terms in a slightly modified PAW total energy functional. The Hamilton operator, the forces, and the stress tensor are derived for this modified PAW functional. A simple way to implement the PAW method in existing plane-wave codes supporting US pseudopotentials is pointed out. In addition, critical tests are presented to compare the accuracy and efficiency of the PAW and the US pseudopotential method with relaxed core all electron methods. These tests include small molecules $({\mathrm{H}}_{2}{,\mathrm{}\mathrm{H}}_{2}{\mathrm{O},\mathrm{}\mathrm{Li}}_{2}{,\mathrm{}\mathrm{N}}_{2}{,\mathrm{}\mathrm{F}}_{2}{,\mathrm{}\mathrm{BF}}_{3}{,\mathrm{}\mathrm{SiF}}_{4})$ and several bulk systems (diamond, Si, V, Li, Ca, ${\mathrm{CaF}}_{2},$ Fe, Co, Ni). Particular attention is paid to the bulk properties and magnetic energies of Fe, Co, and Ni.

From ultrasoft pseudopotentials to the projector augmented-wave method

During the past decade, computer simulations based on a quantum-mechanical description of the interactions between electrons and between electrons and atomic nuclei have developed an increasingly important impact on solid-state physics and chemistry and on materials science—promoting not only a deeper understanding, but also the possibility to contribute significantly to materials design for future technologies. This development is based on two important columns: (i) The improved description of electronic many-body effects within density-functional theory (DFT) and the upcoming post-DFT methods. (ii) The implementation of the new functionals and many-body techniques within highly efficient, stable, and versatile computer codes, which allow to exploit the potential of modern computer architectures. In this review, I discuss the implementation of various DFT functionals [local-density approximation (LDA), generalized gradient approximation (GGA), meta-GGA, hybrid functional mixing DFT, and exact (Hartree-Fock) exchange] and post-DFT approaches [DFT + U for strong electronic correlations in narrow bands, many-body perturbation theory (GW) for quasiparticle spectra, dynamical correlation effects via the adiabatic-connection fluctuation-dissipation theorem (AC-FDT)] in the Vienna ab initio simulation package VASP. VASP is a plane-wave all-electron code using the projector-augmented wave method to describe the electron-core interaction. The code uses fast iterative techniques for the diagonalization of the DFT Hamiltonian and allows to perform total-energy calculations and structural optimizations for systems with thousands of atoms and ab initio molecular dynamics simulations for ensembles with a few hundred atoms extending over several tens of ps. Applications in many different areas (structure and phase stability, mechanical and dynamical properties, liquids, glasses and quasicrystals, magnetism and magnetic nanostructures, semiconductors and insulators, surfaces, interfaces and thin films, chemical reactions, and catalysis) are reviewed. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2008

Ab-initio simulations of materials using VASP: Density-functional theory and beyond.

https://th.fhi-berlin.mpg.de/site/uploads/Publications/CPC-180-2175-2009.pdf

Ab Initio Molecular Simulations with Numeric Atom-Centered Orbitals

Piezoelectric materials, which convert mechanical to electrical energy (and vice versa), are crucial in medical imaging, telecommunication and ultrasonic devices. A new generation of single-crystal materials, such as Pb(Zn1/3Nb2/3)O3-PbTiO3 (PZN-PT) and Pb(Mg1/3Nb2/3)O3-PbTiO3 (PMN-PT), exhibit a piezoelectric effect that is ten times larger than conventional ceramics, and may revolutionize these applications. However, the mechanism underlying the ultrahigh performance of these new materials-and consequently the possibilities for further improvements-are not at present clear. Here we report a first-principles study of the ferroelectric perovskite, BaTiO3, which is similar to single-crystal PZN-PT but is a simpler system to analyse. We show that a large piezoelectric response can be driven by polarization rotation induced by an external electric field. Our computations suggest how to design materials with better performance, and may stimulate further interest in the fundamental theory of dielectric systems in finite electric fields.

Polarization rotation mechanism for ultrahigh electromechanical response in single-crystal piezoelectrics

Electronic structure calculations of solids using the WIEN2k package for material sciences

With its extreme accuracy and reasonable computational efficiency, the linearized augmented planewave (LAPW) method has emerged as the standard by which density functional calculations for transition metal and rare-earth containing materials are judged. This volume presents a thorough and self-conta

David J. Singh

Papers

Planewaves, Pseudopotentials, and the LAPW Method