计算化学公社

标题: 求助，castep收敛性测试的判断标准是总能还是单原子能 [打印本页]

作者
Author: hgkjqx 时间: 2019-3-15 11:33
标题: 求助，castep收敛性测试的判断标准是总能还是单原子能
当时计算黄铁矿晶体时选择fine精度软件预估的截断能是300eV，为了更准确用了350eV,审稿人提出质疑350eV的截断能太低，尤其对于磁性体系。查看文献中做截断能的收敛性测试有两种情况，一种根据两个参数总能的差值，一种根据单原子能的差值（eV/atom）。想请教一下两种方法以哪个为准，具体的判断数值是多少呢？我测试的350和400eV时总能差0.05eV，一共12个原子，单原子能差值应该是0.004eV/atom，不知道这个数值能符合要求吗？有些文献没有写具体数值，只给出了曲线，有的写的是1meV/atom，感觉后者要求很高。

作者
Author: 卡开发发 时间: 2019-3-15 14:48
相对能量收敛要确实比总能收敛更容易，只关注相对能量的情况其实相对能量收敛就可以，但势能面比较复杂的情况建议还是采用参考总能和力的收敛会比较好。1meV/atom间隔50eV的动能截断收敛可能也没那么难，很大的体系可以适当放宽要求。一般1meV/atom是因为VASP手册上那么写，但我觉得这个标准应当按照你SCF的收敛标准和构型优化的标准进行调整。

作者
Author: hgkjqx 时间: 2019-3-15 15:42

卡开发发发表于 2019-3-15 14:48
相对能量收敛要确实比总能收敛更容易，只关注相对能量的情况其实相对能量收敛就可以，但势能面比较复杂的情 ...

谢谢老师。Castep计算时用的fine精度，SCF收敛标准10E-6eV/atom，几何优化收敛标准10E-5eV/atom，默认的fine精度是300eV截断能，我又给加了50eV，感觉应该够了，审稿人问了，我测试发现从350到400之间增加0.004eV/atom，不知道这样审稿人能给过关吗，后面计算量还多着呢，要第一步就不给过就得全部重新算了

作者
Author: 卡开发发 时间: 2019-3-15 16:01

hgkjqx 发表于 2019-3-15 15:42
谢谢老师。Castep计算时用的fine精度，SCF收敛标准10E-6eV/atom，几何优化收敛标准10E-5eV/atom，默认的f ...

你可以再测试几个值看看，然后把收敛的情况列给审稿人看。相对能量如果只是4meV/atom应该比较宽松吧，默认的动能截断不一定靠得住，所以最好还是要做收敛性测试。测试应该做在前面以保证计算的可靠性，不要审稿人问了你再做，返工就很难办了。

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