GW计算光谱

 很久没更文章了，那就更一篇吧

计算光谱（ optical properties）方法很多，PBE最方便因为耗用资源少时间快，可以看到吸收系数的大趋势，缺点就是误差大，不精确。因此HSE/GW就是大家比较热衷的计算方法了。本章分享一下GW方法计算 optical properties（技艺不精，还望大佬看出错误给点建议）。

提醒：以下计算基于VASP5.4来进行分享，版本不同输入文件的参数要相应进行修改

---

gw计算有单步还有两步计算之分，单步容易出错，所以一般使用两步计算（DFT+GW）。

---

1：DFT计算

首先将需要一个高精度优化的结构，然后设置INCAR.这部分分为两小步:

首先是仅包含少量未占用轨道的标准基态计算 

cp -r opt optic_dft1

修改incar文件

INCAR:

SYSTEM=OPTIC 




ISMEAR=0 #zero 
SIGMA=0.05 #For atoms or molecules, use 0.01
EDIFF = 1E-7
LWAVE  = .TRUE.
LCHARG = .TRUE.

KPOINTS POTCAR POSCAR 保持与优化的相同，提交任务。

然后计算大量未占用的轨道

CP -r optic_dft1 optic_dft2

修改incar

SYSTEM=OPTIC_unoccupied

ALGO   = Exact
NELM   = 1
ISMEAR=0 #zero 
SIGMA=0.05 #For atoms or molecules, use 0.01
LWAVE  = .TRUE.
LCHARG = .TRUE.
LOPTICS = .TRUE. 

NBANDS = 516
NCORE = 4

KPOINTS POTCAR POSCAR 保持与不变，提交任务。

注意需要保持LOPTICS = .TRUE. 。计算DFT时也可采用Hybrid，在上述两步的incar中增加Hybrid的参数即可，但由于耗时，我没有计算，需要的参数设置的方法可以自行百度或者私信我。

---

2：GW计算

实际的GW计算是在这一步完成的，这一步可以根据自己需求选择不同的计算风格，如：G0W0,EVGW0,QPGW0等，我采用G0W0来计算……

在上一步基础上计算：

cp -r optic_dft2 optic_gw

修改INCAR:

SYSTEM=OPTIC_GW

ALGO   = GWO
#LSPECTRAL = .TRUE.
NOMEGA    = 50
NELM   = 1
ISMEAR=0 #zero 
SIGMA=0.05 #For atoms or molecules, use 0.01
LWAVE  = .TRUE.
LCHARG = .TRUE.
LOPTICS = .TRUE. 

NBANDS = 516
NCORE = 4

注意该部分的ALGO设置和vasp版本有关（我是vasp5.4），其他版本需要自行查找。

其他保持不变，提交任务

计算结束后，后处理数据，我采用的vaspkit进行处理

./vaspkit

输710即可得到ABSORPTION_2D.dat文件，使用origin进行绘制即可

绘图如下，对比了PBE和GW两种方法得到的光吸收系数