Convergence problem wrt cutoff
Valérie
valerie.... at gmail.com
Tue Mar 24 16:40:25 UTC 2015
Hi all,
I would like to run a geometry optimization of a single molecule with
Quickstep but got stuck at the first step; i.e. get meaningful values of
cutoff and relative cutoffs.
I followed this tutorial: http://www.cp2k.org/howto:converging_cutoff but
the total energy does not converge to great accuracy with my system:
# Grid cutoff vs total energy
# REL_CUTOFF = 100
# Cutoff (Ry) | Total Energy (Ha) | NG on grid 1 | NG on grid 2 | NG on
grid 3 | NG on grid 4
50.00 -68.6591152269 51340 490 0 0
100.00 -71.4960226714 49235 2499 96 0
150.00 -71.3943596348 46461 4879 490 0
200.00 -71.4042102566 44831 6209 790 0
250.00 -71.4043241216 42771 7117 1846 96
300.00 -71.4043468150 41166 8069 2499 96
350.00 -71.4042673160 38450 10341 2847 192
400.00 -71.4043263577 36970 10901 3671 288
450.00 -71.4043454189 35016 11445 4879 490
500.00 -71.4043448309 34404 11749 5187 490
550.00 -71.4043354081 33007 11952 6081 790
600.00 -71.4043430530 32450 12381 6209 790
I was wondering if there was something I could try to fix this or should I
carry on with Cutoff = 400?
Many thanks
Valerie
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20150324/bb4b1217/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Cytosine_E_CRC.inp
Type: chemical/x-gamess-input
Size: 1438 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20150324/bb4b1217/attachment.inp>
More information about the CP2K-user
mailing list