SHG无法通过 VASP 直接计算得出。要获得这些系数，需要计算极化对电场的二阶导数。需要专门编译支持SHG计算的Wannier90版本，通过构建最大局域化Wannier函数，在Wannier表象下计算Berry相位相关的响应，来获得非线性光性质。
本文参考：
https://www.cnblogs.com/jiaming121/p/19454380
支持SHG版本的wannier90下载地址
https://github.com/peio-gg/wannier90/tree/eps_shg_mb
wannier90编译
下载解压后，进入文件夹
cp config/make.inc.ifort make.inc make&& make lib在完成VASP+wannier90的计算并拟合好紧束缚模型的能带结构后，在wannier90.win里加入
berry = trueberry_task=shg #sc：位移电流shift current，shg：二次谐波second harmonic generationberry_kmesh = 15151 #需要测试收敛性，一般是scf密度的数倍#fermi_energy=-3.65 #从OUTCAR读取
kubo_freq_min= 0.0 #光频kubo_freq_max= 2.0kubo_freq_step= 0.02sc_phase_conv=2kubo_adpt_smr=falsesmr_fixed_en_width=0.025#delta函数展宽然后并行执行
mpirun -np xx ~/software/wannier90-3.1.0/wannier90-eps_shg_mb/postw90.x wannier90 |tail -f wannier90.wpout
位移电流
