good news:
with the today commit, (finally !!), there's the possibility to run
PM6 semi-empirical calculation in
CP2K (fully support for all elements) (at the moment only numerical
derivatives.. analytical will come soon).
bad news:
In the original version of PM6 there are some terms that make the
heat of formation a discontinuous function of
the nuclei positions (are just correction terms, constants that
become zero for different topological configurations) that
we didn't implement. Therefore for few situations there could be
discrepancies between the heat of formation of CP2K
and the one of MOPAC (but the electronic energy should anyway be the
same (within numerical errors (see below)).
An example is H-CC-H, that has a constant correction term to the heat
of formation of the triple bond that becomes
zero when the triple bond is broken (it's just a step function).
Why have been these terms omitted?
Simply because there's no trace of them in the PM6 paper.
I'm sure that there could be other hidden corrections/parameters not
explicitly mentioned in the paper (classical correction
like the one for correcting the sp2 pyramidalization, that we didn't
implement as well), so in case of big discrepancies
compared with MOPAC2007 feel free to post here your input file and a
short description of the problem.
What does it mean big discrepancies?
These are the values for the heat of formation (kcalmol) of the
following systems (CP2K/MOPAC2007)
HCl (17.0356/17.0345)
BrCl (2113.534/2113.493)
TiO (202.714/202.633)
If you have some discrepancy larger than the ones above, check always
the value of the electronic energy that
must anyway be VERY similar (the only difference is in the way we
compute integrals.. so just kind of numerical
errors, totally negligible).
Have fun!!
Teo
The analytical derivatives for PM6 and other SE methods using d-
orbitals have been implemented and are
available in the CVS.
Cheers,
Teo
--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+uns...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/e93c6e15-39d4-4b8e-865f-f857da3f9f1b%40googlegroups.com.
To unsubscribe from this group and stop receiving emails from it, send an email to cp...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/e93c6e15-39d4-4b8e-865f-f857da3f9f1b%40googlegroups.com.
is it possible to run
QM / "QM"
eg. DFT/SE like PBE / PM6
via FORCE_MIXING
in cp2k,rather than QM/"MM"?
Best regards
Bernd
&MULTIPLE_FORCE_EVALS FORCE_EVAL_ORDER 2 3 &END &FORCE_EVAL METHOD MIXED &MIXED MIXING_TYPE GENMIX &GENERIC MIXING_FUNCTION X+Y-Z VARIABLES X Y Z &END GENERIC &END &END FORCE_EVAL &FORCE_EVAL METHOD FIST &END FORCE_EVAL &FORCE_EVAL METHOD QS &END FORCE_EVAL &FORCE_EVAL METHOD FIST &END FORCE_EVAL
So the way to set up
a (DFT)PBE-PM6 scheme would be
&MULTIPLE_FORCE_EVALS FORCE_EVAL_ORDER 1 2 &END
&FORCE_EVAL METHOD QS
PM6
....
&END FORCE_EVAL
&FORCE_EVAL METHOD QS
PBE
.....
&END FORCE_EVAL
Would this be correct?
Many thanks
Bernd