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.

Large Scale Plane Wave Pseudopotential Density Functional Theory Calculations on GPU Clusters

The Analysis of a Plane Wave Pseudopotential Density Functional Theory Code on a GPU Machine

Fast plane wave density functional theory molecular dynamics calculations on multi-GPU machines

First-principle nonlocal projector potential calculation on GPU cluster

The Implementation of the Three-Dimensional Unified Gas-Kinetic Wave-Particle Method on Multiple Graphics Processing Units

Implementation of the moving particle semi-implicit method for free-surface flows on GPU clusters.

Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

First Principle Calculation on Multi-GPU Machines

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

DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization

Heterogeneous Programming and Optimization of Gyrokinetic Toroidal Code and Large-Scale Performance Test on TH-1A.

An alternative GPU acceleration for a pseudopotential plane-waves density functional theory code with applications to metallic systems

Accelerating Relativistic Exact-Two-Component Density Functional Theory Calculations with Graphical Processing Units

GPU_PBTE: an Efficient Solver for Three and Four Phonon Scattering Rates on Graphics Processing Units

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

GPGPU Acceleration of All-Electron Electronic Structure Theory Using Localized Numeric Atom-Centered Basis Functions

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

On the Efficient Evaluation of the Exchange Correlation Potential on Graphics Processing Unit Clusters

Very-Large-Scale GPU-Accelerated Nuclear Gradient of Time-Dependent Density Functional Theory with Tamm-Dancoff Approximation and Range-Separated Hybrid Functionals

Parallel Transport Time-Dependent Density Functional Theory Calculations with Hybrid Functional on Summit