fix neb/spin command


fix ID group-ID neb/spin Kspring
  • ID, group-ID are documented in fix command

  • neb/spin = style name of this fix command

Kspring = spring constant for parallel nudging force
(force/distance units or force units, see parallel keyword)


fix 1 active neb/spin 1.0


Add nudging forces to spins in the group for a multi-replica simulation run via the neb/spin command to perform a geodesic nudged elastic band (GNEB) calculation for finding the transition state. Hi-level explanations of GNEB are given with the neb/spin command and on the Howto replica doc page. The fix neb/spin command must be used with the “neb/spin” command and defines how inter-replica nudging forces are computed. A GNEB calculation is divided in two stages. In the first stage n replicas are relaxed toward a MEP until convergence. In the second stage, the climbing image scheme is enabled, so that the replica having the highest energy relaxes toward the saddle point (i.e. the point of highest energy along the MEP), and a second relaxation is performed.

The nudging forces are calculated as explained in (BessarabB)). See this reference for more explanation about their expression.

Restart, fix_modify, output, run start/stop, minimize info:

No information about this fix is written to binary restart files. None of the fix_modify options are relevant to this fix. No global or per-atom quantities are stored by this fix for access by various output commands. No parameter of this fix can be used with the start/stop keywords of the run command.

The forces due to this fix are imposed during an energy minimization, as invoked by the minimize command via the neb/spin command.


This command can only be used if LAMMPS was built with the SPIN package. See the Build package doc page for more info.



(BessarabB) Bessarab, Uzdin, Jonsson, Comp Phys Comm, 196, 335-347 (2015).