Theoretical reference values for the AE6 and BH6 test sets from explicitly correlated coupled-cluster theory

We propose a new computational protocol to obtain highly accurate theoretical reference data. This protocol employs the explicitly correlated coupled-cluster method with iterative single and double excitations as well as perturbative triple excitations, CCSD(T)(F12), using quadruple-z\zeta basis sets. Higher excitations are accounted for by conventional CCSDT(Q) calculations using double-z\zeta basis sets, while core/core-valence correlation effects are estimated by conventional CCSD(T) calculations using quadruple-z\zeta basis sets. Finally, scalar-relativistic effects are accounted for by conventional CCSD(T) calculations using triple-z\zeta basis sets. In the present article, this protocol is applied to the popular test sets AE6 and BH6. An error analysis shows that the new reference values obtained by our computational protocol have an uncertainty of less than 1 kcal/mol (chemical accuracy). Furthermore, concerning the atomization energies, a cancellation of the basis set incompleteness error in the CCSD(T)(F12) perturbative triples contribution with the corresponding error in the contribution from higher excitations is observed. This error cancellation is diminished by the CCSD(T*)(F12) method. Thus, we recommend the use of the CCSD(T*)(F12) method only for small- and medium-sized basis sets, while the CCSD(T)(F12) approach is preferred for high-accuracy calculations in large basis sets.     

Quantum dynamics of vibrationally activated OH-CO reactant complexes

The MP2 complete basis set (CBS) limit for the binding energy of the two low-lying water octamer isomers of D2d and S4 symmetry is estimated at -72.7+/-0.4 kcal/mol using the family of augmented correlation-consistent orbital basis sets of double through quintuple zeta quality. The largest MP2 calculation with the augmented quintuple zeta (aug-cc-pV5Z) basis set produced binding energies of -73.70 (D2d) and -73.67 kcal/mol (S4). The effects of higher correlation, computed at the CCSD(T) level of theory, are estimated at <0.1 kcal/mol. The newly established MP2/CBS limit for the water octamer is reproduced quite accurately by the newly developed all atom polarizable, flexible interaction potential (TTM2-F). The TTM2-F binding energies of -73.21 (D2d) and -73.24 kcal/mol (S4) for the two isomers are just 0.5 kcal/mol (or 0.7%) larger than the MP2/CBS limit.     

Torsional potential of 4,4'-bipyridine: ab initio analysis of dispersion and vibrational effects

Ab initio calculations using restricted Hartree-Fock, second-order M?ller-Plesset perturbation theory (MP2), density-functional theory (DFT), and coupled-cluster methods have been done to obtain the torsional potential-energy profile of the aza-aromatic molecule 4,4'-bipyridine. The torsional potential is evaluated adiabatically by fixing the normalized sum of the dihedral angles through the C-C inter-ring bond at several values along the torsional path and relaxing the remaining degrees of freedom. Previous discrepancies between MP2 and DFT internal rotation barrier heights are removed, and seen to be mostly due to the underestimation of the dispersion energy in the coplanar conformer by MP2 when using relatively small basis sets. The calculations indicate that the barrier height between the twisted global minimum and the 0 degrees conformer is around 1.5-1.8 kcal mol-1 while that corresponding to the 90 degrees one is about 2.0-2.2 kcal mol-1. This same relative energy ordering of the coplanar and perpendicular conformers was experimentally derived from nuclear magnetic resonance (NMR) measurements of 1H dipolar couplings on 4,4'-bipyridine solutions in a nematic liquid crystal, although the barrier heights are much lower than those estimated from NMR experiments in the gas phase. The DFT infrared spectrum and zero-point vibrational energy corrections to the torsional energy profile have also been calculated, the latter having a small influence on the torsional potential-energy profiles.     

Reactions of hydrogen atom with hydrogen peroxide

Rate coefficients are calculated using canonical variational transition state theory with multidimensional tunneling (CVT/SCT) for the reactions H + H2O2 --> H2O + OH (1a) and H + H2O2 --> HO2 + H2 (1b). Reaction barrier heights are determined using two theoretical approaches: (i) comparison of parametrized rate coefficient calculations employing CVT/SCT to experiment and (ii) high-level ab initio methods. The evaluated experimental data reveal considerable variations of the barrier height for the first reaction: although the zero-point-exclusive barrier for (1a) derived from the data by Klemm et al. (First Int. Chem. Kinet. Symposium 1975, 61) is 4.6 kcal/mol, other available measurements result in a higher barrier of 6.2 kcal/mol. The empirically derived zero-point-exclusive barrier for (1b) is 10.4 kcal/mol. The electronic structure of the system at transition state geometries in both reactions was found to have "multireference" character; therefore special care was taken when analyzing electronic structure calculations. Transition state geometries are optimized by multireference perturbation theory (MRMP2) with a variety of one-electron basis sets, and by a multireference coupled cluster (MR-AQCCSD) method. A variety of single-reference benchmark-level calculations have also been carried out; included among them are BMC-CCSD, G3SX(MP3), G3SX, G3, G2, MCG3, CBS-APNO, CBS-Q, CBS-QB3, and CCSD(T). Our data obtained at the MRMP2 level are the most complete; the barrier height for (1a) using MRMP2 at the infinite basis set limit is 4.8 kcal/mol. Results are also obtained with midlevel single-reference multicoefficient correlation methods, such as MC3BB, MC3MPW, MC-QCISD/3, and MC-QCISD-MPWB, and with a variety of hybrid density functional methods, which are compared with high-level theory. On the basis of the evaluated experimental values and the benchmark calculations, two possible recommended values are given for the rate coefficients.     

Heats of formation of krypton fluorides and stability predictions for KrF4 and KrF6 from high level electronic structure calculations

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for KrF+, KrF-, KrF2, KrF3+, KrF4, KrF5+, and KrF6 from coupled-cluster theory (CCSD(T)) calculations with effective core potential correlation-consistent basis sets for krypton. To achieve near chemical accuracy (+/-1 kcal/mol), three corrections were added to the complete basis set binding energies based on frozen core coupled-cluster theory energies: a correction for core-valence effects, a correction for scalar relativistic effects, and a correction for first-order atomic spin-orbit effects. Vibrational zero point energies were computed at the coupled-cluster level of theory. The calculated value for the heat of formation of KrF2 is in excellent agreement with the experimental value. Contrary to the analogous xenon fluorides, KrF2, KrF4, and KrF6 are predicted to be thermodynamically unstable with respect to loss of F2. An analysis of the energetics of KrF4 and KrF6 with respect to fluorine atom loss together with calculations of the transition states for the intramolecular loss of F2 show that fluorine atom loss is the limiting factor determining the kinetic stabilities of these molecules. Whereas KrF4 possesses a marginal energy barrier of 10 kcal/mol toward fluorine atom loss and might be stable at moderately low temperatures, the corresponding barrier in KrF6 is only 0.9 kcal/mol, suggesting that it could exist only at very low temperatures. Although the simultaneous reactions of either two or four fluorine atoms with KrF2 to give KrF4 or KrF6, respectively, are exothermic, they do not represent feasible synthetic approaches because the attack of the fluorine ligands of KrF2 by the fluorine atoms, resulting in F2 abstraction, is thermodynamically favored over oxidative fluorination of the krypton central atom. Therefore, KrF6 could exist only at very low temperatures, and even the preparation of KrF4 will be extremely difficult.     

Accurate ab initio potential energy surfaces for the 3A' and 3A' electronic states of the O(3P)+HBr system

In this work, we report the construction of potential energy surfaces for the (3)A(') and (3)A(') states of the system O((3)P) + HBr. These surfaces are based on extensive ab initio calculations employing the MRCI+Q/CBS+SO level of theory. The complete basis set energies were estimated from extrapolation of MRCI+Q/aug-cc-VnZ(-PP) (n = Q, 5) results and corrections due to spin-orbit effects obtained at the CASSCF/aug-cc-pVTZ(-PP) level of theory. These energies, calculated over a region of the configuration space relevant to the study of the reaction O((3)P) + HBr → OH + Br, were used to generate functions based on the many-body expansion. The three-body potentials were interpolated using the reproducing kernel Hilbert space method. The resulting surface for the (3)A(') electronic state contains van der Waals minima on the entrance and exit channels and a transition state 6.55 kcal/mol higher than the reactants. This barrier height was then scaled to reproduce the value of 5.01 kcal/mol, which was estimated from coupled cluster benchmark calculations performed to include high-order and core-valence correlation, as well as scalar relativistic effects. The (3)A(') surface was also scaled, based on the fact that in the collinear saddle point geometry these two electronic states are degenerate. The vibrationally adiabatic barrier heights are 3.44 kcal/mol for the (3)A(') and 4.16 kcal/mol for the (3)A(') state.     

Coupled-cluster study of the electronic structure and energetics of tetrasulfur, S4

Ab initio electronic structure calculations are reported for S4. Geometric and energetic parameters are calculated using the singles and doubles coupled-cluster method, including a perturbutional correction for connected triple excitation, CCSD(T), together with systematic sequences of correlation consistent basis sets extrapolated to the complete basis set limit. The geometry for the ground state singlet C2v structure of S4 is in good agreement with the microwave structure determined for S4. There is a low-lying D2h transition state at 1.6 kcal/mol which interchanges the long S-S bond. S4 has a low-lying triplet state (3B 1u) in D2h symmetry which is 10.8 kcal/mol above the C2v singlet ground state. The S-S bond dissociation energy for S4 into two S2(3Sigma*g) molecules is predicted to be 22.8 kcal mol(-1). The S-S bond energy to form S3+S(3P) is predicted to be 64 kcal/mol.     

Ab initio calculation of the inversion barrier in the germyl radical

Equilibrium geometry, vibrational frequencies and inversion barrier of the GeH₃ radical have been calculated using ab initio molecular orbital methods at both the SCF and correlated levels. At the SCF level, the effect of several small and medium size basis sets of Ge on these properties are studied systematically. Electron correlation corrections as incorporatedvia fourth-order unrestricted Møller-Plesset perturbation theory in conjunction with large Gaussian basis sets were found to reduce the GeH₃ inversion barrier appreciably relative to the UHF values. The final barrier height of 4.5 kcal/mol can be compared with a recent estimate of 4.4 kcal/mol inferred from the REMPI study.     

Cepa calculations of potential energy surfaces for open-shell systems. I. The reaction of O(3P) with H2(1Σg+)

Ab initio calculations including electron correlation effects (mainly on CEPA-PNO level) have been performed for the potential energy surface of the reaction of ³P oxygen atoms with molecular hydrogen. The collinear abstraction reaction (C_∞v symmetry) and the vertical insertion reaction (C_2v) have been investigated with particular emphasis.The influence of the orbital basis size and of electron correlation both on the reaction energy and on the barrier height and location for the abstraction reaction has been studied in some detail. Our extrapolated value for the barrier for this reaction is 13.4 ± 1.0 kcal/mol, in fair agreement with the experimental activation energy, while the insertion reaction has to pass over a barrier of =48 kcal/mol. The analysis of electron correlation effects reveals that it is compulsory to include all singly and doubly substituted configurations, to correct for unlinked cluster contributions and to use fairly large basis sets if one wants to get accurate ab initio potential surfaces for the reactions of triplet oxygen atoms.     

Ab initio calculations on the barrier height for the hydrogen addition to ethylene and formaldehyde. The importance of spin projection

Heats of reaction and barrier heights have been computed for H + CH₂CH₂ → C₂H₅, H + CH₂O → CH₃O, and H + CH₂O → CH₂OH using unrestricted Hartree-Fock and Møller–Plesset perturbation theory up to fourth order (with and without spin annihilation), using single-reference configuration interaction, and using multiconfiguration self-consistent field methods with 3-21G, 6-31G(d), 6-31G(d,p), and 6-311G(d,p) basis sets. The barrier height in all three reactions appears to be relatively insensitive to the basis sets, but the heats of reaction are affected by p-type polarization functions on hydrogen. Computation of the harmonic vibrational frequencies and infrared intensities with two sets of polarization functions on heavy atoms [6-31G(2d)] improves the agreement with experiment. The experimental barrier height for H + C₂H₄ (2.04 ± 0.08 kcal/mol) is overestimated by 7?9 kcal/mol at the MP2, MP3, and MP4 levels. MCSCF and CISD calculations lower the barrier height by approximately 4 kcal/mol relative to the MP4 calculations but are still almost 4 kcal/mol too high compared to experiment. Annihilation of the largest spin contaminant lowers the MP4SDTQ computed barrier height by 8?9 kcal/mol. For the hydrogen addition to formaldehyde, the same trends are observed. The overestimation of the barrier height with Møller-Plesset perdicted barrier heights for H + C₂H₄ → C₂H₅, H + CH₂O → CH₃O, and H + CH₂O → CH₂OH at the MP4SDTQ /6-31G(d) after spin annihilation are respectively 1.8, 4.6, and 10.5 kcal/mol.     

Predicting new molecular species of potential interest to atmospheric chemistry: the isomers HSBr and HBrS

Coupled cluster singles and doubles with perturbative contributions of connected triples CCSD(T) theory with a series of correlation consistent basis sets was used to predict the existence and characterize for the first time the structures, harmonic frequencies, and energetic quantities of the isomeric species HSBr and HBrS, as well as the transition state connecting them. These calculations consider extrapolation to the complete basis set (CBS) limit, corrections for scalar relativistic effects using the second-order Douglas-Kroll-Hess Hamiltonian, and also correlation of the bromine d electrons in addition to the 14 valence electrons. The species HSBr was found to be more stable than HBrS by 50.93 kcal/mol, with a high barrier height of 60.00 kcal/mol for the interconversion into HBrS. The smaller barrier of 7.90 kcal/mol (ZPE included) for the reverse process, however, should favor a rapid interconversion of HBrS into HSBr if HBrS can also be initially present in a potential synthetic route. If trapped in a matrix, their harmonic frequencies will allow for an unambiguous distinction between the two species. Scalar relativistic corrections and correlation of 24 electrons, although minor for the present purpose of a first time, but accurate, characterization of these species, are needed if chemical accuracy is also pursued. A test of the DFT/B3LYP approach in describing this type of system resulted in good energetic quantities, but geometric parameters and frequencies still lack spectroscopic accuracy. Whether HSBr can act as a temporary bromine reservoir and/or a source of reactive bromine and HS radicals requires further studies that are underway in our group.     

Molecular mechanism for H2 release from BH3NH3, including the catalytic role of the Lewis acid BH3

Electronic structure calculations using various methods, up to the coupled-cluster CCSD(T) level, in conjunction with the aug-cc-pVnZ basis sets with n = D, T, and Q, extrapolated to the complete basis set limit, show that the borane molecule (BH3) can act as an efficient bifunctional acid-base catalyst in the H2 elimination reactions of XHnYHn systems (X, Y = C, B, N). Such a catalyst is needed as the generation of H2 from isoelectronic ethane and borane amine compounds proceeds with an energy barrier much higher than that of the X-Y bond energy. The asymptotic energy barrier for H2 release is reduced from 36.4 kcal/mol in BH3NH3 to 6.0 kcal/mol with the presence of BH3 relative to the molecular asymptote. The NH3 molecule can also participate in a similar catalytic process but induces a smaller reduction of the energy barrier. The kinetics of these processes was analyzed by both transition-state and RRKM theory. The catalytic effect of BH3 has also been probed by an analysis of the electronic densities of the transition structures using the atom-in-molecule (AIM) and electron localization function (ELF) approaches.     

Is HO3 minimum cis or trans? An analytic full-dimensional ab initio isomerization path

The minimum energy path for isomerization of HO(3) has been explored in detail using accurate high-level ab initio methods and techniques for extrapolation to the complete basis set limit. In agreement with other reports, the best estimates from both valence-only and all-electron single-reference methods here utilized predict the minimum of the cis-HO(3) isomer to be deeper than the trans-HO(3) one. They also show that the energy varies by less than 1 kcal mol(-1) or so over the full isomerization path. A similar result is found from valence-only multireference configuration interaction calculations with the size-extensive Davidson correction and a correlation consistent triple-zeta basis, which predict the energy difference between the two isomers to be of only Δ = -0.1 kcal mol(-1). However, single-point multireference calculations carried out at the optimum triple-zeta geometry with basis sets of the correlation consistent family but cardinal numbers up to X = 6 lead upon a dual-level extrapolation to the complete basis set limit of Δ = (0.12 ± 0.05) kcal mol(-1). In turn, extrapolations with the all-electron single-reference coupled-cluster method including the perturbative triples correction yield values of Δ = -0.19 and -0.03 kcal mol(-1) when done from triple-quadruple and quadruple-quintuple zeta pairs with two basis sets of increasing quality, namely cc-cpVXZ and aug-cc-pVXZ. Yet, if added a value of 0.25 kcal mol(-1) that accounts for the effect of triple and perturbative quadruple excitations with the VTZ basis set, one obtains a coupled cluster estimate of Δ = (0.14 ± 0.08) kcal mol(-1). It is then shown for the first time from systematic ab initio calculations that the trans-HO(3) isomer is more stable than the cis one, in agreement with the available experimental evidence. Inclusion of the best reported zero-point energy difference (0.382 kcal mol(-1)) from multireference configuration interaction calculations enhances further the relative stability to ΔE(ZPE) = (0.51 ± 0.08) kcal mol(-1). A scheme is also suggested to model the full-dimensional isomerization potential-energy surface using a quadratic expansion that is parametrically represented by a Fourier analysis in the torsion angle. The method illustrated at the raw and complete basis-set limit coupled-cluster levels can provide a valuable tool for a future analysis of the available (incomplete thus far) experimental rovibrational data.     

Theoretical study for the HNSi–HSiN isomerization: ab initio and DFT calculations

Ab initio and density functional theory (DFT) calculations using the GAUSSIAN 94 program have been performed to investigate the molecular structures of HNSi and HSiN in the ground state as well as the transition state for the HNSi–HSiN isomerization reaction at the 6-311G(d,p), 6-311+G(2d,p) and 6-311+G(2df,p) basis sets. The results show that DFT calculations at higher levels of theory reproduce experimental vibrational frequencies of both HNSi and HSiN better than ab initio methods including electron correlation effects. Those calculated geometries are accurate enough to predict the rotational constant of HNSi. The barrier height for the isomerization reaction is found to be about 10 kcal/mol.     

Benchmark calculations of reaction energies, barrier heights, and transition-state geometries for hydrogen abstraction from methanol by a hydrogen atom

We report benchmark calculations of reaction energies, barrier heights, and transition-state geometries for the reaction of CH(3)OH with H to produce CH(2)OH and H(2). Highly accurate composite methods, such as CBS, G2, G3S, G3X, G3SX, and multi-coefficient correlation methods (MCCMs), are used to calibrate lower-cost methods. We also performed single-level CCSD(T) calculations extrapolated to the infinite-basis limit on the basis of aug-cc-pVXZ (X = 3, 4) correlation consistent basis sets. The benchmark high-level calculations give consensus values of the forward reaction barrier height and the reaction energy of 9.7 kcal/mol and - 6.4 kcal/mol, respectively. To evaluate the accuracy of cost-efficient methods that are potentially useful for dynamics studies of the title reaction, we further include the results obtained by hybrid density functional theory methods and hybrid meta density functional theory methods that have recently been designed for chemical kinetics. Results obtained by popular semiempirical methods are also given for comparison. On the basis of the benchmark gas-phase results, we suggest MC-QCISD/3, MC3BB, and BB1K as reasonably accurate and affordable electronic structure methods for calculating dynamics for the title reaction.     

A systematic ab initio study of the water dimer in hierarchies of basis sets and correlation models

A systematic, high-level ab initio investigation of the water dimer has been performed. The oxygen-oxygen bond distance has been estimated to be around 2.90 ?, about 0.05 ? shorter than the experimentally estimated distance, challenging the accuracy of the latter. The interaction energy has been obtained at −5.0±0.1 kcal/mol, which compares favourably with the experimentally estimated value of −5.4±0.7 kcal/mol. The importance of employing basis sets that include diffuse functions in correlated calculations on hydrogen-bonded systems is confirmed. In correlated calculations on the water dimer and the hydrogen fluoride dimer, the counterpoise-corrected interaction energies converge considerably slower towards the basis set limit than do the uncorrected energies, provided that the correlation-consistent basis sets are augmented with diffuse functions. Received: 12 February 1997 / Accepted: 5 June 1997     

Theoretical study of transition state structure and reaction enthalpy of the F + H2-->HF + H reaction by a diffusion quantum Monte Carlo approach

Ab initio calculations of transition state structure and reaction enthalpy of the F + H2-->HF + H reaction has been carried out by the fixed-node diffusion quantum Monte Carlo method in this study. The Monte Carlo sampling is based on the Ornstein-Uhlenbeck random walks guided by a trial wave function constructed from the floating spherical Gaussian orbitals and spherical Gaussian geminals. The Monte Carlo calculated barrier height of 1.09(16) kcal/mol is consistent with the experimental values, 0.86(10)/1.18(10) kcal/mol, and the calculated value from the multireference-type coupled-cluster (MRCC) calculation with the aug-cc-pVQZ(F)/cc-pVQZ(H) basis set, 1.11 kcal/mol. The Monte Carlo-based calculation also gives a similar value of the reaction enthalpy, -32.00(4) kcal/mol, compared with the experimental value, -32.06(17) kcal/mol, and the calculated value from a MRCC/aug-cc-pVQZ(F)/cc-pVQZ(H) calculation, -31.94 kcal/mol. This study clearly indicates a further application of the random-walk-based approach in the field of quantum chemical calculation.     

Computational study of the release of H2 from ammonia borane dimer (BH3NH3)2 and its ion pair isomers

High-level electronic structure calculations have been used to map out the relevant portions of the potential energy surfaces for the release of H2 from dimers of ammonia borane, BH3NH3 (AB). Using the correlation-consistent aug-cc-pVTZ basis set at the second-order perturbation MP2 level, geometries of stationary points were optimized. Relative energies were computed at these points using coupled-cluster CCSD(T) theory with the correlation-consistent basis sets at least up to the aug-cc-pVTZ level and in some cases extrapolated to the complete basis set limit. The results show that there are a number of possible dimers involving different types of hydrogen-bonded interactions. The most stable gaseous phase (AB)2 dimer results from a head-to-tail cyclic conformation and is stabilized by 14.0 kcal/mol with respect to two AB monomers. (AB)2 can generate one or two H2 molecules via several direct pathways with energy barriers ranging from 44 to 50 kcal/mol. The diammoniate of diborane ion pair isomer, [BH4-][NH3BH2NH3+] (DADB), is 10.6 kcal/mol less stable than (AB)2 and can be formed from two AB monomers by overcoming an energy barrier of approximately 26 kcal/mol. DADB can also be generated from successive additions of two NH3 molecules to B2H6 and from condensation of AB with separated BH3 and NH3 molecules. The pathway for H2 elimination from DADB is characterized by a smaller energy barrier of 20.1 kcal/mol. The alternative ion pair [NH4+][BH3NH2BH3-] is calculated to be 16.4 kcal/mol above (AB)2 and undergoes H2 release with an energy barrier of 17.7 kcal/mol. H2 elimination from both ion pair isomers yields the chain BH3NH2BH2NH3 as product. Our results suggest that the neutral dimer will play a minor role in the release of H2 from ammonia borane, with a dominant role from the ion pairs as observed experimentally in ionic liquids and the solid state.     

An ab initio benchmark study of the H + CO --> HCO reaction

The H + CO --> HCO reaction has been characterized with correlation consistent basis sets at five levels of theory in order to benchmark the sensitivities of the barrier height and reaction ergicity to the one-electron and n-electron expansions of the electronic wave function. Single and multireference methods are compared and contrasted. The coupled cluster method RCCSD(T) was found to be in very good agreement with Davidson-corrected internally-contracted multireference configuration interaction (MRCI+Q). Second-order Moller-Plesset perturbation theory (MP2) was also employed. The estimated complete basis set (CBS) limits for the barrier height (in kcal/mol) for the five methods, including harmonic zero-point energy corrections, are MP2, 4.66; RCCSD, 4.78; RCCSD(T), 4.15; MRCI, 5.10; and MRCI+Q, 4.07. Similarly, the estimated CBS limits for the ergicity of the reaction are: MP2, -17.99; RCCSD, -13.34; RCCSD(T), -13.79; MRCI, -11.46; and MRCI+Q, -13.70. Additional basis set explorations for the RCCSD(T) method demonstrate that aug-cc-pVTZ sets, even with some functions removed, are sufficient to reproduce the CBS limits to within 0.1-0.3 kcal/mol.     

Potential energy scans and vibrational assignments of cyclopropanecarboxylic acid and cyclopropanecarboxamide

The structural stability and internal rotations in cyclopropanecarboxylic acid and cyclopropanecarboxamide were investigated by the DFT-B3LYP and the ab initio MP2 calculations using 6-311G** and 6-311+G** basis sets. The computations were extended to the MP4//MP2/6-311G** and CCSD(T)//MP2/6-311G** single-point calculations. From the calculations the molecules were predicted to exist predominantly in the cis (C=O group eclipses the cyclopropane ring) with a cis-trans barrier of about 4-6kcal/mol. The OCOH torsional barrier in the acid was estimated to be about 12-13kcal/mol while the corresponding OCNH torsional barrier in the amide was calculated to be about 20kcal/mol. The equilibrium constant k for the cis<-->trans interconversion in cyclopropanecarboxylic acid was calculated to be 0.1729 at 298.15K that corresponds to an equilibrium mixture of about 85% cis and 15% trans. The vibrational frequencies were computed at the DFT-B3LYP level. Normal coordinate calculations were carried out and potential energy distributions were calculated for the low energy cis conformer of the molecules. Complete vibrational assignments were made on the basis of normal coordinate calculations and comparison with experimental data of the molecules.