Abstract:In this work, we present a computationally efficient methodology that utilizes a local real-space formulation of the projector augmented wave (PAW) method discretized with a finite-element (FE) basis to enable accurate and large-scale electronic structure calculations. To the best of our knowledge, this is the first real-space approach for DFT calculations, combining the efficiency of PAW formalism involving smooth electronic fields with the ability of systematically improvable higher-order finite-element basis to achieve significant computational gains. In particular, we have developed efficient strategies for solving the underlying FE discretized PAW generalized eigenproblem by employing the Chebyshev filtered subspace iteration approach to compute the desired eigenspace in each self-consistent field iteration. These strategies leverage the low-rank perturbation of the FE basis overlap matrix in conjunction with reduced order quadrature rules to invert the discretized PAW overlap matrix while also exploiting the sparsity of both the local and non-local parts of the discretized PAW Hamiltonian and overlap matrices. Using the proposed approach, we benchmark the accuracy and performance on various representative examples involving periodic and non-periodic systems with plane-wave-based PAW implementations. Furthermore, we also demonstrate a considerable computational advantage ($\sim$ 5$\times$ -- 10$\times$) over state-of-the-art plane-wave methods for medium to large-scale systems ($\sim$ 6,000 -- 35,000 electrons). Finally, we show that our approach (PAW-FE) significantly reduces the degrees of freedom to achieve the desired accuracy, thereby enabling large-scale DFT simulations ($>$ 50,000 electrons) at an order of magnitude lower computational cost compared to norm-conserving pseudopotential calculations.

Real-space electronic-structure calculations with full-potential all-electron precision for transition-metals

Fast and scalable finite-element based approach for density functional theory calculations using projector-augmented wave method

First-principles Calculation of Effective Onsite Coulomb Interaction of 3d Transition Metals: Constrained Local Density Functional Approach with Maximally Localized Wannier Function

Implementation of Density Functional Embedding Theory Within the Projector-Augmented-wave Method and Applications to Semiconductor Defect States.

Improving the precision of forces in real-space pseudopotential density functional theory

An efficient first-principle approach for electronic structures calculations of nanomaterials.

Computing total energies in complex materials using charge self-consistent DFT+DMFT

Accurate Electron-phonon Interactions from Advanced Density Functional Theory

2.5 Million-Atom Ab Initio Electronic-Structure Simulation of Complex Metallic Heterostructures with DGDFT

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

Accessible computational materials design with high fidelity and high throughput

Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set

Adaptive Finite Differencing in High Accuracy Electronic Structure Calculations

ATLAS: A Real-Space Finite-Difference Implementation of Orbital-Free Density Functional Theory

Density Functional Tight-Binding Models for Band Structures of Transition-Metal Alloys and Surfaces across the d-Block

Small Basis Set Density Functional Theory Method for Cost-efficient, Large-scale Condensed Matter Simulations

First-Principles Calculation Methods for Obtaining Scattering Waves to Investigate Transport Properties of Nanostructures

An Elongation Method For First Principle Simulations Of Electronic Structures And Electron Transport Properties Of Finite Nanostructures

Efficient all-electron hybrid density functionals for atomistic simulations beyond 10 000 atoms

All-Electron, Real-Space Perturbation Theory for Homogeneous Electric Fields: Theory, Implementation, and Application within DFT

Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms