Finite-element methods for noncollinear magnetism and spin-orbit coupling in real-space pseudopotential density functional theory

TL;DR

This paper presents an efficient finite-element method for large-scale real-space DFT calculations that incorporate noncollinear magnetism and spin-orbit coupling, offering significant speed-ups over traditional plane-wave methods.

Contribution

The authors develop a finite-element approach within the DFT-FE framework that efficiently handles noncollinear magnetism and spin-orbit coupling, filling a gap in real-space DFT methods.

Findings

Achieves 8x-11x speed-up over plane-wave methods for large systems.

Validates accuracy through comparison with plane-wave implementations.

Demonstrates computational advantages on CPUs and GPUs.

Abstract

We introduce an efficient finite-element approach for large-scale real-space pseudopotential density functional theory (DFT) calculations incorporating noncollinear magnetism and spin-orbit coupling. The approach, implemented within the open-source DFT-FE computational framework, fills a significant gap in real-space DFT calculations using finite element basis sets, which offer several advantages over traditional DFT basis sets. In particular, we leverage the local reformulation of DFT electrostatics to derive the finite-element (FE) discretized governing equations involving two-component spinors. We subsequently utilize an efficient self-consistent field iteration approach based on Chebyshev filtered subspace iteration procedure exploiting the sparsity of local and non-local parts of FE discretized Hamiltonian to solve the underlying nonlinear eigenvalue problem based on a two-grid strategy. Furthermore, we propose using a generalized functional within the framework of noncollinear magnetism and spin-orbit coupling with a stationary point at the minima of the Kohn-Sham DFT energy functional to develop a unified framework for computing atomic forces and periodic unit-cell stresses. Validation studies against plane-wave implementations show excellent agreement in ground-state energetics, vertical ionization potentials, magnetic anisotropy energies, band structures, and spin textures. The proposed method achieves up to 8x-11x speed-ups for semi-periodic and non-periodic systems with $\sim$5000-7000 electrons in terms of minimum wall times compared to widely used plane-wave implementations on CPUs in addition to exhibiting significant computational advantage on GPUs.