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中加入下面这条:

 LEPSILON = .TRUE

按照软件案例,作者可能更喜欢hdf5格式的数据文件,生成force_constants.hdf5文件的命令如下

phonopy --fc vasprun.xml --hdf5

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

phonopy --dim="6 6 3" --readfc --hdf5 --fc-symmetry --mesh="1 1 1" --eigenvectors

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

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

phono3py --dim="2 2 2" --dim_fc2="6 6 3" --fc2 --fc3 -v --br --thm --mesh="48 48 48" --write_gamma --gp=0

Part.4

红外强度计算

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

phonopy-ir --ir_reps --linewidth_hdf5="kappa-m484848-g0.hdf5" --linewidth_temperature=300

Part.5

Raman计算

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

phonopy-raman -d --bands="4 5 6 7 8 9 11 12 13 14 15 17 18 20 21 22 23 25 26 27"

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

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

ALGO = NormalEDIFF = 1E-8ENCUT = 700ISIF = 2ISMEAR = 0LASPH = .TRUE.LCHARG = .FALSE.LEPSILON = .TRUE.LREAL = .FALSE.LWAVE = .FALSE.NSW = 0PREC = AccurateSIGMA = 0.01SYSTEM = SiO2

计算完成后处理OUTCAR文件

phonopy-raman -r OUTCAR.*

最后获得Raman计算结果

phonopy-raman -p --ir-reps --linewidth-hdf5="kappa-m484848-g0.hdf5" --linewidth-temperature=300

愿有所成

IEchoQ

引喻失义   妄自菲薄

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值