[论文解读] Microscopic view of heat capacity of matter: solid, liquid, and gas

Jaeyun Moon, Simon Thébaud+2 more|arXiv (Cornell University)|Oct 12, 2022

Phase Equilibria and Thermodynamics被引用 5

一句话总结

该论文通过从瞬时正则模式（INMs）和速度自相关函数（VACFs）推导出的不稳定性参数，提出了一套统一的微观框架，用于理解固态、液态和气态的热容。该框架表明，液态热容源于类固态声子模式与类气态自由度的叠加，其预测结果与分子动力学模拟的误差在5%以内，为经典物相中所有相态的热力学提供了统一描述。

ABSTRACT

Understanding thermodynamics in liquids at the atomic level is challenging because of strong atomic interactions and lack of symmetry. Recent prior theoretical works have focused on describing heat capacity of liquids in terms of phonon-like excitations but often rely on fitting parameters and ad hoc assumptions. In this work, we perform microscopic analysis on instantaneous normal modes and velocity autocorrelations on molecular dynamics simulations of single element systems over wide ranges of temperature (up to $10^8$ K) and pressure (up to 1 TPa). Our results demonstrate that heat capacity of liquids can be described by a combination of both solid-like and gas-like degrees of freedom, leading to a unified framework to describe heat capacity of all three phases of matter: solid, liquid, and gas.

研究动机与目标

• 解决长期以来因液体具有无序、强关联特性且缺乏对称性而带来的热力学理解难题。
• 通过识别热能载流子的共同微观起源，统一描述固态、液态和气态的热容。
• 开发一种超越以往类声子模型中经验假设或拟合参数的框架。
• 利用涵盖极端热力学条件（高达10⁸ K 和 1 TPa）的宽范围分子动力学模拟，提供液体热容的定量、微观解释。
• 通过揭示液体热载流子的双重固态类与气态类特性，实现对液态金属、热能储存系统及离子导体等材料热性能的工程化与优化。

提出的方法

• 分析瞬时正则模式（INMs），提取势能景观的局域曲率，识别与‘类气态’行为相关的不稳定（虚频）模式。
• 利用分子动力学模拟中的速度自相关函数（VACFs），通过不稳定性参数将光谱贡献分解为类固态（声子类）与类气态成分。
• 定义两个不稳定性参数 IP₁ 和 IP₂，以量化振动光谱中‘类气态’的程度，其中 IP₂ 提供更精确的分解。
• 基于不稳定性参数，将 VACF(ω) 进行光谱分解，划分为类气态（VACF_real(ω)）与类固态（VACF_phonon(ω)）成分。
• 利用基于不稳定性参数的光谱分解，通过公式（4）计算定容热容（C_V），该公式将 C_V 与有效热载流子的积分光谱权重相关联。
• 在多种密度和极端温度（高达10⁸ K）下，对单原子系统（Ar、Si、Fe）进行独立计算的 C_V 值进行验证，结果与预测值高度一致。

实验结果

研究问题

• RQ1 如何统一描述所有物相中热能载流子的物理起源，以解释液体的热容？
• RQ2 类固态声子类与类气态自由度在液体热容中以何种程度共存并共同贡献？
• RQ3 能否从瞬时正则模式和速度自相关函数中提取的不稳定性参数，构建一个无需拟合参数的预测性、定量液体热容框架？
• RQ4 局域势能景观曲率与动力学不稳定性在决定液体中有效热载流子方面起什么作用？
• RQ5 在液体中，类固态与类气态模式的相对贡献如何随温度和密度演化？这一演化能否解释从液态向类气态行为的转变？

主要发现

• 液体的热容可被良好地描述为类固态声子类与类气态自由度的叠加，其中后者源于瞬时正则模式谱中的不稳定（虚频）模式。
• 不稳定性参数 IP₂ 能够稳健地将速度自相关函数分解为类固态与类气态光谱成分，从而实现热容的精确预测。
• 基于不稳定性参数框架预测的热容值，在广泛温度与密度范围内，与独立计算的分子动力学结果相比，误差在5%以内。
• 随着温度升高，虚频模式（表示动力学不稳定性）的数量增加，表明系统向类气态行为转变，而实频模式仍持续存在，反映出残留的类固态特征。
• 径向分布函数（g(r)）随温度升高而逐渐失去结构，证实了从液态向类气态状态的转变，这与热容机制中类气态模式的主导性增强密切相关。
• 该框架成功捕捉了液体的双重特性——保留类固态的结构关联性，同时表现出类气态的动力学无序性，为液体独特的热力学响应提供了微观基础。