| 查看: 5309 | 回复: 10 | |||
[交流]
大家一起总结一下MS中SCF不收敛的解决策略
|
|||
|
1.增大SCF圈数 2.减小density mixing 的 charge值减小震荡 3.加大空键 4.增大SCF收敛容差，在此基础上继续降低 5.更改smearing值大小 请问大家还有没有其他的关键参数，出入mS,一起交流 |
回复此楼» 猜你喜欢
越来越担心硕士毕业后没有竞争力了
已经有14人回复
-
论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章
已经有3人回复
-
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请
已经有3人回复
-
闲聊
已经有6人回复
-
申博求助
已经有4人回复
» 抢金币啦！回帖就可以得到:
-
海南大学邱天宇老师团队招收博士生与招聘博士后
+5/605
-
坐标苏州，诚聘女友
+1/286
-
溶质液晶释放必须用PBS吗？
+1/96
-
南开大学环境科学与工程学院 环境与人体健康方向 招博士研究生、硕士研究生各1名
+1/69
-
溜达闲逛
+1/49
-
如何重现小木虫论坛当年盛况？
+1/43
-
北京理工大学深圳汽车研究院长期招收联培学生/实习生
+1/36
-
大连理工大学化学学院招收2027级理论计算化学方向博士生
+1/35
-
【王宁利院士团队】招聘科研助理10人
+1/20
-
浙大许军课题组招聘
+1/19
-
211大学课题组招收数字公共治理相关的学术博士生1名
+1/10
-
上海交通大学-苏州国家实验室联合培养博士生招生（2027年春/秋季入学·申请考核制）
+1/9
-
开心
+1/7
-
【招聘】荷兰格罗宁根大学CSC公派博士生，材料/高分子/化学相关背景
+1/7
-
【2027博士招生】福州大学姚建年院士课题组招收博士研究生2名
+1/5
-
深圳技术大学招聘天体化学、物理化学领域博士后
+1/4
-
百万单元软组织模型一装配就爆内存：生物力学血流仿真遇上大模型本地部署工作站
+1/3
-
二审Required Reviews Complete了
+1/3
-
全球扩展虚拟生物学计划：18亿美元资源投向何处？
+1/2
-
【招聘】荷兰格罗宁根大学CSC公派博士生，材料/高分子/化学相关背景
+1/1
| 如果计算在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): 谢谢参与
·