| 查看: 5235 | 回复: 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 |
回复此楼» 收录本帖的淘贴专辑推荐
 缺陷计算 |
» 猜你喜欢
化学B02口青基 代表作都是什么水平的？向大佬求助
已经有24人回复
-
看到大家一直在说代表作NC，不明觉厉
已经有15人回复
-
G口yq发了吗
已经有8人回复
-
什么时候会评呢？
已经有6人回复
-
今年香江初审出结果了吗
已经有7人回复
-
函评意见：既然你有中科院****人才项目了，就不要申请面上了。
已经有24人回复
-
河北省基金
已经有8人回复
-
蹲一个6月7号江苏省青基公示
已经有22人回复
-
这样的课题组值得去吗
已经有4人回复
-
年年想打听，年年打听不到。。
已经有27人回复
10楼2020-05-20 15:14:06
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
15