Calculation of molecular energies by atom-bond electronegativity equalization method

The atom-bond electronegativity equalization method (ABEEM), based on the equalization of the “effective electronegativity” of an atom or a bond in a molecule, allows the direct calculation of the charge distribution and the molecular electronegativity of a large molecule. We now demonstrate how the another important quantity, the total molecular energy, can be computed directly and rapidly. It is shown that, based on the ABEEM, the corresponding ab initio values of the total molecular energy can be reproduced with very satisfactory accuracy. In addition, it has been found that in the expression of the energy functional E[ρ] of the ABEEM, the charge-dependent term C correlates quite well with the total molecular bonding energy.     

Molecular electronegativity in density functional theory(VIII)

Based on the density functional theory and the atom-bond electronegativity equalization model (ABEEM), a method is proposed to construct the softness matrix and to obtain the electron population normal modes (PNMs) for a closed system. Using this method the information about the bond charge polarization in a molecule can be obtained easily. The test calculation shows that the PNM obtained by this method includes all the modes about the bond charge polarization explicitly. And the bond charge polarization mode characterized by the biggest eigenvalue, which is the softest one of all modes related with chemical bonds, can describe the charge polarization process in a molecule as exquisitely as the correspondingab initio method.     

Toward a semiempirical density functional theory of chemical binding

The new ideas ofbond electronegativity andbond hardness are introduced, and a semiempirical density functional approach to the theory of molecular electronic structure and chemical binding is outlined. There result effective electronegativity equalization procedures that permit calculation of binding energies as well as partial charges. By a modelling of the bond electronegativity and bond hardness, a density functional interpretation of earlier bond charge models is established. Some numerical results are given for diatomic molecules.Dedicated to Professor J. Koutecký on the occasion of his 65th birthday     

Molecular electronegativity in density functional theory(VIII) Charge polarization modes in a closed system

Based on the density functional theory and the atom-bond electronegativity equalization model (ABEEM), a method is proposed to construct the softness matrix and to obtain the electron population normal modes (PNMs) for a closed system. Using this method the information about the bond charge polarization in a molecule can be obtained easily. The test calculation shows that the PNM obtained by this method includes all the modes about the bond charge polarization explicitly. And the bond charge polarization mode characterized by the biggest eigenvalue, which is the softest one of all modes related with chemical bonds, can describe the charge polarization process in a molecule as exquisitely as the correspondingab initio method.     

Molecular electronegativity in density functional theory(Ⅷ)——Charge polarization modes in a closed system

Based on the density functional theory and the atom-bond electronegativity equalization model (ABEEM), a method is proposed to construct the softness matrix and to obtain the electron population normal modes (PNMs) for a closed system. Using this method the information about the bond charge polarization in a molecule can be obtained easily. The test calculation shows that the PNM obtained by this method includes all the modes about the bond charge polarization explicitly. And the bond charge polarization mode characterized by the biggest eigenvalue, which is the softest one of all modes related with chemical bonds, can describe the charge polarization process in a molecule as exquisitely as the corresponding ab initio method.     

应用ABEEM模型计算铁(Ⅱ)配合物的电荷分布

以密度泛函理论和电负性均衡原理为基础,在原子-键电负性均衡模型中,利用最小二乘法,并结合自编程序,拟合确定了氢、碳、氮、硫以及铁(Ⅱ)等各种类型的原子及相关化学键区域的参数.利用上述参数计算了一些铁(Ⅱ)配合物的电荷分布,计算结果可以和从头算结果很好地相关联.     

分子ABEEM电荷分布的图形表示

利用Matlab软件编写程序, 将基于密度泛函理论和电负性均衡原理发展的原子-键电负性均衡方法σ-π模型(ABEEM σ-π)计算的电荷分布, 特别是键电荷和孤对电子的电荷分布, 用图形表现出来, 对分子的电荷分布给出直观形象的认识, 并以腺嘌呤、腺嘌呤脱氧核苷酸以及丙烯与氯化氢的马氏亲电加成反应为例, 进行应用说明.     

电负性均衡

电负性是分子中一个原子把电子拉向它自身的能力,是化学理论的基本概念之一。继Pauling建立第一个电负性标度后,提出了众多的电负性标度。只是在密度泛函理论的基础上,电负性概念和电负性均衡原理,才被精密地论证。近二十多年来,电负性理论的重要发展是:应用电负性均衡模型或方法,可以快速地计算分子体系的电荷分布,从而确定分子的其他性质,甚至包括分子的结构和反应性指标。通常的电负性均衡方法只把分子划分到原子区域,虽然简单直观,但其精度和应用范围受到限制。原子与键电负性均衡方法,把分子划分到包括原子区域、化学键区域和孤对电子区域,能够较快速精密地计算分子的电荷分布和其他性质,并被应用到构建新一代可极化或浮动电荷力场的探索中,有广阔的应用前景。     

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

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

New Investigation of Potential Acting on an Electron in a Molecule to Draw Molecular Faces

A new approach to calculate the potential acting on an electron in a molecule(PAEM) has been established for drawing the molecular face(MF) of a macromolecule, according to the classic point charge model and the atom-bond electronegativity equalization method(ABEEMσπ) for one electron in a molecule. We introduced a dynamic charge distribution from the view of a local electron movement in a molecule based on the new approach, and as further direct evidence, we calculated some physical quantities using the original ab initio method and the new method to verify the accuracy of the method, such as the boundary distance(BD), molecular face surface area(MFSA) and molecular reactivities indicated by the MFs for a variety of organic molecules. All the results by the new method are in agreement with the results by ab initio method.     

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纯溶液和其水溶液体系进行了分子动力学模拟, 得到的径向分布函数、汽化热和偶极矩等物理量与实验值和其它极化力场方法符合很好, 合理描述了溶质与溶剂之间的静电极化和分子内的电荷转移.     

In-Situ Electronegativity and the Bridging of Chemical Bonding Concepts

One challenge in chemistry is the plethora of often disparate models for rationalizing the electronic structure of molecules. Chemical concepts abound, but their connections are often frail. This work describes a quantum-mechanical framework that enables a combination of ideas from three approaches common for the analysis of chemical bonds: energy decomposition analysis (EDA), quantum chemical topology, and molecular orbital (MO) theory. The glue to our theory is the electron energy density, interpretable as one part electrons and one part electronegativity. We present a three-dimensional analysis of the electron energy density and use it to redefine what constitutes an atom in a molecule. Definitions of atomic partial charge and electronegativity follow in a way that connects these concepts to the total energy of a molecule. The formation of polar bonds is predicted to cause inversion of electronegativity, and a new perspective of bonding in diborane and guanine−cytosine base-pairing is presented. The electronegativity of atoms inside molecules is shown to be predictive of pK_a.     

Charge distribution in homonuclear bonds: A semiempirical modeling

Atoms characterized by nonequivalent electronegativities form chemical bonds by exchanging electrical charge. The fraction of charge exchanged is dictated by the electronegativity differences among the system atoms. In the electronegativity equalization method, the charge distribution is estimated by forcing the system to relax to a common chemical potential, which corresponds to its configuration of energy minimum. By definition, this method cannot be applied to homonuclear bonds. A model is proposed to estimate the charge shared in molecular orbitals of homonuclear molecules. The model expands upon the electronegativity equalization method by adding formalism to describe the spin coupling characteristic of homonuclear bonds. Results are in excellent agreement with other quantum mechanical estimations of the charge distributions. © 2013 Wiley Periodicals, Inc.     

致癌性胺离子反应过程中电荷和化学位移的计算

应用原子-键电负性均衡方法中的σπ模型(ABEEMσπ)计算了由致癌性胺离子所参与的反应过程中的电荷分布,所计算出的电荷分布可以和从头算很好的相关联,并且所需要的时间也大大的缩短;同时应用从头计算程序计算致癌性胺离子反应过程中的NMR化学位移.结果表明,在反应过程中电荷的变化和NMR化学位移的变化有很好的对应关系.     

Development and applications of the ABEEM fluctuating charge molecular force field in the ion-containing systems

The microscopic mechanisms of ion hydration and ion selectivity in biomolecular systems are long-standing research topics,in which the difficulty is how to reasonably and accurately describe the ion-water and ion-biomolecule interactions.This paper summarizes the development and applications of the atom-bond electronegativity equalization fluctuating charge force field model,ABEEM/MM,in the investigations of ion hydration,metalloproteins and ion-DNA bases systems.Based on high-level quantum chemistry calculations,the parameters were optimized and the molecular potential functions were constructed and applied to studies of structures,activities,energetics,and thermodynamic and kinetic properties of these ion-containing systems.The results show that the performance of ABEEM/MM is generally better than that of the common force fields,and its accuracy can reach or approach that of the high-level ab initio MP2 method.These studies provide a solid basis for further investigations of ion selectivity in biomolecular systems,the structures and properties of metalloproteins and other related ion-containing systems.     

Analysis of competing bonding parameters. Part 2. The structure of halosilanes and halogermanes (MH4−nXn,n=1–4; M=Si,Ge; X=F,Cl, Br)

A qualitative rationalization of bonding patterns in halosilanes and halogermanes (MH_4−nX_n, n=1–4; M=Si, Ge; X=F, Cl, Br) is presented. Geometrical and bonding properties in these molecules are discussed on the basis of ab initio molecular orbital calculations employing the natural bond orbital population analysis. The results have been compared with data derived previously for halomethanes. Differences in the n-dependence of the M–Cl and M–Br bond lengths for M=C, Si, Ge are explained by a significant reduction in the closed-shell repulsion between the halide atoms. As M gets larger, a continuous decrease in the X–M–X bond angle is observed. Small bond angles (for n=2, 3) are favoured by the p-rich M orbitals in the M–X bonds. They are opposed, however, by the X⋯X repulsion. As M gets larger, the X⋯X separation for a given bond angle increases. A reduction in the X–M–X bond angle is therefore accomplished without overcompensation due to the X⋯X repulsion energy. The variation in the charge density at M as a function of n has been rationalized by differences in the electronegativity of the terminal atoms H and X. Dipole moments have been computed for the molecules in the series. As in the fluoromethanes, a maximum in the dipole moments at n=2 is explained by a combination of geometric and electronic properties unique to the fluoro-compounds. These are, an n-independent charge density at the F sites and a significant decrease in the M–F bond distance as n increases.     

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)对五种环多肽晶体进行了研究. 与传统力场相比, 该方法中的静电势包含了分子内和分子间的静电极化作用, 以及分子内电荷转移影响, 同时加入了化学键等非原子中心电荷位点, 合理地体现了分子中的电荷分布. 相对其他极化力场模型, 具有计算量较小的特点. 该模型下计算得到的环多肽分子单元相对实验测得的结构的原子位置、氢键长度和二面角的均方根偏差分别为0.009 nm、0.013 nm和5.16°, 能够很好地重复实验结果. 总体上, 其结果优于或相当于其他力场模型, 适用于对实际蛋白质体系的模拟和研究.