| 查看: 1750 | 回复: 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么。 |
回复此楼» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有13人回复
-
2026国自然函评费到账
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
感觉是下周放榜了
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何查看计算精度？
已经有16人回复
- 关于杂化泛函的问题
已经有4人回复
- 用离散相模拟无法收敛，问题何在？
已经有6人回复
- 请教LSDA+U计算的问题
已经有24人回复
- 我这个结构什么时候能够收敛？
已经有12人回复
- neb 收敛精度
已经有7人回复
- 急切请教关于范数求导的问题
已经有10人回复
- NEB计算扩散路径的力收敛精度设置问题
已经有16人回复
- 如何鉴别网络质量的好坏？
已经有32人回复
- 收敛精度
已经有7人回复
- 电子步不收敛，跑了60多电子步第一个离子步还没有完成，这是怎么回事？
已经有28人回复
- 关于985学校自主设置的博士点要不要读……是不是表示专业不太好啊？请轻拍……
已经有15人回复
- 【求助】K点取得越大越不收敛，什么原因！
已经有20人回复
- 【求助】收敛性测试的问题
已经有40人回复
- 【求助】关于精度和k点设置的疑问
已经有13人回复
- 【请教】请教一个关于图片精度的问题
已经有6人回复
- 【求助】请教残差曲线震荡的解决问题
已经有6人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
- 【求助结束】结构优化的时候如何提高力的收敛精度？
已经有10人回复
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★
achuan90: 金币+5, ★★★很有帮助 2013-09-04 18:16:02
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
achuan90
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2254.5
- 帖子: 145
- 在线: 51.2小时
- 虫号: 310682
- 注册: 2006-12-24
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2013-09-04 22:55:12
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
Yes, you are right. It was my mistake~ 4楼2013-09-04 23:06:06
