Abstract:Many systems of great importance in material science, chemistry, solid-state physics, and biophysics require forces generated from an electronic structure calculation, as opposed to an empirically derived force law to describe their properties adequately. The use of such forces as input to Newton's equations of motion forms the basis of the ab initio molecular dynamics method, which is able to treat the dynamics of chemical bond-breaking and -forming events. However, a very large number of electronic structure calculations must be performed to compute an ab initio molecular dynamics trajectory, making the efficiency as well as the accuracy of the electronic structure representation critical issues. One efficient and accurate electronic structure method is the generalized gradient approximation to the Kohn-Sham density functional theory implemented using a plane-wave basis set and atomic pseudopotentials. The marriage of the gradient-corrected density functional approach with molecular dynamics, as pioneered by Car and Parrinello (R. Car and M. Parrinello, Phys Rev Lett 1985, 55, 2471), has been demonstrated to be capable of elucidating the atomic scale structure and dynamics underlying many complex systems at finite temperature. However, despite the relative efficiency of this approach, it has not been possible to obtain parallel scaling of the technique beyond several hundred processors on moderately sized systems using standard approaches. Consequently, the time scales that can be accessed and the degree of phase space sampling are severely limited. To take advantage of next generation computer platforms with thousands of processors such as IBM's BlueGene, a novel scalable parallelization strategy for Car-Parrinello molecular dynamics is developed using the concept of processor virtualization as embodied by the Charm++ parallel programming system. Charm++ allows the diverse elements of a Car-Parrinello molecular dynamics calculation to be interleaved with low latency such that unprecedented scaling is achieved. As a benchmark, a system of 32 water molecules, a common system size employed in the study of the aqueous solvation and chemistry of small molecules, is shown to scale on more than 1500 processors, which is impossible to achieve using standard approaches. This degree of parallel scaling is expected to open new opportunities for scientific inquiry.

CP2K on the road to exascale

CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations

Sparse tensor based nuclear gradients for periodic Hartree-Fock and low-scaling correlated wave function methods in the CP2K software package: A massively parallel and GPU accelerated implementation

Pushing Back the Limit of Ab-initio Quantum Transport Simulations on Hybrid Supercomputers

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

Performance Analysis of CP2K Code for Ab Initio Molecular Dynamics

A Submatrix-Based Method for Approximate Matrix Function Evaluation in the Quantum Chemistry Code CP2K

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

ECP libraries and tools: An overview

Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS

Performance Analysis of CP2K Code for Ab Initio Molecular Dynamics on CPUs and GPUs

Ab initio guided atomistic modelling of nanomaterials on exascale high-performance computing platforms

Enabling particle applications for exascale computing platforms

High Performance Computing for First-Principles Kohn-Sham Density Functional Theory Towards Exascale Supercomputers

An Efficient and Accurate Car-Parrinello-like Approach to Born-Oppenheimer Molecular Dynamics

Automated Multi-platform Testing and Code Coverage Analysis of the CP2K Application

Embracing a new era of highly efficient and productive quantum Monte Carlo simulations

Achieving Performance Portability in Gaussian Basis Set Density Functional Theory on Accelerator Based Architectures in NWChemEx

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

Enabling simulation at the fifth rung of DFT: Large scale RPA calculations with excellent time to solution

Parallelization of Kinetic Theory Simulations