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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算
51/1返回列表
查看: 3118  |  回复: 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的回帖

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的回帖
查看全部 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，你可以设置大一些，计算收敛后会自动停止的。
赞一下(1人) 回复此楼
4楼2010-05-13 16:44:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286分人工智能专业请求调剂愿意跨考！ +3 lemonzzn 2026-03-17 4/200 2026-03-20 11:04 by lemonzzn
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 14/700 2026-03-20 10:46 by EBSD
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +8 步川酷紫123 2026-03-13 8/400 2026-03-20 08:47 by xingguangj
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 本人考085602 化学工程 专硕 +17 不知道叫什么！ 2026-03-15 19/950 2026-03-19 15:06 by 尽舜尧1
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 一志愿华中科技大学，080502，354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
信息提示
关闭
请填处理意见
关闭