Al的声子谱计算

Al的声子谱计算

hhwj340

本文介绍计算Al的声子的计算的过程。使用的软件是ABINIT软件包(version 7.4.2)。使用的赝势是模守恒赝势(TM型),交换关联能使用 GGA近似(PBE型)。计算声子时，应用的理论方法是线性响应理论。计算的具体参数见输入文件。

Step 1:对Al的结构进行驰豫,输入文件如下：

计算得到晶格常数为：7.6500879384(Bohr)

Step2:获得计算声子谱所有所需K点

将Step1中的晶格常数由7.6500879384代替，设置nstep=1，及nline=1,ngkpt12 12 12,简单运行之后即可以得到所有所需K点。K 点可以在生成的out,或者DDB文件中截取，注意如果生成的K点多余50个，out文件中只列出前50个，所以最后在DDB文件中截取。

Step3:计算所有所需Ｋ点的DDB文件，输入文件如下：

计算可以获得所有的DDB文件。Step 4:合并DDB文件

Step 5:使用ABINIT软件包中的anaddb计算所需点的声子谱

Step 6: 使用软件包中的band2eps来绘图

最后得到Al的声子谱：