Thermochemistry of new molecular species: SBr and HSBr

Total energies, optimized geometries, and vibrational frequencies of SBr and HSBr have been evaluated at the coupled cluster level of theory with the correlation consistent basis sets. Extrapolated to the complete basis set limit and with corrections for core-valence, scalar relativistic, and spin-orbit effects, atomization energies were computed and then combined with the experimental heats of formation of the atomic species to generate very accurate heats of formation for the species SBr and HSBr. For SBr, we predict 37.45 and 36.07 kcal/mol for DeltaHf(0 K) and DeltaHf (298.15 K), respectively, in very good agreement with the inferred experimental values of 37.98 and 36.15 kcal/mol. For HSBr, the estimate turns out to be 10.38 and 8.29 kcal/mol for DeltaHf (0 K) and DeltaHf (298.15 K), respectively. Using the more recent HBrO experimental heat of formation at 298.15 K of Lock et al., [J. Phys. Chem. 100, 7972 (1996)] the inferred experimental value for HSBr is predicted to be 8.15 kcal/mol, compared with 8.65 kcal/mol derived from the data of Ruscic and Berkowitz [J. Chem. Phys. 101, 7795 (1994)]. Considering the better agreement of the result with that predicted using the experimental value of DeltaHf(298.15 K) of Lock et al., the author also supports the suggestion made by Denis [J. Phys. Chem. A. 110, 5887 (2006)] that the result of Lock et al. should be preferred over the one of Ruscic and Berkowitz. For DeltaHf(0 K), the author found 10.38 and 10.56 kcal/mol, respectively, for the theoretical and inferred experimental estimates.     

Quantum Monte Carlo calculations of the dissociation energy of the water dimer

We report diffusion quantum Monte Carlo (DMC) calculations of the equilibrium dissociation energy D(e) of the water dimer. The dissociation energy measured experimentally, D(0), can be estimated from D(e) by adding a correction for vibrational effects. Using the measured dissociation energy and the modern value of the vibrational energy Mas et al., [J. Chem. Phys. 113, 6687 (2000)] leads to D(e)=5.00+/-0.7 kcal mol(-1), although the result Curtiss et al., [J. Chem. Phys. 71, 2703 (1979)] D(e)=5.44+/-0.7 kcal mol(-1), which uses an earlier estimate of the vibrational energy, has been widely quoted. High-level coupled cluster calculations Klopper et al., [Phys. Chem. Chem. Phys. 2, 2227 (2000)] have yielded D(e)=5.02+/-0.05 kcal mol(-1). In an attempt to shed new light on this old problem, we have performed all-electron DMC calculations on the water monomer and dimer using Slater-Jastrow wave functions with both Hartree-Fock approximation (HF) and B3LYP density functional theory single-particle orbitals. We obtain equilibrium dissociation energies for the dimer of 5.02+/-0.18 kcal mol(-1) (HF orbitals) and 5.21+/-0.18 kcal mol(-1) (B3LYP orbitals), in good agreement with the coupled cluster results.     

Correlation energy functionals dependent on an effective number of electrons: charged species and equilibrium geometries

Recently proposed spin-dependent and spin-independent correlation energy functionals [Perez-Jimenez et al., J. Chem. Phys. 116, 10571 (2002)] based on an effective number of electrons N are extended to deal with charged systems. By introducing the concept of an effective atomic number Z analogous to N, the spin-dependent functional in combination with Becke's exchange [Becke, Phys. Rev. A 38, 3098 (1988)] yields a mean absolute error (MAE) of 5.4 kcal/mol for the 88 ionization potentials and 58 electron affinities included in the extended G2 set, and a MAE of 4.1 kcal/mol for the 312 data comprising the above plus the 148 enthalpies of formation of the extended G2 set and the 18 total energies of the neutral atoms H through Ar. Geometry optimizations performed on the 53 molecules of the G2-1 test set with the above combination of exchange and correlation functionals yield MAEs of 0.017 A and 1.5 degrees for the 68 bond lengths and 29 angles analyzed as compared with the experimental estimates.     

Closed-shell ring coupled cluster doubles theory with range separation applied on weak intermolecular interactions

We explore different variants of the random phase approximation to the correlation energy derived from closed-shell ring-diagram approximations to coupled cluster doubles theory. We implement these variants in range-separated density-functional theory, i.e., by combining the long-range random phase approximations with short-range density-functional approximations. We perform tests on the rare-gas dimers He(2), Ne(2), and Ar(2), and on the weakly interacting molecular complexes of the S22 set of Jurec?ka et al. [P. Jurec?ka, J. S?poner, J. C?erny?, and P. Hobza, Phys. Chem. Chem. Phys. 8, 1985 (2006)]. The two best variants correspond to the ones originally proposed by Szabo and Ostlund [A. Szabo and N. S. Ostlund, J. Chem. Phys. 67, 4351 (1977)]. With range separation, they reach mean absolute errors on the equilibrium interaction energies of the S22 set of about 0.4 kcal/mol, corresponding to mean absolute percentage errors of about 4%, with the aug-cc-pVDZ basis set.     

Rotationally resolved reactive scattering: imaging detailed Cl+C2H6 reaction dynamics

The hydrogen atom abstraction reaction of Cl (2P3/2) with ethane has been studied using the crossed molecular beam technique with dc slice imaging at collision energies from 3.2 to 10.4 kcal/mol. The products HCl (v,J) (v = 0, J = 0-5) were state-selectively detected using 2+1 resonance enhanced multiphoton ionization. The images were used to obtain the center-of-mass frame product angular distributions and translational energy release distributions. Two general features were found in all probed HCl quantum states at 6.7 kcal/mol collision energy, and these features have distinct translational energy release and angular distributions, as described for HCl (v = 0, J = 2) in a recent preliminary report [Li et al., J. Chem. Phys. 124, 011102 (2006)]. The results for HCl (v = 0, J = 2) at four collision energies were also compared to investigate the energy-dependent dynamics. We discuss the reaction in terms of a variety of models of polyatomic reaction dynamics. The dynamics of this well studied system are more complicated than can be accounted for by a single mechanism, and the results call for further theoretical and experimental investigations.     

Potential energy surface and quantum dynamics study of rovibrational states for HO(3) (X (2)A')

An analytic potential energy surface has been constructed by fitting to about 28 thousand energy points for the electronic ground-state (X (2)A') of HO(3). The energy points are calculated using a hybrid density functional HCTH and a large basis set aug-cc-pVTZ, i.e., a HCTH/aug-cc-pVTZ density functional theory (DFT) method. The DFT calculations show that the trans-HO(3) isomer is the global minimum with a potential well depth of 9.94 kcal mol(-1) with respect to the OH + O(2) asymptote. The equilibrium geometry of the cis-HO(3) conformer is located 1.08 kcal mol(-1) above that of the trans-HO(3) one with an isomerization barrier of 2.41 kcal mol(-1) from trans- to cis-HO(3). By using this surface, a rigorous quantum dynamics (QD) study has been carried out for computing the rovibrational energy levels of HO(3). The calculated results determine a dissociation energy of 6.15 kcal mol(-1), which is in excellent agreement with the experimental value of Lester et al. [J. Phys. Chem. A, 2007, 111, 4727.].     

Estimated MP2 and CCSD(T) interaction energies of n-alkane dimers at the basis set limit: comparison of the methods of Helgaker et al. and Feller

The MP2 (the second-order M?ller-Plesset calculation) and CCSD(T) (coupled cluster calculation with single and double substitutions with noniterative triple excitations) interaction energies of all-trans n-alkane dimers were calculated using Dunning's [J. Chem. Phys. 90, 1007 (1989)] correlation consistent basis sets. The estimated MP2 interaction energies of methane, ethane, and propane dimers at the basis set limit [EMP2(limit)] by the method of Helgaker et al. [J. Chem. Phys. 106, 9639 (1997)] from the MP2/aug-cc-pVXZ (X=D and T) level interaction energies are very close to those estimated from the MP2/aug-cc-pVXZ (X=T and Q) level interaction energies. The estimated EMP2(limit) values of n-butane to n-heptane dimers from the MP2/cc-pVXZ (X=D and T) level interaction energies are very close to those from the MP2/aug-cc-pVXZ (X=D and T) ones. The EMP2(limit) values estimated by Feller's [J. Chem. Phys. 96, 6104 (1992)] method from the MP2/cc-pVXZ (X=D, T, and Q) level interaction energies are close to those estimated by the method of Helgaker et al. from the MP2/cc-pVXZ (X=T and Q) ones. The estimated EMP2(limit) values by the method of Helgaker et al. using the aug-cc-pVXZ (X=D and T) are close to these values. The estimated EMP2(limit) of the methane, ethane, propane, n-butane, n-pentane, n-hexane, n-heptane, n-octane, n-nonane, and n-decane dimers by the method of Helgaker et al. are -0.48, -1.35, -2.08, -2.97, -3.92, -4.91, -5.96, -6.68, -7.75, and -8.75 kcal/mol, respectively. Effects of electron correlation beyond MP2 are not large. The estimated CCSD(T) interaction energies of the methane, ethane, propane, and n-butane dimers at the basis set limit by the method of Helgaker et al. (-0.41, -1.22, -1.87, and -2.74 kcal/mol, respectively) from the CCSD(T)/cc-pVXZ (X=D and T) level interaction energies are close to the EMP2(limit) obtained using the same basis sets. The estimated EMP2(limit) values of the ten dimers were fitted to the form m0+m1X (X is 1 for methane, 2 for ethane, etc.). The obtained m0 and m1 (0.595 and -0.926 kcal/mol) show that the interactions between long n-alkane chains are significant. Analysis of basis set effects shows that cc-pVXZ (X=T, Q, or 5), aug-cc-pVXZ (X=D, T, Q, or 5) basis set, or 6-311G** basis set augmented with diffuse polarization function is necessary for quantitative evaluation of the interaction energies between n-alkane chains.     

Thermochemistry of the hypobromous and hypochlorous acids, HOBr and HOCl

The enthalpies of formation of HOBr and HOCl have been estimated by employing coupled cluster theory in conjunction with the correlation consistent basis sets and corrections for core-valence, relativistic, and anharmonic effects. We have employed three different reactions to estimate the DeltaH(o)(f,298)(HOBr), namely, the atomization reaction and two homodesmic reactions. Our best estimation is DeltaH(o)(f,298) (HOBr) = -15.3 +/- 0.6 kcal/mol and is very likely to lie toward the more negative values. The present value is 1.4 kcal/mol lower than the widely used experimental determination of Ruscic and Berkowitz (J. Chem. Phys. 1994, 101, 7795), DeltaH(o)(f,298)(HOBr) > -13.93 +/- 0.42 kcal/mol. However, it is closer to the more recent measurement of Lock et al. (J. Phys. Chem. 1996, 100, 7972), DeltaH(o)(f,298)(HOBr) = -14.8 +/- 1 kcal/mol. In the case of HOCl we have determined DeltaH(o)(f,298)(HOCl) = -18.1 +/- 0.3 kcal/mol, just in the middle of the two experimental values proposed, -17.8 +/- 0.5 kcal/mol (JANAF), obtained from equilibrium constant measurements, and -18.36 +/- 0.03 kcal/mol (Joens, J. A. J. Phys. Chem. A 2001, 105, 11041), determined from the measurements of the Cl-OH bond energy. If our conclusions are correct, several enthalpies of formation that have been determined by experimental chemists, Orlando and Burholder (J. Phys. Chem. 1995, 99, 1143), and theoretical chemists, Lee (J. Phys. Chem. 1995, 99, 15074), need to be revised, since a larger value was used for DeltaH(o)(f,298)(HOBr). Employing the results obtained by Orlando and Burkholder for Br(2)O we propose DeltaH(o)(f,298)(Br(2)O) = 24.9 +/- 0.6 kcal/mol, and employing Lee's enthalpies of reaction we propose the following DeltaH(o)(f,298): for BrBrO, HBrO, ClOBr, ClBrO, BrClO, BrCN, BrNC, BrNO, BrON, FOBr, and FBrO, 39.5 +/- 1, 41.0 +/- 1, 22.7 +/- 1.5, 34.2 +/- 1.5, 40.9 +/- 1.5, 43.7 +/- 1.5, 80.1 +/- 1.5, 22.3 +/- 1, 46.2 +/- 1, 17.3 +/- 1.5, and 6.3 +/- 1.5 kcal/mol, respectively. We expect that this work will stimulate new experimental measurements of the thermodynamic properties of HOBr and HOCl.     

The role of the excited electronic states in the C+ + H2O reaction

The electronic excited states of the [COH2]+ system have been studied in order to establish their role in the dynamics of the C+ + H2O-->[COH]+ +H reaction, which is a prototypical ion-molecule reaction. The most relevant minima and saddle points of the lowest excited state have been determined and energy profiles for the lowest excited doublet and quartet electronic states have been computed along the fragmentation and isomerization coordinates. Also, nonadiabatic coupling strengths between the ground and the first excited state have been computed where they can be large. Our analysis suggests that the first excited state could play an important role in the generation of the formyl isomer, which has been detected in crossed beam experiments [D. M. Sonnenfroh et al., J. Chem. Phys. 83, 3985 (1985)], but could not be explained in quasiclassical trajectory computations [Y. Ishikawa et al., Chem. Phys. Lett. 370, 490 (2003); J. R. Flores, J. Chem. Phys. 125, 164309 (2006)].     

Infrared and microwave spectra of the acetylene-ammonia and carbonyl sulfide-ammonia complexes: a comparative study of a weak C-H···N hydrogen bond and an S···N bond

We report a combined high resolution infrared and microwave spectroscopic investigation of the acetylene-ammonia and carbonyl sulfide-ammonia complexes using a pulsed slit-nozzle multipass absorption spectrometer based on a quantum cascade laser and a pulsed nozzle beam Fourier transform microwave spectrometer, respectively. The ro-vibrational transitions of the acetylene-ammonia complex have been measured at 6 μm in the vicinity of the ν(4) band of ammonia for the first time. The previously reported pure rotational transitions have been extended to higher J and K values with (14)N nuclear quadrupole hyperfine components detected and analyzed. The spectral analysis reveals that acetylene binds to ammonia through a C-H···N weak hydrogen bond to form a C(3v) symmetric top, consistent with the previous microwave [Fraser et al., J. Chem. Phys., 1984, 80, 1423] and infrared spectroscopic study at 3 μm [Hilpert et al., J. Chem. Phys., 1996, 105, 6183]. A parallel study has also been carried out for the carbonyl sulfide-ammonia complex whose pure rotational and ro-vibrational spectra at 6 μm have been detected and analyzed for the first time. The spectral and the subsequent structural analyses, in conjunction with the corresponding ab initio calculation, indicate that the OCS-NH(3) complex assumes C(3v) symmetry with S pointing to N of NH(3), in contrast to the T-shaped geometries obtained for the isoelectronic N(2)O-NH(3) and CO(2)-NH(3) complexes.     

Improved efficiency of focal point conformational analysis with truncated correlation consistent basis sets

It has been suggested that the computational cost of correlated ab initio calculations could be reduced efficiently by using truncated basis sets on hydrogen atoms (Mintz et al., J Chem Phys 2004, 121, 5629). We now explore this proposal in the context of conformational analysis of small molecules, such as hydrogen peroxide, dimethyl ether, ethyl methyl ether, formic acid, methyl formate, and several small alcohols. It is found that truncated correlation consistent basis sets that lack certain higher angular momentum functions on hydrogen atoms offer accuracy similar to traditional Dunning's basis sets for conformational analysis. Combination of such basis sets with the basis set extrapolation technique to estimate Hartree-Fock and M?ller-Plesset second order energies provides composite extrapolation model chemistries that are significantly more accurate and faster than analogous single point calculations with traditional correlation consistent basis sets. Root mean square errors of best composite extrapolation model chemistries on the used set of molecules are within 0.03 kcal/mol of traditional focal point conformational energies. The applicability of composite extrapolation methods is illustrated by performing conformational analysis of tert-butanol and cyclohexanol. For comparison, conformational energies calculated with popular molecular mechanics force fields are also given.     

Three-dimensional potential energy surface of the Ar-OH(2Pi i) complex

Pure rotational transitions in the ground state for Ar-OH and Ar-OD [Y. Ohshima et al., J. Chem. Phys. 95, 7001 (1991) and Y. Endo et al., Faraday Discuss. 97, 341 (1994)], those in the excited states of the OH vibration, nu(s)=1 and 2, observed by Fourier-transform microwave spectroscopy in the present study, rotation-vibration transitions observed by infrared-ultraviolet double-resonance spectroscopy [K. M. Beck et al., Chem. Phys. Lett. 162, 203 (1989) and R. T. Bonn et al., J. Chem. Phys. 112, 4942 (2000)], and the P-level structure observed by stimulated emission pumping spectroscopy [M. T. Berry et al., Chem. Phys. Lett. 178, 301 (1991)] have been simultaneously analyzed to determine the potential energy surface of Ar-OH in the ground state. A Schrodinger equation, considering all the freedom of motions for an atom-diatom system in the Jacobi coordinate, R, theta, and r, was numerically solved to obtain energies of the rovibrational energy levels using the discrete variable representation method. A three-dimensional potential energy surface is determined by a least-squares fitting. In the analysis the potential parameters, obtained by ab initio calculations at the RCCSD(T) level of theory with a set of basis functions of aug-cc-pVTZ and midbond functions, are used as initial values. The determined intermolecular potential energy surface and its dependence on the OH monomer bond length are compared with those of an isovalent radical complex, Ar-SH.     

Reactivity enhancement of ultracold O(3P)+H2 collisions by van der Waals interactions

The role of van der Waals forces in O((3)P)+H(2)(upsilon=1,j=0) collisions is investigated theoretically at low and ultralow temperatures. Quantum scattering calculations have been performed for zero total angular momentum using the lowest London-Eyring-Polanyi-Sato double-polynomial (3)A(") potential-energy surface reported by [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)] and its recent BMS1 and BMS2 extensions developed by [Brandao et al., J. Chem. Phys. 121, 8861 (2004)] which provide a more accurate treatment of the van der Waals interaction. Our calculations show that van der Waals forces strongly influence chemical reactivity at ultracold translational energies. The presence of a zero-energy resonance for the BMS1 surface is found to enhance reactivity in the ultracold regime and shift the Wigner threshold to lower temperatures.     

Protonated carbonyl sulfide: prospects for the spectroscopic observation of the elusive HSCO+ isomer

Spurred by the apparent conflict between ab initio predictions and infrared spectroscopic evidence regarding the relative stability of isomers of protonated carbonyl sulfide, key stationary points on the isomerization surface of HOCS(+) have been examined via systematic extrapolations of ab initio energies. Electron correlation has been accounted for using second-order M?ller-Plesset perturbation theory and coupled cluster theory through triple excitations [CCSD, CCSD(T), and CCSDT] in conjunction with the correlation consistent hierarchy of basis sets, cc-pVXZ (X=D,T,Q,5,6). HSCO(+) is predicted to lie lower in energy than HOCS(+) by 4.86 kcal mol(-1), computed using the focal point extrapolation scheme of Allen and co-workers [J. Chem. Phys. 99, 4638 (1993)] with corrections for anharmonic zero-point vibrational energy, core correlation, non-Born-Oppenheimer, and scalar relativistic effects. A transition state has been located, constituting the barrier to isomerization of HSCO(+) to HOCS(+), lying 68.9 kcal mol(-1) higher in energy than HSCO(+). This is well above predicted exothermicity [DeltaH(r) (o)(0 K)=48.1 kcal mol(-1), cc-pVQZ CCSD(T)] for the reaction considered in the experiments (HSCO(+)+H(2)-->OCS+H(3) (+)). Though proton tunneling will lead to a lower effective barrier, this prediction is consistent with the lack of HSCO(+) in electrical discharges in H(2)OCS, since the relative populations of HOCS(+) and HSCO(+) will depend on the experimental details of the protonation route rather than the relative thermodynamic stability of the isomers. Anharmonic vibrational frequencies and vibrationally corrected rotational constants from cc-pVTZ CCSD(T) cubic and quartic force constants are provided, to aid in the spectroscopic observation of the energetically favorable but apparently elusive HSCO(+) isomer.     

Alchemical prediction of hydration free energies for SAMPL

Hydration free energy calculations have become important tests of force fields. Alchemical free energy calculations based on molecular dynamics simulations provide a rigorous way to calculate these free energies for a particular force field, given sufficient sampling. Here, we report results of alchemical hydration free energy calculations for the set of small molecules comprising the 2011 Statistical Assessment of Modeling of Proteins and Ligands challenge. Our calculations are largely based on the Generalized Amber Force Field with several different charge models, and we achieved RMS errors in the 1.4-2.2 kcal/mol range depending on charge model, marginally higher than what we typically observed in previous studies (Mobley et al. in J Phys Chem B 111(9):2242-2254, 2007, J Chem Theory Comput 5(2):350-358, 2009, J Phys Chem B 115:1329-1332, 2011; Nicholls et al. in J Med Chem 51:769-779, 2008; Klimovich and Mobley in J Comput Aided Mol Design 24(4):307-316, 2010). The test set consists of ethane, biphenyl, and a dibenzyl dioxin, as well as a series of chlorinated derivatives of each. We found that, for this set, using high-quality partial charges from MP2/cc-PVTZ SCRF RESP fits provided marginally improved agreement with experiment over using AM1-BCC partial charges as we have more typically done, in keeping with our recent findings (Mobley et al. in J Phys Chem B 115:1329-1332, 2011). Switching to OPLS Lennard-Jones parameters with AM1-BCC charges also improves agreement with experiment. We also find a number of chemical trends within each molecular series which we can explain, but there are also some surprises, including some that are captured by the calculations and some that are not.     

Fitting atomic correlation parameters for RECEP (rapid estimation of correlation energy from partial charges) method to estimate molecular correlation energies within chemical accuracy

The accuracy of the RECEP method [Chem Phys 1997, 224, 33 and Chem Phys Lett 1999, 307, 469] has been increased considerably by the use of fitted atomic correlation parameters. This method allows an extremely rapid, practically prompt calculation of the correlation energy of molecules after an HF‐SCF calculation. The G2 level correlation energy and HF‐SCF charge distribution of 41 closed‐shell neutral molecules (composed of H, C, N, O, and F atoms) of the G2 thermochemistry database were used to obtain the fitted RECEP atomic correlation parameters. Four different mathematical definitions of partial charges, as a multiple choice, were used to calculate the molecular correlation energies. The best results were obtained using the natural population analysis, although the other three are also recommended for use. For the 41 molecules, the G2 results were approached within a 1.8 kcal/mol standard deviation (the mean absolute difference was 1.5 kcal/mol). The RECEP atomic correlation parameters were also tested on a different, nonoverlapping set of other 24 molecules from the G2 thermochemistry database. The G2 results of these 24 molecules were approached within a 2.3 kcal/mol standard deviation (the mean absolute difference was 1.9 kcal/mol). This method is recommended to estimate total correlation energies of closed shell ground‐state neutral molecules at stationary (minimums and transition states) points on the potential surface. Extension of the work for charged molecules, radicals, and molecules containing other atoms is straightforward. Numerical example as a recipe is also provided. © 2000 John Wiley & Sons, Inc. J Comput Chem 22: 241–254, 2001     

Aqueous solvation free energies of ions and ion-water clusters based on an accurate value for the absolute aqueous solvation free energy of the proton

Thermochemical cycles that involve pKa, gas-phase acidities, aqueous solvation free energies of neutral species, and gas-phase clustering free energies have been used with the cluster pair approximation to determine the absolute aqueous solvation free energy of the proton. The best value obtained in this work is in good agreement with the value reported by Tissandier et al. (Tissandier, M. D.; Cowen, K. A.; Feng, W. Y.; Gundlach, E.; Cohen, M. J.; Earhart, A. D.; Coe, J. V. J. Phys. Chem. A 1998, 102, 7787), who applied the cluster pair approximation to a less diverse and smaller data set of ions. We agree with previous workers who advocated the value of -265.9 kcal/mol for the absolute aqueous solvation free energy of the proton. Considering the uncertainties associated with the experimental gas-phase free energies of ions that are required to use the cluster pair approximation as well as analyses of various subsets of data, we estimate an uncertainty for the absolute aqueous solvation free energy of the proton of no less than 2 kcal/mol. Using a value of -265.9 kcal/mol for the absolute aqueous solvation free energy of the proton, we expand and update our previous compilation of absolute aqueous solvation free energies; this new data set contains conventional and absolute aqueous solvation free energies for 121 unclustered ions (not including the proton) and 147 conventional and absolute aqueous solvation free energies for 51 clustered ions containing from 1 to 6 water molecules. When tested against the same set of ions that was recently used to develop the SM6 continuum solvation model, SM6 retains its previously determined high accuracy; indeed, in most cases the mean unsigned error improves when it is tested against the more accurate reference data.     

Active-space coupled-cluster methods through connected quadruple excitations

Coupled-cluster methods that include just a subset of all connected triple, quadruple, or both excitation amplitudes, according to the ansatz of and Adamowicz co-workers [Int. Rev. Phys. Chem. 12, 339 (1993); J. Chem. Phys. 99, 1875 (1993); 100, 5792 (1994)] and Piecuch et al. [J. Chem. Phys. 110, 6103 (1999)], have been implemented into parallel execution programs. They are applicable to closed- and open-shell species and they take advantage of real Abelian point-group symmetry. A symbol manipulation program has been invoked to automate the implementation. These methods have been applied to the singlet-triplet separations of five triatomic hydrides (CH2, NH2+, SiH2, PH2+, and AsH2+) with consideration of scalar relativistic effects. They have been shown to be remarkably effective with errors arising from the use of a very small subset of higher-order excitations being no more than a few tenths of 1 kcal/mol.