计算化学公社

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

作者
Author: XYHhhh 时间: 2022-10-22 16:07
标题: 请教VASP截断能设置问题
vasp.at官网有提到,VASP的默认截断能体现在POTCAR中的ENMAX，请问若有多个元素则默认截断能是多个元素的ENMAX的和还是取最大值？若要手动设置截断能，是要从结构优化开始到后面单点能计算以及非自洽计算都设置相同的截断能吗？另外单点能计算是否就是结构优化后的静态计算？小白求各位老师指导

作者
Author: 乐平 时间: 2022-10-22 16:46
本帖最后由乐平于 2022-10-22 16:58 编辑

1）多个元素的情况，程序默认取最大值。

2）最好手动设置，尤其是做晶胞优化（ISIF=3 时），ENCUT 还要取更大一些。网上教程的经验是 1.2 ~ 1.5 倍 ENMAX。当然，最好是做 ENCUT 的收敛测试，即用不同的 ENCUT 先算一组单点能计算，看看是否趋近于某个值，然后取此时对应的 ENCUT。后续计算过程中保持不变。
** 晶胞优化或者结构优化过程中，可以不保存 WAVECAR 和 CHGCAR，也就是 LWAVE = F, LCHARG = F

3）晶胞优化之后建议再做一次单点能计算（也就是静态自洽计算），得到 CHGCAR 和 WAVECAR 文件。后续的计算可以读取这两个文件，有助于加速电子步收敛。

** 所谓 “静态” 就是保持原子核不动，只做电子的自洽场计算。这就是课本上说的 “玻恩-奥本海默近似”。当然，结构优化过程中也用到了 “玻恩-奥本海默近似”，这是每完成一次电子自洽场计算后，需要根据当前固定的原子核势函数梯度的负值（即受力方向）调整原子核之间的距离，然后在新的结构下再进行电子自洽场计算。如此不断循环，直到体系的前后两次计算的能量之差小于程序设定的阈值，也就是 EDIFFG 的值，循环结束。

我看到你的输入文件是 VASPKIT 生成的，非常好。注意发文章时引用 VASPKIT 的文献 ^_^

作者
Author: XYHhhh 时间: 2022-10-22 18:43

乐平发表于 2022-10-22 16:46
1）多个元素的情况，程序默认取最大值。

2）最好手动设置，尤其是做晶胞优化（ISIF=3 时），ENCUT 还要 ...

谢谢老师的耐心解答

作者
Author: @oneone 时间: 2024-7-30 17:04
多个元素取ENMAX最大值，一般是1~1.5倍的ENMAX

作者
Author: 芝麻233 时间: 2025-9-10 16:06

乐平发表于 2022-10-22 16:46
1）多个元素的情况，程序默认取最大值。

2）最好手动设置，尤其是做晶胞优化（ISIF=3 时），ENCUT 还要 ...

老师好，请问如果结构优化和静态自洽时使用了不同的截断能，结构优化时400，静态自洽时600，会产生什么影响吗？

作者
Author: 北大-陶豫 时间: 2025-9-10 20:28

芝麻233 发表于 2025-9-10 16:06
老师好，请问如果结构优化和静态自洽时使用了不同的截断能，结构优化时400，静态自洽时600，会产生什么影 ...

不同 ENCUT 的计算，CHGCAR 无法互相读取。除此之外应该没啥影响了

作者
Author: 芝麻233 时间: 2025-9-12 11:39

北大-陶豫发表于 2025-9-10 20:28
不同 ENCUT 的计算，CHGCAR 无法互相读取。除此之外应该没啥影响了

谢谢老师的解答

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