两端锚定单链聚电解质在交流电场中滞后效应.docVIP

两端锚定单链聚电解质在交流电场中滞后效应 　　【摘 要】 在均匀的交流电场作用下，应用布朗动力学模拟方法，研究了两端锚定的柔性和半柔性的单链聚电解质的动力学行为。研究结果表明，链段的位置随交流电场的变化表现出类似于磁滞回线的滞后效应，随着交流电的频率减小，这种滞后效应逐渐消失，成为一条单一的曲线。本项目针对两端锚定的聚电解质随交流电场变化的动力学行为的研究，文献中未见报道。 　　【关键词】 滞后效应 聚电解质 交流电场 　　【Abstract】 We investigate the dynamics of a single flexible and semiflexible polyelectrolyte （PE）， fixed at both ends under the action of spatially uniform alternating-current （ac） electric fields using Brownian dynamics simulations. The results show that the shapes of the orbits in the position-electric field plane are strongly influenced by the frequencies of ac electric field. These orbits transform gradually from hysteretic cycles to single curves with the decreasing frequencies. 　　【Key words】 Hysteresis effect Polyelectrolyte Ac electric field 　　1 介绍 　　近年来，在使用高频率电场操纵聚电解质分子逐步引起关注，特别是关于DNA方面的研究[1-4]。例如，将长度超过7μm的长肌动蛋白丝两端用两个金属电极吸引，终止交流电压后，观察两个电极之间的情况[5]。前人的工作表明，单链PE在交变电场中的构型转变存在滞后性[6]。在聚电解质的动力学和电学行为的研究中，一个单链PE的操纵是至关重要的。与研究链在静态下的行为相比较[7]，在外部周期力的扰动下的单链的动力学行为更值得研究[8]。最近一项研究表明[9]，通过动态蒙特卡洛模拟，在震荡力下的单链聚合物展现出了周期行和非周期性的动力学行为。在本文中，我们将研究一条两端锚定的柔性或半柔性的PE链，链的两端将处于均匀的交流电场内。通过对链的布朗动力学模拟，我们将研究交流电场的频率对链的位形的影响。 　　2 模型和方法 　　我们的系统包括一条单一的线性聚电解质链，处于一个各项均匀的电场中。链的模型为竹簧链模型[10]。该链由N个粒子（单体）构成，由N-1个距离为的键连接。在我们的模拟中，我们假设聚合物链是均匀带电的，每一个分子的带电量均为αe.连接成链的单体具有FENE势能： 　　（1） 　　其中，r是相邻两个单体之间的键长，是链可以延伸的最大长度，为劲度系数。 　　所有的单体均应满足LJ势： 　　（2） 　　其中，仍然反映键长，表征了LJ势的强度。因此，单体之间的相互作用是可以忽略的。 　　该链还将受到因角度改变而引起的势[11，12]， 　　（3） 　　其中，是两个单体之间的角度。平衡角度被设置为，对于半柔性链，弹性系数定为2。粒子之间也会通过库仑势相互作用，该势被写为 （4） 　　其中是热能，为Bjerrum长度。是真空介电常数，是溶液的相对介电常数。溶剂分子没有明确的表示势的函数，然而存在摩擦力，这是由于粒子在溶液中运动及其受到的随机力造成。 　　由于热浴运动带来的力隐藏在郎之万方程中， 　　（5） 　　其中，m是粒子的质量，是摩擦力系数，是保守力。随机力满足波动耗散定理。该系统被置于形式的交流电场中，是电场的频率，是电场的振幅，电场的方向为z轴方向，如（图1）所示。 　　（图1）两端在x轴上被锚定在原点和x=d处的单PE链。交流电场的方向为z轴方向。实点表示链的中点（中间的单体）。 　　线性链的长度选择为。链两端固定在x轴的点选为（0，0，0）和（d， 0， 0）。我们将研究d/L=0.9的情况。设定如下关系：，，，以及，并且设定为单位时间。电场强度将被用单位描述。 　　接下来，对于柔性和半柔性PE链，将在一个比较大的范围内变动，我们将研究其对链的中点位置z（t）的影响，以及外部电场E（t）。郎之万方程中的Verlet算法每一步的时间要求为。系统需要先进行MD时间以达到平衡，并进入静止系统集合（），紧接着进行MD步骤以得到100个粒子的轨迹。 　　3 结果 　　三维的柔性和半柔性的PE链如（图2）所示。无论是柔