| 查看: 5137 | 回复: 11 | |||
| 【悬赏金币】回答本帖问题，作者dt97511将赠送您 20 个金币 | |||
[求助]
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  TIM截图20190219223303.jpg |
回复此楼» 收录本帖的淘贴专辑推荐
 缺陷计算 |
» 猜你喜欢
10043以后咋出去啊
已经有3人回复
-
两类问题算是白选了~
已经有13人回复
-
25年博士申请
已经有12人回复
-
博士论文被抄袭
已经有43人回复
-
"颜宁:基础研究应顶天立地"能做到基础研究同时顶天立地的才是牛人
已经有6人回复
-
大学直属学院卸任的副院长退休后还享受副处级的养老待遇吗？
已经有4人回复
-
投稿RSC旗下杂志突然间看不到投稿状态，里面投稿的文章信息不见了，有遇到过的吗？
已经有3人回复
-
求论文投稿期刊推荐
已经有4人回复
-
AC审稿
已经有6人回复
-
聊天
已经有13人回复
shadowlisong
铁杆木虫 (正式写手)
- 应助: 39 (小学生)
- 金币: 5572.6
- 散金: 10
- 红花: 11
- 帖子: 562
- 在线: 248.8小时
- 虫号: 3433738
- 注册: 2014-09-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
dt97511
3楼2019-02-20 15:40:50
shadowlisong
铁杆木虫 (正式写手)
- 应助: 39 (小学生)
- 金币: 5572.6
- 散金: 10
- 红花: 11
- 帖子: 562
- 在线: 248.8小时
- 虫号: 3433738
- 注册: 2014-09-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2019-02-21 13:16:07
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +221 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3479.3
- 散金: 14
- 红花: 222
- 帖子: 2464
- 在线: 1247.8小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
9楼2019-08-23 12:19:41
123dingyc
新虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 16121.1
- 红花: 44
- 帖子: 6584
- 在线: 571小时
- 虫号: 2516270
- 注册: 2013-06-21
- 专业: 凝聚态物性I:结构、力学和
2楼2019-02-19 22:52:08
4楼2019-02-20 15:58:25
勿言
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2938.3
- 散金: 30
- 帖子: 250
- 在线: 38.6小时
- 虫号: 4365924
- 注册: 2016-01-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2019-08-18 13:35:14
xiewuhua
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 1613.7
- 散金: 500
- 红花: 12
- 帖子: 585
- 在线: 362.4小时
- 虫号: 1253926
- 注册: 2011-04-03
- 性别: GG
- 专业: 理论和计算化学
7楼2019-08-19 23:31:07
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25617.7
- 散金: 602
- 红花: 166
- 帖子: 3765
- 在线: 873.7小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
8楼2019-08-22 20:42:32
10楼2020-05-20 15:14:06
5