Benchmark studies of variational, unitary and extended coupled cluster methods

Comparative benchmark calculations are presented for coupled cluster theory in its standard formulation, as well as variational, extended, and unitary coupled cluster methods. The systems studied include HF, N(2), and CN, and with cluster operators that for the first time include up to quadruple excitations. In cases where static correlation effects are weak, the differences between the predictions of molecular properties from each theory are negligible. When, however, static correlation is strong, it is demonstrated that variational coupled cluster theory can be significantly more robust than the traditional ansatz and offers a starting point on which to base single-determinant reference methods that can be used beyond the normal domain of applicability. These conclusions hold at all levels of truncation of the cluster operator, with the variational approach showing significantly smaller errors.     

Correction to constrained coupled cluster doubles models based on the second coupled cluster central moment

We develop a correction for the coupled cluster version of the perfect pairing (PP) model. The correction is based on finding modified values of the PP amplitudes such that the second coupled cluster central moment defined in the space of all valence single and double substitutions vanishes and, subject to this constraint, minimizing the deviation between the modified and unmodified PP amplitudes with respect to a chosen metric. We discuss how this correction can be generalized to other constrained doubles models, such as local correlation and active-space models. While the correction is not strictly size consistent and retains some of the deficiencies of the PP model, numerical results indicate that much of the missing active-space coupled cluster singles and doubles correlation energy is recovered.     

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.     

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.     

The relationship between double excitation amplitudes and Z vector components in some post-Hartree-Fock correlation methods

The relationship between Z vector components and excitation amplitudes is analyzed for several post-Hartree-Fock correlation methods limited to double excitation amplitudes. An analytical formula approximating the Z vector for the coupled cluster doubles method is presented and shown to be quite accurate. This approximation is also used to determine the prefactor of the norm of doubly excited states in averaged coupled pair functional-type energy functionals self-consistently leading to better agreement with coupled cluster results.     

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.     

Combustion chemistry: important features of the C3H5 potential energy surface, including allyl radical, propargyl + H2, allene + H, and eight transition states

The C(3)H(5) potential energy surface (PES) encompasses molecules of great significance to hydrocarbon combustion, including the resonantly stabilized free radicals propargyl (plus H(2)) and allyl. In this work, we investigate the interconversions that take place on this PES using high level coupled cluster methodology. Accurate geometries are obtained using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)] combined with Dunning's correlation consistent quadruple-ζ basis set cc-pVQZ. The energies for these stationary points are then refined by a systematic series of computations, within the focal point scheme, using the cc-pVXZ (X = D, T, Q, 5, 6) basis sets and correlation treatments as extensive as coupled cluster with full single, double, and triple excitation and perturbative quadruple excitations [CCSDT(Q)]. Our benchmarks provide a zero-point vibrational energy (ZPVE) corrected barrier of 10.0 kcal mol(-1) for conversion of allene + H to propargyl + H(2). We also find that the barrier for H addition to a terminal carbon atom in allene leading to propenyl is 1.8 kcal mol(-1) lower than that for the addition to a central atom to form the allyl radical.     

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.      

Benchmarking two-photon absorption with CC3 quadratic response theory, and comparison with density-functional response theory

We present a detailed study of the effects of electron correlation on two-photon absorption calculated by coupled cluster quadratic response theory. The hierarchy of coupled cluster models CCS, CC2, CCSD, and CC3 has been used to investigate the effects of electron correlation on the two-photon absorption cross sections of formaldehyde (CH2O), diacetylene (C4H2), and water (H2O). In particular, the effects of triple excitations on two-photon transition cross sections are determined for the first time. In addition, we present a detailed comparison of the coupled cluster results with those obtained from Hartree-Fock and density-functional response theories. We have investigated the local-density approximation, the pure Becke-Lee-Yang-Parr (BLYP) functional, the hybrid Becke-3-parameter-Lee-Yang-Parr (B3LYP), and the Coulomb-attenuated B3LYP (CAM-B3LYP) functionals. Our results show that the CAM-B3LYP functional, when used in conjuction with a one-particle basis-set containing diffuse functions, has much promise; however, care must still be exercised for diffuse Rydberg-type states.     

Toward subchemical accuracy in computational thermochemistry: focal point analysis of the heat of formation of NCO and [H,N,C,O] isomers

In continuing pursuit of thermochemical accuracy to the level of 0.1 kcal mol(-1), the heats of formation of NCO, HNCO, HOCN, HCNO, and HONC have been rigorously determined using state-of-the-art ab initio electronic structure theory, including conventional coupled cluster methods [coupled cluster singles and doubles (CCSD), CCSD with perturbative triples (CCSD(T)), and full coupled cluster through triple excitations (CCSDT)] with large basis sets, conjoined in cases with explicitly correlated MP2-R12/A computations. Limits of valence and all-electron correlation energies were extrapolated via focal point analysis using correlation consistent basis sets of the form cc-pVXZ (X=2-6) and cc-pCVXZ (X=2-5), respectively. In order to reach subchemical accuracy targets, core correlation, spin-orbit coupling, special relativity, the diagonal Born-Oppenheimer correction, and anharmonicity in zero-point vibrational energies were accounted for. Various coupled cluster schemes for partially including connected quadruple excitations were also explored, although none of these approaches gave reliable improvements over CCSDT theory. Based on numerous, independent thermochemical paths, each designed to balance residual ab initio errors, our final proposals are DeltaH(f,0) ( composite function )(NCO)=+30.5, DeltaH(f,0) ( composite function )(HNCO)=-27.6, DeltaH(f,0) ( composite function )(HOCN)=-3.1, DeltaH(f,0) ( composite function )(HCNO)=+40.9, and DeltaH(f,0) ( composite function )(HONC)=+56.3 kcal mol(-1). The internal consistency and convergence behavior of the data suggests accuracies of +/-0.2 kcal mol(-1) in these predictions, except perhaps in the HCNO case. However, the possibility of somewhat larger systematic errors cannot be excluded, and the need for CCSDTQ [full coupled cluster through quadruple excitations] computations to eliminate remaining uncertainties is apparent.     

Development of a relaxation-inducing cluster expansion formalism for treating strong relaxation and correlation effects

We present in this paper a comprehensive account of an explicitly spin-free coupled cluster theory for treating energy differences of open-shell states relative to a closed-shell ground state, where the open-shell states of interest are dominated by a few simple configuration state functions. We develop a valence-universal coupled cluster formalism to achieve this via a novel cluster expansion ansatz for the valence part of the wave operator, where the orbital relaxation and the correlation relaxation accompanying ionization/excitation from the ground state are taken care of to all orders in compact, efficient, and explicitly spin-free manner. The essential difference of our proposed ansatz from the ordinary and the normal-ordered cluster ansatz in vogue is that (a) we allow the valence cluster operators to be connected among themselves with spectator valence lines only and (b) we use suitable combinatoric factors accompanying powers of cluster operators thus connected, which are equal to the number of ways the operators can be joined, leading to the same excitation (the automorphic factor). We emphasize that such an ansatz does not generate terms (diagrams) with chains of cluster operators joined among themselves via spectator lines only. Barring only a few, almost all the terms in the working equations determining the cluster amplitudes involve contraction of the Hamiltonian with the cluster operators via at least one nonspectator line, leading to what we call a "strongly connected" series. The structure of the working equation is remarkably similar to the single-reference closed-shell equation, with a few additional terms. The presence of contractions among cluster operators via spectator lines introduces the additional physical effects of orbital and correlation relaxation using low-body cluster operators. As an illustrative application of the new multireference coupled cluster (CC) theory, we consider in this paper computation of ionization potentials (IPs) of one-valence problem with only one active orbital. The numerical applications are made for both the core- and the inner- and outer-valence IPs for several molecular systems. The numerical values demonstrate the superiority of the relaxation-inducing CC theory, as compared to the normal-ordered ansatz.     

Extracting dominant pair correlations from many-body wave functions

The singular value decomposition of the n-particle excitation operator as determined by coupled cluster or perturbation theory is used to extract the dominant and interesting electron-electron correlations from complex molecular wave functions. As an example of the very general formalism, the decomposition of the T(2) operator obtained from coupled cluster doubles calculations is used to analyze the strength and character of pair correlations in a variety of molecules with interesting electronic structure. The magnitude of the largest singular value(s) determines the strength of the correlation(s), and the corresponding right- and left-hand singular vectors characterize the physical and spatial nature of the correlations. The primary advantage of this tool over natural orbital analysis is that it provides direct associations between the occupied and virtual geminals involved in the correlations.     

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.     

A simple correction to final state energies of doublet radicals described by equation-of-motion coupled cluster theory in the singles and doubles approximation

A computationally inexpensive energy correction is suggested for radicals described by the equation-of-motion coupled cluster method for ionized states in the singles and doubles approximation (EOMIP-CCSD). The approach is primarily intended for doublet states that are qualitatively described by Koopmans' approximation. Following a strategy similar to those used in multireference coupled cluster theory, the proposed correction accounts for all correlation effects through third order in perturbation theory and also includes selected contributions to higher-order energies. As an initial test of the numerical performance of the method, total energies and energy splittings are calculated for some small prototype radicals.     

Reciprocal adjustment of approximate coupled cluster and configuration interaction approaches

The linear version of the externally corrected coupled cluster method with singles and doubles (ecLCCSD), the recently proposed coupled cluster corrections to the multireference configuration interaction (ccMRCI) energies, and the so‐called self‐consistent, size‐consistent [(SC)²] approaches, which are designed to correct for the dynamic correlation effects and the size inconsistency of the MRCI energies, are analyzed and compared using several illustrative examples, including the dissociation of a triple‐zeta (TZ) model of the N₂ molecule. It is emphasized that the exponential cluster ansatz for the wave function is the basis of all these approaches, and appropriate cluster analysis of the MRCI wave function is the key step for both ecLCCSD and ccMRCI. The contributions from the orthogonal complement of the MRCI space, which can be generated by relying on such a cluster analysis, are responsible for a substantial part of the missing correlation energy. The ecLCCSD approach seems to represent a particularly attractive alternative to other highly accurate methods for the calculation of the ground‐state energy in the presence of quasidegeneracy, both due to its efficiency and affordability. It may in fact be regarded as a simple alternative to the iterative reduced multireference (RMR) CCSD method. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 693–703, 2000     

Electronic structure of spherical quantum dots using coupled cluster method

2, 6, 12, and 20 electron quantum dots have been studied using coupled cluster at singles and doubles level and extensive multireference coupled cluster (MRCC) method. A Fock-space version of MRCC (FSMRCC) containing single hole-particle excited determinants has been used to calculate low-lying excited states of the above system. The ionization potential and electron affinity are also calculated. The effect of correlation energy on excitation energy and charge density is shown by calculating them at the high density region (low value of density parameter rs) and at the low density region (high value of density parameter rs).     

Electron-correlated calculations of the nuclear-shielding polarizability and magnetizability polarizability of H2, N2, HF, and CO

The effects of a uniform static electric field on the NMR shielding and magnetizability of H₂, N₂, HF, and CO are reported. The appropriate defining quantities are the shielding and magnetizability polarizabilities and these are calculated in a mixed numeric-analytic scheme with allowance for electron correlation via third-order Møller-Plesset (MP3) theory or linearized coupled cluster double excitation (L-CCD) theory. This is the first time that these theories have been applied to these properties. Our primary conclusion for the shielding polarizabilities is that, with the exception of H₂, the use of coupled cluster theory is necessary to establish the definitive values for these properties because of their extreme sensitivity to electron correlation. For the magnetizability polarizabilities we find the MP3 values to be in better agreement with the L-CCD results than the MP2 values are.