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

24小时热门版块排行榜    

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