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.

Ten-million-atom electronic structure calculations on the K computer with a massively parallel order-N theory

Novel linear algebraic theory and one-hundred-million-atom quantum material simulations on the K computer

One-hundred-nm-scale electronic structure and transport calculations of organic polymers on the K computer

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

Massively Parallel Simulation Of Large-Scale Electromagnetic Problems Using One High-Performance Computing Scheme And Domain Decomposition Method

2.5 Million-Atom Ab Initio Electronic-Structure Simulation of Complex Metallic Heterostructures with DGDFT

Parallel Implementation of Large-Scale Linear Scaling Density Functional Theory Calculations with Numerical Atomic Orbitals in HONPAS.

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

Efficient and Accurate Linear Algebraic Methods for Large-scale Electronic Structure Calculations with Non-orthogonal Atomic Orbitals

Extreme-Scale Density Functional Theory High Performance Computing of DGDFT for Tens of Thousands of Atoms using Millions of Cores on Sunway TaihuLight

Hybrid MPI and OpenMP Parallel Implementation of Large-Scale Linear-Response Time-Dependent Density Functional Theory with Plane-Wave Basis Set

Efficient and Scalable Calculation of Complex Band Structure using Sakurai-Sugiura Method

10-Million Atoms Simulation of First-Principle Package LS3DF

Ultra-large-scale electronic structure theory and numerical algorithm

Million Atom Electronic Structure and Device Calculations on Peta-Scale Computers

High performance computing of DGDFT for tens of thousands of atoms using millions of cores on Sunway TaihuLight

Enhancement of DFT-calculations at Petascale: Nuclear Magnetic Resonance, Hybrid Density Functional Theory and Car-Parrinello calculations

Dgdft: A Massively Parallel Method for Large Scale Density Functional Theory Calculations

Towards Electronic Structure-Based Ab-Initio Molecular Dynamics Simulations with Hundreds of Millions of Atoms

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

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