表面反应动力学机理研究的新进展

介绍了近年来发展起来的研究表面反应动力学机理的新方法－－多和催化微观动力学分析，通过对表面化学键的形成和断裂的估计估算出基元反应的动力学参数，从而描述催化反应的特征。     

基于Monte Carlo模拟的化学反应动力学参数估算

提出并采用基于MonteCarlo模拟与动力学实验相结合的化学反应动力学参数估算方法,由基元反应确定MonteCarlo模拟具体做法,将MonteCarlo模拟结果与动力学实验结果相比较,根据比较结果自动调整和优化动力学参数,从而无需事先确定动力学方程即可有效估算各种化学反应的动力学参数值.采用该方法估算了丙烯氨氧化反应动力学参数,并对估算结果进行了分析与讨论.     

聚合动力学研究方法进展

从构造聚合动力学模型所采用的数学方法出发 ,对聚合动力学模型———概率统计模型、机理模型及基于MonteCarlo模拟模型等的进展进行了综述。     

扩散模型和凝聚模型耦合作用下胶体凝聚动力学的Monte Carlo模拟研究

以扩散模型(Ds(γ)=D0×sγ)和凝聚模型(Pij(σ)=P0×(i×j)σ)为基础,对胶体体系随时间的演变、团簇大小分布及其标度关系、团簇的重均大小S(t)的变化规律以及模型对最终分形维数的影响四个角度进行了比较研究,发现扩散指数γ0和凝聚概率指数σ0对胶体的凝聚动力学过程有相似的影响.本文在较宽的γ和σ取值范围内,对胶体的凝聚动力学进行了模拟研究,对慢速凝聚向快速凝聚的转化机理作了定量分析,并进一步分析了在团簇-团簇凝聚(CCA)模型下,得到类似扩散置限凝聚(DLA)模型的凝聚体的物理意义,结果表明:(1)γ0代表了体系中团簇或单粒做"定向运动"而非无规则的布朗运动的情况.这种"定向运动"的推动力可能来自于大团簇产生的强"长程范德华力"、"电场力"等,或来自于体系边界处的外力场的作用.(2)当σ0时,体系成为先快后慢的慢速凝聚,这可能对应大团簇为一排斥中心,即胶体颗粒存在"排斥力场"的现象.(3)证实了团簇的重均大小在凝聚过程的早期按指数规律增长,而后期按幂函数规律增长的实验现象.模拟研究还表明,胶体体系的凝聚动力学过程,在σ0时是一个存在正反馈机制的非线性动力学过程,而在σ0时则体现出负反馈的特征.     

用化学动力学方法估算颗粒表面的分维

A chemical kinetic method based on the kinetic expression of the second-order reaction on fractals is used to determine the spectral dimension of fractal surfaces. The method can conveniently give the spectral dimension by analyzing the relation between reactant concentration and reaction time (ct～t). The spectral dimension of the Diffusion-Limited Aggregate (DLA) is obtained by Monte Carlo simulation of the desorption process of chemisorbed hydrogen from DLA, and it is in agreement, with that in literature. Also, the spectral dimension of a real catalyst 12% Ni/Al2O3 is obtained by the method.     

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

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

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

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

甲烷部分氧化过程中强制振荡的动力学Monte Carlo模拟

为了有效控制动力学振荡行为以提高反应的转化率, 利用Monte Carlo 方法研究了甲烷催化部分氧化过程中的强制振荡行为, 探讨了原料气周期性变化对振荡动力学和转化率的影响. 研究表明原料气周期性变化不仅可以有效调控振荡的动力学行为, 产生如短周期振荡和双峰振荡等特殊动力学过程, 而且还可以提高反应的转化率. 当强制过程的周期从T/3 增大到2T(T为自发振荡过程的平均周期), 振荡过程从短周期振荡变化为双峰振荡. 对反应过程中CO的转化率进行了计算, 结果表明原料气周期性变化可以有效提高反应的转化率.振荡动力学的改变和转化率的提高主要是因为强制振荡过程使得催化剂表面发生了从氧化态向还原态的转变.     

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

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

受限高分子薄膜界面动力学与缠结的关系

通过Monte Carlo模拟结合原始路径分析(PPA)的方法,阐述了缠结高分子薄膜的界面动力学与缠结程度的关系.研究发现,以吸附势能的临界值εwc≈-0.6 kBT附近为界,当墙壁-高分子作用势能从弱吸引到强排斥变化时,界面层中的链移动快于中心层,只有当墙壁的吸引作用增强到一定程度时,界面层中的链移动才会慢于中心层.界面动力学受到促进或阻碍可能与界面层和中心层的缠结程度直接相关:界面层的缠结程度保持在本体水平上基本不变;中心层的缠结程度在强吸引表面上低于界面层,而在弱吸引和排斥表面上高于界面层.此外,中心层和界面层中高分子链受限程度的变化对薄膜界面动力学行为的转变产生一定影响.对于薄膜中链密度分布情况随墙壁-高分子作用势能变化的分析为相关的物理化学机制提供了理论依据.     

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

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

碘离子催化的H2O2-BAT反应的Monte Carlo模拟

催化动力学分析是痕量分析中检测下限低、选择性较好的一种方法,但目前分析工作者对所用的反应机理常不明瞭,仅能作为“黑箱”应用之。模拟技术可为探讨这类分析反应机理提供有用的信息。前文^[1]报道用溴胺T(BAT)离子电极、跟踪碘离子催化的BAT-H₂O₂反应测定痕量碘效果甚佳。     

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

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

Coarse-graining Calcium Dynamics on Stochastic Reaction-diffusion Lattice Model

We develop a coarse grained (CG) approach for efficiently simulating calcium dynamics in the endoplasmic reticulum membrane based on a fine stochastic lattice gas model. By grouping neighboring microscopic sites together into CG cells and deriving CG reaction rates using local mean field approximation, we perform CG kinetic Monte Carlo (kMC) simulations and find the results of CG-kMC simulations are in excellent agreement with that of the microscopic ones. Strikingly, there is an appropriate range of coarse proportion m, corresponding to the minimal deviation of the phase transition point compared to the microscopic one. For fixed m, the critical point increases monotonously as the system size increases, especially, there exists scaling law between the deviations of the phase transition point and the system size. Moreover, the CG approach provides significantly faster Monte Carlo simulations which are easy to implement and are directly related to the microscopics, so that one can study the system size effects at the cost of reasonable computational time.     

CSTR-BZ复杂化学反应的Monte-Carlo模拟

在Monte-Carlo算法中引入流率项,通过对具有流率分叉结构的CSTR-BZ反应体系SNB(Showalter-Noyes-Bar-Eli)模型的随机分析,讨论了最低反应浓度的意义及流率项连续性的影响.与实验结果及确定性数值积分结果比较表明,对于流动的实际反应体系,基于“质量作用定律”的宏观确定性方法仍具有适用性.     

水煤气变换反应初始微观反应动力学的模拟研究

用蒙特水罗方法对水气转移反应两种典型机理的寝反应动力学进行了研究。结果表明,Ｔｅｍｋｉｎ等人的氧化－还原机理中,无论反应物ＣＯ和Ｈ２Ｏ的相对浓度怎样变化,ＣＯ２和Ｈ２的初始生成速率总是相差很大。     

Calculation of Surface Excitation Parameters by a Monte Carlo Method

Electron inelastic mean free path (IMFP) is an important parameter for surface chemical quantification by surface electron spectroscopy techniques. It can be obtained from analysis of elastic peak electron spectroscopy (EPES) spectra measured on samples and a Monte Carlo simulation method. To obtain IMFP parameters with high accuracy, the surface excitation effect on the measured EPES spectra has to be quantified as a surface excitation parameter (SEP), which can be calculated via a dielectric response theory. However, such calculated SEP does not include influence of elastic scattering of electrons inside samples during their incidence and emission processes, which should not be neglected simply in determining IMFP by an EPES method. In this work a Monte Carlo simulation method is employed to determine surface excitation parameter by taking account of the elastic scattering effect. The simulated SEPs for different primary energies are found to be in good agreement with the experiments particularly for larger incident or emission angles above 60° where the elastic scattering effect plays a more important role than those in smaller incident or emission angles. Based on these new SEPs, the IMFP measurement by EPES technique can provide more accurate data.     

Reaction kinetics of nonideal hyperbranched polymerizations: Influences of chain rigidity and reaction reversibility

To deeply understand the dependence of hyperbranched polymerization on the specialties of molecular chain and chemical reaction, the reaction kinetics of nonideal hyperbranched polymerizations considering the chain rigidity and reaction reversibility are studied using the reactive three‐dimensional bond fluctuation lattice model. It is found that, with the increase of chain rigidity, the formation probability of intramolecular rings decreases and less intramolecular rings with larger size are formed. It results in the increase of the degrees of polymerization and polydispersity index with the rise of chain rigidity at higher conversion. Furthermore, our simulation shows that the reversibility of reaction has strong influences on the equilibrium state and kinetic process of hyperbranched polymerizations. When the ratio of reverse reaction probability to forward reaction probability is larger, the conversion and degrees of polymerization quickly grow to the equilibrium value simultaneously. At smaller reverse reaction probability, however, the weigh‐average degree of polymerization and polydispersity index further increase slowly after approaching the equilibrium conversion. Our results are well‐consistent with the experiments to indicate that the introduction of the specialties of macromolecule and chemical reaction are necessary for quantitative analysis of the realistic hyperbranched polymerizations. © 2012 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem, 2012     

Dynamics of Bond‐Fluctuation Model Chains in Good and Theta Solvents

The dynamic Monte Carlo algorithm is employed to explore the dynamics of flexible linear chains. The chains are represented by the bond‐fluctuation model with and without attractions between non‐bonded units placed at close distances. This mimics the behavior of real chains in the good and poorer solvents. We obtain the chain sizes, diffusion coefficients, Rouse modes, and their relaxation times. We also evaluate the time correlation function of the end‐to‐end vector at different concentrations. Subsequently, we compare the dependence of the simulation results on chain length, solvent quality, concentration, and mode order with the corresponding theoretical predictions. We observe a retardation of diffusion for non‐dilute systems close to the theta state. This retardation is too high to be exclusively attributed to the increase of global friction and can be caused by temporary adherence of the chains to transient clusters.