请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 73|回复: 3

[CASTEP/MS] 求助,castep收敛性测试的判断标准是总能还是单原子能

[复制链接]

27

帖子

0

威望

59

eV
积分
86

Level 2 能力者

发表于 7 天前 | 显示全部楼层 |阅读模式
当时计算黄铁矿晶体时选择fine精度软件预估的截断能是300eV,为了更准确用了350eV,审稿人提出质疑350eV的截断能太低,尤其对于磁性体系。查看文献中做截断能的收敛性测试有两种情况,一种根据两个参数总能的差值,一种根据单原子能的差值(eV/atom)。想请教一下两种方法以哪个为准,具体的判断数值是多少呢?我测试的350和400eV时总能差0.05eV,一共12个原子,单原子能差值应该是0.004eV/atom,不知道这个数值能符合要求吗?有些文献没有写具体数值,只给出了曲线,有的写的是1meV/atom,感觉后者要求很高。

1460

帖子

3

威望

5236

eV
积分
6756

Level 6 (一方通行)

Ab Initio Amateur

发表于 7 天前 | 显示全部楼层
相对能量收敛要确实比总能收敛更容易,只关注相对能量的情况其实相对能量收敛就可以,但势能面比较复杂的情况建议还是采用参考总能和力的收敛会比较好。1meV/atom间隔50eV的动能截断收敛可能也没那么难,很大的体系可以适当放宽要求。一般1meV/atom是因为VASP手册上那么写,但我觉得这个标准应当按照你SCF的收敛标准和构型优化的标准进行调整。

评分

参与人数 1eV +2 收起 理由
sobereva + 2

查看全部评分

满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。

27

帖子

0

威望

59

eV
积分
86

Level 2 能力者

 楼主| 发表于 7 天前 | 显示全部楼层
卡开发发 发表于 2019-3-15 14:48
相对能量收敛要确实比总能收敛更容易,只关注相对能量的情况其实相对能量收敛就可以,但势能面比较复杂的情 ...

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

1460

帖子

3

威望

5236

eV
积分
6756

Level 6 (一方通行)

Ab Initio Amateur

发表于 7 天前 | 显示全部楼层
hgkjqx 发表于 2019-3-15 15:42
谢谢老师。Castep计算时用的fine精度,SCF收敛标准10E-6eV/atom,几何优化收敛标准10E-5eV/atom,默认的f ...

你可以再测试几个值看看,然后把收敛的情况列给审稿人看。相对能量如果只是4meV/atom应该比较宽松吧,默认的动能截断不一定靠得住,所以最好还是要做收敛性测试。测试应该做在前面以保证计算的可靠性,不要审稿人问了你再做,返工就很难办了。
满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-3-22 17:00 , Processed in 0.139069 second(s), 24 queries .

快速回复 返回顶部 返回列表