Project ID: plumID:22.039
Source: urea_plumed.dat
Originally used with PLUMED version: 2.6.1
Stable: zipped raw stdout - zipped raw stderr - stderr
Master: zipped raw stdout - zipped raw stderr - stderr

Click on the labels of the actions for more information on what each action computes
tested onv2.9
tested onmaster
tested on master
# Restart simulation. This command will append results of new simulations to the same output files, e.g. COLVAR.
#RESTART

LOAD
Loads a library, possibly defining new actions. More details
FILE
file to be loaded
=InterFace.cpp
LOAD
Loads a library, possibly defining new actions. More details
FILE
file to be loaded
=PairEntropy.cpp
LOAD
Loads a library, possibly defining new actions. More details
FILE
file to be loaded
=ManyAngle.cpp
# Define groups for the CV
INCLUDE
Includes an external input file, similar to #include in C preprocessor. More details. Show included file
FILE
file to be included
=centers.dat
C:
GROUP
Define a group of atoms so that a particular list of atoms can be referenced with a single label in definitions of CVs or virtual atoms. More details
ATOMS
the numerical indexes for the set of atoms in the group
=1-2400:8 O:
GROUP
Define a group of atoms so that a particular list of atoms can be referenced with a single label in definitions of CVs or virtual atoms. More details
ATOMS
the numerical indexes for the set of atoms in the group
=2-2400:8 N1:
GROUP
Define a group of atoms so that a particular list of atoms can be referenced with a single label in definitions of CVs or virtual atoms. More details
ATOMS
the numerical indexes for the set of atoms in the group
=3-2400:8 N2:
GROUP
Define a group of atoms so that a particular list of atoms can be referenced with a single label in definitions of CVs or virtual atoms. More details
ATOMS
the numerical indexes for the set of atoms in the group
=6-2400:8 W:
GROUP
Define a group of atoms so that a particular list of atoms can be referenced with a single label in definitions of CVs or virtual atoms. More details
ATOMS
the numerical indexes for the set of atoms in the group
=2401-11655:3 # Define the CV
MANY_ANGLE
This action is not part of PLUMED and was included by using a LOAD command More details
... LABEL=ma1 CENTER=C START=Ncenter END=O RCUT=1.0 NOPBC MEAN MOMENT2 ... MANY_ANGLE
MANY_ANGLE
This action is not part of PLUMED and was included by using a LOAD command More details
... LABEL=ma2 CENTER=C START=N1 END=N2 RCUT=1.0 NOPBC MEAN MOMENT2 ... MANY_ANGLE
COORDINATIONNUMBER
Calculate the coordination numbers of atoms so that you can then calculate functions of the distribution of More details
...
LABEL
a label for the action so that its output can be referenced in the input to other actions
=cn
SPECIES
this keyword is used for colvars such as coordination number
=Ncenter
SWITCH
the switching function that it used in the construction of the contact matrix
={RATIONAL R_0=0.6 D_MAX=0.8}
MORE_THAN1
calculate the number of variables that are more than a certain target value
={RATIONAL R_0=8.0 D_MAX=9.0}
MORE_THAN2
calculate the number of variables that are more than a certain target value
={RATIONAL R_0=11.0 D_MAX=12.0}
MOMENTS
the list of moments that you would like to calculate
=2
LOWMEM
this flag does nothing and is present only to ensure back-compatibility
... COORDINATIONNUMBER
INTERFACE
This action is not part of PLUMED and was included by using a LOAD command More details
... LABEL=int CENTER=C Water=W R0CUT=0.8 RWCUT=0.4 R_0=0.5 R_W=0.3 C_0=5 ... INTERFACE
Q6
Calculate sixth order Steinhardt parameters. More details
...
LABEL
a label for the action so that its output can be referenced in the input to other actions
=q6
SPECIES
this keyword is used for colvars such as coordination number
=C
R_0
The r_0 parameter of the switching function
=0.60
MEAN
calculate the mean of all the quantities
LOWMEM
this flag does nothing and is present only to ensure back-compatibility
... Q6
PAIRENTROPY
Calculate the KL Entropy from the radial distribution function More details
...
LABEL
a label for the action so that its output can be referenced in the input to other actions
=s ATOMS=C
MAXR
the maximum distance to use for the rdf
=0.6 SIGMA=0.05 ... PAIRENTROPY
# Construct a bias potential using metadynamics
l1_1:
COMBINE
Calculate a polynomial combination of a set of other variables. More details
ARG
the input to this function
=ma1.mean,ma2.mean,ma1.moment2,ma2.moment2,cn.morethan-1,cn.morethan-2,cn.moment-2,int,q6.mean,s
COEFFICIENTS
the coefficients of the arguments in your function
=1.192779302597,-4.754418849945,0.035175766796,-0.366665005684,0.002653433941,0.011365669779,-0.020645860583,-0.000668735709,-0.072709567845,0.142931178212
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=NO l1_2:
COMBINE
Calculate a polynomial combination of a set of other variables. More details
ARG
the input to this function
=ma1.mean,ma2.mean,ma1.moment2,ma2.moment2,cn.morethan-1,cn.morethan-2,cn.moment-2,int,q6.mean,s
COEFFICIENTS
the coefficients of the arguments in your function
=-4.171144485474,-0.252191036940,0.253239363432,-0.522297024727,-0.006509369239,0.100738972425,-0.119564652443,0.000595316174,1.388457655907,-0.096874527633
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=NO l1r_1:
MATHEVAL
An alias to the ef CUSTOM function. More details
ARG
the input to this function
=l1_1
FUNC
the function you wish to evaluate
=x+2.238606452942
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=NO l1r_2:
MATHEVAL
An alias to the ef CUSTOM function. More details
ARG
the input to this function
=l1_2
FUNC
the function you wish to evaluate
=x+2.579712629318
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=NO
# load metaD
METAD
Used to performed metadynamics on one or more collective variables. More details
...
LABEL
a label for the action so that its output can be referenced in the input to other actions
=metad
ARG
the input for this action is the scalar output from one or more other actions
=l1r_1,l1r_2
SIGMA
the widths of the Gaussian hills
=0.2,0.2
HEIGHT
the heights of the Gaussian hills
=10
BIASFACTOR
use well tempered metadynamics and use this bias factor
=150
TEMP
the system temperature - this is only needed if you are doing well-tempered metadynamics
=300.0
PACE
the frequency for hill addition
=2000
GRID_MIN
the lower bounds for the grid
=-5,-5
GRID_MAX
the upper bounds for the grid
=5,20
GRID_BIN
the number of bins for the grid
=1000,1000
CALC_RCT
calculate the c(t) reweighting factor and use that to obtain the normalized bias [rbias=bias-rct]
... METAD
PRINT
Print quantities to a file. More details
ARG
the input for this action is the scalar output from one or more other actions
=ma1.mean,ma2.mean,ma1.moment2,ma2.moment2,cn.morethan-1,cn.morethan-2,cn.moment-2,int,q6.mean,s,l1r_1,l1r_2,metad.bias,metad.rbias
STRIDE
the frequency with which the quantities of interest should be output
=500
FILE
the name of the file on which to output these quantities
=COLVAR