24小时热门版块排行榜    

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

dhc198611

木虫 (小有名气)

[交流] 【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算 已有8人参与

我用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算,采用nelec=10,nelup =5,neldw =5, 设置反铁磁有序的。但是输出的结果到100步还是不收敛,而且奇怪的是输出结果居然没有显示受力。不知是怎么回事?希望哪位大侠帮忙解决。小弟非常感谢!
回复此楼
科研无止境
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★
dhc198611(金币+2):感谢你的指点 2010-05-13 22:44:39
gavinliu7390(金币+2):谢谢提示! 2010-05-14 07:28:52
需要说明的是nstep不是电子自洽场迭代的步数,电子迭代步数应该在&electron这个模块设置,参数是electron_maxstep,你可以把这个参数设置为任意值。nstep指的是relax,也就是驰誉的步数。
另外,设置反铁磁的例子在程序包里面有,你可以看一下。一般来说不是设置电子的个数,而是设置原子的初始磁矩大小和方向。你可以把input贴出来大家看看,如果你对大家有防御之心,就没必要在这里问了。
7楼2010-05-13 22:29:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答

157806400

金虫 (初入文坛)

至少应该有输入输出文件吧?

dhc198611(金币+2):谢谢参与 2010-05-13 14:11:20
你这样描述很抽象。
2楼2010-05-13 11:16:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dhc198611

木虫 (小有名气)

引用回帖:
Originally posted by 157806400 at 2010-05-13 11:16:54:
你这样描述很抽象。

这是我的INPIUT
    nelec= 10,
    nelup = 5,
    neldw = 5,
    occupations='fixed',  degauss=0.0
    lda_plus_u=.true., Hubbard_U(2)=2.0, Hubbard_U(3)=2.0,
    starting_ns_eigenvalue(3,2,2) = 1.d0
    starting_ns_eigenvalue(3,1,3) = 1.d0

这是加磁性和强关联势,这样设置反铁磁对不对
科研无止境
3楼2010-05-13 14:11:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★
dhc198611(金币+2):感谢你的宝贵意见,但是nstep最大也只能设置成100。 2010-05-13 21:19:05
zzy870720z(金币+2):感谢指导 2010-05-13 22:25:49
我怀疑你的100步不收敛这里指的是scf计算不收敛,这是电子结构计算的问题,有些参数没有设置准确。你的scf计算都没有收敛,当然不会有力的输出。
relax迭代步数的设置是nstep,你可以设置大一些,计算收敛后会自动停止的。
4楼2010-05-13 16:44:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 BMTdnp4OJO9f 2026-11-06 4/200 2026-11-07 19:36 by Xc2FxHsHkBNc
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 BMTdnp4OJO9f 2026-11-06 4/200 2026-11-07 19:03 by Xc2FxHsHkBNc
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 BMTdnp4OJO9f 2026-11-06 4/200 2026-11-07 18:52 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 17:01 by Xc2FxHsHkBNc
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 4/200 2026-11-07 16:28 by Xc2FxHsHkBNc
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 16:17 by Xc2FxHsHkBNc
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 16:06 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 5sJFDKp2kwuI 2026-11-06 6/300 2026-11-07 15:44 by Xc2FxHsHkBNc
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
信息提示
请填处理意见