模拟囊泡形变动力学的新方法离散空间变分法3.pdfVIP

第 54 卷 第 9 期 2005 年 9 月 物 　理 　学 　报 Vol. 54 ,No. 9 ,September ,2005 ( ) 10003290200554 09 400006 ACTA PHYSICA SINICA 2005 Chin. Phys. Soc. 模拟囊泡形变动力学的新方法离散空间变分法 李剑锋 　张红东 　邱 　枫 　杨玉良 ( 复旦大学高分子科学系 ,教育部聚合物分子工程重点实验室 ,上海　200433) ( ) 2004 年 11 月 11 日收到;2004 年 12 月 20 日收到修改稿 　　提出了全新的离散空间变分法研究二维囊泡形变动力学过程 ,克服了原先解析或数值解囊泡形状方程遇到的 困难 ,结果表明动力学终态与已有理论方法的计算结果完全一致 ;说明这一方法正确而且解法稳定有效 ,这一方法 可进一步推广至研究三维无对称性的囊泡形状以及膜与膜之间有长程相互作用的情形 ,这为囊泡形状的研究提供 了新的理论手段. 关键词 : 离散空间变分法 , 囊泡 , 动力学 , Helfrich 方程 PACC : 0420F , 0420J , 8710 , 8720 能在某些特定情况下得到有限的解析解 ,而对方程 1 引 言 离散求数值近似解时又会遇到稳定性和收敛性的限 制 ,它要求离散时间步长和空间步长都要无限小. 为 近年来 ,高分子自组装行为已经成为理论和实 此 ,我们设计了全新的离散空间变分法去重新考虑 验上都很活跃的研究课题[1 —3 ] . 而膜和囊泡就是由 这个问题 ,这种方法解除了空间步长大小的约束 ,只 表面活性剂、生物或合成大分子组成的超分子体系. 要时间步长足够小那么这种方法就是数值稳定的. 双亲性分子溶解在水溶液中能自发形成亲水端指向 本文用该方法分别考察了不可自交叉和自交叉的情 水溶液的双层膜 ,继而形成具有封闭性的囊泡 ,具有 形 ;虽然对于仅由油脂分子构成的双层膜来说膜与 丰富多彩的形态和形状[4 ] ,引起了人们的广泛兴趣. 膜之间的自交叉似乎是不可想象的 ,但也有人认为 自从Helfrich[5 ] 和欧阳钟灿在膜和囊泡的形变开创 如果膜上接上高分子后 ,其接技点也可被视作交叉 性的工作特别是成功解释红血球的形状问题以 点[ 17] ,而且在微分几何中自交叉曲面属于正规子流 后[6 ] ,大量相关的理论工作陆续涌现[7 —11 ] . 这些工作 形(regular surface) ,也是形状方程的一类解 , 因此作 [5 ] 为对新方法的系统验证 ,给出自交叉情形的计算结 在不同程度上都是基于 Helfrich 自由能方程 κ 果也是必要的. ( ) 2 Δ λ ( ) F = H - c0 dA + P dA + dA , 1 2 ∫ ∫ ∫ 由于模拟三维囊泡的变形太繁杂,不适合用来 式中第一项为弯曲弹性能 , κ是膜的弯曲模量 , H 是 描述一种新方法 ;所以为了集中体现离散空间