lammps/doc/pair_dipole.html

217 lines
9.1 KiB
HTML
Raw Normal View History

<HTML>
<CENTER><A HREF = "http://lammps.sandia.gov">LAMMPS WWW Site</A> - <A HREF = "Manual.html">LAMMPS Documentation</A> - <A HREF = "Section_commands.html#comm">LAMMPS Commands</A>
</CENTER>
<HR>
<H3>pair_style dipole/cut command
</H3>
<H3>pair_style dipole/cut/omp command
</H3>
<H3>pair_style dipole/sf command
</H3>
<H3>pair_style dipole/sf/omp command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style dipole/cut cutoff (cutoff2)
</PRE>
<PRE>pair_style dipole/sf cutoff (cutoff2)
</PRE>
<UL><LI>cutoff = global cutoff LJ (and Coulombic if only 1 arg) (distance units)
<LI>cutoff2 = global cutoff for Coulombic (optional) (distance units)
</UL>
<P><B>Examples:</B>
</P>
<PRE>pair_style dipole/cut 10.0
pair_coeff * * 1.0 1.0
pair_coeff 2 3 1.0 1.0 2.5 4.0
</PRE>
<PRE>pair_style dipole/sf 9.0
pair_coeff * * 1.0 1.0
pair_coeff 2 3 1.0 1.0 2.5 4.0
</PRE>
<P><B>Description:</B>
</P>
<P>Style <I>dipole/cut</I> 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.
</P>
<CENTER><IMG SRC = "Eqs/pair_dipole.jpg">
</CENTER>
<P>where qi and qj are the charges on the two particles, pi and pj 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
<A HREF = "#Allen">(Allen)</A> and in <A HREF = "#Toukmaji">(Toukmaji)</A>.
</P>
<P>Style <I>dipole/sf</I> 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 (sligthly) 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 <A HREF = "#Allen">(Allen)</A>. 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:
</P>
<CENTER><IMG SRC = "Eqs/pair_dipole_sf.jpg">
</CENTER>
<CENTER><IMG SRC = "Eqs/pair_dipole_sf2.jpg">
</CENTER>
<P>where epsilon and sigma are the standard LJ parameters, r_c is the
cutoff, qi and qj are the charges on the two particles, pi and pj 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 <A HREF = "#Stoddard">(Stoddard)</A>. The
original (unshifted) formulas for the electrostatic potentials, forces
and torques can be found in <A HREF = "#Price">(Price)</A>. The shifted-force
electrostatic potentials have been obtained by applying equation 5.13
of <A HREF = "#Allen">(Allen)</A>. The formulas for the corresponding forces and
torques have been obtained by applying the 'chain rule' as in appendix
C.3 of <A HREF = "#Allen">(Allen)</A>.
</P>
<P>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.
</P>
<P>Atoms with dipole moments should be integrated using the <A HREF = "fix_nve_sphere.html">fix
nve/sphere update dipole</A> command to rotate the
dipole moments. The <A HREF = "compute_temp_sphere.html">compute temp/sphere</A>
command can be used to monitor the temperature, since it includes
rotational degrees of freedom. The <A HREF = "atom_style.html">atom_style
dipole</A> command should be used since it defines the
point dipoles and their rotational state. The magnitude of the dipole
moment for each type of particle can be defined by the
<A HREF = "dipole.html">dipole</A> command or in the "Dipoles" section of the data
file read in by the <A HREF = "read_data.html">read_data</A> command. Their initial
orientation can be defined by the <A HREF = "set.html">set dipole</A> command or in
the "Atoms" section of the data file.
</P>
<P>The following coefficients must be defined for each pair of atoms
types via the <A HREF = "pair_coeff.html">pair_coeff</A> command as in the examples
above, or in the data file or restart files read by the
<A HREF = "read_data.html">read_data</A> or <A HREF = "read_restart.html">read_restart</A>
commands, or by mixing as described below:
</P>
<UL><LI>epsilon (energy units)
<LI>sigma (distance units)
<LI>cutoff1 (distance units)
<LI>cutoff2 (distance units)
</UL>
<P>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.
</P>
<HR>
<P>Styles with a <I>cuda</I>, <I>gpu</I>, <I>omp</I>, or <I>opt</I> 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 in <A HREF = "Section_accelerate.html">Section_accelerate</A> of the
manual. The accelerated styles take the same arguments and should
produce the same results, except for round-off and precision issues.
</P>
<P>These accelerated styles are part of the USER-CUDA, GPU, USER-OMP and OPT
packages, respectively. They are only enabled if LAMMPS was built with
those packages. See the <A HREF = "Section_start.html#start_3">Making LAMMPS</A>
section for more info.
</P>
<P>You can specify the accelerated styles explicitly in your input script
by including their suffix, or you can use the <A HREF = "Section_start.html#start_7">-suffix command-line
switch</A> when you invoke LAMMPS, or you can
use the <A HREF = "suffix.html">suffix</A> command in your input script.
</P>
<P>See <A HREF = "Section_accelerate.html">Section_accelerate</A> of the manual for
more instructions on how to use the accelerated styles effectively.
</P>
<HR>
<P><B>Mixing, shift, table, tail correction, restart, rRESPA info</B>:
</P>
<P>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 <I>geometric</I>. See the "pair_modify" command for details.
</P>
<P>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
<I>geometric</I>. See the "pair_modify" command for details.
</P>
<P>This pair style does not support the <A HREF = "pair_modify.html">pair_modify</A>
shift option for the energy of the Lennard-Jones portion of the pair
interaction; such energy goes to zero at the cutoff by construction.
</P>
<P>The <A HREF = "pair_modify.html">pair_modify</A> table option is not relevant
for this pair style.
</P>
<P>This pair style does not support the <A HREF = "pair_modify.html">pair_modify</A>
tail option for adding long-range tail corrections to energy and
pressure.
</P>
<P>This pair style writes its information to <A HREF = "restart.html">binary restart
files</A>, so pair_style and pair_coeff commands do not need
to be specified in an input script that reads a restart file.
</P>
<P>This pair style can only be used via the <I>pair</I> keyword of the
<A HREF = "run_style.html">run_style respa</A> command. It does not support the
<I>inner</I>, <I>middle</I>, <I>outer</I> keywords.
</P>
<P><B>Restrictions:</B>
</P>
<P>The <I>dipole/cut</I> style is part of the DIPOLE package. It is only
enabled if LAMMPS was built with that package. See the <A HREF = "Section_start.html#start_3">Making
LAMMPS</A> section for more info.
</P>
<P>The <I>dipole/sf</I> style is part of the USER-MISC package. It is only
enabled if LAMMPS was built with that package. See the <A HREF = "Section_start.html#start_3">Making
LAMMPS</A> section for more info.
</P>
<P><B>Related commands:</B>
</P>
<P><A HREF = "pair_coeff.html">pair_coeff</A>
</P>
<P><B>Default:</B> none
</P>
<HR>
<A NAME = "Allen"></A>
<P><B>(Allen)</B> Allen and Tildesley, Computer Simulation of Liquids,
Clarendon Press, Oxford, 1987.
</P>
<A NAME = "Toukmaji"></A>
<P><B>(Toukmaji)</B> Toukmaji, Sagui, Board, and Darden, J Chem Phys, 113,
10913 (2000).
</P>
<A NAME = "Stoddard"></A>
<P><B>(Stoddard)</B> Stoddard and Ford, Phys Rev A, 8, 1504 (1973).
</P>
<A NAME = "Price"></A>
<P><B>(Price)</B> Price, Stone and Alderton, Mol Phys, 52, 987 (1984).
</P>
</HTML>