Abstract:The Kohn-sham density functional theory (DFT) is a powerful method to describe the electronic structures of molecules and solids in condensed matter physics, computational chemistry and materials science. However, large and accurate DFT calculations within plane waves process a cubic-scaling computational complexity, which is usually limited by expensive computation and communication costs. The rapid development of high performance computing (HPC) on leadership supercomputers brings new opportunities for developing plane-wave DFT calculations for large-scale systems. Here, we implement parallel iterative eigensolvers in large-scale plane-wave DFT calculations, including Davidson, locally optimal block preconditioned conjugate gradient (LOBPCG), projected preconditioned conjugate gradient (PPCG) and the Chebyshev subspace iteration (CheFSI) algorithms, and analyze the performance of these algorithms in massively parallel plane-wave computing tasks. We adopt a two-level parallelization strategy that combines the message passing interface (MPI) with open multi-processing (OpenMP) parallel programming to handle data exchange and matrix operations in the construction and diagonalization of large-scale Hamiltonian matrix within plane waves. Numerical results illustrate that these iterative eigensolvers can scale up to 42,592 processing cores with high peak performance of 30% on leadship supercomputers to study the electronic structures of bulk silicon systems containing 10,648 atoms.Program summaryProgram Title: Plane wave density functional theory (PWDFT)CPC Library link to program files: https://doi.org/10.17632/c8v2mx5vn4.1Developer's repository link: https://bitbucket.org/berkeleylab/scalesLicensing provisions: BSD 3-clauseProgramming language: C++Nature of problem: PWDFT is used for electronic structure calculations based on Kohn-Sham density functional theory. The key challenge to address is a constrained energy minimization problem, which can also be formulated as a nonlinear eigenvalue problem. MPI/OpenMP-based approaches are employed to provide multi-core acceleration for the study of the chemical and material properties of larger-scale molecules and solids.Solution method: PWDFT implements self-consistent field (SCF) iterations and direct constrained minimization algorithms with various acceleration strategies. It is written in C++ and offers parallel acceleration based on MPI/OpenMP.

Ab Initio Electronic Structure Calculations Using a Real-Space Chebyshev-filtered Subspace Iteration Method.

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

Two-Level Chebyshev Filter Based Complementary Subspace Method: Pushing the Envelope of Large-Scale Electronic Structure Calculations.

Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration

Chebyshev Polynomial Filtered Subspace Iteration in the Discontinuous Galerkin Method for Large-Scale Electronic Structure Calculations.

Solving the electronic structure problem for over 100,000 atoms in real-space

Self-consistent Density Functional Calculations without Diagonalization

Adaptive Regularized Self-Consistent Field Iteration with Exact Hessian for Electronic Structure Calculation.

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics I: Theory, Algorithm, and Performance

ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers

Finite-element methods for noncollinear magnetism and spin-orbit coupling in real-space pseudopotential density functional theory

Introduction to First-Principles Simulation Package Abacus Based on Systematically Improvable Atomic Orbitals

Massively Parallel Implementation of Iterative Eigensolvers in Large-Scale Plane-Wave Density Functional Theory

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

Ab initio Calculations in Atoms, Molecules, and Solids, Treating Spin-Orbit Coupling and Electron Interaction on Equal Footing

Adaptive Finite Differencing in High Accuracy Electronic Structure Calculations

Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms

Computing the Self-Consistent Field in Kohn-Sham Density Functional Theory

A General Second Order Complete Active Space Self-Consistent-field Solver for Large-Scale Systems

Large-scale ab initio simulations based on systematically improvable atomic basis

Efficient and Robust Ab Initio Self-Consistent Field Acceleration Algorithm Based on a Semiempirical Model Hamiltonian