Continuum Limit of Nonlocal Electrostatics in Random Media

技術要約：ランダム媒体における非局所静電法の連続体極限

問題提起
古典的な連続体静電理論は、微視的な電荷分布を抑制しつつ、分極のエネルギー的帰結を記述する。しかし、有限温度においては、電荷を持つ原子が平均的な位置の周囲で揺らぐため、ランダムな電荷密度場が生じる。クーロン相互作用の長距離性により、これらのランダムな微視的場からマクロなエネルギーへの遷移は非自明なものとなる。これまでの解析では、クーロン物質の熱力学的極限、確率的格子における均質化、および無秩序結晶の平均場モデルが扱われてきた。本論文は特定の空白領域、すなわち、微視的なゆらぎが局所的な静電相互作用に与える影響と、それがマクロな分極場に与える影響を区別する、所定の定常ランダム電荷密度に対する三次元連続体エネルギー分解の導出に取り組むものである。

手法
著者らは、領域を以下の3つの長さスケールに分離する、二スケール連続体極限アプローチを採用している：微視的スケール $\ell$（原子格子）、材料点スケール $\varepsilon$（多数の原子を含む）、およびマクロスケール $L$（物体）。ここで $\ell \ll \varepsilon \ll L$ である。

1. 数学的枠組み: 微視的な電荷密度 $\rho(x, y, \omega)$ は、格子並進に対して定常かつエルゴード的なランダム場としてモデル化される。この場は、参照セル電荷密度 $\bar{\rho}$ と、測度保存的なエルゴード的格子作用 $\{T_k\}$ によって生成される。
2. エネルギー分解: Marshall and Dayal [24] に従い、全離散静電エネルギーを以下のように分解する：
   • 局所エネルギー ($E_{loc}$): 同一の材料点内にある電荷間の相互作用。
   • 非局所エネルギー ($E_{nloc}$): 異なる材料点間にある電荷間の相互作用。
3. 平均化と収束: 著者らは、格子エルゴード定理を用いて、大きな数の微視的セルに対する空間平均が、決定論的なアンサンブル平均にほとんど確実に収束することを示す。彼らは、有限体積ポテンシャルが無限体積極限へと収束することを分析し、全空間ポテンシャルと、有限の材料点領域に由来する境界補正ポテンシャルの違いを区別する。
4. 仮定: 導出は、以下の3つの主要な仮定に基づいている：
   • (A1) 電荷密度の写像の連続性と、すべての参照セル実現内における厳密な電荷中性。
   • (A2) 微視的クーロンポテンシャルの存在および $L^2$ 可積分性。
   • (A3) 有限体積場補正の、一様に分極した立方体のポテンシャルの勾配への強収束。

主要な貢献と結果
本論文は、上記の仮定の下で、離散的な局所および非局所エネルギーが決定論的な連続体極限にほとんど確実に収束することを確立している。

1. 局所極限: 局所エネルギーは、以下の2つの項の和に収束する：
   • アンサンブル平均された微視的クーロンエネルギー：$\frac{1}{2} \int_\Omega \mathbb{E}[\int_Q \bar{\rho}\bar{h} \, ds] \, dx$。
   • セル脱分極項：$\frac{1}{6} \int_\Omega |p(x)|^2 \, dx$（ここで $p(x)$ は有効分極である）。
2. 非局所極限: 非局所エネルギーは、マクロな静電場エネルギーから脱分極項を差し引いたものに収束する：
   • $\frac{1}{2} \int_{\mathbb{R}^3} |\nabla \phi|^2 \, dx - \frac{1}{6} \int_\Omega |p(x)|^2 \, dx$（ここで $\phi$ は有効分極 $p$ によって生成されるポテンシャルである）。
3. 全エネルギーの相殺: 極めて重要な結果は、セル脱分極項（$+\frac{1}{6}|p|^2$ および $-\frac{1}{6}|p|^2$）の相殺である。その結果、全連続体エネルギーは以下のようになる：
   $$E_{total} = \frac{1}{2} \int_\Omega \mathbb{E}\left[\int_Q \bar{\rho}\bar{h} \, ds\right] dx + \frac{1}{2} \int_{\mathbb{R}^3} |\nabla \phi|^2 dx$$
   これは、アンサンブル平均された微視的セルエネルギーと、マクロな場エネルギーの合成である。
4. 情報の保持: 局所極限と非局所極限は、異なる情報を保持する。非局所極限は決定論的な有効分極 $p(x)$ のみに依存し、微視的なゆらぎの具体的な統計量を失う。局所極限は、電荷密度のアンサンブル平均された統計量を保持する。
5. 検証事例:
   • 周期的な場合: 確率空間が単一の点であるとき、本フレームワークは既知の周期的な二スケール極限を回収する。
   • ランダム変位モデル: 著者らは、独立したセル単位のランダム変位モデルに対して仮定を検証している。彼らは、中心化された微視的なゆらぎが、有効分極や非局所エネルギーを変化させることなく、局所エネルギーを変化させ得る（アンサンブル平均された微視的項を変更することによって）ことを示している。

意義と主張
本論文は、エルゴード平均によって、静電エネルギーの微視的なランダム性への依存性がすべて排除されるわけではないことを主張している。長距離（非局所）相互作用は決定論的な有効分極によって完全に決定される一方で、短距離（局所）相互作用は、微視的な電荷分布の統計的特徴を保持する。

著者らは、有効分極のみに基づく連続体記述は、微視的なゆらぎのエネルギー的影響を捉えるには不十分であることを強調している。そのような記述は、マクロな場エネルギーとその境界への依存性を正しく予測するが、アンサンブル平均された微視的セルエネルギーを説明することはできない。この区別は、分極のゆらぎが脱分極場や表面スクリーニングと相互作用する有限温度の系において特に重要である。本研究は、ランダム媒体におけるこれらの効果を分離するための厳密な数学的基礎を提供し、従来の周期的な解析を、定常かつエルゴード的なランダム設定へと拡張するものである。