A geometrical correction for the inter- and intra-molecular basis set superposition error in Hartree-Fock and density functional theory calculations for large systems

A semi-empirical counterpoise-type correction for basis set superposition error (BSSE) in molecular systems is presented. An atom pair-wise potential corrects for the inter- and intra-molecular BSSE in supermolecular Hartree-Fock (HF) or density functional theory (DFT) calculations. This geometrical counterpoise (gCP) denoted scheme depends only on the molecular geometry, i.e., no input from the electronic wave-function is required and hence is applicable to molecules with ten thousands of atoms. The four necessary parameters have been determined by a fit to standard Boys and Bernadi counterpoise corrections for Hobza's S66×8 set of non-covalently bound complexes (528 data points). The method's target are small basis sets (e.g., minimal, split-valence, 6-31G*), but reliable results are also obtained for larger triple-ζ sets. The intermolecular BSSE is calculated by gCP within a typical error of 10%-30% that proves sufficient in many practical applications. The approach is suggested as a quantitative correction in production work and can also be routinely applied to estimate the magnitude of the BSSE beforehand. The applicability for biomolecules as the primary target is tested for the crambin protein, where gCP removes intramolecular BSSE effectively and yields conformational energies comparable to def2-TZVP basis results. Good mutual agreement is also found with Jensen's ACP(4) scheme, estimating the intramolecular BSSE in the phenylalanine-glycine-phenylalanine tripeptide, for which also a relaxed rotational energy profile is presented. A variety of minimal and double-ζ basis sets combined with gCP and the dispersion corrections DFT-D3 and DFT-NL are successfully benchmarked on the S22 and S66 sets of non-covalent interactions. Outstanding performance with a mean absolute deviation (MAD) of 0.51 kcal/mol (0.38 kcal/mol after D3-refit) is obtained at the gCP-corrected HF-D3/(minimal basis) level for the S66 benchmark. The gCP-corrected B3LYP-D3/6-31G* model chemistry yields MAD=0.68 kcal/mol, which represents a huge improvement over plain B3LYP/6-31G* (MAD=2.3 kcal/mol). Application of gCP-corrected B97-D3 and HF-D3 on a set of large protein-ligand complexes prove the robustness of the method. Analytical gCP gradients make optimizations of large systems feasible with small basis sets, as demonstrated for the inter-ring distances of 9-helicene and most of the complexes in Hobza's S22 test set. The method is implemented in a freely available FORTRAN program obtainable from the author's website.     

MCSCF study of singlet oxygen addition to ethenol—a model of photooxidation reactions of unsaturated and aromatic compounds bearing hydroxy groups

The reaction path of singlet (¹Δ_g) oxygen addition to ethenol (vinyl alcohol)—a model of the reactions of singlet oxygen with aromatic and unsaturated compounds bearing the hydroxy groups—has been studied by means of MCSCF calculations, using various active spaces and basis sets. The effects of dynamic correlation (at the PT2 level) and basis set superposition error (BSSE) on relative energies were also investigated. It was found that including polarization functions is necessary to obtain geometries of the oxygen moiety consistent with the available experimental data. Two possible reaction products were considered: 1-hydroxy-1,2-dioxethane (peroxide) and 2-hydroperoxyethanal-1 (hydroperoxide); their energies are 24.1 and 36.6 kcal/mol (44.1 and 78.2 kcal/mol with the PT2 contribution and BSSE correction) below the dissociation limit, respectively (all energies reported here refer to the 6-31G** basis set and an active space composed of eight orbitals and ten electrons). A common stage of both reactions is the formation of a peralcoxyl intermediate with one of the oxygen atoms attached to the unsubstituted carbon atom; the energies of the respective transition state and that of the intermediate are 30.2 and 18.7 kcal/mol (15.9 and 10.3 kcal/mol with the PT2 contribution and BSSE correction) above the dissociation limit, respectively. The energy of this transition state is the dominant energy barrier to the reaction. The intermediate can then undergo conversion to the dioxethane product, to the perepoxide intermediate, or via a proton transfer, directly to the hydroperoxide, the last route being the most probable one. The perepoxide intermediate, after a proton transfer, also readily gives the hydroperoxide. It was also found that the unimolecular conversion from dioxethane to hydroperoxide via a proton transfer from the hydroxy group accompanied with ring cleavage requires an activation energy of at least 56 kcal/mol, making this reaction path highly improbable. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1668–1681, 1997     

Combined bond-polarization basis sets for accurate determination of dissociation energies. II. Basis set superposition error as a function of the parent basis set

The effect of the parent basis set on the basis set superposition error caused by bond functions is investigated systematically. An important difference between BSSE at the SCF and correlated levels is pointed out. Three new basis sets are defined, denoted 6-311 + G(d,p)B, 6-311 + G(2d,p)B, and 6-311 + G(2df,p)B. BSSE for the first-row hydrides seems to increase uniformly with increasing atomic number of the central atom. Expansion of the valence part of the basis set from 6-31G to 6-311G, as well as adding f functions, has a significant effect on the BSSE. Additional BSSEs incurred by bond functions are less than or equal to 1 kcal/mol for the 6-311 + G(2df,p)B basis set. For the dissociation energies of the first-row hydride species, agreement with experiment within only a few kcal/mol can be obtained even without resorting to isogyric reaction cycles. For high-quality calculations, adding bond functions seems to have definite advantages over expanding the polarization space beyond the [2d1f] level.     

A quantum-chemical study of the LiH + LiH⋅− ⇄ LiH2− + Li⋅ reaction

Parts of the potential energy surface of the title process and related processes have been investigated at the SCF /6-31G **, SCF /6-31++G **, and MP 2/6-31++G ** levels. The investigated reaction is exothermic (?6.23 kcal/mol, MP 4/6-31++G **//MP 2/6–31++G** level, ZPE included): A linear intermediate radical anion, Li? H? Li? H??, is significantly stabilized with respect to LiH + LiH?? (?38.74 kcal/mol, the same level as above). The BSSE at MP 2/6–31++G **//MP 2/6–31++G ** amounts to 1.8 kcal/mol. The title process seems to be suitable for experimental study in molecular beams.     

Ab initio study of hydrogen-bond formation between aliphatic and phenolic hydroxy groups and selected amino acid side chains

Hydrogen bonding was studied in 24 pairs of isopropyl alcohol and phenol as one partner, and water and amino-acid mimics (methanol, acetamide, neutral and protonated imidazole, protonated methylalamine, methyl-guanidium cation, and acetate anion) as the other partner. MP2/6-31+G* and MP2/aug-cc-pvtz calculations were conducted in the gas phase and in a model continuum dielectric environment with dielectric constant of 15.0. Structures were optimized in the gas phase with both basis sets, and zero-point energies were calculated at the MP2/6-31+G* level. At the MP2/aug-cc-pvtz level, the BSSE values from the Boys-Bernardi counterpoise calculations amount to 10-20 and 5-10% of the uncorrected binding energies of the neutral and ionic complexes, respectively. The geometry distortion energy upon hydrogen-bond formation is up to 2 kcal/mol, with the exception of the most strongly bound complexes. The BSSE-corrected MP2/aug-cc-pvtz binding energy of -27.56 kcal/mol for the gas-phase acetate...phenol system has been classified as a short and strong hydrogen bond (SSHB). The CH3NH3+...isopropyl alcohol complex with binding energy of -22.54 kcal/mol approaches this classification. The complete basis set limit (CBS) for the binding energy was calculated for twelve and six complexes on the basis of standard and counterpoise-corrected geometry optimizations, respectively. The X...Y distances of the X-H...Y bridges differ by up to 0.03 A as calculated by the two methods, whereas the corresponding CBS energy values differ by up to 0.03 kcal/mol. Uncorrected MP2/aug-cc-pvtz hydrogen-bonding energies are more negative by up to 0.35 kcal/mol than the MP2/CBS values, and overestimate the CCSD(T)/CBS binding energies generally by up to 5% for the eight studied complexes in the gas phase. The uncorrected MP2/aug-cc-pvtz binding energies decreased (in absolute value) by 11-18 kcal/mol for the ionic species and by up to 5 kcal/mol for the neutral complexes when the electrostatic effect of a polarizable model environment was considered. The DeltaECCSD(T) - DeltaEMP2 corrections still remained close to their gas-phase values for four complexes with 0, +/-1 net charges. Good correlations (R2 = 0.918-0.958) for the in-environment MP2/aug-cc-pvtz and MP2/6-31+G* hydrogen-bonding energies facilitate the high-level prediction of these energies on the basis of relatively simple MP2/6-31+G* calculations.     

Rapid evaluation of the binding energies in hydrogen-bonded amide-thymine and amide-uracil dimers in gas phase

The binding energies and the equilibrium hydrogen bond distances as well as the potential energy curves of 48 hydrogen‐bonded amide–thymine and amide–uracil dimers are evaluated from the analytic potential energy function established in our lab recently. The calculation results show that the potential energy curves obtained from the analytic potential energy function are in good agreement with those obtained from MP2/6‐311+G** calculations by including the BSSE correction. For all the 48 dimers, the analytic potential energy function yields the binding energies of the MP2/6‐311+G** with BSSE correction within the error limits of 0.50 kcal/mol for 46 dimers, only two differences are larger than 0.50 kcal/mol and the largest one is only 0.60 kcal/mol. The analytic potential energy function produces the equilibrium hydrogen bond distances of the MP2/6‐311+G** with BSSE correction within the error limits of 0.050 Å for all the 48 dimers. The analytic potential energy function is further applied to four more complicated hydrogen‐bonded amide–base systems involving amino acid side chain and β‐sheet. The values of the binding energies and equilibrium hydrogen bond distances obtained from the analytic potential energy function are also in good agreement with those obtained from MP2 calculations with the BSSE correction. These results demonstrate that the analytic potential energy function can be used to evaluate the binding energies in hydrogen‐bonded amide–base dimers quickly and accurately. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Adsorption of polycyclic aromatic hydrocarbons and inversion barriers of curved conjugated systems inside the molecular cage ExCage6+

We investigated the hosting of planar and curved π systems by ExCage⁶⁺. The M06‐2X/6‐311G* method and including basis set superposition error (BSSE) corrections showed good agreement with the experimental free energy changes in solution. The mean absolute deviation (MAD) with respect to experiment was 1.3 kcal/mol. The M06‐2X/6‐31G* method exhibited a MAD of 2.1 kcal/mol, so it may be useful to investigate large systems. The good agreement between the M06‐2X/6‐311G*+BSSE results and the M06‐2X/6‐31G* ones was due to a fortuitous error cancelation between basis set incompleteness and BSSE. The interaction energies decreased linearly with the number of CC double bonds present in the PAH, but the shape of the PAH is an important factor in determining the binding strength. Finally, we studied how effective is ExCage⁶⁺ in reducing the inversion barriers of buckybowls. For corannulene, sumanene and dibenzo[a,g]corannulene they are reduced by 2.0, 2.7, and 2.0 kcal/mol, respectively. Although these values indicate an induced fit catalysis by ExCage⁶⁺, they are far from those values calculated inside bilayer graphene. Therefore, much work is necessary to replicate the reduction of inversion barriers observed for graphene.     

BSSE‐corrected geometry and harmonic and anharmonic vibrational frequencies of formamide–water and formamide–formamide dimers

The basis set superposition error (BSSE) influence in the geometry structure, interaction energies, and intermolecular harmonic and anharmonic vibrational frequencies of cyclic formamide–formamide and formamide–water dimers have been studied using different basis sets (6‐31G, 6‐31G**, 6‐31++G**, D95V, D95V**, and D95V++**). The a posteriori “counterpoise” (CP) correction scheme has been compared with the a priori “chemical Hamiltonian approach” (CHA) both at the Hartree–Fock (HF) and second‐order Møller–Plesset many‐body perturbation (MP2) levels of theory. The effect of BSSE on geometrical parameters, interaction energies, and intermolecular harmonic vibrational frequencies are discussed and compared with the existing experimental data. As expected, the BSSE‐free CP and CHA interaction energies usually show less deep minima than those obtained from the uncorrected methods at both the HF and MP2 levels. Focusing on the correlated level, the amount of BSSE in the intermolecular interaction energies is much larger than that at the HF level, and this effect is also conserved in the values of the force constants and harmonic vibrational frequencies. All these results clearly indicate the importance of the proper BSSE‐free correlation treatment with the well‐defined basis functions. At the same time, the results show a good agreement between the a priori CHA and a posteriori CP correction scheme; this agreement is remarkable in the case of large and well‐balanced basis sets. The anharmonic frequency correction values also show an important BSSE dependence, especially for hydrogen bond stretching and for low frequencies belonging to the intermolecular normal modes. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Nitroamine Dimers

Ab initio self-consistent field(SCF) and Mφller-Plesset correlation correction methods employing 6-31G^** basis set have been applied to the optimizations of nitroamine dimers.The binding energies have been corrected for the basis set superposition error (BSSE) and the zero-point energy.Theree optimized dimers have been obtained.The BSSE corrected binding energy of the most stable dimer is predicted to be -31.85kJ/mol at the MP4/6-31G^**//MP2/6-31G^** level.The energy barriers of the Walden conversion for -NH2 group are 19.7kJ/mol and 18.3kJ/mol for monomer and the most stable dimer,respectively.The molecular interaction makes the internal rotation around N1-N2 even more difficult.The thermodynamic properties of nitroamine and its dimers at different temperatures have been calculated on the basis of vibrational analyses.The change of the Gibbs free energy for the aggregation from monomer to the most stable dimer at standard pressure and 298.2 K is predicted to be 14.05kJ/mol.     

Stabilisation energy of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, I, CN) complexes: complete basis set limit calculations at MP2 and CCSD(T) levels

Stabilisation energies of stacked structures of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, CN) complexes were determined at the CCSD(T) complete basis set (CBS) limit level. These energies were constructed from MP2/CBS stabilisation energies and a CCSD(T) correction term determined with a medium basis set (6-31G**). The former energies were extrapolated using the two-point formula of Helgaker et al. from aug-cc-pVDZ and aug-cc-pVTZ Hartree-Fock energies and MP2 correlation energies. The CCSD(T) correction term is systematically repulsive. The final CCSD(T)/CBS stabilisation energies are large, considerably larger than previously calculated and increase in the series as follows: hexafluorobenzene (6.3 kcal mol(-1)), hexachlorobenzene (8.8 kcal mol(-1)), hexabromobenzene (8.1 kcal mol(-1)) and hexacyanobenzene (11.0 kcal mol(-1)). MP2/SDD** relativistic calculations performed for all complexes mentioned and also for benzene[dot dot dot]hexaiodobenzene have clearly shown that due to relativistic effects the stabilisation energy of the hexaiodobenzene complex is lower than that of hexabromobenzene complex. The decomposition of the total interaction energy to physically defined energy components was made by using the symmetry adapted perturbation treatment (SAPT). The main stabilisation contribution for all complexes investigated is due to London dispersion energy, with the induction term being smaller. Electrostatic and induction terms which are attractive are compensated by their exchange counterparts. The stacked motif in the complexes studied is very stable and might thus be valuable as a supramolecular synthon.     

Regioselectivity in cycloaddition reaction between phosphaacetylene and diazomethane: An ab initio study

The [3 + 2] cycloaddition reaction of phosphaacetylene with diazomethane was investigated by means of high level ab initio calculations. It was deduced that the aromatic diazaphosphole is formed via a nonaromatic intermediate. The regiospecificity of the reaction is thus determined by the energy difference between the two transition states that lead to the two possible regioisomeric intermediates. Of the transition states in the concerted pathways, the one leading to the regioisomer with two PC bonds ( 3 ) was found to be more stable at all the levels of theory investigated, including coupled-cluster singles doubles (CCSD)(T)/6-311 + G^*//Møller-Plessett(MP)2/6-311 + G^* (+ basis set superposition, BSSE, correction). The energy difference between the two transition states, however, is always less than 2 kcal/mol. When the free energies in the two reactions are calculated by use of the harmonic frequencies, the energy separation between the two transition structures remains practically unchanged. The free energy of activation ΔG^† was 21 kcal/mol at the CCSD(T)/6-311 + G^* level of theory and use of the MP2/6-31 + G^* frequencies. At the MP2 level, a rather stable complex is obtained in the initial phase of the reaction. However, the stability of the complexes decreases at the CCSD(T) level, and application of the BSSE correction results in unstable complexes. © 1997 by John Wiley & Sons, Inc.     

Rapid evaluation of the binding energies between peptide amide and DNA base

The binding energies and the equilibrium hydrogen bond distances as well as the potential energy curves of 20 hydrogen‐bonded amide–base dimers are evaluated from the analytic potential energy function established in our laboratory recently. The analytic potential energy function is used to calculate the N? H···N, N? H···O?C, C? H···N, and C? H···O?C dipole–dipole attractive interaction energies and C?O···O?C, N? H···H? N, and N? H···H? C dipole–dipole repulsive interaction energies in the 20 dimers composed of DNA bases adenine, guanine, cytosine, or thymine and peptide amide. The calculation results show that the potential energy curves obtained from the analytic potential energy function are in good agreement with those obtained from MP2/6‐311+G** calculations by including the basis set superposition error (BSSE) correction. For all the 20 dimers, the analytic potential energy function yields the binding energies of the MP2/6‐311+G** with BSSE correction within the error limits of 0.50 kcal/mol for 19 dimers, only one difference is larger than 0.50 kcal/mol and the difference is only 0.61 kcal/mol. The analytic potential energy function produces the equilibrium hydrogen bond distances of the MP2/6‐311+G** with BSSE correction within the error limits of 0.030 Å for all the 20 dimers. The analytic potential energy function is further applied to four more complicated DNA base‐peptide amide systems involving amino acid side chain and β‐sheet. The values of the binding energies and equilibrium hydrogen bond distances obtained from the analytic potential energy function are also in good agreement with those obtained from MP2 calculations with the BSSE correction. These results demonstrate that the analytic potential energy function can be used to evaluate the binding energies in hydrogen‐bonded peptide amide–DNA base dimers quickly and accurately. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Infinite-basis calculations of binding energies for the hydrogen bonded and stacked tetramers of formic acid and formamide and their use for validation of hybrid DFT and ab initio methods

Benchmark stabilization energies for planar H-bonded and stacked structures of formic acid tetramers and formamide tetramers were determined as the sum of the infinite basis set limit of MP2 energies and a CCSD(T) correction term evaluated with the 6-31G*(0.25) basis set. The infinite basis (IB) set limit of MP2 energies was determined by two-point extrapolation using the aug-cc-pVXZ basis sets for X = D and T and separate extrapolation of the Hartree-Fock and correlation energies with new IB parameters for augmented basis sets determined here. Final stabilization energies (kcal/mol) for the tetramer studied are in the range of 4.6 to approximately 6.7 kcal/mol and they were used as reference data to test 14 density functionals. Among the tested DFT methods, PWB6K gives the best performance with an average error equal to only 30% of the average binding energy. In contrast, the popular B3LYP functional has an average error of 85%. We recommend the PWB6K method for exploring the potential energy surfaces of organic complexes and clusters and supramolecular assemblies.     

Ab initio study of gas phase and water-assisted tautomerization of maleimide and formamide

Maleimide serves as an important starting material in the synthesis of drugs and enzyme inhibitors. In the present paper, knowing the importance of tautomerization in maleimide for its drug action, potential energy surface of maleimide is studied and its tautomerization has been discussed and compared with tautomerization of formamide. Gas phase tautomerization of maleimide requires large amount of energy (23·21 kcal/mol) in comparison to formamide (15·05 kcal/mol) at HF/6-31+G* level. Thus making the proton transfer reaction a difficult process in gas phase. Water molecule lowers the energy barrier of tautomerization thus facilitating the tautomerization of maleimide to 5-hydroxy-pyrrol-2-one. Water assisted tautomerization of maleimide requires 19·60 kcal/mol energy at HF/6-31+G* and 17·63 kcal/mol energy at B3LYP/6-31+G* level, a decrease of 3·61 and 5·96 kcal/mol over gas phase tautomerization. Whereas, tautomerization of formamide requires 14·16 and 12·84 kcal/mol energy, a decrease of 0·89 and 2·01 kcal/mol energy over gas phase tautomerization at HF/6-31+G* and B3LYP/6-31+G* level, respectively. Water-assisted tautomerization in maleimide and formamide showed that difference in energy barrier reduces to 2·83 kcal/mol from 10·41 kcal/mol (in gas phase) at B3LYP level, which resulted that maleimide readily undergoes tautomerization in water molecule.     

Pi-pi interaction in pyridine

pi-pi Interaction in pyridine dimer and trimer has been investigated in different geometries and orientations at the ab initio (HF, MP2) and DFT (B3LYP) levels of theory using various basis sets (6-31G, 6-31G, 6-311++G) and corrected for basis set superposition error (BSSE). While the HF and DFT calculations show the pyridine dimer and the trimer to be unstable with respect to the monomer, the MP2 calculations show them to be clearly stable, thus emphasizing the need to include electron correlation while determining stacking interaction in such systems. The calculated MP2/6-311++G binding energy (100% BSSE corrected) of the parallel-sandwich, antiparallel-sandwich, parallel-displaced, antiparallel-displaced, T-up and T-down geometries for pyridine dimer are 1.53, 3.05, 2.39, 3.97, 1.91, 1.47 kcal/mol, respectively. The results show the antiparallel-displaced geometry to be the most stable. The binding energies for the trimer in parallel-sandwich, antiparallel-sandwich, and antiparallel-displaced geometry are found to be 3.18, 6.14, and 8.04 kcal/mol, respectively.     

Performance of the VBSCF method for pericyclic and π bond shift reactions

The performance of the valence bond self-consistent field (VBSCF) method was investigated in this paper by predicting the activation barriers and reaction energies in pericyclic and π bond shift reactions for hydrocarbons. The benchmarking set includes 3 electrocyclic reactions, 3 sigmatropic shifts, 3 cycloadditions, 2 cycloreversions, and 7 π bond shift reactions, where the first 11 reactions are taken from Houk's standard set (J. Phys. Chem. A 2003, 107, 11445). Computational results reveal that the VB(CI) method, which performs VBSCF calculations first with full covalent structures and then includes all mono- and di-ionic structures to compute the total energy without further orbital optimization, matches the accuracy of the complete active space SCF (CASSCF) method very well. The mean absolute error values (the deviations from the CASSCF data) are 0.01 kcal/mol for the reaction energy, and 0.26 and 0.32 kcal/mol for the activation energy with the 6-31G and 6-31G(d) basis sets, respectively. © 2018 Wiley Periodicals, Inc.     

Optimal binding site of a methane molecule on the silanol covered (010) surface of silicalite-1: ONIOM calculations

The binding energies and the corresponding structures of a methane molecule on the silanol covered (010) surface of silicalite-1 have been investigated using ab initio methods. Different levels of calculations, HF/6-31G(d), MP2/6-31G(d) and ONIOM (MP2/6-31G(d):HF/6-31G(d)) including the correction of an error due to an unbalance of the basis set, known as basis set super position error (BSSE), as well as the size of the cluster representing the silicalite-1 surface, were systematically examined to validate the model used. The ONIOM method with the BSSE correction was found to be a compromise between accuracy and computer time required. The optimal binding site on the silicalite-1 surface was observed at the configuration where the methane molecule points one H atom toward the O atom of the silanol group. The corresponding binding energy is -1.71 kJ/mol. This value is significantly higher than that of -5.65 kJ/mol when the methane molecule approaches the center of the straight channel. At this configuration, the C atom of methane was observed to locate exactly at the center of the channel. This leads to the conclusion that the methane molecule will relatively seldom be adsorbed on the silanol covered (010) surface of silicalite-1. Instead, the adsorption process will take place directly at the center of the straight channel.     

Basis sets for molecular interactions. 2. Application to H3N?HF,H3N?HOH,H2O?HF, (NH3)2, and H3CH?OH2

Modifications of the standard 6-31G** basis set as recommended in the accompanying paper are found to markedly lower the basis set superposition error (BSSE) in the title complexes, in contrast to enlargement to a triple-ζ scheme or by addition of a diffuse sp shell or a second set of d-functions without prior optimization, all of which lead to BSSE increase. After appropriate correction for correlation and superposition effects, all basis sets (with the exception of the standard 6-31G** and 6-311G** with their very large BSSE) predict the cyclic geometry of NH₃ dimer to be more stable than the linear arrangement. Correlation and BSSE can shift the equilibrium intermolecular distance in H₃CH-OH₂ by up to 0.4 Å. Failure to correct for superposition error leads to a drastic exaggeration of both the SCF and MP2 components of the interaction energy in this complex. Much better estimates are furnished by our recommended basis sets with their smaller superposition errors.     

On the use of electrostatic potential derived charges in molecular mechanics force fields. The relative solvation free energy of cis- and trans-N-methyl-acetamide

We have carried out free energy perturbation calculations on the relative solvation free energy of cis- and trans-N-methyl-acetamide (NMA). Experimentally, the solvation free energy difference has been found to be near zero. Using 6-31G* ab initio electrostatic potential derived charges for both the cis and trans conformations, we calculate a solvation free energy difference of 0.1 ± 0.1 kcal/mol. Using the 6-31G* charges derived for the trans conformation for both the cis and trans models leads to a solvation free energy difference of 0.9 ± 0.1 kcal/mol, compared to the value of 2.2 kcal/mol determined for the OPLS model for trans-NMA.     

Thermodynamic and ab initio analysis of the controversial enthalpy of formation of formaldehyde.

There are two values, -26.0 and -27.7 kcal mol(-1), that are routinely reported in literature evaluations for the standard enthalpy of formation, Delta(f) H(o)(298), of formaldehyde (CH(2)=O), where error limits are less than the difference in values. In this study, we summarize the reported literature for formaldehyde enthalpy values based on evaluated measurements and on computational studies. Using experimental reaction enthalpies for a series of reactions involving formaldehyde, in conjunction with known enthalpies of formation, its enthalpy is determined to be -26.05+/-0.42 kcal mol(-1), which we believe is the most accurate enthalpy currently available. For the same reaction series, the reaction enthalpies are evaluated using six computational methods: CBS-Q, CBS-Q//B3, CBS-APNO, G2, G3, and G3B3 yield Delta(f) H(o)(298)=-25.90+/-1.17 kcal mol(-1), which is in good agreement to our experimentally derived result. Furthermore, the computational chemistry methods G3, G3MP2B3, CCSD/6-311+G(2df,p)//B3LYP/6-31G(d), CCSD(T)/6-311+G(2df,p)//B3LYP/6-31G(d), and CBS-APNO in conjunction with isodesmic and homodesmic reactions are used to determine Delta(f) H(o)(298). Results from a series of five work reactions at the higher levels of calculation are -26.30+/-0.39 kcal mol(-1) with G3, -26.45+/-0.38 kcal mol(-1) with G3MP2B3, -26.09+/-0.37 kcal mol(-1) with CBS-APNO, -26.19+/-0.48 kcal mol(-1) with CCSD, and -26.16+/-0.58 kcal mol(-1) with CCSD(T). Results from heat of atomization calculations using seven accurate ab initio methods yields an enthalpy value of -26.82+/-0.99 kcal mol(-1). The results using isodesmic reactions are found to give enthalpies more accurate than both other computational approaches and are of similar accuracy to atomization enthalpy calculations derived from computationally intensive W1 and CBS-APNO methods. Overall, our most accurate calculations provide an enthalpy of formation in the range of -26.2 to -26.7 kcal mol(-1), which is within computational error of the suggested experimental value. The relative merits of each of the three computational methods are discussed and depend upon the accuracy of experimental enthalpies of formation required in the calculations and the importance of systematic computational errors in the work reaction. Our results also calculate Delta(f) H(o)(298) for the formyl anion (HCO(-)) as 1.28+/-0.43 kcal mol(-1).