| 查看: 5204 | 回复: 10 | |||
[交流]
大家一起总结一下MS中SCF不收敛的解决策略
|
|||
|
1.增大SCF圈数 2.减小density mixing 的 charge值减小震荡 3.加大空键 4.增大SCF收敛容差，在此基础上继续降低 5.更改smearing值大小 请问大家还有没有其他的关键参数，出入mS,一起交流 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
| 如果计算在Max. cycle全部用完之后仍然不能收敛，要先查看一下.castep文件，如果能量比刚开始计算时有明显下降，而且震动幅度也明显减小，只是还没有降到收敛标准以下，那说明Max. cycle还不够，要增加这个数值。 但多数情况下，我们会发现，能量并没有比刚开始计算时的能量有明显下降，而是围绕着某个平均值在反复地震荡，而且震动幅度也没有减小。那说明是参数设定的问题，此时就算再怎么增加Max. cycle的数值，计算最终还是不会收敛的。 这个时候，最好是增加energy cutoff的数值，降低K point set的数值，增加smearing width数值（当然，如果通过调试前两个数值就可以使计算收敛的话，就最好不要动这个数值）。 再有，如果计算是使用的density mixing的话，那要将mix_charge_amp，mix_spin_amp的数值降低，将mix_energy_cutoff的数值增加到 energy cutoff数值的3到4倍，增加mix_history_length的数值。 当然，如果调试上面的参数仍然发现收敛很辛苦的话，就将mixing scheme转成All Band/EDFT,这会解决收敛问题，但是计算时间将会是使用density mixing计算所用时间的3倍以上 |
3楼2019-12-09 09:38:33
2楼2019-12-09 09:36:32
★
douliguang(金币+1): 谢谢参与
douliguang(金币+1): 谢谢参与
|
谢谢，受教了，有点疑问，smearing一般多少合理呢，看help里它是加速收敛的同时也降低精度，如果原来无法收敛，调大smearing却使精度也降低很多，这不等同于调低精度么，那么精度太低结果也不可信，这里精度是scf tolerance，是楼主说的收敛容差。 发自小木虫Android客户端 |
7楼2019-12-15 21:51:19
|
同理charge设置多少合理呢，这个charge看了help我还是不是很清楚含义，所以一般我也不动它。下面有一个preconditioner可以减小电荷密度震荡加速收敛，这个应该也可以选吧，help里也没有提到会有什么不良效应。 发自小木虫Android客户端 |
8楼2019-12-15 21:57:31
9楼2019-12-18 19:48:05
10楼2020-06-14 22:16:16
11楼2020-06-14 22:48:15
简单回复
2019-12-09 11:36
回复
douliguang(金币+1): 谢谢参与
psylhh5楼
2019-12-11 20:55
回复
douliguang(金币+1): 谢谢参与
nono20096楼
2019-12-14 13:08
回复
douliguang(金币+1): 谢谢参与
·
