| 查看: 5221 | 回复: 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 |
回复此楼» 收录本帖的淘贴专辑推荐
 缺陷计算 |
» 猜你喜欢
每天学术时间不能保证，能保证的只有：
已经有9人回复
-
化学B02口青基 代表作都是什么水平的？向大佬求助
已经有7人回复
-
研0二导师分到新来的博士后靠谱吗
已经有7人回复
-
研究生在毕业答辩时挂了，遗憾
已经有13人回复
-
24年博士招生
已经有10人回复
-
2024博士招生
已经有4人回复
-
各位同学能否分享一下实验室的学生劳务发放标准呀？
已经有18人回复
-
每到中夜，情难自抑
已经有56人回复
-
论大家对6070后普通教授导师的看法
已经有10人回复
-
九部门发文：不得将专利授权数量作为人才评价、项目评审、职称评定、高校评价等的条件
已经有19人回复
10楼2020-05-20 15:14:06
6