tip：以下归类的计算都在scf计算以后读取自洽计算的CHGCAR或者WAVCAR进行进一步计算
-
光学性质计算
关键参数
LOPTICS = .TRUE. # 开启光学性质计算
CSHIFT = 0.1 # Kramers-Kronig变换中的虚部偏移（典型值0.1-0.2 eV）
NBANDS = 1000 # 增加能带数量（建议2-3倍于默认值，确保空带足够）
NEDOS = 2000 # 态密度网格点数，控制能谱分辨率
其他
ALGO = Normal # 推荐Normal或Exact
ISMEAR = 0 # 半导体用Gaussian平滑，金属使用1，如果半导体使用0报错则换1
SIGMA = 0.05 # 平滑宽度
LPEAD = .TRUE. # 可选，计算动量矩阵元
相关计算公式后处理
工具：使用vasprun.xml提取介电函数数据，可通过脚本（如p4vasp、VaspKit或Python脚本）绘制吸收光谱、折射率等。
本人使用vaspkit，流程vaspkit------711，会先生成复介电函数实部和虚部文件（real.in和imge.in），然后再从介电函数中提取吸收系数、折射率、反射率、损失曲线的数据。
###################################################################################################
- 使用DFT+U
参数设置
以Gd₂O₃为例，适用性：DFT+U主要用于Gd₂O₃中Gd的4f电子，O的2p电子通常不加U。
LDAU = .TRUE. # 开启DFT+U计算
LDAUTYPE = 2 # 使用Dudarev方法（推荐，Ueff = U - J）
LDAUL = 2 # 应用U的轨道类型（Ti的3d电子，2表示d轨道）
LDAUU = 5.0 0.0 # U值，(Gd: 5.0 eV, O: 0.0 eV，单位eV）
LDAUJ = 0.0 0.0 # J值（Gd: 0.0 eV, O: 0.0 eV，单位eV）
LMAXMIX = 6 # 控制f（d）电子的混合（d轨道用4，f轨道用6）
###################################################################################################
3.其他
LUSE_VDW = .TRUE. # 开启 VASP 内部实现的范德华 (Van der Waals, vdW) 密度泛函修正，考虑长程范德华弱相互作用
BPARAM = 15.7 # 这个参数用于设置 rVV10 范德华修正中的一个经验参数 B
METAGGA=SCAN # SCAN (Strongly Constrained and Appropriately Normed) 泛函是一种先进的 Meta-GGA 泛函
LASPH = .TRUE. # 启用非球形截断，对于过渡金属、稀土元素（具有 d 或 f 电子，其电子云分布通常是非球形的）以及具有共价键合的体系，启用 LASPH = .TRUE. 可以提高计算的准确性