🤖 AI Summary
Parallel computation of three-body interactions—where three molecules from distinct spatial cells participate simultaneously—poses a significant challenge in molecular dynamics simulations.
Method: This paper introduces the first linked-cell framework supporting traversal across three distinct cells. It extends conventional pairwise neighbor searching to triplet-based neighborhood construction and integrates geometric culling with distance-based pruning to drastically reduce redundant computations. The approach is validated using the Lennard-Jones fluid model under both uniform and non-uniform density conditions.
Contribution/Results: The algorithm achieves strong node-level scalability, delivering a substantial increase in molecular updates per second over baseline methods. It is the first systematic solution to the efficiency bottleneck inherent in parallel traversal of three-cell coupled neighborhoods for short-range many-body forces. By enabling scalable, high-fidelity simulation of many-body potentials, this work provides a foundational algorithmic infrastructure for next-generation, high-accuracy molecular dynamics.
📝 Abstract
In this work, algorithms for the parallel computation of three-body interactions in molecular dynamics are developed. While traversals for the computation of pair interactions are readily available in the literature, here, such traversals are extended to allow for the computation between molecules stored across three cells. A general framework for the computation of three-body interactions in linked cells is described, and then used to implement the corresponding traversals. In addition, our analysis is combined with the commonly used cutoff conditions, because they influence the total workload of the computation of interactions. The combinations between traversals and truncation conditions are validated using the well-known Lennard-Jones fluid. Validation case studies are taken from the literature and configured into homogeneous and inhomogeneous scenarios. Finally, strong scalability and performance in terms of molecule updates are measured at node-level.