Phonopy-Spectroscopy 温度线宽实例

温度可通过声子散射改变线宽。这里使用工具作者公开的石英配套数据，在集群上分别运行100、300、600 K的光谱处理，演示读取匹配线宽的实际步骤。

下载完整光谱实例包。在新目录解压后操作；包内包含真实输入、已有DFT响应、25个操作条件、原生输出、数据来源说明和绘图脚本。

1. 哪些数据必须配套#

石英目录包含同一结构的 mesh.yaml、BORN、irreps.yaml、40份真实DFT位移响应及Γ点Phono3py线宽。原始线宽文件提供0–1000 K、每10 K的数据。

本例源码和输入均固定到提交 316fbf4。ZnS和Si是另外的体系，不与石英Born电荷、模式或线宽交叉使用。

2. 先准备接口线宽#

bash prepare.sh

这一步生成明确标注的FWHM接口文件，保留原始Phono3py半峰宽；转换约定见展宽与单位。

3. 比较红外#

bash run.sh quartz-ir-T100
bash run.sh quartz-ir-T300
bash run.sh quartz-ir-T600

底层参数为 --linewidth-hdf5 linewidth-FWHM.hdf5 --linewidth-temperature 300，温度须存在于原始数据中。

4. 比较拉曼#

bash run.sh quartz-raman-T300
bash run.sh quartz-raman-T100
bash run.sh quartz-raman-T600

300 K入口先用原生 -d生成40个位移结构，逐个核对其与原始DFT输出的几何，再执行 -r和 -p。后两个入口复用同一份真实响应，只改变线宽温度。

所有温度都传递同一份 irreps.yaml，以保持简并分组一致。峰表中E类简并模式的强度相加、频率与线宽取组内平均；峰表行数不一定等于全部3N条模式数。

5. 读原始宽度随温度变化#

IR-PeakTable.csv与 Raman-PeakTable.csv给出频率、不可约表示、活动量和FWHM。对应关系来自同一套模式索引，不按“频率看起来差不多”就任意拼接。

这里固定谐波频率和响应张量，只使用给定温度的声子散射线宽。没有额外加入热膨胀、温度频移、完整Stokes玻色因子或激光频率因子；温度对照按这些实际输入范围理解。

原始数据与计算背景见固定上游石英示例。