Partially linearized, fully size-extensive, and reduced multireference coupled-cluster methods. II. Applications and performance

The partially linearized (pl), fully size-extensive multireference (MR) coupled-cluster (CC) method, fully accounting for singles (S) and doubles (D) and approximately for a subset of primary higher than doubles, referred to as plMR CCSD, as well as its plMR CCSD(T) version corrected for secondary triples, as described in Part I of this paper [X. Li and J. Paldus, J. Chem. Phys. 128, 144118 (2008)], are applied to the problem of bond breaking in the HF, F2, H2O, and N2 molecules, as well as to the H4 model, using basis sets of a DZ or a cc-pVDZ quality that enable a comparison with the full configuration interaction (FCI) exact energies for a given ab initio model. A comparison of the performance of the plMR CCSD/CCSD(T) approaches with those of the reduced MR (RMR) CCSD/CCSD(T) methods, as well as with the standard single reference (SR) CCSD and CCSD(T) methods, is made in each case. For the H4 model and N2 we also compare our results with the completely renormalized (CR) CC(2,3) method [P. Piecuch and M. W?och, J. Chem. Phys. 123, 224105 (2005)]. An important role of a proper choice of the model space for the MR-type methods is also addressed. The advantages and shortcomings of all these methods are pointed out and discussed, as well as their size-extensivity characteristics, in which case we distinguish supersystems involving noninteracting SR and MR subsystems from those involving only MR-type subsystems. Although the plMR-type approaches render fully size-extensive results, while the RMR CCSD may slightly violate this property, the latter method yields invariably superior results to the plMR CCSD ones and is more easy to apply in highly demanding cases, such as the triple-bond breaking in the nitrogen molecule.     

Full potential energy curve for N2 by the reduced multireference coupled-cluster method

Relying on a 56-dimensional reference space and using up to the correlation-consistent, polarized, valence-quadruple-zeta (cc-pVQZ) basis sets, the reduced multireference (RMR) coupled-cluster method with singles and doubles (CCSD), as well as its perturbatively corrected version for secondary triples [RMR CCSD(T)], is employed to generate the full potential energy curves for the nitrogen molecule. The resulting potentials are then compared to the recently published accurate analytic potential based on an extensive experimental data analysis [R. J. Le Roy et al., J. Chem. Phys. 125, 164310 (2006)], and the vibrational term values of these potentials are compared over the entire well. A comparison with single-reference CCSD and CCSD(T) results, as well as with earlier obtained eight-reference RMR CC results, is also made. Excellent performance of RMR CCSD, and its systematic improvement with the increasing dimension of the reference space employed, is demonstrated. For the first 19 vibrationally excited levels, which are based on experimentally observed bands, we find an absolute average deviation of 8 cm(-1) from the computed RMR CCSD/cc-pVQZ values. The perturbative correction for triples increases this deviation to 126 cm(-1), but only to 61 cm(-1) when extrapolated to the basis set limit. Both RMR CCSD and RMR CCSD(T) potentials perform well when compared to the experiment-based analytic potential in the entire range of internuclear separations.     

Performance of the general-model-space state-universal coupled-cluster method

The capabilities of the recently developed multireference, general-model-space (GMS), state-universal (SU) coupled-cluster (CC) method have been extended in order to enable the handling of any excited state that represents a single (S) or a double (D) excitation relative to the ground state. A series of calculations concerning the ground and excited states of the CH(+), HF, F(2), H(2)O, NH(2), and CH(2) molecules were carried out so as to assess the performance of the GMS SU CCSD method. For diatomics we have computed the entire potential energy curves, while for triatomics we have focused on vertical excitation energies. We demonstrate how a systematic enlargement of the model space enables a consideration of a larger and larger number of excited states. A comparison of the CC and full configuration interaction or large-scale CI results enables an assessment of the accuracy and reliability of the GMS SU CCSD method within a given basis set. In all cases very good results have been obtained, including highly excited states and those having a doubly-excited character.     

Partially linearized, fully size-extensive, and reduced multireference coupled-cluster methods. I. Formalism and mutual relationship

We describe a fully size-extensive alternative of the reduced multireference (RMR) coupled-cluster (CC) method with singles (S) and doubles (D) that generates a subset of higher-than-pair cluster amplitudes, using linearized CC equations from the full CC chain, projected onto the corresponding higher-than-doubly excited configurations. This approach is referred to as partially linearized (pl) MR CCSD method and characterized by the acronym plMR CCSD. In contrast to a similar CCSDT-1 method [Y. S. Lee et al., J. Chem. Phys. 81, 5906 (1984)] this approach also considers higher than triples (currently up to hexuples), while focusing only on a small subset of such amplitudes, referred to as the primary ones. These amplitudes are selected using similar criteria as in RMR CCSD. An extension considering secondary triples via the standard (T)-type corrections, resulting in the plMR CCSD(T) method, is also considered. The relationship of RMR and plMR CCSD and CCSD(T) approaches is discussed, and their performance and characteristics are the subject of the subsequent Part II of this paper.     

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.     

Benchmarks for electronically excited states: CASPT2, CC2, CCSD, and CC3

A benchmark set of 28 medium-sized organic molecules is assembled that covers the most important classes of chromophores including polyenes and other unsaturated aliphatic compounds, aromatic hydrocarbons, heterocycles, carbonyl compounds, and nucleobases. Vertical excitation energies and one-electron properties are computed for the valence excited states of these molecules using both multiconfigurational second-order perturbation theory, CASPT2, and a hierarchy of coupled cluster methods, CC2, CCSD, and CC3. The calculations are done at identical geometries (MP26-31G*) and with the same basis set (TZVP). In most cases, the CC3 results are very close to the CASPT2 results, whereas there are larger deviations with CC2 and CCSD, especially in singlet excited states that are not dominated by single excitations. Statistical evaluations of the calculated vertical excitation energies for 223 states are presented and discussed in order to assess the relative merits of the applied methods. CC2 reproduces the CC3 reference data for the singlets better than CCSD. On the basis of the current computational results and an extensive survey of the literature, we propose best estimates for the energies of 104 singlet and 63 triplet excited states.     

Coupled cluster methods including triple excitations for excited states of radicals

We report an extension of the coupled cluster iterative-triples model, CC3, to excited states of open-shell molecules, including radicals. We define the method for both spin-unrestricted Hartree-Fock (UHF) and spin-restricted open-shell Hartree-Fock (ROHF) reference determinants and discuss its efficient implementation in the PSI3 program package. The program is streamlined to use at most O(N(7)) computational steps and avoids storage of the triple-excitation amplitudes for both the ground- and excited-state calculations. The excitation-energy program makes use of a Lowdin projection formalism (comparable to that of earlier implementations) that allows computational reduction of the Davidson algorithm to only the single- and double-excitation space, but limits the calculation to only one excited state at a time. However, a root-following algorithm may be used to compute energies for multiple states of the same symmetry. Benchmark applications of the new methods to the lowest valence (2)B(1) state of the allyl radical, low-lying states of the CH and CO(+) diatomics, and the nitromethyl radical show substantial improvement over ROHF- and UHF-based CCSD excitation energies for states with strong double-excitation character or cases suffering from significant spin contamination. For the allyl radical, CC3 adiabatic excitation energies differ from experiment by less than 0.02 eV, while for the (2)Sigma(+) state of CH, significant errors of more than 0.4 eV remain.     

Real or artifactual symmetry breaking in the BNB radical: a multireference coupled cluster viewpoint

The ground state of the linear BNB radical has been examined via the recently developed reduced multireference coupled cluster method with singles and doubles that is perturbatively corrected for triples [RMR CCSD(T)] using the correlation consistent basis sets (cc-pVXZ, X=D, T, and Q). Similar to earlier results that were based on the single reference CCSD(T) and BD(T) approaches, the RMR CCSD(T) method also predicts an asymmetric structure with two BN bonds of unequal length, even though the MR effects significantly reduce the barrier height. The computed frequencies for the symmetric and antisymmetric stretching modes agree reasonably well with the experimental data.     

Quantum-chemical investigation of the structures and electronic spectra of the nucleic acid bases at the coupled cluster CC2 level

We study the ground-state structures and singlet- and triplet-excited states of the nucleic acid bases by applying the coupled cluster model CC2 in combination with a resolution-of-the-identity approximation for electron interaction integrals. Both basis set effects and the influence of dynamic electron correlation on the molecular structures are elucidated; the latter by comparing CC2 with Hartree-Fock and M?ller-Plesset perturbation theory to second order. Furthermore, we investigate basis set and electron correlation effects on the vertical excitation energies and compare our highest-level results with experiment and other theoretical approaches. It is shown that small basis sets are insufficient for obtaining accurate results for excited states of these molecules and that the CC2 approach to dynamic electron correlation is a reliable and efficient tool for electronic structure calculations on medium-sized molecules.     

Full configuration interaction potential energy curves for the X (1)Sigma(g) (+), B (1)Delta(g), and B(') (1)Sigma(g) (+) states of C(2): a challenge for approximate methods

The C(2) molecule exhibits unusual bonding and several low-lying excited electronic states, making the prediction of its potential energy curves a challenging test for quantum chemical methods. We report full configuration interaction results for the X (1)Sigma(g) (+), B (1)Delta(g), and B(') (1)Sigma(g) (+) states of C(2), which exactly solve the electronic Schrodinger equation within the space spanned by a 6-31G( *) basis set. Within the D(2h) subgroup used by most electronic structure programs, these states all have the same symmetry ((1)A(g)), and all three states become energetically close for interatomic distances beyond 1.5 A. The quality of several single-reference ab initio methods is assessed by comparison to the benchmark results. Unfortunately, even coupled-cluster theory through perturbative triples using an unrestricted Hartree-Fock reference exhibits large nonparallelity errors (>20 kcal mol(-1)) for the ground state. The excited states are not accurately modeled by any commonly used single-reference method, nor by configuration interaction including full quadruple substitutions. The present benchmarks will be helpful in assessing theoretical methods designed to break bonds in ground and excited electronic states.     

On the nature of the low-lying singlet states of 4-(Dimethyl-amino)benzonitrile

4-(N,N-Dimethyl-amino)benzonitrile (DMABN) is a prototype molecule for dual fluorescence. The anomalous emission has been attributed to an intramolecular charge-transfer (ICT) state, and the structure of the latter is still subject to some controversy. We applied a recently developed analytical gradient code for the approximate coupled-cluster singles-and-doubles method CC2 in combination with accurate basis sets to address this problem. Fully optimized excited state structures are presented for the ICT state and the so-called locally excited state, and recent transient IR and Raman measurements on the excited states are interpreted by means of calculated harmonic frequencies. Strong evidence is found for an electronic decoupling of the phenyl and the dimethyl-amino moiety, resulting in a minimum structure for the ICT state with a twisted geometry. In contrast to previous findings, the structure of this state is, at least in the gas phase, not C(2v) symmetric but distorted towards C(s) symmetry. The distortion coordinate is a pyramidalization of the phenyl carbon atom carrying the dimethyl-amino group. The results from the CC2 model are supported by single-point calculations using more elaborate coupled-cluster models (CCSD, CCSDR(3)) and by CASSCF calculations.     

Quantum chemical study of the electronic structure of NiCH2 + in its ground state and low-lying electronic excited states

The electronic structure of NiCH(2) (+), representative of transition metal carbene ions, is investigated by means of several methods of quantum chemistry. The relative stabilities of the four low-lying doublet electronic states ((2)A(1), (2)A(2), (2)B(1), and (2)B(2)) are determined at the coupled cluster singles and doubles level (CCSD) and triples level [CCSD(T) and CCSDT-3] with both a Hartree-Fock and density functional theory (Kohn-Sham) reference. The equation-of-motion coupled cluster for treatment of excited states in singles and doubles approximation (EOM-CCSD) is used to characterize the transition energies from the (2)A(1) electronic ground state to the low-lying doublet excited states. The (2)A(2) and (2)B(1) states are nearly degenerate, found to be separated by 940 cm(-1) at the EOM-CCSD level, in agreement with the CASSCF energy ordering. The (2)B(2) state is calculated to be higher in energy by more than 1.0 eV. The spin purity of the low-lying doublet and quadruplet states described by CCSD calculations based on the unrestricted open-shell Hartree-Fock reference is discussed.     

The coupled cluster singles, doubles, and a hybrid treatment of connected triples based on the split virtual orbitals

We have proposed a simple strategy for splitting the virtual orbitals with a large basis set into two subgroups (active and inactive) by taking a smaller basis set as an auxiliary basis set. With the split virtual orbitals (SVOs), triple or higher excitations can be partitioned into active and inactive subgroups (according to the number of active virtual orbitals involved), which can be treated with different electron correlation methods. In this work, the coupled cluster (CC) singles, doubles, and a hybrid treatment of connected triples based on the SVO [denoted as SVO-CCSD(T)-h], has been implemented. The present approach has been applied to study the bond breaking potential energy surfaces in three molecules (HF, F(2), and N(2)), and the equilibrium properties in a number of open-shell diatomic molecules. For all systems under study, the SVO-CCSD(T)-h method based on the unrestricted Hartree-Fock (UHF) reference is an excellent approximation to the corresponding CCSDT (CC singles, doubles, and triples), and much better than the UHF-based CCSD(T) (CC singles, doubles, and perturbative triples). On the other hand, the SVO-CCSD(T)-h method based on the restricted HF (RHF) reference can also provide considerable improvement over the RHF-based CCSD(T).     

Real versus artifactual symmetry-breaking effects in Hartree-Fock, density-functional, and coupled-cluster methods

We have examined the relative abilities of Hartree-Fock, density-functional theory (DFT), and coupled-cluster theory in describing second-order (pseudo) Jahn-Teller (SOJT) effects, perhaps the most commonly encountered form of symmetry breaking in polyatomic molecules. As test cases, we have considered two prototypical systems: the 2Sigmau+ states of D( infinity h) BNB and C3+ for which interaction with a low-lying 2Sigmag+ excited state leads to symmetry breaking of the nuclear framework. We find that the Hartree-Fock and B3LYP methods correctly reproduce the pole structure of quadratic force constants expected from exact SOJT theory, but that both methods appear to underestimate the strength of the coupling between the electronic states. Although the Tamm-Dancoff (CIS) approximation gives excitation energies with no relationship to the SOJT interaction, the random-phase-approximation (RPA) approach to Hartree-Fock and time-dependent DFT excitation energies predicts state crossings coinciding nearly perfectly with the positions of the force constant poles. On the other hand, the RPA excited-state energies exhibit unphysical curvature near their crossings with the ground (reference) state, a problem arising directly from the mathematical structure of the RPA equations. Coupled-cluster methods appear to accurately predict the strength of the SOJT interactions between the 2Sigmau+ and 2Sigmag+ states, assuming that the inclusion of full triple excitations provides a suitable approximation to the exact wave function, and are the only methods examined here which predict symmetry breaking in BNB. However, coupled-cluster methods are plagued by artifactual force constant poles arising from the response of the underlying reference molecular orbitals to the geometric perturbation. Furthermore, the structure of the "true" SOJT force constant poles predicted by coupled-cluster methods, although correctly positioned, has the wrong structure.     

Multi-reference state-universal coupled-cluster approaches to electronically excited states

The multi-reference (MR), general model space (GMS), state-universal (SU), coupled-cluster (CC) method with singles and doubles (GMS-SU-CCSD), as well as its triple-corrected versions GMS-SU-CCSD(T), are employed to assess their ability to describe low-lying excited states of various molecules, with an emphasis on a simultaneous handling of several states of the same symmetry species. A special attention is given to the role of the so-called C-conditions that account for non-vanishing internal cluster amplitudes when relying on an incomplete GMS, as well as to the choice of suitable model spaces and a perturbative account of secondary triples. The ambiguities arising when using large basis sets are also pointed out. To achieve a general assessment of the potential of the GMS-type SU-CC approaches, the vertical excitation energies of several species, including the challenging BN diatomic as well as larger systems, namely formaldehyde, trans-butadiene, formamide, and benzene are considered. These results are compared with those provided by the equation-of-motion EOM-CCSD method and, whenever available, the density functional theory results and experimental data. These comparisons clearly demonstrate the usefulness of GMS-type MR-CC approaches.     

Mechanistic insights into NO3* induced self-terminating radical oxygenations, part 1: a computational study on NO3* and its addition to alkynes

The geometry of the nitrate radical, NO3*, for which unrestricted Hartree-Fock (HF) breaks spatial symmetry of the wave function, was optimized using hybrid density functionals that include varying fractions of Hartree-Fock exchange. Although symmetry breaking was not observed even when the functional with the highest HF exchange (BHandHLYP) was used, only B3LYP correctly describes the D(3h) symmetry of NO3* as ground-state structure with the lowest energy. Further, geometries and energies of the stationary points in the addition of NO3* to ethyne, propyne, and 2-butyne were calculated using ab initio and density functional methods. The reactions proceed through Z-configurated transition states leading to Z-configurated vinyl radicals with the activation barrier decreasing with increasing methyl substitution at the C[triple bond]C by ca. 11 kJ mol(-1) per methyl group. It was found that the results obtained at the BHandHLYP/cc-pVDZ level of theory are in good agreement with the data from single-point QCISD and CCSD(T) calculations.     

Theoretical study of the electronic structure of MCH(2)(+)(M=Fe,Co,Ni)

State of the art coupled cluster (CC) methods are applied to accurately characterize the ground state electronic structure and photoelectron spectra of transition metal carbene ions MCH(2) (+) (M=Fe, Co, and Ni). The geometries and energies of the lowest energy quartet, triplet, and doublet electronic states as well as several low-lying vertical excitation energies of FeCH(2) (+), CoCH(2) (+), and NiCH(2) (+) are reported. The excitation energies are computed using the equation-of-motion CC and for states of different symmetries, by the energy differences of single reference ground and excited states (Delta-CC). The latter employ several reference states; the unrestricted Hartree-Fock, restricted open shell Hartree-Fock, and unrestricted Kohn-Sham. We conclude that the (2)A(1) electronic ground state of NiCH(2) (+) is separated by about 30.0 kJ/mol from the next highest state, and the lowest (4)B(1) and (4)B(2) states of FeCH(2) (+) as well as the (3)A(2) and (3)A(1) states of CoCH(2) (+) are nearly degenerate. The presence of metal-pi*(MCH(2)) charge transfer states with significant oscillator strengths in the visible/near-UV energy domain of the theoretical spectra of FeCH(2) (+) and CoCH(2) (+) are at the origin of the photofragmentation of these compounds observed after irradiation between 310 and 360 nm.     

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

Variational wavefunctions for low-lying excited states

The multiconfiguration method based on the generalized Brillouin theorem is well suited to optimize orbitals in variational wavefunctions for low-lying excited states of a given symmetry. Such wavefunctions are constrained to be orthogonal to and noninteracting with the wavefunctions for all lower states of the same symmetry. Test calculations were performed on the lowest excited ¹S state of Be. For a Hartree-Fock ground state wavefunction, singly excited configurations were insufficient to describe the lowest excited state, and triply excited configurations had to be added. For multiconfiguration ground state wavefunctions, however, singly excited configurations gave good results.