水团簇(H2O)n中多体作用强度与水分子间距的关系

使用高精度从头算方法(含基组重叠误差校正)计算了水团簇(H₂O)_n (n = 8, 10, 16, 20, 22, 24)中的所有二体、三体和四体作用能,分析了水团簇中的多体效应.研究表明,二体作用对体系总作用能的贡献高达70%以上,三体作用对总作用能的贡献可高达25%,四体作用在总作用能中所占比例不超过3%,五体及以上多体作用能在总作用能中所占比例更小,不超过0.5%.本文研究还表明,两个水分子间距小于0.68 nm的二体作用、三个和四个水分子中最近的两个水分子间距小于0.31 nm的三体和四体作用对体系总作用能的贡献高达99.4%.因此,以生物体系为对象的分子模拟方法应该具备准确地模拟两个水分子间距小于0.68 nm的二体作用、三个和四个水分子中最近的两个分子间距小于0.31 nm的三体和四体作用的能力.     

Monte Carlo liquid water simulation with four-body interactions included

A four-body interaction potential for water molecules is derived. The new terms are combined with previously developed two- and three-body potential terms and applied in a Metropolis—Monte Carlo simulation at 298 K for an (N, V, T) ensemble with 512 water molecules. Improvements, relative to the results using only two-body, or two- and three-body interactions, are reported for the correlation function g(OO), for the X-ray and neutron-beam scattering intensities, and for the enthalpy.     

Cooperative effects in extended hydrogen bonded systems involving O?H groups. Ab initio studies of the cyclic S4 water tetramer

The additional energy stabilization due to cooperative effects was calculated in extended hydrogen bonded systems O? H ?O? H ?O? H with unidirectional (homodromic) orientation of the O? H groups. Ab initio restricted Hartree Fock, MP2 and MP3 calculations with geometry optimization and BSSE correction have been performed using the GAUSSIAN 83 program package for the ground states of the linear water dimer with C_s symmetry and the cyclic water tetramer with S₄ symmetry. The latter represents the smallest possible, experimentally observed cooperative structure. A new definition for a cooperativity parameter is proposed. The definition is based on the two-body, non-neighbour interaction energy, plus three- and four-body contributions, including one-body deformation terms in relation to the total interaction energy of the water tetramer. The advantage of this definition is its independence of the reference system, which is necessary in complicated molecular systems with an undefined number of hydrogen bonds, such as disordered or flip-flop systems. According to this definition the energy gain based on cooperativity in the S₄ water tetramer is 29% with the MP3/6-31G** approximation, (30% with HF/4-31G* and 46% with HF/3-21G). The largest contribution of 18% is due to the three-body term on the MP3/6-31G** level, followed by the two-body, non-neighbour term with 11%. The four-body term and the deformation term are in the order of 1% and cancel each other because they have opposite sign.     

A computational study of microsolvation effect on ethylene glycol by density functional method

This study focuses on the conformational analysis of ethylene glycol-(water)n (n=1-3) complex by using density functional theory method and the basis set 6-311++G*. Different conformers are reported and the basis set superposition error corrected total energy is -306.767 5171, -383.221 3135, and -459.694 1528 for lowest energy conformer with 1, 2, and 3 water molecules, respectively, with corresponding binding energy -7.75, -15.43, and -36.28 kcal/mol. On applying many-body analysis it has been found that relaxation energy, two-body, three-body energy have significant contribution to the binding energy for ethylene glycol-(water)3 complex whereas four-body energies are negligible. The most stable conformers of ethylene glycol-(water)n complex are the cyclic structures in which water molecules bridge between the two hydroxyl group of ethylene glycol.     

Computation of some thermodynamic properties of helium–neon and helium–krypton fluid mixtures using molecular dynamics simulation

We have performed molecular dynamics simulations to obtain internal energy and pressure of helium–neon and helium–krypton mixtures at different densities using accurate recently two-body ab initio potentials supplemented by quantum corrections following the Feynman–Hibbs approach. The significance of this work is that the three-body expression of Wang and Sadus [22] was used to improve prediction of the pressures and internal energies of helium + krypton and helium + neon mixtures without requiring an expensive three-body calculation. Our results show a good agreement with the corresponding experimental data.     

Sensitivity of the thermal and acoustic virial coefficients of argon to the argon interaction potential

Second, third, and fourth thermal and acoustic virial coefficients between 100 and 1000 K are computed for different argon interaction models derived from combinations of accurate two- and three-body potentials. Differences between the various interaction models tested mirror the presumed order in the accuracy of these models, but are not well captured at the level of the lowest-order contributions in the virial expansion: While the second- and third-order virial coefficients are found to be rather insensitive to small variations in the two- and three-body potentials, more pronounced differences in higher-order coefficients are currently of limited use in assessing the accuracy of the interaction potential due to difficulties in the unambiguous experimental determination of these higher-order coefficients. In contrast, pressure-volume and speed-of-sound data - both of which are experimentally known to highest accuracies - are found to be insensitive to small variations in the interaction model. All but the least accurate models reproduce experimental pressure-volume and speed-of-sound data near-quantitatively in regions where the (fourth-order) virial expansions apply. All quantities considered are found to be completely unaffected by a non-vanishing quadruple-dipole four-body potential.     

Intermolecular interaction energy of CH4 trimer by symmetry-adapted perturbation theory

The interaction energy for the cyclic CH₄ trimer is studied in terms of symmetry-adapted perturbation theory. The interaction energy around the van der Waals minimum is dominated by attractive dispersion energy, and the repulsive contribution at the smaller angle region is due to the first-order exchange energy. The total interaction energy is approximated by additive two-body components, because of a mutual cancellation between nonadditive three-body ones.     

Effect of three-body interactions on the vapor-liquid phase equilibria of binary fluid mixtures

Gibbs-Duhem Monte Carlo simulations are reported for the vapor-liquid phase coexistence of binary argon+krypton mixtures at different temperatures. The calculations employ accurate two-body potentials in addition to contributions from three-body dispersion interactions resulting from third-order triple-dipole interactions. A comparison is made with experiment that illustrates the role of three-body interactions on the phase envelope. In all cases the simulations represent genuine predictions with input parameters obtained independently from sources other than phase equilibria data. Two-body interactions alone are insufficient to adequately describe vapor-liquid coexistence. In contrast, the addition of three-body interactions results in very good agreement with experiment. In addition to the exact calculation of three-body interactions, calculations are reported with an approximate formula for three-body interactions, which also yields good results.     

基于化学缔合统计理论的链状流体状态方程

基于化学缔合统计理论的链状流体状态方程(EOS)能够反映实际分子的形状、链节成链、缔合等具体信息,在实际流体热力学性质计算中有着广泛应用.一般的链状流体EOS仅考虑相邻链节间的相关性,我们则借助统计力学和计算机模拟结果在模型中纳入了相间链节间的相关性,获得的硬球链流体(HSCF)模型能够更好地预测模型流体的压缩因子和第二维里系数.以HSCF为参考,引入方阱色散微扰项获得了实际方阱链流体(SWCF)EOS;结合根据黏滞球模型导得的缔合项,进一步构建了缔合流体EOS.最近,我们根据微扰理论和积分方程方法又开发了一新的变阱宽方阱链流体(SWCF-VR)模型.SWCF和SWCF-VREOSs可很好地用于计算小分子、聚合物、离子液体等纯流体及混合物的相行为、热焓、表面张力、黏度等热力学及传递性质,显示了模型良好的工程应用价值.本文就本课题组多年来在自由空间范畴内基于化学缔合统计理论开发链状流体EOS及其实际应用作系统的总结.     

Many-body exchange-repulsion in polarizable molecular mechanics. I. Orbital-based approximations and applications to hydrated metal cation complexes

We have quantified the extent of the nonadditivity of the short-range exchange-repulsion energy, E(exch-rep), in several polycoordinated complexes of alkali, alkaline-earth, transition, and metal cations. This was done by performing ab initio energy decomposition analyses of interaction energies in these complexes. The magnitude of E(exch-rep(n-body, n > 2)) was found to be strongly cation-dependent, ranging from close to zero for some alkali metal complexes to about 6 kcal/mol for the hexahydrated Zn(2+) complex. In all cases, the cation-water molecules, E(exch-rep(three-body)), has been found to be the dominant contribution to many-body exchange-repulsion effects, higher order terms being negligible. As the physical basis of this effect is discussed, a three-center exponential term was introduced in the SIBFA (Sum of Interactions Between Fragments Ab initio computed) polarizable molecular mechanics procedure to model such effects. The three-body correction is added to the two-center (two-body) overlap-like formulation of the short-range repulsion contribution, E(rep), which is grounded on simplified integrals obtained from localized molecular orbital theory. The present term is computed on using mostly precomputed two-body terms and, therefore, does not increase significantly the computational cost of the method. It was shown to match closely E(three-body) in a series of test cases bearing on the complexes of Ca(2+), Zn(2+), and Hg(2+). For example, its introduction enabled to restore the correct tetrahedral versus square planar preference found from quantum chemistry calculations on the tetrahydrate of Hg(2+) and [Hg(H(2)O)(4)](2+).     

Towards understanding the effect of electrostatic interactions on the density of ionic liquids

In order to have a better understanding on the electrostatic contribution to the thermodynamic property of ionic liquids (ILs), a two-parameter equation of state (EOS) is developed on the basis of hard sphere perturbation theory by accounting for the dispersion interaction with Cotterman et al.’s EOS for L-J fluid and electrostatic interaction with mean spherical approximation (MSA) approach. The EOS is applicable for the density correlation of molecular liquids, and the resulting parameters, viz. Lennard–Jones dispersive parameter ?/k and soft-core diameter σ, can be used to predict the density of molecular mixtures and the corresponding ILs. The results indicate that the density of IL is always about 10% higher than the corresponding stoichiometric molecular mixture with which the IL is produced as an ionic adduct, for example, IL 1-methyl-3-methylimidazolium dimethylphosphate ([MMIM][DMP]) versus equimolar mixture of 1-methylimidazole (MIM) and trimethylphosphate (TMP). Furthermore, the density enhancement of ILs with respect to their corresponding stoichiometric molecular mixtures can be well represented by the electrostatic contribution among ionic species involved.     

Review of theories on ionization in fast ion-atom collisions with prospects for applications to hadron therapy

This work reviews quantum-mechanical four-body distorted wave theories for double electron capture in collisions between fast heavy multiply charged ions and heliumlike atomic systems. The widely used distorted wave methods of the first- and second-order in the pertinent perturbation series expansions are compared with each other. This tests the presumed importance of double continuum intermediate states of two electrons. Further, the relative performance is evaluated of the second-order theories with and without the eikonalization of the two-electron Coulomb wave functions for double continuum intermediate states. This checks the correctness and usefulness of the eikonalized Coulomb waves when two electrons participate actively to the transition from the initial to the final state of the entire system. We also analyze the significance of the contributions from excited heliumlike states especially in comparison between theory and measurement. The overall goal of the present study is to determine how much of the unprecedented experience gained over several decades in studying high-energy theories of pure three-body charge exchange could be exported directly to four-body double-electron capture without much of additional and essential eleaborations, besides the naturally increased computational demand. In particular, we address the unexpected breakdown of the continuum distorted wave eikonal initial state approximation and the anticipated success of continuum distorted wave theory for double charge exchange in ion-atom collisions at high impact energies.     

醇及脱氧核糖与水分子间三体作用强度的快速准确计算

快速准确预测醇及脱氧核糖分子与水形成的氢键复合物的三体作用强度, 对准确模拟水环境下蛋白质和DNA的结构和功能至关重要. 基于对多体极化作用的理解, 在可极化偶极-偶极作用模型(PBFF)基础上, 将体系中的极性化学键视为化学键偶极, 通过模拟键偶极的极化计算了醇及脱氧核糖与水分子形成的氢键复合物的三体作用能. 通过拟合甲醇与水氢键复合物的三体作用能随分子间距离变化的能量曲线确定了所需的参数. 将模型和所确定的参数应用于计算更多的甲醇、 乙醇及脱氧核糖与水氢键复合物的三体作用能, 检验了模型的准确性和参数的可转移性. 计算结果表明, 可极化偶极-偶极作用模型及所确定的参数能够较好地预测具有不同结构的氢键复合物的三体作用强度, 其精度可与MP2方法的计算精度相当.     

<Emphasis Type="Italic">Ab initio</Emphasis> study of the pressure effect on the structural and electronic properties of crystalline hydrogen azide

Within density functional theory with regard to the dispersion interaction the crystal structure parameters of hydrogen azide are determined. The pressure effect on its structural and electronic properties is studied in the range of 0-10 GPa. By means of the Vinet equation of state the bulk modulus of compression is found to be 9.26 GPa. It is shown that with an increase in the pressure molecules approach each other in molecular layers and this is accompanied by an increase in the total electron density contours, which means the principal possibility for polymerization. The external pressure of 10 GPa leads to the broadening of the upper valence energy bands and a decrease in the band gap from 6.14 eV to 5.51 eV.     

Review of theories on double electron capture in fast ion-atom collisions

This work reviews quantum-mechanical four-body distorted wave theories for double electron capture in collisions between fast heavy multiply charged ions and heliumlike atomic systems. The widely used distorted wave methods of the first- and second-order in the pertinent perturbation series expansions are compared with each other. This tests the presumed importance of double continuum intermediate states of two electrons. Further, the relative performance is evaluated of the second-order theories with and without the eikonalization of the two-electron Coulomb wave functions for double continuum intermediate states. This checks the correctness and usefulness of the eikonalized Coulomb waves when two electrons participate actively to the transition from the initial to the final state of the entire system. We also analyze the significance of the contributions from excited heliumlike states especially in comparison between theory and measurement. The overall goal of the present study is to determine how much of the unprecedented experience gained over several decades in studying high-energy theories of pure three-body charge exchange could be exported directly to four-body double-electron capture without much of additional and essential eleaborations, besides the naturally increased computational demand. In particular, we address the unexpected breakdown of the continuum distorted wave eikonal initial state approximation and the anticipated success of continuum distorted wave theory for double charge exchange in ion-atom collisions at high impact energies.     

Modeling the molecular dynamics of liquid metals at high pressures: Liquid potassium

A way of building potassium models by molecular dynamics using the embedded atom model (EAM) is developed. The contribution from pairwise interaction is presented as power series at an interpartial distance. Embedded potential parameters are determined by the experimental dependence of pressure on volume for a static compression of potassium at 300 K to a pressure of 53 GPa (potential A). By using potential A to describe shock compression and choosing the appropriate temperature at given degree of compression (up to 40000 K at compression to 0.29 of initial volume) it is shown that the model pressure can be made equal to the pressure indicated by the Rankine-Hugoniot relations. The model energy is lower than the actual energy determined by the relations, and the difference in energies increases with temperature almost linearly; such growth corresponds to an excess average heat capacity of about 11.6 J/(mol K), compared to the model heat capacity. It is established that the reasons for this divergence are the inability of the EAM potential to describe the temperature dependency of metal properties precisely, and the appearance of an energy contribution upon heating that is dependent on temperature but not on atom coordinates. Adding another summand to the potential energy (which is dependent on temperature only) allows us to match the heat capacities of real potassium and the models. The dependence of potassium’s melting temperature on pressure is calculated. The calculated melting temperature at 41.2 GPa is 1231 K. Additional data (e.g., the actual temperature on the Rankine-Hugoniot curve and precise quantum mechanics calculations of heat capacity at extreme conditions) is required to eliminate potential ambiguity.     

Group contribution multifluid lattice theory for simultaneous predictions of thermodynamic properties of pure fluids and mixtures

A unified group contribution (GC) lattice equation of state (EOS) was formulated based on the multifluid approximation of the nonrandom lattice fluid theory. The GC-EOS requires segment size and interaction energy parameter from functional group characteristics. The unique feature of the approach is that a single set of group parameters are used for both pure fluids and mixtures. The approach was found to be quantitatively applicable for predicting thermodynamic properties of real pure fluids and mixtures. Its potential utility was demonstrated for vapor pressures, vapor–liquid coexistence densities of pure fluids and phase equilibrium properties of mixtures including polymeric solutions.     

用重整化群理论研究胶体模型体系的相行为

用yukawa势能函数表达胶体颗粒之间的吸引作用。用Duh-Mier-Y-Teran状态方 程表达液相Helmholtz自由能。用一阶微扰理论、固体硬球径向分布函数解析式和 改进的胞腔模型建立固相状态方程，结合建立的状态方程和重整化群理论。研究了 胶体模型体系的液-液相平衡和液-固相平衡。研究表明，颗粒之间色散作用量程参 数的变化对胶休到本世纪末茶杯 系的相行为有特殊需要影响。所得结果与分子模 拟数据吻合良好。     

Conformational and phase transformations of chlorocyclohexane at high pressures by Raman spectroscopy

Pressure induced conformational and phase transformations of chlorocyclohexane (CCH) were investigated in a diamond anvil cell by Raman spectroscopy at room temperature. Pure CCH was compressed up to 20 GPa and then decompressed to ambient pressure. The conformational equilibrium was shifted by pressure from equatorial to axial conformers in the fluid phase below 0.7 GPa, consistent with previous observations. Upon further compression, several solid-to-solid phase transitions were identified by the observation of markedly different Raman patterns as well as different pressure dependences of characteristic Raman modes. The possible structures of these phases were analyzed in correlation with previously observed solid phases at low temperatures. Finally, CCH exhibits pressure hysteresis and partial reversibility upon decompression which result in the formation of the phases with different Raman patterns from those obtained upon compression. The difference can be interpreted as conformational contribution as well as the intrinsic plasticity of CCH crystals.