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.

Large-scale atomic simulation package (LSASP) :program and progress

LASP: Fast Global Potential Energy Surface Exploration

10-Million Atoms Simulation of First-Principle Package LS3DF

ls1 mardyn: The massively parallel molecular dynamics code for large systems

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

${\it Asparagus}$: A Toolkit for Autonomous, User-Guided Construction of Machine-Learned Potential Energy Surfaces

Large-Scale Atomic Simulation via Machine Learning Potentials Constructed by Global Potential Energy Surface Exploration

Lsqc: Low Scaling Quantum Chemistry Program

PASP: Property Analysis and Simulation Package for Materials

Introduction to First-Principles Simulation Package Abacus Based on Systematically Improvable Atomic Orbitals

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

DPA-2: Towards a universal large atomic model for molecular and material simulation

Development of simulation package 'ELSES' for extra-large-scale electronic-structure calculation

DP-GEN: A concurrent learning platform for the generation of reliable deep learning based potential energy models

FSATOOL 2.0: An integrated molecular dynamics simulation and trajectory data analysis program

Long-Short-Range Message-Passing: A Physics-Informed Framework to Capture Non-Local Interaction for Scalable Molecular Dynamics Simulation

Large-scale ab initio simulations based on systematically improvable atomic basis

qvasp: A flexible toolkit for VASP users in materials simulations

Amesp: Atomic and molecular electronic structure program

Atom-centered machine-learning force field package

ESPResSo++ 2.0: Advanced methods for multiscale molecular simulation