LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Commands

### pair_style lj/mdf command

### pair_style buck/mdf command

### pair_style lennard/mdf command

Syntax:

    pair_style style args

- style = lj/mdf or buck/mdf or lennard/mdf
- args = list of arguments for a particular style

        lj/mdf args = cutoff1 cutoff2
          cutoff1 = inner cutoff for the start of the tapering function
          cutoff1 = out cutoff for the end of the tapering function
        buck/mdf args = cutoff1 cutoff2
          cutoff1 = inner cutoff for the start of the tapering function
          cutoff1 = out cutoff for the end of the tapering function
        lennard/mdf args = cutoff1 cutoff2
          cutoff1 = inner cutoff for the start of the tapering function
          cutoff1 = out cutoff for the end of the tapering function

Examples:

    pair_style lj/mdf 2.5 3.0
    pair_coeff * * 1 1
    pair_coeff 1 1 1 1.1 2.8 3.0 3.2

    pair_style buck 2.5 3.0
    pair_coeff * * 100.0 1.5 200.0
    pair_coeff * * 100.0 1.5 200.0 3.0 3.5

    pair_style lennard/mdf 2.5 3.0
    pair_coeff * * 1 1
    pair_coeff 1 1 1 1.1 2.8 3.0 3.2

Description:

The lj/mdf, buck/mdf and lennard/mdf compute the standard 12-6 Lennard-Jones and Buckingham potential with the addition of a taper function that ramps the energy and force smoothly to zero between an inner and outer cutoff.

The tapering, f(r), is done by using the Mei, Davenport, Fernando function (Mei).

where

Here r_m is the inner cutoff radius and r_cut is the outer cutoff radius.

For the lj/mdf pair_style, the potential energy, E(r), is the standard 12-6 Lennard-Jones written in the epsilon/sigma form:

The following coefficients must be defined for each pair of atoms types via the pair_coeff command as in the examples above, or in the data file or restart files read by the read_data or read_restart commands, or by mixing as described below:

- epsilon (energy units)
- sigma (distance units)
- r_m (distance units)
- r_cut (distance units)

For the buck/mdf pair_style, the potential energy, E(r), is the standard Buckingham potential:

- A (energy units)
- \rho (distance units)
- C (energy-distance^6 units)
- r_m (distance units)
- r_cut$ (distance units)

For the lennard/mdf pair_style, the potential energy, E(r), is the standard 12-6 Lennard-Jones written in the $A/B$ form:

The following coefficients must be defined for each pair of atoms types via the pair_coeff command as in the examples above, or in the data file or restart files read by the read_data or read_restart commands, or by mixing as described below:

- A (energy-distance^12 units)
- B (energy-distance^6 units)
- r_m (distance units)
- r_cut (distance units)

Mixing, shift, table, tail correction, restart, rRESPA info:

For atom type pairs I,J and I != J, the epsilon and sigma coefficients and cutoff distance for all of the lj/cut pair styles can be mixed. The default mix value is geometric. See the "pair_modify" command for details.

All of the lj/cut pair styles support the pair_modify shift option for the energy of the Lennard-Jones portion of the pair interaction.

The lj/cut/coul/long and lj/cut/tip4p/long pair styles support the pair_modify table option since they can tabulate the short-range portion of the long-range Coulombic interaction.

All of the lj/cut pair styles support the pair_modify tail option for adding a long-range tail correction to the energy and pressure for the Lennard-Jones portion of the pair interaction.

All of the lj/cut pair styles write their information to binary restart files, so pair_style and pair_coeff commands do not need to be specified in an input script that reads a restart file.

The lj/cut and lj/cut/coul/long pair styles support the use of the inner, middle, and outer keywords of the run_style respa command, meaning the pairwise forces can be partitioned by distance at different levels of the rRESPA hierarchy. The other styles only support the pair keyword of run_style respa. See the run_style command for details.

Restrictions:

These pair styles can only be used if LAMMPS was built with the USER-MISC package. See the Making LAMMPS section for more info on packages.

Related commands:

pair_coeff

Default: none

(Mei) Mei, Davenport, Fernando, Phys Rev B, 43 4653 (1991)