pair_style lj/cut/dipole/cut command

Accelerator Variants: lj/cut/dipole/cut/gpu, lj/cut/dipole/cut/omp

pair_style lj/sf/dipole/sf command

Accelerator Variants: lj/sf/dipole/sf/gpu, lj/sf/dipole/sf/omp

pair_style lj/cut/dipole/long command

Accelerator Variants: lj/cut/dipole/long/gpu

pair_style lj/long/dipole/long command

Syntax

pair_style lj/cut/dipole/cut cutoff (cutoff2)
pair_style lj/sf/dipole/sf cutoff (cutoff2)
pair_style lj/cut/dipole/long cutoff (cutoff2)
pair_style lj/long/dipole/long flag_lj flag_coul cutoff (cutoff2)

  • cutoff = global cutoff LJ (and Coulombic if only 1 arg) (distance units)

  • cutoff2 = global cutoff for Coulombic and dipole (optional) (distance units)

  • flag_lj = long or cut or off

    long = use long-range damping on dispersion 1/r^6 term
cut = use a cutoff on dispersion 1/r^6 term
off = omit disperion 1/r^6 term entirely

  • flag_coul = long or off

    long = use long-range damping on Coulombic 1/r and point-dipole terms
off = omit Coulombic and point-dipole terms entirely

Examples

pair_style lj/cut/dipole/cut 10.0
pair_coeff * * 1.0 1.0
pair_coeff 2 3 1.0 1.0 2.5 4.0

pair_style lj/sf/dipole/sf 9.0
pair_coeff * * 1.0 1.0
pair_coeff 2 3 1.0 1.0 2.5 4.0 scale 0.5
pair_coeff 2 3 1.0 1.0 2.5 4.0

pair_style lj/cut/dipole/long 10.0
pair_coeff * * 1.0 1.0
pair_coeff 2 3 1.0 1.0 2.5 4.0

pair_style lj/long/dipole/long long long 3.5 10.0
pair_coeff * * 1.0 1.0
pair_coeff 2 3 1.0 1.0 2.5 4.0

Description

Style lj/cut/dipole/cut computes interactions between pairs of particles that each have a charge and/or a point dipole moment. In addition to the usual Lennard-Jones interaction between the particles (Elj) the charge-charge (Eqq), charge-dipole (Eqp), and dipole-dipole (Epp) interactions are computed by these formulas for the energy (E), force (F), and torque (T) between particles I and J.

\[\begin{split}E_{LJ} = & 4 \epsilon \left[ \left(\frac{\sigma}{r}\right)^{12} - \left(\frac{\sigma}{r}\right)^6 \right] \\ E_{qq} = & \frac{q_i q_j}{r} \\ E_{qp} = & \frac{q}{r^3} (p \bullet \vec{r}) \\ E_{pp} = & \frac{1}{r^3} (\vec{p_i} \bullet \vec{p_j}) - \frac{3}{r^5} (\vec{p_i} \bullet \vec{r}) (\vec{p_j} \bullet \vec{r}) \\ & \\ F_{qq} = & \frac{q_i q_j}{r^3} \vec{r} \\ F_{qp} = & -\frac{q}{r^3} \vec{p} + \frac{3q}{r^5} (\vec{p} \bullet \vec{r}) \vec{r} \\ F_{pp} = & \frac{3}{r^5} (\vec{p_i} \bullet \vec{p_j}) \vec{r} - \frac{15}{r^7} (\vec{p_i} \bullet \vec{r}) (\vec{p_j} \bullet \vec{r}) \vec{r} + \frac{3}{r^5} \left[ (\vec{p_j} \bullet \vec{r}) \vec{p_i} + (\vec{p_i} \bullet \vec{r}) \vec{p_j} \right] \\ & \\ T_{pq} = T_{ij} = & \frac{q_j}{r^3} (\vec{p_i} \times \vec{r}) \\ T_{qp} = T_{ji} = & - \frac{q_i}{r^3} (\vec{p_j} \times \vec{r}) \\ T_{pp} = T_{ij} = & -\frac{1}{r^3} (\vec{p_i} \times \vec{p_j}) + \frac{3}{r^5} (\vec{p_j} \bullet \vec{r}) (\vec{p_i} \times \vec{r}) \\ T_{pp} = T_{ji} = & -\frac{1}{r^3} (\vec{p_j} \times \vec{p_i}) + \frac{3}{r^5} (\vec{p_i} \bullet \vec{r}) (\vec{p_j} \times \vec{r})\end{split}\]

where \(q_i\) and \(q_j\) are the charges on the two particles, \(\vec{p_i}\) and \(\vec{p_j}\) are the dipole moment vectors of the two particles, r is their separation distance, and the vector r = Ri - Rj is the separation vector between the two particles. Note that Eqq and Fqq are simply Coulombic energy and force, Fij = -Fji as symmetric forces, and Tij != -Tji since the torques do not act symmetrically. These formulas are discussed in (Allen) and in (Toukmaji).

Also note, that in the code, all of these terms (except Elj) have a \(C/\epsilon\) prefactor, the same as the Coulombic term in the LJ + Coulombic pair styles discussed here. C is an energy-conversion constant and epsilon is the dielectric constant which can be set by the dielectric command. The same is true of the equations that follow for other dipole pair styles.

Style lj/sf/dipole/sf computes “shifted-force” interactions between pairs of particles that each have a charge and/or a point dipole moment. In general, a shifted-force potential is a (slightly) modified potential containing extra terms that make both the energy and its derivative go to zero at the cutoff distance; this removes (cutoff-related) problems in energy conservation and any numerical instability in the equations of motion (Allen). Shifted-force interactions for the Lennard-Jones (E_LJ), charge-charge (Eqq), charge-dipole (Eqp), dipole-charge (Epq) and dipole-dipole (Epp) potentials are computed by these formulas for the energy (E), force (F), and torque (T) between particles I and J:

\[ \begin{align}\begin{aligned}\begin{split} E_{LJ} = & 4\epsilon \left\{ \left[ \left( \frac{\sigma}{r} \right)^{\!12} - \left( \frac{\sigma}{r} \right)^{\!6} \right] + \left[ 6\left( \frac{\sigma}{r_c} \right)^{\!12} - 3\left(\frac{\sigma}{r_c}\right)^{\!6}\right]\left(\frac{r}{r_c}\right)^{\!2} - 7\left( \frac{\sigma}{r_c} \right)^{\!12} + 4\left( \frac{\sigma}{r_c} \right)^{\!6}\right\} \\ E_{qq} = & \frac{q_i q_j}{r}\left(1-\frac{r}{r_c}\right)^{\!2} \\ E_{pq} = & E_{ji} = -\frac{q}{r^3} \left[ 1 - 3\left(\frac{r}{r_c}\right)^{\!2} + 2\left(\frac{r}{r_c}\right)^{\!3}\right] (\vec{p}\bullet\vec{r}) \\ E_{qp} = & E_{ij} = \frac{q}{r^3} \left[ 1 - 3\left(\frac{r}{r_c}\right)^{\!2} + 2\left(\frac{r}{r_c}\right)^{\!3}\right] (\vec{p}\bullet\vec{r}) \\ E_{pp} = & \left[1-4\left(\frac{r}{r_c}\right)^{\!3} + 3\left(\frac{r}{r_c}\right)^{\!4}\right]\left[\frac{1}{r^3} (\vec{p_i} \bullet \vec{p_j}) - \frac{3}{r^5} (\vec{p_i} \bullet \vec{r}) (\vec{p_j} \bullet \vec{r})\right] \\ & \\\end{split}\\\begin{split}F_{LJ} = & \left\{\left[48\epsilon \left(\frac{\sigma}{r}\right)^{\!12} - 24\epsilon \left(\frac{\sigma}{r}\right)^{\!6} \right]\frac{1}{r^2} - \left[48\epsilon \left(\frac{\sigma}{r_c}\right)^{\!12} - 24\epsilon \left(\frac{\sigma}{r_c}\right)^{\!6} \right]\frac{1}{r_c^2}\right\}\vec{r}\\ F_{qq} = & \frac{q_i q_j}{r}\left(\frac{1}{r^2} - \frac{1}{r_c^2}\right)\vec{r} \\ F_{pq} = & F_{ij } = -\frac{3q}{r^5} \left[ 1 - \left(\frac{r}{r_c}\right)^{\!2}\right](\vec{p}\bullet\vec{r})\vec{r} + \frac{q}{r^3}\left[1-3\left(\frac{r}{r_c}\right)^{\!2} + 2\left(\frac{r}{r_c}\right)^{\!3}\right] \vec{p} \\ F_{qp} = & F_{ij} = \frac{3q}{r^5} \left[ 1 - \left(\frac{r}{r_c}\right)^{\!2}\right] (\vec{p}\bullet\vec{r})\vec{r} - \frac{q}{r^3}\left[1-3\left(\frac{r}{r_c}\right)^{\!2} + 2\left(\frac{r}{r_c}\right)^{\!3}\right] \vec{p} \\ F_{pp} = &\frac{3}{r^5}\Bigg\{\left[1-\left(\frac{r}{r_c}\right)^{\!4}\right] \left[(\vec{p_i}\bullet\vec{p_j}) - \frac{3}{r^2} (\vec{p_i}\bullet\vec{r}) (\vec{p_j} \bullet \vec{r})\right] \vec{r} + \\ & \left[1 - 4\left(\frac{r}{r_c}\right)^{\!3}+3\left(\frac{r}{r_c}\right)^{\!4}\right] \left[ (\vec{p_j} \bullet \vec{r}) \vec{p_i} + (\vec{p_i} \bullet \vec{r}) \vec{p_j} -\frac{2}{r^2} (\vec{p_i} \bullet \vec{r}) (\vec{p_j} \bullet \vec{r})\vec{r}\right] \Bigg\}\end{split}\end{aligned}\end{align} \]
\[\begin{split} T_{pq} = T_{ij} = & \frac{q_j}{r^3} \left[ 1 - 3\left(\frac{r}{r_c}\right)^{\!2} + 2\left(\frac{r}{r_c}\right)^{\!3}\right] (\vec{p_i}\times\vec{r}) \\ T_{qp} = T_{ji} = & - \frac{q_i}{r^3} \left[ 1 - 3\left(\frac{r}{r_c}\right)^{\!2} + 2\left(\frac{r}{r_c}\right)^{\!3} \right] (\vec{p_j}\times\vec{r}) \\ T_{pp} = T_{ij} = & -\frac{1}{r^3}\left[1-4\left(\frac{r}{r_c}\right)^{\!3} + e3\left(\frac{r}{r_c}\right)^{\!4}\right] (\vec{p_i} \times \vec{p_j}) + \\ & \frac{3}{r^5}\left[1-4\left(\frac{r}{r_c}\right)^{\!3} + 3\left(\frac{r}{r_c}\right)^{\!4}\right] (\vec{p_j}\bullet\vec{r}) (\vec{p_i} \times \vec{r}) \\ T_{pp} = T_{ji} = & -\frac{1}{r^3}\left[1-4\left(\frac{r}{r_c}\right)^{\!3} + 3\left(\frac{r}{r_c}\right)^{\!4}\right](\vec{p_j} \times \vec{p_i}) + \\ & \frac{3}{r^5}\left[1-4\left(\frac{r}{r_c}\right)^{\!3} + 3\left(\frac{r}{r_c}\right)^{\!4}\right] (\vec{p_i} \bullet \vec{r}) (\vec{p_j} \times \vec{r})\end{split}\]

where \(\epsilon\) and \(\sigma\) are the standard LJ parameters, \(r_c\) is the cutoff, \(q_i\) and \(q_j\) are the charges on the two particles, \(\vec{p_i}\) and \(\vec{p_j}\) are the dipole moment vectors of the two particles, r is their separation distance, and the vector r = Ri - Rj is the separation vector between the two particles. Note that Eqq and Fqq are simply Coulombic energy and force, Fij = -Fji as symmetric forces, and Tij != -Tji since the torques do not act symmetrically. The shifted-force formula for the Lennard-Jones potential is reported in (Stoddard). The original (non-shifted) formulas for the electrostatic potentials, forces and torques can be found in (Price). The shifted-force electrostatic potentials have been obtained by applying equation 5.13 of (Allen). The formulas for the corresponding forces and torques have been obtained by applying the ‘chain rule’ as in appendix C.3 of (Allen).

If one cutoff is specified in the pair_style command, it is used for both the LJ and Coulombic (q,p) terms. If two cutoffs are specified, they are used as cutoffs for the LJ and Coulombic (q,p) terms respectively. This pair style also supports an optional scale keyword as part of a pair_coeff statement, where the interactions can be scaled according to this factor. This scale factor is also made available for use with fix adapt.

Style lj/cut/dipole/long computes long-range point-dipole interactions as discussed in (Toukmaji). Dipole-dipole, dipole-charge, and charge-charge interactions are all supported, along with the standard 12/6 Lennard-Jones interactions, which are computed with a cutoff. A kspace_style must be defined to use this pair style. Currently, only kspace_style ewald/disp support long-range point-dipole interactions.

Style lj/long/dipole/long also computes point-dipole interactions as discussed in (Toukmaji). Long-range dipole-dipole, dipole-charge, and charge-charge interactions are all supported, along with the standard 12/6 Lennard-Jones interactions. LJ interactions can be cutoff or long-ranged.

For style lj/long/dipole/long, if flag_lj is set to long, no cutoff is used on the LJ 1/r^6 dispersion term. The long-range portion is calculated by using the kspace_style ewald_disp command. The specified LJ cutoff then determines which portion of the LJ interactions are computed directly by the pair potential versus which part is computed in reciprocal space via the Kspace style. If flag_lj is set to cut, the LJ interactions are simply cutoff, as with pair_style lj/cut. If flag_lj is set to off, LJ interactions are not computed at all.

If flag_coul is set to long, no cutoff is used on the Coulombic or dipole interactions. The long-range portion is calculated by using ewald_disp of the kspace_style command. If flag_coul is set to off, Coulombic and dipole interactions are not computed at all.

Atoms with dipole moments should be integrated using the fix nve/sphere update dipole or the fix nvt/sphere update dipole command to rotate the dipole moments. The omega option on the fix langevin command can be used to thermostat the rotational motion. The compute temp/sphere command can be used to monitor the temperature, since it includes rotational degrees of freedom. The atom_style hybrid dipole sphere command should be used since it defines the point dipoles and their rotational state. The magnitude and orientation of the dipole moment for each particle can be defined by the set command or in the “Atoms” section of the data file read in by the read_data command.

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)

  • cutoff1 (distance units)

  • cutoff2 (distance units)

The latter 2 coefficients are optional. If not specified, the global LJ and Coulombic cutoffs specified in the pair_style command are used. If only one cutoff is specified, it is used as the cutoff for both LJ and Coulombic interactions for this type pair. If both coefficients are specified, they are used as the LJ and Coulombic cutoffs for this type pair.

Styles with a gpu, intel, kk, omp, or opt suffix are functionally the same as the corresponding style without the suffix. They have been optimized to run faster, depending on your available hardware, as discussed on the Speed packages doc page. The accelerated styles take the same arguments and should produce the same results, except for round-off and precision issues.

These accelerated styles are part of the GPU, INTEL, KOKKOS, OPENMP and OPT packages, respectively. They are only enabled if LAMMPS was built with those packages. See the Build package page for more info.

You can specify the accelerated styles explicitly in your input script by including their suffix, or you can use the -suffix command-line switch when you invoke LAMMPS, or you can use the suffix command in your input script.

See the Speed packages page for more instructions on how to use the accelerated styles effectively.

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

For atom type pairs I,J and I != J, the epsilon and sigma coefficients and cutoff distances for this pair style can be mixed. The default mix value is geometric. See the “pair_modify” command for details.

For atom type pairs I,J and I != J, the A, sigma, d1, and d2 coefficients and cutoff distance for this pair style can be mixed. A is an energy value mixed like a LJ epsilon. D1 and d2 are distance values and are mixed like sigma. The default mix value is geometric. See the “pair_modify” command for details.

This pair style does not support the pair_modify shift option for the energy of the Lennard-Jones portion of the pair interaction; such energy goes to zero at the cutoff by construction.

The pair_modify table option is not relevant for this pair style.

This pair style does not support the pair_modify tail option for adding long-range tail corrections to energy and pressure.

This pair style writes its 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.

This pair style can only be used via the pair keyword of the run_style respa command. It does not support the inner, middle, outer keywords.

Restrictions

The lj/cut/dipole/cut, lj/cut/dipole/long, lj/long/dipole/long, and lj/sf/dipole/sf* styles are part of the DIPOLE package. They are only enabled if LAMMPS was built with that package. See the Build package page for more info.

Using dipole pair styles with electron units is not currently supported.

Default

none

(Allen) Allen and Tildesley, Computer Simulation of Liquids, Clarendon Press, Oxford, 1987.

(Toukmaji) Toukmaji, Sagui, Board, and Darden, J Chem Phys, 113, 10913 (2000).

(Stoddard) Stoddard and Ford, Phys Rev A, 8, 1504 (1973).

(Price) Price, Stone and Alderton, Mol Phys, 52, 987 (1984).