高分子动力学模拟：链构象标度律与扩散区制分析

本项目采用粗粒化分子动力学模拟一条 500 珠的高分子链在良溶剂与 θ 溶剂中的动力学。力场为 Kremer-Grest 珠簧模型，截断半径 2.5σ，模拟温度分别为良溶剂条件与 θ 条件。分析回转半径 Rg、均方位移（MSD）与末端距分布。

案例分析

良溶剂下，回转半径随链长标度指数 ν=0.588，与 Flory 理论 0.60 接近，对应溶胀无规线团。θ 溶剂下 ν=0.500，呈理想链行为。MSD 三段区制清晰可见：短时的 Rouse 区指数约 0.52、中间的 entanglement 平台区、长时的 Fick 扩散区。末端距分布峰值对应最可几末端距 3.11 nm，持久长度 0.56 nm。溶剂可及表面积随时间演化显示链逐渐塌缩/溶胀达到稳态。

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

模拟在 GROMACS 2022 中运行，力场参数取自 Kremer-Grest 经典设定，积分步长 0.005τ，Langevin 控温。体系包含单条 500 珠链，盒子大小保证浓度低于重叠浓度。Rg 与末端距统计在 200 个独立快照上完成。验证时标度指数与 Flory-Fisher 理论对比，ν 偏差 <0.02，说明力场与采样充分。

总结教训

第一个教训是标度指数对模拟时长敏感。Rouse 区至少需要 10τ 才能稳定分辨。第二个教训是截断半径 2.5σ 对 KG 模型是标准值，改大会增加计算量但提高范德华精度。第三个教训是良溶剂与 θ 溶剂的区分主要靠 LJ 势中 ε 的调节，需通过第二维里系数标定。

项目经验

1. 高分子 MD 采用 Kremer-Grest 粗粒化珠簧模型。
2. 链长 500 珠，截断半径 2.5σ。
3. 良溶剂下 Rg 标度指数 ν=0.588。
4. θ 溶剂下 Rg 标度指数 ν=0.500。
5. Rouse 区 MSD 指数约 0.52。
6. 末端距最可几值 3.11 nm，持久长度 0.56 nm。
7. MSD 三段区制：Rouse 区、平台区、Fick 扩散区。
8. Rouse 区需 ≥10τ 才能稳定分辨。
9. 溶剂可及表面积演化可判断链构象是否达到稳态。
10. 良溶剂/θ 溶剂通过调节 LJ ε 实现，需用第二维里系数标定。
11. KG 模型截断 2.5σ 是兼顾精度与效率的标准选择。
12. 积分步长 0.005τ，Langevin 热浴控温。
13. 200 个独立快照用于统计 Rg 与末端距。
14. 本项目方法可迁移到聚合物溶液流变、自组装动力学研究。

补充经验

1. 全原子力场可补充侧链细节，但计算成本比粗粒化高 1–2 个量级。
2. 缠结效应在链长 >200 珠时逐渐显现，长链需用滑动弹簧或拓扑约束。
3. 不同溶剂质量可通过调节 LJ ε 参数系统扫描。
4. 温度对标度指数的影响需通过溶剂质量参数系统扫描。

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

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