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.     

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.     

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.     

Ab initio study of the unimolecular pyrolysis mechanisms of monothioformic acid

The mechanisms of the possible unimolecular reactions occurring during the pyrolysis of the four tautomers (two conformers of thiol- and two conformers of thiono-) of monothioformic acid have been proposed and investigated by ab initio methods with STO-3G and 6-31G⁷ basis sets. The effects of valence electron correlation were included by Møller-Plesset (MP) perturbation theory to the fourth order at the 6-31G⁷ level. Our best results of the activation energies are given by MP4/6-31G⁷//HF/ 6-31G⁷ plus scaled zero-point energy. The barrier heights of the dehydrogenation (via a four-centre transition state) and dehydrogensuphidation (via a three-centre transition state) of thiolformic acid pyrolysis are 67.47 and 67.09 kcal mol⁻¹ respectively. The s-cis thionformic acid is dehydrated via a three-centre transition state. The activation energy of the process (81.18 kcal mol⁻¹) is much higher than the activation energy of the dehydrogenation of the s-trans form (68.83 kcal mol⁻¹) which is dehydrogenated via a four-centre transition state. These results suggest that in thionformic acid pyrolysis, the dehyrogenation of the s-cis form is more favourable than the dehydration of the s-trans form.     

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.     

Ab initio study on the mechanism of forming a silapolycyclic compound between silylidene and formaldehyde

The mechanism of the cycloaddition reaction of forming a silapolycyclic compound between singlet silylidene and formaldehyde has been investigated with MP2/6-31G* method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. The energies of the different conformations are calculated by CCSD(T)//MP2/6-31G* method. From the potential energy profile, it can be predicted that the cycloaddition reaction process of forming the silapolycyclic compound (P₂) for this reaction consists of four steps: (I) the two reactants first form a semi-cyclic intermediate INT1a through a barrier-free exothermic reaction of 32.5 kJ mol⁻¹; (II) this intermediate then isomerizes to an active four-membered ring intermediate INT1 via a transition state TS1a with an energy barrier of 30.8 kJ mol⁻¹; (III) INT1 further reacts with formaldehyde to form an intermediate INT2, which is also a barrier-free exothermic reaction of 30.1 kJ mol⁻¹; (IV) INT2 isomerizes to a silapolycyclic compound P₂ via a transition state TS2 with a barrier of 50.6 kJ mol⁻¹. Comparing this reaction path with other competitive reaction paths, we can see that this cycloaddition reaction has an excellent selectivity.     

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

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

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

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.     

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.     

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

Kinetic study of thermal decomposition of CoOOH

The kinetics of the thermal decomposition of CoOOH powder has been studied isothermally in a temperature range of 260—310°C in air. The reaction was found to proceed by the advance of a two-dimensional reaction interface. The kinetics results indicate that there are two phases in the decomposition in this temperature range: up to 280°C with an activation energy E₁ = 34.75 kcal mol⁻¹ and above 280°C with E₂ = 18.91 kcal mol⁻¹. A reaction mechanism is proposed to account for these observations.     

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.     

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 and DFT study on the hydrolysis of carbonyl dichloride

Hydrolysis of carbonyl dichloride or phosgene (Cl₂CO) in gas phase has been investigated at Hartree–Fock, density functional and ab initio levels of theory. The effects of basis sets on the energetics of the reaction have also been explored. Calculations reveal that initially carbonyl dichloride and water form a weak complex and this complex can react further in two ways. In Path 1, water adds on to carbonyl dichloride across carbonyl bond in a concerted fashion to give dichloromethane diol, and this diol decomposes to form chloro formic acid by syn-1,2-elimination of HCl and forms CO₂ and HCl as final products. Path 2 is the concerted addition of water across carbon chlorine bond and elimination of HCl in a single step leading to the formation of chloro formic acid directly. This second path that skips the formation of dichloromethane diol is observed to be very low lying and hence is kinetically favored. Addition of second water molecule to the reacting system is found to catalyze the reaction by stabilizing the complex, intermediate and transition states and reduces the activation energy to 24.6 kcal mol⁻¹ compared to 29.9 kcal mol⁻¹ for a single water molecule.