[论文解读] Thermodynamics and Widom lines in supercritical carbon dioxide

D. Fomin, V. N. Ryzhov+3 more|arXiv (Cornell University)|Nov 25, 2014

Phase Equilibria and Thermodynamics被引用 7

一句话总结

本研究利用分子动力学模拟研究了超临界二氧化碳（CO₂）的热力学响应函数及Widom线，发现不同性质（如比热容、可压缩性）的极大值在密度-温度空间中迅速发散，削弱了单一Widom线的物理意义。研究进一步基于等容比热容提出CO₂的Frenkel线，发现其在log P–T图中与熔化线平行，暗示早期金星可能存在刚性CO₂流体相。

ABSTRACT

Behavior of supercritical fluids attracts a lot of attention nowadays. It is important both from the point of view of fundamental science and technological applications. However, up to now the progress in the field is rather moderate. In this article we report a computational study of supercritical carbon dioxide which is one of the most important fluids for chemical industry. We study the response functions of $CO_2$ in supercritical regime and calculate the locations of their maxima (Widom lines). We also make preliminary calculations of the line of crossover of microscopic dynamics of particles (Frenkel line). The conclusions on the Frenkel line location can be applied to study of the atmosphere of Venus.

研究动机与目标

• 研究超临界CO₂中热力学响应函数（比热容、可压缩性、膨胀系数、密度涨落）的行为。
• 确定Widom线的位置——即响应函数极大值的轨迹——在不同热力学路径（等压线、等温线、等容线）上的分布。
• 探索基于等容比热容的CO₂中Frenkel线的存在性与位置，以区分刚性与非刚性流体动力学。
• 评估这些发现对行星科学的影响，特别是对金星早期大气的影响。
• 将模拟模型与实验数据对比，评估Widom线概念在复杂流体中的物理一致性。

提出的方法

• 采用针对CO₂优化的刚性分子模型势能函数，通过分子动力学（MD）模拟，重现液-气共存及临界点的实验行为。
• 在超临界区域内，沿不同等压线、等温线和等容线计算状态方程及响应函数（cₚ, αₚ, κₜ, χₜ）。
• 通过在ρ–T和P–T平面中定位响应函数极大值，识别Widom线，并比较不同热力学轨迹下的结果。
• 应用等容比热容判据（cᵥ = 3.5）确定Frenkel线，以区分具有振荡动力学的刚性流体与非振荡动力学的非刚性流体区域。
• 利用速度自相关函数支持动态转变的识别，但主要依赖cᵥ进行Frenkel线的估算。
• 将MD结果与实验数据对比，结果显示状态方程和临界参数具有极佳一致性。

实验结果

研究问题

• RQ1 在超临界CO₂中，不同热力学响应函数（如cₚ, κₜ, αₚ）的极大值是否沿单一Widom线对齐，还是迅速发散？
• RQ2 Widom线极大值的位置如何依赖于热力学路径（等温线、等压线、等容线）？这对Widom线的物理意义有何启示？
• RQ3 能否可靠地利用等容比热容判据（cᵥ = 3.5）在CO₂中识别Frenkel线？其与流体微观动力学的关系如何？
• RQ4 Frenkel线对行星科学有何意义，特别是在模拟金星早期大气方面？
• RQ5 CO₂相图中存在多种相（分子相、聚合相、金属相）是否导致多个Widom线？它们如何被区分？

主要发现

• 超临界CO₂中不同热力学响应函数（如cₚ, αₚ, κₜ, χₜ）的极大值并未汇聚成单一Widom线，而是在ρ–T空间中迅速发散，且在T = 320 K（高于临界点4%）时，密度差异超过12%。
• 比热容与等压膨胀系数的极大值随温度升高向更高密度移动，而等温可压缩性与密度涨落的极大值则向更低密度移动，表明其行为具有定性差异。
• 即使同一函数（如cₚ）的极大值，在等压线与等温线路径上也存在显著差异，表明其具有路径依赖性，从而削弱了普遍Widom线概念的合理性。
• 基于cᵥ = 3.5判据估算的CO₂ Frenkel线在低温区域因比热容的超临界极大值而出现弯曲，但在更高压力与温度下，其在log P–T坐标中趋于与熔化线平行。
• Frenkel线计算表明，早期金星（压力达数十MPa）很可能存在刚性CO₂流体，随后发生向非刚性状态的转变，可能对行星表面演化产生影响。
• CO₂的MD模型在状态方程与临界参数方面与实验数据高度一致，验证了该模拟方法在研究超临界流体行为中的有效性。