24小时热门版块排行榜    

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

dhc198611

木虫 (小有名气)

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

我用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算,采用nelec=10,nelup =5,neldw =5, 设置反铁磁有序的。但是输出的结果到100步还是不收敛,而且奇怪的是输出结果居然没有显示受力。不知是怎么回事?希望哪位大侠帮忙解决。小弟非常感谢!
回复此楼
科研无止境
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dhc198611

木虫 (小有名气)

引用回帖:
Originally posted by goldenfisher at 2010-05-13 23:45:30:
ecutwfc       = 35.0可能小了点,另外应该设置ecutrho=8*ecutwfc.另外,你根本没有设置体系的自旋初始值,初始是默认为0的。start_magnetization这个参数你应该设置好。
我觉得你应该设置Fe1的自旋为1,Fe2的自 ...

现在自洽可以收敛,我忘记加start_magnetization这个参数了,主要我在看INPUT_PWSCF时发现start_magnetization的参数说明的时候发现这句话(If you fix the magnetization with "nelup/neldw" or with "multiplicity" or with "tot_magnetization", you should
not specify starting_magnetization.)所以就没加上,现在发现自己以前的想法是错的。

不过现在又出现了新的问题,就是在加反铁磁性计算时,其体系的对称性发生了改变,加磁性之后,体系的对称性更加低,有些原子偏离的非常明显,我在看INPUT_PWSCF参数设置的时候并没有发现可以控制对称性的参数,不知这要怎么设置才能使其对称性保持不变(比如像VASP中ISYM就控制对称性的)!期待大侠的回复。
科研无止境
11楼2010-05-14 11:15:36
已阅   回复此楼   关注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的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见