Theoretical study on the hydration of hydrogen peroxide in terms of ab initio method and atom-bond electronegativity equalization method fused into molecular mechanics

In this paper, the interaction between hydrogen peroxide (HP) and water were systemically studied by atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM) and ab initio method. The results show that the optimized geometries, interaction energies and dipole moments of hydrated HP clusters HP(H₂O) _n (n = 1–6) calculated by ABEEM/MM model are fairly consistent with the MP2/aug-cc-pVTZ//MP2/aug-cc-pVDZ results. The ABEEM/MM results indicate that n = 4 is the transition state structure from 2D planar structure to 3D network structure. The variations of the average hydrogen bond length with the increasing number of water molecules given by ABEEM/MM model agree well with those of ab initio studies. Moreover, the radial distribution functions (RDFs) of water molecule around HP in HP aqueous solution have been analyzed in detail. It can be confirmed that HP is a good proton donor and poor proton acceptor in aqueous solution by analysis of the RDFs.     

Application of the ABEEM/MM model in studying the properties of the water clusters (H2O)n (n= 7- 1O)

ABEEM/MM model has been applied to compute the various properties characterizing water clusters (H2O)n(n = 7-10), such as optimized geometries, the hydrogen bonds number, cluster interaction energies, stabilities, ABEEM charge distributions, dipole moments, structural parameters, and so on, and to describe the transition reflected by the hexamer region from two-dimensional (from dimer to pentamer) to three-dimensional structures (for clusters larger than the hexamer).     

Application of the ABEEM/MM model in studying the properties of the water clusters(H_2O)_n(n=7-10)

ABEEM/MM model has been applied to compute the various properties characterizing water clusters(H2O) n(n = 7-10) ,such as optimized geometries,the hydrogen bonds number,cluster interaction en-ergies,stabilities,ABEEM charge distributions,dipole moments,structural parameters,and so on,and to describe the transition reflected by the hexamer region from two-dimensional(from dimer to pen-tamer) to three-dimensional structures(for clusters larger than the hexamer) .     

应用ABEEM/MM模型研究水分子团簇(H2O)n (n=11～16)的性质

应用ABEEM/MM 模型计算了较大的水分子团簇(H2O)n (n=11~16)的各种性质,如：优化的几何构型, 氢键个数, 结合能, 稳定性, ABEEM 电荷分布, 偶极矩, 以及结构参数、平均氢键个数和强度, 增加的团簇结合能等.结果表明,从立方体结构到笼状结构的过渡出现在n=12的水分子团簇中,随着类似于笼状结构特点的不断增强,五元环的富集程度有所增加.     

Application of the ABEEM/MM model in studying the properties of the water clusters (H2O) n (n=7−10)

ABEEM/MM model has been applied to compute the various properties characterizing water clusters (H₂O) _n (n = 7?10), such as optimized geometries, the hydrogen bonds number, cluster interaction energies, stabilities, ABEEM charge distributions, dipole moments, structural parameters, and so on, and to describe the transition reflected by the hexamer region from two-dimensional (from dimer to pentamer) to three-dimensional structures (for clusters larger than the hexamer).     

Application of the ABEEM/MM model in studying the properties of the water clusters (H2O) n ( n =7?10)

ABEEM/MM model has been applied to compute the various properties characterizing water clusters (H₂O) _n (n = 7−10), such as optimized geometries, the hydrogen bonds number, cluster interaction energies, stabilities, ABEEM charge distributions, dipole moments, structural parameters, and so on, and to describe the transition reflected by the hexamer region from two-dimensional (from dimer to pentamer) to three-dimensional structures (for clusters larger than the hexamer). Supported by the National Natural Science Foundation of China (Grant No. 20373021)     

Sr2+和Ba2+水分子体系结构和结合能的从头算和ABEEM/MM研究

应用从头算方法和ABEEM/MM浮动电荷分子力场, 研究了水合碱土离子团簇Sr²⁺/Ba²⁺(H₂O)_n (n=1-6), 构建了离子-水相互作用的ABEEM/MM势能函数, 获得了水合离子团簇的稳定结构, 计算了结合能. 计算结果表明, ABEEM/MM方法的结果和从头算方法的结果有很好的一致性. 进一步应用ABEEM/MM对Sr²⁺和Ba²⁺水溶液进行了分子动力学模拟. 对Sr²⁺水溶液, 得到的Sr²⁺-水中氧原子的径向分布函数的第一和第二最高峰分别位于0.257和0.464 nm处, 第一和第二水合层的配位水分子数分别为9.2和11.4; 对Ba²⁺水溶液, 得到的Ba²⁺与水中氧原子的径向分布函数的第一和第二最高峰分别位于0.269和0.467 nm处, 第一和第二水合层的配位水分子数分别为9.9和12.4. 这与实验值或其它理论模拟结果有较好的一致性. 对比外层的水分子, 金属离子的极化作用使得溶液中第一水合层中水分子的O―H键长增长, HOH键角减小.     

Study of lithium cation in water clusters: based on atom-bond electronegativity equalization method fused into molecular mechanics

We present a potential model for Li(+)-water clusters based on a combination of the atom-bond electronegativity equalization and molecular mechanics (ABEEM/MM) that is to take ABEEM charges of the cation and all atoms, bonds, and lone pairs of water molecules into the intermolecular electrostatic interaction term in molecular mechanics. The model allows point charges on cationic site and seven sites of an ABEEM-7P water molecule to fluctuate responding to the cluster geometry. The water molecules in the first sphere of Li(+) are strongly structured and there is obvious charge transfer between the cation and the water molecules; therefore, the charge constraint on the ionic cluster includes the charged constraint on the Li(+) and the first-shell water molecules and the charge neutrality constraint on each water molecule in the external hydration shells. The newly constructed potential model based on ABEEM/MM is first applied to ionic clusters and reproduces gas-phase state properties of Li(+)(H(2)O)(n) (n = 1-6 and 8) including optimized geometries, ABEEM charges, binding energies, frequencies, and so on, which are in fair agreement with those measured by available experiments and calculated by ab initio methods. Prospects and benefits introduced by this potential model are pointed out.     

A study of N-methylacetamide in water clusters: based on atom-bond electronegativity equalization method fused into molecular mechanics

N-methylacetamide (NMA) is a very interesting compound and often serves as a model of the peptide bond. The interaction between NMA and water provides a convenient prototype for the solvation of the peptides in aqueous solutions. Here we present NMA-water potential model based on atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM) that is to take ABEEM charges of all atoms, bonds, and lone-pair electrons of NMA and water molecules into the electrostatic interaction term in molecular mechanics. The model has the following characters: (1)it allows the charges in system to fluctuate responding to the ambient environment; (2) for two major types of intermolecular hydrogen bonds, which are the hydrogen bond forming between the lone-pair electron on amide oxygen and the water hydrogen, and the one forming between the lone-pair electron on water oxygen and the amide hydrogen, we take special treatments in describing the electrostatic interaction by the use of the parameters k(lpO=, H) and k(lpO(-), HN(-)), respectively. The newly constructed potential model based on ABEEM/MM is first applied to amide-water clusters and reproduces gas-phase state properties of NMA(H(2)O)(n) (n=1-3) including optimal structures, dipole moments, ABEEM charge distributions, energy difference of the isolated trans- and cis-NMA, interaction energies, hydrogen bonding cooperative effects, and so on, whose results show the good agreement with those measured by available experiments and calculated by ab initio methods. In order to further test the reasonableness of this model and the correctness and transferability of the parameters, many static properties of the larger NMA-water complexes NMA(H(2)O)(n) (n=4-6) are also studied including optimal structures and interaction energies. The results also show fair consistency with those of our quantum chemistry calculations.     

ABEEMσπ/MM力场模型下的多肽构象

在ABEEM/MM蛋白质浮动电荷力场模型的基础上,加入孤对电子和 电荷位点,从而能够体现多肽和蛋白质分子中一些重要的各向异性极化性质,允许非化学键方向的电子转移和极化。利用从头计算数据拟合模型相关参数。计算得到的小分子团簇结合能、偶极矩、氢键键长等性质与从头计算结果符合很好。该经典极化模型力场能够重复量子场下丙氨酸二肽、丙氨酸四肽、甘氨酸三肽的各稳定构象,其稳定性顺序与精密从头计算结果相一致,其结构和能量性质较以往模型有一定提高,并优于其他力场模型。     

Study of water clusters in the n = 2–34 size regime, based on the ABEEM/MM model

Various properties of water clusters in the n = 2–34 size regime with the change of cluster size have been systemically explored based on the newly developed flexible-body and charge-fluctuating ABEEM/MM water potential model. The ABEEM/MM water model is to take ABEEM charges of all atoms, bonds, and lone-pairs of water molecules into the intermolecular electrostatic interaction term in molecular mechanics. The computed correlating properties characterizing water clusters (H₂O) _n (n = 2–34) include optimal structures, structural parameters, ABEEM charge distributions, binding energies, hydrogen bonds, dipole moments, and so on. The study of optimal structures shows that the ABEEM/MM model can correctly predict the following important structural features, such as the transition from two-dimensional (from dimer to pentamer) to three-dimensional (for clusters larger than the hexamer) structures at hexamer region, the transition from cubes to cages at dodecamer (H₂O)₁₂, the transition from all-surface (all water molecules on the surface of the cluster) to one water-centered (one water molecule at the center of the cluster, fully solvated) structures at (H₂O)₁₇, the transition from one to two internal molecules in the cage at (H₂O)₃₃, and so on. The first three structural transitions are in good agreement with those obtained from previous work, while the fourth transition is different from that identified by Hartke. Subsequently, a systematic investigation of structural parameters, ABEEM charges, energetic properties, and dipole moments of water clusters with increasing cluster size can provide important reference for describing the objective trait of hydrogen bonds in water cluster system, and also provide a strong impetus toward understanding how the water clusters approach the bulk limit.     

水合氢离子团簇从头算和ABEEM/MM的理论研究

应用高水平的从头计算方法和ABEEM/MM模型, 研究了水合氢离子团簇H3O＋(H2O)n (n＝1～6), 优化得到了低能构象, 探讨了其结合能和稳定性, 显示出H3O＋(H2O)3局域结构的优势存在．对H3O＋(H2O)6VIa团簇的ABEEM电荷分布进行分析, 表明第一水合层水分子与水合氢离子之间的氢键相互作用要明显强于与第二水层水分子的氢键相互作用. 研究结果表明, ABEEM/MM方法计算的结果和从头算得到的结果存在很好的一致性.     

精密从头算与ABEEM/MM模型对水团簇(H2O)11 9种低能结构的计算

用精密从头算方法研究了(H2O)11的9种低能异构体的性质, 包括优化的几何结构、结合能、偶极矩和氢键个数等, 并且得出了515-a是(H2O)11的全局最低能结构. 同时, 也用ABEEM/MM(atom bond electronegativity equalization method/molecular mechanics) 模型研究了这些性质, 与从头算的结果进行了比较, 得到了相符合的结果. 这显示了ABEEM/MM模型在描述中等大小的水分子团簇结构上是成功的.     

Theoretical Study on Co^3＋ in Aqueous Solution in Terms of ABEEM/MM Model

A detailed theoretical investigation on Co^3＋ hydration in aqueous solution has been carded out by means of molecular dynamics （MD） simulations based on the atom-bond electronegativity equalization method fused into molecular mechanics （ABEEM/MM）. The effective Co^3＋ ion-water potential has been constructed by fitting to ab initio structures and binding energies for ionic clusters. And then the ion-water interaction potential was applied in combination with the ABEEM-7P water model to molecular dynamics simulations of single Co^3＋（aq.） solution, managing to reproduce many experimental structural and dynamical properties of the solution. Here, not only the common properties （radial distribution function, angular distribution function and solvation energy） obtained for Co^3＋ in ABEEM-7P water solution were in good agreement with those from the experimental methods and other molecular dynamics simulations but also very interesting properties of charge distributions, geometries of water molecules, hydrogen bond, diffusion coefficients, vibrational spectra are investigated by ABEEM/MM model.     

原子-键电负性均衡方法融合进分子力场(ABEEM/MM)应用于蛋白质分子(Crambin)模拟

建立应用于多肽和蛋白质模拟的ABEEM/MM浮动电荷力场.利用该模型和参数,对实际蛋白质分子Crambin(植物种子中的一种小的蛋白质)进行模拟,得到了满意的结果,为其更广泛的应用开辟了道路.     

Ab initio and ABEEM/MM fluctuating charge model studies of dimethyl phosphate anion in a microhydrated environment

Dimethyl phosphate (DMP) anion has been used extensively as a model compound to simulate the properties of phosphate group. A 35-point DMP anion potential model is constructed based on the atom-bond electronegativity equalization fluctuating charge molecular force field (ABEEM/MM), and it is employed to study the properties of gas-phase DMP anion and DMP-(H₂O) _n (n = 1–3) clusters. The ABEEM/MM model reproduces well the properties obtained by available experiments and QM calculations, including charge distributions, geometries, and conformational energies of gas-phase DMP-water complexes. Furthermore, molecular dynamics simulation on the DMP anion in aqueous solution based on the ABEEM/MM shows that a remarkable first hydration shell around the nonesterified oxygen atom of DMP anion is formed with a coordination number of 5.2. It is also found that two hydrogen atoms of one water molecule form two hydrogen bonds with two nonesterified oxygen atoms of DMP anion simultaneously. This work could be used as a starting point for us to establish the ABEEM/MM nucleic acid force field.     

利用浮动电荷力场对N-甲基乙酰胺水溶液的分子动力学模拟

利用原子键电负性均衡结合分子力场方法(ABEEM/MM)对N-甲基乙酰胺(NMA)分子的水溶液体系进行了分子动力学模拟. 与经典的力场模型相比, 该方法中的静电势包含了分子内和分子间的静电极化作用, 以及分子内电荷转移影响, 同时加入了化学键等非原子中心电荷位点, 合理体现了分子中的电荷分布. 相对其它极化力场模型, 该模型具有计算量较小的特点. 在该模型下对NMA纯溶液和其水溶液体系进行了分子动力学模拟, 得到的径向分布函数、汽化热和偶极矩等物理量与实验值和其它极化力场方法符合很好, 合理描述了溶质与溶剂之间的静电极化和分子内的电荷转移.