[論文レビュー] Solution of the Metropolis dynamics on a complete graph with application to the Markov chain Mpemba effect

Israel Klich, Marija Vucelja|arXiv (Cornell University)|Dec 31, 2018

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

ひとこと要約

本稿では、完全グラフ上のメトロポリスダイナミクスを解析的に解き、すべての固有値と固有ベクトルを導出することで、系が単調な緩和を示すため、弱いまたは強いメプムバ効果を示さないことが明らかになった。しかし、遷移が遮断される（例：スパースグラフ上）と、より高温の系が速く冷却されるというメプムバ効果が出現する可能性がある。これは、1つの辺が欠落した3状態系の例で示された。

ABSTRACT

We find analytically the complete set of eigenvalues and eigenvectors associated with Metropolis dynamics on a complete graph. As an application, we use this information to study a counter-intuitive relaxation phenomenon, called the Mpemba effect. This effect describes situations when upon performing a thermal quench, a system prepared in equilibrium at high temperatures relaxes faster to the bath temperature than a system prepared at a temperature closer to that of the bath. We show that Metropolis dynamics on a complete graph does not support weak nor strong Mpemba effect, however, when the graph is not complete, the effect is possible.

研究の動機と目的

• エネルギー準位に非解析的依存があるにもかかわらず、完全グラフ上のメトロポリスダイナミクスの全固有系を解析的に決定すること。
• このような系において、より高い初期温度からの速い緩和（メプムバ効果）が生じるかどうかを調査すること。
• メトロポリスダイナミクスに従うマコフ跳躍過程において、メプムバ効果が生じる条件を調査すること。
• 完全グラフと不完全グラフの挙動を対比させ、メプムバ効果を生じさせる構造的要件を同定すること。

提案手法

• ステディステート確率 $ \pi_i $ を降順に再順序付けするための置換行列 $ S_r $ を用い、スペクトル解析を単純化する。
• 再順序付けられた基底において、メトロポリスレート行列 $ R $ の固有値 $ \lambda_k $ と固有ベクトル $ |v_k\rangle $ を導出する。ここで $ \lambda_1 = 0 $ であり、$ k > 1 $ の場合 $ \lambda_k = -\left(\pi_{r_{k-1}}^{-1}Z_{k-1} + k - 2\right) $ である。
• 緩和ダイナミクスの解析を可能にするために、最も緩慢に減衰する固有ベクトル $ |v_2\rangle $ を累積確率 $ Z_k = \sum_{l=k}^n \pi_{r_l} $ の形で表現する。
• LuとRaz（2017）の公式を用い、緩和速度を支配する係数 $ a_2(\beta) $ を、固有ベクトル $ |v_2\rangle $ とバスタームのエネルギー分布を用いて計算する。
• 非単調な振る舞いを示す必要がある弱いまたは強いメプムバ効果を排除するため、$ a_2(\beta) $ の単調性を解析する。
• 1つの遷移が遮断された3状態モデルを構築し、スパースグラフにおけるメプムバ効果の出現を示し、非単調な $ a_2(\beta) $ を示す。

実験結果

リサーチクエスチョン

• RQ1 完全グラフ上のメトロポリスダイナミクスは、非単調な緩和挙動として定義されるメプムバ効果を支持するか？
• RQ2 遷移グラフの構造的特徴（例：完全性）が、メプムバ効果を生じさせるために必要な条件は何か？
• RQ3 完全グラフ上でのメトロポリスレート行列のスペクトル分解が、緩和ダイナミクスに与える影響は何か？
• RQ4 直接の遷移がグラフから削除された場合、メプムバ効果がマコフ跳躍過程で出現する可能性はあるか？
• RQ5 緩和係数 $ a_2(\beta) $ の単調性が、メプムバ効果の有無を決定づける役割を果たすか？

主な発見

• 完全グラフ上でのメトロポリスダイナミクスの全固有系は解析的に解ける。固有値は $ \lambda_k = -\left(\pi_{r_{k-1}}^{-1}Z_{k-1} + k - 2\right) $ であり、固有ベクトル $ |v_k\rangle $ は累積確率 $ Z_k $ の形で表現される。
• 緩和係数 $ a_2(\beta) $ は逆温度 $ \beta $ に関して単調であるため、完全グラフの場合には弱いまたは強いメプムバ効果が生じないことが示唆される。
• 完全グラフでは、メプムバ効果が現れないのは、$ e^{\beta E_1}Z_1(\beta) $ の厳密な単調性に起因し、$ a_2(\beta) $ が $ \beta = \beta_b $ のみで消失し、中間のゼロ点を持たないからである。
• 一方、1つの遷移が遮断された場合（例：3状態系）、緩和係数 $ a_2(\beta) $ は非単調になり、$ a_2(\beta_b) = 0 $ かつ $ \lim_{\beta \to \infty} a_2(\beta) = 0 $ となるため、間接的なメプムバ効果が可能になる。
• スパースグラフでは、対称性の破れと緩和ダイナミクスの非単調性により、メプムバ効果が生じる。完全に接続された場合とは対照的である。
• 結果から、メプムバ効果はメトロポリスダイナミクスの一般的特徴ではなく、下位の遷移グラフのトポロジーに強く依存しており、スパースグラフが非自明な緩和現象を可能にすることが示唆される。