Electronic structure calculations with GPAW: a real-space implementation of the projector augmented-wave method

Jussi Enkovaara, Carsten Rostgaard, Jens J. Mortensen, Juan Chen, Marcin Dułak, Lara Ferrighi, Jeppe Gavnholt, Christian Glinsvad, Ville Haikola, Heine A. Hansen, Henrik H. Kristoffersen, Mikael Kuisma, Ask H. Larsen, Lauri Lehtovaara, Mathias Ljungberg, Olga Lopez-Acevedo, Poul Georg Moses, Jussi Ojanen, Thomas Olsen, Vivien Petzold, Nichols A. Romero, Jess Stausholm-Møller, Mikkel Strange, Georgios A. Tritsaris, Marco Vanin, Michael Walter, Bjørk Hammer, Hannu Häkkinen, Georg K. H. Madsen, Risto M Nieminen, Jens K. Nørskov, Martti Puska, Tapio T. Rantala, Jakob Schiøtz, Kristian S. Thygesen, Karsten W. Jacobsen
Year of publication: 
Journal of Physics: Condensed Matter

Electronic structure calculations have become an indispensable tool in many areas of materials science and quantum chemistry. Even though the Kohn–Sham formulation of the density-functional theory (DFT) simplifies the many-body problem significantly, one is still confronted with several numerical challenges. In this article we present the projector augmented-wave (PAW) method as implemented in the GPAW program package (https://wiki. fysik.dtu.dk/gpaw) using a uniform real-space grid representation of the electronic wavefunctions. Compared to more traditional plane wave or localized basis set approaches, real-space grids offer several advantages, most notably good computational scalability and systematic convergence properties. However, as a unique feature GPAW also facilitates a localized atomic-orbital basis set in addition to the grid. The efficient atomic basis set is complementary to the more accurate grid, and the possibility to seamlessly switch between the two representations provides great flexibility. While DFT allows one to study ground state properties, time-dependent density-functional theory (TDDFT) provides access to the excited states. We have implemented the two common formulations of TDDFT, namely the linear-response and the time propagation schemes. Electron transport calculations under finite-bias conditions can be performed with GPAW using non-equilibrium Green functions and the localized basis set. In addition to the basic features of the real-space PAW method, we also describe the implementation of selected exchange–correlation functionals, parallelization schemes, SCF-method, x-ray absorption spectra, and maximally localized Wannier orbitals.