Some developments in the spin-extended Hartree-Fock method

All the second-order density matrix spin components for the spin-extended Hartree-Fock method are obtained. The coefficients in the final formulae are only ω_sM, ω_sM±1, ω_sM±2, where ω_sM are the weights of pure states of spin s in the initial unprojected determinant with spin projection M. The eigenvalue problem for the best electron density natural orbitals in the spin-extended method is formulated. All the second-order transition density matrix spin components between pure spin basis functions built of orthogonal orbitals and distinguished by different core choice are also found. This basis may be used on CI calculations.     

Application of the unitary group approach to evaluate spin density for configuration interaction calculations in a basis of S2 eigenfunctions

We present an implementation of the spin‐dependent unitary group approach to calculate spin densities for configuration interaction calculations in a basis of spin symmetry‐adapted functions. Using S² eigenfunctions helps to reduce the size of configuration space and is beneficial in studies of the systems where selection of states of specific spin symmetry is crucial. To achieve this, we combine the method to calculate U(n) generator matrix elements developed by Downward and Robb (Theor. Chim. Acta 1977, 46, 129) with the approach of Battle and Gould to calculate U(2n) generator matrix elements (Chem. Phys. Lett. 1993, 201, 284). We also compare and contrast the spin density formulated in terms of the spin‐independent unitary generators arising from the group theory formalism and equivalent formulation of the spin density representation in terms of the one‐ and two‐electron charge densities.     

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.     

Unitary group approach to the many-electron problem. I. Matrix element evaluation and shift operators

This is the first paper in a series of three directed toward the evaluation of spin-dependent Hamiltonians directly in the spin-orbit basis. In this paper we present a new and complete derivation of the matrix elements of the U(n) generators in the electronic Gel'fand basis. The approach employed differs from previous treatments in that the matrix elements of nonelementary generators are obtained directly. A general matrix element formula is derived which explicitly demonstrates the segment level formalism obtained previously by Shavitt using different methods. A simple relationship between the matrix elements of raising and lowering generators is determined which indicates that in CI calculations, only the matrix elements of raising generators need be calculated. Some results on the matrix elements of products of two generators are also presented.     

Graphical techniques in the configuration interaction approach based on pure slater determinants

A graphical approach to the configuration interaction in the basis of pure Slater determinants is presented. The formulation based on the spin-separated two-slope graph (SSTSG ), enabling the selection of determinants with the fixed M_s value, has a direct relation to the well-known concept of the group-function product. The commonly used excitation criterion and the spatial (Abelian) symmetry properties are analyzed in terms of the graph's internal structure. The Slater formulae for the Hamiltonian matrix elements between determinants, in the particle-hole formalism and in the spin-separated form, are related to different classes of loops within graphs. Some aspects of implementation within both the matrix-element-driven (ME ) and integral driven (ID , direct) CI algorithms are discussed. The presented formulation, of a general complete active space (CAS ) CI type, is a basis of the Graphical Determinantal Configuration Interaction (GDCI ) computer program.     

Two‐component calculations of spin–orbit effects for a van der Waals molecule Rn2

Electronic structures of the weakly bound Rn₂ were calculated by the two‐component Møller–Plesset second‐order perturbation and coupled‐cluster methods with relativistic effective core potentials including spin–orbit operators. The calculated spin–orbit effects are small, but depend strongly on the size of basis sets and the amount of electron correlations. Magnitudes of spin–orbit effects on D_e (0.7–3.0 meV) and R_e (−0.4∼−2.2 Å) of Rn₂ are comparable to previously reported values based on configuration interaction calculations. A two‐component approach seems to be a promising tool to investigate spin–orbit effects for the weak‐bonded systems containing heavy elements. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 139–143, 1999     

Ab initio surface hopping excited-state molecular dynamics approach on the basis of spin–orbit coupled states: An application to the A-band photodissociation of CH3I

Ab initio molecular dynamics approach has been extended to multi-state dynamics on the basis of the spin–orbit coupled electronic states that are obtained through diagonalization of the spin–orbit coupling matrix with the multi-state second-order multireference perturbation theory energies in diagonal elements and the spin–orbit coupling terms at the state-averaged complete active space self-consistent field level in off-diagonal elements. Nonadiabatic transitions over the spin–orbit coupled states were taken into account explicitly by a surface hopping scheme with utilizing the nonadiabatic coupling terms calculated by numerical differentiation of the spin–orbit coupled wavefunctions and analytical nonadiabatic coupling terms. The present method was applied to the A-band photodissociation of methyl iodide, CH₃I + hv → CH₃ + I (²P_3/2)/I* (²P_1/2), for which a pioneering theoretical work was reported by Amatatsu, Yabushita, and Morokuma. The present results reproduced well the experimental branching ratio and energy distributions in the dissociative products. © 2018 Wiley Periodicals, Inc.     

Configuration interaction matrix elements for atoms using permutation group algebra

A procedure is described for the efficient evaluation of the energy matrix elements necessary for atomic configuration-interaction calculations. With the orbital configurations of an N electron system in spin state S written as the irreducible representations [2^1/2N?S, 1^2S] of the permutation group S( N ), it is possible to evaluate readily the energy matrix elements of a spin-free Hamiltonian expressed in terms of the generators of the unitary group. We show how the use of angular momentum ladder operators permits the effective generation of a basis of eigenstates of ??², ??_z as well as ??² and ??_z, for which the energy matrix elements may be evaluated with ease.     

Optimization of virtual orbitals in the framework of a multiconfiguration spin-coupled wave function

A new method is introduced for the optimization of nonorthogonal virtual orbitals for use in general multiconfiguration spin-coupled wave functions. The use of a number of highly effective approximations greatly reduces the computational effort involved, the most important being the use of a second-order perturbation expression for the energy and an approximate expression for the elements of the Hessian. As a result, the overall scheme scales very favourably with respect to the numbers of active electrons and of basis functions, making it suitable for the accurate study of large systems. Benchmark calculations are presented for the dissociation of LiH(X¹Σ⁺) and Li₂(X¹Σ⁺ _g ) using a highly compact four-configuration wave function. Standard spin-coupled valence bond expansions in the same virtual space are required to be significantly larger before equivalent results are obtained. The results are shown to compare very favourably with full valence complete active space self-consistent field calculations using an identical basis, and binding energies are within 4% of the values obtained from full configuration interaction calculations in the same basis set. Received: 10 June 1997 / Accepted: 7 October 1997     

Natural spin orbital analysis of diatomic molecular wave functions in terms of generalized diatomic orbitals I. Outline of the method. Results for the ground state of H2

The method of linear combinations of generalized diatomic orbitals (LCGDO) is combined with the method of configuration interaction (CI). CI wave functions obtained in this way are finally submitted to a natural spin orbital analysis; the resulting natural spin orbitals are expansions in terms of generalized diatomic orbitals.For the ground state of H₂, a one-determinantal-approach with a single completely optimized one-electron basis function nearly reproduces the Hartree-Fock-result. The two-determinantal approach with two optimized basis functions of type _g and _u nearly gives the optimized double configuration SCF result.     

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.     

Semiempirical local spin: Theory and implementation of the ZILSH method for predicting Heisenberg exchange constants of polynuclear transition metal complexes

The local spin formalism ( 3 ) for computing expectation values 〈S_A · S_B〉 that appear in the Heisenberg spin model has been extended to semiempirical single determinant wave functions. An alternative derivation of expectation values in restricted and unrestricted cases is given that takes advantage of the zero differential overlap (ZDO) approximation. A formal connection between single determinant wave functions (which are not in general spin eigenfunctions) and the Heisenberg spin model was established by demonstrating that energies of single determinants that are eigenfunctions of the local spin operators with eigenvalues corresponding to high‐spin radical centers are given by the same Heisenberg coupling constants {J_AB} that describe the true spin states of the system. Unrestricted single determinant wave functions for transition metal complexes are good approximations of local spin eigenfunctions when the metal d orbitals are local in character and all unpaired electrons on each metal have the same spin (although spins on different metals might be reversed). Good approximations of the coupling constants can then be extracted from local spin expectation values 〈S_A · S_B〉 energies of the single determinant wave functions. Once the coupling constants are obtained, diagonalization of the Heisenberg spin Hamiltonian provides predictions of the energies and compositions of the spin states. A computational method is presented for obtaining coupling constants and spin‐state energies in this way for polynuclear transition metal complexes using the intermediate neglect of differential overlap Hamiltonian parameterized for optical spectroscopy (INDO/S) in the ZINDO program. This method is referred to as ZILSH, derived from ZINDO, Davidson's local spin formalism, and the Heisenberg spin model. Coupling constants and spin ground states obtained for 10 iron complexes containing from 2 to 6 metals are found to agree well with experimental results in most cases. In the case of the complex [Fe₆O₃(OAc)₉(OEt)₂(bpy)₂]⁺, a priori predictions of the coupling constants yield a ground‐state spin of zero, in agreement with variable‐temperature magnetization data, and corroborate spin alignments proposed earlier on the basis of structural considerations. This demonstrates the potential of the ZILSH method to aid in understanding magnetic interactions in polynuclear transition metal complexes. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

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.     

Singer polymals. I. Evaluation of matrix elements

Hamiltonian matrix elements between pairs of explicitly correlated Singer polymals are derived. A new method for 1/r₁₂ matrix elements is found which avoids the 6 by 6 matrix inversion required at each quadrature point in the old method. The new formula requires about half the computer time of the old one. The electron repulsion and nuclear attraction matrix elements are shown to be related to incomplete elliptic integrals. Simplified formulas are derived for the special cases of cylindrical symmetry, spherical symmetry, and single-center expansions. An upper bound approximation formula for the purpose of neglecting matrix elements is derived. These bounds may be useful for single-center expansions but are found to be too high for the general case.     

The evaluation of molecular electron affinities

The performance of a variety of levels of theory in evaluating molecular electron affinities (EAs) has been systematically examined. Calculations have been carried out for six different basis sets and for nine theoretical procedures including unrestricted (UHF) and restricted (RHF) Hartree-Fock theory, Møler-Plesset perturbation theory (UMP2, UMP3, UMP4), configuration interaction (UCISD, RCISD, RCISD(Q)) and equations-of-motion (EOM) approaches. Electron affinities were evaluated for CH₃, NH₂, OH, F, C₂H, CN, BO, N₃, OCN, and NO₂. Very poor results are generally obtained unless diffuse functions are included in the basis set and electron correlation is incorporated. Even with the largest basis set used in the present study (6-311 + + G(2d, 2p)), there are still residual errors greater than 0.2 eV (UMP4) or 0.6 eV (CISD) in the absolute EAs. However, better results are obtained under certain circumstances for relative EAs. The results appear to be significantly affected by spin contamination in the UHF wave-functions. For those systems for which spin contamination is small, best absolute values of the EAs generally come from the EOM and UMP2 calculations, whereas the most constant errors (thereby allowing systematic correction) are found at the UMP4, CISD, and RCISD(Q) levels. For the systems for which spin contamination is larger, best results are obtained with the CI-based procedures (CISD and RCISD(Q)). The errors in calculated EAs for the molecules with differing electronic characteristics can vary quite widely. Caution must therefore be exercised before applying schemes which rely on a constancy of errors to estimate electron affinities. The UMP procedures appear particularly suspect in this regard if spin contamination is significant. The RCISD(Q) approach is recommended under such circumstances.     

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.     

Unitary group approach to general system partitioning. I. Calculation of U(n = n1 + n2): U(n1) × u(n2) reduced matrix elements and reduced wigner coefficients

This paper is the first in a series of two directed toward a unitary calculus for group-function-type approaches to the many-electron correlation problem. In this paper we present a complete derivation of the matrix elements of the U(n = n₁ + n₂) generators, for the representations approapriate to many-electron systems, in a basis symmetry adapted to the subgroup U(n₁) × U(n₂). Explicit formulae for the fundamental U(n):U(n₁) × U(n₂) reduced Wigner coefficients, which are needed for the general multishell problem, are also obtained. The symmetry properties of the reduced Wigner coefficients and reduced matrix elements are investigated, and a suitable phase convention is given.     

He++N2 Charge Transfer Reaction: An Ab initio Analysis

An ab initio analysis on the involved potential energy surfaces is presented for the investigation of the charge transfer mechanism for the He⁺+N₂ system. At high collision energy, as many as seven low-lying electronic states are observed to be involved in the charge transfer mechanism. Potential energy surfaces for these low-lying electronic states have been computed in the Jacobi scattering coordinates, applying multireference configuration interaction level of theory and aug-cc-pVQZ basis sets. Asymptotes for the ground and various excited states are assigned to mark the entrance (He⁺+N₂) and charge transfer channels (He+N₂⁺). Nonadiabatic coupling matrix elements and quasi-diabatic potential energy surfaces have been computed for all seven states to rationalize the available experimental data on the charge transfer processes and to facilitate dynamics studies.     

Improved treatment for matrix elements of spin-dependent operators

In this paper a general method for the evaluation of the matrix elements of spin-dependent operators is proposed to improve the treatment primitively suggesteed by Cooper and Musher. This approach is largely based on the recent results which the present authors have achieved in the representation theory for the inner- and outer-product reduction of the symmetric group. It is shown that the so-called outer-product coupling coefficients (OPCC ) can be used to generalize the method for constructing the irreducible tensor operators of group S_n. Together with the use of inner-product coupling coefficients (IPCC ), an expression for the matrix elements of spin-dependent operators is presented as the product of a Racah coefficient for S_n and a reduced matrix element which can be expressed in terms of IPCC, OPCC , and the related integrals. The treatment for one- and two-electron spin-dependent operators is discussed in detail.     

Spin-Orbit coupling in aromatic hydrocarbons: A semiempirical evaluation of the radiative lifetimes of naphtalene and anthracene

Singlet-triplet spin-orbit matrix elements, which govern the lowest ³B_1u ← ¹A_g transition in typical aromatic molecules like naphtalene and anthracene, are calculated with INDO molecular orbitals and the conventional spin-orbit one-electron hamiltonian. The correct order of magnitude of the triplet radiative lifetimes is obtained for the two molecules, when INDO MO coefficients are referred to a symmetrically orthogonalized basis. The possibility of using the semiempirical Hamiltonian is explored using ab initio wave-functions for a few test cases. Reasonably accurate doublet-doublet and singlet-triplet matrix elements have been computed.