24小时热门版块排行榜    

查看: 1606  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gemucai

木虫 (正式写手)

[交流] 分享一下,碰到scf计算不收敛已有1人参与

大家好,我最近碰到了scf计算不收敛的问题,进行了100个interaction跳出来。我试了增加K网格、换赝势、增大ecut、改变mix,都不行。

我后来动了动脑子,本来我算的是绝缘体系,但是我对其做了点手脚,会不会这点改动让它变成导体了?于是我改了occupations,结果奇迹发生了,不到20个interaction就收敛了,能带图画出来确实是导体。如果体系改动之后,碰到scf计算不收敛,这种方法也许可以参考一下。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by gemucai at 2012-11-29 15:21:14
occupations默认应该不是smearing,我看了以前没有设置occupation的输出文件,每个k点对应的本征值数目都是电子数的一半。而一旦设置了smearing,默认的就是gauss扩展了。...

对对~脑子短路了,把occupations和smearing给弄混了……
默认是啥来着?不做任何处理么?那样似乎很难收敛啊……
思想重于技巧,内涵重于表象
4楼2012-11-29 18:27:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

souledge

专家顾问 (著名写手)

occupations不是自动就用的么?好象默认是Gaussian扩展吧?
思想重于技巧,内涵重于表象
2楼2012-11-29 10:03:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
2楼: Originally posted by souledge at 2012-11-29 10:03:25
occupations不是自动就用的么?好象默认是Gaussian扩展吧?

occupations默认应该不是smearing,我看了以前没有设置occupation的输出文件,每个k点对应的本征值数目都是电子数的一半。而一旦设置了smearing,默认的就是gauss扩展了。
3楼2012-11-29 15:21:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
4楼: Originally posted by souledge at 2012-11-29 18:27:09
对对~脑子短路了,把occupations和smearing给弄混了……
默认是啥来着?不做任何处理么?那样似乎很难收敛啊……...

我猜想默认是fixed,仅仅是猜想,因为INPUT_PW也没给出来
5楼2012-11-29 19:52:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见