硝酰氯和顺/反亚硝酸氯及其互变异构反应的理论研究

用密度泛函理论（DFT）的B3LYP方法和6－311＋G（3df）基组，计算了气态下硝酰氯和顺／反应硝酸氯的几何构型、电子结构、红外光谱以及热力学性质，并讨论了它们的互变异构反应，分析了过渡态的结构。结果表明，B3LYP／6－311＋G（3df）计算得到的结果与实验值及CCSD（T）方法计算结果吻合，且更适应于研究反应机理，ClNO2转变为cis-ClONO的过渡态（TS1）偏离平面构型；cis-ClONO和trans-ClONO互变反应的过渡态（TS2）属于内旋转位垒；高水平计算表明不存在由ClNO2直接转变为trans-ClONO的过渡态，而是得到了一个十分接近异裂产物的二级马鞍点（2SP）。根据得到的热力学函数计算了气态时各温度下互变异构反应的平衡常数。     

多通道反应O(3P)＋CH2F的理论研究

用密度泛函理论研究了氧原子与氟代甲基自由基的反应.反应中出现的所有物种的平衡构型用B3LYP方法在6-311＋＋G（2d, 2p）基组水平上进行了优化，同时对各物种进行了频率分析；在同一理论水平上计算了各反应通道的势能面变化，分析了反应物、中间体、过渡态、产物的振动模式随反应途径的变化关系，阐明了该多通道反应的反应机理.     

含氧超价化合物OK_n的结构和稳定性

运用Ｂ３ＬＹＰ－ ＤＦＴ 密度函数方法,通过３ －２１Ｇ＊ 基组,对ＯＫｎ（ｎ ＝３ ～６） 系列超价化合物几何进行了优化,计算了各分子优化几何的频率．分析了各分子的稳定构型解离出Ｋ 与Ｋ２ 的反应能,从而预测了它们的热力学稳定性     

乙烯自由基与臭氧反应的DFT计算研究

采用密度泛函B3LYP/6-311G**水平计算研究了O3氧化乙烯基（C2H3）的机理,全参数优化了反应势能面上各驻点的几何构型,用内禀反应坐标（IRC）计算和频率分析方法,对过渡态进行了验证.结果表明，乙烯基（C2H3）与O3之间有很强的反应活性.     

含氧超锂化合物OLi_n的结构和稳定性从头计算研究

使用HF ,MP2 ,MP4和DFT方法 ,采用 6 31G 基组 ,对OLin 系列超价化合物几何进行了优化 .对于各分子优化的稳定构型计算它们解离出Li和Li2 的反应能 ,结果符合已有实验值 ,并预期了OLi6的解离反应能 .同时 ,还计算预言了各OLin 分子的基振动频率 .     

环丙基硅烯的理论研究环丙基硅烯C3H5SiH的重排反应及其机理

本文用限制的Ｈａｒｔｒｅｅ－Ｆｏｃｋ解析梯度方法在３－２１Ｇ和６－３１Ｇ＾＊水平上对环丙基硅烯的重排反应及其机理进行了从头算研究。以６－３１Ｇ＾＊优化构型作了二级微扰计算，并计算了各构型的频率。在此基础上得到了重排反应的热焓ΔＨ，自由能ΔＧ和平衡常数Ｋ，用Ｅｙｒｉｎｇ过渡态理论计算了反应的速度常数ｋ（Ｔ），应用Ｗｏｏｄｗａｒｄ－Ｈｏｆｆｍａｎｎ规则讨论了环丙基硅烯重排反应过程中端基的旋转机理。结果     

F-N=C→F-C≡N的动力学及产物的振动态分布

用数值方案，在RHF／3－21G分子轨道从头算法的水平上，得到了氟化异氰FNC到氟化氰FCN重排反应的反应途径（内禀反应坐标IRC）．沿着IRC；讨论了反应过程中体系几何构型的变化，计算了沿IRC运动与垂直于IRC简正振动之间的耦合常数（BK，F），各振动模式对应的频率（ωK），使用统一的半经典徽扰和无限级突然（SCP－IOS）近似理论计算了在一定能量下产物的振动分配.结果表明，在过渡态后，耦合常数（BK,F）的大小强烈地影响产物的振动态分布,另外用传统过渡态、变分过渡态理论及相关的隧道效应校正计算了该反应的速率常数.     

CH自由基和NO2反应研究Ⅱ.反应的动力学计算

采用MP2(FULL)/6-31G(d)方法从动力学计算上探讨了CH自由基与NO2反应的可能途径,找到了反应物,中间体及产物之间的能量通道和过渡态,报道了它们的构型、电子态及能量.并通过频率分析和IRC方法对所有的过渡态进行了验证.在此基础上求出了各步反应的活化能.在以前热力学研究的基础上,对于可能的反应通道进一步作了动力学分析,找到了反应主产物通道的分支比,与实验得到的分支比基本吻合.     

卤代硅烯异构化为卤代甲硅烷基硅烯的热力学与动力学分析

利用密度泛函理论(DFT) 中的B3LYP方法,全参数优化了卤代硅烯HXSi=SiXH（X=F、Cl、Br、I）异构化反应的反应物、产物及过渡态的几何构型,计算出了它们的振动频率、零点振动能(ZPVE)和总能量,并对它们进行了振动分析,以确定过渡态的真实性。又计算了反应的热力学函数变化,平衡常数及速率常数,比较了不同卤素对反应的影响。热力学与动力学计算结果表明,该异构化反应过程是一个放热的、且在低温下可自发进行的反应,但对于溴代硅烯和碘代硅烯而言,当温度达到1000 K时,反应开始转化为非自发反应。     

NH＋O3→ONH＋O2反应热力学和动力学研究

在量子化学对NH自由基与臭氧O3反应计算的基础上，应用统计热力学方法研究了100～1600 K温度范围内NH和O3反应过程的各热力学量的变化及平衡常数，用经Wigner校正的Eyring过渡态理论计算了不同温度下该反应两不同反应通道的活化热力学量、反应速率常数及频率因子.计算表明，相对于反应通道II，反应通道I不仅有很强的反应自发性，而且在动力学上也是较容易实现的反应.     

羰基氧化物RR^1COO的电子结构

本文采用6-31G基组的abinitio方法对羰基氧化物RR^1COO(R,R^1=H,F,CH~3)进行几何构型优化计算,研究其基态的电子结构。结果表明,RR^1COO的稳定结构为双自由基型,其单重态和双重态的相对稳定性受取代基的影响。H~2COO、H(CH~3)COO和(CH~3)~2COO的基态为单重态(^1A),HFCOO和F~2COO的基态为三重态(^3A),HFCOO和H(CH~3)COO的顺式结构比反式稳定。     

基态N原子与CS2反应机理的理论研究

用密度泛函理论B3LYP/6 311+G 和高级电子相关偶合簇CCSD(T)/6 311+G 方法,计算研究了四重态氮原子与二硫化碳的反应,找到了分别形成CS+NS、NCS+S和CN+S2三个反应通道,优化搜索了各反应的过渡态,并用频率分析和内禀坐标法(IRC)验证了各鞍点构型和反应路径.在三个反应通道中,反应N+CS2→CS+NS由于具有较低的活化能而容易发生,其它两个通道活化能高很难发生,计算结果与实验结果一致.同时,对反应机理进行了详细的讨论.     

Solvent effects on the absorption spectrum and recombination kinetics of diphenylcarbonyl oxide

The absorption spectra and rate constants of diphenylcarbonyl oxide recombination in a series of solvents and their binary mixtures were determined by flash photolysis. An increase in the solvent polarity causes hypsochromic shift of the maximum in the absorption spectrum of Ph₂COO. The analysis of the solvent effect on the recombination rate constant in terms of the four-parameter Koppel—Palm equation shows that the reactivity of carbonyl oxide depends on both specific and non-specific solvations. Quantum chemical B3LYP/6-31G(d) calculations of H₂COO and PhHCOO carbonyl oxides as well as the complexes of H₂COO with acetonitrile and ethylene in different media were performed using a polarized continuum model.     

碳源甲基苯热裂解机理的密度泛函动力学研究

在热力学研究的基础上,用UB3LYP/3-21G^*方法对甲苯热裂解机理进行了动力学研究。计算得到了甲苯的5种热裂解路径的活化能。用过渡状态理论,计算得到了这些路径在298～1223K温度范围内的速率常数。动力学计算结果表明：甲苯在热解温度低于963K时的主反应路径为甲苯热裂解生成苄基自由基的反应,其速控步的活化能△E~0^θ^≠=402.27kJ/mol;当温度高于963K达1223K左右时,主反应路径转为苯环上脱甲基生成苯基和甲基自由基的路径,该路径的活化能△E~0^θ^≠=456.91kJ/mol。以上研究结果与实验结果相一致。     

Dimesitylketone O-oxide: spectroscopic characterization, conformation, and reaction modes: OH formation and OH capture.

Dimesitylketone O-oxide 1b was synthesized by photolysis of dimesityldiazomethane dissolved in an oxygen saturated CCl3F solution at 140 K. Conformation and geometry of 1b were determined by comparing measured NMR chemical shifts with the corresponding chemical shifts calculated at the DFT-IGLO level of theory where it had to be considered that the molecule exists in two enantiomeric forms. Measured and calculated 1H chemical shifts agree within 0.1 ppm while the calculated 13C shift of the COO carbon (210.6 ppm) differs by only 0.4 ppm from the measured shift of 211.0 ppm. The two mesityl rings are perpendicular to each other and enclose angles of 40 and 57 degrees with the COO plane. The preferred rearrangement process of 1b is an H migration from one of the ortho-methyl groups to the terminal O atom of the COO unit. The calculated activation enthalpy of this process is 12.7 kcal/mol (B3LYP/cc-pVTZ). In contrast, the activation enthalpy for isomerization to dioxirane is 5 kcal/mol higher. In CCl3F, the activation barrier for the thermal decay was determined to be 13.8 +/- 0.2 kcal/mol and in acetonitrile 13.1 +/- 0.4 kcal/mol. H migration initiates cleavage of the OO bond and the production of an OH and a benzyl radical. Recombination of the latter in the solvent cage leads to the formation of 2-methylhydroxy-pentamethylbenzophenone, while escape of the OH radical from the solvent cage yields a ketone. These results confirm the possibility of OH production from carbonyl oxides in the solution phase.     

Energetics and structural elucidation of mechanisms for gas phase H/D exchange of protonated peptides

Hydrogen/deuterium exchange reactions involving protonated triglycine and deuterated ammonia (ND(3)) have been examined in the gas phase using a Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. Ab initio and density functional theory (DFT) calculations have been carried out to model the exchanges and to obtain energetics and vibrational frequencies for molecules involved in the proposed exchange mechanisms. Structural optimization and frequency calculations have been performed at the B3LYP level of theory with the 6-311+G(d,p) basis set. Transition states have been calculated at the same level of theory and basis set as above using the QST2 and QST3 methods. Single-point energy calculations have been performed at the MP2/6-311+G(d,p) level. Six labile sites of protonated triglycine were found to undergo H/D exchange. Of these six labile hydrogens, two are amide, three are ammonium, and one is carboxyl. Detailed mechanisms for each of these transfers are proposed. Qualitative onium ion and tautomer mechanisms for the exchanges of ammonium and amide hydrogens, respectively, using semiempirical calculations were suggested in previous studies by Beauchamp et al. As shown by the current ab initio and DFT calculations completed during this study, the mechanisms proposed in that study are notionally correct; however, the tautomer mechanisms are shown here to be the result of the fact that a second stable isomer of protonated triglycine exists in which the amide1 carbonyl oxygen is protonated. The exchange of the carboxyl hydrogen is found to proceed via a transition state resembling an ammonium ion interacting with a carboxylate moiety via two hydrogen bonds. The current work thus provides significant mechanistic and structural detail for a considerably more in-depth understanding of the processes involved in gas phase H/D exchange of peptides.     

Reaction mechanism between carbonyl oxide and hydroxyl radical: a theoretical study

The reaction mechanism of carbonyl oxide with hydroxyl radical was investigated by using CASSCF, B3LYP, QCISD, CASPT2, and CCSD(T) theoretical approaches with the 6-311+G(d,p), 6-311+G(2df, 2p), and aug-cc-pVTZ basis sets. This reaction involves the formation of H2CO + HO2 radical in a process that is computed to be exothermic by 57 kcal/mol. However, the reaction mechanism is very complex and begins with the formation of a pre-reactive hydrogen-bonded complex and follows by the addition of HO radical to the carbon atom of H2COO, forming the intermediate peroxy-radical H2C(OO)OH before producing formaldehyde and hydroperoxy radical. Our calculations predict that both the pre-reactive hydrogen-bonded complex and the transition state of the addition process lie energetically below the enthalpy of the separate reactants (DeltaH(298K) = -6.1 and -2.5 kcal/mol, respectively) and the formation of the H2C(OO)OH adduct is exothermic by about 74 kcal/mol. Beyond this addition process, further reaction mechanisms have also been investigated, which involve the abstraction of a hydrogen of carbonyl oxide by HO radical, but the computed activation barriers suggest that they will not contribute to the gas-phase reaction of H2COO + HO.     

Adsorption of nitrogen oxides on graphene and graphene oxides: insights from density functional calculations

The interactions of nitrogen oxides NO(x) (x = 1,2,3) and N(2)O(4) with graphene and graphene oxides (GOs) were studied by the density functional theory. Optimized geometries, binding energies, and electronic structures of the gas molecule-adsorbed graphene and GO were determined on the basis of first-principles calculations. The adsorption of nitrogen oxides on GO is generally stronger than that on graphene due to the presence of the active defect sites, such as the hydroxyl and carbonyl functional groups and the carbon atom near these groups. These active defect sites increase the binding energies and enhance charge transfers from nitrogen oxides to GO, eventually leading to the chemisorption of gas molecules and the doping character transition from acceptor to donor for NO(2) and NO. The interaction of nitrogen oxides with GO with various functional groups can result in the formation of hydrogen bonds OH???O (N) between -OH and nitrogen oxides and new weak covalent bonds C???N and C???O, as well as the H abstraction to form nitrous acid- and nitric acidlike moieties. The spin-polarized density of states reveals a strong hybridization of frontier orbitals of NO(2) and NO(3) with the electronic states around the Fermi level of GO, and gives rise to the strong acceptor doping by these molecules and remarkable charge transfers from molecules to GO, compared to NO and N(2)O(4) adsorptions on GO. The calculated results show good agreement with experimental observations.     

A computational study of the oxidation of SO2 to SO3 by gas-phase organic oxidants

We have studied the oxidation of SO(2) to SO(3) by four peroxyradicals and two carbonyl oxides (Criegee intermediates) using both density functional theory, B3LYP, and explicitly correlated coupled cluster theory, CCSD(T)-F12. All the studied peroxyradicals react very slowly with SO(2) due to energy barriers (activation energies) of around 10 kcal/mol or more. We find that water molecules are not able to catalyze these reactions. The reaction of stabilized Criegee intermediates with SO(2) is predicted to be fast, as the transition states for these oxidation reactions are below the free reactants in energy. The atmospheric relevance of these reactions depends on the lifetimes of the Criegee intermediates, which, at present, is highly uncertain.     

煤燃烧过程中砷与氮氧化物的反应机理

应用量子化学密度泛函理论B3LYP方法,研究了砷与氮氧化物(N_2O、NO_2和NO)的反应机理。全参数优化了各反应物、中间体、过渡态和产物的几何构型,通过频率分析证实中间体和过渡态的真实性,并通过内禀反应坐标(IRC)计算以进一步确定过渡态。为了得到更精确的能量信息,在B2PLYP水平下计算各结构的单点能,并通过动力学参数深入分析其反应机理。结果表明,砷与三种氮氧化物(N_2O、NO_2和NO)的反应能垒分别为78.45、2.58、155.85 k J/mol。在298-1800 K,各反应速率随温度的升高而增大。由于砷与NO_2的反应能垒较低,其反应速率大于1012cm3/(mol·s),说明该反应容易发生且速率极快。砷与N_2O和NO的反应,在298-900 K,反应速率随温度的升高明显增加;当温度进一步升高,其增加的趋势有所减缓。