正文

“绳子打了结反而更容易断”是日常经验，但分子层面到底为什么、断点会不会随结的复杂程度移动？一篇10月28日发表于《Macromolecules》的反应性分子动力学研究，用200个重复单元的聚乙烯单链、五种结（从三叶结3₁到七叶结7₁）、每种条件100条独立轨迹，给出了迄今最系统的答案：拉断分三个阶段，而断点位置随结复杂度系统地从结的出入口移向结的核心。

背景：要解决什么问题

单分子实验早就发现，打结的DNA在拉力下结会收紧、沿链扩散，而且结越复杂扩散越慢；打结会显著提高拉伸下的断裂率。但此前的模拟大多只做到五交叉以下的结，断点是否总在出入口、高复杂度下会不会出现质变，一直没有定论。更重要的是，断裂力与刚度随结复杂度的标度关系也没有系统数据。

他们做了什么

作者用ReaxFF反应力场做拉伸模拟。ReaxFF是一种允许化学键在模拟中生成与断裂的经验力场——与普通力场“键永远不断”不同，它能直接看到哪根C–C键先断，因此适合研究断裂；代价是绝对断裂力会被系统性高估（本文无结链的峰值断裂力约15 nN，约为量子化学计算的5–7 nN的两倍），所以作者强调只比较相对趋势。模拟流程是：先用OPLS-AA力场在300 K下平衡1纳秒并固定两端，再切换到ReaxFF、以0.25飞秒步长拉伸；每种结类型、每种拉速各跑100条轨迹，统计断点分布并与“均匀随机断裂”的零假设对比。

关键结果

第一，三阶段机制始终成立：熵驱动收紧（链先靠热运动把结收拢）、焓驱动收紧（继续拉伸时键角与键长被迫变形、能量上升）、键断裂。熵是构象混乱程度的度量，焓是键变形储存的能量——前一阶段靠“收拢省构象”，后一阶段靠“硬拉变形”。第二，断裂力随结复杂度呈拉伸指数衰减（结越复杂，断得越容易，但衰减速度比指数慢），而弹性模量随复杂度指数增长——拓扑受限同时让链“更易断、更硬”，这是两个相反方向的效应。第三，断点位置发生 crossover（交叉转变）：低复杂度结（3₁、4₁）主要断在结的出入口；高复杂度结（7₁）则主要断在密集堆积的核心，出入口偏好消失到均匀分布的水平。作者用C–C–C键角分析给出结构解释：随着复杂度上升，出入口与核心的键角畸变差异逐渐消失，几何上的“薄弱点”不再集中在出入口。

为什么重要

它把“结削弱绳子”从一句俗语变成了可调控的设计变量：想让材料在特定位置断、或想提高抗断能力，不能只看结的有无，还要看结的类型与复杂度。断点从出入口到核心的转变，也提示纳米孔转运、凝胶电泳等实验中，结的位置效应可能随复杂度而反转，不能用三叶结的结论外推到所有结。

与MPC＋MD课题的关联

方法论关联度高，体系关联度中等。你的课题是MPC＋MD模拟高分子链，本文的拉伸协议、系综统计（100条轨迹才敢谈分布）与“断点分布 vs 均匀零假设”的检验方式，都可以直接借鉴：例如在MPC流动中统计链的断裂或解缠位置时，同样需要先建立均匀分布的零假设，否则“某处更容易断”可能是链长分布本身造成的假象。但要注意三点局限：一是本文是单链真空拉伸，没有溶剂、没有MPC碰撞、没有多链缠结，结论不能直接外推到溶液剪切；二是ReaxFF的绝对力值有系统偏差，只能谈相对趋势；三是恒温器（Langevin，阻尼时间100飞秒）可能压制断裂瞬间的局部升温，绝对断裂力与断点分布的细节仍需未来用微正则模拟检验。

来源

Macromolecules，2026-10-28，DOI：10.1021/acs.macromol.6c01298