单体慢滴加合成超支化大分子的计算机模拟

1、单体慢滴加合成超支化大分子的计算机模拟第一章 综述 1.1超支化大分子合成及性质概述：在高分子材料的应用中 , 大多采用线形聚合物作为主要材料。 但是线形聚合 物在高分子质量时所具有的高粘度 , 使其无法应用在一些要求低粘度和材料性 能优异的领域。超支化聚合物是这样一类聚合物：由ABX支臂原料和核组分(可 以不加)反应而成，其中A和B是具有反应活性的官能团。最终形成具有众多端 基的聚合物。超支化聚合物与线形聚合物比较具有显著的优越性能 , 正成为人们 所关注的焦点之一。1.1.1 超支化大分子的相关背景：有机物分子的形状是决定其性质的一个重要因素。在过去的 15年, 科学家们 ,尤其是聚合物化

2、学家 , 介绍了一种新的 “树状分子” 。由一系列支化单元组成的 树状支化大分子可分为树枝状大分子 (Dendrimer) 和超支化聚合物 (Hyperbranched polymer) 两大类。树状大分子的合成为了控制分子的尺寸和 形状, 通常需要多步反应 , 并且每步骤都需要采取严格的保护 / 去保护措施和细 致的提纯 , 制备的困难导致价格的昂贵 , 限制了其作为消耗性材料的应用。而超 支化聚合物无需仔细分离提纯 , 可直接由本体聚合制备 , 即一步法合成。超支化聚合物是大分子科学中一个年轻而快速发展的领域, 是一类新崛起的具有特殊性质的聚合物。与传统线型聚合物或星型聚合物不同, 超支化

3、大分子具有三维的结构和大量的端基以及高溶解性、 低粘度、强化学反应活性等性质。 特 别是由于在所有超支化大分子合成方法中聚合物不需要进行纯化或很少需要纯 化, 因而制备简单、成本较低 , 有利于大规模地合成的特点 , 超支化大分子已显示 出诱人的应用前景。从第一次有意识地成功合成超支化聚合物至今已有十多年 , 这一领域取得了突破性进展 , 成为合成化学中的一个新的热点。 毫无疑问 , 了解高 度支化聚合物的制备方法和应用前景对超支化聚合物的发展和开拓其新的应用 领域具有重要的意义。1.1.2 超支化大分子的结构特征：相对于普通的线型高分子及树枝状高分子 , 超支化高分子具有独特的结构 , 这3

4、 种类型的聚合物的结构示意图如图 1 所示：由图 1 可见, 超支化聚合物的分子中只含1个未反应的A基团,而含多个未反应的B基团。超支化聚合物与树枝状 大分子一样，单个分子的形状是球形的，但是树枝状大分子的分子具有完美的分 枝结构,整个分子中无缺陷，因此，树枝状大分子的图1树枝状人分子与超支化聚合物分子结构示意图分子是圆球形，而超支化聚合物的分子中有缺陷，整个分子并不完全对称，所以, 超支化聚合物的单分子形状是椭球形，但这两种结构的分子表面均密布着大量有 反应活性的末端官能团。传统的线型高分子在无外力作用下总是自发地呈蜷曲形 态。当与线型高分子具有相同的端基数目时，超支化高分子的多端基结构决定

5、了 它的无链内缠绕性。超支化高分子具有分子量分散性。而且由于超支化高分子 通常是由ABx型单体自聚合得到，所以超支化高分子的分子量分布比普通的线型 高分子宽。超支化高分子分子量分散性可通过采用多官能度的核，并降低核的浓 度以及采用缓慢滴加单体的方法来进行控制。异构也是超支化高分子与树枝状大分子及线型高分子之间的一个显著差异。因为每个单体的滴加都是随机的，因此，即使对于指定了分子量和分支程度的超 支化高分子，也会有大量异构体出现，并且异构体含量随着单体的复杂程度和高 分子分子量的增大而增加。这种几何异构会影响聚合物的溶解性和固态堆积方式 及其它相关性质。超支化大分子的性质：超支化聚合物的特殊结构

6、决定了它具有与普通线型高分子不同的特殊性质。1.121低粘度超支化聚合物和树枝状大分子的最突出特点是它们表现出惊人的低粘度。童身毅等合成的超支化聚酯2酰胺型两种树脂HAS和HAP ,与市售的聚酯树脂均按60 %二甲苯溶液进行粘度对比，25 C测得的粘度分别为21mPas、12 mPas和6500 mPa - s ,这表明超支化聚合物的溶液粘度低得多。这是由于流体的粘度来 源于聚合物的流体力学体积和分子间的内摩擦，超支化聚合物的分子尺寸小，有 大量短支链存在，以及分子链本身及分子之间无链的缠绕使得分子间相互作用力 小,因而粘度较低。树枝状大分子、超支化聚合物及线型聚合物的特性粘度与相 对分子质量

7、之间的关系如图2所示：图2三种集合物粘度随分子量变化示意图由图2可以看出，树枝状大分子的特性粘度随分子质量增大而增大，经历一个最高值之后下降，超支化聚合物的特性粘度随分子质量增大而增大，但比线型聚合 物的粘度小许多。线型聚合物的特性粘度随分子质量呈线性增大，直到临界分子质量时粘度极度增大。通常在分子量与粘度之间存在着一定的关系，可以表示为Staudinger2krhn 方程：n = K - Ma其中：n 为溶液的特性粘数，K为常数,M为高聚物的分子量，a是与高分子链 在溶液中的形态有关的参数。对于线型高分子而言，其a值为0. 51 ,超支化聚合物的a值一般处于01210. 44之间,这是因为超

8、支化聚合物分子在溶液中呈 球形,无链缠绕,在溶液中虽存在着溶胀现象，但溶胀前后的体积变化较小。良好的溶解性与相对分子量相近的线型大分子相比，超支化聚合物的溶解性有很大的提 高。例如,超支化聚苯和芳香聚酰胺可溶解在有机溶剂中，而对应的线型聚合物则 由于主链的刚性，在有机溶剂中几乎不能溶解。超支化聚合物的分子相对于无规 线团来说,分子结构较紧密，而且在溶液中虽有溶胀现象，但溶胀前后体积的变化 相对较小,分子尺寸虽然相对于无机分子来说是“大分子”，但是比起在溶液中溶 胀的无规线团来说应算是“小分子”了。因此，完全可将超支化聚合物的分子视 为特殊意义上的小分子 , 它所表现出来的牛顿流体行为也可得到解

9、释。1.1.2.3 多官能性超支化聚合物的表面有大量的官能团存在 ,其封端官能度非常大 , 一般为12、 16、32 , 如果保留反应活性基团 ,则反应活性非常高。利用超支化聚合物终端的 多价态羟基基团的氢键作用 , 可在金属表面形成良好的功能化膜。通过对超支化 聚合物端基官能团的改性可以赋予其各种各样的功能。 对第二代超支化聚酯可进 行不同的终端改性 ,当终端接枝甲基丙烯酸时 , 所得聚酯为黄色的粘稠液体 , 可利 用辐射固化 ; 当终端接枝邻苯二甲酸酐时 , 所得的产物是微黄色晶体状固体 , 是一 种非辐射固化的聚酯材料。1.1.2.4 其他性能超支化聚合物具有良好的流动性 , 容易成膜。

10、复合树形分子气体分离膜 , 该膜 是由具有硅氧烷超支化聚合物结构侧基的乙烯基聚合物和热塑性树脂组成, 具有优异的气体渗透性 , 如具有优异的氧气渗透性。此外它还具有阻水性、可模塑加 工性、耐热性和耐候性等。因此 , 在医疗卫生用膜、氧吸收剂、气味吸收剂的填 充材料中都有着广泛的应用。超支化聚合物可以将一种非极性的内层结构与另外一种极性的外层结构结 合在一起 ,比如, 憎水的内层结构和亲水的端基。有研究指出 , 具有酯端基的超支 化聚苯可作为一种单分子胶束 , 酯端基使得该聚合物是水溶性的 , 而憎水性的内 层则可以捕获客体分子。超支化聚合物在溶液中的尺寸受溶液的极性和 pH值等因素的影响。有研

11、究表 明,当溶液的pH值发生改变后,羧端基的树形分子在水中的尺寸可以增加或减少 50 %。另外, 利用超支化聚合物的结构特点 , 通过适当的物理或化学手段 , 还可以 赋予超支化聚合物其他一些特殊性能 , 如光物理及光化学性能、吸附及解吸附性 能等。1.1.3 超支化大分子的合成：缩聚反应是制备超支化聚合物最常用的合成方法 , 通常可分为无控制增长 “一步法”和逐步控制增长“准一步法”。 “一步法”指由 ABx 型单体不加控制 一步反应 , 即一次性将所需的核组分及支臂原料、催化剂投入反应釜合成即可。其优点是合成方法简单 , 一般无需逐步分离提纯 , 且聚合物仍可保持树形大分子 的许多结构特征

12、和性质。但是一步法合成超支化聚合物产率不高 , 支化产物具有 随机支化的特点 ,相对分子质量分布较宽 , 分子质量无法控制。目前已用该方法合 成出一系列超支化大分子。 “准一步法”是指添加 By 型分子为中心“核” , 即将 成核原料与部分支臂原料、 催化剂反应一段时间后 , 再将剩余原料和催化剂加入。 引入By官能团不仅可以控制产物分子量，而且产物分散度也会大大降低。加成反应使合成超支化聚合物的单体数大大增加 , 其缺点为较难控制聚合 度和支化度。在活性自由基聚合中 , 单体既是引发剂也是支化点 , 乙烯基单体在 外激发作用下活化 , 产生多个活性自由基 , 形成新的反应中心 , 引发下一步

13、反应。Hawke等采用含有一个可聚合的苯乙烯和一个具有连接在取代的苯碳原子 的氮氧基团的单体的自缩合自由基聚合合成超支化聚合物。 由于氮氧键的低离解 能，使得热活化自缩合自由基聚合成为可能。Gaynor等采用对一（氯甲基）苯乙烯在铜（I）和,2,2 '-二吡啶存在下通过原子转移自由基聚合“一步法”合成超支化聚苯乙烯。Sakamoto等采用2-2-甲基-1-（3-异基硅氧）-1-丙烯基氧甲基 丙烯酸乙酯为单体 ,在亲核催化剂作用下 , 通过自缩合基团转移聚合合成了超支 化聚甲基丙烯酸酯。以CNM谱表征其结构，该超支化聚合物支化程度取决于用来 激活单体的亲核催化剂。还有人用类似的方法合成了

14、表面嫁接超支化聚合物。 近几年来 , 随着可控的开环聚合技术同支化技术相结合 , 不断出现了许多开环 聚合的新单体,并且涉及到了许多新机理。Cha ng等利用AB型单体,采用质子转移 法开环聚合合成了带有大量末端环氧基的芳香超支化聚醚。合成超支化聚合物的另一种方法是接枝法， 即在已有聚合物的接枝链上再进 行接枝 , 形成发散的扇形树枝状超支化聚合物。合成的扇形树枝状超支化聚合物 是在聚苯乙烯上接枝氯甲基阴离子来制备的 , 这种方法制得的超支化聚合物官 能度大于 5000,分散系数小于 1.3。聚合物的硬度可以通过调整接枝点在聚合物链 上的位置来改变。 静态和动态光散射实验都表明在稀溶液中聚合物

15、分子的表现类 似于硬球。在稀溶液中 , 聚合物的硬度效应随支化度增加而增加。E（封端的趨支化聚合物）图1超支化聚合物的合成超支化大分子的应用1.141在涂料工业中的应用：超支化聚合物在涂料中的应用特别有意义，由于分子中带有许多官能性的端 基,可以有目的地对其进行改性，而且超支化聚合物具有紧凑的分子结构，不易发 生分子间的缠结，当相对分子量增加或浓度增高时，体系仍能保持较低的粘度，从 而使其具有独特的流变性质、很好的成膜性以及极佳的抗化学性、耐候性和力学性能,可用于制作光固化涂料、粉末涂料、无溶剂涂料、高固体组分涂料以及涂 料中的添加剂。如一类改性的甲基丙烯酸酯化超支化聚（胺-酯），该类聚合物可

16、作为光固化涂料的齐聚物，在光引发剂存在时和紫外光辐照下发生快速聚合反应 反应速率和双键转化率随支化分子的结构不同而变化，所得的高度交联的涂膜具 有较好的热学性能。由于其组成的涂料体系粘度低，可以不加或加入少量稀释溶剂,是一种绿色环保材料，且固化膜力学性能优良。聚合物共混改性剂：超支化聚合物具有较低的熔融粘度 , 在传统聚合物加工性能改善方面存在着 巨大的潜力。Hong等人采用超支化大分子作为低密度聚乙烯的加工助剂，发现混 合物的粘度下降 , 在固定挤出速度下流动加快 , 因而挤出所需能量明显降低 , 同时 消除了以前加工过程中经常发生的鲨鱼皮现象。 有人利用所制得的超支化聚合物 作为PMMA勺

17、粘度调节剂，解决了由于超支化聚合物端基为极性基团而与本体不 相容的问题。超支化聚合物与其他聚合物共混时可以使共混物的热稳定性和力学 性能发生很大变化。1.1.4.3 作为环氧树脂的增韧剂： 环氧树脂由于具有优异的力学性能和加工性能而广泛应用于复合材料 , 然而其韧性低而常常限制它的应用。 超支化聚合物具有低粘度、 多官能度和分子链无 链缠结等优良性能 , 有可能作为环氧树脂基复合材料的增韧改性剂。由于超支化 聚合物的极性可以通过对端基的改性进行调整 , 这就有可能使之与其他聚合物相 容或不相容。研究表明 , 在环氧 / 多元胺热固性树脂体系中加入调节了极性的超 支化脂肪族聚酯可以发生诱导相分离

18、 , 所得到的热固性树脂的韧性大幅度提高 , 并保留了原热固性树脂的高模量。1.1.4.4 在药物缓释剂中的应用：超支化聚合物作为药物载体的研究较多 , 可望用于农业、化妆品和医药行业。 使用超支化聚合物作分散剂可使油相在水相中稳定达6个月,粒径在300400nm采用生物相容且低毒性的原料制成超支化药物载体，载体在24h内释放药力, 随着载体相对分子量的增大 , 药力释放更缓慢、更均匀。合成的超支化分子 , 核组分为憎水基团 , 较易与憎水药物相容 , 而支臂的亲水性则促使憎水药物在体 液中的溶解。控制超支化分子尺寸与支化度 , 可生成可控药物缓释剂。超支化聚合物还可作为大分子催化剂和大分子引

19、发剂 , 由于其特殊的结构和 性能 , 超支化聚合物在聚电解质、非线性光学材料、电致发光材料、超支化分子 液晶、固体粒子表面改性等领域也有了广泛的应用 , 并逐渐成为研究的热点。1.2关于Monte Carlo方法的相关介绍：1.2.1Monte Carlo 方法的相关历史Monte Carlo 模拟方法又被称为随机抽样法、随机模拟法、统计实验法。它 起源于 20世纪40年代Los Alamos的物理学家 Von NeumannMetroplis、Ulan和Kahn 对中子扩散问题的模拟 , 即通过计算机对中子行为进行随机抽样模拟 , 然后从大 量的中子样本参数中得到所需的结果。Monte C

20、arlo 方法不仅能用于解决确定性问题，更适合于求解随机性问题， 这时往往对所考虑的问题进行直接模拟。这种情况下， Monte Carlo 方法不仅仅 是一种数值计算方法，而往往是对实际问题的直接模拟。因此，在物理、化学领 域中常称其为“计算机模拟”。近年来，“计算机模拟” 已逐渐成为继“理论科学” 和“实验科学”后科学的第三分支，补充了传统的科学方法理论方法从Monte Carlo方法诞生以来,随着电子计算机的迅速发展,人们开始有意识 地、广泛地、系统地应用随机抽样方法来解决大量数学和物理等方面的问题 , 并 且将Monte Carlo方法作为一种独立的计算方法来进行研究,也使得Monte

21、Carlo 方法向各个学科领域渗透。除了对中子扩散等所谓的粒子输运问题的应用之 外,Monte Carlo方法还被广泛的用来模拟各种平衡态和非平衡态的经典统计力 学问题 , 以及量子统计等问题。Monte Carlo 模拟如此吸引人，并得到如此发展的原因，并非仅仅是因为它 作为理论和实验的补充。必须强调， Monte Carlo 模拟可以提供体系的实验可测 的物理量和通过实验无法测量的物理量。 此外，实验科学家可以改变样品的温度、 压力等外部因素，但是无法改变体系的各相互作用势，而在 Monte Carlo 模拟中 则可以任意地改变体系的相互作用势而不造成任何形式的困难。 实验化学家往往 通过

22、改变化学体系中各组分的初始浓度和相对配比来研究化学反应的宏观动力 学，而这在 Monte Carlo 模拟则可以轻而易举地实现。即便对于实验中难以达到 的某些极限条件，在 Monte Carlo 模拟中也无任何困难而言。1.2.2Monte Carlo 方法的基本思想和特征Monte Carlo 方法能用于求解确定性的数学问题 , 但它更适宜于求解随机性 问题。它的最基本思想是 :为了求解数学、物理或化学等问题 , 首先需要建立一个 概率模型或随机过程 , 使它的参数等于问题的解 ; 当所解的问题本身属随机性问 题时, 则可采用直接模拟法 , 即根据实际物理、化学情况的概率法则来构造 Mont

23、e Carlo 模型; 然后通过对模型或过程的观察或抽样试验来计算所求参数的统计特 征;最后给出所求问题的近似解。 在高分子科学中的 Monte Carlo 模拟主要采用直 接模拟方法 , 因为高分子科学中的很多问题本身就属于随机性问题。Monte Carlo 方法有以下基本特征 : 由于Monte Carlo方法是通过大量简单的重复抽样来实现的,所以Monte Carlo 方法及其程序的结构十分简单 ;与一般数值方法相比，Monte Carlo方法的收敛速度较慢,因此Monte Carlo方 法一般用于求解数值精度要求不太高的一类问题 ;(3) Monte Carlo方法的误差主要取决于样本

24、的容量 N,而对样本中元素所在的空 间无关，即Monte Carlo方法的收敛速度与问题的维数无关，故更适合多维问题的 求解;(4) Monte Carlo 方法对各种问题的适应性很强 , 因为它对问题求解的过程仅取决 于所构造的概率模型 , 这正是 Monte Carlo 方法得以广泛应用的一个原因。关于Monte Carlo模型的概述Monte Carlo 模拟的系统多采用由立方格子整齐地堆积而成的模型 , 即格子 模型。比如50X 50X 50,32 X 32X 256等,每一个格子只可容纳一个高分子链节或 溶剂分子等 , 或代表一个空格。 在模拟过程中 , 可以用以下一些方法来实现高分

25、子 链形状的改变：曲柄运动、摆尾运动、L形翻转、蛇行运动、RRU除一生长法、空 穴扩散法、键长涨落模型等。前三种链节运动方式都只能使高分子链局部运动 , 不能整体移动高分子链 , 故效率较低。蛇行运动算法可以使分子链整体移动 , 但也 不过只改变了首末位置 , 而未改变中间链节的构型。这四种算法适合模拟浓度较 低的系统。 在高浓度系统中 , 运动效率将大大地降低。 故在此基础上 ,Reiter 等提 出了空穴扩散法 , 他将运动的主体定为空穴或溶剂分子 , 而不是高分子的链节。 后 又进一步提出了集体移动法 , 其算法是先将高分子链的一端的链节去除变成一个 空穴, 再应用空穴扩散法改变构型 ,

26、 直至这个空穴又回到作为运动主体的高分子 链的原来的一端。最后又将这个空穴变回原来的链节。在键长涨落模型中 , 链节 间的键长可以伸缩 (比如可以在 1和1/2 之间变化 ), 键角也可以有更多的变化 (比 如90°,180 °,45 °等 )。陆建明提出的中间蛇行运动的算法是结合空穴扩散法和 键长涨落模型法后的一种算法。 其原理是随机选择一个空穴和此空穴邻位上的一 个高分子链节 , 并相互交换位置 , 计算此运动链节与前后相邻链节间的键长 , 若两 个键长的值都不大于 2则实现此构象的改变。 如果只有一个键长大于 2则将此键一 侧所有链节实行蛇行运动 , 即“中

27、间蛇行运动” 。如果两侧键长都大于 2则不能实 现新的构象 ,而维持原构象 ,此算法运动效率较高。Monte Carlo 方法并非仅仅是一种简单的数值计算方法 , 而往往是对实际问 题的直接模拟 , 因此常常被直接称为“计算机实验”。它已经成为与理论方法和 实验方法并立的科学研究的第三分支。 “计算机实验” 之所以占有如此重要的地 位,是因为它可以获得关于所考察体系的任意详尽的信息。虽然Monte Carlo模拟的结果总是伴随着统计误差 , 但从原理上来看 , 只要增加统计样本数就可使误差 减少到所希望的程度。 相反, 解析方法真正可解的实例并不多 , 即使可解也往往带 有某种近似 , 并且所

28、引入的近似的合理性常常事先不能完全知道。 同样, 由实验方 法所获得的关于研究体系的信息也会存在不可靠性。因为在通常情况下 , 人们在 实验前不可能完全了解体系内的各相互作用量 , 因而对实验的设计和实验过程中 各项条件的控制就不可避免地带有盲目性。 甚至有时无法确定实验观察到的现象 究竟是体系的实在特征还是由实验体系中的某些未知杂质或一些未加控制的实 验条件所导致。Monte Carlo方法与一般实验不同，模拟出的各种微观状态足以描 绘出体系的静态和动态行为 , 可以在分子水平上跟踪体系演化过程。 可以说 , 只有 计算机模拟才能达到对一个模型体系进行精确定义 , 才能真正自如的控制体系的

29、各种内部和外部参数 , 从而来考察所感兴趣的因素与模型体系的各种统计性质间 的关系。1.2.4Monte Carlo 方法的意义Monte Carlo 模拟在高分子科学中的应用及其广泛。 首先是因为在高分子科 学中存在着大量可供 Monte Carlo 直接模拟的随机问题。在高分子化学的聚合 物合成中， 存在着大量的不确定性问题。 高分子链有大量的重复单元构成， 分子 量一般在 1,000-1,1000,000 之间。当采用两种和两种以上单体进行共聚时，可 以形成共聚高分子链。 根据聚合方法和单体性质的不同可以生成前段和无规共聚 物等不同的链结构。 由于聚合反应本身的随机性特点， 高分子系统内

30、各成员之间 存在着与其生成机理密切相关的特定分布， 也就是体系中生成的高分子链并非具 有相同的分子量， 而是存在着分子量分布问题。 在多元聚合体系中， 生成的共聚 物不仅具有分子量分布， 而且存在着不同单元在高分子链上的排列问题， 即序列 分布。此外，在多官能团的聚合反应中的枝化问题和凝胶化问题， 高分子链的降 解问题等等，无一不是随机性问题。这使得 Monte Carlo 模拟成为研究高分子反 应的重要手段。对于高分子物理问题， Monte Carlo 模拟有着更深刻的意义。一般高分子可 因热运动而绕其化学键做内旋转， 使高分子链的形状不停的发生变化， 这导致单 个高分子链的构象统计成为十分

31、复杂的问题。 而高分子链的构象或链的形状强烈 依赖于溶剂的性质、 温度等环境因素， 更增加了问题的复杂性。 高分子单链的排 除体积( excluded volume )问题(即自规避行走问题， self-avoiding walks ) 就使高分子科学家困扰多年，直到重整化群理论( renormalization group theory )方法的引进才得到较为圆满的解决。由结构和性质如此复杂的个体所堆砌而成的高分子浓溶液， 乃至高分子本体 的多链体系则具有更复杂、 更深刻的统计内涵， 但同时也给高分子体系的统计理 论研究带来了巨大的困难。这恰恰为 Monte Carlo 方法提供了很好的研究

32、对象。 正是由于高分子科学上的上述复杂性和困难， 几乎在 Monte Carlo 方法刚诞生时， 它就在高分子科学中得到了应用。 Monte Carlo 模拟在高分子科学中的应用的先 驱性工作是 Wall 在 20 世纪 50 年代为研究高分子链的排除体积问题所进行的 Monte Carlo 模拟。从此以后，高分子科学中的 Monte Carlo 的发展与高分子科 学的发展密切相关， Monte Carlo 模拟研究广泛涉及了高分子化学和物理的各个 方面，并取得了丰硕的成果， 这也对现代高分子科学理论的建立和发展起到了十 分重要的推动作用。1.3 键涨落格子模型概述1.3.1 键涨落格子模型的

33、提出Crabb和Downcy等指出，在运用格子链模型对高分子链的动力学行为的模拟研 究中，格子的配位数 z 起着十分重要的作用。一般来说，格子的配位数越大，则 链的局部构象更易松弛。目前大多数的 Monte Carlo 格子模型均是采用修正的 Verdier-Stockmayer 模型。人们提出了各种微松弛模式，例如单链节微松弛模 式、曲柄( C ra n kshaft )微松弛模式等等。 然而，正如 Carmesin 和 Hilhorst40 等指出的那样， 上述链模型的微松弛模式均不能运用到二维格子链体系， 原因在 于，两维空间中的上述微松弛模式存在严重的“阻隔效应” ，它不能将某些局部

34、构象有效地松弛掉而人为导致过慢的松弛特征。另外，即使在三维方格子中， Verdier-Stockmayer 算法所能实现的最大支化度也只有 6，而且支化点无法运 动。同时， Monte Carlo 链模型中的链单元并不真实代表一个链单元，而是以真 实高分子链段的粗粒化近似。 因此， 链单元之间的键长不应是一个常数， 链的轮 廓长度( contour length )应该可以在一定的长度范围内发生涨落，而一般的格 子链模型无法真实地模拟这一点。更重要的是， Verdier-Stockmayer 算法被认 为是“非各态历经”的。这一点对高浓度体系的模拟十分重要。人们对具有键长 涨落特性的 Larson 模型也指出了同样的问题。 针对上述缺点 Carmesin 和 Kremer 提出了二维格子链体系下的“四位置”模型(four site model),后来也称其为“ 键 涨 落 模 型 ( bond fluctuation model )”。 这 个 模 型 克 服 了 Verdier-Stockmayer 模型的上述所有缺点,而且更符合实际情况。在此基础上, Deutsch 和 Binder44 将其应用到三维体系中, 相应的模型称为 “八位置” 模型( eight site model )。1.3.2 三维八位置键涨落