电负性研究I．建立元素电负性标度的有效核电荷数方法

通过价层电离能、价键轨道能量用有效核电荷数法建立了周期表中90种元素的电负性新标度。X＝0。4123√－EV,该式表明电负性值与价键轨道能量的绝对值的平方根成正比,所得数值是一套无量纲的相对参数。元素电负性值随价态的升高与元素非金属性的增强相对应,元素电负性的大小不仅与单个成键电子腾,而且也与参加价键作用的多个电子甚至整个价层都有紧密的联系。氢的元素电负性值不同于Pauling值、等于1．52。用16种氢化物中键的额外离子能△′对（XA－XB）^2作图,两者之间确实具有良好的线性关系。本方法充分体现了目前公认的三大电负性标度的优点,该标度同时也是价层电子在价键状态下的一种能量标度,是对元素周期律的定量描述和反映。     

电负性研究Ⅰ.建立元素电负性标度的有效核电荷数方法

通过价层电离能、价键轨道能量用有效核电荷数法建立了周期表中 90种元素的电负性新标度。χ=0 .41 2 3 -EV,该式表明电负性值与价键轨道能量的绝对值的平方根成正比 ,所得数值是一套无量纲的相对参数。元素电负性值随价态的升高与元素非金属性的增强相对应 ,元素电负性的大小不仅与单个成键电子有关 ,而且也与参加价键作用的多个电子甚至整个价层都有紧密的联系。氢的元素电负性值不同于Pauling值、等于 1 .52。用 1 6种氢化物中键的额外离子能Δ′对 (χA-χB) 2 作图 ,两者之间确实具有良好的线性关系。本方法充分体现了目前公认的三大电负性标度的优点 ,该标度同时也是价层电子在价键状态下的一种能量标度 ,是对元素周期律的定量描述和反映。     

从头算VB-MP2组合方法

提出了一种经济实用的价键计算方法VB-MP2方法.将体系的电子分为非活性电子和活性电子,应用MP2方法计算非活性电子的相关能效应,用VB方法处理活性电子.测试计算表明,该方法保持了价键方法的特点,且计算结果比传统的使用芯-价分离技术的价键方法有较大的改善.     

氯化钕异丙醇配合物的电子结构与络合活化的研究

本文用INDO方法计算了氯化钕异丙醇配合物的电子结构。计算结果指出, 配位后, 与Nd-Cl键对应的分子轨道能量升高, Nd-Cl键级降低, Nd的正电荷降低, Cl的负电荷增多。从而表明。由于配位发生轨道重迭, 使得异丙醇氧的孤对电子移向钕离子的空价轨道。钕离子上的电子云密度增高加强了对Nd—Cl键的价电子的推斥作用, 使Nd-Cl间的价电子密度减小, 导致Nd-Cl键的削弱, 降低了键的能量, 有效地活化了Nd-Cl键,有利于生成活性物种的烷基化反应的进行。同时表明, 对Nd-Cl的络合催化有某种选择性。     

过渡金属羰基簇合物的键价计算和分析

本文运用原子簇化合物键价计算公式 ,对过渡金属羰基簇合物成键情况进行了分析 ,利用金属键轨道数 ,价非键轨道数和金属配体成键轨道数计算簇合物价轨道总数 .计算结果表明 :簇合物价轨道总数与金属键轨道数成线性关系 ,BT=9N -Bn.对于一般簇合物其价轨道总数与Lauher的EHMO计算结果 ,与唐敖庆的结构拓扑规则一致 ,对于反常的高核簇合物其价轨道总数与按化学式计算的 1/ 2VE相吻合     

关于构造合理非正交定域轨道问题的研究

与正交定域轨道相比, 非正交定域轨道更符合经典化学图像、更加紧缩从而更具有可移植性和更符合电子相关能计算和线性标度计算的要求. 非正交定域轨道应该尽可能紧缩而形状及空间分布又与传统化学图像一致. 研究结果表明, 非正交定域轨道组的延伸度是它的正交程度的单调下降函数, 随其正交程度的降低可以缩减到任意程度, 但在超过一定限度以后, 轨道的形状和空间分布不再符合传统化学图像, 最后趋向线性相关. 在没有正交条件约束下, 必须加上另外的限制条件, 优化轨道延伸度泛函才能获得合理的非正交定域轨道. 提出了一种在最小化轨道延伸度泛函的同时最大化轨道重心间距离的方法, 可以得到基本上符合要求的非正交定域轨道.     

化学反应的价键理论研究

根据价键波函数的对称性及基态、激发态化学键与键表的对应 ,可以把轨道对称守恒原理和键对称规律表述成“价键结构对称性匹配” .价键结构对称规则通过对反应物和产物的价键结构进行对称性分析 ,可以直接预测化学反应活性 ,无需依赖理论计算结果 .在价键结构对称性分析的基础上 ,应用多VB结构计算方法 ,建立了反应H HLi→H2 Li和H LiH→HLiH(cyclic)→HLi H的曲线交叉VB图 ,讨论了能垒及过渡态的形成机制 .计算结果表明 ,H原子交换反应H LiH HLi H是一个两步过程 .     

最佳定域分子轨道的新算法

本文提出了最佳群对称定域分子轨道（ＯＳＬＭＯＳ）的概念，报道了ＯＳＬＭＯｓ的生成方法。ＯＳＬＭＯｓ满足分子对称、等价与正交；同时最佳逅近经典的非正交价键型轨道，并可作为ＭＣＳＣＦ与ＣＩ方法中的单电子轨道。     

价键理论中的组态相互作用

近 2 0年来 ,从头计算水平的价键 (VB)方法得到了人们的重视 ,并广泛应用于化学反应等问题的研究[1～ 5] ,然而目前价键理论的计算方法仍然很不完善 .用 VBSCF方法进行计算虽然比较简单 ,能正确地描述化学反应的形成机理 ,但数值结果不理想 ;而用 BOVB方法[4 ] 进行计算虽然可以得到较好的计算结果 ,但存在收敛困难等问题 .分子轨道理论中的组态相互作用是一种简单直接的电子相关能计算方法 ,显然这一方法可以应用于价键方法中 .然而与分子轨道理论方法不同 ,在价键方法中 ,无法直接得到空轨道 ,此外如何选取激发价键函数使得计算结果…     

SN2反应Xl^-＋CH3Xr→XlCH3＋Xr^-（X1=Xr=F，Cl，Br，I）的价键方法研究

应用最近发展的价键组态相互作用(VBCI)方法计算了SN2反应Xl^- CH3Xr→XlCH3 Xr^-(X1=Xr=F，Cl，Br，I)的反应能垒和价键相关参数．计算结果表明．VBCI能垒与采用分子轨道理论的CCSD(T)方法计算的能垒相一致．讨沦了SN2反应的反应参数、     

A new algorithm for inactive orbital optimization in valence bond theory

This paper presents an efficient algorithm for energy gradients in valence bond self-consistent field (VBSCF) method with non-orthogonal orbitals. The frozen core approximation method is extended to the case of non-orthogonal orbitals. The expressions for the total energy and its gradients are presented by introducing auxiliary orbitals, where inactive orbitals are orthogonal, while active orbitals are non-orthogonal themselves but orthogonal to inactive orbitals. It is shown that our new algorithm has a low scaling of (N _a + 1)m ⁴, where N _a and m are the numbers of the active orbitals and basis functions, respectively, and is more efficient than the existing VBSCF algorithms.     

An efficient algorithm for complete active space valence bond self‐consistent field calculation

This article describes a novel algorithm for the optimization of valence bond self‐consistent field (VBSCF) wave function for a complete active space (CAS), so‐called VBSCF(CAS). This was achieved by applying the strategies adopted in the optimization of CASSCF wave functions to VBSCF(CAS) wave functions, using an auxiliary orthogonal orbital set that generates the same configuration space as the original nonorthogonal orbital set. Theoretical analyses and test calculations show that the VBSCF(CAS) method shares the same computational scaling as CASSCF. The test calculations show the current capability of VBSCF method, which involves millions of VB structures. © 2012 Wiley Periodicals, Inc.     

An efficient algorithm for energy gradients and orbital optimization in valence bond theory

An efficient algorithm for energy gradients in valence bond theory with nonorthogonal orbitals is presented. A general Hartree-Fock-like expression for the Hamiltonian matrix element between valence bond (VB) determinants is derived by introducing a transition density matrix. Analytical expressions for the energy gradients with respect to the orbital coefficients are obtained explicitly, whose scaling for computational cost is m(4), where m is the number of basis functions, and is thus approximately the same as in HF method. Compared with other existing approaches, the present algorithm has lower scaling, and thus is much more efficient. Furthermore, the expression for the energy gradient with respect to the nuclear coordinates is also presented, and it provides an effective algorithm for the geometry optimization and the evaluation of various molecular properties in VB theory. Test applications show that our new algorithm runs faster than other methods.     

Evaluation of the charge penetration energy between non-orthogonal molecular orbitals using the Spherical Gaussian Overlap approximation

An overlap dependent formula for evaluating the charge penetration energy between non-orthogonal molecular orbitals is derived using the Spherical Gaussian Overlap approximation. When combined with an accurate multipole representation of the electrostatic energy, such as in the effective fragment potential method, ab initio electrostatic energies are generally reproduced to within 0.2 kcal/mol for a variety of molecular dimers and basis sets. The only larger error is for the DMSO dimer, where the electrostatic energy is overestimated by 0.7 kcal/mol.     

The use of population localised orbitals to interpret molecular orbital calculations

An efficient and general method is derived to calculate population localised molecular orbitals (LMO's) from a given SCF eigenvector matrix, by reduction to an eigenvalue problem. Applications to both localised molecules (NH₃ and C₂H₂) and delocalised ones (B₂H₆, C₆H₆ and butadiene) are discussed in some detail. It is shown that unequal occupation of atomic energy levels leads to non-orthogonal LMO's. The consequences of non-orthogonal atomic hybrid orbitals are discussed, formulas for their overlap in terms of atomic occupation numbers are derived and it is shown that the occupation numbers are connected to LMO atomic orbital coefficients by various sum rules.     

价键理论新进展

概要介绍了现代价键理论的几个主要方法,并讨论了它们各自的特点及其发展现状,并重点介绍了键表方法的基本理论、计算程序及一些应用。     

The HeHe Coulomb and exchange interaction energy

The Coulomb and exchange interaction energy between two helium atoms is calculated, increasing systematically the basis set, by a supermolecular method which excludes the superposition error by using non-orthogonal orbitals. At 5.6 bohr an energy limit of 9.69 K is found at the Hartree-Fock level; for correlated wavefunctions the corresponding value is 11.86 K using only s orbitals, 10.70 K including p orbitals and 10.58 K adding d orbitals.     

MC SCF molecular gradients and hessians: Computational aspects

Molecular gradients and hessians for multiconfigurational self-consistent-field wavefunctions are derived in terms of the generators of the unitary group using exponential unitary operators to describe the response of the energy to a geometrical deformation. Final expressions are cast in forms which contain reference only to the primitive non-orthogonal atomic basis set and to the final orthonormal molecular orbitals; all reference to intermediate orthogonalized orbitals is removed. All of the deformation-dependent terms in the working equations reside in the one- and two-electron integral derivatives involving the atomic basis orbitals. The deformation-independent terms, whose contributions can be partially summed, involve symmetrized density matrix elements which have the same eight-fold index permutational symmetry as te one- and two-electron integral derivatives they multiply. This separation of deformation-dependent and -independent factors allows for single-pass integral-derivative-driven implementation of the gradient and hessian expressions.     

Direct energy functional minimization under orthogonality constraints

The direct energy functional minimization problem in electronic structure theory, where the single-particle orbitals are optimized under the constraint of orthogonality, is explored. We present an orbital transformation based on an efficient expansion of the inverse factorization of the overlap matrix that keeps orbitals orthonormal. The orbital transformation maps the orthogonality constrained energy functional to an approximate unconstrained functional, which is correct to some order in a neighborhood of an orthogonal but approximate solution. A conjugate gradient scheme can then be used to find the ground state orbitals from the minimization of a sequence of transformed unconstrained electronic energy functionals. The technique provides an efficient, robust, and numerically stable approach to direct total energy minimization in first principles electronic structure theory based on tight-binding, Hartree-Fock, or density functional theory. For sparse problems, where both the orbitals and the effective single-particle Hamiltonians have sparse matrix representations, the effort scales linearly with the number of basis functions N in each iteration. For problems where only the overlap and Hamiltonian matrices are sparse the computational cost scales as O(M2N), where M is the number of occupied orbitals. We report a single point density functional energy calculation of a DNA decamer hydrated with 4003 water molecules under periodic boundary conditions. The DNA fragment containing a cis-syn thymine dimer is composed of 634 atoms and the whole system contains a total of 12,661 atoms and 103,333 spherical Gaussian basis functions.