Abstract:Since the classical molecular dynamics simulator LAMMPS was released as an open source code in 2004, it has become a widely-used tool for particle-based modeling of materials at length scales ranging from atomic to mesoscale to continuum. Reasons for its popularity are that it provides a wide variety of particle interaction models for different materials, that it runs on any platform from a single CPU core to the largest supercomputers with accelerators, and that it gives users control over simulation details, either via the input script or by adding code for new interatomic potentials, constraints, diagnostics, or other features needed for their models. As a result, hundreds of people have contributed new capabilities to LAMMPS and it has grown from fifty thousand lines of code in 2004 to a million lines today. In this paper several of the fundamental algorithms used in LAMMPS are described along with the design strategies which have made it flexible for both users and developers. We also highlight some capabilities recently added to the code which were enabled by this flexibility, including dynamic load balancing, on-the-fly visualization, magnetic spin dynamics models, and quantum-accuracy machine learning interatomic potentials.Program SummaryProgram Title: Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS)CPC Library link to program files: https://doi.org/10.17632/cxbxs9btsv.1Developer's repository link: https://github.com/lammps/lammpsLicensing provisions: GPLv2Programming language: C++, Python, C, FortranSupplementary material: https://www.lammps.orgNature of problem: Many science applications in physics, chemistry, materials science, and related fields require parallel, scalable, and efficient generation of long, stable classical particle dynamics trajectories. Within this common problem definition, there lies a great diversity of use cases, distinguished by different particle interaction models, external constraints, as well as timescales and lengthscales ranging from atomic to mesoscale to macroscopic.Solution method: The LAMMPS code uses neighbor lists, parallel spatial decomposition, and parallel FFTs for long-range Coulombic interactions [1]. The time integration algorithm is based on the Størmer-Verlet symplectic integrator [2], which provides better stability than higher-order non-symplectic methods. In addition, LAMMPS supports a wide range of interatomic potentials, constraints, diagnostics, software interfaces, and pre- and post-processing features.Additional comments including restrictions and unusual features: This paper serves as the definitive reference for the LAMMPS code.References[1]S. Plimpton. Fast parallel algorithms for short-range molecular dynamics. J. Comp. Phys., 117:1–19, 1995.[2]L. Verlet. Computer experiments on classical fluids: I. Thermodynamical properties of Lennard–Jones molecules. Phys. Rev., 159:98–103, 1967.

OpenMM 7: Rapid development of high performance algorithms for molecular dynamics

OpenMM 8: Molecular Dynamics Simulation with Machine Learning Potentials

An interoperable implementation of collective‐variable based enhanced sampling methods in extended phase space within the OpenMM package

Molecular dynamics simulations using the drude polarizable force field on GPUs with OpenMM: Implementation, validation, and benchmarks

drMD: Molecular Dynamics for Experimentalists

GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit

GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers

Open‐source molecular modeling software in chemical engineering focusing on the Molecular Simulation Design Framework

NNP/MM: Accelerating molecular dynamics simulations with machine learning potentials and molecular mechanic

Scalable molecular dynamics on CPU and GPU architectures with NAMD

Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters

DMFF: an Open-Source Automatic Differentiable Platform for Molecular Force Field Development and Molecular Dynamics Simulation.

Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS

GROMACS: Fast, flexible, and free

Scalable molecular dynamics with NAMD

Performance modeling of microsecond scale biological molecular dynamics simulations on heterogeneous architectures

NNP/MM: Accelerating Molecular Dynamics Simulations with Machine Learning Potentials and Molecular Mechanics

LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales

GROMACS: A message-passing parallel molecular dynamics implementation

MOIL-opt: Energy-Conserving Molecular Dynamics on a GPU/CPU System