高分子动力学模拟：链长、温度和缠结——三个变量交织成Tg和扩散系数的十度偏差

高分子的动力学是一个时间尺度的阶梯——键振动（飞秒级）→侧基旋转（皮秒级）→链段运动（纳秒级）→整链扩散（微秒到毫秒级）→解缠结和爬行（毫秒到秒级）。全原子MD能覆盖前三级——键振动到链段运动。再往上走，需要对系统做不同程度的”加速”。

聚苯乙烯（PS）是一个经典例子。10条PS短链（每条50个单体，分子量约5200 g/mol）在500K下的全原子模拟，跑200ns能看到整链扩散——链中心均方位移在log-log图上斜率接近1。但分子量5.2万的PS（每条500个单体），在同样的全原子框架下跑200ns，链中心基本在原地——看不到扩散，因为运动的时间尺度已经跳到了微秒以上。

降温速率：Tg不是材料的固有属性

高分子动力学模拟中最常被问的一个量是Tg（玻璃化转变温度）。MD求Tg的标准做法是：从高温（>Tg + 100K）开始，以恒定速率降温，每降一步温度做一段NPT平衡，记录比体积随温度的变化。比体积-温度曲线的两个线性段交点即为Tg。

但MD的降温速率通常在10⁹-10¹⁰ K/s量级——实验DSC的降温速率是10 K/min，即0.17 K/s。差了10⁹倍。在这个量级的速率差下，MD的Tg比实验Tg系统性偏高。

全原子PS（50-mer）以10¹⁰ K/s降温，Tg~420K。实验PS的Tg约373K（对于这个分子量，Tg比无限分子量低约10K）。高了将近50K。

原因是：降温速率越快，链段在某个温度下”冻住”的时间越晚——需要更低的温度才能使弛豫时间超过实验观测窗口。Tg和降温速率的对数关系（Vogel-Fulcher-Tammann行为）决定了每10倍降温速率变化对应约3-5K的Tg移动。从10¹⁰到0.17 K/s的10¹⁰倍差距，理论上的Tg偏移应该在30-50K左右——与MD的结果一致。

Tg不是材料的一个固有常数——它是弛豫动力学在特定时间窗口下的一个截断。MD的Tg偏高不是”算错了”，是时间窗口不同。

链长效应和Rouse模型的适用边界

对于50-mer的短链PS，链运动在500K下遵循Rouse模型——链中心MSD ~ t¹/²（未缠结链的亚扩散行为），扩散系数 D ~ 1/N（反比于链长）。拉伸链长到500-mer（分子量5.2万，超过了PS的缠结分子量Mc≈1.8万），链的运动模式从Rouse变成了管模型（tube model/reptation）。

全原子模拟500-mer链的扩散已不可行。粗粒化方案——用前面那篇文章里的双珠CG模型——把500-mer映射到1000个CG珠子，时间步长10fs，在800K下跑50μs（等效时间，考虑加速因子4.2后约12μs真实时间）。CG模拟能观察到链在缠结网络中的爬行运动——MSD呈现t¹/⁴（管约束下的Rouse运动）到t¹/²（经典爬行，即reptation）的交叉，交叉时间约1.2μs。

这个交叉时间在全原子模拟中根本看不到——全原子能覆盖的最大模拟时间200ns，只能看到链在管中的局部Rouse运动，看不到爬行。

加速因子和标定

CG模型的动力学加速因子随链长也不是常数。50-mer的加速因子~4.2（前面文章已标定），但500-mer的加速因子实测约5.8——更长的链在CG模型中感觉到的有效的”摩擦”更小。部分原因是CG势的软性减少了链之间的粗糙摩擦（atomistic-level的几何互锁在CG珠子层面被抹掉了）。

这意味着不能把短链标定出来的加速因子直接套用给长链。每个链长和每个温度都应该独立标定。标定的基准是：CG和全原子在重叠的时间和空间尺度上（如果有的话）的对应关系。

高分子动力学模拟的终极困境是：要做准Tg，降温速率太快的偏差无法消除（除非换增强采样方法）；要做准链扩散，全原子的时间尺度不够，CG又丢失了化学细节。两种不完美之间的平衡——是高分子模拟永远的现实。