气相离子-分子反应B2H3-+CS2——B2H3S-+CS的理论研究

用密度泛函方法B3LYP/6-311++G(d,p)和高级电子相关的偶合簇法CCSD(T)/6-311++G(d,p)研究了气相离子-分子反应B₂H₃^-+CS₂B₂H₃S^-+CS的机理.结果表明,B₂H₃最可能进攻CS₂中碳原子形成三元环中间体,随后通过氢迁移和最终消除CS的反应步骤形成硫原子转移产物H₃BBS^-+CS,反应大量放热且不需要活化能.B₂H₃直接对CS₂中硫原子进攻夺取硫原子的反应方式存在一定能垒阻碍.计算结果有助于深入了解B₂H₃,B₃H-6和B₄H₇^-等缺电子硼氢负离子的反应行为.     

Theoretical investigation on the protonation reactions and products of the stable [N,C,C,S] isomers

A theoretical study on the protonation system of [N,C,C,S], [H,N,C,C,S]+, was performed at the B3LYP/6-311++G(d,p) and CCSD(T)/6-311++G(2df,2p) (single point) levels of theory. On the doublet [H,N,C,C,S]+ surface, 24 species were located as energy minima and 10 of them were considered as kinetically stable species. The species HNCCS+ with 2A' state and a shallow W-shaped skeleton was predicted to be the global minimum and kinetically the most stable species, being in good agreement with previous experimental findings. Furthermore, the protonation reactions of the stable [N,C,C,S] isomers were investigated in detail. The calculation results indicated that the [N,C,C,S] isomers may be significantly stabilized upon protonation. Finally, the possible covalent structures of the [H,N,C,C,S]+ isomers with considerable stability were briefly discussed.     

Protonated nitro group: structure, energy and conjugation

Structure of protonated nitro compounds was investigated by calculations at the levels MP2(FC)/6-311++G(2d,2p)//MP2(FC)/6-311++G(2d,2p)(nitromethane and reference compounds) or B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p)(nitrobenzene and its 18 meta- and para-substituted derivatives). The group NO2H+ reveals many similarities with the isoelectronic group CO2H as the preferred conformation, conformational equilibrium, and stabilization by interaction (resonance) within the group quantified by means of isodesmic reactions. However, there is a difference in the interaction with donor groups (for instance in 4-nitroaniline) that is much stronger with NO2H+ than with CO2H. This interaction may be called resonance and may be described by standard resonance formulas, but these formulas predict only partially the geometry and cannot explain the great interaction energy.     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

CH3SH与NO2的反应机理及动力学

在G3B3,CCSD(T)/6-311++G(d,p)//B3LYP/6-311++G(d,p)水平上详细研究了CH3SH与基态NO2的微观反应机理.在B3LYP/6-311++G(d,p)水平得到了反应势能面上所有反应物、过渡态和产物的优化构型,通过振动频率分析和内禀反应坐标(IRC)跟踪验证了过渡态与反应物和产物的连接关系.在CCSD(T)/6-311++G(d,p)和G3B3水平计算了各物种的能量,得到了反应势能面.利用经典过渡态理论(TST)与变分过渡态理论(CVT)并结合小曲率隧道效应模型(SCT),分别计算了在200～3000K温度范围内的速率常数kTST,kCVT和kCVT/SCT.研究结果表明,该反应体系共存在5个反应通道,其中N进攻巯基上H原子生成CH3S+HNO2的通道活化势垒较低,为主要反应通道.动力学数据也表明,该通道在200～3000K计算温度范围内占绝对优势,拟合得到的速率常数表达式为k1CVT/SCT=1.93×10-16T0.21exp(-558.2/T)cm3·molecule-1·s-1.     

COS与O2反应的密度泛函理论研究

用量子化学密度泛函理论（DFT）方法，对COS与O2的反应进行了理论研究．在UB3LYP／6—31G^＊，UB3LYP／6—311＋＋G^＊＊水平上，优化了反应势能面上各驻点（反应物、产物、中间体和过渡态）的几何构型，在UB3LYP／6—31G^＊水平上通过内禀反应坐标（IRC）计算和振动分析，对过渡态进行了确认．在CCSD（T）／6—311＋＋G（2d，2p）水平上进行了单点能量计算，并确定了反应机理．研究结果表明，反应主要产物为CO2和SO.     

Theoretical study on the reaction mechanism of NH2(-) with O2 (a1Δg)

A detailed theoretical study of the potential energy surface of poorly understood ion-molecule reaction of NH(2)(-) and O(2) (a(1)Δ(g)) is explored at the density functional theory B3LYP/6-311++G(d,p), ab initio of QCISD/6-311++G(d,p) and CCSD(T)/6-311++G(3df, 2pd) (single-point) theoretical levels for the first time. It is shown that there are six total possible products from P(1) to P(6) on the singlet potential energy surface. Among these, the charge-transfer product P(1) (NH(2) + O(2)(-)) is the most favorable product with predominant abundances, whereas P(4) (NO(-) + H(2)O) and P(2) (HNO + OH(-)) may be the second and third feasible products followed by the almost neglectable P(3) (NO(2)(-) + H(2)), while P(5) (c-NO(2)(-) + H(2)) and P(6) (ONO(-) + H(2)) will not be observed due to their either high barriers or being secondary products. The present theoretical study points out that besides P(1) (NH(2) + O(2)(-)) and P(2) (HNO + OH(-)), P(4) (NO(-) + H(2)O) should be also observed, which is different from the previous experiment study by Anthony Midey et al. in 2008. In addition, almost all of the reaction pathways to products are exothermic and the reaction rate should be very fast since the reaction barriers are very low except for P(5) (c-NO(2)(-) + H(2)) which is in agreement with the measured total reaction rate constant k = 9.0 × 10(-10) cm(3)s(-1) at 300 K in the experiment study. It is expected that the present theoretical study may be helpful for the understanding of the reaction mechanism related to NHX(-), NX(2)(-), PHX(-), and PX(2)(-) (X = H, F, and Cl).     

Computational study on the aminolysis of beta-hydroxy-alpha,beta-unsaturated ester via the favorable path including the formation of alpha-oxo ketene intermediate

The possible mechanisms of the aminolysis of N-methyl-3-(methoxycarbonyl)-4-hydroxy-2-pyridone (beta-hydroxy-alpha,beta-unsaturated ester) with dimethylamine are investigated at the hybrid density functional theory B3LYP/6-31G(d,p) level in the gas phase. Single-point computations at the B3LYP/6-311++G(d,p) and the Becke88-Becke95 1-parameter model BB1K/6-311++G(d,p) levels are performed for more precise energy predictions. Solvent effects are also assessed by single-point calculations at the integral equation formalism polarized continuum model IEFPCM-B3LYP/6-311++G(d,p) and IEFPCM-BB1K/6-311++G(d,p) levels on the gas-phase optimized geometries. Three possible pathways, the concerted pathway (path A), the stepwise pathway involving tetrahedral intermediates (path B), and the stepwise pathway via alpha-oxo ketene intermediate due to the participation of beta-hydroxy (path C), are taken into account for the title reaction. Moreover, path C includes two sequential processes. The first process is to generate alpha-oxo ketene intermediate via the decomposition of N-methyl-3-(methoxycarbonyl)-4-hydroxy-2-pyridone; the second process is the addition of dimethylamine to alpha-oxo ketene intermediate. Our results indicate that path C is more favorable than paths A and B both in the gas phase and in solvent (heptane). In path C, the first process is the rate-determining step, and the second process is revealed to be a [4+2] pseudopericyclic reaction without the energy barrier. Being independent of the concentration of amine, the first process obeys the first-order rate law.     

Theoretical Study on Reaction Mechanism of Aluminum-Water System

A theoretical study on the reaction of aluminum with water in the gas phase was performed using the hybrid density functional B3LYP and QCISD(T) methods with the 6-311+G(d,p) and the 6-311++G(d,p) basis sets. The results show that there are three possible reaction pathways that involve four isomers, seven transition structures, and two possible products for the reaction of aluminum with water. The two most favorable reaction pathways were found, whose intermediates and products agreed quite well with experimental results. The enthalpy and Gibbs free energy change of the reaction between Al and H2O at 298 and 2000 K were calculated. Some results are also in good agreement with the previous calculations or experimental results.     

OH-与CH2ClF反应的阴离子产物通道

理论研究了羟基负离子(OH^-)与氟氯代甲烷(CH₂ClF)反应的阴离子产物通道. 分别在B3LYP/6-31+G(d,p)和B3LYP/6-311++G(2d,p)水平上得到反应势能面上各关键物种的优化构型, 进而计算得到谐振频率和零点能. 基于CCSD(T)/6-311+G(3df,3dp)水平的相对能量, 描述了由质子转移和双分子亲核取代(S_N2)过程生成各阴离子产物的途径. 各阴离子产物途径势垒的计算结果表明质子转移过程是实验中的主要产物通道, 与以往实验测量的结论相符. 此外, 计算还显示双分子亲核取代过程得到了非典型的阴离子产物, 其中动力学效应可能会导致F^-的生成.     

Theoretical study of the reaction mechanism of HCN+ and CH4 of relevance to Titan's ion chemistry

Titan is the largest satellite of Saturn. In its atmosphere, CH4 is the most abundant neutral after nitrogen. In this paper, the complex doublet potential-energy surface related to the reaction between HCN+ and CH4 is investigated at the B3LYP/6-311G(d,p), CCSD(T)/6-311G++(3df,2pd)(single-point), and QCISD/6-311G(d,p) computational levels. A total of seven products are located on the PES. The initial association of HCN+ with CH4 is found to be a prereaction complex 1 (HCNHCH3(+)) without barrier. Starting from 1, the most feasible pathway is the direct H-abstraction process (the internal C-H bond dissociation) leading to the product P1 (HCNH++CH3). By C-C addition, prereaction complex 1 can form intermediate 2 (HNCHCH3(+)) and then lead to the product P2 (CH3CNH++H). The rate-controlling step of this process is only 25.6 kcal/mol. It makes the Path P2 (1) R --> 1 --> TS1/2 --> 2 --> TS2/P2 --> P2 another possible way for the reaction. P3 (HCNCH3(+) + H), P5 (cNCHCH2(+) + H2), and P6 (NCCH3(+) + H2) are exothermic products, but they have higher barriers (more than 40.0 kcal/mol); P4 (H + HCN + CH3(+)) and P7 (H + H2 + HCCNH+) are endothermic products. They should be discovered under different experimental or interstellar conditions. The present study may be helpful for investigating the analogous ion-molecule reaction in Titan's atmosphere.     

OH-与CH2CIF反应的阴离子产物通道

理论研究了羟基负离子(OH-)与氟氯代甲烷(CH2CIF)反应的阴离子产物通道.分别在B3LYP/6-31+G(d,p)和B3LYP/6-311++G(2d,p)水平上得到反应势能面上各关键物种的优化构型,进而计算得到谐振频率和零点能.基于CCSD(T)/6-311+G(3df,3dp)水平的相对能量,描述了由质子转移和双分子亲核取代(SN2)过程生成各阴离子产物的途径.各阴离子产物途径势垒的计算结果表明质子转移过程是实验中的主要产物通道,与以往实验测量的结论相符.此外,计算还显示双分子亲核取代过程得到了非典型的阴离子产物,其中动力学效应可能会导致F-的生成.     

NH0-1+2-3离解能等的高级ab initio计算与评价

A recent experimental determination[1] of the dissociation energies (D0) for H2N-H, H2N+-H and H2N-H+, the ionization energies for NH3 and NH2 resulted in large deviations when compared with those of the earlier values and the QCISD(T)/6-311+G(3df,2p) ab initio calculations. We have performed some higher level ab initio calculations on these data by utilizing the Gaussian 92/DFT and Gaussian 94 pakages of programs and have assessed the available experimental values. Our calculations were carried out at the QCISD (TQ)/aug-cc-pVDZ, G2(QCI), QCISD(T)/6-311 ++G(3df,3pd) and QCISD(T)/aug-cc-pVTZ levels of theory. Geometries were optimized at both of the MP2(full)/6-31G(d) and the MP2(full)/6-31(d,p) levels, and were compared with those of the experiments if available. The MP2(full)/6-31G(d,p) tight-optimized geometries for the neutrals are closer to those of the experiments than those of the MP2 (full)/6-31G(d), and are in excellent agreement with the experimental results as shown in Table 1. In this case, we assumed that the optimized geometries for the cations would be better if p polarization functions are added to the hydrogen atoms. We firstly noted that the symmetry of the NH3+ cation was D3h, other than Cs. as reported in ref.[1]. All of the zero-point energies and the final geometries are calculated at the MP2(full)/6-31G(d,p) level of theory. We have also repeated the QCISD(T )/6-311 + G(3df,2p) calculations of ref. [1], because we could not identify their level of goemetry optimization. It is found that the total energy, -55.244 19 Hartrees, for NH2+(1A1 ) in ref.[1] might be in error. Our result is -55.336 29 Hartrees at the same level of theory. At our highest level [QCISD(T)/aug-cc-pVTZ] of calculations as shown in Table 3, the D0 (temperature at zero Kelvin) values of H2N-H, H2N+-H(3B1for NH2+ ) and H2N- H+ are 4.51, 5.49 and 8.00 eV, respectively. These data reported in re f.[1] were 4.97, 5.59 and 8.41 eV, respectively. Our result on D0(H2N-H) supports the work of ref.[2,3,5,6]. The ionization energies (IE) for NH3 and NH2 (3B1 for NH2+) at our highest level are 10.11 and 11.09 eV while in ref.[1] were 10.16 and 10.78 eV, respectively. For the latter, our result supports the experiment of ref.[3]. Our predicted D0 for HN2+-H and IE for NH2 (1A1 for each NH2+) are 6.80 and 12.39 eV, respectively. These values differ greatly from the predicted values (9.29 and 14.88 eV) of ref.[1] where the total energy of NH2+(1A1) might be in error. The D0 value for HN-H has not been found in ref.[1]. Our result supports the work of ref.[3]. We have also derived all of these values at the temperature of 298K and under the pressure of 101kPa at several levels of thoery as shown in Table 3. On examining the experiment of ref.[1] in detail, it is easy to find that all of the larger deviations might be from a too high value of the appearance potential of proton AP(H+). Indeed, ref.[1] has mentioned that the determintion of AP(H+), due to kinetic shift, would lead to a hihger value for the dissociation energy as has been pointed out by Berkowitz and Ruscic. In this work, we concluded that, besides some mistakes in the theoretical calculations of ref.[1], the dissociation energies for H2N-H and H2N-H+,the IE for NH2 (3B1 for NH2+) might also be unreliable and need to be re-examined.
     

CH2F自由基与HNCO反应机理的理论研究

The reaction mechanism of CH2F radical with HNCO was investigated by density functional theory (DFT)at the B3LYP/6-311++G(d,p) level. The geometries of the reactants, the intermediates, the transition states and the products were optimized. The transition states were verified through the vibration analysis.The relative energies were calculated at the QCISD(T)/6-311++G**//B3LYP/6-311++G(d,p) level. Seven feasible reaction pathways of the reaction were studied. The results indicate that the pathway (5) is the most favorable to occur, so it is the main pathway of the reaction.     

Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

Theoretical studies on the intramolecular hydrogen bond and tautomerism of 8-mercaptoquinoline in the gaseous phase and in solution using modern DFT methods

A theoretical quantum chemical study of the intramolecular hydrogen bonding interactions in 8-mercaptoquinoline has been carried out. Special attention has been paid to the rotation of S-H bond and intramolecular proton-transfer reactions. Therewith, the B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p), MPW1K/6-311++G(d,p), MPW1K/6-31+G(2d,2p), BH&HLYP/6-311++G(d,p), and G96LYP/6-311++G(d,p) methods have been used. By means of the Onsager and PCM reaction field methods, the effects of solvent on hydrogen-bond energies, conformational equilibria, rotational barriers, and tautomerism in aqueous solution have been studied. These simulations were done at the MPW1K/6-311++G(d,p) and B3LYP/6-311++G(d,p) levels. Natural-bond orbital analysis has been performed to study the intramolecular hydrogen bond (IHB) in the gaseous phase and in aqueous medium. The stability of forms under consideration in solution does not coincide with that in the gaseous phase, underlining a great importance of the electrostatic influence of solvent. Double-proton transfer in the prototropic tautomerization of 8-mercaptoquinoline, one water molecule complex in the gaseous phase and in solution, has been systematically studied. The double-proton transfer occurs concertedly and synchronously. The water-assisted tautomerization is kinetically less, but thermodynamically more favorable, compared to that of the single-proton transfer. As in the case with single-proton transfer, for water-assisted reaction, the tautomerization energies and barrier heights decrease with the increase in dielectric constant, which implies faster and more complete tautomerization of 8-mercaptoquinoline in a polar solvent.     

CH3CH2+O(3P)反应机理的理论研究

The reaction for CH3CH2+O(3P) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the QCISD(T)/6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2O＋CH3, CH3CHO+H and CH2CH2+OH in the reaction. For the products CH2O＋CH3 and CH3CHO+H, the major production channels are A1: (R)→IM1→TS3→(A) and B1: (R)→IM1→TS4→(B), respectively. The majority of the products CH2CH2+OH are formed via the direct abstraction channels C1 and C2: (R)→TS1(TS2)→(C). In addition, the results suggest that the barrier heights to form the CO reaction channels are very high, so the CO is not a major product in the reaction.     

H2FCS单分子分解反应的直接动力学研究

采用量子化学QCISD(T)/6-311++G(d,p)//B3LYP/6-311+G(d,p)方法研究了H2FCS单分子分解反应的微观动力学性质, 构建了反应势能剖面. 利用经典过渡态理论(TST)与变分过渡态理论(CVT)并结合小曲率隧道效应模型(SCT), 分别计算了在200～3000 K温度范围内的速率常数kTST、kCVT和kCVT/SCT. 计算结果表明, H2FCS可经过不同的反应通道生成10种小分子产物, 脱H反应和HF消去反应为标题反应的主反应通道, 其中HF消去反应产物HCS可由两条反应通道生成. 在200~3000 K温度区间内得到三条反应通道的表观反应速率常数三参数表达式分别为 , 和 . 速率常数计算结果显示, 量子力学隧道效应在低温区间对反应速率常数的影响显著, 而变分效应在计算温度范围内可以忽略.     

CH3 CH2 +N(4S)反应机理的理论研究

The reaction for CH3CH2+N(4S) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single point calculations for all the stationary points were carried out at the QCISD(T)/ 6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2CH2+3NH and H2CN＋CH3, and the minor products are the CH3CHN+H in the reaction. The majority of the products CH2CH2+3NH are formed via a direct hydrogen abstraction channel. The products H2CN＋CH3 are produced via an addition/dissociation channel. The products CH3CHN+H are produced via an addition/dissociation channel.     

Kinetics of the multichannel reaction of methanethiyl radical (CH3S*) with 3O2

The CH3S* + O2 reaction system is considered an important process in atmospheric chemistry and in combustion as a pathway for the exothermic conversion of methane-thiyl radical, CH3S*. Several density functional and ab initio computational methods are used in this study to determine thermochemical parameters, reaction paths, and kinetic barriers in the CH3S* + O2 reaction system. The data are also used to evaluate feasibility of the DFT methods for higher molecular weight oxy-sulfur hydrocarbons, where sulfur presents added complexity from its many valence states. The methods include: B3LYP/6-311++G(d,p), B3LYP/6-311++G(3df,2p), CCSD(T)/6-311G(d,p)//MP2/6-31G(d,p), B3P86/6-311G(2d,2p)//B3P86/6-31G(d), B3PW91/6-311++G(3df,2p), G3MP2, and CBS-QB3. The well depth for the CH3S* + 3O2 reaction to the syn-CH3SOO* adduct is found to be 9.7 kcal/mol. Low barrier exit channels from the syn-CH3SOO* adduct include: CH2S + HO2, (TS6, E(a) is 12.5 kcal/mol), CH3 + SO2 via CH3SO2 (TS2', E(a) is 17.8) and CH3SO + O (TS17, E(a) is 24.7) where the activation energy is relative to the syn-CH3SOO* stabilized adduct. The transition state (TS5) for formation of the CH3SOO adduct from CH3S* + O2 and the reverse dissociation of CH3SOO to CH3S* + O2 is relatively tight compared to typical association and simple bond dissociation reactions; this is a result of the very weak interaction. Reverse reaction is the dominant dissociation path due to enthalpy and entropy considerations. The rate constants from the chemical activation reaction and from the stabilized adduct to these products are estimated as functions of temperature and pressure. Our forward rate constant and CH3S loss profile are in agreement with the experiments under similar conditions. Of the methods above, the G3MP2 and CBS-QB3 composite methods are recommended for thermochemical determinations on these carbon-sulfur-oxygen systems, when they are feasible.