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

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

定域片断轨道基组的对称化

片断的UHF运算不能保证每个片断轨道具有确切的电子占据数,故Kost定域化是必需的.当片断产生于多键断裂时,在确保目标轨道单占据性的同时,Kost定域破坏了轨道基组的对称性.为此,在Kost定域化后,必须对单占据轨道作2×2对称化旋转后,再作有条件的RHF运算.以乙烯基片断CHCH(波二烯分子中的一个片断)为例,详述了对称化的方法、原理和计算程序.以C-H片断为例,细述三单键片断轨道基组对称化的特殊性.介绍C-HR参考键长选择的判据,探讨键长与选择Gaussian基组大小的关系.     

定域分子轨道中的σ—π分离

对含重键的分子体系,分子轨道定域化会涉及到两种完全不同的重键描述,即等价的重键或“香蕉”键和不等价的σ和π键,文献曾利用杂化轨道法对此进行过讨论。对于定域分子轨道,具体得出哪种描述取决于所采用的定域准则及计算中采用的近似方法。对于从头算法的定域化研究,Boys定域准则强烈地趋向于等价重键描述;Ruedenberg定域准则只是对未共轭的重键体系有较强的等价重键描述倾向,对共轭的重键体系,这种倾向性明显减     

单占据片断轨道的组合式(Kost+Perkin)定域化

为了使开壳层片断轨道具有正确的集居数,正则轨道必须定域化.但是,在确保目标轨道Φ_sdan单占据性的同时,Kost定域化也破坏了Φ_s特定的对称性.实际计算表明,Perkin定域化可以弥补Kost程序的缺陷,将单占据的Φ_s转化成高度定域的、对称的片断轨道.在片断分子中,C-H_R键长r的选择和Kost定域化的方式对Perkin定域化的成败具有重大的影响.当Gaussian基组为STO-3G3-21G和4-31G时,r应为0.1nm;但在6-31G水平下,必须r=0.09nm.     

一种计算量小的构造合理最紧缩定域轨道的方法

提出一种计算量小的构造合理最紧缩定域化轨道的方法.非正交定域轨道(NOLMO)没有“正交化尾巴”,比正交定域轨道(OLMO)更局域化、更紧缩、更具可移植性,从而更适合用于化学问题的理论研究.但若在用变分法确定最大限度局域化NOLMO时单纯取消正交条件而不附加其它的限制,会得到趋于线性相关的不合理结果.提出用强制NOLMO与投影自然键轨道(PNBO)的重心重合代替正交化条件确定合理的最大限度局域化NOLMO(合理最紧缩定域轨道)的方法.对一系列不同类型分子的计算结果表明,用该方法可以得到空间分布合理及线性独立的NOLMO,延伸度与文献已有的最佳结果接近,而计算量大幅度降低.由于构造投影自然键轨道的计算量随分子中原子数目的增加只是线性增长,因此,该方法可以用来构造较大体系的合理最紧缩定域轨道     

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

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

建立在群对称轨道(SMO)及特征组态(CD)基础上的CI方法与程序发展

ISMO－CDCI方法点群理论，特别是不可约张量方法，在量子化学理论方法发展以及简化概念与计算方面，发挥了重要作用．但在国外的量子化学计算程序（如G94）中，在后自洽场（opt－SCF）计算方面，很少用到对称约化．在文献中，只见有关hbeltah群（DZh与它的子群）对称约化用于组态相关计算的报导．由于多重耦合系数的计算复杂，蜕化不可约表示多体问题的对称约化难于得到解决．我们提出了一个统一与普遍的方法，它能解决所有分子体系多体相关的点群对称约化问题［‘－6］，这个方法的核心是群对称轨道（SMO）概念的提出．SMO的基本特…     

EHMO方法的分子轨道定域化

用Foster-Boys的定域化准则讨论了EHMO方法的分子轨道定域化问题,提出用双中心重叠积分近似计算双中心轨道偶极矩积分方法,得到的EHMO定域分子轨道与严格定域化结果接近,与从头计算方法的定域化结果定性一致。     

分子的模糊对称性

将模糊数学方法引入对分子对称点群的研究,建立描述具有不完整分子对称性的模糊点对称群(集合).建立具有模糊对称性分子轨道的模糊表示及其模糊特征标(模糊广义宇称).通过对典型的线状分子、平面分子以及非平面的立体分子等进行分析,展示了一个新的理论化学园地.初步探讨了具有模糊对称性的动态反应体系.从模糊对称性出发,探讨了分子轨道对称守恒原理的半定量特征.     

偶极积分与Boys法分子轨道定域化

定域分子轨道在分子体系的化学图象和物理图象之间充当重要的桥梁作用,它的产生依赖于定域化准则,其中最普遍使用的是Foster-Boys和Edmiston-Ruedenberg(E—R)提出的两种定域化准则。这两种定域化准则是等价的,因而结果也是一致的。但对于E—R定域化来说,由于涉及到大量的多中心积分的计算,计算极为费时,因而远不如Foster-Boys定     

Optimal group symmetric localized molecular orbitals

Summary The concept and generating method of optimum group symmetric localized molecular orbitals (OSLMOs) are proposed. The OSLMOs have strong points of orthogonality, equivalence and symmetry, and they are simultaneously as close to the classical VB structure as possible. By using the OSLMOs as one-electron orbitals the multiconfigurational correlation calculations are reduced. The scheme is also a valuable popularization and development to hybridization theory.     

Study of localized molecular orbitals using group theory methods and its approach to the many-electron correlation problem. IV. The symmetry-adaptation of many-center integrals and hamiltonian matrix elements in MCSCF calculations

Group theoretic methods are presented for the transformations of integrals and the evaluation of matrix elements encountered in multiconfigurational self-consistent field (MCSCF) and configuration interaction (CI) calculations. The method has the advantages of needing only to deal with a symmetry unique set of atomic orbitals (AO) integrals and transformation from unique atomic integrals to unique molecular integrals rather than with all of them. Hamiltonian matrix element is expressed by a linear combination of product terms of many-center unique integrals and geometric factors. The group symmetry localized orbitals as atomic and molecular orbitals are a key feature of this algorithm. The method provides an alternative to traditional method that requires a table of coupling coefficients for products of the irreducible representations of the molecular point group. Geometric factors effectively eliminate these coupling coefficients. The saving of time and space in integral computations and transformations is analyzed. © 1994 by John Wiley & Sons, Inc.     

Performance of dynamically weighted multiconfiguration self-consistent field and spin-orbit coupling calculations of diatomic molecules of Group 14 elements

The efficacy of several multiconfiguration self-consistent field (MCSCF) methods in the subsequent spin-orbit coupling calculations was studied. Three MCSCF schemes to generate molecular orbitals were analyzed: state-specific, state-averaged, and dynamically weighted MCSCF. With Sn(2)(+) as the representative case, we show that the state-specific MCSCF orbitals lead to discontinuities in potential energy curves when avoided crossings of electronic states occur; this problem can be solved using the state-averaged or dynamically weighted MCSCF orbitals. The latter two schemes are found to give similar results when dynamic electron correlation is considered, which we calculated at the level of multiconfigurational quasidegenerate perturbation theory (MCQDPT). We employed the recently developed Douglas-Kroll spin-orbit adapted model core potential, ZFK3-DK3, and the dynamically weighted MCSCF scheme to calculate the spectroscopic constants of the mono-hydrides and compared them to the results obtained using the older set of potentials, MCP-TZP. We also showed that the MCQDPT tends to underestimate the dissociation energies of the hydrides and discussed to what extent coupled-cluster theory can be used to improve results.     

价键理论新进展

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

Pathway analysis of super-exchange electronic couplings in electron transfer reactions using a multi-configuration self-consistent field method

We present a novel pathway analysis of super-exchange electronic couplings in electron transfer reactions using localized molecular orbitals from multi-configuration self-consistent field (MCSCF) calculations. In our analysis, the electronic coupling and the tunneling pathways can be calculated in terms of the configuration interaction (CI) Hamiltonian matrix obtained from the localized MCSCF wave function. Making use of the occupation restricted multiple active spaces (ORMAS) method can effectively produce the donor, acceptor, and intermediate configuration state functions (CSFs) and CIs among these CSFs. In order to express the electronic coupling as a sum of individual tunneling pathways contributions, we employed two perturbative methods: L?wdin projection-iteration method and higher-order super-exchange method. We applied them to anion couplings of butane-1,4-diyl and pentane-1,5-diyl. The results were (1) the electronic couplings calculated from the two perturbative methods were in reasonable agreement with those from a non-perturbative method (one-half value of the energy difference between the ground and first excited states), (2) the main tunneling pathways consisted of a small number of lower-order super-exchange pathways where bonding, anti-bonding, or extra-valence-shell orbitals were used once or twice, and (3) the interference among a huge number of higher-order super-exchange pathways significantly contributed to the overall electronic coupling, whereas each of them contributed only fractionally. Our method can adequately take into account both effects of non-dynamical electron correlation and orbital relaxation. Comparing with the analyses based on the Koopmans' theorem (ignoring both effects) and the ORMAS-CIs from frozen localized reference orbitals (ignoring the effect of orbital relaxation), we discuss these effects.     

Equivalent orbitals for multiconfigurational spin‐tensor electron propagator method (MCSTEP): The vertical ionization potentials of B,NO, CF,and OF

The multiconfigurational spin tensor electron propagator method (MCSTEP) was developed as an implementation of electron propagator/single particle Green's function methods for ionization potentials (IPs) and electron affinities (EAs). MCSTEP was specifically designed for open shell and highly correlated (nondynamically correlated) initial states. For computational efficiency the initial state used in MCSTEP is typically a small complete active space (CAS) multiconfigurational self‐consistent field (MCSCF) state. If in a molecule there are some degenerate orbitals which are not fully or half occupied, usual MCSCF calculations will make these orbitals inequivalent, i.e., the occupied ones will be different from the nonoccupied ones, so that the degeneracy is broken. In this article, we use a state averaged MCSCF method to get equivalent orbitals for the initial state and import the integrals into the subsequent MCSTEP calculations. This gives, in general, more reliable MCSTEP vertical IPs. © 2008 Wiley Periodicals, Inc., 2008     

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

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

Generalized-molecular-orbital theory: Simple multiconfiguration self-consistent-field method

Even after completing a multiconfiguration self-consistent-field (MCSCF ) calculation, one must often include additional configuration interaction (CI ) to obtain quantitative or semiquantitative results. There is some question of whether the prior MCSCF calculation is worthwhile, if additional CI is needed later. We have developed a new MCSCF computational method, which, because of our assumptions about the nature of the configurations, yields one Fock-like operator for all the “filled” orbitals (high occupation numbers) and a second Fock-like operator for all the “virtual” orbitals (low occupation numbers). Since there are only two matrices to build, our method is considerably faster than other MCSCF approaches. Because of these similarities to standard molecular-orbital (MO ) calculations, we have termed our approach generalized-molecular-orbital (GMO ) theory. However, the “virtual” orbitals, unlike those of standard MO theory, are optimized to correlate the “filled” ones and can he used in a subsequent CI calculation. Results are presented for the correlation energy of H₂O, the spectroscopic constants of N₂, the singlet–triplet energy separations in CH₂, and the nature of the chromium–chromium quadruple bond. Although these results are at a very low level of CI , the GMO approach appears to correct for the gross deficiencies of the single-determinant SCF procedure.     

Symmetry adaptation of configuration basis in MCSCF method

Summary A novel approach of space symmetry adaptation is developed for multiconfigurational (MC) functions in fully optimized reaction space and complete active space SCF calculations. The bonded tableau and two box symmetric tableau are basic representations (rep) of configuration functions; the group symmetric localized orbitals are used as one-electron orbitals. The method is proposed for generating a complete and orthonormal set of MC single excited functions. The redundant variable in MCSCF can be eliminated by symmetry adaptation.