高分子玻璃化转变温度的分子动力学模拟：比体积法、MSD法与能量法的对比选择

聚甲基丙烯酸甲酯（PMMA）在320K附近出现了一段令人不安的温度平台——升温曲线上那长达50K的比体积平坦区，让OPLS-AA、PCFF和COMPASS三个力场的结果同时偏离了DSC实验标定的玻璃化转变温度值378K。这个偏差不是”模拟不准”四个字能概括的。它指向一个更根本的问题：当冷却速率被迫维持在10^10 K/s量级时，高分子链段松弛的时间窗口被硬生生压缩了十几个数量级，而比体积法对这种压缩的敏感度远高于预期。

分子动力学模拟中玻璃化转变温度的判定历来有三套判据：比体积-温度曲线拐点、均方位移（MSD）斜率突变、以及非键能与温度的交叉。问题在于，这三套判据在同一组轨迹上给出的Tg值很少重合。Han等人在2016年对聚苯乙烯体系的系统评估中已经注意到这一现象——三种方法的结果离散范围可达25K，但该研究未深入探讨冷却速率与力场选择的交互效应。在PMMA这一侧链更复杂、分子间相互作用更强的体系中，这种离散可能被进一步放大。

比体积法的困境最为直观。从500K以每纳秒降1K的速率冷却至200K的过程中，PMMA的比体积-温度曲线在340K到370K区间呈现出渐进而非陡峭的弯折——这不是一个”拐点”，而是一段模糊的过渡带。模拟盒子的初始尺寸为6.5 nm，含48条聚合度为50的无规立构PMMA链，NPT系综下采用各向同性的Nosé-Hoover控压器。当用双线性拟合强行定位时，分段区间的选取偏差50K能让Tg浮动15K。更棘手的是，当冷却速率从1 K/ns提升到5 K/ns时，比体积法测得的Tg整体上移了约12K，这种速率依赖性比实验中的DSC扫描速率依赖性大两个数量级，根源在于模拟时间尺度内链段无法完成构象弛豫。Williams-Landel-Ferry方程的外推表明，要将Tg偏移控制在2K以内，冷却速率需降至10^7 K/s以下——这在当前原子级MD的计算资源下几乎不可实现。

于是转向MSD法。在温度扫描过程中，计算每个温度窗口内链段质心的均方位移并拟合扩散系数D，期望D-T曲线在Tg处出现斜率转折。实际操作中发现，当温度低于Tg时，MSD在10ns的模拟窗口内仍被限制在0.5 nm²以下——这相当于链段质心仅移动了不到两个单体单元的距离——统计噪声彻底淹没了扩散行为的信号。OPLS-AA力场下D的拟合不确定性在玻璃态区间高达40%，意味着MSD法判定的Tg携带至少±10K的模糊边界。PCFF力场表现稍好：其非键参数中范德华半径的设置对链段堆积的描述更紧致，MSD平台值更低且起伏更小，D-T转折更为清晰，但最终给出的Tg为363K，与实验值仍有15K的差距。这一差距的物理根源在于，PCFF的力场参数是从小分子晶体和液体数据拟合而来，其对高分子长链协同运动的描述天生存在系统性的刚度偏高倾向。

能量法提供了第三条路径。跟踪不同温度下非键相互作用能（范德华+静电）随温度的变化，在Tg附近出现斜率突变——因为玻璃态和橡胶态的热膨胀系数差异导致了平均分子间距的分叉，进而改变了非键能的温度依赖关系。在PMMA体系中，非键能法给出的Tg为370K，与实验值378K最为接近。但这种接近带有运气成分：COMPASS力场的静电参数针对甲基丙烯酸酯基团的拟合恰好补偿了冷却速率带来的系统偏移。换用其他含酯侧链的高分子体系——例如聚甲基丙烯酸乙酯或聚甲基丙烯酸丁酯——这种补偿未必成立。Kraft等人在2018年对聚碳酸酯的分子动力学模拟中同样报告了非键能法与实验值的意外吻合，并指出这可能是力场参数化过程中隐式地吸收了某些实验热力学数据的温度依赖特征所致。

三种方法的交叉比对最终揭示了一套具有操作性的决策框架。比体积法操作最简单，物理图像直接，适合快速筛选力场和冷却方案，但对速率和分段区间的敏感性要求至少三个冷却速率下的重复计算，并通过速率外推来逼近准静态极限。MSD法对力场的堆积参数和链段摩擦特性高度敏感，适合用来检验力场在玻璃态的动力学表现，但低于Tg时的统计收敛需要将模拟窗口延长至50ns以上，计算成本也因此膨胀到比体积法的三到五倍。非键能法在三种方法中与实验值的偏差最小，却对力场的非键参数存在隐式依赖——它的”准确”可能源于误差抵消而非物理描述的完备性，不建议作为唯一的Tg判据使用。

这一轮跑下来的教训比论文里的结论多得多。Tg的分子动力学模拟不是在”算一个数”，而是在测量一个在模拟时间尺度上根本来不及完成的物理过程。三种判据的差异不是方法的优劣之分，而是各自捕捉了玻璃化转变的不同侧面——比体积法捕捉的是宏观自由体积的冻结，MSD法捕捉的是链段扩散动力学的冻结，能量法捕捉的则是分子间相互作用的温度微分。将三者并列、交叉验证，而非择一而从，才是分子动力学模拟中玻璃化转变温度判定的可靠路径。

欢迎访问 keyanxueshu.com 了解分子动力学计算服务。