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.     

Semiempirical molecular dynamics studies of C60/C70 fullerene oxides: C60O, C60O2 and C70O

The spectral analysis indicates that all isomers of C₆₀O, C₇₀O and C₆₀O₂ have an epoxide-like structure (an oxygen atom bridging across a C–C bond). According to the geometrical structure analysis, there are two isomers of fullerene monoxide C₆₀O (the 5,6 bond and the 6,6 bond), eight isomers of fullerene monoxide C₇₀O and eight isomers of fullerene dioxide C₆₀O₂. In order to simulate the real reaction conditions at 300 K, the calculation of the different isomers of C₆₀O, C₆₀O₂ and C₇₀O fullerene oxides was carried out using the semiempirical molecular dynamics method with two different approaches: (a) consideration of the geometries and thermodynamic stabilities, and (b) consideration of the ozonolysis mechanism. According to the semiempirical molecular dynamic calculation analysis, the probable product of this ozonolysis reaction is C₆₀O with oxygen bridging over the 6–6 bond (C_2v). The most probable product in this reaction contains oxygen bridging across in the upper part of C₇₀ (6–6 bond in C₇₀O-2 or C₇₀O-4) an epoxide-like structure. C₆₀O₂-1, C₆₀O₂-3 and C₆₀O₂-5 are the most probable products for the fullerene dioxides. All of these reaction products are consistent with the experimental results. It is confirmed that the calculation results with the semiempirical molecular dynamics method are close to the experimental work. The semiempirical molecular dynamics method can offer both the reaction temperature effect by molecular dynamics and electronic structure, dipole moment by quantum chemistry calculation.     

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.     

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

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.     

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.     

Donor–acceptor interaction of fullerene C60 with triptycene in molecular complex TPC·C60

A molecular complex of fullerene C₆₀ with triptycene, TPC·C₆₀ is obtained. The complex has a three-dimensional packing of C₆₀ molecules. According to the IR spectra, the freezing of free rotation of C₆₀ molecules in the complex is maintained up to 360 K. The XP-spectra of TPC·C₆₀ show the suppression of π–π^* transitions of TPC phenylene rings. The separation of C₆₀ molecules by TPC ones in TPC·C₆₀ results in low intensity of the C₆₀ transitions in the 420–500 nm range in an optical spectrum. This absorption is assumed as that attributed to intermolecular transitions between adjacent C₆₀ molecules.     

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.     

Direct dynamics studies on the hydrogen abstraction reactions CF3O+CH4(CD4)→CF3OH(CF3OD)+CH3(CD3)

The dynamics properties of the hydrogen abstraction reaction CF₃O+CH₄→CF₃OH+CH₃ are studied by dual-level direct dynamics method. Optimization calculations are preformed by B3LYP and MP2 with the 6-311G(d,p) basis set, and the single-point calculations are done at the multi-coefficient correction method based on quadratic configuration interaction with single and double excitations (MC-QCISD) method. The rate constants are evaluated by canonical variational transition-state theory with a small-curvature tunneling correction over a wide range of temperature 200–2000 K. The agreement between theoretical and experimental rate constants is good in the measured temperature range. The calculated results show that the variational effect is small and almost neglected over the whole temperature range, whereas, the tunneling correction plays a role in the lower temperature range. The kinetic isotope effect for the reaction is ‘normal’. The value of k_H/k_D is 2.38 at room temperature and it decreases with the temperature increasing.     

An ab initio study of the Cl(P)+C2H6→C2H5+HCl abstraction reaction

Electronic energies, geometries, and harmonic vibration frequencies for the reactants, products, and transition state for the Cl(³P)+C₂H₆→C₂H₅+HCl abstraction reaction were evaluated at the HF and MP2 levels using several correlation consistent polarized-valence basis sets. Single-point calculations at PMP2, MP4, QCISD(T), and CCSD(T) levels were also carried out. The values of the forward activation energies obtained at the MP4/cc-pVTZ, QCISD(T)/cc-pVTZ, and CCSD(T)/cc-pVTZ levels using the MP2/cc-pVTZ structures are equal to −0.1, −0.4, and −0.3 kcal/mol, respectively. The experimental value is equal to 0.3±0.2 kcal/mol. We found that the MP2/aug-cc-pVTZ adiabatic vibration energy for the reaction (−2.4 kcal/mol) agrees well with the experimental value −(2.2–2.6) kcal/mol. Rate constants calculated with the zeroth-order interpolated variational transition state (IVTST-0) method are in good agreement with experiment. In general, the theoretical rate constants differ from experiment by, at most, a factor of 2.6.     

Reaction of [(C6H6)RuCl2]2 with 7,8-benzoquinoline and 8-hydroxyquinoline

The reaction of [(C₆H₆)RuCl₂]₂ with 7,8-benzoquinoline and 8-hydroxyquinoline in methanol were performed. The obtained complexes have been studied by IR, UV–VIS, ¹H and ¹³C NMR spectroscopy and X-ray crystallography. In the reaction with 8-hydroxyquinoline the arene ruthenium(II) complex oxidized to Ru(III). The electronic spectra of the obtained compounds have been calculated using the TDDFT method. Magnetic properties of [Ru(C₉H₆NO)₃] · CH₃OH complex suggest the antiferromagnetic coupling of the ruthenium centers in the crystal lattice. EPR spectrum of [Ru(C₉H₆NO)₃] · CH₃OH compound indicates single isotropic line only characteristic for Ru³⁺ with spin equal to 1/2.     

A systematic computational study on the reactions of dimethyl sulfoxide with XO (X = NO2, Cl, Br and I)

The reactions of dimethyl sulfoxide (DMSO) with XO (X = NO₂, Cl, Br, I) have been studied at CCSD/6-311G(d,p)//B3LYP/6-311G(d,p) level. Two reaction channels have been considered: (1) the oxygen-atom transfer (OAT) from XO (X = NO₂, Cl, Br, I) to DMSO and (2) the hydrogen-abstraction by XO (X = NO₂, Cl, Br, I). The reaction mechanisms of DMSO with NO₃, ClO and BrO are similar: the OAT channel is the dominant channel and DMSO₂ is the primary product; the hydrogen-abstraction channel is not likely to be competitive with the OAT channel. The DMSO + IO reaction, because two reaction channels have the overall negative reaction activation energies and Gibbs free energies, the reaction could occur on both reaction channels. Furthermore, the formation of the stable complex may vary the yield rate of DMSO₂.     

Reaction of C2HCl2＋O2： Combined TR-FTIR Spectroscopy and Electronic Structure

The product channels and mechanisms of the C2HC12＋O2 reaction are investigated by step-scan time-resolved Fourier transform infrared emission spectroscopy and the G3MP2// B3LYP/6-311G（d,p） level of electronic structure calculations. Vibrationally excited products of HCI, CO, and CO2 are observed in the IR emission spectra and the product vibrational state distribution are determined which shows that HCI and CO are vibrationally excited with the nascent average vibrational energy estimated to be 59.8 and 51.8 kJ/mol respectively. In combination with the G3MP2//B3LYP/6-311G（d,p） calculations, the reaction mechanisms have been characterized and the energetically favorable reaction pathways have been suggested.     

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.     

Theoretical studies on the BN substituted fullerenes C70−2x(BN)x (x=1–3)—isoelectronic equivalents of C70

The equilibrium structures and relative stabilities of BN-doped fullerenes C_70−2x(BN)_x (x=1–3) have been studied at the AM1 and MNDO level. The most stable isomers of C_70−2x(BN)_x have been found out and their electronic properties have been predicted. The calculation results show that the BN substituted fullerenes C_70−2x(BN)_x have considerable stabilities, though they are less stable than their all carbon analog. For C₆₈BN, the isomers whose BN is located in the most chemically active bonds of C₇₀ (namely B and A) are among the most stable species, of which B is predicted to be the ground state. The stabilities of C₆₈BN decrease and the dipole moments increase with increasing the distance between the heteroatoms. For C₆₆(BN)₂, the lowest energy species is the isomer in which the B–N–B–N bond is formed; For C₆₄(BN)₃, the most stable species should have three BN units located in the same hexagon to form B–N–B–N–B–N ring. The ionization potentials and the affinity energies of the most stable species of BN-doped C₇₀ are almost the same as those of C₇₀ because of the isoelectronic relationship. The ionization potentials and affinity energies depend on the relative position of the heteroatoms in C₆₈BN, the chemical reactivities of the isomers whose heteroatoms are well separated should differ significantly from their all carbon analog.     

Homo- and hetero-bridged mixed-valence dinuclear complexes which contain the fragment ‘Rh(C6F5)3’

The reaction of the anionic mononuclear rhodium complex [Rh(C₆F₅)₃Cl(Hpz)]^t- (Hpz = pyrazole, C₃H₄N₂) with methoxo or acetylacetonate complexes of Rh or Ir led to the heterodinuclear anionic compounds [(C₆F₅)₃Rh(μ-Cl)(μ-pz)M(L₂)] [M = Rh, L₂ = cyclo-octa-1,5-diene, COD (1), tetrafluorobenzobarrelene, TFB (2) or (CO)₂ (4); M = Ir, L₂ = COD (3)]. The complex [Rh(C₆F₅)₃(Hbim)]⁻ (5) has been prepared by treating [Rh(C₆F₅)₃(acac)]⁻ with H₂bim (acac = acetylacetonate; H₂bim = 2,2′-biimidazole). Complex 5 also reacts with Rh or Ir methoxo, or with Pd acetylacetonate, complexes affording the heterodinuclear complexes [(C₆F₅)₃Rh(μ-bim)M(L₂)]⁻ [M = Rh, L₂ = COD (6) or TFB (7); M = Ir, L₂ = COD (8); M = Pd, L₂ = η³-C₃H₅ (9)]. With [Rh(acac)(CO)₂], complex 5 yields the tetranuclear complex [{(C₆F₅)₃Rh(μ-bim)Rh(CO)₂}₂]₂₋. Homodinuclear Rh^III derivatives [{Rh(C₆F₅)₃}₂(μ-L)₂]^·- [L₂ = OH, pz (11); OH, S^tBu (12); OH, SPh (13); bim (14)] have been obtained by substitution of one or both hydroxo groups of the dianion [{Rh(C₆F₅)₃(μ-OH)}₂]²⁻ by the corresponding ligands. The reaction of [Rh(C₆F₅)₃(Et₂O)_x] with [PdX₂(COD)] produces neutral heterodinuclear compounds [(C₆F₅)₃Rh(μ-X)₂Pd(COD)] [X = Cl (15); Br (16)]. The anionic complexes 1–14 have been isolated as the benzyltriphenylphosphonium (PBzPh₃⁺) salts.     

Rate coefficients for the reactions of cyclohexadienyl (c-C6H7) radicals with O2 and NO at room temperature

Rate coefficients for the reactions of cyclohexadienyl (c-C₆H₇) radicals with O₂ and NO were measured at 296 ± 2 K. The c-C₆H₇ radicals were detected selectively by laser-induced fluorescence. The rate coefficient for the reaction of c-C₆H₇ with O₂, (4.4 ± 0.5) × 10⁻¹⁴ cm³ molecule⁻¹ s⁻¹, was independent of the bath-gas (He) pressure (13–80 Torr). In the reaction of c-C₆H₇ with NO, thermal equilibrium among c-C₆H₇, NO, and C₆H₇NO was observed. The forward and reverse reactions were in the falloff region, and the equilibrium constant was (1.5 ± 0.6) × 10⁻¹⁵ cm³ molecule⁻¹.     

Synthesis and characterization of hydroxo, pyrazolato and carboxylato derivatives of the PdR(PPh3)moiety (R = C6F5 or C6Cl5)

The hydroxo-complexes [{PdR(PPh₃)(μ-OH)}₂] (R = C₆F₅ or C₆Cl₅) have been obtained by reaction of the corresponding [{PdR(PPh₃)(μ-Cl)}2] complexes with NBu₄OH in acetone. In this solvent, the reaction of the hydroxo-bridged complexes with pyrazole (Hpz) and 3,5-dimethylpyrazole (Hdmpz) in 1:2 molar ratio leads to the formation of the new complexes [{Pd(C₅F₅)(PPh₃)(μ-azolate)}2] and [{Pd(C₆Cl₅)(PPh₃)}₂(μ-OH)(μ-azolate)] (azolate = pz or dmpz). The reaction of the bis(μ-hydroxo) complexes with Hpz and Hdmpz in acetone in 1:1 molar ratio has also been studied, and the resulting product depends on the organic radical (C₆F₅ or C₆Cl₅) as well as the azolate (pz or dmpz). The identity of the isomer obtained has been established in every case by NMR (¹H, ¹⁹F and ³¹P) spectroscopy. The reaction of the bis(μ-hydroxo) complexes with oxalic (H₂Ox) and acetic (HOAc) acids yields the binucle ar complexes [{PdR(PPh₃)}2(μ-Ox)] (R = C₆F₅ or C₆Cl₅) and [{Pd(C₆F₅)(PPh₃)(μ-OAc)}2], respectively. [{Pd(C₆F₅)(PPh₃)(μ-OH)}₂] reacts with PPh₃ in acetone in 1:2 ratio giving the mononuclear complex trans-[Pd(C₆F₅) (OH)(PPh₃)₂], whereas the pentachlorophenylhydroxo complex does not react with PPh₃, even under forcing conditions.     

Theoretical investigation of adenine–BX3 (X=F,Cl) complex

Geometries and binding energies are predicted at B3LYP/6-311+G* level for the adenine–BX₃ (X=F,Cl) systems and four conformers with no imaginary frequencies have been obtained for both adenine–BF₃ and adenine–BCl₃, respectively, and single energy calculations using much larger basis sets (6-311+G(2df,p)) and aug-cc-pVDZ were carried out as well. The most stable conformer is BF₃ or BCl₃ connected to N₃ of adenine and with the stabilization energy of 22.55 or 20.59 kcal/mol at B3LYP/6-311+G* level (BSSE corrected). The analyses for the combining interaction between BX₃ and adenine with natural bond orbital method (NBO) and the atom-in-molecules theory (AIM) have been performed. The results indicate that all the conformers were formed with σ–p type interactions between adenine and BX₃, in which pyridine-type nitrogen or nitrogen atom of amino group offers its lone pair electron to the empty p orbital of boron atom and the concomitances of charge transference from adenine to BX₃ were occurred. Frequency analysis suggested that the stretching vibration of BX₃ underwent a red shift in complexes. Adenine–BF₃ complex was more stable than adenine–BCl₃ although the distance of B–N is shorter in the later.     

The crystal structure of (η-C6Me6)Ti[(μ-Cl)2(AlClEt)]2 and the catalytic activity of the (C6Me6)TiAl2Cl8−xEtx(x = 0---4) complexes towards butadiene

The composition of (C₆Me₆)TiAl₂Cl_8−xEt_x complexes in (C₆Me₆)TiAl₂Cl₈ + n Et₃Al (n = 0.5-6) systems was studied by UV-Vis spectroscopy and the X-ray crystal structure of one of them, (η⁶-C₆Me₆)Ti[(μ-Cl)₂(AlClEt)]₂ (IIa-2), has been determined. The complex crystallizes in the orthorhombic space group Pna2₁ with Z = 4 and lattice parameters a 15.634(3), b 11.355(2), c 14.417(2) Å. The ethyl groups of IIa-2 reside in outer positions of aluminate ligands farther away from the C₆Me₆ ligand. The other part of the complex does not differ remarkably from structures of other (arene)Ti^II complexes. Negligible activity of (C₆Me₆)TiAl₂Cl₈ towards the butadiene cyclotrimerization is considerably increased by addition of 2.5–3.0 equivalents of Et₃Al. As follows from UV-Vis spectra, such systems contain mainly the (C₆Me₆)TiAl₂Cl₅Et₃ complex. It is suggested that the introduction of three Et substituents destabilizes the Ti-(η⁶-C₆Me₆) bond so that the replacement of hexamethylbenzene by butadiene in the first step of a catalytic cycle becomes more feasible.