High-level ab initio potential energy surfaces and vibrational energies of H2CS

Six-dimensional (6D) potential energy surfaces (PESs) of H(2)CS have been generated ab initio using the recently proposed explicitly correlated (F12) singles and doubles coupled cluster method including a perturbational estimate of connected triple excitations, CCSD(T)-F12b [T. B. Adler, G. Knizia, and H.-J. Werner, J. Chem. Phys. 127, 221106 (2007)] in conjunction with F12-optimized correlation consistent basis sets. Core-electron correlation, high-order correlation, scalar relativistic, and diagonal Born-Oppenheimer terms were included as additive high-level (HL) corrections. The resulting 6D PESs were represented by analytical functions which were used in variational calculations of the vibrational term values below 5000 cm(-1). The best PESs obtained with and without the HL corrections, VQZ-F12(*HL) and VQZ-F12?, reproduce the fundamental vibrational wavenumbers with mean absolute deviations of 1.13 and 1.22 cm(-1), respectively. A detailed analysis of the effects of the HL corrections shows how the VQZ-F12 results benefit from error cancellation. The present purely ab initio PESs will be useful as starting points for empirical refinements towards an accurate "spectroscopic" PES of H(2)CS.     

The structures and energies of pentaprismane and hexaprismane - an ab initio study

Ab initio calculations are reported for pentaprismane and hexaprismane using the 3-21G and 6-31G* basis set. The calculated structures and energies are compared with what is known experimentally and with the results of previous molecular mechanics calculations.     

Theoretical studies on thiabenzene and its fused derivatives: DFT and ab initio computations

The molecular structures of thiabenzene ( 1 ), 1‐thianaphthalene ( 2 ), 2‐thianaphthalene ( 3 ), and 9‐thiaanthracene ( 4 ) are studied using HF and DFT methods with 6‐31+G* basis set. The nonplanar boat conformers of 1–4 , with 6π‐electrons in their heterocyclic ring, appear more stable than the corresponding planar conformers with 8π‐electrons in the ring. This study focuses on the stability, the ylide character, the inversion barrier energy of sulfur atoms, and the conformational flexibility of the ring in 1–4 . © 2006 Wiley Periodicals, Inc. Heteroatom Chem 17:376–381, 2006; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20214     

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.     

Furan-formic acid dimers: an ab initio and matrix isolation study

The dimers formed by formic acid (FA) and furan are investigated by ab initio methods and matrix isolation spectroscopy. Nine complexes with binding energies between -3.91 and -0.82 kcal/mol (MP2/6-311++G(d,p) + ZPE + BSSE) are identified. Another five weaker bound complexes are localized at lower level of theory only. The binding in the furan-FA dimers can be described in terms of OH...O, C=O...H, HO...H, CH...O, OH...pi, and CH...pi interactions. Therefore, the furan-FA complexes are classified in two types: (1) the dimers where the OH hydrogen of formic acid interacts with the furan molecule and (2) the dimers where the main interactions of FA with the furan molecule are via the less acidic CH hydrogen. Duning's and Pople's triple and double basis sets were used to study the dependence of the geometries and energies of the complexes from the basis set. BSSE (basis set superposition error) counterpoise corrections (CP) were included during the geometry optimizations of all dimers at the MP2/6-31G(d,p) level of theory. Matrix isolation spectroscopy allowed us to record the IR spectrum of aggregates between FA and furan. By comparison of the experimental IR spectrum with calculated IR spectra of a variety of complexes, it was possible to identify the most stable furan-FA dimer as the major product of the aggregation.     

XMVB: a program for ab initio nonorthogonal valence bond computations

An ab initio nonorthogonal valence bond program, called XMVB, is described in this article. The XMVB package uses Heitler-London-Slater-Pauling (HLSP) functions as state functions, and calculations can be performed with either all independent state functions for a molecule or preferably a few selected important state functions. Both our proposed paired-permanent-determinant approach and conventional Slater determinant expansion algorithm are implemented for the evaluation of the Hamiltonian and overlap matrix elements among VB functions. XMVB contains the capabilities of valence bond self-consistent field (VBSCF), breathing orbital valence bond (BOVB), and valence bond configuration interaction (VBCI) computations. The VB orbitals, used to construct VB functions, can be defined flexibly in the calculations depending on particular applications and focused problems, and they may be strictly localized, delocalized, or bonded-distorted (semidelocalized). The parallel version of XMVB based on MPI (Message Passing Interface) is also available.     

Enhanced aromaticity of the transition structures for the diels-alder reactions of quinodimethanes: evidence from ab initio and DFT computations

The Diels-Alder reactions of various quinodimethanes with ethylene are studied by means of ab initio molecular orbital and density functional theory (DFT) to show the effect of aromaticity on the reaction path. The calculations reveal that these reactions are both kinetically and thermodynamically much more favored than the prototype butadiene-ethylene Diels-Alder reaction due to the aromatization process in the transition state (TS) and product. A progressive aromaticity gain is noticed during the reaction, and hence the partial pi-delocalized peripheral diene ring function is coupled with the six-electron sigma,pi-delocalized cyclic unit resulting in an enhanced aromaticity of the TS. The magnetic criteria such as magnetic susceptibility exaltation and nucleus independent chemical shift provide definitive evidence for and fully support the aromatization process and the aromaticity of the TS. The extent of sigma-pi delocalization and the bond make-break at the TS are consistent with each other, and this is strongly influenced by the adjacent pi-aromatization process. Moreover, the aromaticity trends in the resulting TSs and products parallel the activation and reaction energies; the extent of aromatization increases with increasing reaction rate and exothermicity. This confirms that aromaticity is the driving factor governing cycloadditions involving quinodimethanes.     

Accurate ab initio binding energies of the benzene dimer

Accurate binding energies of the benzene dimer at the T and parallel displaced (PD) configurations were determined using the single- and double-coupled cluster method with perturbative triple correction (CCSD(T)) with correlation-consistent basis sets and an effective basis set extrapolation scheme recently devised. The difference between the estimated CCSD(T) basis set limit electronic binding energies for the T and PD shapes appears to amount to more than 0.3 kcal/mol, indicating the PD shape is a more stable configuration than the T shape for this dimer in the gas phase. This conclusion is further strengthened when a vibrational zero-point correction to the electronic binding energies of this dimer is made, which increases the difference between the two configurations to 0.4-0.5 kcal/mol. The binding energies of 2.4 and 2.8 kcal/mol for the T and PD configurations are in good accord with the previous experimental result from ionization potential measurement.     

Fully ab initio protein‐ligand interaction energies with dispersion corrected density functional theory

Dispersion corrected density functional theory (DFT‐D3) is used for fully ab initio protein‐ligand (PL) interaction energy calculation via molecular fractionation with conjugated caps (MFCC) and applied to PL complexes from the PDB comprising 3680, 1798, and 1060 atoms. Molecular fragments with n amino acids instead of one in the original MFCC approach are considered, thereby allowing for estimating the three‐body and higher many‐body terms. n > 1 is recommended both in terms of accuracy and efficiency of MFCC. For neutral protein side‐chains, the computed PL interaction energy is visibly independent of the fragment length n. The MFCC fractionation error is determined by comparison to a full‐system calculation for the 1060 atoms containing PL complex. For charged amino acid side‐chains, the variation of the MFCC result with n is increased. For these systems, using a continuum solvation model with a dielectricity constant typical for protein environments (? = 4) reduces both the variation with n and improves the stability of the DFT calculations considerably. The PL interaction energies for two typical complexes obtained ab initio for the first time are found to be rather large (?30 and ?54 kcal/mol). © 2012 Wiley Periodicals, Inc.     

An ab initio benchmark study of hydrogen bonded formamide dimers

The five singly and doubly hydrogen bonded dimers of formamide are calculated at the correlated level by using resolution of identity M?ller-Plesset second-order perturbation theory (RIMP2) and the coupled cluster with singles, doubles, and perturbative triples [CCSD(T)] method. All structures are optimized with the Dunning aug-cc-pVTZ and aug-cc-pVQZ basis sets. The binding energies are extrapolated to the complete basis set (CBS) limit by using the aug-cc-pVXZ (X = D, T, Q) basis set series. The effect of extending the basis set to aug-cc-pV5Z on the geometries and binding energies is studied for the centrosymmetric doubly N-H...O bonded dimer FA1 and the doubly C-H...O bonded dimer FA5. The MP2 CBS limits range from -5.19 kcal/mol for FA5 to -14.80 kcal/mol for the FA1 dimer. The DeltaCCSD(T) corrections to the MP2 CBS limit binding energies calculated with the 6-31+G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ basis sets are mutually consistent to within < or =0.03 kcal/mol. The DeltaCCSD(T) correction increases the binding energy of the C-H...O bonded FA5 dimer by 0.4 kcal/mol or approximately 9% over the distance range +/-0.5 Angstrom relative to the potential minimum. This implies that the ubiquitous long-range C-H...O interactions in proteins are stronger than hitherto calculated.     

NIR-FT Raman and infrared spectra and ab initio computations of glycinium oxalate

The single crystals of glycinium oxalate are grown by slow evaporation technique and vibrational spectral analysis is carried out using NIR-FT Raman and FT-IR spectra. The ab initio quantum computations are also performed at HF/6-31 G(d) level to derive the optimized geometry, atomic charges and vibrational frequencies of the glycinium oxalate molecule. Vibrational analysis indicates the presence of peculiar intermolecular C-H...O hydrogen bonding interaction producing "blue shift" of C-H stretching frequency. The vibrational spectra confirm the existence of NH3(+) in glycinium oxalate. Hydroxyl vibrations with different inter and intra molecular H-bonding are analysed, supported by computed results.     

A comment on “Accurate ab initio determination of binding energies for rare-gas dimers by basis set extrapolation”

In this work we comment on an extrapolation scheme presented by Lee in Theoretical Chemistry Accounts, which is based on an extrapolation of energy differences instead of actual energies. In particular, we show that a very similar scheme had been introduced already in 1999, and used to estimate the MP5, CCSDT and FCI complete basis set limits of He₂. Comment to the article “Accurate ab initio determination of binding energies for rare-gas dimers by basis set extrapolation” by J.S. Lee, Theor Chem Acc (2005) 113: 87–94     

Reliable electron affinities of perfluorocyclopropane and perfluorocyclobutane from convergent ab initio computations

To resolve discrepancies concerning the magnitude of the electron affinities of perfluorocyclopropane and perfluorocyclobutane, quantum chemical calculations have been carried out with the MP2 and CCSD(T) methods in conjunction with augmented correlation consistent basis sets (aug-cc-pVX Z, X = D, T, Q). Though no experimental values have been found for perfluorocyclopropane, we estimate its electron affinity to be 0.17 eV (0.00 eV without zero-point vibrational energy corrections). In addition, determination of the electron affinity of perfluorocyclobutane (0.61 and 0.44 eV with and without zero-point vibrational energy corrections, respectively) is in good agreement with experimental values reported by Miller and co-workers (0.63 +/- 0.05 eV). This study also demonstrates that the widely prescribed B3LYP/DZP++ model chemistry for computing electron affinities does not correctly describe these systems.     

The temperature dependence of radiationless transition rates from ab initio computations

The calculation of radiationless transition rates and of their temperature dependence from first principles is addressed by combining reliable electronic computations of the normal modes of the two electronic states with Kubo's generating function approach for the evaluation of the Franck-Condon weighted density of states. The whole sets of normal modes of the involved cofactors have been employed, taking into account the effects of nuclear equilibrium position displacements, of vibrational frequency changes, and of mixing of the normal modes. Application to the case of the elementary electron transfer step between bacteriopheophytin and ubiquinone cofactors of bacterial photosynthetic reaction centers yields a temperature dependence of the electron transfer rates in very good agreement with the experimental data.     

Spectroscopic investigations and ab initio computations of the dye Chromotrope 2R

Detailed analysis of the NIR FT-Raman, FT-IR and UV–visible spectra of the dye Chromotrope 2R (C2R) has been performed. The optimized geometry of the dye is theoretically computed with the HF and DFT levels using the standard 6-31G(d) and LANL2DZ basis sets. Optimized geometry and vibrational spectra indicate that the major species in the solid state are the trans form of hydrogen bonded hydrazone tautomer. The effect of H-bonding in stabilizing a particular type of structure is also discussed. The most preferred trans-configuration for its photochemical activity has been demonstrated on the basis of torsional potential energy surface (PES) scan studies. The optimized geometries and calculated vibrational wavenumbers are evaluated via comparison with experimental values. Electronic spectra are in accordance with the nature of the electronic transitions predicted by time-dependent B3LYP/DZ calculations.     

Accurate ab initio binding energies of alkaline earth metal clusters

The effects of basis set superposition error (BSSE) and core-correlation on the electronic binding energies of alkaline earth metal clusters Y(n) (Y = Be, Mg, Ca; n = 2-4) at the Moller-Plesset second-order perturbation theory (MP2) and the single and double coupled cluster method with perturbative triples correction (CCSD(T)) levels are examined using the correlation consistent basis sets cc-pVXZ and cc-pCVXZ (X = D, T, Q, 5). It is found that, while BSSE has a negligible effect for valence-electron-only-correlated calculations for most basis sets, its magnitude becomes more pronounced for all-electron-correlated calculations, including core electrons. By utilizing the negligible effect of BSSE on the binding energies for valence-electron-only-correlated calculations, in combination with the negligible core-correlation effect at the CCSD(T) level, accurate binding energies of these clusters up to pentamers (octamers in the case of the Be clusters) are estimated via the basis set extrapolation of ab initio CCSD(T) correlation energies of the monomer and cluster with only the cc-pVDZ and cc-pVTZ sets, using the basis set and correlation-dependent extrapolation formula recently devised. A comparison between the CCSD(T) and density functional theory (DFT) binding energies is made to identify the most appropriate DFT method for the study of these clusters.     

Prediction of geometries and interaction energies of complexes formed by small molecules using semiempirical and ab initio methods

The accuracy of the semiempirical quantum mechanics methods (AM1 and PM3), and the ab initio methods (6-31G^** and MP2/6-31G^**) in predicting intermolecular geometries and interaction energies have been evaluated by detailed studies of 17 bimolecular complexes formed by small molecules. Comparisons between calculated and experimental geometries for 12 complexes are presented. It was found that AM1 gave reasonably good predictions of the geometries of complexes such as CH₄ · CH₄, which have very weak interactions, but it is not as good as other methods in predicting intermolecular geometry for complexes where hydrogen bonding interactions play an important role. This is consistent with its inability to reproduce the charge transfer in the formation of hydrogen bonds in these complexes.

PM3 is able to predict intermolecular geometries for most complexes, including those with hydrogen bonding; its major flaw is its tendency to overestimate the strength of the interactions between hydrogen atoms. Care should be taken therefore in using PM3 to study complicated molecular systems with multiple hydrogen atom interactions and the method's weakness in handling complexes in which electrostatic forces are important should also be noted.

Among ab initio methods, both the 6-31G^** and the MP2/6-31G^** were found to outperform AM1 and PM3 in prediction of intermolecular geometry. Both of these ab initio methods showed excellent consistency in geometry prediction for most of the complexes studied, although MP2/6-31G^** is better than 6-31G^**. It is noted that the MP2/6-31G^** did not produce the correct geometry for the CO₂· HF complex.

For 12 complexes for which experimental geometry data are available, AM1, PM3, 6-31G^**, and MP2/6-31G^** successfully predicted the geometry in 10, 12, 12, and 11 cases, respectively. The average errors given by AM1 in the predicted intermolecular distances were 0.264, 0.272, 0.091, and 0.061 Å, respectively. In comparison to the ab initio methods, AM1 and PM3 commonly underestimated the molecular interaction energy in such complexes by ˜ 1–2 kcal mol⁻¹.     

High-level ab initio studies of hydrogen abstraction from prototype hydrocarbon systems

Symmetric and nonsymmetric hydrogen abstraction reactions are studied using state-of-the-art ab initio electronic structure methods. Second-order M?ller-Plesset perturbation theory (MP2) and the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] methods with large correlation consistent basis sets (cc-pVXZ, where X = D,T,Q) are used in determining the transition-state geometries, activation barriers, and thermodynamic properties of several representative hydrogen abstraction reactions. The importance of basis set, electron correlation, and choice of zeroth-order reference wave function in the accurate prediction of activation barriers and reaction enthalpies are also investigated. The ethynyl radical (*CCH), which has a very high affinity for hydrogen atoms, is studied as a prototype hydrogen abstraction agent. Our high-level quantum mechanical computations indicate that hydrogen abstraction using the ethynyl radical has an activation energy of less than 3 kcal mol(-1) for hydrogens bonded to an sp(2) or sp(3) carbon. These low activation barriers further corroborate previous studies suggesting that ethynyl-type radicals would make good tooltips for abstracting hydrogens from diamondoid surfaces during mechanosynthesis. Modeling the diamond C(111) surface with isobutane and treating the ethynyl radical as a tooltip, hydrogen abstraction in this reaction is predicted to be barrierless.