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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算
241/3返回列表123下一页
查看: 3119  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +21 rare12345 2026-03-18 21/1050 2026-03-20 14:31 by 无懈可击111
[考研] 本人考085602 化学工程 专硕 +18 不知道叫什么！ 2026-03-15 20/1000 2026-03-20 13:52 by danney002
[考研] 一志愿中南化学337求调剂 +4 niko- 2026-03-19 5/250 2026-03-20 13:49 by 促天成
[考研] 279分求调剂 一志愿211 +9 chaojifeixia 2026-03-19 10/500 2026-03-20 12:29 by lature00
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +4 西南交通专材355 2026-03-19 4/200 2026-03-20 11:39 by 花开富贵幸福人�
[考研] 317求调剂 +4 申子申申 2026-03-19 9/450 2026-03-20 11:08 by 申子申申
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 328求调剂，英语六级551，有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 招收0805（材料）调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
信息提示
关闭
请填处理意见
关闭