MD Simulation of Hundred-Billion-Metal-Atom Cascade Collision on Sunway Taihulight

TL;DR

This paper introduces MISA-MD, a highly efficient molecular dynamics simulation method optimized for Sunway Taihulight, capable of simulating billion-atom cascade collisions with high accuracy and scalability, advancing nuclear reactor safety research.

Contribution

MISA-MD is a novel MD implementation that employs hash-based data structures and processor-specific optimizations for large-scale, high-fidelity cascade collision simulations.

Findings

Achieves over 79% parallel efficiency on 655-billion-atom systems.

Requires less memory and outperforms LAMMPS in computational speed.

Demonstrates good accuracy and scalability for nuclear material damage modeling.

Abstract

Radiation damage to the steel material of reactor pressure vessels is a major threat to the nuclear reactor safety. It is caused by the metal atom cascade collision, initialized when the atoms are struck by a high-energy neutron. The paper presents MISA-MD, a new implementation of molecular dynamics, to simulate such cascade collision with EAM potential. MISA-MD realizes (1) a hash-based data structure to efficiently store an atom and find its neighbors, and (2) several acceleration and optimization strategies based on SW26010 processor of Sunway Taihulight supercomputer, including an efficient potential table storage and interpolation method, a coloring method to avoid write conflicts, and double-buffer and data reuse strategies. The experimental results demonstrated that MISA-MD has good accuracy and scalability, and obtains a parallel efficiency of over 79% in an 655-billion-atom system. Compared with a state-of-the-art MD program LAMMPS, MISA-MD requires less memory usage and achieves better computational performance.