PLUMED-NEST

Project ID: plumID:20.009
Source: Ub3/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

Click on the labels of the actions for more information on what each action computes

#SETTINGS NREPLICAS=2
MOLINFOThis command is used to provide information on the molecules that are present in your system. More details STRUCTUREa file in pdb format containing a reference structure=template.pdb

The MOLINFO action with label  calculates somethingWHOLEMOLECULESThis action is used to rebuild molecules that can become split by the periodic boundary conditions. More details ENTITY0the atoms that make up a molecule that you wish to align=1-489

#ubi1
a: CENTERCalculate the center for a group of atoms, with arbitrary weights. More details ATOMSthe group of atoms that you are calculating the Gyration Tensor for=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 NOPBC ignore the periodic boundary conditions when calculating distances
The CENTER action with label a calculates the following quantities:
 Quantity   
 Description  
a.value
the position of the center of massa1: CENTERCalculate the center for a group of atoms, with arbitrary weights. More details ATOMSthe group of atoms that you are calculating the Gyration Tensor for=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 NOPBC ignore the periodic boundary conditions when calculating distances
The CENTER action with label a1 calculates the following quantities:
 Quantity   
 Description  
a1.value
the position of the center of massa2: CENTERCalculate the center for a group of atoms, with arbitrary weights. More details ATOMSthe group of atoms that you are calculating the Gyration Tensor for=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 NOPBC ignore the periodic boundary conditions when calculating distances
#ubi2
The CENTER action with label a2 calculates the following quantities:
 Quantity   
 Description  
a2.value
the position of the center of massb: CENTERCalculate the center for a group of atoms, with arbitrary weights. More details ATOMSthe group of atoms that you are calculating the Gyration Tensor for=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 NOPBC ignore the periodic boundary conditions when calculating distances
The CENTER action with label b calculates the following quantities:
 Quantity   
 Description  
b.value
the position of the center of massb1: CENTERCalculate the center for a group of atoms, with arbitrary weights. More details ATOMSthe group of atoms that you are calculating the Gyration Tensor for=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 NOPBC ignore the periodic boundary conditions when calculating distances
The CENTER action with label b1 calculates the following quantities:
 Quantity   
 Description  
b1.value
the position of the center of massb2: CENTERCalculate the center for a group of atoms, with arbitrary weights. More details ATOMSthe group of atoms that you are calculating the Gyration Tensor for=,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 NOPBC ignore the periodic boundary conditions when calculating distances
#ubi3
The CENTER action with label b2 calculates the following quantities:
 Quantity   
 Description  
b2.value
the position of the center of massc: CENTERCalculate the center for a group of atoms, with arbitrary weights. More details ATOMSthe group of atoms that you are calculating the Gyration Tensor for=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 NOPBC ignore the periodic boundary conditions when calculating distances
The CENTER action with label c calculates the following quantities:
 Quantity   
 Description  
c.value
the position of the center of massc1: CENTERCalculate the center for a group of atoms, with arbitrary weights. More details ATOMSthe group of atoms that you are calculating the Gyration Tensor for=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 NOPBC ignore the periodic boundary conditions when calculating distances
The CENTER action with label c1 calculates the following quantities:
 Quantity   
 Description  
c1.value
the position of the center of massc2: CENTERCalculate the center for a group of atoms, with arbitrary weights. More details ATOMSthe group of atoms that you are calculating the Gyration Tensor for=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 NOPBC ignore the periodic boundary conditions when calculating distances

# distance between ubiquitins 
The CENTER action with label c2 calculates the following quantities:
 Quantity   
 Description  
c2.value
the position of the center of massd1: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=a,b NOPBC ignore the periodic boundary conditions when calculating distances
The DISTANCE action with label d1 calculates the following quantities:
 Quantity   
 Description  
d1.value
the DISTANCE between this pair of atomsd2: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=a,c NOPBC ignore the periodic boundary conditions when calculating distances
The DISTANCE action with label d2 calculates the following quantities:
 Quantity   
 Description  
d2.value
the DISTANCE between this pair of atomsd3: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=b,c NOPBC ignore the periodic boundary conditions when calculating distances
#relativ orientation
The DISTANCE action with label d3 calculates the following quantities:
 Quantity   
 Description  
d3.value
the DISTANCE between this pair of atomstor1: TORSIONCalculate one or multiple torsional angles. More details ATOMSthe four atoms involved in the torsional angle=a1,a2,b1,b2 NOPBC ignore the periodic boundary conditions when calculating distances
The TORSION action with label tor1 calculates the following quantities:
 Quantity   
 Description  
tor1.value
the TORSION involving these atomstor2: TORSIONCalculate one or multiple torsional angles. More details ATOMSthe four atoms involved in the torsional angle=a1,a2,c1,c2 NOPBC ignore the periodic boundary conditions when calculating distances
The TORSION action with label tor2 calculates the following quantities:
 Quantity   
 Description  
tor2.value
the TORSION involving these atomstor3: TORSIONCalculate one or multiple torsional angles. More details ATOMSthe four atoms involved in the torsional angle=b1,b2,c1,c2 NOPBC ignore the periodic boundary conditions when calculating distances

#global radius
The TORSION action with label tor3 calculates the following quantities:
 Quantity   
 Description  
tor3.value
the TORSION involving these atomsrg: GYRATIONCalculate the radius of gyration, or other properties related to it. More details ATOMSthe group of atoms that you are calculating the Gyration Tensor for=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 NOPBC ignore the periodic boundary conditions when calculating distances

#linker angle
The GYRATION action with label rg calculates the following quantities:
 Quantity   
 Description  
rg.value
the radius that was computed from the weightsALPHABETAMeasures a distance including pbc between the instantaneous values of a set of torsional angles and set of reference values. More details ...
LABELa label for the action so that its output can be referenced in the input to other actions=linktor1
REFERENCEthe reference values for each of the torsional angles=0
ATOMS1the atoms involved for each of the torsions you wish to calculate=150,152,154,157
ATOMS2the atoms involved for each of the torsions you wish to calculate=152,154,157,159
ATOMS3the atoms involved for each of the torsions you wish to calculate=154,157,159,162
ATOMS4the atoms involved for each of the torsions you wish to calculate=157,159,162,163
ATOMS5the atoms involved for each of the torsions you wish to calculate=159,162,163,164
ATOMS6the atoms involved for each of the torsions you wish to calculate=162,163,164,166
... ALPHABETA

The ALPHABETA action with label linktor1 calculates the following quantities:
 Quantity   
 Description  
linktor1.value
the alpha beta CVALPHABETAMeasures a distance including pbc between the instantaneous values of a set of torsional angles and set of reference values. More details ...
LABELa label for the action so that its output can be referenced in the input to other actions=linktor2
REFERENCEthe reference values for each of the torsional angles=0
ATOMS1the atoms involved for each of the torsions you wish to calculate=313,315,317,320
ATOMS2the atoms involved for each of the torsions you wish to calculate=315,317,320,322
ATOMS3the atoms involved for each of the torsions you wish to calculate=317,320,322,325
ATOMS4the atoms involved for each of the torsions you wish to calculate=320,322,325,326
ATOMS5the atoms involved for each of the torsions you wish to calculate=322,325,326,327
ATOMS6the atoms involved for each of the torsions you wish to calculate=325,326,327,329
... ALPHABETA

The ALPHABETA action with label linktor2 calculates the following quantities:
 Quantity   
 Description  
linktor2.value
the alpha beta CVSAXSCalculates SAXS intensity. More details ...
LABELa label for the action so that its output can be referenced in the input to other actions=test2
ATOMSThe atoms to be included in the calculation, e=1-489
NOPBC Ignore the periodic boundary conditions when calculating distances

MARTINI Calculate SAXS for a Martini model

QVALUE1Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.0178612
QVALUE2Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.029015
QVALUE3Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.0401689
QVALUE4Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.0513227
QVALUE5Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.0624766
QVALUE6Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.0736304
QVALUE7Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.0847843
QVALUE8Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.0959381
QVALUE9Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.107092
QVALUE10Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.118246
QVALUE11Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.129399
QVALUE12Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.140553
QVALUE13Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.151707
QVALUE14Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.162861
QVALUE15Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.174015
QVALUE16Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.185169
QVALUE17Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.196323
QVALUE18Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.207476
QVALUE19Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.21863
QVALUE20Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.229784
QVALUE21Selected scattering lengths in inverse angstroms are given as QVALUE1, QVALUE2, =0.240938


EXPINT1Add an experimental value for each q value=4.68143
EXPINT2Add an experimental value for each q value=4.11699
EXPINT3Add an experimental value for each q value=3.45783
EXPINT4Add an experimental value for each q value=2.76077
EXPINT5Add an experimental value for each q value=2.15094
EXPINT6Add an experimental value for each q value=1.67643
EXPINT7Add an experimental value for each q value=1.32305
EXPINT8Add an experimental value for each q value=1.06339
EXPINT9Add an experimental value for each q value=0.870742
EXPINT10Add an experimental value for each q value=0.721004
EXPINT11Add an experimental value for each q value=0.597613
EXPINT12Add an experimental value for each q value=0.493722
EXPINT13Add an experimental value for each q value=0.407671
EXPINT14Add an experimental value for each q value=0.337263
EXPINT15Add an experimental value for each q value=0.27869
EXPINT16Add an experimental value for each q value=0.228546
EXPINT17Add an experimental value for each q value=0.184443
EXPINT18Add an experimental value for each q value=0.144056
EXPINT19Add an experimental value for each q value=0.106039
EXPINT20Add an experimental value for each q value=0.0725692
EXPINT21Add an experimental value for each q value=0.0480033

... SAXS

The SAXS action with label test2 calculates the following quantities:
 Quantity   
 Description  
test2.score
the Metainference score
test2.sigma
uncertainty parameter
test2.sigmaMean
uncertainty in the mean estimate
test2.neff
effective number of replicas
test2.acceptSigma
MC acceptance for sigma values
test2.q
The # SAXS of qlsaxs: STATSCalculates statistical properties of a set of collective variables with respect to a set of reference values. More details ARGthe labels of the values from which the function is calculated=(test2\.q-.*) PARARGthe input for this action is the scalar output from one or more other actions without derivatives=(test2\.exp-.*)


The STATS action with label lsaxs calculates the following quantities:
 Quantity   
 Description  
lsaxs.sqdevsum
the sum of the squared deviations between arguments and parameters
lsaxs.corr
the correlation between arguments and parameters
lsaxs.slope
the slope of a linear fit between arguments and parameters
lsaxs.intercept
the intercept of a linear fit between arguments and parametersPBMETADUsed to performed Parallel Bias metadynamics. More details ...
LABELa label for the action so that its output can be referenced in the input to other actions=bias
ARGthe labels of the scalars on which the bias will act=d1,d2,d3,tor1,tor2,tor3,rg,linktor1,linktor2
HEIGHTthe height of the Gaussian hills, one for all biases=0.1
PACEthe frequency for hill addition, one for all biases=200
BIASFACTORuse well tempered metadynamics with this bias factor, one for all biases=30
ADAPTIVEuse a geometric (=GEOM) or diffusion (=DIFF) based hills width scheme=DIFF
GRID_WFILESdump grid for the bias, default names are used if GRID_WSTRIDE is used without GRID_WFILES=GRID.0,GRID.1,GRID.2,GRID.3,GRID.4,GRID.5,GRID.6,GRID.7,GRID.8
GRID_WSTRIDEfrequency for dumping the grid=2000
#GRID_RFILES=GRID.0,GRID.1,GRID.2,GRID.3,GRID.4,GRID.5,GRID.6,GRID.7,GRID.8
SIGMAthe widths of the Gaussian hills=10000
SIGMA_MINthe lower bounds for the sigmas (in CV units) when using adaptive hills=0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005
SIGMA_MAXthe upper bounds for the sigmas (in CV units) when using adaptive hills=0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2

GRID_MINthe lower bounds for the grid=0,0,0,-pi,-pi,-pi,0,0,0
GRID_MAXthe upper bounds for the grid=7,14,7,pi,pi,pi,7,6,6

WALKERS_MPI Switch on MPI version of multiple walkers - not compatible with WALKERS_* options other than WALKERS_DIR
UPDATE_FROMOnly update this action from this time=10000
... PBMETAD

The PBMETAD action with label bias calculates the following quantities:
 Quantity   
 Description  
bias.bias
the instantaneous value of the bias potentialMETAINFERENCECalculates the Metainference energy for a set of experimental data. More details ...
SCALEDATA Set to TRUE if you want to sample a scaling factor common to all values and replicas
AVERAGINGStride for calculation of averaged weights and sigma_mean=100
NOISETYPE functional form of the noise (GAUSS,MGAUSS,OUTLIERS,MOUTLIERS,GENERIC)=MOUTLIERS
ARGthe labels of the scalars on which the bias will act=(test2\.q-.*),bias.bias REWEIGHT simple REWEIGHT using the latest ARG as energy PARARGreference values for the experimental data, these can be provided as arguments without derivatives=(test2\.exp-.*)
SCALE_PRIOR either FLAT or GAUSSIAN=FLAT
ADDOFFSET Set to TRUE if you want to sample an offset common to all values and replicas
OFFSET_PRIOR either FLAT or GAUSSIAN=FLAT
OFFSET0 initial value of the offset=0 OFFSET_MINminimum value of the offset=-10 OFFSET_MAXmaximum value of the offset=10 DOFFSETmaximum MC move of the offset=0.005
SCALE0 initial value of the scaling factor=1.00 SCALE_MINminimum value of the scaling factor=0.8 SCALE_MAXmaximum value of the scaling factor=1.2 DSCALEmaximum MC move of the scaling factor=0.001
SIGMA0 initial value of the uncertainty parameter=0.01 SIGMA_MIN minimum value of the uncertainty parameter=0. SIGMA_MAX maximum value of the uncertainty parameter=0.05
OPTSIGMAMEAN Set to NONE/SEM to manually set sigma mean, or to estimate it on the fly=SEM
LABELa label for the action so that its output can be referenced in the input to other actions=bq
WRITE_STRIDE write the status to a file every N steps, this can be used for restart/continuation=10000
STRIDEthe frequency with which the forces due to the bias should be calculated=5
... METAINFERENCE

The METAINFERENCE action with label bq calculates the following quantities:
 Quantity   
 Description  
bq.bias
the instantaneous value of the bias potential
bq.sigma
uncertainty parameter
bq.sigmaMean
uncertainty in the mean estimate
bq.neff
effective number of replicas
bq.acceptSigma
MC acceptance for sigma values
bq.acceptScale
MC acceptance for scale value
bq.weight
weights of the weighted average
bq.biasDer
derivatives with respect to the bias
bq.scale
scale parameter
bq.offset
offset parameterensaxs: ENSEMBLECalculates the replica averaging of a collective variable over multiple replicas. More details ARGthe labels of the values from which the function is calculated=(test2\.q-.*),bias.bias REWEIGHT simple REWEIGHT using the latest ARG as energy
The ENSEMBLE action with label ensaxs calculates the following quantities:
 Quantity   
 Description  
ensaxs..#!custom
the names of the output components for this action depend on the actions input file see the example inputs below for detailsstatsq: STATSCalculates statistical properties of a set of collective variables with respect to a set of reference values. More details ARGthe labels of the values from which the function is calculated=(ensaxs\.test2\.q-.*) PARARGthe input for this action is the scalar output from one or more other actions without derivatives=(test2\.exp-.*)
The STATS action with label statsq calculates the following quantities:
 Quantity   
 Description  
statsq.sqdevsum
the sum of the squared deviations between arguments and parameters
statsq.corr
the correlation between arguments and parameters
statsq.slope
the slope of a linear fit between arguments and parameters
statsq.intercept
the intercept of a linear fit between arguments and parametersenANA: ENSEMBLECalculates the replica averaging of a collective variable over multiple replicas. More details ARGthe labels of the values from which the function is calculated=d1,d2,d3,tor1,tor2,tor3,rg,linktor1,linktor2,bias.bias REWEIGHT simple REWEIGHT using the latest ARG as energy

The ENSEMBLE action with label enANA calculates the following quantities:
 Quantity   
 Description  
enANA..#!custom
the names of the output components for this action depend on the actions input file see the example inputs below for detailsFLUSHThis command instructs plumed to flush all the open files with a user specified frequency. More details STRIDEthe frequency with which all the open files should be flushed=1000
PRINTPrint quantities to a file. More details ARGthe labels of the values that you would like to print to the file=(test2\.q-.*),(test2\.exp-.*) FILEthe name of the file on which to output these quantities=CV.dat STRIDE the frequency with which the quantities of interest should be output=100
PRINTPrint quantities to a file. More details FILEthe name of the file on which to output these quantities=RESTRAINTS ARGthe labels of the values that you would like to print to the file=bq.* STRIDE the frequency with which the quantities of interest should be output=100
PRINTPrint quantities to a file. More details FILEthe name of the file on which to output these quantities=q.dat ARGthe labels of the values that you would like to print to the file=statsq.*,lsaxs.* STRIDE the frequency with which the quantities of interest should be output=100
PRINTPrint quantities to a file. More details FILEthe name of the file on which to output these quantities=Analyse.dat ARGthe labels of the values that you would like to print to the file=d1,d2,d3,tor1,tor2,tor3,rg,linktor1,linktor2 STRIDE the frequency with which the quantities of interest should be output=100
PRINTPrint quantities to a file. More details FILEthe name of the file on which to output these quantities=ENSAXS ARGthe labels of the values that you would like to print to the file=(ensaxs\.test2\.q-.*) STRIDE the frequency with which the quantities of interest should be output=100
PRINTPrint quantities to a file. More details FILEthe name of the file on which to output these quantities=ENANA ARGthe labels of the values that you would like to print to the file=enANA.* STRIDE the frequency with which the quantities of interest should be output=100