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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算
51/1返回列表
查看: 2875  |  回复: 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 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭