计算化学公社

标题: 截断能测试能量走势问题 [打印本页]

作者
Author: 小书童 时间: 2017-2-20 10:11
标题: 截断能测试能量走势问题
本帖最后由小书童于 2017-2-20 11:19 编辑

截断能测试的能量曲线并不是逐渐降低并趋于平稳的，说明测试过程有问题，附上测试结果和测试脚本，麻烦老师们给指点一下，设置的脚本哪里不合理。ps：测试的体系是4*4的石墨烯，Kpoints中G中心，9*9*1的kmesh。

作者
Author: 卡开发发 时间: 2017-2-20 15:01
本帖最后由卡开发发于 2017-2-20 15:04 编辑

不一定会降低，只要平稳了就可以了。4x4石墨烯应该是32原子，500eV的截断能已经足以让能量收敛到1meV/atom左右了，如果再保守一点可以把力的收敛也考虑进去，不过bash脚本写起来可能稍微复杂点。还有就是石墨烯具有金属性，搞不好得考虑ismear≥0.

作者
Author: 小书童 时间: 2017-2-20 17:01

卡开发发发表于 2017-2-20 15:01
不一定会降低，只要平稳了就可以了。4x4石墨烯应该是32原子，500eV的截断能已经足以让能量收敛到1meV/atom ...

嗯嗯，谢谢老师讲解，是我自己理解错了，以为静态计算，截断能越大，能量越低。

作者
Author: denbaiyang 时间: 2021-1-7 22:31
时隔三年，大佬还在吗？请问测试截断能要达到一个什么样的结果才算合理呢？（就比如说石墨烯单胞的测试，我选择400eV到800eV，因为PBE赝势的建议是400eV）是先测试K点还是先测试截断能呢？

作者
Author: 文申的量子化学 时间: 2022-7-15 14:43

卡开发发发表于 2017-2-20 15:01
不一定会降低，只要平稳了就可以了。4x4石墨烯应该是32原子，500eV的截断能已经足以让能量收敛到1meV/atom ...

卡开老师，请问如何判断一个体系的截断能是不是适合这个分子呢？我做的是有机分子，能详细给说说吗

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