VASP 声子寿命、平均自由程与边界散射
从Phono3py实算结果拆解热导贡献，并展示简单漫散射边界近似。
阅读本篇目录
总热导率之外，我们常想知道“哪些声子在传热”“主要平均自由程有多长”。本例继续使用硅的三阶力常数与RTA热导算例，分析300 K的模分辨数据。
1. 先完成三阶力常数与热导计算#
原算例采用16原子超胞、39个实际位移力目录、4 Å成对截断，以及9³q网格。下载包中kappa-m999.hdf5来自该次Phono3py计算，包含frequency、gamma、group_velocity、mode_kappa和weight。
这是小超胞工作流，q网格也尚未充分收敛。以下进一步分析继承相同近似，不因为输出更多列就获得更高精度。
2. 正确换算寿命与平均自由程#
Phono3py文件中的gamma使用普通频率THz。采用τ(ps)=1/[4π·gamma(THz)]，不要漏掉2π。零gamma可能表示未计算的Gamma声学模，不能直接除零。本例剔除3个此类模，并确认它们没有遗漏有限热导贡献。
群速度数值从THz·Å乘100得到m/s，平均自由程Λ=|v|τ。下载脚本按所有有效模计算，并将模热导按Λ从小到大累加。
mode_kappa已经包含对应对称性计数；本例按其约定求和后除以weight总和729，不再额外乘一次weight。模贡献合计得到κxx=96.9848 W/(m·K)，与文件内总热导一致。单位与求和方式见Phono3py输出说明。
3. 本次累计结果#
在这套离散网格数据上，累计约一半κxx所对应的平均自由程约134.5 nm。这是该模型的模分辨结果，不是实验测得的纳米结构临界尺寸；加密q网格会影响长平均自由程部分。
4. 加入一个简单边界散射模型#
用有效漫散射长度L，设1/τ总=1/τ三声子+|v|/L，即Matthiessen近似。每个模的原热导乘1/(1+Λ/L)，得到下面的示范：
| 有效长度L / nm | κxx / W m⁻¹ K⁻¹ |
|---|---|
| 10 | 9.718 |
| 30 | 21.050 |
| 100 | 40.265 |
| 300 | 59.797 |
| 1000 | 77.624 |
| 3000 | 88.249 |
L越小，边界散射越强；L趋于无穷时回到原本体值。这里的L是模型有效长度，不能直接等同于任意薄膜厚度或晶粒尺寸。未考虑镜面反射、取向相关边界、同位素或缺陷散射，也没有真正建立纳米器件原子结构。
5. 如何得到自己的图#
先按原热导教程生成自己的HDF5，将脚本中的文件名和温度改成对应结果，再执行sbatch run-post.sh。核对模贡献能重新加回总热导后，再画寿命分布、累计热导和L依赖；不要只画累计曲线而跳过这一步。
环境与下载#
本页使用鸢算科技集群的实际计算数据。VASP 6.6.1、PBE PAW64、ENCUT=550 eV、EDIFF=1E-6；优化自由原子力目标0.02 eV/Å。GPU程序使用NVHPC 25.3、CUDA 12.8、HPC-X 2.22.1和MKL 2023.2构建。声子后处理使用Phonopy 4.8.3、Phono3py 4.8.2；结构和网格读取使用ASE 3.29.0。
附件README说明各目录依赖。准备有权使用的POTCAR，按POSCAR元素顺序组合，再用目录内sbatch run.sh提交。WAVECAR、CHGCAR需由自己的前序计算产生，不随公开包分发。JSON及图片是本次参考数据。