grep E-fermi OUTCAR 即可提取出来。
2.2 静电势能的计算:
通过在 INCAR 中添加: LVHAR =.TRUE. 这个参数。
加⼊这个参数,计算结束后,VASP 会输出⼀个⽂件: LOCPOT ⽂件。我们可以通过脚本,或者程序对这个⽂件后处理来获取静电势能,怎么处理后⾯会讲的很清楚,⼤家先不要着急⼼慌。
2.3 VASP的计算流程:优化获取稳定的结构
将 CONTCAR 复制成 POSCAR
在 INCAR 中设置单点计算:
⽅法1:设置 NSW=1 或者 0
⽅法2:设置 IBRION = -1
⽅法3:直接删除 INCAR 中 NSW 这⼀⾏,采⽤默认值
function怎么记忆提交任务
注意:
如果第⼀步计算中有保存 WAVECAR,则静态计算的时候可以设置:ISTART = 2 读取⼀下,以便节省计算时间;
如果没有保存,直接算即可。不要在读不读 WAVECAR 这个问题上瞎操⼼;
电⼦步数的默认值是 40,如果你的体系不容易收敛,计算的时候请设置⼀个较⼤的 NELM 值,⽐如 NELM = 500。
3 Cu(111) 表⾯的功函数计算
3.1 准备计算⽂件:
新建⼀个 workfunction 的⽂件夹,将前⾯优化计算的: INCAR KPOINTS POTCAR CONTCAR 和提交任务的脚本复制到这个⽂件夹中。
3.2 修改 INCAR 改为静态(单点)计算
3.3 修改 INCAR,添加 LVHAR = .TRUE. 这⼀⾏1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19System = Cull
ISMEAR = 1
SIMGA = 0.1
ALGO = FAST
ENCUT = 450
EDIFF = 1E-5
IDIPOL =.TRUE.
IDIPOL = 3
NWRITE = 0
LWAVE = .FALSE.
LCHARG = .FALSE.
IBRION = 2
POTIM = 0
EDIFFG = -0.01
NSW = 1
LVHAR = .TRUE.
3.4 将 CONTCAR 重命名为 POSCAR
3.5 提交任务,等待结束。
4 后处理⽅法-1:p4vasp
p4vasp ⼀直是本⼈强烈推荐的VASP计算相关的建模,可视化,查看结构,结果的软件。本节我们讲解⼀下如何通过 p4vasp 对功函数的计算结果可视化。
4.1 导⼊数据Linux 系统下⾯,进⼊计算的⽬录直接敲命令: p4v & 即可 (我们稍微等下 Windows ⽤户的操作,后⾯可视化过程是⼀样的)
Windows 下⾯,⾸先要将计算⽣成的 LOCPOT ⽂件保存到本地电脑上。
注意:有时候这个⽂件可能会很⼤,下载的时候⽐较蛋疼,不想下载的话等待下节的讲解。
打开 p4vasp 并导⼊ LOCPOT ⽂件:file –> loadsystem,然后到 LOCPOT 所在的⽬录,点解导⼊,效果如下
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论