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.

Linear Scaling Time-Dependent Density-Functional Tight-Binding Method for Absorption Spectra of Large Systems

A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

Efficient Algorithms for Estimating the Absorption Spectrum within Linear Response TDDFT.

Fast Algorithms for Estimating the Absorption Spectrum within the Linear Response Time-dependent Density Functional Theory

Turbo charging time-dependent density-functional theory with Lanczos chains

Density Matrix Based Time-Dependent Density Functional Theory and the Solution of Its Linear Response in Real Time Domain.

A density functional computational method based on division of large systems

Sparse-Stochastic Fragmented Exchange for Large-Scale Hybrid Time-Dependent Density Functional Theory Calculations

Sparse-Stochastic Fragmented Exchange for Large-Scale Hybrid TDDFT Calculations

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

Real-Time Propagation of the Reduced One-Electron Density Matrix in Atom-Centered Gaussian Orbitals: Application to Absorption Spectra of Silicon Clusters

Understanding intermolecular interactions of large systems in ground state and excited state by using density functional based tight binding methods

Time-Dependent Quantum Transport: An Efficient Method Based On Liouville-Von-Neumann Equation For Single-Electron Density Matrix

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

Broadband X-Ray Absorption Spectra from Time-Dependent Kohn-Sham Calculations

Accelerating Optical Absorption Spectra and Exciton Energy Computation via Interpolative Separable Density Fitting

Efficient parallel linear scaling method to get the response density matrix in all-electron real-space density-functional perturbation theory

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

Linear Scaling Algorithm for Tight-Binding Molecular Dynamics Simulations.

Broad Yet Narrow: Super-resolution techniques to simulate electronic spectra of large molecular systems

An O(N) Time-Domain Method for Time-Dependent Density Functional Theory