单态卡宾与臭氧反应机理的量子化学研究

为了研究单态卡宾与臭氧反应机理,本文采用密度泛函理论Gaussian-3方法(G3B3)优化了反应物、中间体、过渡态和产物的几何构型。探讨了单态卡宾与臭氧反应可能途径,并通过频率分析对过渡态和中间体进行了验证,研究结果表明：单态卡宾与臭氧反应有两条反应通道,分别具有亲核反应和亲电反应特征,相对而言亲核反应通道较易发生,且为强放热反应。     

五氧化二氮和水反应的从头算研究

在臭氧空洞形成过程中,极地云是一个重要的条件。为了理解极地云的形成和性质,采用不同水平的从头计算方法对N~2O~5+H~2O→2HNO~3反应进行了理论研究。在QCISD(T)//MP2/6-311G(d,p)理论水平上,该气相反应的焓变△H^0=-39.5kJ·mol^-^1,自由能变化△G^0=-36.1kJ/mol,活化能E~a=88.2kJ·mol^-^1(298K)。计算的结果与五氧化二氮易潮解的性质相符。     

Ab initio CI study of the reaction between NH2 and NO

Using ab initio CI calculations we have evaluated the structural, energetic and kinetic parameters of the reaction between NH₂ and NO. In light of the results obtained, it appears that while the formation of molecular nitrogen is highly probable, the reaction pathway leading to N₂H+OH cannot be thermodynamically excluded. The kinetic model based on the RRKM and TST methods leads to a calculated rate constant at 298 K (k = 1.64×10⁻¹¹ cm³ molecule⁻¹ s⁻¹) which is comparable to that determined experimentally and which decreases with temperature in the range 200–700 K.     

Ab initio studies on mechanisms of cycloaddition reaction between dichloro-vinylidene and ethylene

Mechanisms of cycloaddition reaction between singlet dichloro-vinylidene (R1) and ethylene (R2) have been investigated with MP2 and B3LYP /6-31G* methods, including geometry optimization, vibrational analysis, and energy for the involved stationary points on the potential energy surface. CCSD(T)/6-31G* single point calculations are also applied to the geometries from both methods. CCSD(T) relative energies for the stationary points predicted by MP2 and B3LYP are in a very good agreement. Three reaction pathways are located. Along the first one, one intermediate (INT1) is firstly generated, which then rearranges into a three-membered ring compound (P1) via a small barrier of 5.4 kJ/mol; the other two paths share the other intermediate (INT2), which isomerizes into two four-membered ring compounds (P2 and P3) via a chlorine and a hydrogen transfer, respectively. The barriers of the latter two paths are significantly higher by approximately 25 kJ/mol than that of the former (27.2 and 28.8 vs 5.4 kJ/mol), the major reaction is therefore the formation of P1.     

Ab initio investigation of ammonia-borane complexes for hydrogen storage

The structural, electronic, and thermodynamic properties of ammonia-borane complexes with varying amounts of hydrogen have been characterized by first principles calculations within density functional theory. The calculated structural parameters and thermodynamic functions (free energy, enthalpy and entropy) were found to be in good agreement with experimental and quantum chemistry data for the crystals, dimers, and molecules. The authors find that zero-point energies change several H2 release reactions from endothermic to exothermic. Both the ammonia-borane polymeric and borazine-cyclotriborazane cycles show a strong exothermic decomposition character (approximately -10 kcal/mol), implying that rehydrogenation may be difficult to moderate H2 pressures. Hydrogen bonding in these systems has been characterized and they find the N-H bond to be more covalent than the more ionic B-H bond.     

Ab initio studies on the mechanism of the cycloaddition reaction between ketene and allene

The mechanism of the cycloaddition reaction between ketene and allene to form methylene–cyclobutanones has been studied theoretically by HF /3–21G and MP2 /3–21G . These two reactions are believed to be unsynchronous and concerted, taking place through the twisted transition states. Four orbitals are mainly involved in each reaction, which is a “2 × [1 + 1]”-type cycloaddition. The activated barrier for the two reactions are 27.2 and 27.1 kcal/mol, respectively, at the level of MP 2/6–31G * based on the MP 2/3–21G geometries, i.e., these two reactions are compatible. © 1994 John Wiley & Sons, Inc.     

Ab initio study on the mechanism of reaction HNCO+NH_2

Ab initio UMP2 and UQCISD(T) calculations, with 6-311G** basis sets, were performed for the titled reactions. The results show that the reactions have two product channels: NH2+ HNCO→NH3+NCO (1) and NH2+HNCO-N2H3+CO (2), where reaction (1) is a hydrogen abstraction reaction via an H-bonded complex (HBC), lowering the energy by 32.48 kJ/mol relative to reactants. The calculated QCISD(T)//MP2(full) energy barrier is 29.04 kJ/mol, which is in excellent accordance with the experimental value of 29.09 kJ/mol. In the range of reaction temperature 2300-2700 K, transition theory rate constant for reaction (1) is 1.68 × 1011- 3.29 × 1011 mL · mol-1· s-1, which is close to the experimental one of 5.0 ×1011 mL× mol-1· s-1 or less. However, reaction (2) is a stepwise reaction proceeding via two orientation modes, cis and trans, and the energy barriers for the rate-control step at our best calculations are 92.79 kJ/mol (for cis-mode) and 147.43 kJ/mol (for trans-mode), respectively, which is much higher than     

Ab initio chemical kinetics for the OH+HNCN reaction

The kinetics and mechanism of the reaction of the cyanomidyl radical (HNCN) with the hydroxyl radical (OH) have been investigated by ab initio calculations with rate constants prediction. The single and triplet potential energy surfaces of this reaction have been calculated by single-point calculations at the CCSD(T)/6-311+G(3df,2p) level based on geometries optimized at the B3LYP/6-311+G(3df,2p) and CCSD/6-311++G(d,p) levels. The rate constants for various product channels in the temperature range of 300-3000 K are predicted by variational transition-state and Rice-Ramsperger-Kassel-Marcus (RRKM) theories. The predicted total rate constants can be represented by the expressions ktotal=2.66 x 10(+2)xT-4.50 exp(-239/T) in which T=300-1000 K and 1.38x10(-20)xT2.78 exp(1578/T) cm3 molecule(-1) s(-1) where T=1000-3000 K. The branching ratios of primary channels are predicted: k1 for forming singlet HON(H)CN accounts for 0.32-0.28, and k4 for forming singlet HONCNH accounts for 0.68-0.17 in the temperature range of 300-800 K. k2+k7 for producing H2O+NCN accounts for 0.55-0.99 in the high-temperature range of 800-3000 K. The branching ratios of k3 for producing HCN+HNO, k6 for producing H2N+NCO, k8 for forming 3HN(OH)CN, k9 for producing CNOH+3NH, and k5+k10 for producing NH2+NCO are negligible. The rate constants for key individual product channels are provided in a table for different temperature and pressure conditions.     

Ab initio study of dehalohydrogenation reaction of 2-halo-2,3-dihydrophosphinine

Decomposition of 2-fluoro-2,3-dihydrophosphinine (1), 2-chloro-2,3-dihydrophosphinine (3), 2-bromo-2,3-dihydrophosphinine (5) to phosphinine was investigated using Molecular orbital and density functional theory. Study on the B3LYP/6-311+G** level of theory revealed that the required energy for the decomposition of compounds 1, 3, and 5 to phosphinine is 30.56 kcal·mol^?1, 28.23 kcal·mol^?1, and 24.03 kcal·mol^?1, respectively. HF/6-311+G**//B3LYP/6-311+G** calculated barrier height for the decomposition of compound 1, 3, and 5 to phosphinine is 57.56 kcal·mol^?1, 37.26 kcal·mol^?1, and 30.77 kcal·mol^?1, respectively. Also, MP2/6-311+G**//B3LYP/6-311+G** results indicated that the barrier height for the decomposition of compound 1, 3, and 5 to phosphinine is 46.59 kcal·mol^?1, 47.28 kcal·mol^?1, and 42.57 kcal·mol^?1, respectively. Natural bond orbital (NBO) population analysis and nuclear independent chemical shift (NICS) results showed that, reactants are non-aromatic but products of elimination reaction are aromatic, C-H and C-X bonds are broken and H-X bond is appear.     

Ab initio investigation on the reaction path and rate for the gas-phase reaction of HO + H2O <--> H2O + OH

This article describes an ab initio investigation on the potential surfaces for one of the simplest hydrogen atom abstraction reactions, that is, HO + H2O <--> H2O + OH. In accord with the findings in the previously reported theoretical investigations, two types of the hydrogen-bonding complexes [HOH--OH] and [H2O--HO] were located on the potential energy surface. The water molecule acts as a hydrogen donor in the [HOH--OH] complex, while the OH radical acts as a hydrogen donor in the [H2O--HO] complex. The energy evaluations at the MP2(FC) basis set limit, as well as those through the CBS-APNO procedure, have provided estimates for enthalpies of association for these complexes at 298 K as -2.1 approximately -2.3 and -4.1 approximately -4.3 kcal/mol, respectively. The IRC calculations have suggested that the [H2O--HO] complex should be located along the reaction coordinate for the hydrogen abstraction. Our best estimate for the classical barrier height for the hydrogen abstraction is 7.8 kcal/mol, which was obtained from the CBS-APNO energy evaluations. After fitting the CBS-APNO potential energy curve to a symmetrical Eckart function, the rate constants were calculated by using the transition state theory including the tunneling correction. Our estimates for the Arrhenius parameters in the temperature region from 300 to 420 K show quite reasonable agreement with the experimentally derived values.     

Ab initio direct classical trajectory investigation on the SN2 reaction of F- with NH2F: nonstatistical central barrier recrossing dynamics

The bimolecular nucleophilic substitution (S(N)2) reaction of F(a)(-) with NH(2)F(b) has been investigated with the ab initio direct classical trajectory method. According to our trajectory calculations, a dynamic behavior of nonstatistical central barrier recrossing is revealed. Among the 64 trajectories calculated in this work, 45 trajectories follow the dynamic reaction pathways as assumed by statistical theory and other 19 trajectories with central barrier recrossings are nonstatistical. For the nonstatistical trajectories, the central barrier recrossings may originate from the inefficient kinetic energy transfer from the intramolecular modes of the NH(2)F(a) moiety in the dynamic F(b)(-)…H-NH-F(a) complex to the intermolecular modes of the dynamic F(b)(-)…H-NH-F(a) complex on the exit-channel potential energy surface. With respect to the dynamic behavior of the nonstatistical central barrier recrossing, the statistical theories such as the Rice-Ramsperger-Kassel-Marcus and transition state theories without further corrections cannot be used to model the reaction kinetics for this S(N)2 reaction.     

Ab initio investigation of benzene clusters: molecular tailoring approach

An exhaustive study on the clusters of benzene (Bz)(n), n = 2-8, at MP2/6-31++G(??) level of theory is reported. The relative strengths of CH-π and π-π interactions in these aggregates are examined, which eventually govern the pattern of cluster formation. A linear scaling method, viz., molecular tailoring approach (MTA), is efficiently employed for studying the energetics and growth patterns of benzene clusters consisting up to eight benzene (Bz) units. Accuracy of MTA-based calculations is appraised by performing the corresponding standard calculations wherever possible, i.e., up to tetramers. For benzene tetramers, the error introduced in energy is of the order of 0.1 mH (～0.06 kcal/mol). Although for higher clusters the error may build up, further corrections based on many-body interaction energy analysis substantially reduce the error in the MTA-estimate. This is demonstrated for a prototypical case of benzene hexamer. A systematic way of building up a cluster of n monomers (n-mer) which employs molecular electrostatic potential of an (n-1)-mer is illustrated. The trends obtained using MTA method are essentially identical to those of the standard methods in terms of structure and energy. In summary, this study clearly brings out the possibility of effecting such large calculations, which are not possible conventionally, by the use of MTA without a significant loss of accuracy.     

NH2(X^2B1)与C2H4反应机理的从头算研究

用从头算MP2方法,在6-311G^**基组下,对NH2X^2B1)与C2H4的加成和氢迁移反应机理进行了研究,优化得到反应的过渡态,并通过振动分析和内禀反应坐标(IRC)加以证实,计算了两个反应的能垒和1500K～2000K温度范围内的速率常数。结果表明：在1500K～2000K温度范围内加成反应是NH2(XX^2B1)与C2H4的反应的主要通道,同时报道了两个反应沿反应路径变化信息。     

用从头算研究Shilov反应的机理

Shilov反应在CH~4活化中占有中心地位,它有氧化加成和σ迁移两种可能的机理。本文用较大基组的从头算研究了这两种机理的反应过程,认为Shilov反应应按氧化加成机理进行。     

Ab initio study of the F + CH3NHNH2 reaction mechanism

The F + CH(3)NHNH(2) reaction mechanism is studied based on ab initio quantum chemistry methods as follows: the minimum energy paths (MEPs) are computed at the UMP2/6-311++G(d,p) level; the geometries, harmonic vibrational frequencies, and energies of all stationary points are predicted at the same level of theory; further, the energies of stationary points and the points along the MEPs are refined by UCCSD(T)/6-311++g(3df,2p). The ab initio study shows that, when the F atom approaches CH(3)NHNH(2), the heavy atoms, namely N and C atoms, are the favorable combining points. For the two N atoms, two prereaction complexes with C(s) symmetry are generated and there exists seven possible subsequent reaction routes, of which routes 1, 2, 5, and 7 are the main channels. Routes 1, 2, and 5 are associated with HF elimination, with H from the amino group or imido group, and route 7 involves the N-N bond break. Routes 3 and 6 with relation to HF elimination with H from methyl, and route 4 involved the C-N bond break, are all energetically disfavored. For the C atom, the attack of F results in the break of the C-N bond and the products are CH(3)F + NHNH(2). This route is very competitive.     

Ab initio potential energy surfaces of the ion-molecule reaction: C2H2 + O+

High level ab initio calculations using complete active space self-consistent field and multi reference single and double excitation configuration interaction methods with cc-pVDZ (correlation consistent polarized valence double zeta) and cc-pVTZ (triple zeta) basis sets have been performed to elucidate the reaction mechanism of the ion-molecule reaction, C2H2(1Sigmag+) + O+(4S), for which collision experiment has been performed by Chiu et al. [J. Chem. Phys. 109, 5300 (1998)]. The minor low-energy process leading to the weak spin-forbidden product C2H2+ (2Piu) + O(1D) has been studied previously and will not be discussed here. The major pathways to form charge-transfer (CT) products, C2H2+ (2Piu) + O(3P) (CT1) and C2H2+ (4A2) + O(3P) (CT2), and the covalently bound intermediates are investigated. The approach of the oxygen atom cation to acetylene goes over an energy barrier TS1 of 29 kcal/mol (relative to the reactant) and adiabatically leads the CT2 product or a weakly bound intermediate Int1 between CT2 products. This transition state TS1 is caused by the avoided crossing between the reactant and CT2 electronic states. As the C-O distance becomes shorter beyond the above intermediate, the C1 reaction pathway is energetically more favorable than the Cs pathway and goes over the second transition state TS2 of a relative energy of 39 kcal/mol. Although this TS connects diabatically to the covalent intermediate Int2, there are many states that interact adiabatically with this diabatic state and these lead to the other charge-transfer product CT1 via either of several nonadiabatic transitions. These findings are consistent with the experiment, in which charge transfer and chemical reaction products are detected above 35 and 39 kcal/mol collision energies, respectively.     

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

Six product channels have been found in the association reaction of CN + CH₂CO, and a variety of possible complexes and saddle points along the minimum energy reaction paths have been characterized at the UMP2(full)/6‐31G(d) level. The dominant reaction channels are the production of CH₂CN + CO and CH₂NC + CO. The isomerization and dissociation reactions of the major products of CH₂CN and CH₂NC have been investigated using the G2MP2 level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006