Fe~＋（~6D）催化N_2O和CH_4制取甲醇循环反应的理论探究

本文采用密度泛函理论（DFT）B3LYP与耦合簇（CCSD）方法,研究了气相中四重态和六重态势能面上Fe＋催化N2O和CH4制取甲醇的微观机理.运用分子轨道理论和自然键轨道理论（NBO）对反应势能面进行分析,并通过自旋-轨道耦合（SOC）计算,讨论了势能面的交叉情况和自旋翻转的可能性.对Kozuch提出的能量跨度模型引入系间窜越几率加以修正,使其适用于非绝热两态反应.用修正后的能量跨度模型计算了催化剂的转化频率（TOF）,同时确定了整个反应的决速态.     

气相中W+活化CO2分解的自旋禁阻反应机理

用密度泛函理论中的UB3LYP方法, 对W采用相对论校正赝势基组(SDD), 对C、O采用6-311+G(3d)基组, 研究了气相中不同自旋态W+活化CO2分解的反应机理. 计算结果表明, W+活化CO2分解反应以六重态进入反应通道, 经过六重态势能面到四重态势能面的系间窜越(ISC),最后产物WO+和CO以四重态离开反应通道. 运用Harvey方法优化出最低能量交叉点(MECP), 并计算了MECP处的自旋-轨道耦合(SOC)常数(494.95 cm-1), 势能面的交叉和在MECP处较强的自旋-轨道耦合作用降低了自旋禁阻反应能垒, 为反应提供了一条低能反应路径, 反应总放热量为122.33 kJ·mol-1.     

气相中W~+活化CO_2分解的自旋禁阻反应机理

用密度泛函理论中的UB3LYP方法,对W采用相对论校正赝势基组(SDD),对C、O采用6-311+G(3d)基组,研究了气相中不同自旋态W+活化CO2分解的反应机理.计算结果表明,W+活化CO2分解反应以六重态进入反应通道,经过六重态势能面到四重态势能面的系间窜越(ISC),最后产物WO+和CO以四重态离开反应通道.运用Harvey方法优化出最低能量交叉点(MECP),并计算了MECP处的自旋-轨道耦合(SOC)常数(494.95cm-1),势能面的交叉和在MECP处较强的自旋-轨道耦合作用降低了自旋禁阻反应能垒,为反应提供了一条低能反应路径,反应总放热量为122.33kJ.mol-1.     

气相中Ta~+活化分解CO_2的自旋禁阻反应机理

采用密度泛函理论中的UB3LYP方法,对Ta采用相对论校正赝势基组(SDD),C,O采用6-311+G(3d)基组,计算研究了气相中不同自旋态下Ta+活化CO2分解反应的微观机理.计算结果表明,Ta+活化CO2分解反应是典型的自旋禁阻反应,以五重态进入反应通道,经过五重态势能面到三重态势能面的系间窜越(ISC),最终产物TaO+,CO以三重态离开反应通道.运用Harvey方法优化出最低能量交叉点(MECP),并计算了MECP处的自旋-轨道耦合(SOC)常数(204.94cm-1),较大的SOC值说明了在MECP处的自旋-轨道耦合作用较强,势能面的交叉降低了自旋禁阻反应能垒,为反应提供了一条低能反应路径.运用Landau-Zener跃迁几率公式计算了MECP处系间窜越几率,较大的系间窜越几率说明了该自旋禁阻反应速率较快.     

1CH2＋N2O反应的势能面

利用密度泛函理论（B3LYP）计算了1CH2＋N2O反应的反应物、中间体、过渡态及产物 的几何构型.进而用从头算方法（QCISD（T））计算了单点能量.由此描绘了反应的势能面, 确定了反应的最终产物通道为N2＋H2CO和NO＋HCN＋H.后者比前者有更大的分支比.N2、H2CO 、NO、HCN的存在有待于实验检测.作者认为,反应在室温下是加成-消除机理,而在高温下 可以通过直接取代的机理获得N2＋H2CO.     

CH2SH与NO2双自由基反应机理的理论研究

在B3LYP/6-311＋＋G(2df,p)水平上优化了标题反应驻点物种的几何构型, 并在相同水平上通过频率计算和内禀反应坐标(IRC)分析对过渡态结构及连接性进行了验证. 采用双水平计算方法HL//B3LYP/6-311＋＋G(2df,p)对所有驻点及部分选择点进行了单点能校正, 构建了CH2SH＋NO2反应体系的单重态反应势能剖面. 研究结果表明, CH2SH与NO2反应体系存在4条主要反应通道, 两个自由基中的C与N首先进行单重态耦合, 形成稳定的中间体HSCH2NO2 (a). 中间体a经过C—N键断裂和H(1)—O(2)形成过程生成主要产物P1 (CH2S＋trans-HONO), 此过程需克服124.1 kJ•mol－1的能垒. 中间体a也可以经过C—N键断裂及C—O键形成转化为中间体HSCH2ONO (b), 此过程的能垒高达238.34 kJ•mol－1. b再经过一系列的重排异构转化得到产物P2 (CH2S＋cis-HONO), P3 (CH2S＋HNO2)和P4 (SCH2OH＋NO). 所有通道均为放热反应, 反应能分别为－150.37, －148.53, －114.42和－131.56 kJ•mol－1. 标题反应主通道R→a→TSa/P1→P1的表观活化能为－91.82 kJ•mol－1, 此通道在200～3000 K温度区间内表观反应速率常数三参数表达式为kCVT/SCT＝8.3×10－40T4.4 exp(12789.3/T) cm3•molecule－1•s－1.     

金催化CO氧化湿度增强效应的理论研究

实验发现纳米金催化的CO氧化有良好的湿度增强效应,但有关机制仍不清楚.我们应用密度泛函理论研究了湿度增强效应的微观机制,以Au4团簇为例,研究了金催化CO氧化的微观机理,考察了H2O在反应中的角色和作用.计算结果表明,H2O与Au4团簇一样,在反应中扮演催化剂的角色,参与反应的进行、改变反应历程、降低反应能垒.催化循环包含4个基元步骤：O2＋H2O→OOH＋OH,CO＋OOH→CO2＋OH,CO＋OH→COOH,和COOH＋OH→CO2＋H2O,其中自由基OOH和OH的形成是催化循环的速控步骤,其能垒为100.31kJ/mol,明显低于非水参与反应的能垒（161.41kJ/mol）.目前的结果合理地解释了实验观测的CO催化氧化的湿度增强效应,给出了其微观反应机制.     

D-苯丙氨酸与Cu＋(1S0, 3d10)气相反应理论研究

用量子化学密度泛函(DFT)方法研究D-苯丙氨酸与一价基态金属阳离子Cu^＋在气相中反应的机理. 在B3LYP/6-31G*水平上, 优化了反应包含的4个反应通道的反应物、中间体、过渡态和产物的几何构型, 并采用B3LYP/DZVP, B3LYP/[6-311＋G**(C,H,O)＋Lanl2dz(Cu)], B3LYP/6-311＋G**, MP2/6-311＋G** 等方法对各驻点进行了单点能计算. 通过对计算结果的分析, 获得了其单重态反应势能面的一般轮廓、各驻点几何构型优化参数, 明确了其反应机理.     

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～3000 K温度范围内的速率常数kTST, kCVT和kCVT/SCT. 研究结果表明, 该反应体系共存在5个反应通道, 其中N进攻巯基上H原子生成CH3S＋HNO2的通道活化势垒较低, 为主要反应通道. 动力学数据也表明, 该通道在200～3000 K计算温度范围内占绝对优势, 拟合得到的速率常数表达式为k1CVT/SCT＝1.93×10－16T0.21exp(－558.2/T) cm3&;#8226;molecule－1&;#8226;s－1.     

气相中Y＋, Zr＋, Nb＋与CO2反应的理论研究

以Y^＋, Zr^＋, Nb^＋与CO₂反应作为第二前过渡金属离子与CO₂反应的范例体系. 采用密度泛函UB3LYP方法, 对于Y, Zr, Nb采用Stuttgart赝势基组, 对于CO₂采用6-311＋G(2d)基组, 计算研究了三种金属离子在基态和激发态时与CO₂气相反应的机理. 结果表明三种金属离子与CO₂反应以高自旋进入反应通道, 在反应过程中发生系间窜越, 以低自旋中间体和最终产物离开反应通道. 用内禀坐标单点垂直激发计算的方法找到了势能面交叉点, 并作了相应的讨论. 因为有金属离子的参与, 使单分子CO₂的强吸热分解反应变为生成CO和MO^＋的放热过程.     

On the catalytic role of Re+ in the oxygen transport activation of N2O by CO

The mechanism of the cyclic reaction N₂O(¹∑) + CO(¹∑⁺) → N₂(¹∑ _g ⁺ ) + CO₂(¹∑ _g ⁺ ) catalyzed by Re⁺ has been investigated on quintet and septet potential energy surfaces (PES). The reactions were studied by the B3LYP density functional method and the CCSD(T) theory. The calculated results of different PES show that the reaction proceeds in a two-step manner and spin crossing between different PES occurs. The involving crossing points (CPs) between the quintet and septet PES have been discussed by means of the intrinsic reaction coordinate approach. And the O-atom affinities testified that Re⁺ can capture O from N₂O and transfer O atom to CO in the two spin state, which are thermodynamically allowed. Furthermore, the spin–orbit coupling (SOC) is calculated between electronic states of different multiplicities at the CPs. For CP1 and CP2, the computed SOC constants are 8.34 and 10.09 cm^?1, respectively, obtained by using one-electron spin–orbit Hamiltonian in GAMESS. Therefore, the intersystem crossing at CP1 and CP2 occurs with a little probability because of the small SOC involved.     

Ti^＋离子和C2H4分子自旋禁阻反应中C—H键活化机理的理论研究

用密度泛函B3LYP方法,研究了二重态和四重态势能面自旋禁阻反应Ti^＋（^4F,3d^24s^1）＋C2H4→TiC2H2^＋（^2A2）＋H2的微观机理．通过自旋．轨道耦合的计算讨论了势能面交叉点和可能的自旋翻转过程．中间体IM1-^4B2处,四重态和二重态间的旋-轨耦合值为59．3cm^-1．自旋多重度必将发生变化,从四重态系间穿越到二重态势能面形成共价型复合物IM1-^2A1,同时导致四重态势能面的势垒明显降低．到插入中间体IM2后,二重态势能面上有两条不同的反应路径,即分步和协同路径,后者在二重态势能面上得到放热产物TiC2H2^＋（^2A2）＋H2具有较低的活化势垒,4．52kcal／mol,其主反应路径为：Ti^＋＋C2H4→^4IC→IM1—^4B2→4.2ISC→IM1—^2A1→[^2TSins]→IM2-^2A”→[^2TSMCTS]→IM5→TiC2H2^＋（^2A2）＋H2．     

氧负离子与乙烯自由基反应的理论研究

在G3MP2B3理论水平下研究了氧负离子与乙烯自由基的反应机理. 反应入口势能面的刚性扫描显示: 对于不同的初始反应取向, 体系存在3种不同的反应机理, 分别对应直接脱水、插入反应和直接键合成中间体通道. 其中, 通过插入反应形成的富能中间体[CH2＝C—OH]－及键合中间体[CH2＝CHO]－都可以进一步经异构化和解离生成其它各种可能产物, 如C2H－＋H2O, OH－＋CH2C和 ＋CO产物通道. 基于计算得到的反应势垒的相对高度, 直接脱水反应显然是该反应体系最主要的产物通道, 同时我们还结合Mulliken电荷布居分析研究了其中涉及的电子交换过程. 由此, 计算结果证实了以往OH－与C2H2反应的实验研究结果. 此外, 还对比了该反应体系、氧原子与乙烯自由基、氧负离子与乙烯分子三个反应的不同机理.     

Theoretical study of the reaction of Ni with OCS

The reactivity of Ni⁺ with OCS on both doublet and quartet potential energy surfaces (PES) has been investigated at the B3LYP/6-311+G(d) level. The object of this investigation was the elucidation of the reaction mechanism. The calculated results indicated that both the CS and CO bond activations proceed via an insertion–elimination mechanism. Intersystem crossing between the doublet and quartet surfaces may occur along both the CS and CO bond activation branches. The ground states of NiS⁺ and NiO⁺ were found to be quartets, whereas NiCO⁺ and NiCS⁺ have doublet ground states. The CS bond activation is energetically much more favorable than the CO bond activation. All theoretical results are in line with early experiments.     

A two-state computational investigation of methane C--H and ethane C--C oxidative addition to [CpM(PH3)]n+ (M = Co, Rh, Ir; n = 0, 1)

Reductive elimination of methane from methyl hydride half-sandwich phosphane complexes of the Group 9 metals has been investigated by DFT calculations on the model system [CpM(PH(3))(CH(3))(H)] (M = Co, Rh, Ir). For each metal, the unsaturated product has a triplet ground state; thus, spin crossover occurs during the reaction. All relevant stationary points on the two potential energy surfaces (PES) and the minimum energy crossing point (MECP) were optimized. Spin crossover occurs very near the sigma-CH(4) complex local minimum for the Co system, whereas the heavier Rh and Ir systems remain in the singlet state until the CH(4) molecule is almost completely expelled from the metal coordination sphere. No local sigma-CH(4) minimum was found for the Ir system. The energetic profiles agree with the nonexistence of the Co(III) methyl hydride complex and with the greater thermal stability of the Ir complex relative to the Rh complex. Reductive elimination of methane from the related oxidized complexes [CpM(PH(3))(CH(3))(H)](+) (M = Rh, Ir) proceeds entirely on the spin doublet PES, because the 15-electron [CpM(PH(3))](+) products have a doublet ground state. This process is thermodynamically favored by about 25 kcal mol(-1) relative to the corresponding neutral system. It is essentially barrierless for the Rh system and has a relatively small barrier (ca. 7.5 kcal mol(-1)) for the Ir system. In both cases, the reaction involves a sigma-CH(4) intermediate. Reductive elimination of ethane from [CpM(PH(3))(CH(3))(2)](+) (M = Rh, Ir) shows a similar thermodynamic profile, but is kinetically quite different from methane elimination from [CpM(PH(3))(CH(3))(H)](+): the reductive elimination barrier is much greater and does not involve a sigma-complex intermediate. The large difference in the calculated activation barriers (ca. 12.0 and ca. 30.5 kcal mol(-1) for the Rh and Ir systems, respectively) agrees with the experimental observation, for related systems, of oxidatively induced ethane elimination when M = Rh, whereas the related Ir systems prefer to decompose by alternative pathways.     

Theoretical investigation of the reactions of La atom and La<Superscript>+</Superscript> cation with carbonyl sulfide

The potential energy surfaces for the La+SCO and La⁺+ SCO reactions have been theoretically investigated by using the DFT (B3LYP/ECP/6-311+G(2d)) level of theory. Both ground and excited state potential energy surfaces (PES) are discussed. The present results show that the reaction mechanism is insertion mechanism both along the C-S and C-O bond activation branches, but the C-S bond activation is much more favorable in energy than the C-O bond activation. The reaction of La atom with SCO was shown to occur preferentially on the ground state (doublet) PES throughout the reaction process, and the experimentally observed species, have been explained according to the mechanisms revealed in this work. While for the reaction between La⁺ cation with SCO, it involves potential energy curve-crossing which dramatically affects reaction mechanism, and the crossing points (CPs) have been localized by the approach suggested by Yoshizawa et al. Due to the intersystem crossing existing in the reaction process of La⁺ with SCO, the products SLa⁺(η²CO) and OLa⁺(η²CS) may not form. This mechanism is different from that of La + SCO system. All our theoretical results not only support the existing conclusions inferred from early experiment, but also complement the pathway and mechanism for this reaction.     

气相VO（∑＋）与CH3OH（1^A′）分子自旋禁阻反应C—H，O—H键活化机理的理论研究

采用密度泛函理论（DFT）B3LYP与cCsD方法研究了二重态和四重态势能面自旋禁阻反应VO（∑’）活化cH30H（1^A′）分子c—H,0—H键的微观机理．通过自旋一轨道耦合的计算讨论了势能面交叉点和可能的自旋翻转过程．在MEcP处,四重态和二重态问的旋轨耦合常数为131．14cm^-1．自旋多重度发生改变,从四重态系间穿越到二重态势能面形成中间体2^IM1,导致反应势能面的势垒明显降低．     

Spin-inversion mechanisms in the reactions of transition metal cations (Sc+, Ti+, V+, Cr+, Mn+, Fe+, Co+, Ni+, and Cu+) with OCS in the gas phase: A perspective from automated reaction path search calculations

The mechanisms including spin-inversion have been systematically studied for the M⁺ + OCS → MS⁺ + CO/MO⁺ + CS (M denotes a transition metal from Sc to Cu) ion-molecule reactions using the automated reaction path search method. We used the lowest mixed-spin potential energy surface obtained from the diagonalization of the spin-coupled Hamiltonian matrix, whose diagonal elements are taken to be the lowest two spin states. This scheme can effectively locate approximate minimum energy crossing points between the two potential energy surfaces with different spin multiplicities. The spin-orbit couplings at spin-inversion points have been calculated to understand the efficiencies of nonadiabatic transitions. The obtained reaction pathways and the calculated spin-orbit couplings are employed to interpret previous experimental studies.