<p>Perform constant NVE updates of position and velocity for atoms in the
group each timestep. A limit is imposed on the maximum distance an
atom can move in one timestep. This is useful when starting a
simulation with a configuration containing highly overlapped atoms.
Normally this would generate huge forces which would blow atoms out of
the simulation box, causing LAMMPS to stop with an error.</p>
<p>Using this fix can overcome that problem. Forces on atoms must still
be computable (which typically means 2 atoms must have a separation
distance > 0.0). But large velocities generated by large forces are
reset to a value that corresponds to a displacement of length <em>xmax</em>
in a single timestep. <em>Xmax</em> is specified in distance units; see the
<aclass="reference internal"href="units.html"><spanclass="doc">units</span></a> command for details. The value of <em>xmax</em> should be
consistent with the neighbor skin distance and the frequency of
neighbor list re-building, so that pairwise interactions are not
missed on successive timesteps as atoms move. See the
<aclass="reference internal"href="neighbor.html"><spanclass="doc">neighbor</span></a> and <aclass="reference internal"href="neigh_modify.html"><spanclass="doc">neigh_modify</span></a> commands
for details.</p>
<p>Note that if a velocity reset occurs the integrator will not conserve
energy. On steps where no velocity resets occur, this integrator is
exactly like the <aclass="reference internal"href="fix_nve.html"><spanclass="doc">fix nve</span></a> command. Since forces are
unaltered, pressures computed by thermodynamic output will still be
very large for overlapped configurations.</p>
<divclass="admonition note">
<pclass="first admonition-title">Note</p>
<pclass="last">You should not use <aclass="reference internal"href="fix_shake.html"><spanclass="doc">fix shake</span></a> in conjunction
with this fix. That is because fix shake applies contraint forces
based on the predicted postitions of atoms after the next timestep.
It has no way of knowing the timestep may change due to this fix,
which will cause the constraint forces to be invalid. A better
strategy is to turn off fix shake when performing initial dynamics
that need this fix, then turn fix shake on when doing normal dynamics
with a fixed-size timestep.</p>
</div>
<p><strong>Restart, fix_modify, output, run start/stop, minimize info:</strong></p>
<p>No information about this fix is written to <aclass="reference internal"href="restart.html"><spanclass="doc">binary restart files</span></a>. None of the <aclass="reference internal"href="fix_modify.html"><spanclass="doc">fix_modify</span></a> options
are relevant to this fix.</p>
<p>This fix computes a global scalar which can be accessed by various
<aclass="reference internal"href="Section_howto.html#howto-15"><spanclass="std std-ref">output commands</span></a>. The scalar is the
count of how many updates of atom’s velocity/position were limited by
the maximum distance criterion. This should be roughly the number of
atoms so affected, except that updates occur at both the beginning and
end of a timestep in a velocity Verlet timestepping algorithm. This
is a cumulative quantity for the current run, but is re-initialized to
zero each time a run is performed. The scalar value calculated by
this fix is “extensive”.</p>
<p>No parameter of this fix can be used with the <em>start/stop</em> keywords of
the <aclass="reference internal"href="run.html"><spanclass="doc">run</span></a> command. This fix is not invoked during <aclass="reference internal"href="minimize.html"><spanclass="doc">energy minimization</span></a>.</p>
Built with <ahref="http://sphinx-doc.org/">Sphinx</a> using a <ahref="https://github.com/snide/sphinx_rtd_theme">theme</a> provided by <ahref="https://readthedocs.org">Read the Docs</a>.