forked from lijiext/lammps
171 lines
6.1 KiB
HTML
171 lines
6.1 KiB
HTML
<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 edip command
|
|
</H3>
|
|
<P><B>Syntax:</B>
|
|
</P>
|
|
<PRE>pair_style edip
|
|
</PRE>
|
|
<PRE>pair_style edip/omp
|
|
</PRE>
|
|
<P><B>Examples:</B>
|
|
</P>
|
|
<P>pair_style edip
|
|
pair_coeff * * Si.edip Si
|
|
</P>
|
|
<P><B>Description:</B>
|
|
</P>
|
|
<P>The <I>edip</I> style computes a 3-body <A HREF = "#EDIP">EDIP</A> potential which is
|
|
popular for modeling silicon materials where it can have advantages
|
|
over other models such as the <A HREF = "pair_sw.html">Stillinger-Weber</A> or
|
|
<A HREF = "pair_tersoff.html">Tersoff</A> potentials. In EDIP, the energy E of a
|
|
system of atoms is
|
|
</P>
|
|
<CENTER><IMG SRC = "Eqs/pair_edip.jpg">
|
|
</CENTER>
|
|
<P>where phi2 is a two-body term and phi3 is a three-body term. The
|
|
summations in the formula are over all neighbors J and K of atom I
|
|
within a cutoff distance = a.
|
|
Both terms depend on the local environment of atom I through its
|
|
effective coordination number defined by Z, which is unity for a
|
|
cutoff distance < c and gently goes to 0 at distance = a.
|
|
</P>
|
|
<P>Only a single pair_coeff command is used with the <I>edip</I> style which
|
|
specifies a EDIP potential file with parameters for all
|
|
needed elements. These are mapped to LAMMPS atom types by specifying
|
|
N additional arguments after the filename in the pair_coeff command,
|
|
where N is the number of LAMMPS atom types:
|
|
</P>
|
|
<UL><LI>filename
|
|
<LI>N element names = mapping of EDIP elements to atom types
|
|
</UL>
|
|
<P>As an example, imagine a file Si.edip has EDIP values for Si.
|
|
</P>
|
|
<P>EDIP files in the <I>potentials</I> directory of the LAMMPS
|
|
distribution have a ".edip" suffix. Lines that are not blank or
|
|
comments (starting with #) define parameters for a triplet of
|
|
elements. The parameters in a single entry correspond to the two-body
|
|
and three-body coefficients in the formula above:
|
|
</P>
|
|
<UL><LI>element 1 (the center atom in a 3-body interaction)
|
|
<LI>element 2
|
|
<LI>element 3
|
|
<LI>A (energy units)
|
|
<LI>B (distance units)
|
|
<LI>cutoffA (distance units)
|
|
<LI>cutoffC (distance units)
|
|
<LI>alpha
|
|
<LI>beta
|
|
<LI>eta
|
|
<LI>gamma (distance units)
|
|
<LI>lambda (energy units)
|
|
<LI>mu
|
|
<LI>tho
|
|
<LI>sigma (distance units)
|
|
<LI>Q0
|
|
<LI>u1
|
|
<LI>u2
|
|
<LI>u3
|
|
<LI>u4
|
|
</UL>
|
|
<P>The A, B, beta, sigma parameters are used only for two-body interactions.
|
|
The eta, gamma, lambda, mu, Q0 and all u1 to u4 parameters are used only
|
|
for three-body interactions. The alpha and cutoffC parameters are used
|
|
for the coordination environment function only.
|
|
</P>
|
|
<P>The EDIP potential file must contain entries for all the
|
|
elements listed in the pair_coeff command. It can also contain
|
|
entries for additional elements not being used in a particular
|
|
simulation; LAMMPS ignores those entries.
|
|
</P>
|
|
<P>For a single-element simulation, only a single entry is required
|
|
(e.g. SiSiSi). For a two-element simulation, the file must contain 8
|
|
entries (for SiSiSi, SiSiC, SiCSi, SiCC, CSiSi, CSiC, CCSi, CCC), that
|
|
specify EDIP parameters for all permutations of the two elements
|
|
interacting in three-body configurations. Thus for 3 elements, 27
|
|
entries would be required, etc.
|
|
</P>
|
|
<P>At the moment, only a single element parametrization is
|
|
implemented. However, the author is not aware of other
|
|
multi-element EDIP parametrizations. If you know any and
|
|
you are interest in that, please contact the author of
|
|
the EDIP package.
|
|
</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>This pair style does not support the <A HREF = "pair_modify.html">pair_modify</A>
|
|
shift, table, and tail options.
|
|
</P>
|
|
<P>This pair style does not write its information to <A HREF = "restart.html">binary restart
|
|
files</A>, since it is stored in potential files. Thus, you
|
|
need to re-specify the pair_style and pair_coeff commands 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>
|
|
<HR>
|
|
|
|
<P><B>Restrictions:</B>
|
|
</P>
|
|
<P>This angle style can only be used if LAMMPS was built with the
|
|
USER-MISC package. See the <A HREF = "Section_start.html#start_3">Making LAMMPS</A>
|
|
section for more info on packages.
|
|
</P>
|
|
<P>This pair style requires the <A HREF = "newton.html">newton</A> setting to be "on"
|
|
for pair interactions.
|
|
</P>
|
|
<P>The EDIP potential files provided with LAMMPS (see the potentials directory)
|
|
are parameterized for metal <A HREF = "units.html">units</A>.
|
|
You can use the SW potential with any LAMMPS units, but you would need
|
|
to create your own EDIP potential file with coefficients listed in the
|
|
appropriate units if your simulation doesn't use "metal" units.
|
|
</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 = "EDIP"></A>
|
|
|
|
<P><B>(EDIP)</B> J. F. Justo et al., Phys. Rev. B 58, 2539 (1998).
|
|
</P>
|
|
</HTML>
|