嵌段高分子尾形链构象性质的Monte Carlo研究

基于简立方格点模型对AB两嵌段高分子尾形链的构象性质及其链节的空间分布进行了MonteCarlo模拟.结果表明,链的尺寸、形状和链节的空间分布等统计性质和B链节与平面壁之间的相互吸引能有关.随着B链节与壁之间的吸引能的增加,链的尺寸和形状均呈现出先下降后升高的变化趋势,而且B链节的比例越大,这种变化越明显.     

高分子链构象性质和溶解性的Monte Carlo模拟

运用Monte Carlo方法对线型高分子链格点模型的构象进行了模拟,研究了单链体系的构象的尺寸（采用均方末端距,均方回转半径来表征）、形状（采用平均非球形因子来表征）和溶解性随溶剂与链段间相互作用能的变化情况.结果表明,、和随着εPS的增大而减小,具有相同的变化规律;随着Δε增大,溶剂变得越来越不利于溶解,高分子链的形状蜷曲为椭球形,高分子线团相互穿插交叠.     

金刚石格点上自避行走尾形链的构象统计理论

研究了金刚石格点上自避随机行走（ＳＡＷ）尾形链，采用精确计数和ＭｏｎｔｅＣａｒｌｏ模拟方法求得该ＳＡＷ尾形链的构象数Ｃ（Ｄ）１（Ｎ）和均方末端距［ｈ（Ｄ）１（Ｎ）］２及其分量随链长Ｎ的变化关系．发现它们与自由ＳＡＷ链一样都服从标度律，从这些量的计算机实验数据拟合求出了金刚石格点上ＳＡＷ尾形链的临界指数和格点指数．计算结果还表明短链ＳＡＷ在壁的法向与ＮＲＷ尾形链一样有所伸展，均方末端距的法向分量几乎是平行分量的２倍；但随Ｎ→∞，链自回避效应对壁的作用有所屏蔽．这些都与简立方格子模型上得到的结果一致．     

平行板间高分子链构象熵的动力学Monte Carlo模拟

用动力学Monte Carlo方法模拟了受限于两平行板之间的高分子链,并用扫描法计算了链的构象熵S,研究了构象熵相对于自由链的减小量(S0-S)与平行板间距D和高分子链长n的关系.结果证实了de Gennes的自由能标度关系,并给出了标度关系适用的范围.当D非常小时,高分子链受到强烈限制,S0-S与n成正比,表明单链节受到平行板的平均排斥作用力与链长无关.随着D增大,平行板对构象熵的影响越来越弱,单链节受到平行板的平均排斥作用力随链长的增长而增大.当D比较大时,平行板对构象熵的影响近似可以忽略,高分子链构象熵与自由空间中的结果一致.     

端基附壁高分子链的构象统计理论——Ⅱ.SAW尾形链

研究了SAW尾形链,采用精确计数和Monte Carlo模拟方法求得d维(d=2,3)SAW尾形链的构象数C_1~d和均方末端距[h_1~d(N)]~2及其分量随链长N的变化关系.发现它们与自由SAW链一样都服从标度律,从模拟数据拟合求出了这些SAW尾形链的临界指数.还表明短链SAW在壁的法向与NRW尾形链一样有所伸展,均方末端距的法向分量几乎是平行分量的2倍;但随N→∞,链自回避效应对壁的作用有所屏蔽,这与NRW尾形链有明显不同.     

磁性高分子链的磁性质和构象性质的Monte Carlo模拟

采用退火 (Annealing)MonteCarlo方法 ,从高温到低温顺序模拟了简立方格点上考虑最近邻Ising相互作用的磁性高分子链在不同温度的磁性质和构象性质 .磁性高分子链在低温下存在自发磁矩 ,无限长链的临界温度Tc=1 77± 0 0 5J kB.在临界温度附近 ,高分子链经历了从伸展的无规线团到紧缩球体的塌缩相变 .对链的尺寸、形状、近邻数及能量的分析表明 ,高分子链的构象性质从温度Tc=1 77开始发生较明显的变化 ,这表明高分子Ising链的相变是Ising相互作用和链节运动协同作用的结果 .     

聚硅烷链形状的Monte Carlo研究

用ＭｏｎｔｅＣａｒｌｏ方法对θ溶液中考虑二级相互作用的聚硅烷链的形状进行了研究。结果表明聚硅烷链的形状明显偏离球形，它与链长有关，长链极限的比值＜Ｌ２１＞∶＜Ｌ２２＞∶＜Ｌ２３＞约为１∶２．７∶１２．１。还发现聚硅烷链转动惯量的最长主轴与末端距矢量之间夹角的统计平均值＜θ＞的极限值约为２７°。     

基于格子链的缩聚反应的动态Monte Carlo模拟

采用描述自回避格子链的键长涨落模型, 以动态Monte Carlo方法对AB型单体的线型缩聚反应动力学过程进行了模拟. 通过该方法可以得到反应过程中链的瞬时构象, 还可以得到反应程度、聚合度、分子量分布及其随时间的演化. 模拟得到了合理的结果, 同时验证了无规线团尺寸与平均链长的标度关系, 表明该方法用于研究逐步聚合反应过程是可行的, 并且与一般的研究聚合反应的Monte Carlo方法相比, 还能够同时得到构象等空间信息. 还比较了不同大小的模拟体系所得到的分子量和多分散系数的异同, 讨论了有限元胞效应.     

端点附壁的高分子链形状的Monte Carlo模拟

分别基于简立方格点和四面体格点模型对一端吸附在无限大平面的高分子链（平面接枝高分子链）的形状进行了Ｍｏｎｔｅ Ｃａｒｌｏ模拟,结果表明,接枝高分子链的形状更偏离球形,〈Ｌ＾２１〉：〈Ｌ＾２２〉：〈Ｌ＾２３〉的极限值约为１：２．７５：１２．５,其中〈Ｌ＾２１〉,〈Ｌ＾２２〉和〈Ｌ＾２３〉分别为回转半径张量的本征值Ｌ＾２１,Ｌ＾２２和Ｌ＾２３（Ｌ＾２１〈Ｌ＾２２〈Ｌ＾２３）的统计平均;链长相同时,接枝     

模型高分子尾型链的构象统计理论

基于金刚石格点上受限于半空间的随机行走，借概率论中的反射原理推导出了模型高分子尾型链的构象分布．得到了链长为Ｎ的模型尾型链的允许构象数Ｃ１（Ｎ）和末端距分布函数．在四选择模型中，发现Ｃ１（Ｎ）／４Ｎ与Ｎ－１／２成正比．模型尾型链在与壁正交方向上的均方末端距分布与平行方向比较扩张至２倍，而后者与自由链对应的分量相等．这些解析结果得到了精确计数和ＭｏｎｔｅＣａｒｌｏ模拟计算机“实验”的支持．     

聚合物链构象的Monte Carlo模拟:单链的构象熵

本工作提出了一种新的直接用Monte Carlo模拟方法计算聚合物链构象熵,我们称这种方法为构象计数法。在本文中首先叙述了这个方法的计算原理,然后将它应用到简单立方格子中的无规自避行走问题上。具体计算了短链(链步数n=7～19)和长链(n～2100)时的熵在短链下所得到的熵数据与准确值相比偏差在0.04％以内(2000个样本)。在长链下所得到的熵数据与重整化群理论的结果相比,误差在0.8％左右(样本数为300)。同时本文还计论了限制在边长为2的立方体内,最大链步数n=26的受限链的构象熵。由本方法所求得的数据和用直接计算构象数得到的熵准确值也符合得很好。偏差在0.6％以内(2000个样本)。本文还对权重因子进行了讨论,提出改进数据准确度的方法。     

高分子链坍塌转变动力学过程的动态蒙特卡罗模拟

基于键长涨落格子模型和动态蒙特卡罗模拟方法,引入疏水相互作用,对均聚高分子链坍塌转变动力学过程进行了模拟研究.模拟发现,其坍塌时间呈高斯分布,而平均坍塌时间随淬火深度的变化类似于蛋白质折叠,但难以发现局部最小;另外,平均坍塌时间随链长呈指数形式变化.在其坍塌动力学过程中,高分子链构象先由橄榄球状演变为项链状,进而演变为香肠状,最后形成近球状的熔融球;基于团簇数目、团簇大小和非球面参数等参量,对前人提出的动力学过程四阶段划分进行了更为清晰的界定.     

DYNAMIC MONTE CARLO STUDY ON THE CORRELATION BETWEEN SHAPE AND SIZE OF LINEAR POLYMER CHAINS

The correlation between shape and size of linear chains on the simple cubic lattice is investigated using a dynamicMonte Carlo technique. A positive correlation between the asphericity parameter A and the square of the end-to-end distanceR~2, as well as that between A and the square of the radius of gyration S~2, is found for both RW and SAW chains, indicatingthat a chain conformation of small size is usually more spherical than one of large size. The result can explain why the shapeof the SAW chain deviate much more from a sphere than that of the RW chain, and can also explain the similar dependenceof size and shape on chain stiffness and on the distance of the first bead of a chain from an infinitely large flat surface.     

MONTE CARLO STUDY ON THE CRITICAL ADSORPTION POINT OF BOND-FLUCTUATED POLYMER CHAINS TETHERED ON ADSORBING SURFACES

 The behavior of three-dimensional bond fluctuation model chains tethered on an adsorbing flat surface was simulated by the Monte Carlo method. The dependence of the number of surface contacts M on the interaction strength e and the chain length N was investigated by a finite-size scaling law M = N[a₀ + a₁N^1/v_k + O((N ^1/v_k)²)] for e near the critical adsorption point e_c, i.e., k ≡(e-e_c)/e_c closes to 0. The critical adsorption point was estimated to be e_c = 0.93, and the exponents  = 0.49 and 1/v= 0.57.     

MONTE CARLO STUDY ON THE CRITICAL ADSORPTION POINT OF BONDFLUCTUATED POLYMER CHAINS TETHERED ON ADSORBING SURFACES

The behavior of three-dimensional bond fluctuation model chains tethered on an adsorbing flat surface was simulated by the Monte Carlo method.The dependence of the number of surface contacts M on the interaction strengthεand the chain length N was investigated by a finite-size scaling law M = Nφ[a0 + a1N1/V/κ+ O((N1/vκ)2)] forεnear the critical exponentsφ= 0.49 and 1/v= 0.57.     

Monte Carlo模拟方法研究管道受限对高分子链缠结的影响

采用基于格子模型的Monte Carlo模拟方法,考察了链长分别为30和250的线形高分子在管道内受限状态下的尺寸及扩散性质.结果表明,当线形高分子的链长低于其缠结链长时,其受限条件下高分子的性质随受限程度的变化而单调变化;而当线形高分子的链长大于缠结链长时,高分子的性质表现出非单调性,进而说明受限具有一定程度的解缠结...     

甲烷化反应机理的蒙特卡罗模拟

甲烷化反应机理的蒙特卡罗模拟郭向云，李永旺，王琴，罗庆云，钟炳（中国科学院山西煤炭化学研究所太原０３０００１）关键词蒙特卡罗方法，模拟，甲烷化机理催化反应机理可用普通动力学模型研究，也可用蒙特卡罗（ＭＣ）模型研究。对于一些比较复杂的催化反应，其动力学...