[論文レビュー] Thermal response of a Fermi-Pasta-Ulam chain with Andersen thermostats

Federico D’Ambrosio, Marco Baiesi|arXiv (Cornell University)|Jun 21, 2017

Advanced Thermodynamics and Statistical Mechanics参考文献 46被引用数 2

ひとこと要約

本論文は、アンデルセン温度制御を用いて、慣性的で非平衡状態の系に対するフラクチュエーション・レスポンス関係を構築し、微小な摂動に依存せずに熱的感受率（比熱や熱膨張率など）を正確に計算可能にする。この手法は、特に強い非平衡状態において、直接的な摂動シミュレーションを上回り、境界に熱浴を持つフェルミ・パスタ・ユリム鎖の系において滑らかで誤差の小さい感受率推定を可能にする。

ABSTRACT

The linear response to temperature variations is well characterised for equilibrium systems but a similar theory is not available, for example, for inertial heat conducting systems, whose paradigm is the Fermi-Pasta-Ulam (FPU) model driven by two different boundary temperatures. For models of inertial systems out of equilibrium, including relaxing systems, we show that Andersen thermostats are a natural tool for studying the thermal response. We derive a fluctuation-response relation that allows to predict thermal expansion coefficients or the heat capacitance in nonequilibrium regimes. Simulations of the FPU chain of oscillators suggest that estimates of susceptibilities obtained with our relation are better than those obtained via a small perturbation.

研究の動機と目的

• 標準的な平衡状態のフラクチュエーション・レスポンス関係が失敗する、慣性的で非平衡状態の系における熱的感受率の線形応答理論を構築すること。
• 駆動されたフェルミ・パスタ・ユリム（FPU）鎖のような系に熱フラックスが存在する際、比熱や熱膨張率といった熱的応答係数を信頼性高く計算するためのフレームワークの欠如に応えること。
• アンデルセン温度制御が、リッジ・ノイズ系に代わる数学的に取り扱いやすく物理的に意味のある代替手段を提供することを示すこと。
• 数値的に、提案されたフラクチュエーション・レスポンス関係が、シミュレーションにおける微小摂動による推定と比較して、より正確で安定した感受率推定を可能にすることを示すこと。

提案手法

• 確率的経路の重み比に基づき、アンデルセン温度制御で駆動される慣性系における熱的感受率のフラクチュエーション・レスポンス関係を導出する。
• アンデルセン温度制御における速度再スケーリングの確率的性質を活用し、リッジ・ノイズに基づくアプローチで一般的に生じる特異性を回避する。
• 感受率を観測量とエントロピー生成率の関数として表現し、主要な結果は式 (12) から (14) に示され、経路重み比の対称的および反対称的成分を含む。
• FPU鎖の長さに対する境界温度の変化に伴う感受率を、平衡状態および非平衡定常状態の両方を用いて計算する。
• 異なる温度に設定された2つの境界に温度制御を施したFPU鎖の数値的シミュレーションを実施し、導出された式を用いて感受率を計算し、直接的摂動結果と比較する。
• 1/t における線形補外を用いて、t → ∞ のときの静的感受率 χs_L を推定し、収束性を向上させ、統計的ノイズを低減する。

実験結果

リサーチクエスチョン

• RQ1 慣性的で非平衡状態の系、例えば境界に熱浴を持つFPU鎖に対して、熱的感受率のためのフラクチュエーション・レスポンス関係を導出可能か？
• RQ2 導出されたフラクチュエーション・レスポンス式の性能は、熱的感受率推定における直接的微小摂動法と比較してどうか？
• RQ3 系の長さに対する境界温度変化への感受率は、平衡状態と非平衡定常状態で異なるか？ その差は非平衡効果をどのように示唆するか？
• RQ4 非平衡状態において発散する可能性がある経路重み比の対称的および反対称的成分（C+ と C−）を用いて、感受率を正確に再構成できるか？
• RQ5 静的感受率 χs_L は、非平衡定常状態において平均温度と熱フラックスの方向に依存するか？

主な発見

• 導出されたフラクチュエーション・レスポンス関係（式 12–14）は、微小摂動に依存せず、慣性的で非平衡状態の系における熱的感受率を堅牢かつ正確に計算する手法を提供する。
• 数値的結果から、フラクチュエーション・レスポンス式による感受率推定は、直接的摂動シミュレーションによるものよりも滑らかでノイズが少ないことが示された。
• 平衡状態では、応答の対称的および反対称的成分（C+ と C−）が等しくなることが期待通りであり、感受率は標準的なクーブォー公式 χ = 2C− と一致する。
• 非平衡定常状態では、C+ と C− が顕著に乖離し、強い非平衡効果が示唆されるが、依然として完全な式により正しい感受率推定が得られる。
• 静的感受率 χs_L は温度に依存して増加する：平衡状態では T = 0.1 で χs_L = 0.524(3)、T = 0.3 で 0.601(2)、T = 0.5 で 0.625(2) となる。
• 非平衡状態では、熱フラックスの方向に依存する：T1 = 0.1 かつ TN = 0.5 のとき χs_L = 0.625(2)、T1 = 0.5 かつ TN = 0.1 のときには 0.553(3) にまで低下する。これは、平衡平均では捉えきれない真正の非平衡効果を示している。