Project ID: plumID:23.046
Source: github/MccJ25/plumed.dat
Originally used with PLUMED version: 2.7.1
Stable: zipped raw stdout - zipped raw stderr - stderr
Master: zipped raw stdout - zipped raw stderr - stderr
# This is the plumed file to use on linear native MJ25 - using 2 CVs #RESTART
# unitsUNITSThis command sets the internal units for the code. More detailsLENGTH=Athe units of lengths
# MoleculeWHOLEMOLECULESThis action is used to rebuild molecules that can become split by the periodic boundary conditions. More detailsENTITY0=1-292the atoms that make up a molecule that you wish to align
# distances # first is 1N-8CD and the others are for the CA of resid19 to resids 1to8DISTANCECalculate the distance between a pair of atoms. More detailsATOMS=5,246the pair of atom that we are calculating the distance betweenLABEL=d1a label for the action so that its output can be referenced in the input to other actionsDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=12,246the pair of atom that we are calculating the distance betweenLABEL=d2a label for the action so that its output can be referenced in the input to other actionsDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=19,246the pair of atom that we are calculating the distance betweenLABEL=d3a label for the action so that its output can be referenced in the input to other actionsDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=29,246the pair of atom that we are calculating the distance betweenLABEL=d4a label for the action so that its output can be referenced in the input to other actionsDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=36,246the pair of atom that we are calculating the distance betweenLABEL=d5a label for the action so that its output can be referenced in the input to other actionsDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=53,246the pair of atom that we are calculating the distance betweenLABEL=d6a label for the action so that its output can be referenced in the input to other actionsDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=77,246the pair of atom that we are calculating the distance betweenLABEL=d7a label for the action so that its output can be referenced in the input to other actionsDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=83,246the pair of atom that we are calculating the distance betweenLABEL=d8a label for the action so that its output can be referenced in the input to other actions
CV1 :DISTANCECalculate the distance between a pair of atoms. More detailsATOMS=1,91 CV2 :the pair of atom that we are calculating the distance betweenCOMBINE...Calculate a polynomial combination of a set of other variables. More detailsARG=d1,d2,d3,d4,d5,d6,d7,d8the input to this functionCOEFFICIENTS=-0.3221,-0.0361,-0.1063,-0.1538,-0.0913,-0.5176,-0.6589,-0.3867the coefficients of the arguments in your functionPERIODIC=NO ...if the output of your function is periodic then you should specify the periodicity of the function
# activade well-tempered metadynamics metad :METADUsed to performed metadynamics on one or more collective variables. More detailsARG=CV1,CV2 ... # Deposit a Gaussian every 500 time steps (= 1ps), with initial height # equal to 5 kJ/mol and bias factor equal to 80 # This version we use a double value for the height and bias factorthe input for this action is the scalar output from one or more other actionsPACE=500the frequency for hill additionHEIGHT=5the heights of the Gaussian hillsBIASFACTOR=80use well tempered metadynamics and use this bias factorTEMP=300 # Gaussian width (sigma) should be chosen based on the CV fluctuations in unbiased run # try 1/2 or 1/3 of the estimated fluctuations # standard deviation for distance 1N - 8N was 0.91 # standard deviation for CA distances was between 6.84-7.06the system temperature - this is only needed if you are doing well-tempered metadynamicsSIGMA=0.5,3.5 # Gaussians will be written to file and also stored on gridthe widths of the Gaussian hillsFILE=HILLSa file in which the list of added hills is storedGRID_WSTRIDE=100000 #GRID_RFILE=grid1.datwrite the grid to a file every N stepsGRID_WFILE=grid1.datthe file on which to write the gridGRID_NOSPLINEdon't use spline interpolation with gridsGRID_MIN=0,-200the lower bounds for the gridGRID_MAX=35,0the upper bounds for the gridGRID_BIN=350,500 ... # Print both collective variables on COLVAR file every 100 stepsthe number of bins for the gridPrint quantities to a file. More detailsARG=d1,d2,d3,d4,d5,d6,d7,d8,CV1,CV2,metad.biasthe input for this action is the scalar output from one or more other actionsFILE=COLVARthe name of the file on which to output these quantitiesSTRIDE=100the frequency with which the quantities of interest should be output