Project ID: plumID:20.009
Source: Ub4/plumed-cv.dat
Originally used with PLUMED version: 2.6
Stable: zipped raw stdout - zipped raw stderr - stderr
Master: zipped raw stdout - zipped raw stderr - stderr
#SETTINGS NREPLICAS=2MOLINFOThis command is used to provide information on the molecules that are present in your system. More detailsSTRUCTURE=template.pdba file in pdb format containing a reference structureWHOLEMOLECULESThis action is used to rebuild molecules that can become split by the periodic boundary conditions. More detailsENTITY0=1-652the atoms that make up a molecule that you wish to align
#ubi1 a:CENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=1,3,5,7,11,13,16,18,20,22,23,26,28,30,32,34,36,38,40,42,44,46,48,50,52,54,56,59,60,63,65,67,69,72,74,75,77,79,81,83,85,87,90,92,94,98,99,100,103,105,107,109,111,112,115,117,119,121,123,127,129,131,133,136,138,140,142,144,148,150the group of atoms that you are calculating the Gyration Tensor forNOPBCa1:ignore the periodic boundary conditions when calculating distancesCENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=1,3,5,7,11,13,16,18,20,22,23,26,28,30,32,34,36,38,40,42,44,46,48,50,52,54,56,59,60,63,65,67,69,72,74,75the group of atoms that you are calculating the Gyration Tensor forNOPBCa2:ignore the periodic boundary conditions when calculating distancesCENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=77,79,81,83,85,87,90,92,94,98,99,100,103,105,107,109,111,112,115,117,119,121,123,127,129,131,133,136,138,140,142,144,148,150the group of atoms that you are calculating the Gyration Tensor forNOPBC#ubi2 b:ignore the periodic boundary conditions when calculating distancesCENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=164,166,168,170,174,176,179,181,183,185,186,189,191,193,195,197,199,201,203,205,207,209,211,213,215,217,219,222,223,226,228,230,232,235,237,238,240,242,244,246,248,250,253,255,257,261,262,263,266,268,270,272,274,275,278,280,282,284,286,290,292,294,296,299,301,303,305,307,311,313the group of atoms that you are calculating the Gyration Tensor forNOPBCb1:ignore the periodic boundary conditions when calculating distancesCENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=164,166,168,170,174,176,179,181,183,185,186,189,191,193,195,197,199,201,203,205,207,209,211,213,215,217,219,222,223,226,228,230,232,235,237,238the group of atoms that you are calculating the Gyration Tensor forNOPBCb2:ignore the periodic boundary conditions when calculating distancesCENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=,240,242,244,246,248,250,253,255,257,261,262,263,266,268,270,272,274,275,278,280,282,284,286,290,292,294,296,299,301,303,305,307,311,313the group of atoms that you are calculating the Gyration Tensor forNOPBC#ubi3 c:ignore the periodic boundary conditions when calculating distancesCENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=327,329,331,333,337,339,342,344,346,348,349,352,354,356,358,360,362,364,366,368,370,372,374,376,378,380,382,385,386,389,391,393,395,398,400,401,403,405,407,409,411,413,416,418,420,424,425,426,429,431,433,435,437,438,441,443,445,447,449,453,455,457,459,462,464,466,468,470,474,476the group of atoms that you are calculating the Gyration Tensor forNOPBCc1:ignore the periodic boundary conditions when calculating distancesCENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=327,329,331,333,337,339,342,344,346,348,349,352,354,356,358,360,362,364,366,368,370,372,374,376,378,380,382,385,386,389,391,393,395,398,400,401the group of atoms that you are calculating the Gyration Tensor forNOPBCc2:ignore the periodic boundary conditions when calculating distancesCENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=403,405,407,409,411,413,416,418,420,424,425,426,429,431,433,435,437,438,441,443,445,447,449,453,455,457,459,462,464,466,468,470,474,476the group of atoms that you are calculating the Gyration Tensor forNOPBC#ubi4 e:ignore the periodic boundary conditions when calculating distancesCENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=490,492,494,496,500,502,505,507,509,511,512,515,517,519,521,523,525,527,529,531,533,535,537,539,541,543,545,548,549,552,554,556,558,561,563,564,566,568,570,572,574,576,579,581,583,587,588,589,592,594,596,598,600,601,604,606,608,610,612,616,618,620,622,625,627,629,631,633,637,639the group of atoms that you are calculating the Gyration Tensor forNOPBCe1:ignore the periodic boundary conditions when calculating distancesCENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=490,492,494,496,500,502,505,507,509,511,512,515,517,519,521,523,525,527,529,531,533,535,537,539,541,543,545,548,549,552,554,556,558,561,563,564the group of atoms that you are calculating the Gyration Tensor forNOPBCe2:ignore the periodic boundary conditions when calculating distancesCENTERCalculate the center for a group of atoms, with arbitrary weights. More detailsATOMS=566,568,570,572,574,576,579,581,583,587,588,589,592,594,596,598,600,601,604,606,608,610,612,616,618,620,622,625,627,629,631,633,637,639the group of atoms that you are calculating the Gyration Tensor forNOPBCignore the periodic boundary conditions when calculating distances
# distance between ubiquitins d1:DISTANCECalculate the distance between a pair of atoms. More detailsATOMS=a,bthe pair of atom that we are calculating the distance betweenNOPBCd2:ignore the periodic boundary conditions when calculating distancesDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=a,cthe pair of atom that we are calculating the distance betweenNOPBCd3:ignore the periodic boundary conditions when calculating distancesDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=a,ethe pair of atom that we are calculating the distance betweenNOPBCd4:ignore the periodic boundary conditions when calculating distancesDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=b,cthe pair of atom that we are calculating the distance betweenNOPBCd5:ignore the periodic boundary conditions when calculating distancesDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=b,ethe pair of atom that we are calculating the distance betweenNOPBCd6:ignore the periodic boundary conditions when calculating distancesDISTANCECalculate the distance between a pair of atoms. More detailsATOMS=c,ethe pair of atom that we are calculating the distance betweenNOPBCignore the periodic boundary conditions when calculating distances
#relativ orientation tor1:TORSIONCalculate a torsional angle. More detailsATOMS=a1,a2,b1,b2the four atoms involved in the torsional angleNOPBCtor2:ignore the periodic boundary conditions when calculating distancesTORSIONCalculate a torsional angle. More detailsATOMS=a1,a2,c1,c2the four atoms involved in the torsional angleNOPBCtor3:ignore the periodic boundary conditions when calculating distancesTORSIONCalculate a torsional angle. More detailsATOMS=a1,a2,e1,e2the four atoms involved in the torsional angleNOPBCtor4:ignore the periodic boundary conditions when calculating distancesTORSIONCalculate a torsional angle. More detailsATOMS=b1,b2,c1,c2the four atoms involved in the torsional angleNOPBCtor5:ignore the periodic boundary conditions when calculating distancesTORSIONCalculate a torsional angle. More detailsATOMS=b1,b2,e1,e2the four atoms involved in the torsional angleNOPBCtor6:ignore the periodic boundary conditions when calculating distancesTORSIONCalculate a torsional angle. More detailsATOMS=c1,c2,e1,e2the four atoms involved in the torsional angleNOPBCignore the periodic boundary conditions when calculating distances
#global radius
rg:GYRATIONCalculate the radius of gyration, or other properties related to it. More detailsATOMS=1,3,5,7,11,13,16,18,20,22,23,26,28,30,32,34,36,38,40,42,44,46,48,50,52,54,56,59,60,63,65,67,69,72,74,75,77,79,81,83,85,87,90,92,94,98,99,100,103,105,107,109,111,112,115,117,119,121,123,127,129,131,133,136,138,140,142,144,148,150,152,154,157,159,162,163,164,166,168,170,174,176,179,181,183,185,186,189,191,193,195,197,199,201,203,205,207,209,211,213,215,217,219,222,223,226,228,230,232,235,237,238,240,242,244,246,248,250,253,255,257,261,262,263,266,268,270,272,274,275,278,280,282,284,286,290,292,294,296,299,301,303,305,307,311,313,315,317,320,322,325,326,327,329,331,333,337,339,342,344,346,348,349,352,354,356,358,360,362,364,366,368,370,372,374,376,378,380,382,385,386,389,391,393,395,398,400,401,403,405,407,409,411,413,416,418,420,424,425,426,429,431,433,435,437,438,441,443,445,447,449,453,455,457,459,462,464,466,468,470,474,476,478,480,483,485,488,489,490,492,494,496,500,502,505,507,509,511,512,515,517,519,521,523,525,527,529,531,533,535,537,539,541,543,545,548,549,552,554,556,558,561,563,564,566,568,570,572,574,576,579,581,583,587,588,589,592,594,596,598,600,601,604,606,608,610,612,616,618,620,622,625,627,629,631,633,637,639,641,643,646,648,651,652the group of atoms that you are calculating the Gyration Tensor forNOPBCignore the periodic boundary conditions when calculating distances
#linker angleALPHABETA...Calculate the alpha beta CV More detailsLABEL=linktor1a label for the action so that its output can be referenced in the input to other actionsREFERENCE=0the reference values for each of the torsional anglesATOMS1=150,152,154,157the atoms involved for each of the torsions you wish to calculateATOMS2=152,154,157,159the atoms involved for each of the torsions you wish to calculateATOMS3=154,157,159,162the atoms involved for each of the torsions you wish to calculateATOMS4=157,159,162,163the atoms involved for each of the torsions you wish to calculateATOMS5=159,162,163,164the atoms involved for each of the torsions you wish to calculateATOMS6=162,163,164,166 ... ALPHABETAthe atoms involved for each of the torsions you wish to calculateALPHABETA...Calculate the alpha beta CV More detailsLABEL=linktor2a label for the action so that its output can be referenced in the input to other actionsREFERENCE=0the reference values for each of the torsional anglesATOMS1=313,315,317,320the atoms involved for each of the torsions you wish to calculateATOMS2=315,317,320,322the atoms involved for each of the torsions you wish to calculateATOMS3=317,320,322,325the atoms involved for each of the torsions you wish to calculateATOMS4=320,322,325,326the atoms involved for each of the torsions you wish to calculateATOMS5=322,325,326,327the atoms involved for each of the torsions you wish to calculateATOMS6=325,326,327,329 ... ALPHABETAthe atoms involved for each of the torsions you wish to calculateALPHABETA...Calculate the alpha beta CV More detailsLABEL=linktor3a label for the action so that its output can be referenced in the input to other actionsREFERENCE=0the reference values for each of the torsional anglesATOMS1=476,478,480,483the atoms involved for each of the torsions you wish to calculateATOMS2=478,480,483,485the atoms involved for each of the torsions you wish to calculateATOMS3=480,483,485,488the atoms involved for each of the torsions you wish to calculateATOMS4=483,485,488,489the atoms involved for each of the torsions you wish to calculateATOMS5=485,488,489,490the atoms involved for each of the torsions you wish to calculateATOMS6=488,489,490,492 ... ALPHABETAthe atoms involved for each of the torsions you wish to calculateSAXS...Calculates SAXS intensity. More detailsLABEL=test2a label for the action so that its output can be referenced in the input to other actionsATOMS=1-652The atoms to be included in the calculation, eNOPBCIgnore the periodic boundary conditions when calculating distancesMARTINICalculate SAXS for a Martini modelQVALUE1=0.02Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE2=0.03Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE3=0.04Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE4=0.05Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE5=0.06Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE6=0.07Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE7=0.08Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE8=0.09Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE9=0.1Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE10=0.11Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE11=0.12Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE12=0.13Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE13=0.14Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE14=0.15Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE15=0.16Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE16=0.17Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE17=0.18Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE18=0.19Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE19=0.2Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,QVALUE20=0.21Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2,EXPINT1=5.24833Add an experimental value for each q valueEXPINT2=4.39027Add an experimental value for each q valueEXPINT3=3.50363Add an experimental value for each q valueEXPINT4=2.70379Add an experimental value for each q valueEXPINT5=2.07676Add an experimental value for each q valueEXPINT6=1.62809Add an experimental value for each q valueEXPINT7=1.31061Add an experimental value for each q valueEXPINT8=1.07411Add an experimental value for each q valueEXPINT9=0.887992Add an experimental value for each q valueEXPINT10=0.738255Add an experimental value for each q valueEXPINT11=0.617356Add an experimental value for each q valueEXPINT12=0.518766Add an experimental value for each q valueEXPINT13=0.436834Add an experimental value for each q valueEXPINT14=0.367909Add an experimental value for each q valueEXPINT15=0.310039Add an experimental value for each q valueEXPINT16=0.261567Add an experimental value for each q valueEXPINT17=0.220234Add an experimental value for each q valueEXPINT18=0.183526Add an experimental value for each q valueEXPINT19=0.149622Add an experimental value for each q valueEXPINT20=0.117966 ... SAXSAdd an experimental value for each q value
lsaxs:STATSCalculates statistical properties of a set of collective variables with respect to a set of reference values. More detailsARG=(test2\.qthe input for this action is the scalar output from one or more other actionsPARARG=(test2\.expthe input for this action is the scalar output from one or more other actions without derivativesPBMETAD...Used to performed Parallel Bias metadynamics. More detailsLABEL=biasa label for the action so that its output can be referenced in the input to other actionsARG=d1,d2,d3,d4,d5,d6,linktor1,linktor2,linktor3,rg,tor1,tor2,tor3,tor4,tor5,tor6the input for this action is the scalar output from one or more other actionsHEIGHT=0.075the height of the Gaussian hills, one for all biasesPACE=200the frequency for hill addition, one for all biasesBIASFACTOR=40use well tempered metadynamics with this bias factor, one for all biasesADAPTIVE=DIFFuse a geometric (=GEOM) or diffusion (=DIFF) based hills width schemeGRID_WFILES=GRID.0,GRID.1,GRID.2,GRID.3,GRID.4,GRID.5,GRID.6,GRID.7,GRID.8,GRID.9,GRID.10,GRID.11,GRID.12,GRID.13,GRID.14,GRID.15dump grid for the bias, default names are used if GRID_WSTRIDE is used without GRID_WFILESGRID_WSTRIDE=2000 #GRID_RFILES=GRID.0,GRID.1,GRID.2,GRID.3,GRID.4,GRID.5,GRID.6,GRID.7,GRID.8,GRID.9,GRID.10,GRID.11,GRID.12,GRID.13,GRID.14,GRID.15frequency for dumping the gridSIGMA=1000the widths of the Gaussian hillsSIGMA_MIN=0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005the lower bounds for the sigmas (in CV units) when using adaptive hillsSIGMA_MAX=0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2the upper bounds for the sigmas (in CV units) when using adaptive hillsGRID_MIN=0,0,0,0,0,0,0,0,0,0,-pi,-pi,-pi,-pi,-pi,-pithe lower bounds for the gridGRID_MAX=7,14,18,7,14,7,6,6,6,7,pi,pi,pi,pi,pi,pithe upper bounds for the gridWALKERS_MPISwitch on MPI version of multiple walkers - not compatible with WALKERS_* options other than WALKERS_DIRUPDATE_FROM=2000 ... PBMETADOnly update this action from this timeMETAINFERENCE...Calculates the Metainference energy for a set of experimental data. More detailsSCALEDATASet to TRUE if you want to sample a scaling factor common to all values and replicasNOISETYPE=MOUTLIERSfunctional form of the noise (GAUSS,MGAUSS,OUTLIERS,MOUTLIERS,GENERIC)AVERAGING=100Stride for calculation of averaged weights and sigma_meanARG=(test2\.q-.*),bias.biasthe input for this action is the scalar output from one or more other actionsREWEIGHTsimple REWEIGHT using the latest ARG as energyPARARG=(test2\.expreference values for the experimental data, these can be provided as arguments without derivativesSCALE_PRIOR=FLATeither FLAT or GAUSSIANADDOFFSETSet to TRUE if you want to sample an offset common to all values and replicasOFFSET_PRIOR=FLATeither FLAT or GAUSSIANOFFSET0=0initial value of the offsetOFFSET_MIN=-10minimum value of the offsetOFFSET_MAX=10maximum value of the offsetDOFFSET=0.01maximum MC move of the offsetSCALE0=1.00initial value of the scaling factorSCALE_MIN=0.8minimum value of the scaling factorSCALE_MAX=1.2maximum value of the scaling factorDSCALE=0.001maximum MC move of the scaling factorSIGMA0=0.01initial value of the uncertainty parameterSIGMA_MIN=0minimum value of the uncertainty parameterSIGMA_MAX=0.05maximum value of the uncertainty parameterOPTSIGMAMEAN=SEMSet to NONE/SEM to manually set sigma mean, or to estimate it on the flyLABEL=bqa label for the action so that its output can be referenced in the input to other actionsWRITE_STRIDE=10000write the status to a file every N steps, this can be used for restart/continuationSTRIDE=5 ... METAINFERENCEthe frequency with which the forces due to the bias should be calculated
ensaxs:ENSEMBLECalculates the replica averaging of a collective variable over multiple replicas. More detailsARG=(test2\.q-.*),bias.biasthe input for this action is the scalar output from one or more other actionsREWEIGHTstatsq:simple REWEIGHT using the latest ARG as energySTATSCalculates statistical properties of a set of collective variables with respect to a set of reference values. More detailsARG=(ensaxs\.test2\.qthe input for this action is the scalar output from one or more other actionsPARARG=(test2\.expthe input for this action is the scalar output from one or more other actions without derivativesFLUSHThis command instructs plumed to flush all the open files with a user specified frequency. More detailsSTRIDE=1000the frequency with which all the open files should be flushedPrint quantities to a file. More detailsARG=(test2\.q-.*),(test2\.expthe input for this action is the scalar output from one or more other actionsFILE=CV.datthe name of the file on which to output these quantitiesSTRIDE=100the frequency with which the quantities of interest should be outputPrint quantities to a file. More detailsFILE=RESTRAINTSthe name of the file on which to output these quantitiesARG=bqthe input for this action is the scalar output from one or more other actionsSTRIDE=100the frequency with which the quantities of interest should be outputPrint quantities to a file. More detailsFILE=q.datthe name of the file on which to output these quantitiesARG=statsq.*,lsaxsthe input for this action is the scalar output from one or more other actionsSTRIDE=100the frequency with which the quantities of interest should be outputPrint quantities to a file. More detailsFILE=Analyse.datthe name of the file on which to output these quantitiesARG=d1,d2,d3,d4,d5,d6,linktor1,linktor2,linktor3,rg,tor1,tor2,tor3,tor4,tor5,tor6the input for this action is the scalar output from one or more other actionsSTRIDE=100the frequency with which the quantities of interest should be outputPrint quantities to a file. More detailsFILE=ENSAXSthe name of the file on which to output these quantitiesARG=(ensaxs\.test2\.qthe input for this action is the scalar output from one or more other actionsSTRIDE=100the frequency with which the quantities of interest should be output
enANA:ENSEMBLECalculates the replica averaging of a collective variable over multiple replicas. More detailsARG=d1,d2,d3,d4,d5,d6,linktor1,linktor2,linktor3,rg,tor1,tor2,tor3,tor4,tor5,tor6,bias.biasthe input for this action is the scalar output from one or more other actionsREWEIGHTsimple REWEIGHT using the latest ARG as energyPrint quantities to a file. More detailsFILE=ENANAthe name of the file on which to output these quantitiesARG=enANAthe input for this action is the scalar output from one or more other actionsSTRIDE=100the frequency with which the quantities of interest should be output