[論文レビュー] 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）を用い、不安定性パラメータを介して、固体的（フォノン的）および気体的成分へのスペクトル的分解を実施する。
• 2つの不安定性パラメータ、IP₁およびIP₂を定義し、振動スペクトル内の『気体的性』の度合いを定量化する。IP₂はより正確な分解を提供する。
• 不安定性パラメータに基づき、VACF(ω) を気体的（VACF_real(ω)）および固体的（VACF_phonon(ω)）成分にスペクトル的分解する。
• 式4に示される式を用いて、有効な熱輸送粒子のスペクトル重みの積分と関連づけ、定積比熱（C_V）を不安定性パラメータに基づくスペクトル分解から計算する。
• アルゴン（Ar）、ケイ素（Si）、鉄（Fe）の単原子系を用い、さまざまな密度および極端な温度（最大10⁸ K）の条件下で、独立に計算されたC_Vと比較して予測を検証する。

実験結果

リサーチクエスチョン

• RQ1 すべての物質相（固体・液体・気体）における熱エネルギー輸送粒子の物理的起源を統一的に記述するには、液体の比熱をどのように記述すればよいか？
• RQ2 固体的フォノン的および気体的自由度が、液体の比熱にどの程度共存し、寄与しているか？
• RQ3 瞬間正規モードおよび速度自己相関関数から導出される不安定性パラメータは、フィッティングパラメータを必要とせず、予測的かつ定量的な液体比熱のフレームワークを提供できるか？
• RQ4 局所的なポテンシャルエネルギー障壁の曲率および力学的不安定性が、液体における有効な熱輸送粒子を決定する役割を果たすか？
• RQ5 液体において、温度および密度の変化に伴い、固体的および気体的モードの寄与比はどのように変化するか？この変化は、液体から気体的挙動への遷移を説明できるか？

主な発見

• 液体の比熱は、固体的フォノン的および気体的自由度の重ね合わせとして良好に記述され、後者は瞬間正規モードスペクトル内の不安定（虚数振動数）モードに起因する。
• 不安定性パラメータIP₂は、速度自己相関関数を固体的および気体的スペクトル成分に堅牢に分解でき、比熱の正確な予測を可能にする。
• 不安定性パラメータフレームワークに基づく予測された比熱値は、広範な温度および密度範囲において、独立に計算された分子動力学シミュレーション結果と5%未満の誤差で一致する。
• 温度が上昇するに従い、虚数振動数モード（力学的不安定性を示す）の割合が増加し、気体的挙動への移行が示唆される一方で、実数振動数モードは依然として存在し、残留する固体的特徴を反映している。
• 対分布関数（g(r)）は、温度上昇に伴い構造の喪失が進行することを示しており、液体から気体的状態への遷移が確認されており、これは比熱機構における気体的モードの優位性増加と相関している。
• 本フレームワークは、固体のように構造的相関を保持しながらも、気体のように動的不規則性を示す液体の特異な熱力学的応答の微視的基盤を的確に捉えている。