A route to the thermodynamics of colloid-polymer mixtures from structural information

技术摘要：从结构信息建立胶体-聚合物混合物热力学路径的研究

问题陈述
液态理论传统上依赖于微观结构（例如径向分布函数 $g(r)$ 和静态结构因子 $S(k)$）与宏观热力学性质之间的严密联系。然而，预测复杂流体（如胶体-聚合物混合物）的热力学和相行为面临着重大挑战。传统方法通常需要明确知道底层的相互作用势。在许多实验和复杂的软物质系统中，这些有效相互作用是未知的、取决于状态点的，或者受到多体效应的强烈影响，从而使得简单的对势近似（如平均力势）在有限密度下不再准确。此外，获得准确的热力学性质通常需要对广泛的状态点进行热力学积分（TI），这在计算上是非常耗时的。现有的逆向方法（如逆向玻尔兹曼迭代，IBI）可以从结构数据中重构势能，但生成的势能往往保留了对特定热力学状态的依赖性，限制了其对更广泛相行为的预测能力。

方法论
作者开发了一个框架，旨在直接从单个热力学状态点获得的结构相关性来预测热力学和相行为，从而绕过了对显式相互作用势的需求。该方法利用径向分布函数（$g(r)$）或静态结构因子（$S(k)$）作为输入来构建自由能泛函。

该形式区分了两类相互作用：

1. 排斥体积相互作用： 具有硬核排斥的系统（其中当 $r < \sigma_u$ 时，$g(r)=0$）。
2. 软相互作用： 粒子可以相互穿透的系统（其中 $g(r)$ 保持有限值）。

对于具有排斥体积相互作用的系统，总相互作用势被分解为一个短程排斥核（$\phi_{ref}$）和一个剩余尾部（$\phi_{tail}$）。该工作的创新之处在于确定了有效硬球直径（$\sigma_{opt}$）。作者并非使用基于已知势能的标准判据（如 Barker-Henderson 方法），而是通过最小化参考硬球 RDF 与由 IBI 重构势（仅作为隔离排斥贡献的中间步骤）导出的 RDF 之间的差异来确定 $\sigma_{opt}$。

随后，过剩自由能（$F^{exc}$）被拆分为参考贡献（$F^{exc}_{ref}$）和尾部贡献（$F^{exc}_{tail}$）：

• 参考系统： $F^{exc}_{ref}$ 使用基本测度理论（FMT）或 Carnahan-Starling 状态方程计算有效硬球流体。
• 尾部贡献： 作者提出了两种互补的途径来确定代表尾部的积分相互作用强度（$A$）：
  1. 热力学积分（TI）途径： 假设结构相关性由排斥体积效应主导，通过耦合平均的配对相关函数来近似全系统 $g(r)$。
  2. 直接相关函数（DCF）途径： 利用通过 Ornstein-Zernike (OZ) 关系获得的配对直接相关函数 $c^{(2)}(r)$。过剩自由能表示为全系统的 DCF 与参考系统 DCF 之差（$\Delta c^{(2)}$）。

对于软相互作用（例如聚合物），该框架利用随机相位近似（RPA）关系，将积分强度直接与 DCF 或第二维里系数联系起来。

该框架通过构建所有配对类型的积分相互作用强度矩阵（$A_{ij}$）扩展到了多组分混合物（胶体-聚合物系统）。参考自由能考虑了硬球相关性、耗尽效应（通过自由体积理论）以及软聚合物-聚合物相互作用，并将尾部贡献以微扰方式加入。

主要贡献

• 单状态点预测： 主要贡献在于提出了一种方法，仅通过来自单个状态点的结构数据即可推断出状态方程和相行为，而无需进行广泛的热力学积分或预先知晓相互作用势。
• 无需势能的逆向框架： 该方法避免了重构状态依赖的有效势，而是直接从结构相关性构建自由能泛函数，仅将 IBI 作为定义参考系统几何结构的过渡工具。
• 路径的鲁棒性： 研究表明，DCF 途径和 TI 途径在确定积分相互作用强度方面能够得出一致的结果，验证了使用源自 $c^{(2)}(r)$ 或 $g(r)$ 的结构核进行自由能计算的有效性。
• 对混合物的适用性： 该方法成功处理了二元胶体-聚合物混合物，涵盖了自相关和交叉相关，包括具有非理想（软）聚合物相互作用的系统。

结果
该框架针对以下分子模拟数据进行了验证：

1. 单组分系统： 对 Lennard-Jones (12-6) 和广义 Mie (24-12) 势的测试表明，TI 途径和 DCF 途径预测的状态方程在低密度到高密度范围内与模拟数据高度吻合。基于稀释极限近似的维里途径在较高密度下表现较差。预测的 12-6 系统双节线曲线与模拟趋势一致，尽管临界点略有低估，这可能是由于势能截断差异所致。
2. 胶体-聚合物混合物： 应用于具有高斯聚合物势和 Mie 胶体势的二元混合物。该框架准确预测了固定胶体组分下的状态方程（压力 vs 密度）。
3. 相行为预测： 该方法成功预测了胶体-聚合物混合物的流体-流体共存（双节线）。由此得出的相图（仅由单个输入状态的结构相关性导出）正确识别了向聚合物富集/胶体贫乏相和胶体富集/聚合物贫乏相的偏析现象，即使在聚合物-聚合物相互作用为非理想的情况下也是如此。

意义与主张
本文声称建立了一条直接从实验测得的结构信息中推断相行为的有效途径，即使在微观相互作用并非预先已知的情况下亦是如此。这对于有效相互作用难以精确确定且随系统状态变化的实验软物质系统尤为重要。

作者强调，其公式化方案为构建无需显式相互作用势或特定系统热力学模型的自由能泛函数提供了一条通用路径。作者也谦逊地指出，目前的表述限于各向同性相互作用，并未明确考虑取向自由度或强相关系统（例如高电荷胶体），在这些系统中仅靠配对相关性可能是不够的。然而，他们认为该框架仍适用于结构相关性可由各向同性分量描述的各向异性系统。这项工作表明，单个状态点的结构信息包含足够的信息来预测复杂流体的更广泛热力学响应，为计算昂贵的模拟或不确定的势模型提供了一种实用的替代方案。