vasp计算自由能：振动熵贡献与准谐近似的精度路径

vasp计算自由能的需求，往往出现在那些需要预测材料热力学稳定性的课题中。这个项目在研究钛酸锶（SrTiO3）的高温相变行为时，自由能随温度的变化曲线是判断相变温度的核心依据。项目组需要用VASP计算不同温度下的亥姆霍兹自由能，而声子频率的准确获取成为了整条计算链路的咽喉。

vasp计算自由能的基础，是通过声子频率获取振动熵对自由能的贡献。项目组采用VASP结合Phonopy软件的有限位移法，在SrTiO3原胞中施加0.01 Å的原子位移，计算力常数矩阵。2×2×2超胞的声子色散曲线在Gamma点附近出现了明显的虚频，最低频率为-1.2 THz。系统在这里卡了将近三周，项目组反复排查后认定，超胞尺寸不足导致长程力截断，将超胞扩展至4×4×4后，虚频消失，最低声子频率稳定在0.3 THz。

差距不会说谎。2×2×2超胞给出的室温自由能比4×4×4结果偏高0.15 eV/原子，这个偏差足以使相变温度的预测偏移200 K以上。项目组最终确认，对于钙钛矿结构，超胞尺寸至少需要4×4×4才能保证声子频率的收敛，这与Togo在Computer Physics Communications上发表的Phonopy方法论文中的建议一致。

振动自由能的计算公式将声子频率与温度耦合：F_vib = kT Σ ln[2sinh(ℏω/2kT)]。项目组在300 K下计算得到的SrTiO3振动自由能为-0.23 eV/原子，与中子散射实验推导的-0.21 eV/原子偏差控制在10%以内。这个精度在后续的相图构建中发挥了决定性作用。

vasp计算自由能在有限温度下的核心方法，是准谐近似（QHA）。项目组在7个体积点（从58 ų到66 ų）上分别计算声子频率，构建振动自由能对体积的函数F(V,T)。在每个温度下对F(V,T)进行Birch-Murnaghan状态方程拟合，得到该温度下的平衡体积和自由能。

QHA的结果揭示了SrTiO3从立方相到四方相的转变温度。项目组计算的相变温度为105 K，实验值为106 K，吻合度令人振奋。然而，项目组也注意到，在500 K以上，QHA预测的自由能与实验值开始出现系统性偏差，偏差量在700 K时达到0.08 eV/原子。根源在于QHA忽略了非谐声子-声子相互作用，在高温下这一效应不可忽略。

项目组随后引入了热力学积分方法来修正高温区的自由能。通过VASP的分子动力学模块，在NVT系综下采集不同温度的系综平均势能，再以谐振子参考系统为基准进行积分。在700 K下，热力学积分给出的自由能修正量为-0.06 eV/原子，修正后的值与实验数据的偏差缩小至0.02 eV以内。这一方法与Kresse课题组在Physical Review Letters上的工作思路一致。

声子频率计算对VASP参数的敏感性极高。项目组在参数测试中发现，ENCUT从400 eV提升至520 eV后，声子频率的最大变化达到0.15 THz；而EDIFF从1×10⁻⁶降至1×10⁻⁸后，力常数的收敛精度改善了两个数量级。项目组认定，声子计算必须使用`EDIFF=1E-8`和`PREC=Accurate`的组合，任何降低精度的做法都会在声子频率中引入不可控的噪声。

K点密度同样不可忽视。项目组在Gamma中心K点从4×4×4到12×12×12的测试中，发现自由能在8×8×8后收敛至0.005 eV/原子以内。

vasp计算自由能必须正视准谐近似的固有局限。QHA假设声子频率仅依赖于体积，不显式包含温度对声子线宽的影响。对于具有强非谐性的材料（如PbTiO3、GeTe等铁电或相变存储材料），QHA在相变温度附近的预测误差可达0.1-0.2 eV/原子。

此外，电子热激发对自由能的贡献在金属体系中不可忽略。项目组在补充计算铜的自由能时发现，电子熵在1000 K时贡献了-0.03 eV/原子的自由能修正，忽略这一项会使高温自由能系统性偏高。被证明有效的做法是，对金属体系始终开启`ISMEAR=-1`并设置合理的展宽参数，同时通过Mermin自由能公式纳入电子熵贡献。回过头看，将声子QHA与电子热激发分别计算再线性叠加，是在计算效率和物理准确性之间的最优平衡——这套方案在项目组后续的三个金属间化合物热力学研究中均得到了验证。