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

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

高分子链形状与尺寸关联的Monte Carlo模拟

运用MonteCarlo方法对线型高分子链格点模型的构型进行了模拟,研究了构型的尺寸(采用平方末端距R²,平方回转半径S²来表征)和形状(由非球形因子A表征)之间的关联.对任何长度的高分子链,其关联系数CA,R²和CA,S²均为正值,表明高分子链的形状与尺寸之间存在正关联,即尺寸小的构型其非球形因子A一般也小,反之尺寸大的构型其非球形因子A一般也大.关联系数CA,R²和CA,S²均随链长的增大而减小,近似地与链长的倒数(n^-1)成正比.研究还表明,关联系数的极限值(链长n很大时)与格点的类型无关,与链样本产生的方式也无关,但与链是否考虑排斥体积有关,考虑了排斥体积后,关联系数增大.     

Monte Carlo模拟研究高分子单链在基体中扩散的拓扑效应

高分子流体的性质同高分子的链结构及其动力学行为密切相关. 在高分子共混物中,共混组分的性质不仅依赖于自身的拓扑结构,还受到其它组分分子拓扑的影响. 本文中采用基于格子模型的Monte Carlo模拟方法研究了在高分子基体中扩散的4种不同拓扑组合（环形或线形高分子链）体系中目标单链的静态和动态性质. 结果发现,环形目标单链的性质受基体分子拓扑结构的影响要大于线形目标单链;其中环/线这一拓扑组合中目标单链的扩散机理相比其本体已经发生了较大改变,链末端在其中起了重要作用. 此外,我们对引起这一现象出现的可能原因做了分析.     

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

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

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

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

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

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

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

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

铝酸钠熔体Monte Carlo法计算机模拟研究

用计算机Monte Carlo 法研究铝酸钠熔体结构.计算了径向分布函数和局部结构.表明NaAlO_2熔体是由畸变的共隅AlO_4四面体离子集团、Na~+离子和xNa~+·yO~(2-)离子集团所组成.     

高分子单链热力学量的Monte Carlo模拟

提出了一种Monte Carlo模拟方法,可以计算高分子链的热力学量.该方法采用简单立方格子上的自避行走链和紧邻相互作用模型,计算了链长200以内的链在无量纲相互作用参数△_ε/(?)T不同数值下的构象熵和Helmholtz自由能,其中△_ε为一次接触的能量,(?)为Boltzmann常数,T为绝对温度.模拟结果和直接计数法的准确结果显示出较好的一致性 并且得到排除体积和相互作用效应互相抵消的(?)条件下的参数△_ε/(?)T之值大约为-0.6且与链长无关.     

铝酸钠熔体Monte Carlo法计算机模拟研究

用计算机Monte Carlo法研究铝酸钠熔体结构. 计算了径向分布函数和局部结构. 表明NaAlO_2熔体是由畸变的共隅AlO_4四面体离子集团、Na~+离子和xNa~+·yO~(2-)离子集团所组成.     

高分子动态Monte Carlo模拟的协同运动算法

在高分子动力学的研究中 ,动态 Monte Carlo模拟发挥了重要的作用 [1] .动态 Monte Carlo模拟的关键是选择具有物理真实性的高分子运动算法 .目前广为采用的算法是经 Hilhorst和 Deutch[2 ] 修正的 Verdier- Stockmayer算法 [3] 以及 Carmesin和 Kremer等的键长涨落算法 [4 ] ,陆建明和杨玉良曾提出一种高分子动态算法 [5] ,在他们的算法中很长一段链节可能作蛇行运动 ,但是冯捷等 [6 ] 指出这种运动模式不满足微观可逆性条件 .本文对该运动模式进行修正 ,得到一种协同运动算法 ,并对其动力学行为进行检验 .1　算　法在平面正方形格子…     

Monte Carlo Simulation of Thin Film Polymer Melts

We present Monte Carlo simulation data on conformations and dynamics of polymer melts confined in narrow slits of different widths and compare with data of bulk systems. We find that in confined geometries the chains swell laterally; they retain and even expand their spatially long-range correlations compared to bulk polymers and in contrast to the assumption of a complete screening of excluded volume. Long chains in bulk melts show entangled dynamics with a clear signature of a t^1/4-power law for the mean square displacements of innermost monomers at intermediate time scales. This behavior is gradually lost by confining the melts in slits with decreasing width. For ultra-thin films, the dynamics appears to follow a Rouse-like behavior over the entire subdiffusive regime. However, the terminal relaxation time is significantly increased compared to Rouse relaxation. This interesting observation was not reported previously and is the focus of our ongoing research.     

Monte Carlo Simulation of Coil-to-Globule Transition of Compact Polymer Chains: Role of Monomer Interacting

Coil-to-globule transitions are fundamental problems existing in polymer science for several decades; however, some features are still unclear, such as the effect of chain monomer interaction. Herein, we use Monte Carlo simulation to study the coil-to-globule transition of simple compact polymer chains. We first consider the finite-size effects for a given monomer interaction, where the short chain exhibits a one-step collapse while long chains demonstrate a two-step collapse, indicated by the specific heat. More interestingly, with the decrease of chain monomer interaction, the critical temperatures marked by the peaks of heat capacity shift to low values. A closer examination from the energy, mean-squared radius of gyration and shape factor also suggests the lower temperature of coil-to-globule transition.     

表面活性剂与高分子链混合体系的模拟

计算机模拟了高分子链对表面活性剂胶束形成过程的影响，以及高分子链构象性质随胶束化过程的变化.结果表明,当高分子链与表面活性剂之间的相互作用强度超过临界值后，高分子链的存在有利于表面活性剂胶束的形成.临界聚集浓度（CAC）与临界胶束浓度（CMC）的比值CAC/CMC随高分子链长的增大和相互吸引作用的增强而减小.在CAC之前，高分子链与表面活性剂分子只有动态的聚集；但在CAC之后，表面活性剂胶束随表面活性剂浓度X的增加而增大，并静态地吸附在高分子链上，形成表面活性剂/高分子聚集体.随着表面活性剂分子的加入，高分子链的均方末端距和平均非球形因子先保持恒定；从X略小于CAC开始， 和快速减小，至极小值后又逐渐增大.模拟结果支持高分子链包裹在胶束表面的实验模型.     

Monte Carlo simulation of the structure of mono- and bidisperse polyethylene nanocomposites

The structure of bidisperse polyethylene(PE) nanocomposite mixtures of 50:50(by mole) of long and short chains of C160H322/C80H162 and C160H322/C40H82 filled with spherical nanoparticles were investigated by a coarse-grained, on lattice Monte Carlo method using rotational isomeric state theory for short-range and Lennard-Jones for long-range energetic interactions. Simulations were performed to evaluate the effect of wall-to-wall distance between fillers(D), polymer-filler interaction(w) and polydispersity(number of short chains in the mixture) on the behavior of the long PE chains. The results indicate that long chain conformation statistics remain Gaussian regardless of the effects of confinement, interaction strength and polydispersity. The various long PE subchain structures(bridges, dangling ends, trains, and loops) are influenced strongly by confinement whereas monomer-filler interaction and polydispersity did not have any impact. In addition, the average number of subchain segments per filler in bidisperse PE nanocomposites decreased by about 50% compared to the nanocomposite system with monodisperse PE chains. The presence of short PE chains in the polymer matrix leads to a reduction of the repeat unit density of long PE chains at the interface suggesting that the interface is preferentially populated by short chains.     

Monte Carlo Simulation of Polymer Reactions at Interfaces

Adhesion of immiscible polymers during two‐component injection moulding may be improved by transreactions of properly functionalised components. We performed MC simulations based on the three‐dimensional coarse‐grained bond fluctuation model (BFM) including a thermal interaction potential in with energy to characterise the behaviour of several selected types of chemical reactions, which are governed by activation energies of E_A = 0, 1, 3 and 5 k_BT. The consumption of reactive monomers for all the reactions in the time interval below the Rouse time τ_R exhibits a typical crossover from a kinetic‐controlled to a diffusion‐controlled behaviour and can be described by a bimolecular kinetic ansatz.

     

Monte Carlo simulation of miscibility of polymer blends with repulsive interactions: effect of chain structure

The effects of the chain structure and the intramolecular interaction energy of an A/B copolymer on the miscibility of the binary blends of the copolymer and homopolymer C have been studied by means of a Monte Carlo simulation. In the system, the interactions between segments A, B and C are more repulsive than those between themselves. In order to study the effect of the chain structure of the A/B copolymer on the miscibility, the alternating, random and block copolymers were introduced in the simulations, respectively. The simulation results show that the miscibility of the binary blends strongly depends on the intramolecular interaction energy () between segments A and B within the A/B copolymers. The higher the repulsive interaction energy, the more miscible the A/B copolymer and homopolymer C are. For the diblock copolymer/homopolymer blends, they tend to form micro phase domains. However, the phase domains become so small that the blend can be considered as a homogeneous phase for the alternating copolymer/homopolymer blends. Furthermore, the investigation of the average end-to-end distance () in different systems indicates that the copolymer chains tend to coil with the decrease of whereas the of the homopolymer chains depends on the chain structure of the copolymers. As for the system containing the alternating or the random copolymers, the homopolymer chains also tend to coil with the decrease of . However, for the systems including the block copolymers, there is a slight difference in the of the homopolymer chains with the variation of .     

Monte Carlo simulation on layered polymeric films

Thin films of polymer blends composed of alternating copolymer, diblock copolymer and/or homopolymer are studied using Monte Carlo simulation. A multilayer morphology is observed in the film, that is, the blended polymers assemble into individual domains arranged from interior to the surfaces of the film. The coexisting components residing throughout the neighboring domains in the film make no distinguishable interface between any neighboring domains. By this means, it forms a vertical composition gradient in the polymeric film. Being different from layer-by-layer deposition of polyelectrolyte or hydrogen bonding approach etc., the layered structure in this study is formed by polymer blending in one step. Alternating copolymers are found to be essential components to form vertical composition gradient （layered structure） in thin films.     

高分子链超声裂解新的Monte Carlo算法

设计了一种新的Monte Carlo算法，在计算机上模拟了高分子链的超声裂解过程。首先根据Ovenall模型用遍历法对高分子链进行模拟断裂，随后进行倒接，进而反演了整个裂解过程。在获得大量模拟结果的基础上，分别研究了极限断裂次数的频率分布，最可几极限断裂次数与聚合度之间的关系，极限断片长度的平均分布，以及分子量多分散系数与裂解程度之间的关系。