A Gaussian-2 ab initio study of [C2H5S] ions: III. H2 and CH4 eliminations from CH3CHSH and CH3SCH2

Gaussian-2 ab initio calculations were performed to examine the six modes of unimolecular dissociation of cis-CH₃CHSH⁺ (1⁺), trans-CH₃CHSH⁺ (2⁺), and CH₃SCH₂⁺ (3⁺): 1⁺→CH₃⁺+trans-HCSH (1); 1⁺→CH₃+trans-HCSH⁺ (2); 1⁺→CH₄+HCS⁺ (3); 1⁺→H₂+c-CH₂CHS⁺ (4); 2⁺→H₂+CH₃CS⁺ (5); and 3⁺→H₂+c-CH₂CHS⁺ (6). Reactions (1) and (2) have endothermicities of 584 and 496 kJ mol⁻¹, respectively. Loss of CH₄ from 1⁺ (reaction (3)) proceeds through proton transfer from the S atom to the methyl group, followed by cleavage of the C–C bond. The reaction pathway has an energy barrier of 292 kJ mol⁻¹ and a transition state with a wide spectrum of nonclassical structures. Reaction (4) has a critical energy of 296 kJ mol⁻¹ and it also proceeds through the same proton transfer step as reaction (3), followed by elimination of H_2. Formation of CH₃CS⁺ from 2⁺ (reaction (5)) by loss of H₂ proceeds through protonation of the methine (CH) group, followed by dissociation of the H₂ moiety. Its energy barrier is 276 kJ mol⁻¹. On both the MP2/6-31G* and QCISD/6-31G* potential-energy surfaces, the H₂ 1,1-elimination from 3⁺ (reaction (6)) proceeds via a nonclassical intermediate resembling c-CH₃SCH₂⁺ and has a critical energy of 269 kJ 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.     

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.     

Thermochemistry of the higher chlorine oxides ClOx (x=3, 4) and Cl2Ox (x=3–7)†

Heats of formation for ClO₃, ClO₄, Cl₂O₃, Cl₂O₄, Cl₂O₅, Cl₂O₆ and Cl₂O₇ molecules are determined at the B3LYP, B3PW91, mPW1PW91 and B1LYP levels of the density functional theory employing a series of extended basis sets, and using Gaussian-3 model chemistries. Modified Gaussian-3 calculations, which employ accurate B3LYP/6-311+G(3d2f) molecular geometries and vibrational frequencies, were also performed. Heats of formation were calculated from both total atomization energies and isodesmic reaction schemes. The latter method in conjunction with Gaussian-3 models leads to the most reliable results. The best values at 298 K for ClO₃, ClO₄, Cl₂O₃ and Cl₂O₄ as derived from an average of G3//B3LYP and G3//B3LYP/6-311+G(3d2f) calculations are 43.1, 54.8, 31.7 and 37.4 kcal mol⁻¹. From calculations carried out at the G3(MP2)//B3LYP and G3(MP2)//B3LYP/6-311+G(3d2f) levels, heats of formation for Cl₂O₅, Cl₂O₆ and Cl₂O₇ are predicted to be 53.2, 52.2 and 61.5 kcal mol⁻¹. All best values are reproduced within 1 kcal mol⁻¹ by using mPW1PW91/6-311+G(3d2f) isodesmic energies. Enthalpy changes for relevant Cl–O bond fission reactions are reported. Comparisons with previous thermodynamics data are made.     

Quantum chemical study of the hydrogen-bonded HXeOH–H2O complex

A combined MP2 and DFT/B3LYP study of the HXeOH–H₂O complex is presented. These computational methods have been used to extract information on the structural, energetical and vibrational properties of the complex. Additionally, we have applied anharmonic vibrational calculations based on the MP2-computed intermolecular potential energy surface. Large perturbations both on the subunit structures and their fundamental vibrational modes are found upon complexation. Large changes of anharmonicity of the HXeOH subunit reflects the perturbation of the molecule's electronic structure. The computed BSSE-corrected interaction energies are −40.23 and −38.94 kJ mol⁻¹ at the CCSD(T)//MP2 and CCSD(T)//B3LYP levels of theory, respectively. The estimated deformation energy contribution to the interaction energy is about 5%, which is very large compared with classical hydrogen-bonded complexes. The topological analysis of the Electron Localization Function (ELF) was applied to study further the hydrogen-bonded interaction between the two complex partners. The obtained interaction pattern suggests that the interaction between HXeOH and H₂O is a typical hydrogen bond interaction driven mainly by electrostatic interactions.     

Structure and conformational properties of carbonylbisimidosulfuryl fluoride, O=C(N=S(O)F2)2

The molecular structure and conformational properties of O=C(N=S(O)F₂)₂ (carbonylbisimidosulfuryl fluoride) were determined by gas electron diffraction (GED) and quantumchemical calculations (HF/3-21G* and B3LYP/6-31G*). The analysis of the GED intensities resulted in a mixture of 76(12)% syn–syn and 24(12)% syn–anti conformer (ΔH⁰=H⁰(syn–anti)−H⁰(syn–syn)=1.11(32) kcal mol⁻¹) which is in agreement with the interpretation of the IR spectra (68(5)% syn–syn and 32(5)% syn–anti, ΔH⁰=0.87(11) kcal mol⁻¹). syn and anti describe the orientation of the S=N bonds relative to the C=O bond. In both conformers the S=O bonds of the two N=S(O)F₂ groups are trans to the C–N bonds. According to the theoretical calculations, structures with cis orientation of an S=O bond with respect to a C–N bond do not correspond to minima on the energy hyperface. The HF/3-21G* approximation predicts preference of the syn–anti structure (ΔE=−0.11 kcal mol⁻¹) and the B3LYP/6-31G* method results in an energy difference (ΔE=1.85 kcal mol⁻¹) which is slightly larger than the experimental values. The following geometric parameters for the O=C(N=S)₂ skeleton were derived (r_a values with 3σ uncertainties): C=O 1.193 (9) Å, C–N 1.365 (9) Å, S=N 1.466 (5) Å, O=C–N 125.1 (6)° and C–N=S 125.3 (10)°. The geometric parameters are reproduced satisfactorily by the HF/3-21G* approximation, except for the C–N=S angle which is too large by ca. 6°. The B3LYP method predicts all bonds to be too long by 0.02–0.05 Å and the C–N=S angle to be too small by ca. 4°.     

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

Direct dynamic study on the hydrogen abstraction reaction C2(Πu)+H2 → C2H+H

The ab initio direct dynamics method at the G2//UQCISD/6-311 + G(d,p) level is employed to study the hydrogen abstraction reaction C₂(³Π_u)+H₂ → C₂H+H over a wide temperature range 100–4650 K. The barrier heights obtained for the forward and reverse reactions are 7.78 and 17.53 kcal/mol, respectively. Comparing with one recent experiment, the calculated forward rate constants over the temperature range 2580–4650 K are about 4.4–13.5 times greater and show a steeper temperature-dependent effect. This indicates that further experimental investigation on this simple radical reaction may still be desired. Finally, G2//UQCISD/6-311 + G(2df,2p) calculations are performed to test the reliability of the G2//UQCISD/6-311 + G(d,p) results.     

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.     

Synthesis and enthalpy of formation of SrB4O7·3H2O

Hydrated strontium borate, SrB₄O₇·3H₂O, has been synthesized and characterized by XRD, FT-IR, DTA-TG and chemical analysis. The molar enthalpy of solution of SrB₄O₇·3H₂O in 1 mol dm⁻³ HCl(aq) was measured to be (21.15 ± 0.29) kJ mol⁻¹. With incorporation of the previously determined enthalpies of solution of Sr(OH)₂·8H₂O(s) in [HCl(aq) + H₃BO₃(aq)] and H₃BO₃ in HCl(aq), and the enthalpies of formation of H₂O(l), Sr(OH)₂·8H₂O(s) and H₃BO₃(s), the enthalpy of formation of SrB₄O₇·3H₂O was found to be −(4286.7 ± 3.3) kJ mol⁻¹.     

An ab initio potential energy surface for the υ4 vibrational mode of hydrogen peroxide

In this paper, the levels and the torsional microstates of hydrogen peroxide are determined from fully optimized ab initio calculations using a nuclear model in one dimension. Calculations have been performed at the MP2 level with the 6-311 G(2df,2pd), 6-31 1+G(2df,2pd), cc-pVTZ and AUG-cc-pVTZ basis sets including polarization orbitals and diffuse functions. The most stable conformation, calculated with the MP2/AUG-cc-pVTZ approach, is a trans–gauche conformer lying at 67.5° from the trans structure. By using the same level of calculations, the heights of the trans and cis barriers have been determined to be 386.5 and 2643.8 cm⁻¹ in a good agreement with the experimental data. The variational torsional levels split into four components by the tunnelling effect of the barriers. The splitting of the fundamental level caused by the trans barrier has been found to be 11.8683 cm⁻¹, whereas the splitting caused by the cis barrier is insignificant under n=2. Current ab initio energies confirm the experimental assignments and verify the separability of the torsion from the rest of the vibrations. However, the experimental relation of dependence on the torsion of the rotational constants cannot be reproduced in one-dimension and depends on several additional vibrational effects.     

An ab initio and density functional theory study of the dimethylzinc-hydrogen selenide adduct: (CH3)2Zn:SeH2

The molecular structure (equilibrium geometry) and binding energy of the dimethylzinc (DMZn)-hydrogen selenide (H₂Se) adduct, (CH₃)₂Zn:SeH₂, have been computed with ab initio molecular orbital and density functional theory (DFT) methods and, where possible, compared with experimental results. The structure of the precursors DMZn and H₂Se are perturbed to only a small extent upon adduct formation. (CH₃)₂Zn:SeH₂ was found to be 3 kcal mol⁻¹ less stable than the precursors at the B3LYP/6-311 + G(2d,p)//B3LYP/6-311 + G(2d,p) level of computation, indicating that the (CH₃)₂Zn:SeH₂ adduct is unlikely to be a stable gas-phase species under chemical vapour deposition conditions. Further calculations at the B3LYP/6-311 + G(2d,p)//B3LYP/6-311 + G(2d,p) level of computation suggest that the 1:2 adduct species, (CH₃)₂Zn:(SeH₂)₂, is much less stable than the 1:1 adduct and consequently the precursors by 19 kcal mol⁻¹.     

Matrix isolation and ab initio studies of the H2S---CO complex

The structure, energetics, and vibrational properties of complexes formed between H₂S and CO have been investigated by matrix isolation FTIR spectroscopy and ab initio molecular orbital theory. Two stable computational minima were found representing nearly linear hydrogen bonds between the subunits. The H₂S---CO and H₂S---OC species were calculated to be bound by 5.22 and 1.54 kJ mol⁻¹, respectively. The computational results were reproduced by experimental assignments for the carbon attached complex. The stretching vibrations of the complex subunits were found to be similarly perturbed upon complexation both experimentally and computationally.     

Low-temperature heat capacity and thermodynamic properties of crystalline [Re2(Ala)4(H2O)8](ClO4)6 (Re = Eu, Er; Ala = alanine)

[Re₂(Ala)₄(H₂O)₈](ClO₄)₆ (Re=Eu, Er; Ala=alanine) were synthesized, and the low-temperature heat capacities of the two complexes were measured with a high-precision adiabatic calorimeter over the temperature range from 80 to 370 K. For [Eu₂(Ala)₄(H₂O)₈](ClO₄)₆, two solid–solid phase transitions were found, one in the temperature range from 234.403 to 249.960 K, with peak temperature 243.050 K, the other in the range from 249.960 to 278.881 K, with peak temperature 270.155 K. For [Er₂(Ala)₄(H₂O)₈](ClO₄)₆, one solid–solid phase transition was observed in the range from 270.696 to 282.156 K, with peak temperature 278.970 K. The molar enthalpy increments, ΔH_m, and entropy increments,ΔS_m, of these phase transitions, were determined to be 455.6 J mol⁻¹, 1.87 J K⁻¹ mol⁻¹ at 243.050 K; 2277 J mol⁻¹, 8.43 J K⁻¹ mol⁻¹ at 270.155 K for [Eu₂(Ala)₄(H₂O)₈](ClO₄)₆; and 4442 J mol⁻¹, 15.92 J K⁻¹ mol⁻¹ at 278.970 K for [Er₂(Ala)₄(H₂O)₈](ClO₄)₆. Thermal decompositions of the two complexes were investigated by use of the thermogravimetric (TG) analysis. A possible mechanism for the thermal decomposition is suggested.     

Ab initio quantum chemical studies of reaction mechanism for CH2CO with NCO

An extensive quantum chemical study of the potential energy surface (PES) for all possible isomerization and dissociation reactions of CH₂CO with NCO is reported at the DFT (B3LYP/6-311++G(d,p)) and CCSD(T)/cc-pVDZ//B3LYP/6-311++G(d,p) levels of theory. For the CH₂CO+NCO reaction, the formation of CO+CH₂NCO via an addition–elimination mechanism is the dominant channel on the doublet surface. While the formation of CO+CH₂OCN via bimolecular substitution reaction is in the secondary. Meanwhile, the isomerization and dissociation reactions of the products, CH₂NCO and CH₂OCN, also have been investigated using the same theoretical approach. It can be concluded that these reaction channels are not feasible kinetically at low or fairly high temperatures. On the basis of the ab initio data, the total rate constants for the CH₂CO+NCO reaction in the T=296–560 K range have been computed using conventional transition state theory with Wigner tunneling correction and fitted by a rate expression as k=2.14×10⁻¹² (cm³ molecule⁻¹ s⁻¹) exp(654.29/T). The calculated total rate constants with Wigner tunneling correction for the CH₂CO+NCO reaction are in good agreement with the available experimental values.     

Ab initio computation of the potential energy surfaces of the water·hydrocarbon complexes H2O·C2H2, H2O·C2H4 and H2O·CH4: minimum energy structures, vibrational frequencies and hydrogen bond energies

On the basis of ab initio MP2/6–31 + + G(2d,2p) calculations, we examined the potential energy surfaces of the water·hydrocarbon complexes H₂O·CH₄, H₂O·C₂H₂ and H₂O·C₂H₂ to locate all the minimum energy structures and estimate the hydrogen bond energies and vibrational frequencies associated with the C(spⁿ)---H·O and the O---H·C(spⁿ) bonds (n = 1−3). Our calculations show that H₂O·C₂H₂, H₂O·C₂H₄ and H₂O·CH₄ have two minimum energy structures (i.e., the C---H·O and O---H·C hydrogen bond forms), but H₂O·C₂H₄ has only one when the vibrational motion is taken into account, the O---H·C hydrogen bond form. We have also computed the barrier for the interconversion from one minimum to the other. The fully optimized geometries of H₂O·CH₄, H₂O·C₂H₄ and H₂O·C₂H₂ as well as the vibrational shifts of the C---H stretching frequencies in their C---H·O hydrogen-bonded forms are in good agreement with the available experimental data. The calculated hydrogen bond energies show that the C(spⁿ---H·O bond strengths decrease in the order C(sp)---H·O>C(sp²)---H·O>C(sp³)---O>C(sp³---H·O, which is also consistent with the available experimental data.     

Reaction path dynamics and theoretical rate constants for the SiH3Cl+H→SiH2Cl+H2 reaction by ab initio direct dynamics method

Ab initio direct dynamics method has been used to study the title reaction. Electronic structure information including geometries, gradients and force constants (Hessians) are calculated at the UQCISD/6-311+G^** level. Energies along the minimum energy path are improved by a series of single-point G2//QCISD calculations. The changes of the geometries, vibratioanal frequencies, potential energies and total curvature along the reaction path are discussed. The rate constants in the temperature range 200–3000 K are calculated by canonical variational transition state theory with small-curvature tunneling correction (CVT/SCT) method. The results show that the variational effect is small and in the lower temperature range, the small curvature tunneling effect is important for the reaction.     

Some ab initio valence bond studies of the H5 anion

The results of some ab initio valence bond calculations with 1s basis sets are reported for the anion H⁻₅ as H₂H⁻H₂. Both single-zeta and double-zeta calculations have been performed, with attention given primarily to the results of the single-zeta calculations. For the latter calculations, it is shown that both polarization of the H₂ bonds, and charge transfer from the H⁻ to the H₂ contribute to the binding of the linear anion. In order that bending may occur, it is necessary in these calculations to include a 1s midbond function, whose exponent is the same as that for the 1s atomic orbital of the H⁻. The midbond function is located midway between the hydrogen atoms that are adjacent to the H⁻. With a single-zeta basis, the calculated bond angle is 165°.     

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

Hydrothermal synthesis, crystal structure and characterization of [Zn(H2O)4]2 [H2As6V15O42(H2O)]·2H2O

The compound [Zn(H₂O)₄]₂[H₂As₆V₁₅O₄₂(H₂O)]·2H₂O (1) has been synthesized and characterized by elemental analysis, IR, ESR, magnetic measurement, third-order nonlinear property study and single crystal X-ray diffraction analysis. The compound 1 crystallizes in trigonal space group R3, a=b=12.0601(17) Å, c=33.970(7) Å, γ=120°, V=4278.8(12) Å³, Z=3 and R1(wR2)=0.0512 (0.1171). The crystal structure is constructed from [H₂As₆V₁₅O₄₂(H₂O)]⁴⁻ anions and [Zn(H₂O)₄]²⁺ cations linked through hydrogen bonds into a network. The [H₂As₆V₁₅O₄₂(H₂O)]⁶⁻ cluster consists of 15 VO₅ square pyramids linked by three As₂O₅ handle-like units.