A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization

Abstract

We present a subspace projection technique to conduct large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization. The proposed method treats both metallic and insulating materials in a single framework, and is applicable to both pseudopotential as well as all-electron calculations. The key ideas involved in the method include: (i) employing a higher-order spectral finite-element basis that is amenable to mesh adaption; (ii) using a Chebyshev filter to construct a subspace which is an approximation to the occupied eigenspace in a given self-consistent field iteration; (iii) using a localization procedure to construct a non-orthogonal localized basis spanning the Chebyshev filtered subspace; (iv) using a Fermi-operator expansion in terms of the subspace-projected Hamiltonian represented in the non-orthogonal localized basis to compute relevant quantities like the density matrix, electron density and band energy. We demonstrate the accuracy and efficiency of the approach on benchmark systems involving pseudopotential calculations on metallic aluminum nano-clusters up to 3430 atoms and on insulating alkane chains up to 7052 atoms, as well as all-electron calculations on silicon nano-clusters up to 3920 electrons. The benchmark studies revealed that accuracies commensurate with chemical accuracy can be obtained, and a subquadratic-scaling with system size was observed for the range of materials systems studied. In particular, for the alkane chains---close to linear-scaling is observed, whereas, for aluminum nano-clusters---the scaling is observed to be $\mathcal{O}(N^{1.46})$. For all-electron calculations on silicon nano-clusters, the scaling with the number of electrons is computed to be $\mathcal{O}(N^{1.75})$. Furthermore, significant computational savings have been realized with the proposed approach with respect to reference calculations.

25 pages, 18 figures

References in corpus (5)

• Compressed Modes for Variational Problems in Mathematics and Physics
• Recent progress with large-scale ab initio calculations: the CONQUEST code
• Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration
• All-electron density functional theory and time-dependent density functional theory with high-order finite elements
• An efficient and accurate decomposition of the Fermi operator

Cited by in corpus (21)

• DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations
• RESCU: a Real Space Electronic Structure Method
• DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization
• Large-scale all-electron density functional theory calculations using an enriched finite element basis
• Cyclic Density Functional Theory : A route to the first principles simulation of bending in nanostructures
• Electronic-structure study of an edge dislocation in Aluminum and the role of macroscopic deformations on its energetics
• Configurational forces in electronic structure calculations using Kohn-Sham density functional theory
• Electronic structure study of screw dislocation core energetics in Aluminum and core energetics informed forces in a dislocation aggregate
• Real-time time-dependent density functional theory using higher order finite element methods
• Orbital-free DFT study of the energetics of vacancy clustering and prismatic dislocation loop nucleation in aluminum
• Isogeometric analysis in electronic structure calculations
• Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements
• Spectrum-splitting approach for Fermi-operator expansion in all-electron Kohn-Sham DFT calculations
• Finite element method and isogeometric analysis in electronic structure calculations: convergence study
• Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low rank approximation of the dielectric matrix
• Treecode-accelerated Green Iteration for Kohn-Sham Density Functional Theory
• Ionic forces and stress tensor in all-electron DFT calculations using enriched finite element basis
• Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures
• Globally Constructed Adaptive Local Basis Set for Spectral Projectors of Second Order Differential Operators
• GPU Acceleration of Collinear and Noncollinear DFT Using a Numerical Atomic Orbital-Based DFT Code
• Tensor-structured algorithm for reduced-order scaling large-scale Kohn-Sham density functional theory calculations