A relativistic 4-component general-order multi-reference coupled cluster method: initial implementation and application to HBr

We present the initial implementation of a determinant-based general-order coupled cluster method which fully accounts for relativistic effects within the four-component framework. The method opens the way for the treatment of multi-reference problems through a state-selective expansion of the model space. The evaluation of the coupled cluster vector function is carried out via relativistic configuration interaction expansions. The implementation is based on a large-scale configuration interaction technique, which may efficiently treat long determinant expansions of more than 10⁸ terms. We demonstrate the capabilities of the new method in calculations of complete potential energy curves of the HBr molecule. The inclusion of spin–orbit interaction and higher excitations than coupled cluster double excitations, either by multi-reference model spaces or the inclusion of full iterative triple excitations, lead to highly accurate results for spectral constants of HBr. An erratum to this article can be found at     

Two- and four-component relativistic generalized-active-space coupled cluster method: implementation and application to BiH

A string-based coupled-cluster method of general excitation rank and with optimal scaling which accounts for special relativity within the four-component framework is presented. The method opens the way for the treatment of multi-reference problems through an active-space inspired single-reference based state-selective expansion of the model space. The evaluation of the coupled-cluster vector function is implemented by considering contractions of elementary second-quantized operators without setting up the amplitude equations explicitly. The capabilities of the new method are demonstrated in application to the electronic ground state of the bismuth monohydride molecule. In these calculations simulated multi-reference expansions with both doubles and triples excitations into the external space as well as the regular coupled-cluster hierarchy up to full quadruples excitations are compared. The importance of atomic outer core-correlation for obtaining accurate results is shown. Comparison to the non-relativistic framework is performed throughout to illustrate the additional work of the transition to the four-component relativistic framework both in implementation and application. Furthermore, an evaluation of the highest order scaling for general-order expansions is presented.     

The structure of the Si9H12 cluster: A coupled cluster and multi-reference perturbation theory study

Full geometry optimizations using both singles and doubles coupled cluster theory with perturbative triple excitations, CCSD(T), and second order multi-reference perturbation theory, MRMP2, have been employed to predict the structure of Si9H12, a cluster commonly used in calculations to represent the Si(100) surface. Both levels of theory predict the structure of this cluster to be symmetric (not buckled), and no evidence for a buckled (asymmetric) structure is found at either level of theory.     

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     

Fock space multi-reference coupled cluster study of transition moment and oscillator strength

Summary A method of calculating transition moment and oscillator strength within the framework of the Fock space multi-reference coupled cluster method is described. Diagrammatic technique is used to obtain coupled cluster equations. The general form of equations for the transition moment betweenN-electron ground and excited states is obtained. MBPT analysis of the final equations is done. The excitation energies, dipole transition moments and oscillator strengths for theCH ⁺ molecule are calculated.     

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.     

Vibrational coupled cluster response theory: a general implementation

The calculation of vibrational contributions to molecular properties using vibrational coupled cluster (VCC) response theory is discussed. General expressions are given for expectation values, linear response functions, and transition moments. It is shown how these expressions can be evaluated for arbitrary levels of excitation in the wave function parameterization as well as for arbitrary coupling levels in the potential and property surfaces. The convergence of the method is assessed by benchmark calculations on formaldehyde. Furthermore, excitation energies and infrared intensities are calculated for the fundamental vibrations of furan using VCC limited to up to two-mode and up to three-mode excitations, VCC[2] and VCC[3], as well as VCC with full two-mode and approximate three-mode couplings, VCC[2pt3].     

A CSF-based multi-reference coupled pair approximation III. An application to F2, As2 and As+2

Using the newly developed multi-reference coupled pair approximation program code, the adiabatic potential curves of the ground states of F ₂, As ₂ and As ₂ ⁺ were calculated. Computed spectroscopic constants of these molecules were found to be in good agreement with experimental values. The resulting binding energy of As ₂ (3.86 eV) was compared with the experimental value of 3.99 eV [15] and the best multi-reference configuration interaction value (3.58 eV) reported previously by the present authors. The calculated first adiabatic ionization potential of As ₂ (9.67 eV) was found to be in good agreement with the experimental result. Received: 5 July 1997 / Accepted: 27 August 1997     

A sequential transformation approach to the internally contracted multireference coupled cluster method

The internally contracted multireference coupled cluster (ic-MRCC) approach is formulated using a new wave function ansatz based on a sequential transformation of the reference function (sqic-MRCC). This alternative wave function simplifies the formulation of computationally viable methods while preserving the accuracy of the ic-MRCC approach. The structure of the sqic-MRCC wave function allows folding the effect of the single excitations into a similarity-transformed Hamiltonian whose particle rank is equal to the one of the Hamiltonian. Consequently, we formulate an approximation to the sqic-MRCC method with singles and doubles (included respectively up to fourfold and twofold commutators, sqic-MRCCSD[2]) that contains all terms present in the corresponding single-reference coupled cluster scheme. Computations of the potential energy curves for the dissociation of BeH(2) show that the untruncated sqic-MRCCSD scheme yields results that are almost indistinguishable from the ordinary ic-MRCCSD method. The energy obtained from the computationally less expensive sqic-MRCCSD[2] approximation is found to deviate from the full ic-MRCCSD method by less than 0.2 mE(h) for BeH(2), while, in the case of water, the harmonic vibrational frequencies of ozone, the singlet-triplet splitting of p-benzyne, and the dissociation curve of N(2), sqic-MRCCSD[2] faithfully reproduces the results obtained via the ic-MRCCSD scheme truncated to two commutators. A formal proof is given of the equivalence of the ic-MRCC and sqic-MRCC methods with the internally contracted and full configuration interaction approaches.     

Analytic gradient for the multireference Brillouin-Wigner coupled cluster method and for the state-universal multireference coupled cluster method

We present the analytic gradient theory and its pilot implementation for the multireference Brillouin-Wigner coupled cluster (BWCC) method and for the state-universal multireference coupled cluster method. The analytic gradient has been derived for three cases: (i) BWCC method without a size-extensivity correction, (ii) BWCC method with the iterative size-extensivity correction, and (iii) for the rigorously size-extensive state-universal method. The pilot implementation is based on full-configuration interaction expansions and is presently limited to single and double excitation levels; however, the resulting equations are general. For BWCC methods, they also do not contain terms explicitly mixing amplitudes of different reference configurations and can thus be implemented in an efficient way. The analytic gradients have been verified with respect to numerically computed ones on the example of CH2 molecule, and geometry optimizations of CH2 and SiH2 have been carried out.     

A new method to generate spin-orbit coupled potential energy surfaces: effective relativistic coupling by asymptotic representation

A new method has been developed to generate fully coupled potential energy surfaces including derivative and spin-orbit coupling. The method is based on an asymptotic (atomic) representation of the molecular fine structure states and a corresponding diabatization. The effective relativistic coupling is described by a constant spin-orbit coupling matrix and the geometry dependence of the coupling is accounted for by the diabatization. This approach is very efficient, particularly for certain systems containing a very heavy atom, and yields consistent results throughout nuclear configuration space. A first application to a diatomic system is presented as proof of principle and is compared to accurate ab initio calculations. However, the method is widely applicable to general polyatomic systems in full dimensionality, containing several relativistic atoms and treating higher order relativistic couplings as well.     

Multireference Brillouin-Wigner coupled cluster method with singles, doubles, and triples: efficient implementation and comparison with approximate approaches

We have developed an efficient implementation of the multireference Brillouin-Wigner coupled cluster method with full iterative treatment of connected singles, doubles, and triples. Its computational costs are too high for applications to larger molecules; however, it can be used as a useful tool for benchmarking approximate methods. Performance of the method has been tested on the ground and low-lying excited states of the oxygen molecule and on the singlet-triplet gap in methylene; the results are in good agreement with experimental data.     

An exponential multi-reference wavefunction ansatz: connectivity analysis and application to N2

A connectivity analysis for the exponential multi-reference wavefunction ansatz (MRexpT) (J Chem Phys 123:84102, 2005) is carried out. Assuming a complete model space and separating interactions carrying active labels the cluster operator carrying no active labels is fully connected. The valence (active) part of the MRexpT cluster operator, however remains disconnected. Consequently, the MRexpT correlation energy scales linearly with the number of non-active electrons as single reference coupled cluster does while MRexpT additionally can treat multi reference cases. Therefore, MRexpT should be well suited to be applied to a large number of molecular applications. Its applicability to periodic systems with multi-reference unit cells however seems to be limited. An application to the triple bond breaking of the N₂ molecule is presented.     

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.     

The orbital-specific-virtual local coupled cluster singles and doubles method

We extend the orbital-specific-virtual tensor factorization, introduced for local M?ller-Plesset perturbation theory in Ref. [J. Yang, Y. Kurashige, F. R. Manby and G. K. L. Chan, J. Chem. Phys. 134, 044123 (2011)], to local coupled cluster singles and doubles theory (OSV-LCCSD). The method is implemented by modifying an efficient projected-atomic-orbital local coupled cluster program (PAO-LCCSD) described recently, [H.-J. Werner and M. Schu?tz, J. Chem. Phys. 135, 144116 (2011)]. By comparison of both methods we find that the compact representation of the amplitudes in the OSV approach affords various advantages, including smaller computational time requirements (for comparable accuracy), as well as a more systematic control of the error through a single energy threshold. Overall, the OSV-LCCSD approach together with an MP2 correction yields small domain errors in practical calculations. The applicability of the OSV-LCCSD is demonstrated for molecules with up to 73 atoms and realistic basis sets (up to 2334 basis functions).     

Aspects of size extensivity in unitary group adapted multi-reference coupled cluster theories: the role of cumulant decomposition of spin-free reduced density matrices

We present in this paper a comprehensive study of the various aspects of size extensivity of a set of unitary group adapted multi-reference coupled cluster (UGA-MRCC) theories recently developed by us. All these theories utilize a Jeziorski–Monkhorst (JM) inspired spin-free cluster Ansatz of the form \(|\varPsi \rangle = \sum\nolimits_\mu \varOmega _\mu |\phi _\mu \rangle c_\mu\) with \(\varOmega _\mu =\{\exp ({T_\mu })\}\) , where \(T_\mu\) is expressed in terms of spin-free generators of the unitary group \(U(n)\) for n-orbitals with the associated cluster amplitudes. \(\{...\}\) indicates normal ordering with respect to the common closed shell \(core\) part, \(|0\rangle\) , of the model functions, \(\{\phi _\mu \}\) which is taken as the vacuum. We argue and emphasize in the paper that maintaining size extensivity of the associated theories is consequent upon (a) connectivity of the composites, \(G_\mu\) , containing the Hamiltonian \(H\) and the various powers of \(T\) connected to it, (b) proving the connectivity of the MRCC equations which involve not only \(G_\mu\) s but also the associated connected components of the spin-free reduced density matrices (RDMs) obtained via their cumulant decomposition and (c) showing the extensivity of the cluster amplitudes for non-interacting groups of orbitals and eventually of the size-consistency of the theories in the fragmentation limits. While we will discuss the aspect (a) above rather briefly, since this was amply covered in our earlier papers, the aspect (b) and (c), not covered in detail hitherto, will be covered extensively in this paper. The UGA-MRCC theories dealt with in this paper are the spin-free analogs of the state-specific and state-universal MRCC developed and applied by us recently.We will explain the unfolding of the proof of extensivity by analyzing the algebraic structure of the working equations, decomposed into two factors, one containing the composite \(G_\mu\) that is connected with the products of cumulants arising out of the cumulant decomposition of the RDMs and the second term containing some RDMs which is disconnected from the first and can be factored out and removed. This factorization ultimately leads to a set of connected MRCC equations. Establishing the extensivity and size-consistency of the theories requires careful separation of truly extensive cumulants from the ones which are a measure of spin correlation and are thus connected but not extensive. We have discussed in detail, using diagrams, the factorization procedure and have used suitable example diagrams to amplify the meanings of the various algebraic quantities of any diagram. We conclude the paper by summarizing our findings and commenting on further developments in the future.     

Charge-transfer separability and size-extensivity in the equation-of-motion coupled cluster method: EOM-CCx

We study the charge-transfer separability (CTS) property of the Fock space (FS) and equation-of-motion (EOM) coupled cluster (CC) methods by analysing the charge-transfer (CT) excitation energy versus the donor-acceptor (D-A) distance. All FS-CC approaches fulfill the CT separability condition which is not the case for the standard EOM-CC approaches. This defect of the EOM-CC scheme can be fixed by slight modification of the H matrix's diagrammatic structure, namely by adding some "dressing" composed of disconnected terms. The latter guarantee CTS of the respective EOM-CC scheme and marginally improve local excitations. The newly proposed variant of the EOM-CCSD approach is termed EOM-CCSDx (size-extensive EOM-CCSD).     

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.     

Metalloenzyme mimic: diironhexacarbonyl cluster coupled to redox-active 4-mercapto-1,8-naphthalic anhydride ligands

Transition Metal Chemistry - The non-innocent redox-active ligand, 4-mercapto-1,8-naphthalic anhydride (HS-NAH), has been used in the design and synthesis of a diironhexacarbonyl complex,...