Abstract:We present a computationally efficient approach to perform large-scale all-electron density functional theory calculations by enriching the classical finite element basis with compactly supported atom-centered numerical basis functions that are constructed from the solution of the Kohn-Sham (KS) problem for single atoms. We term these numerical basis functions as enrichment functions, and the resultant basis as the enriched finite element basis. The enrichment functions are compactly supported through the use of smooth cutoff functions, which enhances the conditioning and maintains the locality of the basis. The integrals involved in the evaluation of the discrete KS Hamiltonian and overlap matrix in the enriched finite element basis are computed using an adaptive quadrature grid based on the characteristics of enrichment functions. Further, we propose an efficient scheme to invert the overlap matrix by using a block-wise matrix inversion in conjunction with special reduced-order quadrature rules to transform the discrete Kohn-Sham problem to a standard eigenvalue problem. Finally, we solve the resulting standard eigenvalue problem by using a Chebyshev polynomial based filtering technique to compute the relevant eigenspectrum. We demonstrate the accuracy, efficiency and parallel scalability of the proposed method on semiconducting and heavy-metallic systems of various sizes, with the largest system containing 8694 electrons. We obtain accuracies in the ground-state energies that are within $\sim 1$mHa with ground-state energies obtained from classical finite element as well as gaussian basis sets. We observe a $50-300$ and $\sim 8$ fold reduction in the overall computational time when compared to classical finite element and gaussian basis, respectively. We also observe good parallel scalability up to 384 processors for a benchmark system comprising of 280-atom silicon nano-cluster.

Comparison Study of Finite Element and Basis Set Methods for Finite Size Scaling

Trace Distance and Scaling Behavior of A Coupled Cavity Lattice at Finite Temperature

Linear Scaling Electronic Structure Calculations With Numerical Atomic Basis Set

High-precision finite-size scaling analysis of the quantum-critical point of S=1/2 Heisenberg antiferromagnetic bilayers

Inverse Volume Scaling of Finite-Size Error in Periodic Coupled Cluster Theory

Accelerating the Convergence of Auxiliary-Field Quantum Monte Carlo inSolids with Optimized Gaussian Basis Sets

Finite Size Scaling, Fisher Zeroes and Super Yang-Mills

Understanding finite size effects in quasi-long-range orders for exactly solvable chain models

Exact scaling of geometric phase and fidelity susceptibility and their breakdown across the critical points

Scaling of Geometric Phase and Fidelity Susceptibility Across the Critical Points and Their Relations

Finite-size versus finite-temperature effects in the critical long-range $O(N)$ model

Finite-size scaling of string order parameters characterizing the Haldane phase

Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System

An evaluation of scaling factors for multiparameter scaling procedures based on DFT force fields

Systematic study of finite-size effects in quantum Monte Carlo calculations of real metallic systems

A Nonempirical Scaling Correction Approach for Density Functional Methods Involving Substantial Amount of Hartree-Fock Exchange.

Elective-momentum Ultrasize Quantum Monte Carlo Method

Auxiliary basis expansions for large-scale electronic structure calculations

Finite-size versus finite-temperature effects in the critical long-range O(N) model

Large-scale all-electron density functional theory calculations using an enriched finite element basis