Theoretical study of structural properties, infrared spectra, and energetic properties of C4H4 isomers

A high level of ab initio and several density functional theory studies were performed with the target being to elucidate structural, spectroscopic, and energetic properties of triafulvene, cyclobutadiene, and tetrahedrane as three C₄H₄ isomers. For triafulvene all of the experimental values are available, but for cyclobutadiene only the IR spectra is known, while for tetrahedrane none of these physical values have been estimated experimentally. An excellent agreement between the computed and the experimental values were obtained where data are available, which makes it confident that the computed data should be very accurate. New structural data, IR spectra, and enthalpies of formation for both cyclobutadiene and tetrahedrane are offered.     

Complete basis set, the MP2 ab initio, and hybrid density functional theory evaluation of ionization potential and electron affinity for PH, PH2, PHF, PF, and PF2

A systematic study was performed on the small molecular systems built from phosphor, hydrogen and fluorine with the target being to evaluate accurately their ionization potentials and electron affinities, as well as influence fluorine on the ionization potential of phosphor as a central atom. To determine the accuracy of hybrid density functional methods for computing those energies, ionization energies for hydrogen, fluorine and phosphor were calculated and compared with the experimental and CBSQ values. To demonstrate the accuracy of this method, both the ionization potential and the electron affinity for phosphorus and fluorine atoms were calculated and compared with the experimental data. For both PF and PF₂, an identical electron affinity of 0.72 eV and for PH and PHF 1.0 eV were suggested.     

Relative energies of isomers of C36F2

The symmetry unrestricted C₃₆F₂ isomers formed from fullerene C₃₆, the initial symmetry of which is C_6v, C_6h, or D_2d, have been extensively studied with semi-empirical (AM1 and PM3) calculations. Based on the relationship between the isomer's stability and the adding positions, three patterns of the adding sites of F₂ moiety in the additive reactions have been deducted. The results of the π-orbital axis vector (POAV) analysis indicate that the chemical reactivity of C₃₆ is the result of the high strain in the C₃₆ cage. But, in order to form stable compounds, the effects, which guide the F₂ moiety to select carbon atoms in the C₃₆ cage, are dominated by the conjugate effect in C₃₆F₂ system rather than the strain release in the C₃₆ cage.     

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.     

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

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.     

Exploring the lowest energy triplet potential energy surface for cyclic C4H4 isomers with the complete basis set ab initio method. Is the transformation of triafulvene into cyclobutadiene possible in their excited states?

High level computational studies were performed with the aim being to explore the possibility of converting the experimentally available triafulvene into the hard-to-detect cyclobutadiene. The method is based on the simple approach used to excite triafulvene into the triplet state, and then through various reaction channels, come to the aromatic triplet cyclobutadiene. Triplet cyclobutadiene is only a few kcal/mol higher in energy than singlet cyclobutadiene and should be easily relaxed into the latter. Several reaction pathways that include only a concerted mechanism, as well as reaction pathways that include the radical formation–recombination were also explored. Some possible approaches for experimentally obtaining the singlet or triplet cyclobutadiene were suggested.     

Singlet and triplet potential energy surfaces of C3H2

The detailed singlet and triplet potential energy surfaces of C₃H₂ involving nine isomers and 13 transition structures are studied at the G3 level of theory. The rearrangement mechanisms and the electronic properties of various isomers in a broad energy range have been studied in both singlet and triplet states. Cyclopropenylidene and propargylene are found to be the most stable isomers in the singlet and triplet states, respectively. The singlet isomers are found to be more kinetically stable species as a result of high conversion barriers through which they pass. The calculations indicate that cyclopropyne in its triplet state is the least kinetically stable isomer. It is realized that the G3 method comprises both computational cost and accuracy and thus can be applied to investigation of potential energy surface of small molecules.     

Ab initio molecular orbital calculations of nuclear spin-spin coupling constants in PH2, PH3, PH4 and P2H4

Ab initio molecular orbital calculations of nuclear spin-spin coupling constants in PH⁻₂, PH₃, PH⁺₄ and P₂H₄ have been carried out employing SCF perturbation theory. Basis set dependence of all the four contributing terms has been studied in order to find the criterion for the selection of basis sets to be employed for computing this property. The dependence of the coupling constants of PH⁻₂ on its geometry has also been found. This study also discusses the requirement for satisfactory computation of couplings in cases where none of the coupling nuclei is a proton. It is found that bond-centred functions along with at least double zeta basis sets reproduce coupling constants quite satisfactorily. In all the cases studied, uncontracted core basis functions yield couplings which are in better agreement with experimental couplings than those obtained with contracted core functions.     

时间分辨红外发射光谱研究C2H3+NO反应体系

The elementary reaction of C2H3＋ NO has been reported for the first time in this paper.C2H3 radical was produced by laser photolysis of vinyl bromide at 248 nm. Vibrationally excited reaction products H2CO,NCO and HCN were observed. Two exothermic reaction channels leading to HCN＋ H2CO and CH3＋ NCO are identified.     

Ni_(x)Cu_(y)-B_(24)N_(28)催化CO_(2)氧化C_(3)H_(8)速控步骤的反应机理研究北大核心CSCD

采用密度泛函理论(DFT)研究了C_(3)H_(8)和CO_(2)在Ni_(x)Cu_(y)-B_(24)N_(28)(x+y=4,x=1、2、3、4)表面吸附及速控步骤反应机理.计算了C_(3)H_(8)、CO_(2)和相应中间体在Ni_(x)Cu_(y)-B_(24)N_(28)表面的吸附能以及6条可能路径下的反应热和活化能.计算结果表明,C_(3)H_(8)和CO_(2)在Ni_(x)Cu_(y)-B_(24)N_(28)表面是物理吸附,C_(3)H_(8)+CO_(2)→CH_(3)CHCH_(3)+OCOH是最有利的路径,其在不同催化剂表面的活化能顺序是NiCu_(3)-B_(24)N_(28)(1.42 eV)、Ni_(2)Cu_(2)-B_(24)N_(28)(1.57 eV)、Ni_(3)Cu-B_(24)N_(28)(1.62 eV)、Ni_(4)-B_(24)N_(28)(1.75 eV).由此可知,在Ni_(x)Cu_(y)-B_(24)N_(28)催化CO_(2)氧化C_(3)H_(8)的体系中,Cu含量直接影响其催化活性,即NiCu_(3)-B_(24)N_(28)用于催化CO_(2)氧化C_(3)H_(8)有一定优势.     

NO reduction with C3H8 on Co/Al2O3 catalyst. effect of O2 concentration

The selective catalytic reduction (SCR) of NO by propane in excess oxygen-containing gas mixture was studied on Co/Al₂O₃ catalyst. The oxygen concentration is very important for the reaction. The NO conversion to N₂ without oxygen is 3% at 800 K and when the O₂ concentration is raised up to 8% the NO conversion reaches its maximum value of 60% at 800 K. Characterization results by TPR and UV-Vis spectroscopy show that in the catalyst, species strongly interacting with tetrahedral and octahedral Co²⁺ ions in the support are present. This revised version was published online in June 2006 with corrections to the Cover Date.     

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.     

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.     

Dipole polarizability pseudospectra and C6 dispersion coefficients for H2–H2

Pseudostate decomposition of static dipole polarizabilities for ground state H₂⁺ from a Givens-Householder diagonalization of the excitation operator (

Full-size image

H" height="17" width="20">₀−E₀) over an N-term basis of appropriate symmetry allows for a rapidly convergent evaluation of C₆ dispersion coefficients for H₂⁺–H₂⁺. 27-term pseudospectra of hypergeneralized James functions with a 30-term GGJ⁺ unperturbed wavefunction with an optimized scale factor δ=0.918 at R=2₀ give C₆ and γ₆ values that are accurate to no less than nine significant figures.     

CCl2自由基与H2O分子反应动力学研究

用213 nm激光光解CCl4产生CCl2自由基，用LP LIF技术测定了室温下基态CCl2自由基与H2O分子的反应速率常数为(5.45±0.95)×10－14 cm3•molecule－1•s－1.在G2MP2理论水平上计算了CCl2＋H2O反应的最低单重态势能面，揭示了插入与加成 消除两种反应机理，得到了三个可能的产物通道:HCl＋HClCO、HCl＋trans ClCOH以及HCl＋cis ClCOH.并用RRKM TST和传统过渡态理论计算了这三个通道的分支比及其温度效应.结果说明在低温下（273 K），插入机理的产物通道的分支比远大于加成 消除机理的产物通道， HCl＋HClCO是主要产物，分支比为77.4％，其次是HCl＋cis ClCOH，分支比为22.6％.而在高温下（3000 K），加成 消除机理的反应通道大于插入机理， HCl＋trans ClCOH分支比为82.3％.     

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.     

Reactivity of TpIr(C2H4)(DMAD) with carboxylic acids. A DFT study on geometrical isomers and structural characterization

The thermally unstable adduct Tp^Me2Ir(C₂H₄)(DMAD), which was generated “in situ” by the reaction of DMAD with Tp^Me2Ir(C₂H₄)₂ (1) at low temperature, reacted with different carboxylic acids to produce the following compounds: Tp^Me2Ir(E-C(CO₂Me)CH(CO₂Me))(H₂O)(OC(O)C₆H₄R), (R = H, 2a; o-OH, 2b; o-Cl, 2c; m-Cl, 2d; o-NO₂, 2e; m-NO₂, 2f;o-Me, 2g;p-Me, 2h) and Tp^Me2Ir(E-C(CO₂Me)CH(CO₂Me))(H₂O)(OC(O)Me) 3. In the reaction of derivative 2a with Lewis bases, Tp^Me2Ir(E-C(CO₂Me)CH(CO₂Me))(L)(OC(O)C₆H₅), (L = Py, 4a; m-Br-Py, 4b; m-Cl-Py, 4c; NCMe, 5) were obtained, of which 4b and 4c were isolated as a mixture of two isomers in which the substituted pyridine ring was present at different rotational orientations. All new compounds prepared were characterized by ¹H and ¹³C{¹H} NMR spectroscopy, the structure of compounds 2d, 2h and 4a being determined by X-ray diffraction analysis. DFT was used to analyze the relative stability and the structural orientation of the isomers.     

丙烷在Ce-Fe/Al2O3/cordierite上选择性催化还原NO

利用凝胶溶胶法和浸渍法制备了Ce-Fe/Al₂O₃/cordierite催化剂,实验研究了其催化丙烷选择性还原NO的特性。结果表明,当铈的负载量从1%增加至5%时,Ce-Fe/Al₂O₃/cordierite的C₃H₈-SCR性能先增强后减弱,3.5Ce-Fe/Al₂O₃/cordierite具有最佳的脱硝性能,在有氧条件下,600℃时可实现96.5%的脱硝效率。Ce的加入能够提升Fe/Al₂O₃/cordierite催化剂的抗硫性能。烟气中通入0.02%的SO₂后,3.5Ce-Fe/Al₂O₃/cordierite催化丙烷还原NO的转化率始终维持在93%,而没有经过Ce修饰的Fe/Al₂O₃/cordierite的NO转化率从88%下降为80%左右。利用XRD、N₂吸附-脱附、SEM、H₂-TPR、吡啶吸附红外光谱等手段研究了催化剂的物理化学性质。结果表明,加入助剂铈能与Fe形成了固溶体,增加催化剂表面Lewis酸浓度和氧化还原能力,从而提高了催化丙烷还原NO的性能。过多的铈引入会减少Fe₂O₃结晶体的形成,不利于在C₃H₈-SCR反应中形成NO₂/NO₃^-物种,从而导致NO还原效率下降。     

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.