机械-静电-构象耦合下聚电解质链的非线性力-伸长响应与亚稳异质性

《Small Structures》：Nonlinear Force-Extension Response of a Polyelectrolyte Chain Undergoing Mechano–Electrostatic–Conformational Coupling and Metastable Heterogeneity

　　

拉伸柔性聚电解质会将机械形变、电荷调控与构象异质性耦合在一起。研究人员开发了一个单链热力学框架（single-chain thermodynamic framework），该框架综合了有限链伸展性（finite chain extensibility）、头-尾构象分配（head–tail conformational partitioning）、区域电离（regional ionization）、反离子熵（counterion entropy）、静电屏蔽（electrostatic screening）和介电效应（dielectric effects）。链被表示为致密且近似各向同性的头部与伸展且高度取向的尾部共存。自由能最小化预测了从致密球状（compact globule）经蝌蚪状共存态（tadpole-like coexistence state）到高度伸展构象的多阶段转变。计算得到的力-伸长曲线再现了分子模拟中报道的主要非线性趋势。因此，该模型为解释拉伸柔性聚电解质的机械与电荷调控耦合行为提供了一个物理透明的框架。

论文解读

刺激响应性高分子材料可在温度、pH、离子强度、化学配体或外加形变等刺激下发生相态、体积、力学性质与输运行为的可逆变化。单分子力谱（SMFS）等技术能够在纳米尺度和皮牛力分辨率下测量单个大分子的形变，但仅靠力学测量难以直接获得局部电荷与构象信息。中性柔性聚合物在良溶剂中的力-伸长行为可由朗之万或蠕虫链模型（WLC）描述，而聚电解质的拉伸行为因长程静电作用与反离子效应复杂得多。经典Manning反离子凝聚模型假定无限长刚性线电荷，不适用于柔性链；近年发展的自洽热力学方法虽纳入链柔性与电荷调控，但多数平均场模型仍假定空间均匀电荷分布，并常采用高斯链统计，因而不能描述强拉伸下有限链伸展性以及致密球体与伸展尾部共存等异质构象。因此，建立一个能同时考虑构象不对称、有限伸展性、静电屏蔽和局部电荷调控的单链热力学框架，具有重要意义。

研究人员构建了一个单链热力学模型，用于描述柔性聚电解质在外力作用下的机械-静电-构象耦合行为。模型将链划分为致密且近似各向同性的头部和高度取向的拉伸尾部，用取向参数χ=N_tail/N描述构象分配；头、尾区域分别赋予独立的电离度α_head和α_tail，以体现局部电荷异质性。头部静电自由能采用体积依赖的Debye–Hückel近似，尾部则采用一维屏蔽静电势描述。总自由能包含非高斯弹性熵、反离子结合组合熵、自由离子平动熵、静电相互作用能以及介电失配贡献，并通过数值求解?F/?χ=0、?F/?α_head=0和?F/?α_tail=0获得平衡构象，再由f=dF_eq/dR计算力-伸长曲线。该研究为理论建模，无实验样本队列；主要计算参数为N=1000、b=1.0、ξ=1.0、φ=1.0、c_s=0.01、l_B^head=3.0、l_B^tail=0.6、δ_head=1.5、δ_tail=1.0、Ω=5×10⁵。

研究结果包括三个部分。

3.1 宏观力学响应与构象相变规律。通过数值求解平衡方程并计算力-伸长关系、弹性模量、弹性自由能F_en和取向参数χ，发现力-伸长行为不同于中性柔性聚合物的连续变形。小伸长时链保持致密态并抵抗形变，机械功主要用于结构重排和破坏局部静电稳定；随后进入头-尾共存区，弹性模量下降，柔顺性增强；当大部分链段进入尾区后，模量再次上升。F_en随伸长先升后降再迅速上升，分别在R_e≈275 nm和R_e≈625 nm附近出现特征转变点，对应准球状到蝌蚪状、再到高度伸展态的转变。χ在低伸长时接近零，随后快速增大，说明链段由头部转移至取向尾部，以释放头部过度拉伸造成的熵惩罚。

3.2 区域电离度演化与电荷调控机制。通过计算全局与头/尾区域电离度以及静电能量贡献，发现尾部保持较高电离度，因为其伸展、暴露的几何构型降低局部电荷拥挤，有利于反离子释放；头部则在收缩过程中逐渐被中和，以降低致密区域保留净电荷的静电代价。尾部静电贡献随其增长而增加，头部静电贡献随其收缩与中和而减小，头-尾静电相互作用也因头区电荷减少而减弱。这表明构象演化与电荷调控相互耦合：结构不对称创造不同的静电环境，电荷再分布又稳定了非等效区域。

3.3 局部静电环境效应。通过改变头、尾区域的局部Bjerrum长度（l_B，两电荷静电作用能与热能相等时的特征距离）并计算力-伸长曲线、取向参数和头区电离/静电能，发现头区Bjerrum长度变化对早期构象演化影响较弱：在力平台结束前χ曲线几乎重合，表明头-尾共存结构的初期形成和早期生长受头区静电条件影响有限。模型同时考察了尾区Bjerrum长度对力-伸长响应和取向参数的影响，进一步说明局部静电环境可调控构象转变与电荷再分布。

综合引言中对现有理论局限的分析和数值结果，可以得出以下结论：柔性聚电解质在外部拉伸下可自发发生构象分配，形成紧凑头区与取向尾区；自由能最小化预测链经历从致密球状经蝌蚪状共存态到高度伸展态的多阶段转变；计算得到的力-伸长曲线再现了分子模拟中报告的主要非线性趋势。区域电离度呈现头部逐渐中和、尾部保持较高电离的异质性，且局部静电环境会调节构象转变与电荷再分布。该模型为解释拉伸柔性聚电解质的机械-电荷调控耦合行为提供了一个物理透明的理论框架，有助于连接单链局部电荷调控与宏观力学响应。