TDDFT计算软件Octopus学习笔记(五)：介电函数和吸收谱(ZnO)

然后进行时间传播的计算

计算结束后，相应信息储存在 ./td.general 文件夹下的 gauge_field 文件中。 Octopus自带小程序来处理这个文件，只需要在含有inp文件的目录下(即 ./ 目录)运行 oct-dielectric-function 命令，小程序自动运行，就会在 ./td.general 文件夹中生成 “dielectric_function”， “inverse_dielectric_function” 和 “chi” 这3个文件。 其中 “dielectric_function” 文件即包含有介电函数的信息。其文件内部结构如下：

共7列，分别是：能量，介电函数在x方向的实部，虚部，在y方向实部，虚部，z方向实部，虚部。
使用画图小软件gnuplot，简单画图，即可看x方向虚部。

还可以通过介电函数实部和虚部来计算吸收谱：

注1：时间步长的确定

时间步长的选择十分重要，太大的步长会导致计算发散，太小的步长又会浪费计算时间。合理的选择思路如下：在不施加外界微扰的情况下进行时间传播计算，传播两百步左右，由于未施加扰动，系统始终处于基态，所以在第一步的总能量和第二百步的总能量应该相同，但是实际计算中，如果时间步长太大，起点和终点的总能量会略有不同，更大的时间步长甚至会导致总能量发散，出现NaN字样，即 Not A Number。但是，起点和终点的总能量调节至完全相同需要很小的步长，会浪费计算时间，所以，选择一个好的步长是十分重要的。