Algebraic Approach to Coupled Cluster Theory

REDUCED MULTIREFERENCE COUPLED-CLUSTER METHOD AND ITS APPLICATION TO THE PYRIDYNE DIRADICALS

The reduced multireference (RMR) coupled-cluster (CC) method with singles and doubles (RMR CCSD) that employs a modest-size MR CISD wave function as an external source for the most important (primary) triples and quadruples in order to account for the nondynamic correlation effects in the presence of quasidegeneracy, and which is further perturbatively corrected for the remaining (secondary) triples, RMR CCSD(T), is employed to compute the molecular geometry and the energy of the lowest-lying singlet and triplet states, as well as the corresponding singlet–triplet splitting, for all possible isomers of the m, n-pyridyne diradicals. A comparison is made with earlier results that were obtained by other authors, and the role of the multireference effects for both the geometry and the spin multiplicity of the lowest state, as described by the RMR-type methods, is demonstrated on the example of 2,6- and 3,5-pyridynes.

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.

Binding in transition metal complexes: Reduced multireference coupled-cluster study of the MCH2+ (M=Sc to Cu) compounds

The recently developed reduced multireference coupled-cluster method with singles and doubles (RMR CCSD), which is perturbatively corrected for triples [RMR CCSD(T)], is employed to compute binding energies of nine transition metal ions with CH2. Unlike analogous compounds involving main-group elements, the MCH2+ (M=Sc to Cu) transition metal complexes often exhibit a non-negligible multireference character. The authors thus employ the RMR CCSD(T) method, which represents an extension of the standard single-reference (SR) CCSD(T) method and can account for multireference effects, while employing only small reference spaces. In this way the role of quasidegeneracy effects on the binding energies of these complexes can be assessed at a higher SD(T) level than is possible with the widely used ab initio methods, namely, with the standard SR CCSD(T) approach, and provide a new benchmark for these quantities. The difference between the RMR and the standard CCSD(T) methods becomes particularly evident when considering nonequilibrium geometries.

Second order explicitly correlated R12 theory revisited: A second quantization framework for treatment of the operators’ partitionings

Second order R12 theory is presented and derived alternatively using the second quantized hole-particle formalism. We have shown that in order to ensure the strong orthogonality between the R12 and the conventional part of the wave function, the explicit use of projection operators can be easily avoided by an appropriate partitioning of the involved operators to parts which are fully describable within the computational orbital basis and complementary parts that involve imaginary orbitals from the complete orbital basis. Various Hamiltonian splittings are discussed and computationally investigated for a set of nine molecules and their atomization energies. If no generalized Brillouin condition is assumed, with all relevant partitionings the one-particle contribution arising in the explicitly correlated part of the first order wave function has to be considered and has a significant role when smaller atomic orbital basis sets are used. The most appropriate Hamiltonian splitting results if one follows the conventional perturbation theory for a general non-Hartree-Fock reference. Then, no couplings between the R12 part and the conventional part arise within the first order wave function. The computationally most favorable splitting when the whole complementary part of the Hamiltonian is treated as a perturbation fails badly. These conclusions also apply to MP2-F12 approaches with different correlation factors.

Correction for Triples in Reduced Multireference Coupled-Cluster Approaches

The performance of the recently proposed version of the reduced multireference (RMR) coupled-cluster (CC) method with singles and doubles (SD), which employs a modest-size configuration interaction wave function as an external source for a small subset of approximate connected three- and four-body cluster amplitudes that are primarily responsible for the nondynamic correlation effects, and which has been perturbatively corrected for the remaining triples along the same line as in the standard CCSD(T) method (Li X., Paldus J.: J. Chem. Phys. 2006, 124, 174101), referred to by the acronym RMR CCSD(T), is being tested by evaluating equilibrium spectroscopic constants for a demanding system of the beryllium dimer, as well as by computing atomization energies for several di- and triatomics. The focus is on the equilibrium properties, since it has been demonstrated earlier that the RMR CCSD method corrects well for the nondynamic correlation in bond-breaking situations. We find that in all the cases we have examined, the RMR CCSD(T) method does in fact improve the performance of CCSD(T) even in the vicinity of the equilibrium geometry. For states possessing a moderate multireference character, the improvement in computed thermochemical properties relative to CCSD(T) amounts to a few kJ/mol, a meaningful amount when striving for chemical accuracy.

A truncated version of reduced multireference coupled-cluster method with singles and doubles and noniterative triples: Application to F2 and Ni(CO)n (n=1, 2, and 4)

A perturbatively truncated version of the reduced multireference coupled-cluster method with singles and doubles and noniterative triples RMR CCSD(T) is described. In the standard RMR CCSD method, the effect of all triples and quadruples that are singles or doubles relative to references spanning a chosen multireference (MR) model space is accounted for via the external corrections based on the MR CISD wave function. In the full version of RMR CCSD(T), the remaining triples are then handled via perturbative corrections as in the standard, single-reference (SR) CCSD(T) method. By using a perturbative threshold in the selection of MR CISD configuration space, we arrive at the truncated version of RMR CCSD(T), in which the dimension of the MR CISD problem is significantly reduced, thus leaving more triples to be treated perturbatively. This significantly reduces the computational cost. We illustrate this approach on the F2 molecule, in which case the computational cost of the truncated version of RMR CCSD(T) is only about 10%-20% higher than that of the standard CCSD(T), while still eliminating the failure of CCSD(T) in the bond breaking region of geometries. To demonstrate the capabilities of the method, we have also used it to examine the structure and binding energy of transition metal complexes Ni(CO)n with n=1, 2, and 4. In particular, Ni(CO)2 is shown to be bent rather than linear, as implied by some earlier studies. The RMR CCSD(T) binding energy differs from the SR CCSD(T) one by 1-2 kcal/mol, while the energy barrier separating the linear and bent structures of Ni(CO)2 is smaller than 1 kcal/mol.

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.

Multireference Brillouin-Wigner coupled clusters method with noniterative perturbative connected triples

We developed and implemented the state-specific Brillouin-Wigner coupled cluster method with singles, doubles, and noniterative perturbative triples, called MR BWCCSD(T), for a general number of closed- and open-shell reference configurations. To assess the accuracy of the method, we performed calculations of the three lowest electronic states of the oxygen molecule and of the automerization barrier of cyclobutadiene. For the oxygen molecule, the results were in a good agreement in comparison with those of the iterative MR BWCCSDTalpha method. For cyclobutadiene, the effect of connected triples was found to be minor, which is in agreement with the previous study by and Balková and Bartlett [J. Chem. Phys. 101, 8972 (1994)].

General-model-space state-universal coupled-cluster methods for excited states: Diagonal noniterative triple corrections

The recently developed multireference, general-model-space, state-universal coupled-cluster approach considering singles and doubles (GMS SU CCSD) has been extended to account perturbatively for triples, similar to the ubiquitous single-reference CCSD(T) method. The effectiveness of this extension in handling of excited states and its ability to account for the static and nondynamic correlation effects when considering spin- and/or space-symmetry degenerate levels within the spin-orbital formalism is examined on the example of low-lying excitation energies of the C2, N2, and CO molecules and a comparison is made with the (N,N)-CCSD method used for the same purpose. It is shown that while the triple corrections are very effective in improving the absolute energies, they have only a modest effect on the corresponding excitation energies, which may be even detrimental if both the ground- and excited-state levels cannot be given a balanced treatment. While the triple corrections help to avoid the symmetry-breaking effects arising due to the use of the spin-orbital formalism, they are much less effective in this regard than the (N,N)-CCSD approach.

Accurate quantum-chemical prediction of enthalpies of formation of small molecules in the gas phase

The coupled-cluster approach, including single and double excitations and perturbative corrections for triple excitations, is capable of predicting molecular electronic energies and enthalpies of formation of small molecules in the gas phase with very high accuracy (specifically, with error bars less than 5 kJmol-1), provided that the electronic wavefunction is dominated by the Hartree-Fock configuration. This capability is illustrated by calculations on molecules containing O-H and O-F bonds, namely OH, FO, H2O, HOF, and F2O. To achieve this very high accuracy, it is imperative to account for electron-correlation effects in a quantitative manner, either by using explicitly correlated two-particle basis functions (R12 functions) or by extrapolating to the limit of a complete basis. Besides taking into account harmonic zero-point vibrational energies, it is also necessary to account for anharmonic corrections to the zero-point vibrational energies, to include the core orbitals into the coupled-cluster calculations, and to account for spin-orbit corrections and scalar relativistic effects. These additional corrections constitute small but significant contributions in the range of 1-4 kJmol-1 to the enthalpies of formation of the aforementioned molecules. The highly accurate coupled-cluster results, obtained by employing R12 functions and by including various corrections, are compared with standard Kohn-Sham density-functional calculations as well as with the Gaussian-2 and complete-basis-set model chemistries.