Semiempirical and ab initio calculations on the alkaline hydrolysis of the β-lactam ring Influence of the solvent

We used semiempirical and ab initio calculations to investigate the nucleophilic attack of the hydroxyl ion on the β-lactam carbonyl group. Both allowed us to detect reaction intermediates pertaining to proton-transfer reactions. We also used ab initio calculations and the PM3 semiempirical method to investigate the influence of the solvent on the process. The AMSOL method predicts the occurrence of a potential energy barrier of 20.7 kcal mol⁻¹ due to the desolvation of the hydroxyl ion in approaching the β-lactam carbonyl group. Using the supermolecular approach and a water solvation sphere of 20 molecules around the solute, the potential energy barrier is lowered to 17.5 kcal mol⁻¹. Ab initio calculations using the SCRF method predict a potential energy barrier of 13.6 kcal mol⁻¹. These three values, especially the last two, are very close to the experimental value of 16.7 kcal mol⁻¹.     

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.     

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.     

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.     

Ab initio calculations of the vibrational spectra of the ammonia and fluoramide complexes with HF

The structures, energetics, vibrational frequencies and IR intensities of the H₃N HF, H₃N F₂ and NH₂FHF (three isomers) complexes were examined using the self-consistent field method within the 6-311G^** basis set. The interaction energies were calculated using the MP2 approach. The results are compared with monomer calculations and experimental data. The complex NH₂FHF was found to exist in three forms: one with the HF molecule hydrogen bonded to the nitrogen lone pair of NH₂F (D₀ =7.403 kcal mol⁻¹), another a complex formed through the F atom lone pair (D₀=4.698 kcal mol⁻¹) and third a cyclic structure (D₀=5.644 kcal mol⁻¹).     

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.     

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⁻¹.     

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⁻¹.     

MNDO study of the proton affinity of fluorinated formaldehydes and acetones

The MNDO molecular orbital method is employed to calculate the proton affinities of fluorinated formaldehydes and acetones. Agreement with experimentally reported proton affinities is good. In the acetone series a decrease in proton affinity is calculated for each successive fluorine substituent. The calculated strength of the intramolecular hydrogen bond in the protonated fluoro-formaldehydes and acetones is 0.6—2.7 kcal mol⁻¹, in good agreement with the experimental value of 2—3 kcal mol⁻¹ in the protonated fluoroacetones. Examination of the calculated charge distribution shows that the trends in proton affinity can be understood qualitatively both in terms of initial-state and final-state effects caused by the fluorine substituents. Protonation at the fluorine atom is less stable by about 25 kcal mol⁻¹ than protonation at the oxygen atom for monofluoroacetone.     

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.     

Ab initio study of the two competitive channels HO2 + H → H2 + O2 and HO2 + H → H2O + O

Saddle point geometries and barrier heights have been calculated for the H abstraction reaction HO₂(²A″)+H(²S) → H₂(¹Σ⁺_g)+O₂(³Σ⁻_g) and the concerted H approach-O removing reaction HO₂ (²A″)+H(²S) → H₂O(¹A₁)+O(³P) by using SDCI wavefunctions with a valence double-zeta plus polarization basis set. The saddle points are found to be of C_s symmetry and the barrier heights are respectively 5.3 and 19.8 kcal by including size consistent correction. Moreoever kinetic parameters have been evaluated within the framework of the TST theory. So activation energies and the rate constants are estimated to be respectively 2.3 kcal and 0.4×10⁹ ℓ mol⁻¹ s⁻¹ for the first reaction, 20.0 kcal and 5.4.10⁻⁵ ℓ mol⁻¹ s⁻¹ for the second. Comparison of these results with experimental determinations shows that hydrogen abstraction on HO₂ is an efficient mechanism for the formation of H₂ + O₂, while the concerted mechanism envisaged for the formation of H₂O + O is highly unlikely.     

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.     

Kinetic study of trifluoromethylation with s-(trifluoromethyl) dibenzothiophenium salts

The kinetic parameters were determined for C-trifluoromethylation of aniline with S-(trifluoromethyl)dibenzothiophenium triflate (1), its 3,7-dinitro derivative (2) and S-(trifluoromethyl)diphenylsulfonium triflate (3) in DMF-d₇. The higher reactivity of heterocyclic 1 compared with non-heterocyclic 3 could be explained on the basis of its greatly enhanced activation entropy (ΔS^≠: −11.2 cal mol ⁻¹ K⁻¹ for 1; −47.1 for 3), but not its enhanced activation enthalpy (ΔH^≠: 21.2 kcal mol⁻¹ for 1; 12.1 for 3). The aromatic delocalization of the heterocyclic ring may thus be only slightly disturbed by the S-trifluoromethyl substituent. The high reactivity of 2 was attributed to the great electron deficiency caused by two nitro groups in addition to the heterocyclic salt system (ΔH^≠ 17.0 kcal mol⁻¹, ΔS^≠ −9.1 cal mol⁻¹ K⁻¹ for 2). The reaction mechanism is discussed; the conventional S_N2 attack mechanism was ruled out and a mechanism for a side-on attack to the CF₃-S bond may possibly be applicable.     

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.     

Theoretical investigations of structure and enzymatic mechanisms of aspartyl proteinases : Part I. Ab-initio calculations on an active site model: hydrogen diformiate with H2O and H3O

The active site of aspartyl proteinases (Asp) was modelled as two formiates connected with a proton and set in geometry corresponding to Asp 32 and Asp 215 side chain carboxylate groups of endothiapepsin. The shared solvent molecule was alternatively H₂O and H₃O⁺. Their positions and those of hydrogen-bonded protons were optimized using the STO-3G basis set. Full geometry optimizations were made of the hydrogen diformiate complexes with H₂O and H₃O⁺. Asymmetric hydrogen-bonded structures resulted from these calculations, except for the fully optimized complex with H₂O. In the complexes with H₃O⁺, one proton moved consistently to the proximate carboxylic oxygen yielding a neutral, hydrated formic acid dimer. Interaction energies and proton potential energy curves were calculated using the 4-31G basis set. The interaction energy with H₂O was found to be 20.49 kcal mol⁻¹ and 202.75 kcal mol⁻¹ with H₃O⁺.     

Ab initio molecular orbital study on three feasible mechanisms for substitution of the vinylic carbon in F2C=C(OMs)BMe3−

A substitution on 2,2-difluorovinylic carbon was investigated by using ab initio molecular orbital calculations. Three feasible mechanisms, which are the S_N1-like, the S_N2-type and the addition-elimination mechanisms, were ex- amined for a model borate, 2,2-difluoro-1-mesyloxyvinyl(trimethyl)borate. Four TSs were obtained depending on the position of Li⁺ around the vinylborate although activation energies in the gas phase are rather high (ca. 30–40 kcal mol⁻¹) in comparison with that expected from the experimental conditions. It was confirmed at the SCRF-IPCM calculations that the solvent effect reduces the acti- vation energy of one S_N2-type mechanism very much (4. l kcal mol⁻¹ at the B3LYP/6-31+G^*//RHF/6-31+G^/s* level of theory) while those for the other mechanisms do not change very much. Therefore, the S_N2-type mechanism is applicable to the substitution reaction observed for the vinylborate.     

Correlation between hydration and deprotonation in hexaaqua metal complexes

The hydration energy of metallic cations determined with density functional calculations using a double-numerical plus p-polarization basis set, related to the acidity constants of hexaaqua metal complexes, was investigated in the present study. From the results calculated by Vosko-Wilk-Nusair (VWN), Becke-Perdew (BP) and Becke-Lee-Yang-Parr (BLYP) density functionals, a global linear correlation with the observed acidity constants in both main group [Mg(II), Ca(II) and Al(III)] and (post-)transition group [Mn(II), Zn(II), Cd(II), Sc(III), Cr(III), Fe(III), Ga(III) and In(III)] hexaaqua metal complexes has been established:

VWN density functional: pK_a = 16.5760 + 0.0173E_hydr kcal mol⁻¹

BP density functional: pK_a = 15.7329 + 0.0182E_hydr kcal mol⁻¹

BLYP density functional: pK_a = 15.9448 + 0.0185E_hydr kcal mol⁻¹     

Protonation of group 14 oxides: the proton affinity of SnO

Molecular orbital calculations are reported for the monoxides, XO, of group 14 elements (X = C, Si, Ge and Sn) and for both isomers, XOH⁺ and HXO⁺, of the protonated monoxides. Structure optimisation has been carried out using the Density Functional Theory employing the B3LYP procedure and at both Hartree-Fock and MP2 (full) levels, all with a variety of medium-sized Gaussian basis sets. In all XO molecules the oxygen atom is the preferred site for protonation, except when X = C where HCO⁺ is the lower energy isomer. Barriers to interconversion between the two isomers XOH⁺ and HXO⁺ are over-estimated by the Hartree-Fock calculations, but with wave functions that include electron correlations they generally fall into the range 27-44 kcal mol⁻¹. Proton affinities increase as the atomic number of X increases, and values calculated by averaging over all wave functions that include electron correlation, give the following proton affinities: for CO, 141.5; for SiO, 189.3; for GeO, 196.1; and for SnO, 215.6 (all in kcal mol⁻¹).     

Centrifugal distortion and the ring puckering vibration in the microwave spectrum of 2,3-dihydrofuran

The microwave spectrum of 2,3-dihydrofuran has been reinvestigated and measurements for the ground and first five excited states of the ring puckering vibration have been extended to higher frequencies and rotational quantum numbers in order to study the vibrational dependence of the rotational and centrifugal distortion constants. The ring puckering potential function derived by Green from the far infrared spectrum does not reproduce the vibrational dependence of the rotational constants well. A slightly different potential function is derived which gives a reasonable fit both to the far infrared spectrum and the rotational constants. This changes the barrier to ring inversion from 1.00 kJ mol⁻¹ to 1.12 kJ mol⁻¹. The vibrational dependence of the centrifugal distortion constants is accounted for satisfactorily by the theory developed by Creswell and Mills. An attempt to reproduce the vibrational dependence of the rotational and centrifugal distortion constants using the ring puckering potential function and a simple model for this vibration has very limited success.