24小时热门版块排行榜    

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

dhc198611

木虫 (小有名气)

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

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

goldenfisher

金虫 (著名写手)

dhc198611(金币+2): 2010-05-13 21:17:54
另外,我步明白你为什么要设置starting_ns_eigenvalue这个参数,你确定你的设置是对的吗?如果不确定,建议你用默认值
5楼2010-05-13 16:46:41
已阅   回复此楼   关注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贴出来大家看看,如果你对大家有防御之心,就没必要在这里问了。
7楼2010-05-13 22:29:34
已阅   回复此楼   关注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的回帖

goldenfisher

金虫 (著名写手)


zzy870720z(金币+1):感谢专家指导 2010-05-15 21:14:42
dhc198611(金币+2):感谢指点!对于这种加磁性和强关联势的体系我要好好再研究一下了。 2010-05-15 21:30:41
PWSCF里面没有专门控制对称性的参数,如果你可以不考虑对称性,可以考虑使用nosym这个参数。
12楼2010-05-15 21:11:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dhc198611 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见