+
+ +
+

fix nvk command

+
+

Syntax

+
fix ID group-ID nvk
+
+
+
    +
  • ID, group-ID are documented in fix command
  • +
  • nvk = style name of this fix command
  • +
+
+
+

Examples

+
fix 1 all nvk
+
+
+
+
+

Description

+

Perform constant kinetic energy integration using the Gaussian +thermostat to update position and velocity for atoms in the group each +timestep. V is volume; K is kinetic energy. This creates a system +trajectory consistent with the isokinetic ensemble.

+

The equations of motion used are those of Minary et al in +(Minary), a variant of those initially given by Zhang in +(Zhang).

+

The kinetic energy will be held constant at its value given when fix +nvk is initiated. If a different kinetic energy is desired, the +velocity command should be used to change the kinetic +energy prior to this fix.

+
+

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. +This fix is not invoked during energy minimization.

+
+
+

Restrictions

+

The Gaussian thermostat only works when it is applied to all atoms in +the simulation box. Therefore, the group must be set to all.

+

This fix has not yet been implemented to work with the RESPA integrator.

+

This fix is part of the USER-MISC package. It is only enabled if LAMMPS +was built with that package. See the Making LAMMPS section for more info.

+

Related commands: none

+

Default: none

+
+

(Minary) Minary, Martyna, and Tuckerman, J Chem Phys, 18, 2510 (2003).

+

(Zhang) Zhang, J Chem Phys, 106, 6102 (1997).

+
+
+ + +
+