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 = true
berry_task=shg #sc：位移电流shift current，shg：二次谐波second harmonic generation
berry_kmesh = 15151 #需要测试收敛性，一般是scf密度的数倍
#fermi_energy=-3.65 #从OUTCAR读取

kubo_freq_min= 0.0 #光频
kubo_freq_max= 2.0
kubo_freq_step= 0.02
sc_phase_conv=2
kubo_adpt_smr=false
smr_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
图片

 

图片

位移电流

图片