  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.
