QE (PWSCF) 计算光学性质

用DFT计算固体光学的方法，比较先进的是在GW，BSE方程的基础上加入自能、准粒子相互作用的修正，以及含时密度泛函等方法，在QE（pwscf）中有gwl，yambo，tddfpt 等模块发展相关的方法。今天不准备讨论这些方法，而是介绍一种基于LDA/GGA的通过Kohn-Sham轨道跃迁矩阵元计算介电常数的QE后处理程序epsilon.x，编译了QE的PP模块（make pp）之后，目录e-6.1/PP/src下会生成这个可执行文件。

源文件qe-6.1/PP/src/epsilon.f90给出了输入参数

以及参数的默认值。

输入文件的例子，
&inputpp
outdir='OUTepsilon'
prefix='pwscf'
calculation='eps'
/
&energy_grid
smeartype='gauss'
smear=0.15d0
wmax=30.0d0
nw=1000
shift=0.0d0
/
命名为eps.in，运行epsilon.x < eps.in > eps.out

k点选取时要关掉对称性，用Monkhorst-Pack自动生成k点时要设置
nosym=.true. 和 noinv=.true. 这时所有的k点的权重相等。
epsilon.x的计算现在只支持模守恒赝势，输出4个文件eels_pwscf.dat epsi_pwscf.dat epsr_pwscf.dat ieps_pwscf.dat，分别是eels谱（电子能量损失谱），介电函数的虚部、实部，虚轴上的介电函数。

输出格式如下，例如epsi_pwscf.dat 的内容如下，意义见本帖8楼