GPAW: An open Python package for electronic structure calculations.

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

Plane wave implementation of the magnetic force theorem for magnetic exchange constants: application to bulk Fe, Co and Ni.

We present a plane wave implementation of the magnetic force theorem, which provides a first principles framework for extracting exchange constants parameterizing a classical Heisenberg model description of magnetic materials. It is shown that the full microscopic exchange tensor may be expressed in terms of the static Kohn-Sham susceptibility tensor and the exchange-correlation magnetic field. This formulation allows one to define arbitrary magnetic sites localized to predefined spatial regions, hence rendering the problem of finding Heisenberg parameters independent of any orbital decomposition of the problem. The susceptibility is calculated in a plane wave basis, which allows for systematic convergence with respect to unoccupied bands and spatial representation. We then apply the method to the well-studied problem of calculating adiabatic spin wave spectra for bulk Fe, Co and Ni, finding good agreement with previous calculations. In particular, we utilize the freedom of defining magnetic sites to show that the calculated Heisenberg parameters are robust towards changes in the definition of magnetic sites. This demonstrates that the magnetic sites can be regarded as well-defined and thus asserts the relevance of the Heisenberg model description despite the itinerant nature of the magnetic state.