用原子-键电负性均衡方法计算镉（Ⅱ）配合物的电荷分布

在以密度泛函理论和电负性均衡原理为基础的原子-键电负性均衡方法的σ-π模型（ABEEMσ-π）中,利用最小二乘法并结合自编程序,通过大量的量子化学计算拟合确定了氢、氮、氧、硫、氯以及镉（Ⅱ）等各种类型的原子及相关化学键区域的ABEEMσ-π参数.将这些参数应用到研究镉（Ⅱ）配合物的电荷分布,结果显示由ABEEMσ-π模型计算的电荷分布与量子化学计算结果有很好的一致性,线性相关系数在0.93~0.99之间.进一步验证了ABEEMσ-π模型的合理性和可靠性,发展了该模型在第五周期过渡金属中的应用.     

稠杂环化合物ABEEM参数的拟合及简单应用

应用以密度泛函理论和电负性均衡原理为基础的原子-键电负性均衡方法中的σ-π模型（ABEEMσ-π）,通过大量的量子化学计算,利用最小二乘法,并结合我们自编的程序,拟合确定了稠杂环化合物的ABEEMσ-π各参数.进而,将这些参数应用到研究叶酸及维生素B2等稠杂环药物分子体系的电荷分布,结果显示由ABEEMσ-π模型获得的电荷分布与量子化学方法的计算结果有很好的一致性,线性相关系数均达到0.98以上.这说明我们拟合的参数是正确的,也进一步验证了ABEEMσ-π模型的合理性和可靠性.     

应用ABEEM σπ模型研究含氮杂环化合物的电荷分布

应用在密度泛函理论和电负性均衡原理基础上发展的原子-键电负性均衡方法的σπ模型(即ABEEM σπ模型),并结合自编程序,对5种典型的含氮杂环化合物的电荷分布进行了计算. 计算结果表明,所得的电荷分布可以和从头计算结果有非常好的一致性,进一步验证了ABEEM σπ模型的合理性和可靠性     

应用ab initio方法和ABEEMσπ模型研究含铝金属酶

应用原子-键电负性均衡方法中的σπ模型(ABEEMσπ模型),通过大量量子化学计算,拟合确定了含铝金属酶体系的ABEEMσπ参数.将这些参数应用到含铝金属酶大分子体系的电荷分布及Fukui函数的计算,结果显示,ABEEMσπ模型计算得到的电荷分布及Fukui函数与从头算和实验结论均有很好的一致性.还进一步计算分析了1L3R酶与丝氨酸结合前后的分子各区域的电荷分布,结果表明,Al 3+是1L3R酶的活性中心,根据结合后分子的Fukui函数可以得出丝氨酸会使1L3R酶的活性降低.另外,通过比较两者结合前后Al 3+的广义Fukui函数,证明了广义Fukui函数可用于该体系分子间反应活性的比较,同时也说明利用ABEEMσπ模型来预测含铝金属酶的抑制剂是可行的.     

应用ABEEM_(σπ)模型研究非金属硼化物的电荷分布

通过大量量子化学计算,拟合确定了非金属硼化物体系模型分子的ABEEMσπ参数.将这些参数应用到ABEEMσπ模型中计算非金属硼化物模型分子的电荷分布,计算结果显示,ABEEMσπ模型计算得到的电荷分布与从头算计算的电荷分布都有很好的一致性.还计算了测试分子的电荷分布从而验证了ABEEMσπ参数的可转移性,电荷分布变化的规律是在BF3,BCl3,BBr3,BI3分子中B原子电荷逐渐减小,F,Cl,Br,I原子电荷逐渐增大.ABEEMσπ模型可以很好地应用于非金属硼化物的结构和性质的分析.     

水合碱金属离子Na＋（H2O）n和K＋（H2O）n（n=1～6）体系电荷分布的理论研究

应用量子化学方法和原子-键电负性均衡方法(ABEEM)对Na+(H2O)n和K+(H2O)n(n=1～6)体系进行研究.采用MP2/6-31++G(d,p)方法进行几何构型的优化及频率的计算,在MP2/6-311++G(2d,2p)方法下计算了能量.以HF/STO-3G所计算的体系电荷为基准,拟合确定了Na+和K+的ABEEM参数,应用其计算所得的Na+(H2O)n和K+(H2O)n(n=1～6)的电荷分布,与从头计算的结果一致.     

应用 ABEEMσπ模型研究鸟嘌呤异构体的电荷分布和偶极矩

以密度泛函和电负性均衡原理为基础发展的原子-键电负性均衡方法的σπ模型（ABEECσπ）.并结合自编程序,计算了鸟嘌呤五种异构体的电荷分布和偶极矩,计算结果表明,所计算出的电荷分布可以和从头算结果很好地关联.所计算的偶极矩与其它从头计算方法和力场方法相比总体趋势相同,进一步说明了ABEEMσπ模型的合理性.     

磷脂分子DPPC的浮动电荷分子力场

基于量子化学(QM)计算和原子-键电负性均衡浮动电荷分子力场(ABEEMσπ/MM),构建了1,2-二棕榈酰磷脂酰胆碱(1,2-dipalmitoyl-sn-phosphatidylcholine,DPPC)分子的可极化势能函数.首先将DPPC分子分成乙基三甲基铵(ETMA)、乙酸甲酯(MAS)、磷酸二甲酯(DMP)和烷烃等官能团模型小分子片段,依据这些模型分子的QM结果,优选并确定相关电荷参数及力场参数.应用ABEEMσπ/MM对模型小分子的计算结果与QM计算结果符合很好,其中,稳定结构的键长、键角、二面角的绝对平均偏差(AAD)分别为0.000 5nm、1.59°和1.33°,两者的电荷分布的线性相关系数为0.968 8.进一步将上述函数应用到研究DPPC分子的结构和电荷分布,并与QM结果相比,结果表明:键长、键角、二面角的AAD和均方根偏差(RMSD)分别仅为0.000 01和0.000 02nm、0.006°和0.009°、0.007°和0.011°;两者计算得到的电荷分布的线性相关系数为0.965 9.     

应用ABEEMσπ确定大分子体系的电荷分布

以密度泛函理论（DFT）和电负性均衡原理为基础,明确处理了双键的结构,发展建立原子一键电负性均衡方法中的σπ模型（ABEEMσπ）,本模型将双键划分为一个σ键区域和四个π键区域（每个双键原子各有2个π键区域）,其中,σ电荷中心位于两成键原子之间共价半径之比处;π电荷中心垂直于双键所在平面,置于双键原子上下两侧。本文给出分子中各部分有效电负性的精密公式,以及参数确定方法,并应用该方法简捷快速地计算了C18N3O3H27C15NO2H21,以促甲状腺激素放激素TRH等有机和生物大分子体系的电荷分布,得到的结束可以和相应的从头算项媲美。     

并行程序实现ABEEM σπ模型电荷分布计算

对以密度泛函理论和电负性均衡原理为基础发展的原子-键电负性均衡方法中的σπ模型(ABEEM σπ模型)进行分析得出,利用该模型计算分子体系的电荷分布时,最耗费时间的部分是求解线性方程组.根据解线性方程组的串行程序,我们提出在并行环境下不带平方根的Cholesky分解方法.结果表明,利用改编后的并行程序能快速而准确地计算分子的电荷分布,算法随着矩阵规模的增大,并行效率也随之增高,即分子体系越大结果越理想.因此本算法适用于大规模问题的计算.     

(Urea)n(n=2～5)内分子间相互作用的理论化学研究

用MP2/6-311++G~(**)//B3LYP/6-311++G~(**)方法对(Urea)_n(n=2～5)复合物进行结构优化、频率和结合能的计算,并考虑基组重叠误差.研究发现,(Urea)_n(n=2～5)中存在N—H…O和N—H…N两种类型的氢键,其中,N—H…O型氢键键长较短,氢键角较大,作用强度较大.通过Mulliken电荷布居分析,发现复合物中参与形成氢键的原子所带电荷变化明显,而没有参与形成氢键的原子电荷几乎不变.基于原子-键电负性均衡方法,确定尿素团簇的相关ABEEM电荷参数,计算得到的电荷分布与量子化学结果有很好的一致性.     

应用ABEEMσπ模型研究氨基酸和多肽分子的电荷分布

在密度泛函理论和电负性均衡原理的理论框架下 ,发展了直接用于计算原子及键电荷分布的原子 键电负性均衡方法的σπ模型 ,本文应用此模型计算了氨基酸和一些较大多肽分子的电荷分布 .     

HIV-1病毒中NCp7锌指区域和配体的电荷分布研究

采用B3LYP/6-311++G(d,p)方法,优化了4个锌指蛋白分子和31个模型分子的几何构型,使用HF/STO-3G计算了优化后体系的Mulliken电荷分布.使用线性回归和最小二乘法,拟合确定了ABEEMσπ参数(包括参考电荷,价态电负性和价态硬度).使用ABEEMσπ参数,计算获得了4个锌指蛋白和11个配体的电荷分布.ABEEMσπ方法所计算的电荷分布与从头计算方法的线性相关方程斜率k接近1.000,截距b接近于0.000,线性相关系数R在0.943 4以上.由此可见,ABEEMσπ方法与从头计算方法相比有很好的一致性.验证了笔者所拟合的ABEEMσπ参数是合理的和可转移的,可应用于类似体系的电荷分布的计算.     

应用ABEEM方法计算含金属离子Ga~(3+)蛋白的电荷分布

首次开发ABEEM方法应用于含金属离子Ga~(3+)蛋白体系的研究.ABEEM方法将分子电荷分解到了原子区域、σ键区域、π键区域和孤对电子区域.对金属离子Ga~(3+)与蛋白之间的相互作用采用成键模型,Ga~(3+)与配体原子之间有键电荷分布.通过蛋白晶体数据库的搜索,总结出Ga~(3+)离子和蛋白相互作用的模型分子,确定了相关的新的电荷参数.应用ABEEM方法对模型分子的电荷分布、电荷转移和Ga~(3+)离子的电荷进行了计算和分析.结果表明,ABEEM方法计算的电荷可以和从头算的HF/STO-3G方法的结果相比拟,可以快速给出所有的模型分子的电荷分布.并且通过金属离子Ga~(3+)蛋白大分子体系的电荷计算验证了ABEEM方法以及电荷参数的正确性和可转移性.高价态的Ga~(3+)离子和蛋白的相互作用的理论研究,为其动力学模拟研究奠定了基础.     

应用ABEEMσπ模型确定大分子体系的电荷分布

以密度泛函理论(DFT)和电负性均衡原理为基础，明确处理了双键的结构，发展建立原子一键电负性均衡方法中的σπ模型(ABEFMσπ).本模型将双键划分为一个σ键区域和四个π键区域(每个双键原子各有2个π键区域)，其中，σ电荷中心位于两成键原子之间共价半径之比处；π电荷中心垂直于双键所在平面，置于双键原子上下两侧.本文给出分子中各部分有效电负性的精密公式，以及参数确定方法.并应用该方法简捷快速地计算了C     

应用ABEEM／MM模型研究甘氨酸与水的氢键作用

应用原子-键电负性均衡方法中的浮动电荷分子力场（ABEEM／MM）,我们构建了1个新的甘氨酸（Glycine）-水势能函数,并将其应用到甘氨酸与水组成的二聚体的研究中．首先研究了甘氨酸单体的几何构型,构型结果与可获得的实验数据显示了很好的一致性;根据对甘氨酸中氮原子和氧原子的孤对电子电荷分布的ABEEM／MM方法分析,可以解释水与作用点形成氢键的强弱;进而研究了甘氨酸与1个水作用的结合能,ABEEM／MM模型与从头计算方法得到的结果显示了很好的一致性．表明我们势能函数的合理性以及参数的正确性．     

分子中电荷分布和分子电负性的预测

使用原子－键电负性均衡模型（ＡＢＥＥＭ），计算了一些任意选取的有机大分子的电荷分布和分子电负性。计算结果很好地重复了相应的从头计算结果，但却更加快捷、省时，展示了原子－键电负性均衡模型的广阔应用前景。     

应用ABEEMσπ/MM方法计算小分子肽的总能量

把分子力场(MM)和原子-键电负性均衡方法(ABEEMσπ)协调融合在一起,采用ABEEMσπ/MM方法的分子能量计算方法,选取直链烷烃(n=1~10)、丙酮、丁酮、3-戊酮、N-甲基乙酰胺(NMA)、丙氨酸二肽和若干种氨基酸作为模型分子,计算模型分子能量,从而得到不同类型原子价态能量参数.应用这些原子价态能量参数,计算小分子肽的能量.把所得结果与从头算(MP2/6-311++G(d,p))方法计算的结果相比,得出小分子肽能量的绝对偏差小于9.000 0 kcal/mol,相对偏差小于19.000 0×10-6,且均方根偏差3.450 0 kcal/mol,相对均方根偏差6.919 0×10-6.以上结果表明,ABEEMσπ/MM方法计算的分子能量结果与从头算(MP2/6-311++G(d,p))方法计算的结果具有较好的一致性,此外,用ABEEMσπ/MM方法计算分子能量要远远快于从头算(MP2/6-311++G(d,p))方法.     

含镍超氧化物歧化酶反应活性的原子-键电负性均衡方法预测研究

应用原子-键电负性均衡方法,计算了含镍超氧化物歧化酶活性中心原子的电荷分布和Fukui函数．结果表明,超氧化物歧化酶活性中心与超氧阴离予自由基作用时,金属离子电荷转移约为0．15e,而配体原子的电荷转移却很小;同时金属离子的Fukui函数大于配位原子的Fukui函数．电荷转移和Fukui函数计算结果表明,金属离子是含镍超氧化物歧化酶活性中心,该预测与量子化学理论计算一致．     

用DFT和ABEEMσπ对不同形状单层石墨烯结构研究

选取完整的单层矩形(41个碳环)石墨烯构型和完整单层正六边形(37个碳环)石墨烯构型,采用B3LYP/6-31G(d,p)方法进行几何优化,计算完整石墨烯构型的HOMO和LUMO分布、能隙及电荷分布.研究结果显示,完整的矩形石墨烯的HOMO和LUMO只分布在锯齿形边缘,能隙约0.15 eV;正六边形石墨烯的HOMO和LUMO广泛地分布在整个片层上,能隙约2.0 eV.石墨烯的几何构型包含扶手椅形边缘和锯齿形边缘,然后定义了2种描述C原子的ABEEM参数,即石墨烯边缘碳原子、内部碳原子,以及与这2种碳原子相关的键.使用ABEEM参数计算出的电荷分布与HF/STO-3G方法所计算的电荷分布相比较,得到的线性相关系数都达到了0.99以上.认为所获得的碳原子和π键ABEEM参数是合理的,可以应用于大石墨烯分子体系.