FENE链分子相图的分子动力学研究.pdfVIP

FENE 链分子相图的分子动力学研究 1 1 2 1 欧阳文泽 孙昭艳 吕中元 安立佳 1 （中国科学院长春应用化学研究所高分子物理与化学国家重点实验室，长春， 130022） 2 吉林大学理论化学研究所，理论化学计算国家重点实验室，长春，130023） （ 利用计算机模拟方法研究流体的气液相图，在科学及工业上都有着重要的意 义。模拟相图一般有几种方法：传统的模拟方法，它要求体系包含大量粒子，并 且对每个状态点都要进行单独的模拟，因此从算法的角度来看，非常不经济；吉 布斯系综蒙特卡罗方法（GEMC），该方法中某一给定温度下的气液相平衡可以利 用一次模拟直接完成，使计算时间大大缩短；与GEMC 相对应的分子动力学模拟 方法，使之具有GEMC 的优点并且可以描述体系的动力学。 本文中我们采用了一种链分子输运的分子动力学方法来模拟计算链分子的 相图。在模拟中，我们的链分子模型为珠簧模型 （bead-spring model）,即把链 分子看成一个个被弹簧连接的接点，键作用即同一条链内相邻的链节点间的作用 [1] 用FENE （finitely extensible nonlinear elastic）来描述 ，  k 2 r 2 − R ln(1− ) r R  0 2 0 UFENE ( r)  2 R0 ∞ else 而任意两个链节点的作用则为12-6 的lennard-Jones 相互作用 σ 12 σ 6 φ(r ) 4ε[( ) =−( ) ] ij r r ij ij 为体现相图的普适性，采用了约化单位(reduced units)，并以粒子半径σ，粒子 间相互作用能量ε和粒子质量m 作为基本单位。 模拟方法与过程简述如下：首先，两个箱中分别均匀放置N 个链分子，链的 各节点初始速度由Maxwell-Bolzmann 分布给出。然后每个箱进行 105 步NVT 模 拟（无链分子交换），以消除模型建构中形成的局部不平衡。接着允许链分子交 换，进行MD 模拟。模拟箱中的粒子 （各个链节点）位移由作用的运动方程决定， [2,3] 在三维周期性边界条件下，力与势能的计算采用最小镜像约定 。粒子速度及 位置的改变由差分法求解牛顿方程来实现，为保证模拟的稳定性，采用了跳蛙算 [3] [4] 法 (leap-frog algorithm)及Berendsen 热浴及压力浴弱耦合方法 。每2000 步模拟，被允许可以进行一次链分子的交换，以确保在链交换前模拟箱各自已达 到平衡。链的交换由两箱中链分子的化学势决定 （每次交换，从化学势大的箱中 随机选择一条链移入另一箱中）。这里，我们用Rosenbluth 方法来计算链分子的 [5] 化学势 。 模拟计算了不同链长n = 2 到n = 10 下的链分子相图。Fig.1 给出了n = 4 下的链