Project ID: plumID:25.005
Source: dimer_octahedral/plumed.dat
Originally used with PLUMED version: 2.8.2
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.10
tested onmaster
#RESTART
UNITSThis command sets the internal units for the code. More details LENGTHthe units of lengths=A ENERGYthe units of energy=kcal/mol TIMEthe units of time=fs
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

# define groups of different atoms
AG: 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=1-3,12-14
S: 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=4-6,15-17
C: 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=10,11,21-24
N: 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=7-9,18-20
H: 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-36
O_S: 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=41-726:6

# define center of mass of Ag atoms
c: COMCalculate the center of mass for a group of atoms. More details ATOMSthe list of atoms which are involved the virtual atom's definition=1-3,12-14

# define fixed atom near center of the box
c_box: FIXEDATOMAdd a virtual atom in a fixed position. More details ATcoordinates of the virtual atom=10.63933,9.75917,9.83800

# distance between COM of Ag atoms and center of the box
d: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c,c_box

# define coordination numbers for Ag and S atoms
cn_s1: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=16 GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=14 R_0The r_0 parameter of the switching function=3.25 NN The n parameter of the switching function =6 MM The m parameter of the switching function; 0 implies 2*NN=8 
cn_s2: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=4  GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=1  R_0The r_0 parameter of the switching function=3.25 NN The n parameter of the switching function =6 MM The m parameter of the switching function; 0 implies 2*NN=8
cn_s3: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=6  GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=2  R_0The r_0 parameter of the switching function=3.25 NN The n parameter of the switching function =6 MM The m parameter of the switching function; 0 implies 2*NN=8
cn_s4: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=17 GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=3  R_0The r_0 parameter of the switching function=3.25 NN The n parameter of the switching function =6 MM The m parameter of the switching function; 0 implies 2*NN=8
cn_s5: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=5  GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=12 R_0The r_0 parameter of the switching function=3.25 NN The n parameter of the switching function =6 MM The m parameter of the switching function; 0 implies 2*NN=8
cn_s6: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=15 GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=13 R_0The r_0 parameter of the switching function=3.25 NN The n parameter of the switching function =6 MM The m parameter of the switching function; 0 implies 2*NN=8
cn_ag_box: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=AG GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=c_box R_0The r_0 parameter of the switching function=3.0 NN The n parameter of the switching function =6 MM The m parameter of the switching function; 0 implies 2*NN=8

# CV2 solvent coordination around center of the box
cn: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=c_box GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=O_S R_0The r_0 parameter of the switching function=3.0 NN The n parameter of the switching function =6 MM The m parameter of the switching function; 0 implies 2*NN=10

# CV1 solute coordination
s1: COMBINECalculate a polynomial combination of a set of other variables. More details ARGthe values input to this function=cn_ag_box,cn_s1,cn_s2,cn_s3,cn_s4,cn_s5,cn_s6 COEFFICIENTS the coefficients of the arguments in your function=1,1,1,1,1,1,1 PERIODICif the output of your function is periodic then you should specify the periodicity of the function=NO

# WT-Metad block
metad: METADUsed to performed metadynamics on one or more collective variables. More details ARGthe labels of the scalars on which the bias will act=s1,cn PACEthe frequency for hill addition=100 FILE a file in which the list of added hills is stored=HILLS TEMPthe system temperature - this is only needed if you are doing well-tempered metadynamics=300 BIASFACTORuse well tempered metadynamics and use this bias factor=15 SIGMAthe widths of the Gaussian hills=0.1,0.025 HEIGHTthe heights of the Gaussian hills=2

# define lowest coordination of different atoms to monitor ML potential stability
cn_H: COORDINATIONNUMBERCalculate the coordination numbers of atoms so that you can then calculate functions of the distribution of More details SPECIESAthe list of atoms for which the symmetry function is being calculated=25-36,42-726:6,38-726:6,39-726:6,40-726:6 SPECIESBthe list of atoms that can be in the environments of each of the atoms for which the symmetry function is being calculated=1-726 R_0The r_0 parameter of the switching function=1.2 LOWEST this flag allows you to recover the lowest of these variables LOWMEM this flag does nothing and is present only to ensure back-compatibility
cn_CH3: COORDINATIONNUMBERCalculate the coordination numbers of atoms so that you can then calculate functions of the distribution of More details SPECIESAthe list of atoms for which the symmetry function is being calculated=37-726:6 SPECIESBthe list of atoms that can be in the environments of each of the atoms for which the symmetry function is being calculated=42-726:6,38-726:6,39-726:6,40-726:6 R_0The r_0 parameter of the switching function=1.5 LOWEST this flag allows you to recover the lowest of these variables LOWMEM this flag does nothing and is present only to ensure back-compatibility
cn_CO: COORDINATIONNUMBERCalculate the coordination numbers of atoms so that you can then calculate functions of the distribution of More details SPECIESAthe list of atoms for which the symmetry function is being calculated=41-726:6 SPECIESBthe list of atoms that can be in the environments of each of the atoms for which the symmetry function is being calculated=37-726:6 R_0The r_0 parameter of the switching function=1.5 LOWEST this flag allows you to recover the lowest of these variables LOWMEM this flag does nothing and is present only to ensure back-compatibility

# define minimum distance between all hydrogen atoms Hgrp: 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-36,42-726:6,38-726:6,39-726:6,40-726:6 dHH: DISTANCESCalculate the distances between multiple piars of atoms More details GROUPACalculate the distances between all the atoms in GROUPA and all the atoms in GROUPB=Hgrp GROUPBCalculate the distances between all the atoms in GROUPA and all the atoms in GROUPB=Hgrp MINcalculate the minimum value={BETA=20} NOPBC ignore the periodic boundary conditions when calculating distances
# define wall for S and non bonding Ag S1: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=16 GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=2,3,12,13 R_0The r_0 parameter of the switching function=2.75 S2: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=4 GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=3,12,13,14 R_0The r_0 parameter of the switching function=2.75 S3: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=6 GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=1,3,12,13 R_0The r_0 parameter of the switching function=2.75 S4: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=17 GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=1,2,12,14 R_0The r_0 parameter of the switching function=2.75 S5: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=5 GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=1,2,13,14 R_0The r_0 parameter of the switching function=2.75 S6: COORDINATIONCalculate coordination numbers. More details GROUPAFirst list of atoms=15 GROUPBSecond list of atoms (if empty, N*(N-1)/2 pairs in GROUPA are counted)=1,2,3,14 R_0The r_0 parameter of the switching function=2.75 # define walls to keep system in the region of interest a1: ANGLECalculate one or multiple angle/s. More details ATOMSthe list of atoms involved in this collective variable (either 3 or 4 atoms)=1,13,727 a2: ANGLECalculate one or multiple angle/s. More details ATOMSthe list of atoms involved in this collective variable (either 3 or 4 atoms)=12,14,727 a3: ANGLECalculate one or multiple angle/s. More details ATOMSthe list of atoms involved in this collective variable (either 3 or 4 atoms)=2,3,727 a4: ANGLECalculate one or multiple angle/s. More details ATOMSthe list of atoms involved in this collective variable (either 3 or 4 atoms)=1,727,12 a5: ANGLECalculate one or multiple angle/s. More details ATOMSthe list of atoms involved in this collective variable (either 3 or 4 atoms)=2,727,13 a6: ANGLECalculate one or multiple angle/s. More details ATOMSthe list of atoms involved in this collective variable (either 3 or 4 atoms)=3,727,14 d1: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,1 d2: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,2 d3: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,3 d4: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,12 d5: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,13 d6: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,14 d1cs: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,4 d2cs: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,5 d3cs: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,6 d4cs: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,15 d5cs: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,16 d6cs: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,17 d1cn: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,7 d2cn: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,8 d3cn: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,9 d4cn: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,18 d5cn: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,19 d6cn: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=c_box,20 d1n: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=20,1 d2n: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=7,2 d3n: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=8,3 d4n: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=9,12 d5n: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=18,13 d6n: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=19,14 d1s: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=16,1 d2s: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=4,2 d3s: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=5,3 d4s: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=15,12 d5s: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=17,13 d6s: DISTANCECalculate the distance/s between pairs of atoms. More details ATOMSthe pair of atom that we are calculating the distance between=6,14 uwall_s_ag: UPPER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=S1,S2,S3,S4,S5,S6 ATthe positions of the wall=0.5,0.5,0.5,0.5,0.5,0.5 KAPPAthe force constant for the wall=100,100,100,100,100,100 uwall_a: UPPER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=a1,a2,a3 ATthe positions of the wall=1.5,1.5,1.5 KAPPAthe force constant for the wall=100,100,100 lwall_a: LOWER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=a1,a2,a3 ATthe positions of the wall=0.65,0.65,0.65 KAPPAthe force constant for the wall=100,100,100 lwall_a2: LOWER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=a4,a5,a6 ATthe positions of the wall=2.4,2.4,2.4 KAPPAthe force constant for the wall=100,100,100 lwall: LOWER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=d1,d2,d3,d4,d5,d6 ATthe positions of the wall=1.5,1.5,1.5,1.5,1.5,1.5 KAPPAthe force constant for the wall=100,100,100,100,100,100 uwall: UPPER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=d1,d2,d3,d4,d5,d6 ATthe positions of the wall=5,5,5,5,5,5 KAPPAthe force constant for the wall=100,100,100,100,100,100 lwallcn: LOWER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=d1cn,d2cn,d3cn,d4cn,d5cn,d6cn ATthe positions of the wall=3,3,3,3,3,3 KAPPAthe force constant for the wall=100,100,100,100,100,100 lwallcs: LOWER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=d1cs,d2cs,d3cs,d4cs,d5cs,d6cs ATthe positions of the wall=3,3,3,3,3,3 KAPPAthe force constant for the wall=100,100,100,100,100,100 uwalln: UPPER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=d1n,d2n,d3n,d4n,d5n,d6n ATthe positions of the wall=2.7,2.7,2.7,2.7,2.7,2.7 KAPPAthe force constant for the wall=100,100,100,100,100,100 uwalls: UPPER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=d1s,d2s,d3s,d4s,d5s,d6s ATthe positions of the wall=2.7,2.7,2.7,2.7,2.7,2.7 KAPPAthe force constant for the wall=100,100,100,100,100,100 uwall_cn: UPPER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=cn ATthe positions of the wall=5.0 KAPPAthe force constant for the wall=100 cn_wall: LOWER_WALLSDefines a wall for the value of one or more collective variables, More details ARGthe arguments on which the bias is acting=cn_CH3.lowest,cn_CO.lowest ATthe positions of the wall=2.55,0.47 KAPPAthe force constant for the wall=100,100 res: RESTRAINTAdds harmonic and/or linear restraints on one or more variables. More details ARGthe values the harmonic restraint acts upon=d ATthe position of the restraint=0.0 KAPPA specifies that the restraint is harmonic and what the values of the force constants on each of the variables are=100 # define print statement PRINTPrint quantities to a file. More details STRIDE the frequency with which the quantities of interest should be output=100 FILEthe name of the file on which to output these quantities=colvar ARGthe labels of the values that you would like to print to the file=*