forked from lijiext/lammps
73 lines
2.2 KiB
Plaintext
73 lines
2.2 KiB
Plaintext
"LAMMPS WWW Site"_lws - "LAMMPS Documentation"_ld - "LAMMPS Commands"_lc :c
|
|
|
|
:link(lws,http://lammps.sandia.gov)
|
|
:link(ld,Manual.html)
|
|
:link(lc,Section_commands.html#comm)
|
|
|
|
:line
|
|
|
|
improper_style ring command :h3
|
|
|
|
[Syntax:]
|
|
|
|
improper_style ring :pre
|
|
|
|
[Examples:]
|
|
|
|
improper_style ring
|
|
improper_coeff 1 8000 70.5 :pre
|
|
|
|
[Description:]
|
|
|
|
The {ring} improper style uses the potential
|
|
|
|
:c,image(Eqs/improper_ring.jpg)
|
|
|
|
where K is a prefactor, theta is the angle formed by the atoms
|
|
specified by (i,j,k,l) indices and theta0 its equilibrium value.
|
|
|
|
If the 4 atoms in an improper quadruplet (listed in the data file read
|
|
by the "read_data"_read_data.html command) are ordered i,j,k,l then
|
|
theta_{ijl} is the angle between atoms i,j and l, theta_{ijk} is the
|
|
angle between atoms i,j and k, theta_{kjl} is the angle between atoms
|
|
j,k, and l.
|
|
|
|
The "ring" improper style implements the improper potential introduced
|
|
by Destree et al., in Equation (9) of "(Destree)"_#Destree. This
|
|
potential does not affect small amplitude vibrations but is used in an
|
|
ad-hoc way to prevent the onset of accidentially large amplitude
|
|
fluctuations leading to the occurrence of a planar conformation of the
|
|
three bonds i-j, j-k and j-l, an intermediate conformation toward the
|
|
chiral inversion of a methine carbon. In the "Impropers" section of
|
|
data file four atoms: i, j, k and l are specified with i,j and l lying
|
|
on the backbone of the chain and k specifying the chirality of j.
|
|
|
|
The following coefficients must be defined for each improper type via
|
|
the "improper_coeff"_improper_coeff.html command as in the example
|
|
above, or in the data file or restart files read by the
|
|
"read_data"_read_data.html or "read_restart"_read_restart.html
|
|
commands:
|
|
|
|
K (energy/radian^2)
|
|
theta0 (degrees) :ul
|
|
|
|
theta0 is specified in degrees, but LAMMPS converts it to radians
|
|
internally; hence the units of K are in energy/radian^2.
|
|
|
|
:line
|
|
|
|
[Restrictions:]
|
|
|
|
This improper style can only be used if LAMMPS was built with the
|
|
USER-MISC package. See the "Making LAMMPS"_Section_start.html#start_3
|
|
section for more info on packages.
|
|
|
|
[Related commands:]
|
|
|
|
"improper_coeff"_improper_coeff.html
|
|
|
|
:link(Destree)
|
|
[(Destree)] M. Destree, F. Laupretre, A. Lyulin, and J.-P. Ryckaert,
|
|
J Chem Phys, 112, 9632 (2000).
|
|
|