应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 强关联磁性体系SCF不收敛
121/2返回列表12下一页
查看: 3702  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

浔荆

木虫 (正式写手)

[求助] 强关联磁性体系SCF不收敛 已有1人参与

请教各位老师和研友，我在用CASTEP计算一个稀土氧化物(UO2)时，采用DFT+U的方法来处理f电子的强关联效应，体系呈反铁磁性。
问题是，试了很多参数，包括增大max scf cycles、调小mixing amplitude、增加smearing、增加empty bands、增大k点和截断能、更换贋势和泛函等，电子自洽SCF均不收敛，一直在-3到-5量级之间震荡。
使用all bands/edft倒是收敛性不错，但是计算代价太高，机器内存不足跑不动。
所以想请教一下各位，对于这类稀土磁性元素，有没有什么比较好的方法来促进收敛呢？
回复此楼

» 收录本帖的淘帖专辑推荐

稀土 第一性原理计算

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-04-27 21:32:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

【答案】应助回帖

★ ★ ★ ★ ★
franch: 金币+5, 应助指数+1, 谢谢回帖交流。。。 2014-05-05 22:20:16
franch: 回帖置顶 2014-05-05 22:20:19
引用回帖:
4楼: Originally posted by 浔荆 at 2014-04-27 22:29:42
谢谢你的回复，为何增加spin mixing amplitude呢？我之前一直都在减少spin mix......all bands/edft虽然能收敛，但是只能在小体系上算，放到大的超胞上算会超出内存，所以不是精度的问题。...

恩，似乎是我打错了，应该减小spin mix和charge mix，而增加对应cutoff vector。。。突然想起来前两天和别人说的时候也写错了。。

我指的是你先采用allbands/edft的coarse精度进行优化，经过优化构型有所调整，在这个新的结构上再使用density mix收敛可能会容易收敛一些（这样的情况我碰到过）。计算超胞的话先在小体系优化收敛后建立超胞在重新优化。

目前来看，即便算得到结果结果应该也不会好，这样的体系确实头疼。。。。
赞一下(3人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
11楼2014-04-27 23:36:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-05-05 22:17:43
empty bands要充分，k和Ecut也要合理。如果呈现不断的震荡没有收敛的趋势，即便采用很大的max scf cyc也是无意义的，可以适当减小charge mix及适当增加spin mix。

如果allbands/edft可以收敛，那么不妨试试用coarse精度下的allbands/edft做优化，达到几何优化收敛，换成更高的精度转为density mixing的方法进行优化。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
2楼2014-04-27 22:03:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
浔荆: 金币+2, ★有帮助 2014-04-27 22:29:53
4f加U是很难收敛，castep这方面默认参数基本上是算不出结果的。主要还是要加大空带条数。

castep现在完全可以加U做结果优化，从你的信息来看，你的结构优化并没有加U。其实+U优化和没有加结构还是有区别的。
赞一下 回复此楼
3楼2014-04-27 22:13:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

浔荆

木虫 (正式写手)
引用回帖:
2楼: Originally posted by 卡开发发 at 2014-04-27 22:03:27
empty bands要充分，k和Ecut也要合理。如果呈现不断的震荡没有收敛的趋势，即便采用很大的max scf cyc也是无意义的，可以适当减小charge mix及适当增加spin mix。

如果allbands/edft可以收敛，那么不妨试试用coa ...

谢谢你的回复，为何增加spin mixing amplitude呢？我之前一直都在减少spin mix......all bands/edft虽然能收敛，但是只能在小体系上算，放到大的超胞上算会超出内存，所以不是精度的问题。
赞一下 回复此楼
4楼2014-04-27 22:29:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

浔荆

木虫 (正式写手)
引用回帖:
3楼: Originally posted by mywai520 at 2014-04-27 22:13:18
4f加U是很难收敛，castep这方面默认参数基本上是算不出结果的。主要还是要加大空带条数。

castep现在完全可以加U做结果优化，从你的信息来看，你的结构优化并没有加U。其实+U优化和没有加结构还是有区别的。

empty bands 已经加到120%了，收敛性还是很差。现在我的问题不是结构问题，而是一次scf自洽循环都完成不了。。。
赞一下 回复此楼
5楼2014-04-27 22:31:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-05-05 22:18:26
你用的方法都不同了，精度设置上肯定不同。所以你若是不用all bands/edft，必须要把精度调高。这个MS说明书中就讲的很清楚。

不要比较两个不同方法下同一体系的参数，这样么有意义。
赞一下 回复此楼
6楼2014-04-27 22:32:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

浔荆

木虫 (正式写手)
引用回帖:
6楼: Originally posted by mywai520 at 2014-04-27 22:32:39
你用的方法都不同了，精度设置上肯定不同。所以你若是不用all bands/edft，必须要把精度调高。这个MS说明书中就讲的很清楚。

不要比较两个不同方法下同一体系的参数，这样么有意义。

没太懂你的意思。。。
你是说需要调高截断能和k点才能是density mixing收敛么？
赞一下 回复此楼
7楼2014-04-27 23:12:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)
引用回帖:
7楼: Originally posted by 浔荆 at 2014-04-27 23:12:39
没太懂你的意思。。。
你是说需要调高截断能和k点才能是density mixing收敛么？...

是的
回复此楼
8楼2014-04-27 23:16:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

浔荆

木虫 (正式写手)
引用回帖:
8楼: Originally posted by mywai520 at 2014-04-27 23:16:22
是的...

你好，请问这个在帮助文件的哪里呢？关于all bands/edft 帮助里提的不多，对他还没有太多理解
赞一下 回复此楼
9楼2014-04-27 23:25:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)
★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-05-05 22:18:48
引用回帖:
9楼: Originally posted by 浔荆 at 2014-04-27 23:25:29
你好，请问这个在帮助文件的哪里呢？关于all bands/edft 帮助里提的不多，对他还没有太多理解...

在MS 帮助直接搜all bands/edft，这个版主文件的确不多。

上面说的是+U计算，参数设置的要求还是很详细的。
赞一下 回复此楼
10楼2014-04-27 23:35:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 浔荆 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 申请科研项目纯靠实力能成功吗？ 200+7 David1982 2026-10-13 11/550 2026-10-15 10:30 by zyqchem
[职场人生] 40+中年人失业后，还能找到好工作吗？ +5 xqy2009 2026-10-13 5/250 2026-10-15 09:49 by 淘淘塑物性通
[硕博家园] 极限拉扯，如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么？ +4 好运来～9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭