[論文レビュー] Exact Analytic Second Virial Coefficient for the Lennard-Jones Fluid

Byung Chan Eu|ArXiv.org|Sep 17, 2009

Aerodynamics and Fluid Dynamics Research参考文献 3被引用数 4

ひとこと要約

本稿では、Lennard-Jones流体の第二ビリアル係数の正確な解析的表現を、すべての温度範囲にわたって提示する。変数変換により、元来は調和振動子に類似した力学的構造が明らかになる。解は放物型円柱関数および発展型超幾何関数を用いて表現され、低温および高温極限における正確な挙動を示し、無次元単位でのボイル温度は3.41793であり、文献値と非常に良好に一致する。

ABSTRACT

An exact analytic form for the second virial coefficient, valid for the entire range of temperature, is presented for the Lennard-Jones fluid in this paper. It is derived by making variable transformation that gives rise to the Hamiltonian mimicking a harmonic oscillator-like dynamics for negative energy. It is given in terms of parabolic cylinder functions or confluent hypergeometric functions. Exact limiting laws for the second virial coefficient in the limits of low and high temperatuers are also deduced for the Lennard-Jones fluid.

研究の動機と目的

• 全温度範囲にわたって有効なLennard-Jones流体の第二ビリアル係数の正確な解析的表現を導出すること。
• 第二ビリアル係数の積分形の背後にある力学的構造を解明し、隠れた調和振動子に類似した系を明らかにすること。
• 温度がゼロおよび無限大に近づく際の第二ビリアル係数の正確な極限挙動を導出すること。
• 正確な解からボイル温度を解析的に決定し、既存の文献値と比較すること。

提案手法

• 第二ビリアル係数の積分に対して変数変換を適用し、∫ y^(-p) e^(-y) e^(2α√y) dy 型の積分に変換する。
• 部分積分を用いて、B₁/₄(α) および B₃/₄(α) を変数とする1階の連立微分方程式に還元する。
• これらの微分方程式を1つの2階微分方程式に統合：d²ψ/dα² - 2α dψ/dα - 3ψ = 0。これは放物型円柱方程式の一種として認識される。
• 変数変換 z = √2 α および ψ(z) = e^(z²/4) φ(z) を施すことで、標準的な放物型円柱方程式に変換され、発展型超幾何関数を用いて解ける。
• 解はクンマーの発展型超幾何関数 M(a,b,z) を用いて表現され、B₂ の閉形式解析的表現が得られる。
• 発展型超幾何関数の漸近展開を用いて、T → 0 および T → ∞ の極限におけるB₂の正確な極限法則を導出する。

実験結果

リサーチクエスチョン

• RQ1 Lennard-Jones流体の第二ビリアル係数の正確な解析的形は、すべての温度で有効であるか？
• RQ2 導出に用いられた積分変換によって、どのような隠れた力学的系が明らかになるか？
• RQ3 温度がゼロおよび無限大に近づく際の第二ビリアル係数の正確な極限挙動は何か？
• RQ4 正確な解析的解から得られるLennard-Jones流体のボイル温度は何か？

主な発見

• 第二ビリアル係数は、放物型円柱関数および発展型超幾何関数を用いて、すべての温度で有効な正確な解析的形で導出された。
• T → ∞ の極限挙動は B₂ → 4√2 Γ(3/4) v₀ (εβ)^(1/4) であり、高温で正の、かつ非常に小さな値を示すことが確認された。
• T → 0 の極限挙動は B₂ → -16√2π v₀ e^(εβ) (εβ)^(3/2) であり、低温で発散する負の値を示す。
• ボイル温度は解析的に T_B* = 3.41793（無次元単位）として決定され、文献値の3.42と非常に高い精度で一致した。
• 正確な解は、積分表法よりも物理的洞察を深めるために、積分の背後にある隠れた調和振動子に類似した力学的構造を明らかにした。
• 本手法は、12および6の指数の組み合わせを持つポテンシャルにのみ適用可能であり、2階微分方程式が閉形式で得られる場合に限る。