The following CPC paper is the canonical reference to use for citing LAMMPS. It gives an overview of the code including its parallel algorithms, design features, performance, and brief highlights of many of its materials modeling capabilities. If you wish, you can also mention the URL of the LAMMPS website in your paper, namely https://www.lammps.org.
LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales, A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, S. J. Plimpton, Comp Phys Comm, 271 (2022) 10817.
Here is the DOI which has a link to the PDF: https://doi.org/10.1016/j.cpc.2021.108171. And here is a BibTeX entry for the paper.
This earlier JCP paper was the original citation for LAMMPS. You can cite it if you want to refer to the parallel spatial-decomposition strategy LAMMPS still uses:
Fast Parallel Algorithms for Short-Range Molecular Dynamics, S. Plimpton, J Comp Phys, 117, 1-19 (1995). (abstract)
Below are other papers that describe specific algorithms used within LAMMPS.
Note that several capabilities in LAMMPS were implemented by folks who also wrote papers about their work. If you use their feature in your simulations and write your own paper, it would be appropriate to cite their work specifically. The list of these features is as follows. If your feature and paper is missing from this list, please let us know, and we'll add it.
|Command or package||Author||Paper to cite (see details below)|
|DIFFRACTION package||Shawn Coleman||MSMSE 2013|
|pair srp||Tim Sirk||J Chem Phys 2012|
|long-range per-atom stress||Tim Sirk||J Chem Phys 2013|
|REAXFF package and fix qeq/reax||Metin Aktulga||Parallel Computing 2012|
|PHONON package||Ling-Ti Kong||Comp Phys Comm 2011|
|EFF package||Andres Jaramillo-Botero||J Comp Chem 2011|
|GPU package||Mike Brown||Comp Phys Comm 2011 or 2012|
|fix srd||Matt Peterson||J Chem Phys 2010|
|pair gayberne||Mike Brown||J Chem Phys 2009|
|POEMS package||Rudra Mukherjee||Intl J of Non-Linear Mechanics 2008|
|PERI package for Peridynamics||Mike Parks||Comp Phys Comm 2008|
|neighbor multi||Pieter in 't Veld||Comm Phys Comp 2008|
|fix orient/fcc||Koen Janssens||Nature Materials 2006|
|fix bond/swap||Rolf Auhl||J Chem Phys 2003|
When you run LAMMPS and use one of these commands, a "log.cite" file is output, which gives a BibTeX form of the corresponding citation. You can use the -nocite or -nc command-line switch to disable output of the file.
This paper describes the virutal x-ray and electron diffraction computation as implemented in the compute xrd and compute saed commands respectively.
S. P. Coleman, D. E. Spearot, L. Capolungo, MSMSE, 21, 055020 (2013).
This paper describes the soft segmental repulsive potential (SRP) that can prevent DPD polymer chains from crossing each other, as implemented by the pair_style srp command.
An enhanced entangled polymer model for dissipative particle dynamics, T. Sirk, Y. Sliozberg, J. Brennan, M. Lisal, J. Andzelm, J Chem Phys, 136, 134903 (2012).
This paper describes the per-atom stress calculation that can now be used in LAMMPS with the long-range PPPM solver.
Characteristics of thermal conductivity in classical water models, T. W. Sirk, S. Moore, E. F. Brown, J Chem Phys, 138, 064505 (2013). (abstract)
This paper gives an overview of the many-body potentials available in LAMMPS.
Computational Aspects of Many-body Potentials, S. J. Plimpton and A. P. Thompson, MRS Bulletin, 37, 513-521 (2012). (abstract)
This paper describes work by Metin Aktulga (now at LBNL) and collaborators at Purdue to implement the ReaxFF force field in LAMMPS, available as the REAXFF package, and the fix/qeq/reaxff command.
Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques, H. M. Aktulga, J. C. Fogarty, S. A. Pandit, A. Y. Grama, Parallel Computing, 38, 245-259 (2012). (abstract)
This paper describes how the fix phonon command in LAMMPS works to compute dynamical matrices and phonon dispersion relations:
Phonon dispersion measured directly from molecular dynamics simulations, L. T. Kong, Comp Phys Comm, 182, 2201-2207 (2011). (abstract)
This paper describes work by Andres Jaramillo-Botero and collaborators at Caltech to implement their electron force field (eFF) in LAMMPS as the EFF package.
Large-Scale, Long-Term Nonadiabatic Electron Molecular Dynamics for Describing Material Properties and Phenomena in Extreme Environments, A. Jaramillo-Botero, J. Su, A. Qi, W. A. Goddard III, J Comp Chem, 32, 497-512 (2011). (abstract)
These papers describe work by Mike Brown (ORNL) to adapt several of the LAMMPS kernels (force calculation, neighbor list construction) for "hybrid" processors, meaning ones that hvae both multicore CPUs and GPUs.
Implementing Molecular Dynamics on Hybrid High Performance Computers - Short Range Forces, W. M. Brown, P. Wang, S. J. Plimpton, A. N. Tharrington, Comp Phys Comm, 182, 898-911, (2011). (abstract)
Implementing Molecular Dynamics on Hybrid High Performance Computers - Particle-Particle Particle-Mesh, W. M. Brown, A. Kohlmeyer, S. J. Plimpton, A. N. Tharrington, Comp Phys Comm, 183, 449-459 (2012). (abstract)
This paper describes the stochastic rotation dynamics (SRD) model as implemented in LAMMPS for use as a coarse-grained cheap solvent. This is for a pure SRD fluid model; a paper describing mixture systems (coarse-grained solute particles in SRD fluid) will be forthcoming.
Mesoscale Hydrodynamics via Stochastic Rotation Dynamics: Comparison with Lennard-Jones Fluid, M. K. Petersen, J. B. Lechman, S. J. Plimpton, G. S. Grest, P. J. in't Veld, P. R. Schunk, J Chem Phys, 132, 174106 (2010). (abstract)
This paper describes different ways of formulating per-atom and global virial and stress calculations, including how it is done in LAMMPS.
General formulation of pressure and stress tensor for arbitrary many-body interaction potentials under periodic boundary conditions, A. P. Thompson, S. J. Plimpton, W. Mattson, J Chem Phys, 131, 154107 (2009). (abstract)
This paper describes the ellipsoidal Gay-Berne potential in LAMMPS and its application to mixture systems.
Liquid crystal nanodroplets in solution, W. M. Brown, M. K. Petersen, S. J. Plimpton, and G. S. Grest, J Chem Phys, 130, 044901 (2009). (abstract)
This paper describes the coupling of LAMMPS to the POEMS multi-body dynamics solver, available in LAMMPS via the fix poems command.
Substructured molecular dynamics using multibody dynamics algorithms, R. M. Mukherjee, P. S. Crozier, S. J. Plimpton, K. S. Anderson, Intl J of Non-Linear Mechanics, 43, 1045-1055 (2008). (abstract)
This paper describes the implementation of Peridynamics in LAMMPS, which is a particle-based model of continuum mechanics, suitable for modeling materials at the mesoscale and macroscale.
Implementing peridynamics within a molecular dynamics code, M. L. Parks, R. B. Lehoucq, S. J. Plimpton, S. A. Silling, Comp Phys Comm, 179, 777-783 (2008). (abstract)
This paper describes the neighboring and communication algorithms developed for efficient simulation of mixtures with particles of widely varying size as implemented in the neighbor multi and communicate multi commands:
Accurate and Efficient Methods for Modeling Colloidal Mixtures in an Explicit Solvent using Molecular Dynamics, P. J. in 't Veld, S. J. Plimpton, G. S. Grest, Comp Phys Comm, 179, 320-329 (2008). (abstract)
This paper describes the grain boundary driving force methodology that allows more rapid computation of grain boundary mobility, as implemented in the fix orient/fcc command:
Computing the Mobility of Grain Boundaries, K. G. F. Janssens, D. Olmsted, E.A. Holm, S. M. Foiles, S. J. Plimpton, and P. M. Derlet, Nature Materials, 5, 124-127 (2006). (abstract)
This paper discusses the Monte Carlo bond-swapping algorithm used in LAMMPS as implemented in the fix bond/swap command:
Equilibration of long chain polymer melts in computer simulations, R. Auhl, R. Everaers, G. S. Grest, K. Kremer, S. J. Plimpton, J Chem Phys, 119, 12718-12728 (2003). (abstract)
This paper discusses the implementation of the PPPM solver (variant of PME) and rRESPA within LAMMPS, including a discussion of how the parallel FFTs work.
Particle-Mesh Ewald and rRESPA for Parallel Molecular Dynamics Simulations, S. J. Plimpton, R. Pollock, M. Stevens, in Proc of the Eighth SIAM Conference on Parallel Processing for Scientific Computing, Minneapolis, MN (March 1997). (abstract)