二氢卟吩e6赖氨酸酰胺稳定性及电子光谱

二氢卟吩e6赖氨酸酰胺稳定性及电子光谱 二氢卟吩e6赖氨酸酰胺稳定性及电子光谱 32 物理化学(WuliHuaxueXuebao) ActaP.一Chim.Sin.,2011,27(1):32-38January [Article]www.whxb.pku.edu.cn 二氢卟吩e6赖氨酸酰胺稳定性及电子光谱 万景柏'刘淑贞'潘树英'许旋,李旭'叶冰 1华南师范大学化学与环境学院,广州510006;广东省高校电化学储能与发电技术重点实验室,广州510006; 华南师范大学,电化学储能材料与技术教育部工程研究中心,广州510006) 摘要:采用密度泛函理论在B3LYP/6.31G水上对二氢卟吩e6(简称e6)及设计的6个e6赖氨酸酰胺进行几 何优化,对优化构型用B3LYP/6-31G"法进行单点计算,并用含时密度泛函理论LSDA/6—31G"方法计算电子吸 收光谱.结果表明,赖氨酸的,一NH与e6连接的酰胺更稳定,其中,15位的乙酰胺Y最稳定.形成赖氨酸酰胺改 善了e6的水溶性,有利j药物吸收.各e6赖氨酸酰胺的前线轨道集中于=氢卟吩环,由于连接酰胺基侧链的 二氢卟吩环平面性有所下降,前线轨道能隙略为升高,最大电子吸收波长相对于e6蓝移16—39nm,但仍处于 光动力治疗窗口"600—900nm".酰胺链的构象对吸收波长影响较大,Yj个较稳定构象中,酰胺链垂直于二氢 卟吩环的YE和Y的二氢卟吩环血性较好,最大吸收波长比酰胺链与二氢卟吩环近似平面的Y红移53,5O nm,三者平均值较e6红移18nm. 关键词:密度泛函理论:二氢卟吩e6赖氨酸酰胺:电子吸收光谱;稳定性 中图分类号:O641 StabilityandElectronicSpectraofChlorine6LysineAmides WANJing—Bai'LIUShu—Zhen'PANShu-YingXUXuan'_-.,LIXu'YEBing' fSchoolofChemistry&Environment,SouthChinaNormalUniversity,Guangzhou,51000~R.China;2KeyLaboratoryof ElectrochemicalTechnologyonEnergyStorageandPowerGenerationinGuangdongUniversities,Guangzhou51000~PR. China;EngineeringResearchCenterofMaterialsandTechnologyforElectrochemicalEnergyStorage, MinistryofEducationofChina.SouthChinaNormalUniversity,Guangzhou510006,PR.China) Abstract:Thegeometricconfigurationsofchlorine6andsix usingdensityfunctionaltheoryattheB3LYP/6—31Gleve designede6lysineamideswereoptimized Basedontheobtainedminimumenergy structure.asinglepointcalculationwascarriedoutattheB3LYP/6— 31GIeve1.Theelectronicspectraof thesecompoundswerecalculatedusingtimedependentdensityfunctionaltheoryattheLSDA/6—31G" leve1.Theresultsindicatethatthee6lysineamides.inwhichthecarboxylofe6connectswiththe,一NH2of lysine,aremorestable.Amongthem,the15-acetamideYhasthehigheststability.Theformationoflysine amidesimprovesthewater— solubilityandleadstoaslightlypoorcoplanarityofthechlorinmacrocyclein e6.Therefore,thefrontierorbitalgapsofthee6lysineamidesareslightlyhigherthanthatofe6,causing theIongwavelengthabsorptionloasmallblue— shiftof16-39nm.TheadsorptionwavelengthisstilIjnthe rangeof600— 900amforphotodynamictherapyFu~hermore,thelongwavelengthabsorptionisstrongly affectedbytheconformationofthemolecule.BycomparisonwithY.inwhichthelysineamidegroupand thechlorinmacrocyclearealmostcoplanar.theplanarityofthechlorinringofY1andY,2,inwhichthe lysineamideisalmostperpendiculartothechlorinring.jsimprovedandresultsinredshiftsof53nmand 50nmfortheirlongwavelengthabsorptions,respectively,theaverageadsorptionwavelengthofYe,Y,1, andY口is18nmlargerthanthatofe6. Received:July24201O:Revised:October26.20l0:PublishedoilWeb:December9,2010. 'Correspondingauthor.Email:xuxuan@scnu.edu.cn;Tel:+86.20—852l6889. TheprojcotwassupportedbytheProjectofMinistryofEducationofGuangdongProvince,China(2ol0B090400184)andNaturalScience FoundationofGuangdongProvince.Chinar9l510631O10000371. 广东省教育部产学研项目(2010B090400184)和广东省自然科学资金项目 (91510631O1000037)资助 (0EditorialofficcofActaPhysico—ChimicaSinica NO.1万景柏等:二氢卟吩e6赖氨酸酰胺稳定性及电子光谱33 KeyWords:Densityfunctionaltheory;Chlorine6lysineamide;Electronicabsorptionspectrum Stability 二氢卟吩e6(chlorine6,简称e6)是光动力学疗 法(PDT)治疗恶性肿瘤的光敏剂之一.其机理是: 光敏剂或其代谢产物聚集于靶组织,在适当波长光 的作用下,将恶性组织周围的三重态o分子激发成 单重态氧而诱导肿瘤细胞坏死.e6被认为是第二 代新型光敏剂,与第一代光动力抗肿瘤药物如血卟 啉衍生物(HPD),光敏素II(photofrinII)等复杂的混 合卟啉制剂相比,具有结构明确,红光区(>600am) 吸收系数高等优点,且肿瘤光动力疗效更理想, 对肿瘤的选择性较高.此外,e6还具有亲脂性和一 定的亲水性,稳定性有所提高,但亲水性还不足以提 高药物的吸收效果. 为得到更稳定,选择性更高,疗效更好的e6衍 生物,已合成了一系列酯类,酰胺类,位取代e6衍 生物及其Zn,pd2,Sn配合物],测定了化合物的 药理活性.如方瑛等研究了e6的C15位单甲酯的 抗肿瘤活性,O1shevskaya等报道了e6酰胺的 Pd(II),Sn(IV)~HZn(II)等金属配合物可提高光敏剂 在靶细胞内的积聚,增加单重态氧量子产率而改善 光敏剂的光动力疗效.由于氨基酸是人体必需的营 养物质,与肿瘤组织有较强的亲和性,氨基酸与e6 结合的酰胺有望提高化合物的稳定性和对肿瘤组织 的识别能力.据报道,e6单天冬氨酸酰胺(NPe6)" 为有希望的光动力敏化剂,其光谱吸收波长为664 nm,对小鼠纤维肉瘤具有明显的光动力治疗效果, 对皮肤癌及早期肺癌具有一定的疗效,已进入II期 临床试验.因此,设计和开发其他e6氨基酸酰胺有 望得到更理想的光敏剂. 目前对共轭程度和对称性较高的卟啉类衍生物 光敏剂的结构和电子光谱的理论研究[14-161较多,但关 于二氢卟吩衍生物的理论研究还很少.Kuzmitskii 等用CNDO/S方法结合实验考察了二氢卟吩P6 及其衍生物的电子结构和荧光淬灭的原因:徐志广 等报道了位取代基对e6吸收波长的影响.本文设 计了6个e6赖氨酸酰胺分子结构,用密度泛函理论 对其稳定性和光谱性质进行理论研究. 1计算方法 1.1化合物分子模型 目前尚未有e6晶体结构的报道.将e6三甲酯A 晶体结构的三个酯基改为羧基,建立e6分子初始 L何结构模型.对e6的三个羧基分别连接赖氨酸, 设计了6个e6赖氨酸酰胺分子模型,如图1所示:赖 氨酸的一NH2分别与e6的甲羧基,乙羧基,丙羧基连 接形成的酰胺J,Y,B:赖氨酸的s—NH:分别与e6的 甲羧基,乙羧基,丙羧基形成的酰胺J,Y,,B,. 1.2方法 采用密度泛函理论(DFT)的B3LYP方法,分 别用6—31G和6.31G基组对e6,e6三甲酯进行几何 优化,发现两种基组优化的几何差异很小.因此,为 节省计算机时用B3LYP/6.31G法优化所有分子的 几何结构.为提高计算精度对各优化几何结构用 B3LYP/631G"法进行了振动频率和自然键轨道分 析(NBO)的单点计算;并采用含时密度泛函理论 (TD—DFT)LSDA[n1/6—31G"方法计算各分子的电子 吸收光谱,考察前30个激发态的能量,吸收波长和 振子强度等.所有计算采用Gaussian03软件包f231.此 外,用HyperChem程序的原子参数法计算了各 分子的油水分配系数1值. HcH3 H CH. A:R1=R2=l=OCH3e6:RI:R2=R3=OH B.:R1=NHCH(COOH)CH2(CH2)3NH2B:RjNI-I(CH2),CH(NH)COOH R2==OHR2=R3=OH Yn:R2=NHcH(cooH)cH:(cHz)NHzY:R2=NH(CH2)4CH(NH:)COOH R1=R3=OHRl=Rs=OH L:RsNHCH(COOH)CH2(CH2)3NH2J:RsNH(CH2),CHrNH)COOH R1=R2—OHRj=R2=OH 图1二氢卟盼e6,e6三甲酯和e6赖氨酸酰胺的分子结构 与原子编号 Fig.1Molecularstructuresandatomiclabelsofchlorin e6.e6trimethylande6lysineamides 34ActaPhys.一Chim.Sin.2011Vo1.27 2结果与讨论 2.1e6赖氨酸酰胺的几何结构 ,6.31G和6-31G"基组优化的e6三 由表1可见 甲酯A的结构参数均与实验值A..l】的较接近,其 中,6.31G"基组计算的二氢卟吩环的N—N距离更 接近实验值,但两种基组优化的A和e6的结构没有 显着差异,故采用6.31G基组优化所有分子的结构. 用6—31G'基组优化的e6二氢卟吩环的主要键 长与e6三甲酯A的偏差为0.0000—0.0053am,主要 键角偏差0.2.1.4.,环上桥头碳与邻近原子问的二 面角及四个N原子的二面角偏差0.6.一1.9..由于e6 的13,l5,17位侧链不含甲酯基,空间位阻较小,故 二氢卟吩环平面性较高,二面角N23C14C15C16仅 为1.1.,低于A的12.6..因此,用6.31G基组计算的 e6二氢卟吩环结构较合理,与e6三甲酯结构也很 接近. B3LYP/6—3lG方法优化的各分子构型用 B3LYP/6—31G"方法计算的振动频率均无虚频.与e6 相比,计算的6个e6赖氨酸酰胺中,环中央的N—N 距离偏差为0.0001—0.0042rim,二氢卟吩环主要原 子问键长偏差(0.0000—0.0007rim)较小,桥头碳与临 近两个原子构成的键角偏差也只有0.1.1.0.,桥头 碳与邻近原子的二面角及四个N原子的二面角偏差 仅02.一3.0..但因13,15,17位侧链分别形成了赖氨 酸酰胺,空间位阻增大,使二面角N23C14C15C16为 10.5o13.9o,高于e6的1.1.,但与A的12.6.较接近, 表明形成赖氨酸酰胺使二氢卟吩环的平面性略为降 低.由于Y的酰胺链与环近乎垂直,与17,13位侧链 的排斥较小,所以其二面角N23C14C15C16较小与 e6相差不大.各e6赖氨酸酰胺异构体二氢卟吩环的 结构差异主要在酰胺链上. 2.2e6赖氨酸酰胺的稳定性和亲水性 由表2中B3LYP/6.31G"方法计算的零点振动 能校正的分子能量可见,赖氨酸的s.NH形成的e6 酰胺J,Y,B,比一NH:形成的e6酰胺L,Y,B能量 低0.002—0.012a.U.,表明以s—NH:结合的e6赖氨酸 酰胺更稳定.赖氨酸一c0OH的吸电性虽使6c—NH 中N原子的负电荷密度一0.899e略小于s—NH 的一0.904e,但其负电荷密度差异不大.s.NH:形成的 酰胺更稳定的原因可能是空问位阻较小的缘故.e6 的l3,15,17位连接甲羧基,乙羧基和丙羧基,有一 定的空问位阻,赖氨酸6c.C除NH外还有COOH,使 .[.NH:与e6形成酰胺时空间位阻较大,因此.NH:形 成的酰胺的稳定性较低. s—NH与e6结合的酰胺J,,B中,乙羧基形成 的酰胺能量最低,最稳定.这可能是由于l5位的 酰胺处于l3位的甲羧基和17位的丙羧基之间,酰胺 基上的N,0,H与13,17位羧基的H,0形成分子内 氢键,使分子能量降低. 由表2可见,考虑了原子的杂化形式,电荷密度 表1B3LYP/6.31G,B3LYP/6.31G"水平上优化的e6,e6三甲酯及e6赖氨酸酰胺几何 参数 Table1Selectedgeometricalparametersofchlorine6,e6trimethylande6lysineamidesopti mizedattheB3LYP/6-31G andB3LYP/6—31Glevels calculatedattheB3LYP/6-31G"level;calculatedattheB3LYP/6.31G'level No.1万景柏等:二氢卟吩e6赖氨酸酰胺稳定性及电子光谱35 表2B3LYP/6-31G"水平上计算的e6~e6赖氨酸酰胺的分子能量(a.II.)和油水分配 系数 Table2Energy(a.n.)andoil-waterpartitioncoefficientofchlorine6ande6lysineamidescalc ulatedat theB3LYP/6.31G"level 和原子的溶剂可及性等因素的原子参数法[24-25计算 的各e6赖氨酸酰胺异构体的油水分配系数l均小 于e6的,表明引入含,C00H和一NH:等亲水基团 的赖氨酸使e6赖氨酸酰胺水溶性比e6有所改善,有 利于药物的吸收. 2.3电子光谱 2.3.1电子结构与分子轨道 化合物的光谱性质与分子电子结构密切相关, 为探讨电子跃迁实质,将e6赖氨酸酰胺分成4个分 子片,各片原子对前线分子轨道的贡献:二氢 卟吩环C,N,H原子(D),乙烯基(0),环上的甲基和乙 基(丫),13,15,17位侧链(6),结果列于表3. 从表3的前线分子轨道组成和图2的轨道图 可见,e6及e6赖氨酸酰胺的前线轨道(HOMO, LUMO,HOMO—l,LUMO+1)的电子云与卟啉类化 合物的前线轨道一致,均集中在二氢卟吩环上, 13,15和17位侧链对前线轨道的贡献不大,这表明 形成赖氨酸酰胺对e6的前线轨道分布没有明显影 响. 与e6相比,因二氢卟吩环的平面性略为降低, 离域效应略为减小,由图3可见,6个最稳定的e6 赖氨酸酰胺的HOMO,LUMO能级均稍高于e6. LUMO能级升高的幅度稍大于HOMO,因此LUMO 与HOMO能隙?略为增大(Y除外).这表明e6赖 氨酸酰胺在热力学上较e6稳定,但其最大吸收波 长将略小于e6.由于Y的15位酰胺链近乎垂直于 二氢卟吩环,降低了空间位阻,使二氢卟吩环的平 面性较好,离域强,故其HOMO与LUMO能级差 ?小于e6. 2.3.2电子吸收光谱 与卟啉类化合物相似,e6有两个强吸收带刮: 约400rLrn短波处的强吸收峰和640—660D1TI长波处 吸收峰.由于光动力治疗采用600—900am范围(称 为光动力治疗窗口)的吸收光,故长波处吸收峰对光 动力治疗和诊断有重要作用.用LSDA/6—31G"方法 计算的e6及其赖氨酸酰胺的电子吸收光谱列于表4. 根据B3LYP/6.31G方法优化的几何构型,用含 时密度泛函理论B3LYP/6—31G"和LSDA/6—31G"方 表3e6及e6赖氨酸酰胺的前线分子轨道组成(%) Table3Frontiermolecularorbitalcomponents(%)ofchlorine6ande6lysineamides Thefrontiermolecularorbitalcomponentsofchlorinering(p),vinyl(0),sumofthemethyland ethylonchlorinering,andsunlofthethree sidechainsin13,15,and17positions(6).HOMO:thehighestoccupiedmolecularorbital;LU MO:thelowestunoccupiedmolecularorbital; M0:molecularorbita1 ActaPhys.一Chim.Sin.2011Vo1.27 鼙.攀罄. 攀攀 YHOMO一1HOMOLUMO 图2Y,Y和e6的前线分子轨道图 Fig.2FrontiermolecularorbitalsofY.,Yqande6 法计算e6的电子吸收光谱.B3LYP/6—31G方法计算 的e6短波处的吸收波长392am与实验值402nm较 接近,但其最大吸收波长584nm与实验值643nm 相差较大.而由表4可见,LSDA/6—3lG"方法计算的 e6的长波处和短波处的吸收波长644和386nm均 与实验值很接近.因此,选用LSDA/6—31G"方法计 算所有分子的电子吸收光谱.计算的e6赖氨酸酰胺 电子光谱与e6的相似,保留了e6的特征吸收峰.6 个最稳定的e6赖氨酸酰胺的最大吸收波长位于 605—627am,相对于e6略为蓝移,但仍处于600一 表4LSDA/6.31G"水平上计算和实验得到的e6及e6赖氨酸酰胺的电子光谱数据 Table4ExDerimentaiandcalculatedelectronicspectrumdataofchlorine6ande6lysine— amidesatthe— LSDA/6-31G'"level theoscillatorstrength;H:HOMO,L=LUMO NO.1万景柏等:二氢卟吩e6赖氨酸酰胺稳定性及电子光谱37 图3LSDA/6.31G"水平上计算的e6及e6赖氨酸酰胺分子 轨道能级图 Fig.3Energylevelsoffrontiermolecularorbitalsof chlorine6ande6lysineamidesattheLSDA/6-31G'level 900nlTl范围,不影响光动力治疗.由于e6赖氨酸酰 胺的LUMO与HOMO能隙?较e6的略大,故其 长波处的吸收波长略小于e6.虽然Y的略小于 e6,但因其主要跃迁不含HOM0一LUMO跃迁,故 其吸收波长仍小于e6.e6及e6赖氨酸酰胺的最大吸 收主要由HOMO—LUM0跃迁产生.由表3和图2 可见,e6赖氨酸酰胺的HOMO—l,H0MO,LUMO 和LUMO+1轨道主要集中在二氢卟吩环,所以长波 处吸收主要由二氢卟吩环的跃迁产生.e6强度 最大的386nm吸收峰与实验值402nm接近,其赖 氨酸酰胺在该处的吸收峰377—427nlll也位于400 nlTl附近.该处吸收涉及多组跃迁轨道. 2.3.3构象对吸收波长的影响 e6赖氨酸酰胺中酰胺链的构象对电子吸收波 长影响较大.以最稳定的Y为例加以分析:以酰胺 链与二氢卟吩环共平面为初始构型优化得到最稳定 的Y;以酰胺链与二氢卟吩环垂直,且酰胺链在二 氢卟吩环的上下侧的两个初始构型优化得到Y.和 Y,所有结构的振动频率均没有虚频.由表2可见, Y最稳定,经零点振动能校正的能量比Y,.和Y低 0.011和0.008a_u..从图4可知,中酰胺链与二氢 卟吩环近似平面且伸向17位的丙羧基方向,-NH 的N8与丙羧基的H9形成较强的氢键,N8一H9距 离仅O.1646illn,且N8H9010键角为166.0..分子内 较强氢键的形成使Y的能量最低.优化后Y,和 的结构很相似,酰胺链均基本垂直于二氢卟吩环且 位于与丙羧基相反的一侧,两者能量略高于Y也可 以稳定存在,H7与甲羧基的O12距离为0.1870和 0.1872nm,其N6H7012键角均为159.0.,也能形成 分子内氢键,但比Y,中的弱,因此Y.和Y也较稳 定.此外,由于Y和Y的酰胺链垂直于二氢卟吩环 且与丙羧基位于二氢卟吩环的相反方向,使C13, C15,C17三个侧链的空间位阻较小,二氢卟吩环的 平面性较高,其二面角N23C14C15C16仅为一1.0o 和一0.7.,远低于Y的11.2.. Y,Y.和Y三个构象的HOMO,HOMO一1, LUMO和LUMO+1前线分子轨道均主要集中于二 氢卟吩环,由于YE.和Y.二氢卟吩环的平面性较高, 图4e6赖氨酸酰胺Y(a)和Y(b)的结构 Fig.4Structuresofe6lyslneamideYE(a)andYn(b) 38ActaPhys.一Chim.Sin.2011Vo1.27 离域效应较强,故前线轨道能级及能级差较低.因 此,由表4可见,酰胺链垂直于二氢卟吩环的Y.和 Y的最大吸收波长为681和678nm,比Y的628nm 大53和50nlF1.三个构象的最大吸收波长平均值为 662nm,略大于e6最大吸收波长644nm.而构象对 400m的吸收峰影响则不大. 3结论 用密度泛函方法对二氢卟吩e6及设计的6个 e6赖氨酸酰胺的稳定性,前线分子轨道能和电子吸 收光谱等性质研究发现,赖氨酸的s.NH:比.NH:与 e6形成的酰胺更稳定,e6乙羧基形成的酰胺分子Y 最稳定.e6赖氨酸酰胺的水溶性较e6有所改善,前 线轨道与e6相似,基本集中于二氢I:1吩环,但二氢 卟吩环平面性略为下降,使前线轨道能隙略为升高, 最大电子吸收波长605—627nm相对于e6略为蓝 移,但仍处于600—900nm光动力治疗窗口.构象对 最大吸收波长影响较大,e6乙羧基形成的酰胺的较 稳定构象,Y.和Y:最大吸收波长平均值为662 nm,比e6略为红移:酰胺链垂直于二氢卟吩环的Y 和Yn的最大吸收波长比酰胺链与二氢卟吩环近似 平面的Y约长53,50rim. Referenees 1Bisland,S.K.;Singh,D.;Gafi6py,JBioconjugateChem.,1999, 10:982 2Uzdenskya,A.B.;Dergachevaa,O.Y;Zhavoronkova,A.A; Reshemikovb,A.V;Ponomarevc,G.VSciences,2004,74: 2185 3Chen,J.B.;Ou,M.R.;Li,Z.Q;Xu,x.PStraitPharmaceutical Journal,2003,15(6):8[陈杰波,欧敏锐,李忠琴,许小平.海峡 药学,2003,15(6):8】 4Spikes,J.D.Photochem.Photobio1.B:Bio1.,1990,6:259 5Zhang,L.;Xu,D.YChin.Org.Chem.,1999,l9(4):424 [章玲,许德余.有机化学,1999,19(4):424】 6Pellerito,C.;DAgati,P;Fiore,T.;Mansueto,C.;Mansueto,V; Stocco,G.;Nagy,L.;Pellerito,LInorg.Biochem.,2005,99(6): 1294 7Vermathen,M.;Vermathen,P;Simonis,U.;Bigler,RLangmuir, 2008,24(21):12521 8Mennenga,A.;Gartner,,;Lubitz,W.;Gt~mer,HChem Chem.尸,2006,8:5444 9Fang,;Xu,DChin.LaserMed.勋曙.,1999,8(4):235 防瑛,许德余中国激光医学杂志,1999,8):235] 10O1shevskaya,VA.;Savchenko,AN.;Zaitsev,AV;Kononova, EG.;Pe~ovskii,PV;Ramonova,A.A.;Tatarskiy,VV; Uvarov,O.V;Moisenovich,M.M.;Kalinin,VN.;Shtil,A.A. OrganometChem.,2009,694:1632 IIHargus,J.A;Fronczek,FR.;Vicente,M.G.H.;Smith,K.M. PhotochemPhotobio1.,2007,83:1006 12Harubumi,K.;Kinya,F;Masami,O.LungCance~2003,42(1): 1O3 13Nelson,JS.;Robe~,G.;Bems,M.W.CancerRes.,1987,47: 4681 14Walsh,RJ.;Gordon,K.C.;Office~D.LMoLStruct.,2006, 759:i7 15Koseki,J.;Maezono,R.;Tachikawa,M.Chem.Phys.,2008, 129:085103/1 16Kuzmitsky,VA.;Volkovich,D.I. _,App1.Spectrosc.,2008,75 (1):27 17Kuzmitskii,VA.;Knyukshto,VN.;Gael,VI.;Zenkevich,E. I.;Sagun,E.1_;Pukhlikova,N.A.;Lebedeva,VS.;Mironov,A. F』AppI.Spectrosc.,2003,70(1):43 18Liu,X.Q.:Xu,Z.G.ComputApp1.Chem.,2009,26(3):339 [刘学琴,徐志广计算机与应用化学,2009,26(3):339】 19Pajunen,A.;Stapelbroek—Mollmarm,M.E,;Hyrminen,PH.Aota Cryste1996,52:743 20Lee,C.;Yang,W;PaR.GPhys.Rev.B,1988,37:785 21Reed,A.E;Weinhold,F.Chem.Phys.,1983,78:4066 22Vosko,S.H.;Wilk,L.;Nusair,M.Chem.P.,1980,58:1200 23Frisch,MJ.;Trucks,G.W;Schlegel,H.B.;eta1.Gaussian03, RevisionB.03.Pi~sburgh,PA:GaussianInc.,2003 24Ghose,A.K.;Pfitchea,A.;Cfippen,G.M.Comput.Chem., 1988,9(1):80 25Viswanadhan,VN.;Ghose,A.K.;RevankaqG.N.;Robins,R. K.ChemInf.Comput.Sci.,1989,29(3):163 26Ren,X.F.;Ren,A.M.;Wang,Q.;Feng,J.K.ActaPhys.一C. Sin.,2010,26(1):llO[任雪峰,任爱民,王钦,封继康.物理 化学,2010,26(1):ll0] 27Isakau,HA.;TruJ(11acheva,V;Pego~RPharm.Biomed. Ana1.,2007,45:20 28Oladkova,0.L;Parkhats,M.V;Oorbachova,A.N.;Terekhov, S.N.SpectrochimicaActaPartA,2010,76:388 29Isakau,H.A.;Parkhats,M.V;Knyukshto,VN.;Dzhagarov,B. MPe~ov,E.P.Pe~ov,P.T.J.Photochem.Photobio1.B:Bi0I.. 2008,92:165