Matrix element factorization in the unitary group approach for configuration interaction calculations

Techniques of diagrammatic spin algebra are employed to derive segment factorization formulas for spin-adapted matrix elements of one- and two-electron excitation operators. The spin-adapted basis is formed by the Yamanouchi–;Kotani geneological coupling method, and therefore constitutes an irreducible basis of the unitary group U(N), as prescribed by Gel'fand and Tsetlin. Several features distinguish this paper from similar work that has recently been published. First, intermediate steps in the derivation of each segment factor are fully documented. Comprehensive tables list the spin diagrams and phases that contribute to the possible segment factors. Second, a special effort has been made to distinguish between those parts of a segment factor that can be ascribed to a spin diagram and those parts which arise from the orbitals. The results of this paper should thus be useful for those who wish to extend diagrammatic spin algebra to evaluation of matrix elements for states built from nonorthogonal orbitals. Third, a novel graphical method has been introduced to keep track of phase changes that are induced by line up permutations of creation and annihilation operators. This technique may be useful for extension of our analysis to higher excitations. The necessary concepts of second quantization and diagrammatic spin algebra are developed in situ, so the present derivation should be accessible to those who have little prior knowledge of such methods.     

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.     

The modified Hartree–Fock self-consistent field equation

A one-electron correlation operator is introduced into the Hartree–Fock self-consistent field equation. The correlation operator is derived from the second-order perturbation theory. Energies of atomic and molecular systems calculated from this modified Hartree–Fock equation are equal to that from second-order perturbation of Hartree–Fock equation. The modified equation can also be solved self-consistently by the LCAO approximation. We also presented the modified expressions for other operators.     

Molecular structural formulas as one-electron density and hamiltonian operators: the VIF method extended

The valency interaction formula (VIF) method is given a broader and more general interpretation in which these simple molecular structural formulas implicitly include all overlaps between valence atomic orbitals even for interactions not drawn in the VIF picture. This applies for VIF pictures as one-electron Hamiltonian operators as well as VIF pictures as one-electron density operators that constitute a new implementation of the VIF method simpler in its application and more accurate in its results than previous approaches. A procedure for estimating elements of the effective charge density-bond order matrix, Pmunu, from electron configurations in atoms is presented, and it is shown how these lead to loop and line constants in the VIF picture. From these structural formulas, one finds the number of singly, doubly, and unoccupied molecular orbitals, as well as the number of molecular orbitals with energy lower, equal, and higher than -1/2Eh, the negative of the hydrogen atom's ionization energy. The VIF results for water are in qualitative agreement with MP2/6311++G3df3pd, MO energy levels where the simple VIF for water presented in the earlier literature does not agree with computed energy levels. The method presented here gives the simplest accurate VIF pictures for hydrocarbons. It is shown how VIF can be used to predict thermal barriers to chemical reactions. Insertion of singlet carbene into H2 is given as an example. VIF pictures as one-electron density operators describe the ground-state multiplicities of B2, N2, and O2 molecules and as one-electron Hamiltonian operators give the correct electronegativity trend across period two. Previous implementations of VIF do not indicate singly occupied molecular orbitals directly from the pictorial VIF rules for these examples. The direct comparison between structural formulas that represent electron density and those that represent energy is supported by comparison of a simple electronegativity scale, chiD=N/n2, with well-known electronegativity scales of Pauling, Mulliken, and Allen. This scale comes from the method used to calculate Pmumu for sp3 hybridized period-two elements and is comparable to electronegativity because it has the same form as <1/r> for hydrogenic orbitals. It therefore provides a physical basis for the representation of one electron density and Hamiltonian operators by the same VIF picture.     

Consistent modifications of SINDO1: I. Approximations and parameters

A consistent modification, MSINDO, of the semiempirical MO method SINDO1 is presented. Different basis sets are used for one- and two-center interactions. The treatment of the core matrix elements in the nonorthogonal basis is retained with changes only for hydrogen and 3d orbitals. Orthogonalization corrections are now restricted to nonvanishing core matrix elements in the INDO approximation. The set of atomic parameters is increased, but bond parameters are no longer used. An automatic nonlinear least-squares algorithm with a restricted step constraint is used for the optimization of parameters. Heats of formation are adjusted with inclusion of zero-point energies obtained by a scaling procedure of the force constant matrix. The present version MSINDO provides significant improvements over previous versions. A brief comparison for ground-state properties of the elements H, C, N, O, F, and Na to Cl is given. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 563–571, 1999     

Comparative study of the Hartree–Fock and the Hartree–Fock–Slater method

The Hartree–Fock method (standard Roothaan closed-shell HF –LCAO theory) and the Hartree–Fock–Slater method (restricted HFS –LCAO –DV method developed by Baerends and Ros) have been compared with emphasis on the respective one-electron equations and on the matrix elements of the respective Fock operators. Using the same STO basis in the two cases, the matrix elements of the Fock operators and of their separate one-electron, Coulomb, and exchange contributions have been calculated for the same orbitals and density of the ground state of the diatomic molecule ZnO. The effects of methodical (exchange potential) and numerical (DV method, density fit) differences between the HF and HFS methods on the various matrix elements have been analyzed. As expected the methodical effect prevails and is responsible for the higher (less negative) values of the matrix elements of the HFS Fock operator compared to those of the HF Fock operator. Numerical effects are observable also and are caused by the difference in integration procedures (DV method), not by the density fit.     

Maximal orbital analysis of molecular wavefunctions

We describe a new way to decompose one-electron orbitals of a molecule into atom-centered or fragment-centered orbitals by an approach that we call “maximal orbital analysis” (MOA). The MOA analysis is based on the corresponding orbital transformation (COT) that has the unique mathematical property of maximizing any sub-trace of the overlap matrix, in Hilbert metric sense, between two sets of nonorthogonal orbitals. Here, one set comprises the molecule orbitals (Hartree–Fock, Kohn–Sham, complete-active-space, or any set of orthonormal molecular orbitals), the other set comprises the basis functions associated with an atom or a group of atoms. We show in prototypical molecular systems such as a water dimer, metal carbonyl complexes, and a mixed-valent transition metal complex, that the MOA orbitals capture very well key aspects of wavefunctions and the ensuing chemical concepts that govern electronic interactions in molecules. © 2019 Wiley Periodicals, Inc.     

Theory of two shells of atomic electrons using non-orthogonal radial orbitals

A theory for handling non-orthogonal radial orbitals of two shells of atomic electrons based on the mathematical apparatus of irreducible tensor operators is presented. The general expressions for one- and two-electron operator matrix elements are given.     

Valence bond approach exploiting Clifford algebra realization of Rumer–Weyl basis

A detailed algorithm is described that enables an implementation of a general valence bond (VB ) method using the Clifford algebra unitary group approach (CAUGA ). In particular, a convenient scheme for the generation and labeling of classical Rumer–Weyl basis (up to a phase) is formulated, and simple rules are given for the evaluation of matrix elements of unitary group generators, and thus of any spin-independent operator, in this basis. The case of both orthogonal and nonrothogonal atomic orbital bases is considered, so that the proposed algorithm can also be exploited in molecular orbital configuration interaction calculations, if desired, enabling a greater flexibility for N-electron basis-set truncation than is possible with the standard Gel'fand–Tsetlin basis. Finally, an exploitation of this formalism for the VB method, based on semiempirical Pariser–Parr–Pople (PPP )-type Hamiltonian and nonorthogonal overlap-enhanced atomic orbital basis, and its computer implementation, enabling us to carry out arbitrarily truncated or full VB calculations, is described in detail.     

Scale‐adaptive tensor algebra for local many‐body methods of electronic structure theory

While the formalism of multiresolution analysis, based on wavelets and adaptive integral representations of operators, is actively progressing in electronic structure theory (mostly on the independent‐particle level and, recently, second‐order perturbation theory), the concepts of multiresolution and adaptivity can also be utilized within the traditional formulation of correlated (many‐particle) theory based on second quantization and the corresponding (generally nonorthogonal) tensor algebra. In this article, we present a formalism called scale‐adaptive tensor algebra, which introduces an adaptive representation of tensors of many‐body operators via the local adjustment of the basis set quality. Given a series of locally supported fragment bases of a progressively lower quality, we formulate the explicit rules for tensor algebra operations dealing with adaptively resolved tensor operands. The formalism suggested is expected to enhance the applicability of certain local correlated many‐body methods of electronic structure theory, for example, those directly based on atomic orbitals (or any other localized basis functions in general). © 2014 Wiley Periodicals, Inc.     

Valence-bond calculations on ZNO and HGO using integrals computed through the semiempirical AM1 method

Valence-bond calculations have been carried out on ZnO and HgO using a basis set of Slatertype atomic orbitals and the one- and two-electron integrals as computed in the semiempirical AM 1 molecular orbital method. The zero differential overlap approximation has been used to calculate integrals between atomic orbital Slater determinants using the rules for matrix elements between determinants formed by orthogonal orbitals. Diabatic and adiabatic curves have been analyzed for the two systems, and results compared with molecular orbital AM 1 results. © 1992 John Wiley & Sons, Inc.     

Calculation of transition matrix elements by nonsingular orbital transformations

A general strategy is described for the evaluation of transition matrix elements between pairs of full class CI wave functions built up from mutually nonorthogonal molecular orbitals. A new method is proposed for the counter‐transformation of the linear expansion coefficients of a full CI wave function under a nonsingular transformation of the molecular‐orbital basis. The method, which consists in a straightforward application of the Cauchy–Binet formula to the definition of a Slater determinant, is shown to be simple and suitable for efficient implementation on current high‐performance computers. The new method appears mainly beneficial to the calculation of miscellaneous transition matrix elements among individually optimized CASSCF states and to the re‐evaluation of the CASCI expansion coefficients in Slater‐determinant bases formed from arbitrarily rotated (e.g., localized or, conversely, delocalized) active molecular orbitals. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Unitary group tensor operator algebras for many-electron systems: II. One- and two-body matrix elements

Relying on our earlier results in the unitary group Racah-Wigner algebra, specifically designed to facilitate quantum chemical calculations of molecular electronic structure, the tensor operator formalism required for an efficient evaluation of one- and two-body matrix elements of molecular electronic Hamiltonians within the spin-adapted Gel'fand-Tsetlin basis is developed. Introducing the second quantization-like creation and annihilation vector operators at the unitary group [U(n)] level, appropriate two-box symmetric and antisymmetric irreducible tensor operators as well as adjoint tensors are defined and their matrix elements evaluated in the electronic Gel'fand-Tsetlin basis as single products of segment values. Using these tensor operators, the matrix elements of one- and two-body components of a general electronic Hamiltonian are found. Explicit expressions for all relevant quantities pertaining to at most two-column irreducible representations that are required in molecular electronic structure calculations are given. Relationships with other approaches and possible future extensions of the formalism to partitioned bases or spin-dependent Hamiltonians are discussed.On leave from: Department of Chemistry, Xiamen University, Xiamen, Fujian, PR China.     

Orthogonal natural atomic orbitals form an appropriate one-electron basis for expanding CASSCF wave functions into localized bonding schemes and their weights

Localized bonding schemes and their weights have been obtained for the pi-electron system of nitrone by expanding complete active space self-consistent field wave functions into a set of Slater determinants composed of orthogonal natural atomic orbitals (NAOs) of Weinhold and Landis (Valency and Bonding: A Natural Bond Orbital Donor-Acceptor Perspective, 2005). Thus, the derived bonding schemes are close to orthogonal valence bond structures. The calculated sequence of bonding scheme weights accords with the sequence of genuine resonance structure weights derived previously by Ohanessian and Hiberty (Chem Phys Lett 1987, 137, 437), who employed nonorthogonal atomic orbitals. This accord supports the notion that NAOs form an appropriate orthogonal one-electron basis for expanding complete active space self-consistent field wave functions into meaningful bonding schemes and their weights.     

An efficient self-consistent field method for large systems of weakly interacting components

An efficient method for removing the self-consistent field (SCF) diagonalization bottleneck is proposed for systems of weakly interacting components. The method is based on the equations of the locally projected SCF for molecular interactions (SCF MI) which utilize absolutely localized nonorthogonal molecular orbitals expanded in local subsets of the atomic basis set. A generalization of direct inversion in the iterative subspace for nonorthogonal molecular orbitals is formulated to increase the rate of convergence of the SCF MI equations. Single Roothaan step perturbative corrections are developed to improve the accuracy of the SCF MI energies. The resulting energies closely reproduce the conventional SCF energy. Extensive test calculations are performed on water clusters up to several hundred molecules. Compared to conventional SCF, speedups of the order of (N/O)2 have been achieved for the diagonalization step, where N is the size of the atomic orbital basis, and O is the number of occupied molecular orbitals.     

Big data analysis of ab Initio molecular integrals in the neglect of diatomic differential overlap approximation

Most modern semiempirical quantum-chemical (SQC) methods are based on the neglect of diatomic differential overlap (NDDO) approximation to ab initio molecular integrals. Here, we check the validity of this approximation by computing all relevant integrals for 32 typical organic molecules using Gaussian-type orbitals and various basis sets (from valence-only minimal to all-electron triple-ζ basis sets) covering in total more than 15.6 million one-electron (1-e) and 10.3 billion two-electron (2-e) integrals. The integrals are calculated in the nonorthogonal atomic basis and then transformed by symmetric orthogonalization to the Löwdin basis. In the case of the 1-e integrals, we find strong orthogonalization effects that need to be included in SQC models, for example, by strategies such as those adopted in the available OMx methods. For the valence-only minimal basis, we confirm that the 2-e Coulomb integrals in the Löwdin basis are quantitatively close to their counterparts in the atomic basis and that the 2-e exchange integrals can be safely neglected in line with the NDDO approximation. For larger all-electron basis sets, there are strong multishell orthogonalization effects that lead to more irregular patterns in the transformed 2-e integrals and thus cast doubt on the validity of the NDDO approximation for extended basis sets. Focusing on the valence-only minimal basis, we find that some of the NDDO-neglected integrals are reduced but remain sizable after the transformation to the Löwdin basis; this is true for the two-center 2-e hybrid integrals, the three-center 1-e nuclear attraction integrals, and the corresponding three-center 2-e hybrid integrals. We consider a scheme with a valence-only minimal basis that includes such terms as a possible strategy to go beyond the NDDO integral approximation in attempts to improve SQC methods. © 2018 Wiley Periodicals, Inc.     

Calculation of molecular systems via coordinate wave functions

The system of charges is in a state with a given total spin S, which is described by a configuration of one-electron orbitals with arbitrary filling (subject to the Pauli principle). Expressions are derived for the matrix elements of operators F and G that are independent of the spin. The energy of the interaction between the completely filled orbitals and the singly filled ones is found to be independent of the spin of the latter. The formulas may be used with the tables of [2] to derive directly the expressions for the matrix elements of a configuration having an arbitrary number of completely filled orbitals and up to six singly filled ones.