CH自由基和NO~2反应研究: I. 反应的热力学计算

探讨了CH自由基与NO~2反应的可能路径,通过计算确定了反应物,产物和稳定中间体的电子状态和平衡构型,并运用Gaussian-3方法和MRCISD方法对可能的反应路径进行了热力学计算。在多数情况下与实验值符合较好。对于个别与理论计算差别较大的实验值进行了评述。     

CH2ClO与NO反应机理的理论研究

采用B3LYP,MP2方法在6-31 (d,p)和6-311 G(d,p)水平研究了CH2ClO自由基与NO反应的微观机理,找到了三个可能的反应通道.并得到了各反应通道的反应物、中间体、过渡态和产物的优化构型、谐振频率.成功地解释了Wu等的实验结论.从电子密度拓扑分析的角度,讨论了化学反应过程中化学键的变化规律,为实验研究大气化学反应提供理论依据.找到了该反应的结构过渡态(结构过渡区)和能量过渡态,发现了反应热与结构过渡区之间的关系.     

Direct ab initio dynamics calculation of the reaction rates of CH3OCl with OH

A direct ab initio dynamics method was carried out for the reaction CH3OCl + OH --> products. Three abstraction channels from chlorine atom, in-plane hydrogen, and out-of-plane hydrogen atoms at the CH3 group have been found. The optimized geometries and frequencies of the stationary points and the minimum-energy paths (MEPs) were calculated at the MP2/6-311G(d,p) level. To improve the reaction enthalpy and potential barrier, single-point calculations were made at three higher levels of theory, the approximate QCISD(T)/6-311++G(3df,2pd), G3, and G3(MP2) levels. Furthermore, the rate constants for three abstraction channels were evaluated using canonical variational transition state theory (CVT) with the small-curvature tunneling correction (SCT) over a wide temperature range of 220-2000 K at above three higher theory levels, respectively. The calculated rate constants as well as branching rates are in reasonable agreement with the experimental values in the temperature region 250-341 K. The present results indicate H-abstraction especially from out-of-plane hydrogen is the main reaction pathway, while Cl-abstraction is much less competitive.     

One- and three-parameter density functional theory and high level ab initio theory study of the CH+CH disproportionation reaction

Complete basis set (CBS) ab initio computational studies were performed with the target being to explore the CH+CH potential energy surface. Several closed and open shell intermediates were located on the potential energy surface. Computed enthalpies for the branching reactions, as well as heats of formation are in excellent agreement. Although CBS computed energies are of high quality, this computational study is not capable of predicting the branching product ratio due to fact that neither the MP2 nor the 6-311G(2d,2p) basis set are sufficient to locate the reactant complexes and the transition state structures for the hydrogen and carbon transfer reactions in the reaction complexes. To properly explore the CH+CH potential energy surface a much higher ab initio theory level is required.     

über die Reaktionen von CH3OCl,CF3OCl,CF3OF und CF3OH mit dem supersauren System HF/MF5 (M = As,Sb) Darstellung und Charakterisierung von CH3OCl(H)+MF6− und CF3OCl(H)+MF6−

On the Reactions of CH₃OCl, CF₃OCl, CF₃OF, and CF₃OH with the Superacid System HF/MF₅ (M = As, Sb). Preparation and Characterization of CH₃OCl(H)⁺MF₆^? and CF₃OCl(H)⁺MF₆^? The preparation of the chlorine oxoniumsalts CH₃OCl(H)⁺MF₆^? and CF₃OCl(H)⁺MF₆^? by protonation of CH₃OCl and CF₃OCl in the superacid solution of HF/MF₅ (M = As, Sb) is described. However CF₃OF and CF₃OH have not been protonated under the same conditions. In the case of CF₃OH the formation of F₂CO · MF₅ is observed. The novel compounds are characterized by nmr- and vibrational spectroscopy.     

Vibrational mode analysis for the multichannel reaction of CH3Cl + OH

The recently presented ab initio calculations for the reaction system of CH₃Cl + OH (Dehestani and Shojaie, Int J Quantum Chem, in press) are applied to the vibrational mode analysis. Extending previous work, we use the vibrational mode analysis to elucidate the relationships of the reactants, the transition state, the intermediates (IM), and the products. The extensive investigation shows that the reaction mechanism is reliable. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Structures,Vibrational Spectra,and Relative Energetics of CH3COIO3 Isomers

The structures, vibrational spectra, relative energetics, and enthalpies of formation of CH₃COIO₃ isomers have been investigated with B3LYP, B3P86 and B3PW91 methods in conjugation with the 6‐31+G(d), 6‐311+G(d,p) and 6‐311++G(3df,3pd) basis sets. The CH₃COOIO₂ structure was found to be the most stable form among the isomers with an estimated enthalpy of formation of ?314.6 kJ·mol^?1. The enthalpies of formation for CH₃COOOOI, CH₃COOOIO and CH₃COIO₃ are ?180.7, ?184.9 and ?50.6 kJ·mol^?1, respectively. The implication of the formation of CH₃COIO₃ isomers from the atmospheric cross‐reactions of the acetylperoxy (CH₃COO₂) and iodine monoxide (IO) radicals was examined and the possible dissociation products of the most likely CH₃COIO₃ isomers were determined.     

Reaction of CH3O2 and HO2: Ab initio characterization of dimer structure and vibrational mode analysis for reaction mechanisms

All species involved in the multi‐channel reaction of CH₃O₂ with HO₂ have been investigated using density functional theory (DFT). The molecular geometries for various species are optimized employing the B3LYP method implementing the 6‐311++G** basis set. The relative energies of all species are calculated at the same level theory. The results show that there are two kinds of channels: singlet and triplet. The singlet channel involves four intermediates and six transition states. The triplet channel includes two intermediates and two transition states. There are four kinds of reaction products: CH₃OOH + ¹O₂, CH₃OH + O₃, CH₄ + 2O₂, and CH₃OOH + ³O₂. The vibrational mode analysis is used to elucidate the relationships of the intermediates, the transition states, and the products. The extensive investigation shows that the reaction mechanism is reliable. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Ab initio study of the potential energy surface and product branching ratios for the reaction of O(1D) with CH3CH2Br

The potential energy surface of O(¹D) + CH₃CH₂Br reaction has been studied using QCISD(T)/6‐311++G(d,p)//MP2/6‐311G(d,p) method. The calculations reveal an insertion‐elimination reaction mechanism of the title reaction. The insertion process has two possibilities: one is the O(¹D) inserting into C? Br bond of CH₃CH₂Br producing one energy‐rich intermediate CH₃CH₂OBr and another is the O(¹D) inserting into one of the C? H bonds of CH₃CH₂Br producing two energy‐rich intermediates, IM1 and IM2. The three intermediates subsequently decompose to various products. The calculations of the branching ratios of various products formed though the three intermediates have been carried out using RRKM theory at the collision energies of 0, 5, 10, 15, 20, 25, and 30 kcal/mol. CH₃CH₂O + Br are the main decomposition products of CH₃CH₂OBr. CH₃COH + HBr and CH₂CHOH + HBr are the main decomposition products for IM1; CH₂CHOH + HBr are the main decomposition products for IM2. As IM1 is more stable and more likely to form than CH₃CH₂OBr and IM2, CH₃COH + HBr and CH₂CHOH + HBr are probably the main products of the O(¹D) + CH₃CH₂Br reaction. Our computational results can give insight into reaction mechanism and provide probable explanations for future experiments. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

C2BH3稳定性与芳香性的从头算研究

在密度泛函理论B3LYP/6-311 G(2df)水平上,对C2BH3异构体进行结构优化和简谐振动频率计算。结果表明C2BH3基态为平面环状结构(1A1,C2V)。分子轨道分析显示基态有一个垂直于分子平面、双电子占据的π分子轨道;其三元环几何中心核独立化学位移(NICS)为较大负值,这些表明基态分子具有较强的芳香性。在相同的理论水平上,本文最后详细地分析了基态的红外振动光谱。     

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.     

An ab initio study of the fcc and hcp structures of helium

The hexagonal close packed (hcp) and face centered cubic (fcc) structures of helium are studied by using a new ab initio computational model for large complexes comprising small subsystems. The new model is formulated within the framework of the energy incremental scheme. In the calculation of intra- and intersystem energies, model systems are introduced. To each subsystem associated is a set of partner subsystems defined by a vicinity criterion. In the independent calculations of intra- and intersystem energies, the calculations are performed on model subsystems defined by the subsystems considered and their partner subsystems. A small and a large basis set are associated with each subsystem. For partner subsystems in a model system, the small basis set is adopted. By introducing a particular decomposition scheme, the intermolecular potential is written as a sum of effective one-body potentials. The binding energy per atom in an infinite crystal of atoms is the negative value of this one-body potential. The one-body potentials for hcp and fcc structures are calculated for the following nearest neighbor distances (d0): 4.6, 5.1, 5.4, 5.435, 5.5, 5.61, and 6.1 a.u. The equilibrium distance is 5.44 a.u. for both structures. The equilibrium dimer distance is 5.61 a.u. For the larger distances, i.e., d0 > 5.4 a.u., the difference of the effective one-body potentials for the two structures is less than 0.2 microE(h). However, the hcp structure has the lowest effective one-body potential for all the distances considered. For the smallest distance the difference in the effective one-body potential is 3.9 microE(h). Hence, for solid helium, i.e., helium under high pressure, the hcp structure is the preferred one. The error in the calculated effective one-body potential for the distance d0 = 5.61 a.u. is of the order of 1 microE(h) (approximately 0.5%).     

An ab initio study of the CH3I photodissociation. I. Potential energy surfaces

The multireference spin-orbit (SO) configuration interaction (CI) method in its Lambda-S contracted SO-CI version is employed to calculate two-dimensional potential energy surfaces for the ground and low-lying excited states of CH3I relevant to the photodissociation process in its A absorption band. The computed equilibrium geometry for the X A1 ground state, as well as vibrational frequencies for the nu2 umbrella and nu3 symmetric stretch modes, are found to be in good agreement with available experimental data. The 3Q0+ state converging to the excited I(2P1/2o) limit is found to possess a shallow minimum of 850 cm(-1) strongly shifted to larger internuclear distances (RC-I approximately 6.5a0) relative to the ground state. This makes a commonly employed single-exponent approximation for analysis of the CH3I fragmentation dynamics unsuitable. The 4E(3A1) state dissociating to the same atomic limit is calculated to lie too high in the Franck-Condon region to have any significant impact on the A-band absorption. The computed vertical excitation energies for the 3Q1, 3Q0+, and 1Q states indicate that the A-band spectrum must lie approximately between 33,000 and 44,300 cm(-1), i.e., between 225 and 300 nm. This result is in very good agreement with the experimental findings. The lowest Rydberg states are computed to lie at >or=49,000 cm(-1) and correspond to the ...a(1)2n3a1(6sI) leading configuration. They are responsible for the vacuum ultraviolet absorption lines found experimentally beyond the A-band spectrum at 201.1 nm (49,722 cm(-1)) and higher.     

Reaction Mechanism Investigation Using Vibrational Mode Anal—ysis for the Multichannel Reacition of CH3O＋CO

On the basis of the computed results got by the Gaussian 94 package at B3LYP/6-311 G** level,the reaction mechanism of CH3O radical with CO has been investiagted thoroughly via the vibrational model analysis ,And the relationships among the reactants,eight transition states,four intermediates and various products involved this multichannel reation are eluci-dated,The vibrational mode anaysis shows that the reaction mechanism is relialbe.     

An ab initio study of formamide

Calculated energy and molecular properties of the ground and low-energy excited states of formamide are presented at the ground state geometry. Satisfactory results are obtained except for the ¹* energy which remains too high by 1 eV (which is nevertheless a large improvement over previous calculations). The predicted triplet energies lie at 5.4 eV (³ n*) and 5.8 eV (³*).     

An ab initio study of possible pathways in the thermal decomposition of NaAlH4

Density functional theory (DFT) has been used to study the structural stability of possible intermediate alanate structures, Na₅Al₃H₁₄ and Na₂AlH₅, in the thermal decomposition of NaAlH₄. Na₅Al₃H₁₄ crystallizes in the space group P4/mnc with lattice constants , and c/a=1.52. It is shown that both Na₅Al₃H₁₄ and Na₂AlH₅ have the right thermodynamics and can fit in as an intermediate state during the thermal decomposition process of NaAlH₄. The heat of formation of Na₅Al₃H₁₄ is −60 kJ/mol H₂, which is intermediate between that of NaAlH₄ (−51 kJ/mol H₂) and Na₃AlH₆ (−69.7 kJ/mol H₂). An alternative decomposition pathway based on Na₂AlH₅ has also been discussed. Frequency analysis showed that the least energetic Na₂AlH₅ structure has imaginary frequencies, implying that it is unstable. The presence of soft phonon modes also shows that Na₅Al₃H₁₄ is mechanically metastable. These results are consistent with the notion that they are the intermediate states that lead to the formation of AlH₃. This facilitates the mass transport of aluminum atoms in the decomposition pathway of NaAlH₄.     

Interconversion behavior of the C-H bond in the CH4+ radical cation: ab initio molecular dynamics study

Thermal motion of CH4+ is investigated by performing an ab initio molecular dynamics method with the second-order M?ller-Plesset (MP2)/6-311G** force field. In the trajectories obtained at 400 K, we have observed rapid interconversion behavior of the geometrical parameters of CH4+ with the frequency of 0.6/ps, where the C-H pair forming the small angle around 55 degrees is switched to another pair on subpicosecond time scale. The switching patterns are found to be classified into the following two types. Type 1: one C-H of the small angled C-H pair is switched to one C-H of the other C-H pair. Type 2: the small angled C-H pair is switched to the other C-H pair, which has been newly observed in the present ab initio MD calculation. The four C-H bonds of CH4+ are characterized by the long and short C-H bonds in a time region of the trajectories, and also for the time-evolution of C-H bonds such interconversion behavior is observed. The switching patterns of the geometrical parameters are compared with those in the interconversion scheme between six equivalent C2v symmetry structures of CH4+ [Paddon-Row, M. N. et al., J Am Chem Soc 1985, 107, 7696]. We have also investigated the electronic energy fluctuation due to thermal motion of CH4+. The standard deviation of total electronic energy at 400 K is evaluated to be 1.2 kcal/mol.     

An ab initio molecular orbital study of the trimethylgallium-arsine adduct: (CH3)3Ga:AsH3

The structure, harmonic frequencies, and binding energy of the trimethylgallium-arsine adduct, (CH₃)₃Ga: AsH₃, have been computed using ab initio molecular orbital methods, and, where possible, compared with experimental results. The structures and frequencies of the precursors trimethylgallium and arsine are perturbed to only a small extent upon adduct formation. The binding energy of (CH₃)₃Ga: AsH₃ is found to be 5.2 kcal/mol lower than that for H₃Ga:AsH₃ at the MP2/HUZSP^*//RHF/HUZSP^* level of computation.     

An ab initio molecular orbital study of the grignard reagents CH3MgCl and [CH3MgCl]2: The Schlenk equilibrium

Ab initio molecular orbital calculations are used to study the modified Schlenk equilibrium: 2RMgCl (RMgCl)² MgR₂ + MgCl₂ Mg(Cl₂)MgR₂ with R=H and CH₃. In the absence of any solvents, calculations indicate that the formation of the various possible bridged dimers (RMgCl)₂ is substantially exothermic. However, using dimethylether as a model solvent, we show that the formation of the dimer (Me₂O)(CH₃)Mg(Cl₂)Mg(CH₃)(OMe₂) is exothermic only when entropic effects are included.     

硝基甲烷异构化反应势能面的ab initio研究

在B3LYP/6-311+ +G(2d,2p)水平上,优化得到硝基甲烷CH3NO2的10种异构体和23个异构化反应过渡态,并用G2MP2方法进行了单点能计算.根据计算得到的G2MP2相对能量,探讨了CH3NO2势能面上异构化反应的微观机理.研究表明,反应初始阶段的CH3NO2异构化过程具有较高的能垒,其中CH3NO2的两个主要异构化反应通道,即CH3NO2→CH3ONO和CH3NO2→CH2N(O)OH的活化能分别为270.3和267.8 kJ/mol,均高于CH3NO2的C-N键离解能.因而,从动力学角度考虑, CH3NO2的异构化反应较为不利.