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

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

D_4/N-β-氨乙基-γ-氨丙基甲基二甲氧基硅烷本体共聚反应聚合度变化的Monte Carlo模拟

用MonteCarlo方法模拟八甲基环四硅氧烷 (D4 )与N β 氨乙基 γ 氨丙基甲基二甲氧基硅烷(APAEDMS)的本体开环共聚动力学 .模拟过程采用自由体积理论简化处理扩散效应并与本征反应动力学耦合 .本征动力学常数通过模拟主要共聚基元反应得到 ,基于优化的动力学常数通过模拟从分子水平揭示D4 APAEDMS本体开环共聚反应过程主要表现为活性阴离子聚合行为 ,同时又伴有部分逐步聚合特征 .     

平行反应的热动力学研究法

根据热动力学基本理论, 推导了平行准一级反应和平行准一二级反应的热动力学方程, 建立了热谱解析平行反应动力学参数的热动力学研究法, 并利用该法研究了两个模拟平行反应体系的热动力学, 实验结果验证了方法的正确性.     

工业Fe-Mn催化剂上基于详细反应机理的F-T合成动力学模型Ⅰ.烯烃再吸附动力学研究

 利用转篮式无梯度反应器,在工业Fe-Mn催化剂上,在较宽的工业操作相关的反应条件下进行了F-T合成反应动力学研究. 首次提出了利用转篮式无梯度反应器反应气氛和反应温度均一的优势,将烃生成反应动力学的估算从传统的对烃生成和水煤气变换这两类发生在不同活性中心的反应同时进行估算的方法中分离出来,简化了烃生成动力学模型的计算. 在基于亚甲基插入的亚烷基机理动力学模型基础上,考虑到乙烯与催化剂表面强的相互作用,将乙烯和乙烷的生成动力学参数单独计算. 动力学模型计算的链增长、烷烃和烯烃生成的活化能均与文献报道值具有较好的一致性. 由F-T合成动力学模型计算的合成气消耗速率、甲烷生成速率和C5+的生成速率较好地与实验值吻合. 通过动力学模型并结合实验结果分析发现,未考虑除化学反应之外的非本征因素的烯烃再吸附动力学模型不能够正确预测烃产物分布偏离ASF规律及烯烷比随碳数增加而下降的现象.     

硅橡胶硫化反应场的数值模拟

以化学反应动力学原理为基础,实验确定了硅氢加成反应的动力学模型;引入 了硫化率增量的概念,获得了全量硫化率和硫化率增量的数值计算式;根据统计理 论,获得了起始分子无规分布时的交联结构参数的数值计算式;介绍了非稳态硫化 反应场的有限元模拟步骤。在此基础上,设计了硅橡胶硫化反应过程的有限元模拟 软件;应用该软件,根据输入的初边值条件,可计算任一空间点在任一时刻的交联 反应参数和交联结构参数,进而可预测制品性能、优化反应参数,或优化设计高分 子材料。最后给出了一个典型算例,并实验验证了该模拟理论及算法的合理性。     

Kalman滤波用于酶反应动力学的研究

采用卡尔曼滤波研究酶反应动力学及分析，依据酶存在下过氧化氢放出氧的反应建立了扩展Ｋａｌｍａｎ滤波模型，用滤波处理确定了反应速率与起始速率浓度的关系，并估算出酶反应动力学参数Ｋｍ。     

高分子浓溶液相分离动力学的Monte Carlo模拟

运用链动力学MonteCarlo方法模拟了高分子浓溶液的spinodal相分离动力学过程 .结果表明 :分相早期 ,散射峰的位置向左移动 ,不符合经典的Cahn Hilliard线性理论 ;分相后期 ,聚合物相形态处于逾渗状态 ,结构因子可以被标度化 ,且基本符合相应的Furukawa标度律 .模拟结果与相关实验相符 .揭示了线团尺寸在深度淬冷过程中先急剧收缩、然后随相分离的进行而逐步上升的复杂变化 ,表明链动力学MonteCarlo方法能同时考察高分子链的构象变化和多链体系的分相过程以及两者的关联     

绝热法研究己内酰胺阴离子聚合尼龙动力学

采用己内酰胺钠盐、N-75缩二脲作为反应催化体系, 确定反应温度在145-160 ℃之间, 通过计算得到动力学参数: 反应级数为准一级、活化能在73.2-77.1 kJ·mol-1之间、指前因子在2.9×1011-3.6×1011 mol1-n·s-1范围内. 本实验条件下测定并计算的反应热为134.5-137.3 J·g-1, 与文献值(138.6 J·g-1)吻合. 并在前人基础上修正并建构了己内酰胺阴离子绝热反应动力学模型, 对反应过程的模拟结果与实验数据基本吻合, 从而证明了本模型的正确合理.     

不同金属催化水煤气变换反应活性的Monte Carlo模拟研究

运用BOC-MP方法对Cu(110),Cu(111),Pd(111)和Au(111)等过渡金属催化的WGS反应的可能微观动力学步骤进行了详尽的能学数据计算,并结合MonteCarlo方法对WGS反应的表面氧化还原机理进行了计算机模拟。结果表明,Cu的催化活性优于Pd,Au的催化活性,并获得了相应金属上WGS反应的表观活化能及动力学指前因子(相对值);在此基础上,对该反应的结构敏感性进行了研究,发现该反应为一结构敏感反应,与实验结果相符。     

合成甲醇Cu基催化剂结构敏感性的模拟

运用BOC-MP方法对Cu（111）、Cu（110）上CO2＋H2合成甲醇的甲酸基机理中的可能基元反应步骤的动力学参数进行了计算．在此基础上，结合MonteCarlo模拟技术对该反应过程进行了随机模拟．模拟结果表明，该反应在钢的两种不同单晶表面上的反应活化能值非常接近（其比值为1：1.2，基本吻合于实验值），表明该反应活性与催化剂表面原子排布的关系不甚紧密，即该反应为一非结构敏感反应．     

O—酰基—α—酮肟光分解反应的Monte Carlo处理

本文首次对光化学反应体系用Monte Carlo方法进行模拟处理。通过5个O-酰基-α-酮肟光分解反应的Monte Carlo模拟,可避免解析解中由于对吸收光强须采用一级近似求解动力学微分方程组,而造成拟合反应在后期产生与实验结果的偏差。     

苯乙烯在顺1,4-聚丁二烯上的接枝共聚反应动力学的Monte Carlo模拟

用ＭｏｎｔｅＣａｒｌｏ模拟方法对苯乙烯在顺１，４－聚丁二烯上的接枝反应动力学进行了研究。模拟结果表明，对于不同的橡胶浓度，接枝的和均聚的聚苯乙烯的数均聚合度、接枝效率与苯乙烯在顺１，４－聚丁二烯上的接枝共聚反应的动力学的解析解十分吻合，证明本方法能够有效地应用于自由基型接枝共聚合反应体系。     

Permeation and gating in proteins: kinetic Monte Carlo reaction path following

We present a new Monte Carlo technique, kinetic Monte Carlo reaction path following (kMCRPF), for the computer simulation of permeation and large-scale gating transitions in protein channels. It combines ideas from Metropolis Monte Carlo (MMC) and kinetic Monte Carlo (kMC) algorithms, and is particularly suitable when a reaction coordinate is well defined. Evolution of transition proceeds on the reaction coordinate by small jumps (kMC technique) toward the nearest lowest-energy uphill or downhill states, with the jumps thermally activated (constrained MMC). This approach permits navigation among potential minima on an energy surface, finding the minimum-energy paths and determining their associated free-energy profiles. The methodological and algorithmic strategies underlying the kMCRPF method are described. We have tested it using an analytical model and applied it to study permeation through the curvilinear ClC chloride and aquaporin pores and to gating in the gramicidin A channel. These studies of permeation and gating in real proteins provide extensive procedural tests of the method.     

Incorporation of inhomogeneous ion diffusion coefficients into kinetic lattice grand canonical monte carlo simulations and application to ion current calculations in a simple model ion channel

To deal with inhomogeneous diffusion coefficients of ions without altering the lattice spacing in the kinetic lattice grand canonical Monte Carlo (KLGCMC) simulation, an algorithm that incorporates diffusion coefficient variation into move probabilities is proposed and implemented into KLGCMC calculations. Using this algorithm, the KLGCMC simulation method is applied to the calculation of ion currents in a simple model ion channel system. Comparisons of ion currents and ion concentrations from these simulations with Poisson-Nernst-Planck (PNP) results show good agreement between the two methods for parameters where the latter method is expected to be accurate.     

Monte Carlo estimation of kinetic parameters in polymerization reactions

A general method for estimating kinetic parameters in polymerization reactions using Monte Carlo simulation to represent the models of the reactions is developed. From a statistical point of view, the procedure is a Bayesian one in which a posterior probability density surface (PPDS) is calculated for points on a grid in the parameter space. A smoothing function is fitted to the PPDS, then a posterior probability region, which is similar to a confidence region, is calculated for the parameters. An application to a relatively trivial example, the Mayo–Lewis copolymerization model is shown in detail. Many other potential applications are suggested.     

Quantum effect induced kinetic molecular sieving of hydrogen and deuterium in microporous materials

We report here our investigations using Monte Carlo and molecular dynamics (MD) simulations, as well as quasi-elastic neutron scattering experiments, to study the adsorption and diffusion of H₂ and D₂ in zeolite Rho. In the simulations, quantum effects are incorporated via the Feynman-Hibbs variational approach. At low temperatures, we observe a reversal of kinetic molecular sieving in which D₂ diffuses faster than H₂. Based on fits of bulk data, we suggest new set of potential parameters for hydrogen, with the Feynman-Hibbs variational approach used for quantum corrections. The transport properties obtained from MD simulations are in excellent agreement with the experimental results, with both showing significant quantum effects on the transport at low temperature. The MD simulation results on two different structures of zeolite Rho clearly demonstrate that the quantum effect is very sensitive to pore size. High transport flux selectivity is noted at low temperatures, suggesting feasibility of kinetic isotope separation.     

Molecule-based kinetic modeling by Monte Carlo methods for heavy petroleum conversion

A methodology for kinetic modeling of conversion processes is presented.The proposed approach allows to overcome the lack of molecular detail of the petroleum fractions and to simulate the reactions of the process by means of a two-step procedure.In the first step,a synthetic mixture of molecules representing the feedstock is generated via a molecular reconstruction method,termed SR-REM molecular reconstruction.In the second step,a kinetic Monte Carlo method,termed stochastic simulation algorithm(SSA),is used to simulate the effect of the conversion reactions on the mixture of molecules.The resulting methodology is applied to the Athabasca vacuum residue hydrocracking.An adequate molecular representation of the vacuum residue is obtained using the SR-REM algorithm.The reaction simulations present a good agreement with the laboratory data for Athabasca vacuum residue conversion.In addition,the proposed methodology provides the molecular detail of the vacuum residue conversion throughout the reactions simulations.     

Monte Carlo source simulation technique for solution of interference reactions in INAA experiments: a preliminary report

A new method using Monte Carlo source simulation of interference reactions in neutron activation analysis experiments has been developed. The neutron spectrum at the sample location has been simulated using the Monte Carlo code MCNP and the contributions of different elements to produce a specified gamma line have been determined. The produced response matrix has been used to measure peak areas and the sample masses of the elements of interest. A number of benchmark experiments have been performed and the calculated results verified against known values. The good agreement obtained between the calculated and known values suggests that this technique may be useful for the elimination of interference reactions in neutron activation analysis.