A Generalized Model for Disordered Random Sequential Adsorption with Charge-Dependent Deposition

技術要約：電荷依存型堆積を伴う無秩序なランダム逐次吸着の一般化モデル

問題提起
ランダム逐次吸着（RSA）は、レニーの一次元カー・パーキング問題に代表される、不可逆的な粒子堆積を記述する基本的な非平衡モデルである。標準的な連続系RSAでは、既存の構成に関わらず、幾何学的に許容されるすべての堆積座標が一様にサンプリングされる。しかし、コロイドの結合、タンパク質の吸着、表面反応など、多くの物理的な吸着プロセスは、表面の不均一性、外部電場、または粒子と基板の相互作用によって、到着位置や付着位置にバイアスがかかる影響を受ける。

従来のRSAの拡張では、多くの場合、「クエンチされた（凍結された）」ランドスケープ（例：ランダムなサイトやパターン化された基板）を介して無秩序性を導入するか、試行レートや受理確率を修正する。これに対し、本研究が対象とするのは、堆積測度自体が系の状態に基づいて動的に進化するという、より独特なメカニズムである。具体的には、利用可能なギャップ内における堆積座標が、そのギャップを挟む電荷と、入射粒子の端点電荷に依存する条件付き密度によって決定されるシナリオをモデル化することを目的としている。これは、幾何学的性質と電荷状態が多種型のランダム断片化プロセスとして共に進化する「動的に生成された無秩序」を生み出す。

手法
著者らは、電荷を持つ単位セグメントに関する一般化された一次元RSAモデルを提案している。核心となる革新は、一様な堆積カーネルを電荷依存型のベータ・カーネルに置き換えた点にある。

1. モデル:

   • 空いているギャップの長さを $x$ とし、その両端の電荷を $\alpha, \beta \in \{+, -\}$ とする。
   • 端点電荷 $\sigma, \tau$ を持つ入射単位セグメントが、位置 $u$（左端点）に堆積される。
   • 条件付き確率密度 $P_{\alpha\beta}^{\sigma\tau}(u|x)$ は、正規化されたベータ分布によって定義される：
     $$P_{\alpha\beta}^{\sigma\tau}(u|x) \propto u^{a_{\alpha\sigma}-1}(x-1-u)^{b_{\tau\beta}-1}$$
   • パラメータ $a_{\alpha\sigma} = 1 + \lambda q_\alpha q_\sigma$ および $b_{\tau\beta} = 1 + \lambda q_\tau q_\beta$ は、相互作用パラメータ $0 \le \lambda < 1$ と電荷の積（$q_+ = 1, q_- = -1$）に依存する。
   • 物理的解釈: $\lambda$ は実効的な逆温度（$\lambda = \varepsilon/k_B T$）として機能する。$\lambda=0$ のとき、モデルは一様なカー・パーキングへと回帰する。$\lambda \to 1$ に近づくにつれ、堆積は反対符号の境界付近で局在化し（引力）、同符号の境界付近では抑制される（斥力）。$\lambda \to 1$ の極限は、ゼロ温度極限ではなく、可積分性の閾値を表している。

2. 解析的枠組み:

   • 著者らは、ギャップの長さ $x$ と境界電荷 $\alpha, \beta$ に対する完全な占有統計をエンコードする、**4成分確率母関数（PGF）**方程式 $G_{\alpha\beta}(z, x)$ を定式化している。
   • このPGFは、堆積イベントによって生成される2つの娘ギャップの条件付き独立性に由来する、積分再帰式を満たす。
   • このPGFから、平均占有数（$f_{\alpha\beta}$）および分散（$s_{\alpha\beta}$）に関する閉じた再帰方程式が導出される。

3. 数値解および厳密解:

   • 厳密な有限シェル解: 区間 $x \in [n, n+1)$ （シェル）に対する解析解が導出されている。最初の非自明なシェル（$2 \le x < 3$）は、分布、平均、および分散の厳密な表現を与える。
   • ガウス・ヤコビ求積法: ベータ・カーネルの端点の特異性を効率的に扱うため（特に $\lambda \to 1$ の場合）、著者らはガウス・ヤコビ求積法を採用している。
   • モンテカルロ・シミュレーション: 条件付きカーネルから最終的なジャム状態（充填状態）を直接サンプリングするための、棄却なしモンテカルロ・アルゴリズムが開発され、求積法の結果に対する独立した検証として用いられている。

主な結果

1. 厳密な再帰式とベンチマーク: 本論文は、完全な占有分布に関する厳密な再帰式を提供している。最初の数シェルの解析解は厳密なベンチマークとして機能し、ガウス・ヤコビ法およびモンテカルロ法が、平均および分散を高い精度（$10^{-4}$ 以内）で厳密解と一致させることを確認している。
2. 漸近的な自己平均性: 中心的な理論的結果として、ジャム占有密度が自己平均的であることの証明が示されている。
   • 初期境界電荷は有限サイズ挙動に影響を与えるが、漸近的な平均密度（$\rho$）および分散密度（$\chi$）は、境界状態（$\alpha, \beta$）に依存しない。
   • これは、電荷状態の遷移行列にパロン・フロベニウスの定理を適用することで証明される。反復的な堆積プロセスが初期の電荷情報を粒子の分布統計へと置き換えるため、系はバルクにおいて初期の境界条件を「忘却」する。
3. 電荷選択性（$\lambda$）と組成の影響:
   • 被覆率: 電荷選択的な堆積は、一様なRSAと比較してジャム密度を一般に増加させる。これは、境界に好まれるイベント（反対電荷）が、セグメントをギャップの端付近に配置することで、より大きな利用可能な娘区間を残すためである。
   • 粒子混合物: この効果は、反対の端点電荷（$p_{+-}, p_{-+}$）が支配的な混合物において最も顕著である。$\lambda = 0.99$ において、この混合物は密度 $\rho \approx 0.9451$ を達成し、これは一様なRSAの極限（$\approx 0.7476$）や、同符号電荷が支配的な場合（$\approx 0.8670$）よりも大幅に高い。
   • ゆらぎ: 分散密度（$\chi$）は明確な挙動を示す。反対電荷が支配的な混合物では、高密度に伴ってゆらぎが強く抑制される（$\chi \approx 0.0243$）ことが示され、これはより再現性の高いジャム状態であることを示唆している。逆に、同符号電荷が支配的な混合物では、抑制されるチャネルと促進されるチャネルの競合により、ゆらぎが増大する。

意義および主張
本論文は、以下の3つの異なる領域間の扱いやすい接続を提供すると主張している：

1. 一様なRSA: $\lambda = 0$ で回帰する。
2. 動的に生成された無秩序: 堆積測度がシステムの状態と共に進化するモデルであり、クエンチされた無秩序とは異なる。
3. 特異な電荷セグメント動力学: $\lambda \to 1$ において、著者らの先行研究で扱われた決定論的かつ電荷選択的なルールへと連続的に接続される。

著者らは、この枠組みが、状態依存的な堆積を伴う連続系RSAモデルにおいて、完全な占有統計（単なる平均だけでなく）の解析的研究を可能にすることを強調している。複雑で履歴依存的な電荷の進化にもかかわらず、ジャム密度が自己平均的であるという実証は、重要な理論的貢献である。これらの結果は、電荷の組成と相互作用の強さが、無秩序なジャム層の効率（充填密度）と再現性（ゆらぎ）の両方を制御するために使用できることを示唆している。

本論文は、特定の実験セットアップや将来の応用を提案するものではなく、その貢献は純粋に統計力学および非平衡確率過程の領域内にある。