fix nve/dotc/langevin command
Syntax
fix ID group-ID nve/dotc/langevin Tstart Tstop damp seed keyword value
ID, group-ID are documented in fix command
nve/dotc/langevin = style name of this fix command
Tstart,Tstop = desired temperature at start/end of run (temperature units)
damp = damping parameter (time units)
seed = random number seed to use for white noise (positive integer)
keyword = angmom
angmom value = factor factor = do thermostat rotational degrees of freedom via the angular momentum and apply numeric scale factor as discussed below
Examples
fix 1 all nve/dotc/langevin 1.0 1.0 0.03 457145 angmom 10
fix 1 all nve/dotc/langevin 0.1 0.1 78.9375 457145 angmom 10
Description
Apply a rigid-body Langevin-type integrator of the kind “Langevin C” as described in (Davidchack) to a group of atoms, which models an interaction with an implicit background solvent. This command performs Brownian dynamics (BD) via a technique that splits the integration into a deterministic Hamiltonian part and the Ornstein-Uhlenbeck process for noise and damping. The quaternion degrees of freedom are updated though an evolution operator which performs a rotation in quaternion space, preserves the quaternion norm and is akin to (Miller).
In terms of syntax this command has been closely modelled on the fix langevin and its angmom option. But it combines the fix nve and the fix langevin in one single command. The main feature is improved stability over the standard integrator, permitting slightly larger timestep sizes.
Note
Unlike the fix langevin this command performs also time integration of the translational and quaternion degrees of freedom.
The total force on each atom will have the form:
Tstart and Tstop have to be constant values, i.e. they cannot be variables. If used together with the oxDNA force field for coarse-grained simulation of DNA please note that T = 0.1 in oxDNA units corresponds to T = 300 K.
The damp parameter is specified in time units and determines how
rapidly the temperature is relaxed. For example, a value of 0.03
means to relax the temperature in a timespan of (roughly) 0.03 time
units
The random # seed must be a positive integer. A Marsaglia random number generator is used. Each processor uses the input seed to generate its own unique seed and its own stream of random numbers. Thus the dynamics of the system will not be identical on two runs on different numbers of processors.
The keyword/value option has to be used in the following way:
This fix has to be used together with the angmom keyword. The particles are always considered to have a finite size. The keyword angmom enables thermostatting of the rotational degrees of freedom in addition to the usual translational degrees of freedom.
The scale factor after the angmom keyword gives the ratio of the rotational to the translational friction coefficient.
An example input file can be found in examples/PACKAGES/cgdna/examples/duplex2/. Further details of the implementation and stability of the integrators are contained in (Henrich). The preprint version of the article can be found here.
Restrictions
These pair styles can only be used if LAMMPS was built with the CG-DNA package and the MOLECULE and ASPHERE package. See the Build package page for more info.
Default
none
(Davidchack) R.L Davidchack, T.E. Ouldridge, M.V. Tretyakov. J. Chem. Phys. 142, 144114 (2015).
(Miller) T. F. Miller III, M. Eleftheriou, P. Pattnaik, A. Ndirango, G. J. Martyna, J. Chem. Phys., 116, 8649-8659 (2002).
(Dunweg) B. Dunweg, W. Paul, Int. J. Mod. Phys. C, 2, 817-27 (1991).
(Henrich) O. Henrich, Y. A. Gutierrez-Fosado, T. Curk, T. E. Ouldridge, Eur. Phys. J. E 41, 57 (2018).