24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2879  |  回复: 23
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dhc198611

木虫 (小有名气)

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

我用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算,采用nelec=10,nelup =5,neldw =5, 设置反铁磁有序的。但是输出的结果到100步还是不收敛,而且奇怪的是输出结果居然没有显示受力。不知是怎么回事?希望哪位大侠帮忙解决。小弟非常感谢!
回复此楼
科研无止境
已阅   回复此楼   关注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的,和你这个很相似。
9楼2010-05-13 23:45:30
已阅   回复此楼   关注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,你可以设置大一些,计算收敛后会自动停止的。
4楼2010-05-13 16:44:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见