2. 聚合物链结构与力场基础

聚合物分子动力学（MD）模拟的第一步，是把一条真实的高分子链翻译成“力场能读懂的语言”。这一步决定了模拟能否复现实验中的链尺寸、玻璃化转变温度、密度、模量乃至自组装形貌。力场本质上是势能函数 U(r) 及其参数集合，总能量通常拆成键合项与非键合项两部分：

U_total = U_bond + U_angle + U_dihedral + U_improper
        + U_vdw + U_coulomb

键合项描述沿链骨架的成键、键角与二面角扭转；非键合项描述原子对之间的范德华（LJ 12-6）与静电作用。聚合物与普通小分子的关键差别在于：链由成千上万个重复单元通过共价键连接，构象熵主导行为，且同一体系往往同时存在局部键合细节与长程链段松弛，因此力场的“分辨率”选择成为首要工程决策。

三种分辨率力场的适用边界

全原子（All-Atom, AA）力场为每个原子（含氢）显式建模，能描述氢键、侧基取向、手性、局部堆积等细节。典型代表有 PCFF、COMPASS、OPLS-AA、CHARMM CGenFF、AMBER GAFF。它适合研究：小到中等链长（聚合度几十到几百）的局部结构、界面吸附、氢键网络、玻璃态局部动力学。代价是自由度大、时间步长受 C–H 伸缩振动限制（约 1 fs），难以覆盖微秒以上的链松弛。

联合原子（United-Atom, UA）力场把 CH₂、CH₃ 等基团合并为单个作用位点，氢原子被隐式处理。代表有 TraPPE-UA、NERD、PYS。它保留了链的拓扑与二面角构象统计，但省去了高频 C–H 振动，时间步长可放宽到 2–5 fs，适合研究聚乙烯、聚丙烯等烯烃类聚合物的熔体密度、链尺寸与相行为。其边界在于：涉及极性侧基、氢键或精确静电的体系，UA 的简化会带来明显偏差。

粗粒化（Coarse-Grained, CG）力场把若干单体映射为一个珠子，如 Martini、SDK、MARTINI 聚合物模型。它牺牲化学细节换取时间与空间尺度，可模拟自组装、相分离、膜穿透、纳米颗粒-聚合物界面等介观现象。CG 的适用边界是：不能直接给出原子级结构信息，且参数依赖映射方案与目标性质，需通过“自下而上”或“自上而下”策略标定。

键合项参数化：从量子化学到构象统计

键合项的核心是让模型链在局部构象上与真实分子一致。以聚乙烯为例，C–C 键伸缩可用简谐势 k_b(r−r₀)²，r₀≈1.53 Å，k_b 由 DFT 势能面曲率得到；键角弯曲同理，C–C–C 平衡角约 112°。真正决定链柔性的是二面角势：反式（t）、左旁式（g⁺）、右旁式（g⁻）三个能量极小值之间的能垒与能量差。TraPPE-UA 中聚乙烯的二面角采用 OPLS 形式：

U_dihedral = c1(1+cosφ) + c2(1−cos2φ) + c3(1+cos3φ)
c1 = 1.19 kJ/mol, c2 = −0.40 kJ/mol, c3 = 2.06 kJ/mol

这组参数使 t 态最稳定，g 态能量高约 2.1 kJ/mol，与量子化学结果吻合。参数化流程通常是：先用 DFT（如 B3LYP/6-31G*）扫描 φ 从 0° 到 360°，得到相对能量曲线；再用最小二乘拟合上述余弦级数；最后在 MD 中统计 t/g 比例，与实验或高精度计算对比。若 g 态比例偏高，链会过度柔性，导致回转半径 R_g 偏小。

非键合项参数化：密度与内聚能的平衡

非键合项决定链间堆积与热力学性质。LJ 参数 ε、σ 通常通过拟合实验密度和蒸发焓获得。以聚苯乙烯为例，全原子 OPLS-AA 中苯环碳的 ε≈0.29 kJ/mol、σ≈3.55 Å，需与电荷分布配合。一个常见陷阱是：只拟合密度而忽略内聚能密度（CED），会导致链间作用过强或过弱，表现为玻璃化转变温度 T_g 偏移。验证时应同时检查 298 K 下的密度（PS 约 1.05 g/cm³）与 CED（约 300–400 J/cm³）。

静电项在含极性侧基的聚合物中不可忽略。电荷可通过 RESP 拟合 HF/6-31G* 静电势得到，并保持重复单元总电荷为零。若使用 UA 或 CG 模型，电荷常被合并或隐式处理，此时需重新标定介电屏蔽，否则会出现虚假的链聚集。

聚合物专用力场的选择与验证

选择力场时，先问三个问题：体系是否含氢键或强极性？目标性质是局部结构还是介观形貌？可接受的计算成本是多少？若研究聚乙烯熔体的链尺寸，TraPPE-UA 是性价比最高的选择；若研究聚电解质与水的相互作用，必须用全原子力场并显式水；若研究嵌段共聚物微相分离，CG Martini 可在微秒尺度给出有序结构。

验证方法应分层进行。第一层是单链验证：在稀溶液或真空下计算 R_g、持续长度 l_p、特征比 C∞，与实验或文献值对比。第二层是熔体验证：检查密度、CED、T_g、链间径向分布函数 g(r)。第三层是动力学验证：比较扩散系数、Rouse 松弛时间或缠结模量。若某一层偏差超过 10%–15%，应回溯参数化环节，而不是盲目延长模拟时间。

最后需要强调：力场没有“绝对正确”，只有“对目标性质足够好”。同一聚合物在不同力场下的 R_g 可能相差 10%–20%，这并非错误，而是分辨率与参数化目标的差异。工程上应明确模拟目标，选择匹配的力场，并用实验或高精度计算做闭环验证，而不是追求单一力场的普适性。