水煤气变换反应的微观动力学分析──Cu基催化剂上H_2S中毒失活原因的 Monte Carlo 模拟研究

用ＢＯＣＭＰ方法及ＭｏｎｔｅＣａｒｌｏ模拟技术，对Ｈ２Ｓ导致Ｃｕ基催化剂失活的原因进行了计算分析。研究结果表明，当原料气中存在Ｈ２Ｓ时，ＷＧＳ反应的活化能明显高于无Ｈ２Ｓ时的活化能，随着表面Ｈ２Ｓ浓度的增大（θ＝０，０．１０，０．２５），反应的活化能也逐渐变大（其大小比为１∶１．３４∶３．３），究其原因可归结为Ｈ２Ｓ的存在使得反应物分子的吸附热减小，从而使Ｈ２Ｏ的解离吸附（ＷＧＳ反应的速控步骤）活化能增大。     

Using molecular dynamics simulations to identify the key factors responsible for chiral recognition by an amino acid-based molecular micelle

Molecular dynamics (MD) simulations were used to investigate the binding of six chiral compounds to the amino acid-based molecular micelle (MM) poly-(sodium undecyl-(L)-leucine-leucine) or poly(SULL). The MM investigated is used as a chiral selector in capillary electrophoresis. The project goal was to characterize the chiral recognition mechanism in these separations and to move toward predictive models to identify the best amino acid-based MM for a given separation. Poly(SULL) was found to contain six binding sites into which chiral compounds could insert. Four sites had similar sizes, shapes, and electrostatic properties. Enantiomers of alprenolol, propranolol, 1,1′-bi-2-naphthyl-2,2′-diyl hydrogen phosphate, 1,1′-bi-2-naphthol, chlorthalidone, or lorazepam were separately docked into each binding pocket and MD simulations with the resulting intermolecular complexes were performed. Solvent-accessible surface area calculations showed the compounds preferentially associated with binding sites where they penetrated into the MM core and shielded their non-polar atoms from solvent. Furthermore, with five of the six compounds the enantiomer with the most favorable free energy of MM association also experienced the most favorable intermolecular hydrogen bonding interactions with the MM. This result suggests that stereoselective intermolecular hydrogen bonds play an important role in chiral discrimination in separations using amino acid-based MMs.GRAPHICAL ABSTRACT     

Numerical Study of Pore Sizes Distribution in Gels

We introduce a new numerical technique for the calculation of the pore size distribution in two-dimensional disordered systems. Our method is based on a triangulation technique which allows a closer measurement of pores surface without any morphological hypothesis.In this work, we focus our calculations in simulated gels. Such materials are modeled in two different conditions: by means of the Diffusion-Limited and Reaction-Limited Cluster-cluster Aggregation algorithms, DLCA and RLCA, respectively. In both situations, when the particles concentration decreases, the average pores size increases. The more compact cluster in RLCA, compared with DLCA, is consistent with the pore size distribution we have calculated. The simulated mean pore size is quantitatively in agreement with experimental data from literature.     

ZnCl2熔盐的分子动力学模拟

近年来，若干作者根据几种粒子间函数，对ZnCl2熔盐结构做过分子动力学模拟［1－3］其出的Zn－Cl和Cl－Cl离子间的偏径向分布函数与中子衍射实测值符合较好，但Zn－Zn离子间距以及Zn－Zn间配位数计算值多偏高．鉴于Busing势函数在多价卤化物馆盐的分子动力学计算中应用效果较好[4],我们试用Busing势函数为基础对ZnCl2熔盐结构和能量做分子动力学计算．1研究方法计算所用粒子势为Busing势函数此处，Zi为离子的电行数（ZZn。＋＝2，Zcl－＝1），几；是离子有效半径，人为＊离子的“硬度”参数·据文献问，f二0．005071，尸zn。十二0．…     

Polymer crystallization dynamics,as reflected by differential scanning calorimetry. Part 2: Numerical simulations

With the aid of a model for the kinetics of polymer crystallization, as put forward in previous publications, the shape of DSC-curves and their position on the temperature scale were simulated for various conditions of heat transfer in the apparatus. It turns out that the outcome is very dependent on the assumptions made with respect to these heat transfer conditions. For the ideal condition — no temperature differences between sample, pan and furnace — an invariable shape is predicted for the DSC-curves. They only shift to lower temperatures with increasing cooling rates. For more realistic conditions, the curves not only shift but become broader and their maxima decrease. They show a more familiar appearance. These calculations are very involved, however, A simple balance equation is shown to yield equivalent results, if a dimensionless characteristic number like the Nusselt number remains considerably smaller than one. This number contains an effective heat transfer coefficient between sample and furnace which, surprisingly, should not be too high. Apparently, the heat capacity of the pan does not play an important role under these conditions. This is investigated in Appendix II. Appendix I describes the procedure of the numerical simulations.     

Effect of enzyme concentration on the dynamic behavior of a membrane-bound enzyme system

The dynamic behavior of the reaction-diffusion system, composed of glucose oxidase (EC 1.1.3.4) immobilized at a uniform concentration in a membrane, used as a glucose electrode is represented by a diffusion equation with a nonlinear reaction-term in one-dimensional space. The mathematical model is analyzed by computer simulation, that is, numerical integration of the equation under various initial and boundary conditions, to examine the effect of enzyme concentration on the response characteristics (responsiveness and linearity in response) of the electrode. The analysis of the responses of the system to stepwise changes in the boundary value (glucose concentration in simple solution) infers that the enzyme concentration governs the patterns of the spatial distributions of the substrates (glucose and dissolved oxygen) in steady states and transient responses. It is also revealed that the response characteristics of the electrode are optimized with concentration of immobilized enzyme and that the system establishes the steady states at the same spatial distributions of the substrates, regardless of the boundary value. The diffusion of the substrates and the oxygen concentration also have significant effects on the response characteristics of the electrode.     

金属铜升温熔化过程的分子动力学模拟

采用分子动力学方法模拟了金属铜的升温熔化过程.原子间作用势采用FS (Finnis-Sinclair)势，结构分析采用双体分布函数(PCF)、均方位移(MSD)等方法.计算结果表明,在连续升温过程中，金属铜在1444 K熔化，在该熔化点的扩散系数为4.31×10－9 m2•s－1.上述结论与实验值相当接近，并且比之采用EAM镶嵌原子势所作模拟得到的结果更佳，说明FS势可以用来处理象液铜这样较复杂的无序体系.本文指出了升温速率在金属熔化过程中所起的作用.     

Dynamics of the formation of two-dimensional ordered structures

The growth of ordered domains in lattice gas models, which occurs after the system is quenched from infinite temperature to a state below the critical temperatureT _c, is studied by Monte Carlo simulation. For a square lattice with repulsion between nearest and next-nearest neighbors, which in equilibrium exhibits fourfold degenerate (2×1) superstructures, the time-dependent energy E(t), domain size L(t), and structure functionS(q, t) are obtained, both for Glauber dynamics (no conservation law) and the case with conserved density (Kawasaki dynamics). At late times the energy excess and halfwidth of the structure factor decrease proportional tot ^–x, whileL(t) t ^x, where the exponent x=1/2 for Glauber dynamics and x1/3 for Kawasaki dynamics. In addition, the structure factor satisfies a scaling lawS(k,t)=t ^2xS(kt^x). The smaller exponent for the conserved density case is traced back to the excess density contained in the walls between ordered domains which must be redistributed during growth. Quenches toT>T _c, T=T_c (where we estimate dynamic critical exponents) andT=0 are also considered. In the latter case, the system becomes frozen in a glasslike domain pattern far from equilibrium when using Kawasaki dynamics. The generalization of our results to other lattices and structures also is briefly discussed.