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

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

Spectroscopic Investigation of Asphaltene Aggregation in Carbon Tetrachloride Medium Containing Microquantities of Water

By fluorescence spectroscopy it has been established that coal derived asphaltene (CDA) forms aggregates in CCl₄ medium. The new feature of this compound is that even before aggregate formation it can incorporate water molecules and the critical aggregation concentration (CAC) in nonaqueous medium increases with increase in ω (= water:CDA mole ratio). The fluorescence spectroscopic determination of the equilibrium constant (K) of the molecular complex of [60] fullerenes with CDA in carbon tetrachloride medium has revealed that K increases on addition of microquantities of water. This observation has been rationalized by assuming that the nitrogenatoms at one end of CDA molecules attach themselves to the water micropool and thus producing a different conformation with a wider distance at the other end. This pre-organization further facilitates inclusion. By utilizing fluorescence variation of [60]fullerene/CDA complex in CCl₄ as continuous medium with change in ω a critical ω (=5) has been found. The formation of CDA/[60]fullerene complex has also been established by NMR spectroscopy.     

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

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

Monte Carlo Simulation in Electron Probe Microanalysis. Comparison of Different Simulation Algorithms

Practical aspects of Monte Carlo simulation of EPMA experiments are considered. Simulations are performed using the general-purpose Monte Carlo code system Penelope, which is briefly described. This code includes geometry tools and variance reduction methods that allow the practical simulation of x-ray spectra from samples with complex geometries in moderate computing times. The reliability of simple interaction models and approximations, which have been frequently used in EPMA studies, is analyzed by studying their effects on the simulated x-ray spectra.     

胶粒分形粒子簇反应控制聚集动态行为位垒影响的Monte Carlo模拟

通过对反应控制聚集过程的MonteCarlo模拟,从微观及介观层次上探讨了胶粒间相互作用位能曲线上位垒高度的变化对胶粒分形粒子簇大小分布和动态标度函数及聚集动力学行为的影响规律。     

布朗扩散系数对胶体分形粒子簇扩散控制聚集动力学影响的Monte Carlo模拟

使用表征粒子簇结构的几何形状因子,通过对扩散控制聚集过程的模拟,从微观或介观层次研究了粒子簇结构对粒子簇增长速率和速率常数的影响规律,并与实验结果进行了对比分析.     

自由基聚合反应的Monte Carlo模拟方法

本文将化学反应动力学的MonteCarlo模拟方法运用到引发剂引发的自由基聚合反应的非稳态动力学,针对自由基聚合反应动力学数值模拟所特有的"无伸缩问题",采用"偏倚抽样法"解决了MonteCarlo模拟中的"无伸缩问题",模拟结果与非稳态动力学解的结果完全一致,此算法易推广到研究更复杂的自由基聚合反应体系。     

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

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

Comprehensive Study of Free Radical Copolymerization Using a Monte Carlo Simulation Method, 1

Summary: In order to investigate the influence of reactivity ratios and initial feed composition on the microstructure of macromolecules in free radical copolymerization, a comprehensive study was carried out using a Monte Carlo simulation method. As a result, a new procedure was introduced to modify the works of others on the initiation step. The variation of the copolymer composition and the fashion of the arrangement of monomers in simulated chains were evaluated as a function of copolymerization parameters. The model was capable of monitoring any change in azeotropy as well as the magnitude and direction of composition drift from the azeotrope point. The maximum reachable conversion (MRC) was predicted for different combinations of initial feed compositions and reactivity ratios. According to the simulation results, a critical conversion where the macromolecules produced inherited the maximum allowed alternation was obtained for the reactivity ratios given.

Change of sequence distribution of simulated copolymer chains with conversion for various initial feed compositions on a triangular graph (r_A = 0.5, r_B = 0.9).     

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

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

分形表面上CO氧化反应的蒙特卡罗模拟

众所周知,人们在理想完好的单晶表面上的动力学研究已取得了很大的成功[1],然而对许多实际催化剂来说,传统催化理论和实验之间不符的情况仍然俯拾皆是．现在人们对造成这种情况的原因已经有了更加深入的认识：大量的实验事实表明,在表面催化反应体系中,不但催化剂表面具有复杂的分形结构,而且催化剂表面上的活性中心分布也具有复杂的分形特性[2]．“分形”是指那些具有分数维数的几何对象．这些对象往往是不规则的,不能用通常的欧氏几何来描述[3]．将分形引入催化科学中最早的是Pfeifer和Avnir等[4]人,至今已有近二十年的历史,现…     

Monte Carlo Simulation of Ethylene Copolymers: A Look into the Monomer Sequence Distribution

A comprehensive mathematical model was developed using Monte Carlo simulation to describe the mechanism of ethylene and α-olefin copolymerization. The model studies the polymerization mechanism using coordination catalysts and is able to predict molecular weight and detailed chemical composition distributions. This work is considered to be a useful tool that enables us to understand and described the monomer sequence distribution as a function of chain length in semi-batch polymerization reactors.     

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.

     

Carboxyl Terminated Polymer Chain Extension Using a Bisoxazoline Coupling Agent: Monte Carlo Simulation

Summary: The kinetics of the chain extension reaction of a carboxyl‐terminated polymer using a bisoxazoline coupling agent was simulated by the Monte Carlo method based on the master equation. The effects of temperature and stepwise addition of the chain extender was examined. A comparison between simulated results and those calculated by an amended kinetic model was made. The results show that the highest coupling efficiency and the highest can be obtained when the initial concentrations of carboxyl and oxazoline groups are equal, which is in good agreement with the experiments. It is found that a higher reaction temperature could lead to a bigger coupling efficiency, a higher and narrower MWDs. The stepwise addition of the chain extender can only postpone the chain extension reaction, but cannot affect the final coupling efficiency and the MWDs when the concentration of the oxazoline group is lower than that of the carboxyl group. However, stepwise addition of the chain extender favors bigger coupling efficiency and narrower MWDs when the concentration of the oxazoline group exceeds that of the carboxyl group.

The relationship among three group of polymers, r and ratios of k₂ to k₁. The dots are the data simulated by MC and solid lines are relative data calculated by the kinetic model.     

组合取样的误差理论研究及其Monte Carlo模拟

基于多项分布理论,建立了组合取样方差与取样总体中各子区的大小及其组分含量之间的关系和计算公式,探讨了组合取样常数的物理意义.以颗粒物质为例,探讨了组合取样的逻辑质量单元的概念及其意义,为确定组合取样中的份样质量提供了理论依据.本文对于完善组合取样误差理论,保证组合取样的质量均具有重要意义.     

Monte Carlo模拟表面上吸附态物种的迁移对反应行为的影响

用MonteCarlo方法研究了非均相催化剂表面吸附态氢原子的迁移对催化反应活性的影响,模拟结果表明,吸附态氢原子扩散很慢时,表面活性位很快被氢原子饱和,转换频率TOF增大到一定程度时很快下降;而当表面吸附态氢原子的扩散速率达到足够大的程度时,TOF将不再受氢原子扩散的影响.     

凝并和成核机理下颗粒尺度分布的Monte Carlo求解

颗粒的凝并和成核现象影响其尺度分布,现有的MonteCarlo方法描述颗粒尺度分布的时间演变过程存在若干困难.提出了一种新的多重MonteCarlo(MMC)算法,基于时间驱动,利用加权的虚拟颗粒的思想,在模拟过程中保持虚拟颗粒总数不变和计算区域体积不变.利用该算法对“常凝并核,一阶成核”的情况下颗粒尺度分布的时间演变过程进行了数值求解,所得结果与数值解相符,表明MMC算法具有高且稳定的计算精度.另外,MMC算法由于跟踪比实际颗粒数目少得多的虚拟颗粒而具有较低的计算代价.     

蒙特卡洛模拟在分子模拟教学中的应用*

以有机化合物的模板诱导沉积实验为例,介绍了蒙特卡洛模拟在分子模拟教学中的应用。本实验在有机化合物诱导沉积实验的基础上,将诱导沉积实验、蒙特卡洛模拟方法系统综合在一个实验教学中,让学生在学习蒙特卡洛模拟方法的同时,更直观、形象地理解分子间相互作用对有机分子聚集形貌的影响,加深学生对有机分子自组装概念的理解。通过实际上机操作,学生还能初步了解Dev-C++,RasMol等相关软件的使用,并加深学生对结构化学及计算化学课程相关知识的理解。     

Gibbs系综Monte Carlo模拟甲烷的吸附平衡

在263、298和313 K下,对甲烷在1.91 nm的活性炭孔中的吸附平衡进行了Gibbs系综Monte Carlo(GEMC)模拟的研究.改进了GEMC方法,使之可用于模拟指定压力下的吸附平衡.通过改进的GEMC模拟,得到了在1.91 nm的活性炭中甲烷在263、298和313 K时的吸附等温线;发现263 K时的超额吸附量要大于298 K、313 K时的超额吸附量; 且不同温度下的超额吸附等温线均存在一最大超额吸附.263 K时,超额吸附量在5.0 MPa时出现最大值;而298 K、313 K时超额吸附量则在7.0 MPa时出现最大值.此工作为不同温度下天然气吸附存贮过程的开发及设计提供了依据.     

巨正则系综Monte Carlo模拟方法确定活性炭的微孔尺寸

根据299K下甲烷在活性炭中的吸附实验数据,通过调节狭缝微孔的孔宽参数,利用巨正则系综MonteCarlo(GCEMC)方法得到不同孔宽下流体的微观结构以及吸附等温线.比较并拟合模拟结果和实验数据,确定了活性炭微孔的平均孔宽,为下一步求解微孔尺寸分布以及为预测吸附剂在不同温度下吸附不同吸附质分子时的吸附性能提供了基础与指导.模拟中,甲烷分子采用单点Lennard-Jones球型分子模型,活性炭用狭缝孔来近似表征,流体分子与单个狭缝墙的相互作用采用著名的Steele的10－4-3势能模型.模拟表明,此方法为考察介孔材料的微孔分布以及微孔平均孔宽提供了新的思路.