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

技术摘要：非均匀系统中的均匀分布

问题陈述
本研究探讨了经典密度泛函理论（DFT）在非均匀流体中的逆问题。标准的 DFT 应用旨在解决“正问题”——即在给定外部势能 $V(r)$ 的情况下确定平衡密度分布 $\rho(r)$——而本研究寻求其逆过程：构建特定的外部壁面势能 $V(r)$，以在靠近基底的流体区域内产生完全平坦且恒定的平衡密度分布 $\rho(r) = \rho_b$。

其动机源于观察到现实中的基底通常会诱导堆积振荡、耗尽层或吸附现象。作者旨在识别能够显式抑制这些非均匀性的“结构抵消”壁面场。该问题与润湿和干燥的界面理论相关，在这些理论中，当壁-流体相互作用与流-流相互作用满足特殊的匹配条件时，会出现平坦的密度分布。本研究调查了所需的补偿场如何取决于基底几何形状（平面、球形和柱形）以及流体相互作用。

方法论
作者采用 Rosenfeld 的基本度量理论（FMT）来描述硬球参考系统，并使用平均场处理吸引相互作用。研究分析了两种单组分流体模型：

1. 纯排斥硬球流体。
2. 具有硬球排斥及截断型 Lennard-Jones 类吸引尾部的流体。

研究考虑了三种基本几何形状：

• 平面： 位于 $z=R$ 处的平坦壁面。
• 球形： 半径为 $\varrho$ 的实心球体。
• 柱形： 半径为 $\varrho$ 的无限长实心圆柱体。

核心方法论是将预设的恒定密度分布代入由巨正则势泛函导出的 Euler-Lagrange 方程。外部势能通过下式恢复：
$$\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}$），验证了解析与数值构建的有效性。

意义与主张
本文声称提供了一种显式的微观实现方式来构建结构抵消壁面场，超越了形式化或纯数值的反演方法。其主要贡献在于：

• 显式解析构建： 推导了用于抵消界面结构的平面和球形几何下壁面势能的闭合形式表达式。
• 几何依赖性： 清晰地展示了基底曲率如何修正补偿场，特别强调了相比于纯排斥系统，曲率对具有吸引相互作用系统的影响更为强烈。
• 参考资料： 本工作提供了一套关于三种基本几何形状下加权密度和单体直接相关函数的紧凑公式集，可作为 DFT 实现的有用参考表达式。
• 理论背景： 研究表明，在显式的非局部密度泛函中，可以实现通常在润湿理论中被视为渐近或唯象匹配条件的平坦分布条件。

作者指出，补偿势能的范围是有限的，这取决于 FMT 核函数和截断吸引相互作用的范围。他们承认，在平均场理论之外，临界涨落会引入长程贡献，而这超出了当前有限程处理的范围。该工作建议未来的扩展方向包括多组分流体、系统的曲率展开，以及将这种逆向视角应用于动力学和流变学背景。