[論文レビュー] Derivation of the non-equilibrium generalized Langevin equation from a generic time-dependent Hamiltonian

Roland R. Netz|arXiv (Cornell University)|Oct 1, 2023

Advanced Thermodynamics and Statistical Mechanics被引用数 8

ひとこと要約

この論文は時変ハミルトニアンを用いたモリ投影により厳密な非平衡一般化 Langevin 方程式（GLE）を導出し、ガウス観測量は平衡の場合と同じGLE形式に従うことを示し、非ガウス補正およびデータ駆動のパラメータ抽出手法を提供します。

ABSTRACT

It has been become standard practice to describe steady-state non-equilibrium phenomena by Langevin equations with colored noise and time-dependent friction kernels that do not obey the fluctuation-dissipation theorem, but since these Langevin equations are typically not derived from first-principle Hamiltonian dynamics it is not clear whether they correspond to physically realizable scenarios. By exact Mori projection in phase space we derive the non-equilibrium generalized Langevin equation (GLE) from a generic many-body Hamiltonian with a time-dependent force h(t) acting on an arbitrary phase-space dependent observable $A$. The GLE is obtained in explicit form to all orders in $h(t)$. For non-equilibrium observables that correspond to a Gaussian process, the resultant GLE has the same form as the equilibrium Mori GLE, in particular the memory kernel is proportional to the total force autocorrelation function. This means that the extraction and simulation methods developed for equilibrium GLEs can be used also for non-equilibrium Gaussian variables. This is a non-trivial and very useful result, as many observables that characterize non-equilibrium systems display Gaussian statistics. For non-Gaussian non-equilibrium variables correction terms appear in the GLE and in the relation between the complementary force autocorrelations and the memory kernels, which are explicitly given in term of cubic correlation functions of $A$. Interpreting the time-dependent force h(t) as a stochastic process, we derive non-equilibrium corrections to the fluctuation-dissipation theorem and methods to extract all GLE parameters from experimental or simulation data, thus making our non-equilibrium GLE a practical tool to study and model general non-equilibrium systems.

研究の動機と目的

• 時変力 h(t) が観測量 A_S に作用する、一般的な時変ハミルトニアンから厳密な非平衡GLEを導出する。
• モリ投影が得るGLE形が、ガウス観測量に対して平衡の場合と一致することを示し、非ガウス観測量に対する補正を特定する。
• 非平衡GLEを一般化された非平衡フラクドディシペーション関係につなぎ、データからの実用的なパラメータ抽出法を概説する。

提案手法

• 時間依存ハミルトニアン H(ω,t) = H0(ω) - h(t) A_S(ω) を定義する。
• リューヴィル算子と時間順序付き演算子指数を用いてシュレーディンガー描像とハイゼンベルク描像を関連付ける。
• モリ投影を適用して、明示的な記憶カーネルと補完力項を持つ非平衡GLEを得る。
• 硬さ K(t)、記憶カーネル Γ(s,t)、結合 Γ_A(s,t) の厳密な表現を導出し、それらが h(t) に依存することを示す。
• ガウス観測量と非ガウス観測量を議論し、ガウス観測量は改変されたカーネルとともに平衡GLE形式を回復することを示す。
• h(t) を確率過程として扱い、非平衡補正と非平衡フラクドディシペーション関係を導出する。

実験結果

リサーチクエスチョン

• RQ1 時変ハミルトニアン下で、ガウス観測量に対して非平衡GLEは平衡モリGLE構造を保持するか。
• RQ2 非ガウス観測量のゆらぎはGLEと記憶カーネルと力の自己相関との関係をどう修正するか。
• RQ3 h(t) が時変または確率過程のとき、フルクドディシペーション定理の explicit な非平衡補正は何か。
• RQ4 実務的に非平衡時系列データからGLEパラメータをすべて抽出できるか。
• RQ5 h(t) を確率過程と解釈すると、GLEの均質性と解けやすさにどう影響するか。

主な発見

• ガウス観測量の場合、非平衡GLEは平衡モリGLEと同じ形を保ち、記憶カーネルは全力の自己相関に結びつく。
• 非ガウス観測量は、A の三次以上の相関関数に依存する明示的な補正項をGLEに導入する。
• h(t) を確率的とみなすと、時相同性を持つ平均化済みのカーネルと一般化された非平衡フラクドディシペーション関係を得る。
• この枠組みは実験データやシミュレーションデータからGLEパラメータを抽出する明示的な式を提供する。
• このアプローチは h(t) に対して非摂動的で全次に対して厳密であり、非平衡系の現実的なデータベース化モデル化を可能にする。