Nonequilibrium nucleation theory for nonconserved fields: from active matter to population dynamics

本論文は、界面プロファイルの偏差を考慮した反応座標を定義することにより、非保存の秩序変数を持つ系に対する非平衡核生成理論を展開し、その枠組みをアクティブマターや個体群動態のモデルに適用して数値結果と良好な一致を示すことに成功した。

要約

穏やかな湖を眺めているところを想像してみてください。突然、小さな波紋が生まれます。物理学の世界では、これは「メタステーブル（準安定）」状態、つまり一見安定しているように見えるが、実は別の状態へと急変するのを待っている状態のようなものです。通常、新しい相（水の中に氷が形成されることや、新しい種がエコシステムを支配することなど）が現れるためには、「臨界ドロップレット（臨界液滴）」が形成されなければなりません。もしドロップレットが小さすぎれば、それは縮んで消えてしまいます。しかし、もし十分に大きければ、制御不能なほどに成長し、システム全体を飲み込んでいきます。

このプロセスは**核生成（ニュークリエーション）**と呼ばれます。

数十年にわたり、科学者たちは、穏やかでバランスの取れたシステム（平衡状態）において、これがどのように起こるかについてのルールブックを持ってきました。それは、ルールが公平で予測可能なゲームのようなものです。ドロップレットの表面によるエネルギーコストと、その内部によるエネルギー利得を計算すれば、その事象が起こる確率を正確に予測できます。

問題：物事が「アクティブ」になるとルールが変わる
しかし、現実の世界はもっと混沌としており、「アクティブ（能動的）」です。一緒に泳ぐ魚の群れ、自律的に動く細菌、あるいは資源を奪い合う動物の個体群などを考えてみてください。これらのシステムでは、通常のバランス（「詳細釣り合い」と呼ばれます）が崩れています。古いルールブックはもはや通用しません。

以前、科学者は、物（粒子の数など）の総量が一定に保たれるシステムに対してのみ、新しいルールブックを書くことができました。しかし、個体数の増減や、粒子がエネルギーを消費するアクティブ・フルイド（能動流体）のように、総数が一定ではない重要なシステムが多く存在します。これらは「非保存系」です。これまで、このような混沌とした非保存系における核生成を予測する信頼できる方法は存在しませんでした。

解決策： 「丘」を測る新しい方法
この論文は、これらの中乱雑な非保存系に特化した新しい理論（非平衡核生成理論、または NNT）を紹介しています。

ここでの核心となるアイデアを、比喩を使って説明します。

重い岩を丘の上から反対側の谷へと押し上げる場面を想像してください。

• 従来の方法（平衡状態）： 岩が滑らかで対称的な経路を転がり落ちると仮定します。登っていく経路と、転がり落ちる経路は、正確に鏡合わせの形をしています。あなたはただ丘の高さを測れば、答えを知ることができます。
• 新しい現実（アクティブ／非保存系）： 丘は滑りやすく、地面は常に動いています。ランダムなノイズによって岩が押し上げられる「登る」経路は、岩が転がり落ちる「下る」経路とは一致しません。
  • もし「下る」経路を使って「登る」経路を推測しようとすると、間違った答えを導き出してしまいます。丘が実際よりも2倍高いと判断してしまい、本来なら起こりうる事象を「不可能である」と予測してしまうかもしれません。

ブレイクスルー： 「反応座標」の定義
著者らは、ドロップレットのサイズを測定する非常に特殊な方法を考案することで、この問題を解決しました。

ドロップレットを単なる「物質の塊」としてではなく、「境界がぼやけていて、ゆらゆらと動く球体」として捉えます。アクティブなシステムでは、成長するドロップレットの端と、縮小するドロップレットの端は異なって見えます。

• 著者らは、もしドロップレットの「サイズ」を非常に巧妙な数学的手法で定義すれば、その境界の乱雑なゆらぎを無視できることに気づきました。
• これにより、システム全体の複雑で混沌とした動きを、「ドロップレットの大きさはどれくらいか？」という単純な一次元の物語へと投影することができました。
• これによって、成長する際と縮小する際でドロップレットが取る異なる形状に惑わされることなく、エネルギー障壁（丘）の真の高さを計算することが可能になったのです。

検証対象
彼らは単に数学的な計算を行っただけでなく、それを2つの現実的なシナリオでテストしました。

1. 個体群力学： 森の中にいる動物の種を想像してください。時として、新しい、より強い種が森を支配しようと試みます。この論文は、個体数が激しく変動し、「ノイズ（ランダムな出生と死）」がシステムの一部となっている状況下でも、その新しい種が森を制圧できる確率をどのように計算するかを示しました。
2. アクティブ・マター（アクティブ・モデルA）： これは、自律的に動く粒子（例えば細菌や合成ロボット）をモデル化したものです。彼らは、これらのシステムにおいて、ドロップメントが乗り越えるべき「丘」が、従来の理論が予測していたものよりも大幅に異なることを発見しました。

大きな教訓
この論文は、アクティブな非保存系においては、成長の経路が減衰の経路の逆であると単純に仮定することはできず、もしそうすれば間違いを犯す、ということを証明しています。

「核（新しい相の種）」のサイズを測定する方法を慎重に定義することで、彼らはコンピュータ・シミュレーションと完璧に一致する、新しい正確な理論を作り上げました。これは、科学者が、平衡状態に基づいた古い仮定に頼ることなく、細菌コロニーの形成から新しい種の侵入に至るまで、稀に起こる劇的な出来事をより正確に予測できるようになったことを意味しています。

技術概要

技術要約：非保存場における非平衡核形成理論

問題提起
古典的核形成理論（CNT）は、平衡系における安定相から準安定相への形成を、バルクの利得と界面コストの間の自由エネルギーの競合を定量化することによって、成功裏に記述している。平衡状態において、核形成障壁は自由エネルギーのサドルポイントによって決定され、最適軌道（インスタントン）は決定論的な緩和経路の時間反転である。しかし、詳細釣り合いが存在しない（非平衡）系では、これらの単純化は成立しない。すなわち、自由エネルギーに基づく議論では稀なイベントの確率を推定できず、またインスタントン軌道もあらかじめ既知ではない。

保存される秩序パラメータを持つ系（例：Active Model B+）に対しては非平衡核形成理論（NNT）が開発されているが、非保存の秩序パラメータを持つ系への厳密な拡張は依然として困難な課題であった。このような非保存のダイナミクスは、個体群動態、バイオ分子コンデンセート、およびアクティブ物質（例：Active Model A）を含む多くの反応拡散系を支配している。中心となる課題は、ノイズに駆動される液滴の界面プロファイルが、決定論的な緩和過程で見られるプロファイルから逸脱するため、従来の現象論的なCNTの適用では捉えきれない形で核形成障壁を変化させる可能性があることである。

手法
著者らは、以下の確率微分方程式に従う単一のスカラー非保存秩序パラメータ $\phi(x,t)$ を持つ場の方程式に対する、厳密なNNTを構築している：
$$\partial_t \phi = -\mu[\phi] + \sqrt{2T}D(\phi)\xi$$
ここで、$\mu$ は決定論的なドリフト（必ずしも自由エネルギーの汎関数微分ではない）であり、$\xi$ はガウス型ホワイトノイズである。

核心となる手法上の革新は、液滴半径 $R(t)$ の確率論的ダイナミクスを、以下の1次元ランジュバン方程式へと系統的に展開することである：
$$\dot{R} = M(R)[- \partial_R U(R)] + \sqrt{2T M(R)}\zeta$$
ここで、$U(R)$ は準ポテンシャル、$M(R)$ は創発的な移動度、$\zeta$ はノイズである。

導出における主要なステップは以下の通りである：

1. 反応座標の定義： 著者らは、液滴半径 $R(t)$ を恣意的に定義するのではなく、補助関数 $\psi'(r-R)$ を含む積分条件の一意な解として定義した。この特定の選択は、ランジュバン方程式から界面プロファイルの偏差（$\epsilon$）を投影除去するように設計されている。
2. 界面偏差の処理： 平衡CNTや保存系NNTとは異なり、成長する液滴の界面プロファイル（インスタントン軌道）は、緩和プロファイルとは異なることを著者らは示している。彼らは、これらの偏差（$\epsilon_A$ と表記）は小さい（$O(v_0, 1/R, \sqrt{T})$）ものの、適切に処理されない限りエネルギー障壁に影響を与えることを示している。
3. $\psi'$ の選択： $\psi'$ を随伴線形作用素（$L^\dagger$）の核（カーネル）に含まれるように選択することで、ドリフト項における未知の場偏差 $\epsilon$ への依存性を排除し、$M(R)$ および $U(R)$ の解析的な導出を可能にしている。
4. 検証： 解析的な結果は、幾何学的最小作用法（GMAM）とリッツ法を組み合わせた、元の場の方程式のフレイドリン・ウェントセル作用量を数値的に最小化することによって検証されている。

主要な貢献と結果

1. NNTの解析的公式： 本論文は、非保存の非平衡系に対する、臨界半径 $R_c$、移動度 $M(R)$、および準ポテンシャル障壁 $U(R)$ の明示的な式を導出している。
   • 臨界半径は $R_c = \frac{(d-1)\sigma}{\int \psi' \mu_0[\phi_0]}$ で与えられる。
   • 準ポテンシャル障壁 $U(R_c)$ は、ノイズ構造 $D(\phi)$ および $\psi'$ の具体的な選択に依存するため、無ノイズのダイナミクスや界面張力 $\sigma$ のみから推論することはできない。
2. 時間反転緩和（TRR）仮説の破綻： 著者らは、一般的な非平衡系において、インスタントンが時間反転された緩和軌道を辿ると仮定することは誤った結果を導くことを示した。具体的には、TRR仮定は核形成障壁 $U(R_c)$ を $O(1)$ の因子（例：Active Model A では2倍以上）で過大評価し、低ノイズ極限において核形成時間を無限に長く見積もってしまう。
3. 特定モデルへの適用：
   • 個体群動脈： 強力なアリー効果を持つフィッシャー型方程式に適用した。その結果、核形成障壁は乗法的ノイズ $D(\phi)$ の具体的な形式に依存し、平衡系の予測とは異なることを示した。
   • Active Model A (AMA)： 非平衡相形成の最小場理論に適用した。本理論は、活性パラメータ $\lambda$ および外部場 $h$ に対する障壁の依存性を正しく予測する。作用量の数値的最小化は、$R_c$ と $U(R_c)$ の両方について解析的予測を裏付けており、同時にTRR仮定が著しく失敗することも示している。
4. 界面プロファイルの歪み： インスタントン・プロファイルと緩和プロファイルの間の偏差を定量化し、それが $\epsilon_A \sim O(v_0)$ とスケールすることを示し、解析的な仮定と一致することを証明した。

意義
本論文は、非保存の非平衡場に対する、根拠に基づいた初の核形成理論の拡張を提供すると主張している。その主要な意義は、非平衡系における核形成障壁は緩和経路を時間反転させることで近似できるという誤解を正した点にある。著者らは、界面プロファイルにおける偏差が、非保存ダイナミクスにおける障壁計算において主要な効果であることを示している。

反応座標を慎重に定義することで、偏差を投影除去する閉じたランジュバン方程式を復元し、正しい準ポテンシャルを得ることに成功した。この枠組みにより、詳細釣り合いが破れているアクティブ物質から生態学的な個体群動態に至るまでのシステムにおいて、核形成率の解析的な予測が可能となる。本研究は、アクティブ物質におけるCNTの現象論的な適用と厳密な統計力学との間の溝を埋めるものであり、数値的な稀なイベントアルゴリズムのベンチマークを提供するとともに、同様のアプローチがバイオ分子コンデンセートや多成分秩序パラメータ系にも適用可能であることを示唆している。