RESCU: a Real Space Electronic Structure Method

Abstract

In this work we present RESCU, a powerful MATLAB-based Kohn-Sham density functional theory (KS-DFT) solver. We demonstrate that RESCU can compute the electronic structure properties of systems comprising many thousands of atoms using modest computer resources, e.g. 16 to 256 cores. Its computational efficiency is achieved from exploiting four routes. First, we use numerical atomic orbital (NAO) techniques to efficiently generate a good quality initial subspace which is crucially required by Chebyshev filtering methods. Second, we exploit the fact that only a subspace spanning the occupied Kohn-Sham states is required, and solving accurately the KS equation using eigensolvers can generally be avoided. Third, by judiciously analyzing and optimizing various parts of the procedure in RESCU, we delay the $O(N^3)$ scaling to large $N$, and our tests show that RESCU scales consistently as $O(N^{2.3})$ from a few hundred atoms to more than 5,000 atoms when using a real space grid discretization. The scaling is better or comparable in a NAO basis up to the 14,000 atoms level. Fourth, we exploit various numerical algorithms and, in particular, we introduce a partial Rayleigh-Ritz algorithm to achieve efficiency gains for systems comprising more than 10,000 electrons. We demonstrate the power of RESCU in solving KS-DFT problems using many examples running on 16, 64 and/or 256 cores: a 5,832 Si atoms supercell; a 8,788 Al atoms supercell; a 5,324 Cu atoms supercell and a small DNA molecule submerged in 1,713 water molecules for a total 5,399 atoms. The KS-DFT is entirely converged in a few hours in all cases. Our results suggest that the RESCU method has reached a milestone of solving thousands of atoms by KS-DFT on a modest computer cluster.

19 pages, 11 figures

References in corpus (6)

• Libxc: a library of exchange and correlation functionals for density functional theory
• Recent progress with large-scale ab initio calculations: the CONQUEST code
• Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration
• A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization
• Parallel eigensolvers in plane-wave Density Functional Theory
• SIESTA-PEXSI: Massively parallel method for efficient and accurate \textit{ab initio} materials simulation without matrix diagonalization

Cited by in corpus (16)

• DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations
• DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization
• Intriguing Magnetoelectric Effect in Two-dimensional Ferromagnetic/Perovskite Oxide Ferroelectric Heterostructure
• Chebyshev polynomial filtered subspace iteration in the Discontinuous Galerkin method for large-scale electronic structure calculations
• Meta-GGA density functional calculations on atoms with spherically symmetric densities in the finite element formalism
• Valley filtering effect of phonons in graphene with a grain boundary
• Electronic Structure and Band Gap Engineering of Two-Dimensional Octagon-Nitrogene
• Thermal engineering in low-dimensional quantum devices: a tutorial review of nonequilibrium Green's function methods
• Stabilities and novel electronic structures of three carbon nitride bilayers
• Solving the electronic structure problem for over 100,000 atoms in real-space
• Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low rank approximation of the dielectric matrix
• A Time-Dependent Random State Approach for Large-scale Density Functional Calculations
• A two-density approach to the general many-body problem and a proof of principle for small atoms and molecules
• Strain-induced topological phase transition in two-dimensional platinum ditelluride
• Finite-element methods for noncollinear magnetism and spin-orbit coupling in real-space pseudopotential density functional theory
• Three-dimensional flat bands and possible interlayer triplet pairing superconductivity in the alternating twisted NbSe${}_2$ moiré bulk