Effective interactions in active Brownian particles

問題の背景
能動ブラウン粒子（ABP）は、自己推進力を駆動するためにエネルギーを散逸させる非平衡系であり、運動誘起相分離（MIPS）、クラスター形成、フレーミングといった集団的挙動を示す。この分野における大きな課題は、能動的物質に対する直接的な熱力学的枠組みが欠如していることである。具体的には、圧力や化学ポテンルといった概念がこれらの系に対して厳密に定義できるのか、また、その非平衡構造が有効平衡ポテンシャルによって記述できるのかどうかは依然として不明である。先行研究では、希薄系に対する有効ポテンシャル（平均力ポテンシャルによるもの）や、密度に依存しない多体相互作用が導出されているが、様々な密度や相互作用タイプにわたって能動懸濁液の構造を正確に再現できる、単一の密度依存的な有効対ポテンシャルを得る手法は欠けていた。

手法
著者らは、二次元ABP懸濁液の有効対ポテンシャル（$u_{eff}(r)$）を決定するために、テスト粒子挿入（TPI）法に基づく逆手法を用いている。このアプローチは、動径分布関数 $g(r)$ を計算するための2つの異なるスキームを一致させることを目的としている。

1. 距離ヒストグラム（DH）法： 能動粒子の定常状態のスナップショット（LAMMPSを用いて実施）から $g_{DH}(r)$ を直接計算する。
2. テスト粒子挿入（TPI）法： 仮説上の有効対ポテンシャルから導かれた有効ポテンシャルエネルギー $\Psi_{eff}$ を用いて $g_{TPI}(r)$ を計算する。

手法の中核は、初期のポテンシャルの推定値（平均力ポテンシャル $w(r)$）を、$g_{TPI}(r)$ が $g_{DH}(r)$ に収束するまで更新していく反復的な予測子・修正子スキーム（シュモマーズのアルゴリズム）である。得られる $u_{eff}(r)$ は、モンテカルロ法を用いてシミュレーションを行った際に、元の能動系の構造相関（$g(r)$）を再現するような「受動的」な系のポテンシャルとして定義される。

$u_{eff}(r)$ が得られた後、著者らは系を平衡状態にある受動系として扱い、以下の値を算出する：

• 有効化学ポテンシャル ($\mu_{eff}$): 元のABPの座標を用い、過剰化学ポテンシャルのためのTPI公式を用いて導出される。
• 有効圧力 ($P_{eff}$): 体積の変化に伴う有効自由エネルギーの変化を調べるために、テストボリューム法（ビリアル方程式のノイズの多い微分を避ける手法）を用いて算出される。

本研究では、3つのパラメータ（基礎となる受動的相互作用ポテンシャル：レナード・ジョーンズ、ウィークス・チャンドラー・アンダーソン、および二段階長さスケール・ショルダーポテンシャル）、ペクレ数（$Pe$、活性の強さを表す）、および粒子数密度（$\rho$）を系統的に変化させている。

主な結果

• 構造の再現： 逆手法は、受動的モンテカルロシミュレーションにおいて元の能動系の $g(r)$ を正確に再現する有効対ポテンシャル $u_{eff}(r)$ を生成することに成功した。これは様々なポテンシャル型や密度において成立しており、一定の温度と密度において、与えられた定常状態の構造に対して一意の有効対ポテンシャルが存在することを裏付けている。
• 創発的な引力： 純粋に反発的な受動ポテンシャル（WCAおよびショルダー）に対して、導出された $u_{eff}(r)$ は引力の井戸を示す。これは、活性が実効的な引力を誘起していることを示しており、この現象がクラスター形成やMIPSを駆動している。井戸の深さと形状は、ペクレ数と密度に依存する。
• 密度依存性： 真の平衡ポテンシャルとは異なり、導出された $u_{eff}(r)$ は明示的に密度依存的である。密度が増加するにつれて有効相互作用は変化し、これは能動系における粒子の速度低下を反映している。
• 平均力ポテンシャルとの比較： 低密度において、有効ポテンシャル $u_{eff}(r)$ は平均力ポテンシャル $w(r)$ とよく一致する。しかし、密度が高くなるにつれて両者は乖離し、$u_{eff}(r)$ は $w(r)$ が高密度領域で捉えきれない高次の構造相関を考慮している。
• 熱力学的性質： 著者らは、能動系を等価な受動系として扱うことで、有効化学ポテンシャルと有効圧力を測定できることを報告している。ショルダーポテンシャルの場合、$\mu_{eff}$ と $P_{eff}$ は密度および活性に対して非単調な挙動を示し、構造転移と相関する極大値を示す。

意義と主張
本論文は、非平衡な能動系の構造を、平衡論的な有効対ポテンシャルを用いて記述するための実践的なアプローチを提供すると主張している。主な意義は、ABPの本質的な非平衡性にもかかわらず、その定常構造が、密度依存的な有効ポテンシャルを持つ等価な受動系へと写像できることを示した点にある。

著者らは、自らの貢献を「能動的物質を記述する統計力学的枠組みの構築に関する問いに光を当てる」ためのツールとして控えめに位置づけている。彼らは、これらの有効ポテンシャルが能動系の真の熱力学的状態変数であるとは主張していない（能動系における機械的圧力は議論の分かれる状態変数であることに触れている）。代わりに、これらの有効ポテンシャルを用いることで「有効な」化学ポテンシャルや圧力を測定することが可能であり、能動的物質の動力学と平衡統計力学との間の架け橋となることを提唱している。本研究は、活性が複雑な多体効果をもたらす一方で、それらが構造解析のための、二体かつ密度依存的なポテンシャルによって効果的に捉えられる可能性があることを示唆している。