Unitary group approach to the many-electron problem. III. Matrix elements of spin-dependent Hamiltonians

This is the final paper in a series of three directed toward the evaluation of spin-dependent Hamiltonians. In this paper we derive the reduced matrix elements of the U(2n) generators in a basis symmetry adapted to the subgroup U(n) × U(2) (i.e., spin-orbit basis), for the representations appropriate to many-electron systems. This enables a direct evaluation of the matrix elements of spin-dependent Hamiltonians in the spin-orbit basis. An alternative (indirect) method, which employs the use of U(2n) ↓ U(n) × U(2) subduction coefficients, is also discussed.     

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.     

A spin-dependent unitary group approach to many-electron systems

In this article we derive a segment-level formula for the matrix elements of the U(2n) generators in a basis symmetry adapted to the subgroup U(n) × U(2) (i.e., spin-orbit basis), for the representations appropriate to many-electron systems. This enables the direct evaluation of the matrix elements of spin-dependent Hamiltonians.     

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 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.     

Subduction coefficients for 〈2N/2–S, 12S〉 ↓ 〈2, 1〉 ⊗ 〈2, 1〉 of U(n) ↓ U(n1) ⊗ U(n2)

The determination of the subduction coefficients for states of the unitary group U(n) under the restrictions U(n) ↓ U(n₁) ? U(n₂) have been considered for the spin free states of many electron systems. Using the transformation properties of the tensor basis spanning the irreducible representation 〈2^N/2–S, 1^2S〉 of U(n) under the permutations of electron coordinates, a simple programmable procedure has been developed for the determination of these coefficients. The procedure has been illustrated using a simple example.     

Matrix elements of the U(2n) generators in the spin-orbit basis

An algorithm for the evaluation of the matrix elements of the U(2n) generators in the spin-orbit basis [induced by the U(n) ? U(2) subgroup] is proposed. This algorithm can be used for calculations with spin-dependent Hamiltonians for many-electron problems.     

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.     

The evaluation of matrix elements for non-canonical Weyl tableau basis states adapted toU(n 1+n 2)⊃U(n 1)×U(n 2)

Summary This is the first paper of a series of two, which enables the evaluation ofU(n) generator matrix elements in the non-canonical Weyl tableau basis adapted to subgroupU(n ₁)×U(n ₂). In this paper the explicit closed formulae for subduction coefficients are presented. These formulae will become useful through an inductive method to be presented in the second paper.     

Unitary group tensor operator algebras for many-electron systems. III. Matrix elements in U(n 1 +n 2) ⊃ U(n 1) × U(n 2) partitioned basis

Exploiting our earlier results [J. Math. Chem. 4 (1990) 295–353 and 13 (1993) 273–316] on the unitary group U(n) Racah-Wigner algebra, specifically designed for quantum chemical calculations of molecular electronic structure, and the related tensor operator formalism that enabled us to introduce spin-free orbital equivalents of the second quantization-like creation and annihilation operators as well as higher rank symmetric, antisymmetric and adjoint tensors, we consider the problem of U(n) basis partitioning that is required for group-function type approaches to the many-electron problem. Using the U(n) U(n ₁) × U(n ₂),n =n ₁ +n ₂ adapted basis, we evaluate all required matrix elements of U(n) generators and their products that arise in one- and two-body components of non-relativistic electronic Hamiltonians. The formalism employed naturally leads to a segmented form of these matrix elements, with many of the segments being identical to those of the standard unitary group approach. Relationship with similar approaches described earlier is briefly pointed out.     

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.     

Unitary group approach to general system partitioning. II. U(n) matrix element evaluation in a composite basis

An explicit segment level formalism is derived for the matrix elements of the U(n) generators in an arbitrary (multishell) composite basis. The results of this paper, which contain the usual (spin-independent) unitary calculus approach as a limiting case, yield a more powerful and versatile algorithm than the traditional (spin-independent) unitary group formalism.     

Some simplification in spin-free configuration interaction studies

A simplification has been attempted in the procedures for determining the matrix elements of the generators of the unitary group U(n) over a tensor basis spanning the irreducible representation 2 ^N/2–S , 1 ^2S for an N-electron system. It has been shown that these matrix elements require, for their determination, only the corresponding representation matrices of cyclic permutations of the group S _N . A viable algorithm has been obtained for determining these representation matrices.     

Spin-free quantum chemistry. XXIII. The generator-state approach

In the unitary-group formulation of quantum chemistry, the spin-projected, configuration-state spaces of quantum chemistry are realized by the irreducible representation spaces (IRS ) of the freeon unitary group U_(n), where n is the number of freeon orbitals. The Pauli-allowed IRS are labeled by the partitions [λ] = [2^(N/2)?s, 1^2S], where N and S are the particle number and the spin, respectively. The generator-state approach (GSA ) to the unitary-group formulation consists of (1) the construction of the overcomplete, nonorthonormal generator basis for each IRS ; (2) the Lie-algebraic computation of matrix elements over generator states; (3) the Moshinsky–Nagel construction of the complete, orthonormal Gel'fand basis in terms of the generator basis; and (4) the computation of matrix elements over Gel'fand states in terms of matrix elements over generator states.     

Calculation of SN 1 + N 2 ⊃ SN1 ⊗ sn 2 and U(n1 + n2) ⊃ u (n1) ⊗ u (n2) subduction coefficients by using spin graph

By use of the graphical method of spin algebra, simple, and closed expressions for S_N1+_N2 ? S_N1 ? S_N2 and U(n₁ + n₂) ? U(n₁) ? U(n₂), subduction coefficients are derived.     

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.     

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.     

Method of recurrent construction of Löuwdin spin-adapted wave functions. II. Local representation of creation–annihilation operators

Some compositions of the addition and subtraction operators and recurrence relations for the Sanibel-type coefficients c_{u, v} (n, s, M) generated by these compositions are studied. A local representation of the fermion creation–annihilation operators via the addition and subtraction operators is obtained. Operators of single excitations, coupling, and decoupling operators, in terms of which the unitary group generators can be expressed are defined. The resulting representation of the nonelementary unitary group generators is much more simple than in the Gelfand–Tzetlin basis and in the most general case contains only six logically different terms, each of them possessing quite transparent physical significance.     

A unitary group formulation of many-body theory: Diagram systematics and use of the spin shifts

This paper shows that the spin-shift formalism developed in B. T. Pickup and A. Mukhopadhyay [Int. J. Quantum Chem. 26 , 101 (1984)] supports a one-component diagrammatics which has a systematics akin to that in the spin-orbital many-body theory. The diagrams are neither Goldstone nor Yutsis type, and characterize the chain U(2R) ? U(R)?SU(2) on which the spin-shift formalism is based. Accordingly, while the lines in such diagrams are labeled by the orbital indices, the diagram structure adequately reflects the irreducible representation of the group U(R). In this sense the paper presents a unitary group approach to the natural generalization of the usual many-body theory for the spin-adapted cases. A set of very simple rules is derived; their similarity with the corresponding rules in the ordinary many-body theory and practical utility are discussed in connection with (a) matrix elements over many-electron spin states and (b) closed- and open-shell many-body perturbation theory. A possibility of integral-driven many-body perturbation theory for open-shells is indicated. Connections of this formalism with others are also discussed.     

Unitary group approach to spin-adapted open-shell coupled cluster theory

We show that the irreducible tensor operators of the unitary group provide a natural operator basis for the exponential Ansatz which preserves the spin symmetry of the reference state, requires a minimal number of independent cluster amplitudes for each substitution order, and guarantees the invariance of the correlation energy under unitary transformations of core, open-shell, and virtual orbitals. When acting on the closed-shell reference state with n_c doubly occupied and n_v unoccupied (virtual) orbitals, the irreducible tensor operators of the group U(n_c) ? U(n_V) generate all Gelfand-Tsetlin (GT) states corresponding to appropriate irreducible representation of U(n_c + n_v). The tensor operators generating the M-tuply excited states are easily constructed by symmetrizing products of M unitary group generators with the Wigner operators of the symmetric group S_M. This provides an alternative to the Nagel-Moshinsky construction of the GT basis. Since the corresponding cluster amplitudes, which are also U(n_c) ? U(n_s) tensors, can be shown to be connected, the irreducible tensor operators of U(n_c) ? U(n_v) represent a convenient basis for a spin-adapted full coupled cluster calculation for closed-shell systems. For a high-spin reference determinant with n, singly occupied open-shell orbitals, the corresponding representation of U(n), n=n_c + n_v + n_s is not simply reducible under the group U(n_c) ? U(n_s) ? U(n_v). The multiplicity problem is resolved using the group chain U(n) ? U(n_c + n_v) ? U(n_s) ? U(n_c) ?U(n_s)? U(n_v) ? U(n_v). The labeling of the resulting configuration-state functions (which, in general, are not GT states when n_c > 1) by the irreducible representations of the intermediate group U(n_c + n_v) ?U(n_s) turns out to be equivalent to the classification based on the order of interaction with the reference state. The irreducible tensor operators defined by the above chain and corresponding to single, double, and triple substitutions from the first-, second-, and third-order interacting spaces are explicitly constructed from the U(n) generators. The connectedness of the corresponding cluster amplitudes and, consequently, the size extensivity of the resulting spin-adapted open-shell coupled cluster theory are proved using group theoretical arguments. The perturbation expansion of the resulting coupled cluster equations leads to an explicitly connected form of the spin-restricted open-shell many-body perturbation theory. Approximation schemes leading to manageable computational procedures are proposed and their relation to perturbation theory is discussed. © 1995 John Wiley & Sons, Inc.