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.

An Efficient Algorithm for Calculating the Heat Capacity of A Large-Scale Molecular System

A Method for Calculating the Heats of Formation of Medium-Sized and Large-Sized Molecules

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

Calculating the Density of States for Large-Scale Molecular Systems

Accelerated molecular dynamics force evaluation on graphics processing units for thermal conductivity calculations.

The Heat Capacity of Triatomic Gases: An Analytical Approach

A quantum algorithm for heat conduction with symmetrization

Spectral Decomposition of Thermal Conductivity: Comparing Velocity Decomposition Methods in Homogeneous Molecular Dynamics Simulations

An enhanced version of the heat exchange algorithm with excellent energy conservation properties

Calculation of absolute molecular entropies and heat capacities made simple

Cluster-Based parallel simulation for large scale molecular dynamics in microscale thermophysics

Molecular Simulation of Caloric Properties of Fluids Modelled by Force Fields with Intramolecular Contributions: Application to Heat Capacities

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

A Concurrent Multiscale Method Based on Smoothed Molecular Dynamics for Large-Scale Parallel Computation at Finite Temperature

Homogeneous nonequilibrium molecular dynamics method for heat transport and spectral decomposition with many-body potentials

Towards Large-Scale AFQMC Calculations: Large Time Step Auxiliary-Field Quantum Monte Carlo

Finite Difference Approach for Multiscale Computations of Atomic Chain at Finite Temperature

Peridynamic modeling of transient heat conduction in solids using a highly efficient algorithm

Dgdft: A Massively Parallel Method for Large Scale Density Functional Theory Calculations

Large-Scale Normal Coordinate Analysis On Distributed Memory Parallel Systems

Implementation of the Density-functional Theory on Quantum Computers with Linear Scaling with respect to the Number of Atoms