A divide-and-conquer/cellular-decomposition framework for million-to-billion atom simulations of chemical reactions

To enable large-scale atomistic simulations of material processes involving chemical reactions, we have designed linear-scaling molecular dynamics (MD) algorithms based on an embedded divide-and-conquer (EDC) framework: first principles-based fast reactive force-field (F-ReaxFF) MD; and quantum-mechanical MD in the framework of the density functional theory (DFT) on adaptive multigrids. To map these O(N) algorithms onto parallel computers with deep memory hierarchies, we have developed a tunable hierarchical cellular-decomposition (THCD) framework, which achieves performance tunability through a hierarchy of parameterized cell data/computation structures and adaptive load balancing through wavelet-based computational-space decomposition. Benchmark tests on 1920 Itanium2 processors of the NASA Columbia supercomputer have achieved unprecedented scales of quantum-mechanically accurate and well validated, chemically reactive atomistic simulations—0.56 billion-atom F-ReaxFF MD and 1.4 million-atom (0.12 trillion grid points) EDC–DFT MD—in addition to 18.9 billion-atom non reactive space–time multiresolution MD. The EDC and THCD frameworks expose maximal data localities, and consequently the isogranular parallel efficiency on 1920 processors is as high as 0.953. Chemically reactive MD simulations have been applied to shock-initiated detonation of energetic materials and stress-induced bond breaking in ceramics in corrosive environments.