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

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.

27 Pages, 10 Tables, 6 Figures

References in corpus (26)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Advanced capabilities for materials modelling with Quantum ESPRESSO
• Optimized norm-conserving Vanderbilt pseudopotentials
• The PseudoDojo: Training and grading a 85 element optimized norm-conserving pseudopotential table
• Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems
• Daubechies wavelets as a basis set for density functional pseudopotential calculations
• DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations
• Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration
• Daubechies Wavelets for Linear Scaling Density Functional Theory
• Higher-order adaptive finite-element methods for Kohn-Sham density functional theory
• SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Extended systems
• RESCU: a Real Space Electronic Structure Method
• SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Isolated clusters
• DGDFT: A Massively Parallel Method for Large Scale Density Functional Theory Calculations
• Magnetocrystalline anisotropy of FePt: a detailed view
• Large-scale all-electron density functional theory calculations using an enriched finite element basis
• Magnetocrystalline Anisotropy of Fe-based $L{1}_0$ Alloys: Validity of Approximate Methods to Treat the Spin-Orbit Interaction
• Non-collinear density functional theory
• Evaluation of Exchange-Correlation Energy, Potential, and Stress
• Configurational forces in electronic structure calculations using Kohn-Sham density functional theory
• Density functional theory study of skyrmion pinning by atomic defects in MnSi
• All-electron density functional calculations for electron and nuclear spin interactions in molecules and solids
• Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements
• High-order finite element method for atomic structure calculations
• Ionic forces and stress tensor in all-electron DFT calculations using enriched finite element basis
• Fast hardware-aware matrix-free algorithm for higher-order finite-element discretized matrix multivector products on distributed systems