[论文解读] First principles phonon calculations in materials science

Atsushi Togo, Isao Tanaka|arXiv (Cornell University)|Jun 29, 2015

High-pressure geophysics and materials参考文献 32被引用 11

一句话总结

本论文提出了一套全面的框架，用于在材料科学中基于密度泛函理论（DFT）和有限位移法，在简谐近似下进行从头算声子计算。该研究展示了使用phonopy预测金属Al和Cu的声子色散关系、动力学稳定性和相变路径的应用，结果表明特定波矢（如Γ、M、X）处的声子不稳定性会驱动结构转变，且计算结果与实验测得的晶格热导率和相变能垒高度一致。

ABSTRACT

Phonon plays essential roles in dynamical behaviors and thermal properties, which are central topics in fundamental issues of materials science. The importance of first principles phonon calculations cannot be overly emphasized. Phonopy is an open source code for such calculations launched by the present authors, which has been world-widely used. Here we demonstrate phonon properties with fundamental equations and show examples how the phonon calculations are applied in materials science.

研究动机与目标

• 建立一个稳健、开源的从头算声子计算框架，适用于多种材料。
• 实现对晶体材料中声子色散关系、动力学稳定性和相变路径的精确预测。
• 通过群论和能量势垒分析，将声子不稳定性与结构转变（特别是Al和Cu等金属）联系起来。
• 通过非谐声子-声子相互作用计算晶格热导率，并与实验数据进行验证。
• 阐明声子在理解超导性、力学性能和固体中热输运行为中的作用。

提出的方法

• 采用投影缀加平面波（PAW）方法和广义梯度近似（PBE）在VASP中进行电子结构计算。
• 使用30×30×30的k点网格和300 eV的平面波截断能，采用0.2 eV的Methfessel-Paxton展宽以实现收敛。
• 采用有限位移法，构建3×3×3超胞（含108个原子），位移量为0.01 Å，以计算力常数。
• 通过求解由二阶力常数和原子质量定义的动力学矩阵D(q)的本征值问题，计算简谐声子模式。
• 利用晶体对称性降低计算成本并提高力常数矩阵构建的数值精度。
• 通过分析高对称点（Γ、M、X、K）处的声子不稳定性，映射相变路径，并在受空间群约束下进行几何优化。

实验结果

研究问题

• RQ1 如何通过从头算声子计算预测Al和Cu等材料中晶体结构的动力学稳定性？
• RQ2 特定波矢（如Γ、M、X）处的声子不稳定性在驱动结构相变中起什么作用？
• RQ3 声子-声子相互作用如何影响晶格热导率？这些效应能否从头算出发精确计算？
• RQ4 过渡金属中BCC、HCP和FCC结构之间的相变在能量上最有利的路径是什么？
• RQ5 声子能带结构与对称性分析能否揭示如Cu中的ω相等动态不稳定但能量上可及的中间相？

主要发现

• 声子计算表明，Cu的简单立方（SC）相在Γ点处动力学不稳定，导致通过对称性降低的中间相向体心立方（BCC）相转变。
• SC-Cu在M点的不稳定性导致向面心立方（FCC）结构转变，该路径在P4/nmm空间群下稳定。
• X点不稳定性导致两条可能的对称性畸变路径通往FCC结构，其最高对称性空间群类型分别为Pmma和Cmcm。
• 此前被忽视的Cu中ω相被识别为关键中间态，具有较低的能垒，解释了BCC–FCC转变的可行性。
• 基于从头算声子-声子相互作用强度计算的晶格热导率与多种化合物的实验测量结果高度一致。
• 研究发现，经典Bain路径并非唯一可行路径；相反，声子不稳定性与对称性约束共同定义了一个更精确的、多路径的相变网络，其中ω相作为关键中间态起核心作用。