气相中HRnCCH和X(X=H_2O,NH_3)反应机理的理论研究

在CCSD(T)//MP2/aug-cc-pVTZ-pp理论水平上,研究了HRnCCH与大气中H2O及NH3分子反应的机理,反应主要包括HRnCCH与HRnOH及HRnNH2之间的转化、H2O和NH3在HRnCCH中的碳碳三键上的加成反应以及HRnCCH与双分子水反应等.结果表明,HRnCCH与H2O反应生成HCCH和HRnOH及HRnCCH与NH3反应生成HCCH和HRnNH2的能垒分别为54.1和75.2 kJ/mol,而生成HRnCHC(OH)H,HRnC(OH)CH2,HRnCHC(NH2)H和HRnC(NH2)CH2的活化能分别为219.6,220.5,174.4和182.4kJ/mol,此结果表明HRnCCH反应性较弱且是稳态存在的.此外,在HRnCCH与H2O反应中加入单个水分子,仍然生成HRnCHC(OH)H,但反应活化能却降低了96.4 kJ/mol,说明水分子对该反应有明显的催化作用.     

HNCO+OH-→NH2+CO2反应理论研究

用从头算UHF/6-31G基组研究了异氰酸和羟基生成氨基和二氧化碳即HNCO+OH→NH2+CO2的反应机理.优化得到了反应途径上的过渡态和中间体,并通过振动分析对过渡态和中间体进行了确认.在UMP4/6-31G水平上计算了它们的能量,同时对零点能进行了较正.计算结果表明:此反应是多步反应,先后通过3个过渡态(TS1,TS2,TS3),2个内旋转位垒(PSI,TSII),4个中间体(IM1,IM2,IM3,IM4),其中,IM3→TS2这一步为整个反应的决速步骤,速控步的活化能为202.388lJ/mol.与异氰酸和羟基作用的另一反应通道(即HNCO+OH→H2O+NCO)的活化能(69.038kJ/mol)比较,可看出所研究反应通道为次要反应通道,这与实验结果是一致的.     

金催化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催化氧化的湿度增强效应,给出了其微观反应机制.     

HNCO + CN反应途径的从头计算

在6-311G(d, p)基组水平上, 采用全电子的UMP2方法对HNCO和CN自由基的反应途径进行了研究, 结果表明, 反应存在如下两条反应通道: HNCO + CN→NCO + HCN (1); HNCO + CN→HNCN + CO (2). 其中反应(1)是吸氢反应, 氢转移过程在分子间以协同方式进行, UMP2(full)/6-311G(d, p)计算位垒是 20.80 kJ/mol, 在反应的温度区间(1000~2100 K), 传统过渡态理论得出的速率常数值范围是0.32×10^−11~6.9×10^−11 cm³·mol^−1·s^−1, 支持了实验预测结果. 反应(2)是分步反应, 在反应途径上存在分子中间体, UMP2(full)/6-311G(d, p)计算位垒为83.42 kJ/mol. 反应通道(2)的位垒大于反应通道(1).     

HCNO分子裂解机理的理论研究

实验提示应用248 nm UV波长对HCNO分子进行直接光解, 该分子可能发生裂解, 得到某些产物. 为了揭示HCNO分子的裂解机理, 选择HCNO分子的一组相对能级作为理论研究的起始点, 即¹A' (0.00 kJ/mol), ³A' (255.01 kJ/mol), ³A" (282.37 kJ/mol)和¹A" (341.59 kJ/mol), 进而找到了合理的反应路径, 阐明了相应的裂解机理, 得到的主要产物为H＋NCO, HCN＋O和NH＋CO, 与实验提示的结果相符合.     

CH3＋HNCO反应机理的理论研究

在6-311＋＋G**基组水平上,采用UMP2方法对自由基CH3与HNCO反应机理进行了研究，全参数优化了反应通道上各驻点的几何构型.结果表明, 自由基CH3与HNCO分子间反应有三条反应通道,第一为CH3与HNCO分子间经过生成一个稳定化能为4.56 kJ•mol－1的含氢键的分子复合物M后，经过渡态TS生成另一个产物复合物M′，然后分解为甲烷和NCO自由基；第二是CH3与HNCO分子间通过生成稳定反式中间体trans-int，其经过渡态trans-ts分解成产物CH3NH和CO；第三是CH3与HNCO分子间通过生成稳定顺式中间体cis-int，其经过渡态cis-ts分解成产物CH3NH和CO.比较三条反应通道的反应活化能，表明CH3与HNCO反应较易生成CH4＋NCO.     

Au4团簇上甲酸分解反应机理的理论研究

采用密度泛函理论方法对Au4团簇上甲酸分解反应的反应机理进行了研究,并考察了Au4团簇的两个催化活性位点。在路径Ⅰ和路径Ⅱ中,HCOOH分解的产物是CO2和H2。在路径Ⅲ和路径Ⅳ中,HCOOH分解的最终产物为CO和H2O。此外,本文研究了CO2、H2和CO、H2O两种产物的相互转化,即路径Ⅴ和路径Ⅵ。研究结果表明,路径Ⅰ和路径Ⅱ的活化自由能垒较低,即在Au4团簇上HCOOH更易分解得到CO2和H2,此外两种产物之间不容易转化。进一步研究发现团簇的大小及CeO2载体对HCOOH分解脱氢路径的活化自由能垒有一定的影响。     

NH2 + HNCO反应机理的从头计算

在6-311G(d,p)基组水平上, 采用全电子的UMP2和UQCISD(T)方法对自由基NH₂和HNCO反应机理进行了研究, 结果表明, 反应存在如下两条反应通道: NH₂ + HNCO→NH₃ + NCO (1)和NH₂ + HNCO→N₂H₃ + CO (2). 反应(1)是吸氢反应, QCISD(T,full)// MP2(full)/6-311G(d,p) 计算位垒是29.04 kJ/mol. 与实验估计值29.09 kJ/mol一致. 在反应的温度区间(2300～2700 K),传统过渡态理论得出的速率常数值的范～围是1.68×1011～3.29×10¹¹cm³·mol^-1·s^-1, 支持了反应速率常数应小于等于5.0×10¹¹cm³·mol^-1·s^-1的实验结论. 对反应(1), 理论研究还得出反应物分子可通过分子间作用生成氢键复合物(HBC), 其能量相对于反应物降低32.41 kJ/mol. 反应(2)是一个可经过顺式或反式方式进行的分步反应, 在反应分子间第1步生成N—N键, 再经过一个C—N键断裂过渡态生成产物. 反应(2)控速步骤的位垒为92.79 kJ/mol(顺式)或147.43 kJ/mol(反式). 反应(2)位垒高于反应(1).     

HNCO+OH——NH~2+CO~2反应理论研究

用从头算UHF/6-31G基组研究了异氰酸和羟基生成氨基和二氧化碳即HNCO+OH--NH~2+CO~2的反应机理.优化得到了反应途径上的过渡态和中间体,并通过振动分析对过渡态和中间体进行了确认.在UMP4/6-31G水平上计算了它们的能量,同时对零点能进行了较正.计算结果表明:此反应是多步反应,先后通过3个过渡态(TS1,TS2,TS3),2个内旋转位垒(TSI,TSII),4个中间体(IM1,IM2,IM3,IM4),其中,IM3--TS2这一步为整个反应的决速步骤,速控步的活化能为202.388kJ/mol.与异氰酸和羟基作用的另一反应通道(即HNCO+OH--H~2O+NCO)的活化能(69.038kJ/mol)比较,可看出所研究反应通道为次要反应通道,这与实验结果是一致的。     

HNCO+OH——NH~2+CO~2反应理论研究

用从头算UHF/6-31G基组研究了异氰酸和羟基生成氨基和二氧化碳即HNCO+OH--NH~2+CO~2的反应机理.优化得到了反应途径上的过渡态和中间体,并通过振动分析对过渡态和中间体进行了确认.在UMP4/6-31G水平上计算了它们的能量,同时对零点能进行了较正.计算结果表明:此反应是多步反应,先后通过3个过渡态(TS1,TS2,TS3),2个内旋转位垒(TSI,TSII),4个中间体(IM1,IM2,IM3,IM4),其中,IM3--TS2这一步为整个反应的决速步骤,速控步的活化能为202.388kJ/mol.与异氰酸和羟基作用的另一反应通道(即HNCO+OH--H~2O+NCO)的活化能(69.038kJ/mol)比较,可看出所研究反应通道为次要反应通道,这与实验结果是一致的。     

Theoretical study of formamide decomposition pathways

The chemical transformations of formamide (NH(2)CHO), a molecule of prebiotic interest as a precursor for biomolecules, are investigated using methods of electronic structure computations and Rice-Rampserger-Kassel-Marcus (RRKM) theory. Specifically, quantum chemical calculations applying the coupled-cluster theory CCSD(T), whose energies are extrapolated to the complete basis set limit (CBS), are carried out to construct the [CH(3)NO] potential energy surface. RRKM theory is then used to systematically examine decomposition channels leading to the formation of small molecules including CO, NH(3), H(2)O, HCN, HNC, H(2), HNCO, and HOCN. The energy barriers for the decarboxylation, dehydrogenation, and dehydration processes are found to be in the range of 73-78 kcal/mol. H(2) loss is predicted to be a one-step process although a two-step process is competitive. CO elimination is found to prefer a two-step pathway involving the carbene isomer NH(2)CHO (aminohydroxymethylene) as an intermediate. This CO-elimination channel is also favored over the one-step H(2) loss, in agreement with experiment. The H(2)O loss is a multistep process passing through a formimic acid conformer, which subsequently undergoes a rate-limiting dehydration. The dehydration appears to be particularly favored in the low-temperature regime. The new feature identifies aminohydroxymethylene as a transient but crucial intermediate in the decarboxylation of formamide.     

Ab initio study of the reaction path of the reaction HNCO+CN

HNCO is a convenient photolytic source of NCO and NH radicals for laboratory kinetics studies of elementary reaction[1] and plays an important role in the combustion and atmosphere chemistry. It can re- move deleterious compounds rapidly from exhausted ga…     

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     

Decomposition of protonated noradrenaline and normetanephrine assisted by NH2 migration studied by electrospray tandem mass spectrometry and molecular orbital calculations

As part of a research program on neurotransmitters in a biological fluid, the fragmentations characterising catecholamines protonated under electrospray ionisation (ESI) conditions, under low collision energy in a triple-quadrupole mass spectrometer, were investigated. The decompositions of protonated noradrenaline (VH) and normetanephrine (VIH) were studied. Both precursor ions eliminate first H2O at very low collision energy, and the fragmentations of [MH-H2O]+ occur at higher collision energy. The breakdown graphs of [MH-H2O]+ ions, with collision energy varying from 0-40 eV in the laboratory frame, are presented. [VIH-H2O]+ ions lose competitively NH3 and CH3OH. For [VH-H2O]+ the loss of NH3 is dominant while H2O is eliminated at very low abundance at all collision energies. All of these secondary fragmentations are followed at higher collision energies by elimination of CO. These fragmentations are interpreted by means of ab initio calculations up to the B3LYP/6-311+G(2d,2p) level of theory. The elimination of H2O requires first the isomerisation of N-protonated forms, chosen as energy references, to O-protonated forms. The isomerisation barriers are calculated to be lower than 81 kJ/mol above the N-protonated forms. The elimination of NH3 from [MH-H2O]+ requires first the migration, via a cyclisation, of the amine function from the linear chain to the aromatic ring in order to prevent the formation of unstable disubstituted carbocations in the ring. The barriers associated with the loss of NH3 are located 220 and 233 kJ/mol above VH and 219 kJ/mol above VIH. The energy barrier for the loss of ROH is located 236 and 228 kJ/mol above VH and VIH, respectively. The absence of ions corresponding to [VH-2H2O]+ is due to a parasitic mechanism with an activation barrier lower than 236 kJ/mol that leads to a stable species unable to fragment, thus preventing the second loss of H2O. Losses of CO following the secondary fragmentations involve activation barriers higher than 330 kJ/mol.     

Theoretical study of mechanism of forming a silapolycyclic compound between methylenesilylene and acetone

The mechanism of the cycloaddition reaction of forming a silapolycyclic compound between singlet methylenesilylene and acetone has been investigated with MP2/6‐31G* method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. The energies of the different conformations are calculated by CCSD(T)//MP2/6‐31G* method. From the potential energy profile, we predict that the cycloaddition reaction of forming a silapolycyclic compound between singlet methylenesilylene and acetone has two competitive dominant reaction pathways. First dominant reaction pathway consists of four steps: (I) the two reactants (R1, R2) first form an intermediate (INT1) through a barrier‐free exothermic reaction of 46.2 kJ/mol; (II) intermediate (INT1) then isomerizes to a planar four‐membered ring product (P3) via transition state (TS3) with an energy barrier of 47.1 kJ/mol; (III) planar four‐membered ring product (P3) further reacts with acetone (R2) to form an intermediate (INT4), which is also a barrier‐free exothermic reaction of 40.0 kJ/mol; (IV) intermediate (INT4) isomerizes to a silapolycyclic compound (P4) via transition state (TS4) with an energy barrier of 57.0 kJ/mol. Second dominant reaction pathway consists of three steps: (I) the two reactants (R1, R2) first form a four‐membered ring intermediate (INT2) through a barrier‐free exothermic reaction of 0.5 kJ/mol; (II) INT2 further reacts with acetone (R2) to form an intermediate (INT5), which is also a barrier‐free exothermic reaction of 45.4 kJ/mol; (III) intermediate (INT5) isomerizes to a silapolycyclic compound (P5) via transition state (TS5) with an energy barrier of 49.3 kJ/mol. P4 and P5 are isomeric compounds. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Isomerization of azomethine- and azo-compounds and their protonated species

Barriers of rotation and inversion, respectively, have been calculated for the species H₂CNH (I), H₂CNCH₃ (II), NHNH (III), NHNCH₃ (IV) and their protonated species. For any unprotonated molecule the barrier of inversion is consistently lower than the barrier of rotation. Tile inversion barriers are: 27.8 (I), 23.8 (II), 51.9 (III) and 46.1 (IV) kcal/mole. In the case of azomethine species, protonation results in an increased rotational barrier (from 50.8 to 74.7 kcal/mole for II). In the case of azo species barriers of inversion are lowered on protonation (from 51.9 to 30.1 for II and from 46.1 to 24.4 kcal/mole for IV). All barriers are given with reference to the trans-isomer (azo). Proton affinities for the azomethine species are higher than those of the corresponding azo species (223.3 for 1, 199.9 kcal/mole for II).     

<Emphasis Type="Italic">Ab initio</Emphasis> study of the reaction path of the reaction HNCO+CN

We report ab initio UMP2 calculations of the reaction of CN with HNCO using 6-311G(d,p) basis sets. The obtained results show that the reaction has two product channels: HNCO+CN→HCN+NCO (1) and HNCO+CN→HNCN+CO (2). Channel (1) is a hydrogen abstraction reaction, which is a concerted process. The calculated potential energy barrier is 20.80 kJ/mol at UMP2(full)/6-311G(d,p) level. In the range of reaction temperature (1000-2100 K), the conventional transition theory rate constant for channel (1) ranges from 0.32×10⁻¹¹ to 6.9×10⁻¹¹cm³· mol⁻¹· s⁻¹, which is close to the experimental value. Channel (2) is a stepwise reaction involving an intermediate during the process of reaction. The UMP2(full)/6-311G(d,p) potential energy barrier is 83.42 kJ/mol for the rate-controlling step, which is much higher than that of channel (1).