水凝胶分子动力学模拟

我做过一个水凝胶网络的分子动力学项目，客户想定量看交联密度怎么影响溶胀比和水分子扩散。我这个项目把高分子交联网络、溶胀行为和输运性质串在了一起。

左侧图是我构建的水凝胶细观模型。粉色波浪线代表聚合物链，橙色方块是交联点，浅蓝色小点是水分子。水凝胶的本质是亲水性高分子链通过化学或物理交联形成三维网络，水分子填充在网格里。我做 MD 时，先搭一条聚乙烯醇或聚丙烯酸链，按目标交联密度加入交联剂，然后把体系放到水盒子里溶胀平衡。

右侧图是溶胀比 Q 和水扩散系数 D 随交联密度变化的关系。横轴是交联密度（mol%），纵轴是相对值。红色虚线 Q 随交联密度升高而下降，因为交联点越多，网络越紧密，能容纳的自由水越少；蓝色实线 D 也下降，说明水分子在致密网络里迁移更慢。图中数据点显示：交联密度从 1 mol% 升到 8 mol%，Q 从约 6.2 降到 2.2，D 从约 3.5×10⁻⁶ cm²/s 降到 0.5×10⁻⁶ cm²/s。这个趋势和客户做的重量法溶胀实验定性一致。

模拟流程：我用 LAMMPS 建初始构型，聚合物用 OPLS-AA 或 GAFF 力场，水用 TIP3P。先做 NPT 下溶胀：固定聚合物网络，让水分子自发填充并平衡；然后放开全部自由度跑 NPT 生产 50 ns。溶胀比通过对比干网络和溶胀后的体积得到，水扩散系数通过 MSD 拟合。

踩坑经验：第一个是交联反应的实现。MD 不能直接做化学交联，我通常先按目标交联度人工指定交联点对，再用约束或特殊键参数模拟化学键。第二个是溶胀平衡时间：高分子网络弛豫很慢，50 ns 可能还不够，我通过监测密度和 Rg 判断，必要时候跑到 200 ns。第三个是水含量：初始水分子数要足够多，让体系能充分溶胀；太少会人为限制 Q。

工程经验：水凝胶设计是溶胀和力学性能的折中。交联密度低，溶胀大、水输运快，但强度差；交联密度高，强度高，但溶胀小、药物释放慢。做 MD 时，我建议同时输出杨氏模量和溶胀比，给客户一个完整的性能窗口。另外，如果客户研究的是温敏或 pH 敏水凝胶，还需要在力场里加入响应性单体参数。

更多 GROMACS/LAMMPS 水凝胶和软物质模拟的案例，我整理在 [水凝胶分子动力学模拟](https://www.keyanxueshu.com/category/md/) 栏目里。[水凝胶分子动力学模拟](https://www.keyanxueshu.com/)

溶胀行为还可以从“自由水”和“结合水”的角度再解读。通过分析水分子与聚合物链的氢键寿命，可以把水分为紧密结合水（寿命 > 100 ps）、松散结合水（10-100 ps）和自由水（<10 ps）。高交联密度下，紧密结合水比例上升，自由水比例下降，宏观上表现为溶胀比降低和扩散系数下降。我们做报告时，会把水分子分类比例和右侧 Q、D 曲线一起呈现。另一个工程上有用的量是溶胀前后的网络孔径分布。孔径越大，药物或离子输运越快，但机械强度越差。这个信息可以从 RDF 或 Voronoi 分析得到，适合给客户做设计权衡。