用于直接探究核模型的 KS-DFT 分子动力学的 Sesame 风格分解

构建状态方程的一种常见范式是将热力学分解为三项的叠加：静态晶格冷曲线、原子核热运动的贡献以及电子热激发的贡献。虽然晶体的统计力学模型为固相中的核贡献提供了可处理的框架，但对于熔化温度以上（$C_v^{(\text{nuc})}\approx 3R$）的核贡献及其如何过渡到高温极限（$C_v^{(\text{nuc})} \sim \frac{3}{2}R$）知之甚少。在这项工作中，我们描述了一种从量子分子动力学 (QMD) 中提取热核贡献和热电子贡献的算法。随后，我们使用 VASP QMD 软件包探测正常密度下液态铝的热核行为，并将结果与半经验模型——Johnson 通用模型、Chisolm 高温液体模型和 CRIS 模型——进行比较。