Si笼形结构声子与热导率DFT计算
第一性原理；声子谱；晶格热导率；Si笼形结构；Phonopy
一、文章简要介绍
热电材料可将废热直接转化为电能，其性能由无量纲优值ZT决定，降低晶格热导率kappa_latt是提升ZT的关键途径。Si笼形(clathrate)结构是“声子玻璃-电子晶体”概念的典型代表：其开放的框架结构能大幅抑制热输运，同时保持半导体材料的电学性能，且Si元素廉价无毒、储量丰富。英国曼彻斯特大学Skelton团队在Molecules期刊上发表了系统性的第一性原理研究(DOI:10.3390/molecules27196431)，利用VASP密度泛函理论结合Phonopy/Phono3py，全面计算了金刚石Si及五种框架结构的声子动力学与热输运性质。
该研究构建的六种结构分别为：金刚石Si(d-Si)、Type-II笼形(C-II)、oC24框架、Type-I笼形(K-II/C-I)、Type-VI笼形(K-V/C-VI)和Type-V笼形(K-VII/C-V)。采用PBEsol泛函、PAW赝势、500 eV截断能，所有结构均优化至力收敛10^{-2} eV/A。计算了完整声子色散、红外/拉曼光谱、热力学稳定性以及三阶力常数下的晶格热导率，并与实验数据进行了严格对比验证。
二、仿真步骤
第一步：结构优化。从Materials Project和ICSD数据库获取六种Si结构的初始构型，使用spglib进行对称性分析后，在VASP中进行全优化。采用PBEsol泛函描述电子交换关联能，PAW赝势处理Si的3s/3p价电子，截断能500 eV。能量收敛至10^{-8} eV，原子力收敛至10^{-2} eV/A，k点网格密度从金刚石Si的6x6x6到C-II的4x4x4不等。
第二步：声子谱计算。利用Phonopy通过有限位移法构建超胞（金刚石Si为3x3x3共216原子，oC24为4x4x2共384原子），计算二阶力常数Phi(2)。在规则q点网格上插值得到声子色散曲线和态密度曲线，Gaussian展宽0.064 THz。第三步：光谱模拟。使用DFPT方法计算Gamma点的Born有效电荷张量和极化率导数，分别得到红外活性和平均拉曼活性，生成模拟红外/拉曼光谱，峰宽由MBPT三声子自能虚部计算。
第四步：热导率计算。利用Phono3py计算三阶力常数Phi(3)（超胞位移步长3x10^{-2} A），在单模弛豫时间近似(SM-RTA)下计算声子线宽和晶格热导率kappa_latt。金刚石Si的q点网格为36x36x36，其他结构为9x9x9至16x16x16。第五步：定量分解。应用CRTA模型将kappa_latt分解为谐波项（群速度平方加权）和平均寿命项，定量归因各结构热导率差异的微观来源。
三、关键结果与发现
声子谱分析揭示，C-II结构的谱密度显著向低频偏移（图3累积曲线），其化学键明显弱于其他结构。红外与拉曼模拟光谱为实验物相鉴定提供了参考标准：金刚石Si仅有一个拉曼峰在508 cm^{-1}（与实验值520 cm^{-1}误差约2%），而C-II和K-II/C-I因原胞原子数多（102~138个），光谱特征丰富得多。
热力学分析表明五种框架结构均为亚稳相，相对金刚石Si的晶格能差从K-VII/C-V的6.82 kJ/mol/atom到C-II的32.6 kJ/mol/atom。声子熵对框架结构的稳定化贡献仅0.12~0.92 kJ/mol/atom，远不足以补偿晶格能差距。热导率计算呈现出清晰排序：金刚石Si(136.2) > oC24(40.9) ~ K-II/C-I(43.5) > K-V/C-VI(36.3) > K-VII/C-V(31.2) > C-II(6.33 W/m/K)。C-II的kappa_latt不到金刚石Si的5%——CRTA分解揭示其群速度项(0.46)和平均寿命(13.8 ps)均最低，弱化学键降低了声子群速度，强三声子散射缩短了寿命。金刚石Si的计算值与两组独立实验在150~1000 K范围内的平均误差仅8.5%，验证了方法的可靠性。
