正文

薄膜里的高分子链打结概率不是随厚度单调变化的，而是在狭缝厚度约等于链的本体尺寸时达到峰值——这篇2026年10月19日在线发表于《Macromolecules》的研究，用大规模分子动力学模拟把“受限熔体里的打结与缠结”这一此前只有零星结果的问题系统地算清楚了。说明：近7天内未检索到细节完整的MPC/MD新论文，本栏放宽至10月15–19日发表的论文，日期均如实标注。

背景：为什么要在薄膜里数“结”

高分子熔体可以看成一大锅互相缠绕的长链，链与链之间的缠结决定了材料是像液体一样流动还是像橡胶一样回弹；链自身打结（像耳机线自己绕出的结）则会改变局部运动与力学响应。把熔体压进纳米薄膜（比如涂层、光刻胶、柔性电子里的聚合物层）时，受限会改变链的形状与缠结，于是打结概率也跟着变。问题是：变薄一定让结变少吗？薄膜里的局部性质能不能从厚膜推出来？刚度（链有多硬）会不会改变答案？此前只有单链的格子模拟给过线索，真实熔体里的数据一直缺。

他们做了什么：Kremer-Grest链 + LAMMPS，在不同厚度的狭缝里算结

作者用的是高分子模拟的标准模型——Kremer-Grest珠簧模型：把链看成一串珠子，相邻珠子用谐振弹簧连着（键长约0.967σ，弹簧刚度400 kBT/σ²），不相邻的珠子之间只有纯排斥作用（WCA势，即只推开、不吸引，相当于只算体积排除），再加Kratky-Porod弯曲势调刚度（参数B=0为柔性链，B=2为半柔性链）。所谓σ、τLJ是模拟里的约化长度与时间单位，相当于把真实单位换算成模型单位。体系是两面平壁夹着的熔体，链长N=256、512、1024，链数48–768条；先用Langevin动力学（一种带随机碰撞与摩擦的恒温方法）平衡至少2个链松弛时间，再用Nosé-Hoover恒温方法生产至少10个链松弛时间，N=1024时典型生产时长达5×10⁷ τLJ，步长0.01 τLJ；另用单链reptation蒙特卡洛（一种让链像蛇一样蠕动换构象的抽样方法）作对照。

关键结果：打结概率在h≈Rg处见顶，缠结却单调下降

第一，打结概率随狭缝厚度h呈非单调变化：在h约等于本体回转半径Rg,bulk（链在自由状态下大概占多大的尺寸）附近达到最大值，h再小就趋向零（二维极限下链几乎铺平、打不出结），h再大则回归本体值。第二，缠结行为不同：缠结长度随受限单调增大、构象愈发扁平——也就是说“结多了”不等于“缠结多了”，两者要分开看。第三，分层解析：厚狭缝里靠近壁面的局部性质（回转半径调制、壁面附近打结增强）积分起来，可以重建薄狭缝的整体厚度依赖，说明薄膜行为可以从局部信息拼出来。第四，刚度效应：半柔性链的打结最大值略向更小的h移动，刚度只是微调峰位、不改变非单调趋势。

为什么重要

它给薄膜与受限体系提供了一个可操作的标尺：判断受限强弱，先看h与Rg的比值——h≈Rg附近是打结最强的区域，而不是越薄越弱。这对涂层、薄膜器件里链拓扑导致的老化、断裂与输运行为有直接参考价值，也澄清了一个常见误解：把“打结”与“缠结”混为一谈会得出错误趋势。

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

关联直接但需区分：本文的链模型、平衡与生产流程、松弛时间标度与缠结长度分析都可复用，若课题涉及受限、薄膜或界面链，h≈Rg是判断受限效应强弱的基准。但本文是熔体MD、没有显式溶剂，流体力学耦合部分不能直接照搬，需要另参考MPCD溶剂的文献；把本文的拓扑统计方法嫁接到MPC+MD的剪切或受限模拟中，是自然的下一步。

来源

Macromolecules（ACS），2026-10-19，DOI：10.1021/acs.macromol.6c01338