[論文レビュー] 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₂）の熱力学的応答関数およびウィドム線を調査し、異なる性質（例：比熱、圧縮率）の極大値が密度-温度空間で急速に乖離することを明らかにした。これは、単一のウィドム線が物理的意味を持つとは限らないことを示唆している。さらに、定体比熱に基づきCO₂のフレンケル線を提案し、それが対数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₂における熱力学的応答関数（比熱、圧縮率、体積膨張率、密度フラクチュエーション）の挙動を調査すること。
• 等圧線、等温線、定体積線などの異なる熱力学的経路に沿って、ウィドム線（応答関数の極大値の軌跡）の位置を特定すること。
• 定体積比熱に基づくフレンケル線の存在と位置を検討し、剛体的・非剛体的流体ダイナミクスを区別すること。
• これらの発見が、特に金星の初期大気に関する惑星科学に与える影響を評価すること。
• 実験データとの比較を通じてシミュレーションモデルの妥当性を検証し、複雑な流体におけるウィドム線概念の物理的整合性を評価すること。

提案手法

• 実験的液体-気相共存状態および臨界点を再現するように最適化された剛体分子モデルポテンシャルを用いた分子動力学（MD）シミュレーションを実施。
• 超臨界領域における等圧線、等温線、定体積線に沿って、状態方程式および応答関数（cₚ、αₚ、κₜ、χₜ）を計算。
• ρ–T平面およびP–T平面において、応答関数の極大値を特定することでウィドム線を同定し、異なる熱力学的経路間での結果を比較。
• 定体積比熱基準（cᵥ = 3.5）を用いてフレンケル線の位置を特定し、振動的ダイナミクスを示す剛体的流体と非振動的ダイナミクスを示す非剛体的流体の領域を区別。
• 速度自己相関関数を用いて動的転移の同定を支援したが、主にcᵥを用いてフレンケル線を推定。
• 実験データとの比較によりMD結果を妥当性検証し、状態方程式および臨界パラメータにおいて優れた一致を示した。

実験結果

リサーチクエスチョン

• RQ1 超臨界CO₂における異なる熱力学的応答関数（例：cₚ、κₜ、αₚ）の極大値は、単一のウィドム線に一致するのか、それとも乖離するのか？
• RQ2 ウィドム線極大値の位置は、熱力学的経路（等温線、等圧線、定体積線）に依存するのか？その結果、ウィドム線の物理的意味にどのような含意があるか？
• RQ3 定体積比熱基準（cᵥ = 3.5）を用いてCO₂にフレンケル線を信頼性高く同定できるか？また、その流体の微視的ダイナミクスとはどのように関係するか？
• RQ4 フレンケル線は、特に金星の初期大気をモデル化するにあたり、どのような意義を有するか？
• RQ5 CO₂の相図に見られる複数の相（分子相、ポリマー相、金属相）が、複数のウィドム線を生じるのか？それらはどのように区別されるか？

主な発見

• 超臨界CO₂における異なる熱力学的応答関数（例：cₚ、αₚ、κₜ、χₜ）の極大値は、単一のウィドム線に集中せず、ρ–T空間で急速に乖離し、T = 320 K（臨界点より4%上）で密度差が12%を超える。
• 比熱と等圧体積膨張率の極大値は温度上昇に伴い高密度側にシフトするが、等温圧縮率と密度フラクチュエーションの極大値は低密度側にシフトするため、定性的に異なる挙動を示す。
• 同じ関数（例：c₅）の極大値でも、等圧線と等温線に沿って計算すると顕著に異なるため、経路依存性が強く、普遍的なウィドム線の概念に疑問を呈する。
• cᵥ = 3.5基準によるCO₂のフレンケル線は、低温域で比熱の超臨界極大値に起因する曲がりを示すが、高圧・高温領域では対数P–T座標において融点線にほぼ平行となる。
• フレンケル線の計算結果から、初期の金星（最大数10MPaの圧力）では剛体的CO₂流体が存在した可能性があり、その後非剛体的状態に移行したとされ、惑星表面の進化に影響を与えた可能性がある。
• CO₂のMDモデルは、状態方程式および臨界パラメータにおいて実験データと優れた一致を示し、超臨界流体の挙動を研究するにあたりシミュレーション手法の妥当性が裏付けられた。