北京理工大学博士学位论文synthesis of Ns_ salts. Here we provide a theoretical investigation on the decompositionand isomerization pathways for some complexes with general formula MN5, where M arethe alkali metal atoms Li, Na, K, and also Rb, with the intention to support experimentalstudies on the pentazole derivatives. The results indicate that the pyramidal MN5 are alllocal minima on their respective potential energy surfaces. However, the low barrier heightssuggest that the pyramidal structures are kinetically unstable and they will rapidly rearrangeto the most stable planar structures. The calculated barriers for the planar sturctures indicatethat the stability increases going from lighter to heavier metals (from Li to Rb) with themost stable of RbN5. In additional, rate constants using TST, CVT, and CVT/SCTcalculations for the decomposiiton reactions MNs-B-MN}+N2 weer performed. Thecalculated results imply that the small-curvature tunneling corrections are necessayr incomputing the rate constants of these decomposition reactions. Theerfore, we predicted thatthe theoretical CVT/SCT rate constants are probably more reliableKeywords: high-energy-density material (HEDM), nitrogen cluster, metal-nitrogen cluster,ab initio, density functional theory(DFT), decomposition and isomerizationbar rier, eraction rate constantiv研究成果声明本人郑重声明:所提交的学位论文是我本人在指导教师的指导下进行的研究工作获得 的研究成果。尽我所知,文中除特别标注和致谢的地方外,学位论文中不包含其他人已经发表或撰写过的研究成果,也不包含为获得北京理工大学或其它教育机构的学位或证书所使用过的材料。与我一同工作的合作者对此研究工作所做的任何贡献均已在学位论文中作了明确的说明并表示了谢意。特此声明。 签名:勺浅营日期:2D0于7.多关于学位论文使用权的说明 本人完全了解北京理工大学有关保管、使用学位论文的规定,其中包括:①学校有权保管、并向有关部门送交学位论文的原件与复印件;②学校可以采用影印、缩印或其它复制手段复制并保存学位论文;③学校可允许学位论文被查阅或借阅;④学校可以学术交流为目的,复制赠送和交换学位论文:⑤学校可以公布学位论文的全部或部分内容(保密学位论文在解密后遵守此规定)。签名:今1豁育导师签名日期:2-v378日期:今/尸2.又E北京理工大学博十学位论文第一章绪论.,氮原子簇研究的意义近年来,随着碳的同素异形体C6。的发现和由此产生的丰富而有魅力 的化学、物理、材料应用等方面的研究,激发起人们研究寻找新型的氮同素异形体的热情【1-40]。因为在元素周期表中,第一周期的C, N和。是最相近元素,并且现在已经得到了由纯碳和纯氧元素构成的分子,如C60,Coo... [41-45]和02, 03二等,这就使人们联想到是否可以有由纯氮构成的相似构型的稳定分子。众所周知,氮与CH基团是同瓣对应体,在一定程度上成键性质是彼此相似的,现在我们己知大量的碳氢化合物(CH)。是客观存在的,其C原子骨架构成了点(甲烷)、线(链烷、链烯、链炔)、面(苯、蔡、稠苯)、体(多面体C烷CH。和C笼原子簇)等多种结构,这启发我们具有相似构型的氮原子簇也可能稳定存在。但是N原子毕竟与CH原子团不同,它存在着孤对电子,孤对电子之间的相互排斥作用一方面会影响分子的几何结构,另一方面也会降低分子的稳定性,使一些分子事实上并不存在或是以亚稳态形式存在,因此N系列各种原子簇的结构研究引起了人们广泛的兴趣。除了结构意义之外,氮原子簇吸引科学家广泛关注的另一个原因是 它们作为一种新型环保型高能量密度材料(HEDM)的潜在应用价值【8,9]0我们已经知道N2的离解能(N三N键能)是225 kcal/mol,而标准的N=N双键是100 kcal/mol, N-N单键大约是40 kcal/mol。显然,N=N键能比三倍的N-N键能(3X40=120 kcal/mol)大得多,这表明对任何可能形成的包含氮氮双键和单键的N。原子簇(n>3),当它分解成N:分子时,将会释放出大量热量,使其成为很好的储能材料。Engelke[14」曾将N8立方烷的能量性质与另一种高能材料HMX(奥克托金)进行了比较(见下表):北京理工大学博士学位论文N8(0h)HMX密度g/cm'3爆速m m/ps爆压Gpa 爆热KJ/mol>2.2514.75 1372203191993︸所得的结果非常诱人,如果N8立方烷能够被合成,将会是一种性能优异的高能量密度材料。氮原子簇在分解时不仅能释放出巨大的能量,可以作为炸药,而且如果作为推进剂,通常可以通过氮原子簇分解为氮气而自身产生做功的工质,并不需要另外携带供氧物质,特别适用于火箭的能源,这对火箭推进剂来说是一个革命性的进步,而且分解生成的N2是无公害分子,因而是一个环保型的HEDM,具有潜在和广阔的高科技及军事应用前景。1.2氮原子簇研究现状当前国际军事装备的较量,本质上是武器和装备的高科技水平的竞 争。作为武器系统的能源材料,HEDM的性能有所改进,都会对武器系统的性能产生重大的甚至决定性的影响。从上世纪八十年代中期开始,高能亚稳态的氮原子簇作为一种新型的HEDM一直是世界上众多科学家关注和探索的目标之一。i.2.1氮原子簇的实验研究由于N。原子簇合成过程中的高吸热性,使其实验合成和使用一直面 临极大的挑战。事实上,到目前为止,已知的N。原子簇非常有限。从分子氮于1772年被分离出来后,N3一是在1890年由Curtius等【341第一个被合成出来的离子氮原子簇。此后的一百多年,尽管人们也进行了很多实验合成的研究工作【46-531,但大都没有取得成功。直至1998年末,Christe小组[[351才取得了重大突破。他们用NZF+AsF6和HN3在无水的HF溶液中,在一78℃下合成了NS+的AsFb-盐,而且证明了NS+是有史以来制造出北京理_L大学博士学位论文的最有破坏力的炸药,接着,他们「54]又报告合成出了另外一种心1.19稳定的N5+盐, N5+[SbF6r,这种NZN3 新的盐在70℃分解,并且抗冲击性较强。同时Christe等人也通过14N115N NMR和振动光谱(IR,N4+Raman)研究确定了NS+离子(16)是具有C2V对称性的线型结构。N5+盐的成功合成引起了人们极大的关注和很高的评价,这不仅是由于它是一种新的化合物,而且更重要的是为合成其他更多的聚氮化合物打开了通道。2001年初,Olah等人[55]试图在一78℃由氨基重氮盐制取N42+离子,然而遗憾的是,他们没有得到N4的二价阳离子。2002年7月,Vij等人[[56]通过断裂苯基五吟衍生-sN物的C-N键,第一次实验报道了82八长期寻求的cyclo-N5一的存在。此﹀外,实验上还用质谱方法检测到干太一些微量氮原子簇。如国外以激丫光溅射方法,采用飞行质谱仪检|测出N4-N30等偶数氮原子簇的会峰。用紫外检测到N2n+1(n=1-4)的光谱[33,39-42]。最近F.CacaceN6等人[[36]采用单电子还原气态N4+图1.1实验合成和检测到的氮的中性化一再离子化质谱,实验制原子簇类型(图中键长单位为埃)备和检测到了气态亚稳态的N4分子有超过1微秒的寿命。这些实验上的成果都极大地鼓舞了科学家的信心,为合成和检测其它更多的聚氮化合物北京理_r大学博士学位论文展开了光明的前景。2.2氮原子簇的理论研究如上所述,由于实验合成工作的巨大困难,因此大部分对氮原子簇 的研究工作局限在理论预测上。近年来,随着计算机技术的不断进步和计算方法的发展,有大量的理论研究被报道。Pyykko[26]和M atsunaga[33]早在1991年就预测了NS+的结构,而且预见如能实现氮原子簇的合成,其性能将是异常优越。Engelke等[12-15]曾对氮原子簇进行了理论研究,研究表明有可能合成出亚稳态的氮原子簇。最近,B artlett等「57」系统计算了N:到NI。原子簇及其阴、阳离子原子簇的结构、能量和光谱性质后指出聚氮化合物的结构是迷人的并具有潜在的爆炸性,他预测大量的聚氮化合物正等着人们去发现,关键的问题是合成的方法和分解的势垒。近来,美国LLNL的Mailhiot等人提出T“聚合氮" [Polymeric Nitrogen]这个全新的概念,其设想己进入初步实验阶段。他们认为在超高压,高温下氮原子能以共价键三维地连结起来,而成为网状的新物质。将氮原子均以共价键连结起来,则同一体积中所贮存的能量大得多,同一体积中氮原子的密度也大得多,因此,估计能量密度比起通常的含能材料要高达三倍之多。计算表明在常温常压下这样的聚合氮是亚稳态的,估算当这种聚合氮起爆时,将发出巨大的能量。这是关于氮的高能量密度材料非常有希望的新进展。在寻找新型氮原子簇的时候,常用的方法是和其它元素组成的已知 化合物对比。例如,磷的最稳定形态为Pa,因此相应结构的正四面体N4分子((1) (图1.2)也被期望是可以存在的,并且N;的能量比2N:高约(1)(2)(3)北京理工大学博士学位论文1.3.1热力学稳定性的研究按照原子簇结构与化学键之间的依据关系,并根据氮原子簇的结构 规则,首先对可能的氮原子簇进行结构设计,再利用当前的计算条件,尽可能的在高精度计算下,对所研究的氮原子簇所有可能的异构体进行全部优化和频率计算,并进行零点能校正,以期得到全势能面上的各种稳定异构体。1.3.2动力学稳定性的研究如前所述,研究的氮原子簇类型是否是有价值的HEDM侯选者,仅 仅有热力学稳定性是不够的,在动力学上也必须有足够的稳定性以保证其在使用和运输过程的安全性,也就是说这些异构体要有足够高的异构化和分解反应势垒。因此,在热力学研究的基础上,进一步对异构体所有可能的异构化和分解反应途径进行研究,找出反应的过渡态并利用内察反应坐标理论进行复核确认后,得到反应的中间产物、过渡态等表征指标及反应势垒高度的计算值,进而讨论其稳定性,从而为氮原子簇实验合成和检测提供更好的理论依据。本文在研究对象上不仅计算了一些纯氮原子簇,还计算了一系列金 属氮原子簇。从研究内容看,不单进行通常的几何结构优化和能量计算,而且还进行了分解、合成途径、位垒高度和反应速率常数等动力学数据的计算。从研究方法来说,从分子结构设计到分子反应设计均采用高精度量子化学方法进行理论计算和模拟,并结合谱学表征计算,形成了一套比较完整的计算体系。北京理工大学博士学位论文今考文献1. Trinquier G, Malrieu J, Daudey J, Chem. Phys. Left., 1981, 80: 5522. Franc] M M, Chesick J P, J. Phys. Chem., 1990, 94: 5263. Lee T J, Rice J E, J. Chem. Phys., 1991, 94: 12154. Glukhovtsev M N, Schleyer P v R, Int. J. of Quantum Chem., 1993, 46: 1195. Korkin A A, Balkova A, BartlettR J, et al.,J. Phys. Chem., 1996, 100: 57026. Lauderdale W J, Stanton J F, BartlettR J, J. Phys. Chem,1992, 96: 11737. Gimarc B M, Zhao M, Inorg. Chem., 1996,35:32898. Mathew L, Leiginger M. L., Schaefer H F, J. Phys. Chem. A, 1997,101:44609. Glukhovtsev M N, Jiao H, Schleyer P v R, Inorg. Chem., 1996,35: 712410. Saxe P, Schaefer H F, J. Am. Chem. Sec., 1983, 105: 176111. Glukhovtsev M N, Schleyer P v R, Chem. Phys. Let., 1992, 198: 54712. Engelke凡Stine J R, J. Phys. Chem., 1990, 94: 568913. Engelke R, J. Phys. Chem., 1992, 96: 1078914. Engelke R, J. Am. Chem. Soc., 1993, 115: 296115. Engelke凡J. Org. Chem,1992,57:484116. Li Q S, Hu X G, Xu W G, Chem. Phys. Lett., 1998, 287: 9417. Leininger M L, Sherrill C D, Schaefer H F, J. Phys. Chem,1995,99:232418. Tian A M, Ding F J, Zhang L F, et al., J. Phys. Chem. A, 1997,101: 194619. Gagliardi L, Evangelisti S, WidmarkP O, et al., Theor. Chem. Ace,1997,97:13620. Gagliardi L, Evangelisti S, RoosB O, J. Mol. Struc. (Theochem), 1998, 428, 121. Bliznyuk A A, Shen M, Schaefer H F, Chem. Phys. Let,1992, 198, 24922. Chen C, Lu L H, Yang Y W, J. Mol. Struct., 1992, 253, 123. Michels H H, Montgomery J A Jr., Christe K O, J. Phys. Chem,1995,99:18724. Gagliardi L, Evangelisti S, Bernhardsson A, Lindh R, Roos B O, Int. J. Quant. Chem.,2000,77:311 25. Nguyen M T, Ha T K, Chem. P冲s. Let., 2000, 317: 13526. Pyykk6 P, Runeberg N, J. Mol. Struct. (Theochem), 1991, 234: 279北京理工大学博士学位论文27. Li Q, Hu X, Xu W, Chem. Phys. Lett., 1998, 287: 949928. Benjamin M G, Zhao M, Inorg. Chem., 1996,35: 328929. Li J, Liu C W, Lu J X, J. Mol. Struct. (Theochem), 1993, 99: 22330. Mailhiot C, et al, Phys. Rev. B, 1992, 46: 1441931. Horenzana H E, Mailhiot C, et al., UCRL-ID-123644, DE96010714, 199632. Gonzalez C, Schlegel H B, J. Phys. Chem., 1990, 94: 552333. Matsunaga T, Kayaku Gukkuishi, 1998, 4: 20134. Curtius T, Ber. Dtsch Chem. Ges., 1890, 23: 302335. Christe K O, Wilson W W, Sheehy J A, et al., Angew. Chem., 1999, 38: 200436. Cacace F, Petris G, Troiani A, Science, 2002, 295: 480内︐j矛丹0﹂0Engelke R, Blais N C, Sander R K, J. Phys. Chem. A, 1999,103: 5611GonzalezC, Schlegel H B, J. Chem. Phys., 1989, 90: 2154Stamatovic A, Mark T D, J. Chem. Phys,1988,88:428939. Scheier P,40. Jonkman H T, Michl J, J. Am. Chem. Soc., 1981, 103: 73341. Kroto H W, Heath J R, et al, Nature, 1985, 318: 142. Curl R F, Smalley R E, Science, 1988, 242: 101743. Balasubramanian K, Chem. Phys. Lett., 1993, 204: 48644. Smalley R E, The Science, March/April 1991(New York Academy of Science, NewYor k) pp.22-2845. Thompson W E, Jacox M E, J. Chem. Phys. 1990, 93: 385646. Nguyen M T, Ha, T K, Chem. Ber., 1996,129: 115747. Huisgen R, Ugi I, Angew. Chem., 1956, 68: 70548. Wallis J D, Dunita J D, J. Chem. Soc. Chem. Commun 1983,910: 35849. Casewit C, Wenninger J, Roberts J D, J. Am. Chem. Soc., 1981, 103: 642850. Workentin M S, Wagner B D, Ne颤F, Zgierski, M Z, Lusztky J, Siberand W, WaynerD D M, J. Phys. Chem., 1995, 99: 9451. Workentin M S, Wagner B D, Lusztky J, Wayner D D M, J. Am. Chem. Soc., 1995,117:119 52. Vogel A, Wright R E, Kenkley H, Angew. Chem. Int. Ed. Engl. 1980,19: 71753. Gorini J A C Farras J, Feliz M, Olivella S, Sole A, Vilarrasa J, J. Chem. Soc. Chem,北京理工大学博士学位论文Commun. 1986,95:954. Vij A, Wilson W W Vij V, Tham FS, Sheehy J A, Christe K O, J. Am. Chem. Soc.,2001,1 23:630855. Olah G A, Prakash G KS , RasulqJ. Am. Chem. Soc. 2001, 123: 330856. Vij A, Pavlovich JqWilson W W, Vij V, Christe K O, Angew. Chem. Int.Ed. 2002,41:3051 57. Bartlet R J, Chem Indust. 2000,4: 14058. Gagliardi L, Evangelisti S, Barone V, Roos B O, Chem. Phys Lett. 2000, 320: 51859. Nguyen M T, Ha T -K, Chem. Phys. Let. 2001, 335: 31160. ChungQSchmidt M W Gordon M S, J. Phys. Chem. A 2000,104: 564761. Li Q S,Wang L J, J. Phys. Chem., 2001, 105: 120362. Korkin A A, Balkova A, BartletR J, Boyd R J, Schleyer P v R, J. Phys. Chem., 1996,100:5702 63. Yarkony D R, J. Am. Chem. Soc., 1992, 114: 540664. Schmidt M W, Gordon M S, Boatz J A, Int. J. Quantum. Chem., 2000, 76: 43465. Wilson K J, Perera S A, BatrletR J, J. Phys. Chem. A 2001, 105: 769366. Straka M, Chem. Phys. Lett. 2002, 358: 53167. Gagliardi L, Pyykk6 P, J. Am. Chem. Soc. 2001, 123: 970068. Lein M, Frunzke J, Timoshkin A, Frenking G, Chem. Eur. J. 2001, 7: 415569. Hammerl A, Klapbtke T M, Inorg. Chem. 2002, 41: 90670. Kwon O, Almond P M, Mckee M L, J. Phys. Chem., 2002,106: 6864北京理工大学博十学位论文第二章理论原理和计算方法量子化学是应用量子力学的基本原理来讨论分子的结构、性能以及分 子间相互作用与化学反应的一门科学,在众多的学科领域中有着广泛的应用,其对相关体系的研究主要是通过求解体系的Schr6dinger方程来实现的。然而,对于实际上重要的体系,不引入一些近似,就不可能求解其Schrodinger方程,到目前为止,仅对简单的H2`. He和H:进行了精确的求解。依据采用的近似程度,可以把量子化学计算方法划分为两大类:从 头计算法和半经验性近似方法。从头计算法只限于引入一些基本近似,非相对论近似、Born-Oppenheimer近似[1-2]和单电子近似[3]等,而不引用任何原子的实验资料。采用这些近似之后,使分子的Schrbdinger方程化为可解的形式。然而,在实践中这类从头算方法仍存在较大的困难,为了将量子化学方法广泛地用于复杂化合物或者化学反应体系,不可避免地要借助于半经验量子化学计算方法。 本论文主要应用量子化学从头算和密度泛函理论(DFT)方法计算研究了几种聚氮原子簇的结构和稳定性方面的性质。下面将涉及的基本原理、计算方法和实施方案简述如下:2.1势能面势能面( potential energy surface)概念是由著名的Born-Oppenheimer近似【1]引入的,它是该构型的函数.在描述核运动的薛定愕方程中,势能面相当于势能项。反应体系核的运动可以形象地描述为在势能面上的运动。有了势能面的解析或数值网格,解含时或非含时的薛定is方程,就可以由此获得核运动及化学反应的重要信息。然而,利用从头算方法构造完整的势能面是一件非常困难的事,即使不考虑构造过程中的技术难度,仅从计算量上看就是惊人的,即使选用最简单的Hartree-Fock方法也是很困难的。13 北京理工大学博士学位论文pv, =2F,CwC;S,是原子轨道中;和中,的重叠积分:(2.6)S u,一执(l从(l》:(2.7)。,为相应的轨道能量。HFR方程也可用矩阵形式表示为 通常把F, h和G矩阵分别称为Fock矩阵,Hamilton矩阵和电子排FC=SCF F=h+G (2.8)斥矩阵。矩阵P=几,称为单电子密度矩阵,S=-SN,称为重叠矩阵。HFR方程在形式上是求解本征值问题,C相当于算符F的本征值,C相当于算符F属于本征值£的本征向量。但是,它与一般的本征值问题不同,因为算符F本身与轨道组合系数{C砂有关·于是,在求解HFR方程时只能用迭代的方法,即所谓的自洽场(Self-Consistent-Field, SCF)方法。迭代是否收敛的判据有两种,一种是本征值判据,一种是本征向量判据。在Gaussian程序中,本征值判据为10-1,本征向量判据为10-1.简单的分子轨道理论可用两种方式扩展到开壳层体系,一种是自旋限 制性的Hartree-Fock(RHF)理论,另一种是自旋非限制的Hartree-Fock(UHF)理论。在UHF理论中,两组分子轨道分别由两套展开系数来定义 平卜Icp,4uWp=YCwcN(2.9)保持a ,1 i两组轨道各自正交归一的条件,经过与闭壳层相类似的数学处理,可得Lr }Fa一:S4v f":一“全} F1P1一。DSHV }pi一“这时,Fock矩阵元为(2.10)(2.11)北京理工大学博士学位论文F,}'}=H"""十E L-. P",[(PI1 +P,kw,入︑口 ︑︑ 了了一PL }pk I Vaj(2.12)Fu一H"、艺艺P, [(pz、P,kvv ︑义︑︑仃尹尹一P.\1}xI V4(2.13)密度矩阵元为P0,=艺Cw C}ps=Y Cw Cn.(2.14)这就是开壳层分子的Hartree-Fock-Roothaan方程。2. 3从头算基组的选择【5-7]如上所述,分子轨道可以表示成原子轨道或任一组基函数的线性组 合。在从头算方法中,广泛使用的基函数有两种:Slater型基函数((STO)和Gauss型基函数(GTO) eSl ater型基函数下的原子轨道在描述电子云分布方面较其它基函数具有一定的优越性。然而在计算包含Slate:原子轨道的三中心和四中心积分时,由于用到1/r12的无穷级数展开式,从而使积分变得十分复杂。如果用若干个Gauss函数的线性组合去拟合Slater原子轨道,将会使难以处理的多中心积分大大简化。Slater函数反映分子中的电子运动远比Gauss函数好,但计算时远不及Gauss函数方便。为把这两个函数的优点结合起来,可以用Slater函数向Gauss函数集合展开的方法来选择基函数。在Slater函数向球Gauss函数集合展开的过程中,根据计算需要仅近似地选择k项,系数和参数用最小二乘法进行选择。这个基函数集合通常用STO-KG表示。较STO-KG更好的基组是“分裂价基”(Split Valence BasisSets)。分裂价基是通过每个价层Slater原子轨道用两个或多个新的具有不同指数的STO-KG来逼近,从而改善计算结果与实验结果的偏差。例如6-311G,其中每个内层轨道用一个STO-6G逼近而每个价层轨道用一个STO-3G和两个STO-1 G逼近。 对于元素周期表中第三周期以下的元素(主族及过渡金属元素),随着d轨道的出现,其化学性质及结构特点都将发生变化,此时含有极化函数北京理工大学博士学位论文的Gauss基对于计算结果处理变得更合理。此外,对于含有孤对电子、带电离子和弱相互作用体系可以在已有的基函数的基础上引入扩散函数,如6-311+G*基组。2.4后自治场(Post-SCF)方法[1113]自洽场( (SCF)方法中假定电子在原子核及其它电子形成的平均势场中独立地运动,考虑了电子间平均作用,但没考虑电子间瞬时相关,即平均场中独立运动的两个自旋反平行电子可能在某一瞬间在空间某点同时出现。由于电子间Coulomb排斥,这是不可能的。单组态自洽场(SCF)计算未考虑这种电子相关作用,导致计算误差。处理这一电子相关问题的方法被称为电子相关方法或后自洽场(Post-SCF)方法,其中包括微扰理论(MP)、组态相互作用方法(Cl)等。2.4.1 Meller-Plesset微扰理论[14171电子相关能与体系的总能量相比是个小量,通常只有双重激发组态占 重要地位,因此可以用多体微扰理论(Many-body Perturbation Theory,MBPT)计算电子相关能(1-3)0将Hamilton算符表示为H=Ho+V(2.15)式中n。为无微扰的Hamilton算符,伞为微扰量。薛定愕方程n`P = E'P可表示为位一n。)甲卜VI甲)将甲按H。的本征函数系。。展开HncF;=E;D(; (i=0, 1, 2,…,‘)(2.16)(2.17)(2.18)T一艺a,D(, (ao=1)北京理工大学博士学位论文一般EKE,·设正交归一化条件为(,,}电卜气,(0o I T)一,,则由(“·,“,式可得( E一E加、=(。,1V1 T) (2.19)用投影算符P0 = Imo)(D,0}把函数甲投影到。。子空间内,PoJT)一(J D(OLPTI if0 ))一a0ID(0)则(2.20)甲=叉 a,D(,二。。+ -aY;D(;1一P六._、_久_, .、_A= (P+一一-=VJT)=W+UVI`Y)E一Ho (2.21)用叠代法求解上式,得甲一j(GVI lD(0)(2.22)E一E0 +()q0IVIT)-一<D(o{ Vj-0卜>(GVr结果得到(2.23)甲一1D(o)+I甲‘,,)+一甲‘,,)+…E=Eo+SO)+E (Z)+E(')+二GVGV】其中,TO)一GVl(po),甲(,)=。。),…(2.24)E( ')=)(。。{,},。),。‘”一(,。一,。。‘”_(}QIVGVGVIc}o),...VGV(2.25)将E按微扰参数x的级数进一步展开,并将((15)式改为n=Ho +2V(2.19)和(2.24)式整理得:代入 (E。一E;}‘一()a,I}v+E。一EIT)(2.26)北京理丁大学博十学位论文`P=j(G.XVrl(UQ)(2.27)一(}o1 },V}.-0、}Go一,一>其中xV=xV一(E一Ea),行=一Eb一Ho _I一P(2.28)(2.29)无微扰算符宜。的最简单的选法是取Fock算符之和(这种Hamilton量的划分称Moller-Plesset划分)。基于这种划分的MP微扰方法,根据选取微扰项的不同,分为二级微扰(MP2)理论、三级微扰(MP3)理论和四级微扰(MP4)理论,其中含四级微扰项(MP4SDTQ)的理论包含了四电子跃迁的组态,有很高的计算精度。2.4.2组态相互作用(Configuration interaction, CI) [13, 18-20]组态相互作用( CI)是最早提出的计算电子相关能的方法之一。从一组Fock空间完备的单电子基函数{)x(K`数集合{`1'k)}出发,可造出一个完备的行列式函D(k=(N!)一,/2}jqk (X I )T12 (X 2) ... TkN (X N) , (2.30)任何多电子波函数都可以用它来展开。一般{ `Pw(x)}称为空间轨道, l`i'ki称为组态空间。在组态相互作用( C工)方法中,将多电子波函数近似展开为有限个行列式波函数的线性组合(CI展开)T一艺CaDa=。。+YEca)q;+}}c LrY C aijbD(ia b +}}L r Z -,c C aijbk+…Qa力ij(2.31).,6p i,j'I用变分法确定系数CS,即选取CS使体系能量取极小值,得到广义本征方程He=ScE(2.32)北京理r大学博士学位论文其中H=\D(SIHIV-`), S =恤lD(, ),。为系数矩阵,满足以下条件Co Sc;二艺C.,psasry-8m(2.33)若中:是正交归一集合,则上两式变为(2.34)=6因(2.35)组态相互作用( CI)方法中。s称为组态函数,简称组态。它是一种行列式函数。为提高计算效率,一般它满足一定的对称条件,如自旋匹配条件、对称匹配条件等。完全的CI计算能给出精确的能量上界,而且计算出的能量具有广延量的性质即“大小一致性”。然而,由于CI展开式收敛慢且考虑多电子激发时组态数增加很快,通常只能考虑有限的激发,如CISD表示考虑了单、双激发。这种截断的CI计算不具有大小一致性。Pople等人通过在CI方程中引入新的项从而使非完全CI计算具有大小一致性,新项以二次项出现,该方法就称为QCI(Quadratic Configuration Interaction)方法。QCISD方法除避免了CISD中的大小不一致性外,还包含了更高级别的电子相关能。2. 5密度泛函理论(Density Functional Theory, DFT) [21-331DFT理论的发展经过了一个漫长的过程,直到1964年才有一个确定 的模型。该理论的开始可以追溯到1927年前后Thomas和Fermi的工作。他们认识到可以从统计学的角度(statistical considerations)来近似处理原子中电子的分布。Thomas将这一假设陈述为:“对于每个h3体积中2倍速度(At the rate of two)的电子运动,电子均匀(uniformly)分布于六维空N之中,并且存在一个核电荷,由电子分布所决定。”Thomas-Fermi的电子密度公式即由上述假设推得。 我们将Thomas-Fermi建立的模型称为均匀电子气模型,假定三维势箱中电子分布是均匀的,该模型中动能表述为:T" H=C,- f P,Y(r*,CF =希(3n21一2.871(2.36)能量为:北京理工大学博士学位论文E_一C} IpX(rpr一Z四+ 1 ffp(n)P(rz}dr,dr,(2.37)2“I r,一r,由上式可以看出,在T-F模型下,体系的能量是其电子密度的泛函,这便是T-F理论蕴藏着密度泛函的思想。对于T-F电子气模型,若用变分法和自洽场手段求电子密度和体系的能量,则这种方法称为密度泛函理论(Density Functional Theory, DFT)< 1964年,Hohenberg和Kohn首先证明:对于电子数不变的基态体系,其Hamilton量是电子在空间密度分布P(r)唯一确定的泛函,反之亦然.与此同时,对于一个近似的(或试探性的)电子密度PW,如果有P(r) ?“和介(rk1r=N(N为电子总数),则变分原理成立:Eo、E}] (2.38)此称为Hohenberg-Kohn定理或Hohenberg-Kohn变分定理。其中,最具吸引力的是对上式右端进行变分处理。变分时,只涉及到电子的空间坐标三个变量,待优化参数仅仅是包含了三个变量的电子密度(相应HF方法中一般有n2个待优化的参数)。为此工作,Kohn获得了1998年Nobel化学奖。Hohe nberg和Kohn未给出泛函E卜01的具体形式,1965年,Kohn和Sham将能量表达式写成:E=ET +E0 +E'+Ex` (2.39)其中,E7是电子运动的动能项,E“包含原子核与电子吸引以及原子核间的排斥势能,E“为电子排斥项(或理解为电子密度本身的库仑相互作用),Exc为交换相关项并包含其余的电子一电子相互作用。ET十Ev +E“对应一定电荷分布下的经典能量,而E xc包括了量子力学波函数要求的反对称化的交换能量和单个电子运动的动态相关两部分能量。因此,势能泛函Ex`'b( ,刀的具体形式是解决问题的关键。Hohenberg和Kohn证明Ex`7完全由电子密度决定。在实际应用中E xc被近似为只和自旋密度以及梯度有关的积分:EX` (P)= ,f.f(P a(r),Pp(FjVPa(FjVPg(r)}3r(2.40)其中P., P。分别是。、p自旋密度。P表示总电子密度。一般将EX`7写成:Ex "(p)=EX (p)+E`,扣)(2. 41)北京理工人学博十学位论文上式右端第一项称为交换泛函(Exchange Functional),第二项称为相关泛函(Correlation Functional),两者又可进一步写成局域(local)泛函部分(只与电子密度有关),非局域((non-local)泛函部分(只与密度梯度有关)和杂化(hybrid)交换和相关函数。定域交换泛函一般定义为: 二/:、月_、 ELxDA一亏(云{J p"d'r(2.42)这种形式被用做计算归一化电子气的交换能,但是,不能很好地描述分子体系。Becke于1988年在LDA交换泛函的基础上发展了梯度校正泛函:r''R-ke88一“LDA一‘J 1 + bysinh-' x‘P SXZ(2.43)其中二=P场PI。参数爆为了拟和惰性气体原子的交换能,Beck。取它的值为0.0042Hartrees。如上式所示,Becke的泛函是定域LDA交换泛函的校正,它成功地弥补了LDA的许多不足。同样,相关泛函也有定域相关泛函和梯度校正泛函。如Per dew和Wang1991年提出的相关泛函的定域部分:E c一Jps,(,(p(F)),}d)'+' , r一L4c7p],毛二Pa +Pa (2.44)纯的DFT方法定义为一种交换泛函和一种相关泛函的组合。例如「3飞 3_PY。一P著名的BLYP泛函,组合了Becke的梯度校正交换泛函和Lee, Yang, Parr的梯度校正相关泛函。Hartree-Fork理论的公式中确实包含一个交换项。Becke提出了将x`Hartree-Fock交换项和DFT的交换泛函棍合,和DFT相关一起定义E 为:45)E Eh xyc=chrtd = c二E众+cxr EexF + e。D二E器Fr ExDhc Y( 2.其中c为常数。例如,Becke的3参数泛函B3LYP可以通过下式定义E}R3f.YP=:X+cELDA + c。a ( XEHF一:LDA)+CxAE88+:赢,Ec VWN 3 +c+ c。r (/\ cELYP一:V WN3) (2.46)=072,cc=0.81。其中,c, = 0.20, cX现在,DFT己成为电子结构理论中解决许多难题的强有力的而且有效北京理工大学博士学位论文的工具,并成功地扩展到激发态以及与时间有关的基态性质研究。现代DFT如B3LYP方法计算精度大致与MP2方法相似,与MP2相比DFT计算的基组要求也小得多。更重要的是,在给定的基组下,DFT方法的计算资源要求与HF相似而比一般的相关方法小得多。2. 6速率常数理论「34-37]传统过渡态理论是建立在以下四个假设的基础上:①核运动绝热假定 即核运动和电子运动分开处理,把核看成是在一定电子态的势能面上运动;②核运动可以做经典处理,忽略隧道效应和相对论效应;③反应体系服从Boltzmann统计平衡分布规律;④全部反应轨线都一次穿过分割面,不再返回反应物区。据上述基本思想可推得反应速率表达式为:k(T)=Q Q'(T)脚QR(T)exp[(_flAE(sx ))(2.47)式中a表示反应物到产物的反应通道数,也称对称数;,li= (ka7) ', k。是Boltzmann常数;Qx表示过渡态的单位体积的配分函数;Q”表示反应物单位体积的配分函数;AE(S")经典势垒;Sm是过渡态点(S=O)。可以看出,传统过渡态理论是把分割面选在S=O的过渡态处用统计力学方法计算反应速率常数。但是这种基于轨线不返回假设得出的结果往往高于真实的反应速率常数。正则变分过渡态理论( CVT)是在传统过渡态的基础上发展起来的,它通过变分的方法选择不同的分隔面,通过改变划分过渡态分隔面的方法来寻求最小的速率常数。以k(s, 7)对反应坐标s变分求km,.(s, 7)等价于以4 G(s, 7)对反应坐标s变分求△Gm.(s, 7)ak(ass:,T)二、=“即aAGa(s一s, T)5一、一。(2.48)北京理工大学博士学位论文就是将分割面s选在△G(s, 7)最大值处s(通常不再为0)。因此,正则变分过渡态理论速率常数表达式为:k(s, T,一6,hK}exAG(s,T)R丁= - 。_。「 n exq-PV(so)+ln KQ'Q(sR,(T)一_,8fh QQF(Q s(oT,T) )exp[ /iv‘一”(2.49)显然,变分过渡态理论允许轨线有返回。2.7计算方案的实施和结果分析2.7.1分子平衡几何优化和过渡态(TS)搜索电子结构计算的最一般用处或许是预言分子几何。在 Born-Oppenheimer近似下,作为核坐标函数的分子总能量定义了位能面。原子核在这种位能面上运动定义了可能的化学反应、分子振动等。位能面上最有意义的点是一些驻点:鞍点定义了过渡态(TS),而极小点给出了分子平衡几何。G98程序中几何优化从初始构型开始,计算能量和梯度,然后决定下 一步的方向和步长,其方向总是能量下降最快的方向进行。大多数的优化也计算能量的二阶导数来修正矩阵,从而表明在该点的曲度。有一些系统的优化很难进行,采用Gaussian98默认的方法得不到结果,产生的原因往往是计算出的力矩阵与实际的相差太远。当默认方法得不到结果时,就要采取一些其他的方法。Gaussian提供很多的选择,比如选用关键词OPT=ReadFC, OPT=CaleFC等。在进行 TS搜索和构型优化时先要选定反应坐标(Reactions北京理下大学博十学位论文Coordinates, RC),固定RC值并进行优化得到该RC处的最小能量Emin逐渐增大RC使之从反应物变到产物构成了反应途径,在RC变化过程中Em;。出现一个最大值,过渡态位于该点附近,进一步进行过渡态优化得到所要寻找的TS. G98程序为实际计算提供了两类可选择的方案:1)二次同步搜索(Quadratic Synchronous Transit, QST)方法,包括QST2和QST3两个可选项,前者需要输入反应物和产物分子构型,TS的初始构型由程序根据产物和反应物构型自动进行估算;后者需要输入反应物、产物和猜测过渡态的初始构型。这类方法是用Quasi -Newton-Rhapson(利用一阶微分对二阶微分进行估算)方法或者本征值跟踪法(Eigenvector Following)对TS进行几何优化,简称STQN (Synchronous Transit-Guided Quasi-Newton); 2)过渡态优化,OPT=TS。这一方法一般与关键词Calcfc(计算力常数)连用,即用Newton-Rhapson方法对力常数矩阵进行数值计算(默认的方法为Quasi-Newton-Rhapson),这样虽然加大了计算量,同时也增加了TS优化成功的机会.为确保用以上方法得到的是所讨论反应的过渡态,须进一步通过内 察反应坐标(Intrinsic Reaction Coordinates IRC)计算来检验。IRC计算从过渡态出发沿着虚频所对应的振动模式的两个相反方向逐点进行几何构型优化,得到的最终几何构型分别对应于反应物和产物。2.7.2振动频率的计算和分析一个由N个原子组成的分子在其平衡构型附近的总能量可写为:E一T·v一2 1Y 3Nq?·evV`q9·+ _ 2 1E3N ( (aaq1,avq; ),q;q;(2.50)对于经典的动力学运动方程,有:9;=一艺f;q;,.1=1,2,3,...,3N(2.51)式中f;i称为二次力常数(quadratic force constants),它是在平衡构型时热能相对于加权质量笛卡尔位移qj的二阶导数,北京理工大学博士学位论文2了les︐.︑ 即f 与 a=v一 一 aq丙;(2.52)对方程( (2.48)可通过标准的方法解得3N个正则模式的振动频率,其中振动模式(对于线性分子是五个)的频率为零,对应于平动和转动的自由度。使用Fr eq关键词,G98程序在自洽场解HF方程的基础上可进行分子的红外/拉曼振动频率、振动强度计算和振动模式分析,还通过对振动、转动、平移和电子配分函数的计算预测化合物零点能和各种热力学函数,如系统的嫡和焙。2.7.3布居数分析【5]电子密度分布图能很好地表达分子中电荷的空间分布,但是作图方 法比较麻烦,并且没有把分子与其组成原子联系起来。为了表明电荷在各组成原子之间的分布情况,Mulliken在1955年提出了“布居数分析”的方法。按照分子轨道理论,分子的电荷密度是由各分子轨道上的电子贡献 到一个分子轨道上,或者没有电子或者有整数个电子。但是,我们也可以把分子电荷看成是由组成分子的各原子的原子轨道上的电子组成的。从原子轨道角度看,每个轨道上的电子数不一定为整数,电子在各原子轨道上有一个分布或者“布居”。分析这种分布,对于了解分子中原子的成键情况有帮助,与化学家传统的认为分子是由原子组成的观点也比较一致。设分子有N个原子( 用A, B, C,…标记),n个电子,分子轨道为中,=Y-1c,OXAA卜 (2.53)我们有 其中x AN表示A原子的,轨道。(2.54) 咖f=艺艺艺艺c沁CBaiXApXBXA Bp两边积分,并乘以ni可得北京理工大学博十学位论文n; =n,艺艺!c"i{,-F 2n;艺艺艺艺CA,uCBxiSAU.aaA>日px(2.55)利用这个式子可以做出电荷分布的分析如一F:(1)分子轨道中j中的电子分布在原子轨道xA,上的电荷为n(i,AR)=ni ICA,U12 (2)对i求和,得到原子轨道xA,上的电荷为(2.56)(2.57)n(Ag)一艺n(i,Aix)一艺n; I CAW},一tP".(3)对A原子的所有原子轨道求和,得到A原子上的电荷为(2.58) n(A)一艺n(AP) =ZPG.Nt.(4)分子轨道i中,属于A原子的轨道和B原子的轨道的重叠区电荷为n (i,AU,BX)=2n;c沁CBMSAa.Bk (2.54;(5)对分子轨道i求和,得到的A原子的u轨道和B原子的x轨道的总重叠电荷为艺2nic沁CBxisAN,BX=2P-,声,、=2pxpn(AR, B?.)=艺n(i, AP, Bk.)二通常把(Pxw)称为按原子轨道的布居矩阵。(2.60)(6)对A和B原子的所有轨道求和,得到的A原子的轨道和B原子的轨道的总重叠电荷为n(A, B) = 2艺艺氏S}, = 2艺艺Pz,t(7)最后,对分子中所有原子求和,自然得到(2.61)-Y-YPL-s-=E D-。二艺n(A)十艺艺n(A,B)一艺艺瑞十2艺艺凡SXA-A卜人>B轰( 2.62)其中求和遍及所有原子轨道。1Mul liken把重叠区电荷平均分配给有关原子轨道。对布居矩阵的1行求和,得到的A原子p轨道上的总电荷为n 、一肠+1t2艺凡、十112艺巧一肠+艺心,(2.63)北京理工大学博士学位论文因为丙一几,,引入112是由于只有一半重叠电荷分给原子A的。轨道。对属于原子A的所有原子轨道求和,得到的A原子上的总电荷为n*一InA ,二Y-A-YQuv(2.64)ppV于是,A原子上的净电荷为qA=ZA-nA-本论文的全部计算工作均是应用Gaus sian98W量子化学软件包,在PC机上完成的。论文中的几何构型图采用Chem3D和ChemDraw程序绘制完成,应用POLYRATE 8.2程序计算了不同温度下的分解反应速率常数。北京理工大学博士学位论文今考文献1. Born M, Oppenheimer J R, Quantum theory of molecules. Ann Physik, 1927,84:4572. Born M, Huang K, Dynamical theory of cyrstal lattices. Oxford University Press, NewYor k, 19543. Hartree D, Calculations of atomic structure. Wiley. 19574.唐敖庆,杨忠志,李前树,量子化学.北京:科学出版社.19825徐光宪,黎乐民,王德民.量子化学基本原理和计算方法(上、中册).北京:科学出版社.1994 6.丁益宏,新型原子簇一一两类重要星际分子势能面的量子化学研究:[博士学位论 文].长春:吉林大学理论化学研究所,19987.肖鹤鸣,硝基化合物的分子轨道理论,北京:国防工业出版社19938. Hehre W J, Radom L, Schleyer P v P, People J A. Ab initio molecular orbital theory.New York: Wiley&Sons, Inc., 19869. Dykstra C E. Ab initio calculation of the structures and properties of molecules.NewYor k: Elsevier Science Publishers, 198810. Szabo A, Ostlund N S. Modern quantum chemistyr. New York: McGraw-Hill, 198911. Lowdin P O, Correlation problem in many-electron quantum mechanics. I. Review ofdi ferent approaches and discussion of some current idesa. Adv. Chem. Phys., 1959,2: 20732212. Wilson S, Electron correlation in molecules. Oxfore: Clarendon Press, 198413. Bauschlicher C W Jr, LanghoffS R, Adv. Chem. Phys., 1990, 77: 10314. Moller C, Plesset M S, Phys. Rev., 1934,46: 61815. Head-Gordon M, Pople J A, Frisch M J, Chem. Phys. Lett., 1988, 153: 50316. Frisch M J, Head-Gordon M, Pople J A, Chem. Phys. Let., 1990,166: 27517. Frisch M J, Head-Gordon M, Pople J A, Chem. Phys. Let., 1990, 166:28118. Siegbahn P E M, J. Chem. Phys., 1980, 72: 164719. Knowles P J, Handy N C, Chem. Phys. Let., 1984, Ill: 31520. Pople J A, Binkley J S, Seeger R, Int. J. Quant. Chem. Symp,1976,10:130北京理工大学博十学位论文21. Hohenberg P, Kohn W, Phys. Rev. B, 1964,136: 86422. Kohn砚Sham LJ , Phys. Rev. A, 1965, 140: 113323. Parr RqYang W, Density functional theory of atoms and moleculesOxf ord University Press, 1989New York:k:24. Labanowski J K, Andzelm J, Density functional methods in chemistry.New YorS pringer-Verlag, 199125. Ziegler T, Chem. Rev., 1991, 91: 65126. Johnson BqGill PM砚Pople J A, J. Chem. Phys.,1993, 98: 561227. Siegbahn P E M, Adv. Chem. Phys,1996,93:33328. Vosko S H, Wilk L, Nusair M, Can. J. Phys., 1980, 58: 120029. Becke A D, Phys. Rev. A, 1988, 38: 309830. Perdew JP, Phys. Rev. B, 1986, 33: 882231. Lee C, Yang砚Parr RqPhys. Rev. B, 1988,37: 78532. Becke AD, J. Chem. Phys., 1993, 98: 137233. Becke A D, J. Chem. Phys., 1993, 98: 564834. Garrett B C, Truhlar DqJ. Phys. Chem,1979, 83, 105235. Truhlar l)qGarrettB C, Ace. Chem. Res,1980, 13, 440Dynamics,36. Truhlar, DG;Isaacson, A. D., Garrett, B. C. Theory of Chemical ReactionBaer , M. Eds., Boca raton, CRC Press, FL. 198537. Garrett B C, Truhlar DqJ. Am. Chem. Soc., 1980, 1胆,2559北京理工大学博十学位论文第三章3. 1引言N,原子簇结构和动力学稳定性研究最近一段时间以来,聚氮和富氮化合物的研究吸引了众多实验和理 论化学家的广泛关注,一方面这不仅是由于氮原子簇本身的结构意义,而且更重要的是由于它们作为环境友好高能量密度材料(HEDM)的潜在应用价值。虽然N和CH是等电子体,并且己有很多稳定的(CH)。化合物,比如苯和聚链烷,然而,大多数相应的氮原子簇化合物仍没有制备出来。在1772年首次分离出N2分子、1890年合成出叠氮阴离子N厂以后,直到1999年,Cbriste和合作者「1]才成功地合成出第三种稳定的聚氮化合物N5+盐(NS+AsF6-).最近,Cbriste等人[2]也报告了另外一种稳定的N5+盐,N犷SbF6,这种新的盐在70℃分解,并且抗冲击性较强。另外,也有光谱探测到N3, N3十,N4+和N6一等自由基和离子的研究报道。有没有更多种类的聚氮化合物的存在?在过去几年里,理论上已经 预测了很多稳定的聚氮化合物结构[3-24],诸如N4[5-8], N6[5-6, 10-121和NS[5-6, 13-16]. Bartlett等人[17]对NZ到N,。原子簇以及它们的正、负离子的结构、能量和光谱性质进行了系统研究,指出大量的聚氮化合物正等待人们去发现。Gagliardi等人「18]在他们最近的一篇论中报告了N5-负离子有足够的稳定性且进一步研究了从N31 Ns一和NS十合成N:和Nio的可能性。在氮原子簇的合成研究过程中,科学家们对偶数氮原子簇寄予了更 多的希望。然而,有趣的是,到目前为止,只有N3一和N5十被成功地合成。这启示我们奇数氮原子簇的合成也许是更有希望的。显然,N7和它的离子原子簇应该是下一个很有希望的合成目标。然而,到目前为止,只有Q.S.Li等[21]和Wang等人[[22]对N7原子簇进行了一些理论研究,Li等[211的计算结果表明开链Cs和C:结构((2,见图3.1)是N7势能面((PES)上的最小点,而Wang等人[22]采用Gaussian-3方法确定了12个低自旋异构体为势能面上的局域极小点,其中包括一个新的最稳定的异构体(1I32 北京理工大学博士学位论文见图3.1)a鉴于以上的计算结果,我们采用Ab initio和DFT方法对以前研究过的三种异构体((1, 2和5)和两种新的异构体((3和4)进行了研究。为了进一步评价这些异构体的动力学稳定性,我们也进行了分解和异构化反应机理的研究。这是首次从动力学角度对N:原子簇的稳定性进行探讨。我们希望目前的研究工作能对将来寻找稳定的氮原子簇提供一些理论的数据。3. 2计算方法对自由基( 奇数电子)体系而言,可以选用自旋限制性开壳HF(ROHF)或自旋非限制性HF(UHF)方法。Byrd等人[251最近对33个小自由基分子的结构和频率进行了系统的计算研究,结果表明ROHF方法比UHF方法能更好地描述自由基的性质。另外他们的研究结果也显示B3LYP方法的计算结果可与CCSD相比。所以从机时和计算准确性来看,B3LYP是对自由基分子进行常规计算的一个很好方法。在这部分工作中,我们分别采用限制性、非限制性的二级微扰理论(ROMP2和UMP2)和密度泛函理论B3LYP方法对分子的几何和振动频率进行计算。在有关波函数的计算中,核轨道被冻结。整个计算采用6-31G*和6-311G*基组。频率的计算可以判定所得结构是为全实频的极小点还是具有一个虚频的过渡态。在B3LYP和UMP2水平下进行了内察反应坐标(IRC)[26, 27]的计算,验证所得的过渡态确实是连接反应物和生成物的最小能量路径。为使势能面上的能量更精确,我们还用CCSD(T)/6-31 I G*方法在B3 LYP/6-311 G*优化出的几何的基础上进行了单点能的校正。本文中所有的计算均利用Gaussian98 [28]量子化学软件包完成。除特别说明外,键长单位为埃〔 A),键角为度((degree),全部能量为Hartrees,相对和零点振动能为kcal/molo北京理_r.大学博十学位论文3. 3结果与讨论3. 3. 1 N:原子簇稳定异构体图3. 1显示了研究的五种异构体1一5的优化几何,图中的“T”代表二面角。它们的点群、基态的对称性、全部能量和零点能都列在表3.1中。表3.2列出了五种异构体的相对能量(包括零点能修正)。两种新的开链线性异构体((3和4)的振动频率和红外强度列于表3.3中。109.24B3LYP/6-31 G'UMP2/6-31G'1.317103.1}}{}1 (Cs)3 B3LYP/6.31G'170.6 B3LYP/6-311G'172.4 UMP216-31G'170.2 ROMP2/6-31 G'1. 26810302 (C2v)1.1441.1361.3921.2541721.93139 B3LYP/6-31 G'131 B3LYP/6-311G'115 UMP2/6-31G'157 ROMP2/6-31G"1.3421.3391.3801.3273 (CS)34北京理工大学博士学位论文3 1.3525166.0166.9171.51653T}4}=106 B3LYP16-370'15. 7 B3LYP/6-311054. 4 UMP2/6-310*00 ROM P2/6-3 10 -4 (C2 or C2)105.98106.8B3LYP16-31G*1067.}7UMP2/56-3 I G*B3LYP16-3IIG*20117.8T312407.4T645{_;;;1.48675 (CS)118. 3 7.51133 68图3.1五个N7(1-5)异构体的优化几何具有Cs对称性的结构1由一个平面五元环和一个N2侧链组成。王 等[221预测异构体1是最低能量的N,异构体。我们的计算结果表明1在B3LYPl6-31 G*和UMP2r6-31 G*水平下具有全正的振动频率,是势能面上的局域极小点。但是当采用较高的6-311G*基组时,发现在优化过程中1离解为一个N3和两个N2碎片。另外,在ROMP216-31 G*计算水平下也没有找到稳定的异构体1。因此我们判断异构体1并不是N:势能面上真正的局域极小点。线性N7有三种异构体(具有C2,对称性的异构体2, Cs对北京理工大学博十学位论文称性的异构体3和C2或C2v对称性的异构体4),这三种构型都是N7势能面上的极小点。从结构上看,3相应于2中N4-N5-N6-N7沿着N3-N4轴旋转得到,并且键长和键角只有轻微的改变。异构体2和3几乎具有一致的能量:在CCSD(T)/6-311 G*//B3LYP/6-311 G*+ZPE(B3LYP/6-311 G*)水平下,3仅仅比2的能量高0.7kcal/mol。在B3LYP和UMP2计算水平上,4具有C2对称性,然而,在ROMP2水平上,构型4是一个平面分子,具有C2,对称性。在CCSD(T)/6-31 1G*//B3LYP/6-311G*+ZPE(B3LYP/6-311 G*)水平下,4的能量比2高10.9 kcal/mol。如图3.1所示,在构型4中,端键N1-N2(N6-N7)的键长是1.146人(B3LYP/6-31G*), 1.138人表3.1 1-5异构体总能量E (Hartrees)和零点能ZPE (kcal/mol)CCSD(T)a/B3LYP/6-31G*B3LYP/6-311G*E ZPEUMP2/6-31G*ZPEE ROMP2/6.31G*ZPE6-311G*吟思思)vzcE ZPE一83A5404-383.0303317.917.6一82.03926一83.12725一83.12509-383.1072617.517.316.519.220.420.320.9174一81.99309-383.02845-383.0112217.416.6-381.99041-381.99381一82.02187一82.02156-382.0100417.917.716.6-382.20728-382.20592-382.18827- 382.14338-382.95789 17.8 -383.04359 17.6 -381.90899(Cs) 单点能计算在CCSD(T)/6-31 1G*//B3LYP/6-311G*水平上表3.2 1-5异构体的相对能量(kcal/mol)IsomersB3LYFJ6-31G* B3LYP16-311G* UMP2/6-31G* RGMP2/6-31G* CCSD(T)&/6-311G*-14.60乃1 (Cs)2 (C2.-)3 (Cs)4 (C2 or C205 (Cs)1.011.045‘7︸-曰众八L︐乙11.552.6柳叻6.1挑众o认n6J0.00.710.940.2.相对能量在CCSD(T)16-311G*//B3LYP/6-31 IG*+ZPE(B3LYP/6-31 IG*)水平(B3LYP/6一3111G*), 1.150A(UMP2/6-31G*), 1.174A(ROMP2/6一31G*),比N3-N4(N4-N5)键长和N2-N3(N5-N6)键长短很多。从结构看,构型4相当于构型3中N5-N6-N7沿着N4-N5轴旋转得到,同时对称性由C¥改变为C2(或C2v)。从表3.3显示的频率可以看到,UMP2计算给出的构型3和4的最大振动频率分别为3328cm“和2962cm-1,这些数值显然太大了。这表明UMP2方法对于计算这种类型自由基的振动频率有较大的误差。Wang等[[22]报道两环结构5是N:势能面上的局域极小点。但是我们的计36北京理工大学博士学位论文算结果显示,在ROMP2/6-3I G*优化水平下,5离解为多个碎片。下面动力学稳定性分析也表明构型5的分解势垒很低,说明构型5确实并不是一个稳定的异构体。表33somers异构体3和4的振动频率((crir')和红外强度(km/mol)(括号中)3 (Cs)IIG* UMP2/6-31G* ROMP2/6-31G*mode B3LYP/6一IG*B3LYP/6-311311 4(0)29(1 )A113( 0)117( 0)128人125( 0)126( 0)231人140( 0)233(1 )234(1 )299A281( 0)236( 2)284( 0)445(1 )358A370( 5)366( 5)441645( 3)A424( 3)428( 2)A人AAAAAAA 7)482(622(11)773(12)954(11)1003(47)1141(37)1232(127)2166(408)2235(552) 6)480(629(11 )770(1 4)943(1 3)995( 48)1137( 41)1231(140)2154(446)2223(579)715(31)757(49)791(28)907(160)923(46)469647775940979108312251052(509)1159(442)2950(604)3328(22)231623614 (C2)a A 29(0) 48(0) 98(2) C2-CN A2,121A 155(0) 151(0) 168(0) A,, 167B 198(3) 186(3) 192(0) B2, 178B 272(0) 247(0) 369(37) B1,236B 362(9) 387(7) 495(13) A2,266A 377(0) 393(0) 531(4) A1,417A 421(1) 435(1) 533(7) B1,437B 614(11) 617(14) 724(24) B2,636A 867(7) 862(2) 816(11) B2,759A 878(45) 871(50) 1030(39) A,,821B 969(12) 948(20) 1040(458) Al, 931A 1081(5) 1083(12) 1272(18) Al, 980B 1111(58) 1100(66) 1568(17) B2,991B 2100(200) 2093(230) 2834(73) A1,2286A. .,} 2,166,呱211)___2154(245)一一2962(99)一一一B2, 2389一一a在ROMP2/6-31G*水平下,4具有C2,对称性。”ROMP2计算没有红外强度结果。3.3.2异构化和分解反应及其过渡态我们采用B3LYP和MP2方法对三个线型异构体2( C2,), 3(Cs)和4(C2或C2V)的异构化反应机理进行了研究。然而,我们只在¥3LYP/6-31 G*和B3LYP/6-31IG*水平下找到两个过渡态,TS213和TS3/4,如图3.2所示。IRC计算证实TS2/3确实连接N,势能面上的两个极小点2和3。在37 北京理1:大学博士学位论文CCSD(T)/6-31 1G*//B3LYP/6-311G* + ZPE(B3LYP/6-311G*)水平下,TS2/3的能量比构型3高6.7 kcal/mol。这个反应势垒很低的原因是由于从3异构化反应到2仅仅是N4-N5-N6-N7沿着N3-N4键旋转即可。构型4可以通过过渡态TS3/4到达3。基于TS3/4和4之间的能量差,这个异构化反应的势垒高在B3LYP/6-31 G*和B3LYP/6-311 G*(包括零点能修正)水平下分别为1.3 kcal/mol和0.8 kal/mol。然而,在CCSD(T)/6-311G*//B3LYP/6-311 G*+ZPE(B3LYP/6-311 G*)水平下,TS3/4的能量比构型4低。如此低的垒高说明了构型4是动力学非常不稳定的。2441.413106.611234=168.4168. 2 87.7 175.7178. 6B3LYP/6-311 G*12345=88.113456=175.86-31G*14567=178刀B3LYP/TS2/3 (Cl)\︑门/了1};{」:2551.250T5432=3.5 13456=63.0 B3LYP/6-31G*4. 5 61.0 B3LYP/6-311G*14321=178314567=166.4176. 8 165.0;:;TS3/4 (C1)图3.2过渡态TS2/3和TS3/4的优化几何结构38北京理工大学博十学位论文表3.4过渡态总能量E (Harteres),零点能ZPE (kcal/mol)和相对能量(E,1) (kcal/mol)B3LYP/6-31G*-383.03033 2 (C20383.01517TS2/3 (Cl)-TS3/4 (C,)-383.01009鲡ZPEI7.617.2I6.8B3LYP/6-311G*ZPE一83.12725-383.11212-383.1064117.517.116.717.0一.01912.3Iah一0妇一82.20728-382.19488-382.1898545.9512.37心一047.01I IG*单点能在CCSD(T)/6-31 IG*//B3LYP/6-311G*水平上。相对能量(E.0在CCSD(T)16-31 IG*//B3LYP/6-311G*+ZPE(B3LYP/6-311G*)水平上。 另外,我们也对最稳定异构体2的分解机理进行了研究。Nguyen和Ha[19]曾报道线型N5(C20是不稳定的,极易分解为N2+N3。因此我们认为线型异构体2分解放出两个N2的反应应该是一个协同分解反应。对于两个N2的消除反应,用B3LYP和RQMP2方法我们没有找到这个过渡态,一旦将对称性改变,便直接分解为多个N2碎片。但在UMP2/6-31 G*水平上,找到了过渡态TS-1(见图3.3)。在CCSD(T)/6-311 G*//UMP2/6-31 G*水平下,分解反应的能垒为34.8 kcal/mol,但是由于UMP2方法在频率计算中的误差,我们并不能得到2是动力学稳定的结论。Fan等[29]在计算N5自由基的分解时也发现基于UHF方法的MBPT(2)和MP2也有同样的缺点。Lauderdale和Bartlett[30]的研究表明ROHF-MBPT(2)方法对某些自由基显示了较好的计算结果。通过以上讨论,基于异构体2在B3LYP和ROMP2理论下无势垒分解,我们得出2是动力学非常不稳定的结论是合理的。一般说来,包括开壳层分子的过渡态应该采用多参考态的方法,比如MR-CISDo2刀86;2尹z产 户产UMP216-3TS-1(Cj图3.3 UMP2水平下TS-1的优化几何结构北京理工大学博士学位论文表3.5 TS-1总能量E (¥atrrees),零点能ZPE (kcal/mol)和相对能量(Ere,) (kcal/mol)UMP2/6-31G*CCSD(T)/6-311 G*E,e12 (C20Eat-382.202200.0TS-1(Cs)-381.99309-381.908380.050.6-382.1426334.8*相对能量(E.1)在CCSD(T)/6-311 G*//UMP2/6-31 G*+ZPE(UMP2/6-31 G)水平上。如前所述,构型5在B3LYP和UMP2水平下是势能面上的极小点, 但在ROMP2优化过程中分解为多个碎片。我们在B3LYP水平下对5的分解机理进行了研究,找到一个过渡态TS-2,如图3.4所示,从5转变到过渡态TS-2的过程中,随着N5-N7键伸长,闭合四元环打开。这一分解反应几乎不需要活化能,其势垒仅为0.5 kcal/molo IRC复核计算表明过渡态TS-2确实连接反应物5和产物碎片(一个N3和两个N2),三个键(N1-N3, N4-N5, N4-N6)断裂得到这些分解产物。5 L642=116.4L312=1 13.7Z421=106.7L542=105.105. 4 116.9113. 9 106.897.8 2325Z764=97.9924213=7.4T5421=-5.67. 5-5. 526421=-103.8 276420106.1 B3LYP/6-31G*-103. 2 106.0 B3LYP/6-311G*};;{TS-2 (CO图3.4 TS-2的优化几何结构3.4结论我们用UMP2, ROMP2和B3LYP方法对三个以前研究的N:异构体和两个新的N:异构体的结构和动力学稳定性进行了研究。用简谐振动频率分析了驻点是否是极小点或是过渡态,用内票反应坐标IRC解析了过渡态是否是连接反应物和产物的最低能量路径,得到以下结论:北京理工大学博十学位论文(1)以前报告的最稳定的异构体1在B3LYP/6-31 I G*和ROMP2/6-31G*水平下是不稳定的,在优化过程中即分解为多个碎片。 (2)三个线型异构体((2, 3和4)都是N7势能面上的极小点。构型3经过过渡态TS2/ 3到2的异构化反应势垒是6.7kcal/mol。构型4到3的异构化反应几乎不需要活化能。在B3LYP和ROMP2理论水平下, 最稳定的直线构型2无势垒分解,虽然UMP2的分解势垒很高,但 是UMP2在频率计算上的较大误差说明UMP2的计算结果并不可 靠。 (3)动力学数据显示,以前报告的两环构型5也不可能稳定存在。(4)我们的计算结果也显示,ROMP2方法对振动频率的计算结果要比UMP2方法好,似乎ROMP2方法更适合于计算这种类型的自由基。 北京理工大学博十学位论文参考文献1. Christe K O, Wilson W W, Sheehy J A, Boatz J A, Angew. Chem. Int. Ed. Eng., 1999,38: 20042. Christe K O, Dixon D A, McLemore D, Wilson W W, Sheely JA, Boatz J A, J.Fl uorine. Chem. 2000,101: 1513. Olah GA, Prakash G KS , Rasul G J, Am. Chem. Soc. 2001, 123:33084. Pyykko P, Runeberg N, J. Mol. Struct. (THEOCHEM) 1991, 234: 2795. Glukhovtsev M N, Jiao H, Schleyer P v R, Inorg. Chem. 1996, 35: 71246. Lauderdale W J, Stanton J F, BatrletR J, J. Phys. Chem. 1992, 96: 11737. Dunn K M, Morokuma K, J. Chem. Phys. 1995, 102: 49048. Lee T J, Rice J E, J. Chem. Phys. 1991, 94: 12159. Bitererova M, Brinck T, Ostmark H, J. Phys. Chem. A 2000,104: 1199910. Engelke R, J. Phys. Chem. 1992, 96: 1078911. Gagliardi L, Evangelisti S, Barone V, Roos B O, Chem. Phys. Let. 2000, 320: 51812. Tobita M, BartlettR J, J. Phys. Chem. A 2001, 105: 410713. Leininger M L, Sherrill C D, Schaefer H F, J. Phys. Chem. 1995, 99: 232614. Gagliardi L, Evangelisti S, Bemhardsson A, Lindh R, Roos B. O, Inter. J. Quant. Chem.2000,77:311 15. Tian A-M, Ding F-J, Zhang L, Me Y-M, Schaefer H F, J. Phys. Chem. A 1997, 101:1946 16. Chung G, Schmidt M W, Gordon M S, J. Phys. Chem. A 2000,104: 564717. http://www.qtp.ufl.edu/-batrlett/polyhnitrogen.doc.gz18. Gagliardi L, Orlandi G, Evangelisti S, Roos B O, J. Chem. Phys. 2001, 22: 1073319. Nguyen M T, Ha T. K, Chem. Phys. Let. 2001, 335: 31120. Wang X, Hu H.-R, Tian A-M, Wong N. B, Chien S-H, Li W-K, Chem. Phys. Lett.2000,329:483 21. Q-S Li, X-G Hu, W-G Xu, Chem. Phys. Lett. 1998,287: 9422. Wang X, Tian A-M, Wong N B, Law C-K, Li W-K, Chem. Phys. Lett. 2001, 338: 367北京理工大学博十学位论文23. Y-D Liu, J-F Zhao, Q-S Li, Theo. Chem. Acc. 2002,107:14024. Michels H H, Montgomery J A Jr, Christe K O, Dixon D A, J. Phys. Chem. A 1995, 99:187 25. Byrd E F C, Sherrill D, Head-Gordon M, J. Phys. Chem. A 2001,105: 973626. Gonzalez C, Schlegel H B, J. Chem. Phys. 1989, 90: 215427. Gonzalez C, Schlegel H B, J. Phys. Chem. 1990,94: 552328. Frisch M J, et al., GAUSSIAN 98, Revision A.9. Gaussian, Inc., Pitsburgh, PA, 199829. Fau S, Wilson K J, BartlettR J, J. Phys. Chem. A 2002,106:463930. Lauderdale W J, Stanton J F, Gauss J, Watts J D, BartlettR J,187:21 Phys. Lett. 1991,43北京理工大学博士学位论文第四章N7土原子簇分解机理的理论研究4.1引言C6 。等富勒烯发现之后,人们一方面继续研究原子数目更多、结构更新颖的碳原子簇,另一方面也努力开展其它的原子簇的研究。其中,氮原子簇的研究吸引了众多实验和理论化学家的广泛关注,一方面这不仅是由于氮原子簇本身的结构意义,而且更重要的是由于它们作为环境友好高能量密度材料(HEDM)的潜在应用价值。另外,氮原子簇分解放热时形成稳定的N:气体,如果用它作为火箭燃料,则可以自身产生工质,并且由于能量释放是通过分解反应而不是氧化反应,不需要携带氧化剂,因而若能使氮原子簇作为能源,将是火箭技术的一场革命,因此氮原子簇的研究不仅具有重要的学科意义,而且有广阔的高科技和军事应用前景,为各国科学家和政府所重视。然而,氮原子簇的合成一直是一个极大的挑战。在1772年首次分离 出N:分子、1890年合成出叠氮阴离子N3一以后,直到1999年,Christe等[1]才成功地合成出第三种稳定的聚氮化合物N5+盐(N5气AsF6r)。最近,Christe和合作者[2]也报告了另外一种稳定的NS+P-t', NS+AsF6-,这种新的盐在70℃分解,并且抗冲击性较强。Olah[3]等人试图在一78℃由氨基重氮盐制取N42+离子,然而,遗憾的是,他们没有得到N;的二价阳离子。在过去几年里,理论上己经预测了很多稳定的聚氮化合物结构 [4-21]o Schleyer等人[7]用Ab intio和DFT方法计算了中性偶数氮原子簇(N2, N4, N6, N8, Nio和N12)的结构和能量。Bartlett等人[22]对N2到Nio氮原子簇以及它们的正、负离子的结构、能量和光谱性质进行了系统的研究,并且指出大量的聚氮化合物正等待人们去发现。Gagliardi等人【23]在他们最近的一篇论文中报告了NS一负离子有足够的稳定性并且进一步研究了从N3-, N5一和NS+合成N:和N,。的可能性。在氮原子簇的合成研究过程中,科学家们对偶数氮原子簇寄予了更多 的希望。然而,有趣的是,到目前为止,只有N3一和N5+被成功地合成。北京理工大学博士学位论文这启示我们奇数氮原子簇的合成也许是更有希望的。显然,N:原子簇应该是下一个合成的目标。最近,我们【25]采用从头算和密度泛函理论对四种N7+和五种N7一异构体的几何、能量、振动频率和成键进行了研究。振动频率分析的结果表明九种N7t离子异构体在HF16-31 G* . 3LYP/6-31 G*.MP2/6-31 G*和B3LYP/6-311(+)G*水平下都是稳定点或接近稳定点。其中,最稳定的N7+异构体为开链线型结构,具有C2}对称性,与己知的奇数氮原子离子(N犷和N9十原子簇)的最稳定构型相似。最稳定的N7一异构体也为开链的线型结构,具有C2}或C2的对称性,而且N7土离子原子簇相对于7/2 N:的能量表明,N7士离子可能被合成并可能是潜在的高能量密度材料。但是,作为高能量密度材料候选者,对氮原子簇来说,最关键的问 题是分解的势垒和合成的方法。分解势垒必须足够高才能提供燃料正常储存的稳定性。一般情况下分解势垒的高度应大于20 kcal/mol,最好是大于30 kcal/mol。因此,对任何可能的高能密度材料候选者都必须要探讨它们可能的分解途径进而确定分解势垒。以前只有一些偶数氮原子簇如N4[8-11], N6[12-14]和N8[17-18]的分解机理的研究报道,本研究的目的就是采用Ab initio和DFT方法,在前面工作的基础上,进一步对N7+和N7一原子簇的分解机理进行研究进而确定这些反应的势垒,即确定它们的动力学稳定性。我们希望我们的研究能对将来氮原子簇的寻找和合成提供一些理论基础。4.2计算方法所有的计算都应用Gauss ian 98程序包[31]完成。反应物、产物和过渡态的几何用两种不同的方法,即MP2和DFT进行了全优化。本文应用的DFT方法是一种HF/DFT杂化的方法B3LYP[26-27], MP2方法采用固核近似。对于阳离子N犷,应用的是6-31 G(d)基组,而对于阴离子N7-,应用是弥散基组6-31+G(d)。同时在以上的理论水平下进行了振动频率的计算,用来表征稳定点和进行零点能校正。所有稳定点已经被确定为极小点(虚频个数=0)或过渡态(虚频个数=1)。连接反应物和产物的最小能量途径用内察反应坐标方法[28-29]进行了证实。北京理I几大学博士学位论文为了获得更可靠的分解势垒,相对能量计算应用了与G2( MP2)[30]方法类似的基组外推法,在这里标记为G2*(MP2)。对MP2/6-31(+)G(d)几何,全部能量采用如下公式计算:E [G2*(MP2)] =E [MP2/6-311G(d, p)]+)MP2+AEQC + ZPE+AE HLC+eLC这里,QMP2=E [MP2/6-311+G(3df, 2p)]-E[(MP2/6-311G(d, p)],)EQCI=E [QCISD(T)/6-311G(d, p)] -E[MP2/6-311G(d, p)],DEHLC=-0.00019rla 0.00595rlp, and OHLC=0.00114lip,11n, Ilp分别是a和R价电子的数目。显然,G2*(MP2)能量和直接在QCISD(T)/6-31 1+G(3df, 2p)//MP2/6-31(+)G(d)水平下计算的能量近似。为了比较,对B3LYP/6-31(+)G(d)优化几何也用较大的基组6-311+G(3df,2p)进行了单点能计算。本文所有的的能量均包括零点能修正。4.3结果和讨论每个异构体的分解途径将分别在下面的各部分进行讨论。为了方便, 字母C和A分别代表阳离子和阴离子。优化的N犷和N7-稳定点和分解过渡态的几何结构分别显示在图4.1、图4.2和图4.3,最低振动频率列于表4.1,相对能量和零点能列于表4.2.表4.1异构体i (SCs)各异构体和过渡态的最低振动频率(单位:cm')MP2/ 6-31(+)G(d) B3LYP/6-31(+)G(d)-50944-030169-565N,一A-1(C2)TS2 (C0A- 2(Cs)TS3-a (C,)A一( C,)TS3-b (C,)TS3 (Cs)_-71358210-75567北京理工大学博士学位论文表4.2相对能量(Ere3) (kcallmol)和零点能ZPE (kcalmol)异构体MP216-31(+)G(d)ZPEEr<I八U月护G2"MP2E.10.0B3LYP16-31(+)G(d)ZPEErelOO18316.716.216.4B3LYPI6-311+G(3df,2p)00.04.8-44.4N7' C-1(C2v)aNs, (C2v) + N2N,一^-1 (C2)' (TS2 CI)A- 2 (Cs)TS3 (Cs)A- 3 (CI)TS1( Cs)-56.2000.317.816.015石一-353.1.00.04.5-40.816.615.97219.417.516.41.210.715.016.0l7.9TS3- a (C I)4. 317石16.416.612.524.116.829.0 TS3-b (C I)A-Complex (Cs)N3 +N2+N213.1-5.9-161.2 5.51498. 6-140.912.223石16.517.617.113.715.226.520.8一131.618.2-127.6a: N7-和N?一的最稳定异构体的能量分别设为零能量。b' G2*(MP2)//MP2/6-31(+)G(d)单点能用MP2/6-31(+)G(d)的零点能修正.‘B3LYP16-31 I+G(3df,2p)//B3LYP16-31(+)G(d)单点能用B3LYP/6-31(+)G(d)的零点能修正a4.3. t阳离子N广的分解机理MP2J6-31G(d)31G(d)c-1 (G,,)戮562M于,216-31G(d),矛I.石351 597151.9 B3LYPf6一(d)TSl (CS)图4.1优化的N7+ C-1和TS1的几何结构和参数从前面的研究[24]可知,C2、对称的开链线性结构C-1是所有N7十异构体中最稳定的。从图4.1我们可以看到,在MP2/6-31 G(d)和B3LYP/6-31 G(d)47北京理工大学博士学位论文水平下,Nl-N2(N1-N3)和N2-N4(N3-N5)的键长分别为1.311A, 1.302A和1.327A, 1.321A。这些键长位于HN=NH中的N=N双键(1.311.)和芳香NN键(1.350人)之间,表明这些键在某种程度上有一定的共辘性。终端键角(L 246和L 357=166.80)接近1800,而其它的键角接近1200e在MP2/6-31 G(d)和B3LYP/6-31 G(d)水平下,我们找到了一个过渡态TS1(见图4.1)。如表4.1所示,在上面两种水平下,过渡态TS1虚频的大小分别为一595和一509。从图4.1可以看到,在由C-1转变为TS1的过程中,N3-N5键拉长并且对称性由C2,变为Cs,同时其它的键长和键角只有很轻微的改变。这些暗示N3-N5键很容易断裂,从而N厂离解为NS++N2。在以上两种水平下完成的IRC计算证实,在N广势能面上,TS1确实连接C-1和NS++N2。在B3LYP/6-311+G(3df,2p)//B3LYP/6-31 G(d)水平下,TS1和N7+(C2})的能量差为4.8 kcal/mol, G2"(MP2)降低这个差值到3.1 kcal/mol(见表4.2)。如此小的能垒高度表明C-1是动力学不稳定的,不适合作为高能密度材料候选者。Nguyen等[14]己经讨论了阳离子N犷分解为N3++N2的离解机理并且指出开链N¥+具有一个相对高的分解能垒((170-220kJ/mol) e4.3.2阴离子N7一的分解机理Mi chels等〔19〕研究了C2和CS对称的N(N3)2一的结构和稳定性。在前面的研究中〔24],我们己经指出在MP2和DFT水平下,C2构型的线性异构体A-1是N7一势能面上的最小点。从图4.2我们可以看到,对异构体A-1来说,最长和最弱的键是中心Nl-N2(Nl-N3)键,它的长度接近在N214中N-N单键的键长(rN.N=1.45劫o Mulliken电荷分布分析表明,负电荷主要集中在中心氮原子上,这暗示A-1的离解机理是N-N的断裂而不是N2的直接消除。我们已经找到的一个过渡态TS2 (C,)(见图4.2),在MP2/6-31+G(d)和B3LYP/6-31+0(d)水平下,虚频分别为一403和-138(见表4.1)。在B3LYP/6-311+G(3df,2p)//B3LYP/6-31+G(d)水平下,TS2的能量仅仅比A-1高5.0 kcal/mol,而02*MP2降低这个能垒值到1.2 kcal/mol(见表4.2)。和稳定的异构体A-1相比,TS2中Nl-N2, N3-N5键长增长,而北京理工大学博士学位论文N1-N3键长缩短。从TS2开始,IRC计算显示,NI-N2和N3-N5两个键同时断裂形成N3一和两个N2分子,NI-N2和N3-N5这两个键同时断裂也许是由于线性N;不稳定直接分解为两个N2分子造成的。然而,如此低的分解势垒表明,实验能否观测到这个异构体是非常值得怀疑的。174.32/6-31+G(d)LYP/6-31+G(d)A-I {C2)1748厂入,‘3I+G(dB3LYP/6刁卜G(d)丸213=164.5T5312=135名173力 86. 2TS2 (c1)T7531=110.5T6421=180.2171. 5144.3 图4.2优化的N7一 A-1和TS2的几何结构和参数我们前面的研究[ [24] E经显示,具有C:对称性的环状异构体A-2和笼状异构体A-Complex都是N7_势能面上的局域极小点(见图4.3)。从表2可以看到,包括零点能修正,在B3LYP/6-31I+G(3df,2p)//B3 LYP/6-31+G(d)和G2*(MP2)水平下,A-2的能量t匕A-1分别高23.5 kcal/mol和10.7 kcal/mot。上面两种理论水平都预测A-Complex和A-2异构体具有相似的能量,A-Complex的能量比A一分别低2.7 kcal/mol和2.1kcal/mola Mullik。电荷分布分析表明,多余的负电荷主要集中在N2和N3上。从结构观点看,A-Complex可以看作是N。环和叠氮阴离子N3-之间的一个弱的分子一离子络合物。我们也在MP2/ 6-31+G(d)和B3LYP/6-31+G(d)水平下考查T环状异构北京理工大学博十学位论文体A-2的分解机理。从计算的结果我们发现,对于A-2而言,它的分解机理强烈地依赖于所使用的计算方法。通常情况下,对氮原子簇来说,MP2方法由于过多地考虑电子相关能的影响从而导致一些计算的键长比实际键长长。这似乎也是在MP2和B3LYP两种计算方法下,A-2有不同的分解机理的原因。在MP2水平下,我们找到了一个过渡态TS3。从图4.3中我们可以看到,与稳定的异构体A-2比较,TS3仍然具有CS对称,只是N2-N4和N3-N5键长增长了大约0.3人,同时N4-N5的键长大大的缩短。在MP2/6-31+G(d)水平下完成的IRC计算证实,通过A-Complex,A一分解为N3-+2N2。这个分解过程可以表示为:A-2(Cs)--TS3(CS)-A-Complex-N3-+2N2。然而,在B3LYP/6-31+G(d)水平下,A-2转化为A-Complex是通过一个异构化反应完成的。如图4.3所示,在B3LYP势能面上我们找到了两个过渡态TS3-a和TS3-b. IRC计算证实,过渡态TS3-a连接势能面上的两个稳定点A-2和A-3,过渡态TS3-b连接A-3和A-Com川ex。从表4.1可以看出,具有C.对称性的异构体A-3在MP2和B3LYP水平下具有全正的振动频率,也是N7-势能面的极小点。在B3LYP/6-311+G(3df,2p)//B3LYP/6-31+G(d)和G2*MP2水平下,A-3的能量比A-2分别低8.3 kcal/mol和5.2 kcal/mol(见表4.2)。在B3LYP/6-311+G(3df,2p)// B3LYP/6-31+G(d)水平下,TS3-a和A-2的能量差为5.5kcal/mol,这个较低的势垒高度表明这个异构化反应很容易发生。另一方面,TS3-b和A-3的能量差为11.3 kcal/mol,这表明在特殊的实验条件下可能能探测到异构体A一的存在,但是作为亚稳态的高能密度材料,显然它不是合适的候选者。总之,MP2和B3LYP方法都预测环状阴离子A-2不具有足够的动力学稳定性。北京理工大学博士学位论文MP2/6-31+G(d)B3L)甲/6.31+G(d)叼︐1.744︐A-2(Cs)TS3 (Cs) [MP2/6-31+G(d))MP216-31+0(d)B3LYP/6.31+G(d)1357325TS3-a (CI) [B3LYP/6-31+G(d)]1.325A-3 (C 1)842.231.1韵1.215 MP2/6-31 +G(d)184田LYP/6-31+G(d)69.4 MP216.31+G(d)冲173.4 B3LYP16.31+G闭)‘︑63:{拿1}102.0︐淤2沸︑︑.︑白﹂弓‘2.404;﹄·︑2539898石·︑2.587;·﹄八︐飞划‘2.559·-1.49838叮284的A-Complex (Cs)1.285TS3- b (C 1)图4.3优化的h17一 A-2, A-3和过渡态TS3, TS3-a, TS3一的几何结构和参数5I北京理工大学博士学位论文4.4结论在以前热力学稳定结构研究的基础上,本文进一步用从头算和密度 泛函方法对低自旋N7+和N厂原子簇的动力学稳定性进行了研究。计算结果表明:(1)热力学最稳定的线性N7-(A-1)异构体离解时,发生协同分解反应,同时消除两个N2分子( N7-(CZ)-TS2--N3一 +N2+N2)。在我们最高的 计算水平下,G2*(MP2)计算的这个分解反应势垒是1.2 kcal/mole而对同样是线性结构,N7 +(C-1)异构体消除两个N2分子是一个分步过程。对N, +(C2V)-" TS1--NS++N2分解反应,G2*(MP2)计算的势垒是3. 1 kcal/mol。如此低的势垒表明热力学最稳定的线型N7+, N7-异构体都不具有足够的动力学稳定性。 (2)对具有CS对称性的七元环状异构体N7-(A-2)来说,MP2和DFT计算的分解势垒都小于20 kcal/mot,表明它也不是一个具有实际意义的亚稳态类型。 北京理工大学博士学位论文参考文献1. Christe K O, Wilson W W, Sheehy J A, Boatz J A, Angew. Chem. Int. Ed. 1999, 38:2004 2. Vij A, Wilson W W, Vij V, ThanFr S, Sheehy J A, Christe K O, J. Am. Chem. Soc. 2001,123: 63083. Olah G A, Prakash G KS , Rasul仅J. Am. Chem. Soc. 2001, 123: 33084. Glukhovtsev M N, Schleyer P二凡Chem. Phys. Lett 1992, 198: 5475. Vogler A, Wright RE, Kunkley H, Angew. Chem. Int. Ed. 1980,19: 7176. Engelke R, J. Chem. Phys. 1992,96:107897. Glukhovtsev M N,Jiao凡Schleyer P v R, Inorg.Chem. 1996, 35: 71248. Lee T J, Rice J E, J. Chem. Phys. 1991, 94: 12159. Dunn K M, Morokuma K J, J. Chem. Phys. 1995,102:490410. Lauderdale W J, Stanton J F, BartlettR J, J. Phys .Chem. 1992, 96: 117311. Bitererova M, Brinck T, t'jstmark H, J. Phys.Chem. A 2000, 104: 1199912. Korkin A, Balkova A, BartletR J, Boyd R J, Schleyer P v R, J. Phys.Chem. 1996, 100:5702 13. Gagliardi L, Evangelisti S, Barone V, Roos B O, Chem. Phys. Letf. 2000, 320: 51814. Nguyen M T, Ha T.-K, Chem. Phys. Lett. 2001, 335: 31115. Gagliardi L, Evangelisti S, Widmark P O, Roos B O, Theor. Chem. Acc. 1997, 97: 13616. Gagliardi L, Evangelisti S, Roos B O, Widmark P O, J. Mol. Struct. (THEOCHEM)1998,428:1. 17. Gagliardi L, Evangelisti S, Bemhardsson A, Lindh R, Roos B O, Int. J.Qutantum.Chem.2000,77:311 18. ChungqSchmidt M砚Gordon M S, J. Phys.Chem. A 2000,104: 564719. Michels H H, Montgomery J A Jr, Christe K O, Dixon D, J. Phys. Chem 1995, 99: 18720. Nguyen M T, Ha T.-K, Chem.Ber. 1996,129:115421. Leininger M L, Sherrill C D, Schaefer H F, J. Phys. Chem 1995,99: 232422. R J Bartlet, http: //www.gtp.ufl.edu/-bartlett/ polynirtogen .doc.gz北京理工大学博十学位论文23. Gagliardi L, OrlandiqEvangelisti S, Roos B O, J. Chem. Phys. 2001, 114: 1073324. Y Dong Liu, J-Fang Zhao, Q-Shu Li, Theor. Chem. Acc. 2001, 107: 14025. Q-Shu Li, L-Jie Wang, J. Phys.Chem. A 2001, 105: 120326. Becke A D, J. Chem. Phys. 1993,98: 564827. Lee C, Yang W, Parr RqPhys. Rev. B 1988, 37: 78528. Gonzalez C, Schlegel H B, J. Chem. Phys. 1989, 90: 215429. Gonzalez C, Schlegel H B, J. Phys. Chem. 1990, 94: 5 52330. Curtiss L A, Raghavachari K, Pople J A, J. Chem. Phys. 1993, 98: 129331. Frisch M J, et al. Gaussian 98, Revision A.9.; Gaussian, Inc., Pitsburgh, PA, 199854北京理工大学博士学位论文第五章5.1引言Nt:原子簇势能面的理论研究作为一种新型环保型高能量密度材料,氮原子簇在近年来吸引了众 多实验化学家和理论化学家的广泛关注「1-29]。然而,它们的合成一直是一个极大的挑战。到目前为止,已知的聚氮化合物只有Net N3一和N犷,其中第三种稳定的聚氮化合物一N5+盐是1999年由Christe小组[1]合成。这个新的发现为合成更多稳定的氮原子簇打开了光明的前景。Gl ukhovtsev等人[8]采用从头算和DFT方法计算了中性氮原子簇N4,Na, Na, Nro和Nit结构和稳定性·最近,Bartlett[9]在一篇综述中指出聚氮化合物的结构是迷人的并具有潜在的爆炸性,他预测大量的聚氮化合物正等着人们去发现,关键的问题是合成的方法和分解的势垒。Lauderdale和Stanton[10]采用SCF, CC(相关簇)和MBPT(多体微扰)方法计算了四种N4构型的能量和振动频率,指出正四面体N4(Td)在势能面上有一局域极小,是一种亚稳态分子,但遗憾的是正四面体烷分子至今仍没有被实验观察到,同时他们也分析了由立方体N8(O砂分解为4N:的能级相关图,指出N¥--4N2是Woodward-Hoffmann禁阻反应,表明立方烷N。应该有一个明显的分解能垒高度,如果能被合成,它将是稳定的。Gagliardi等人[11]和Nguyen等[121分别研究了N。的分解机理,发现最稳定线性二叠氮Ns异构体的分解反应是协同反应并且中等大小的活化能表明在室温下只具有很短的寿命。Chung等人[171报道了三个低能量的N:异构体的分解机理。他们的计算结果显示这三个低能量的N。异构体的分解垒高大都约在20 eom/lack2001年春,Ol ah等人[24]在J. Am. Chem. Soc.杂志上撰文预测二价阳离子N6 2+ . N4 2+与叠氮阴离子(N3_)可发生放热反应分别生成N12和N1oa前面的理论研究己经报道了四个热力学稳定的N12异构体,但迄今为止,还没有对它们分解机理和合成途径的研究报道,而这对于决定N12原子簇是否是有价值的HEDM是非常有意义的。在本章的研究中我们第一次报55 北京理工大学博士学位论文道了两个新的N12异构体。另外,为了更进一步评价它们的动力学稳定性,采用DFT方法对这六个异构体可能的分解途径进行了研究。希望目前的计算结果能为将来更进一步的实验合成提供一些理论数据。5.2计算方法所有计算采用Gaus sian 98量子化学软件包完成【30]。几何优化是在B3LYP[31-32]/6-31 G*水平上完成的,在同样水平下完成的振动频率计算用于表征驻点的性质:极小点(虚频个数=0)或过渡态(虚频个数=1),同时频率计算给出零点振动能(ZPE),这些数值可以用来对体系的总能量进行校正。在同样水平上对各过渡态进行了内票反应坐标(intrinsic reactioncoordinate, IRC)[33-34]复核计算,以确认过渡态结构与相应的异构体是否正确连接。在B3LYP/6-31 G*几何基础上进行了B3LYP/6-311+G(3df,2p)单点计算,并加上B3LYP/6-31 G*零点能校正。5.3结果和讨论5.3.1 N12异构体图5. 1显示了六个异构体的几何结构和计算的几何参数。全部能量、零点能(ZPE)和最低振动频率列于表5.1。从表5.1可以看到这六个异构体有全正的振动频率,都是N12势能面上的局域极小点。表5.1构型1-6全部能量(hartrees),零点能((kcal/mol)和最低振动频率((cm 1)ZPE一一一一一-一-一-一一一B3LYP/6-31 G*一一-一.)B3j访PLYP厉二/6-盯砰面〕f311+G(3R,. 2p)//B3LYP/6-31 G*内︐︐矛O门jU1 (Cm)3 (Cn)4 (Dad)5 (Dad)2 (Cs )-656.70563-656.6815335.333.932.331.632.831.2-656.91946-656.65036-656.58891-656.38734勺︵]1‘-656.90909-656.88599-656.80932-656.59997-656.36743294208北京理工大学博十学位论文表5.2构型1-6与6N:分子的相对能量“(kcal/mol)6NI 1 2 3 4 5 60. 0 316.4 320.5 334.4 381.8 5144658.7a.单点能计算在B3 LYP/6-311+G(3df,2p)//B3LYP/6-31 G*+ZPE(B3 LYP/6-31 G*)1.3821 (C2n)111.012852 (C办3 (C2n57北京理工大学博士学位论文1.421z2-1-10.9-111. 44 (D3砂5 (Did)‘(D6h)图5.1六个异构体(1-6)的优化几何Gl ukhovtsev等[8]曾报道构型1和3是两个低能量的N12异构体,Li等也报道了笼状N12异构体5和6,异构体2和4是我们报告的两个新构型。从表5.2可以看到,这六个异构体的稳定性顺序是1>2>3>4>5>6。以前的研究【8,291显示,对一些较大的偶数氮原子簇来说,芳香五元环是一个稳定的结构单元。对N12原子簇来说,构型1含有两个五元环结构表现了最大的稳定性,含有一个五元环和一个N7侧链的构型2是第二稳定的异构体,在B3LYP/6-311+G(3df,2p)//B3LYP/6-31G*水平下(包括零点能修正),它的能量仅仅比1高4.1 kcal/mol。如图5.1所示,除了端键N11-N12的键长为1.131A,比氮气中的氮氮三键的实验值(1.o98A)稍长,端键键角接近180。外,其余的键长和键角大都接近于芳香氮氮键的键长((1.35A)和58 北京理工大学博士学位论文120。的键角。与小氮原子簇直线构型最稳定的规律不同,C2 h对称的直线构型3的能量比1和2都高,在B3LYP/6-311+G(3df,2p)水平上,3的能量比最稳定的构型1高18.0 kcal/mol,如图5.1显示,构型3中,除端键N I-N2(或NI I-N12)的键长比氮气中的氮氮三键的键长长一些,端键键角接近1800外,其余的键长基本位于单双键键长之间,键角接近10900构型4是一个具有Da d对称性的十二元环的环状异构体,能量比1高65.4 kcal/mol。从图5.1显示的键长数据可以看出,构型4中的十二根键为单双键交替排列,键角非常接近于理想的¥p2杂化角的值(1200)a除了线型和平面型结构的异构体外,我们也对体状即多面体N原子 簇给予了更多的关注。这是因为多面体N原子簇全部由三键连的N-N单键构成,当它们分解放出N:时将会释放出比线型和面型异构体更多的能量。如果它们是稳定的并被合成,将会是非常有价值的高能量密度材料。除了图1中的两个笼状异构体外,我们也设计了更多由三元环、四元环、五元环以及六、七元环组成的笼状结构,如图5.2所示。但是很遗憾,这些结构仅在HF方法下是稳定的,如果采用更高的算法,比如DFT方法,在优化过程中构型发生分解,产生多个N2碎片。北京理一E大学博士学位论文图5.2设计的一些笼状N12异构体在B3LYP水平下,如图5. 1所示的构型5和6是两个稳定的高对称性笼状异构体,从表5.2可以看到,它们的能量大大高于前面四个线型和面型异构体。显然一方面环张力降低了构型5和6的稳定性,另一方面,正如Ha等人[[181研究笼状N2。时指出的那样,氮原子上的孤对电子的相互作用也是笼状结构高能的一个原因。相对稳定性也可从几何参数上反映,笼状异构体全部由弱的单键组成,而其余能量较低的四个异构体中既有单键也有双键。在B3LYP/ 6-31 1+G(3df,2p)水平下,这六个异构体与6个N2分子的相对能量列于表5.2,显然都是高能量分子。5.3.2 N,2异构体的分解途径从结构分析,N1 2异构体的离解有两种途径:一是直接消除一个或多个N2分子,二是先异构化为其它的构型再发生离解反应。我们采用DFT方法在B3LYP/6-31 G*水平上研究了以上六个稳定异构体的分解机理。过渡态(TSl-TS11)的全部能量、零点能((ZPE)和最低振动频率列于表5.3。表5.4列出了极小点(1, 2, 3, 4, 5, 6和两个Ni o异构体)与它们各自相应的过渡态的能量差。北京理工大学博士学位论文表5.3过渡态(TS1-TS21)和分解产物(Nlo, NS, N6)的全部能量[E(hartrees)],零点能[ ZPE(kca/lmol)]和最低振动频率[”1 (Cm 1)]B3LYP/ 6-31G` B3LYP/6-31 I+G(3df,2p)I somersZPEUl//B3LYP/6-31 G'TSI (Cs)-656.6892933.2-494-656.90973TS2 (Cs)-656.6583432.1721-656.88498TS3 (Cs)-656.6626031.7一98-656.89086TS4 (Cs)-656.6251530.5-265-656.86004TS5 (CZ,,)-656.61734287一44-656.85238TS6 (Cs)-547.2190626石一11-547.40719TS7 (Cs)-547.2238526.0-504-547.41337TS8 (Cs)-547.1885024.8刀47-547.38448TS9 (Cl)一56.5859729.8-442-656.80636TS10 (C,)-656.3815731.2-415-656.59344TSIl (C,)-656.1537330.3-351-656.36685-547.2459601NNia (Cs)45-547.43498-54721614-547.41322翰(C2n)(Cs)一37.79857招价24-437.94960437.777752省.‘肠-437.93599眺(D(C22hn))-437.776402凡j‘U1l642-437.92416乃凡﹂-328.44399表5.4过渡态与对应的稳定异构体的能量差趣〔包括ZPE修正)IsomersB3LYP/6-31 G*B3LYP/6-31 1+G(3df,2p)a/lB3LYP/6-31G*1(C2n)nn﹄U0.0TSI po勺妇4.02 (Cs)︸U们0.0 TS212名15.1TS3 9.711.43 (C2n)0.00.0 TS414.014.5TS5 17.117.54 (D3d)0.00.0TS9 0.050.065 (D3a)0.00.0TS10 2.0256 (D6h)0.00.0TS11 刁7-0.5Nio (Cs)0.00.0TS6 15.115.6TS7 11.511.2Nto (C2h)18_712.1TS吕34.028.1厂蔺蘸弄毓疏丽不藏下花石云2p)/B3LYP/6-31G*+ZPE(B3LYP/6-31G*)水平上61北京理工大学博士学位论文5.3.2.1构型1的分解机理:l (Cs) -. TSI -. N1o (Cs) +N21.370TSI (C劝图5.3过渡态TSI的优化几何从上面的讨论中我们已看到,构型1是最稳定的Ni t异构体。对1来说,可能的分解途径有两种:NS环断裂或NS环先打开、再断裂。我们没有找到N,环打开的过渡态,而找到了过渡态TS1相应于NS环断裂直接放出一个N2a TS1具有C¥对称性而不是Klapotke[35]报道的C:对称性。如图5.3所示,在TSI中,N4-N12和N8-N 10键长为1.685和1.676A,暗示N4和N12(N8和Nl0)之间的静电作用力很弱,因此这两个键很容易断裂放出N2分子.包括零点能修正,TS1的能量比1高8.2 kcal/mol,B3LYP/6-311+G(3df,2p)单点计算将这个值降低到4.0 kcal/mol。如此低的垒高表明制备这个分子可能性非常小。5.3.2.2构型2的分解机理:2 (Cs) -" TS2 -- N2 o (CS) +从2 (Cs) ̄TS3 ̄Nio (C2id十怂平面构型的异构体2与最稳定的N8几何构型很相似,由一个NS环和 一个N:侧链组成。从结构上分析,其分解途径包括两种方式:NS环断裂放出N:或侧链N-N键断裂放出N:分子。我们找到了两个过渡态相对于这两种不同的离解方式。一方面,构型2可以通过过渡态TS2断裂北京理工大学博士学位论文TS2 (Cs)月.才︐;1.693\-.TS3 (Cs)图5.4过渡态TS2和TS3的优化几何侧链上的N-N键放出一个NZ分子,同时生成含有一个Ns环和一个N5侧链的平面N,o(Cs)分子,这个分解反应的势垒高度是12.8kcal/mol,B3LYP/6-311+G(3df,2p)单点计算给出一个稍高的值,15.1 kcal/mol。如图5.4所示,与稳定的构型2相比,TS2中N11-N12键长缩短0.014A,断裂的N10-Nll键长增长0.3A.另一方面,构型2五元环上的N3- N4和N1-N5键可以通过过渡态TS3断裂五元环上的N-N键,经过9.7 kcal/mol(B3LYP/6-31G*+ZPE)的势垒产生N2和一个线型NFo,而B3LYP/6-311+G(3df,2p)单点计算势垒高度值是11.4kcal/mol。如图5.4所示,N3-N4和N1-N5键长分别从1.376和1.347A增加到过渡态TS3的1.693和1.697A,而其他的键长和键角改变很小,暗示N3-N4和N1-NS这两个键很容易断裂从而放出N:分子。我们注意到,上述离解产物之一,平面Nl o(Cs)分子构型与2相似,由一个N,环和一个NS侧链构成,据我们所知,以前没有对这个异构体的北京理工大学博士学位论文报道。我们的DFT计算结果显示,平面N1o(Gs)是NIo势能面上的一个局域极小点。对其分解途径的研究也与2相似,可以从环上断裂放出N2和从侧链断裂产生Nz两个分解机理来进行。OiN (CS)的分解:NJO (Cs)--+TS6-N8(Cs)+N2 Nto(Cs)-TS7-Na(C2)h+N2、1.575TS6 (Cs)1.679723T57 (CS)图5.5产物N1武C好和过渡态TS6-7的几何结构如图5. 5所示,TS‘和T57分别为侧链N-N键断裂和五元环断裂的过渡态,在B3LYP/6-31G*水平上,垒高分别为15.1和11.5 kcal/mol,提64 北京理工大学博士学位论文高的B3LYPJ6-31 1+G(3df,2p)单点计算值为15.6和11.2 kcallmol。在Nlo(Cs)经TS6离解到N8(Cs)过程中,N8-N9键增长了0.3.E,而端键N9-N10键缩短了大约0.04人。显然,N8-N9键会很容易断裂放出一个N2分子同时生成最稳定的N武Cs)。另一方面,N1。分子中的五元环也可以断裂经过渡态TS7离解为直线N以C2h)和一个N2分子。分析构型2和具有相似结构的Nl o(Cs)的离解反应我们发现,含有五元环和N-N侧链的异构体的离解通道有两条:五元环-N-N键断裂或侧链N-N键断裂。我们的计算结果表明:五元环中N-N键断裂放出N2比侧链断裂放出N2的势垒稍低。所以含有N,环结构的氮原子簇异构体在发生离解反应时,应该优先发生环断裂反应。N1o (C2h)的分解:Nio (C`zid - TS8 --> Ns(c2d + Nz另外,我们也对产物,线型Nl o (C2h)分子的分解途径进行了理论计算,找到了过渡态TS8,如图5.6所示。包括零点能修正,TS8的能量比线型Nlo(C2h)高15.3 kcallmol。以前的论文[171己对NS(Cs)和线型N8(C2h)的分解途径进行了研究并指出它们的分解势垒都低于20 kcallmol,这里我们就不再对其进行研究。N1o (C2h)、一2、1.5681103107.0118.7‘、、之又44.7TS8 (Cs)图5.6产物Nlo (C2h)和过渡态TS8的几何结构北京理工大学博士学位论文5.3.2.3构型3的分解机理:3(C2h) --> TS4 - Njo(C2h) + N23(C2h)*TS5*3N2+2N3106.3105.11103TS4 (Cs)TS5 (C2h)图5.7过渡态TS4和TS5的优化几何对开链构型3来说,我们找到了两个过渡态TS4和TS5,如图5. 7所示。与稳定构型3相比,在TS4中,终端N10-NII键长增长了大约0.3A,而N11-N12键长缩短了0.014A,暗示N10-Nll键很容易断裂放出一个N2分子,同时生成线型NIo(C2h)分子。在B3LYP/6-311+G(3df,2p)//B3LYP/6-31G*水平上此离解反应的势垒是14.5 kcal/mol,与线型NIO(C2h)放出一个N:分子的势垒接近。另外,开链构型3也可以经过渡态TS5直接分解为3个N2和2个N;碎片,活化势垒为17.5 kcal/mol。如图5.7所示,在TS5中,链中的两个键N6-N7和N7-N8同时从1.407A增长到1.845A,而对称性保持不变。IRC计算证实N3-N4, N5-N6, N7-N8, N9-N10同时破裂产生多个产物碎片。其中N3-N4和N 9-N 10键的断裂也许是由于线型NS自由基不稳定造成的。我们曾试图寻找线型构型3发生均裂,分解为两个最稳定线型N6分子的可能分解途径,但遗憾的是,我们没有找到这北京理丁大学博士学位论文个分解反应的过渡态。5.3.2.4构型4的分解机理:4 (D3d) - TS9 - N6 (C2h) + 3N2在B3LYP/ 6-31 G*水平上我们对十二元环异构体4的分解途径进行了研究,找到了过渡态TS9,如图5.8所示。除了N2-N11. N1-N10键显著增长到1.594A, 1.607A外,总的来说,TS,的几何与稳定构型4的几何很相似。IRC复核计算证实经过TS9,构型4中N1-N10, N2-N11, N3-N12和N6-N9键同时断裂放出3个N2分子,同时生成一个线型N6分子。这一分解反应IL乎不需要活化能,在B3LYP/6-311+G(3df,2p)水平上,活化势垒仅仅为0.06 kcal/mol。如此低的势垒表明环状构型4是不可能存在的。1.213、君11.3787 .247飞 1 IN1.211一刀}}1、.594’卜’ 1.607工 .‘”TS9 (Cl)图过渡态TS,的优化几何5.3.2.5笼状构型5和6的分解机理:5 (D3d) }TS10}N8 (Dan)+2N26(D6n)}TSII-4 (D3d)}TS9}N6(C28)+3N2从结构特点分析,我们认为由三元环、四元环、五元环或六元环组成 的笼状分子稳定性应该比较差,更易分解。在B3LYP/6-31 G*水平上我们找到了笼型结构5的过渡态TS10。如图5.9所示,除了某些键的键长有所增长之外,TS10的几何非常接近反应物5的几何,暗示这个分解反应的活化能不会很高。IRC解析表明,5经过过渡态TS10分解放出两个N2分子,同时生成由两个五元环稠合而成稳定的具有Den对称性的N8分子。与我们的最初猜测相同,这个分解反应的势垒非常低,在B3LYP/6-311+G(3df,2p)//B3LYP+ZPE(B3LYP/6-31G*)计算水平上仅仅为2.5 kcal/mole67 北京理工大学博十学位论文1.537-.、1.551户尹沪口TS11 (C2h)TS10 (CO图5.9过渡态TS10和TS11的优化几何Ql l等[28]通过对N6(D6h)-N12(D6h)-N6(D6h)和Ni2(D6h)- 6N2的能级相关图分析指出,不论由N12(D6h)分子分解成二个N。分子或六个N2分子都是基态下对称性禁阻的反应,表明Ni2(D6h)分子是具有一定稳定性的结构。但我们对六棱柱型的异构体6的分解途径研究表明,它的分解首先经过一个异构化反应,即经过过渡态TS11转化为构型4再分解。如图5.9所示,与稳定构型6相比,TS11的N-N键长和键角都有些轻微的改变并且对称性从D6h变为C2h。同时我们注意到这个异构化反应几乎不需要活化能,在¥3LYP/6-31G*水平下(不包括零点能修正)异构化势垒仅仅为0.2kcal/mot,当进行零点能修正或在更高的计算水平下我们发现这个势垒成为一个负值,因此,六棱柱样的笼状构型‘可能并不存在。北京理工大学博士学位论文Relative energy aF1少C)1m0 刘咖02期巧均心明︐创2翻引幼︒比心油禅叫心呱1-咖J肋J加刁孤小撕‘ 胜 份r卜L甲受朴卜介f.卜r..片蛋.‘r..卜r炙解1.叶r.L份‘IL司肠r介.性.人而...L介r通‘护百L叮r.一6迎玛TS14、 ‘TS9么一介:----,··--一、 ‘亏一.‘二__一,二一、广一“’一”“‘’-一、TS5 、 、 亏、.,r峨汀龟- J一.又N6+3N23N2+2N3图5.10在B3LYP/6-31 G*水平上计算的N12分解反应的势能面示意图54结论我们用B3LYP/ 6-31 G*方法对四个以前研究过的Nit异构体((1, 3, 5和6)和两种新的Ni 2异构体((2和4)的分解和异构化机理进行了理论研究,得到了如下结论:(1) N-N三键连的笼状异构体((5和‘)由于环张力和氮原子上孤对电子的相互作用,使它们具有比线型和面型异构体高得多的能量,十二元环构型4的能量也较高。动力学数据显示环和笼型异构体的异构化或分解势垒极低,稳定性差,不是有实际应用价值的亚稳态分子。(2)含有芳香N5环结构的分子(1和2)比直线型分子(3)的能量低,但其分解时却发生五元环中N-N优先断裂,放出一个N:分子。相对而言,直线链状构型具有较高的动力学稳定性。但要在实验中被观察到,还需要特殊的条件。北京理工大学博士学位论文今考文献1. Christe K O, Wilson, W. W.; Sheehy, J. A.; Boatz, J. A. Angew. Chem. Int. Ed. Engl.1999,38:2004 2. See, for example: Chem. Eng. News 1999, 77, Jan 25, 7; 1999, 77, Nov 29, 38; 2000, 78,Aug 14, 413. Vij, A.; Wilson, W. W.; Vij, V.; Tham, F. S.; Sheehy, J. A.; Christe, K. O. J. Am. Chem.Soc. 2001, 123: 63084. Cacace, F.; Petris, G; Troiani, A. Science 2002, 295: 4805. Hammerl, A.; Klap6tke, T. M. Inorg. Chem. 2002,41: 9066. Pyykk6, P.; Runeberg, J. J. Mol. Struct. (THEOCHEM) 1991, 234: 2797. Lee, T. J.; Rice, J. E. J. Chem. Phys. 1991, 94: 12158. Glukhovtsev, M. N.; Jiao, H.; Schleyer, P. v. R. Inorg. Chem. 1996, 35: 71249. Bartlett, R. J. Chem Indust. 2000,4: 140IO.Lauderdale, W. J.; Stanton, R. F.; Bartlet, R. J. J. Phys. Chem. 1992, 96: 117311.Gagliardi, L.; Evangelisti, S.; Barone, V.; Roos, B. O. Chem. Phys. Let. 2000, 320: 51812.Nguyen, M. T.; Ha, T -K. Chem. Phys. Lett. 2001, 335: 31113. Wilson, K. J.; Perera, S. A.; Bartlet, R. J. J. Phys. Chem. A 2001, 105: 769314.Schmidt, M. W.; Gordon, M. S. Boatz, J. A. Int. J. Quantum. Chem. 2000, 76: 43415.Gagliardi, L.; Evangelisti, S.; Bemhardsson, A.; Lindh, R.; Roos, B. O. Int. J. Quantum.Chem. 2000, 277: 31116.Gagliardi, L.; Orlandi, G; Evangelisti, S;Roos, B. O. J. Chem, Phys. 2001, 114: 1073317.Chung,G;Schmidt, M. W.; Gordon, M. S. J. Phys. Chem. A 2000, 104: 564718.Ha, T -K;Suleimenov, O.; Nguyen, M.工Chem. Phys. Left. 1999,315: 32719.Korkin, A;Balkova, A.; Bartlet, R. J.; Boyd, R. J.; Schleyer, P. v. R. J. Phys. Chem.1996,100:5702 20.Dunn, K.; Morokuma, K. J. Chem. Phys. 1995, 102: 490421.Tobita, M.; Batrlett, R. J. J. Phys. Chem. 2001, 105: 410722.Gimarc, B. M.; M, Zhao. Inorg. Chem. 1996,35: 328923.Xin, W.; H-R, Hu;An-M, Tian.; Wong, N. B.; Siu-H, Chien.; W-K, Li. Chem. Phys.北京理工大学博士学位论文Let t. 2000, 329: 48324.Olah.GA.; Prakash, G K. S.; Rasul, G. J. Am. Chem. Soc. 2001,123: 330825.Bliznyuk, A. A.; Shen, M;Schaefer, H. F. Chem. Phys. Let. 1992,198: 24926.An-M, Tian.; F-J. Ding.; L-F, Zhang.; Y -M, Xie.; Schaefer, H. F. J. Phys. Chem. 1997,101:1946 27.Q-Shu, Li.; L-J. Wang. J. Phys. Chem. A 2001, 105: 197928.Qu. H.; Q-Shu, Li.; H-Sun, Zhu. Chin. Sci. Bull. 1997, 42: 46229.Manaa, M. R. Chem. Phys. Let. 2000, 331: 26230.Frisch, M. J;et al., GAUSSIAN 98, Revision A.9. Gaussian, Inc. Pitsburgh, PA, 199831.Becke, A. D. J. Chem. Phys. 1993, 98: 564832.Lee, C;Yang, W.; Parr, R. G Phys. Rev. B 1988, 37: 785内j几︐.Gonzalez, C.; Schlegel, H. B. J. Chem. Phys. 1989, 90: 2154月jd.Gonzalez, C.; Schlegel, H. B. J. Phys. Chem. 1990, 94: 552335.Klap8tke, T. M.; Harcourt, R. D. J. Mol. Struct. (THEOCHEM) 2001, 541: 237北京理工大学博士学位论文第六章MN5(M=Li, Na, K, Rb)动力学稳定性的研究6.1引言 聚氮和富氮化合物一直是理论和实验化学家关注的焦点[1-101。除了结构的重要性外,更是由于它们用作环境友好高能量密度材料的潜在意义。然而,由于合成过程的高放热性,它们的合成对实验化学家一直是一个巨大的挑战。N犷盐[11]的成功合成,极大鼓舞了人们对合成其它稳定的多氮化合物的信心。最近,Cacace等人t12]用单电子还原气态N4+的中性化一再离子化质谱,检测到气态亚稳态的N4分子在298K下有超过1微秒的寿命。2002年7月,Vij等人[[131通过断裂苯基五哇衍生物的C-N键,第一次实验报道了长期寻求的cyclo-N5一的存在。这些新的发现为实验合成其它更多的聚氮化合物打开了光明的前景。N5一阴离子的成功发现为实验合成它的稳定的盐提供了极大可能。对 于取代的N5一已有一些理论的研究。Bartlett等[14]用偶合簇(CC)和多体微扰理论(MBPT)研究T HN,的结构和稳定性。Schleyer等人[151采用abinitio方法对N5Li进行3研究。Nguyen等人[16]对TM(COW r1 5-N5)仃M=Fe2+, Mn,和Cr-)完成了Mickel理论计算(EHF),结果表明这些络合物能够稳定存在。最近,Gagliardi和即ykk6[17-18]理论预测了两类新的化合物的存在,Sc( n 7-N7)和n 5-N5--M- IT7- N73-。Lein等人【191预测二茂铁样结构的Fe( 15- N5)是一种很稳定的络合物。Burke等[[201采用从头算和密度泛函的方法研究了MN5(M=Na, K, Mg, Ca,和Zn)的性质,结果显示金属-N5环络合物之间静电相互作用占主要地位,并且他们也指出,金属阳离子对N犷环的电子分布、振动和NMR光谱都产生了一定的影响。然而,他们没有研究这些化合物的动力学稳定性。为了更进一步判断这些化合物的动力学稳定性,我们对MNs(M= Li, Na, K, Rb)的异构化和分解途径进行了理论研究,目的是为Ns_金属衍生物的实验合成提供一些理论的支持。北京理工大学博士学位论文6.2计算方法全部计算都利用Gaus sian 98程序完成〔21]a采用密度泛函理论中的B3LYP[22]和B3PW91 [23-24]方法计算了各种异构体的结构、能量和振动频率。6-311+G*基组用于氮原子和金属原子Li, Na, K,对于第四周期重金属元素Rb,使用的是有效核势基组(Eeffective Core Potential,ECP )[25-27].振动频率的计算用于表征稳定点(虚频个数=0)或过渡态(虚频个数=1)。用内IK反应坐标((IRC)[28-29]计算对过渡态所连接反应物和生成物的最小能量路径进行了复核。为使势能面上的能量更精确,我们还用QCISD(T)/6-311+G*方法在B3LYP/6-31 I +G*优化出的几何的基础上进行了单点能的校正.另外,为了获得更多的的电子结构信息,我们在B3LYP/6-311+G*理论水平下进行T NBO分析〔30, 31]fl 我们用POLYRATE 8.2程序[32]计算了不同温度下的反应速率常数。为了比较,分别用传统过渡态理论(TST)和正则变分过渡态理论(CVT)计算了反应的速率常数【33-34]。还基于离心为主的小曲率半经典绝热基态(CD-SCSAG)方法,用小曲率隧道发射系数对CVT进行了校正。63结果和讨论从结构的观点看,金属原子Li, Na, K和Rb与N,环如下图所示的三种方式成键:金字塔结构A,平面结构B and C。︐了矛NA B C (X=Li, Na, K和Rb)图6. 1-6.4分别展示了LiN5, NaNS, KNs,和RbN,的稳定构型、过渡态和离解产物的几何结构和计算的参数。全部能量、零点能((ZPE)和最低振动频率分别列在表6.1, 6.3, 6.5和6.7中。表6.2, 6.4, 6.6和6.8列出了相对能量。北京理工大学博士学位论文6.3.1 LiN。的结构和动力学稳定性与( CH)5Li相似,具有C5,对称性的金字塔结构LiNs-A是势能面上的极小点口但它的能量比平面结构LiNs-B(C2,)高8.2 kcal/mol。与NS环以单键相连的平面结构LiNs-TS-C不是势能面的极小点。在B3LYP和B3PW91水平上完成的IRC计算证实它是金属Li沿着N5环迁移的过渡态。在QCISD(T)水平上,LiNs -TS-C的能量比LiNs-B高4.8 kcal/mole从图6.1可以看到,LiNs-A中N-N键长相等并且接近孤立的N5阴离子的N-N键长(在B3LYP水平下为1.326.),表明在LiN,一^络合物中N5阴离子的几何结构受金属的影响很小。然而,在构型LiNs-B和LiNs-TS-C中,金属阳离子的存在使NS一环中N-N键长发生了变化,使N5环芳香性降低。我们进一步对上面这三种构型的异构化和分解反应机理进行了研究。金字塔构型LiNs-A经过过渡态TSl可重排到能量最低的异构体LiNs-B,在QCISD(T)水平上,这个异构化反应的势垒仅为2.0 kcal/mol,表明金字塔构型LiNs(C5v)的动力学稳定性很低。另一方面,我们也在DFT水平上对异构体LiNs-B的单分子分解反应机理进行了研究。IRC计算表明LiNs-B” TS2(Cs) - LiN3(C m v)+N:的分解反应经过线性络合物LiNs-complex(C-,)完成。LiNs-complex(见图6.1)具有全正的振动频率,是势能面上的一个极小点。能量上,LiNs-complex比B低35.6 kcal/mole基于LiNs-B和TS2的能量差,LiNs-B的分解垒高预测是19.9 kcal/mol(表6.5)。以前的研究【35]表明弥散函数可以充分地减少BSSE影响,并且我们的计算也显示BSSE修正的影响相对较小(小于键能的5%)。因此,在我们目前的工作中,键能的计算忽略了BSSE。如表6.5所示,LiN5-complex的键能在QCISD(T)水平上预测是6.2 kcal/mol o Glukhovtsev等[15]指出线型LiN3是最稳定的异构体。如图6.1所示,在B3LYP/6-311+G*水平下,Li-N键长是1.7o9Aa NBO分析表明,Li-N的Wiberg键序(WBI)接近于零,而且,正电荷主要在Li原子上。因此,LiN3本质上是离子型化合物且静电相互作用占优势。北京理工大学博士学位论文B3LYP/6-311+G*886 B3 P W91 /6-31 I+G"1.305LiNs-B (C2v)LiNS-A (C5.)Li-TS-C (C2.)1.219TSl (C¥)1.1501.1971.70193LiN3TS2 (Cs)2.149281.094LiNs-Complex图6.1优化的LiNS的稳定构型和过渡态的几何结构75北京理工大学博士学位论文表6.1 LiN5全部能量[E(hartrees)],零点能[ZPE(kcal/mol)]和最低振动频率〔。!(cm')]B3LYP/6-311+G*B3PW91/6-311+0*QCISD(T)/6-Species311+G* aE ZPE u,ZPE ofLi Ns-A-281.3101914.6220-281.1930614.8 220一80.61544TS1 -281.3068814.3-179-281.18977814.6 -180一80.61183UNs- B-281.3272815.0158-281.2096115.3 160-280.62916Li N5-TS-C-281.31864147-72-281.2012915.0 -72-280.62096TS2 -281.2969913.2-588一81.1731813.3 -601-280.59453LiNs-Complex-2813829812.822一81.2521712.8 25-280.68247Li N3(C"v)-171.813098名159-171.733378.9 140一171.35645UN尹Na 一81.3727912.3-281.2433912.4一80.67169.单点能计算在QCISD(T)/6-31 1+G*/IB3LYP/6-311+G.水平上。表6. 2 LiNs异构体的相对能量(kcal/mol)SpeciesB3LYP/6-311+G* B3PW91/6-311+G*QCISD(TX6-311+0*aUN. ,-A 10.39.982TS1 12.111.810.2Li Ns习0OO00.0Li 凡-TS-C5.14.94.8TS2 17.220.919.9LiNg -Complex-37.1一9.2-35.6Li N3+Nx一13一41-29.4a相对能量在QCISD(T)/6-31 1+G*1/B3LYP16-3I I+G*+ZPE(B3LYPl6-311+G*)水平上。6. 3. 2 NaN,和KN。的结构和动力学稳定性 与LiN5情况相似,在B3LYP和B3PW91水平上,具有C2,对称性的NaN5-B和KNS-B·分别是他们各自势能面上能量最低的异构体。在QCISD(T)计算水平下,金字塔结构分别比平面构型B高5.7和晰kcal/mot,如表6.4, 6.6所示。TS3(Cs)和T55(Cs)是找到的由A异构化到B的过渡态,异构化垒高在QCISD(T)水平上分别为0.7和2.7kcal/mol,如此低的垒高表明金字塔构型NaN3-A和KN¥-A似乎都是动力学不稳定的,他们将很快异构化为NaNs-B和KNS-B。与LiNs相似,NoNs-TS-C76 北京理工大学博十学位论文和KN5-TS-C也都是对应于金属Na和K沿N3环迁移的过渡态。在QCISD(T)水平下,T S-C和B的能量差为4.0和5.1kcal/mot。我们也对能量最低异构体NaNs-B和KN5-B的分解机理进行了调查,找到了两个具有CS对称性的过渡态TS4和TS6,如图6.2, 6.3所示。分解反应NaNs-B-TS4-NaN5-complex-YNaN3(C一,)+N2和KN5-B一TS6- KN5-complex-KN3(C-,)+N2的垒高在QCISD(T)水平上预测是22.0和22.5 kcal/mol,显示出比LiN5-B高的动力学稳定性。另外由表6.4, 6.6也可以看到,NaN5-complex和KN5-complex的键能分别是4.1和3.2 kcal/mol(不考虑BSSE修正)。表6.3 NaN,全部能量[E(hartrees)]>零点能[ZPE(kcal/mol)]和最低振动频率[o I (cm I)lB3LYP/6-31 l+G*B3PW9 1/6-311+G*QCISD(T)/6-Species311+G*'EZP已EZPE ulNaNs魂 一36.0854414.0 100 -435.9334314.3t00一35.01437TS3 -436.0842913.9一2 -435.9322014.2一2-435.01311NaNs -B436.0976514.4 101 -435.9449814.7101-435.02414NaN5一5- C-436.0903414.1 -80 -435.9378114.4-80-435.01723TS4 436.0627712.4 -626 -435.9038712.6-629434.98595NaNs-Complex礴36.14128II名12 -435.9763411乡13一35.06606NaNAC动 -326.575017.9 82 -326.460778.093-325.74366NaN3 斗NJ-436.1347011.4一-435.9707811.5-435.05889,单点能计算在QCISD(T)/6-311十G*/B3LYP/6-311+6*水平上。表6. 4 NaN,异构体的相对能量(kcal/mol)SpeciesB3LYP/6-311+G* B3PW9116-31 1+G*QCISD(T)/6-311+G*'NaNs-A7.36吕5.7TS3 7.97.56.4NaNsB 0.00.000NaNs -TS-C434.24.0TS4 19.923.722.0NaNs-Complex-30.0一2.5-28.9N皿Nl+N, -26.3-19.4-24.8a相对能量在QCISD(f)/6-311+G*J/B3LYP/6-31 i+G*+2PE(B3LYP/6-3I I +G*)水平上。77北京理工大学博士学位论文表6.5 KN,全部能量[E(hartrees)],零点能[ZPE(kcaUmol)]和最低振动频率[u 1 (cm')]B3LYP/6-311十G*B3PW91/6-311+G*QCISD(T)/6-Species311+G*a E ZPE u[ZPE uiKNS ,A-873.7429714.0108-873.57777143113-872.51979TS5 -873.7401913.9-81-873.5748514.2-82-872.51539KNs一 B-873.7489014.281-873.5834014.580-872.52268KNs -TS-C-873.7409714.079 -873.5755714.2-79-872.51417TS6 -873.7130312.2-626-873.5412812.4-630-872.48367KNS-Complex-873.7872311.613-873.6103311.715-872.55985KN3 (C--764.224407.679-764.097467.774-763.23874KN3 +N2-873.7841011.1-873.6074711.2-872.55398单点能计算在QCISD(T)/6-31 I+G*/B3LYP/6-3I I +G*水平上。表6.6 KNS异构体的相对能量(kcal/mol)SpeciesB3LYP/6-311十G* B3PW91/6-311+G*QCISD(T)/6-311+G*KNs-A3.53.31. 6TS5 5.25.143KNs -B0.00.00.0KN5 -TS-C4.84.65.1TS6 20.524.322.5KNs-Complex-26.7-19.7-25.9KN3 +NZ-25.2-18.4-22.7a相对能量在QCISD(T)/6-31 1+G*/B3LYP/6-311+G*+ZPE(B3LYP/6-311+G*)水平上。78北京理工大学博士学位论文’︐’︐︒“口已口从口绍l曰︑︑{宝翌NaN犷B(CZv)NaNrA(Csv)11111.327NaN宁TS一C(C知)322}黑1·戮TS3(Cs)进翌︑︑}粼全}琵︑︑︐︑︑}蹂︑掇了嘟︑︑︑︑︑{黑里且子盆NaN3TS4(Cs)1. 1561」952刀76垂翌鹅些场一里塑一铆圣_撰1. 0950上塑曝NaN犷ComPIel图6.2优化的NaNS的稳定构型和过渡态的几何结构79北京理工大学博士学位论文1.315B3LYP/6-31 l+G"281.323KN-5-B (C2v)KN5-A (C5.)108.231.319KN5-TS-C (C2v)TS5 (C户1.1591.234豁1564.1.30, ,, 3.033砌1.097 '3.102 1.095TS6 (Cs)KN5-Complex1_.15751972.380KN3 图6.3优化的KN5的稳定构型和过渡态的几何结构80北京理工大学博士学位论文1316RbNs-B (C20RbNs-A (C5,,)圣6_.4犯08.223RbNS-TS-C (C2v)108.0.1219T57 (Cs)1.162 1.193153.41,3.护2砂,72沪3.2691.176TSS (Cs)1_098RbN5-Complex1.1601.158:一{9522.59846RbN3图6.4优化的RbN,的稳定构型和过渡态的几何结构81北京理工大学博士学位论文6. 3. 3 RbN。的结构和动力学稳定性鉴于从轻金属都重金属(由Li到K) MN5的稳定性逐渐加强,我们也对RbNs的可能结构和稳定性进行了研究。就我们所知,没有重金属与6n电子芳香环之间的键的研究被报道。如图6.4和表6.8所示,以两点与环相连的平面构型RbNs-B和金字塔构型RbNs-A都是势能面上的局域极小点。在QCISD(T)水平上,RbNs-B的能量LL RbNs-A低1.3 kcal/mol eRbNs-TS-C也是Rb原子沿N5环迁移的过渡态,在QCISD(T)水平下,其异构化势垒是4.2 kcal/mol。与NaNs和KN 5相似,我们也研究了RbNs-A与RbNs-B之间的异构化和RbNs-B的分解反应。从表6.8可以看到,这个很低的异构化垒高(在QCISD(T)水平下为3.1 kcal/mol)表明金字塔型的异构体RbNs-A动力学也是不稳定的。在同样水平下,RbNs-B的分解势垒高预测为23.0 kcal/mol,表现出比NaN5和KN,稍高的动力学稳定性.NBO分析表明,金属原子的电正性越大,对阴离子N5一环的电子结构影响越小。QISD(T)/6-311+G*//B3LYP/6-311+G*+ZPE(B3LYP/6-311+G*)水平上,环状N5一 (Dsn)的分解势垒是25.8 kcal/mol o表6.7 RbN。全部能量[E(hartrees)],零点能[ZPE(kcal/mol)]和最低振动频率[o I (cm 1)]B3LYP /6-311+6*'B3PW91/6-311+G*'QCISDSpecies311+G*'bE ZPE v,ZPE.1RbN,A -297.9087113.987-297.5878314.189一97.00782TS7 -297.9065713.8-63-297.8187214.1-64-297.00474RbN5 -B-297.9130514.168一97.8250114.466-297.01002RbN5 -TS-C-297.9060813.8-72一97名180611.1-71-297.00287TS8 -297.8763312刀-626一97.7820212.2-630-296.96998RbNS-Complea-297.9482911.658-297.8482411.755-297.04459RbN3 (C刁-188.385187.454-188.335747554-187.72430RbN2 +N2-297.9448810.9-297.8457511.0-297.03953对N原子使用是6-311+G*,而对重金属Rb使用的是ECP基组。单点能计算在QCISD(T)/6-31 1+G*/B3LYP/6-311+G*水平上。82北京理工大学博士学位论文表6.8 RbN5异构体的相对能量(kcal/mol)SpeciesB3LYP/6-311+G*B3P W91 /6-31 l+G*QCISD(T)/6-311+G* aRbNs-A2.11. 9 1.3TS? 3.4333.1RbNs -B0.00.0 0.0RbNs -TS-C4.04.0 4.2TS8 21.024.8 23.0RbN5 -Complex一4.6-17.3 -24.2RbN3 +Nz一2.7-15.9 -21.1对N原子使用是6-311+G*,而对重金属Rb使用的是ECP基组。相对能量在QCISD(T)/6-311 +G*//B3 LYP/6-31 l+G*+ZPE(B3LYP/6-311+G*)水平上。6.3.4分解反应.Ns一 B--MN,+Nz的速率常数利用POLYRATE 8.2程序,我们计算了分解反应MN5-B- MN3十N2(M=Li, Na, K,和Rb)的速率常数。采用TST, CVT和CVT/SCT方法计算的速率常数值分别列于表6.9-6.12。从表中可以看到,TST和CVT值非常接近,表明变分的影响很小,而CVT/SCT方法计算的速率常数值比TST和CVT值都稍大,表明小曲率隧道修正对这种分解反应是必要的。从上面的讨论可以看出,CVT/SCT理论速率常数似乎更可靠。表6.9LiNs-B-TS2--LiN3+N2的分解反应的速率常数(cm 3 mol-I S-1)T (K)TSTCVTCVT/SCT250. 01.5729E-021.4047E-021.7099E-01298. 05.0754〔十004.2282E+002.3679E十01350. 04.6188E+023.6062E+021.2313E+03400. 01.2012E+048.9040E+032.2434E十04450. 01.5413E十051.0942E+052.2379E+05500. 01.2014E+068.2259E+051.4489E+06600. 02.6678E+071.7266E十072.5035E+07700. 02.4793E十081.5396E+081.9868E+08800. 01.3297E+098.0023E+089.5847E十081000.01.4102E+101.3610E十101.4978E+101250.09.3938E+108.7236E+109.1422E+101500.03.3358E+113.0188E+I13.0867E+1183北京理工大学博士学位论文表6.10NaNs-B--TS4--NaN3+N2的分解反应的速率常数(cm 3 mol"I S")T (K)TSTCVTCVT/ SCT250. 05.6030E-055.3578E-053.9320E-03298. 04.5213E-024.3563E-021.0328E+00350. 08.3911E+008.1337E+009.2528E+0 l400. 03.6427E十023.5466E+022.5125E+03450. 06.9700E十036.8100E+033.4265E+04500. 07.4834E+047.3325E+042.8511 E十05600. 02.6893E十062.6463E+067.1862E十06700. 03.5267E+073.4808E+077.4810E十07800. 02.4495E+082.4229E十084.4322E+081000. 03.7313E+093.7018E十095.5354E+091250. 03.3211 E十103.3021E+104.3037E+101500. 01.4307E+111.4244E十111.7173E+11表6.11KN5-B-"TS6--KN3+N:的分解反应的速率常数((cm3 m01'1 S'I)T (K)TSTCVTCVT/ SCT250. 01.5128E-053.8066E-061.3014E-02298. 01.4822E-023.6230E-033.2355E-02350. 03.1979E十007.6445E-011.0884E+00400. 01.5468E+023.6415E+013.9523E十01450. 03.2197E十037.4962E+027.0945E十02500. 03.6984E+048.5404E十037.3959E+03600. 01.4710E+063.3611E十052.6244E+05700. 02.0745E十074.7123E+063.4956E+06800. 01.5218E+083.4459E+072.4944E+071000. 02.5030E+095.6562E十084.0502E+081250. 02.3693E+105.3599E十093.8925E十091500. 01.0635E+112.4121 E十101.7912E+1084北京理工大学博士学位论文表6.12RbNs-B-TS8-RbN3十NZ的分解反应的速率常数((cm 3 mol"I S-1)T (K)TSTCVTC VT/SCT250. 05.4546E-065.4210E-063.1626E+00298. 06.3493E一36.3235E-035.6253E+00350. 01.5658E+001.5618E+002.2702E+01400. 08.3358E+018.3236E+013.9822E+02450. 01.8698E+031.8684E+036.1754E+03500. 02.2803E+042.2795E+046.0056E+04600. 09.9224E+059.9223E十051.9599E+06700. 01.4921E+071.4919E+072.4709E+07800. 01.1486E+081.1479E斗081.6915E+081000. 02.0209E+092.0172E+092.5806E+091250. 02.0189E+102.0113E十102.3425E十101500. 09.3936E十108.2784E+105.4935E+106.4结论我们采用密度泛函理论中的B3LYP和B3PW91方法系统研究了N5 -阴离子环与碱金属(Li, Na, K和Rb)络合而成的一系列化合物的结构性质和分解途径,得到如下结论:(1)与己知的金字塔结构的(CH)5Li相似,金字塔结构的LiN5-A,NaN5 -A, KN5-A和RbN5-A也都是它们各自势能面的极小点。然而,非常低的异构化势垒表明这种金字塔结构的异构体并不稳定,它们 将很快重排为最稳定的平面构型B. (2)从轻金属都重金属元素(由Li到Rb), MN5的分解稳定性逐渐加强。这可能是由于金属原子的电正性越大,对阴离子N5 一环的电子结构影响越小的缘故。 (3)我们也采用TST, CVT和CVT/SCT方法计算T分解反应MN5-B--MN3 十N2的速率常数,计算结果显示,对这种类型的分解反应来说,小曲率隧道修正方法计算的速率常数似乎更可靠。 85北京理工大学博士学位论文今考文献1. Pyykko P, Runeberg J, J. Mol. Struct. (THEOCHEM) 1991, 234: 2792. Yarkony D R, J. Am. Chem. Soc. 1992,114: 54063. Glukhovtsev M N, Jiao H, SchleyerP v R, Inorg. Chem. 1996, 35, 71244. BartlettR J, Chem Indust. 2000, 4:1405. Nguyen M T, Ha T -K, Chem. Phys. Lett. 2001, 335: 3116. Workentin M S, Wagner B D, Lusztyk J, Wayner D M, J. Am. Chem. Soc. 1995. 117:119 7. Gagliardi L, OrlandiqEvangelisti S, Roos B O J, Chem. Phys. 2001, 114: 107338. Michels H H, Montgomery JA Jr, Christe K O, Dixon D A, J. Phys. Chem. 1995, 99187 9. Hammerl A, Klapotke T M, Inorg. Chem. 2002, 41: 90610. Olah G A, Prakash G KS , RasulqJ. Am. Chem. Soc. 2001,123:330811. Christe K O, Wilson W W, Sheehy J A, Boatz J A, Angew. Chem. Int. Ed. 1999, 38:2004 12. Cacace F, PetrisQTroiani A, Science. 2002, 295: 48013. Vij A, Pavlovich J Q Wilson W W, Vij V, Christe K O, Angew. Chem. Int.Ed. 2002,41:3051 14. Ferris K F, BartlettR J, J. Am. Chem. Soc. 1992, 114: 830215. Glukhovtsev M N, Schleyer P v R, Maerker C J. Phys. Chem. 1993, 97: 820016. Nguyen M T, McGinn M A, Hegarty A F, Elguero J, Polyhedron. 1985, 4: 172117. Gagliardi L, Pyykk6 P, J Am. Chem. Soc. 2001, 123: 970018. Gagliardi L, Pyykko P, J. Phys. Chem. A 2002, 106: 469019. Lein M, Frunzke J, Timoshkin A, FrenkingqChem. Eur. J. 2001, 7: 415520. Burke L A, Butler R N, Stephens J C J. Chem. Soc., Perkin Trans. 2001, 2: 167921. Frisch M J, et al., GAUSSIAN 98, Revision A.9. Gaussian, Inc., Pitsburgh, PA, 199822. Becke A D, J. Chem. Phys. 1993, 98: 564823. Lee C, Yang W, Parr RG Phys. Rev. B 1988, 37: 78524. Perdew J, Wang P Y, Phys. Rev. B 1996, 54, 16533北京理工大学博士学位论文25. Andrae D, Haeussermann U, Dolg M, Stoll H, Preuss H, Theor. Chim. Acta. 1990, 77:123 26. Nicklsas A, Dolg M, Stoll H, Preuss H, J. Chem. Phys. 1995, 102: 894227. Leininger T, Nicklsas A, Stoll H, Dolg M, Schwerdtfeger P, J. Chem. Phys. 1996, 105:1052 28. Gonzalez C, Schlegel H B, J. Chem. Phys. 1989, 90: 215429. Gonzalez C, Schlegel H B, J. Phys. Chem. 1990, 94: 552330. Carpenter J E, Weinhold F, J. Mol. Struct. (THEOCHEM) 1988, 169: 4131. Carpenter J E, Ph. D. Thesis, University of Wisconsin, Madison. 198732. Stickler R, HuW一P, et al. Comput. Phys. Commum 1995, 88: 34133. Lide D R, (Ed) CRC Handbook of Chemistyr and Physics; CRC perss, Boca Raton,Fl oroda, 73TH, 1992-1993; p.9-18; p.9-3134. King C M, Nixon E R, J. Phys. Chem. 1968, 48: 168535. Feller D, J. Chem. Phys. 1992, 96: 6104北京理工大学博士学位论文致 谢谨向我的导师李前树教授致以诚挚的感谢和崇高的敬意。本论文是 在他的悉心指导和严格要求下完成的。他在科研中开拓进取、孜孜不倦的精神和他敏捷的思维、渊博学识都给我留下了深刻的印象,他传授给我的科学知识和研究方法将使我终生受益。感谢乔治亚大学的谢尧明博士、北京理工大学的徐文国教授、张锡 宏教授和张绍文博士,与他们在学术方面的讨论使我获益非浅;感谢北京理工大学的严沛然教授给予的帮助;同研究小组的高凤新老师、殷鹏刚、吕瑞花、关君、程丽萍、吕立强、李亚萍、张越、金俏同学以及已经毕业的李国良博士、陆军博士、李欢军博士、朱锐博士、徐卫国博士、王朝阳博士、张星辰博士、武振羽博士、金宏威博士、刘永东博士在共同学习和生活中也给予我许多帮助,在此表示感谢。感谢我的父母、爱人和远在美国学习的弟弟,他们的艰苦奋斗精神 是我的精神支柱,他们始终如一的支持和鼓励是我学业和研究工作得以顺利完成的保障,还要感谢我的女儿,她可爱的笑声为我的学习生活增添了色彩和动力。感谢所有曾经和正在给我帮助的人们。 北京理工大学博十学位论文攻读博士期间所撰写的论文1. Qian Shu Li, Jun士ants Z红ao. "Theoretical Study of Potential Energy Surfacesf or Nit Clusters" Joumal of Physical ChemistyrA 2002, 106, 5367-5372.2. Qian Shu Li, Jun Fang Zhao. "A Theoretical Study on Decomposition Pathwaysof N7' and Ni Clusters" Joumal of Physical ChemistyrA 2002, 106, 5928-5931.3. Jun士ang Zhao, Qian Shu Li. "Structures and Kinetic Stability of N7 Cluster"Chemi cal Physics Letter 2003, 36811-2, 12-194.迎匹旦迎眨塑丝,Nan Li, Qian Shu Li. "A Kinetic Stability Study of MN5 (M=Li, Na,K, Rb)" Theoretical ChemistyrAccounts(印刷中)·5. Yong Dong Liu, Jun Fang Zhao, Qian Shu Li. "Structures and stability of N7'and Ni clusters" Theoretical Chemistry Accounts 2002, 107, 140-146.6. Qian Shu Li, Jun- Fang Zhao, Yaoming Xie, Schaefer, H. F. "Structures,Ther mochemistyr and Electron Affinities of the CX3 (X=H, CI, F, Br)" MolecularPhysi cs (invited article) 2002, 100, 3615.7.壑鱼羞,李楠等,“氮笼N,:的量子化学计算”淳等学友钻学学蔺2002, 23,1944-1947. 8. Jun Fang Zhao, Qian Shu Li. "A Kinetic Stability Study of MN5' (M=Be, Mg, Ca,Sr , and Ba)", lntemational Journal of Quantum Chemistyr(审理中).9.壑煌羞,刘永东等,N12原子簇结构和稳定性研究第八次全国量子化学学术会议论文集长春2002, 7089北京理工大学博士学位论文作者简介1969年1 0月出生于河北省石家庄市。1987年考入河北轻化工学院(现河北科技大学)化T-系有机化工专业,1991年毕业并获得工学学士学位,同年分配到河北省石油化工规划设计院任助理工程师。1995年调入石家庄铁道学院基础部任教,1997年评为讲师。1996年考入石油大学(北京)化学工程系化学工艺专业,1999年获得工学硕士学位。2000年考入北京理工大学化工与材料学院(现属理学院),攻读应用化学溥士学位,指导老师为李前树教授。在读博士期间,参加了国家自然科学基金和国防科技预研基金项目的部分研究工作。