Project ID: plumID:21.002
Source: Crystallization/IceIh/96molecules/Multithermal/plumed.dat
Originally used with PLUMED version: 2.8
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
#SETTINGS NREPLICAS=2
# vim:ft=plumed
Enables syntax highlighting for PLUMED files in vim. See here for more details.

#RESTART
vol:
VOLUME
Calculate the volume of the simulation box. More details

ENVIRONMENTSIMILARITY
Measure how similar the environment around atoms is to that found in some reference crystal structure. This action is a shortcut and it has hidden defaults. More details
...
SPECIES
this keyword is used for colvars such as coordination number
=1-288:3
SIGMA
the width to use for the gaussian kernels
=0.052
CRYSTAL_STRUCTURE
Targeted crystal structure
=CUSTOM
LABEL
a label for the action so that its output can be referenced in the input to other actions
=refcv
REFERENCE_1
PDB files with relative distances from central atom
=env1h.pdb
REFERENCE_2
PDB files with relative distances from central atom
=env2h.pdb
REFERENCE_3
PDB files with relative distances from central atom
=env3h.pdb
REFERENCE_4
PDB files with relative distances from central atom
=env4h.pdb
MORE_THAN
calculate the number of variables that are more than a certain target value
={RATIONAL R_0=0.5 NN=8 MM=16}
MEAN
calculate the mean of all the quantities
... ENVIRONMENTSIMILARITY

# Construct a bias potential
energy:
ENERGY
Calculate the total potential energy of the simulation box. More details

# Construct bias potential with OPS
ecv:
ECV_MULTITHERMAL_MULTIBARIC
Expand a simulation to sample multiple temperatures and pressures. More details
...
ARG
the labels of the potential energy and of the volume of the system
=energy,vol
TEMP
temperature
=330
TEMP_MIN
the minimum of the temperature range
=300
TEMP_MAX
the maximum of the temperature range
=350
PRESSURE
pressure
=0.06022140857
PRESSURE_MIN
the minimum of the pressure range
=0.06022140857
PRESSURE_MAX
the maximum of the pressure range
=0.06022140857 ...
ecv2:
ECV_UMBRELLAS_LINE
Target a multiumbrella ensemble, by combining systems each with a parabolic bias potential at a different location. This action has hidden defaults. More details
ARG
the input for this action is the scalar output from one or more other actions
=refcv.morethan
CV_MIN
the minimum of the CV range to be explored
=0.0
CV_MAX
the maximum of the CV range to be explored
=96.0
SIGMA
sigma of the umbrella Gaussians
=1.0
BARRIER
a guess of the free energy barrier to be overcome (better to stay higher than lower)
=150
TEMP
temperature
=330

opes:
OPES_EXPANDED
On-the-fly probability enhanced sampling with expanded ensembles for the target distribution. This action has hidden defaults. More details
ARG
the label of the ECVs that define the expansion
=ecv.*,ecv2
FILE
a file with the estimate of the relative Delta F for each component of the target and of the global c(t)
=DeltaF.data
PACE
how often the bias is updated
=100
WALKERS_MPI
switch on MPI version of multiple walkers

ENVIRONMENTSIMILARITY
Measure how similar the environment around atoms is to that found in some reference crystal structure. This action is a shortcut and it has hidden defaults. More details
...
SPECIES
this keyword is used for colvars such as coordination number
=1-288:3
SIGMA
the width to use for the gaussian kernels
=0.052
CRYSTAL_STRUCTURE
Targeted crystal structure
=CUSTOM
LABEL
a label for the action so that its output can be referenced in the input to other actions
=refcv2
REFERENCE_1
PDB files with relative distances from central atom
=env1c.pdb
REFERENCE_2
PDB files with relative distances from central atom
=env2c.pdb
MORE_THAN
calculate the number of variables that are more than a certain target value
={RATIONAL R_0=0.5 NN=10 MM=20}
MEAN
calculate the mean of all the quantities
... ENVIRONMENTSIMILARITY

diff:
MATHEVAL
An alias to the CUSTOM function that can also be used to calaculate combinations of variables using a custom expression. More details
ARG
the input to this function
=refcv2.mean,refcv.mean
FUNC
the function you wish to evaluate
=((x-0.230)/(0.585-0.230)-(y-0.254)/(0.812-0.254
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=NO
UPPER_WALLS
Defines a wall for the value of one or more collective variables, More details
ARG
the arguments on which the bias is acting
=diff
AT
the positions of the wall
=0.1
KAPPA
the force constant for the wall
=100000
EXP
the powers for the walls
=2
LABEL
a label for the action so that its output can be referenced in the input to other actions
=uwall
# Avoid other structures
Q6
Calculate sixth order Steinhardt parameters. This action is a shortcut. More details
SPECIES
this keyword is used for colvars such as coordination number
=1-288:3
SWITCH
the switching function that it used in the construction of the contact matrix
={CUBIC D_0=0.3 D_MAX=0.35}
VMEAN
calculate the norm of the mean vector
LABEL
a label for the action so that its output can be referenced in the input to other actions
=q6
diff2:
MATHEVAL
An alias to the CUSTOM function that can also be used to calaculate combinations of variables using a custom expression. More details
ARG
the input to this function
=q6.vmean,refcv.mean
FUNC
the function you wish to evaluate
=((x-0.0702)/(0.429-0.0702)-(y-0.254)/(0.812-0.254
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=NO
UPPER_WALLS
Defines a wall for the value of one or more collective variables, More details
ARG
the arguments on which the bias is acting
=diff2
AT
the positions of the wall
=0.1
KAPPA
the force constant for the wall
=100000
EXP
the powers for the walls
=2
LABEL
a label for the action so that its output can be referenced in the input to other actions
=uwall2
PRINT
Print quantities to a file. More details
STRIDE
the frequency with which the quantities of interest should be output
=500
ARG
the input for this action is the scalar output from one or more other actions
FILE
the name of the file on which to output these quantities
=COLVAR