空间调制的驱动外场下活性聚合物的动力学行为*

高艺雯 王影 田文得 陈康

(苏州大学物理科学与技术学院,软凝聚态物理及交叉研究中心,苏州 215006)

活性聚合物由于其可形变的长链构象而表现出更丰富的动力学行为.本文通过郎之万动力学模拟研究了在周期性外场调制下活性聚合物的运动行为.我们发现柔性聚合物链有4 种跨区域运动模式.外场的周期性区域宽度在活性聚合物行为的调控方面起到了关键作用.比如随着区域宽度的改变活性聚合物的扩散系数可以发生高达两个数量级的改变,其余如链的统计空间密度分布、整体受力情况、特征尺寸和取向等都表现出明显的非单调变化.这些变化反映了不同区域宽度下,聚合物链的构象特征及运动模式的改变.我们的研究为设计和调控活性聚合物的动力学行为提供了新思路,有望为可形变链状或长条状活性物质在生物、医学和材料领域的潜在应用提供参考.

1 引言

活性物质是由可以自主运动的个体组成的非平衡体系.活性物质个体种类繁多,形状各异,尺寸可跨越微观和宏观.细胞内,吸附在微管和微丝上的马达蛋白可以通过ATP 的水解把化学能转化为机械运动,从而实现一系列生物功能,例如细胞内物质输运和协助有丝分裂.正因为有了这些微观机器的驱动,各种生物,小到细菌大到鲸鱼,才表现出生生不息的活力.体外实验表明,在嫁接有马达蛋白的衬底上,微管微丝如同自驱动的聚合物链(活性聚合物链)一样,可以自发组装成团簇、旋涡、条带[1]以及向列相涡旋阵列等动态结构[2].另外,许多游动的微观生物体[3]呈现出细长的链状结构,还有细菌[4]往往通过鞭毛或纤毛等链状结构来推动自身.这些都是活性聚合物的例子,它们表现出令人着迷的集体行为.除了生物活性聚合物,实验上成功合成了由两面神胶体粒子组装成的活性聚合物链[5],为活性聚合物的研究提供了更可控的平台.

活性聚合物的动力学与自组装行为涉及聚合物物理与非平衡物理的交叉,因而吸引了一系列理论与模拟的研究[6,7].例如,聚合物链在活性粒子浴中的运动[8]、聚合物链在活性短链浴中的运动[9]、活性链在多孔介质中的运动[10],还有活性链在通道中的输运[11]等.对单根活性蠕虫状链的研究表明[12],在不同驱动力强度和链刚性下,活性聚合物在二维可呈现出普通的聚合物线团态、弱螺旋态和螺旋态.若聚合物的前端存在负载,则活性聚合物会形成圆、摆动或螺旋结构[13].

对于很多活性物质体系,个体的活性可以通过外部条件调控,如光照条件[14].目前对活性物质的理论模拟研究大部分都假设个体的活性不随空间位置变化.少数针对活性粒子气的理论与模拟研究表明活性粒子的空间浓度分布与其局域活性驱动速率成反比[15];粒子速率的概率分布存在非高斯性,这种非高斯特征随着活性驱动速率在空间变化的加强而加强;而且在短时区域,活性粒子的速度-速度关联以及均方位移受活性驱动速率在空间变化的影响很大[16].活性聚合物与活性胶体粒子不同,它的尺寸与构象是不可忽视的因素.聚合物的长链结构,使它的不同部分可以同时处于不同空间活性区域,从而在链上产生不均匀的驱动.一个有趣的问题是活性聚合物在驱动强度受到空间位置调制的情况下,它的动力学行为会表现出怎样的特征.这方面的研究将为外场调控活性聚合物的群体行为提供基础,有望为它们在生物、医学和工业领域的潜在应用提供参考.

本文主要研究了活性聚合物在有外场调制情况下的运动,考虑了最简单的情形,即在周期性外场下,空间形成活性与非活性交替的区域.发现柔性链在大区域宽度下(d ≫2Rg0,Rg0为无规线团态的回旋半径),链长时间受限于单一区域中,并进行缓慢的布朗平移扩散;在中等区域宽度下(2RL

2 体系与模型

我们使用布朗动力学模拟研究了空间周期调制场下的单条活性聚合物链.该链由N=50 个直径为σ,质量为m的活性粒子线性排列而成.考虑最简单的一维周期性调制场,即在x方向形成周期性交替的活性区域(active 区域)和非活性区域(passive 区域).x方向周期数为n,每个区域宽度为d,因此体系x方向尺寸Lx=2×n×d.体 系y方向是均匀的,尺寸Ly=100σ,如图1 所示.周期之间以黑色实线划分,每个周期内点线左侧为active 区,右侧为passive 区.区域宽度d是本文主要调控参数之一.

图1 空间调制场与活性聚合物示意图(d=32σ ,Lx=192σ)Fig.1.Schematic diagram of spatial modulation field and an active polymer (d=32σ ,Lx=192σ).

活性聚合物链上每个粒子的运动遵循郎之万动力学方程:

其中,ri是第i个粒子的坐标,i为粒子的编号,m为粒子的质量,t为粒子运动的时间,−∇i是势函数的负梯度,表示粒子受到的力.方程(1)右边第一项中的U({ri}) 是相互作用能,由三部分组成:

Ubond是键连接的弹簧谐振势:

其中r0=1σ是键的平衡长度;Ka=1000kBT/σ2为弹簧系数,kB为玻尔兹曼常数,T为温度.Ubend是链的弯曲势:

其中,K为刚性系数,θi是由粒子i −1 ,i和i+1组成的相邻两个键之间的夹角.UWCA是非相邻粒子间的体积排斥相互作用势[17]:

其中,ε为相互作用强度,rij=|ri −rj|为粒子i和j之间的距离,作用势在rc=21/6σ处截断.

其中,单位矢量=(ri+1−ri−1)/|ri+1−ri−1| 表示自驱动力方向沿活性链切线方向,如图2 所示,Fa为自驱动力大小.

图2 活性力示意图Fig.2.Schematic diagram of active force.

我们利用LAMMPS 软件进行模拟[18],模拟中使用约化单位.设置m=1,σ=1和ε=1,因而时间单位τ=固定体系的约化温度kBT=1,摩擦系数γ=10 (摩擦系数足够大以保证聚合物链的运动满足过阻尼条件).当聚合物链上的粒子位于活性区域时,它会受到驱动力Fa=10 .如果整条K=0 的柔性链处于活性区域,它会形成直径2RL≈8的稳定螺旋(图3(a)和图3(b)).而当粒子位于非活性区域时,设置Fa=0,即粒子做普通的布朗运动.如果整条K=0 的柔性链都处于非活性区域,它呈无规线团状,模拟得其尺寸2Rg0≈14(Rg0表示均方根回旋半径;符合二维自回避链公式Rg0≈0.35σN3/4≈6.6[19],N为活性链的聚合度).在正式计算前,让聚合物链在无自驱动力条件下进行充分的热平衡.在之后的模拟中,选取时间步长为 0.005τ,每个轨迹模拟总时长为 50000τ.对所有的参数条件都做了10 次独立的模拟.

3 结果与讨论

3.1 柔性链(K=0)的跨区域运动

以K=0 的柔性聚合物链为例来探讨柔性聚合物链的运动.也模拟了其他小K值的情况,如K=3或 5 等,结果都定性一致.柔性聚合物链的特点是在无外场(非活性区)情况下(Fa=0)做无规布朗运动,整条链呈无规线团构象,而整条链在有外场情况下(活性区,Fa=10),会由于弯曲及自身的碰撞,形成稳定的螺旋结构,但整个螺旋的迁移几乎不受自主驱动的影响,表现出无规热运动的特征.在本文设置的周期性外场中,在大区域宽度时(d ≫2Rg0),聚合物链有很大的概率位于活性区,做螺旋运动,当然也有很大概率位于非活性区做无规布朗运动(图3(a)).当区域宽度逐渐减小,处于中等区域宽度时(2RL

图3 柔性链在不同区域宽度下的构象快照 (黄色粒子代表聚合物链头部,蓝色粒子代表聚合物链尾端) (a) d=32,活性区域的稳定螺旋态以及非活性区域的无规线团状态;(b)—(f) d=8,活性区的螺旋态、非活性区的拉伸的无规线团以及各种跨区域瞬时构象;(g),(h) d=3,跨区域的拉伸和卷曲塌缩构象Fig.3.Snapshots of the conformation of the flexible chain at different region widths (Yellow particles represent the head of polymer chain,while blue particles represent the tail of polymer chain): (a) d=32,stable helical state in the active region and random coil state in the inactive region;(b)–(f) d=8,helical state in the active region,stretched random coil state and various cross-regional instantaneous conformations in the inactive region;(g),(h) d=3,crossregional stretched and curled collapsed conformations.

如前所述,在中等及窄区域宽度(d<2Rg0)时,聚合物链的跨区域运动频繁,构象变化强烈.我们发现跨区域运动存在4 种典型的过程.第一种是聚合物链的头部从非活性区进入活性区.由于头部链段获得驱动力形成定向加速,从而试图拖拽后面链段也进入活性区,实质上这时处于非活性区的链段成为了整条链自驱运动的“负载”.由于柔性,进入活性区的头部会形成卷曲,往螺旋态发展.当卷曲很早发生且螺旋的形成速度大于链被拖入活性区的速度,就会看到形成的不完整螺旋一边长大一边向区域边界靠近,最终返回非活性区(图4(a)—(c)).卷回非活性区的聚合物螺旋失去驱动力,其紧致的螺旋结构无法维持,向无规线团转变.第二种也是聚合物链的头部从非活性区进入活性区.但头部链段没有很快形成卷曲螺旋,而是有相当部分链段已经进入活性区后螺旋才形成(在活性区狭窄时会出现跨多区域的情况,图4(d)).此情况下,非活性区域链段对链整体自驱运动的阻碍比较弱,所以链段被拖入活性区的速度不慢于螺旋的生长,最终我们看到活性区的螺旋在形成过程中没有向边界靠近,可以在活性区形成完整的螺旋态(图4(d)—(f)).第三种同样是聚合物链的头部从非活性区进入活性区.进入活性区后,链始终没有发生卷曲(图4(g)),最终头部穿越活性区进入下一个非活性区(图4(h)).此时,由于失去自驱动力,头部链段的定向运动消失,成为了后续链段进入非活性区的阻碍.在活性区域链段的不断进入及推挤作用下,聚合物链在非活性区靠近边界处形成主要沿y方向延展的堆积.整体上看,在此跨越过程中,聚合物链两端处于非活性区呈无规的线团或堆积状态,而中间由拉伸的处于活性区的链段连接.第四种是聚合物链的中间部分从非活性区进入活性区.由于驱动力是沿链切线方向,所以链上粒子的自主运动产生“跟随”效应.进入活性区的链段,其前端粒子会返回非活性区,而后端会拖拽其后的链段跟随进入活性区,从而使得其后链段产生明显的(纵向即y方向)拉伸效应(图4(j)—(l)).第四种跨区域运动出现的最频繁.

图4 柔性链跨区域运动的快照(d=8) (a)—(c) 第一种跨区域运动模式,即头部从非活性区进入活性区,很快形成螺旋并最终返回非活性区;(d)—(f) 第二种跨区域运动模式,即头部从非活性区进入活性区,在大部分链段进入活性区后形成螺旋并最终在活性区形成完整螺旋;(g)—(i) 第三种跨区域运动模式,即头部从非活性区进入活性区,并跨越整个活性区进入下一个非活性区;(j)—(l) 第四种跨区域运动模式,即中间部分从非活性区进入活性区造成其后面的链段产生明显的(沿纵向即y 方向的)拉伸Fig.4.Snapshots of cross region motion of the flexible chain(d=8).(a)–(c) The first mode of cross region motion: the head enters active region from the passive region,rapidly forms spiral and eventually returns to the passive region.(d)–(f) The second mode of cross region motion: the head enters active region from the passive region,forms spiral after most of the chain segments have entered active region and eventually forms complete spiral in the active region.(g)–(i) The third cross regional motion mode: the head enters active region from the passive region and crosses the entire active region to the next passive region.(j)–(l) The fourth mode of crosss region motion: the middle segments enter active region from the passive zone causing significant stretch (longitudinal,i.e.y -direction) of the chain segments behind it.

我们对聚合物链的整体运动进行了统计,计算了聚合物链质心的均方位移(mean square displacement,MSD):

其中,rcm(t′)是t′时刻聚合物链质心的位置,〈·〉代表时间和10 次独立轨迹的平均.图5 给出了不同区域宽度条件下聚合物质心的MSD,看到区域宽度对聚合物链的运动影响巨大.例如在长时间,窄区域宽度(d=3)和大区域宽度(d=50)之间MSD 的值差了两个数量级,曲线的特征彼此差别也很大.具体来看,大区域宽度(d ≫2Rg0)时,如d=32和d=50,由于聚合物链跨区域运动很少,聚合物链的运动主要是在非活性区做无规布朗运动或在活性区做螺旋运动,这两种情况下,质心的运动以普通热扩散机制为主,所以两条曲线整体上满足普通的扩散规律 ∆(t)∝tα且α=1 .由模拟曲线可得到长时的扩散系数D=∆(t)/(4t)≈0.0018,与Rouse 模型下预测的聚合物链整体扩散系数D=kBT/(Nγ)=0.002符合[20].另外,图5 在短时部分的两个曲线存在振荡,其原因是MSD 中有很大一部分贡献来自于在活性区做稳定螺旋运动的聚合物链.而螺旋的中心和聚合物链的质心并不重叠,所以除了平移热扩散,质心还环绕中心运动,这就导致了MSD 曲线的振荡.

图5 不同区域宽度下柔性链的质心均方位移 (不同颜色曲线对应不同区域宽度,断线表示不同的扩散系数)Fig.5.Mean square displacement (center-of-mass) of the flexible chain under different regional widths (Different color curves correspond to different regional widths,and the broken line indicates different diffusion coefficients).

在中等区域宽度(2RL

窄区域宽度(d<2RL≈8)时,发现MSD 曲线出现了明显的从早期超扩散到后期普通扩散的转变.这说明跨区域运动已经超越聚合物链在活性与非活性区的普通扩散运动而居于主导.频繁的跨区域运动特别是第三种跨区域运动模式,链在活性区受驱动力拉伸定向横跨整个活性区域,极大促进了链的迁移.由模拟曲线可以得到d=3 情况下,长时的扩散系数D=∆(t)/(4t)≈0.12 .如前所述,对比d=50,扩散系数相差了两个数量级.

MSD 给出了聚合物链运动的统计信息,而更直观的运动特点蕴含在所有轨迹当中.图6 展示了聚合物轨迹的一些例子,给出的是聚合物质心在x和y方向的演化.需要指出的是同样参数条件下,各独立轨迹差别很大,但表现出的特征比较类似.比如在窄区域宽度d=3,4时,x方向质心的运动轨迹表现出平台和跳跃频繁交替.其中平台对应于聚合物链被暂时局域在某非活性区域,而跳跃对应于跨越整个活性区域的运动,也就是图4(g)—(i)所指第三种跨区域运动模式.d=5 时这种跨越整个活性区域的运动已经减少很多.在给出的中等区域宽度d=10 的例子中已经没有看到这种跨越.当然它是可以发生的,只是出现概率很小.在d=20的例子中发生了一次这种跨越.对比之下,y方向质心的运动轨迹表现出正常的涨落行为,没有明显的平台和跳跃现象.

图6 柔性链在不同区域宽度的周期性外场中,质心沿(a) x 方向和(b) y 方向随时间运动的位置变化.在d=3,4,5 时,x 方向的位置变化出现平台与跳跃,平台对应链被限制在非活性区域,跳跃则对应链跨越活性区域;y 方向的位置变化出现正常的涨落.在d=10,20 时,链的运动幅度相对较小Fig.6.In the periodic external field of the flexible chain with different regional widths,the position change of the center of the mass moving along (a) x direction and (b) y direction with time.Different color curves represent different regional widths d respectively.At d=3,4,5,the position changes in the x direction appear platform and jump.The chain corresponding to platform is restricted to the passive region,while the chain corresponding to jump crosses the active region.The change of position in the y direction shows normal fluctuation.At d=10,20,the movement of the chain is relatively small.

聚合物链在不同区域构象与运动形式的差异以及各向异性的跨区域运动,这些因素导致聚合物链在空间的统计分布不均匀.为此,统计计算了不同条件下聚合物链在活性区域与非活性区域出现的概率情况.定义聚合物链在活性区出现的概率为φa=〈Na(t)/N〉,其中Na(t)是t时刻处于活性区域的聚合物链粒子数目.聚合物链在非活性区出现的概率则为φp=1−φa.图7(a)给出了不同区域宽度下聚合物链在两种区域中的概率分布.在大区域宽度d=50 情况下,聚合物链几乎有同等概率(φp≈φa≈0.5)位于活性区与非活性区,分别呈紧致的螺旋态构象和无规线团状.但在中等区域宽度情况(8≤d≤20) 下,链更多停留在非活性区域,ϕa几乎接近于0.这是因为此时区域宽度与链的构象尺寸相当或更小,螺旋态在活性区域无法长时间稳定存在,会频繁接触区域边界进入非活性区.而链处于非活性区时,同样会频繁跨越边界,但从图4 的跨区域运动模式看到,在跨区域运动过程中大部分结果(第一、三、四种情况)都是链返回或进入另一片非活性区域.第三种情况在活性区形成的螺旋也会很快接触边界重新进入非活性区.所以聚合物链主要逗留在非活性区,而活性区聚合物链的概率主要来自于跨区域运动过程中部分链在活性区的逗留以及偶尔形成的完整螺旋在活性区的短暂逗留.窄区域宽度(d≤8)时,区域尺寸远小于链的构象尺寸,跨区域运动几乎不间断,所以φa升高.甚至在d=3 时,聚合物链基本上始终以跨越多个区域的形式存在(图4(g),(h)),而φa也接近0.4.

从前面的叙述,我们看到柔性聚合物链进入活性区域后倾向于弯曲形成螺旋,而活性区的聚合物螺旋进入非活性区后会一定程度上保留螺旋式的弯曲.为了定量化这种链沿一个方向的持续弯曲程度(螺旋程度),定义平均螺旋角θ为

其中,由叉乘 (ri −ri−1)×(ri+1−ri) 的方向决定符号函数 sgn[θi(t)] 的正负(正负方向的选取不影响结果).(8)式是链整体的平均螺旋角.同时分别计算了在活性区和非活性区链段的平均螺旋角,

当d=50,大区域宽度下,聚合物链在活性区形成稳定的多圈螺旋所以螺旋角θa值较高接近0.17π(图7(b)).非活性区的链段则处于无规线团状,所以螺旋角θp非常小接近 0.01π .由于聚合物链在两种区域出现概率相当,所以总的平均螺旋角θ介于两个值之间.随着区域宽度变窄,在中等区域宽度情况下(8≤d≤20),活性区的链段螺旋角θa保持较大的值,甚至比d=50 时还略有上升,超过0.2π.这是因为在此宽度下,活性区链段常形成不完整的螺旋,即只有内圈部分(图4(a)和图4(b)).而内圈部分对应更高的螺旋角.但此时,聚合物主要分布于非活性区,所以总的平均螺旋角θ由θp的贡献为主.非活性区螺旋角θp相较d=50 时也有明显上升,达到约 0.02π,这源于跨区域运动过程中返回非活性区的链段部分保留了在活性区时的螺旋性(图4(a)—(c)和图4(j)—(l)).当区域宽度很窄时(d<8),聚合物链整体跨越多个区域,活性区链段的螺旋性因而明显减弱(图3(f)—(h)),但带动了非活性区螺旋性的进一步上升.

图7 柔性链的区域概率分布 φ (a)和平均螺旋角 θ (b)与区域宽度d 的变化关系Fig.7.Regional probability distribution φ of the flexible chain (a) and the average spiral angle θ (b) in relation to the variation of the region width d.

下面分析聚合物链的整体受力情况.关注作用在链上的自驱动力,为此统计了作用在整条链上的自驱动力沿x方向和y方向的平均分力以及合力:

另外,分别计算了聚合物链上粒子与聚合物链质心的短时平均运动速度,定义为

类似地,可以定义速度沿x方向和y方向分量(vx和vcx,vy和vcy)的平均大小.图8(b)是得到的速度变化的曲线,看到粒子运动的平均速度要远大于聚合物链质心运动的平均速度,说明链相当一部分运动是原地徘徊.特别是在大d情况下,虽然螺旋中粒子在自驱动下运动很快,但链作为整体平移很慢.在d<8时,粒子和整条链的y方向运动被显著加强,和图8(a)对应,原因也一致.

图8 (a) 柔性链的活性力 Fa与区域宽度 d 的非单调变化关系;(b) 柔性链的粒子和质心的短时平均速度(v 和 vc)与区域宽度 d 的变化关系Fig.8.(a) Nonmonotonic variation relationship between the active force Faand the region width d of the flexible chain;(b) the short-term average velocity of particles and center of mass of the flexible chain (vand vc) in relation to the variation of the region width d.

在周期外场下,聚合物链的构象变化剧烈,但总体上区域宽度狭窄时(d<2Rg0),x方向的受限会导致聚合物链沿y方向拉伸.为此,我们计算了链的回转半径Rg和末端距Re,并分别给出了它们x方向和y方向贡献.两个特征尺寸的定义式为

从图9 可看到,在区域宽度很大(d≥20)时,Rg和Re都有一个下降,这是活性区聚合物链形成稳定紧致螺旋导致的.d=50时,Rgx(Rex)和Rgy(Rey)相差不大,反映了聚合物链基本不受x方向边界的约束.些许差别来自于链跨越边界时沿x方向的拉伸行为.在区域宽度中等或很窄(d<20)时,链会显著地沿y方向拉伸,即Rgy(Rey)≫Rgx(Rex).有趣的是,在d=6左右发现了Rg(Re)和Rgy(Rey)的峰值.d>6时,区域宽度d的减小促进链沿y方向的拉伸.而d<6 时,由于链常处于跨越多个区域的状态(图3(g),(h)),区域宽度d的减小,反而使链感受到更均匀的环境,受限效应减弱,链沿y方向的拉伸也减弱.

链各向异性的构象还可以通过键的取向信息来获得.这里定义了聚合物链的沿y方向的向列相序参数S:

类似(9)式和(10)式,可分别定义活性区与非活性区链段沿y方向的向列相序参数:

从图10 可看到,聚合物链在活性区更倾向于沿y方向取向.随着区域宽度减小,受限程度的加强,无论活性区还是非活性区的链段沿y方向取向都得到加强.而峰值同样出现在d=5 左右,和图9 的Rg,Re曲线类似,其背后的物理机制也一致.

图9 柔性链的回转半径 Rg(a)和均方末端距 Re (b)与区域宽度d 的变化关系,在d=6 附近 Rg,Re 出现峰值,两者走势相似Fig.9.Radius of gyration of the flexible chain Rg (a) and the mean square end-to-end distance Re (b) in relation to the variation of the region width d.Peaks of RgandRe occur around d=6,both are similar.

图10 柔性链的向列序参数 S 与区域宽度 d 的非单调性变化关系,在 d=5左右存在峰值,Sa 相对变化较大Fig.10.Nonmonotonic relationship of the nematic order parameter S and the region width d of the flexible chain.There is a peak around d=5,and the change of Sa is large relatively.

3.2 半刚性链(K=10)的跨区域运动

3.1 节讨论了柔性链在周期性外场中的行为.链的刚柔性对此有很重要的影响,故这里以K=10为例来探讨半刚性链在跨区域运动过程中的表现.还模拟了其他大K值的情况,如K=20,50,100结果定性一致.它们共同的特点是在活性区时不能形成卷曲螺旋态,但链可以呈现明显弯曲.此外,还模拟了K=1000 的情况,此时链刚性很强,如细棒几乎不弯曲,其运动方式简单,这里省略不讨论.研究发现半刚性链的一个比较突出的行为特征是当它从活性区进入非活性区时会呈现屈曲拱起的构象.原因类似图4(g)—(i)的第三种跨越过程,即其头部一进入非活性区就失去了自驱力而慢下来,而后续链段在自驱力下往前推挤,但由于链有一定刚性不会发生不规则堆积,而是沿一个方向拱起,如图11 所示.

图11 半刚性链跨区域运动的快照(d=8),头部从非活性区进入活性区失去自驱力,链段往前推挤形成屈曲拱起Fig.11.Snapshot of the across region motion of the semirigid chain (d=8).The head enters the active region from passive and loses self-propelled force,and the chain segment is pushed forward to form buckling.

半刚性链虽然特征尺寸很大,但由于上述的拱起运动模式,使其在区域宽度不是很狭窄时很容易在非活性区形成沿y方向的伸展构象,并长时间逗留在非活性区做无规布朗运动,如图12(a)所示.随着区域宽度的减小(d<6),链形成的拱起很容易跨越边界导致链持续不稳定的跨区域运动,如图12(b)所示.这种持续不稳定的运动会导致半刚性链形成屈曲折叠构象(图12(c)).当区域宽度进一步减小(d= 3),半刚性链会始终跨越多个区域,产生多处拱起并在空间呈蛇行运动(图12(d)).

图12 半刚性链在不同区域宽度下的构象快照 (a) d=8,非活性区的沿y 方向伸展状态;(b),(c) d=6,跨区域的屈曲和屈曲折叠构象;(d) d=3,跨越多个区域的蛇行构象Fig.12.Snapshots of the conformation of semi-rigid chains at different domain widths: (a) d=8,stretching along the y-direction in the inactive region;(b),(c) d=6,buckling and buckling-folded conformation across regions;(d) d=3,serpentine conformation spanning multiple regions.

参照柔性链部分的分析,我们对半刚性链的区域密度分布、活性力与短时平均速度随着区域宽度的变化进行了统计(图13),得到了单调变化的曲线.因为在区域宽度较大(d>8)时,除了偶尔的跨区域运动,半刚性链基本上都逗留在非活性区,所以活性区的统计密度接近 0,链所受驱动力趋于0,短时平均速度也很小.当区域宽度变窄d≤6 时,半刚性链的跨区域运动开始频繁,活性区的统计密度、链所受驱动力、短时平均速度都有所上升.在d=3时,半刚性链在空间几乎是连续光滑的运动,活性区与非活性区链粒子密度几乎相当.特别值得注意的是在此条件下,链的平均横向受力超越纵向受力,即,表明链更趋于横向取向跨越活性区(图13(b)).这点可以和下面计算的向列相序参数互相印证,如图14(c)所示,这个序参数在d=3变为负值.

图13 (a)半刚性链的区域概率分布 ϕ 与区域宽度 d 的单调变化关系;(b)半刚性链的活性力 Fa与区域宽度 d 的变化关系;(c)半刚性链的粒子和质心的短时平均速度(v 和 vc)与区域宽度 d 的变化关系Fig.13.(a) Monotonous variation relationship with the regional probability distribution ϕ and the region width d of the semi-rigid chain;(b) relationship between the active force Fa and the variation of the region width d of the semi-rigid chain; (c) the short-term average velocity of particles and center of mass of the semi-rigid chain (v and vc) in relation to the region width d.

类似地,统计计算了半刚性链的特征尺寸Rg,Re以及向列相序参数,如图14 所示.发现和前面所述链跨越区域运动频繁相对应,d=6 成为转折点.当d<6,半刚性链的x方向特征尺寸急剧增加,而链整体沿y方向向列相序急剧减小,在d=3时甚至变为负值.

图14 半刚性链的回转半径 Rg(a)和均方末端距 Re (b)随区域宽度 d 的变化关系,在 d=6 时链的整体尺寸出现最低点;(c)半刚性链的向列序参数随区域宽度 d 的变化关系,在 d=6时,Sa 出现 峰值Fig.14.Radius of gyration of the semi-rigid chain Rg (a)and the mean square end-to-end distance Re (b) in relation to the variation of the region width d.The overall size of the chain at d=6 is the lowest point.(c) The relationship between the nematic order parameter of the semi-rigid chain with the width of the region d.At d=6,theSa shows peak.

4 结论

关于活性聚合物的研究,近几年才刚刚起步,有很多值得探索的问题.本文着眼于活性聚合物单链运动行为的调控,设想了简单的周期性外场对聚合物的活性进行空间调制,发现柔性聚合物链有四种跨区域运动模式.与此相对应,随着外场周期性区域宽度的变化,活性聚合物的迁移运动MSD、统计空间密度分布、整体受力情况、特征尺寸和取向等表现出三种不同的行为.柔性链在活性区域的螺旋态尺寸 2RL和在非活性区域的无规线团尺寸2Rg0共同标定了区域的宽窄.半刚性链的运动模式相对简单.在一般区域宽度(≳2RL)下倾向于在非活性区逗留并沿y方向伸展.只当区域宽度很窄(≲2RL)时,半刚性链的运动会发生向全空间连续运动的转变.相应地,链的统计尺寸、极性和空间密度分布等都发生显著的转变.

本文只是对外场调控活性聚合物运动的初步探索,其实体系可以改变很多值得探索的物理参数.例如,本文中链的长度没有改变;活性区域和非活性区域的宽度可以不一致;外场的空间变化可以不是方波型而是光滑缓变的;外场的调制可以不是单一方向(x方向)的变化,而是采用复杂的空间图案.当然更有趣且更值得研究的是多体活性聚合物在外场调制下的自组装,这也是我们正在开展的研究课题.最后,关于活性聚合物行为调控方面的研究不仅在基础科学方面是有意义的课题,同时它有望为链状活性物质的设计及其在生物、医学、材料等领域的潜在应用提供参考.