24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4090  |  回复: 10

douliguang

银虫 (正式写手)


[交流] 大家一起总结一下MS中SCF不收敛的解决策略

1.增大SCF圈数
2.减小density mixing 的 charge值减小震荡
3.加大空键
4.增大SCF收敛容差,在此基础上继续降低
5.更改smearing值大小


请问大家还有没有其他的关键参数,出入mS,一起交流
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

douliguang

银虫 (正式写手)


如果计算在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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

douliguang

银虫 (正式写手)


不收敛的时候,首先是增加SCF的迭代次数,其次是减少density mixing的值,先用大的smearing,后减少,你试试,这样还不收敛的话,就是你的结构或者软件或者计算机的问题。
2楼2019-12-09 09:36:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

丁风缘

铁虫 (小有名气)



douliguang(金币+1): 谢谢参与
谢谢,受教了,有点疑问,smearing一般多少合理呢,看help里它是加速收敛的同时也降低精度,如果原来无法收敛,调大smearing却使精度也降低很多,这不等同于调低精度么,那么精度太低结果也不可信,这里精度是scf tolerance,是楼主说的收敛容差。

发自小木虫Android客户端
7楼2019-12-15 21:51:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

丁风缘

铁虫 (小有名气)


同理charge设置多少合理呢,这个charge看了help我还是不是很清楚含义,所以一般我也不动它。下面有一个preconditioner可以减小电荷密度震荡加速收敛,这个应该也可以选吧,help里也没有提到会有什么不良效应。

发自小木虫Android客户端
8楼2019-12-15 21:57:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

丁风缘

铁虫 (小有名气)


说错了,smearing不降精度,harris approximation降精度

发自小木虫Android客户端
9楼2019-12-18 19:48:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

DR_Zhou

新虫 (正式写手)



douliguang(金币+1): 谢谢参与
请教一下,如果换成All band/EDFT之后也无法收敛该怎么办呢?
10楼2020-06-14 22:16:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fgmsyutong

木虫 (文坛精英)



douliguang(金币+1): 谢谢参与
把density mixing改为all bands/edft

发自小木虫Android客户端
11楼2020-06-14 22:48:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2019-12-09 11:36   回复  
douliguang(金币+1): 谢谢参与
psylhh5楼
2019-12-11 20:55   回复  
douliguang(金币+1): 谢谢参与
nono20096楼
2019-12-14 13:08   回复  
douliguang(金币+1): 谢谢参与
·
相关版块跳转 我要订阅楼主 douliguang 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见