[論文レビュー] Exactly solvable model illustrating far-from-equilibrium predictions

O. Mazonka, Christopher Jarzynski|ArXiv.org|Dec 7, 1999

Advanced Thermodynamics and Statistical Mechanics被引用数 43

ひとこと要約

本稿では、熱的バスタイムに依存する調和的ばねによって引きずられる粒子の正確に解けるモデルを提示し、任意の引き出しレートにおいて、非平衡統計力学の主要な結果の妥当性を示している。これには、定常状態、一時的、詳細なフラクチュアーション定理、および自由エネルギー差に対する非平衡仕事関係が含まれる。モデルは、平均仕事と分散の明示的表現を通じて、非線形応答を示し、非平衡状態でさえも、仕事分布がガウス分布に従い、Jarzynski等式を満たすことを解析的に確認している。

ABSTRACT

We describe an exactly solvable model which illustrates the Fluctuation Theorem and other predictions for systems evolving far from equilibrium. Our model describes a particle dragged by a spring through a thermal environment. The rate at which the spring is pulled is arbitrary.

研究の動機と目的

• 非平衡統計力学の予測を明確に示す、取り扱いやすく解析的に解けるモデルを提供すること。
• 制御された設定において、フラクチュエーション定理（定常状態、一時的、詳細）の妥当性を検証すること。
• 有限レートで駆動される系において、非平衡仕事関係が自由エネルギー差に対して成り立つかを検証すること。
• 線形応答理論を超えた非線形応答効果が、仕事統計にどのように現れるかを示すこと。

提案手法

• 時間に依存する調和ポテンシャル Uλ(x) = (k/2)(x − λΔx)² + αx における粒子をモデル化し、λ(t) = t/τ により一定レートでばねが引き伸ばされる状況を表す。
• 位置 y と仕事 w を追跡するため、連合確率分布 f(y, w, t | y₀) の時間発展をFokker-Planck方程式で記述する。
• ガウス型のアンザッツと行列形式を用いてFokker-Planck方程式を正確に解き、時間に依存する確率密度関数の明示的表現を得る。
• 初期条件を λ = 0 におけるカノニカル集合からサンプリングすることで、仕事分布 η(W) を積分して計算する。
• 平均仕事 ⟨W⟩ と分散 σ²W を導出し、それらが引き出しレート τ と系のパラメータにどのように依存するかを示す。
• 非平衡仕事関係 ⟨e^−βW⟩ = e^−βΔF が正確に成り立つことを確認し、ゆっくり駆動の極限でも σ²W がゼロでないことを示す。

実験結果

リサーチクエスチョン

• RQ1 線形応答領域を超えて、任意に速いか遅いレートで駆動される系において、一時的フラクチュエーション定理が成立するか？
• RQ2 明示的な時間反転対称性を持つモデルにおいて、詳細フラクチュエーション定理を解析的に検証できるか？
• RQ3 有限レートで駆動される系において、非平衡的散逸を伴っても、非平衡仕事関係 ⟨e^−βW⟩ = e^−βΔF が満たされるか？
• RQ4 非線形応答が存在する場合、仕事分布は平衡状態からどのようにずれるか？ そのずれは解析的に定量化可能か？
• RQ5 駆動される調和系における仕事の平均と分散の正確な形は何か？ それらは自由エネルギー差とどのように関係するか？

主な発見

• 仕事分布 η(W) は正確にガウス分布であり、平均 ⟨W⟩ = αΔx + μγ(Δx)²/τ と分散 σ²W = 2μγ(Δx)²/(βτ) で与えられる。ここで μ(τ) は非線形応答を記述する。
• 平均仕事は自由エネルギー差 ΔF = αΔx を 1/τ の比例定数で上回り、不可逆的散逸が確認される。
• フラクチュエーション-散逸に類似した関係 ⟨W⟩ = ΔF + βσ²W/2 が成り立つが、分散 σ²W は平衡揺らぎではなく、非平衡的ダイナミクスに起因する。
• 非平衡仕事関係 ⟨e^−βW⟩ = e^−βΔF は正確に満たされ、このモデルにおけるJarzynski等式の妥当性が裏付けられる。
• 準静的極限（τ → ∞）では ⟨W⟩ → ΔF かつ σ²W → 0 となり、可逆的・決定論的な自由エネルギー変化が回復される。
• 詳細フラクチュエーション定理は、P₊ と P₋ の正確な連合確率分布の形を通じて検証され、必要な指数関数的比 exp(ΔS/kB) が得られる。