Cp2k convergence test: is it possible that cutoff of 250 Ry has fairly lower total energy than 350 Ry (with same rel_cutoff 60 Ry)?

anyu... at gmail.com anyu... at gmail.com
Fri Aug 24 14:29:24 UTC 2018


Dear CP2K developers and users,

I did the convergence tests with cp2k 3.0 for hBN system, following the 
instructions of https://www.cp2k.org/howto:converging_cutoff . The energy 
values and the NG on grids are as follows. I wonder why higher cutoff 
values give fairly higher energies, for example, the total energy for 
cutoff 350 ry is higher than 300 ry, I thought in principle, higher cutoff 
values could have more stable energy (lower value), at least not that 
higher than the number with lower cutoffs. But here in my case, it's not, 
can anybody help me to see why?  Is 250Ry and rel_cutoff of 70 Ry here a 
good choice? 

Also, when one starts a calculation,  should the convergence tests always 
be taken with different values of cut off and rel_cutoff? or it's safe to 
just use larger values such as (cutoff 400 ry, rel_cutoff 80ry)?  

Thanks a lot in advance for all your suggestions.

Best regards,
Yun

# Grid cutoff vs total energy
# Date: Fri Aug 24 13:38:48 CEST 2018
# PWD: /home/Convergence_test
# REL_CUTOFF = 60
# Cutoff (Ry) | Total Energy (Ha) | NG on grid 1 | NG on grid 2 | NG on 
grid 3 | NG on grid 4
     50.00  -833.3666229006     559488   81536       0       0
    100.00  -824.4110369450     457024  180160    3840       0
    150.00  -824.3925595993     424576  134912   81536       0
    200.00  -824.3899074807     345088  210560   81536    3840
    250.00  -824.3979520649     294272  162752  180160    3840
    300.00  -824.3920162797   294272  162752  180160    3840
    350.00  -824.3858376968    271744  185280  119488   64512
    400.00  -824.3841314539    271744  152832  134912   81536
    450.00  -824.3831328479     271744  152832  134912   81536
    500.00  -824.3833590835     236416  184320  138752   81536


# Rel Grid cutoff vs total energy
# Date: Fri Aug 24 14:54:23 CEST 2018
# PWD: /home/Convergence_test/Rel_cutoff
# CUTOFF = 250
# Rel Cutoff (Ry) | Total Energy (Ha) | NG on grid 1 | NG on grid 2 | NG on 
grid 3 | NG on grid 4
     10.00  -824.9353843726     24384  212032  184320  220288
     20.00  -824.4674485372     117312  176960  162752  184000
     30.00  -824.3980791451    236416  184320  138752   81536
     40.00  -824.3980825574     271744  152832  134912   81536
     50.00  -824.3979520649     294272  162752  180160    3840
     60.00  -824.3979520649     294272  162752  180160    3840
     70.00  -824.3979509738     345088  210560   81536    3840
     80.00  -824.3979509738     345088  210560   81536    3840
     90.00  -824.3979509738     420736  138752   79616    1920
    100.00  -824.3979509738    424576  134912   81536       0


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20180824/ec8be3db/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: convergence_test.inp
Type: chemical/x-gamess-input
Size: 9336 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20180824/ec8be3db/attachment.inp>


More information about the CP2K-user mailing list