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.

Solving the electronic structure problem for over 100,000 atoms in real-space

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

Two-Level Chebyshev Filter Based Complementary Subspace Method: Pushing the Envelope of Large-Scale Electronic Structure Calculations.

An Efficient DFT Solver for Nanoscale Simulations and Beyond

10-Million Atoms Simulation of First-Principle Package LS3DF

Chemical bonding in large systems using projected population analysis from real-space density functional theory calculations

Large-scale all-electron density functional theory calculations using an enriched finite element basis

Towards Very Large-Scale Electronic-Structure Calculations

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

2.5 Million-Atom Ab Initio Electronic-Structure Simulation of Complex Metallic Heterostructures with DGDFT

Towards Electronic Structure-Based Ab-Initio Molecular Dynamics Simulations with Hundreds of Millions of Atoms

Self-consistent Density Functional Calculations without Diagonalization

Parallel implementation of electronic structure eigensolver using a partitioned folded spectrum method

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

An efficient first-principle approach for electronic structures calculations of nanomaterials.

A Comprehensive Ab Initio Study of Electronic, Optical and Cohesive Properties of Silicon Quantum Dots of Various Morphologies and Sizes up to Infinity

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

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

Million Atom Electronic Structure and Device Calculations on Peta-Scale Computers

Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms

Finite-element methods for noncollinear magnetism and spin-orbit coupling in real-space pseudopotential density functional theory