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.

LUNAR: Automated Input Generation and Analysis for Reactive LAMMPS Simulations

Amber Free Energy Tools: Interoperable Software for Free Energy Simulations Using Generalized Quantum Mechanical/molecular Mechanical and Machine Learning Potentials.

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

CURATOR: Building Robust Machine Learning Potentials for Atomistic Simulations Autonomously with Batch Active Learning

AENET-LAMMPS and AENET-TINKER: Interfaces for Accurate and Efficient Molecular Dynamics Simulations with Machine Learning Potentials

Molecular Modeling of Reactive Systems with REACTER

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Molecular dynamics simulation of complex reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) soft- ware package

Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques

An Atomistic Machine Learning Package for Surface Science and Catalysis

Modeling chemical reactions in classical molecular dynamics simulations

Atomate: A high-level interface to generate, execute, and analyze computational materials science workflows

Active learning meets metadynamics: Automated workflow for reactive machine learning potentials

Implementing Reactivity in Molecular Dynamics Simulations with Harmonic Force Fields

PyLAT: Python LAMMPS Analysis Tools

sGDML: Constructing Accurate and Data Efficient Molecular Force Fields Using Machine Learning

Molecular Dynamics Simulations of Active Matter using LAMMPS

Catalytic Large Atomic Model (CLAM): A Machine-Learning-Based Interatomic Potential Universal Model

DLPGEN: Preparing Molecular Dynamics Simulations with Support for Polarizable Force Fields

Atom-centered machine-learning force field package