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.

Optimization and Parallelization of a force field for silicon using OpenMP

General Multiobjective Force Field Optimization Framework, with Application to Reactive Force Fields for Silicon Carbide

LMFF: Efficient and Scalable Layered Materials Force Field on Heterogeneous Many-Core Processors.

10-Million Atoms Simulation of First-Principle Package LS3DF

Force field optimization by end-to-end differentiable atomistic simulation

Convergence Analysis and Performance of an Extended Central Force Optimization Algorithm

Optimization of Reactive Force Field Simulation: Refactor, Parallelization, and Vectorization for Interactions

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

Interpolative Separable Density Fitting Decomposition for Accelerating Hybrid Density Functional Calculations with Applications to Defects in Silicon

Parallel solution, adaptivity, computational convergence, and open-source code of 2d and 3d pressurized phase-field fracture problems

Global optimization of parameters in the reactive force field ReaxFF for SiOH

Searching for Globally Optimal Functional Forms for Inter-Atomic Potentials Using Parallel Tempering and Genetic Programming

The Vectorization of the Tersoff Multi-Body Potential: An Exercise in Performance Portability

Reactive force field simulations of silicon clusters

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

Structure Optimization with Stochastic Density Functional Theory

Machine Learning Force Field for Thermal Oxidation of Silicon

Parallel-in-time integration of the Landau-Lifshitz-Gilbert equation with the parallel full approximation scheme in space and time

Student Cluster Competition 2017, Team Peking University: Reproducing Vectorization of the Tersoff Multi-Body Potential on the Intel Broadwell Architecture

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