24小时热门版块排行榜    

查看: 1480  |  回复: 3

achuan90

木虫 (小有名气)

[求助] 请教关于收敛精度

在文章中经常看到迭代过程收敛精度为2×10- 5 eV·atom- 1,为什么就说“即作用在每个原子上的力不大于0.5 eV·nm- 1, 内应力不大于0.1 GPa,”
在abinit中收敛精度设置是不是就是toldfe和tolvrs,两者又有什么区别啊。
如果迭代过程收敛精度为2×10- 5 eV·atom- 1,在abinit是这样设置么toldfe  7.4×10- 8么。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
achuan90: 金币+5, ★★★很有帮助 2013-09-04 18:16:02
我觉得LZ多加了个“即”字,或者写文章的人就没搞清楚自己在说什么。之前的收敛精度是指SCF计算的收敛阀值,而后面的力和内应力则分别对应了原子位置计算的力收敛阀值和晶格形貌计算的应力阀值。在Abinit中,这三个则是分别对应了toldfe,toldff,和strtarget。
SCF的收敛阀值可以使用toldfe,也可以使用tolvrs,两者的区别在Abinit官方网站说得很明白,前者是指前后两步的总能差阀值,后者是计算所得势与上步所得势的差别阀值。
如果需要2*10^(-5) eV/atom,那么对于计算有N个原子的体系,toldfe需要设置为:2*10^(-5) / ( 27.2113845 * N ),当然,也可以直接把toldfe设置为eV单位的,这样只需要除原子个数就可以了。
思想重于技巧,内涵重于表象
2楼2013-09-04 10:13:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

achuan90

木虫 (小有名气)

引用回帖:
2楼: Originally posted by souledge at 2013-09-04 10:13:34
我觉得LZ多加了个“即”字,或者写文章的人就没搞清楚自己在说什么。之前的收敛精度是指SCF计算的收敛阀值,而后面的力和内应力则分别对应了原子位置计算的力收敛阀值和晶格形貌计算的应力阀值。在Abinit中,这三个 ...

你好!感谢你的答复。另外:
如果需要2*10^(-5) eV/atom,那么对于计算有N个原子的体系,toldfe应该是上面乘以N吧。也就是toldfe需要设置为2*10^(-5) / ( 27.2113845 )* N而不是除以吧。
3楼2013-09-04 22:55:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

引用回帖:
3楼: Originally posted by achuan90 at 2013-09-04 22:55:12
你好!感谢你的答复。另外:
如果需要2*10^(-5) eV/atom,那么对于计算有N个原子的体系,toldfe应该是上面乘以N吧。也就是toldfe需要设置为2*10^(-5) / ( 27.2113845 )* N而不是除以吧。...

Yes, you are right. It was my mistake~
思想重于技巧,内涵重于表象
4楼2013-09-04 23:06:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 achuan90 的主题更新
信息提示
请填处理意见