A theoretical study of the reaction of SiH2 with C2H2 and C2D2

Potential energy surfaces of the reaction of SiH₂ and C₂H₂ (and C₂D₂) have been calculated by means of ab initio molecular orbital theory at the QCISD/6-311G++(2df, 2p)//MP2/6-31G(d, p) level with corrections for the triple excitations to the QCISD energies. The barrier heights for the two reaction channels of the adduct, thus calculated, were further utilized for the dynamical calculation of the rate constants in the framework of quantum statistical Rice-Ramsperger-Kassel theory. Contributions of the rate constants of the various pathways to the total rate constant (K_T) for the disappearance of the reactants are critically examined and compared with experiment. The pressure dependence of K_T(C₂H₂) is primarily due to the formation of silirene. K_T(C₂D₂) is consistently higher than K_T(C₂H₂). The standard heat of formation of silirene is predicted to be 72.1 ± 3 kcal/mol. Rearrangement of silirene to vinylsilylene requires an activation energy smaller than that to silylacetylene.     

Ab initio MO study on the global minimum structure of C3H3 anion: propynl-1-yl versus allenyl anions

The geometrical structures of the C₃H⁻₃ anion are surveyed at the coupled-cluster doubles (CCD) level of theory with the aug-cc-pVDZ basis set. To clarify the CCD geometries, the stable two isomers -- propynl-l-yl 1 and allenyl 2 anions -- are further optimized at the coupled-cluster singles, doubles (triples) (CCSD(T)) level of theory both with the aug-cc-pVDZ and aug-cc-pVTZ basis sets. The final energies are calculated at the CCSD(T) and the complete active space self-consistent field (CASSCF) multi-reference internally contracted CI (MRCI) levels of theory with the aug-cc-pVTZ basis set. At the MRCI level of theory including both the corrections due to the cluster energies (MRCI+Q) and the zero-point vibrational energies, the allenyl anion 2 is about 1.3 kcal mol⁻¹ lower in energy than the propynl-l-yl anion 1. These results contrast with the previous theoretical estimates, where the propynl-l-yl anion 1 is 2-3 kcal mol⁻¹ lower in energy than the allenyl anion 2. The activation energies of the intramolecular hydrogen transfer in the 1 → 2 conversion reactions are 63.5 kcal mol⁻¹ at the MRCI+Q level of theory with the aug-cc-pVTZ basis set including the zero-point energy corrections. The adiabatic electron affinity of the planer propargyl (H₂CCCH) radical, which is the global minimum of the C₃H₃ radical, is calculated to be 0.976 eV (after correction for the zero-point energy changes) at the CCSD(T) level of theory with the aug-cc-pVTZ basis set. The present electron affinity is in fairly good agreement with the experimental one (0.893 eV) observed by Oakes and Ellison.     

Theoretical calculations of energetics, structures, and rate constants for the H + CH3OH hydrogen abstraction reactions

Barrier heights, structures, reaction energies, and rate constants are calculated with the DFT, MP2, and CCSD(T) methods for the first two channels of the H + CH₃OH reactions. The best estimate results based on CCSD(T) calculations give reaction enthalpies () for the first (−8.4 kcal/mol) and second (0.8 kcal/mol) reactive channels, which are comparable to the experimental values, −8.8 ± 0.9 and −0.3 ± 0.9 kcal/mol, respectively. Rate constants and activation energies calculated with the IVTST-0 method using CCSD(T)/cc-pVTZ geometries and frequencies are also in excellent agreement with experiment data.     

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.     

Synthesis and crystal structures of (C8h8)Nd(C5H9C5H4) (THF)2 and [(C8H8)Gd(C5H9C4H4(THF)][(C8H8)GD(C5H9C5H4(THF)2]

LnCl₃ (Ln=Nd, Gd) reacts with C₅H₉C₅H₄Na (or K₂C₈H₈) in THF (C₅H₉C₅H₄ = cyclopentylcyclopentadienyl) in the ratio of 1 : to give (C₅H₉C₅H₄)LnCl₂(THF)_n (orC₈H₈)LnCl₂(THF)_n], which further reacts with K₂C₈H₈ (or C₅H₉C₅H₄Na) in THF to form the litle complexes. If Ln=Nd the complex (C₈H₈)Nd(C₅H₉C₅H₄)(THF)₂ (a) was obtained: when Ln=Gd the 1 : 1 complex [(C₈H₈)Gd(C_%H₉)(THF)][(C₈H₈)Gd(C₅H₉H₄)(THF)₂] (b) was obtained in crystalline form.

The crystal structure analysis shows that in (C₈H₈)Ln(C₅H₉C₅H₄)(THF)₂ (Ln=Nd or Gd), the Cyclopentylcyclopentadieny (η⁵), cyclooctatetraenyl (η⁸) and two oxygen atoms from THF are coordinated to Nd³⁺ (or Gd³⁺) with coordination number 10.

The centroid of the cyclopentadienyl ring (Cp′) in C₅H₉C₅H₄ group, cyclooctatetraenyl centroid (COTL) and two oxygens (THF) form a twisted tetrahedron around Nd³⁺ (or Gd³⁺). In (C₈H₈)Gd(C₅H₉C₅H₄)(THF), the cyclopentyl-cyclopentadienyl (η⁵), cyclooctatetraenyl (η⁸) and one oxygen atom are coordinated to Gd³⁺ with the coordination number of 9 and Cp′, COT and oxygen atom form a triangular plane around Gd³⁺, which is almost in the plane (dev. -0.0144 Å).     

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.     

Pathways for the thermally induced dehydrogenation of C60H2

The thermolysis of C₆₀H₂ to yield C₆₀ and H₂ was studied by hybrid density functional theory (B3LYP/6-311G**//B3LYP/3-21G). The concerted loss of dihydrogen requires an activation energy of 92 kcalmol⁻¹ atT=452 K. An alternative radical mechanism, which is first order in the C₆₀H₂ concentration, has an activation energy at 452 K of only 61 kcalmol⁻¹. Monitoring of the C₆₀H₂ decomposition in 1,2-dichloro-[D₄]-benzene solution by NMR spectroscopy indicates a pseudo first-order reaction with an activation energy of 61.38±2.35 kcalmol⁻¹.     

Geometry and nature of the binding of the pre-reactive complex C2H4…ClF from its rotational spectrum

The encounter complex C₂H₄…ClF was isolated by using a fast-mixing nozzle before chemical reaction could occur between the components and was characterised by Fourier-transform microwave spectroscopy. Rotational constants, centrifugal distortion constants, Cl nuclear quadrupole constants and Cl spin-rotation constants were determined for the isotopomers C₂H₄…³⁵ClF and C₂H₄…³⁷ClF. The complex has C_2v symmetry with the ClF subunit perpendicular to the plane of C₂H₄ and oriented so that Cl is closer to C₂H₄. Both the centrifugal distortion constant Δ_J and the Cl nuclear quadrupole coupling constants indicate that the complex is relatively weakly bound and it is concluded that the interaction between the subunits is largely electrostatic in origin.     

Spectres infrarouges et structure des anions 2-methylallyle, CH3-C(CH2)2- et triméthylèneméthane, C(CH2)3

The infrared spectra of solid samples of C₄H₇K and C₄D₇K have been investigated in the 4000 to 30 cm⁻¹ range. A complete assignment of intramolecular fundamentals of C₄H₇⁻ and C₄D₇⁻ ions and of potassium-allyl vibrations is proposed and the intramolecular force constants are calculated. The C(CH₂)₃²⁻ anion has been identified spectroscopically. Structures of C₃H₅⁻, C₄H₇⁻ and C(CH₃)₃²⁻ are discussed and compared with those optimised by the MINDO/3 method.     

Theoretical studies and rate constant calculations of the reactions C2F5CHO with OH radicals and Cl atoms

The hydrogen abstraction reactions of C₂F₅CHO with OH radicals and Cl atoms have been investigated theoretically by a dual-level direct dynamics method. In this study, the optimized geometries and frequencies of the stationary points are calculated at the MP2/cc-pVDZ level of theory. The energies of the stationery points and the selected points along the minimum energy paths are further refined at the MC-QCISD level using the MP2 geometries. Complexes with energies less than those of the reactants or products are located at the entrance or the exit channels of the two reactions. This result indicates that both of reactions proceed via indirect reaction mechanisms. The enthalpies of formation for the reactant C₂F₅CHO and the product radical C₂F₅CO are estimated by isodesmic reactions at the MC-QCISD//MP2/cc-pVDZ level. At the same level, the rate constants are calculated by canonical variational transition state theory (CVT) incorporating with the small-curvature tunneling correction (SCT) in the temperature range 200–1000 K. Good agreement between the calculated and experimental rate constants is obtained at the room temperature. Due to the lack of the kinetic data of these reactions, the fitted three-parameter expressions based on the CVT/SCT rate constants within 200–1000 K are k₁ = 1.64 × 10⁻²⁴ T^4.33 exp (−566.1/T) and k₂ = 6.33 × 10⁻¹⁵ T^1.35 exp (550.3/T) cm³ molecule⁻¹ s⁻¹, respectively.     

Theoretical study on the hydrogen abstraction reaction of CH3CH2F (HFC-161) with Cl atom

A direct dynamics method is employed to study the hydrogen abstraction reaction of CH₃CH₂F+Cl. Three distinct transition states are located, one for -H abstraction and two for β-H abstraction. The potential energy surface (PES) information is obtained at the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) and G2//MP2/6-311G(d,p) level. Based on the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) results, the rate constants of the three reaction channels are evaluated by using the canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions over the temperature range of 220–2800 K. The calculated results indicate that -H abstraction dominates the total reaction almost over the whole temperature range.     

Dissociation reactions of the vinyl radical in the AA″ state

In the present theoretical work we have explored mechanisms of dissociation reactions of the vinyl radical in the A²A″ state (C₂H₃ (A²A″)) and examined possible pathways for nonadiabatic dissociation of C₂H₃ (A²A″) into C₂H₂ (X¹Σ_g⁺). In the calculations we used the complete active space self-consistent field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) methods in conjunction with the cc-pVDZ and cc-pVTZ basis sets. Mechanisms for the following three dissociation channels of C₂H₃ in the A²A″ state were explored: (1) C₂H₃ (A²A″) → C₂H₂ (trans, ³A_u) + H, (2) C₂H₃ (A²A″) → C₂H₂ (cis, ³A₂) + H, and (3) C₂H₃ (A²A″) → H₂CC (³A₂) + H. The CASSCF and CASPT2 potential energy curve calculations for the C₂H₃ (A²A″) dissociation channels (1)–(3) indicate that there is neither transition state nor intermediate for each of the channels. At the CASPT2//CASSCF/cc-pVTZ level, the dissociation energies for channels (1)–(3) are predicted to be 84.3, 91.1, and 86.9 kcal/mol, respectively. For a recently observed nonadiabatic dissociation of C₂H₃ (A²A″) into C₂H₂ (X¹Σ_g⁺) + H [J. Chem. Phys. 111 (1999) 3783], two previously suggested internal conversion (IC) pathways were examined based on our CASSCF and CASPT2 calculations. Our preliminary CASSCF and CASPT2 calculations indicate that the assumed IC pathway via the twisted C₂H₃ (A²A^′) structure might be feasible. The CASSCF/cc-pVTZ geometry optimization and frequency analysis calculations were performed for the four C_2v bridge structures in the ²B₂, ²A₂, ²B₁, and ²A₁ states along the pathways of the 1²A^′ (X²A^′), 1²A″ (A²A″), 2²A″, and 2²A^′ states of C₂H₃, respectively, and the CASPT2//CASSCF/cc-pVTZ energetic results indicate that the assumed IC pathway, via a C_2v (²A₂) structure and then ²A₂/²A₁ surface crossing, be not feasible since at their excitation wavelengths (327.4 and 366.2 nm) the C_2v (²A₂) structure could not be accessed.     

Synthesis and structural characterization of [Sn(1,2−O2C6H3NO2−4)(C4H8O)n---A Tin(II) catechol complex

[Sn(1,2−O₂C₆H₃NO₂−4)(C₄H₈O)]_n has a chain polymeric solid-state structure based on very distorted square-pyramidal SnO₅ moieties linked by oxygen bridges from the catechol ligands. It is, however, monomeric in dimethyl sulphoxide and soluble in many organic solvents.     

An ab initio and matrix isolation infrared study of the 1:1 C2H2–CHCl3 adduct

The details of weak C–Hπ interactions that control several inter and intramolecular structures have been studied experimentally and theoretically for the 1:1 C₂H₂–CHCl₃ adduct. The adduct was generated by depositing acetylene and chloroform in an argon matrix and a 1:1 complex of these species was identified using infrared spectroscopy. Formation of the adduct was evidenced by shifts in the vibrational frequencies compared to C₂H₂ and CHCl₃ species. The molecular structure, vibrational frequencies and stabilization energies of the complex were predicted at the MP2/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels. Both the computational and experimental data indicate that the C₂H₂–CHCl₃ complex has a weak hydrogen bond involving a C–Hπ interaction, where the C₂H₂ acts as a proton acceptor and the CHCl₃ as the proton donor. In addition, there also appears to be a secondary interaction between one of the chlorine atoms of CHCl₃ and a hydrogen in C₂H₂. The combination of the C–Hπ interaction and the secondary ClH interaction determines the structure and the energetics of the C₂H₂–CHCl₃ complex. In addition to the vibrational assignments for the C₂H₂–CHCl₃ complex we have also observed and assigned features owing to the proton accepting C₂H₂ submolecule in the acetylene dimer.     

The synthesis, X-ray and DFT structures of the free ansa–cyclopentadiene ligand C5H5CMe2CMe2C5H5

The title compound 2,3-dicyclopentadiene-2,3-dimethylbutane (C₅H₅CMe₂CMe₂C₅H₅) 1 shows the typical staggered conformation of a highly substituted ethane derivative with the two largest substituents (C₅H₅) adopting a trans position. The molecule shows C₂ symmetry about the central C–C bond. Due to the high substitution, the central bond of the ethane is elongated to 160.0 pm (X-ray structure analysis) while the DFT calculation finds a value of 159.2 pm.     

Peculiar basis set dependence of the energetics of C2S2H2 isomers. In search of adequate and affordable basis set for routine calculations

Quantum mechanical calculations using the Hartree-Fock (HF), post-Hartree-Fock {MP2 and CCSD(T)}and gradient-corrected hybrid density functional variant, B3LYP methodology in conjugation with basis sets like 6-31G*, 6-311G** and cc-pVTZ, were employed to critically probe the right quality of basis set and the appropriate level of theory that can be applied in modeling the organo-sulfur compounds. The relative stability ordering of 1,2-dithiete and dithioglyoxal reveal the delicate basis set dependence especially on ‘S’ atom. Single point calculations at the B3LYP, MP2 and CCSD(T) levels using a series of basis set ranging from double ζ quality all the way up to the quintuple ζ quality were done to illustrate the effect of electron correlation and the basis set dependency for these compounds. Basis set requirements are much higher than affordable for medium sized molecules as very slow convergence is seen even when the calculations are carried out with basis set as high as cc-pV5Z at the B3LYP and MP2 level of theories.     

Homochiral cyclopentane-based C1-symmetric P,P ligands C5H8(PPh2)(PR2) from C2-symmetric C5H8(PCl2)2

Treatment of 1,2-trans-C₅H₈(PCl₂)₂ with 1,2-C₂H₄(NHPr-i)₂ gave the C₂-symmetric perhydro-1,6,2,5-diazaphosphocine C₅H₈{P(Cl)N(Pr-i)CH₂}₂-cyclo, which produced dissymmetric C₅H₈(PPh₂){P[N(Pr-i)CH₂]₂-cyclo} on further reaction with PhMgBr. Cleavage of the P---N bonds with gaseous HCl afforded C₅H₈(PPh₂)(PCl₂), which was converted to C₅H₈(PPh₂){P(OPh)₂}₂ by reaction with phenol. All chiral P,P derivatives were obtained as racemates as well as resolved (1R,2R)- and (1S,2S)-enantiomers.     

Determination of the first ionization energy of furan (C4H4O) from an extrapolation of two nd Rydberg series observed in the mass-resolved (2 + 1) resonance enhanced multiphoton ionization spectrum

The first ionization energy of furan (C₄H₄O) has been determined from a short extrapolation of two nd (n=6–22) Rydberg series observed in the mass-resolved (2 + 1) resonance enhanced multiphoton ionization spectrum as IE=71673 ± 3 cm⁻¹. This value confirms the higher of the two values in the literature.     

Structures and relative stabilities of [C2H6N] ions: A non-empirical and MNDO study

Non-empirical molecular orbital calculations have been performed on the seven possible [C₂H₆N]⁺ isomers at different levels of sophisitcation. The results suggest that the ions possessing a resonance contribution with a filled octet are the most stable ones. Barriers to CH₃ — and NH₃ —rotation have also been computed and are generally quite low ( 5 kJ mol⁻¹).     

Structures and energetics of C3H6S isomers: a Gaussian-3 ab initio study

Twenty-two isomers/conformers of C₃H₆S^+√ radical cations have been identified and their heats of formation (ΔH_f) at 0 and 298 K have been calculated using the Gaussian-3 (G3) method. Seven of these isomers are known and their ΔH_f data are available in the literature for comparison. The least energy isomer is found to be the thioacetone radical cation (4⁺) with C_2v symmetry. In contrast, the least energy C₃H₆O^+√ isomer is the 1-propen-2-ol radical cation. The G3 ΔH_f298 of 4⁺ is calculated to be 859.4 kJ mol⁻¹, ca. 38 kJ mol⁻¹ higher than the literature value, ≤821 kJ mol⁻¹. For allyl mercaptan radical cation (7⁺), the G3 ΔH_f298 is calculated to be 927.8 kJ mol⁻¹, also not in good agreement with the experimental estimate, 956 kJ mol⁻¹. Upon examining the experimental data and carrying out further calculations, it is shown that the G3 ΔH_f298 values for 4⁺ and 7⁺ should be more reliable than the compiled values. For the five remaining cations with available experimental thermal data, the agreement between the experimental and G3 results ranges from fair to excellent.

Cation CH₃CHSCH₂^+√ (10⁺) has the least energy among the eleven distonic radical cations identified. Their ΔH_f298 values range from 918 to 1151 kJ mol⁻¹. Nevertheless, only one of them, CH₂=SCH₂CH₂^+√ (12⁺), has been observed. Its G3 ΔH_f298 value is 980.9 kJ mol⁻¹, in fair agreement with the experimental result, 990 kJ mol⁻¹.

A couple of reactions involving C₃H₆S^+√ isomers CH₂=SCH₂CH₂^+√ (12⁺) and trimethylene sulfide radical cation (13⁺) have also been studied with the G3 method and the results are consistent with experimental findings.