git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@1705 f3b2605a-c512-4ea7-a41b-209d697bcdaa

This commit is contained in:
sjplimp 2008-04-02 23:05:25 +00:00
parent 67f28d776b
commit b84693f8ed
4 changed files with 72 additions and 26 deletions

View File

@ -13,10 +13,10 @@
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style lj/smooth Rin cutoff
<PRE>pair_style lj/smooth Rin Rc
</PRE>
<UL><LI>Rin = global inner cutoff beyond which force smoothing will be applied (distance units)
<LI>cutoff = global cutoff for lj/smooth interactions (distance units)
<UL><LI>Rin = inner cutoff beyond which force smoothing will be applied (distance units)
<LI>Rc = outer cutoff for lj/smooth interactions (distance units)
</UL>
<P><B>Examples:</B>
</P>
@ -32,9 +32,12 @@ applied between the inner and outer cutoff.
<CENTER><IMG SRC = "Eqs/pair_lj_smooth.jpg">
</CENTER>
<P>The polynomial coefficients C1, C2, C3, C4 are computed by LAMMPS to
cause the force to vary smoothly from Rin to Rc. At Rin the force and
its 1st derivative will match the unsmoothed LJ formula. At Rc the
force and its 1st derivative will be 0.0.
cause the force to vary smoothly from the inner cutoff Rin to the
outer cutoff Rc.
</P>
<P>At the inner cutoff the force and its 1st derivative
will match the unsmoothed LJ formula. At the outer cutoff the force
and its 1st derivative will be 0.0. The inner cutoff cannot be 0.0.
</P>
<P>IMPORTANT NOTE: this force smoothing causes the energy to be
discontinuous both in its values and 1st derivative. This can lead to
@ -50,13 +53,11 @@ commands, or by mixing as described below:
</P>
<UL><LI>epsilon (energy units)
<LI>sigma (distance units)
<LI>Rin (distance units)
<LI>cutoff (distance units)
<LI>innner (distance units)
<LI>outer (distance units)
</UL>
<P>The last 2 coefficients are optional. If not specified, the global
Rin and cutoff are used. Rin cannot be 0.0. If Rin = cutoff, then no
force smoothing is performed for this type pair; the standard LJ
formula is used.
<P>The last 2 coefficients are optional inner and outer cutoffs. If not
specified, the global values for Rin and Rc are used.
</P>
<HR>

View File

@ -10,10 +10,10 @@ pair_style lj/smooth command :h3
[Syntax:]
pair_style lj/smooth Rin cutoff :pre
pair_style lj/smooth Rin Rc :pre
Rin = global inner cutoff beyond which force smoothing will be applied (distance units)
cutoff = global cutoff for lj/smooth interactions (distance units) :ul
Rin = inner cutoff beyond which force smoothing will be applied (distance units)
Rc = outer cutoff for lj/smooth interactions (distance units) :ul
[Examples:]
@ -29,9 +29,12 @@ applied between the inner and outer cutoff.
:c,image(Eqs/pair_lj_smooth.jpg)
The polynomial coefficients C1, C2, C3, C4 are computed by LAMMPS to
cause the force to vary smoothly from Rin to Rc. At Rin the force and
its 1st derivative will match the unsmoothed LJ formula. At Rc the
force and its 1st derivative will be 0.0.
cause the force to vary smoothly from the inner cutoff Rin to the
outer cutoff Rc.
At the inner cutoff the force and its 1st derivative
will match the unsmoothed LJ formula. At the outer cutoff the force
and its 1st derivative will be 0.0. The inner cutoff cannot be 0.0.
IMPORTANT NOTE: this force smoothing causes the energy to be
discontinuous both in its values and 1st derivative. This can lead to
@ -47,13 +50,11 @@ commands, or by mixing as described below:
epsilon (energy units)
sigma (distance units)
Rin (distance units)
cutoff (distance units) :ul
innner (distance units)
outer (distance units) :ul
The last 2 coefficients are optional. If not specified, the global
Rin and cutoff are used. Rin cannot be 0.0. If Rin = cutoff, then no
force smoothing is performed for this type pair; the standard LJ
formula is used.
The last 2 coefficients are optional inner and outer cutoffs. If not
specified, the global values for Rin and Rc are used.
:line

View File

@ -72,8 +72,30 @@ an optional coefficient. If not specified, the outer cutoff in the
table itself (see below) will be used to build an interpolation table
that extend to the largest tabulated distance. If specified, only
file values up to the cutoff are used to create the interpolation
table.
table. The format of this file is described below.
</P>
<HR>
<LI>Here are some guidelines for using the pair_style table command to
best effect:
</UL>
<LI>Vary the number of table points; you may need to use more than you think
to get good resolution.
<LI>Always use the <A HREF = "pair_write.html">pair_write</A> command to produce a plot
of what the final interpolated potential looks like. This can show up
interpolations "features" you may not like.
<LI>Start with the linear style; it's the style least likely to have problems.
<LI>Use as large an inner cutoff as possible. This avoids fitting splines
to very steep parts of the potential.
<LI>Make sure the potential tabulation you are feeding to LAMMPS is not
pathological.
<HR>
<P>The format of a tabulated file is as follows (without the
parenthesized comments):
</P>

View File

@ -69,7 +69,29 @@ an optional coefficient. If not specified, the outer cutoff in the
table itself (see below) will be used to build an interpolation table
that extend to the largest tabulated distance. If specified, only
file values up to the cutoff are used to create the interpolation
table.
table. The format of this file is described below.
:line
Here are some guidelines for using the pair_style table command to
best effect: :ule,l
Vary the number of table points; you may need to use more than you think
to get good resolution. :l
Always use the "pair_write"_pair_write.html command to produce a plot
of what the final interpolated potential looks like. This can show up
interpolations "features" you may not like. :l
Start with the linear style; it's the style least likely to have problems. :l
Use as large an inner cutoff as possible. This avoids fitting splines
to very steep parts of the potential. :l
Make sure the potential tabulation you are feeding to LAMMPS is not
pathological. :l
:line
The format of a tabulated file is as follows (without the
parenthesized comments):