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.

Massively Parallel Computational Chemistry with the Super Instruction Architecture and ACES4

Leveraging the Super Instruction Architecture to Develop Massively Parallel Computational and Environmental Chemistry Applications

MASSIVE: A Multidisciplinary Applications-Oriented Simulation and Visualization Environment

A Parallel, Distributed Memory Implementation of the Adaptive Sampling Configuration Interaction Method

Redesign and Accelerate the AIREBO Bond-Order Potential on the New Sunway Supercomputer

Solving the Bethe-Salpeter equation on massively parallel architectures

Massive-Parallel Implementation of the Resolution-of-Identity Coupled-Cluster Approaches in the Numeric Atom-Centered Orbital Framework for Molecular Systems

A parallel algorithm for the concurrent atomistic-continuum methodology

Combining low-cost electronic structure theory and low-cost parallel computing architecture

10-Million Atoms Simulation of First-Principle Package LS3DF

Scalable fine-grained parallelization of plane-wave-based ab initio molecular dynamics for large supercomputers

High-Scalable Collaborated Parallel Framework for Large-Scale Molecular Dynamic Simulation on Tianhe-2 Supercomputer

mAMBER: A CPU/MIC Collaborated Parallel Framework for AMBER on Tianhe-2 Supercomputer

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

Robust and effective ab initio molecular dynamics simulations on the GPU cloud infrastructure using the Schrödinger Materials Science Suite

Software-Hardware Co-Optimization for Computational Chemistry on Superconducting Quantum Processors

Massively parallel simulation of multi-physics effects in sip using high-performance computing scheme

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

SIESTA-PEXSI: Massively Parallel Method for Efficient and Accurate Ab Initio Materials Simulation Without Matrix Diagonalization.

Atomistic simulation of low-dimensional nanostructures toward extreme-scale supercomputing