阿司匹林的热解机理及热动力学研究

物理化学学报( Wuli huaxue Xuebao)Acta Phys. -Chim. Sin., 2004, 20(1): 50-54January阿司匹林的热解机理及热动力学研究徐芬孙立贤谭志诚梁建国1周丹红2邸友莹兰孝征张涛(中国科学院大连化学物理研究所,材料热化学实验室,大连116023;1湖南省药品检验所,长沙410000辽宁师范大学化学系,大连116029)摘要在用热重法研究了阿司匹林的热稳定性实验的基础上,通过量子化学方法( ab initio dFt)计算了阿司匹林分子的键级,据此计算结果提岀了阿司匹林的热解机理,按此机理得到的理论计算值与实验结果一致;运用Freeman-Carroll Kissinger和 Ozawa三种方法分别计算了阿司匹林的热解动力学参数:活化能(E)、反应级数(n)和指前因子(A),其热解动力学方程为:da/dt=4.74×10lexp(-(100.34±5.18)×103/RT)(1-a)280用差示扫描量热法测定的该物质的熔点、摩尔熔化焓和摩尔熔化熵分别为(409.19±0.22)K、(29,.17±0.41kJ·mol-1和(71.09±1,06)J·mol-.K-1关键词:阿司匹林,量化计算,热解杋理和热动力学,热重法,差示扫描量热法中图分类号:阿司匹林(乙酰苯甲酸)是临床上应用的广谱1.1试剂与仪器药物之一.除作为解热、镇痛、消炎、抗风湿药物阿司匹林由湖南省药品检验所提供,其纯度符外,还在治疗心脑血管疾病、糖尿病、缓解晚期癌瘤合中国药典2000版标准症状等方面得到广泛应用12.关于阿司匹林及其Setsys16/18型综合分析仪,DSC-141型差式衍生物在体内的溶解动力学、释放机理及吸收过程扫描量热计,均为法国 Setaram公司制造等方面的研究工作目前有很多报导.但阿司匹1.2热重(TG)·微分热重(DTG)和差式扫描量林热稳定性方面的研究报导不多.阿司匹林的热热(DSC)曲线测量性能对于该药物的生产、稳定性的评价和临床应用热重曲线以AL2O3为参比物,在N2气氛中(N2都具有十分重要的意义.张健等人从热动力学方流速为30mL·min-1),控制升温速率分别为5面对阿司匹林的热稳定性及热解动力学进行了研10、15、20、25K·min-.样品池系容积为100μL究,但尚未见有关阿司匹林热解机理方面的报导.的Al2O3坩埚近些年来,人们从量子化学的角度来预测有机物的DSC曲线:N2流速为50mL·min-1,升温速率热解机理,得到了令人满意的结果.本实验室已为10K·min-;样品池为带盖的铝坩(容积为30对多种药物及材料的热力学特性进行了深入μL)的研究.本文在已有工作的基础上,通过量化计算,得到了有关阿司匹林分子的键级,用该计算结2结果与讨论果预测了它的热解机理,同时运用三种不同方法系2.1阿司匹林的热解机理统地研究了阿司匹林的热分解动力学参数,使用差测试了不同升温速率对阿司匹林的热重曲线式扫描量热法(DSC, differential scanning calorime-的影响.实验结果表明,其热分解过程包括了两个try)进一步测定了该物质熔化过程的热力学性质参阶段,升温速率改变,每个阶段所对应的温度范围数略有变化随养涅违潔的坦高,温度范围略向高温方中国煤化率基本保持不变1实验部分现仅NMH实率为例(见图1),其200306-13收到初稿,2003-09-15收到修改稿.联系人:孙立贤(E-mail: Ixsun@dicp.ac,cn;Tel&Fax:0411-4379213)家自然科学基金(20073047)资助项目No. I孙立贤等阿司匹林的热解机理及热动力学研究51DIG088°0973C100(用40050060070080090图2阿斯匹林分子中的化学键键级Fig 2 Bond orders of aspirin molecule图1阿斯匹林的TG-DTG曲线Fig 1 Plot of TG-DTG for aspirin匹林第一步热分解的动力学参数.设阿司匹林热分中TG为热重曲线,DTG为微分热重曲线,第一阶解的动力学机理函数为f(a)=(1-a)",由 Freeman段的温度范围为420~550K,失重率为33.67%;第Caro法可得方程式(1)、(2)二阶段的温度范围为550~720K,失重率为In(da/dt)+E/RT=In A+ nIn(1-a)(1)和ld dw dt为了了解阿司匹林的热解机理,我们首先试图△1-a4.575△|K1-a)从量子化学的角度进行理论探讨,并着重讨论第其中,a为t时刻物质已反应的质量分数;R为气个阶段的热解过程.本文应用 Cerius2软件包中ab体常数;T为绝对温度;E为活化能;A为指前因tio dft方法对阿司匹林体系的几何结构进行优子;n为反应级数.以△lg(da/dt)/△g(1-a)对化,并计算出分子的 Maver键级.优化计算采用局(1/7)/Alg(1-a)作图,应为一直线,由直线的域密度近似的方法( local density approximation,斜率可算出E;截距为见图3)LDA)基组选择DNP( double-numerical basis func然后用 Kissinger和 Ozawa方法再分别计tions together with polarization functions),该基组相算它的活化能E.根据 Kissinger法当于631· Gauss轨道基组.所有计算采用hn(B/T2)=n(AR/E)+ln[-df(a)/dalp-Cerius2软件中的基于密度泛函理论的DMol3模E/RTP块,在SGIO2工作站完成,结果见图2以ln(B/T2)对1/T(其中T为DTG的峰温,β为由图2可知在阿司匹林分子中C1-O3键的键升温速率)作图,由直线的斜率可求出E(见图4)级最小.因此,阿司匹林热解时,首先应从此断裂最后根据 Ozawa法失去CO2和0.5CHCH,生成苯甲酸.按此历程热解的理论失重率为32.21%,这基本上与热重曲线的第一阶段失重率33.67%相对应.因此可以推断阿司匹林的第一阶段的热解反应机理如下COOHCoOH20}y-240-546982XCrackingE=104. 6 k]. molcocH3420~550K25R2-093CO2,+0.5 CH, CH,104△(1/T)/△lgl-a)2.2阿司匹林热解反应的动力学研究图3中国煤化工4(1/)/△g(1-a)2.2.1热动力学计算的理论方法CNMHGFig3 root uruigluc/ut)ang 1-a)us采用了三种不同方法来计算该物质热分解的活△(1/T)/△lg(1-a) during decomposition化能(E).首先根据 Freeman-Carroll法计算阿司process of aspirinActa Phys. -Chim. Sin. Wuli Huaxue Xuebao ) 2004Vol 20-102}r=1523-1179228X9804R2=099949106图41n(B/T2)对T1的关系曲线图5ln(da/dt)+E/RT对ln(1-a)曲线Fig 4 Curve of In(B/T)us TEFig 5 Curve of In(da/dt)+E/RTus In(1-a)d(gB)/d(1/T)=-0.4567E/R(4)阿司匹林热分解的活化能(E),结果见图4.由图中式中T为在相等质量损失时,与升温速率β、B、的斜率可求E为98.04 k]. molB……相对应的温度然后,以lgB再用 Ozawa法计算E值,结果见表2.当a在对1/T作图得一直线,由其斜率可计算出E0.2~0.8范围内时,E值在90.11~107.32kJ而指前因子(A)可根据式(5舶关系mol-1范围内变动,其均值为96.71kJ·mol-,它与ln(da/dt)+E/RT=lnA+mln(1-α)(5)前两种方法得到的E值很接近,这表明阿司匹林热以mn(da/dt)+E/RT对mn(1-a)作图,由截距可分解动力学机理函数f(a)=(1-a)"的假设是成求得A,斜率为立的.由三种方法得到的E的平均值为(100.34±2.2.2热解反应动力学参数的计算5. 18)k]. mol-I以10K·min-的升温速率测得阿司匹林的图5为以ln(da/dT)+E/RT对ln(1-a)绘TG-DTG曲线用 Freeman-Carroll法计算其热分解制的曲线(E=100.34kJ·mol-),由截距求得A活化能(E)和反应级数(n)(见图3和表1).由表1为474×101;由斜率可得n也为2.8.将求得的可知,E为(106.28±1.7)kJ·mol-;n为2.8±A、E、n值代入动力学方程0.3.da/dt= Alex(- E/Rt)I(分别以5、10、15、20、25K·min-的升温速率得测得阿司匹林的 TG-DTG曲线,用 Kissinger法计算da/dt=4.74×10exp(-(100.34±5.18)103/RT)](1-a)2303(6)表1用 Freeman-Carroll法测得的活化能E和反应级数nTable 1 Determination of e and n based on the freeman-Dynamic parameter value Mean value Standard deviEndotherm2.80.3表2用0zawa法测得的TH中国煤化工SC曲线CNMHGSCfor aspirinTable 2 Determination of e based on the Ozawa methodE/kJ·ml91.0990.1191.8910.7898.31107.32No. I孙立贤等阿司匹林的热解机理及热动力学研究表3阿斯匹林熔化过程的热力学参数的测定值Table 3 Thermodynamic parameters in fusion process of aspirin△Hm/ k]. mol409.44414.67408.98414.7970.9409.11414.52Mean value409.19414.6729.1771.092.3阿司匹林熔化过程热力学参数的测定35(6):461阿司匹林的熔点为409K;但NST中未见hany0.MN: dssas.R.s. larn. Acta He熔化焓和熔化熵报导.本文用差式扫描量热法2000,74(4):3516 Chronos, N, A. F: Wilson, D. J,: Janes, S. L: Hutton, R. H.(DSC)测量了阿司匹林的熔化过程(见图6)参照Buller, N. P; Goodall, A. H. Clin. Sci., 1994, 87: 575热重曲线图1可知,图6中的吸热峰出现在失重之72h1n,;Ww1,s..;Cha,D.H.1 South central unit前,故该峰为熔化吸热峰.从DSC曲线准确测得其Nationalities( Nat. Sci. Edition ) 2002, 21(4): 21 [5t t相应的熔点(Tn)熔化焓(ΔHn);并根据热力学公土筠,陈栋华.中南民族大学学报( Zhongnan Mingzu Daxue式ΔSn=ΔHn/T计算出熔化熵(ΔSa),详细结果Xuebao),2002,21(4):21]见表3Chinese J. Inorganic Chem,2001,17:538[王惠,杨海峰冉新权,文振翼,史启祯.无机化学学报( Wuji Huaxue Xue3结论0),2001,17:538]通过量化计算对阿司匹林的热解机理进行了9Wang,x:m.xc.:1.X: uyang,:uo,MD理论探讨,并运用热力学方法对其动力学过程进行S.S. J. Wuhan Univ. Nat Sci. Ed.), 2001, 47: 399 [f4宪,颜肖慈,黎新,欧阳礼,罗明道,屈松生.武汉大学学报了基础研究,得出了阿司匹林第一步热解的动力学Wuhan Daxue Xuebao ) 2001, 47:3991方程为0 Wang, J; Gu, J D,; Tian, A M. Acta Chimica sinica, 2002, 60da/dt=4.74×l0lexp(-(100.34±5.18)×103/(4):590[汪敬,顾健德,田安民.化学学报( Huaxue xueRT)](1-a)2bao),2002,60(4):590]测定了阿司匹林的熔点、熔化焓和熔化熵分别为11 Xu, F: Sun, L, X; Tan, Z C,; Nan, Z D: Yu, P,: Zhang, T.J. Thermal Analysis and Calorimetry, 2003, 74: 335(409.19±0.22)K、(29.17±0.41)kJ·mol-和12 Nan, Z. D: Tan, Z.C.: Sun, L.X. Acta Phys. -Chim. Sin, 2002(71.09±1.06)J·mol-….K-18(10):947南照东,谭志诚,孙立贤.物理化学学报(WReferences13 Freeman, E S; Carroll, B. J. Phys. Chem., 1958, 62: 394I Altman, R: Scazziota, A: Cordero Fuenes, J. Thromb Res.14 Liu, Z H: Hatayama, L Z. Handbook of analytical chemistry1988,51:22 Zou, H Q: Gu, W.X. Guang Dong Pharma, 1998,(1)16[刘振海,龟山立子编,分析化学手册,热分析.北邹恒琴,古维新.广东药学( Guangdong Yaoxue),1988,(1):11京:化学工业出版社,2000:4613 Gilian, P. M: Mary, T K. Analytical Chemistry, 1998, 70(2)Science press,2001:108胡荣祖,史启祯.热分析动力学t Wang, J; Zhang, R. H. Acta Pharmaceutica Sinica, 2000, 35(6)北京科学出版社,2001:108461[王晋,张汝华.药学学报( Yaoxne xuebao)2000,THemEActa Phys. -Chim. Sin. Wuli Huaxue Xuebao ) 2004Vol 20Studies on Thermal decomposition mechanism and Kinetics of aspirinXu Fen Sun Li-Xian Tan Zhi-Cheng Liang Jian-Guo Zhou Dan-Hong Di You-Ying Lan Xiao-Zheng Zhang TaoMaterial Thermochemistry Laboratory, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023;Hunan Institute of Drug Detection, Changsha 410000; Department of Chemistry, Liaoning Normal University,116029)Abstract The mechanism of thermal decomposition of aspirin was studied by both thermogravimetry andMayer bond orders calculated by Cerius2 software. The parameters of thermal decomposition kinetics for aspirinsuch as activation energy( E), reaction order ( n) and frequency factor (A) were obtained by thermogravime-try. The kinetic equation of thermal decomposition of aspirin is expressed asda/dt=4.74×10texp(-(100.34±5.18)×103/R)1(1-a)23=3Melting point, molar enthalpy and entropy of fusion of aspirin were examined to be(409. 19+0. 22)K,(29. 17+0. 41)kJ. mol- and(71. 09+1.06)J.mol-I. K, respectively, by means of DSc with a rising temperatureprogram. The results obtained show that theoretical mass-loss derived from Mayer bond orders of aspirin is in verygood agreement with that obtained from thermogravimetry experimentKeywords: Aspirin, Quantum chemistry calculation, Mechanism and kinetics of thermal decompositionThermogravimetry, Differential scanning calorimetry(dSc)Received: June 13, 2003; Revised: September 15, 2003. Correspondent: Sun Li-Xian(E-mail: lxsun@ dicp ac, cn; Tel& Fax: 0411-4379213)The Project Supported by NSFC (20073047)中国化学会第24届学术年会征文启事中国化学会第24届学术年会定于2004年4月24-27日在湖南省长沙市举行.会议由中国化学会主办湖南大学承办湖南省化学化工学会协办大会主题为新世纪的中国化学——机会与挑战会议主席为白春礼院士,王柯敏任会议执行主席;姚建年、周梦君仼会议秘书长.会议共设21个分会绿色化学、纳米化学、化学生物学、化学信息学、有序超薄膜与超分子化学、分子磁体、有机固体电子过程与分子电子学、高分子化学与物理、生物有机药物化学、生命科学中的分析化学、电化学、光化学、复杂体系的分离分析、晶体化学、有机和高分子光电子材料、催化化学、理论化学和计箅与模拟化学、配位化学与生物无机化学、固体化学与物理无机化学、21世纪的应用化学、国防科技中的特种化学问题3个专题论坛牝学教育论坛、化学科普论坛、化学与社会论坛会议征集前瞻性综述、学科现状、前沿及展望、系统性研究成果及原创性硏究工作等论文.会议交流形式为大会特邀报告、分会邀请报告、专题报告与讨论、墙报展讲等,报名表及应征论文请直接寄给各分会联系人,有关会议详细内容、报名表、各分会地址、联系方式请从中国化学会网页下载网址为htp:!//www.cs/2004nh htm会议期间还将组织成果、产品以及仪器展示,包括实用技术交流产品、专业图书及试剂展示实验仪器、设备演示顶目洽谈狺息发布等.会议将设置展台提供信息服务.需要者请从速联系索取有关通知内容中国煤化工欢迎广大科技工作者踊跃投稿积极参加,欢迎相关企业高校科研院所利用会议注册费会员600元非会员850元,学生凭证减半不含博士后)年会CNMHG理入会手续中国化学会联系人刘钧地址北京2709信箱邮编:100080电话/传真10-62568157;电子信箱csj@icas