Molecular-Level Understanding of Polyethylene Glycol Adsorption on TiO2(001) at the Solid–Water Interfaces

技术摘要：聚乙二醇在 TiO₂(001) 固体–水界面上的分子级吸附理解

问题陈述
聚乙二经（PEG）改性是增强二氧化钛（TiO₂）纳米颗粒的胶体稳定性、水兼容性和界面功能性的标准策略。然而，预测氧化物/液态水界面的饱和 PEG 覆盖度仍然是一个重大挑战。这种困难源于 PEG 吸附受复杂的热力学平衡支配，涉及表面结合、链间堆积、水合作用、水分置换以及构象熵。虽然以往的研究已经探讨了 PEG 在锐钛矿 TiO₂(101) 晶面上的吸附，但锐钛矿 TiO₂(001) 晶面独特的表面原子结构和更高的反应活性可能会导致不同的界面行为。因此，需要对 TiO₂(001) 上 PEG 吸附进行分子层面的、热力学层面的理解，以便在水性环境中理性控制饱和覆盖度。

研究方法
本研究采用结合密度泛函理论（DFT）、从头算分子动力学（AIMD）和高斯过程回归（GPR）的计算框架，研究了水合锐钛矿 TiO₂(001) 上的 PEG-5 (n=5) 吸附。

• 模拟设置： 使用 VASP 软件包进行自旋极化 DFT 计算，采用 vdW-DF2 泛函以考虑色散相互作用，并对 Ti 3d 轨道使用 Hubbard U 校正（3.0 eV）。TiO₂(001) 表面使用具有三个 Ti–O 层的 (2 × 2) 原胞进行建模。
• AIMD 方案： 模拟在 298.15 K 的正则系综下进行，液态水密度约为 1.00 g/cm³。PEG 分子以垂直于表面的刷状构象进行初始化。系统经过平衡后，进行 10 ps 的 AIMD 轨迹模拟。每隔 1 ps 提取一个快照，在 0 K 下进行几何优化，并进行平均以获得代表性的能量。
• 覆盖度分析： 研究了覆盖度为 0.25、0.50、0.75 和 1.00 单层（ML）的 PEG-5 吸附情况。
• 自由能分解： 吸附自由能 ($\Delta G_{ads}$) 被分解为：
  1. 焓贡献： 端基吸附、链间相互作用以及 PEG–水溶剂化。
  2. 零点能 (ZPE)。
  3. 熵贡献： 采用先前关于 TiO₂(101) 研究中的熵值作为一阶近似，假设熵损失主要由构象限制而非特定晶面结合引起。
  4. H₂O 脱附： 基于 AIMD 中观察到的分子吸附水物种进行计算。
• 统计分析： 使用带有径向基函数核的高斯过程回归（GPR）来从离散的 AIMD 数据中插值得到连续的自由能曲线。采用自助重采样法（1,000 次迭代）来估计插值曲线的 95% 置信区间。

关键结果

• 结构演变： 在所有覆盖度下，PEG 链均呈现刷状构象（$R_F/D_m > 1.0$）。随着覆盖度从 0.25 增加到 1.00 ML，链间重叠增加，水分子被逐步从表面置换。在 0.75 和 1.00 ML 时，致密的刷状结构显著降低了水对 PEG 链的可及性。
• 焓贡献：
  • 端基吸附： 单位分子的归一化能量在 0.25–0.50 ML 时最为有利，但由于立体拥挤和 Ti 位点局部电子环境的变化，在高覆盖度下会减弱。
  • 链间相互作用： 吸引作用随覆盖度增加而增加，并在 0.75–1.00 ML 附近趋于饱和。
  • PEG–水溶剂化： 该项在高覆盖度下变得越来越不利。PEG–水氢键的总数在 0.50 ML 时达到峰值，并在 0.75 和 1.00 ML 时急剧下降，表明水被排除在致密刷状结构之外。
• 热力学稳定性：
  • 总吸附焓 ($\Delta H_{PEG}$) 在所有覆盖度下均保持有利，但在中间覆盖度（0.50 ML）处达到最大值，此时链间相互作用与水合作用达到了最佳平衡。
  • 总吸附自由能 ($\Delta G_{ads}$) 呈非单调变化。最有利的数值（$\Delta G_{ads} = -7.85$ eV）出现在 0.50 ML 处。
  • 在 0.75 和 1.00 ML 时，$\Delta G_{ads}$ 显著变得不再有利（在 1.00 ML 时上升至 -1.30 eV），这主要是由于与构象限制和溶剂化减少相关的熵损失增加所致。
• 微分吸附： 微分吸附自由能 ($\Delta G_{diff}$) 在 0.50 ML 以下为负值，表明增加更多 PEG 是热力学有利的。然而，$\Delta G_{diff}$ 在超过 0.50 ML 后变为正值（在 0.75 ML 时为 2.48 eV），表明进一步的吸附是不利的。
• GPR 插值： GPR 插值的自由能曲线证实了在 0.51 ML 附近存在极小值（插值值为 -8.11 eV），且在 0.50 ML 数据点周围具有较窄的 95% 置信区间，支持了离散 AIMD 的发现。

意义与主张
本文声称提供了在水性环境中理解 TiO₂(001) 上 PEG 饱和度的分子级热力学基础。主要发现是，饱和 PEG-5 覆盖度并非仅由 PEG–表面结合强度决定，而是由有利相互作用（表面结合、链间堆积、水合作用）与不利惩罚（水置换、构象熵损失）之间的竞争所决定的。

具体而言，研究预测在水合 TiO₂(001) 上的饱和 PEG-5 覆盖度为 0.50 ML（对应于 3.49 分子/nm² 的表面密度）。达到这一极限的原因是，超过此覆盖度后，由于 PEG–水氢键的丧失、溶剂化程度降低以及与致密刷状堆积相关的熵惩罚增加，热力学驱动力被抵消。这项工作表明，理性的 PEG 覆盖度控制需要平衡这些相互竞争的焓效应和熵效应，而不仅仅是最大化表面结合位点。