Matrix‐covariant representation of high‐order configuration interaction and coupled cluster theories

We present the closed form of the reduced density matrices (RDMs) of arbitrary order for configuration interaction (CI) wave functions at any excitation level, up to the full CI. A special operator technique due to Bogoliubov is applied and extended. It focuses on constructions of matrix‐covariant expressions independent of the basis set used. The corresponding variational CI equations are given in an explicit form containing the matrices related to conventional excitation operators. A subsequent transformation of the latter to an irreducible form makes it possible to generate the matrix‐covariant representation for coupled cluster (CC) models. Here this transformation is performed for a simplified high‐order CC scheme somewhat reminiscent of the quadratic CI model. A generalized spin‐flip approximation closely related to high‐order CI and CC models is presented, stressing on a possible inclusion of nondynamical and dynamical correlation effects for multiple bond breaking. A derivation of the full CI and simple CC models for systems involving effective three‐electron interactions is also given, thereby demonstrating the capability of the proposed method to deal with complicated many‐body problem. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

The spin-flip extended single excitation configuration interaction method

An extension of the spin-flip single excitation configuration interaction (SF-CIS) method is introduced. The extension, abbreviated as SF-XCIS, includes all configurations in which no more than one virtual level of the high spin triplet reference becomes occupied and no more than one doubly occupied level becomes vacant. The number of such configurations is quadratic with molecule size, and the method is implemented in a direct algorithm whose cost scales in the same way with molecule size as CIS itself, thus permitting applications to large systems. Starting from a spin restricted triplet determinant, SF-XCIS yields spin-pure singlet, triplet, and quintet states, and treats both half-occupied reference orbitals in a fully balanced way to allow application to strongly correlated problems. Tests on bond dissociation in the HF molecule, the torsional potential of ethylene, and excited states of polyenes show encouraging improvements using SF-XCIS compared to SF-CIS and a previously suggested extension, the spin-complete CIS model.     

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     

Explicitly correlated coupled cluster F12 theory with single and double excitations

Full explicitly correlated F12 coupled cluster theory with single and double excitations and with Slater-type geminal as a correlation factor is introduced and implemented within the standard approximation. The variant "C" that does not require integrals over the commutator between the kinetic operator and the correlation factor has been used. All the necessary integrals are analytically calculated. With variant C also, first results are reported for the correlation factor being the interelectronic distance coordinate, i.e., for original R12 method. Calculations have been performed for a set of eight molecules including CH(2)((1)A(1)), CH(4), NH(3), H(2)O, HF, CO, N(2), and F(2), as well as for the constituting atoms. Atomization energies are reported too.     

The convergence of complete active space self-consistent-field configuration interaction including all single and double excitation energies to the complete basis set limit

Examination of the convergence of full valence complete active space self-consistent-field configuration interaction including all single and double excitation (CASSCF-CISD) energies with expansion of the one-electron basis set reveals a pattern very similar to the convergence of single determinant energies. Calculations on the lowest four singlet states and the lowest four triplet states of N(2) with the sequence of n-tuple-zeta augmented polarized (nZaP) basis sets (n=2, 3, 4, 5, and 6) are used to establish the complete basis set limits. Full configuration-interaction (CI) and core electron contributions must be included for very accurate potential energy surfaces. However, a simple extrapolation scheme that has no adjustable parameters and requires nothing more demanding than CAS(10e(-),8orb)-CISD/3ZaP calculations gives the R(e), omega(e), omega(e)X(e), T(e), and D(e) for these eight states with rms errors of 0.0006 Angstrom, 4.43 cm(-1), 0.35 cm(-1), 0.063 eV, and 0.018 eV, respectively.     

Quasidegenerate scaled opposite spin second order perturbation corrections to single excitation configuration interaction

Scaled opposite spin (SOS) second order perturbative corrections to single excitation configuration interaction (CIS) are extended to correctly treat quasidegeneracies between excited states. Two viable methods, termed as SOS-CIS(D(0)) and SOS-CIS(D(1)), are defined, implemented, and tested. Each involves one empirical parameter (plus a second for the SOS-MP2 ground state), has computational cost that scales with the fourth power of molecule size, and has storage requirements that are cubic, with only quantities of the rank of single excitations produced and stored during iterations. Tests on a set of low-lying adiabatic valence excitation energies and vertical Rydberg excitations of organic and inorganic molecules show that the empirical parameter can be acceptably transferred from the corresponding nondegenerate perturbation theories without any further fitting. Further tests on higher excited states show that the new methods correctly perform for surface crossings for which nondegenerate approaches fail. Numerical results show that SOS-CIS(D(0)) appears to treat Rydberg excitations in a more balanced way than SOS-CIS(D(1)) and is, therefore, likely to be the preferred approach. It should be useful for exploring excited state geometries, transition structures, and conical intersections for states of medium to large organic molecules that are dominated by single excitations.     

Coupled cluster studies. III. Comparison of the numerical behaviour of coupled cluster doubles with configuration interaction and perturbation theory. Basis set and geometry optimizations

In the beryllium atom most of the valence electron correlation energy can be obtained with only three equivalent excitations. It is shown that this is possible also with only one primitive set of p-type gaussians, provided their exponent is carefully optimized. Further it is demonstrated that this is due to the near degeneracy of the 1s²2s² and 1s²2p² states of Be, because it does not occur in He and H₂. Møller—Plesset (MP) perturbation theory gives less correlation energy than coupled cluster doubles (CCD) in these cases but the same optimum exponent for the set of p functions. For the molecules CO and CO₂ where MP2 is known to overestimate the correlation energy in comparison to CCD the convergence properties of MP are studied up to fourth order. It is shown that the results of MP4(DQ) and CCD are very similar also in these cases. For comparison results for ethylene are also discussed. Due to size consistency CCD is superior over configuration interaction with double excitations (CID). Equilibrium properties computed with CCD and MP4(DQ) are shown to be as reliable as those given by CID, although the variational theorem does not hold for CCD and MP4(DQ). The effect of single excitations turned out to be negligible in these cases. In basis set optimizations all correlation methods studied worked out similar.     

Quasiclassical trajectory study of the Cl+CH4 reaction dynamics on a quadratic configuration interaction with single and double excitation interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the Cl+CH(4) reactive system has been constructed using the interpolation method of Collins and co-workers [J. Chem. Phys. 102, 5647 (1995); 108, 8302 (1998); 111, 816 (1999); Theor. Chem. Acc. 108, 313 (2002)]. The ab initio calculations have been performed using quadratic configuration interaction with single and double excitation theory to build the PES. A simple scaling all correlation technique has been used to obtain a PES which yields a barrier height and reaction energy in good agreement with high level ab initio calculations and experimental measurements. Using these interpolated PESs, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations, and internal energy distributions has been carried out for the Cl+CH(4) and Cl+CD(4) reactions, and the theoretical results have been compared with the available experimental data. It has been shown that the calculated total reaction cross sections versus collision energy for the Cl+CH(4) and Cl+CD(4) reactions is very sensitive to the barrier height. Besides, due to the zero-point energy (ZPE) leakage of the CH(4) molecule to the reaction coordinate in the quasiclassical trajectory (QCT) calculations, the reaction threshold falls below the barrier height of the PES. The ZPE leakage leads to CH(3) and HCl coproducts with internal energy below its corresponding ZPEs. We have shown that a Gaussian binning (GB) analysis of the trajectories yields excitation functions in somehow better agreement with the experimental determinations. The HCl(v'=0) and DCl(v'=0) rotational distributions are as well very sensitive to the ZPE problem. The GB correction narrows and shifts the rotational distributions to lower values of the rotational quantum numbers. However, the present QCT rotational distributions are still hotter than the experimental distributions. In both reactions the angular distributions shift from backward peaked to sideways peaked as collision energy increases, as seen in the experiments and other theoretical calculations.     

Review of biorthogonal coupled cluster representations for electronic excitation

Single-reference coupled-cluster (CC) methods for electronic excitation are based on a biorthogonal representation (bCC) of the (shifted) Hamiltonian in terms of excited CC states, also referred to as correlated excited (CE) states, and an associated set of states biorthogonal to the CE states, the latter being essentially configuration interaction (CI) configurations. The bCC representation generates a non-hermitian secular matrix, the eigenvalues representing excitation energies, while the corresponding spectral intensities are to be derived from both the left and right eigenvectors. Using the perspective of the bCC representation, a systematic and comprehensive analysis of the excited-state CC methods is given, extending and generalizing previous such studies. Here, the essential topics are the truncation error characteristics and the separability properties, the latter being crucial for designing size-consistent approximation schemes. Based on the general order relations for the bCC secular matrix and the (left and right) eigenvector matrices, formulas for the perturbation-theoretical order of the truncation errors (TEO) are derived for energies, transition moments, and property matrix elements of arbitrary excitation classes and truncation levels. In the analysis of the separability properties of the transition moments, the decisive role of the so-called dual ground state is revealed. Due to the use of CE states, the bCC approach can be compared to so-called intermediate state representation (ISR) methods based exclusively on suitably orthonormalized CE states. As the present analysis shows, the bCC approach has decisive advantages over the conventional CI treatment, but also distinctly weaker TEO and separability properties in comparison to a full (and hermitian) ISR method.     

Ground and valence-excited states of SF: a multireference configuration interaction with single and double excitations+Q study

Ab initio calculations on the ground and valence-excited states of the sulfur monofluoride radical have been performed using entirely uncontracted all-electron augmented correlation consistent polarized valence quintuple zeta basis sets and the internally contracted multireference configuration interaction with single and double excitations method and Davidson correction (+Q). Potential-energy curves of all valence electronic states and the spectroscopic constants of several bound states are fitted. It is the first time that the entire 27-omega states generated from the 12 valence lambda-S states which come from the S(3P(g)) and F(2P(u)) atomic states of SF radical have been studied theoretically. The effects of spin-orbit coupling and the avoided crossing rule between omega states of the same symmetry are analyzed. The calculated results reproduce well the available experimental values and predict the properties of several bound excited states that have never been observed in experiment. The transition properties of the dipole-allowed transitions from bound excited states to the ground state are predicted for the first time, including the transition dipole moments, the Franck-Condon factors, and the radiative lifetimes.     

Ab initio coupled cluster determination of the equilibrium structures of cis- and trans-1,2-difluoroethylene and 1,1-difluoroethylene

The equilibrium structures of cis- and trans-1,2-difluoroethylene and 1,1-difluoroethylene, C(2)H(2)F(2), have been determined with high-level coupled cluster techniques combined with large basis sets, explicit consideration of core/valence, and scalar relativistic and higher order correlation effects. Excellent agreement was found with new semiexperimental structures, increasing the level of confidence in both approaches. Differences in bond lengths among ethylene and the fluoroethylenes are discussed.     

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.     

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 lowest singlet-triplet excitation energy of BN: A converged coupled cluster perspective

The notoriously small X 3Pi-a 1Sigma+ excitation energy of the BN diatomic has been calculated using high-order coupled cluster methods. Convergence has been established in both the one-particle basis set and the coupled cluster expansion. Explicit inclusion of connected quadruple excitations T4 is required for even semiquantitative agreement with the limit value, while connected quintuple excitations T5 still have an effect of about 60 cm(-1). Still higher excitations only account for about 10 cm(-1). Inclusion of inner-shell correlation further reduces Te by about 60 cm(-1) at the CCSDT, and 85 cm(-1) at the CCSDTQ level. Our best estimate, Te = 183+/-40 cm(-1), is in excellent agreement with earlier calculations and experiment, albeit with a smaller (and conservative) uncertainty. The dissociation energy of BN(X 3Pi) is De = 105.74+/-0.16 kcal/mol and D0 = 103.57+/-0.16 kcal/mol.     

Electronic spectra and photodissociation of vinyl chloride: a symmetry-adapted cluster configuration interaction study

The vertical absorption spectrum and photodissociation mechanism of vinyl chloride (VC) were studied by using symmetry-adapted cluster configuration interaction theory. The important vertical pi --> pi* excitation was intensively examined with various basis sets up to aug-cc-pVTZ augmented with appropriate Rydberg functions. The excitation energy for pi --> pi* transition obtained in the present study, 6.96 eV, agrees well with the experimental value, 6.7-6.9 eV. Calculated excitation energies along with the oscillator strengths clarify that the main excitation in VC is the pi --> pi* excitation. Contrary to the earlier theoretical reports, the results obtained here support that the C-Cl bond dissociation takes place through the n(Cl-)sigma(C-Cl)* state.     

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     

Viscosity and structure configuration properties of equilibrium and supercooled liquid cobalt

The shear viscosity of liquid cobalt at the pressure p = 1.5 bar and at the temperatures corresponding to equilibrium liquid and supercooled liquid states is measured experimentally and evaluated by means of molecular dynamics simulations. Further, the shear viscosity is also calculated within the microscopic theoretical model. Comparison of our experimental, simulation and theoretical results with other available data allows one to examine the issue about the correct temperature dependence of the shear viscosity of liquid cobalt. A strong correlation between the viscosity and the configuration entropy of liquid cobalt over the considered temperature range is found, which can be taken into account by the Rosenfeld’s model.     

An efficient and accurate approximation to double substitution coupled cluster wavefunctions

An approximate CCD theory in an efficient electron pair operator form has been applied to He₂ with very accurate results.     

Third- and fourth-order perturbation corrections to excitation energies from configuration interaction singles

Complete third-order and partial fourth-order Rayleigh-Schrodinger perturbation corrections to excitation energies from configuration interaction singles (CIS) have been derived and termed CIS(3) and CIS(4)(P). They have been implemented by the automated system TENSOR CONTRACTION ENGINE into parallel-execution programs taking advantage of spin, spatial, and index permutation symmetries and applicable to closed- and open-shell molecules. The consistent use of factorization, first introduced by Head-Gordon et al. in the second-order correction to CIS denoted CIS(D), has reduced the computational cost of CIS(3) and CIS(4)(P) from O(n(8)) and O(n(6)) to O(n(6)) and O(n(5)), respectively, with n being the number of orbitals. It has also guaranteed the size extensivity of excited-state energies of these methods, which are in turn the sum of size-intensive excitation energies and the ground-state energies from the standard M?ller-Plesset perturbation theory at the respective orders. The series CIS(D), CIS(3), and CIS(4)(P) are usually monotonically convergent at values close to the accurate results predicted by coupled-cluster singles and doubles (CCSD) with a small fraction of computational costs of CCSD for predominantly singly excited states characterized by a 90%-100% overlap between the CIS and CCSD wave functions. When the overlap is smaller, the perturbation theory is incapable of adequately accounting for the mixing of the CIS states through higher-than-singles sectors of the Hamiltonian matrix, resulting in wildly oscillating series with often very large errors in CIS(4)(P). Hence, CIS(3) and CIS(4)(P) have a rather small radius of convergence and a limited range of applicability, but within that range they can be an inexpensive alternative to CCSD.     

Thermodynamics and equilibrium structure of Ne38 cluster: quantum mechanics versus classical

The equilibrium properties of classical Lennard-Jones (LJ38) versus quantum Ne38 Lennard-Jones clusters are investigated. The quantum simulations use both the path-integral Monte Carlo (PIMC) and the recently developed variational-Gaussian wave packet Monte Carlo (VGW-MC) methods. The PIMC and the classical MC simulations are implemented in the parallel tempering framework. The classical heat capacity Cv(T) curve agrees well with that of Neirotti et al. [J. Chem. Phys. 112, 10340 (2000)], although a much larger confining sphere is used in the present work. The classical Cv(T) shows a peak at about 6 K, interpreted as a solid-liquid transition, and a shoulder at approximately 4 K, attributed to a solid-solid transition involving structures from the global octahedral (Oh) minimum and the main icosahedral (C5v) minimum. The VGW method is used to locate and characterize the low energy states of Ne38, which are then further refined by PIMC calculations. Unlike the classical case, the ground state of Ne38 is a liquidlike structure. Among the several liquidlike states with energies below the two symmetric states (Oh and C5v), the lowest two exhibit strong delocalization over basins associated with at least two classical local minima. Because the symmetric structures do not play an essential role in the thermodynamics of Ne38, the quantum heat capacity is a featureless curve indicative of the absence of any structural transformations. Good agreement between the two methods, VGW and PIMC, is obtained. The present results are also consistent with the predictions by Calvo et al. [J. Chem. Phys. 114, 7312 (2001)] based on the quantum superposition method within the harmonic approximation. However, because of its approximate nature, the latter method leads to an incorrect assignment of the Ne38 ground state as well as to a significant underestimation of the heat capacity.