Phono3py 计算累积晶格热导率及绘图

Phono3py是一个主要用python写的声子-声子相互作用相关性质的模拟包，可以基于有限位移算法实现三阶力常数和晶格热导率的计算过程，同时输出包括声速，格林奈森常数，声子寿命和累积晶格热导率等参量。

相关介绍和安装请参考往期推荐。

phono3py快速安装教程

QE+phono3py 热导率计算

VASP+phono3py:快速计算晶格热导率

理论到实践：VASP+Phono3py计算Phonon Lifetime

phonopy 和 phono3py 安装教程Phonopy-Spectroscopy计算材料红外和Raman光谱

本文以最新版phono3py（版本3.10）的example为基础，介绍如何基于力常数计算晶格热导率后，获得累积晶格热导率。

参考链接：

https://github.com/phonopy/phono3py

https://github.com/phonopy/phono3py/tree/develop/example/Si-PBEsol

以硅计算的文件为例。采用

2×2×2超晶胞。

生成有限位移文件后，并通过VASP计算获得力常数fc3.hdf5和fc2.hdf5后，进行热导率计算，采用11×11×11的网格，

phono3py-load --mesh 11 11 11 --br

有关计算结果会被写入到 kappa-m111111.hdf5中，通过读取该文件并生成累积晶格热导率

 phono3py-kaccum kappa-m111111.hdf5 |tee kaccum.dat

kaccum.dat可以

使用gnuplot绘制

gnuplot> p 'kaccum.dat' i 30 u 1:2 w l, 'kaccum.dat' i 30 u 1:8 w l

发现大部分晶格热导率来自6 THz以下的声子模。

而对应的Si的声子谱和声子态密度为：