The barrier height of the F+H2 reaction revisited: coupled-cluster and multireference configuration-interaction benchmark calculations

Large scale coupled-cluster benchmark calculations have been carried out to determine the barrier height of the F+H2 reaction as accurately as possible. The best estimates for the barrier height of the linear and bent transition states amount to 2.16 and 1.63 kcal/mol, respectively. These values include corrections for core correlation, scalar-relativistic effects, spin-orbit effects, as well as the diagonal Born-Oppenheimer correction. The CCSD(T) basis-set limits are estimated using extrapolation techniques with augmented quintuple and sextuple-zeta basis sets, and remaining N-electron errors are determined using coupled-cluster singles, doubles, triples, quadruples calculations with up to augmented quintuple-zeta basis sets. The remaining uncertainty is estimated to be less than 0.1 kcal/mol. The coupled-cluster results are used to calibrate multireference configuration-interaction calculations with empirical scaling of the correlation energy.     

Extrapolating to the one-electron basis set limit in polarizability calculations

We report calculations of polarizabilities using total energies extrapolated to the complete basis set limit. A dual-level scheme has been employed, with the complete basis set limit of the correlation energy determined by the recently reported uniform singlet- and triplet-pair extrapolation method. The finite field approach has been employed, with tensors and averaged polarizabilities for the ground electronic states of H 2, N 2, CO, and H 2O reported and compared with available experimental data in the literature. Exploratory results are also presented for C 6H 4NO 2NH 2.     

Linear-response theory for Mukherjee's multireference coupled-cluster method: Excitation energies

The recently presented linear-response function for Mukherjee's multireference coupled-cluster method (Mk-MRCC) [T.-C. Jagau and J. Gauss, J. Chem. Phys. 137, 044115 (2012)] is employed to determine vertical excitation energies within the singles and doubles approximation (Mk-MRCCSD-LR) for ozone as well as for o-benzyne, m-benzyne, and p-benzyne, which display increasing multireference character in their ground states. In order to assess the impact of a multireference ground-state wavefunction on excitation energies, we compare all our results to those obtained at the single-reference coupled-cluster level of theory within the singles and doubles as well as within the singles, doubles, and triples approximation. Special attention is paid to the artificial splitting of certain excited states which arises from the redundancy intrinsic to Mk-MRCC theory and hinders the straightforward application of the Mk-MRCC-LR method.     

Coupled-cluster methods with perturbative inclusion of explicitly correlated terms: a preliminary investigation

We propose to account for the large basis-set error of a conventional coupled-cluster energy and wave function by a simple perturbative correction. The perturbation expansion is constructed by L?wdin partitioning of the similarity-transformed Hamiltonian in a space that includes explicitly correlated basis functions. To test this idea, we investigate the second-order explicitly correlated correction to the coupled-cluster singles and doubles (CCSD) energy, denoted here as the CCSD(2)(R12) method. The proposed perturbation expansion presents a systematic and easy-to-interpret picture of the "interference" between the basis-set and correlation hierarchies in the many-body electronic-structure theory. The leading-order term in the energy correction is the amplitude-independent R12 correction from the standard second-order M?ller-Plesset R12 method. The cluster amplitudes appear in the higher-order terms and their effect is to decrease the basis-set correction, in accordance with the usual experience. In addition to the use of the standard R12 technology which simplifies all matrix elements to at most two-electron integrals, we propose several optional approximations to select only the most important terms in the energy correction. For a limited test set, the valence CCSD energies computed with the approximate method, termed , are on average precise to (1.9, 1.4, 0.5 and 0.1%) when computed with Dunning's aug-cc-pVXZ basis sets [X = (D, T, Q, 5)] accompanied by a single Slater-type correlation factor. This precision is a roughly an order of magnitude improvement over the standard CCSD method, whose respective average basis-set errors are (28.2, 10.6, 4.4 and 2.1%). Performance of the method is almost identical to that of the more complex iterative counterpart, CCSD(R12). The proposed approach to explicitly correlated coupled-cluster methods is technically appealing since no modification of the coupled-cluster equations is necessary and the standard M?ller-Plesset R12 machinery can be reused.     

Fermi resonance in CO2: a combined electronic coupled-cluster and vibrational configuration-interaction prediction

The authors present a first-principles prediction of the energies of the eight lowest-lying anharmonic vibrational states of CO(2), including the fundamental symmetric stretching mode and the first overtone of the fundamental bending mode, which undergo a strong coupling known as Fermi resonance. They employ coupled-cluster singles, doubles, and (perturbative) triples [CCSD(T) and CCSDT] in conjunction with a range of Gaussian basis sets (up to cc-pV5Z, aug-cc-pVQZ, and aug-cc-pCVTZ) to calculate the potential energy surfaces (PESs) of the molecule, with the errors arising from the finite basis-set sizes eliminated by extrapolation. The resulting vibrational many-body problem is solved by the vibrational self-consistent-field and vibrational configuration-interaction (VCI) methods with the PESs represented by a fourth-order Taylor expansion or by numerical values on a Gauss-Hermite quadrature grid. With the VCI, the best theoretical estimates of the anharmonic energy levels agree excellently with experimental values within 3.5 cm(-1) (the mean absolute deviation). The theoretical (experimental) anharmonic frequencies of the Fermi doublet are 1288.9 (1285.4) and 1389.3 (1388.2) cm(-1).     

Coupled-cluster studies of the electronic excitation spectra of silanes

The electronic excitation spectra of unsubstituted linear silanes (n-Si(m)H(2m+2), m = 1-6), cyclopentasilane (c-Si5H10), and neopentasilane (neo-Si5H12) have been studied at the coupled-cluster approximate singles and doubles (CC2) level using Dunning's quadruple-zeta basis sets augmented with diffuse functions (aug-cc-pVQZ). Comparisons with measured ultraviolet spectra for Si2H6 and n-Si3H8 show that CC2 calculations using these basis sets yield excitation energies in good agreement with experiment. The calculated excitation thresholds for Si2H6 and n-Si3H8 of 7.61 and 6.68 eV are only 0.05 eV larger than the gas-phase values of 7.56 and 6.63 eV, respectively. For n-Si4H10, n-Si5H12, and neo-Si5H12, the calculated excitation thresholds of 6.51, 6.14, and 6.87 eV for the lowest dipole-allowed transitions are about 0.4 eV larger than the corresponding liquid-phase data of 6.05, 5.77, and 6.53 eV; the discrepancy can mainly be attributed to solvent effects. The obtained excitation thresholds for n-Si6H14 is 5.85 eV, whereas no experimental data are available for its optical gap. Calculations using the Karlsruhe triple-zeta valence basis sets augmented with single and double sets of polarization functions show that very large basis sets augmented with diffuse functions are needed for obtaining accurate excitation energies. The optical gaps for silanes obtained using the triple-zeta polarization basis sets were found to be 0.4 and 0.2 eV larger than those obtained using Dunning's quadruple-zeta basis sets. Excitation thresholds calculated at density functional theory levels using generalized gradient approximation are 0.7-1.0 eV smaller than the experimental values and by employing hybrid functionals they are 0.3-0.4 eV below the experimental thresholds. By adding the present basis-set correction and environmental effects to the previously calculated CC2 value for the excitation threshold of the Si29H36 silicon nanocluster, the extrapolated absorption threshold is 4.0 eV as compared to the recently reported experimental value of 3.7 eV.     

Generalized uniform singlet- and triplet-pair extrapolation of the correlation energy to the one electron basis set limit

The relationship between the triplet- and singlet-pair interaction coefficients in the uniform singlet- and triplet-pair extrapolation method recently suggested for extrapolating ab initio energies to the one-electron basis set limit is analyzed. Based on the premise that such a ratio is invariant over the configuration space of the molecule, generalizations of the method are investigated and their performance tested on extrapolations with MP2, CCD, CCSD, and MRCI(Q) energies. The best variant requires raw energies calculated using augmented correlated consistent basis sets with cardinal numbers up to X = 6 at a single geometry. A scheme is also suggested that performs better than the traditional X(-3) law and possibly the original uniform singlet- and triplet-pair extrapolation method but requires data only up to X = Q.     

Is HO3 minimum cis or trans? An analytic full-dimensional ab initio isomerization path

The minimum energy path for isomerization of HO(3) has been explored in detail using accurate high-level ab initio methods and techniques for extrapolation to the complete basis set limit. In agreement with other reports, the best estimates from both valence-only and all-electron single-reference methods here utilized predict the minimum of the cis-HO(3) isomer to be deeper than the trans-HO(3) one. They also show that the energy varies by less than 1 kcal mol(-1) or so over the full isomerization path. A similar result is found from valence-only multireference configuration interaction calculations with the size-extensive Davidson correction and a correlation consistent triple-zeta basis, which predict the energy difference between the two isomers to be of only Δ = -0.1 kcal mol(-1). However, single-point multireference calculations carried out at the optimum triple-zeta geometry with basis sets of the correlation consistent family but cardinal numbers up to X = 6 lead upon a dual-level extrapolation to the complete basis set limit of Δ = (0.12 ± 0.05) kcal mol(-1). In turn, extrapolations with the all-electron single-reference coupled-cluster method including the perturbative triples correction yield values of Δ = -0.19 and -0.03 kcal mol(-1) when done from triple-quadruple and quadruple-quintuple zeta pairs with two basis sets of increasing quality, namely cc-cpVXZ and aug-cc-pVXZ. Yet, if added a value of 0.25 kcal mol(-1) that accounts for the effect of triple and perturbative quadruple excitations with the VTZ basis set, one obtains a coupled cluster estimate of Δ = (0.14 ± 0.08) kcal mol(-1). It is then shown for the first time from systematic ab initio calculations that the trans-HO(3) isomer is more stable than the cis one, in agreement with the available experimental evidence. Inclusion of the best reported zero-point energy difference (0.382 kcal mol(-1)) from multireference configuration interaction calculations enhances further the relative stability to ΔE(ZPE) = (0.51 ± 0.08) kcal mol(-1). A scheme is also suggested to model the full-dimensional isomerization potential-energy surface using a quadratic expansion that is parametrically represented by a Fourier analysis in the torsion angle. The method illustrated at the raw and complete basis-set limit coupled-cluster levels can provide a valuable tool for a future analysis of the available (incomplete thus far) experimental rovibrational data.     

Calculated electronic transitions of the water ammonia complex

We have calculated vertical excitation energies and oscillator strengths of the low lying electronic transitions in H2O, NH3, and H2ONH3 using a hierarchy of coupled cluster response functions [coupled cluster singles (CCS), second order approximate coupled cluster singles and doubles (CC2), coupled cluster singles and doubles (CCSD), and third order approximate coupled cluster singles, doubles, and triples (CC3)] and correlation consistent basis functions (n-aug-cc-pVXZ, where n=s,d,t and X=D,T,Q). Our calculations indicate that significant changes in the absorption spectra of the photodissociative states of H2O and NH3 monomers occur upon complexation. In particular, we find that the electronic transitions originating from NH3 are blueshifted, whereas the electronic transitions originating from H2O are redshifted.     

Accurate benchmark calculation of the reaction barrier height for hydrogen abstraction by the hydroperoxyl radical from methane. Implications for C(n)H(2n+2) where n = 2 --> 4

The CH4 + HO2(*) reaction is studied by using explicitly correlated coupled-cluster theory with singles and doubles (CCSD-R12) in a large 19s14p8d6f4g3h basis (9s6p4d3f for H) to approach the basis-set limit at the coupled-cluster singles-doubles level. A correction for connected triple excitations is obtained from the conventional CCSD(T) coupled-cluster approach in the correlation-consistent quintuple-zeta basis (cc-pV5Z). The highly accurate results for the methane reaction are used to calibrate the calculations of the hydroperoxyl-radical hydrogen abstraction from other alkanes. For the alkanes C(n)H(2n+2) with n = 2 --> 4, the reactions are investigated at the CCSD(T) level in the correlation-consistent triple-zeta (cc-pVTZ) basis. The results are adjusted to the benchmark methane reaction and compared with those from other approaches that are commonly used in the field such as CBS-QB3, CBS-APNO, and density functional theory. Rate constants are computed in the framework of transition state theory, and the results are compared with previous values available.     

Extrapolation to the limit of a complete basis set for electronic structure calculations on the N2 molecule

Results obtained from nonrelativistic electronic structure calculations using finite Gaussian basis sets are extrapolated to the limit of a complete basis set, employing the results of explicitly correlated coupled-cluster calculations including singles and doubles substitutions (CCSD). For N₂, the basis-set limits for the electronic binding energy, equilibrium bond length and harmonic vibrational wave number are established for the CCSD model including a perturbative correction for triples substitutions and for the internally contracted multireference configuration interaction method. The resulting numbers are in good agreement with experimental values. Received: 2 December 1997 / Accepted: 3 February 1998 / Published online: 17 June 1998     

Ab initio calculation of the low-lying vibrational states of C2H2(A) in full dimensionality

We report full-dimensional calculations of vibrational energies of trans-C2H2(A) using the code MULTIMODE and with a full-dimensional potential energy surface obtained by fitting singles and doubles coupled-cluster equations-of-motion (EOM-CCSD) energies using a [3s 2p 1d] atomic natural orbital basis. The EOM-CCSD calculations were done with the code "ACES II". We compare the properties of the potential surface to previous calculations at the trans minimum and also compare the vibrational energies to experimental ones.     

Basis set limits of the second order Moller-Plesset correlation energies of water, methane, acetylene, ethylene, and benzene

We report second order Moller-Plesset (MP2) and MP2-F12 total energies on He, Ne, Ar, H(2)O, CH(4), C(2)H(2), C(2)H(4), and C(6)H(6), using the correlation consistent basis sets, aug-cc-pVXZ (X=D-7). Basis set extrapolation techniques are applied to the MP2 and MP2-F12/B methods. The performance of the methods is tested in the calculations of the atoms, He, Ne, and Ar. It is indicated that the two-point extrapolation of MP2-F12/B with the basis sets (X=5,6) is the most reliable. Similar accuracy is obtained using two-point extrapolated conventional MP2 with the basis sets (X=6,7). For the molecules investigated the valence MP2 correlation energy is estimated within 1 mE(h).     

Basis-set extrapolation techniques for the accurate calculation of molecular equilibrium geometries using coupled-cluster theory

To reduce remaining basis-set errors in the determination of molecular equilibrium geometries, a basis-set extrapolation (BSE) scheme is suggested for the forces used in geometry optimizations. The proposed BSE scheme is based on separating the Hartree-Fock and electron-correlation contributions and uses expressions obtained by straightforward differentiation of well established extrapolation formulas for energies when using basis sets from Dunning's hierarchy of correlation-consistent basis sets. Comparison with reference data obtained at the R12 coupled-cluster level [CCSD(T)-R12] demonstrates that BSE significantly accelerates the convergence to the basis-set limit, thus leading to improvements comparable to or even better than those obtained by increasing the cardinal number in the used basis set by one. However, BSE alone is insufficient to improve agreement with experiment, even after additional consideration of inner-shell correlation and quadruple-excitation effects (mean error and standard deviation with extrapolation are -0.014 and 0.047 pm in comparison with mean error and standard deviation of -0.002 and 0.036 pm without extrapolation). Improvement is obtained only when other contributions of similar magnitude as the BSE contributions (e.g., pentuple-excitation effects and relativistic effects) are also considered. A rather large discrepancy (of the order of a few tenths of a picometer) is observed for the F(2) molecule indicating an enhanced basis-set requirement for the various contributions in this case.     

Extrapolation of electron correlation energies to finite and complete basis set targets

The electron correlation energy of two-electron atoms is known to converge asymptotically as approximately (L+1)(-3) to the complete basis set limit, where L is the maximum angular momentum quantum number included in the basis set. Numerical evidence has established a similar asymptotic convergence approximately X(-3) with the cardinal number X of correlation-consistent basis sets cc-pVXZ for coupled cluster singles and doubles (CCSD) and second order perturbation theory (MP2) calculations of molecules. The main focus of this article is to probe for deviations from asymptotic convergence behavior for practical values of X by defining a trial function X(-beta) that for an effective exponent beta=beta(eff)(X,X+1,X+N) provides the correct energy E(X+N), when extrapolating from results for two smaller basis sets, E(X) and E(X+1). This analysis is first applied to "model" expansions available from analytical theory, and then to a large body of finite basis set results (X=D,T,Q,5,6) for 105 molecules containing H, C, N, O, and F, complemented by a smaller set of 14 molecules for which accurate complete basis set limits are available from MP2-R12 and CCSD-R12 calculations. beta(eff) is generally found to vary monotonically with the target of extrapolation, X+N, making results for large but finite basis sets a useful addition to the limited number of cases where complete basis set limits are available. Significant differences in effective convergence behavior are observed between MP2 and CCSD (valence) correlation energies, between hydrogen-rich and hydrogen-free molecules, and, for He, between partial-wave expansions and correlation-consistent basis sets. Deviations from asymptotic convergence behavior tend to get smaller as X increases, but not always monotonically, and are still quite noticeable even for X=5. Finally, correlation contributions to atomization energies (rather than total energies) exhibit a much larger variation of effective convergence behavior, and extrapolations from small basis sets are found to be particularly erratic for molecules containing several electronegative atoms. Observed effects are discussed in the light of results known from analytical theory. A carefully calibrated protocol for extrapolations to the complete basis set limit is presented, based on a single "optimal" exponent beta(opt)(X,X+1,infinity) for the entire set of molecules, and compared to similar approaches reported in the literature.     

Coupled-cluster study of isomers of H2SO2

A theoretical study has been made on six isomers of H2SO2 using coupled-cluster singles and doubles with noniterative triple excitations (CCSD(T)). The isomers studied are sulfoxylic acid (S(OH)2; C2 and Cs conformers), sulfinic acid (HS(=O)OH; 2 C1 conformers), dihydrogen sulfone (H2SO2; C2v), sulfhydryl hydroperoxide (HSOOH; C1), thiadioxirane (Cs), and dihydrogen persulfoxide (H2SOO; Cs). Molecular geometries, harmonic vibrational frequencies, and infrared intensities of all species were obtained using the CCSD(T) method and the 6-311++G(2d,2p) basis set. All aforementioned species were found to be local minima, with the exception of thiadioxirane, which has one imaginary frequency. A prior possible infrared observation of sulfinic acid was reassessed on the basis of the present data. In agreement with previous MP2 results, the present CCSD(T) data provide support for at most 4 of the 8 observed frequencies. The CCSD(T) frequencies and intensities should be of assistance in future identification of H2SO2 isomers by vibrational spectroscopy. Relative energies were calculated using the CCSD(T) method and several larger basis sets. As found previously, the lowest energy species is C2 S(OH)2, followed by Cs S(OH)2, HS(=O)OH, H2SO2, HSOOH, thiadioxirane, and H2SOO. Expanding the basis set significantly reduces the relative energies of HS(=O)OH and H2SO2. The CCSD(T) method was used with extended basis sets (up to aug-cc-pV(Q+d)Z) and basis set extrapolation in two reaction schemes to calculate the DeltaH degrees t (25 degrees C) of C2 S(OH)2. The two reaction schemes gave -285.8 and -282.7 kJ mol-1, which are quite close to a prior theoretical estimate (-290 kJ mol-1).     

Coupled-cluster studies of the lowest excited states of the 11-cis-retinal chromophore

The first few excited states of the 11-cis-retinal (PSB11) chromophore have been studied at the coupled-cluster approximative singles and doubles (CC2) level using triple-zeta quality basis sets augmented with double sets of polarisation functions. The two lowest vertical excitation energies of 2.14 and 3.21 eV are in good agreement with recently reported experimental values of 2.03 and 3.18 eV obtained in molecular beam measurements. Calculations at the time-dependent density functional theory (TDDFT) level using the B3LYP hybrid functional yield vertical excitation energies of 2.34 and 3.10 eV for the two lowest states. Zero-point vibrational energy (ZPVE) corrections of -0.09 and -0.17 eV were deduced from the harmonic vibrational frequencies for the ground and excited states calculated at the density functional theory (DFT) and TDDFT level, respectively, using the B3LYP hybrid functional.     

The sensitivity of B3LYP atomization energies to the basis set and a comparison of basis set requirements for CCSD(T) and B3LYP

The atomization energies of the 55 G2 molecules are computed using the B3LYP approach with a variety of basis sets. The 6–311 + G(3df) basis set is found to yield superior results to those obtained using the augumented-correlation-consistent valence-polarized triple-zeta set. The atomization energy of SO₂ is found to be the most sensitive to basis set and is studied in detail. Including tight d functions is found to be important for obtaining good atomization energies. The results for SO₂ are compared with those obtained using the coupled-cluster singles and doubles approach including a perturbational estimate of the triple excitations.     

Simple coupled-cluster singles and doubles method with perturbative inclusion of triples and explicitly correlated geminals: The CCSD(T)R12 model

To approach the complete basis set limit of the "gold-standard" coupled-cluster singles and doubles plus perturbative triples [CCSD(T)] method, we extend the recently proposed perturbative explicitly correlated coupled-cluster singles and doubles method, CCSD(2)(R12) [E. F. Valeev, Phys. Chem. Chem. Phys. 8, 106 (2008)], to account for the effect of connected three-electron correlations. The natural choice of the zeroth-order Hamiltonian produces a perturbation expansion with rigorously separable second-order energy corrections due to the explicitly correlated geminals and conventional triple and higher excitations. The resulting CCSD(T)(R12) energy is defined as a sum of the standard CCSD(T) energy and an amplitude-dependent geminal correction. The method is technically very simple: Its implementation requires no modification of the standard CCSD(T) program and the formal cost of the geminal correction is small. We investigate the performance of the open-shell version of the CCSD(T)(R12) method as a possible replacement of the standard complete-basis-set CCSD(T) energies in the high accuracy extrapolated ab initio thermochemistry model of Stanton et al. [J. Chem. Phys. 121, 11599 (2004)]. Correlation contributions to the heat of formation computed with the new method in an aug-cc-pCVXZ basis set have mean absolute basis set errors of 2.8 and 1.0 kJmol when X is T and Q, respectively. The corresponding errors of the standard CCSD(T) method are 9.1, 4.0, and 2.1 kJmol when X=T, Q, and 5. Simple two-point basis set extrapolations of standard CCSD(T) energies perform better than the explicitly correlated method for absolute correlation energies and atomization energies, but no such advantage found when computing heats of formation. A simple Schwenke-type two-point extrapolation of the CCSD(T)(R12)aug-cc-pCVXZ energies with X=T,Q yields the most accurate heats of formation found in this work, in error on average by 0.5 kJmol and at most by 1.7 kJmol.