CH3O与ClO双自由基反应机理的量子化学研究

在QCISD(T)/6-311＋G(d,p)//B3LYP/6-311＋G(3df,3pd)水平上, 对CH₃O与ClO双自由基反应进行了理论研究. 结果表明, 该反应共有三个反应通道, 产物分别为HOCl＋CH₂O, CH₂O₂＋HCl和CH₃Cl＋O₂(¹Δ). 不论从动力学角度, 还是从热力学角度看, 形成产物HOCl＋CH₂O的通道均是最有利的, 因此为主要反应通道, 这与实验观察到的结果是一致的.     

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** 等方法对各驻点进行了单点能计算. 通过对计算结果的分析, 获得了其单重态反应势能面的一般轮廓、各驻点几何构型优化参数, 明确了其反应机理.     

激发态Ti＋(3d14s2)与丙炔醇气相反应理论研究

用量子化学密度泛函(DFT)方法研究了激发态Ti^＋(3d¹4s²)与丙炔醇(PPA)气相反应的机理. 在B3LYP/DZVP水平上, 优化了反应的两个通道的反应物、中间体、过渡态和产物的几何构型, 并在MP4/[6-311＋G**(C,H,O)＋Lanl2dz (Ti)]水平上计算了各驻点的单点能量. 为了确证过渡态的真实性, 在B3LYP/DZVP水平上进行了内禀坐标(IRC)计算和频率分析, 获得了二重态反应势能面, 确定了反应机理. 研究结果表明生成产物为[C₃H₃O]^＋和Ti—H的通道是主要反应途径.     

激发态Ti＋(3d14s2)与丙炔醇气相反应理论研究

用量子化学密度泛函(DFT)方法研究了激发态Ti^＋(3d¹4s²)与丙炔醇(PPA)气相反应的机理. 在B3LYP/DZVP水平上, 优化了反应的两个通道的反应物、中间体、过渡态和产物的几何构型, 并在MP4/[6-311＋G**(C,H,O)＋Lanl2dz (Ti)]水平上计算了各驻点的单点能量. 为了确证过渡态的真实性, 在B3LYP/DZVP水平上进行了内禀坐标(IRC)计算和频率分析, 获得了二重态反应势能面, 确定了反应机理. 研究结果表明生成产物为[C₃H₃O]^＋和Ti—H的通道是主要反应途径.     

CH3S与HCS双自由基反应的密度泛函理论研究

应用量子化学从头算和密度泛函理论(DFT)对CH₃S与HCS双自由基单重态反应进行了研究. 在MPW1PW91/ 6-311G(d,p)水平上优化了反应通道上各驻点(反应物、中间体、过渡态和产物)的几何构型, 用内禀反应坐标(IRC)计算和频率分析方法对过渡态进行了验证. 在QCISD(t)/6-311＋＋G(d,p)水平上计算各物种的单点能, 并对总能量进行了零点能校正. 研究结果表明, CH₃S与HCS反应为多通道反应, 有4条可能的反应通道, 反应物首先通过S…S弱相互作用形成具有竞争反应机理的五元环硫-硫偶合中间体a和链状硫-硫偶合中间体c, 再由此经过氢迁移、离解、异构化等不同机理得到主要产物P1 (2CH₂S), 次要产物P2 (CH₃SH＋CS), P3 (CH₄＋CS₂)和P4 [CH₂(SH)CSH]. 根据势能面分析, 所有反应均为放热反应, 生成P1的反应热为－165.55 kJ•mol^－1. 通道R→a→TSa/b→b→P1为标题反应的主通道, 其速控步骤a→TSa/b→b在200～2000 K温度区间内的速率常数可以表示为k₁^CVT/SCT＝1.75×10¹⁰T^0.65exp(－907.6/T) s^－1. P3及P4的生成需要越过很高的活化能垒, 是动力学禁阻步骤, 但在反应体系中加入合适催化剂, 改变其反应机理, 有可能使生成CH₂(SH)CSH, CH₄及CS₂的反应易于进行.     

CH3S自由基H迁移异构化及脱H2反应的直接动力学研究

采用密度泛函方法(MPW1PW91)在6-311G(d,p)基组水平上研究了CH₃S自由基H迁移反应CH₃S→CH₂SH (R1), 脱H₂反应CH₃S→HCS＋H₂ (R2)以及脱H₂产物HCS异构化反应HCS→CSH (R3)的微观动力学机理. 在QCISD(t)/6- 311＋＋G(d,p)//MPW1PW91/6-311G(d,p)＋ZPE水平上进行了单点能校正. 利用经典过渡态理论(TST)与变分过渡态理论(CVT)分别计算了各反应在200～2000 K温度区间内的速率常数k^TST和k^CVT, 同时获得了经小曲率隧道效应模型(SCT)校正后的速率常数k^CVT/SCT. 结果表明, 反应 R1, R2 和R3的势垒△E^≠分别为160.69, 266.61和241.63 kJ/mol, R1为反应的主通道. 低温下CH₃S比CH₂SH稳定, 高温时CH₂SH比CH₃S更稳定. 另外, 速率常数计算结果显示, 量子力学隧道效应在低温段对速率常数的计算有显著影响, 而变分效应在计算温度段内对速率常数的影响可以忽略.     

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

采用密度泛函B3LYP/6-311G**和高级电子相关耦合簇CCSD(T)/6-311G**方法计算研究了CH₃与NO反应机理, 全参数优化了反应势能面上各驻点的几何构型, 用内禀反应坐标(IRC)计算和频率分析方法, 对过渡态进行了验证. 研究结果表明: CH₃与NO是一多通道多步骤的复杂反应, 可以分别在单重态和三重态势能面上进行. 经过缔合, 氢转移和离解等复杂过程, 最终得到8种产物(P1～P8).     

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

采用密度泛函B3LYP/6-311G**和高级电子相关耦合簇CCSD(T)/6-311G**方法计算研究了CH₃与NO反应机理, 全参数优化了反应势能面上各驻点的几何构型, 用内禀反应坐标(IRC)计算和频率分析方法, 对过渡态进行了验证. 研究结果表明: CH₃与NO是一多通道多步骤的复杂反应, 可以分别在单重态和三重态势能面上进行. 经过缔合, 氢转移和离解等复杂过程, 最终得到8种产物(P1～P8).     

CH2=CHCl与O(3P)反应的理论研究

用量子化学密度泛函理论和QCISD(Quadratic configuration interaction calculation)方法,对0(^3P)与CH2CHCl的反应进行了理论研究．在UB3LYP／6—311 G(d,p),UB3LYP／6—31 (3df,3pd)计算水平上,优化了反应物、产物、中间体和过渡态的几何构型,并在UQCISD(T)／6—311 G(2df,2pO)水平上计算了单点能量．为了确证过渡态的真实性,在UB3LYP／6—311 G(3df,3pd)水平上进行了内禀坐标(IRC)计算和频率分析,并确定了反应机理．研究结果表明,反应主要产物为CH2CHO和Cl．     

CH3S与HCS双自由基反应的密度泛函理论研究

应用量子化学从头算和密度泛函理论(DFT)对CH₃S与HCS双自由基单重态反应进行了研究. 在MPW1PW91/ 6-311G(d,p)水平上优化了反应通道上各驻点(反应物、中间体、过渡态和产物)的几何构型, 用内禀反应坐标(IRC)计算和频率分析方法对过渡态进行了验证. 在QCISD(t)/6-311＋＋G(d,p)水平上计算各物种的单点能, 并对总能量进行了零点能校正. 研究结果表明, CH₃S与HCS反应为多通道反应, 有4条可能的反应通道, 反应物首先通过S…S弱相互作用形成具有竞争反应机理的五元环硫-硫偶合中间体a和链状硫-硫偶合中间体c, 再由此经过氢迁移、离解、异构化等不同机理得到主要产物P1 (2CH₂S), 次要产物P2 (CH₃SH＋CS), P3 (CH₄＋CS₂)和P4 [CH₂(SH)CSH]. 根据势能面分析, 所有反应均为放热反应, 生成P1的反应热为－165.55 kJ•mol^－1. 通道R→a→TSa/b→b→P1为标题反应的主通道, 其速控步骤a→TSa/b→b在200～2000 K温度区间内的速率常数可以表示为k₁^CVT/SCT＝1.75×10¹⁰T^0.65exp(－907.6/T) s^－1. P3及P4的生成需要越过很高的活化能垒, 是动力学禁阻步骤, 但在反应体系中加入合适催化剂, 改变其反应机理, 有可能使生成CH₂(SH)CSH, CH₄及CS₂的反应易于进行.     

CH2CO＋CH2的多通道反应的理论研究

利用密度泛函(DFT)和自然键轨道理论(NBO)及高级电子耦合簇[CCSD(T)]和电子密度拓扑(AIM)方法, 对单重态和三重态CH₂与CH₂CO反应的微观机理进行了研究. 在B3LYP/6-311＋G(d,p)水平上优化了反应通道各驻点的几何构型. 在CCSD(T)/6-311＋G(d,p)水平上计算了各物种的单点能量, 并对总能量进行了校正. 计算表明, 单重态CH₂与CH₂CO的C—H键可发生插入反应, 与C＝C、C＝O可发生加成反应, 存在三条反应通道, 产物为CO和C₂H₄, 从能量变化和反应速控步骤能垒两方面考虑, 反应II更容易发生. 对反应通道中的关键点进行了自然键轨道及电子密度拓扑分析. 三重态CH₂与CH₂CO的反应存在三条反应通道, 一条是与C-H键的插入反应, 另一条是三重态CH₂与C＝C发生加成反应, 产物为CO和三重态C₂H₄, 通道II势垒较低, 更容易发生. 最后一条涉及双自由基的反应活化能最大, 最难发生.     

HO2自由基与NO2反应通道的理论研究

应用密度泛函理论DFT/B3LYP对HO₂＋NO₂反应进行了研究, 在B3LYP/6-311G^**和CCSD(T)/6-311G^**水平上计算了HO₂自由基与NO₂分子反应的单重态和三重态反应势能面, 计算结果表明, 单重态反应势能面中的直接氢抽提反应机理是此反应的主要反应通道, 即HO₂自由基的氢原子转移到NO₂分子的氮原子上形成产物P₁ (HNO₂＋³O₂), 另一个可能的反应通道是单重态反应势能面上HO₂中的端位氧原子进攻NO₂分子中的氮原子形成中间体1 (HOONO₂), 接着中间体1 (HOONO₂)经过氢转移形成产物P₂ (trans-HONO＋³O₂), 以上两个反应通道都是放热反应通道, 分别放热90.14和132.52 kJ•mol^－1.     

HO2自由基与NO2反应通道的理论研究

应用密度泛函理论DFT/B3LYP对HO₂＋NO₂反应进行了研究, 在B3LYP/6-311G^**和CCSD(T)/6-311G^**水平上计算了HO₂自由基与NO₂分子反应的单重态和三重态反应势能面, 计算结果表明, 单重态反应势能面中的直接氢抽提反应机理是此反应的主要反应通道, 即HO₂自由基的氢原子转移到NO₂分子的氮原子上形成产物P₁ (HNO₂＋³O₂), 另一个可能的反应通道是单重态反应势能面上HO₂中的端位氧原子进攻NO₂分子中的氮原子形成中间体1 (HOONO₂), 接着中间体1 (HOONO₂)经过氢转移形成产物P₂ (trans-HONO＋³O₂), 以上两个反应通道都是放热反应通道, 分别放热90.14和132.52 kJ•mol^－1.     

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

采用密度泛函理论B3LYP方法研究了SiH₂自由基与HNCO的反应机理, 并在B3LYP/6-311＋＋G**水平上对反应物、中间体、过渡态进行了全几何参数优化, 通过频率分析和内禀反应坐标(IRC)确定了中间体和过渡态. 为了得到更精确的能量值, 又用QCISD(T)/6-311＋＋G**方法计算了在B3LYP/6-311＋＋G**水平优化后的各个驻点的相对能量. 计算结果表明SiH₂自由基与HNCO的反应有五条反应通道, 其中顺式反应通道SiH₂＋HNCO→IM3→ TS4→IM5→TS5→IM6→SiH₂NH＋CO反应能垒最低, 为主反应通道.     

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

采用密度泛函理论B3LYP方法研究了SiH₂自由基与HNCO的反应机理, 并在B3LYP/6-311＋＋G**水平上对反应物、中间体、过渡态进行了全几何参数优化, 通过频率分析和内禀反应坐标(IRC)确定了中间体和过渡态. 为了得到更精确的能量值, 又用QCISD(T)/6-311＋＋G**方法计算了在B3LYP/6-311＋＋G**水平优化后的各个驻点的相对能量. 计算结果表明SiH₂自由基与HNCO的反应有五条反应通道, 其中顺式反应通道SiH₂＋HNCO→IM3→ TS4→IM5→TS5→IM6→SiH₂NH＋CO反应能垒最低, 为主反应通道.     

Criegee自由基CH2O2与H2O反应机理及动力学的理论研究

在6-311＋G(2d,2p)水平下, 采用密度泛函理论(DFT)的B3LYP方法, 研究了Criegee 自由基CH₂O₂与H₂O的反应. 结果表明反应存在三个通道: CH₂O₂＋H₂O®HOCH₂OOH (R1); CH₂O₂＋H₂O®HCO＋OH＋H₂O (R2); CH₂O₂＋H₂O®HCHO＋H₂O₂ (R3), 各通道的势垒高度分别为43.35, 85.30和125.85 kJ/mol. 298 K下主反应通道(R1)的经典过渡态理论(TST)与变分过渡态理论(CVT)的速率常数k^TST与k^CVT均为2.47×10^－17 cm³•molecule^－1•s^－1, 而经小曲率隧道效应模型(SCT)校正后的速率常数k^CVT/SCT为 5.22×10^－17 cm³•molecule^－1•s^－1. 另外, 还给出了200～2000 K 温度范围内拟合得到的速率常数随温度变化的三参数Arrhenius方程.     

三接触弯曲构型的蓝移氢键CH3…Y的理论研究

运用量子化学从头算方法, 在MP2/6-311＋＋G(d,p), MP2/6-311＋＋G(2df,2p), MP2/6-311＋＋G(3df,3pd)和QCISD/6-311＋＋G(d,p)水平上, 研究了CH₃F, CH₃Cl和CH₃Br作为质子给体与Cl^－, Br^－作为质子接受体形成的氢键CH₃…Y. 计算结果表明: 6种复合物中C—H键收缩, 伸缩振动频率增大, 形成蓝移氢键. 分子中原子(Atoms in Molecules, AIM)分析表明, 这些复合物的电子密度拓扑性质与普通氢键有着本质的不同, 在Y…H之间不存在键临界点, 而在Y与C之间存在键临界点, 因此这些相互作用严格地不能称为氢键. 自然键轨道(Natural bond orbital, NBO)分析表明, 在这些复合物中弯曲的CH…Y的特殊结构使得分子间超共轭n(Y)®σ*(C—H)减小到可以忽略; 质子接受体的电子密度没有转移到σ*(C—H)上, 而是转移到了σ*(C—X) (X＝F, Cl, Br)上; 存在一定程度的重杂化; 分子内超共轭相互作用减小使得σ*(C—H)的电子密度减少. 这些因素共同导致C—H伸缩振动频率的蓝移.     

CCl2与CH2O插入反应机理及热力学与动力学特性的理论研究

采用密度泛函B3LYP/6-311G*和高级电子相关耦合簇[CCSD(T)/6-311G*]方法计算研究了CCl₂与CH₂O的插入反应机理, 全参数优化了反应势能面各驻点的几何构型, 用内禀反应坐标(IRC)和频率分析方法, 对过渡态进行了验证. 研究结果表明: 反应(1)是单重态二氯卡宾与甲醛插入反应的主反应通道. 该反应由两步组成: (i)两反应物首先经一无能垒的放热反应, 放出9.73 kJ•mol^－1的热量, 生成一中间体IM1, (ii)中间体IM1经一过渡态TS1, 发生H的转移, 生成产物P1, 其势垒为47.32 kJ•mol^－1. 用RRKM-TST理论计算了300～1900 K温度范围内反应(1)的压力效应. 用经Wigner校正的Eyring过渡态理论研究了不同温度下该反应的热力学和动力学性质. 从热力学和动力学角度综合分析, 在高压限101325 Pa下, 该反应进行的适宜温度范围为400～1800 K, 如此, 反应既有较大的自发趋势和平衡常数, 又具有较快的反应速率.