2010-10-01 04:38:58 +08:00
|
|
|
<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>fix neb command
|
|
|
|
|
</H3>
|
|
|
|
|
<P><B>Syntax:</B>
|
|
|
|
|
</P>
|
2010-10-01 05:22:06 +08:00
|
|
|
<PRE>fix ID group-ID neb Kspring
|
2010-10-01 04:38:58 +08:00
|
|
|
</PRE>
|
|
|
|
|
<UL><LI>ID, group-ID are documented in <A HREF = "fix.html">fix</A> command
|
|
|
|
|
<LI>neb = style name of this fix command
|
2010-10-01 05:22:06 +08:00
|
|
|
<LI>Kspring = inter-replica spring constant (force/distance units)
|
2010-10-01 04:38:58 +08:00
|
|
|
</UL>
|
|
|
|
|
<P><B>Examples:</B>
|
|
|
|
|
</P>
|
|
|
|
|
<PRE>fix 1 active neb 10.0
|
|
|
|
|
</PRE>
|
|
|
|
|
<P><B>Description:</B>
|
|
|
|
|
</P>
|
2010-10-01 05:22:06 +08:00
|
|
|
<P>Add inter-replica forces to atoms in the group for a multi-replica
|
2010-10-01 04:38:58 +08:00
|
|
|
simulation run via the <A HREF = "neb.html">neb</A> command to perform a nudged
|
|
|
|
|
elastic band (NEB) calculation for transition state finding. Hi-level
|
|
|
|
|
explanations of NEB are given with the <A HREF = "neb.html">neb</A> command and in
|
2011-08-26 01:01:01 +08:00
|
|
|
<A HREF = "Section_howto.html#howto_5">this section</A> of the manual. The fix neb
|
2010-10-01 05:22:06 +08:00
|
|
|
command must be used with the "neb" command to define how
|
|
|
|
|
inter-replica forces are computed.
|
2010-10-01 04:38:58 +08:00
|
|
|
</P>
|
2010-10-01 05:22:06 +08:00
|
|
|
<P>Only the N atoms in the fix group experience inter-replica forces.
|
|
|
|
|
Atoms in the two end-point replicas do not experience these forces,
|
|
|
|
|
but those in intermediate replicas do. During the initial stage of
|
|
|
|
|
NEB, the 3N-length vector of interatomic forces Fi = -Grad(V) acting
|
|
|
|
|
on the atoms of each intermediate replica I is altered, as described
|
|
|
|
|
in the <A HREF = "#Henkelman1">(Henkelman1)</A> paper, to become:
|
2010-10-01 04:38:58 +08:00
|
|
|
</P>
|
|
|
|
|
<PRE>Fi = -Grad(V) + (Grad(V) dot That) That + Kspring (|Ri+i - Ri| - |Ri - Ri-1|) That
|
|
|
|
|
</PRE>
|
|
|
|
|
<P>Ri are the atomic coordinates of replica I; Ri-1 and Ri+1 are the
|
|
|
|
|
coordinates of its neighbor replicas. That (t with a hat over it) is
|
|
|
|
|
the unit "tangent" vector for replica I which is a function of Ri,
|
|
|
|
|
Ri-1, Ri+1, and the potential energy of the 3 replicas; it points
|
2010-10-01 05:22:06 +08:00
|
|
|
roughly in the direction of (Ri+i - Ri-1); see the
|
|
|
|
|
<A HREF = "#Henkelman1">(Henkelman1)</A> paper for details.
|
2010-10-01 04:38:58 +08:00
|
|
|
</P>
|
2010-10-01 05:22:06 +08:00
|
|
|
<P>The first two terms in the above equation are the component of the
|
|
|
|
|
interatomic forces perpendicular to the tangent vector. The last term
|
|
|
|
|
is a spring force between replica I and its neighbors, parallel to the
|
|
|
|
|
tangent vector direction with the specified spring constant <I>Kspring</I>.
|
2010-10-01 04:38:58 +08:00
|
|
|
</P>
|
|
|
|
|
<P>The effect of the first two terms is to push the atoms of each replica
|
|
|
|
|
toward the minimum energy path (MEP) of conformational states that
|
|
|
|
|
transition over the energy barrier. The MEP for an energy barrier is
|
|
|
|
|
defined as a sequence of 3N-dimensional states which cross the barrier
|
|
|
|
|
at its saddle point, each of which has a potential energy gradient
|
|
|
|
|
parallel to the MEP itself.
|
|
|
|
|
</P>
|
|
|
|
|
<P>The effect of the last term is to push each replica away from its two
|
|
|
|
|
neighbors in a direction along the MEP, so that the final set of
|
|
|
|
|
states are equidistant from each other.
|
|
|
|
|
</P>
|
2010-10-01 05:22:06 +08:00
|
|
|
<P>During the second stage of NEB, the forces on the N atoms in the
|
|
|
|
|
replica nearest the top of the energy barrier are altered so that it
|
|
|
|
|
climbs to the top of the barrier and finds the saddle point. The
|
|
|
|
|
forces on atoms in this replica are described in the
|
|
|
|
|
<A HREF = "#Henkelman2">(Henkelman2)</A> paper, and become:
|
2010-10-01 04:38:58 +08:00
|
|
|
</P>
|
2010-10-01 05:22:06 +08:00
|
|
|
<PRE>Fi = -Grad(V) + 2 (Grad(V) dot That) That
|
2010-10-01 04:38:58 +08:00
|
|
|
</PRE>
|
|
|
|
|
<P>The inter-replica forces for the other replicas are unchanged from the
|
2010-10-01 05:22:06 +08:00
|
|
|
first equation.
|
2010-10-01 04:38:58 +08:00
|
|
|
</P>
|
|
|
|
|
<P><B>Restart, fix_modify, output, run start/stop, minimize info:</B>
|
|
|
|
|
</P>
|
|
|
|
|
<P>No information about this fix is written to <A HREF = "restart.html">binary restart
|
|
|
|
|
files</A>. None of the <A HREF = "fix_modify.html">fix_modify</A> options
|
|
|
|
|
are relevant to this fix. No global or per-atom quantities are stored
|
2011-08-26 01:01:01 +08:00
|
|
|
by this fix for access by various <A HREF = "Section_howto.html#howto_15">output
|
|
|
|
|
commands</A>. No parameter of this fix can
|
|
|
|
|
be used with the <I>start/stop</I> keywords of the <A HREF = "run.html">run</A> command.
|
2010-10-01 05:22:06 +08:00
|
|
|
</P>
|
|
|
|
|
<P>The forces due to this fix are imposed during an energy minimization,
|
|
|
|
|
as invoked by the <A HREF = "minimize.html">minimize</A> command via the
|
|
|
|
|
<A HREF = "neb.html">neb</A> command.
|
2010-10-01 04:38:58 +08:00
|
|
|
</P>
|
|
|
|
|
<P><B>Restrictions:</B>
|
|
|
|
|
</P>
|
2011-08-27 02:53:00 +08:00
|
|
|
<P>This command can only be used if LAMMPS was built with the REPLICA
|
2011-08-26 00:46:23 +08:00
|
|
|
package. See the <A HREF = "Section_start.html#start_3">Making LAMMPS</A> section
|
|
|
|
|
for more info on packages.
|
2010-10-01 04:38:58 +08:00
|
|
|
</P>
|
|
|
|
|
<P><B>Related commands:</B>
|
|
|
|
|
</P>
|
|
|
|
|
<P><A HREF = "neb.html">neb</A>
|
|
|
|
|
</P>
|
|
|
|
|
<P><B>Default:</B> none
|
|
|
|
|
</P>
|
2010-10-01 05:22:06 +08:00
|
|
|
<A NAME = "Henkelman"></A>
|
2010-10-01 04:38:58 +08:00
|
|
|
|
2010-10-01 05:22:06 +08:00
|
|
|
<P><B>(Henkelman1)</B> Henkelman and Jonsson, J Chem Phys, 113, 9978-9985 (2000).
|
2010-10-01 04:38:58 +08:00
|
|
|
</P>
|
2010-10-01 05:22:06 +08:00
|
|
|
<A NAME = "Henkelman"></A>
|
2010-10-01 04:38:58 +08:00
|
|
|
|
2010-10-01 05:22:06 +08:00
|
|
|
<P><B>(Henkelman2)</B> Henkelman, Uberuaga, Jonsson, J Chem Phys, 113,
|
|
|
|
|
9901-9904 (2000).
|
2010-10-01 04:38:58 +08:00
|
|
|
</P>
|
|
|
|
|
</HTML>
|
