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.

A multi-scale code for flexible hybrid simulations

Performing highly efficient Minima Hopping structure predictions using the Atomic Simulation Environment (ASE)

Coupling finite element method with large scale atomic/molecular massively parallel simulator (LAMMPS) for hierarchical multiscale simulations

Multi-physics simulations using a hierarchical interchangeable software interface

A Hybrid Approach for Solving Coupled Maxwell and Schrödinger Equations Arising in the Simulation of Nano-Devices

Scalable and fast heterogeneous molecular simulation with predictive parallelization schemes

Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms

Efficient all-electron hybrid density functionals for atomistic simulations beyond 10 000 atoms

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

Multiscale quantum algorithms for quantum chemistry

Code modernization strategies for short-range non-bonded molecular dynamics simulations

Overcoming the chemical complexity bottleneck in on-the-fly machine learned molecular dynamics simulations

10-Million Atoms Simulation of First-Principle Package LS3DF

Accurate interatomic force field for molecular dynamics simulation by hybridizing classical and machine learning potentials

Accelerating atomistic simulations with piecewise machine-learned ab Initio potentials at a classical force field-like cost

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

Adaptive-precision potentials for large-scale atomistic simulations

Atomistic simulation of low-dimensional nanostructures toward extreme-scale supercomputing

Leveraging the Super Instruction Architecture to Develop Massively Parallel Computational and Environmental Chemistry Applications

FINETUNA: Fine-tuning Accelerated Molecular Simulations