电负性均衡

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

原子-键电负性均衡方法中的σπ模型及应用

带有双键或共轭双键的化合物可进行加成、氧化和聚合等反应,在有机、生物和制药等领域用途广泛．密度泛函理论下的电负性均衡方法在研讨分子电荷分布和反应性等方面有独特优势［１～４］．在电负性均衡方法中,Ｍｏｒｔｉｅｒ［１］的电负性均衡方法引人注目,但其未考虑分子中化学键的存在．Ｇｈｏｓｈ［２］的半经验电负性均衡方法考虑了化学键电荷,但键电荷取值假定太简单,只能处理双原子分子．Ｙａｎｇ［３,４］等的原子－键电负性均衡方法同时考虑了分子中原子和化学键的存在,取得令人满意的结果．至今尚无人明确地考虑双键的结构…     

确定大分子体系电荷分布的新方法

近年来，应用密度泛函理论探求电负性和硬度与体系中电荷分布之间的关系，成为令人关注的课题[1~5]．Mortier[1,2]的电负性均衡方法（EEM）没有考虑化学键电荷．Ghosh[4]的半经验电负性均衡方法考虑了化学键电荷，但键电荷取值假设太简单．本文同时考虑了分子中的原子电荷和键电荷的作用，给出了分子中原子和化学键有效电负性的精密公式，这些公式为发展一个系统的精密电负性均衡方法及应用奠定了基础．1理论方法以密度泛函理论[6]为基础，将分子的单电子密度ρmol（r）按式（1）分割：式中ρα（r）是分子中α原子区域的单电子密度，ρα…     

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

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

原子-键电负性均衡方法参数的确定与探讨

原子-键电负性均衡方法(ABEEM)是以密度泛函理论(DFT)和电负性均衡原理为基础发展而来，它明确地考虑了化学键是不引入任何实验数据的带纯理论性和计算的方法．使用统一标准并具有代表性和全面性地选择了200多个模型分子，利用可得到较准确结构的MP2／6-31G*优化结构，心／STO-3G单点计算得到Mulliken重叠布居，再用最小二乘法拟合得到许多主族元素在分子体系中的诸原子(包括单、双和叁键等不同成键状态)和化学键的ABEEM参数．所得到的原子的价态电负性可与已提出的其他电负性标度相比拟，计算CO得到的电荷负端为C(与从头计算的结果相反)，结果与实验相符，且原子电荷的正负不完全由原子电负性决定．     

密度泛函理论下的分子电负性——Ⅵ.原子键电负性均衡模型

以密度泛函理论和电负性均衡原理为基础 ,将体系的单电子密度分割为原子的单电子密度和键的单电子密度 ,得到了分子总能量、分子中原子以及键的有效电负性的表达式 .基于这些表达式 ,提出了直接计算体系的总能量和体系中电荷分布的新方法 .此模型比其他电负性均衡方法更合理 .对大量分子的总能量和电荷分布所进行的计算结果表明所提出的原子 键电负性均衡模型是可行的 ,可以应用于确定各类生物和有机大分子的总能量和电荷分布 ,给出较EEM和MEEM模型更近于从头计算的结果     

密度泛函理论下的分子电负性--Ⅷ. ABEEM方法中封闭体系的电荷极化模式

在密度泛函理论（ＤＦＴ）和原子－键电负性均衡模型（ＡＢＥＥＭ）基础上,给出了建立ＡＢＥＥＭ软度矩阵和封闭体系的电子布居正则模式（ＰＮＭ）的方法．该方法可用于计算和讨论分子体系内部的电荷极化过程,计算结果表明,考虑了孤对电子（ｌｐ）后,原子－键电负性均衡模型下的电子布居正则模式明确地指出了分子内各个部分电荷间的极化．从ＡＢＥＥＭ可以得到与所有化学键电荷极化有关的电子布居正则模式,而这些与化学键电荷极化有关的电子布居正则模式中的最软模式可以正确地预测相应的从头计算结果,但是,方法更简单,结果更明了。     

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

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

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

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

环多肽晶体的浮动电荷极化力场模拟

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

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.     

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(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.     

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.     

Local HSAB rationalization of diels–alder reactions by means of ab initio and ABEEMσπ methods: Stereoselectivity and reaction rate

In this article, four series of Diels–Alder (DA) reactions (totally, 68 reactions) have been investigated, including their transition states, by means of ab initio method at MP2/6‐311+G(d,p) level and atom‐bond electronegativity equalization method (ABEEMσπ) method. The rationalization of stereoselectivities of these DA reactions has been done in the light of the local hard and soft acids and bases (HSAB) criterion. The results show that the finite difference approximation with ab initio method can not always be used to predict the stereoselectivities of these reactions, while ABEEMσπ method can properly be applied to predict the stereoselectivities of these DA reactions. Moreover, we have proposed the generalized Fukui function and local softness which involve the number of the atoms of a molecule. For the first time, we here demonstrate that local HSAB criterion with the generalized local softness can be utilized not only to predict the main products of these DA reactions but also to rationalize the relative magnitudes of the reaction rate constants of these DA reactions. © 2013 Wiley Periodicals, Inc.     

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

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

Valence state parameters of all transition metal atoms in metalloproteins—development of ABEEMσπ fluctuating charge force field

To promote accuracy of the atom‐bond electronegativity equalization method (ABEEMσπ) fluctuating charge polarizable force fields, and extend it to include all transition metal atoms, a new parameter, the reference charge is set up in the expression of the total energy potential function. We select over 700 model molecules most of which model metalloprotein molecules that come from Protein Data Bank. We set reference charges for different apparent valence states of transition metals and calibrate the parameters of reference charges, valence state electronegativities, and valence state hardnesses for ABEEMσπ through linear regression and least square method. These parameters can be used to calculate charge distributions of metalloproteins containing transition metal atoms (Sc‐Zn, Y‐Cd, and Lu‐Hg). Compared the results of ABEEMσπ charge distributions with those obtained by ab initio method, the quite good linear correlations of the two kinds of charge distributions are shown. The reason why the STO‐3G basis set in Mulliken population analysis for the parameter calibration is specially explained in detail. Furthermore, ABEEMσπ method can also quickly and quite accurately calculate dipole moments of molecules. Molecular dynamics optimizations of five metalloproteins as the examples show that their structures obtained by ABEEMσπ fluctuating charge polarizable force field are very close to the structures optimized by the ab initio MP2/6–311G method. This means that the ABEEMσπ/MM can now be applied to molecular dynamics simulations of systems that contain metalloproteins with good accuracy. © 2014 Wiley Periodicals, Inc.     

从头算和ABEEMσπ/MM对(CH_3OH)_n(n=3～12)和[Na(CH_3OH)_n]~+(n=3～6)体系的研究

采用从头算方法(ab initio)和原子-键电负性均衡浮动电荷分子力场方法(ABEEMσπ/MM),对甲醇团簇(CH_3OH)n(n=3~12)和[Na(CH_3OH)_n]~+(n=3~6)体系的结构、电荷分布和结合能进行研究.依据从头算结果构建上述体系的ABEEMσπ/MM浮动电荷势能函数,并确定相关参数.结果表明,ABEEMσπ/MM所获得的结构和结合能等均优于OPLS/AA力场,并与从头算结果相符,其中键长的平均绝对偏差(AAD)小于0.004 nm,键长、键角和结合能的相对均方根偏差(RRMSD)分别小于3.8%,1.7%和6.8%;电荷分布与从头算结果的线性相关系数均大于0.99.