聚2-吡咯烷酮的非等温结晶动力学

用差示量热扫描热分析仪(DSC)测试了不同降温速率下聚2-吡咯烷酮(PPD)样品的温度-热焓曲线,样品黏均分子量为2.2×10~4,熔点为272℃。采用Jeziorny法、Ozawa法和莫志深法分析了PPD的非等温结晶动力学。结果表明,在给定降温速率范围内,Ozawa法不适用于描述PPD的非等温结晶动力学过程,Jeziorny法只适用于描述PPD的主结晶阶段,而莫志深法能很好地描述整个结晶过程。Jeziorny法处理结果表明,PPD主结晶阶段的Avrami指数(n)为1.68~1.78,晶体生长为准二维生长。莫志深法处理结果表明,在单位结晶时间里达到某一相对结晶度所需的降温速率随相对结晶度的增加而增大。用Kissinger方程求得PPD的非等温结晶活化能为-31.9kJ/mol。     

CHARACTERIZATIONS ON THE THIXOTROPY-LOOP TESTS USING UCM MODEL WITH A RATE-TYPE KINETIC EQUATION

The theoretical characterizations on the triangular-form thixotropy-loop tests of an LDPE melt (PE-FSB-23D022/Q200) were conducted in the present paper by using a new thixotropy model, which is constituted by the upper convected Maxwell model and a rate-type kinetic equation. The new thixotropic Maxwell model can partially describe well three reported thixotropy-loop experiments by comparison with the previous calculations of the variant form of the thixotropy-type Huang model. It is noted that the stress deviations between the experiments and the predictions of the new thixotropic Maxwell model are much slighter than those deviations obtained by using the variant Huang model at the same condition, although both models include five parameters. The constitution of the new thixotropic Maxwell model is more reasonable than that of the variant Huang model.     

Poisson-Boltzmann与Donnan模型计算压实膨润土孔隙水与外部溶液间离子平衡的差异性比较

压实膨润土孔隙水与外部溶液之间的离子平衡是影响离子在压实膨润土中扩散的影响因素之一, 表征这一平衡的离子平衡系数可用压实膨润土的宏观属性参数通过Donnan 模型计算得到. 通过对膨润土主体矿物蒙脱石的TOT层结构单元进行简化, 构建了一个压实膨润土的单类孔隙结构模型, 辅以一个尺度变量H, 用Poisson-Boltzmann (PB)理论模型计算上述离子平衡系数. 对比计算结果, 发现PB模型计算的离子平衡系数总是大于Donnan模型的结果, 而参数H是联系这两种模型之间的桥梁. 通过对参数H取极限H→0, 实现了从PB模型到Donnan 模型的数学变换, 并从机理上讨论了两种模型之间的差异及其在实际扩散问题中的应用.分析表明PB模型更符合离子在压实膨润土中扩散的实际情况, 更适于处理实际扩散问题.     

A Simplified Model for the Etch Rate of Novolac-Based Photoresist

The novolac-based resins used as positive-tone photoresists are frequently etched in an oxygen plasma. It is desirable to have a predictive model of the photoresist etch rate but, for process improvement, control, and analysis, the development of a rigorous mechanistic model is impractical. Instead, a simplified mechanistic model is derived, here, according to the method proposed by Hougen and Watson for the study of fluid–solid interactions. This model derivation method is employed in order to arrive at a functional form that represents chemical etching of the resist by oxygen radicals, assisted by the plasma ion flux. Values for model parameters are determined from process data by nonlinear regression. The quality of the model fit to the data is tested statistically.     

动力学方程控制表面反应的模拟模型和方法

提出了动力学方程控制表面反应的模拟模型和方法.该模型从最基本的质量作用定律出发,获得表面反应的动力学方程.而表面反应通过格子模拟反应器进行.通过表面催化样板反应"CO表面催化氧化"检验了该模拟模型和方法,与实验结果吻合.该模型可在其它复杂的表面催化反应体系中推广应用.     

A microscopic model for calculations of chemical processes in aqueous solutions

A new microscopic model for calculations of chemical processes in aqueous solutions is presented. The model, referred to here as the “surface constrained soft sphere dipoles” (SCSSD) model avoids the problems of the continuum models by explicitly including the solvent molecules. Each solvent molecule is represented as a point dipole attached to the center of a soft sphere. The solvation energy is evaluated by minimizing the solute-solvent energy with respect to the orientations and positions of those dipoles while constraining the surface dipoles to have the orientations and positions of the bulk solvent. The model is demonstrated by calculating the energetic of charge separation in aqueous solution and evaluating the corresponding dielectric constant. The SCSSD model can be used for quantitative studies of ionic reactions in solutions. This is demonstrated by calculation of the potential surface for the dissociation of formic acid in aqueous solution.     

煤焦在燃烧过程中孔隙结构变化的模拟

煤焦在燃烧过程中的物理特性，如比表面积和孔径分布会发生连续变化，直接测量煤焦在燃烧过程中的孔隙结构变化很困难，但可以通过合适的数学模型来观察，二维的圆柱孔模型已大量用来对煤焦气化与燃烧过程中表面积和孔隙结构的变化进行模拟，这个模型把孔隙分成两大部分－－大孔与小孔，因为小孔构成比表面的绝大部分，所以在反应过程中比表面积的变化可以由单一小孔模型来拟合，本文采用了用Ｔｓｅｎｇ和Ｅｄｇａｒ提出的孔模型对几     

Solvation effect on conformations of 1,2:dimethoxyethane: charge-dependent nonlinear response in implicit solvent models

The physical content of and, in particular, the nonlinear contributions from the Langevin-Debye model are illustrated using two applications. First, we provide an improvement in the Langevin-Debye model currently used in some implicit solvent models for computer simulations of solvation free energies of small organic molecules, as well as of biomolecular folding and binding. The analysis is based on the implementation of a charge-dependent Langevin-Debye (qLD) model that is modified by subsequent corrections due to Onsager and Kirkwood. Second, the physical content of the model is elucidated by discussing the general treatment within the LD model of the self-energy of a charge submerged in a dielectric medium for three different limiting conditions and by considering the nonlinear response of the medium. The modified qLD model is used to refine an implicit solvent model (previously applied to protein dynamics). The predictions of the modified implicit solvent model are compared with those from explicit solvent molecular dynamics simulations for the equilibrium conformational populations of 1,2-dimethoxyethane (DME), which is the shortest ether molecule to reproduce the local conformational properties of polyethylene oxide, a polymer with tremendous technological importance and a wide variety of applications. Because the conformational population preferences of DME change dramatically upon solvation, DME is a good test case to validate our modified qLD model. The present analysis of the modified qLD model provides the motivation and tools for studying a wide variety of other interesting systems with heterogeneous dielectric properties and spatial anisotropy.     

Two-dimensional model of molecule—surface scattering

Rotational excitation in molecule—surface scattering is analyzed in a simple quantum two-dimensional model. The true molecule—surface potential is replaced by a hard-core ellipsoid. The model accurately describes the dominant features of the transition probability distribution. This is demonstrated by comparing rotational transition probabilities calculated from our model with the results of IOS calculations of Voges and Schinke using realistic potential. The model predicts a new scaling property concerning the anisotropy parameter of the potential surface.     

Damped torsional oscillator model for polymer molecules

At present the widely used model for explaining viscoelastic and dielectric properties of polymer solutions is that of Rouse and Bueche. Here the polymer molecule is considered as an array of Gaussian subunits, each of which acts as an entropy spring. The motion of these segments is damped by the viscous drag of the surrounding solvent (RB model). An alternative model is presented, in which the segments are torsional oscillators consisting of two or three backbone links, and the damping is due to hindered internal rotation (DTO model). The mathematical treatment of these two models is essentially identical, but the physical interpretation of the constants used is very different. The DTO model has previously been applied by one of us to the interpretation of viscoelastic data. It is here applied to the interpretation of dielectric loss data. It is shown that dielectric measurements in dilute solution should very readily discriminate between the two approaches. Finally it is shown that the relaxation time computed from the DTO model is in closer agreement with published NMR data on poly(propylene oxide) 2025, than is the RB relaxation time. The postulates of the DTO model appear to be confirmed for this low molecular weight polymer. An even more sensitive distinction should be available by studies of the relaxation time as a function of polymer concentration.     

Mathematical modeling of the relation between myosin phosphorylation and stress development in smooth muscles

In this paper the 4-state latch bridge model proposed by Rembold and Murphy is expanded; first by incorporation of the analytical expression of Ca2+ dependent MLCK activation from the work of Kato et al. and second, by inclusion of the myosin dephosphorylation based on the Michaelis-Menten kinetics. The analysis of the proposed model and the comparison with the original model results as well as with the experimental data is presented. The model is able to predict the steady-state isometric stress and the myosin phosphorylation in dependence on steady cytosolic [Ca2+] as well as the temporal evolution of the system in dependence on the input Ca2+ signal in the form of biphasic transient, whereby our model results are in several aspects in better agreement with experimental observations.     

The titration of clay minerals II. Structure-based model and implications for clay reactivity

The potentiometric titration and CEC data presented in part I are modeled in this paper, part II. Two models are compared: the two pK, three complexation sites plus exchange sites nonelectrostatic model developed by Baeyens and Bradbury and a model based on the MUSIC approach developed by Hiemstra and Van Riemsdijk. Both morphological and structural information is used to develop this new model. Morphological information is taken from the literature, while structural information is taken from the literature and constrained by supporting FTIR experiments. The Baeyens and Bradbury model is found to reproduce the general tendency of the titration curve, whereas the model based on the Hiemstra and Van Riemsdijk MUSIC approach provides a better fit to the experimental data. The former uses only 3 edge reaction sites, whereas the latter uses at least 27 edge reaction sites. Five main reactive sites are sufficient to fit the MUSIC model curve, but the model allows us to derive the properties of 22 other reactive sites. Logically, the greater the number of sites, the better the fit. Nevertheless, fewer adjustable parameters are necessary for the Hiemstra and Van Riemsdijk MUSIC model than for the Baeyens and Bradbury model, thanks to structural and morphological constraints. The precision of the potentiometric titration curve is insufficient to verify that the properties of the 27 sites given by the MUSIC model are effective. Thus, we coupled some properties of clay minerals, such as dissolution, to the modeled acid-base properties of these sites to assess our model. We then questioned the ability of simplified models such as the Baeyens and Bradbury model to predict the interactions between clay minerals and solutions in natural environments. In addition, we derived the cation exchange selectivity coefficients for CaCl+ ionic pairs and H+ from our CEC data and gave an estimate for the CaOH+ selectivity coefficient.     

Analysis and Simulation of the Dynamics of a Catalyzed Model Reaction: CO Oxidation on Zeolite Supported Palladium

The particle size effect on the oscillatory behavior during CO oxidation over zeolite-supported Pd catalysts is simulated with the help of a deterministic point model and a stochastic mesoscopic model. The point model is developed on the basis of the well-known Sales–Turner–Maple model, which is modified to consider the slow processes of oxidation and reduction of the Pd bulk as well as the effects of the bulk oxidation on the catalyst activity. It is demonstrated that the point model developed can simulate many experimental trends, e.g., the dependence of the catalytic activity and the waveform of the oscillations on the particle size and the pretreatment of the catalyst, as well as the counterclockwise hysteresis, depending on the reaction rate during the cyclic variation of the CO inlet concentration. The higher activity of the smaller particles can be explained by the attainment of a more reduced state of Pd in smaller particles in the course of the reaction. The stochastic model simulates the reaction by a Markovian chain of elementary stages of the reaction. The model variables are the numbers of reagent atoms. Transition probabilities of the stochastic model are chosen in accordance with the rates of the developed point model. It is shown that intrinsic fluctuations and correlations of stochastic variables can significantly change the reaction dynamics on nm-sized particles by extending the oscillatory region in the parameter space.     

Analysis of the small-angle X-ray scattering of linear polyethylene by a paracrystalline lattice model

Three paracrystalline lattice models for the interpretation of the small angle scattering of polyethylene are discussed: The “lattice model”, the “stapel model” (often referred to as the lamellar stack model) and the “proportional model”. While the applicability of the first model is restricted, the latter models differ in the statistical assumptions of lamellar and interlamellar thickness distributions. The principal advantage of the proportional model over the stapel model is its applicability through the adjustment of only three parameters: long period, crystallinity and one statistical parameter. Small angle X-ray curves of linear polyethylene are interpreted by the proportional model. The results are in good agreement with stapel model calculations.     

Improved partition equilibrium model for predicting analyte response in electrospray ionization mass spectrometry

A previously proposed partition equilibrium model for quantitative prediction of analyte response in electrospray ionization mass spectrometry is modified to yield an improved linear relationship. Analyte mass spectrometer response is modeled by a competition mechanism between analyte and background electrolytes that is based on partition equilibrium considerations. The correlation between analyte response and solution composition is described by the linear model over a wide concentration range and the improved model is shown to be valid for a wide range of experimental conditions. The behavior of an analyte in a salt solution, which could not be explained by the original model, is correctly predicted. The ion suppression effects of 16 : 0 lysophosphatidylcholine (LPC) on analyte signals are attributed to a combination of competition for excess charge and reduction of total charge due to surface tension effects. In contrast to the complicated mathematical forms that comprise the original model, the simplified model described here can more easily be employed to predict analyte mass spectrometer responses for solutions containing multiple components. Copyright © 2008 John Wiley & Sons, Ltd.     

Scaled one-electron hamiltonian model for open-shell LCAO–MO–SCF calculations: Comparison with the restricted open-shell method of roothaan

The performance of a recently proposed scaled one-electron Hamiltonian (SOEH ) model is tested against parallel sets of restricted open-shell calculations by the method of Roothaan. It is found that the energy calculated by SOEH model, in general, lies slightly higher than the energy computed by the restricted open-shell method of Roothaan lending credibility to the application of variational argument to the scaled pseudoenergy functional (E^av) for deriving the SOEH model. The numerical stability of the converged SOEH energy with respect to changes in trial vectors indicates the reliability of the method. The SOEH model is shown to perform well in the calculation of geometries of radicals and ions. The convergence behavior of the SOEH model is compared with that of the restricted open-shell method of Roothaan.     

正癸烷燃烧机理及航空煤油点火延时动力学模拟

以单一正癸烷作为国产航空煤油的单组分替代模型, 应用自有的碳氢燃料反应机理生成程序ReaxGen-Combustion构建了燃烧反应的详细机理. 以国产煤油在加热型激波管上的燃烧实验为参考, 对比研究了文献报道的3组分替代模型(模型Ⅰ)、2组分替代模型(模型Ⅱ)以及本文的单组分替代燃烧反应机理(模型Ⅲ)在预测我国航空煤油点火延时特性方面的实用性. 结果表明, 温度在1052~1538 K时, 模型Ⅰ预测的点火延时与实验值相差较大; 模型Ⅲ在温度高于1176 K时的预测值与实验值符合较好, 在1052~1176 K之间时则相差较大; 模型Ⅱ与模型Ⅲ预测值符合很好, 由于前者考虑了低温反应机理, 因而对1052~1176 K区间的预测精度与模型Ⅲ相比有所改善. 计算还发现, 模型Ⅱ中添加的20%(质量分数)1,2,4-三甲基苯对高温段点火延时未产生明显影响.     

The femtosecond dynamics of electron transfer in a modified photosynthesis reaction center: Quantum, classical, and kinetic analyses

The dynamics of electron transfer in a modified photosynthesis reaction center in which electron transfer from the bridge to the acceptor is blocked is considered. A microscopic model of the process is suggested. Within this model, the diabatic electronic states of the donor and bridge are described by one-dimensional displaced harmonic oscillators. The dynamics of the population of electronic states is calculated by the quantum method of wave packets and classical and kinetic modeling. The suggested model is used to study the qualitative dependence of the dynamics of electron transfer on the nonadiabatic interaction potential. The parameters of the model are determined by comparing the experimental and calculated absorption spectra of the product of electron transfer. It is shown that kinetic models can be used to approximately describe the dynamics of electron transfer in reaction centers. The boundaries of the applicability of the kinetic method are considered.     

Crystallographic B-factors highlight energetic frustration in aldolase folding

Kinetic simulations of the folding and unfolding of the mammalian TIM barrel protein aldolase were conducted to determine if a minimalist monomeric Gō model, using the native structure to determine attractive energies in the protein chain, could capture the experimentally determined folding pathway. The folding order, that is, the order in which different secondary structures fold, between the Gō model simulations and that from hydrogen-deuterium exchange experiments, did not agree. To explain this discrepancy, two alternate variant of the basic Gō model were simulated: (1) a monomer Gō model with native contact energies weighted by a statistical potential (SP model) and (2) a monomer Gō model with native contact energies inversely weighted by crystallographic B factors (B model). The B model demonstrated the best agreement between simulation and experiments. The success of the B model is attributed to the ability of B factors to highlight local energetic frustration in the aldolase structure which results in weaker native contacts in these frustrated regions. The predictive success of the B model also reveals the potential use of B factor information for energetic weighting in general protein modeling.