【摘要】 往期推文中我们已经详细介绍了运用直接法来计算声子谱，具体可点击链接进行查看第一性原理||VASP计算声子谱（以SiO2-HP为例）。那么本期将为大家介绍另一种方法，除了直接法还可以应用密度泛函微扰理论（DFPT）来进行计算声子谱。

 

往期推文中我们已经详细介绍了运用直接法来计算声子谱，具体可点击链接进行查看第一性原理||VASP计算声子谱（以SiO2-HP为例）。那么本期将为大家介绍另一种方法，除了直接法还可以应用密度泛函微扰理论（DFPT）来进行计算声子谱。

 

DFPT是直接计算得到force constants文件，在经过后处理得到对应的声子谱。今天，今天我们简单描述通过VASP 结合DFPT怎么计算声子谱，以NaCl为例，详细参考官网（VASP-DFPT & phonopy calculation — Phonopy v.2.12.0），具体操作如下：

 

第一步：对元胞进行高精度的结构优化

 

EDIFF = 1E-8，甚至更高

 

构优化第二步：准备计算的POSACR

 

从POSCAR-unitcell得到完美超胞（2*2*2），在phonopy环境下，运行命令：

 

phonopy -d --dim=”2 2 2” -c POSCAR-unitcell

 

会得到一系列文件，包括SPOSCAR, phonopy_disp.yaml, and POSCAR-{number}，如下所示：

 

 

我们需要用到的只有SPOSCAR，用作接下来计算的POSCAR文件。

 

第三步：设置计算的INCAR

 

设置IBRION=8，其他参数根据实际情况修改。

 

 

第四步：提交计算任务

 

提交计算，计算结束后，会产生vasprun.xml，其包含了Hessian，从而得到力常数。

 

第五步：后处理，得到力常数

 

使用命令：

 

 

对上一步的计算结果进行处理，得到FORCE_CONSTANTS文件。

 

第六步：后处理，得到声子谱

 

准备band.conf文件，如下所示：

 

 

再使用命令，即可得到：

 

 

本文所有内容文字、图片和音视频资料，版权均属科学指南针网站所有，任何媒体、网站或个人未经本网协议授权不得以链接、转贴、截图等任何方式转载。