| 查看: 5218 | 回复: 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 |
回复此楼» 猜你喜欢
研究生在毕业答辩时挂了，遗憾
已经有13人回复
-
24年博士招生
已经有10人回复
-
2024博士招生
已经有4人回复
-
各位同学能否分享一下实验室的学生劳务发放标准呀？
已经有18人回复
-
每到中夜，情难自抑
已经有56人回复
-
论大家对6070后普通教授导师的看法
已经有10人回复
-
九部门发文：不得将专利授权数量作为人才评价、项目评审、职称评定、高校评价等的条件
已经有19人回复
-
没收到邮件
已经有7人回复
-
24or25材料专业申博
已经有11人回复
-
B口人才项目
已经有19人回复
wfbgc
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1802.1
- 散金: 10
- 红花: 2
- 帖子: 220
- 在线: 172小时
- 虫号: 448198
- 注册: 2007-11-01
- 性别: GG
- 专业: 原子和分子物理
12楼2023-04-20 12:35:19
123dingyc
新虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 16121.1
- 红花: 44
- 帖子: 6584
- 在线: 571小时
- 虫号: 2516270
- 注册: 2013-06-21
- 专业: 凝聚态物性I:结构、力学和
2楼2019-02-19 22:52:08
shadowlisong
铁杆木虫 (正式写手)
- 应助: 39 (小学生)
- 金币: 5572.6
- 散金: 10
- 红花: 11
- 帖子: 562
- 在线: 248.9小时
- 虫号: 3433738
- 注册: 2014-09-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
dt97511
3楼2019-02-20 15:40:50
4楼2019-02-20 15:58:25
5