pair_style command
Syntax
pair_style style args
- style = one of the styles from the list below
- args = arguments used by a particular style
Examples
pair_style lj/cut 2.5
pair_style eam/alloy
pair_style hybrid lj/charmm/coul/long 10.0 eam
pair_style table linear 1000
pair_style none
Description
Set the formula(s) LAMMPS uses to compute pairwise interactions. In LAMMPS, pair potentials are defined between pairs of atoms that are within a cutoff distance and the set of active interactions typically changes over time. See the bond_style command to define potentials between pairs of bonded atoms, which typically remain in place for the duration of a simulation.
In LAMMPS, pairwise force fields encompass a variety of interactions, some of which include many-body effects, e.g. EAM, Stillinger-Weber, Tersoff, REBO potentials. They are still classified as “pairwise” potentials because the set of interacting atoms changes with time (unlike molecular bonds) and thus a neighbor list is used to find nearby interacting atoms.
Hybrid models where specified pairs of atom types interact via different pair potentials can be setup using the hybrid pair style.
The coefficients associated with a pair style are typically set for each pair of atom types, and are specified by the pair_coeff command or read from a file by the read_data or read_restart commands.
The pair_modify command sets options for mixing of type I-J interaction coefficients and adding energy offsets or tail corrections to Lennard-Jones potentials. Details on these options as they pertain to individual potentials are described on the doc page for the potential. Likewise, info on whether the potential information is stored in a restart file is listed on the potential doc page.
In the formulas listed for each pair style, E is the energy of a pairwise interaction between two atoms separated by a distance r. The force between the atoms is the negative derivative of this expression.
If the pair_style command has a cutoff argument, it sets global cutoffs for all pairs of atom types. The distance(s) can be smaller or larger than the dimensions of the simulation box.
Typically, the global cutoff value can be overridden for a specific pair of atom types by the pair_coeff command. The pair style settings (including global cutoffs) can be changed by a subsequent pair_style command using the same style. This will reset the cutoffs for all atom type pairs, including those previously set explicitly by a pair_coeff command. The exceptions to this are that pair_style table and hybrid settings cannot be reset. A new pair_style command for these styles will wipe out all previously specified pair_coeff values.
Here is an alphabetic list of pair styles defined in LAMMPS. They are also listed in more compact form on the Commands pair doc page.
Click on the style to display the formula it computes, any additional arguments specified in the pair_style command, and coefficients specified by the associated pair_coeff command.
There are also additional accelerated pair styles included in the LAMMPS distribution for faster performance on CPUs, GPUs, and KNLs. The individual style names on the Commands pair doc page are followed by one or more of (g,i,k,o,t) to indicate which accelerated styles exist.
- none - turn off pairwise interactions
- hybrid - multiple styles of pairwise interactions
- hybrid/overlay - multiple styles of superposed pairwise interactions
- zero - neighbor list but no interactions
- adp - angular dependent potential (ADP) of Mishin
- agni - machine learned potential mapping atomic environment to forces
- airebo - AIREBO potential of Stuart
- airebo/morse - AIREBO with Morse instead of LJ
- atm - Axilrod-Teller-Muto potential
- awpmd/cut - Antisymmetrized Wave Packet MD potential for atoms and electrons
- beck - Beck potential
- body/nparticle - interactions between body particles
- body/rounded/polygon - granular-style 2d polygon potential
- body/rounded/polyhedron - granular-style 3d polyhedron potential
- bop - BOP potential of Pettifor
- born - Born-Mayer-Huggins potential
- born/coul/dsf - Born with damped-shifted-force model
- born/coul/dsf/cs - Born with damped-shifted-force and core/shell model
- born/coul/long - Born with long-range Coulombics
- born/coul/long/cs - Born with long-range Coulombics and core/shell
- born/coul/msm - Born with long-range MSM Coulombics
- born/coul/wolf - Born with Wolf potential for Coulombics
- born/coul/wolf/cs - Born with Wolf potential for Coulombics and core/shell model
- brownian - Brownian potential for Fast Lubrication Dynamics
- brownian/poly - Brownian potential for Fast Lubrication Dynamics with polydispersity
- buck - Buckingham potential
- buck/coul/cut - Buckingham with cutoff Coulomb
- buck/coul/long - Buckingham with long-range Coulombics
- buck/coul/long/cs - Buckingham with long-range Coulombics and core/shell
- buck/coul/msm - Buckingham with long-range MSM Coulombics
- buck/long/coul/long - long-range Buckingham with long-range Coulombics
- buck/mdf - Buckingham with a taper function
- buck6d/coul/gauss/dsf - dispersion-damped Buckingham with damped-shift-force model
- buck6d/coul/gauss/long - dispersion-damped Buckingham with long-range Coulombics
- colloid - integrated colloidal potential
- comb - charge-optimized many-body (COMB) potential
- comb3 - charge-optimized many-body (COMB3) potential
- cosine/squared - Cooke-Kremer-Deserno membrane model potential
- coul/cut - cutoff Coulombic potential
- coul/cut/soft - Coulombic potential with a soft core
- coul/debye - cutoff Coulombic potential with Debye screening
- coul/diel - Coulomb potential with dielectric permittivity
- coul/dsf - Coulombics with damped-shifted-force model
- coul/long - long-range Coulombic potential
- coul/long/cs - long-range Coulombic potential and core/shell
- coul/long/soft - long-range Coulombic potential with a soft core
- coul/msm - long-range MSM Coulombics
- coul/shield - Coulombics for boron nitride for use with ilp/graphene/hbn potential
- coul/streitz - Coulombics via Streitz/Mintmire Slater orbitals
- coul/wolf - Coulombics via Wolf potential
- coul/wolf/cs - ditto with core/shell adjustments
- dpd - dissipative particle dynamics (DPD)
- dpd/fdt - DPD for constant temperature and pressure
- dpd/fdt/energy - DPD for constant energy and enthalpy
- dpd/tstat - pair-wise DPD thermostatting
- dsmc - Direct Simulation Monte Carlo (DSMC)
- e3b - Explicit-three body (E3B) water model
- drip - Dihedral-angle-corrected registry-dependent interlayer potential (DRIP)
- eam - embedded atom method (EAM)
- eam/alloy - alloy EAM
- eam/cd - concentration-dependent EAM
- eam/cd/old - older two-site model for concentration-dependent EAM
- eam/fs - Finnis-Sinclair EAM
- edip - three-body EDIP potential
- edip/multi - multi-element EDIP potential
- edpd - eDPD particle interactions
- eff/cut - electron force field with a cutoff
- eim - embedded ion method (EIM)
- exp6/rx - reactive DPD potential
- extep - extended Tersoff potential
- gauss - Gaussian potential
- gauss/cut - generalized Gaussian potential
- gayberne - Gay-Berne ellipsoidal potential
- granular - Generalized granular potential
- gran/hertz/history - granular potential with Hertzian interactions
- gran/hooke - granular potential with history effects
- gran/hooke/history - granular potential without history effects
- gw - Gao-Weber potential
- gw/zbl - Gao-Weber potential with a repulsive ZBL core
- hbond/dreiding/lj - DREIDING hydrogen bonding LJ potential
- hbond/dreiding/morse - DREIDING hydrogen bonding Morse potential
- ilp/graphene/hbn - registry-dependent interlayer potential (ILP)
- kim - interface to potentials provided by KIM project
- kolmogorov/crespi/full - Kolmogorov-Crespi (KC) potential with no simplifications
- kolmogorov/crespi/z - Kolmogorov-Crespi (KC) potential with normals along z-axis
- lcbop - long-range bond-order potential (LCBOP)
- lebedeva/z - Lebedeva interlayer potential for graphene with normals along z-axis
- lennard/mdf - LJ potential in A/B form with a taper function
- line/lj - LJ potential between line segments
- list - potential between pairs of atoms explicitly listed in an input file
- lj/charmm/coul/charmm - CHARMM potential with cutoff Coulomb
- lj/charmm/coul/charmm/implicit - CHARMM for implicit solvent
- lj/charmm/coul/long - CHARMM with long-range Coulomb
- lj/charmm/coul/long/soft - CHARMM with long-range Coulomb and a soft core
- lj/charmm/coul/msm - CHARMM with long-range MSM Coulombics
- lj/charmmfsw/coul/charmmfsh - CHARMM with force switching and shifting
- lj/charmmfsw/coul/long - CHARMM with force switching and long-rnage Coulombics
- lj/class2 - COMPASS (class 2) force field with no Coulomb
- lj/class2/coul/cut - COMPASS with cutoff Coulomb
- lj/class2/coul/cut/soft - COMPASS with cutoff Coulomb with a soft core
- lj/class2/coul/long - COMPASS with long-range Coulomb
- lj/class2/coul/long/soft - COMPASS with long-range Coulomb with a soft core
- lj/class2/soft - COMPASS (class 2) force field with no Coulomb with a soft core
- lj/cubic - LJ with cubic after inflection point
- lj/cut - cutoff Lennard-Jones potential with no Coulomb
- lj/cut/coul/cut - LJ with cutoff Coulomb
- lj/cut/coul/cut/soft - LJ with cutoff Coulomb with a soft core
- lj/cut/coul/debye - LJ with Debye screening added to Coulomb
- lj/cut/coul/dsf - LJ with Coulombics via damped shifted forces
- lj/cut/coul/long - LJ with long-range Coulombics
- lj/cut/coul/long/cs - ditto with core/shell adjustments
- lj/cut/coul/long/soft - LJ with long-range Coulombics with a soft core
- lj/cut/coul/msm - LJ with long-range MSM Coulombics
- lj/cut/coul/wolf - LJ with Coulombics via Wolf potential
- lj/cut/dipole/cut - point dipoles with cutoff
- lj/cut/dipole/long - point dipoles with long-range Ewald
- lj/cut/soft - LJ with a soft core
- lj/cut/thole/long - LJ with Coulombics with thole damping
- lj/cut/tip4p/cut - LJ with cutoff Coulomb for TIP4P water
- lj/cut/tip4p/long - LJ with long-range Coulomb for TIP4P water
- lj/cut/tip4p/long/soft - LJ with cutoff Coulomb for TIP4P water with a soft core
- lj/expand - Lennard-Jones for variable size particles
- lj/expand/coul/long - Lennard-Jones for variable size particles with long-range Coulombics
- lj/gromacs - GROMACS-style Lennard-Jones potential
- lj/gromacs/coul/gromacs - GROMACS-style LJ and Coulombic potential
- lj/long/coul/long - long-range LJ and long-range Coulombics
- lj/long/dipole/long - long-range LJ and long-range point dipoles
- lj/long/tip4p/long - long-range LJ and long-range Coulombics for TIP4P water
- lj/mdf - LJ potential with a taper function
- lj/sdk - LJ for SDK coarse-graining
- lj/sdk/coul/long - LJ for SDK coarse-graining with long-range Coulombics
- lj/sdk/coul/msm - LJ for SDK coarse-graining with long-range Coulombics via MSM
- lj/sf/dipole/sf - LJ with dipole interaction with shifted forces
- lj/smooth - smoothed Lennard-Jones potential
- lj/smooth/linear - linear smoothed LJ potential
- lj/switch3/coulgauss/long - smoothed LJ vdW potential with Gaussian electrostatics
- lj96/cut - Lennard-Jones 9/6 potential
- local/density - generalized basic local density potential
- lubricate - hydrodynamic lubrication forces
- lubricate/poly - hydrodynamic lubrication forces with polydispersity
- lubricateU - hydrodynamic lubrication forces for Fast Lubrication Dynamics
- lubricateU/poly - hydrodynamic lubrication forces for Fast Lubrication with polydispersity
- mdpd - mDPD particle interactions
- mdpd/rhosum - mDPD particle interactions for mass density
- meam/c - modified embedded atom method (MEAM) in C
- meam/spline - splined version of MEAM
- meam/sw/spline - splined version of MEAM with a Stillinger-Weber term
- mesocnt - mesoscale model for (carbon) nanotubes
- mgpt - simplified model generalized pseudopotential theory (MGPT) potential
- mie/cut - Mie potential
- mm3/switch3/coulgauss/long - smoothed MM3 vdW potential with Gaussian electrostatics
- momb - Many-Body Metal-Organic (MOMB) force field
- morse - Morse potential
- morse/smooth/linear - linear smoothed Morse potential
- morse/soft - Morse potential with a soft core
- multi/lucy - DPD potential with density-dependent force
- multi/lucy/rx - reactive DPD potential with density-dependent force
- nb3b/harmonic - non-bonded 3-body harmonic potential
- nm/cut - N-M potential
- nm/cut/coul/cut - N-M potential with cutoff Coulomb
- nm/cut/coul/long - N-M potential with long-range Coulombics
- oxdna/coaxstk -
- oxdna/excv -
- oxdna/hbond -
- oxdna/stk -
- oxdna/xstk -
- oxdna2/coaxstk -
- oxdna2/dh -
- oxdna2/excv -
- oxdna2/hbond -
- oxdna2/stk -
- oxdna2/xstk -
- oxrna2/coaxstk -
- oxrna2/dh -
- oxrna2/excv -
- oxrna2/hbond -
- oxrna2/stk -
- oxrna2/xstk -
- peri/eps - peridynamic EPS potential
- peri/lps - peridynamic LPS potential
- peri/pmb - peridynamic PMB potential
- peri/ves - peridynamic VES potential
- polymorphic - polymorphic 3-body potential
- python -
- quip -
- reax/c - ReaxFF potential in C
- rebo - 2nd generation REBO potential of Brenner
- resquared - Everaers RE-Squared ellipsoidal potential
- sdpd/taitwater/isothermal - smoothed dissipative particle dynamics for water at isothermal conditions
- smd/hertz -
- smd/tlsph -
- smd/tri_surface -
- smd/ulsph -
- smtbq -
- snap - SNAP quantum-accurate potential
- soft - Soft (cosine) potential
- sph/heatconduction -
- sph/idealgas -
- sph/lj -
- sph/rhosum -
- sph/taitwater -
- sph/taitwater/morris -
- spin/dipole/cut -
- spin/dipole/long -
- spin/dmi -
- spin/exchange -
- spin/magelec -
- spin/neel -
- srp -
- sw - Stillinger-Weber 3-body potential
- table - tabulated pair potential
- table/rx -
- tdpd - tDPD particle interactions
- tersoff - Tersoff 3-body potential
- tersoff/mod - modified Tersoff 3-body potential
- tersoff/mod/c -
- tersoff/table -
- tersoff/zbl - Tersoff/ZBL 3-body potential
- thole - Coulomb interactions with thole damping
- tip4p/cut - Coulomb for TIP4P water w/out LJ
- tip4p/long - long-range Coulombics for TIP4P water w/out LJ
- tip4p/long/soft -
- tri/lj - LJ potential between triangles
- ufm -
- vashishta - Vashishta 2-body and 3-body potential
- vashishta/table -
- yukawa - Yukawa potential
- yukawa/colloid - screened Yukawa potential for finite-size particles
- zbl - Ziegler-Biersack-Littmark potential
Restrictions
This command must be used before any coefficients are set by the pair_coeff, read_data, or read_restart commands.
Some pair styles are part of specific packages. They are only enabled if LAMMPS was built with that package. See the Build package doc page for more info. The doc pages for individual pair potentials tell if it is part of a package.
Default
pair_style none