Phonopy-Spectroscopy计算材料红外和Raman光谱

经验分享

        Phonopy- spectroscopy是一个为Phonopy代码添加模拟振动光谱能力的项目。

        它可以从Phonopy或VASP计算红外(IR)强度，和计算Raman活度张量和远离共振近似的标量平均强度。

        参考该软件examples的计算流程，整理了这篇教程。

准备内容

软件安装

声子谱&三阶力常数计算

红外强度计算

Raman计算

Part.1

软件安装

软件安装包可从github下载

https://github.com/JMSkelton/Phonopy-Spectroscopy/

解压后，主要文件如下

首先要将执行命令所在的文件夹Scripts添加到环境变量中，或者根据所在目录直接调用文件夹内命令。

然后将文件夹SpectroscoPy添加到Python环境可调用的库中。因为我这里使用的是Anaconda，所以将该文件夹移动到site-packages文件夹中。

由于在我的使用环境中出现了一些 错

在phonopy.py中添加下图中阴影内容后即可正常使用

Part.2

声子谱&三阶力常数计算

声子谱的计算主要目的是获得二阶力常数文件FORCE_CONSTANTS，而为了计算光学性质，同时需要计算波恩电荷，写出BORN文件。需要在计算声子谱的INCAR中加入下面这条：

生成mash.yaml或mash.hdf5文件：

使用phono3py软件进行三阶力常数文件的计算。具体计算流程可以参考VASP+phono3py:快速计算晶格热导率

生成&Gamma点模式的声子线宽:

Part.4

红外强度计算

这时需要使用之前准备好的BORN文件，生成室温(300 K)线宽的模拟红外光谱和峰值表,命令如下

Part.5

Raman计算

通过irreps.yaml的振动信息，并参考D3点组的字符表，判断可能显示Raman活动的模式，并生成计算文件。

bands的值为需要计算的模式。

计算的INCAR 参考软件的example，如下

计算完成后处理OUTCAR文件

最后获得Raman计算结果

愿有所成

引喻失义   妄自菲薄

文章知识点与官方知识档案匹配，可进一步学习相关知识Python入门技能树首页概览212224 人正在系统学习中