小分子自组装动力学模拟

结合案例图分析，本文以这个两亲分子胶束形成项目为例，讲自组装模拟的采样与力场坑。这是一个两亲性小分子胶束形成研究项目。客户想预测临界胶束浓度和胶束尺寸，以指导实验合成。第一张案例图是最大簇尺寸随时间变化，约 60 ns 后趋于平稳达约 50 个分子。第二张案例图是平衡后聚集数分布直方图，峰值在 50 附近。

项目过程：两亲分子胶束形成项目的执行链路

方法上粗粒化（CG）能跑更长但丢原子细节、且力场参数需额外验证；客户要的是具体胶束尺寸和形貌，我选了全原子显式水 MD，虽然慢但结果可直接和实验散射对比。力场参数对两亲分子至关重要——很多没现成参数，我先用 RESP 从量化单点拟合电荷、再用 CGenFF 补二面，验证过液相密度和实验偏差 < 3%。 我的模拟：构建含 150 个小分子的溶液盒，加水和离子，2 fs 时间步、300 K、长时间 MD，通过团簇尺寸和分布判断自组装。客户据此调整亲水/疏水链长。

案例分析：从团簇曲线到胶束尺寸

图里读出的数要落地。第一步确认已到平衡——最大簇尺寸平台且分布直方图单峰。我对平衡后 20 ns 的团簇尺寸算均值和标准差，50±5 个分子且标准差小，说明胶束尺寸均一；若标准差大说明在球/棒间动态翻转。 第二步把形貌（球/棒/层）从快照判定，本例球状、直径约 4 nm。误差量级上，力场未拟合时胶束尺寸偏大 30%、形貌从球变成棒，拟合后回到 50±5 的球状。第三步结合 CMC：我只模拟了单一浓度，在报告里标注”本浓度下胶束尺寸 50±5，CMC 需另行计算”。

总结教训：力场参数与采样时间的坑

做的过程中我踩过实打实的坑。第一坑是采样时间：我画过最大簇尺寸随模拟时间的曲线——20 ns 时只聚到 15 个分子（看似没自组装），40 ns 到 35，60 ns 才稳定在 50±5，证明只跑 20 ns 会误判”不自组装”，采样必须到平衡。 第二坑是力场参数，初版直接用通用力场没拟合，胶束尺寸和形貌都错。第三坑是团簇判定 cutoff（0.5 nm），cutoff 取太大把溶剂分子也算进团簇，分布峰被拉宽。

项目经验：自组装模拟交付的可复用清单

作为可复用的项目经验，第一，力场参数对小分子至关重要，没现成必须拟合验证（液相密度偏差 < 3%）。第二，除了胶束尺寸，把形貌和形成动力学也描述清楚，对实验才有指导意义。第三，采样必须到平衡（本例 60 ns），否则尺寸只是运气数字。客户据此调整亲水/疏水链长。说到底，自组装模拟最怕”跑了就交”——没到平衡、参数没验证，给出的胶束尺寸只是运气数字。

小分子自组装动力学模拟

计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验

需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。