《试验十七高分子链构象的计算机模拟试验目的对高分子链构象》由会员分享,可在线阅读,更多相关《试验十七高分子链构象的计算机模拟试验目的对高分子链构象(5页珍藏版)》请在金锄头文库上搜索。
1、实验十七 高分子链构象的计算机模拟一、实验目的1对高分子链构象有直观、形象的了解。2了解 Monte Carlo 方法的原理。二、实验原理l 高分子的链构象 链构象就是分子链在空间中的形状和尺寸,高分子溶液和本体的许多性质,诸如热力学、光学、电学、声学、流体动力学和力学等性质,都与链构象有关。而链的柔性, 使构象具有统计性。 可以说, 链构象理论是研究高分子的基础, 我们在高 分子物理课程中曾学到,高分子链构象可以分成两类,即理想链和真实链。理想链即无干扰状态下的高分子链, 它可以用无规行走来描述, 故又称无规链, 其均方末端距(h2)可以表示为22 (h2) = NL2 (1)式中 N 是链
2、节数; L 是链节长度。又M = MoN (2)式中M为高分子链的分子量; Mo为链节分子量。这样,由(1 )、(2)两式,我们 就可以得到均方末端距同链分子量 M 的关系。真实链即有排除体积效应的高分子链。由于原子之间相互作用的存在。两个原 子不能在空间占据同一位置。换言之,在分子链中,组成链的原子相互被排除于同 一体积。这种效应即称为排除体积效应。在链构象的统计理论和数学模拟中又称自 回避链。排除体积效应引起分子链扩张,即分子链有较大的均方末端距,表示为:(h2) = NvL v 6/5一般来讲,高分子链是真实链,但高分子链若处于0状态(在0温度、0溶剂中)则成为理想链。这时,从无规行走的
3、角度来看,无规行走的迹线可以前后重叠。本 实验的主要目的就是直观地从计算机上“看到”这两种链、同时,通过观察高分子 的链构象,找出高分子线团尺寸与分子链节数的关系,亦即均方末端距与分子量的 关系。2. 链构象的模拟 从高分子链的二次结构上看。链分子是由链节依照一定的键角和旋转角,一个 一个地连接而成的。一个具有 N 个链节数的高分子链,从计算机模拟的角度来看, 相当于一步一步走 N 步而成的。 不同的链构象是不同行走方式的必然结果。 高分子 链在某一时刻为何种构象完全是随机的,这就要求行走方式的选择必须是随机的。计算机对这种随机过程的描述,是采用Monte Carlo方法。Monte Carl
4、o方法又称计算机随机模拟方法,统计试验方法。简单地说,当我们欲计算某物理量时,先建立某种合适的物理模型,然后大量地随机取样,(通过某种随机过程获得样本,该样本同此物理量有直接关系,且可以得到此物理量的一个 数值结果)将这些样本的结果做统计平均。这个统计平均值就是我们欲求的物理量。 当然,样本数越多,这个平均值越接近真实值。一般地,我们要选取成千上万个样,才能获得较好的模拟结果。所以Monte Carlo方法是离不开计算机的,其详细的原理可见附录。在本实验中,所模拟的对象是柔性高分子链。其构象可以用高分子线团来描述。我们利用无规飞行模型和晶格模型两种物理模型,来模拟理想链和自回避链的构 象。(1
5、)无规飞行模型为便于观察,我们选取二维平面的无规飞行模型。在这个模型中,高分子链可 以朝任何方向行走,飞行迹线可以重叠。形象地说,这种飞行迹线就是布朗粒子的 行走迹线。这个过程,并将这成千上万个末端距相加,模拟的基本思想是:在平面上任选一点 为原点,以1为步长,利用计算机的随机数, 随机取定第一步的行走方向。记录下这次飞 行迹线;再以第二点为起始点, 利用随机数, 随机取定第二步的行走方向,记录下这次飞 行迹线,如此类推,我们就可以得到N步的飞行迹线,这就相当于一条链节数为 N的高 分子链的构象,记录下原点到最后一步的末 点的距离,即为这种构象的末端距。然后, 我们再重复上述过程,就可以获得链
6、节数N的高分子链的另一种构象,其末端距也会不 同;如此类推,一般我们要重复成千上万次 再除以重复过程数目,就得到链节数为N的高分子链的均方末端距。显然,这种模拟的结果,相当于键角可取任意值的自 由连接链,图17-1是一条无规飞行链的分子链构象。(2)晶格模型单链构象的Monte Carlo模拟,常常是采用晶格模型。这里我们就简单介绍二 维平面方格模型上的无规链和自回避链的模拟。无规链的模拟如图17-2所示的方格模型,其显示的为一次无规行走所产生的高分子链。为简单计,该方格边长为 1,行走迹线是沿正方形的边,每一步(即长度为 1的边长) 代表一个链节,链节长度就为1,不同于无规飞行模型,方格上的
7、每一步只有四种方向可以选择,见图 17-3。不同于无规飞行模型的模拟之处在于:行走方向只有种选择,至于每一步究竟取何种行走方向,仍由随机数随机取定,显然这相当于链 角受限的高分子链。自回避链的模拟图17-2“-:图17-4是一次自回避行走的链构象。X图17-3 辻心沁 :自回避链同无规链相比,主要区别就在于模型中链是不能相互交叉重叠的,即 每一步的坐标不能有相同者。因此,在编程序时,只要排除那些坐标与原先走过的 坐标相同的选择走向。其余的模拟思想仍同前。 具体程序参见附录。三、实验步骤1 熟悉计算机的使用(1) 打开电源;(2) 在A驱动器中插入启动盘;(3) 开启计算机开关(先开启屏幕开关)
8、;(4) 计算机处于“ A ”状态,输入程序名“TEST”再按“RETURN ”键,即可开始实验。(5) 实验完毕后,先关电源,再关屏幕。2.观察无规链和自回避链的行走(1) 选择(1)是无规飞行链,步数自己设定,(但不要超过200步)*(2) 选择(2)是在平面方格上的无规行走链,步数自己设定(但不要超过200步)(3) 选择(3)是在平面方格上的自回避链,步数自己设定(但不要超过100步)以上每项内容,自己设定五个不同的步数N,而每种N值又要重复20次,并记录每次的末端距 h2 ( N),数据记录、处理类似于下表:N =前123i57g1011-20hl3. 无规链的链末端距的计算50步)
9、(1) 选择(4)是在立方格子上的无规行走链,步数自己设定(但不要超过(2) 记录步数N与相应的链末端距(h2)、N(3) 做In ( h2) NI nN的双对数图,并用最小二乘法拟合求出斜率。思考题1 对最小二乘法拟合的结果作一讨论,实验误差的主要来源是什么? 2你能举出一两个例子说明Monte Carlo方法在高分子中的应用吗?3你能否编写一个在平面三角格子之上无规行走链的BASIC程序吗?参考文献1 吴大诚,高分子构象统计理论导引,成都,四川教育出版社,1985。2 B. Alder, S. Fembach, M. Rotenberg. Methods in Computational
10、Physics, 1963, 217242。*两个原因:一是屏幕限制,步数太大,迹线会超出屏幕。二是时间限制,步数太 大,时间会较长。当然,我们也可通过缩小步长以使屏幕上显示出较多步数。附录1 Monte Carlo (蒙特卡罗)方法Monte carlo是摩纳哥的一个城市,以赌博闻名于世,赌博是有随机性的。因此,在数学上就借此来代表一种用人工抽样试验来估计数值并求解未知量的计算方法; 即按照某一统计量在样本中的平均值求解所需的估计值,故也称统计试验方法,或 计算机随机模拟方法。这种方法可用来解决那些无法求分析解或数值解的复杂而困 难的问题。我们先通过下面的例子,来了解Monte Carlo方
11、法的基本精神。考虑一个射击运动员的射击成绩G,令X表示弹着点到靶心的距离,g (X)表示得分,f (X)表示该运动员的弹着点分布的密度函数,则G (X) = g (X) f (X) dX另一方面,如果该运动员进行了实弹射击,弹着点依次为兄X1, X2Xn,则平均得分为- 1 NGNg(Xn)N心很明显,Gn是G的一个近似估计,Monte carlo方法计算积分 G,正是用Gn作为G的近似估计。由此例可以看出,Monte carlo方法解题的一般过程是,首先构成一个概率空间,然后在该概率空间中确定一个依赖随机变量X (可以为任意维)的统计量g (X),其数学期望Eg =g X dF X正好等于所
12、要求的值 G,其中F(X)是X的分布函数;最后,产生随机变量X的简单子样Xi, X2Xn,用其相应的统计量 g (Xi),g (Xn)的算术平均值1 NGNg XnN n4作为G的近似估计。换句法说,Monte Carlo方法的最低要求是,能确定这样一个与计算步数N有关的统计估计量 Gn,当N ;匚时,Gn依概率收敛于所要求的值 G,亦即,对于 任意的 & 0,应有lim P |GN -G|: ; =11 N 1 f其中P(A)表示事件A的概率。用Monte Carlo方法解决具体问题的一般程序是:先将此问题用适当的物理模型描述,再选取适当的数学模型,然后通过大量地随机取样,求解此数学问题,从
13、而就解决了这一物理问题。这种方法的特征是采用人工随机取样试验,在取样中最关键的是试验的随机性,所以此法的第一要求就是提供“高质量”的随机数。由具有单位均匀分布的总体中所产生的简单子样称为随机数序列,其中的每一个体称为随机数。在计算机中一般是一种余数法来产生随机数,这种方法简述如下:令:R=(kR 1+C hd N 0兰 RcN式中符号mod意思是用 k -1 C除以N,然后只取余数,按上式此余数即为 Ni; 于是,由某一随机数 Ri一 1就可产生下一个随机数 Ri,这就是递推公式的涵义。例:令N = 8, k= 5, C= 3及R0= 4则得到一个随机数序列为4, 7, 6, l, 0, 3,
14、 2, 5, 4, 7继续算下去,我们会看到这个随机序到实际上是4, 7, 6, 1 , 0 , 3 , 2 , 5这8个数为周期的无限循环。随机数,顾名思义是要求要有相互独立性和均匀分布性的。显然,这个随机数序列的质量不高。在实际应用中N参数是很大的数,以保证产生的随机数序列在所有随机取样的样本中,能较好地满足相互独立性和均匀分布性。一般计算机中的随机数 (函数名为RND )能适合我们一般计算的要求, 它给出的是 一个0, 1区间均匀分布的随机序列。(只要将Ri除以N,即可由随机数 Ri成为0,1 区间的随机数)Monte Carlo方法已经在自然科学的许多领域得到了广泛的应用,在高分子化学和物理上,也取得了许多成果。