Project ID: plumID:25.004
Source: H2_flooding/coverage_33.3/plumed.dat
Originally used with PLUMED version: 2.9
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
# vim: ft=plumedEnables syntax highlighting for PLUMED files in vim. See here for more details.
UNITSThis command sets the internal units for the code. More details LENGTHthe units of lengths=A TIMEthe units of time=0.001
The UNITS action with label calculates somethingHtotThe GROUP action with label Htot calculates the following quantities:| Quantity | Type | Description |
| Htot | atoms | indices of atoms specified in GROUP |
: GROUPDefine 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 ATOMSthe numerical indexes for the set of atoms in the group=25,26,27,28,37,38,39,40,49,50,51,52,61,62,63,64,72,73,74,75,108,109,110,111,119,120,121,122,131,132,133,134,143,144,145,146,155,156,157,158,191,192,193,194,203,204,205,206,215,216,217,218,227,228,229,230,238,239,240,241,273,274,275,276,285,286,287,288,297,298,299,300,309,310,311,312,320,321,322,323,356,357,358,359,368,369,370,371,380,381,382,383,392,393,394,395,404,405,406,407,440,441,442,443,452,453,454,455,464,465,466,467,476,477,478,479,488,489,490,491,524,525,526,527,536,537,538,539,548,549,550,551,560,561,562,563,572,573,574,575,608,609,610,611,619,620,621,622,631,632,633,634,642,643,644,645,654,655,656,657,690,691,692,693,702,703,704,705,714,715,716,717,725,726,727,728,737,738,739,740
HThe GROUP action with label H calculates the following quantities:| Quantity | Type | Description |
| H | atoms | indices of atoms specified in GROUP |
: GROUPDefine 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 ATOMSthe numerical indexes for the set of atoms in the group=749,752,755,756,754,747,748,753,750,751
aThe FIXEDATOM action with label a calculates the following quantities:| Quantity | Type | Description |
| a | atoms | virtual atom calculated by FIXEDATOM action |
: FIXEDATOMAdd a virtual atom in a fixed position. This action has hidden defaults. More details ATcoordinates of the virtual atom=0,0,47.5
a: FIXEDATOMAdd a virtual atom in a fixed position. This action uses the defaults shown here. More details ATcoordinates of the virtual atom=0,0,47.5 SET_MASS mass of the virtual atom=1 SET_CHARGE charge of the virtual atom=0
dHHThe DISTANCES action with label dHH calculates the following quantities:| Quantity | Type | Description |
| dHH | vector | the DISTANCES between the each pair of atoms that were specified |
| dHH_min | scalar | the minimum colvar |
: DISTANCESCalculate the distances between multiple piars of atoms This action is a shortcut. More details GROUPCalculate the distance between each distinct pair of atoms in the group=H MINcalculate the minimum value={BETA=20.0}
# dHH: DISTANCES GROUP=H MIN={BETA=20.0}
dHHThe DISTANCE action with label dHH calculates the following quantities:| Quantity | Type | Description |
| dHH | vector | the DISTANCE for each set of specified atoms |
: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMS1the pair of atom that we are calculating the distance between=752,749 ATOMS2the pair of atom that we are calculating the distance between=755,749 ATOMS3the pair of atom that we are calculating the distance between=755,752 ATOMS4the pair of atom that we are calculating the distance between=756,749 ATOMS5the pair of atom that we are calculating the distance between=756,752
dHH_me_minThe CUSTOM action with label dHH_me_min calculates the following quantities:| Quantity | Type | Description |
| dHH_me_min | vector | the vector obtained by doing an element-wise application of an arbitrary function to the input vectors |
: CUSTOMCalculate a combination of variables using a custom expression. More details ARGthe values input to this function=dHH FUNCthe function you wish to evaluate=exp(20.0/x) PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO
dHH_mec_minThe SUM action with label dHH_mec_min calculates the following quantities:| Quantity | Type | Description |
| dHH_mec_min | scalar | the sum of all the elements in the input vector |
: SUMCalculate the sum of the arguments More details ARGthe values input to this function=dHH_me_min PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO
dHH_minThe CUSTOM action with label dHH_min calculates the following quantities:| Quantity | Type | Description |
| dHH_min | scalar | an arbitrary function |
: CUSTOMCalculate a combination of variables using a custom expression. More details ARGthe values input to this function=dHH_mec_min FUNCthe function you wish to evaluate=20.0/log(x) PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO
# --- End of included input --- cHHThe COORDINATIONNUMBER action with label cHH calculates the following quantities:| Quantity | Type | Description |
| cHH | vector | the coordination numbers of the specified atoms |
| cHH_max | scalar | the maximum colvar |
: COORDINATIONNUMBERCalculate the coordination numbers of atoms so that you can then calculate functions of the distribution of This action is a shortcut and it has hidden defaults. More details SPECIESthe list of atoms for which the symmetry function is being calculated and the atoms that can be in the environments=H R_0The r_0 parameter of the switching function=1.5 NN The n parameter of the switching function =6 MM The m parameter of the switching function; 0 implies 2*NN=12 MAXcalculate the maximum value={BETA=0.05}
cHH: COORDINATIONNUMBERCalculate the coordination numbers of atoms so that you can then calculate functions of the distribution of This action is a shortcut and uses the defaults shown here. More details SPECIESthe list of atoms for which the symmetry function is being calculated and the atoms that can be in the environments=H R_0The r_0 parameter of the switching function=1.5 NN The n parameter of the switching function =6 MM The m parameter of the switching function; 0 implies 2*NN=12 MAXcalculate the maximum value={BETA=0.05} D_0 The d_0 parameter of the switching function=0.0
# cHH: COORDINATIONNUMBER SPECIES=H R_0=1.5 NN=6 MM=12 MAX={BETA=0.05}
cHH_grpThe GROUP action with label cHH_grp calculates the following quantities:| Quantity | Type | Description |
| cHH_grp | atoms | indices of atoms specified in GROUP |
: GROUPDefine 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 ATOMSthe numerical indexes for the set of atoms in the group=H
cHH_matThe CONTACT_MATRIX action with label cHH_mat calculates the following quantities:| Quantity | Type | Description |
| cHH_mat | matrix | a matrix containing the weights for the bonds between each pair of atoms |
: CONTACT_MATRIXAdjacency matrix in which two atoms are adjacent if they are within a certain cutoff. More details GROUPspecifies the list of atoms that should be assumed indistinguishable=H R_0The r_0 parameter of the switching function=1.5 D_0 The d_0 parameter of the switching function=0.0 NN The n parameter of the switching function =6 MM The m parameter of the switching function; 0 implies 2*NN=12
cHH_onesThe CONSTANT action with label cHH_ones calculates the following quantities:| Quantity | Type | Description |
| cHH_ones | vector | the constant value that was read from the plumed input |
: ONESCreate a constant vector with all elements equal to one More details SIZEthe number of ones that you would like to create=10
cHHThe MATRIX_VECTOR_PRODUCT action with label cHH calculates the following quantities:| Quantity | Type | Description |
| cHH | vector | the vector that is obtained by taking the product between the matrix and the vector that were input |
: MATRIX_VECTOR_PRODUCTCalculate the product of the matrix and the vector More details ARGthe label for the matrix and the vector/scalar that are being multiplied=cHH_mat,cHH_ones
cHH_caverageThe MEAN action with label cHH_caverage calculates the following quantities:| Quantity | Type | Description |
| cHH_caverage | scalar | the mean of all the elements in the input vector |
: MEANCalculate the arithmetic mean of the elements in a vector More details ARGthe values input to this function=cHH PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO
cHH_me_maxThe CUSTOM action with label cHH_me_max calculates the following quantities:| Quantity | Type | Description |
| cHH_me_max | vector | the vector obtained by doing an element-wise application of an arbitrary function to the input vectors |
: CUSTOMCalculate a combination of variables using a custom expression. More details ARGthe values input to this function=cHH FUNCthe function you wish to evaluate=exp(x/0.05) PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO
cHH_mec_maxThe SUM action with label cHH_mec_max calculates the following quantities:| Quantity | Type | Description |
| cHH_mec_max | scalar | the sum of all the elements in the input vector |
: SUMCalculate the sum of the arguments More details ARGthe values input to this function=cHH_me_max PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO
cHH_maxThe CUSTOM action with label cHH_max calculates the following quantities:| Quantity | Type | Description |
| cHH_max | scalar | an arbitrary function |
: CUSTOMCalculate a combination of variables using a custom expression. More details ARGthe values input to this function=cHH_mec_max FUNCthe function you wish to evaluate=0.05*log(x) PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO
# --- End of included input ---
dzThe ZDISTANCES action with label dz calculates the following quantities:| Quantity | Type | Description |
| dz_max | scalar | the maximum colvar |
: ZDISTANCESCalculate the z components of the vectors connecting one or many pairs of atoms. This action is a shortcut. More details GROUPACalculate the distances between all the atoms in GROUPA and all the atoms in GROUPB=a GROUPBCalculate the distances between all the atoms in GROUPA and all the atoms in GROUPB=H MAXcalculate the maximum value={BETA=0.05}
# dz: ZDISTANCES GROUPA=a GROUPB=H MAX={BETA=0.05}
dzThe DISTANCE action with label dz calculates the following quantities:| Quantity | Type | Description |
| dz.x | vector | the x-component of the vector connecting the two atoms |
| dz.y | vector | the y-component of the vector connecting the two atoms |
| dz.z | vector | the z-component of the vector connecting the two atoms |
: DISTANCECalculate the distance/s between pairs of atoms. More details COMPONENTS calculate the x, y and z components of the distance separately and store them as label ATOMS1the pair of atom that we are calculating the distance between=a,749 ATOMS2the pair of atom that we are calculating the distance between=a,752 ATOMS3the pair of atom that we are calculating the distance between=a,755 ATOMS4the pair of atom that we are calculating the distance between=a,756 ATOMS5the pair of atom that we are calculating the distance between=a,754
dz_me_maxThe CUSTOM action with label dz_me_max calculates the following quantities:| Quantity | Type | Description |
| dz_me_max | vector | the vector obtained by doing an element-wise application of an arbitrary function to the input vectors |
: CUSTOMCalculate a combination of variables using a custom expression. More details ARGthe values input to this function=dz.z FUNCthe function you wish to evaluate=exp(x/0.05) PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO
dz_mec_maxThe SUM action with label dz_mec_max calculates the following quantities:| Quantity | Type | Description |
| dz_mec_max | scalar | the sum of all the elements in the input vector |
: SUMCalculate the sum of the arguments More details ARGthe values input to this function=dz_me_max PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO
dz_maxThe CUSTOM action with label dz_max calculates the following quantities:| Quantity | Type | Description |
| dz_max | scalar | an arbitrary function |
: CUSTOMCalculate a combination of variables using a custom expression. More details ARGthe values input to this function=dz_mec_max FUNCthe function you wish to evaluate=0.05*log(x) PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO
# --- End of included input ---
target1The CUSTOM action with label target1 calculates the following quantities:| Quantity | Type | Description |
| target1 | scalar | an arbitrary function |
: CUSTOMCalculate a combination of variables using a custom expression. More details ARGthe values input to this function=cHH.max FUNCthe function you wish to evaluate=step(x-0.8) PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO
OPES_METADOn-the-fly probability enhanced sampling with metadynamics-like target distribution. This action has hidden defaults. More details ...
LABELa label for the action so that its output can be referenced in the input to other actions=opes1The OPES_METAD action with label opes1 calculates the following quantities:| Quantity | Type | Description |
| opes1.bias | scalar | the instantaneous value of the bias potential |
| opes1.rct | scalar | estimate of c(t). log(exp(beta V)/beta, should become flat as the simulation converges. Do NOT use for reweighting |
| opes1.zed | scalar | estimate of Z_n. should become flat once no new CV-space region is explored |
| opes1.neff | scalar | effective sample size |
| opes1.nker | scalar | total number of compressed kernels used to represent the bias |
ARGthe labels of the scalars on which the bias will act=cHH.max
FILE a file in which the list of all deposited kernels is stored=Kernels1.data
TEMP temperature=700
PACEthe frequency for kernel deposition=250
BARRIERthe free energy barrier to be overcome=60
RESTARTallows per-action setting of restart (YES/NO/AUTO)=NO
EXCLUDED_REGIONkernels are not deposited when the action provided here has a nonzero value, see example above=target1
... OPES_METAD
OPES_METADOn-the-fly probability enhanced sampling with metadynamics-like target distribution. This action uses the defaults shown here. More details ...
LABELa label for the action so that its output can be referenced in the input to other actions=opes1
ARGthe labels of the scalars on which the bias will act=cHH.max
FILE a file in which the list of all deposited kernels is stored=Kernels1.data
TEMP temperature=700
PACEthe frequency for kernel deposition=250
BARRIERthe free energy barrier to be overcome=60
RESTARTallows per-action setting of restart (YES/NO/AUTO)=NO
EXCLUDED_REGIONkernels are not deposited when the action provided here has a nonzero value, see example above=target1
SIGMA the initial widths of the kernels=ADAPTIVE COMPRESSION_THRESHOLD merge kernels if closer than this threshold, in units of sigma=1
... OPES_METAD
COMMITTORDoes a committor analysis. More details ...
ARGthe labels of the values which is being used to define the committor surface=dz.max
STRIDE the frequency with which the CVs are analyzed=1000
BASIN_LL1List of lower limits for basin #=6.0
BASIN_UL1List of upper limits for basin #=15.0
FILEthe name of the file on which to output the reached basin=FORMATION
... COMMITTOR
PRINTPrint quantities to a file. More details ARGthe labels of the values that you would like to print to the file=* STRIDE the frequency with which the quantities of interest should be output=10 FILEthe name of the file on which to output these quantities=COLVAR
FLUSHThis 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=1
