24小时热门版块排行榜    

查看: 5220  |  回复: 11
【悬赏金币】回答本帖问题,作者dt97511将赠送您 20 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dt97511

银虫 (小有名气)

[求助] vasp中+U进行结构优化并计算分析能带,结果没有改善已有2人参与

大家好,本人最近在学习使用vasp计算钙钛矿的结构优化,由于想分析一下电子结构对结构变化的影响,需要计算能带,但DFT中能带会因为自相互作用误差,能带偏小,所以考虑使用GGA+U的方法对能带进行一定修正(自己服务器计算,无法支持杂化泛函等高计算量方法),但发现我自己在网上摸索学习出来的方法,在能带计算结果上并没有太大改善,现将自己的操作步骤和INCAR等文件写出,希望大家看看哪里有问题。

1.结构优化
首先通过MS建模,导出为cif文件,然后用VESTA打开,保存为vasp格式文件,重命名后作为POSCAR
使用cat命令合并各原子的赝势库文件,得到POTCAR
使用0自动生成Monkhorst-Pack法的4*4*4的k点,得到KPOINTS
参考文献和网络上的设置,调整INCAR如下所示:
SYSTEM=convergence test
ENCUT=450
GGA=PE
ISIF=3
ISMEAR=1
SIGMA=0.05
ISPIN=2
EDIFF=1E-5
EDIFFG=1E-3
ISTART=0
ICHARG=2
IBRION=2
NSW=100
#LCALPOL=.TRUE.

LDAU=.TRUE.
LDAUTYPE=2
LDAUL=3 -1 -1
LDAUU=3.5 0 0
LDAUJ=0.5 0 0
LDAUPRINT=2
在终端输入mpirun -np 4 vasp开始计算,得到CONTCAR

2.能带计算
将结构优化后的CONTCAR改为POSCAR,复制KPOINTS,POTCAR,INCAR,并将INCAR中的ISMEAR=1改为-5,其余均未改变,计算得到CHGCAR
将CHGCAR,POTCAR,POSCAR,INCAR复制,新建一个line模式的K点文件(由于教材中没有详细介绍line模式下各行的参数设置,此处直接使用的教材里的设置,没有根据自己的模型进行改动,不知道有没有影响),如下:
Automatic mesh
10 ! 10 intersections
Line-mode
cart
0 0 0 ! gamma
0 0 1 ! X

0 0 1 ! X
0.5 0 1 ! W

0.5 0 1 ! W
0 0 0 ! gamma
同时,注释掉ISMEAR=-5,ISIF=3,NSW=100。更改ISTART=0,变成1,ICHARG=2变成11,计算后得到vasprun.xml
使用P4VASP打开vasprun.xml文件,得到如下的能带图(为什么这里只有能带图,没有DOS图,我看别人有的是DOS和能带图都有的),再export为dat格式,在origin中打开,得到详细数据,分析发现不管INCAR中有没有下面这段+U的设置(已与POTCAR的原子顺序对应),能带均没有太大改变,请问各位我是不是有哪些步骤没有作对呀
LDAU=.TRUE.
LDAUTYPE=2
LDAUL=3 -1 -1
LDAUU=3.5 0 0
LDAUJ=0.5 0 0
LDAUPRINT=2

vasp中+U进行结构优化并计算分析能带,结果没有改善
TIM截图20190219223303.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shadowlisong

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
+U 的数值要参考之前的文献,如果没有要测试。

还有就是对非强关联体系,没有d f 轨道电子的元素种类,+U效果并不明显。

» 本帖已获得的红花(最新10朵)

Hello,world
3楼2019-02-20 15:40:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shadowlisong

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
dt97511: 金币+10, 有帮助 2019-02-21 16:48:46
引用回帖:
4楼: Originally posted by dt97511 at 2019-02-20 15:58:25
我试了两个常见的U值,都没有太大改善。

我加的是澜(La)元素,我看一本参考书上说:强关联体系是锕系和过渡金属氧化物,并没有提到澜系,不知道是澜系是非强关联体系还是写漏了,但我还是看到有部分文献是+U计 ...

步骤看起来没什么大问题,你用P4VASP做PDOS 可以看分波态密度,加不加U 应该对 d f 轨道的分布形状有影响
Hello,world
5楼2019-02-21 13:16:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dt97511 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 研究生在毕业答辩时挂了,遗憾 +10 spacr2022 2024-05-27 13/650 2024-06-01 17:45 by keyaner23
[基金申请] 九部门发文:不得将专利授权数量作为人才评价、项目评审、职称评定、高校评价等的条件 +16 sjtu2012 2024-05-28 19/950 2024-06-01 15:24 by 天外飞去来
[硕博家园] 实验室太吵闹,无法安静学习,怎么办? +6 utahh 2024-05-31 11/550 2024-06-01 13:27 by utahh
[教师之家] 选择 +6 459582015 2024-05-31 6/300 2024-06-01 07:32 by otani
[教师之家] 中年 (金币+3) +18 459582015 2024-05-28 19/950 2024-06-01 00:41 by 沈婉婷.Girl
[找工作] 找工作如此之难 +7 探123 2024-05-25 7/350 2024-05-31 15:58 by 平南小刀
[教师之家] 在大地上我们只过一生---看完我的阿勒泰上头了好几天,完结那天晚上几乎失眠 +10 瞬息宇宙 2024-05-27 11/550 2024-05-31 15:38 by 烟 火
[论文投稿] 纠结选哪一个期刊,电化学领域 50+8 Freya163 2024-05-28 10/500 2024-05-31 15:09 by wzykobe
[有机交流] 苯酚和炔溴反应成醚 6+3 TONGMEIMEI 2024-05-26 12/600 2024-05-31 14:41 by TONGMEIMEI
[基金申请] 离职后国自然项目不能变更单位,在新单位还能申请新的国自然项目吗 5+3 宜兰 2024-05-31 6/300 2024-05-31 12:05 by felicity6056
[文学芳草园] 对对子啊 +5 天若孤独 2024-05-29 7/350 2024-05-31 09:00 by wjykycg
[基金申请] 入职高校3年发表10+SCI,尽人事听天命 +29 kaoyan250 2024-05-27 40/2000 2024-05-31 08:44 by Xiaolin81
[材料综合] 真空封石英管 北京 +4 dessha 2024-05-29 5/250 2024-05-30 16:40 by mpdfwxgui
[文学芳草园] 物是人非 +4 myrtle 2024-05-30 4/200 2024-05-30 15:05 by mapenggao
[有机交流] 液相纯度高,但产品析不太出来 10+4 cui19236 2024-05-27 9/450 2024-05-30 07:45 by yuanjijoy
[硕博家园] 又想换工作 +15 brightmj 2024-05-27 20/1000 2024-05-29 23:25 by zyqchem
[有机交流] 机理求助 120+4 15147165026 2024-05-26 10/500 2024-05-28 14:42 by 江东闲人
[有机交流] 奇怪的物质 100+4 桃桃PXS 2024-05-27 7/350 2024-05-28 10:22 by 091602
[基金申请] 面上基金会评专家,有回避机制吗? +4 huang1991js 2024-05-27 4/200 2024-05-27 19:08 by 星火12
[硕博家园] 课题做到一半担心做不出来怎么办 +4 mkpooooo 2024-05-27 7/350 2024-05-27 11:56 by kanmiaolucky
信息提示
请填处理意见