The performance of the rapid estimation of basis set error and correlation energy from partial charges method on new molecules of the G3/99 test set

Experimental enthalpies of formation have been approximated using single-point Hartree–Fock (HF)–self-consistent-field (SCF) total energies plus the rapid estimation of basis set error and correlation energy from partial charges (REBECEP) energy corrections. The energy corrections are calculated from the HF–SCF partial atomic charges and optimized atomic energy parameters. The performance of the method was tested on 51 closed-shell neutral molecules (50 molecules from the G3/99 thermochemistry database plus urea, composed of H, C, N, O, and F atoms). The predictive force of the method is demonstrated, because these larger molecules were not used for the optimization of the atomic parameters. We used the earlier RECEP-3 [HF/6-311+G(2d,p)] and REBECEP [HF/6-31G(d)] atomic parameter sets obtained from the G2/97 thermochemistry database (containing small molecules) together with natural population analysis and Mulliken partial charges. The best results were obtained using the natural population analysis charges, although the Mulliken charges also provide useful results. The root-mean-square deviations from the experimental enthalpies of formation for the selected 51 molecules are 1.15, 3.96, and 2.92 kcal/mol for Gaussian-3, B3LYP/6-11+G(3df,2p), and REBECEP (natural population analysis) enthalpies of formation, respectively (the corresponding average absolute deviations are 0.94, 7.09, and 2.27 kcal/mol, respectively). The REBECEP method performs considerably better for the 51 test molecules with a moderate 6-31G(d) basis set than the B3LYP method with a large 6-311+G(3df,2p) basis set. Received: 10 March 2001 / Accepted: 5 July 2001 / Published online: 11 October 2001     

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.     

Simple energy corrections for precise atomization energies of CHON molecules

We test a few ways to improve atomization energies of CHON molecules and found that the best way to do it is simply by correcting the atom energies of the participating atoms. Extraordinary improvement on the average errors is obtained. For the HF/6-31G** level of theory an average error of 271.2 kcal/mol on 115 molecules is improved to 6.7 kcal/mol simply by correcting the atomic energy of the four chon atoms. The corrections to density functional methods allow us to reach a chemical accuracy of 2 kcal/mol.     

Rapid Estimation of Basis Set Error and Correlation Energy Based on Mulliken Charges and Mulliken Matrix with the Small 6-31g* Basis Set

Good, density functional quality (B3LYP/6-31G*) ground state total electronic energies have been approximated using single point Hartree–Fock-self consistent field (HF-SCF/6-31G*) total energies and Mulliken partial charges versus. Mulliken matrix (electrons assigned to atoms and atoms pairs from Mulliken population analysis). This is a development of our rapid estimation of basis set error and correlation energy from partial charges (REBECEP) method, published earlier (see references [21,22,30]. The development is as follows: (1) A larger set of atoms (H, C, N, O, F, Si, P and S) are considered as building blocks for closed shell, neutral, ground state molecules at their equlibrium geometry; (2) geometries near equilibrium geometry are also considered; (3) A larger set, containing 115 molecules, was used to fit REBECEP parameters; (4) most importantly, electrons belonging to chemical bonds (between atom pairs) are also considered (Mulliken matrix) in addition to the atoms (Mulliken charges), using more REBECEP parameters to fit and yielding a more flexible algorithm. With these parameters a rather accurate closed shell ground state electronic total energy can be obtained from a small basis set HF-SCF calculation in the vicinity of optimal geometry. The 3.3 kcal/mol root mean square deviation of REBECEP improves to 1.5 kcal/mol when using Mulliken matrix instead of Mulliken charges.     

Multiple bonds and excited states from the Hartree-Fock-Heitler-London method

The recently proposed Hartree-Fock-Heitler-London, HF-HL, method (Corongiu, G. J. Phys. Chem. A 2006, 110, 11584) previously tested for single bond molecules is validated by potential energy computations for open and closed shells, single and multiple bonds, in ground and excited states of homopolar diatomic molecules of the first and second period. The simple HF-HL function, including the configurations for 2s/2p near degeneracy and avoiding state crossing, yields correct dissociation products, qualitatively correct binding, and accounts for non-dynamical correlation. Addition of ionic structures improves the ab initio HF-HL function and yields about 95% of the experimental binding energies on average. Computed excitation energies are also in agreement with laboratory values as verified for the 3 Pi u and 3 Zeta g- excited states of the C2 molecule. Computation of the remaining dynamical correlation using a semiempirical functional yields binding energies with an average deviation of 1.5 kcal/mol from laboratory values, and total energies with an average deviation of 0.7 kcal/mol from exact nonrelativistic dissociation energies.     

Accurate thermochemistry from corrected Hartree–Fock results: rapid estimation of nearly experimental quality total energy using the small 6-31G(d) basis set

Gaussian-3 ground-state total electronic energies have been approximated using single point 6-31G(d) basis set Harteee–Fock self-consistent-field (HF-SCF) total energies and partial charges based on our earlier rapid estimation of correlation energy from partial charges method. Sixty-five closed-shell neutral molecules (composed of H, C, N, O, and F atoms) of the G2/97 thermochemistry database were selected for the present study. The main feature in this work is that the␣basis set error has been treated by the least squares fit of rapid estimation of basis set error and correlation energy from partial charges (REBECEP) atomic parameters. With these parameters a rather accurate closed-shell ground-state electronic total energy can be obtained from a small basis set HF-SCF calculation in the vicinity of stationary points. The average absolute deviation of the best REBECEP enthalpies of formation from the experimental enthalpies of formation is 1.39 kcal/mol for the test set of 65 enthalpies of neutral molecules. Received: 11 December 2000 / Accepted: 6 February 2001/Published online: 11 October 2001     

Fast,efficient generation of high-quality atomic charges. AM1-BCC model: II. Parameterization and validation

We present the first global parameterization and validation of a novel charge model, called AM1-BCC, which quickly and efficiently generates high-quality atomic charges for computer simulations of organic molecules in polar media. The goal of the charge model is to produce atomic charges that emulate the HF/6-31G* electrostatic potential (ESP) of a molecule. Underlying electronic structure features, including formal charge and electron delocalization, are first captured by AM1 population charges; simple additive bond charge corrections (BCCs) are then applied to these AM1 atomic charges to produce the AM1-BCC charges. The parameterization of BCCs was carried out by fitting to the HF/6-31G* ESP of a training set of >2700 molecules. Most organic functional groups and their combinations were sampled, as well as an extensive variety of cyclic and fused bicyclic heteroaryl systems. The resulting BCC parameters allow the AM1-BCC charging scheme to handle virtually all types of organic compounds listed in The Merck Index and the NCI Database. Validation of the model was done through comparisons of hydrogen-bonded dimer energies and relative free energies of solvation using AM1-BCC charges in conjunction with the 1994 Cornell et al. forcefield for AMBER.(13) Homo- and hetero-dimer hydrogen-bond energies of a diverse set of organic molecules were reproduced to within 0.95 kcal/mol RMS deviation from the ab initio values, and for DNA dimers the energies were within 0.9 kcal/mol RMS deviation from ab initio values. The calculated relative free energies of solvation for a diverse set of monofunctional isosteres were reproduced to within 0.69 kcal/mol of experiment. In all these validation tests, AMBER with the AM1-BCC charge model maintained a correlation coefficient above 0.96. Thus, the parameters presented here for use with the AM1-BCC method present a fast, accurate, and robust alternative to HF/6-31G* ESP-fit charges for general use with the AMBER force field in computer simulations involving organic small molecules.     

A new force field (ECEPP-05) for peptides, proteins, and organic molecules

Parametrization and testing of a new all-atom force field for organic molecules and peptides with fixed bond lengths and bond angles are described. The van der Waals parameters for both the organic molecules and the peptides were taken from J. Phys. Chem. B 2003, 107, 7143 and J. Phys. Chem. B 2004, 108, 12181. First, the values of the 1-4 nonbonded and electrostatic scale factors appropriate to the new force field were determined by computing the conformational energies of six model molecules, namely, ethanol, ethylamine, propanol, propylamine, 1,2-ethanediol, and 1,3-propanediol with different values of these factors. The partial atomic charges of these molecules were obtained by fitting to the electrostatic potentials calculated with the HF/6-31G quantum-mechanical method. Two different charge models (single- and multiple-conformation-derived) were also considered. We demonstrated that the charge model has a stronger effect on the conformational energies than the 1-4 scaling. The choice of a charge model affected the conformational energies of even the smallest molecules considered, whereas the effect of the 1-4 electrostatic or nonbonded scaling was apparent only for 1,3-propanediol. The best agreement with high-level ab initio data was obtained with the multiple-conformation-derived charges and with no scaling of the 1-4 nonbonded or electrostatic interactions (scale factors of 1.0). Next, the torsional parameters of a large number of neutral and charged organic molecules, assumed to be models of the side chains of the 20 naturally occurring amino acids, were computed by fitting to rotational energy profiles obtained from ab initio MP2/6-31G calculations. The quality of the fits was high with average errors for torsional profiles of less than 0.2 kcal/mol. To derive the torsional parameters for the peptide backbone, the partial atomic charges of the 20 neutral and charged amino acids were obtained by fitting to the electrostatic potentials of terminally blocked amino acids using the HF/6-31G quantum-mechanical method. Then, the phi-psi energy maps of Ac-Ala-NMe and Ac-Gly-NMe were computed using MP2/6-31G//HF/6-31G quantum-mechanical methods. The phi-psi energy map of Ac-Ala-NMe was used for refinement of the nonbonded parameters for the backbone nitrogen and hydrogen bonded to it. Subsequently, the main-chain torsional parameters were obtained by fitting the molecular mechanics energies to the phi-psi energy maps of Ac-Ala-NMe and Ac-Gly-NMe. The transferability of the entire force field was demonstrated by reproducing the main energy minima of terminally blocked Ala3 from the literature. The performance of the force field was also evaluated by simulating crystal structures of small peptides. By comparison of simulated and experimental data, examination of the torsional-angle and atom-positional root-mean-square deviations of the energy-minimized crystal structures from the corresponding X-ray model structures demonstrated high accuracy of the force field.     

ESCAPE: A novel approach for a fast estimation of dynamic correlation energies: Application to large organic molecules

Advanced wave function-based quantum chemical ab initio methods, such as CCSD(T), are able to calculate the energies of small- to medium-sized molecules with chemical accuracy. Unfortunately, these methods scale quite unfavorably with the size of the system and are getting too time consuming—and too expensive—for larger molecules. In order to be able to treat larger organic molecules, we propose a novel scheme for a quick and reliable estimate of molecular correlation energies, which we call ESCAPE (ES timation of C orrelA tion energies by P air E nergies). It is based on the pair correlation energies for localized molecular orbitals that have been generated by CCSD[T] and fitted to suitable functional forms. All fit parameters are stored in a large parameter file. Aiming at chemical accuracy (±1 kcal/mol), we have first limited our approach to aliphatic hydrocarbons. The total molecular CCSD[T] correlation energies of a training set of 41 aliphatic hydrocarbons could be reproduced with a mean absolute error (MAE) of 0.56 kcal/mol or 0.11%. A similar accuracy could be obtained for a test set of 11 additional hydrocarbons with up to eight carbon atoms (MAE of 0.65 kcal/mol or 0.09%). In a more critical test, we checked the small energy differences for a set of 13 isomerization reactions. The comparison with experimental data showed that we could reach chemical accuracy as well. Our estimate (MAE of 0.55 kcal/mol) is slightly inferior to the CCSD[T] result (MAE of 0.17 kcal/mol), but superior to SCF, DFT/B3LYP, and DFT/B3LYP + D3. Moreover, in all cases, we obtained the correct sign, that is, the correct equilibrium structure. A similar accuracy could be reached in an application to the three lowest isomers of the C₆₀ molecule. Using the example of a set of eight alcohols, we were able to proof the method's ability for molecules including heteroatoms. Three fast steps are necessary for the application to any aliphatic hydrocarbon or alcohol: (1) An SCF calculation at the selected molecular geometry; it can be fast since a medium size basis set is generally sufficient. (2) The localization of the occupied molecular orbitals and determination of their properties (center of charge and spatial extent). (3) Estimate of the correlation energy using the existing parameter file. © 2019 Wiley Periodicals, Inc.     

Accurate calculation of the heats of formation for large main group compounds with spin-component scaled MP2 methods

Three MP2-type electron correlation treatments and standard density functional theory (DFT) approaches are used to predict the heats of formation for a wide variety of different molecules. The SCF and MP2 calculations are performed efficiently using the resolution-of-the-identity (RI) approximation such that large basis set (i.e., polarized valence quadruple-zeta quality) treatments become routinely possible for systems with 50-100 atoms. An atom equivalent scheme that corrects the calculated atomic energies is applied to extract the "real" accuracy of the methods for chemically relevant problems. It is found that the spin-component-scaled MP2 method (SCS-MP2, J. Chem. Phys, 2003, 118, 9095) performs best and provides chemical accuracy (MAD of 1.18 kcal/mol) for a G2/97 test set of molecules. The computationally more economical SOS-MP2 variant, which retains only the opposite-spin part of the correlation energy, is slightly less accurate (MAD of 1.36 kcal/mol) than SCS-MP2. Both spin-component-scaled MP2 treatments perform significantly better than standard MP2 (MAD of 1.77 kcal/mol) and DFT-B3LYP (MAD of 2.12 kcal/mol). These conclusions are supported by results obtained for a second test set of complex systems containing 70 molecules, including charged, strained, polyhalogenated, hypervalent, and large unsaturated species (e.g. C60). For this set, DFT-B3LYP performs badly (MAD of 8.6 kcal/mol) with many errors >10-20 kcal/mol while the spin-component-scaled MP2 methods are still very accurate (MAD of 2.8 and 3.7 kcal/mol, respectively). DFT-B3LYP shows an obvious tendency to underestimate molecular stability as the system size increases. Out of six density functionals tested, the hybrid functional PBE0 performs best. All in all, the SCS-MP2 method, together with large AO basis sets, clearly outperforms current DFT approaches and seems to be the most accurate quantum chemical model that routinely can predict the thermodynamic properties of large main group compounds.     

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.     

The correlation consistent composite approach (ccCA): an alternative to the Gaussian-n methods

An alternative to the Gaussian-n (G1, G2, and G3) composite methods of computing molecular energies is proposed and is named the "correlation consistent composite approach" (ccCA, ccCA-CBS-1, ccCA-CBS-2). This approach uses the correlation consistent polarized valence (cc-pVXZ) basis sets. The G2-1 test set of 48 enthalpies of formation (DeltaHf), 38 adiabatic ionization potentials (IPs), 25 adiabatic electron affinities (EAs), and 8 adiabatic proton affinities (PAs) are computed using this approach, as well as the DeltaHf values of 30 more systems. Equilibrium molecular geometries and vibrational frequencies are obtained using B3LYP density functional theory. When applying the ccCA-CBS method with the cc-pVXZ series of basis sets augmented with diffuse functions, mean absolute deviations within the G2-1 test set compared to experiment are 1.33 kcal mol(-1) for DeltaHf,0.81 kcal mol(-1) for IPs, 1.02 kcal mol(-1) for EAs, and 1.51 kcal mol(-1) for PAs, without including the "high-level correction" (HLC) contained in the original Gn methods. Whereas the HLC originated in the Gaussian-1 method as an isogyric correction, it evolved into a fitted parameter that minimized the error of the composite methods, eliminating its physical meaning. Recomputing the G1 and G3 enthalpies of formation without the HLC reveals a systematic trend where most DeltaHf values are significantly higher than experimental values. By extrapolating electronic energies to the complete basis set (CBS) limit and adding G3-like corrections for the core-valence and infinite-order electron correlation effects, ccCA-CBS-2 often underestimates the experimental DeltaHf, especially for larger systems. This is desired as inclusion of relativistic and atomic spin-orbit effects subsequently improves theoretical DeltaHf values to give a 0.81 kcal mol(-1) mean absolute deviation with ccCA-CBS-2. The ccCA-CBS method is a viable "black box" method that can be used on systems with at least 10-15 heavy atoms.     

The limitations of Slater's element-dependent exchange functional from analytic density-functional theory

Our recent formulation of the analytic and variational Slater-Roothaan (SR) method, which uses Gaussian basis sets to variationally express the molecular orbitals, electron density, and the one-body effective potential of density-functional theory, is reviewed. Variational fitting can be extended to the resolution of identity method, where variationality then refers to the error in each two-electron integral and not to the total energy. However, a Taylor-series analysis shows that all analytic ab initio energies calculated with variational fits to two-electron integrals are stationary. It is proposed that the appropriate fitting functions be charge neutral and that all ab initio energies be evaluated using two-center fits of the two-electron integrals. The SR method has its root in Slater's Xalpha method and permits an arbitrary scaling of the Slater-Gàspàr-Kohn-Sham exchange-correlation potential around each atom in the system. The scaling factors are Slater's exchange parameters alpha. Of several ways of choosing these parameters, two most obvious are the Hartree-Fock (HF) alpha(HF) values and the exact atomic alpha(EA) values. The former are obtained by equating the self-consistent Xalpha energy and the HF energies, while the latter set reproduces exact atomic energies. In this work, we examine the performance of the SR method for predicting atomization energies, bond distances, and ionization potentials using the two sets of alpha parameters. The atomization energies are calculated for the extended G2 set of 148 molecules for different basis-set combinations. The mean error (ME) and mean absolute error (MAE) in atomization energies are about 25 and 33 kcal/mol, respectively, for the exact atomic alpha(EA) values. The HF values of exchange parameters alpha(HF) give somewhat better performance for the atomization energies with ME and MAE being about 15 and 26 kcal/mol, respectively. While both sets give performance better than the local-density approximation or the HF theory, the errors in atomization energy are larger than the target chemical accuracy. To further improve the performance of the SR method for atomization energies, a new set of alpha values is determined by minimizing the MAE in atomization energies of 148 molecules. This new set gives atomization energies half as large (MAE approximately 14.5 kcal/mol) and that are slightly better than those obtained by one of the most widely used generalized-gradient approximations. Further improvements in atomization energies require going beyond Slater's functional form for exchange employed in this work to allow exchange-correlation interactions between electrons of different spins. The MAE in ionization potentials of 49 atoms and molecules is about 0.5 eV and that in bond distances of 27 molecules is about 0.02 A. The overall good performance of the computationally efficient SR method using any reasonable set of alpha values makes it a promising method for study of large systems.     

The quest for best suited references for configuration interaction singles calculations of core excited states

Near edge X‐ray absorption fine structure (NEXAFS) simulations based on the conventional configuration interaction singles (CIS) lead to excitation energies, which are systematically blue shifted. Using a (restricted) open shell core hole reference instead of the Hartree Fock (HF) ground state orbitals improves (Decleva et al., Chem. Phys., 1992, 168, 51) excitation energies and the shape of the spectra significantly. In this work, we systematically vary the underlying SCF approaches, that is, based on HF or density functional theory, to identify best suited reference orbitals using a series of small test molecules. We compare the energies of the K edges and NEXAFS spectra to experimental data. The main improvement compared to conventional CIS, that is, using HF ground state orbitals, is due to the electrostatic influence of the core hole. Different SCF approaches, density functionals, or the use of fractional occupations lead only to comparably small changes. Furthermore, to account for bigger systems, we adapt the core‐valence separation for our approach. We demonstrate that the good quality of the spectrum is not influenced by this approximation when used together with the non‐separated ground state wave function. Simultaneously, the computational demands are reduced remarkably. © 2016 Wiley Periodicals, Inc.     

Investigation of the use of B3LYP zero-point energies and geometries in the calculation of enthalpies of formation

The use of B3LYP/6–31G^* zero-point energies and geometries in the calculation of enthalpies of formation has been investigated for the enlarged G2 test set of 148 molecules [J. Chem. Phys. 106 (1997) 1063]. A scale factor of 0.96 for the B3LYP zero-point energies gives an average absolute deviation nearly the same as scaled HF/6–31G^* zero-point energies for G2, G2(MP2), and B3LYP/6–311 + G(3df,2p) enthalpies. A scale factor of 0.98, which has been recommended in some studies, increases the average absolute deviation by about 0.2 kcal/mol. Geometries from B3LYP/6–31G^* are found to do as well as MP2/6–31G^* geometries in the calculation of the enthalpies of formation.     

Gaussian-4 theory

The Gaussian-4 theory (G4 theory) for the calculation of energies of compounds containing first- (Li-F), second- (Na-Cl), and third-row main group (K, Ca, and Ga-Kr) atoms is presented. This theoretical procedure is the fourth in the Gaussian-n series of quantum chemical methods based on a sequence of single point energy calculations. The G4 theory modifies the Gaussian-3 (G3) theory in five ways. First, an extrapolation procedure is used to obtain the Hartree-Fock limit for inclusion in the total energy calculation. Second, the d-polarization sets are increased to 3d on the first-row atoms and to 4d on the second-row atoms, with reoptimization of the exponents for the latter. Third, the QCISD(T) method is replaced by the CCSD(T) method for the highest level of correlation treatment. Fourth, optimized geometries and zero-point energies are obtained with the B3LYP density functional. Fifth, two new higher level corrections are added to account for deficiencies in the energy calculations. The new method is assessed on the 454 experimental energies in the G305 test set [L. A. Curtiss, P. C. Redfern, and K. Raghavachari, J. Chem. Phys. 123, 124107 (2005)], and the average absolute deviation from experiment shows significant improvement from 1.13 kcal/mol (G3 theory) to 0.83 kcal/mol (G4 theory). The largest improvement is found for 79 nonhydrogen systems (2.10 kcal/mol for G3 versus 1.13 kcal/mol for G4). The contributions of the new features to this improvement are analyzed and the performance on different types of energies is discussed.     

Accurate thermochemistry and spectroscopy of the oxygen-protonated sulfur dioxide isomers

Despite the promising relevance of protonated sulfur dioxide in astrophysical and atmospheric fields, its thermochemical and spectroscopic characterization is very limited. High-level quantum-chemical calculations have shown that the most stable isomer is the cis oxygen-protonated sulfur dioxide, HOSO(+), while the trans form is about 2 kcal mol(-1) less stable; even less stable (by about 42 kcal mol(-1)) is the S-protonated isomer [V. Lattanzi et al., J. Chem. Phys., 2010, 133, 194305]. The enthalpy of formation for the cis- and trans-HOSO(+) is presented, based on the well tested HEAT protocol [A. Tajti et al., J. Chem. Phys., 2004, 121, 11599]. Systematically extrapolated ab initio energies, accounting for electron correlation through coupled cluster theory, including up to single, double, triple and quadruple excitations, have been corrected for core-electron correlation, anharmonic zero-point vibrational energy, diagonal Born-Oppenheimer and scalar relativistic effects. As a byproduct, proton affinity of sulfur dioxide and atomization energies have also been obtained at the same levels of theory. Vibrational and rotational spectroscopic properties have been investigated by means of composite schemes that allow us to account for truncation of basis set as well as core correlation. Where available, for both thermochemistry and spectroscopy, very good agreement with experimental data has been observed.     

Predicting hydration free energies of neutral compounds by a parametrization of the polarizable continuum model

A parametrization of the polarizable continuum model (PCM) is presented having the experimental hydration free energies of 215 neutral molecules as target. The cavitation and dispersion contributions were based on the Tu?on-Silla-Pascual-Ahuir (Tu?on; et al. Chem. Phys. Lett. 1993, 203, 289) and Floris-Tomasi (Floris, F.; Tomasi, J. J. Comput. Chem. 1989, 10, 616) expressions, respectively. Both the polar and nonpolar contributions were evaluated on the same solvent-excluding molecular surface that used unscaled Bondi atomic radii. The parametrization was provided for the HF, Xalpha, LSDA, B3LYP, and mPW1PW91 methods at the 6-31G(d) basis set, and the results are in fair agreement with the experimental data. For the sake of comparison, the PCM(UAHF) and our parametrization (PCM2), both at HF level, have produced DeltaG(PCM(UAHF)) = aDeltaGexp (a = 1.02 +/- 0.02, r = 0.945, sd = 0.987, Ftest = 1778) and DeltaG(PCM2) = aDeltaGexp (a = 0.95 +/- 0.02, r = 0.952, sd = 0.843, Ftest = 2070), respectively. The mean absolute deviations from experimental data were 0.67 and 0.68 kcal/mol for PCM(UAHF) and PCM2, respectively.     

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.     

Quantum monte carlo study of heats of formation and bond dissociation energies of small hydrocarbons

A quantum Monte Carlo (QMC) benchmark study of heats of formation at 298 K and bond dissociation energies (BDEs) of 22 small hydrocarbons is reported. Diffusion Monte Carlo (DMC) results, obtained using a simple product trial wavefunctions consisting of a single determinant and correlation function, are compared to experiment and to other theory including a version of complete basis set theory (CBS‐Q) and density functional theory (DFT) with the B3LYP functional. For heats of formation, the findings are a mean absolute deviation from experiment of 1.2 kcal/mol for CBS‐Q, 2.0 kcal/mol for B3LYP, and 2.2 kcal/mol for DMC. The mean absolute deviation of 31 BDEs is 2.0 kcal/mol for CBS‐Q, 4.2 kcal/mol for B3LYP, and 2.5 kcal/mol for DMC. These findings are for 17 BDEs of closed‐shell molecules that have mean absolute deviations from experiment of 1.7 kcal/mol (CBS‐Q), 4.0 kcal/mol (B3LYP), and 2.2 kcal/mol (DMC). The corresponding results for the 14 BDEs of open‐shell molecules studied are 2.4 kcal/mol (CBS‐Q), 4.3 kcal/mol (B3LYP), and 2.9 kcal/mol (DMC). The DMC results provide a baseline from which improvement using multideterminant trial functions can be measured. © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 583–592, 2005