[論文レビュー] Two theorems that relate discrete stochastic processes to microscopic mechanics

Hal Tasaki|ArXiv.org|Jun 7, 2007

Stochastic processes and statistical mechanics参考文献 1被引用数 8

ひとこと要約

本稿は、確率的で離散時間のマルコフ過程において、深いつながりのないポテンシャル井戸を粒子が飛び移動する場合の物理的に正しい遷移確率を、微視的ハミルトニアン力学と有効な確率的挙動とを厳密に結びつけることで導出する。一般条件下で詳細つり合いが成立することを証明し、遷移率が exp(−β(b_{x,y} − F_x)) に比例する明示的形を導出する。ここで b_{x,y} は障壁高さ、F_x は制限自由エネルギーであり、非平衡モデルに物理的根拠を持つルールを提供する。

ABSTRACT

Starting from a classical mechanics of a ``colloid particle'' and $N$ ``water molecules'', we study effective stochastic dynamics of the particle which jumps between deep potential wells. We prove that the effective transition probability satisfies (local) detailed balance condition. This enables us to rigorously determine precise form of the transition probability when barrier potentials have certain regularity and symmetry.

研究の動機と目的

• 微視的ハミルトニアン力学から遷移確率を導出することにより、非平衡系における物理的に一貫した確率的ダイナミクスの枠組みを確立すること。
• 非平衡確率的モデルにおける強いルール依存性問題を解消するために、物理的根拠を持つ遷移確率ルールを同定すること。
• 駆動力が働く系であっても、一般条件下で有効な遷移確率が詳細つり合いを満たすことを証明すること。
• 障壁高さと制限自由エネルギーの観点から、遷移確率の明示的で厳密な形を導出し、平衡系および駆動系の両方に適用可能なものとする。
• 物理的に意味のあるダイナミクスを持つ離散確率的モデルを用いて、非平衡定常状態を研究する基盤を提供すること。

提案手法

• コロイド粒子と N 個の水分子を含む周期的ボックス内の古典的力学系を定式化し、粒子および分子の運動エネルギー、外部ポテンシャル、粒子-分子相互作用を含むハミルトニアンを導入する。
• 粒子が局所的に捕らえられる空間的領域 W_x を用いてポテンシャル井戸を定義し、それらの間には高い障壁が存在する。
• 各井戸内で準平衡状態をモデル化するため、制限された正準集合を用い、分配関数 Z_x と特性関数 χ_x を定義する。
• 有効遷移確率 P^{(T)}_{x→y} を、ハミルトニアン力学に従って W_x に初期位置を持つ粒子が時間 T 後に W_y に見つかる確率として定義する。
• T_rel ≪ T ≪ T_wait を満たす時間 T を設定することで、高々1回の遷移しか起こらないようにし、マルコフ近似を正当化する。
• 時間反転対称性を用いて詳細つり合いを証明し、W_y で一定ポテンシャルを持つ修正モデルを導入することで、P^{(T)}_{x→y} が障壁高さと自由エネルギーの観点から評価可能であることを示す。

実験結果

リサーチクエスチョン

• RQ1 確率的モデルにおける深いポテンシャル井戸間を粒子が飛び移動する場合の、物理的に正しい遷移確率の形は何か？
• RQ2 微視的ハミルトニアン力学から導かれる有効な確率的ダイナミクスは、詳細つり合いを満たすか？
• RQ3 外部場の存在が遷移確率に与える影響は何か？また、詳細つり合い条件は依然として保持されるか？
• RQ4 現実的な仮定のもとで、障壁高さと制限自由エネルギーの観点から、遷移確率を厳密に表現できるか？
• RQ5 ポテンシャルの対称性および正則性が、遷移確率の形に果たす役割は何か？

主な発見

• 有効遷移確率 P^{(T)}_{x→y} は局所的詳細つり合い条件を満たしており、平衡統計力学と整合的であることが保証される。
• 遷移確率が厳密に exp(−β(b_{x,y} − F_x)) に比例することが示され、ここで b_{x,y} は障壁高さ、F_x は井戸の制限自由エネルギーである。
• 導出された形 (21) は、平衡系および駆動系の両方で有効であり、外部場は自由エネルギー項に ±E/2 の寄与を、場の方向に与える。
• 係数 ρ(T) または ˜ρ(T) は障壁の詳細に依存するが、井戸の識別番号とは無関係であり、普遍的なスケーリング行動を示す。
• 結果として、非平衡モデルにおける確率的ダイナミクスに物理的根拠を持つ、恣意的でないルールが得られ、強いルール依存性問題が解決される。
• ガウス型ポテンシャルおよび弱い結合の仮定のもとで明示的な計算が可能となり、P^{(T)}_{x→x±e₁} ∝ ρ(T) exp[−β(b_{x,x±e₁} − F_x^{(0)} ∓ E/2)] が得られる。ここで F_x^{(0)} は明示的に与えられる。