Ab initio calculations of structures and interaction energies of toluene dimers including CCSD(T) level electron correlation correction

The intermolecular interaction energy of the toluene dimer has been calculated with the ARS-F model (a model chemistry for the evaluation of intermolecular interaction energy between ARomatic Systems using Feller's method), which was formerly called as the AIMI model III. The CCSD(T) (coupled cluster calculations with single and double substitutions with noniterative triple excitations) interaction energy at the basis set limit has been estimated from the second-order Moller-Plesset perturbation interaction energy at the basis set limit obtained by Feller's method and the CCSD(T) correction term obtained using a medium-size basis set. The cross (C(2)) dimer has the largest (most negative) interaction energy (-4.08 kcal/mol). The antiparallel (C(2h)) and parallel (C(S)) dimers (-3.77 and -3.41 kcal/mol, respectively) are slightly less stable. The dispersion interaction is found to be the major source of attraction in the toluene dimer. The dispersion interaction mainly determines the relative stability of the stacked three dimers. The electrostatic interaction of the stacked three dimers is repulsive. Although the T-shaped and slipped-parallel benzene dimers are nearly isoenergetic, the stacked toluene dimers are substantially more stable than the T-shaped toluene dimer (-2.62 kcal/mol). The large dispersion interaction in the stacked toluene dimers is the cause of their enhanced stability.     

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.     

Electron affinities of small uracil-water complexes: a comparison of benchmark CCSD(T) calculations with DFT

We present benchmark CCSD(T) calculations of the adiabatic electron affinities (AEA) and the vertical detachment energies (VDE) of the uracil molecule interacting with one to three water molecules. Calculations with rather large aug-cc-pVTZ basis set were only tractable when the space of virtual orbitals was reduced to about 60% of the full space employing the OVOS (Optimized Virtual Orbital Space) technique. Because of the microhydration, the valence-bound uracil anion is stabilized leading to gradually more positive values of both AEA and VDE with increasing number of participating water molecules. This agrees with experimental findings. Upon hydration by three water molecules, the electron affinity of uracil increased in comparison with AEA of the isolated molecule by about 250 up to 570 meV, depending on the geometry of the complex. CCSD(T) results confirm trends determined by DFT calculations of the microhydrated uracil and its anion, even if electron affinities of the free and hydrated uracil molecule are overestimated by DFT by up to 300 meV.     

Stacking of Benzene with Metal Chelates: Calculated CCSD(T)/CBS Interaction Energies and Potential‐Energy Curves

Accurate values for the energies of stacking interactions of nickel‐ and copper‐based six‐membered chelate rings with benzene are calculated at the CCSD(T)/CBS level. The results show that calculations made at the ωB97xD/def2‐TZVP level are in excellent agreement with CCSD(T)/CBS values. The energies of [Cu(C₃H₃O₂)(HCO₂)] and [Ni(C₃H₃O₂)(HCO₂)] chelates stacking with benzene are ?6.39 and ?4.77 kcal mol^?1, respectively. Understanding these interactions might be important for materials with properties that are dependent on stacking interactions.     

Binding free energies in the SAMPL5 octa-acid host–guest challenge calculated with DFT-D3 and CCSD(T)

We have tried to calculate the free energy for the binding of six small ligands to two variants of the octa-acid deep cavitand host in the SAMPL5 blind challenge. We employed structures minimised with dispersion-corrected density-functional theory with small basis sets and energies were calculated using large basis sets. Solvation energies were calculated with continuum methods and thermostatistical corrections were obtained from frequencies calculated at the HF-3c level. Care was taken to minimise the effects of the flexibility of the host by keeping the complexes as symmetric and similar as possible. In some calculations, the large net charge of the host was reduced by removing the propionate and benzoate groups. In addition, the effect of a restricted molecular dynamics sampling of structures was tested. Finally, we tried to improve the energies by using the DLPNO–CCSD(T) approach. Unfortunately, results of quite poor quality were obtained, with no correlation to the experimental data, systematically too positive affinities (by ~50 kJ/mol) and a mean absolute error (after removal of the systematic error) of 11–16 kJ/mol. DLPNO–CCSD(T) did not improve the results, so the accuracy is not limited by the energy function. Instead, four likely sources of errors were identified: first, the minimised structures were often incorrect, owing to the omission of explicit solvent. They could be partly improved by performing the minimisations in a continuum solvent with four water molecules around the charged groups of the ligands. Second, some ligands could bind in several different conformations, requiring sampling of reasonable structures. Third, there is an indication the continuum-solvation model has problems to accurately describe the binding of both the negatively and positively charged guest molecules. Fourth, different methods to calculate the thermostatistical corrections gave results that differed by up to 30 kJ/mol and there is an indication that HF-3c overestimates the entropy term. In conclusion, it is a challenge to calculate binding affinities for this octa-acid system with quantum–mechanical methods.     

Intermolecular interaction between hexafluorobenzene and benzene: Ab initio calculations including CCSD(T) level electron correlation correction

The intermolecular interaction energy of hexafluorobenzene-benzene has been calculated with the ARS-E model (a model chemistry for the evaluation of the intermolecular interaction energy between aromatic systems using extrapolation), which was formerly called the AIMI model. The CCSD(T) interaction energy at the basis-set limit has been estimated from the MP2 interaction energy at the basis-set limit and the CCSD(T) correction term obtained using a medium-sized basis set. The slipped-parallel (Cs) complex has the largest (most negative) interaction energy (-5.38 kcal/mol). The sandwich (C6v) complex is slightly less stable (-5.07 kcal/mol). The interaction energies of two T-shaped (C2v) complexes are very small (-1.74 and -0.88 kcal/mol). The calculated interaction energy of the slipped-parallel complex is about twice as large as that of the benzene dimer. The dispersion interaction is found to be the major source of attraction in the complex, although electrostatic interaction also contributes to the attraction. The dispersion interaction increases the relative stability of the slipped-parallel benzene dimer and the hexafluorobenzene-benzene complex compared to T-shaped ones. The electrostatic interaction is repulsive in the slipped-parallel benzene dimer, whereas it stabilizes the slipped-parallel hexafluorobenzene-benzene complex. Both electrostatic and dispersion interactions stabilize the slipped-parallel hexafluorobenzene-benzene complex, which is the cause of the preference of the slipped-parallel orientation and the larger interaction energy of the complex compared to the benzene dimer.     

Magnitude and nature of carbohydrate-aromatic interactions in fucose-phenol and fucose-indole complexes: CCSD(T) level interaction energy calculations

The CH/π contact structures of the fucose-phenol and fucose-indole complexes and the stabilization energies by formation of the complexes (E(form)) were studied by ab initio molecular orbital calculations. The three types of interactions (CH/π and OH/π interactions and OH/O hydrogen bonds) were compared and evaluated in a single molecular system and at the same level of theory. The E(form) calculated for the most stable CH/π contact structure of the fucose-phenol complex at the CCSD(T) level (-4.9 kcal/mol) is close to that for the most stable CH/π contact structure of the fucose-benzene complex (-4.5 kcal/mol). On the other hand the most stable CH/π contact structure of the fucose-indole complex has substantially larger E(form) (-6.5 kcal/mol). The dispersion interaction is the major source of the attraction in the CH/π contact structures of the fucose-phenol and fucose-indole complexes as in the case of the fucose-benzene complex. The electrostatic interactions in the CH/π contact structures are small (less than 1.5 kcal/mol). The nature of the interactions between the nonpolar surface of the carbohydrate and aromatic rings is completely different from that of the conventional hydrogen bonds where the electrostatic interaction is the major source of the attraction. The distributed multipole analysis and DFT-SATP analysis show that the dispersion interactions in the CH/π contact structure of fucose-indole complex are substantially larger than those in the CH/π contact structures of fucose-benzene and fucose-phenol complexes. The large dispersion interactions are responsible for the large E(form) for the fucose-indole complex.     

Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs

MP2 and CCSD(T) complete basis set (CBS) limit interaction energies and geometries for more than 100 DNA base pairs, amino acid pairs and model complexes are for the first time presented together. Extrapolation to the CBS limit is done by using two-point extrapolation methods and different basis sets (aug-cc-pVDZ - aug-cc-pVTZ, aug-cc-pVTZ - aug-cc-pVQZ, cc-pVTZ - cc-pVQZ) are utilized. The CCSD(T) correction term, determined as a difference between CCSD(T) and MP2 interaction energies, is evaluated with smaller basis sets (6-31G** and cc-pVDZ). Two sets of complex geometries were used, optimized or experimental ones. The JSCH-2005 benchmark set, which is now available to the chemical community, can be used for testing lower-level computational methods. For the first screening the smaller training set (S22) containing 22 model complexes can be recommended. In this case larger basis sets were used for extrapolation to the CBS limit and also CCSD(T) and counterpoise-corrected MP2 optimized geometries were sometimes adopted.     

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

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

Parameterization of a B3LYP specific correction for non-covalent interactions and basis set superposition error on a gigantic dataset of CCSD(T) quality non-covalent interaction energies

A vast number of non-covalent interaction energies at the counterpoise corrected CCSD(T) level have been collected from the literature to build a diverse new dataset. The whole dataset, which consists of 2027 CCSD(T) energies, includes most of the published data at this level. A large subset of the data was then used to train a novel, B3LYP specific, empirical correction scheme for non-covalent interactions and basis set superposition error (abbreviated as B3LYP-MM). Results obtained with our new correction scheme were directly compared to benchmark results obtained with B3LYP-D3(1) and M06-2X(2) (two popular density functions designed specifically to accurately model non-covalent interactions). For non-covalent complexes dominated by dispersion or dipole-dipole interactions all three tested methods give accurate results with the medium size aug-cc-pVDZ(3-6) basis set with MUE's of 0.27 (B3LYP-MM), 0.32 (B3LYP-D3) and 0.47 kcal/mol (M06-2X) (with explicit counterpoise corrections). These results validate both B3LYP-D3 and M06-2X for interactions of this type using a much larger data set than was presented in prior work. However, our new dispersion correction scheme shows some clear advantages for dispersion and dipole-dipole dominated complexes with the small LACVP* basis set, which is very popular in use due to its low associated computational cost: The MUE for B3LYP-MM with the LACVP* basis set for this subset of complexes (without explicit counterpoise corrections) is only 0.28 kcal/mol, compared to 0.65 kcal/mol for M06-2X or 1.16 kcal/mol for B3LYP-D3. Additionally, our new correction scheme also shows major improvements in accuracy for hydrogen-bonded systems and for systems involving ionic interactions, for example cation-π interactions. Compared to B3LYP-D3 and M06-2X, we also find that our new B3LYP-MM correction scheme gives results of higher or equal accuracy for a large dataset of conformer energies of di- and tripeptides, sugars, and cysteine.     

Explicit correlation and intermolecular interactions: investigating carbon dioxide complexes with the CCSD(T)-F12 method

We have optimized the lowest energy structures and calculated interaction energies for the CO(2)-Ar, CO(2)-N(2), CO(2)-CO, CO(2)-H(2)O, and CO(2)-NH(3) dimers with the recently developed explicitly correlated coupled cluster singles doubles and perturbative triples [CCSD(T)]-F12 methods and the associated VXZ-F12 (where X = D,T,Q) basis sets. For a given cardinal number, we find that results obtained with the CCSD(T)-F12 methods are much closer to the CCSD(T) complete basis set limit than the conventional CCSD(T) results. The relatively modest increase in the computational cost between explicit and conventional CCSD(T) is more than compensated for by the impressive accuracy of the CCSD(T)-F12 method. We recommend use of the CCSD(T)-F12 methods in combination with the VXZ-F12 basis sets for the accurate determination of equilibrium geometries and interaction energies of weakly bound electron donor acceptor complexes.     

Accurate quantum chemical energies for the interaction of hydrocarbons with oxide surfaces: CH(4)/MgO(001)

We examine the adsorption of CH(4) on the MgO(001) surface by a hybrid approach. It combines MP2 calculations with extrapolation to the complete basis set limit for the adsorption site and the CH(4)-CH(4) pair interactions in the adsorbate layer, with DFT+dispersion calculations under periodic boundary conditions for the whole system. To the total binding energy of 10.7 kJ mol(-1), the DFT+D(ispersion) correction contributes 0.7 kJ mol(-1) only, showing that the Mg(9)O(9) two-layer surface model is an excellent choice and that the interaction between the CH(4) molecules in the adsorbate layer is dominated by pair interactions. Contributions due to relaxation of the atom positions of 0.6 kJ mol(-1) (evaluated at DFT+dispersion) and of higher order correlation effects of 2.0 kJ mol(-1) (evaluated by CCSD(T)) yield a final estimate of 13.3 kJ mol(-1). To this total adsorption energy, the lateral interactions between the CH(4) molecules in the adsorbate layer contribute substantially, 4.1 kJ mol(-1)."Observed" desorption energies of 15.3 and 16.0 kJ mol(-1) have been derived from the observed Arrhenius desorption barriers (12.6 and 13.1 kJ mol(-1)) using thermal enthalpy contributions and a substantial zero-point energy (4.2 kJ mol(-1)) calculated from DFT+D vibrational frequencies. The comparison shows that our final hybrid MP2?:?PBE+D+ΔCCSD(T) estimate has reached chemical accuracy. It misses 2-3 kJ mol(-1) of binding only, which is most likely due to missing higher order correlation effects.PBE+D(ispersion) itself yields an adsorption energy that agrees within 1 kJ mol(-1) with our final hybrid MP2?:?PBE+D+ΔCCSD(T) estimate.     

Two methanes are better than one: a density functional theory study of the reactions of Mo2Oy- (y = 2-5) with methane

The mechanisms of chemical reactions of molybdenum suboxide clusters Mo(2)O(n)- (n = 2-5) with methane are investigated using B3LYP hybrid density functional theory and polarized basis sets. In particular, we focus on the reactions of the most stable structural isomers of Mo(2)O(2,3,4,5)- that lead to single molybdenum species such as HMoO(2)CH(3)-, as seen in the recent experimental study of Jarrold and co-workers. We find that, while all experimentally observed products are unfavorable due to the high amount of energy required to cleave the metal oxide, the formation of HMoO(2)CH(3)- is least endothermic. Even in this case, the thermodynamics of these reactions is very unfavorable when a single methane is reacted with the metal oxide. However, we find that the sequential addition of two methanes produces HMoO(2)CH(3)- (and another neutral molecule whose identity depends on the number of oxygens in the metal oxide) at a much lower thermodynamic cost. Further, the overall reaction barriers are much lower when the second methane adds prior to the Mo(2)O(2,3,4,5)- cleavage. The methane addition at each metal center oxidizes the metals to produce a species that is then stable enough to afford the Mo-Mo cleavage.     

Calculated conformer energies for organic molecules with multiple polar functionalities are method dependent: Taxol (case study)

Background  

Molecular mechanics (MM) and quantum chemical (QM) calculations are widely applied and powerful tools for the stereochemical and conformational investigations of molecules. The same methods have been extensively used to probe the conformational profile of Taxol (Figure 1) both in solution and at the β-tubulin protein binding site.     

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

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

True stabilization energies for the optimal planar hydrogen-bonded and stacked structures of guanine...cytosine, adenine...thymine, and their 9- and 1-methyl derivatives: complete basis set calculations at the MP2 and CCSD(T) levels and comparison with experiment

Planar H-bonded and stacked structures of guanine...cytosine (G.C), adenine...thymine (A...T), 9-methylguanine...1-methylcytosine (mG...mC), and 9-methyladenine...1-methylthymine (mA...mT) were optimized at the RI-MP2 level using the TZVPP ([5s3p2d1f/3s2p1d]) basis set. Planar H-bonded structures of G...C, mG...mC, and A...T correspond to the Watson-Crick (WC) arrangement, in contrast to mA...mT for which the Hoogsteen (H) structure is found. Stabilization energies for all structures were determined as the sum of the complete basis set limit of MP2 energies and a (DeltaE(CCSD(T)) - DeltaE(MP2)) correction term evaluated with the cc-pVDZ(0.25,0.15) basis set. The complete basis set limit of MP2 energies was determined by two-point extrapolation using the aug-cc-pVXZ basis sets for X = D and T and X = T and Q. This procedure is required since the convergency of the MP2 interaction energy for the present complexes is rather slow, and it is thus important to include the extrapolation to the complete basis set limit. For the MP2/aug-cc-pVQZ level of theory, stabilization energies for all complexes studied are already very close to the complete basis set limit. The much cheaper D-->T extrapolation provided a complete basis set limit close (by less than 0.7 kcal/mol) to the more accurate T-->Q term, and the D-->T extrapolation can be recommended for evaluation of complete basis set limits of more extended complexes (e.g. larger motifs of DNA). The convergency of the (DeltaE(CCSD(T)) - DeltaE(MP2)) term is known to be faster than that of the MP2 or CCSD(T) correlation energy itself, and the cc-pVDZ(0.25,0.15) basis set provides reasonable values for planar H-bonded as well as stacked structures. Inclusion of the CCSD(T) correction is essential for obtaining reliable relative values for planar H-bonding and stacking interactions; neglecting the CCSD(T) correction results in very considerable errors between 2.5 and 3.4 kcal/mol. Final stabilization energies (kcal/mol) for the base pairs studied are very substantial (A...T WC, 15.4; mA...mT H, 16.3; A...T stacked, 11.6; mA...mT stacked, 13.1; G...C WC, 28.8; mG...mC WC, 28.5; G...C stacked, 16.9; mG...mC stacked, 18.0), much larger than published previously. On the basis of comparison with experimental data, we conclude that our values represent the lower boundary of the true stabilization energies. On the basis of error analysis, we expect the present H-bonding energies to be fairly close to the true values, while stacked energies are still expected to be about 10% too low. The stacking energy for the mG...mC pair is considerably lower than the respective H-bonding energy, but it is larger than the mA...mT H-bonding energy. This conclusion could significantly change the present view on the importance of specific H-bonding interactions and nonspecific stacking interactions in nature, for instance, in DNA. Present stabilization energies for H-bonding and stacking energies represent the most accurate and reliable values and can be considered as new reference data.     

Activity coefficients in concentrated electrolytes: a comparison of the Brunauer-Emmett-Teller (BET) model with experimental values

Electrolyte mean ionic activity coefficients, γ_±, are calculated from the Stokes-Robinson application of the Brunauer-Emmett-Teller (BET) adsorption model for seven electrolytes (NaOH, HCl, KOH, CaCl₂, LiCl, LiBr and Ca(NO₂)₃). Only two model parameters are needed, which are derived from water vapor pressure measurements. Results were compared with experimental mean ionic activity coefficients (γ_±). Because the standard state for the BET model is the anhydrous electrolyte rather than the infinitely dilute solution, it is necessary to adjust for the differing standard states by comparing the BET and experimental γ_± values at one ‘anchoring’ concentration. The higher this ‘anchoring’ concentration, the better is the agreement between the BET and experimental mean ionic activity coefficients over the entire concentration (molality, mol kg⁻¹) range except at the two extremes of nearly pure water and nearly pure electrolyte. This is because the BET model lacks a finite limit at infinite dilution, and the experimental data are referenced to the infinite dilute solution. The BET model is well-behaved and in good agreement with experimental data on γ_±.     

High-level ab initio calculations for the four low-lying families of minima of (H2O)20. I. Estimates of MP2/CBS binding energies and comparison with empirical potentials

We report estimates of complete basis set (CBS) limits at the second-order M?ller-Plesset perturbation level of theory (MP2) for the binding energies of the lowest-lying isomers within each of the four major families of minima of (H(2)O)(20). These were obtained by performing MP2 calculations with the family of correlation-consistent basis sets up to quadruple zeta quality, augmented with additional diffuse functions (aug-cc-pVnZ, n=D, T, Q). The MP2/CPS estimates are -200.1 (dodecahedron, 30 hydrogen bonds), -212.6 (fused cubes, 36 hydrogen bonds), -215.0 (face-sharing pentagonal prisms, 35 hydrogen bonds), and -217.9 kcal/mol (edge-sharing pentagonal prisms, 34 hydrogen bonds). The energetic ordering of the various (H(2)O)(20) isomers does not follow monotonically the number of hydrogen bonds as in the case of smaller clusters such as the different isomers of the water hexamer. The dodecahedron lies ca. 18 kcal/mol higher in energy than the most stable edge-sharing pentagonal prism isomer. The TIP4P, ASP-W4, TTM2-R, AMOEBA, and TTM2-F empirical potentials also predict the energetic stabilization of the edge-sharing pentagonal prisms with respect to the dodecahedron, albeit they universally underestimate the cluster binding energies with respect to the MP2/CBS result. Among them, the TTM2-F potential was found to predict the absolute cluster binding energies to within <1% from the corresponding MP2/CBS values, whereas the error for the rest of the potentials considered in this study ranges from 3% to 5%.     

A new non-symmetric N(OH)3 species: Comparison with the C3 species and thermochemistry at the HF,DFT, MP2, MP4 and CCSD(T) levels of theory

A new non-symmetric N(OH)₃ species more stable than the C₃ structure is found at Density Functional Theory (B3PW91, B3LYP), MP2, MP4 and CCSD(T) levels with extended basis sets. C₁ and C₃ structures are qualitatively different from those of the As(OH)₃ molecule. Energy differences and interconversion barriers become smaller with increasing inclusion of electronic correlation. However at the MP2, MP4 and CCSD(T) levels, these differences increase with basis set quality. ZPE corrections lead to barrier collapse but only at the CCSD(T)/AVTZ level. The C₁ and C₃ MP2/AVTZ infrared spectra are given for future studies.