计算化学公社

标题: ENCUT设置问题讨论 [打印本页]

作者
Author:
Stars    时间: 2023-12-20 00:01
标题: ENCUT设置问题讨论
在Gaussian中计算单点时会用高一个档次的标准计算能量,而在优化结构时用较低的标准。那么VASP有没有这个说法?因为平面波的“基组”大小是由截断能ENCUT决定的,但结构优化时用较小的ENCUT而能量计算时用较大的ENCUT的做法貌似很少采用(具体是文献中经常写ENCUT was set to ***eV,似乎是整个计算中用相同的ENCUT),这是什么原因?

作者
Author:
wsz    时间: 2023-12-20 00:44
对ENCUT做过收敛性测试的话,就可以一直用这个ENCUT,再增加意义不大。
作者
Author:
Stars    时间: 2023-12-20 09:29
本帖最后由 Stars 于 2023-12-20 09:31 编辑
wsz 发表于 2023-12-20 00:44
对ENCUT做过收敛性测试的话,就可以一直用这个ENCUT,再增加意义不大。

所以说,在VASP里面可以通过一个合理的ENCUT截断达到能量收敛,因此可以从头到尾用相同的ENCUT值?但在Gaussian中为减轻计算压力,可以用小一点的基组完成优化?如果按照VASP的思路,平常使用Gaussian时也应当做基组收敛性检查?或者说VASP在进行比如结构优化计算时可以用粗糙的截断标准?
作者
Author:
Stars    时间: 2023-12-20 09:34
另外像这种问题,回复时同时涉及了量子化学与第一性原理的内容的应该如何分区?如有不当,麻烦管理员了!
作者
Author:
ahxb    时间: 2023-12-20 10:09
主要是少有人做这样的benchmark,vasp算个Hartree-Fock计算量得加两三个数量级,更不用说post-HF,缺乏用作标准的数据。
作者
Author:
白皛白    时间: 2023-12-20 10:53
VASP里提高精度貌似更多的是通过使用更严格的SCF收敛限与更多的K点,调整ENCUT没有怎么听说过。
贴一段vaspwiki中ENCUT条目的话
We strongly recommend specifying the energy cutoff ENCUT always manually in the INCAR file to ensure the same accuracy between calculations. Otherwise, the default ENCUT may differ among the different calculations (e.g., for the calculation of the cohesive energy), with the consequence that the total energies, for instance, can not be compared.
似乎也建议始终使用一致的ENCUT进行计算。但是不知道比如结构优化时使用ENCUT=300,单点时使用ENCUT=500,然后进行能量比较时只在单点计算得到的能量之间进行比较这样是否可行。
作者
Author:
Stars    时间: 2023-12-20 11:04
白皛白 发表于 2023-12-20 10:53
VASP里提高精度貌似更多的是通过使用更严格的SCF收敛限与更多的K点,调整ENCUT没有怎么听说过。
贴一段vas ...

我好像意识到问题了。就是在VASP计算时k点对计算时间影响更大,截断不会带来这么大的时间影响,所以就都统一用一样的截断了
作者
Author:
Weldingspock    时间: 2023-12-20 15:05
VASP里面需要测试ENCUT对体系能量的收敛影响,此外,对于含有一些元素的体系(例如C,N),不够大的ENCUT会引入Pulay stress,对体系的计算造成影响。https://www.vasp.at/wiki/index.p ... ns_and_Pulay_stress
一般建议设置ENCUT在ENMAX的1.3倍以上,设置1.5倍基本就是一个稳妥的选择了。但是很可惜这个标准在杂化钙钛矿里面并没有很好的执行,大多数计算杂化钙钛矿的文章ENCUT设置在450eV-500 eV.很难保证能量的收敛
作者
Author:
snljty2    时间: 2023-12-21 10:33
原理上不同任务应该做不同的截断能测试的,比如几何优化任务实际是要求随着截断能增加,几何结构的变化基本可以忽略不计。比如反应势垒的要求是随着截断能增加,计算的反应势垒几乎不变,等等。实际机时有限,一般经验性测一下单点能,选一个估计基本足够的截断能也是比较常见的做法。
作者
Author:
snljty2    时间: 2023-12-21 10:33
本帖最后由 snljty2 于 2023-12-21 10:42 编辑

[手抖发重了,麻烦管理员删除,谢谢]




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3