24小时热门版块排行榜&&&&
(小有名气)
在线: 573.9小时
虫号: 780288
注册: 性别: GG专业: 凝聚态物性 II ：电子结构
【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算已有8人参与
我用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算，采用nelec=10,nelup =5,neldw =5, 设置反铁磁有序的。但是输出的结果到100步还是不收敛，而且奇怪的是输出结果居然没有显示受力。不知是怎么回事？希望哪位大侠帮忙解决。小弟非常感谢！
& 猜你喜欢
已经有8人回复
已经有7人回复
已经有19人回复
已经有45人回复
已经有10人回复
已经有26人回复
已经有41人回复
已经有35人回复
已经有76人回复
已经有24人回复
& 本主题相关价值贴推荐，对您同样有帮助:
已经有56人回复
已经有3人回复
已经有4人回复
已经有31人回复
已经有6人回复
已经有6人回复
已经有10人回复
已经有13人回复
已经有9人回复
已经有4人回复
科研无止境
(初入文坛)
在线: 6.5小时
虫号: 796630
注册: 性别: GG专业: 功能材料计算
至少应该有输入输出文件吧？
dhc198611(金币+2):谢谢参与
你这样描述很抽象。
(小有名气)
在线: 573.9小时
虫号: 780288
注册: 性别: GG专业: 凝聚态物性 II ：电子结构
引用回帖:Originally posted by
你这样描述很抽象。 这是我的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
这是加磁性和强关联势，这样设置反铁磁对不对
科研无止境
(著名写手)
散金: 2790
在线: 346.9小时
虫号: 117764
注册: 性别: GG专业: 原子和分子物理
★ ★ dhc198611(金币+2):感谢你的宝贵意见，但是nstep最大也只能设置成100。
21:19:05zzy870720z(金币+2):感谢指导
我怀疑你的100步不收敛这里指的是scf计算不收敛，这是电子结构计算的问题，有些参数没有设置准确。你的scf计算都没有收敛，当然不会有力的输出。
relax迭代步数的设置是nstep，你可以设置大一些，计算收敛后会自动停止的。
(著名写手)
散金: 2790
在线: 346.9小时
虫号: 117764
注册: 性别: GG专业: 原子和分子物理
dhc198611(金币+2):
另外，我步明白你为什么要设置starting_ns_eigenvalue这个参数，你确定你的设置是对的吗？如果不确定，建议你用默认值
(小有名气)
在线: 573.9小时
虫号: 780288
注册: 性别: GG专业: 凝聚态物性 II ：电子结构
引用回帖:Originally posted by goldenfisher at
另外，我步明白你为什么要设置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个参数我不太了解，但是按计算的提示是要加的。
科研无止境
(著名写手)
散金: 2790
在线: 346.9小时
虫号: 117764
注册: 性别: GG专业: 原子和分子物理
★ ★ dhc198611(金币+2):感谢你的指点
22:44:39gavinliu7390(金币+2):谢谢提示！
需要说明的是nstep不是电子自洽场迭代的步数，电子迭代步数应该在&electron这个模块设置，参数是electron_maxstep,你可以把这个参数设置为任意值。nstep指的是relax，也就是驰誉的步数。
另外，设置反铁磁的例子在程序包里面有，你可以看一下。一般来说不是设置电子的个数，而是设置原子的初始磁矩大小和方向。你可以把input贴出来大家看看，如果你对大家有防御之心，就没必要在这里问了。
(小有名气)
在线: 573.9小时
虫号: 780288
注册: 性别: GG专业: 凝聚态物性 II ：电子结构
引用回帖:Originally posted by goldenfisher at
需要说明的是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,
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.0&&0.0000& &0&&0&&0
Bi&&0.0&&0.0000& &0&&0&&0
Fe1&&0.2&&0.9172
Fe2&&0.9&&0.9169&&
O&&0.8&&0.2429& &
O -0.0&&0.2763& &
O&&0.3 -0.8008& &
O&&0.0&&0.1991& &
O&&0.2&&0.2762& &
O&&0.1&&0.2430& &
K_POINTS {automatic}
&&4 4 4 0 0 0
这是我的INPUT，我想收敛的话，应该能在100内能搞定的，但是电子步就是不收敛。我的初始位置是经过VASP弛豫的，在PWSCF里应该不会差很多的。
科研无止境
(著名写手)
散金: 2790
在线: 346.9小时
虫号: 117764
注册: 性别: GG专业: 原子和分子物理
★ ★ gavinliu7390(金币+2):谢谢提示！
07:29:06dhc198611(金币+4):非常感谢
07:44:26dhc198611(金币+1):期待你的回复（关于如何控制PWSCF弛豫时的对称性）
ecutwfc& && & = 35.0可能小了点，另外应该设置ecutrho=8*ecutwfc.另外，你根本没有设置体系的自旋初始值，初始是默认为0的。start_magnetization这个参数你应该设置好。
我觉得你应该设置Fe1的自旋为1，Fe2的自旋为-1，把体系的总磁矩置零。
我觉得你应该仔细看input_pwscf文件，然后再看例子。例子里有计算FeO的，和你这个很相似。
(著名写手)
在线: 529.1小时
虫号: 932458
注册: 专业: 凝聚态物性 II ：电子结构
★ dhc198611(金币+1):谢谢参与
[ Last edited by lono75 on
at 10:49 ]
由天命,成人事！
相关版块跳转
第一性原理
我要订阅楼主
的主题更新
小木虫,学术科研互动社区,为中国学术科研免费提供动力
违规贴举报删除请联系客服电话： 邮箱：（全天候） 或者 QQ:
广告投放与宣传请联系 李想 QQ：
QQ：&&邮箱：
Copyright & 2001-, All Rights Reserved. 小木虫 版权所有
浏览器进程
打开微信扫一扫
随时随地聊科研24小时热门版块排行榜&&&&
【悬赏金币】回答本帖问题，作者ayikongjian将赠送您 10 个金币
(初入文坛)
在线: 32.1小时
虫号: 535631
注册: 专业: 凝聚态物性 II ：电子结构
pwscf能否做含有磁性体系的电声耦合计算
各位pwscf高手，
& && &&&刚接触pwscf不久，想问一个也许比较白痴的问题：pwscf计算超导临界温度，能否做含有磁性体系的电声耦合计算呢？好像听人说过超导和磁性不能共存之类的，那理论上讲pwscf能否计算磁性体系的电声耦合呢？具体怎么操作的呢？和一般非磁性体系计算过程一致吗？
& 猜你喜欢
已经有8人回复
已经有7人回复
已经有19人回复
已经有45人回复
已经有10人回复
已经有26人回复
已经有41人回复
已经有35人回复
已经有76人回复
已经有24人回复
& 本主题相关价值贴推荐，对您同样有帮助:
已经有26人回复
已经有11人回复
已经有34人回复
已经有7人回复
已经有13人回复
已经有16人回复
已经有5人回复
已经有25人回复
已经有7人回复
已经有31人回复
已经有7人回复
已经有12人回复
已经有14人回复
已经有28人回复
已经有7人回复
已经有23人回复
已经有13人回复
已经有21人回复
(著名写手)
散金: 2790
在线: 346.9小时
虫号: 117764
注册: 性别: GG专业: 原子和分子物理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ayikongjian: 金币+10, ★★★很有帮助
引用回帖:: Originally posted by 刘仕晨 at
是的，但是这个必须使得occupation为smearing，我是绝缘体系，而且要计算波恩有效电荷，按照手册需要设置occupation为fixed，这个时候starting_magnetization就不能使用了，只能使用tot_magnetization，但是这样的 ... 原则上，没有一定要限制occupation为fixed。对绝缘体，同样可以设置smearing，但是这个时候你要注意degauss的数值必须足够小。这样你可以算出与fixed一样的结果。
(著名写手)
散金: 2790
在线: 346.9小时
虫号: 117764
注册: 性别: GG专业: 原子和分子物理
引用回帖:: Originally posted by 刘仕晨 at
谢谢goldenfisher老师，早就听说老师的大名了，没想到这次得到了老师的指导，那我计算声子的时候，因为要考虑LO-TO的spliting，所以要设置epsil=.true.来计算波恩有效电荷和介电张量，这个时候设置为smearing也没有 ... 首先，如果你的带隙足够大，那么degauss没有必要多小。原则上，小于带隙宽度即可，当然为了准确，还是小点好。我一般取0.01左右即可，太小了会不收敛。算声子谱应该是和smearing与否没有本质关系。
(著名写手)
散金: 2790
在线: 346.9小时
虫号: 117764
注册: 性别: GG专业: 原子和分子物理
报错是有可能的，如果是这样，你只能分开计算了
(著名写手)
散金: 2790
在线: 346.9小时
虫号: 117764
注册: 性别: GG专业: 原子和分子物理
引用回帖:: Originally posted by 刘仕晨 at
谢谢老师，我刚试了下， 用smearing计算声子，同时设置epsil=.true.确实报错了，显示no elec. field with metals，还是这个意思。
老师您说分开算是什么意思？分开计算波恩有效电荷和介电常数么，如果是这样，那就 ... 我有点不懂啊！反铁磁一定要金属设置吗？我记得不用啊！你把你的设置贴？来我看看
(著名写手)
散金: 2790
在线: 346.9小时
虫号: 117764
注册: 性别: GG专业: 原子和分子物理
【答案】应助回帖
引用回帖:: Originally posted by 刘仕晨 at
不是的，老师，可能我表达的不清楚，我的目的是要计算一个绝缘体系的声子谱，其磁结构为G-type的反铁磁，并且需要计算波恩有效电荷和介电张量来进行对声子谱修正，这样的情况下，设置occupation为smearing的话，ph ... 金属和磁性本质上是两个独立的维度！对绝缘体同样可以设置反铁磁！我记得例子中有关于FeO的例子，你找出来看看应该就明白了！
相关版块跳转
第一性原理
我要订阅楼主
的主题更新
小木虫,学术科研互动社区,为中国学术科研免费提供动力
违规贴举报删除请联系客服电话： 邮箱：（全天候） 或者 QQ:
广告投放与宣传请联系 李想 QQ：
QQ：&&邮箱：
Copyright & 2001-, All Rights Reserved. 小木虫 版权所有
浏览器进程
打开微信扫一扫
随时随地聊科研 上传我的文档
 下载
 收藏
免责声明：1，本账号发布文档部分来源于互联网，仅用于技术分享交流用，版权为原作者所有。 2，文档内容部分来自网络意见，与本账号立场无关。涉及政治言论一律相应删除，请大家监督。 3，如有侵犯原您的版权,请提出指正，我们将立即删除相关资料有其它问题也欢迎与本人联系，谢谢！
 下载此文档
反铁磁性和亚铁磁性
下载积分：2000
内容提示：反铁磁性和亚铁磁性
文档格式：PDF|
浏览次数：52|
上传日期： 20:11:11|
文档星级：
全文阅读已结束，如果下载本文需要使用
 2000 积分
下载此文档
该用户还上传了这些文档
反铁磁性和亚铁磁性
官方公共微信反铁磁性_中国百科网
    
铁磁性ferromagnetism过渡族金属（如铁）及它们的合金和化合物所具有的磁性。研究简史铁磁理论的奠基者，法国物理学家P.-E.外斯于1907年提出了铁磁现象的唯象理论。他假定铁磁体内部存在强大的“分子场”，即使无外磁场，也能使内部...
。 　　BiFeO3(BFO)是一种具有菱形扭曲的钙钛矿结构的多铁性材料，室温下同时具有两种结构有序，即铁电有序（Tc=1103 K）和G 型反铁磁有序（TN = 643 K），是少数室温下同时具有铁电性和铁磁性的多铁性材料之一。研究结果显示薄膜形态的...
收录时间:日 12:53:48 来 源：未知作者:反铁磁性,物质,自旋,
上一篇： &(&&)
创建分享人
Copyright by ；All rights reserved.