计算环烷烃沸点的新方法——基团键贡献法2∞2第22卷第3期.194~199有机化学Chine~Joamalof(…0…v0l22.2∞2I3.J94~199研究简报?计算环烷烃沸点的新方法——基团键贡献法王克强,王捷冯瑞英(.洛阳师范学院化学系洛阳471022)(许昌职业技术学院许昌461000)洛阳市第一高级中学洛阳471002)摘要根据分子中基团的特性和连接性.发展了一种直接根据分子结构信息计算环烷烃沸点的新方法基团键贡献法,该方法既考虑分子中基团的特性,又考虑基团之问的连接性(化学键),具有基团贡献法和化学键贡献法的特点对256种环烷烃沸点的计算结果表明,计算值与宴验值的一致性令人满意,平均误差071啦.关键词结构性能关系,基团键,基团键贡献法,沸点,烷烃GroupBondContribution'__.?.___-——ANewMethodforCalculating BomngPointsofCycloalkanesW ANG,Ke-QiangwA№,Jie.FENG,Rui—Yiag(mofChoTaist~,Normalc.,h~/ang471022)(~XuchangI.'ocatiormlTerhnical0f,461000)('.增sHigh&hodNoJ,4~002) AbslractOilcharacteristicsandconnectivi~,ofthegroupsinmolecules.groupbondscanbea ppliedtocharacterizemolecularstructure.Groupbondscontaininformationofgrouppropertyandco nnectivityinmoleculcs.Groupbonds,obtaineddirectlyfr(1堋molecularstructure,COIlbeusedtocNcuhWthepropertiesofmolecules.Anewmethod,thegroupbondcontributionmeth~,WItSdevelopedtocalculatetb eboilingpoints0fcycloalkmmsfrommolecularstructurelrhecalculatndresultsshowedthatthecaleutatedb oilingpointsof cycloalkaneswereingoodagreementwiththeexperimentaldata,withilloanrelativedeviatio nbeing0.71%for256cyeloalkanes(includingmonoeycloalkanesandpolyeycloalkmaes)egroupbondcontri butionmethodhasadvantagesoverthegroupcontributionmethodKeywordsatrueture-prepet~-relationship,groupbond,groupbendcontributhmmethod.bo ilingpoints,CVC】0日Ikane有机化台物结构性能定量关系的研究,是化学中'个十分活跃的研究领域.多年来,国内外许多研究者作了大量研究的工作,发展了许多根据分子结构预测化合物性能的方法.在有机化台物性能的预测方法中,应用较广的是基团贡献法和化学键贡献法,由于这些方法对同分异构体的区分能力相对较差,预测精度不够理想,作者之?4—6根据分子结构的特点,将基团贡献法和化学键贡献法结合在一起,发展了一种预测化合物性能的新方法一基团键贡献法,该方法具有基团贡献法和化学键贡献法的特点.可有效地预测化合物的性能.在有机化合物结构性能研究中,环状化台物的研究相对较为E-nHn:~'3558101@sdauO00t.~3558IOI@37InelRemlwdApril28,2001;re~sedSeptetnber17,2001;~ptedOctober22.ⅫNo3王克强等:计算:烷烃沸点的新方法摹囝键贡献法困难,无论是基团贡献法.还是化学键贡献法,都需要刘同的环进行修正,拓扑方法也需要对环状化台物分别处理.例如在环烷烃的研究中,拓扑方法对链烷烃和环烷烃的研究是分成陌类进行的,Estrade【提出计算环烷烃拓扑指数的方法.Schultz对单环烷烃的沸点与分子结构之间的关系作了研究,取得了较好的结果.本文用基团键贡献法发展一种计算环皖烃沸点的方法.对环烷烃沸点的计算结果表明.计算值与实验值的一致性令人满意1基本原理和方法结构决定性能是化学中的一条基本规律.化合物的性能不仅与分子中原子的种类和数El有关,而日更与原子的连接性有关根据前文的研究结果町知,分子结构可用基团键进行表征.显然,分子中基团的种类不同,基团键就不同;分子中基团的连接性不同,基团键也不同为简化起见,我们暂不考虑分子中环的大小,将环烷烃分子中的基团分为烷烃基团和环烷烃基团两大类,据此确定分子中基团键的数日对于烷烃,分子中有4种基团:分子中存在着1O种基团键.对于环烷烃,分子中有7种基团c其中包括烷烃分子的3种基团).其中三种新的基团为:为简化起,我们将环烷烃基团与烷烃基团之问形成的基团键作为烷烃基团键处理,仅考虑环烷烃基团之间形成的新基团键.同时根据环的大小进行适,的修正.环烷烃分子中有3种新基团,根据排列组台原理可以确定环烷烃分子中存在着C{+C;=6种新的基团键.表1列出了烷烃和环烷烃分子中存在的16种基团键,其中前l0种(C~C139)为烷烃基团键,后6种(GB『1_~GB』)为环烷烃基团键.分子中基团键的种类和数日可以直接根据分子结构确定,也可以根据分子结构的邻接矩阵和染色矩阵计算得到,具体计算方法见前文6.环烷烃与烷烃的主要差别在于前者分子中存在着环状结构,我们分两步计算环烷烃的沸点.第一步暂不考虑环烷烃与烷烃的差别,将环烷烃摹团烷烃基团处理,计算得到环烷烃沸点的近似值rlI(烷烃);第二步考虑环烷烃与烷烃的差别,计算环烷烃基团和环对沸点的影响,计算沸点的修JE量△r¨1这样即可计算得到环烷烃的沸点:(环烷烃)=Tb(烷烃J+AT1.(1)因此,只要根据分子结构计算,(烷烃)和△,即可得到环烷烃的沸点r(环烷烃).自先.我们探讨烷烃沸点T(烷烃)的计算.沸点是与分子结构密切相关的物理量,分了结构不同,沸点就不同既然分子结构可用基团键表征,那么,沸电必然为基团键的函数,因此.,(烷烃)满足:rlIc烷烃)=fl(GB)].显然,从理论上确定这函数是相当困难的,我们用数学方法进行近似处理我们发现,烷烃的沸点r(烷烃)随着分子巾碳原子数(碳原子数"与基团键之间存在着定量关系[TJ)的增加而增大,烷烃同分异构体之问的沸差异可用基团键定量描述,即碳原子数反映烷烃沸点的体变化趋势,基团键则从细节上反映r分子结构差异对沸点的影响.对烷烃沸点的拟台结果表明,沸点随着分子中碳原子的增加而增大,随着基团键的变化而改变,且烷烃的沸点1(烷烃)/n与基团键数n(GB)的线性函数和n(GB)/n显着相关,因此.我们选择的函数为b(烷烃)/=144.723+∑(+/)n(GB)…(2)式中为分子中的碳原子数,12.(GB,)为基团键GB. 的数目,日和C是与基团键特性有关的常数(基团键参数).表l列出了本文用753种烷烃沸点数据J 拟合得到基团键参数.其次,我们探讨环对沸点的影响.环烷烃中环对沸点的影响△r与环的大小及环上取代基的位置有关,前者对沸点的影响可根据环的人小进j修正, 后者对沸点的影响可根据环烷烃的基团键修正我们发现,环对△.的影响不仪与环的大小和数日有关,而且与分子中的碳原子数有关,在环的数目和大小一定的条件下,环的影响随着碳原子数的增加而增大;环烷烃基团键对沸点的影响同样与基团键的数目和分子中的碳原于数有关,且基团键的影响随着分子中碳原子数的增加而减小因此,我们选择的△函数为:l96有机化学V ol22.2002Bj)e-O.frS~c/n+哺㈥Tal~~靴2ParametersB.雌ofth团e式中(GBj)为环烷烃基团键GBj的数目,日是与基团键特性有关的常数(基团键参数),为环的数目.A为环的修正量.表2列出了本文用256种环烷烃沸点数据[拟合得到基团键参数和环修正量.表3列出了确定基团键参数所用的基团键数目. 表1计算Tb(烷烃)的基团键参数日和l'al~e1ParametemBi,Cofthegroupbondsforcalculating (alkane)2计算结果及讨论应用表2提供的基团键参数日,,c和环修正量△,可直接根据分子结构计算环烷烃的沸点.值得指出的是,在确定分子中环的大小和数目时,应注意各环应尽可能独立,不能存在环包含环的情况,即以取小环为主,所取的大环不能包含小环.例如,四环[2, 2,I,02一,os,]庚烷(图1)分子中的环为:1个四元环,2个三元环,1个五元环;三环r5,2,1,o2]癸烷分子中有3个五元环计算沸点的基本方法是:第一ringsforcalculatingATb下标C表示环烷烃基团之间形成的基团键表3确定基团键参数所用基团键的数目NTalde3Numberofthegroupbondstrsedincorre]atingthe paratt~etem丑andC代码Ⅳ代码Ⅳ代码Ⅳ代码N代码GB112492187G355GB】0737GB"9lGB2l5l91337G252GBl_494GB】d14G璐l015G57027GB】246GB】51由于G玮仅存在于乙烷中.其基团键的数目为H川H2圈1P耳环[2,2.1,o2~,os:庚烷和--~[5.2.1,0=]癸烷Figure1Quadrlcyclo[2.2,1.~,.]heptmr~andlricyclo[52,1.02,]deeane步暂不考虑环烷烃与烷烃的差别,将环烷烃基团作为烷烃基团处理,根据分子结构确定分子中基团键的种类,数目,代人式(2)计算(烷烃);第二步考虑环烷烃与烷烃的差别,根据分子中环烷烃基团之03王克强等:计算环烷烃沸点的新方法一基团键贡献法197 间形成的基团键及环的大小和数曰,代人式(3)计算△Th,将计算得到的I(烷烃)和△I代人式(1)即呵计算得到环烷烃的沸点.例如四环一2,2,l,.",03']庚烷,"=7,第一步,将环烷烃基团作为烷烃基团处理,分子中有2种基团键:n(G)=2,n(GB7)=8,代人式(2)得到Th(烷烃)=3O4.24K;第2步,考虑环烷烃基团的与烷烃基团的差别.环烷烃基团之间形成2种基团键:n(GB】1):2,n(GB】3)=8,分子中有1个四元环,2个三元环和1个五元环,代人式(3)可得△=82.5lK,代人式(1)得(环烷烃)=(烷烃)+△】,=386.75K;对于三环[5,2,1,02,]癸烷,,fr=10,第一步,将环烷烃基团作为烷烃基团处理,分子中有3种基团键:n(G)=3,(G)=6,n(GB7)=3,代人式(2)得到Th(烷烃)=339,05K;第2步.考虑环烷烃基团的与烷烃基团的羞别,环烷烃基团之间形成3种基团键:n(GB)=3,n(GB】1)=6,n(GBl)=3,分子中有3个五元环,代人式(3)可得△Tb=120.41K,代人式(1)得I(环烷烃)=(烷烃)+△TI=45946K.为确定本文方法的可靠程度和适用范围,我们对环烷烃(包括环丙烷到环二十烷)进行计算,并与实验值I进行比较.限于篇幅,表4列出了部分环烷烃沸点的计算结果,全体256种环烷烃的计算结果绘于图2.表4部分环烷烃沸点it-算值(ca1.)和实验值Tb(exl~)的比较Table4CⅥ)aI1s0nofthec~leulatedboilingpointTbmt}Iexperime.t.1dataofsoⅡcyclOalka~es环烷烃Tk,(exp)Th(ca1)E1%环烷烃(唧)Tb(eat)£/%甲基环丙烷27388271.88—0.73戊基环烷4768247587—0201一顺一2-二甲基环丙烷31018310.81020辛基环己烷5367553303—0.69异丙基环丙烷154340.57272十一烷基环己烷58635581.60—0.811一甲基一反2-£基环丙烷331.8133714161十四烷基环己烷6281562366—071 1顺一2一反一3三甲基环丙烷332.853H430.48十七烷基环己烷664.1566o58—054己基环丁烷34:375344.300.16=十烷基环己烷6951569332—0261一反一2-二甲基环丁烷33315337.30l24二十三烷基环己烷72215722.560.06秫-1'1AIA1nm+"1…AA…丙基环戊烷己基环戊烷壬基环戊烷十二烷基环戊烷五烷基环戊烷十八烷基环戊烷=十一烷基环1萱烷二四烷基环戊烷=十七烷基环戊烷三}烷基环戊烷三十三烷基环戊烷三十六烷基环戊烷己摹环戊烷】一反_2'二甲基环戊烷丙基环戊烷l一甲基一顺一2一乙基环戊烷】一甲基一反一3一乙基环戊烷1一』匝一2一顺一3一三甲基环戊烷l一顺一2顺一4-i甲基环戊烷丁基环戊烷特丁基戊烷1一甲基一反一2一丙基环戊烷41)4.1O476055351558435626.15662.1569'3157201574515767】57s7.15805.1537662 365024O4.10 401.20 39425 39615390.15 429754l800419.5241)447O.09 473】8—0.60 529.82一1.00 5_77.91—1.10 61952—106 65601—093 68834—069 717.】9—041 743.09—028 7054.5—009 787610068o6阻0.21 377.86033 37083】59 404.47009 39856—0.66 394.980l9 39150一l17 38793—0.57429.03—0.17422.5811O422.55072=十九烷基环己烷三十=烷基环己烷三十五烷基环己烷甲基环己烷1一顺一2一=甲基环己烷1一反一3-二甲基环己烷l雨基环己烷1甲基一顺2-乙基环己烷l一甲基一反一3一己基环己烷1,1,2一三甲基环己烷】一顺0丽一3一三甲基环已烷1一顺2一顺4三甲基环己烷l一反2一反4三甲基环己烷二环己甲烷1,】一二环己丁烷l_l一二环己庚烷1.1一二环己癸烷1.1一二环己十三烷】,】一二环己十六烷环千烷环十二烷环十五烷68.15788l5妍15374084O2939760 42987 42925424】5 418.35426l54211541437525.95 566.12605l564.O15 670.15 697.15 451555l2l5559】5772O58 794080758】3.640.8O 37806106 39953—084 39596—041 43】350.34 424.87一1.o2 421.63一O59 4】575—062 4l779一】96 414.55一1.57 414.5500452489一O205690305l608.9l0.62643.42O516736】()527002】04444896—0.57510.75—02755537一O68l98有机化学V o1.22,2g02续表环烷烃TbE/%环烷烃To7E/%(exp)(ca1)(exp.)(ca1)1一甲基一1一异而基环戊烷42l1542319048环十八烷5981560246()721一甲基一顺一3一异丙基环戊烷4151541717049联环庚烷56415565870311圳2一二己基环戊烷4267l42315~083环戊基环己烷4882548244—119I一反一3一二乙基环戊烷423154199l~0762,4,6-三甲基二环[3,1.1]庚烷4426543597—15】l一顺一2一二甲基一1一乙基环戊烷41715417740142一甲基二环[2,2.2]辛烷43l15435541O21一顺一2一一甲基一反一3一乙基环戊烷418.15416.07~050二环[2.1.O]戊烷318.653155l一0981一顺一2一二甲基一顺4乙基环戊烷420.1541284~17,4环『1n]_丁辟1'11钔一n11一顺一3一_I甲基一反一2一乙基环戊烷410.15416.07I44四环[2,2.I,.03]庚烷381.15386.75I.471一顺一3一二甲基反4乙基环戊烷410.15412.84065顺一一环[3.3,0]辛烷410.15412300.521.1,2,2_四甲基环戊烷496.15405.27~0.22反.=环[3,3.O]辛烷405.15412301.76 1.1-顺一2_顺一4一四甲基环戊烷403.15403.720.142一甲基二环l2,2,2辛烷43115435541.o21,1一顺一3一顺一4一四甲基环戊烷496.154.9O.89~I.313二环4,2.0一辛烷4091540892-一0.嘶1,2,2-反一3一四甲基环虎烷4l1.15409.79~0.33二环4.1,0.庚烷3896538017—2.43I一顺一2一反一3一婀4四甲基环戊烷4O415408.991.20三环:5,2,l,02癸烷4661545946—1.44I一反2一反一3一顺4四甲基环戊烷4O4154O8991.20螺[2.2]戊烷3131530746一12己基环己烷4O493406.5904l顺一六氢茚44O1543886—0.29表中环烷烃是用计算机选取的,方法是:对256种环烷烃编号为No1~256,选取No. 2,N05.8,….约占总数的三分之一!II2环烷烃沸点计算值Tb(c)与实验值(.xp)的比较Figure2(_^】m90nofthecalctrlatedboihngI~intTh(ca1)htheexperimentaldatars(e~p)ofcycl0alkane对256环烷烃的计算结果表明,无论是对含有一个环的简单环烷烃,还是对含有多个环的复杂环烷烃,沸点的计算值与实验值".]的一致性令人满意,计算结果的平均绝对误差3.29K,平均相对误差0.71%在计算的256种环烷烃中,计算误差小于0.5%的106种,计算误差在0.5%~l0%之问的90种,计算误差1.0%一1.5%之间的38种,计算误差1.5%~2.0%之间的l2种,计算误差2.0%~2.5%之间的6种,计算误差2.5%~3.0%之间的2种.计算误差大于3.0%的2种,最大计算误差3.13%.显然,绝大多数化合物的计算误差在l0%以内.文献上应用较广的计算烷烃沸的方法是Jobaek方法",该方法对上述256种烷烃计算结果的平均绝对误差46.89K,平均相对误差7.1】%.显然,与Joback方法相比,基团键贡献法将计算精度提高了一个数量级.因此,基团键贡献法具有适用范围,..,计算结果可靠的特点.此外,基团键贡献法是根据分子中基团的特性和连接性发展的,易丁实现计算机程序化设计,我们开发了相应的计算机程序,只需输A少量分子结构参数即可计算得到沸点预测值3结论本文根据分子中基团的特性和连接性,发展了一种根据计算环烷烃沸点的新方法——基团键贡献法,该方法既考虑了分子中基团的特性,又考虑了基团的连接性,同时具有基团贡献法和化学键贡献法的特点,对环烷烃(包括单环和多环烷烃)的计算结果表明,沸点计算值与实验值的一致性令人满意.No3王克强等:计算环烷烃沸点的新方法一基团键贡献法l99 1Reid,R.C;Prausrdtz,J.M;Poling,BE.Properties‰andt/qt,i&,4thEd.,McGraw-HillBook Company.NewY0rk.1987.2Xu,ZH;Mao,Z.X.;Wang,LS;Pang,YLHandbo&Chemi*'alProperey"EstimationMeth,~,Huaxue Got~,gyeChiCle.Bing,I99l(inCtfine~)(许志宏.毛卓雄,王乐珊,鹿瑶琳译,亿学性质估算方法手册,化学工业出版社,北京,1991.)3Wang,K.Q;Hu.YChin,.忡Ctwm.1997,17,230(inChinese)(王克强.胡英,有机化学,1997.t7.230.)4Wang,KQJ~.jiHngHuax,w1996,13.182 (inChisese)(土克强.计算机与应用化学,1996,17,182.)5Gutman,I;Estrade,E.JChemComputSci1996.36.5416schuhz.HP:Schultz1993.33,2407Wang,K.Q.Ch/nJChit~se)PjCheml《0唧nOrgCtwm.1999,19,304(in(于克强,有机化学,1997,19,304)8:~ericanPetroleumhmfituteResearchProject44,Seleaed pⅡofl#opertyofH)~warbonsandRelatedCom~unds1lrmo[hIacResearchCenterdtheAmericanPetroleum hrstitute,Texas,1,Tablea9Weast,RCCRCHandbookChen~tO,andm.66thFA.CRCPress.h,e..Horida,l986,04210Dean,J.Ahinge'HandbookofChemic-.13thEd. McGraw-HillBookCompany,NewY ork.1赠5.PP.7~82. 】lMa,PS踊uHuagong肛^u卧咖S,%uce,Huaxue Gon~,eChubami~e,l~ijing,1993(inChinese).(马肺生,石油化工数据手册,化学工业出版社,北京, 1993,(Y0104281LI,L.T:F,YY)。