units in GENPOT section
Noam Bernstein
no... at bollweevil.gdbg.org
Wed Jun 2 15:23:39 UTC 2010
Hi all - I'm trying to set up a DFT calculation with vdW corrections,
in the form of a non-bonded potential. I've been looking at the
example in tests/QS/regtest-dft-vdw-corr, where I found the
following section.
&GENPOT
ATOMS C C
RCUT [angstrom] 14.9
RMIN [angstrom] 0.2
RMAX [angstrom] 15.0
VARIABLES R
PARAMETERS CSIX RZERO D
VALUES -216.9850 6.57936 2.9688
FUNCTION (CSIX/(R**6))*((1.0 - EXP(-1.0*D*((R/RZERO)**7)))**4)
&END
I think I understand what it's doing, and the functional form is consistent
with what's in, for example, Williams and Malhotra (Chem Phys v 327).
However, I haven't been able to figure out what are the units. The
PARAMETERS must implicit have units (CSIX must be in units of
energy*length^6, e.g.), but the documentation says that PARAMETERS are
in cp2k_internal units, and I couldn't find those defined on the units page of
the documentation. Is my guess that they're Hartree and Bohr correct?
thanks,
Noam
P.S. Is there another good source of parameters that I should know about,
in addition to that Williams and Malhotra paper?
More information about the CP2K-user
mailing list