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.     

Orbital-optimized coupled-cluster theory does not reproduce the full configuration-interaction limit

It is shown that due to the mixing of the usual projection approach of coupled cluster with variational orbital optimization, orbital-optimized coupled cluster (OCC) fails to reproduce the full configuration-interaction (full CI) limit when the cluster operator becomes complete. It is pointed out that the fulfillment of the projected singles equations, which define the orbital gradient in Brueckner coupled cluster (BCC), is mandatory for a correct behavior. As numerical examples we present general OCC and BCC calculations up to the full CI limit on CH(2) and an active-space model of ozone. The observed deviations of OCC from full CI are of the order of the correlation error obtained in calculations with up to quadruples excitations. Thus the failure of OCC may be considered tolerable in more approximate calculations but clearly prohibitive for any benchmark application. For applications to active-space models a hybrid approach for OCC is suggested in which for active particle-hole rotations the Brueckner orbital gradient is employed, whereas for the remaining orbital rotations the variational orbital gradient is retained.     

Calculated electronic transitions of the water ammonia complex

We have calculated vertical excitation energies and oscillator strengths of the low lying electronic transitions in H2O, NH3, and H2ONH3 using a hierarchy of coupled cluster response functions [coupled cluster singles (CCS), second order approximate coupled cluster singles and doubles (CC2), coupled cluster singles and doubles (CCSD), and third order approximate coupled cluster singles, doubles, and triples (CC3)] and correlation consistent basis functions (n-aug-cc-pVXZ, where n=s,d,t and X=D,T,Q). Our calculations indicate that significant changes in the absorption spectra of the photodissociative states of H2O and NH3 monomers occur upon complexation. In particular, we find that the electronic transitions originating from NH3 are blueshifted, whereas the electronic transitions originating from H2O are redshifted.     

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.     

Analytical energy gradients for a unitary coupled-cluster theory

Our recently developed and tested unitary multiconfigurational coupled-cluster electronic wavefunction method is extended to permit, for the first time, the analytical evaluation of energy derivatives. The unitary nature of this method admits a variational energy functional whose stationary nature plays a key role in simplifying our derivation. Explicit expressions are given for the gradient (first energy derivative) for both the full unitary coupled cluster and its coupled electron pair approximation (CEPA).     

The quadratic coupled cluster doubles model

Beginning from the bi-variational expression for the standard coupled cluster doubles (CCD) energy, we propose including a term quadratic in the left operator, . As this makes the left-hand and right-hand wave functions more similar, the resulting functional should better approximate the variational expression. The energy given by the quadratic functional is extensive and the stationary equations may be solved in O(N⁶) time. These equations have been implemented and the potential energy surfaces for HF dissociation, H₂O double dissociation and N₂ dissociation are examined. It is found that the quadratic functional effectively reproduces the fully variational results in all cases.     

Some non-relativistic many electron theories: single reference case

There has been substantial progress in the development of electron correlation methods with some benefits and drawbacks. In this paper, we give a review of electron correlation effects on many body particles system. We focus more on atoms and molecules rather than solid state and take into account single reference case only. We mainly discuss perturbation theory, coupled electron many electron theory and a few versions of coupled electron pair approximations with comparison configuration interaction and some coupled cluster methods in which coupled cluster method is crucially important for the crystalline solids, the electron gas and the heat of the reaction. We also show some results, reported by several authors, to fairly compare and judge the methods’ feasibility mentioned above.      

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     

A universal state-selective approach to multireference coupled-cluster non-iterative corrections

A new form of the asymmetric energy functional for multireference coupled cluster (MRCC) theories is discussed from the point of view of an energy expansion in a quasidegenerate situation. The resulting expansion for the exact electronic energy can be used to define the non-iterative corrections to approximate MRCC approaches. In particular, we show that in the proposed framework the essential part of dynamic correlation can be encapsulated in the so-called correlation Hamiltonian, which in analogy to the effective Hamiltonian, is defined in the model space (M(0)). The proper parametrization of the exact/trial wavefunctions leads to the cancellation of the overlap-type numerators and to a connected form of the correlation Hamiltonian and size-extensive energies. Within this parametrization, when the trial wavefunctions are determined without invoking a specific form of the MRCC sufficiency conditions, the ensuing correction can be universally applied to any type of the approximate MRCC method. The analogies with other MRCC triples corrections to MRCC theories with singles and doubles (MRCCSD) are outlined. In particular, we discuss the approach, which in analogy to the Λ-Mk-MRCCSD(T) method [F. A. Evangelista, E. Prochnow, J. Gauss, H. F. Schaefer III, J. Chem. Phys. 132, 074107 (2010)], introduces an approximate form of the triply-excited clusters into the effective and correlation Hamiltonians. Since the discussed corrections can be calculated as a sum of independent reference-related contributions, possible parallel algorithms are also outlined.     

A nine-dimensional calculation of the vibrational OH stretching and HOH bending spectrum of the water trimer

We have studied the vibrational high-frequency spectrum of the water trimer computationally. We expand an earlier study [J. Chem. Phys. A 2009, 113, 9124-9132] where we approximated the water trimer as three individually vibrating water monomer units. Some intramolecular potential energy coupling terms are now included in the previous model. The six OH bond lengths and the three HOH bending angles are used as the internal coordinates. The kinetic energy operator is a sum of the kinetic energy operators of the monomer units. We use the coupled cluster method with single, double, and perturbative triple excitations method [CCSD(T)] with augmented correlation consistent polarized valence triple-ζ (aug-cc-pVTZ) basis set to calculate the potential energy surface (PES). The counterpoise correction is included in the one-dimensional part of the PES. We calculate the vibrational energy eigenvalues using the variational method. The corresponding eigenfunctions are used to obtain the absorption intensities.     

Alternative single-reference coupled cluster approaches for multireference problems: the simpler, the better

We report a general implementation of alternative formulations of single-reference coupled cluster theory (extended, unitary, and variational) with arbitrary-order truncation of the cluster operator. These methods are applied to compute the energy of Ne and the equilibrium properties of HF and C(2). Potential energy curves for the dissociation of HF and the BeH(2) model computed with the extended, variational, and unitary coupled cluster approaches are compared to those obtained from the multireference coupled cluster approach of Mukherjee et al. [J. Chem. Phys. 110, 6171 (1999)] and the internally contracted multireference coupled cluster approach [F. A. Evangelista and J. Gauss, J. Chem. Phys. 134, 114102 (2011)]. In the case of Ne, HF, and C(2), the alternative coupled cluster approaches yield almost identical bond length, harmonic vibrational frequency, and anharmonic constant, which are more accurate than those from traditional coupled cluster theory. For potential energy curves, the alternative coupled cluster methods are found to be more accurate than traditional coupled cluster theory, but are three to ten times less accurate than multireference coupled cluster approaches. The most challenging benchmark, the BeH(2) model, highlights the strong dependence of the alternative coupled cluster theories on the choice of the Fermi vacuum. When evaluated by the accuracy to cost ratio, the alternative coupled cluster methods are not competitive with respect to traditional CC theory, in other words, the simplest theory is found to be the most effective one.     

Explicitly correlated treatment of the Ar-NO+ cation

We present an application of the recently developed explicitly correlated coupled cluster method to the generation of the three-dimensional potential energy surface (PES) of the Ar-NO(+) cationic complex. A good overall agreement is found with the standard coupled clusters techniques employing correlation consistent atomic basis sets (aug-cc-pVnZ, n= D, T, Q) of Wright et al. This PES is then used in quantum close-coupling scattering and variational calculations to treat the nuclear motions. The bound states energies of the Ar-NO(+) complex obtained by both approaches are in good agreement with the available experimental results. The analysis of the vibrational wavefunctions shows strong anharmonic resonances between the low frequency modes (intermonomer bending and stretching modes) and the wavefunctions exhibit large amplitude motions.     

Singlet�Ctriplet energy gap for trimethylenemethane, oxyallyl diradical, and related species: single- and multireference computational results

We have investigated the through-bond exchange interactions in three non-Kekulé hydrocarbon diradicals on the basis of single- and multireference coupled cluster and related broken-symmetry (BS) methods. The singlet?Ctriplet energy gap (S-T gap) and diradical characters for these species are evaluated. It is found that the spin contamination involved in the BS solutions is non-negligible and the approximate spin-projection method greatly improves the usual BS solutions. As for Mukherjee??s state-specific multireference coupled cluster (MkMRCC) computations, the size-consistent correction with the UHF localized natural orbitals (ULO) is useful to obtain the qualitatively correct 2J values.     

In search of definitive signatures of the elusive NCCO radical

Previous experimental assignments of the fundamental vibrational frequencies of NCCO have been brought into question by subsequent unsuccessful attempts to observe IR signatures of this radical at these frequencies. Here we compute the fundamental vibrational frequencies by applying second-order vibrational perturbation theory to the complete quartic force field computed at the all-electron (AE) coupled cluster singles, doubles, and perturbative triples level [CCSD(T)] with the correlation-consistent, polarized core-valence quadruple-zeta (cc-pCVQZ) basis set, which has tight functions to correctly describe core correlation. The AE-CCSD(T)/cc-pCVQZ geometric parameters are r(e)(N-C)=1.1623 A, r(e)(C-C)=1.4370 A, r(e)(C-O)=1.1758 A, theta(e)(N-C-C)=168.55 degrees , and theta(e)(C-C-O)=132.22 degrees . Our CCSD(T)/cc-pCVQZ values of the characteristic stretching frequencies nu(1) and nu(2) are 2171 and 1898 cm(-1), respectively, in stark contrast to the experimentally derived values of 2093 and 1774 cm(-1). Finally, focal-point extrapolations using correlation-consistent basis sets cc-pVXZ (X=D,T,Q,5,6) and electron correlation treatments as extensive as full coupled cluster singles, doubles, and triples (CCSDT) with perturbative accounting of quadruple excitations [CCSDT(Q)] determine the vibrationless barrier to linearity of NCCO and the dissociation energy (D(0)) of NCCO-->NC+CO to be 8.4 and 26.5 kcal mol(-1), respectively. Using our precisely determined dissociation energy, we recommend a new 0 K enthalpy of formation for NCCO of 50.9+/-0.3 kcal mol(-1).     

Explicitly correlated treatment of H2NSi and H2SiN radicals: electronic structure calculations and rovibrational spectra

Various ab initio methods are used to compute the six dimensional potential energy surfaces (6D-PESs) of the ground states of the H(2)NSi and H(2)SiN radicals. They include standard coupled cluster (RCCSD(T)) techniques and the newly developed explicitly correlated RCCSD(T)-F12 methods. For H(2)NSi, the explicitly correlated techniques are viewed to provide data as accurate as the standard coupled cluster techniques, whereas small differences are noticed for H(2)SiN. These PESs are found to be very flat along the out-of-plane and some in-plane bending coordinates. Then, the analytic representations of these PESs are used to solve the nuclear motions by standard perturbation theory and variational calculations. For both isomers, a set of accurate spectroscopic parameters and the vibrational spectrum up to 4000 cm(-1) are predicted. In particular, the analysis of our results shows the occurrence of anharmonic resonances for H(2)SiN even at low energies.     

Coupled cluster and DFT calculations of 14N nuclear quadrupole coupling constants

The dependence of ¹⁴N quadrupole coupling constants calculated using coupled cluster theory on the level of approximation is examined for a series of small molecules. For HCN, HNC, CH₃CN, and CH₃NC, we use the coupled cluster singles‐and‐doubles with a noniterative perturbative triples correction—CCSD(T)—approach, and we analyze the basis set dependence of the results. For aziridine, diazirine, and cyclopropyl cyanide, we use the CCSD(T) approach, but smaller basis sets, and for the largest studied molecules—quinuclidine and hexamine—we present CCSD results. The differences between computed and experimental values for the best basis sets used are ≈ 5% at the CCSD level and decrease noticeably at the CCSD(T) level. The ‐ N≡C bonds are an exception—in this case the quadrupole coupling constants are very small, hence the differences between theory and experiment become larger (up to 9%). We also consider the performance of density functional theory, comparing the results for different density functionals with the coupled cluster values of the same constants. Most of the functionals provide results systematically improved with respect to the Hartree–Fock values, with ¹⁴N coupling constants in ‐ N≡C bonds being again an exception. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Bond breaking with auxiliary-field quantum Monte Carlo

Bond stretching mimics different levels of electron correlation and provides a challenging test bed for approximate many-body computational methods. Using the recently developed phaseless auxiliary-field quantum Monte Carlo (AF QMC) method, we examine bond stretching in the well-studied molecules BH and N(2) and in the H(50) chain. To control the sign/phase problem, the phaseless AF QMC method constrains the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. With single Slater determinants from unrestricted Hartree-Fock as trial wave function, the phaseless AF QMC method generally gives better overall accuracy and a more uniform behavior than the coupled cluster CCSD(T) method in mapping the potential-energy curve. In both BH and N(2), we also study the use of multiple-determinant trial wave functions from multiconfiguration self-consistent-field calculations. The increase in computational cost versus the gain in statistical and systematic accuracy are examined. With such trial wave functions, excellent results are obtained across the entire region between equilibrium and the dissociation limit.     

A method of incorporating quadruple correction in the scheme of multi-reference singly and doubly excited configuration interaction — a CSF based coupled pair approximation

Summary We develop an approximate size consistent method within a framework of the multi-reference configuration interaction scheme. The Rayleigh-Schrödinger perturbation theory is employed with a specific selection of the unperturbed part of the electronic Hamiltonian. The second order energy is obtained by a set of equations similar to the quasidegenerate variational perturbation theory of Cave and Davidson. The approximate fourth order energy is obtained by solving a set of equations similar to the coupled electron pair approximation (CEPA). The method has been tested for two simple systems, BeH₂ and N₂, and the results are quite encouraging.     

Treating molecules in arbitrary spin states using the parametric two-electron reduced-density-matrix method

Minimizing the electronic energy with respect to a parameterized two-electron reduced density matrix (2-RDM) is known as a parametric variational 2-RDM method. The parametric 2-RDM method with the M 2-RDM parametrization [D. A. Mazziotti, Phys. Rev. Lett. 101, 253002 (2008)] is extended to treat molecules in arbitrary spin states. Like its singlet counterpart, the M parametric 2-RDM method for arbitrary spin states is derived using approximate N-representability conditions, which allow it to capture more correlation energy than coupled cluster with single and double excitations at a lower computational cost. We present energies, optimized bond lengths, potential energy curves, and occupation numbers for a set of molecules in a variety of spin states using the M and K parametric 2-RDM methods as well as several wavefunction methods. We show that the M parametric 2-RDM method can describe bond breaking of open-shell molecules like triplet B(2) and singlet and triplet OH(+) even in the presence of strong correlation. Finally, the computed 2-RDMs are shown to be nearly N-representable at both equilibrium and non-equilibrium geometries.