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

Evangelista, Francesco A.; Allen, Wesley D.; Schaefer, Henry F.
July 2007
Journal of Chemical Physics;7/14/2007, Vol. 127 Issue 2, p024102
Academic Journal
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 (F2) 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 F2 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 a1 vibrational frequencies of ozone, the antisymmetric stretching frequency ω3(b2) 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.


Related Articles

  • Vibrational excitation energies from vibrational coupled cluster response theory. Seidler, Peter; Christiansen, Ove // Journal of Chemical Physics;5/28/2007, Vol. 126 Issue 20, p204101 

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

  • Toward the observation of quartet states of the ozone radical cation: Insights from coupled cluster theory. Speakman, Lucas D.; Turney, Justin M.; Schaefer, Henry F. // Journal of Chemical Physics;6/7/2008, Vol. 128 Issue 21, p214302 

    Since the discovery of ozone depletion, the doublet electronic states of the ozone radical cation have received much attention in experimental and theoretical investigations, while the low-lying quartet states have not. In the present research, viable pathways to the quartet states from the...

  • Improved statistical quark model for the nucleon structure function. Mírez, Carlos // AIP Conference Proceedings;5/22/2010, Vol. 1245 Issue 1, p137 

    An improved statistical quark model, with quark energy levels given by a central linear confining potential, is used to obtain the light sea-quark asymmetry, d/ū, and also for the corresponding difference d–ū, inside the nucleon. In the model, a temperature parameter is adjusted by...

  • Interaction of CO with Kr: Potential energy surface and bound states. Eryin Feng; Zhongquan Wang; Mingyan Gong; Zhifeng Cui // Journal of Chemical Physics;11/7/2007, Vol. 127 Issue 17, p174301 

    The first ab initio potential energy surface of the Kr–CO complex is developed using single and double excitation coupled-cluster theory with noniterative treatment of triple excitations. Mixed basis sets, aug-cc-pVQZ for the C and O atoms and aug-cc-pVQZ-PP for the Kr atom, with an...

  • A full-configuration interaction “nuclear orbital” method to study doped 3HeN clusters (N<=4). de Lara-Castells, M. P.; Delgado-Barrio, G.; Villarreal, P.; Mitrushchenkov, A. O. // Journal of Chemical Physics;12/14/2006, Vol. 125 Issue 22, p221101 

    An efficient full configuration interaction (FCI) treatment, based on the Jacobi-Davidson algorithm, is developed in order to study small doped 3HeN clusters. The state of each He atom in a given cluster is described by a set of wave-functions which by extention of the quantum-chemistry notation...

  • New calculations of α-decay half-lives by the Viola-Seaborg formula. Tiekuang Dong; Zhongzhou Ren // European Physical Journal A -- Hadrons & Nuclei;Oct2005, Vol. 26 Issue 1, p69 

    The Viola-Seaborg formula is a well-known formula for α-decay half-lives of heavy nuclei. In this work we obtain new parameters of this formula through a least-square fit to even-even nuclei between Z = 84 and Z = 110 with N greater than 126. On average, the formula can reproduce the...

  • Ab initio study of the energy difference between the cyclic and linear forms of the C6 molecule. Pless, V.; Suter, H. U.; Engels, B. // Journal of Chemical Physics;9/1/1994, Vol. 101 Issue 5, p4042 

    The energy difference between the three lowest-lying isomers of C6, the linear 3Σ-g state, and the two ring forms, the benzene structure (1A1g) possessing D6h symmetry and a distorted cyclic form (1A1’, D3h symmetry) have been calculated using various ab initio methods. Variational...

  • Quantum dynamics of vibrationally activated OH–CO reactant complexes. He, Yong; Goldfield, Evelyn M.; Gray, Stephen K. // Journal of Chemical Physics;7/8/2004, Vol. 121 Issue 2, p823 

    A six-dimensional wave packet study of the unimolecular decay of vibrationally activated OH–CO reactant channel complexes is presented. The ab initio based Lakin–Troya–Schatz–Harding potential energy functions for the A′ and A″ states are employed. Good...

  • Reaction of chlorine atom with trichlorosilane from 296 to 473 K. DeSain, John D.; Valachovic, Linda; Jusinski, Leonard E.; Taatjes, Craig A. // Journal of Chemical Physics;12/14/2006, Vol. 125 Issue 22, p224308 

    The reaction of trichlorosilane (HSiCl3) with atomic chlorine (Cl) has been investigated by using infrared kinetic spectroscopy of the HCl product. The overall second order rate constant for the reaction has been determined as a function of temperature by using pseudo-first-order kinetic...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics