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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算
241/3返回列表123下一页
查看: 3150  |  回复: 23
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dhc198611

木虫 (小有名气)

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

我用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算，采用nelec=10,nelup =5,neldw =5, 设置反铁磁有序的。但是输出的结果到100步还是不收敛，而且奇怪的是输出结果居然没有显示受力。不知是怎么回事？希望哪位大侠帮忙解决。小弟非常感谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
科研无止境
1楼 2010-05-13 08:19:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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，你可以设置大一些，计算收敛后会自动停止的。
赞一下(1人) 回复此楼
4楼2010-05-13 16:44:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
dhc198611(金币+2): 2010-05-13 21:17:54
另外，我步明白你为什么要设置starting_ns_eigenvalue这个参数，你确定你的设置是对的吗？如果不确定，建议你用默认值
赞一下 回复此楼
5楼2010-05-13 16:46:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dhc198611

木虫 (小有名气)
引用回帖:
Originally posted by goldenfisher at 2010-05-13 16:46:41:
另外，我步明白你为什么要设置starting_ns_eigenvalue这个参数，你确定你的设置是对的吗？如果不确定，建议你用默认值

starting_ns_eigenvalue这个是+U的设置，即使默认页还是不收敛。迭代100步的时候，不收敛就不会输出力，还有PWSCF的默认迭代次数为100，即使你设置nstep=200，程序也会在100步的时候自动停掉。
我的主要问题在于，当加入
  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
这几句的时候，其体系就是不收敛，去掉的话就能收敛了，所以可能是我的反铁磁有序设置有问题。
请问有没有人知道怎么设置  nelec，nelup ， neldw 这3个参数，这3个参数怎么样设置成反铁磁有序。对于这3个参数我不太了解，但是按计算的提示是要加的。
赞一下 回复此楼
科研无止境
6楼2010-05-13 21:17:43
已阅   回复此楼   关注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贴出来大家看看，如果你对大家有防御之心，就没必要在这里问了。
赞一下(1人) 回复此楼
7楼2010-05-13 22:29:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dhc198611

木虫 (小有名气)
引用回帖:
Originally posted by goldenfisher at 2010-05-13 22:29:34:
需要说明的是nstep不是电子自洽场迭代的步数，电子迭代步数应该在&electron这个模块设置，参数是electron_maxstep,你可以把这个参数设置为任意值。nstep指的是relax，也就是驰誉的步数。
另外，设置反铁磁的 ...

&control
    calculation='relax',
    restart_mode='from_scratch',
    pseudo_dir = '~/pseudo/',
    outdir='.',
    prefix='BFO',
    tprnfor = .true.
    tstress = .true.
    etot_conv_thr = 1e-5
    forc_conv_thr = 5e-5
    nstep=200
/
&system   
    ibrav=  5,
    celldm(1) =10.7241882,
    celldm(4) =0.51099,
    nat=10,
    ntyp= 4,
    ecutwfc       = 35.0
    nspin = 2
    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,
/
&electrons
    mixing_beta   = 0.3,
/
&IONS
/
ATOMIC_SPECIES
Bi  208.98  Bi.pbe-d-mt.UPF
Fe1  55.845   Fe.pbe-sp-van_ak.UPF
Fe2  55.845   Fe.pbe-sp-van_ak.UPF
O  15.9994  O.pbe-van_ak.UPF
ATOMIC_POSITIONS {crystal}
Bi  0.0000000000000000  0.0000000000000000  0.0000000000000000   0  0  0
Bi  0.5000000000000000  0.5000000000000000  0.5000000000000000   0  0  0
Fe1  0.2225464345999172  0.2225464345999172  0.2225464345999172
Fe2  0.7225464345999169  0.7225464345999169  0.7225464345999169  
O  0.5367314264732763 -0.0624437787918008  0.3877084883892429   
O -0.0624437787918008  0.3877084883892430  0.5367314264732763   
O  0.3877084883892429  0.5367314264732763 -0.0624437787918008   
O  0.0367314264732762  0.8877084883892430  0.4375562212081991   
O  0.8877084883892431  0.4375562212081992  0.0367314264732762   
O  0.4375562212081991  0.0367314264732761  0.8877084883892430   
K_POINTS {automatic}
  4 4 4 0 0 0
   
这是我的INPUT，我想收敛的话，应该能在100内能搞定的，但是电子步就是不收敛。我的初始位置是经过VASP弛豫的，在PWSCF里应该不会差很多的。
赞一下 回复此楼
科研无止境
8楼2010-05-13 22:44:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★
gavinliu7390(金币+2):谢谢提示！ 2010-05-14 07:29:06
dhc198611(金币+4):非常感谢 2010-05-14 07:44:26
dhc198611(金币+1):期待你的回复（关于如何控制PWSCF弛豫时的对称性） 2010-05-14 11:17:08
ecutwfc       = 35.0可能小了点，另外应该设置ecutrho=8*ecutwfc.另外，你根本没有设置体系的自旋初始值，初始是默认为0的。start_magnetization这个参数你应该设置好。
我觉得你应该设置Fe1的自旋为1，Fe2的自旋为-1，把体系的总磁矩置零。
我觉得你应该仔细看input_pwscf文件，然后再看例子。例子里有计算FeO的，和你这个很相似。
赞一下(1人) 回复此楼
9楼2010-05-13 23:45:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)
★
dhc198611(金币+1):谢谢参与
http://valenhou.blog.edu.cn/2007/133214.html

[ Last edited by lono75 on 2010-5-14 at 10:49 ]
赞一下(1人) 回复此楼
由天命,成人事！
10楼2010-05-14 10:26:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dhc198611 的主题更新
241/3返回列表123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 还有化工二轮调剂的学校吗 5+4 化工人999 2026-04-09 18/900 2026-04-09 15:49 by 紫翼精灵
[考研] 070300化学279求调剂 +17 哈哈哈^_^ 2026-04-08 18/900 2026-04-09 10:49 by 三七七想上岸
[考研] 材料工程322 +18 哈哈哈吼吼吼哈 2026-04-07 19/950 2026-04-09 10:44 by cymywx
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 化学308分求调剂 +21 你好明天你好 2026-04-07 23/1150 2026-04-08 22:32 by 凯凯要变帅
[考研] 一志愿山东大学化学与化工学院材料与化工专硕，360分求调剂 +8 不愿透露姓名的� 2026-04-02 8/400 2026-04-08 20:12 by luoyongfeng
[考研] 275 求调剂 +8 Lei812514 2026-04-07 8/400 2026-04-08 12:46 by chemisry
[考研] 调剂 +3 电气300求调剂不 2026-04-08 6/300 2026-04-08 09:39 by 电气300求调剂不
[考研] 计算机408｜在校多次国家级竞赛获奖｜申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 生物调剂 +5 橙子橙子橙子啊 2026-04-05 9/450 2026-04-07 15:31 by 上岸快快
[考研] 调剂 +4 mcbbc 2026-04-06 5/250 2026-04-07 12:33 by upczlm1989
[考研] 材料调剂 +17 小刘同学吖吖 2026-04-06 18/900 2026-04-07 11:41 by 诗与自由
[考研] 293求调剂 +3 勇远库爱314 2026-04-06 3/150 2026-04-07 11:15 by hugr
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 一志愿武汉理工大学080200机械工程308分，求调剂 +4 终不似从前 2026-04-05 4/200 2026-04-06 11:46 by 考研学校招点人
[考研] 生物与医药调剂 +4 十七sa 2026-04-05 4/200 2026-04-05 20:05 by lys0704
[考研] 计算机11408，286分求调剂 +7 木子念晞 2026-04-05 7/350 2026-04-05 19:02 by chy09050039
[考研] 材料调剂 +12 一样YWY 2026-04-04 12/600 2026-04-05 08:24 by 544594351
[考研] 309求调剂 +4 快乐的小白鸽 2026-04-04 5/250 2026-04-04 15:55 by cql1109
[考研] 求调剂 +3 农业工程与信息� 2026-04-04 3/150 2026-04-04 12:19 by 舍而后得
信息提示
关闭
请填处理意见
关闭