Theoretical study of the N---H tautomerism in free base porphyrin

The N---H tautomerism of free base porphyrin is investigated at the semiempirical spin-unrestricted AM1 (UAM1) and ab initio RHF/3-21G levels. The UAM1 method provides delocalized geometries for all stationary structures without imposing any symmetry constraint. RHF/3-21G geometry optimizations have to be performed under symmetry restrictions to ensure that realistic delocalized structures are obtained. Both the semiempirical and the ab initio calculations predict that the interconversion between trans tautomers proceeds in an asynchronous two-step process via intermediate cis tautomers. The cis tautomers are characterized as minima in the potential energy surface and are 8–10 kcal mol⁻¹ higher in energy. The activation energy for the trans → cis interconversion is calculated to be approximately 23 kcal mol⁻¹ at the 3-21G level. The activation energy for the synchronous trans → trans interconversion is higher and has a value of 30.5 kcal mol⁻¹. The activation energies obtained at the semiempirical UAM1 level are twice as large as the ab initio values.     

Electronic structures and electronic spectra of the linkage isomers NSO and SNO

The relative stabilities and electronic structures of the linkage isomers NSO⁻ and SNO⁻ have been determined by the MNDO and ab initio Hartree—Fock—Slater methods. Both approaches predict a higher stability for SNO⁻ by ca. 100 kcal mol⁻¹, but an overlap population analysis indicates substantially higher bond orders for NSO⁻ compared to SNO⁻. The calculations also reveal a low energy pathway with a barrier of ca. 6 kcal mol⁻¹ for the isomerization process NSO⁻ → SNO⁻. Good agreement was found between the observed UV-visible absorption bands for NSO⁻ (λ_max 379 nm) and SNO⁻ (λ_max 340 nm) and calculated values of the electronic transition energies.     

Electronic structure and conformational properties of the amide linkage Part 5. Internal rotation and inversion in 1-formylaziridine

Internal rotation and nitrogen inversion in 1-formylaziridine (1) have been investigated by quantum mechanical (ab initio and MNDO) calculations, especially with respect to the variation of the geometry of the aziridine ring. While conformational stability is mainly determined by the n(N)/π(CO) interaction, the bond lengths within the ring are affected by the amount of interaction between the π(CO) orbital and the Walsh orbital ω_A. To separate the two types of interaction, calculations were also performed on formylcyclopropane (9). The torsional potential of 1 has a minimum close to the perpendicular conformation 1b. The two bisected conformations, 1a and 1c, are transition states for internal rotation. For nitrogen inversion, a barrier of 1.44 kcal mol⁻¹ (ab initio) was calculated. Calculations on 1-cyanoaziridine (7) gave inversion barriers of 5.81 (ab initio) and 12.31 kcal mol⁻¹ (MNDO). Probably due to methodical reasons the ab initio values seem to be too low, as calculations with different basis sets for aziridine indicate.     

An ab initio molecular orbital study on the gas-phase rapid reaction of the silicon cation Si with ammonia

The gas-phase rapid ion-molecule reaction Si⁺ (²P) + NH₃→ SiNH₂⁺ + H is theoretically investigated by the ab initio molecular orbital methods. Several possible pathways (A, B, C) on its potential energy surface have been examined, discussed and compared. Theoretical calculations indicate that pathway A is favourable in energy and that the reaction begins by forming a collision complex of the ion-dipole molecule Si-NH⁺₃, which forms with no barrier into the first energy well of the reaction coordinate. Migration of an H atom from an N atom to a Si atom forms the intermediate HSi-NH⁺₂, which corresponds to the second energy well and can fragment to the observed product SiNH⁺₂ by losing an H atom from the Si atom. The barriers for migration and fragmentation are 52.5 and 38.6 kcal mol⁻¹ respectively. Pathway A has a negative activation energy of −42.1 kcal mol⁻¹.     

Ab initio calculations of the rotational barriers in formamide and acetamide: The effects of polarization functions and correlation

Ab initio calculations have been used to determine the gas-phase rotational barrier about the CN bond in formamide and acetamide. The results indicate that the inclusion of polarization functions in the basis set leads to a substantial decrease (ca. 5 kcal mol⁻¹) in the calculated barrier height at the SCF level. Electron correlation effects decrease the barrier by less than 1 kcal mol⁻¹, while the addition of zero point energy corrections changes the barrier height only slightly. Based upon the current calculations, the 0 K rotational barriers for isolated formamide and acetamide are predicted to be 14.2 and 12.5 kcal mol⁻¹, respectively.     

Can chlorine anion catalyze the reaction of HOCl with HCl?

The reaction of HOCl + HCl → Cl₂ + H₂O in the presence of chlorine anion Cl⁻ has been studied using ab initio methods. The overall exothermicity is 15.5 kcal mol⁻¹ and this reaction has been shown to have a high activation barrier of 46.5 kcal mol⁻¹. Cl⁻ is found to catalyze the reaction via the formation of HOCl·Cl⁻, ClH·HOCl·Cl⁻ and Cl⁻·H₂) intermediate ion-molecule complexes or by interacting with a concerted four-center transition state of the reaction of HOCl + HCl.     

Theoretical studies on the potential energy surface and rovibrational states for the electronic ground state of carbonyl sulfide

A potential energy surface for the electronic ground state of carbonyl sulfide was optimized by using a self-consistent field-configuration interaction method and involving the recent observed vibrational band origins up to 8000 cm⁻¹ for the Σ state. The root mean square error for this refinement was found to be 0.27 cm⁻¹. The calculated quartic force constants from the refined potential are very close to the recent high level ab initio calculations. The vibrational energy levels for the Π and Δ states and for some isotopomers of carbonyl sulfide molecule were calculated to test the refined potential. The calculated energy levels are in good agreement with the experimental values.     

Ab initio quantum-chemical study of the unimolecular pyrolysis mechanisms of acetic acid

The mechanisms of unimolecular dehydration and decarboxylation reactions occurring during the pyrolysis of acetic acid above 700°C have been investigated by ab initio methods. The atomic basis set influence as well as the electron correlation effects are analyzed by using a variety of basis sets, ranging from minimal to polarized split-valence, and by introducing the Møller-Plesset (MP) perturbation theory. With an activation barrier of 76.0 kcal mol⁻¹, the concerted dehydration process occurs via a four-centre transition state. On the other hand, the decarboxylation process could be described by two different mechanisms depending on the nature of the kinetic experiments. While in flow systems, the decarboxylation of acetic acid takes place by a concerted mechanism via a four-centre transition state with an activation energy of 77.3 kcal mol⁻, the results suggest rather a water-catalyzed concerted mechanism via a six-membered transition state for the reaction carried out in batch systems, the activation barrier amounting to 64.0 kcal mol⁻¹.     

A model potential for the internal rotation of neighbouring rings of bithiophene and bipyrrole

Based upon ab initio Hartree-Fock calculations we propose a model potential for simplifying the study of the internal rotation of neighbouring rings. The quality of the proposed analytic potential was numerically verified in the case of bithiophene and bipyrrole by determining the deviations of the predicted energy values along the torsional angle with respect to 19 independent reference values. It is found that our model potential describes correctly the torsional motion the average deviations of the predicted energies along the torsional angle with respect to the ab initio reference energies are 0.20 and 0.22 kJ mol⁻¹ for bithiophene and bipyrrole, respectively. In both molecules, when going from a coplanar to a perpendicular conformation, two stable isomers have been detected. These isomers are separated by a barrier height of about 4 kJ mol⁻¹. Full optimization of the structural parameters along the torsional angle shows that they remain almost constant along the isomerization process.     

Thermochemistry of organosilicon compounds VII. Permethylcyclosilazanes and 1,1,3,3-tetramethyldisilazane

The enthalpy of formation (ΔH_f⁰), enthalpy of evaporation (ΔH_v⁰) and enthalpy of atomization (ΔH_a) of permethylcyclosilazanes (Me₂SiNH)_n (n = 3, 4) and 1,1,3,3-tetramethyldisilazane (Me₂SiH)₂NH have been determined. The enthalpies of formation of these compounds were compared with those calculated by the Benson-Buss-Franklin and Tatevskii additive schemes. In higher permethylcyclosilazanes the energy of the endocyclic Si---N bond is 306 ± 2 kJ mol⁻¹ (73 kcal mol⁻¹), that is 12 ± 2 kJ mol⁻¹ (3 kcal mol⁻¹) lower than the energy of the acyclic Si---N bond. The strain energy of the cyclotrisilazane ring is estimated to be 10.5 kJ mol⁻¹ (2.5 kcal mol⁻¹), whereas the energy of the ring Si---N bond is 295 kJ mol⁻¹ (70.5 kcal mol⁻¹).

The thermochemical data for permethylcyclosilazanes were compared with the corresponding values for permethylcyclosiloxanes calculated from the results of previously reported studies.     

Kinetic study of the reaction H2O2+H→H2O+OH by ab initio and density functional theory calculations

Theoretical investigations on the kinetics of the elementary reaction H₂O₂+H→H₂O+OH were performed using the transition state theory (TST). Ab initio (MP2//CASSCF) and density functional theory (B3LYP) methods were used with large basis set to predict the kinetic parameters; the classical barrier height and the pre-exponential factor. The ZPE and BSSE corrected value of the classical barrier height was predicted to be 4.1 kcal mol⁻¹ for MP2//CASSCF and 4.3 kcal mol⁻¹ for B3LYP calculations. The experimental value fitted from Arrhenius expressions ranges from 3.6 to 3.9 kcal mol⁻¹. Thermal rate constants of the title reaction, based on the ab initio and DFT calculations, was evaluated for temperature ranging from 200 to 2500 K assuming a direct reaction mechanism. The modeled ab initio-TST and DFT–TST rate constants calculated without tunneling were found to be in reasonable agreement with the observed ones indicating that the contribution of the tunneling effect to the reaction was predicted to be unimportant at ambient temperature.     

Barrier to internal rotation in 2,3- and 3,3-difluoropropenal and their methyl derivatives

The conformational stability and structure of 2,3-dimethylpropenal, 2,3-difluoropropenal and their 3,3-dimethyl and 3,3-difluoro derivatives were investigated utilizing ab initio calculations with 3-21G and 6-31G basis sets. For 2,3-dimethylpropenal and 3,3-difluoropropenal the s-trans was predicted to be the low-energy form. In the case of 3,3-dimethylpropenal and 2,3-difluoropropenal the s-cis was predicted by both levels of calculation to be the more stable conformer. Full optimization was performed at the transition states and the barriers to internal rotation were calculated. Methyl and fluorine substitution were found to significantly increase the barrier to interconversion in propenal. The relative change in the barrier depends on the position and the type of the substituent. The trans to cis barrier in 2,3-dimethylpropenal was calculated to be about 3 kcal mol⁻¹ greater than that in 3,3-dimethylpropenal, while the cis to trans barrier in 2,3-difluoropropenal was predicted to be about 7 kcal mol⁻¹ higher than the corresponding one in 3-3-difluoropropenal.     

An ab initio MO study of allene episulfide, cyclopropanethione and thioxyallyl

The electronic structures of allene episulfide, cyclopropanethione and thioxyallyl were examined by ab initio MO calculations and were compared with those of the corresponding oxygen compounds, allene oxide, cyclopropanone and oxyallyl. The difference in reactivities of allene episulfide and allene oxide was also speculatively estimated from the calculated electronic structures. The lowest singlet state of thioxyallyl was predicted to be the B₂ state, which corresponds to the σ, π-diradical. A small activation energy is required for the cyclization of the B₂ state to give allene episulfide. The A₁ singlet state lies 11 kcal mol⁻¹ higher than the B₂ singlet state and undergoes the disrotatory rotation of methylene groups to give cyclopropanethione with no activation energy.     

On the vibrational spectra and conformational stability of 1,1-dimethylhydrazine from temperature dependent FT-IR spectra of krypton solutions and ab initio calculations

Variable temperature (−105 to −150 °C) studies of the infrared spectra (3500–400 cm⁻¹) of 1,1-dimethylhydrazine, (CH₃)₂NNH₂, in liquid krypton have been carried out. No convincing spectral evidence could be found for the trans conformer which is expected to be at least 600 cm⁻¹ less stable than the gauche form. The structural parameters, dipole moments, conformational stability, vibrational frequencies, and infrared and Raman intensities have been predicted from MP2/6-31G(d) ab initio calculations. The predicted infrared and Raman spectra are compared to the experimental ones. The adjusted r₀ parameters from MP2/6-311+G(d,p) calculations are compared to those reported from an electron diffraction study. The energy differences between the gauche and trans conformers have been obtained from MP2 ab initio calculations as well as from density functional theory by the B3LYP method calculations from a variety of basis sets. All of these calculations indicate an energy difference of 650–900 cm⁻¹ with the B3LYP calculations predicted the larger values. The potential function governing the conformational interchange has been predicting from both types of calculations and comparisons have been made. The barrier to internal rotation by the independent rotor model of the inner methyl group is predicted to have a value of 1812 cm⁻¹ and that of the outer one of 1662 cm⁻¹ from ab initio MP2/6-31G(d) calculations. These values agree well with the experimentally determined values of 1852±16 and 1558±12 cm⁻¹, respectively, from a fit of the torsional transitions with the coupled rotor model. For the coupled rotor model the predicted V′₃₃ (sin 3τ₀ sin 3τ₁ term) value which ranged from 190 to 232 cm⁻¹ is in reasonable agreement with the experimental value of 268±3 cm⁻¹ but the predicted V₃₃ (cos 3τ₀ cos 3τ₁ term) value of −73 to −139 cm⁻¹ is 25% smaller and of the opposite sign of the experimental value of 333±22 cm⁻¹. These theoretical and spectroscopy results are compared to similar quantities of some corresponding molecules.     

Effect of small cations on the hydrogen bond between an N-aromatic heterocycle and amine

The influence of an Li⁺ ion on the structure and bonding of the H-bond interaction of an N-aromatic heterocycle and a singly bonded N-amine is studied by ab initio SCF calculations with imidazole (Imh) and NH₃ serving as models of the two families. Full optimization of structures have been carried out for DZ, DZP and TZP basis sets. The computed H-bond energy for Imh/ NH₃ of −6.4 kcal/mol is in close agreement with a recent experiment. An appreciable three-body interaction of −3.8 kcal/mol is found for the complex Imh/NH₃/Li⁺.     

Equilibrium structure of the carbon dioxide-water complex in the gas phase: an ab initio and density functional study

High-level ab initio (MP2/6-311++G(2d,2p) geometry, Gaussian-2, MP4(SDTQ) and QCISD(T) binding energies) and density-functional (Becke3LYP/6-311++G(2df,2pd)) calculations have been performed on the charge-transfer complex between water and carbon dioxide. The complex appears to have two equivalent non-planar minima of C_s symmetry. Minima are separated by transition states with C₁ symmetry, whereas the totally planar structure with C_2v symmetry is a second-order transition state. All the critical points lie at approximately the same energy (less than 0.05 Kj mol⁻¹ difference). Therefore, the experimentally observable structure should be planar. The best equilibrium intermolecular distance for this complex calculated at the MP2/6-311++G(2d,2p) level is 2.800 Å. Our best estimate of the observable intermolecular distance (corrected for anharmonicity) is 2.84 Å, in agreement with the experimentally derived value of 2.836 Å. Our best estimate of the binding energy at the QCISD(T) level, taking into account the variation of the distance owing to anharmonicity and the use of more sophisticated theoretical treatments, is −12.0 ± 0.2 kJ mol⁻¹. Our best estimate of the barrier to internal rotation, also at the MP2/6-311++G(2d,2p) level, is 4.0 kJ mol⁻¹, outside the error limits of the experimental determination (3.64 ± 0.04 kJ mol⁻¹). Density functional theory at the level employed here gives an equilibrium intermolecular distance that is too large (2.857 Å), a binding energy that is too small (8.1 kJ mol⁻¹), attributable neither to geometry nor to the basis set, and also a barrier to internal rotation that is slightly too small (3.39 kJ mol⁻¹). The overall picture is, however, reasonably good.     

An ab initio study on the equilibrium structure and torsional potential energy function of dinitrogen tetroxide

The molecular parameters of dinitrogen tetroxide, N₂O₄, have been determined in large-scale ab initio calculations using the multiconfigurational second-order perturbation method, CASSCF/CASPT2, and basis sets of double- to quadruple-zeta quality. With the largest basis set employed, cc-pVQZ for nitrogen and cc-pVTZ for oxygen, the structural equilibrium parameters are determined to be r(NN) = 1.7940 Å, r(NO) = 1.1906 Å and (NNO) = 112.55°. The potential energy barrier at the staggered conformation of the molecule is found to be 2313 cm⁻¹, and the binding energy of the NN bond is calculated to be 4616 cm⁻¹ (13.2 kcal/mol).     

Ab initio and density function theory computational studies of the CH4 + H → CH3 + H2 reaction

The activation barrier for the CH₄ + H → CH₃ + H₂ reaction was evaluated with traditional ab initio and Density Functional Theory (DFT) methods. None of the applied ab initio and DFT methods was able to reproduce the experimental activation barrier of 11.0-12.0 kcal/mol. All ab initio methods (HF, MP2, MP3, MP4, QCISD, QCISD(T), G1, G2, and G2MP2) overestimated the activation energy. The best results were obtained with the G2 and G2MP2 ab initio computational approaches. The zero-point corrected energy was 14.4 kcal mol⁻¹. Some of the exchange DFT methods (HFB) computed energies which were similar to the highly accurate ab initio methods, while the B3LYP hybrid DFT methods underestimated the activation barrier by 3 kcal mol⁻¹. Gradient-corrected DFT methods underestimated the barrier even more. The gradient-corrected DFT method that incorporated the PW91 correlational functional even generated a negative reaction barrier. The suitability of some computational methods for accurately predicting the potential energy surface for this hydrogen radical abstraction reaction was discussed.     

Structure and conformation of cyclobutylsilane: an electron diffraction and ab initio study

A gas electron diffraction study of cyclobutylsilane results in a mixture of equatorial and axial conformers, with the equatorial confomer slightly more stable (Δ G = 0.8 ± 0.4 kJ mol⁻¹). The cyclobutyl ring is distorted with the adjacent bonds longer (C₁---C₂ = 1.573 (4) Å) than the opposite bonds (C₂---C₃ = 1.557 (4) Å). The experimental values for the energy difference between the two conformers and for the geometric parameters are reproduced very well by ab initio calculations. The importance of silicon 3d orbitals in the interpretation of ring distortion is ambiguous, but on the basis of the ab initio calculations the participation of silicon 3d functions is negligible.     

Ab initio study of the hydration of carbon dioxide: Additional comments based on refined calculations

Ab initio calculations on the formation of carbonic acid from the hydration of carbon dioxide with water dimer are re-examined. Fully optimized geometries of the three stationary points (minima and transition state) with the 3-21G basis set are reported. They possess non-planar structures. The inclusion of polarization (with the 6-31G* basis) and electron correlation (via Møller-Plesset perturbation theory to second through to fourth-order using the 6-31G basis) tends to enlarge the energy barrier (35–40 kcal mol⁻¹) for the double hydrogen transfer. This suggests that the neutral hydrolysis of CO₂ could require more water molecules (an oligomer) in an autocatalytic process rather than a dimer.