PLUMED-NEST

Project ID: plumID:24.008
Source: 5FU/Volume_MetaD_Input/plumed_run.dat
Originally used with PLUMED version: 2.4
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

MOLINFO
This command is used to provide information on the molecules that are present in your system. More details
 
STRUCTURE
a file in pdb format containing a reference structure
=../Structure/yCD_5FU.pdb
The MOLINFO action with label  calculates something
WHOLEMOLECULES
This action is used to rebuild molecules that can become split by the periodic boundary conditions. More details
 
ENTITY0
the atoms that make up a molecule that you wish to align
=1-2423 
ENTITY1
the atoms that make up a molecule that you wish to align
=2424-4846 
ENTITY2
the atoms that make up a molecule that you wish to align
=4847-4858 
ENTITY3
the atoms that make up a molecule that you wish to align
=4859-4870
FIT_TO_TEMPLATE
This action is used to align a molecule to a template. More details
 
REFERENCE
a file in pdb format containing the reference structure and the atoms involved in the CV
=../Structure/heavy_atoms.pdb 
TYPE
 the manner in which RMSD alignment is performed
=OPTIMAL

# Group definitions
sph: 
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
 
NDX_FILE
the name of index file (gromacs syntax)
=../Structure/index.ndx 
NDX_GROUP
the name of the group to be imported (gromacs syntax) - first group found is used by default
=ChainA-H_ZN-ChainA 
The GROUP action with label sph calculates somethingurl: 
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
 
NDX_FILE
the name of index file (gromacs syntax)
=../Structure/index.ndx 
NDX_GROUP
the name of the group to be imported (gromacs syntax) - first group found is used by default
=UR1_noH
The GROUP action with label url calculates somethinglig: 
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
 
NDX_FILE
the name of index file (gromacs syntax)
=../Structure/index.ndx 
NDX_GROUP
the name of the group to be imported (gromacs syntax) - first group found is used by default
=UR1
The GROUP action with label lig calculates somethingref: 
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
 
NDX_FILE
the name of index file (gromacs syntax)
=../Structure/index.ndx 
NDX_GROUP
the name of the group to be imported (gromacs syntax) - first group found is used by default
=Protein-H_ZN
# Group sys 
The GROUP action with label ref calculates somethingsys: 
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
 
NDX_FILE
the name of index file (gromacs syntax)
=../Structure/index.ndx 
NDX_GROUP
the name of the group to be imported (gromacs syntax) - first group found is used by default
=Protein_ZN_UR1_UR2

The GROUP action with label sys calculates something
WRAPAROUND
Rebuild periodic boundary conditions around chosen atoms. More details
 
ATOMS
wrapped atoms
=lig 
AROUND
reference atoms
=sph

# Now we need to remake whole molecules
WHOLEMOLECULES
This action is used to rebuild molecules that can become split by the periodic boundary conditions. More details
 
ENTITY0
the atoms that make up a molecule that you wish to align
=1-2423 
ENTITY1
the atoms that make up a molecule that you wish to align
=2424-4846 
ENTITY2
the atoms that make up a molecule that you wish to align
=4847-4858 
ENTITY3
the atoms that make up a molecule that you wish to align
=4859-4870 
ENTITY4
the atoms that make up a molecule that you wish to align
=4957-4959 
ENTITY5
the atoms that make up a molecule that you wish to align
=4960-4962
FIT_TO_TEMPLATE
This action is used to align a molecule to a template. More details
 
REFERENCE
a file in pdb format containing the reference structure and the atoms involved in the CV
=../Structure/chainA_template.pdb 
TYPE
 the manner in which RMSD alignment is performed
=OPTIMAL

# Fix
fix: 
FIXEDATOM
Add a virtual atom in a fixed position. More details
 
AT
coordinates of the virtual atom
=-1.101616,-0.678271,0.017793

The FIXEDATOM action with label fix calculates somethingsph_center: 
COM
Calculate the center of mass for a group of atoms. More details
 
ATOMS
the list of atoms which are involved the virtual atom's definition
=sph
The COM action with label sph_center calculates somethingurl_center: 
COM
Calculate the center of mass for a group of atoms. More details
 
ATOMS
the list of atoms which are involved the virtual atom's definition
=url

The COM action with label url_center calculates somethingdis: 
DISTANCE
Calculate the distance between a pair of atoms. More details
 
ATOMS
the pair of atom that we are calculating the distance between
=sph_center,fix


# Define volume based CVs 
The DISTANCE action with label dis calculates the following quantities:
 Quantity   
 Description  
dis.value
the DISTANCE between this pair of atomsabs: 
DISTANCE
Calculate the distance between a pair of atoms. More details
 
ATOMS
the pair of atom that we are calculating the distance between
=url_center,sph_center 
COMPONENTS
 calculate the x, y and z components of the distance separately and store them as label
The DISTANCE action with label abs calculates the following quantities:
 Quantity   
 Description  
abs.x
the x-component of the vector connecting the two atoms
abs.y
the y-component of the vector connecting the two atoms
abs.z
the z-component of the vector connecting the two atoms
abs.value
the DISTANCE between this pair of atomsrho: 
DISTANCE
Calculate the distance between a pair of atoms. More details
 
ATOMS
the pair of atom that we are calculating the distance between
=url_center,sph_center
The DISTANCE action with label rho calculates the following quantities:
 Quantity   
 Description  
rho.value
the DISTANCE between this pair of atomsrel_x: 
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
=abs.z,abs.y 
FUNC
the function you wish to evaluate
=(x-y)/sqrt(2 
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=NO
The MATHEVAL action with label rel_x calculates somethingrel_y: 
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
=abs.x 
FUNC
the function you wish to evaluate
=x 
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=NO
The MATHEVAL action with label rel_y calculates somethingrel_z: 
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
=abs.z,abs.y 
FUNC
the function you wish to evaluate
=(x+y)/sqrt(2 
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=NO

The MATHEVAL action with label rel_z calculates somethingtheta: 
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
=rel_z,rho 
FUNC
the function you wish to evaluate
=acos(x/y 
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=0.,pi
The MATHEVAL action with label theta calculates somethingphi: 
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
=rel_x,rel_y 
FUNC
the function you wish to evaluate
=atan2(y,x 
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=-pi,pi

The MATHEVAL action with label phi calculates somethingrmsd: 
RMSD
Calculate the RMSD with respect to a reference structure. More details
 
REFERENCE
a file in pdb format containing the reference structure and the atoms involved in the CV
=../Structure/heavy_atoms_5FU1.pdb 
TYPE
 the manner in which RMSD alignment is performed
=OPTIMAL

# Restraining potential of the sphere
The RMSD action with label rmsd calculates somethingrestr: 
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
=rho 
AT
the positions of the wall
=2.8 
KAPPA
the force constant for the wall
=200
The UPPER_WALLS action with label restr calculates somethingrestr2: 
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
=rmsd 
AT
the positions of the wall
=0.1 
KAPPA
the force constant for the wall
=200

# Coordination number
The UPPER_WALLS action with label restr2 calculates somethingc: 
COORDINATION
Calculate coordination numbers. More details
 
GROUPA
First list of atoms
=url 
GROUPB
Second list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)
=sph 
R_0
The r_0 parameter of the switching function
=0.45

# Metadynamics
The COORDINATION action with label c calculates the following quantities:
 Quantity   
 Description  
c.value
the value of the coordination
METAD
Used to performed metadynamics on one or more collective variables. More details
 ...
ARG
the input for this action is the scalar output from one or more other actions
=rho,theta,phi
GRID_MIN
the lower bounds for the grid
=0,0.,-pi
GRID_MAX
the upper bounds for the grid
=5,pi,pi
SIGMA
the widths of the Gaussian hills
=0.1,pi/16.,pi/8
HEIGHT
the heights of the Gaussian hills
=1.2
PACE
the frequency for hill addition
=500
BIASFACTOR
use well tempered metadynamics and use this bias factor
=20
TEMP
the system temperature - this is only needed if you are doing well-tempered metadynamics
=300
LABEL
a label for the action so that its output can be referenced in the input to other actions
=metad
CALC_RCT
 calculate the c(t) reweighting factor and use that to obtain the normalized bias [rbias=bias-rct]
... METAD

The METAD action with label metad calculates the following quantities:
 Quantity   
 Description  
metad.bias
the instantaneous value of the bias potential
metad.rbias
the instantaneous value of the bias normalized using the c(t) reweighting factor [rbias=bias-rct]
metad.rct
the reweighting factor c(t)
PRINT
Print quantities to a file. More details
 
ARG
the input for this action is the scalar output from one or more other actions
=rho,dis 
FILE
the name of the file on which to output these quantities
=distance.dat 
STRIDE
 the frequency with which the quantities of interest should be output
=500
PRINT
Print quantities to a file. More details
 
ARG
the input for this action is the scalar output from one or more other actions
=metad 
FILE
the name of the file on which to output these quantities
=metad_data.dat 
STRIDE
 the frequency with which the quantities of interest should be output
=500
PRINT
Print quantities to a file. More details
 
ARG
the input for this action is the scalar output from one or more other actions
=restr 
FILE
the name of the file on which to output these quantities
=sphere_restraint.dat 
STRIDE
 the frequency with which the quantities of interest should be output
=500
PRINT
Print quantities to a file. More details
 
ARG
the input for this action is the scalar output from one or more other actions
=abs.x,abs.y,abs.z,rel_x,rel_y,rel_z 
FILE
the name of the file on which to output these quantities
=xyz_coord.dat 
STRIDE
 the frequency with which the quantities of interest should be output
=500
PRINT
Print quantities to a file. More details
 
ARG
the input for this action is the scalar output from one or more other actions
=rho,theta,phi 
FILE
the name of the file on which to output these quantities
=rtp_coord.dat 
STRIDE
 the frequency with which the quantities of interest should be output
=500
PRINT
Print quantities to a file. More details
 
ARG
the input for this action is the scalar output from one or more other actions
=c 
FILE
the name of the file on which to output these quantities
=all_coordination_45.dat 
STRIDE
 the frequency with which the quantities of interest should be output
=500
PRINT
Print quantities to a file. More details
 
ARG
the input for this action is the scalar output from one or more other actions
=restr2 
FILE
the name of the file on which to output these quantities
=rmsd_restraint.dat 
STRIDE
 the frequency with which the quantities of interest should be output
=500


FLUSH
This command instructs plumed to flush all the open files with a user specified frequency. More details
 
STRIDE
the frequency with which all the open files should be flushed
=500