An orbital-invariant internally contracted multireference coupled cluster approach

We have formulated and implemented an internally contracted multireference coupled cluster (ic-MRCC) approach aimed at solving two of the problems encountered in methods based on the Jeziorski-Monkhorst ansatz: (i) the scaling of the computational and memory costs with respect to the number of references, and (ii) the lack of invariance of the energy with respect to rotations among active orbitals. The ic-MRCC approach is based on a straightforward generalization of the single-reference coupled cluster ansatz in which an exponential operator is applied to a multiconfigurational wave function. The ic-MRCC method truncated to single and double excitations (ic-MRCCSD) yields very accurate potential energy curves in benchmark computations on the Be + H(2) insertion reaction, the dissociation of hydrogen fluoride, and the symmetric double dissociation of water. Approximations of the ic-MRCC theory in which the Baker-Campbell-Hausdorff expansion is truncated up to a given number of commutators are found to converge quickly to the full theory. In our tests, two commutators are sufficient to recover a total energy within 0.5 mE(h) of the full ic-MRCCSD method along the entire potential energy curve. A formal analysis shows that the ic-MRCC method is invariant with respect to rotation among active orbitals, and that the orthogonalization procedure used to produce the set of linearly independent excitation operators plays a crucial role in guaranteeing the invariance properties. The orbital invariance was confirmed in numerical tests. Moreover, approximated versions of the ic-MRCC theory based on a truncated Baker-Campbell-Hausdorff expansion, preserve the orbital invariance properties of the full theory.     

A state-specific partially internally contracted multireference coupled cluster approach

A state-specific partially internally contracted multireference coupled cluster approach is presented for general complete active spaces with arbitrary number of active electrons. The dominant dynamical correlation is included via an exponential parametrization of internally contracted cluster operators ( ?T) which excite electrons from a multideterminantal reference function. The remaining dynamical correlation and relaxation effects are included via a diagonalization of the transformed Hamiltonian ?? =e(- ?T)H?e( ?T) in the multireference configuration interaction singles space in an uncontracted fashion. A new set of residual equations for determining the internally contracted cluster amplitudes is proposed. The second quantized matrix elements of ?? , expressed using the extended normal ordering of Kutzelnigg and Mukherjee, are used as the residual equations without projection onto the excited configurations. These residual equations, referred to as the many-body residuals, do not have any near-singularity and thus, should allow one to solve all the amplitudes without discarding any. There are some relatively minor remaining convergence issues that may arise from an attempt to solve all the amplitudes and an initial analysis is provided in this paper. Applications to the bond-stretching potential energy surfaces for N(2), CO, and the low-lying electronic states of C(2) indicate clear improvements of the results using the many-body residuals over the conventional projected residual equations.     

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.     

Orbitally invariant internally contracted multireference unitary coupled cluster theory and its perturbative approximation: theory and test calculations of second order approximation

A unitary wave operator, exp (G), G(+) = -G, is considered to transform a multiconfigurational reference wave function Φ to the potentially exact, within basis set limit, wave function Ψ = exp (G)Φ. To obtain a useful approximation, the Hausdorff expansion of the similarity transformed effective Hamiltonian, exp (-G)Hexp (G), is truncated at second order and the excitation manifold is limited; an additional separate perturbation approximation can also be made. In the perturbation approximation, which we refer to as multireference unitary second-order perturbation theory (MRUPT2), the Hamiltonian operator in the highest order commutator is approximated by a Mo?ller-Plesset-type one-body zero-order Hamiltonian. If a complete active space self-consistent field wave function is used as reference, then the energy is invariant under orbital rotations within the inactive, active, and virtual orbital subspaces for both the second-order unitary coupled cluster method and its perturbative approximation. Furthermore, the redundancies of the excitation operators are addressed in a novel way, which is potentially more efficient compared to the usual full diagonalization of the metric of the excited configurations. Despite the loss of rigorous size-extensivity possibly due to the use of a variational approach rather than a projective one in the solution of the amplitudes, test calculations show that the size-extensivity errors are very small. Compared to other internally contracted multireference perturbation theories, MRUPT2 only needs reduced density matrices up to three-body even with a non-complete active space reference wave function when two-body excitations within the active orbital subspace are involved in the wave operator, exp (G). Both the coupled cluster and perturbation theory variants are amenable to large, incomplete model spaces. Applications to some widely studied model systems that can be problematic because of geometry dependent quasidegeneracy, H4, P4, and BeH(2), are performed in order to test the new methods on problems where full configuration interaction results are available.     

Development and applications of a unitary group adapted state specific multi-reference coupled cluster theory with internally contracted treatment of inactive double excitations

Any multi-reference coupled cluster (MRCC) development based on the Jeziorski-Monkhorst (JM) multi-exponential ansatz for the wave-operator Ω suffers from spin-contamination problem for non-singlet states. We have very recently proposed a spin-free unitary group adapted (UGA) analogue of the JM ansatz, where the cluster operators are defined in terms of spin-free unitary generators and a normal ordered, rather than ordinary, exponential parametrization of Ω is used. A consequence of the latter choice is the emergence of the "direct?term" of the MRCC equations that terminates at exactly the quartic power of the cluster amplitudes. Our UGA-MRCC ansatz has been utilized to generate both the spin-free state specific (SS) and the state universal MRCC formalisms. It is well-known that the SSMRCC theory requires suitable sufficiency conditions to resolve the redundancy of the cluster amplitudes. In this paper, we propose an alternative variant of the UGA-SSMRCC theory, where the sufficiency conditions are used for all cluster operators containing active orbitals and the single excitations with inactive orbitals, while the inactive double excitations are assumed to be independent of the model functions they act upon. The working equations for the inactive double excitations are thus derived in an internally contracted (IC) manner in the sense that the matrix elements entering the MRCC equations involve excitations from an entire combination of the model functions. We call this theory as UGA-ICID-MRCC, where ICID is the acronym for "Internally Contracted treatment of Inactive Double excitations." Since the number of such excitations are the most numerous, choosing them to be independent of the model functions will lead to very significant reduction in the number of cluster amplitudes for large active spaces, and is worth exploring. Moreover, unlike for the excitations involving active orbitals, where there is inadequate coupling between the model and the virtual functions in the SSMRCC equations generated from sufficiency conditions, our internally contracted treatment of inactive double excitations involves much more complete couplings. Numerical implementation of our formalism amply demonstrates the efficacy of the formalism.     

Parallel internally contracted multireference configuration interaction

A parallel implementation of the internally contracted (IC) multireference configuration (MRCI) module of the MOLPRO quantum chemistry program is described. The global array (GA) toolkit has been used in order to map an existing disk-paging small-memory algorithm onto a massively parallel supercomputer, where disk storage is replaced by the combined memory of all processors. This model has enabled a rather complicated code to be ported to the parallel environment without the need for the wholesale redesign of algorithms and data structures. Examples show that the parallel ICMRCI program can deliver results in a fraction of the time needed for equivalent uncontracted MRCI computations. Further examples demonstrate that ICMRCI computations with up to 10⁷ variational parameters, and equivalent to uncontracted MRCI with 10⁹ configurations, are feasible. The largest calculation demonstrates a parallel efficiency of about 80% on 128 nodes of a Cray T3E-300. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1215–1228, 1998     

Implementation of the multireference Brillouin-Wigner and Mukherjee's coupled cluster methods with non-iterative triple excitations utilizing reference-level parallelism

In this paper we discuss the performance of the non-iterative state-specific multireference coupled cluster (SS-MRCC) methods accounting for the effect of triply excited cluster amplitudes. The corrections to the Brillouin-Wigner and Mukherjee's MRCC models based on the manifold of singly and doubly excited cluster amplitudes (BW-MRCCSD and Mk-MRCCSD, respectively) are tested and compared with exact full configuration interaction results for small systems (H(2)O, N(2), and Be(3)). For the larger systems (naphthyne isomers) the BW-MRCC and Mk-MRCC methods with iterative singles, doubles, and non-iterative triples (BW-MRCCSD(T) and Mk-MRCCSD(T)) are compared against the results obtained with single reference coupled cluster methods. We also report on the parallel performance of the non-iterative implementations based on the use of processor groups.     

A convenient decontraction procedure of internally contracted state-specific multireference algorithms

Internally contracted state-specific multireference (MR) algorithms, either perturbative such as CASPT2 or NEVPT2, or nonperturbative such as contracted MR configuration interaction or MR coupled cluster, are computationally efficient but they may suffer from the internal contraction of the wave function in the reference space. The use of a low dimensional multistate model space only offers limited flexibility and is not always practicable. The present paper suggests a convenient state-specific procedure to decontract the reference part of the wave function from a series of state-specific calculations using slightly perturbed zero-order wave functions. The method provides an orthogonal valence bond reading of the ground state and an effective valence Hamiltonian, the excited roots of which are shown to be relevant. The orthogonal valence bond functions can be considered quasidiabatic states and the effective valence Hamiltonian gives therefore the quasidiabatic energies and the electronic coupling among the quasidiabatic states. The efficiency of the method is illustrated in two case problems where the dynamical correlation plays a crucial role, namely, the LiF neutral/ionic avoided crossing and the F(2) ground state wave function.     

Triple excitations in state-specific multireference coupled cluster theory: application of Mk-MRCCSDT and Mk-MRCCSDT-n methods to model systems

We report the first implementation with correct scaling of the Mukherjee multireference coupled cluster method with singles, doubles, and approximate iterative triples (Mk-MRCCSDT-n, n=1a,1b,2,3) as well as full triples (Mk-MRCCSDT). These methods were applied to the classic H4, P4, BeH(2), and H8 model systems to assess the ability of the Mk-MRCCSDT-n schemes to accurately account for triple excitations. In all model systems the inclusion of triples via the various Mk-MRCCSDT-n approaches greatly reduces the nonparallelism error (NPE) and the mean nonparallelism derivative diagnostics for the potential energy curves, recovering between 59% and 73% of the full triples effect on average. The most complete triples approximation, Mk-MRCCSDT-3, exhibits the best average performance, reducing the mean NPE to below 0.6 mE(h), compared to 1.4 mE(h) for Mk-MRCCSD. Both linear and quadratic truncations of the Mk-MRCC triples coupling terms are viable simplifications producing no significant errors. If the off-diagonal parts of the occupied-occupied and virtual-virtual blocks of the Fock matrices are ignored, the storage of the triples amplitudes is no longer required for the Mk-MRCCSDT-n methods introduced here. This proves to be an effective approximation that gives results almost indistinguishable from those derived from full consideration of the Fock matrices.     

High-order excitations in state-universal and state-specific multireference coupled cluster theories: model systems

For the first time high-order excitations (n>2) have been studied in three multireference couple cluster (MRCC) theories built on the wave operator formalism: (1) the state-universal (SU) method of Jeziorski and Monkhorst (JM) (2) the state-specific Brillouin-Wigner (BW) coupled cluster method, and (3) the state-specific MRCC approach of Mukherjee (Mk). For the H4, P4, BeH(2), and H8 models, multireference coupled cluster wave functions, with complete excitations ranging from doubles to hextuples, have been computed with a new arbitrary-order string-based code. Comparison is then made to corresponding single-reference coupled cluster and full configuration interaction (FCI) results. For the ground states the BW and Mk methods are found, in general, to provide more accurate results than the SU approach at all levels of truncation of the cluster operator. The inclusion of connected triple excitations reduces the nonparallelism error in singles and doubles MRCC energies by a factor of 2-10. In the BeH(2) and H8 models, the inclusion of all quadruple excitations yields absolute energies within 1 kcal mol(-1) of the FCI limit. While the MRCC methods are very effective in multireference regions of the potential energy surfaces, they are outperformed by single-reference CC when one electronic configuration dominates.     

Reduced multireference coupled cluster method with singles and doubles: Perturbative corrections for triples

The reduced multireference coupled-cluster method with singles and doubles (RMR CCSD) that employs multireference configuration interaction wave function as an external source for a small subset of approximate connected triples and quadruples, is perturbatively corrected for the remaining triples along the same lines as in the standard CCSD(T) method. The performance of the resulting RMR CCSD(T) method is tested on four molecular systems, namely, the HF and F(2) molecules, the NO radical, and the F(2) (+) cation, representing distinct types of molecular structure, using up to and including a cc-pVQZ basis set. The results are compared with those obtained with the standard CCSD(T), UCCSD(T), CCSD(2), and CR CCSD(T) methods, wherever applicable or available. An emphasis is made on the quality of the computed potentials in a broad range of internuclear separations and on the computed equilibrium spectroscopic properties, in particular, harmonic frequencies omega(e). It is shown that RMR CCSD(T) outperforms other triply corrected methods and is widely applicable.     

On the internally contracted multireference CI method with full contraction

The configuration interaction (CI ) method where the efficiency of the generators of the unitary group is most fully exploited is the internally contracted multireference CI method. In the most recent version of this method the semi-internal configurations were kept uncontracted, which means that the number of configurations can still be quite large. In the present study the necessary formulas are derived for the case where the semi-internal states are also contracted. The highest density matrix that appears in these formulas is of order 5, and the computational treatment of this large matrix is discussed in detail.     

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.     

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.     

Fock space multireference coupled cluster theory: noniterative inclusion of triples for excitation energies

In this paper we propose and numerically implement a specific scheme for calculating the excitation energies (EEs) within the Fock space multireference coupled cluster framework, which includes the contributions from noniterative triples cluster amplitudes. These contribute to the EEs at the third order. We present results for CH⁺ and N2, and study the effects of these noniterative triples on EEs. Received: 28 July 1997 / Accepted: 8 December 1997     

The coupled cluster approach with a hybrid treatment of connected triple excitations based on the restricted Hartree-Fock reference

A generalization of the coupled cluster (CC) singles, doubles, and a hybrid treatment of connected triples [denoted as CCSD(T)-h] [Shen et al., J. Chem. Phys. 132, 114115 (2010)] to the restricted Hartree-Fock (RHF) reference is presented. In this approach, active (or pseudoactive) RHF orbitals are constructed automatically by performing unitary transformations of canonical RHF orbitals so that they spatially mimic the natural orbitals of the unrestricted Hartree-Fock reference. The present RHF-based CCSD(T)-h approach has been applied to study the potential energy surfaces in several typical bond breaking processes and the singlet-triplet gaps in a diradical (HFH)(-1). For all systems under study, the overall performance of CCSD(T)-h is very close to that of the corresponding CCSD(T) (CC singles, doubles, and triples), and much better than that of CCSD(T) (CC singles, doubles, and perturbative triples).