Uniform distributions in nonuniform systems: Wall potentials generating constant density profiles in classical density functional theory

技術要約：非一様系における一様分布

問題提起
本研究は、不均一流体における古典密度汎関数理論（DFT）における逆問題を扱う。標準的なDFTの応用は、与えられた外部ポテンシャル $V(r)$ に対して平衡密度プロファイル $\rho(r)$ を決定するという「直接問題」を解くものであるが、本研究ではその逆、すなわち、基板に隣接する流体領域において、完全に平坦で一定の平衡密度プロファイル $\rho(r) = \rho_b$ を生成するために必要な、特定の外部壁ポテンシャル $V(r)$ を構築することを目的とする。

この動機は、現実的な基材が通常、パッキング振動、枯渇層、あるいは吸着を誘起するという観察に基づいている。著者らは、これらの不均一性を明示的に抑制する「構造相殺型」の壁場（wall field）を特定することを目指している。この問題は、壁－流体間および流体－流体間の相互作用が特別な整合条件を満たす場合に平坦な密度プロファイルが生じるという、濡れや乾燥の界面理論とも関連している。本研究では、要求される補償場が基材の幾何学的形状（平面、球状、円筒状）および流体の相互作用にどのように依存するかを調査する。

手法
著者らは、硬球参照系を記述するためにRosenfeldの基本測度理論（FMT）を用い、引力的相互作用については平均場近似による扱いを採用している。2種類の単成分流体モデルを解析している：

1. 純粋に反発的な硬球流体。
2. 硬球の反発に加えて、切断されたレナード・ジョーンズ型の引力的テイルを持つ流体。

研究では、以下の3つの基本的な幾何学的形状を考慮している：

• 平面（Planar）: $z=R$ にある平らな壁。
• 球状（Spherical）: 半径 $\varrho$ の固体の球。
• 円筒状（Cylindrical）: 半径 $\varrho$ の無限に長い固体の円筒。

核心となる手法は、グランドポテンシャル汎関数から導出されるオイラー＝ラグランジュ方程式に、規定された一定の密度プロファイルを代入することである。外部ポテンシャルは以下のように回収される：
$$\beta V(\xi) = \beta\mu + c^{(1)}(\xi) - \ln[\Lambda^3 \rho(\xi)]$$
ここで、$c^{(1)}$ は一次の直接相関関数である。著者らは、各幾何学的形状における重み付き密度と $c^{(1)}$ の明示的な式を導出している。

• 解析的処理: 平面および球状の幾何学においては、FMTの畳み込み積分が一次元の形式に還元され、閉じた解析解を可能にする。
• 数値的処理: 円筒状の幾何学においては、FMTカーネルの還元により完全楕円積分が生じるため、単純な閉じた形式を得ることができず、したがってこのケースは数値積分によって評価される。

主な結果

1. 平面幾何学: 逆問題は、壁ポテンシャルのための明示的かつコンパクトな解析的表現を与える。硬球の場合、ポテンシャルは純粋に反発的で短距離であり、その振幅は充填率 $\eta$ とともに増加する。引力的相互作用が含まれる場合、ポテンシャルは引力成分を持つようになり、これはパッキングと凝集性の間の競争を反映している。接触値のポテンシャルは、厳密な総和則 $\beta p_{hs}/\rho_b = 1 + \beta V_{hs}^{pl}(R)$ を満たす。
2. 球状幾何学: 解析解は存在するが、$\eta$ と基材半径 $\varrho$ の両方に依存するため、大幅に複雑である。結果は、純粋な反発的硬球ポテンシャルの場合、曲率の影響は弱いものの、引力的相互作用を持つ流体においてはより顕著な影響を与えることを示している。具体的には、球の半径が減少するにつれて、引力の井戸が浅くなる。
3. 円筒状幾何学: 壁ポテンシャルは数値的に決定された。結果は、同じ半径に対して、円筒のポテンシャルが平面と球状の極限の中間に位置することを示しており、これは円柱の中間的な曲率を反映している。
4. 検証: 導出されたポテンシャルを用いた独立したDFT計算により、規定された平坦な密度プロファイルが数値的な精度内（相対偏差は $10^{-4}$ のオーダー）で再現されることが確認され、解析的および数値的な構成が検証された。

意義と主張
本論文は、形式的または純粋に数値的な反転を超えて、構造相殺型壁場の明示的な微視的実現を提供すると主張している。主な貢献は以下の通りである：

• 明示的な解析的構成: 界面構造を相殺する壁ポテンシャルの閉じた形式の表現を、平面および球状の幾何学において導出した。
• 幾何学的依存性: 基材の曲率がいかに補償場を変化させるかを明確に示しており、特に、純粋な反発系と比較して、引力的相互作用を持つ系において曲率依存性が強いことを強調している。
• リファレンス資料: 本研究は、3つの基本幾何学的形状における重み付き密度および一次直接相関関数のコンパクトな公式集を提供しており、DFT実装における有用な参照式として機能する。
• 理論的文脈: 本研究は、濡れの理論において漸近的または現象論的な整合条件として扱われることが多い平坦プロファイル条件が、非局所的な密度汎関数論の中で明示的に実現可能であることを確立している。

著者らは、補償ポテンシャルの範囲は有限であり、FMTカーネルおよび切断された引力的相互作用に結びついていることを注記している。平均場理論を超えると、臨界ゆらぎが長距離の寄与を導入する可能性があるが、これらは現在の有限範囲の取り扱いの範囲外である。本研究は、多成分流体、系統的な曲率展開への拡張、およびこの逆の視点を動力学やレオロジーの文脈へ適用することを示唆している。