Vibrational coupled cluster theory

The theory and first implementation of a vibrational coupled cluster (VCC) method for calculations of the vibrational structure of molecules is presented. Different methods for introducing approximate VCC methods are discussed including truncation according to a maximum number of simultaneous mode excitations as well as an interaction space order concept is introduced. The theory is tested on calculation of anharmonic frequencies for a three-mode model system and a formaldehyde quartic force field. The VCC method is compared to vibrational self-consistent-field, vibrational M?ller-Plesset perturbation theory, and vibrational configuration interaction (VCI). A VCC calculation typically gives higher accuracy than a corresponding VCI calculation with the same number of parameters and the same formal operation count.     

Vibrational excitation energies from vibrational coupled cluster response theory

Response theory in the context of vibrational coupled cluster (VCC) theory is introduced and used to obtain vibrational excitation energies. The relation to the vibrational configuration interaction (VCI) approach is described, and the increase in accuracy of VCC response energies relative to VCI energies is discussed theoretically in terms of a perturbational order expansion and demonstrated numerically. To illustrate the theory, a pilot implementation is used to obtain anharmonic vibrational frequencies for fundamental, first overtone and combination excitations of formaldehyde as well as for the fundamental transitions of ethylene.     

Quasi-variational coupled cluster theory

We extend our previous work on the construction of new approximations of the variational coupled cluster method. By combining several linked pair functional transformations in such a way as to give appropriately balanced infinite-order contributions, in order to approximate (L) well at all orders, we formulate a new quantum chemical method, which we name quasi-variational coupled cluster. We demonstrate this method to be particularly robust in the regime of strong static electron correlation, improving significantly on our earlier approximate variational coupled cluster approach.     

Coupling term derivation and general implementation of state-specific multireference coupled cluster theories

Simple closed-form expressions are derived for the "same vacuum" renormalization terms that arise in state-specific multireference coupled cluster (MRCC) theories. Explicit equations are provided for these coupling terms through the triple excitation level of MRCC theory, and a general expression is included for arbitrary-order excitations. The first production-level code (PSIMRCC) for state-specific and rigorously size-extensive Mukherjee multireference coupled cluster singles and doubles (MkCCSD) computations has been written. This code is also capable of evaluating analogous Brillouin-Wigner multireference energies (BWCCSD), including a posteriori size-extensivity corrections. Using correlation-consistent basis sets (cc-pVXZ, X=D,T,Q), MkCCSD and BWCCSD were tested and compared on two classic multireference problems: (1) the dissociation potential curve of molecular fluorine (F(2)) and (2) the structure and vibrational frequencies of ozone. Comparison with experimental data shows that the Mukherjee method is generally superior to the Brillouin-Wigner theory in predicting energies, structures, and vibrational frequencies. Particularly accurate results for F(2) are obtained by applying the MkCCSD method with localized molecular orbitals. Although the MkCCSD theory greatly improves upon single-reference CCSD for the geometric parameters and a(1) vibrational frequencies of ozone, the antisymmetric stretching frequency omega(3)(b(2)) remains pathological and cannot be properly treated without the inclusion of connected triple excitations. Finally, preliminary multireference MkCCSD results are reported for the singlet-triplet splittings in ortho-, meta-, and para-benzyne, coming within 1.5 kcal mol(-1) of experiment in all cases.     

无自旋混杂的耦簇理论

本文用APCCSD(T)，自旋投影方法消除CCSD中的自旋混杂，自旋投影算符^P用自旋湮灭算符^As 1和^As 2的乘积近似表出，提高了计算精度。     

Approximate variational coupled cluster theory

We show that it is possible to construct an accurate approximation to the variational coupled cluster method, limited to double substitutions, from the minimization of a functional that is rigorously extensive, exact for isolated two-electron subsystems and invariant to transformations of the underlying orbital basis. This approximate variational coupled cluster theory is a modification and enhancement of our earlier linked pair functional theory. It is first motivated by the constraint that the inverse square root of the matrix that transforms the cluster amplitudes must exist. Low-order corrections are then included to enhance the accuracy of the approximation of variational coupled cluster, while ensuring that the computational complexity of the method never exceeds that of the standard traditional coupled cluster method. The effects of single excitations are included by energy minimization with respect to the orbitals defining the reference wavefunction. The resulting quantum chemical method is demonstrated to be a robust approach to the calculation of molecular electronic structure and performs well when static correlation effects are strong.     

Dispersion coefficients for first hyperpolarizabilities using coupled cluster quadratic response theory

The frequency dependence of third-order properties can in the normal dispersion region be expanded in a Taylor series in the frequency arguments. The dispersion coefficients thus obtained provide an efficient way of expressing the dispersion of frequency-dependent properties and are transferable between different optical processes. We derive analytic expressions for the dispersion coefficients of third-order properties in coupled cluster quadratic response theory and report an implementation for the three coupled cluster models CCS, CC2, and CCSD. Calculations are performed for the first hyperpolarizability of the NH₃ molecule. The convergence of the dispersion expansion with the order of the coefficients is examined and we find good convergence up to about half the frequency at which the first pole in the hyperpolarizability occurs. Padé approximants improve the convergence dramatically and extend the application range of the dispersion expansion to frequencies close to the first pole. The sensitivity of the dispersion coefficients on the dynamic correlation treatment and on the choice of the one-electron basis set is investigated. The results demonstrate that, contrary to presumptions in the literature, the dispersion coefficients are sensitive to basis set effects and correlation treatment similar to the static hyperpolarizabilities. Received: 26 March 1998 / Accepted: 21 July 1998 / Published online: 19 October 1998     

Brueckner coupled cluster response functions

A derivation of the linear response function for the Brueckner coupled cluster method is presented that enables the calculation of second-order molecular properties such as frequency-dependent polarizabilities. By using the Brueckner orbital variant of coupled cluster theory, the spurious pole structure inherent in the standard coupled cluster approach with orbital relaxation is avoided. © 1994 John Wiley & Sons, Inc.     

Localized optimized orbitals, coupled cluster theory, and chiroptical response properties

The impact of orbital localization on the efficiency and accuracy of the optimized-orbital coupled cluster model is examined for the prediction of chiroptical properties, in particular optical rotation. The specific rotations of several test cases-(P)-[4]triangulane, (S)-1-phenylethanol, and chiral conformers of 1-fluoropentane, heptane, and nonane-were computed using an approach in which localization is enforced throughout the orbital optimization and subsequent linear response computation. This method provides a robust local-correlation scheme for future production-level implementation. Although the cross-over point between the canonical and localized coupled cluster approach lies at larger molecules than for ground-state energies, the scheme presented should still provide reduced scaling sufficient to investigate much larger molecules than are presently accessible.     

Block-correlated coupled cluster theory: the general formulation and its application to the antiferromagnetic Heisenberg model

The general formalism of the block-correlated coupled cluster (BCCC) method, an alternative multireference coupled cluster method for calculating the ground-state electronic structures of molecular systems, has been presented. The BCCC theory is constructed in terms of a complete set of many-electron states of individual blocks, assumed that the whole system could be partitioned into a set of blocks. The reference state in the BCCC is selected as a tensor product of the most important many-electron state of each system block. By truncating the cluster operator to a certain n-block correlation level, an approximate but size-extensive BCCC method, denoted as BCCCn, is defined. For reducing the computational effort but without much loss of accuracy, the reduced density matrix is introduced to generate an optimal subset of many-electron states for each block. I have implemented the BCCCn (n=2,3) methods within the S=1/2 Heisenberg Hamiltonian, and applied them to calculate the ground-state energies of one-dimensional spin chains and quasi-one-dimensional two-leg spin ladders. The calculated results show that with the appropriate partition of the studied systems the BCCC3 method can yield quite satisfactory ground-state energies for these spin systems.     

Vibrational absorption spectra from vibrational coupled cluster damped linear response functions calculated using an asymmetric Lanczos algorithm

We report the theory and implementation of vibrational coupled cluster (VCC) damped response functions. From the imaginary part of the damped VCC response function the absorption as function of frequency can be obtained, requiring formally the solution of the now complex VCC response equations. The absorption spectrum can in this formulation be seen as a matrix function of the characteristic VCC Jacobian response matrix. The asymmetric matrix version of the Lanczos method is used to generate a tridiagonal representation of the VCC response Jacobian. Solving the complex response equations in the relevant Lanczos space provides a method for calculating the VCC damped response functions and thereby subsequently the absorption spectra. The convergence behaviour of the algorithm is discussed theoretically and tested for different levels of completeness of the VCC expansion. Comparison is made with results from the recently reported [P. Seidler, M. B. Hansen, W. Gyo?rffy, D. Toffoli, and O. Christiansen, J. Chem. Phys. 132, 164105 (2010)] vibrational configuration interaction damped response function calculated using a symmetric Lanczos algorithm. Calculations of IR spectra of oxazole, cyclopropene, and uracil illustrate the usefulness of the new VCC based method.     

Multireference coupled cluster theory in Fock space

Summary A Fock space multireference coupled cluster method based on incomplete model spaces is described. Some of the essential computational aspects of the theory are discussed with the aid of the diagrammatic representation of the equations. An application to the calculation of ionization potentials and excitation energies ofs-tetrazine is presented along with comparisons with conventionalab initio calculations and experimental results.     

The Fock space coupled cluster method: theory and application

Summary The Fock space coupled cluster method and its application to atomic and molecular systems are described. The importance of conserving size extensivity is demonstrated by the electron affinities of the alkali atoms. Two types of intruder states are discussed, one attributable to the orbital energy spectrum and the other caused by two-electron interactions. They are illustrated by the excited states of Li₂ and by¹ S states of Be, respectively. It is shown how both problems may be solved using incomplete model spaces. The selection of the model space in a moderately dense spectrum is discussed in connection with N₂ excited states.Supported in part by the U.S.-Israel Binational Science Foundation     

Method of moments approach and coupled cluster theory

Summary The single reference coupled cluster (CC) approach to the many-electron correlation problem is examined from the viewpoint of the method of moments (MM). This yields generally an inconsistent (overcomplete) set of equations for cluster amplitudes, which can be solved either in the least squares sense or by selective projection process restricting the number of equations to that of the unknowns. These resulting generalized MM-CC equations always contain the standard CC equations as a special case. Since, in the MM-CC formalism, the Schrödinger equation will be approximately satisfied on a subspace spanned by non-canonical configurations, this procedure may be helpful in extending the standard single reference CC theory to quasi-degenerate situations. To examine the potential usefulness of this idea, we explore the linear version of the CC approach for systems with a quasi-degenerate reference, in which case the standard linear theory is plagued with singularities due to the intruder states. Implications of this analysis for the structure of the wavefunction are also briefly discussed.Killam Research Fellow 1987–89     

Two-particle density matrix cumulant of coupled cluster theory

A two-particle density matrix obtained from the expectation value expression for the coupled cluster wave function is separated into the antisymmetrized product of the one-particle density matrices and the remaining cumulant part. It is demonstrated that the proposed formula for the coupled cluster cumulant is a valid cumulant expression, since it is a connected, and therefore size-extensive quantity. It is also shown that the density matrices from coupled cluster gradient theory, when used to define a cumulant, result in the expression containing disconnected terms. The proposed formulation of the coupled cluster cumulant makes it easy to develop size-extensive truncation schemes. As an example, explicit equations for the cumulant at the coupled cluster single, double and triple excitation level are presented.     

Thermochemistry of the fluoroformyloxyl radical: a computational study based on coupled cluster theory

The standard enthalpy of formation of FCO(2) (X (2)B(2)) was determined by a computational approach based on coupled cluster theory [CCSD(T)] with energies extrapolated to the basis-set limit, with additional corrections accounting for core-valence correlation, scalar relativity, spin-orbit coupling, and zero-point vibrational motions. Utilizing a variety of independent reaction schemes, our best estimate is Delta(f)H(o)(0)(FCO(2)) = -86.0 +/- 0.6 kcal mol(-1) [Delta(f)H(o)(298) )(FCO(2)) = -86.7 +/- 0.6 kcal mol(-1)], which is shown to be more accurate than previous theoretical and experimental values. The chosen computational procedure was also applied to HCO (X (2)A'), where we find excellent agreement with experiment, and to FCO (X (2)A'), where we recommend an improved value of Delta(f)H(o)(0)(FCO) = -42.1 +/- 0.5 kcal mol(-1) [ Delta(f)H(o)(298)(FCO) = -42.0 +/- 0.5 kcal mol(-1)]. Further theoretical results concern the C-F bond dissociation energy, electron affinity, ionization energy, first and second excitation energies in FCO(2), fluoride ion affinity of CO(2), and equilibrium geometries of the molecules treated presently. For FCO (X (2)A') we propose an improved equilibrium structure: r(e)(CF) = 132.5(2) pm, r(e)(CO) = 116.7(2) pm, and theta(e)(FCO) = 127.8(2)(o).     

Fock space multireference coupled cluster theory: Study of shape resonance

The complex absorbing potential along with correlated independent particle potential (CIP) Fock space multireference coupled cluster method is used for the study of resonances. We have studied shape resonance of e^?‐ F₂, e^?‐ N₂O and e^?‐CO molecules. In particular, we have studied e^?‐ F₂ scattering at different bond lengths to know whether is bound at the equilibrium bond length of F₂. © 2013 Wiley Periodicals, Inc.     

Formulation and implementation of a unitary group adapted state universal multi-reference coupled cluster (UGA-SUMRCC) theory: Excited and ionized state energies

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.     

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.     

Coupled cluster response theory in parameter subspaces

We introduce a response function formalism that enables smaller number of parameters than that defined in standard coupled cluster response theory. This is essential in the development of reduced scaling methods. The formalism is general and it applies to all parameterizations at all levels of the coupled cluster hierarchy. We show that to achieve physically reasonable results the parameterization must fulfill certain criteria. The linear response functions are derived and discussed in the context of optimized virtual orbitals and Cholesky decomposition of the cluster amplitudes. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009