正文

主链与支链的相对刚度，决定了非对称的机械互锁环与化学键合环在剪切流里是翻跟头还是“履带转动”——这篇 2026 年 10 月 15 日发表于《RSC Advances》的研究，用多粒子碰撞动力学（MPCD）耦合分子动力学（MD）的混合模拟，第一次系统地比较了非对称与对称两种情形。对做 MPC+MD 模拟的人来说，它最直接的价值是：同一套模拟框架里，环的大小与连接方式这两个参数，就能调出完全不同的运动模式。

背景：为什么要比较两种“双环”

环形高分子就是把线性链的首尾接起来得到的闭环，在自然界对应质粒 DNA、动质体 DNA 这类结构。两个环可以有两种连接方式：一种是像钥匙环一样机械互锁、没有化学键把它们拴死的“索烃”（catenane，文中记为 PC），另一种是两个环之间有一根不能弯折、不能扭转的化学键直接连着的“键合环”（bonded rings，记为 BR）。此前对柔性双环的研究发现两者都会翻滚，但机械互锁的是断断续续的“滑移翻滚”，化学键合的是连续的“梯度翻滚”；而半柔性、非对称（两个环大小不一样）的情况在剪切流里会怎样，此前并不清楚。

他们做了什么：MPCD 做溶剂，MD 做链

作者用粗粒化的“珠簧模型”（bead-spring，也就是 Kremer-Grest 模型）来搭链：把每个环看成一串由弹簧连起来的小珠子，忽略具体的化学组成与氢键等细节，专注于环大小与刚度带来的宏观效应。溶剂不是用显式的水分子一个个算，而是用 MPCD 处理：溶剂粒子先自由运动一步，再按空间格子分组、把同一格子内粒子的速度绕质心旋转一次——这一步以很低的计算代价给出了真实的流体力学相互作用（即一个珠子运动会拖动周围流体、进而影响别的珠子），这正是剪切流模拟里不能省略的部分；而链本身的运动则用常规 MD 积分。所谓“半柔性”，是指链有一定的弯曲刚度，介于完全柔软的面条与完全僵硬的棍子之间。两个环大小不同时记为“非对称”，大小相同时记为“对称”，以此对照。

关键结果：连接方式决定运动模式

第一，形状：通过分析环的“扁长程度”（prolateness，越接近长条数值越大）、“非球形度”（asphericity，偏离球形越远数值越大）与三维表面积，作者发现非对称的 PC 与 BR 在剪切下都呈现相同的盘状结构；而对称情形下两者不同，对称的 PC 要先经过“剪切变稀”（shear thinning，即剪切越强、表观黏度越低、链被拉得越扁）之后才形成盘状。第二，运动模式：非对称的 PC 做翻滚运动（tumbling，像车轮一样整体翻转），非对称的 BR 则做“履带转动”（tank-treading，像坦克履带那样整体位置不动、表面物质循环滚动）；对称的 PC 是“滑移翻滚”与翻滚交替，对称的 BR 是“履带转动”与滑移交替。第三，时间演化：通过追踪两个环质心之间的相对位置矢量与环面法向量的分量随时间的变化，作者把两类体系的运动归为“连续运动”与“间歇运动”两大类——前者平滑持续，后者时动时停。

为什么重要

这项工作的意义在于把“拓扑连接方式”从一个静态的结构标签，变成了可以预测动力学行为的变量：同样是两个环，机械互锁还是化学键合、大小是否对称，就决定了它在流动中是翻滚还是履带转动、是连续还是间歇。这对设计分子机器、拓扑软物质材料有直接启发——如果你希望某个双环器件在流动中稳定地“履带转动”而不翻滚，选键合环、调非对称可能是更稳的路线；反之，想要间歇性的翻滚响应，机械互锁的对称结构更合适。

与 MPC+MD 模拟高分子链课题的关联

关联非常直接：这篇用的正是 MPCD 溶剂 + MD 珠簧链的混合方案，与课题的技术路线同源，可以直接借鉴其“用 MPCD 保证流体力学相互作用、用珠簧模型控制刚度与拓扑”的做法。更具体地说，文中用扁长程度、非球形度、质心相对矢量与法向量这套分析量来区分翻滚与履带转动，这套“形状—运动”判据完全可以移植到自己体系的剪切模拟中，用来判断不同刚度组合下链的运动模式是否发生了转变。需要注意的局限是：该模型忽略了氢键与溶剂化等化学细节，结论是拓扑与刚度层面的普适规律，不能直接对应某一种具体化学组成的实验数值。

来源

RSC Advances，2026-10-15，DOI：10.1039/d6ra05054a