Project ID: plumID:21.039
Source: data/chignolin/1_opes_multi/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
# vim: ft=plumed
Enables syntax highlighting for PLUMED files in vim. See here for more details. 


####################################
#  >> Chignolin <<
#  OPES Multicanonical simulation
####################################

#generated with gmx editconf -f topol.tpr -o reference.pdb, see https://www.plumed.org/doc-v2.7/user-doc/html/_m_o_l_i_n_f_o.html

MOLINFO
This command is used to provide information on the molecules that are present in your system. More details
 
MOLTYPE
 what kind of molecule is contained in the pdb file - usually not needed since protein/RNA/DNA are compatible
=protein 
STRUCTURE
a file in pdb format containing a reference structure
=chignolin-ref.pdb

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-166

# Define CVs

# Select Calpha
PROTEIN: 
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-166
CA: 
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
=5,26,47,61,73,88,102,109,123,147
# RMSD
rmsd_ca: 
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
=chignolin-ca.pdb 
TYPE
 the manner in which RMSD alignment is performed
=OPTIMAL
# END-TO-END DISTANCE
end: 
DISTANCE
Calculate the distance between a pair of atoms. More details
 
ATOMS
the pair of atom that we are calculating the distance between
=5,147
# HBONDS
hbonds: 
CONTACTMAP
Calculate the distances between a number of pairs of atoms and transform each distance by a switching function. More details
 
ATOMS1
the atoms involved in each of the contacts you wish to calculate
=23,145 
ATOMS2
the atoms involved in each of the contacts you wish to calculate
=45,120 
ATOMS3
the atoms involved in each of the contacts you wish to calculate
=56,100 
ATOMS4
the atoms involved in each of the contacts you wish to calculate
=56,107 
SWITCH
The switching functions to use for each of the contacts in your map
={RATIONAL R_0=0.4 NN=6 MM=8} 
SUM
 calculate the sum of all the contacts in the input

ene: 
ENERGY
Calculate the total potential energy of the simulation box. More details


# Define Bias
mc: 
ECV_MULTITHERMAL
Expand a simulation to sample multiple temperatures simultaneously. More details
 
ARG
the label of the internal energy of the system
=ene 
TEMP_MIN
the minimum of the temperature range
=270 
TEMP_MAX
the maximum of the temperature range
=700 
opes: 
OPES_EXPANDED
On-the-fly probability enhanced sampling with expanded ensembles for the target distribution. More details
 
PACE
how often the bias is updated
=500 
ARG
the label of the ECVs that define the expansion
=mc 
WALKERS_MPI
 switch on MPI version of multiple walkers
 

# PRINT

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

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
=5000

ENDPLUMED
Terminate plumed input. More details