计算化学公社

标题: 请教CP2K截断设置问题 [打印本页]

作者
Author: 一条君 时间: 2021-8-23 22:27
标题: 请教CP2K截断设置问题
本帖最后由一条君于 2021-8-24 15:27 编辑

MD体系包含N、O、F元素，报错类别pw/pw_grids.F:1612。从最开始cut off=400改到560，rel_cut off=40，仍然报错，(1)请问还是截断设置的有问题吗？（参考https://www.cp2k.org/faq:cutoff，取所含元素中最大的指数项F元素的14，*40=560），还是应该设置成
a “对于包含Na, N, O, F, Ni, Ga等元素的计算需要设置到1000 Ry确保计算精度”呢 ;  b 如果设置成1000，64核是不是也跑不动了
还是另有原因呢？
(2)对于下面grid  1>4，是否也反应出应该继续调大截断？
(3)cp2k报错类别有官方说明吗，谢谢~

-------------------------------------------------------------------------------
----                            MULTIGRID INFO                         ----
-------------------------------------------------------------------------------
count for grid       1:       3651075       cutoff [a.u.]       260.00
count for grid       2:       2784864       cutoff [a.u.]          86.67
count for grid       3:       2520111       cutoff [a.u.]          28.89
count for grid       4:       1482248       cutoff [a.u.]          9.63
total gridlevel count  :    10438298

CPASSERT failed
pw/pw_grids.F:1612

9 pw_grid_assign
8 pw_grid_setup_internal
7 pw_grid_setup
6 pw_env_rebuild
5 qs_env_rebuild_pw_env
4 qs_forces
3 velocity_verlet
2 qs_mol_dyn_low
1 CP2K

作者
Author: sobereva 时间: 2021-8-24 11:29
那是平面波的动能截断，完全不是“截断半径”，标题有严重误导性
截断值取多大是最便宜的够用的值，可以做收敛性测试得到

贴完整报错信息，或者上传输出文件便于别人判断

可以把报错的关键提示放Google上搜，没有官网的全面的Q&A，不过官方的Google group里对很多常见的问题都有人回答过怎么解决

作者
Author: 一条君 时间: 2021-8-24 15:31
本帖最后由一条君于 2021-8-24 15:58 编辑

sobereva 发表于 2021-8-24 11:29
那是平面波的动能截断，完全不是“截断半径”，标题有严重误导性
截断值取多大是最便宜的够用的值，可以做 ...

谢谢老师，已修改标题，并上传了输出文件，查阅group后发现可能是内存问题，由于锡元素只有一个，尝试不了减小“晶胞”，还请老师再指教一下，谢谢（大概需要多大内存呢；这个量级波函数，换替他理论方法可以吗）

作者
Author: sobereva 时间: 2021-8-25 09:09
别用那么老CP2K版本

欢迎光临计算化学公社 (http://bbs.keinsci.com/)