We present a general-order spin-free formulation of the single-reference closed-shell coupled-cluster method. We show that the working equations of a fully biorthogonal contravariant projection formulation of the residual equations, as near-universally used in closed-shell CCSD, can also be defined at the CCSDT and CCSDTQ levels, despite singularities in the spin projection manifolds. We describe permutation-group based techniques for obtaining and simplifying the equations encountered in general second-quantization-based methods; this includes a permutation group based approach of evaluating second-quantized matrix elements into tensor contraction networks, and the use of Portugal's double coset canonical representation technique [Int. J. Mod. Phys. C 13, 859 (2002)] for eliminating redundant terms. A computer implementation of our techniques is simple, because no operator-valued symbolic algebra is required. Explicit working equation lists for closed-shell CCSD, CCSDT, and CCSDTQ in the semi-biorthogonal formulation are provided. We also release open-source computer programs for both deriving and numerically evaluating these equations.
We present a general-order spin-free formulation of the single-reference closed-shell coupled-cluster method. We show that the working equations of a fully biorthogonal contravariant projection formulation of the residual equations, as near-universally used in closed-shell CCSD, can also be defined at the CCSDT and CCSDTQ levels, despite singularities in the spin projection manifolds. We describe permutation-group based techniques for obtaining and simplifying the equations encountered in general second-quantization-based methods; this includes a permutation group based approach of evaluating second-quantized matrix elements into tensor contraction networks, and the use of Portugal's double coset canonical representation technique [Int. J. Mod. Phys. C 13, 859 (2002)] for eliminating redundant terms. A computer implementation of our techniques is simple, because no operator-valued symbolic algebra is required. Explicit working equation lists for closed-shell CCSD, CCSDT, and CCSDTQ in the semi-biorthogonal formulation are provided. We also release open-source computer programs for both deriving and numerically evaluating these equations.
We present a symbolic manipulation algorithm for the efficient automated implementation of rigorously spin-free coupled cluster (CC) theories based on a unitary group parametrization. Due to the lack of antisymmetry of the unitary group generators under index permutations, all quantities involved in the equations are expressed in terms of non-antisymmetric tensors. Given two tensors, all possible contractions are first generated by applying Wick's theorem. Each term is then put down in the form of a non-antisymmetric Goldstone diagram by assigning its contraction topology. The subsequent simplification of the equations by summing up equivalent terms and their factorization by identifying common intermediates is performed via comparison of these contraction topologies. The definition of the contraction topology is completely general for non-antisymmetric Goldstone diagrams, which enables our algorithm to deal with noncommuting excitations in the cluster operator that arises in the unitary group based CC formulation for open-shell systems. The resulting equations are implemented in a new code, in which tensor contractions are performed by successive application of matrix-matrix multiplications. Implementation of the unitary group adapted CC equations for closed-shell systems and for the simplest open-shell case, i.e., doublets, is discussed, and representative calculations are presented in order to assess the efficiency of the generated codes.
This thesis deals with selected aspects of a new implementation of many-body methods which can be formulated in the framework of second quantization, in particular the coupled-cluster (CC) method with arbitrary excitation level. Coupled-cluster is one of the most successful and widely used quantum chemical methods for accurate calculations on small to medium-sized molecules. Since it employs a nonlinear parametrization of the wave function, its implementation is rather difficult, in particular if higher (i.e. more than double) excitations are to be included. The latter is necessary to obtain highly accurate results and also as a prerequisite for the generalization to multi-reference cases. The implementation described here has a twofold focus. One is on generality and flexibility regarding the method to be implemented, the other is on efficiency. To achieve flexibility, it is useful to have a machinery which automatically derives working equations for a given method. We realize this by applying techniques of second quantization. This work treats in particular the last step of this procedure, namely the simplification of the resulting equations by the identification of equivalent terms. The algorithm used here is based on the interpretation of algebraic terms as graphs. The derived CC equations then have to be solved iteratively. The efficiency of the program is mainly determined by the evaluation of the occurring expressions, which has to be done in each iteration step. The evaluation is split up in a sequence of tensor contractions. Their generic implementation is complicated by the particular structure of the involved tensors. We reduce each contraction to a sequence of matrix multiplications, which requires a previous data rearranging. But since matrix multiplication is the most efficient operation on modern computers, this additional effort pays off. Preliminary tests show that our program is at least as fast as the most efficient general coupled-cluster implementation so far, and the relation is expected to improve for calculations with larger basis sets where the matrix multiplication becomes the time-determining step. Finally, we give an outlook to possible further developments. In particular, the preparation of the equations before the actual evaluation (factorization) offers much potential for optimizations which we do not exploit at the moment, in contrast to the program with which we compare, which employs at least a partial optimization at this point.
The Coupled-Cluster (CC) method is one of the most popular and efficient correlation methods in quantum chemistry. \nEspecially the CCSD(T) approximation, which includes single and double excitations by means of the application of the cluster operator \nto a reference determinant and triple excitations via a perturbative treatment, has become a standard tool \nin quantum chemical applications. However, the method is restricted to relatively small system sizes due to its unfavourable scaling (N^7, where N is the number of basis functions applied). \nFor the full treatment of triple excitations the scaling advances to N^8 and every further excitation level increases the exponent by two. \n \nThe goal of this work is to reduce the calculation time for closed shell systems at least by a factor growing with the excitation level for arbitrary truncation levels. \nThis is done by restricting the spatial parts of the spin orbitals and thus treat pairs of spin orbitals on the same footing. The restrictions can be easily constructed \nfor the CCSD model by spin integration. \nThe derivation of the restrictions arising in higher excited case will be done employing \nthe spatial orbital excitation operators E. \n \nIn the first part of this work an algorithm is derived that is capable of the derivation of the energy and amplitude equations for arbitrary excitation levels. In the second part \nan implementation of this algorithm is presented.
The traditional state universal multi-reference coupled cluster (SUMRCC) theory uses the Jeziorski-Monkhorst (JM) based Ansatz of the wave operator: Ω = Σ(μ)Ω(μ)|φ(μ)><φ(μ)|, where Ω(μ) = exp(T(μ)) is the cluster representation of the component of Ω inducing virtual excitations from the model function φ(μ). In the first formulations, φ(μ)s were chosen to be single determinants and T(μ)s were defined in terms of spinorbitals. This leads to spin-contamination for the non-singlet cases. In this paper, we propose and implement an explicitly spin-free realization of the SUMRCC theory. This method uses spin-free unitary generators in defining the cluster operators, {T(μ)}, which even at singles-doubles truncation, generates non-commuting cluster operators. We propose the use of normal-ordered exponential parameterization for Ω:Σ(μ){exp(T(μ))}|φ(μ)><φ(μ)|, where {} denotes the normal ordering with respect to a common closed shell vacuum which makes the "direct term" of the SUMRCC equations terminate at the quartic power. We choose our model functions {φ(μ)} as unitary group adapted (UGA) Gel'fand states which is why we call our theory UGA-SUMRCC. In the spirit of the original SUMRCC, we choose exactly the right number of linearly independent cluster operators in {T(μ)} such that no redundancies in the virtual functions {χ(μ)(l)} are involved. Using example applications for electron detached/attached and h-p excited states relative to a closed shell ground state we discuss how to choose the most compact and non-redundant cluster operators. Although there exists a more elaborate spin-adapted JM-like ansatz of Datta and Mukherjee (known as combinatoric open-shell CC (COS-CC), its working equations are more complex. Results are compared with those from COS-CC, equation of motion coupled cluster methods, restricted open-shell Hartree-Fock coupled cluster, and full configuration interaction. We observe that our results are more accurate with respect to most other theories as a result of the use of the cluster expansion structure for our wave operator. Our results are comparable to those from the more involved COS-CC, indicating that our theory captures the most important aspects of physics with a considerably simpler scheme.
Coupled-cluster many-electron theory (CCMET) has been developed in the second-quantized formulation with and without spin projection. In this article CCMET is developed in the unitary-group formulation where the group is U(p) and p is the number of orbitals in the basis set. The zero-order ground state is the highest-weight state of an irreducibly invariant subspace of U(p) and the excitation operators are infinitesimal generators. These irreducible spaces are uniquely labeled by Young diagrams which the Pauli principle limits to no more than one column for fermion orbitals and no more than two columns for freeon (spin-free) orbitals. In the latter case the spin is one-half the difference in the lengths of the two columns. We employ both the Gel'fand and the generator bases for the irreducible invariant vector spaces. Matrix elements for the former are evaluated by techniques due to Gel'fand, Biedenharn, Louck, Paldus, and Shavitt and for the latter directly by Lie algebraic and diagrammatic techniques as a simple function of the weight components of the highest-weight state. For freeon or fermion orbitals the results are equivalent to those obtained by the second-quantized formulation with or without spin projection.