一种模拟沸石分子筛扩散反应的Monte Carlo方法

在考虑了沸石分子筛中吸附、脱附、扩散和反应的特点的基础上，发展了一种模拟沸石分子筛中扩菜反应的ＭｏｎｔｅＣａｒｌｏ方法。用Ｌａｎｇｍｕｉｒ吸附等温式验证了该方法的合理性，并证明了轰击频率λ代表了反应物气相压力，而吸附、脱附几率之比即为吸附平衡常数。同时，通过改变轰击频率λ和反应几率Ｐｔ，得到了吸附等温线和有效因子随Ｔｈｉｅｌｅ模数变化的曲线。     

中性载体阳离子电极膜电化学的计算机数字模拟

对中性载体阳离子电极膜3种可能的传输模式进行了数字模拟,研究了淌度、浓度、配合物稳定常数以及阴离子亲脂性对电极电位响应的影响。模拟结果对于廓清这类电极响应机理,解释有关实验现象及指导新电极设计提供了有用的信息。     

Voltammetric sizing of a sphere.

The size of a glass sphere positioned in the center of a microdisk electrode is determined by using a simple electrochemical procedure and is confirmed, additionally, by a microscopical measurement of the sphere at the time of the electrochemical measurement. The cyclic voltammetric response of the naked electrode and of the electrode with the sphere positioned in its center is recorded over a wide range of scan rates (0.002-1.5 V s(-1)). The size of the sphere is then determined by comparison of the experimental voltammogram with simulations for each individual scan rate.     

On the physical adsorption of gases on carbon materials from molecular simulation

In this paper we present a series of work covering a range of aspects relating molecular simulation to experiment. The importance of surface mediation type effects to the adsorption of simple and complex gases is demonstrated. Coupled with the adsorption of simple gases is their projection area when used for surface area determination. The pressure dependence of a projection area is demonstrated for argon at 77 and 87.3 K. A simple model is used to account for the degree of graphitisation of a surface is demonstrated and used to account for the isosteric heat behaviour of non-graphitised carbon blacks. Turning from surfaces to porous solids, an alternative treatment of experiment data (either sub or super critical) is presented that avoids the ambiguity of excess amounts adsorbed. Using this method one is able to obtain pore size distributions and amounts adsorbed without relying on such things as helium expansion volumes. Since this type of method is usually applied to composite solids we also demonstrate the correct method for calculating the heat of adsorption using independent sets of simulations. The final topic covered in this paper is an example of the information that can be gained from the heat capacity of an adsorbed phase.     

钙钛矿型CaSiO_3解压缩和重压缩过程的分子动力学模拟

用分子动力学模拟方法 (MD)研究了 3 0 0K时钙钛矿型CaSiO3 ,从高压到负压的解压缩过程 .MD模拟获得的P V关系与实验数据相近 ,与已报道的MD模拟数据基本一致 ,所得体积模量也在实验数据分布范围内 .减压缩和重压缩的MD模拟数据与实验结果相似 .钙钛矿型CaSiO3 解压缩成非晶态时 ,存在两个结构破坏阶段 :破坏硅氧八面体和破坏钙氧二十面体 .当钙氧二十面体被破坏后 ,重压缩不能得到钙钛矿型结构 .只要钙氧二十面体未被破坏 ,重压缩可恢复钙钛矿型结构 .本研究得到的结果尚未见相关报道 .由MD模拟数据计算了CaSiO3 系统的红外光谱 ,分析这些数据可知钙钛矿型CaSiO3 解压缩非晶化是一个二阶软模相变 .研究表明钙钛矿型CaSiO3 结构存在一个等容的亚稳极限 ,其解压缩非晶化是一个受动力学控制的亚稳状态 .     

Dynamic Study of Gemini Surfactant and Single-chain Surfactant at Air／Water Interface

Molecular dynamics (MD) simulation are used to study the properties of gemini surfactant of ethyl-α,ω-bis(dodecyldimethylammonium bromide) (C12C2C12) and dodecyl-trimethylammonium bromide (DTAB) at the air/water interface, respectively. In the two systems,t he surfactant concentrations are both 28 wt.%, and other conditions are also the same. After reaching the thermodynamic equilibrium, the concentration profiles, the radial distributions functions (RDF) and the mean squared displacement (MSD) are investigated. The results reveal that the surface activity of C12C2C12 surfactant is higher than DTAB surfactant.     

分子筛模型中非骨架离子位置的研究

从理论上研究了非骨架阳离子可能的分布位置．计算了分子筛晶体基本结构单元在整个空间的静电势分布，找出空间中静电势的各个极小点．在静电势最低点逐个放置钠离子，优化，重新计算模型分子筛的电势分布．研究了阳离子放置过程中静电势变化的趋势．     

Free energy perturbation approach to the critical assessment of selective cyclooxygenase-2 inhibitors

The discovery of selective cyclooxygenase-2 (COX-2) inhibitors represents a major achievement of the efforts over the past few decades to develop therapeutic treatments for inflammation. To gain insights into designing new COX-2-selective inhibitors, we address the energetic and structural basis for the selective inhibition of COX isozymes by means of a combined computational protocol involving docking experiment, force field design for the heme prothetic group, and free energy perturbation (FEP) simulation. We consider both COX-2- and COX-1-selective inhibitors taking the V523I mutant of COX-2 to be a relevant structural model for COX-1 as confirmed by a variety of experimental and theoretical evidences. For all COX-2-selective inhibitors under consideration, we find that free energies of binding become less favorable as the receptor changes from COX-2 to COX-1, due to the weakening and/or loss of hydrogen bond and hydrophobic interactions that stabilize the inhibitors in the COX-2 active site. On the other hand, COX-1-selective oxicam inhibitors gain extra stabilization energy with the change of residue 523 from valine to isoleucine because of the formations of new hydrogen bonds in the enzyme-inhibitor complexes. The utility of the combined computational approach, as a valuable tool for in silico screening of COX-2-selective inhibitors, is further exemplified by identifying the physicochemical origins of the enantiospecific selective inhibition of COX-2 by -substituted indomethacin ethanolamide inhibitors.