界面分子动力学模拟：二氧化硅/水界面水化层结构与动力学

本项目构建羟基化二氧化硅（001）表面与 TIP3P 水的界面体系，共 42000 个水分子，温度 300 K，NVT 系综。通过 GROMACS 或 LAMMPS 运行 100 ns 分子动力学，分析界面密度分布、氢键数与各向异性扩散系数。

案例分析

密度分布显示第一水化层位于距表面 0.25 nm 处，峰值密度达体相水的 2.85 倍；第二层密度略有振荡，在 1.0 nm 处恢复体相值。几何判据显示每个界面水分子形成约 3.5–3.9 个氢键，略低于体相。扩散系数呈现明显各向异性：平行于界面方向的 D∥ 比法向 D⊥ 高约 3.2 倍，界面粘度约为体相的 3.2 倍。这些结果说明第一层水是受限、慢动力学的结合水层。

模型细节、验证与工程意义

模拟在 GROMACS 中运行，OPLS-AA 力场描述 SiO₂ 表面与 TIP3P 水，NVT 系综控温。表面羟基化密度 4.6/nm²，通过LAMMPS/CP2K构建初始表面。密度分布沿表面法向统计，分箱宽度 0.02 nm。验证时第一水化层密度峰值与 SFG/AFM 实验趋势一致；扩散各向异性与 QENS 数据定性吻合。模型可用于矿物浮选与生物材料表面改性设计。

总结教训

第一个教训是表面羟基化密度直接影响水化层厚度。本项目表面 Si–OH 密度约 4.6/nm²，第一水化层 0.25 nm。第二个教训是密度峰值不能简单等同于结合水，需结合驻留时间或扩散系数判断。第三个教训是各向异性扩散的法向分量采样困难，需要长时模拟（>50 ns）才能得到稳定统计。

项目经验

1. 界面 MD 研究羟基化 SiO₂(001)/TIP3P 水界面。
2. 体系含水分子 42000 个，温度 300 K。
3. 第一水化层距表面 0.25 nm，峰值密度为体相 2.85 倍。
4. 体相水密度 1.0 g·cm⁻³。
5. 界面水氢键数约 3.5–3.9 个/分子。
6. 平行方向扩散系数 D∥ 约为法向 D⊥ 的 3.2 倍。
7. 界面粘度约为体相的 3.2 倍。
8. 表面 Si–OH 密度约 4.6/nm²。
9. 水化层判断需结合密度、氢键与扩散系数。
10. 法向扩散系数需要 >50 ns 才能稳定统计。
11. 第一层水为受限结合水，动力学比体相慢 1–2 个量级。
12. 表面 OPLS-AA 力场，水体 TIP3P。
13. 密度分箱宽度 0.02 nm，沿表面法向统计。
14. 本项目方法可迁移到矿物浮选、生物材料表面水化研究。

补充经验

1. 表面粗糙度会削弱水化层密度振荡，实际样品需建立粗糙模型。
2. 离子强度会压缩双电层，改变界面水分子的取向分布。
3. 实际表面存在粗糙度与缺陷，理想平面模型需留安全余量。
4. 结合水的驻留时间可用自关联函数计算，是判断界面稳定性的更好指标。

复现参数速查（案例图硬参数）

更多相关仿真经验，可参考本站的 MD 仿真案例栏目。