Ergodicity Breaking in Active Run-and-Tumble Particles in a Double-Well Potential

技術要約：二重ウェルポテンシャルにおける能動的ラン・アンド・タンブル粒子のエルゴード性破れ

問題提起
本研究は、一次元対称二重ウェルポテンシャル $U(x) = \frac{1}{4}(x^2 - a^2)^2$ 内を移動するラン・アンド・タンブル粒子（RTP）のダイナミクスを調査している。障壁横断や受動的なブラウン粒子における定常状態は、よく理解されている現象であるが（クラマースの法則および一意なボルツマン分布に従う）、熱雑音が存在しない状況下での能動的粒子の振る舞いは、依然としてほとんど解明されていない。具体的には、著者らは、RTPの固有の能動的性質（速度 $v_0$ と方向転換率 $\gamma$ による持続的な運動）が、有限の障壁を持つランドスケープにおける位置のゆらぎ、エルゴード性、および障壁横断時間にどのように影響するかという問題に取り組んでいる。

手法
著者らは、粒子の位置 $x(t)$ を以下のランジュバン方程式を用いてモデル化している：
$$\dot{x}(t) = f(x) + v_0\sigma(t)$$
ここで、$f(x) = -U'(x)$ は決定論的な力であり、$\sigma(t) = \pm 1$ は粒子の内部方位を表す二値ノイズ（転換率 $\gamma$ で切り替わる）である。

解析は、ポテンシャル力 $f(x)$ と能動速度 $v_0$ の関係を検討することから進められる。著者らは、臨界障壁高さ（パラメータ $a$ によって制御される）を特定している。そして、以下の2つのレジームを区別している：

1. エルゴード・レジーム ($a < a_c$): 力の大きさ $|f(x)|$ が関連領域内で $v_0$ を決して超えないため、粒子は転換点間の全領域にアクセスできる。
2. 非エルゴード・レジーム ($a > a_c$): 力の大きさが特定の領域で $v_0$ を超えるため、「実効的な壁」が形成され、粒子は方位に関わらずその力を克服することができない。

定常状態を特徴付けるために、著者らは以下を用いている：

• 解析的導出: 閉じ込められたポテンシャル中のRTPに関する既知の結果を用いて、定常確率密度関数 $P_{st}(x)$ を導出する。
• ヒット確率（到達確率）: 非エルゴード・レジームにおいて、中間領域 $|x_0| < \ell_1$ から出発した粒子が、最終的に左または右のウェルにトラップされる「分裂確率」（またはヒット確率）$\pi(x_0)$ を計算するために、後退フォッカー・プランク方程式を解く。
• 数値シミュレーション: 解析的な予測を、ランジュバン方程式の数値積分によって検証する。

主要な貢献および結果

1. 強いエルゴード性破れ:
   本論文は、障壁高さが臨界閾値（$a > a_c$）を超える場合、系が強いエルゴード性破れを示すことを証明している。初期条件に依存せず一意な長時間分布を持つ受動的システムとは異なり、RTPの定常分布は初期位置 $x_0$ に決定的に依存する：

   • $x_0$ が一方のウェルの吸引圏内（具体的には $|x_0| > \ell_1$）にある場合、粒子はそのウェルに永久にトラップされる。
   • $x_0$ が中間領域（$|x_0| < \ell_1$）にある場合、粒子は有限の確率 $\pi(x_0)$ でいずれかのウェルに到達できる。長時間分布は、$\pi(x_0)$ によって決定される、2つのウェルにおける分布の加重和となる。
   • 著者らは、$\pi(x_0)$（式12）および結果としての非エルゴード定常分布（式13）の厳密な解析的表現を提供している。

2. 臨界障壁高さと転換点:
   臨界値は $a_c = \left(\frac{3\sqrt{3}v_0}{2}\right)^{1/3}$ として導出される。

   • $a < a_c$ の場合、2つの転換点 ($\pm \ell$) が存在し、系は $[-\ell, \ell]$ 上に支持された一意の定常分布を持つエルゴード的な系となる。
   • $a > a_c$ の場合、$f(x) = \pm v_0$ という方程式は6つの転換点 ($\pm \ell_1, \pm \ell_2, \pm \ell_3$) を与える。これらの点は位相空間を断片化された領域に分割し、中間ゾーンの外側から出発した場合、粒子が中央の障壁を越えることを妨げる。

3. 障壁横断時間とクラマースの法則の破れ:
   エルゴード・レジーム ($a < a_c$) において、著者らは障壁頂上での定常密度を介して定義される特性障壁横断時間 $\tau$ を推定している：$\tau = 1/[v_0 P_{st}(0)]$。

   • 本研究は、$\tau$ が標準的なクラマース・アレニウス則 ($\tau \sim e^{\beta \Delta V}$) を満たさないことを示している。
   • 障壁高さが臨界閾値に下から接近するにつれ（$a \to a_c^-$）、横断時間は発散する。
   • この発散は、異常な指数を持つヴォーゲル・フルマー・タンメルマン型の形式に従う：$\tau \sim (a_c - a)^{-1/2}$。

意義および主張
著者らは、本研究が、熱雑音のない二重ウェルポテンシャル内における能動的粒子のエルゴード性破れに関する初の解析的特性付けを提供するものであると主張している。主な意義は、能動的なダイナミクスが障壁横断の統計力学を根本的に変え得ることを示した点にある：

• 非エルゴード性: 系が全位相空間を探索できない非エルゴード相の存在は、能動的なダイナミクスの直接的な帰結であり、これは受動的なブラウン系には見られない特徴である。
• 初期条件依存性: 定常状態は一意ではなく、初期条件に依存する。これは、能動物質に適用される標準的な平衡統計力学の仮定に挑戦するものである。
• 異常スケーリング: 臨界点付近での障壁横断時間の指数 $1/2$ を伴う発散は、熱的活性化と比較して、能動的な障壁横断の明確に異なる普遍性クラスを浮き彫りにしている。

結論として、有限の温度は熱ゆらぎによる障壁横断を可能にするが、「純粋な能動的性質」のダイナミクスは、平衡現象とは質的に異なるトラッピングとエルゴード性破れのメカニズムを導入することを、本論文は示している。