甲基苯热裂解机理的理论研究

第30卷第1期燃料化学学报ol. 30 No I2002年2月JOURNAL OF FUEL CHEMISTRY AND TECHNOLOGY文章编号:0253-24092002m1-0049-05甲基苯热裂解机理的理论研究乔占平2,王惠,赵文立,冉新权(1.西北大学化学系,陕西西安710069;2.南阳师院化学系,河南南阳473069)摘要:在半经验UAM方法及从头算UHF/3-21G*水平的基础上用考虑了电子相关效应的密度泛函UB3LYP321G※方法对碳前驱体甲基苯的热裂解反应机理进行了研究并对反应物、产物自由基及过渡态的结构进行了能量梯度法全优化。计算结果表明甲基苯的热裂解温度较低时,首先支持苯环甲基上的C-H键断裂继而发生生成联二甲苯的反应。随着温度的提高生成苯自由基的反应比例将增加同时发生少量的苯环上C-H键断裂的反应。该反应机理与实验结果相一致用密度泛函方法UB3LYP321G兴计算的反应标准焓变ΔH与实验值吻合较好。关键词:碳前驱体;甲基苯;热裂解机理;密度泛函UB3LYP3-21G兴中图分类号:0611,0641文献标识码:A低分子有机化合物的热裂解炭化机理及热解炭化初期的反应过程,目前尚未阐明清楚。它的研究T H对于有效利用有机质能源、开发各种新型碳材料如石墨电极、碳纤维、原子反应堆用石墨材料、医学用碳材料、及碳/碳复合材料等具有重要的意义。CH CH,由于有机化合物热裂解反应的复杂性利用实验手段进行上述两方面的研究显得非常有限对反但随着热解温度的提高(800℃左右)甲苯脱甲基的应历程只能根据不同温度段的取样分析进行大致的比例将要多些同时也有少量苯环上的脱氩反应发生。为了证实这一实验结果并荻得理论上的支持判断。量子化学理论计算在研究反应机理和预测反应方向方面有许多成功的范例1-61且量子化学用我们对甲基苯可能裂解的反应路径进行了如图1的于低级芳香烃作为碳碳复合材料碳前驱体热解机模拟设计。为了方便我们对不同反应路径、反应物理方面的研究国内外尚不多见。因此本文将以碳及自由基进行了编号见图1下面的叙述均以编前驱体甲苯为例首先用半经验UAM方法及从头号代之。算UHF/3-21G兴方法对其可能裂解的初期反应路径进行模拟设计及热力学计算。并在此基础上用考虑CH. +H了电子相关效应的密度泛函B3LYP7方法( Becksthree parameters hybrid method using the LYP correlationfunctional在3-21G水平上进行了非限制性计算同cHc④时将计算结果与热力学实验结果相比较以摸索一套适合低级芳香烃热裂解机理研究的有效计算方法。反应路径设计S Deutsch等8在研究甲苯的热裂解反应时发图苯热裂解反应路径设计现甲苯中断裂能最小的键是甲基C—H键其次是Figure 1 Design of toluene pyrolysis reaction path苯环与甲基相联的C—C键。经实验证实500℃左右最初断裂的键是甲基C—H键朕二甲苯是主要2的初期产物其反应机理如下TYH中国煤化工CNMHG体系时用半经验的收稿日期:201-06-26;修回日期:2001-1221基金项目:陕西省自然科学基金2001H03和陕西省教育厅重大产业化培育基金017004)联系作者:王惠西北大学化学系,Te(0298303270E-mil:xb@pub,kaoline.com作者简介万病|5)男河南南阳人硕士生物理无机化学专业。燃料化学学报30卷SCF-MO-AM方法进行非限制性 Hartree-Fock( UHF)计反应的△HF、△E、△E、△H、△C和△S。对热计算虽耗时短但结果有时不太可靠。若计算时考力学支持的主反应路径用UAM法进行了过渡虑电子相关效应则计算大分子体系的UHF/3-21G*态的寻找和优化并通过振动频率分析和RC方法方法又有一定的缺陷。而密度泛函理论在考虑电子对过渡态进行了验证得到了主反应的活化能。具相关效应的同时不仅可进行大分子体系的计算而体计算过程如下:于设计的热裂解反应为单分子且计算结果准确、可靠”。为了摸索适合低级芳香反应因此使单分子热解可能得到的自由基之间保烃类化合物尤其是分子结构较复杂的化合物热分解持一个大的距离约为正常键长的3倍~5倍作为反应机理研究的理论计算方法我们首先以〔s-产物构型用Qsn方法优化并经振动频率分析确sian98程序包中的半经验方法UAMI及从头算UHF定有唯一的虚频最后用RC法验证当ⅢRC法得到方法对各反应物及产物的几何构型进行优化。将优的反应物和产物与优化的构型和热力学量一致时化好的平衡几何构型作为初始输入构型再用考虑反应路径得到确认。全部计算在西北大学物理无机了电子相关效应的密度泛函方法UB3LYP3-21*研究所 Dell Dimension pll-550计算机上完成。进行进一步计算。对化合物中可能产生热裂解的化学键作了模拟均裂的热力学量及振动频率计算对3结果与讨论热力学量进行了振动零点能校正,计算得到了各设3.1优化的平衡几何构型参数图2为用密度泛HHJAO3436D723DasHIAtomic numberHgD, 254双7D.445)中国煤化工D910.14.ICNMHG图2反应物、中间体及产物优化的几何构型及部分几何构型参数键长m键角、二面角度)products at UB3LYP/3-21G* level( bond lengths, nm. Bond angles and Dihedral, deg1期乔占平等:甲基苯热裂解机理的理论研究函方法UB3LP/3-21G兴计算的A、B、C….H等9种比较反应路径1和2的难易即可。从前面表1的热反应物、产物及产物自由基的几何构型图中数据为力学量变参数Δ及ΔG的比较可知热力学结果部分键长、键角及二面魚D)数据。由图可看出反更支持反应路径1为主反应路径这是因为无论是应产物A的R28R7,15)的键长分别为0.1085从能量的角度还是从平衡的角度反应路径1的值m和0.1099m与相应的实验值o.1085m和均最低。因此,路径6比路径7和8更容易发生。0.1100m极为接近,说明密度泛函方法UB3LYP/3-这是因为路径产物H联二苯厢和路径&产物I2lG兴计算结果的准确性。将反应物A与自由基B、二苯甲烷的速控步是路径2这一步反应竞争不过C、D、E、F相应的键长比较发现由于自由基的生成生成(〔联二甲苯路径6舶反应。若温度升高路其键长均略有变化且变化方向变大或变小与预想径2的反应比例将增加即路径7和8的产物比例亦的变化方向一致将增加。这一结果与联二甲苯是主要的初期产物自由基B、C、DF、F是由反应物A热解断裂键且随着热解温度的提高(800℃左右)甲苯脱甲基的后生成的自由基因此自由基的对称性与反应物A比例将要多些同时也有少量苯环上的脱氢反应发相比亦有变化这可由键角和二面角的比较看出。生的实验结果相一致。3.2热裂解反应的热力学量变及反应机理表1,表1由 UB3LYP/321G法计算的反应路径1、2、3、4和2为UB3LYP方法在3-21G兴水平上计算的8种热裂5的标准热力学参数解反应路径的有关热力学数据。由表1可看出,计 Table 1 The standard thermodynamic parameters of the reaction算得到的△E。、ΔE、△H及ΔC其热力学量由小process of five types by UB3LYP/3-21G到大的顺序为路径1<路径2<路径3~路径4Reaction path路径5。由于ΔH为设计热裂解路径的反应标准焓△E0"/kmol)371.16414.19462.37462.49464.55变从能量的角度考虑,△P愈小反应愈容易进行;△E(m)1373.30417.26463.67460468.12ΔG为标准自由能变从平衡的角度考虑,△C愈小△Hk〔mol)375.78419.74466.14468.50470.60反应达平衡时的产物愈多。因此,无论是从能量上△G/k〔mol)1348.31376.0243972435.2743677看还是从平衡上热力学均支持热裂解路径1即甲/从ml)K10.09200.14650.0850.11020.134基苯的热裂解反应首先从反应路径1开始接下来△E0=△E+ZPE;△E=△E0+△Ea+△En+是反应路径2。这一计算结果与大谷杉郎等的甲苯△E1-;△H=△E+nRTG=△H-7s中甲基C—H键的断裂能小于苯环与甲基相联的表2反应路径67和8的标准热力学参数C—C键的断裂能的实验结果相一致。该计算结果也证实了随着温度的提高脱甲基的比例要多些的结Table 2 The standard thermodynamameters of the论。reaction path 6 ,7 and 8表2为路径67和8的热力学量变参数。路径K生成G)为路径1的接续反应潞径飞生成H)为Reaction path△E0△E△H△G/kK mol K)路径2的接续反应路径&生成Ⅰ)为路径1和2的6255.86-255.18-257.66-207.17-0.1734接续反应即苯甲基自由基和苯自由基的交叉反应732-410.13-0.1949由表2可看出3条路径的反应焓变△H均为负值91-298.34即自由基结合生成化合物G、H和I的一步都是放表3计算的热裂解路径1的活化能、活化焓、活化热反应,该热力学计算结果表明这3条接续反应极自由能和活化熵易发生。经仔细优化寻找3条反应路径的过渡态,Table 3 Calculated activity energy activity enthalpy发现其反应过程是一无过渡态的无能垒反应动力中国煤化工 activity entropy of the学计算结果亦表明3条路径的反应极易发生。由此CNMHG)可知从甲苯热解开始经路径1和2再至路径6、7 Activity ene和8生成产物G、H和I的3条反应的速控步应为k〔mol)E△H△G/k〔molK)路径1和路径2。也就是说,比较生成物G、H和IPath1283.9928776287.75267.630.0674的3条反、7和8的反应难易只需要通过燃料化学学报30卷表4反应物和部分自由基的实验标准生成焓Table 4 Experimental standard formation enthalpy of the reactant and free radicals 1-131Reactant andC Hs CH2C Hs CH, C Hs CHCH,CHs CHCHsCHa△H/mol)50.17±0.42204325.1±8.4143.2182.10142.3±5218.0°芳at300K由于本文设计的热裂解路徴(如路径1)为一个由表3可看岀,路径6速控步(路径1)的活化能分子裂解成两个自由基在这样的基元反应中无需△E0=283.9kJ/mol再形成新的化学键活化能即为被裂解的化学键的3.3反应焓变的计算值与实验值的比较我们以键能"。因此在这类反应中热力学量如反应的焓可查到的反应物及自由基的标准生成炫(见表4)变△H可当作键能。由于这一特点本文作者在研对反应路径12、67的标准反应焓变的实验值进行究低级芳香烃系列化合物的热解机理时发现没计了计算并将实验值与计算值进行了比较见表5的这一类热裂解反应路径过渡态是经频率分析确发现考虑了相关能的密度泛函方法 UB3LYP321G*认且再用BC认定的反应路径舶活化能的大小与的△H计算结果与实验值吻合较好其中最大误差反应的热力学量相近似所以往往热力学支持的路为7k/ml左右最小误差仅为0.22kJ/mdl左右说径亦是动力学支持的路径m。正因为此我们只对明密度泛函方法UB3YP3-2G*较适合低级芳香反应路径1作了动力学计算汁算方法为UAM法,烃热裂解反应的热力学量计算且数据比较准确可QST2优化过渡态相应的过渡态热力学量见表3靠表5反应路径1、2、6和7反应焓变的计算值和实验值的比较able 5 Comparing standard enthalpy change k, I/mol )between calculated and experimental data△H/kmol)(cal)△H/k(mol)(exp)UB3LYP/3-21G#375.78419.74-257.64-468.32371.83417.23-264.80-468.14结论结果首先支持反应路径1及其后续反应路径6为主〔1)将标准反应焓变的计算值与实验值比较,反应路径产物为联二甲苯。随着温度的提高按路发现考虑了相关能的密度泛函方法UB3LP/3-21G*径2的反应比例将增加即反应路径7和8的反应比的计算结果与实验值吻合较好说明密度泛函方法例将增加同时发生少量路径3、45的反应。该结UB3LYP/3-21G*较适合低级芳香烃热裂解反应的热果与实验结果相一致。力学量计算(3)AM1计算的主反应路径速控(路径1舶(2)甲基苯热裂解反应机理如下热力学计算活化能为283.99kJ/mod参考文献[1]石土金李宗和刘若庄.HNCO+OH=H2O+NCO的反应机理J]物理化学学报19913)247-252[2 Hui Wang Haifeng Yang Xinquan Ran et al. Pyrolysis mechanism of carbon matrix precursor cyclohexane II J ]. MOlecularStructurd Throchem),2001, 571: 115-131[3]杨海峰王惠段锦霞等.有机配体 Schiff碱42-羟基苯基)亚中国煤化T1无机化学学报,200172)215-22[4]王惠杨海峰翟高红等.碳源甲基苯热裂解机理的密度泛函εCNMHG1,591):17-21[5]朱卡克甘正汀苏克和等.CH的自由基和NO2反应研究J]化学学报200055)515-518[6]王惠翟高红杨海峰等·碳前驱体αHλHN热解杋理的热力学和动力学υr研究J]高等学校仳学学报2001,2x5)800-804[7] Becke bBH5Ey-functional thermochemistry. Il. The role of exact exchange J]. J Chem Phys 1993, 98(5): 648-6561期乔占平等:甲基苯热裂解机理的理论研究[8] Deutsch S Keieger K A. Free radical from the decomposition of toluene and butene -[ J ]. J Phy Chem 1962, 66( 19)1573[ 9] Bicout D, Field M. Quantum mechanical simulation method for studing biological system( les houches workship IM]nger- Verlag, 1995,1-2[10]傅献彩陈瑞华.物理化敩M]北京:人民教育出版社,1982.214-21511 Eugene S Domalski, Elizalzeth D Hearing. Estimation of the thermodynamic properties of hydrocarbons at 298. 15 KI J ]. J Phyem Ref Data988,174):1637-167812 J Sharon G Lias, Joel F Liebman, Rhoda D Levin. Evaluated gas phase basicities and proton affinities of molecules heats of formation of protonated molecule J ] J Phys Chem Ref Data, 1984,13 3)195-80813 David R I( Ed. ) Handbook of Chemistry and Physic M]. 63ed, 1982-1983, 55 ed, 1974-1975THEORETICAL STUDY ON PYROLYSIS MECHANISM OF TOLUENEQIAO Zhan-ping, WANG Hui', ZHAO Wen-Ii',RAN Xin-quan(1. Chemistry Department, Northwest University, Xi'an 710069,Chin2. Chemistry Department, Nan- yang normal College, Nan- yang 473061, ChinaAbstract: On the basis of the calculating results by UAMI and UHF/3-21G* methods , the pyrolysis mechanism oftoluene has been studied at UB3LYP/3-21G* level. The geometries of the reactant producted radicals and TransitionState were optimized by energy gradient method. The results showed that C-H bond of methyl on benzene will splitfirst at low temperature followd by the formation of bibenzil. The rate of benzene radical formation increases with theincrease of the temperature, and a few of C-H bonds on benzene will be broken at the same time. This mechanism isin accord with experimental result. In addition, the data of AH calculated at UB3LYP/3-21G* level coincide with theKey words: carbon matrix precursor toluene i pyrolysis mechanism UB3LYP/3-21G* computationFoundation item This work was supported by the Natural Science Foundation of Shaan-Xi province( 2001 H03 )and the breeding industri-alized fund of the education committee of Shaan-Xi Province 01Z004 ), which is gratefully acknowledCorresponding auther: WANG Hui, Department of Chemistry Northwest University E-mail: xb@ pub. xaonline comAuthor introduction: QIAO Zhan-ping( 1965.), male, Graduate, physical chemistry中国煤化工CNMHG