Bubble formation in active binary mixture model

技術要約：活性二成分混合モデルにおけるバブル形成

問題提起
能動的相分離、あるいは運動誘起相分離（MIPS）は、自己駆動粒子が自発的に高密度相と低密度相に分離する集団現象である。能動的相分離の顕著な特徴は、平衡状態の相分離とは異なり、高密度相の中に低密度相のバブル（泡）が持続的に形成されることである。Active Model B+ (AMB+) のような理論モデルは、この挙動を引き起こす特定の非平衡項（$\lambda$ 項および $\zeta$ 項）を特定しているが、微視的な粒子モデルにおいてバブル形成を系統的に制御する基礎的なパラメータは依然として不明である。さらに、バブルの存在は、臨界挙動や普遍性クラスの調査を困難にする。なぜなら、バブルが臨界指数を平衡状態のIsing普遍性クラスから変化させる可能性があるためである。

手法
著者らは、これらの現象を調査するために、長方形格子上の活性二成分混合モデルを導入する。このシステムは、溶質（solute）と溶媒（solvent）の2種類の種で構成され、両者は極性（自己推進の方向）を持つ。モデルは以下の確率的ダイナミクスに従う：

1. 極性の回転： 粒子は、回転率 $\gamma$ で極性を $\pm 90^\circ$ 回転させる。
2. 交換ダイナミクス： 隣接する溶質・溶媒ペアが位置を交換する。
   • 受動的交換： 極性に関わらず、単位レート（1）で行われる。
   • 能動的交換： 各粒子の極性が運動方向に一致している場合、それぞれのレート $\varepsilon_{\text{solute}}$ および $\varepsilon_{\text{solvent}}$ で発生する。
   • 制約： 同一種間の交換（溶質同士または溶媒同士）は禁止されている。

本研究では、主に3つの解析手法を用いている：

• 数値シミュレーション： モンテカルロ（MC）シミュレーションを実行し、活動性の非対称性（$\varepsilon_{\text{solute}}$ 対 $\varepsilon_{\text{solvent}}$）および密度（$\rho$）の変化に伴う定常状態の相図、バブルサイズ分布、および時間発展を観察する。
• 平均場理論（MFT）： 本質的な相挙動を捉え、線形安定性解析を行うために、単一サイトおよび4サイトの平均場近似を導出する。4サイト近似は、単一サイトモデルで無視される局所的な相関を取り入れている。
• 有限サイズスケーリング解析： バブル形成が抑制される対称極限（$\varepsilon_{\text{solute}} = \varepsilon_{\text{solvent}}$）において、臨界指数（$\nu, \beta, \gamma$）を分析し、その普遍性クラスを決定する。

主な結果

1. 活動性の非対称性がバブル形成を制御する：

   • 中程度の非対称性： 溶質の活動性に対して溶媒の活動性が中程度である場合（例：$\varepsilon_{\text{solvent}} \approx 0.2-0.6 \varepsilon_{\text{solute}}$）、バブル形成が著しく促進される。深く相分離した状態において、バブルサイズ分布はべき乗則 $N(A) \sim A^{-\alpha}$ （$\alpha \approx 1.75$）に従う。その結果、バブル分率 $f_{\text{bubble}}$ はシステムサイズに対して $L_x^{0.5}$ でスケーリングし、バブル相分離またはマイクロ相分離への傾向を示す。
   • 対称な活動性： $\varepsilon_{\text{solute}} = \varepsilon_{\text{solvent}}$ のとき、バブル形成は抑制される。システムは、どちらの相の中にもバブルを成長させる駆動力となる対称性が失われるため、持続的なバブルを伴わないマクロな相分離を示す。
   • 高い溶媒活動性： 溶媒の活動性が溶質の活動性に匹着するかそれを上回ると、二つの相の間のマクロな対称性が回復するため、バブルは抑制される。

2. 増強メカニズム：

   • 平均場理論と線形安定性解析により、溶質の活動性と溶媒の活動性が競合していることが明らかになった。溶質の活動性は、均一状態をマクロな相分離へと不安定化させる傾向がある一方、溶媒の活動性は、密度に応じて均一状態を安定化させるか、あるいはマイクロ相分離を誘起する。
   • バブルが増強される領域では、界面付近において溶媒の極性密度が溶質に富む相に向かって整列する。この溶媒の活動性に駆動された整列が、溶質に富むバルク内での溶媒バブルの拡大を促進する。

3. 臨界挙動と普遍性：

   • バブルが抑制される対称極限（$\varepsilon_{\text{solute}} = \varepsilon_{\text{solvent}}$）において、著者らは能動的相分離の臨界指数を推定している。
   • 有限サイズスケーリング解析により、$\nu \approx 0.89$、$\beta/\nu \approx 0.119$、$\gamma/\nu \approx 1.77$ という指数が得られた。
   • これらの比率は、二次元Ising普遍性クラスの値（$\nu=1, \beta/\nu=0.125, \gamma/\nu=1.75$）にかなり近いが、推定された $\nu$ はわずかに低い。著者らは、漸近的な指数がIsingクラスと一致するか、あるいは逸脱するかを確定的に確認するには、より大規模なシミュレーションが必要であると述べている。

意義と主張
本論文は、活動性の非対称性を、能動物質の相転移における主要な制御パラメータとして確立している。溶媒に自己推進能力を持たせることで、標準的な能動粒子モデルでは制御が困難であった「バブル形成の強化または抑制」という調整可能なメカニズムを提供している。

著者らは、以下の点で自らの知見が非平衡系の普遍性に新たな洞察を与えると考えている：

• 本モデルは、ドライな能動物質（溶媒が実質的に空の空間であるもの）と、多成分能動系との間のギャップを埋めるものである。
• 対称性によってバブルを抑制できることにより、能動的相分離がIsing普遍性クラスに属するかどうかをテストするための、よりクリーンな基準となる臨界現象の調査が可能になる。
• 平均場解析は、微視的な交換規則とマクロな相挙動を結びつけ、バブル増強および抑制の定性的なメカニズムをうまく再現している。

結論として、現在の結果はバブル抑制下でのIsing普遍性クラスへの強い関連性を示唆しているものの、決定的な分類にはさらなる大規模シミュレーションが必要である。また、活動性の非対称性によるバブルの制御可能性は、活性ブラウン粒子などの他の二成分能動系においても観察可能であることが示唆されている。