水-甲醇体系的Monte Carlo分子膜拟

本文应用Monte Carlo分子膜拟方法对水、甲醇水-甲醇的1:1混合物、甲醇无限稀释时的水溶液和水无限稀释时的甲醇溶液等五个体系进行了研究。采用TIP分子位能函数, 得到了上述体系的热力学性质、原子径向分布函数、分子氢键配位数分布。并以Monte Carlo分子模拟获得的结构函数与X射线衍射实验结果进行了比较。     

单官能团单体C存在下A_2-B_2型缩聚产物的分子量分布

分别以概率统计方法和Monte Carlo模拟获得了单官能团单体C存在下A_2-B_2型线型缩聚产物的分子量分布。Monte Carlo方法与统计计算得出一致的结果,表明建立的Monte Carlo模型是合理的。研究了单官能团单体C的反应活性对缩聚产物分子量分布的影响,对反应体系中存在的各种不同的分子类型的消长情况分别进行了动态模拟。     

流体熵相关性质的Monte Carlo模拟新方法

提出了用于计算实际体系熵相关性质的Monte Carlo 多级取样分子模拟方法.应用这一方法,对硬球流体的化学势及Helmholtz 自由能进行了估算,得到了满意的结果.计算化学势时,不存在通常试验粒子方法所遇到的高密度问题.该方法特别适合规律性的系统研究,较之普通模拟方法要有效得多.模拟得到的硬球体系无限稀释组份的超额化学势与对比直径的关系,在相变区域为一条双凹曲线;无论是在相变区还是在单相区,Carnahan-Starling 公式对这一关系的描述均有较大偏差.     

流体熵相关性质的Monte Carlo模拟新方法

提出了用于计算实际体系熵相关性质的Monte Carlo多级取样分子模拟方法. 应用这一方法, 对硬球流体的化学势及Helmholtz自由能进行了估算, 得到了满意的结果. 计算化学势时, 不存在通常试验粒子方法所遇到的高密度问题. 该方法特别适合规律性的系统研究, 较之普通模拟方法要有效得多. 模拟得到的硬球体系无限稀释组份的超额化学势与对比直径的关系, 在相变区域为一条双凹曲线; 无论是在相变区还是在单相区, Carnahan-Starling公式对这一关系的描述均有较大偏差.     

连锁聚合反应的动态Monte Carlo模拟及其浓度效应

用动态Monte Carlo方法和键长涨落模型相结合分别对自由基聚合、离子聚合动力学过程进行了计算机模拟, 在三维空间中统计了转化率、聚合度、分子量的数量和重量分布, 并将模拟结果和理论预测进行了对比分析, 同时得到不同浓度下均方回转半径与平均链长的标度关系, 说明三维空间中标度指数与浓度的相关性. 还初步考察了扩散对反应体系分子量的影响. 说明动态Monte Carlo方法用于研究连锁聚合反应动力学过程是有效的, 而且能够同时得到经典的聚合反应Monte Carlo模拟方法所难以得到的链构象、分子扩散等空间相关的信息.     

超临界正己烷—甲醇溶液体系的Monte Carlo模拟：浓度对溶液构型？…

用Ｍｏｎｔｅ Ｃａｒｌｏ分子模拟方法对正己烷－甲醇溶液体系处于临界状态时微观结构随浓度变化的情况进行了研究。模拟采用随机边界条件，得到了不同浓度条件下体系各基团的径向分布函数。结果表明，在溶液临界状态，正己烷分子有较强的聚集行为，随着溶液中甲醇浓度的变大，正己烷分子的聚集程度逐渐下降；甲醇分子周围正己烷分子的分布有一个最佳的浓度范围，在甲醇浓度为２２．５％时，配位数达到最大；超临界正己烷－甲醇溶液     

超临界正己烷甲醇溶液体系的 Monte Carlo 模拟──浓度对溶液构型影响的考察

用ＭｏｎｔｅＣａｒｌｏ分子模拟方法对正己烷－甲醇溶液体系处于临界状态时微观结构随浓度变化的情况进行了研究。模拟采用随机边界条件，得到了不同浓度条件下体系各基团的径向分布函数。结果表明，在溶液临界状态，正己烷分子有较强的聚集行为，随着溶液中甲醇浓度的变大，正己烷分子的聚集程度逐渐下降；甲醇分子周围正己烷分子的分布有一个最佳的浓度范围，在甲醇浓度为２２．５％时，配位数达到最大；超临界正己烷－甲醇溶液中甲醇分子的作用在不同浓度条件下具有不同的特点，当甲醇浓度较低时，甲醇分子之间有氢键的作用，随着甲醇浓度的提高，甲醇分子之间氢键的作用变得越来越弱，而且分子的取向变得越来越无序。     

La—LaCl3—KCl熔体的计算机模拟研究

用Monte Carlo方法对La-LaCl_3-KC1体系在1223K时的结构进行了计算机模拟研究,得到了熔体中诸离子对的偏径向分布函数和体系在1223K的势能和内能。结果表明,在熔体中La(Ⅲ)有相当一部分以LaCl_6~3六配位形式存在,而La(Ⅱ)则主要以LaCl_4~2四配位形式存在。结果还表明,熔体中的自由体积分布不均匀,存在许多不规则的空孔和缝隙,其中XK~+·YC1~-集团内的缝隙比纯KC1熔体中的缝隙明显增多。     

电解质模型流体的MonteCarlo分子模拟

采用正则系综Monte Carlo分子模拟方法，模拟得到了一系列状态条件下1：1 价和2：2价对称电解质模型流体的径向分布函数、构型能和压缩因子的“机器实验 数据”，这些数据对构筑电解质溶液的新参考流体具有重要意义。模拟过程中，电 解质溶液被简化为硬球离子和硬球点偶扳子的混合物，离子—离子、离子—偶极和 偶极—偶极长程位能计算采用了Ewald求和方法。最后，对离子—偶极子混合物的 结构和热力学性质与体系温度、密度和浓度的关系进行了分析和讨论。     

芳氧基稀土化合物催化ε-己内酯开环聚合的链转移反应及Monte Carlo模拟

设计了三(2,6-二叔丁基-4-甲基苯氧基)稀土配合物[Ln(OAr)3]催化ε-己内酯(CL)开环聚合反应体系中的链转移反应.证明存在分子内/分子间酯交换反应,并验证聚合物在链转移反应过程中发生明显的分子量下降和分子量分布变宽.用Monte Carlo方法对聚ε-己内酯(PCL)链转移过程进行了计算机模拟,定量解释了链转移过程中环链分子数量比值的增加、环和链分子的数均分子量下降以及分子量分布变宽等现象,并发现随着链转移反应的进行,体系的各项参数先后趋于定值,最终体系将进入"环链平衡"状态.     

Hydration structure around a non-electrolyte molecule as revealed by monte carlo calculation: have you seen the “iceberg”?

Monte Carlo calculations have been carried out for pure water and an infinitely dilute aqueous solution of methanol at 298.15 K at ordinary density by the Metropolis method in NTV ensemble. The hydration structure around methanol revealed by a graphic display technique indicates a structure-forming effect near the hydrophobic group.     

Solubility of carbon dioxide in aqueous solutions of methanol. Predictions by molecular simulation and comparison with experimental data

The solubility of carbon dioxide in pure methanol, and in aqueous solutions of methanol, was computed using the Gibbs ensemble Monte Carlo (GEMC) technique for 313, 354, and 395 K at pressures up to 9 MPa. Three solvent mixtures (of methanol and water) with methanol mole fractions of 10, 50, and 75 mole percent (in the gas-free solvent mixture) were studied. The Monte Carlo simulations were conducted in an isothermal-isobaric ensemble applying effective pair potentials for the pure components from literature. Common mixing rules without any adjustable binary interaction parameters were used to describe the interactions between the mixture components. Overall, a good agreement between simulation results and recently published experimental data is achieved.     

Test of athermal terms of activity coefficient models by Monte Carlo simulation with hard-core models

The infinite dilution activity coefficients of exactly athermal fluids were calculated by Monte Carlo simulation with hard-core models. The hard-core models used in this work were hard-sphere and hard-spherocylinder models. The Widom test particle method was adopted to calculate the residual chemical potentials of solutes in pure solvent and in pure solute solutions. The infinite dilution activity coefficients of solutes were obtained from the residual chemical potentials of solutes. The infinite dilution activity coefficients calculated by Monte Carlo simulation were compared with those of athermal terms in activity coefficient equations. Staverman–Guggenheim equation overestimates the activity coefficients. The deviations of activity coefficients increase with increasing the hard-core volume of solute. Flory–Huggins equation based on molar volume gives good results for the hard-spherocylinder systems. Elbro-FV equation gives good results for both the hard-sphere and hard-spherocylinder systems.     

Molecular simulation study on adsorption of methanol/water mixtures in mesoporous silicas modified pore surface silylation

Two types of molecular simulation techniques have been utilized to investigate adsorption of methanol/water mixtures in a mesoporous silica with a hydrophobic pore surface: the NVT-ensemble Molecular Dynamics method with the melt-quench algorithm for modeling a fully-silylated mesoporous silica and the μVT-ensemble Orientaional-Biased Monte Carlo method for calculating adsorption isotherms. Adsorption isotherms of methanol and water at 333 K are calculated for an equi-relative-pressure mixture (each component has the same relative pressure which is defined as the ratio of the partial pressure to the saturation pressure of the pure gas) together with pure gases. In the case of the pure gas, water hardly adsorb even at elevated pressures, while the adsorption isotherm for methanol shows the condensable adsorption. On the other hand, in the case of the mixture, water molecules are substantially adsorbed along with methanol molecules, showing an isotherm representing the condensation mechanism. In addition, it is found that the separation factor of methanol to water is the highest in the case of monolayer adsorption from a liquid mixture.     

Four-m calculations on methanol-water solutions

The structural characteristics of methanol in aqueous solutions, on a molecular level, can be elucidated by four types of calculation: molecular orbital, multiparametric optimization of intermolecular potential function, Monte Carlo, and molecular dynamics. As a first step, the potential between water and methanol was determined by ab initio LCAO SCF molecular orbital calculations with the STO—3G basis set and subsequent multiparametric fitting. This and water—water potentials were used for Monte Carlo calculation on an aqueous methanol solution containing a 1:216 mole ratio of methanol to water. Hydration around methanol is briefly discussed.     

Chain conformation and solvent partitioning in reversed-phase liquid chromatography: Monte Carlo simulations for various water/methanol concentrations

Many structural models for the stationary phase in reversed-phase liquid chromatography (RPLC) systems have been suggested from thermodynamic and spectroscopic measurements and theoretical considerations. To provide a molecular picture of chain conformation and solvent partitioning in a typical RPLC system, a particle-based Monte Carlo simulation study is undertaken for a dimethyl octadecyl (C(18)) bonded stationary phase on a model siliceous substrate in contact with mobile phases having different methanol/water concentrations. Following upon previous simulations for gas-liquid chromatography and liquid-liquid phase equilibria, the simulations are conducted using the configurational-bias Monte Carlo method in the Gibbs ensemble and the transferable potentials for phase equilibria force field. The simulations are performed for a chain surface density of 2.9 micromol/m(2), which is a typical bonded-phase coverage for mono-functional alkyl silanes. The solvent concentrations used here are pure water, approximately 33 and 67% mole fraction of methanol and pure methanol. The simulations show that the chain conformation depends only weakly on the solvent composition. Most chains are conformationally disordered and tilt away from the substrate normal. The interfacial width increases with increasing methanol content and, for mixtures, the solvent shows an enhancement of the methanol concentration in a 10 Angstrom region outside the Gibbs dividing surface. Residual surface silanol groups are found to provide hydrogen bonding sites that lead to the formation of substrate bound water and methanol clusters, including bridging clusters that penetrate from the solvent/chain interfacial region all the way to the silica surface.     

超临界CO2中甲醇和乙醇无限稀释扩散系数的分子动力学模拟与实验测定

采用分子动力学模拟(MD)方法对甲醇和乙醇在超临界二氧化碳中的无限稀释扩散系数进行了模拟计算, 并应用泰勒分散理论, 采用超临界色谱仪对模拟结果进行了实验验证. 模拟计算值与实验值较吻合, 且变化规律基本一致, 表明采用这种新方法可以准确有效地预测超临界体系的扩散性质, 能够方便地应用于工程设计.     

Difference spatial distribution function analysis of aqueous solutions. II. Hydration structures of dimethyl ether, 180 degrees ethyl methyl ether and 0 degree ethyl methyl ether solutions

Monte Carlo simulations are systematically presented to demonstrate the influence of the hydrophobic group's steric bulk on hydration structure. We have simulated a dimethyl ether (DME), two conformations for ethyl methyl ether (0 degree EME and 180 degrees EME), and 0 degree ethanol solutions. Spatial distribution function (SDF), goo(x,y,z) and difference SDF (DSDF), delta goo(x,y,z), obtained from MC simulation in an infinitely dilute aqueous solution of ether show the three-dimensional probability of an atom-atom pair distribution between solute and solvent atoms. Based on the results of SDF in an infinitely dilute aqueous solution of ether, the distribution of hydration water molecules can be divided into hydrogen acceptor (HA) and hydrophobic hydration (HH), regions, and the spatial orientation of the hydrogen-bonded water in the HA region is found to form a triple-layer structure, as it does in alcohol solutions. From the results of an analysis of the DSDF delta goo(x,y,z) between the SDFs of EME and DME, it is apparent that the distribution changes of hydration water molecules in ether solutions are essentially similar to those in the alcohol solutions. Further, we show that the hydration water molecules are distributed mainly in the stable area in the binding energy's (BE) contour maps for each region.     

Structure of the liquid-vapor interface of water-methanol mixtures as seen from Monte Carlo simulations

Monte Carlo simulation of the vapor-liquid interface of water-methanol mixtures of different compositions, ranging from pure water to pure methanol, have been performed on the canonical (N, V, T) ensemble at 298 K. The analysis of the systems simulated has revealed that the interface is characterized by a double layer structure: methanol is strongly adsorbed at the vapor side of the interface, whereas this adsorption layer is followed at its liquid side by a depletion layer of methanol of lower concentration than in the bulk liquid phase of the system. The dominant feature of the interface has been found to be the adsorption layer in systems of methanol mole fractions below 0.2, and the depletion layer in systems of methanol mole fractions between 0.25 and 0.5. The orientation of the molecules located at the depletion layer is found to be already uncorrelated with the interface, whereas the methanol molecules of the adsorption layer prefer to align perpendicular to the interface, pointing straight toward the vapor phase by their methyl group. Although both the preference of the molecular plane for a perpendicular alignment with the interface and the preference of the methyl group for pointing straight to the vapor phase are found to be rather weak, the preference of the methyl group for pointing as straight toward the vapor phase as possible within the constraint imposed by the orientation of the molecular plane is found to be fairly strong. One of the two preferred orientations of the interfacial water molecules present in the neat system is found to disappear in the presence of methanol, because methanol molecules aligned in their preferred orientation can replace these water molecules in the hydrogen-bonding pattern of the interface.