1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
|
"LAMMPS WWW Site"_lws - "LAMMPS Documentation"_ld - "LAMMPS Commands"_lc :c
:link(lws,http://lammps.sandia.gov)
:link(ld,Manual.html)
:link(lc,Section_commands.html#comm)
:line
fix colvars command :h3
[Syntax:]
fix ID group-ID colvars configfile keyword values ... :pre
ID, group-ID are documented in "fix"_fix.html command :ulb,l
colvars = style name of this fix command :l
configfile = the configuration file for the colvars module :l
keyword = {input} or {output} or {seed} or {tstat} :l
{input} arg = colvars.state file name or prefix or NULL (default: NULL)
{output} arg = output filename prefix (default: out)
{seed} arg = seed for random number generator (default: 1966)
{tstat} arg = fix id of a thermostat or NULL (default: NULL) :pre
:ule
[Examples:]
fix colvars peptide peptide.colvars.inp seed 2122 input peptide.colvars.state output peptide
fix colvars all colvars.inp :pre
[Description:]
This fix interfaces LAMMPS to a "collective variables" or "colvars"
module library which allows to calculate potentials of mean force
(PMFs) for any set of colvars, using different sampling methods:
currently implemented are the Adaptive Biasing Force (ABF) method,
metadynamics, Steered Molecular Dynamics (SMD) and Umbrella Sampling
(US) via a flexible harmonic restraint bias. This documentation
describes only the colvars fix itself and LAMMPS specific parts of the
code. The documentation of the colvars implementation itself is
available as part of the "NAMD online
documentation"_http://www.ks.uiuc.edu/Research/namd/2.8/ug/node53.html
There are example scripts for using this package with LAMMPS in
examples/USER/colvars.
The implementation of the portable collective variable library is also
documented in "(Henin)"_#Henin
:line
The only mandatory argument to the fix is the filename to the colvars
input file that contains all input that is independent from the MD
program in which the colvars library has been integrated.
The {group-ID} entry is ignored. The collective variable module will
always apply to the entire system, i.e. use the group {all}.
The {input} keyword allows to specify a state file that would contain
the information required in order to continue a calculation, e.g.
from a restart. Setting it to NULL will start a new colvars run.
The {output} keyword allows to specify the output prefix. All output
files generated will use this prefix followed by the ".colvars." and
a word like "state" or "traj".
The {seed} keyword contains the seed for the random number generator
that will be used in the colvars module.
The {tstat} keyword can be either NULL or the label of a thermostating
fix that thermostats all atoms in the fix colvars group. This will be
used to provide the colvars module with the current thermostat target
temperature.
[Restart, fix_modify, output, run start/stop, minimize info:]
No information about this fix is written to "binary restart
files"_restart.html.
The "fix_modify"_fix_modify.html {energy} option is supported by this
fix to add the energy change from the biasing force added by the fix
to the system's potential energy as part of "thermodynamic
output"_thermo_style.html.
This fix computes a global scalar which can be accessed by various
"output commands"_Section_howto.html#howto_15. The scalar is the
cummulative energy change due to this fix. The scalar value
calculated by this fix is "extensive".
[Restrictions:]
This fix is part of the USER-COLVARS package. It is only enabled if
LAMMPS was built with that package. See the "Making
LAMMPS"_Section_start.html#start_3 section for more info.
There can only be one colvars fix active at a time. Since the
colvars module itself can handle an arbitrary number of collective
variables and always applies to the entire system, this is not
really a deficit in practice.
[Related commands:]
"fix smd"_fix_smd.html
[Default:]
The default options are input = NULL, output = out, seed = 1966, and
tstat = NULL.
:line
:link(Henin)
[(Henin)] Hénin, Fiorin, Chipot, Klein, J. Chem. Theory Comput., 6,
35-47 (2010)
|
