Metadynamics variables,
Lucas Lodeiro
eluni... at gmail.com
Wed Aug 22 22:02:46 UTC 2018
Hi all!
I want to run a MTD for a water system, and i want to select the values of
some variables as i want: Gaussian height 0.0005 Ha (0.3 kcal/mol
approximately), width 0.1 bohr and deposition time step of 40 fs.
The colvar definition between the atoms uses a Rc=1.5 Angstrom. The
timestep MTD is 0.5 fs.
I prepared an input file:
&FREE_ENERGY
METHOD METADYN
&METADYN
DO_HILLS
NT_HILLS 80
WW [kcalmol] 0.3
&METAVAR
COLVAR 1
SCALE 0.1
&END METAVAR
&PRINT
&COLVAR
COMMON_ITERATION_LEVELS 3
&EACH
MD 1
&END
&END
&HILLS
COMMON_ITERATION_LEVELS 3
&EACH
MD 1
&END
&END
&END
&END METADYN
&END FREE_ENERGY
&SUBSYS
&COLVAR
&COORDINATION
ATOMS_FROM 1
ATOMS_TO 2
R0 [angstrom] 1.5
ND 12
NN 6
&END COORDINATION
&END COLVAR
&END SUBSYS
The WW is the heicht of the gaussian in kcal/mol, and NT_HILL 80 is the
number of steps to enter Gaussians, 80*0.5 fs = 40 fs, all ok. But my doubt
is, SCALE is the width of the gaussian, in the manual the default unity
say, [internal_cp2k], and i dont understand if i put 0.1 in that variable
like as in the input, i use a 0.1 bohr width...
How can i use safely 0.1 bohr in the width of gaussian?
Regards
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20180822/97b84637/attachment.htm>
More information about the CP2K-user
mailing list