Homec4science
Diffusion lammps (lammps-icms) Browse

lammps/doc/pair_charmm.html
aba05c2cf991lammps-icms

Last Change
Raw File

pair_charmm.html
View Options

<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 lj/charmm/coul/charmm command
</H3>
<H3>pair_style lj/charmm/coul/charmm/implicit command
</H3>
<H3>pair_style lj/charmm/coul/long command
</H3>
<H3>pair_style lj/charmm/coul/long/opt command
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style style args
</PRE>
<UL><LI>style = <I>lj/charmm/coul/charmm</I> or <I>lj/charmm/coul/charmm/implicit</I> or <I>lj/charmm/coul/long</I> or <I>lj/charmm/coul/long/opt</I>
<LI>args = list of arguments for a particular style
</UL>
<PRE> <I>lj/charmm/coul/charmm</I> args = inner outer (inner2) (outer2)
inner, outer = global switching cutoffs for Lennard Jones (and Coulombic if only 2 args)
inner2, outer2 = global switching cutoffs for Coulombic (optional)
<I>lj/charmm/coul/charmm/implicit</I> args = inner outer (inner2) (outer2)
inner, outer = global switching cutoffs for LJ (and Coulombic if only 2 args)
inner2, outer2 = global switching cutoffs for Coulombic (optional)
<I>lj/charmm/coul/long</I> args = inner outer (cutoff)
inner, outer = global switching cutoffs for LJ (and Coulombic if only 2 args)
cutoff = global cutoff for Coulombic (optional, outer is Coulombic cutoff if only 2 args)
</PRE>
<P><B>Examples:</B>
</P>
<PRE>pair_style lj/charmm/coul/charmm 8.0 10.0
pair_style lj/charmm/coul/charmm 8.0 10.0 7.0 9.0
pair_coeff * * 100.0 2.0
pair_coeff 1 1 100.0 2.0 150.0 3.5
</PRE>
<PRE>pair_style lj/charmm/coul/charmm/implicit 8.0 10.0
pair_style lj/charmm/coul/charmm/implicit 8.0 10.0 7.0 9.0
pair_coeff * * 100.0 2.0
pair_coeff 1 1 100.0 2.0 150.0 3.5
</PRE>
<PRE>pair_style lj/charmm/coul/long 8.0 10.0
pair_style lj/charmm/coul/long/opt 8.0 10.0
pair_style lj/charmm/coul/long 8.0 10.0 9.0
pair_coeff * * 100.0 2.0
pair_coeff 1 1 100.0 2.0 150.0 3.5
</PRE>
<P><B>Description:</B>
</P>
<P>The <I>lj/charmm</I> styles compute LJ and Coulombic interactions with an
additional switching function S(r) that ramps the energy and force
smoothly to zero between an inner and outer cuoff. It is a widely
used potential in the <A HREF = "http://www.scripps.edu/brooks">CHARMM</A> MD code.
See <A HREF = "#MacKerell">(MacKerell)</A> for a description of the CHARMM force
field.
</P>
<CENTER><IMG SRC = "Eqs/pair_charmm.jpg">
</CENTER>
<P>Both the LJ and Coulombic terms require an inner and outer cutoff.
They can be the same for both formulas or different depending on
whether 2 or 4 arguments are used in the pair_style command. In each
case, the inner cutoff distance must be less than the outer cutoff.
It it typical to make the difference between the 2 cutoffs about 1.0
Angstrom.
</P>
<P>Style <I>lj/charmm/coul/charmm/implicit</I> computes the same formulas as
style <I>lj/charmm/coul/charmm</I> except that an additional 1/r term is
included in the Coulombic formula. The Coulombic energy thus varies
as 1/r^2. This is effectively a distance-dependent dielectric term
which is a simple model for an implicit solvent with additional
screening. It is designed for use in a simulation of an unsolvated
biomolecule (no explicit water molecules).
</P>
<P>Style <I>lj/charmm/coul/long</I> computes the same formulas as style
<I>lj/charmm/coul/charmm</I> except that an additional damping factor is
applied to the Coulombic term, as in the discussion for pair style
<I>lj/cut/coul/long</I>. Only one Coulombic cutoff is specified for
<I>lj/charmm/coul/long</I>; if only 2 arguments are used in the pair_style
command, then the outer LJ cutoff is used as the single Coulombic
cutoff.
</P>
<P>Style <I>lj/charmm/coul/long/opt</I> is an optimized version of style
<I>lj/charmm/coul/long</I> that should give identical answers. Depending
on system size and the processor you are running on, it may be 5-25%
faster (for the pairwise portion of the run time).
</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:
</P>
<UL><LI>epsilon (energy units)
<LI>sigma (distance units)
<LI>epsilon_14 (energy units)
<LI>sigma_14 (distance units)
</UL>
<P>Note that sigma is defined in the LJ formula as the zero-crossing
distance for the potential, not as the energy minimum at 2^(1/6)
sigma.
</P>
<P>The latter 2 coefficients are optional. If they are specified, they
are used in the LJ formula between 2 atoms of these types which are
also first and fourth atoms in any dihedral. No cutoffs are specified
because this CHARMM force field does not allow varying cutoffs for
individual atom pairs; all pairs use the global cutoff(s) specified in
the pair_style command.
</P>
<P>If the pair_coeff command is not used to define coefficients for a
particular I != J type pair, the mixing rule for epsilon and sigma for
all CHARMM potentials is to use the <I>arithmetic</I> formulas documented
by the <A HREF = "pair_modify.html">pair_modify</A> command. The <A HREF = "pair_modify.html">pair_modify
mix</A> setting is thus ignored for CHARMM potentials.
</P>
<P><B>Restrictions:</B>
</P>
<P>The <I>lj/charmm/coul/charmm</I> and <I>lj/charmm/coul/charmm/implicit</I>
styles are part of the "molecule" package. The <I>lj/charmm/coul/long</I>
style is part of the "kspace" package. The <I>lj/charmm/coul/long/opt</I>
style is part of the "opt" package and also requires the "kspace"
package. They are only enabled if LAMMPS was built with those
package(s). See the <A HREF = "Section_start.html#2_3">Making LAMMPS</A> section
for more info.
</P>
<P>On some 64-bit machines, compiling with -O3 appears to break the
Coulombic tabling option used by the <I>lj/charmm/coul/long</I> style. See
the "Additional build tips" section of the Making LAMMPS documentation
pages for workarounds on this issue.
</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 = "MacKerell"></A>
<P><B>(MacKerell)</B> MacKerell, Bashford, Bellott, Dunbrack, Evanseck, Field,
Fischer, Gao, Guo, Ha, et al, J Phys Chem, 102, 3586 (1998).
</P>
</HTML>
c4science · Help