SIESTA-PEXSI: Massively parallel method for efficient and accurate \textit{ab initio} materials simulation without matrix diagonalization

Abstract

We describe a scheme for efficient large-scale electronic-structure calculations based on the combination of the pole expansion and selected inversion (PEXSI) technique with the SIESTA method, which uses numerical atomic orbitals within the Kohn-Sham density functional theory (KSDFT) framework. The PEXSI technique can efficiently utilize the sparsity pattern of the Hamiltonian and overlap matrices generated in SIESTA, and for large systems has a much lower computational complexity than that associated with the matrix diagonalization procedure. The PEXSI technique can be used to evaluate the electron density, free energy, atomic forces, density of states and local density of states without computing any eigenvalue or eigenvector of the Kohn-Sham Hamiltonian. It can achieve accuracy fully comparable to that obtained from a matrix diagonalization procedure for general systems, including metallic systems at low temperature. The PEXSI method is also highly scalable. With the recently developed massively parallel PEXSI technique, we can make efficient use of more than $10,000$ processors on high performance machines. We demonstrate the performance and accuracy of the SIESTA-PEXSI method using several examples of large scale electronic structure calculations, including 1D, 2D and bulk problems with insulating, semi-metallic, and metallic character.

References in corpus (8)

• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• Numerical atomic orbitals for linear scaling
• Nearsightedness of Electronic Matter
• O(N) methods in electronic structure calculations
• PT-Scotch: A tool for efficient parallel graph ordering
• Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework I: Total energy calculation
• Pole-based approximation of Fermi-Dirac function

Cited by in corpus (15)

• SIESTA: recent developments and applications
• Challenges in Large Scale Quantum Mechanical Calculations
• RESCU: a Real Space Electronic Structure Method
• DGDFT: A Massively Parallel Method for Large Scale Density Functional Theory Calculations
• Electronic Structure of Large-Scale Graphene Nanoflakes
• Reconstruction of moiré lattices in twisted transition metal dichalcogenide bilayers
• Chebyshev polynomial filtered subspace iteration in the Discontinuous Galerkin method for large-scale electronic structure calculations
• Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms
• Gradient-based stochastic estimation of the density matrix
• Optimizations of the Eigensolvers in the ELPA Library
• Efficient Computation of Sparse Matrix Functions for Large-Scale Electronic Structure Calculations: The CheSS Library
• Manipulation of Spin Transport in Graphene/Transition Metal Dichalcogenide Heterobilayers upon Twisting
• Influence of point defects on the electronic and topological properties of monolayer WTe${}_2$
• Robust Determination of the Chemical Potential in the Pole Expansion and Selected Inversion Method for Solving Kohn-Sham density functional theory
• Atomic relaxation and flat bands in strain-engineered transition metal dichalcogenide bilayer moiré systems