Configuration interaction matrix elements. II. Graphical approach to the relationship between unitary group generators and permutations

The explicit expressions for the matrix elements of unitary group generators between geminally antisymmetric spin-adapted N-electron configurations in terms of the orbital occupancies and spin factors, given as spin function matrix elements of appropriate orbital permutations, are derived using the many-body time-independent diagrammatic techniques. It is also shown how this approach can be conveniently combined with graphical methods of spin algebras to obtain explicit expressions for the spin factors, once a definite coupling scheme is chosen. This method yields explicit expressions for the orbital permutations defining the spin factors. However, if desired, the explicit determination of line-up permutations can be avoided in this approach, since they are implicitly contained in the orbital diagrams. It also clearly indicates why the geminally antisymmetric spin functions have to be used when a simple formalism is desired.     

Role of coupling operators on open shell RHF equations

We have derived the expressions for the extremum condition of E, corresponding to any wave function. These expressions are given as a function of the spin orbitals. We have carried out the derivation considering the spin orbitals as vectors belonging to an orthonormal basis. The corresponding variational equations have been derived introducing the condition that the norm of the wave function is constant, as the only additional constraint.From the expression obtained for the first variation of the matrix elements of H, as a function of the spin orbitals, we have derived the RHF equations for a simple case.In the present procedure, the couplings between orbitals of different shells appear directly, being defined explicitly, and they may be taken as corresponding with the elements of a Hermitian matrix.The calculations that we have carried out show that the coupling operators defined in the paper give results which are variationally correct.     

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.     

Coupled cluster approaches with an approximate account of triexcitations and the optimized inner projection technique

Summary A general orthogonally spin-adapted formalism for coupled cluster (CC) approaches, with an approximate account of triexcited configurations, and for optimized inner projection (OIP) technique is described. Modifying the linear part of the CC equations for pair clusters (CCD) we obtain the orthogonally spin-adapted, non-iterative version of the CCDT-1 method of Bartlett et al. [J. Chem. Phys. 80, 4371 (1984), 81, 5906 (1984), 82, 5761 (1985)]. Similar modification of an approximate coupled pair theory corrected for connected quadruply excited clusters (ACPQ) yields a new approach called ACPTQ. Both the CCDT-1 and ACPTQ methods can be formulated in terms of effective interaction matrix elements between the orthogonally spin-adapted biexcited singlet configurations. The same matrix elements also appear in the orthogonally spin-adapted form of the CCD + T(CCD) perturbative estimate of triply excited contributions due to Raghavachari [J. Chem. Phys. 82, 4607 (1985)] and Urban et al. [J. Chem. Phys. 83, 4041 (1985)], and in the OIP method when applied to the Pariser-Parr-Pople (PPP) model Hamiltonians. We use the diagrammatic approach based on the graphical methods of spin algebras to derive the explicit form of these interaction matrix elements. Finally, the relationship between different diagrammatic spin-adaptation procedures and their relative advantages are discussed in detail.Also affiliated with the Department of Chemistry, and Guelph-Waterloo Center for Graduate Work in Chemistry, Waterloo Campus, University of Waterloo, Waterloo, Ontario, Canada. Killam Research Fellow 1987–89     

Spin-free quantum theory. XXV. The unitary-group formulation of fermion many-body theory

The fermion unitary group formulation (UGF ) of many-body theory is based on the unitary group U(2n) where n is the number of freeon orbitals. This formulation, which conserves particle-number but not spin, is isomorphic to the particle-number-conserving, second-quantized formulation (SQF ). In UGF we derive the familiar diagrammatic algorithm for matrix elements, M(Y) = (?1)^H+L where H and L denote the numbers of hole lines and loops in the diagram D(Y) of M(Y). The unitary group derivation is considerably simpler than is the conventional, second-quantized derivation that employs time-dependence, Wick's theorem, normal-order, and contractions. In neither fermion UGF nor SQF is spin conserved. We carry out in UGF the spin-projection (symmetry adaptation to SU (2)) of the fermion vectors and obtain with a spin-free Hamiltonian the same matrix elements as with the freeon UGF (part 24 of this series). The fermion unitary group formulation for a spin-free Hamiltonian should be regarded as an alternate path to spin-free quantum chemistry.     

Configuration interaction matrix elements. I. Algebraic approach to the relationship between unitary group generators and permutations

Matrix elements of unitary group generators between spin-adapted antisymmetric states are shown to be proportional to spin matrix elements of so-called “line-up” permutations. The proportionality factor is given explicitly as a simple function of the orbital occupation numbers. If one bases the theory on ordered orbital products, the line-up permutations are given a priori. The final formulas have a very simple structure; this is a direct consequence of the fact that the spin functions have been taken to be geminally antisymmetric.     

Unitary group approach to the many-electron problem. II. Adjoint tensor operators for U(n)

In this paper we present a derivation of the U(n) adjoint coupling coefficients for the representations appropriate to many-electron systems. Since the states of a many-fermion system are to comprise the totally antisymmetric Nth rank tensor representation of U(2n), the work of this paper enables the matrix elements of the U(2n) generators to be evaluated directly in the U(n) × U(2) (i.e., spin orbit) basis using their transformation properties as adjoint tensor operators. A connection between the adjoint coupling coefficients, as derived in this paper, and the matrix elements of certain (spin independent) two-body operators is also presented. This indicates that in CI calculations, one may obtain the matrix elements of spin-dependent operators from the known matrix elements of certain spin-independent two-body operators. In particular this implies a segment-level formula for the matrix elements of the U(2n) generators in the spin-orbit basis.     

Spin-free approach for evaluation of electronic matrix elements using character operators of ℒN

A spin-free method is presented for evaluating electronic matrix elements over a spin-independent many-electron Hamiltonian. The spin-adapted basis of configuration state functions is obtained using a nonorthogonal spin basis consisting of projected spin eigenfunctions. The general expressions for the matrix elements are given explicitly, and it is demonstrated how the matrix elements may be obtained simply from the knowledge of the irreducible characters of the permutation group ℒ_N. The presented formulas are very general and may be applied in connection with both spin-coupled valence bond studies and in conventional configuration interaction (CI) methods based on an orthonormal orbital basis. © 1996 John Wiley & Sons, Inc.     

Symmetric group approach to relativistic CI. I. General formalism

General formulas for matrix elements of spin-dependent operators in a basis of spin-adapted antisymmetrized products of orthonormal orbitals are derived. The resulting formalism may be applied to construction of the Hamiltonian matrices both for Pauli and for projected no-pair relativistic configuration interaction methods. From a formal point of view, it is a generalization of the symmetric group approach to the CI method for the case of spin-dependent Hamiltonians. © 1997 John Wiley & Sons, 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.     

New derivation of the Waller–Hartree–Fock spatial wave function

A new derivation is given for the Waller–Hartree–Fock double-determinantal spatial wave function. One starts from the single-determinant wave function in which a orbitals are doubly occupied, and decomposes it into a sum of products of spatial and spin functions. The spatial product of the first genealogical spin eigenfunction is a double-determinantal function. The derivation is based on the simple form of U_1?(P) when the representation matrix is obtained from the genealogical spin eigenfunction.     

Efficient generation of Heisenberg Hamiltonian matrices for VB calculations of potential energy surfaces

The spin‐Hamiltonian valence bond theory relies upon covalent configurations formed by singly occupied orbitals differing by their spin counterparts. This theory has been proven to be successful in studying potential energy surfaces of the ground and lowest excited states in organic molecules when used as a part of the hybrid molecular mechanics—valence bond method. The method allows one to consider systems with large active spaces formed by n electrons in n orbitals and relies upon a specially proposed graphical unitary group approach. At the same time, the restriction of the equality of the numbers of electrons and orbitals in the active space is too severe: it excludes from the consideration a lot of interesting applications. We can mention here carbocations and systems with heteroatoms. Moreover, the structure of the method makes it difficult to study charge‐transfer excited states because they are formed by ionic configurations. In the present work we tackle these problems by significant extension of the spin‐Hamiltonian approach. We consider (i) more general active space formed by n ± m electrons in n orbitals and (ii) states with the charge transfer. The main problem addressed is the generation of Hamiltonian matrices for these general cases. We propose a scheme combining operators of electron exchange and hopping, generating all nonzero matrix elements step‐by‐step. This scheme provides a very efficient way to generate the Hamiltonians, thus extending the applicability of spin‐Hamiltonian valence bond theory. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Orthogonally spin-adapted coupled-cluster equations involving singly and doubly excited clusters. Comparison of different procedures for spin-adaptation

A particularly compact form of the orthogonally spin-adapted coupled-cluster equations involving all singly and doubly excited clusters is derived for the general case of a non-Hartree–Fock closed-shell reference determinant. The diagrammatic approach based on the graphical methods of spin algebras is applied. The relationship of different diagrammatic procedures for spin-adaptation, employing both bare and spin-adapted two-electron interaction vertices, is discussed. A comparison with the results obtained with algebraic spin-adaption approaches is also given.     

Second quantization for systems with a constant number of particles in the Dirac notation

The relation between the completeness condition for an appropriate one-particle basis set and the occupation number representation (second quantization) is shown for the time-independent case. The explicit expressions for the basic symmetric operators are derived in the Dirac bra–ket notation. The physical meaning of these operators, the algebra as well as the connections with the one-electron density matrix and with the projector on the Fermi sea in the one-electron approximation, follow directly from these expressions. The generalization for a nonorthogonal basis and the algebra for corresponding basic operators are formulated. The connection with the notion of the molecular diagrams of different kinds for the nonorthogonal atomic orbitals is shown. The Mulliken populations and the Chirgwin–Coulson bond orders are equal to the diagonal and offdiagonal elements of the molecular diagram 1, respectively. The matrix elements of the projector on the Fermi sea in the one-electron approximation in the representation of nonorthogonal atomic orbitals are elements of the molecular diagram 2.     

Application of graphical methods of spin algebras to limited CI approaches. I. Closed shell case

The time independent diagrammatic technique based on the mathematical methods of quantum electrodynamics (second quantization, Wick's theorem, Feynman-like diagrams) is combined with graphical techniques of spin algebras to derive general expressions for the matrix elements of spin independent one- and two-particle operators between spin symmetry adapted ground, mono- and bi-excited configurations of a closed shell system. Two coupling schemes are considered for bi-excited states and their relative merits are discussed. Finally, the results are used to derive compact expressions for the coupling coefficients of the direct configuration interaction from molecular integrals (CIMI ) method.     

Explicit formulas for the Hamiltonian matrix elements

The explicit formulas for the evaluation of the Hamiltonian matrix elements are presented. The calculation of the integral coefficients is independent of both the nature of the orbitals and th spin coupling schemes. It is fully automatic and only dependent on the number of doubly and singly occupied orbitals. Further-more, the symmetric group representation matrices are not needed, and the N! problem can be avoided.     

On the spin and symmetry adaptation of the density matrix renormalization group method

We present a spin-adapted density matrix renormalization group (DMRG) algorithm designed to target spin and spatial symmetry states that can be difficult to obtain while using a non-spin-adapted algorithm. The algorithmic modifications that have to be introduced into the usual density matrix renormalization group scheme in order to spin adapt it are discussed, and it is demonstrated that the introduced modifications do not change the overall scaling of the method. The new approach is tested on HNCO, a model system, that has a singlet-triplet curve crossing between states of the same symmetry. The advantages of the spin-adapted DMRG scheme are discussed, and it is concluded that the spin-adapted DMRG method converges better in almost all cases and gives more parallel curves to the full configuration interaction result than the non-spin-adapted method. It is shown that the spin-adapted DMRG energies can be lower than the ones obtained from the non-spin-adapted scheme. Such a counterintuitive result is explained by noting that the spin-adapted method is not a special case of the non-spin-adapted one; consequently, the spin-adapted result is not an upper bound for the non-spin-adapted energy.     

Harmonic oscillator tensors. I. The nondegenerate case

It is shown that the Heisenberg Lie algebra of the nondegenerate harmonic oscillator leads to a basis {J₊, J₀, J_?} of LASU (2). The Hamiltonian of the system is proportional to J₀, and the basis elements give rise to irreducible tensors in the associative enveloping algebra of the Heisenberg Lie algebra. The construction of these irreducible tensors is studied with special attention being paid to the case in which they act upon a single vector space spanned by the harmonic oscillator basis functions. A tensor coupling rule is developed, and useful application is made of it in the calculation of general expressions for vibrational operators and their matrix elements. Throughout, the value of the additional algebraic quantum numbers (l, m) is emphasized.