Accelerating Atomic Orbital-based Electronic Structure Calculation via Pole Expansion and Selected Inversion

Abstract

We describe how to apply the recently developed pole expansion and selected inversion (PEXSI) technique to Kohn-Sham density function theory (DFT) electronic structure calculations that are based on atomic orbital discretization. We give analytic expressions for evaluating the charge density, the total energy, the Helmholtz free energy and the atomic forces (including both the Hellman-Feynman force and the Pulay force) without using the eigenvalues and eigenvectors of the Kohn-Sham Hamiltonian. We also show how to update the chemical potential without using Kohn-Sham eigenvalues. The advantage of using PEXSI is that it has a much lower computational complexity than that associated with the matrix diagonalization procedure. We demonstrate the performance gain by comparing the timing of PEXSI with that of diagonalization on insulating and metallic nanotubes. For these quasi-1D systems, the complexity of PEXSI is linear with respect to the number of atoms. This linear scaling can be observed in our computational experiments when the number of atoms in a nanotube is larger than a few hundreds. Both the wall clock time and the memory requirement of PEXSI is modest. This makes it even possible to perform Kohn-Sham DFT calculations for 10,000-atom nanotubes with a sequential implementation of the selected inversion algorithm. We also perform an accurate geometry optimization calculation on a truncated (8,0) boron-nitride nanotube system containing 1024 atoms. Numerical results indicate that the use of PEXSI does not lead to loss of accuracy required in a practical DFT calculation.

References in corpus (3)

• Nearsightedness of Electronic Matter
• O(N) methods in electronic structure calculations
• An efficient and accurate decomposition of the Fermi operator

Cited by in corpus (35)

• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• SIESTA: recent developments and applications
• First-principles Green's-function method for surface calculations: a pseudopotential localized basis set approach
• Challenges in Large Scale Quantum Mechanical Calculations
• ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers
• RESCU: a Real Space Electronic Structure Method
• DGDFT: A Massively Parallel Method for Large Scale Density Functional Theory Calculations
• ELSI -- An Open Infrastructure for Electronic Structure Solvers
• Electronic Structure of Large-Scale Graphene Nanoflakes
• GPU-Acceleration of the ELPA2 Distributed Eigensolver for Dense Symmetric and Hermitian Eigenproblems
• A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization
• GPGPU Acceleration of All-Electron Electronic Structure Theory Using Localized Numeric Atom-Centered Basis Functions
• Chebyshev polynomial filtered subspace iteration in the Discontinuous Galerkin method for large-scale electronic structure calculations
• Edge reconstruction in armchair phosphorene nanoribbons revealed by discontinuous Galerkin density functional theory
• Linear scaling DFT calculations for large Tungsten systems using an optimized local basis
• Efficient Computation of Sparse Matrix Functions for Large-Scale Electronic Structure Calculations: The CheSS Library
• The CECAM Electronic Structure Library and the modular software development paradigm
• Minimax rational approximation of the Fermi-Dirac distribution
• Influence of point defects on the electronic and topological properties of monolayer WTe${}_2$
• Mixed Precision Fermi-Operator Expansion on Tensor Cores From a Machine Learning Perspective
• Performance analysis of electronic structure codes on HPC systems: A case study of SIESTA
• Reduced Scaling Hilbert Space Variational Monte Carlo
• PSelInv -- A Distributed Memory Parallel Algorithm for Selected Inversion : the Symmetric Case
• Spectrum-splitting approach for Fermi-operator expansion in all-electron Kohn-Sham DFT calculations
• Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low rank approximation of the dielectric matrix
• Enhancing the scalability and load balancing of the parallel selected inversion algorithm via tree-based asynchronous communication
• Robust Determination of the Chemical Potential in the Pole Expansion and Selected Inversion Method for Solving Kohn-Sham density functional theory
• The Enduring Relevance of Semiempirical Quantum Mechanics
• An algorithm for real and complex rational minimax approximation
• Globally Constructed Adaptive Local Basis Set for Spectral Projectors of Second Order Differential Operators
• Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework II: Force, vibration, and molecular dynamics calculations
• PEXSI-$\Sigma$: A Green's function embedding method for Kohn-Sham density functional theory
• Minimax separation of the Cauchy kernel
• Incomplete selected inversion for linear-scaling electronic structure calculations
• Numerical solution of large scale Hartree-Fock-Bogoliubov equations