CaSin(n=1-10)团簇结构_稳定性与光谱性质的理论研究
- 格式:pdf
- 大小:668.13 KB
- 文档页数:5
团簇和纳米体系物理发展团簇和纳米体系是20世纪末发展起来的崭新领域,它所研究的对象是既不同于原子、分子,又不同于宏观物体的中间体系,尺度范围大约在0.1-100nm,这是人们过去从未进行研究的新领域,是人们认识物质世界的新层次。
它的丰富物理内涵,对物理学提出了新的挑战,也是当前物理学与其它学科交叉最富有活力的热点领域。
团簇简介原子和分子团簇,简称团簇(Cluster)或微团簇(microclusters),是几个乃至上千个原子、分子或离子通过物理或化学结合力组成相对稳定的微观和亚微观聚集体,其物理和化学性质随所包含的原子数目而变化。
团簇的空间尺度是几个至几百个纳米的范围,用分子描述显得太大,用小块固体描述又显得太小,许多性质既不同于单个原子分子,又不同于固体或液体,也不能用两者性质作简单线性外延和内插得到。
因此,人们把团簇看作是介于原子分子和宏观固体之间物质结构的新层次,有人称之为物质的“第五态”。
正因为如此,团簇可作为各种物质由原子分子向大块物质转变过程中的特殊物相,或者说它代表了凝聚态物质的初始状态,团簇的研究有助于我们认识大块凝聚物质的某些性质和规律。
团簇科学是研究团簇的原子组态和电子结构、物理和化学性质及其向大块物质演化过程中与尺寸的关联,团簇同外界环境的相互作用规律等。
团簇科学处于多学科交叉的范畴。
从原子分子物理、凝聚态物理、量子化学、表面科学、材料科学甚至核物理学引入的概念和方法交织在一起,构成当前团簇究的中心议题,并逐渐发展成一门介于原子分子物理和固体物理之间的新型学科。
团簇研究的基本问题是:弄清团簇如何由原子、分子一步步发展而成,以及随着这种发展、团簇的性质将如何变化,当尺寸多大时,团簇发展成宏观固体人们知道,由若干原子构成的分子,可在气相、液相和固相中稳定的单元存在,而团簇作为原子聚集体往往产生于非平衡条件,很难在平衡的气相中产生。
当团簇尺寸较小时,每增加一个原子,团簇的结构发生变化,即所谓重构。
团簇独特性质所谓团簇,就是由相距较远的小团簇合并成的大团簇。
团簇大小取决于尺度因子,可以用它们来衡量分子尺寸的大小。
在团簇化学中,团簇常表现为无限小,其原因是因为原子尺度的限制。
从原子的角度讲,分子或团簇都是一个整体,只有当这些粒子处于较大的空间范围时才有可能实现。
团簇是研究基本粒子结构与物理性质的一种重要方法,在原子和分子物理中得到了广泛应用,也是探索核反应过程中粒子的产生、衰变和反应机理的有力工具。
一般情况下,由几个或几十个成员组成的团簇叫做链状团簇,简称团簇。
最简单的一种链状团簇—— 1+1型团簇是由一个金原子与一个磷原子组成的,而最简单的12+1型团簇则由13个磷原子组成,因此12+1型团簇又被称为“核素”。
虽然目前已经发现了数百种的团簇,但是像“ 12+1型”、“ 13+1型”这样规整的形式还未出现过。
由于团簇有着丰富的特征,所以人们不断地提出各种假说,试图解释这些特殊形状。
4、团簇具有“体积效应”。
团簇之所以呈现球状,主要是由于其原子在外界影响下形成团簇后,由于受到附近电子云的影响,原子间的距离大大增加,相当于同周期原子体积之和。
这种微观结构上的改变,使得团簇对磁场有很强的敏感性。
团簇虽小,却对周围的环境非常敏感,表现出比晶体更为复杂的内部运动状态。
由于原子之间存在相互作用,因此通过研究其中每一个原子的运动轨迹可以准确地推断出团簇的结构。
5、团簇具有“吸附效应”。
团簇形成过程中会发生大量粒子间的碰撞,同时会不断向团簇周围辐射一定波长范围的光波。
这种现象称为光的“吸附效应”,即吸附光子效应。
团簇不仅可以吸收可见光,也可以吸收X射线和γ射线等。
许多团簇还能吸收紫外光。
因此,利用这一性质可以研究分子的振动性质。
6、团簇具有“量子效应”。
在团簇的形成过程中,由于电子轨道的不断重新排列,团簇的平均能量不断升高。
根据能级交错原理,高能态粒子数密度将超过低能态粒子数密度。
而且,随着团簇的稳定性逐渐增加,能级交错也越来越明显,从而形成明显的激发态和各种亚稳态,出现各种带色的特征光谱。
团簇构造、功能及多级演化重大研究计划2021年度项目指南团簇是介于原子/分子与宏观物质之间的多核聚集体,具有确定的原子组成和化学结构,代表了凝聚态物质的初生态,是关联宏观性质和物质微观结构的理想模型,对深刻认识和理解物质转化的规律具有重大意义。
一、科学目标通过化学、物理、生命、材料、环境、信息等多学科交叉,发展新型团簇及其多级结构构筑的新概念、新策略、新方法和新反应,建立团簇高精度和高分辨表征的新技术,在原子水平上揭示团簇特殊性质的结构基础与演变规律,理解团簇结构与功能的关联,制备功能团簇基材料与器件,解决基于团簇的变革性技术中的关键科学问题,促进相关学科的发展。
二、核心科学问题本重大研究计划将聚焦团簇构效关系,探索物质结构与性能随团簇尺寸变化的规律,揭示团簇稳定性机制,理解多级团簇体系中主体与环境的作用机制,实现功能导向的多级团簇结构的精准构筑和宏量制备。
(一)团簇的稳定性机制。
具有特殊结构与独特性能的新型团簇的发现、团簇形成机理和稳定化机制的理解、各种化学键及弱相互作用的认知。
(二)团簇电子结构的规律。
团簇结构及稳定性随团簇尺寸的演变规律、团簇的“幻数”特性、团簇的构效关系。
(三)多级团簇功能的调控原理。
多级团簇功能与团簇内聚集态、簇际相互作用、团簇与环境耦合的关系。
三、2021年度重点资助研究方向针对原子团簇、分子团簇及其多级体系(特别是生物体系)中的重要科学问题,发展团簇研究的新方法和新理论,揭示团簇的形成机制、稳定性规律和构效关系,构造具有独特功能的团簇材料与器件。
本重大研究计划2021年度重点资助以下研究方向:(一)团簇及其多级结构的精准构筑。
建立团簇的定向设计、高效合成策略和宏量制备方法,认识和理解团簇稳定性机制,揭示团簇制备的调控规律。
理解团簇内聚集态、簇际相互作用、团簇与环境及外场间相互作用的规律,揭示团簇基功能材料的构效关系,进而实现功能复合。
重点支持新型团簇体系的合成方法和多级组装结构的研究,特别注重模拟酶团簇、主族元素团簇、稀土和锕系元素团簇、中性团簇及团簇组装的功能材料体系。
CoBn(n=1~7)团簇的结构和磁性张文庆;张先富;胡志刚;侯丹峰【期刊名称】《磁性材料及器件》【年(卷),期】2009(040)003【摘要】利用密度泛函理论中的广义梯度近似对CoBn(n=1~7) 团簇的结构、稳定性和磁性进行了分析.综合CoBn团簇基态的平均结合能、二阶能量差分、垂直电离势、最高占据分子轨道与最低未占据分子轨道之间的能隙分析可知,n=3和7的三维团簇的稳定性表现突出.Mulliken布局分析表明体系的磁矩主要来自局域d 电子的贡献.CoBn团簇中Co原子的磁矩和总磁矩分别在1<n≤3和4≤n<7振荡变化,在n=3和7发生"猝灭".CoB3和CoB7团簇较高的稳定性和"磁矩猝灭"是由于较多的电荷转移和在费米能级强烈的spd杂化.【总页数】3页(P12-14)【作者】张文庆;张先富;胡志刚;侯丹峰【作者单位】河南科技学院,机电学院,河南新乡,453003;河南大学,物理与电子学院,河南开封,475004;永城职业学院,电子工程系,河南永城,476600;河南科技学院,机电学院,河南新乡,453003;河南科技学院,机电学院,河南新乡,453003【正文语种】中文【中图分类】O641;O482.5【相关文献】1.Gen Gd(n=1~10)团簇结构稳定性与磁性研究 [J], 岳莉;吴位巍;张颂2.Mn掺杂(ZnTe)12团簇结构和磁性质 [J], 谢建明;汪恒;刘成林3.TM1~2Mn@(BN)48笼团簇的结构与磁性研究 [J], 梁文娟; 刘永文; 荆补琴; 卢珍4.GanTi(n=1~10)团簇结构稳定性与磁性研究 [J], 岳莉; 吴位巍; 张颂5.NixPt(x=-9)合金团簇结构稳定性与磁性的第一性原理研究 [J], 岳莉因版权原因,仅展示原文概要,查看原文内容请购买。
WmBn(m+n≤7)团簇电子结构与光谱性质的计算研究张秀荣;尹琳;陈晨【摘要】采用密度泛函理论中的B3LYP方法,在LANL2DZ基组水平上优化了WmBn(m+n≤7)团簇的几何结构,得到了它们的基态构型,并对其自然键轨道( natural bond orbital ,NBO)以及振动光谱进行了计算研究。
结果表明:团簇中NBO电荷分布状况与团簇的对称性相关;W原子比B原子电荷调节能力强,易与其他原子形成化学键;最大频率以及最强峰值对应的振动模式大部分是B 原子或W原子的伸缩振动以及摇摆振动。
%Geometric structures of WmBn (m+n≤7) clusters are optimized by using the density functional theory at the B3LYP/LANL2DZ level.For the ground state configurations , the polarizability, magnetism,NBO(natural bond orbital ) and vibrational spectrum are analyzed .The calculated results show that the charge distribution of NBO is related to the symmetry of cluster;W charge regulation capability is better than B , which is easy to form chemical bonds;the strongest peak vibration mode corresponds to stretching or wagging vibration of atom B or W .【期刊名称】《江苏科技大学学报(自然科学版)》【年(卷),期】2013(000)004【总页数】10页(P399-408)【关键词】WmBn(m+n≤7)团簇;NBO;振动光谱【作者】张秀荣;尹琳;陈晨【作者单位】江苏科技大学数理学院,江苏镇江212003;江苏科技大学数理学院,江苏镇江212003;江苏科技大学材料科学与工程学院,江苏镇江212003【正文语种】中文【中图分类】O641过渡金属混合团簇具有特殊的电子结构和物理化学性质,为此成为目前团簇科学的热门课题[1-8].W合金材料,因其有很高的硬度、较好的耐磨、耐蚀性和延展性、抗高温氧化、易脱模、不粘着等特殊性质,引起了科技工作者的广泛关注[9-15].文献[16]利用电子散射衍射和映像分析技术研究了在烧结过程中WC颗粒的生长行为;文献[17]利用电沉积方法获得的Ni-W合金具有很强的延展性和抗拉强度,可以沿一个角度旋转180°而不使材料断裂.文献[18]对BmN (m=2~9)团簇的几何构型、电子结构、振动频率、自然键轨道(natural bond orbital,NBO)等性质进行了理论研究,得到了BmN (m=2~9)团簇结构的稳定性信息.文献[19]对AlBn+(n=2~10)团簇几何结构、稳定性、电子结构和成键特性进行了系统理论的研究,得到了AlBn+(n=2~10)团簇的最稳定结构,并且通过对其红外振动光谱的研究得出了硼原子间更容易成键的结论.文献[20]计算得出FeBn(n≤15)团簇基态结构中Fe的d轨道和B原子的P轨道存在着明显的杂化现象,研究表明FeB3、FeB5、FeB12和FeB15团簇较相邻团簇稳定.文献[21]对FeBn (n≤6)团簇的磁性做了系统研究,发现除了FeB5团簇外,FeBn (n≤6)团簇的总磁矩和Fe原子磁矩随团簇尺寸的增大而减小.文献[22]计算出了W6Sin0, ±(n=1,2)团簇基态及亚稳态结构,发现Si-Si之间不成键,并且分析了其芳香性及磁性.文中将对WmBn团簇的电子结构和光谱性质进行研究.1 计算方法文中采用量子化学程序Gaussian03,在B3LYP/LANL2DZ水平上对WmBn混合团簇基态结构的电子结构和光谱性质进行了计算研究.为了寻找WmBn的基态结构,设计了WmBn(m+n≤7)团簇多种可能的几何结构,进行几何参数优化,把优化之后无虚频的结构定为稳定结构,把能量最低且没有虚频的结构定为基态结构.然后对基态结构的电子结构和光谱性质进行了计算研究.2 结果与讨论2.1 基态结构图1给出了WmBn(m+n≤7)团簇的基态构型.图中颜色较深的大球为W原子,较浅的小球为B原子.通过图1可知,当只有一个W原子时,团簇的构型是平面结构,当m≥2,且m+n≥4时,除W3B团簇外,其余团簇的基态结构均为立体结构.对称性最高的是W2B3团簇,其对称性为D3h.图1 WmBn(m+n≤7)团簇的基态结构Fig.1 Ground state structures of WmBn(m+n≤7) clusters2.2 自然键轨道分析文中采用自然键轨道(NBO)方法分析了WmBn (m+n≤7)团簇的自然电荷布局以及成键性质.处于稳定状态的原子,核外电子将尽可能地按能量最低原理排布,另外它们还要遵守泡利不相容原理和洪特规则.处于基态的原子中电子以1s,2s,2p,3s,3p,4s,3d,4p,5s,4d 次序排布,从而得出电子组态.自由W原子最外电子层排布为5s25p65d46s2.B原子最外电子层排布为1s22s22p1.表1列出了WmBn (m+n≤7)团簇基态结构的各轨道上的NBO电荷分布.表1 WmBn (m+n≤7)团簇基态结构的自然电子组态和电荷Table 1 Natural electron configuration and atomic charge of the ground state structures of WmBn(m+n≤7) clusters团族自然电子组态电荷团族自然电子组态电荷WBWB21B2S(1.58)2p(1.44)-0.0261B2S(1.26)2p(1.72)3p(0.01)0.0053W6S(0.97)5d(4.96)6p(0.12)-0.009 WB3WB41B2S(1.12)2p(2.06)3p(0.01)-0.1851B2S(1.02)2p(2.09)-0.113 2B2S(1.12)2p(2.06)3p(0.01)-0.1852B2S(1.01)2p(2.21)3p(0.01)-0.2333B2S(1.11)2p(1.80)3p(0.01)0.0823B2S(1.02)2p(1.91)3p(0.01)0.0654W6S(0.82)5d(4.87)6p(0.06)7S(0.01)0.2884B2S(1.06)2p(2.02)3p(0.01)-0.096 5W6S(0.78)5d(4.83)6p(0.05)0.377 WB5WB61B2S(0.97)2p(2.12)3p(0.01)-0.1011B2S(0.95)2p(2.01)3p(0.01)0.032 2B2S(1.04)2p(2.00)3p(0.01)-0.0582B2S(1.00)2p(2.06)3p(0.01)-0.068 3B2S(1.04)2p(2.00)3p(0.01)-0.0583B2S(0.96)2p(1.99)3p(0.01)0.044 4B2S(0.97)2p(2.12)3p(0.01)-0.1014B2S(1.00)2p(2.06)3p(0.01)-0.0685B2S(1.05)2p(1.93)3p(0.01)0.0165B2S(0.95)2p(2.01)3p(0.01)0.0326W6S(0.80)5d(4.88)6p(0.05)0.3026B2S(0.66)2p(2.65)3p(0.02)-0.3317W6S(0.59)5d(4.98)6p(0.11)6d(0.01)0.359W2BW2B21W6S(1.00)5d(4.82)6p(0.15)0.0841W6S(0.34)5d(2.55)6p(0.14)-0.001 2W6S(1.00)5d(4.82)6p(0.15)0.0842B2S(0.68)2p(0.81)0.0013B2S(1.37)2p(1.79)3p(0.01)-0.1693B2S(0.68)2p(0.81)0.0014W6S(0.34)5d(2.55)6p(0.14)-0.001 W2B3W2B41B2S(1.10)2p(1.91)3p(0.01)-0.0161W6S(0.26)5d(2.51)6p(0.15)0.111 2B2S(1.10)2p(1.91)3p(0.01)-0.0162B2S(0.46)2p(1.11)-0.068 3B2S(1.10)2p(1.91)3p(0.01)-0.0163B2S(0.46)2p(1.11)-0.0684W6S(0.57)5d(5.35)6p(0.10)6d(0.01)0.0244B2S(0.53)2p(1.01)-0.0435W6S(0.57)5d(5.35)6p(0.10)6d(0.01)0.0245B2S(0.53)2p(1.01)-0.0436W6S(0.26)5d(2.51)6p(0.15)0.111 W2B5W3B1B2S(0.96)2p(2.04)3p(0.01)-0.2092W6S(1.05)5d(4.81)6p(0.10)6d(0.01)0.0903B2S(1.03)2p(1.90)3p(0.02)0.0603W6S(0.78)5d(4.97)6p(0.27)6d(0.01)0.029 4B2S(0.96)2p(2.04)3p(0.01)-0.0154B2S(1.44)2p(1.70)3p(0.01)-0.1495B2S(1.03)2p(1.90)3p(0.02)0.0606W6S(0.68)5d(5.17)6p(0.17)6d(0.01)0.0287 W6S(0.47)5d(5.05)6p(0.42)6d(0.02)0.091因为每个原子处在不等价的空间位置,所受到的势场不同,当一部分原子失去电荷,另一部分得到电荷就出现了电荷转移现象.从表1发现,钨原子的6s轨道上的NBO电荷分布在0.22至1.05,钨原子NBO电荷主要集中在5d轨道,分布在2.42至5.35,同时6p轨道上的NBO电荷在0.05至0.52,分布相对较少.当m<3时,部分钨原子含有6d轨道,当m≥3时,钨原子都含有6d轨道,但是电荷分布均极少,在0.01至0.02范围内.只有WB和WB3团簇中的钨原子含有7s轨道,电荷分布也较少.由上述分析以及表1看出,除W2B4、W3B2、W3B4和W4B2团簇外,钨原子的5d轨道和6p轨道得到电子,6s轨道失去电子,说明钨原子内部出现了轨道杂化现象,结论与WnNim (n+m=8)团簇[23]自然轨道电荷转移现象相吻合.硼原子中2s轨道上的NBO电荷分布在0.45至1.82,硼原子NBO电荷主要集中在2p轨道(除了WB 团簇),此时NBO电荷分布在0.81至2.65,大部分硼原子含有3s轨道,但是电荷分布极少.据此以及表1可知,除W2B2,W3B4和W4B2团簇外,硼原子的2p轨道得到电子,2s轨道失去电子,硼原子内部也发生了轨道杂化,结论与WnNim (n+m=8)团簇[23]自然轨道电荷转移现象也相吻合.深入探讨可以发现,WmBn (m+n≤7)团簇中,除了WB2和W2B2团簇外,与钨原子相邻的硼原子上的2p轨道得到电荷数目大于2s轨道失去电荷数目,而钨原子6s 轨道失去电子数目大于5d,6d轨道得到的电荷数目,说明钨原子中的部分电荷转移到了硼原子,所以在两个轨道电子的相互作用的过程中,电荷从钨原子转移到硼原子,形成了W-B键.而对于WB2和W2B2团簇,W原子得到了电荷,B原子得到了电荷,即B原子2p轨道上的部分电荷转移到了W原子上,W原子变成了电荷的受体,使得W原子的6s轨道和B原子的2p轨道形成了复杂的化学键.综上所述,在W原子与B原子内部轨道之间发生了电荷转移现象,即轨道杂化现象,这些杂化轨道在原子之间相互作用形成化学键,决定了团簇的稳定性和特殊的物理化学性质.结合表1和基态结构图发现WmBn (m+n≤7)团簇中NBO电荷分布状况与团簇的对称性相关,在对称性较高的团簇中位置相同的原子,其NBO电荷分布状况也相同.如在W2B3团簇(对称性为D3h)中,1B,2B和3B的电荷分布相同均为2S(1.10)2p(1.91)3p(0.01),4W和5W原子上的NBO电荷分布均为6S(0.57)5d(5.35)6p(0.10)6d(0.01);W4B3团簇(对称性为C3V)中,对称位置上的2W和3W的NBO电荷分布状况相同,4W和5W,6B和7B的NBO电荷分布状况也分别相同.表1同时也列出了WmBn (m+n≤7)团簇基态结构中各个原子上的净电荷分布情况.在WmBn (m+n≤7)团簇中,W原子的净电荷分布在-0.065e至0.377e,B原子的净电荷分布在-0.331e至0.083e,W原子的净电荷分布比B原子的净电荷分布范围较大一些,说明W原子比B原子对电荷调节能力要强,易与其他原子形成化学键.在W原子和B原子相互作用形的过程中,原子间发生了电荷转移,大部分W原子呈正电性,大部分B原子呈负电性,B原子从W原子得到电子,这部分W原子提供电荷,在形成团簇的过程中起主要作用.2.3 光谱分析一个多原子的化合物分子可能存在很多振动方式,但并不是所有的分子振动都能吸收红外光.当分子的振动不致改变分子的偶极矩时,它就不能吸收红外辐射,不具有红外活性,即红外光谱的吸收强度由振动中的偶极距变化大小决定的.分子振动(和点阵振动)与转动引起分子极化率发生变化,则产生拉曼光谱.图2给出了WmBn(m+n≤7)团簇基态结构的红外(IR)光谱图和拉曼(Raman)光谱图.通过GaussView来判定各团簇光谱峰值所对应频率的振动方式的归属情况.由图2可以看出,WB2团簇红外光谱只有1个振动峰位于频率695 cm-1处,该处的振动模式为2个B原子沿W-B键的方向做不对称伸缩振动;拉曼光谱中有3个振动峰,其最强振动峰位于频率847 cm-1处,该处的振动模式为2个B原子沿W-B键的方向做对称伸缩振动.WB3团簇红外光谱有4个振动峰,其最强振动峰所对应频率1 010 cm-1处,振动方式为3个B原子沿着B-B键方向往返伸缩振动,在振动过程中化学键的键长发生变化,次强峰位于频率877cm-1处,该处的振动模式为2个B原子沿W-B键的方向做不对称伸缩振动,从而引起化学键键长的变化;拉曼光谱中有只有1个振动峰,位于频率827 cm-1处,其振动模式为W原子和其中1个B原子固定不动,其他2个B原子做对称伸缩振动;WB3团簇的IR和Raman主峰都在高频段,表明在高频段该团簇红外和拉曼活性都很强.WB4团簇红外光谱有多个振动峰,其中最强振动峰位于频率443 cm-1处,其振动模式为B原子做面内弯曲振动;拉曼光谱和红外光谱一样,几乎整个频段上都有峰值分散出现,其最强振动峰位于频率995 cm-1处,其振动模式为不对称伸缩振动;次强峰位于波数672 cm-1处,此时2个B原子做不对称伸缩振动,而引起的结构变形以及化学键键长的改变.WB5团簇红外光谱有4个振动峰,最强峰位于频率529 cm-1处,其振动模式为平面剪式振动;拉曼光谱中有多个振动峰,最强振动峰位于频率1 205 cm-1处,振动模式为平面对称伸缩振动.a)WB2b) WB3d) WB5e) WB6f) W2Bg) W2B2h) W2B3i) W2B4 j) W2B5 k) W3B l) W3B2 m) W3B3 n) W3B4p) W4B2q) W4B3r) W5Bs) W5B2t) W6B图2 WmBn (m+n≤7)团簇IR光谱和Raman光谱Fig.2 IR spectra and Raman spectra of WmBn(m+n≤7) clustersWB6团簇红外光谱在1 000 ~1 200 cm-1段内有连续多个峰值,说明在这个频段内该团簇红外活性很好,且主峰值相差较小,最强振动峰位于频率1 112 cm-1处,振动模式为位于几何中心的B原子做不对称伸缩振动;拉曼光谱中有多个振动峰,最强振动峰位于频率805cm-1处,振动模式为平面剪式振动.W2B团簇红外光谱有3个振动峰,最强峰的振动模式是B原子的对称伸缩振动,位于频率727 cm-1处;拉曼光谱有3个振动峰,最小峰强度小的几乎看不出来,最强峰的振动模式为平面对称伸缩振动,位于727 cm-1处.W2B2团簇红外光谱有4个振动峰,最强振动峰位于频率为584 cm-1处,振动模式为2个B原子在各自的平面内做摇摆振动;拉曼光谱中有多个振动峰,最强峰位于频率为797 cm-1处,振动模式为2个B原子分别沿W-B键做对称伸缩振动.W2B3团簇红外光谱有3个振动峰,在频率288 cm-1处的是最强峰,振动模式为下方2个B原子的上下摇摆振动; 拉曼光谱中只有2个振动峰,最强峰位于频率为700 cm-1处,振动模式为呼吸振动.W2B4团簇红外光谱有多个振动峰,最强振动峰位于467 cm-1处,振动模式为位于两翼的B原子做面外摇摆振动;拉曼光谱的振动峰几乎分布在整个频段上,最强峰位于913 cm-1处,振动模式为各个B原子在同一平面内做剪式振动.W2B5团簇红外光谱在200~450 cm-1段内有多个连续的峰值,但强度都很小,其中最强振动峰峰位于339 cm-1出,振动模式为面外弯曲振动;拉曼光谱的峰值几乎分散整个频段,在1 179 cm-1处的振动最强,振动模式为沿B-B键方向的对称伸缩振动.W3B团簇的红外光谱在258 cm-1处出现了最强谱峰,是由于B原子做面外摇摆振动产生的;拉曼光谱和红外光谱类似,分布在相对较高和较低的频段,在718 cm-1处出现了最大谱峰,是由于W原子之间的相互伸缩振动而产生的.W3B2团簇红外光谱和拉曼光谱都有多个振动峰,红外光谱最强振动峰位于545 cm-1处,振动模式为面内弯曲振动;拉曼光谱最强振动峰位于545 cm-1处,同红外光谱相同,振动模式为面内弯曲振动.W3B2团簇与WB4团簇的分布情况类似,IR和Raman的峰值比较分散,几乎整个频段上都有.W3B3团簇红外光谱有多个振动峰,分布在整个频段上,最强振动峰位于621 cm-1处,位于左右下方的2个B原子做前后面外摇摆振动;拉曼光谱虽然有多个振动峰,但是只有一个显著的峰,位于735 cm-1处,该振动模式为3个B原子沿W-B键做对称伸缩振动.W3B4团簇红外光谱分布在整个频段上,最强振动峰位于739 cm-1处,振动模式为扭曲振动;拉曼光谱也分散在整个频段,最强振动峰位于269 cm-1处,振动模式为剪式振动.W4B团簇红外光谱有多个较强振动峰,最强振动峰位于646 cm-1处,振动模式为B 原子的伸缩振动;拉曼光谱的最强振动峰位于646 cm-1处,振动模式为B原子的伸缩振动.W4B2团簇的红外光谱和拉曼光谱峰都比较多,而且都是分段集中;红外光谱最强振动峰位于627 cm-1处,振动模式为2个B原子的前后不对称摇摆振动;拉曼光谱的最强振动峰位于703 cm-1处,振动模式为2个B原子的前后对称伸缩振动.W4B3团簇红外光谱和拉曼光谱都有多个振动峰,而且都比较分散;红外光谱最强振动峰位于456 cm-1处,振动模式为位于顶端的B原子的摇摆振动;拉曼光谱最强振动峰位位于822 cm-1处,振动模式为3个B原子沿B-B键做伸缩振动.W5B团簇红外光谱在100~260 cm-1频段上有多个连续的峰值,最强振动峰位于118 cm-1处,振动模式为W-W键和W-B键的伸缩振动;拉曼光谱有多个振动峰,最强振动峰位于118 cm-1处,振动模式为W-W键和W-B键的伸缩振动.W5B2团簇红外光谱只有3个振动峰,193 cm-1处的最强振动峰振动模式为伸缩振动;拉曼光谱有多个振动峰,最强振动峰位于595 cm-1处,振动模式表现为2个B原子的剪式振动.W6B团簇红外光谱在70~230 cm-1范围内集中了许多小峰,最强振动峰位于549 cm-1处,振动模式为B原子的摇摆振动;拉曼光谱主要集中分布在70~300 cm-1频段间,最强峰位于76 cm-1处,振动模式为剪式振动.通过对WmBn(m+n≤7)团簇基态结构的光谱分析可知,WmBn团簇的红外和拉曼光谱的振动峰都较多,主要分布在30至1 288 cm-1处,最大频率以及最强峰值对应的振动模式大部分是B原子或W原子的伸缩振动以及摇摆振动.周围环境的改变对频率的变化影响较小,力常数较大致使体系的电偶极矩变化大,所以基团伸缩振动常在高频处出现吸收,W原子之间化学键的力常数较大,容易产生较高的振动频率,所以一般来说频率值较大的几乎都是W原子之间的伸缩振动.WB2,WB3,W2B,W2B2和W2B3团簇红外光谱的吸收峰相对较少,振动模式大部分为伸缩振动.从红外光谱强度图上可以看出,WB5,WB6和W4B2团簇的最强峰分别达到了58.69 km.mol-1,42.20 km.mol-1和40.53 km.mol-1,和其他团簇的最强峰相差较多,说明这两个团簇的振荡模式大大改变了其结构中的电子云分布,偶极距的急剧变化导致了这一现象.3 结论采用密度泛函理论中的B3LYP方法,在LANL2DZ基组水平上优化了WmBn(m+n≤7)团簇的几何结构,得到了它们的基态构型,并对其自然键轨道(NBO)以及振动光谱进行了计算研究.对NBO的分析可知,团簇中NBO电荷分布状况与团簇的对称性相关,在对称性较高的团簇中位置相同的原子,其NBO电荷分布状况也相同,W原子比B原子对电荷调节能力强,易和其它原子形成化学键;振荡光谱分析表明,振动频率主要分布在30~1 288 cm-1处,最大频率以及最强峰值对应的振动模式大部分是B原子或W原子的伸缩以及摇摆振动.参考文献[1] Suetin D V, Shein I R, Ivanovskii A L. Structural, electronic properties and stability of tungsten mono-and semi-carbides: a first principles investigation [J].Journal of Physics and Chemistry of Solids,2009, 70(1): 64-71.[2] Ding F, Larsson P, Larsson J A, et al. The importance of strong carbon-metal adhesion for catalytic nucleation of single-walled carbon nanotubes[J].Nano Lett,2008, 8(2): 463-468.[3] 张秀荣,高从花,吴礼清,等, WnNim (n+m≤7; m=1,2) 团簇电子结构与光谱性质的理论研究[J].物理学报,2010,59(8):5429-5438.Zhang Xiurong, Gao Conghua, Wu Liqing,et al.The theory study of electronic structures and spectram properties of WnNim (n+m≤7; m=1,2) clusters[J].Acta Physica Sinica,2010,59(8):5429-5438.(in Chinese)[4] Ludwig A, Cao J, Dam B, et al. Opto-mechanical characterization of hydrogen storage properties of Mg-Ni thin film compositionspreads[J].Applied Surface Science,2007, 254(3): 682-686.[5] 仇毅翔,李佳,王曙光.配体稳定的二元过渡金属团簇[PdAu8(PR3)8]2+ (R=Me, OMe, H, F, Cl, CN)的量子化学理论研究[J].化学学报, 2010,68(7):611-616.Qiu Yixiang,Li Jia,Wang Shuguang. Theoretical investigations on ligand-stabilized binary transition-metal cluster [PdAu8(PR3)8]2+ (R=Me, OMe, H, F, Cl, CN)[J].Acta Chimica Sinica,2010,68(7):611-616.(in Chinese)[6] Zhang Lin, Zhang Caibei, Qi Yang. Local structure changes of Cu55 cluster during heating[J].Chinese Physics,2007, 16(1): 77-82.[7] 张秀荣,张伟,高从花,等.(OsB)n(n=1-6)团簇结构与稳定性的第一性原理研究[J]. 计算机与应用化学,2010,27(2):221-224.Zhang Xiurong,Zhang Wei,Gao Conghua,et al.First-principles study of structures and stable properties of (OsB)n(n=1-6) clusters[J].Computers and Applied Chemistry,2010,27(2):221-224.(in Chinese)[8] 华英杰,王崇太,孟长功.NiTi合金表面镍离子释放的电子结构理论计算[J]. 计算机与应用化学,2005,22(12):1115-1118.Hua Yingjie,Wang Chongtai,Meng Changgong.The theoretical calculation of electronic structure for the release of nickel ion on NiTi alloysurface[J].Computers and Applied Chemistry,2005,22(12):1115-1118.(in Chinese)[9] 朱佳,金华,李奕,等.TiO2(110)表面负载W3O10团簇构型和电子结构的理论研究[J].化学学报,2011,69(8):905-911.Zhu Jia,Jin Hua,Li Yi,et al.Theoretical study on the geometry and electronic structure of W3O10 clusters supported on the TiO2 (110) surface[J].Acta Chimica Sinica,2011,69(8):905-911.(in Chinese)[10] Yamaguchi W, Murakami J. Geometries of small tungstenclusters[J].Chemical Physics,2005, 316(1/3): 45-52.[11] 姚明珍,顾牡,梁玲,等. PbWO4晶体空位型缺陷电子结构的研究[J].物理学报,2002,51(1):125-128.Yao Mingzhen,Gu Mu,Liang Ling, et al.Electronic structures of defects associated with intrinsic vacancies in PbWO4 crystals[J].Acta Physica Sinica,2002,51(1):125-128.(in Chinese)[12] Suetin D V, Shein I R, Ivanovskii A L. Structural, elastic and electronic properties and formation energies for hexagonal (W0.5Al0.5)C in comparison with binary carbides WC and Al4C3 from first-principles calculations[J].Physica B, 2008, 403:2654-2661.[13] Wei A, Heath T C. Electronic structure calculations of gas adsorption on boron-doped carbon nanotubes sensitized with tungsten [J].Chemical Physics Letters,2009, 482(4/6): 274-280.[14] Shane M S, Adam W S, Michael D M. Optical spectroscopy of tungsten carbide WC[J].Journal of Chemical Physics,2002,116: 993-1002.[15] 邵泽旭,张启仁,刘延禹,等.CaWO4晶体中F型色心电子结构的研究[J].物理学报,2007,56(7):4089-4093.Shao Zexu,Zhang Qiren,Liu Yanyu,et al.Electronic structures of CaWO4 crystal with F type color center[J]. Acta Physica Sinica,2007,56(7):4089-4093.(in Chinese)[16] Mannesson K, Elfwing M. Analysis of WC grain growth during sintering using electron backscatter diffraction and image analysis[J].International Journal of Refractory Metals & Hard Materials,2008,26:449-455.[17] Tohru Yamasaki. High-strength nanocrystalline Ni-W alloys produced by electrodeposition[J].Mater Phys Mech,2000,1:127-132.[18] 马文瑾,王艳宾,张静,等.BmN (m=2-9)团簇的结构特征与稳定性[J]. 物理化学学报,2007, 23(2): 169-172.Ma Wenjin,Wang Yanbin,Zhang Jing,et al.Structure characteristics and stability of BmN(m=2-9) clusters[J].Acta Phys-Chim Sin,2007, 23(2): 169-172. (in Chinese)[19] 李莉莎,刘甫,孙久雨,等,AlBn+(n=2~10)团簇结构和红外振动光谱研究[J].光子学报,2011,40(2):321-326.Li Lisha,Liu Fu, Sun Jiuyu,et al.Theoretical study of structure and infrared vibration spectra about AlBn+(n=2~10) clusters[J]. Acta Photonica Sinica,2011,40(2):321-326. (in Chinese)[20] 刘火雁,雷雪玲,陈杭,等,FeBN(N≤15)团簇结构、电子性质和磁性的密度泛函理论研究[J].原子与分子物理学报,2011,28(2):258-266.Liu Huoyan,Lei Xueling,Chen Hang,et al.Structures,electronic properties and magnetisms of FeBN(N≤15) clusters:density functional theory investigations[J].Journal of atomic and molecular physics,2011,28(2):258-266. (in Chinese)[21] 杨致,闰玉丽,赵文杰,等,FeBN(N≤6)团簇的结构与磁性[J].物理学报,2007,56(5):2590-2595.Yang Zhi,Yan Yuli,Zhao Wenjie,et al.Structures and magnetism ofFeBN(N≤6) clusters[J].Acta Physica Sinica,2007,56(5):2590-2595. (in Chinese)[22] 张秀荣,刘小芳,康张李.W6Sin0, ±(n=1,2)团簇结构与电子性质的密度泛函理论研究[J].江苏科技大学学报:自然科学版,2010,24(6):619-624.Zhang Xiurong,Liu Xiaofang,Kang Zhangli.Density functional theory study on the structure and electronic properties of W6Sin0, ±(n=1,2)clusters[J].Journal of Jiangsu University of Science and Technology:Natural Science Edition,2010,24(6):619-624. (in Chinese)[23] 张秀荣,李扬,杨星. WnNim(n+m=8)团簇结构与电子性质的理论研究[J]. 物理学报, 2011,10:236-245.Zhang Xiurong,Li Yang,Yang Xing.Theoretical study on structural and electronic properties of WnNim(n+m=8) clusters[J].Acta Physica Sinica,2011,10:236-245. (in Chinese)。
团簇和纳米体系物理发展团簇和纳米体系是20世纪末发展起来的崭新领域,它所研究的对象是既不同于原子、分子,又不同于宏观物体的中间体系,尺度范围大约在0.1-100nm,这是人们过去从未进行研究的新领域,是人们认识物质世界的新层次。
它的丰富物理内涵,对物理学提出了新的挑战,也是当前物理学与其它学科交叉最富有活力的热点领域。
团簇简介原子和分子团簇,简称团簇(Cluster)或微团簇(microclusters),是几个乃至上千个原子、分子或离子通过物理或化学结合力组成相对稳定的微观和亚微观聚集体,其物理和化学性质随所包含的原子数目而变化。
团簇的空间尺度是几个至几百个纳米的范围,用分子描述显得太大,用小块固体描述又显得太小,许多性质既不同于单个原子分子,又不同于固体或液体,也不能用两者性质作简单线性外延和内插得到。
因此,人们把团簇看作是介于原子分子和宏观固体之间物质结构的新层次,有人称之为物质的“第五态”。
正因为如此,团簇可作为各种物质由原子分子向大块物质转变过程中的特殊物相,或者说它代表了凝聚态物质的初始状态,团簇的研究有助于我们认识大块凝聚物质的某些性质和规律。
团簇科学是研究团簇的原子组态和电子结构、物理和化学性质及其向大块物质演化过程中与尺寸的关联,团簇同外界环境的相互作用规律等。
团簇科学处于多学科交叉的范畴。
从原子分子物理、凝聚态物理、量子化学、表面科学、材料科学甚至核物理学引入的概念和方法交织在一起,构成当前团簇究的中心议题,并逐渐发展成一门介于原子分子物理和固体物理之间的新型学科。
团簇研究的基本问题是:弄清团簇如何由原子、分子一步步发展而成,以及随着这种发展、团簇的性质将如何变化,当尺寸多大时,团簇发展成宏观固体人们知道,由若干原子构成的分子,可在气相、液相和固相中稳定的单元存在,而团簇作为原子聚集体往往产生于非平衡条件,很难在平衡的气相中产生。
当团簇尺寸较小时,每增加一个原子,团簇的结构发生变化,即所谓重构。
第30卷第6期2014年12月 分 子 科 学 学 报 JOURNAL OF MOLECULAR SCIENCE Vol.30No.6 Decomber 2014收稿日期:2014-05-30基金项目:河南省科技厅基础与前沿资助项目(132300410209,132300410290);河南省教育厅自然科学基金资助项目(2010B140012,2011B140015).联系人简介:李根全(1961-),男,教授,博士,主要从事团簇电子结构研究.E-mail:genquan_li1961@163.com[文章编号]1000-9035(2014)06-0481-05CaSin(n=1~10)团簇结构、稳定性与光谱性质的理论研究白音布和1,张 宇2,尉乔南3,李根全*,2,张 帅2(1.通辽职业学院机电工程学院,内蒙古通辽028000;2.南阳师范学院物理与电子工程学院,河南南阳473061;3.南阳理工学院电子与电气工程学院,河南南阳473061)[摘 要] 利用密度泛函理论(DFT)中的B3LYP方法,在6-311G(d)基组的水平上系统研究了CaSin(n=1~10)团簇的几何构型、稳定性与光谱(红外与拉曼)性质.研究结果表明,Ca-Sin团簇构型是在CaSin-1构型上戴帽1个原子而形成的;当n≥4,CaSin团簇的最低能量结构均为立体构型;Ca原子的掺杂降低了体系的化学稳定性;CaSi3与CaSi5是幻数结构;在相同的观察频段内,CaSi3团簇的红外与拉曼活性在低频段均表现较好,而在高频段拉曼活性则表现较差,与之不同的是CaSi5团簇的红外与拉曼活性在整个频段内都表现的较好.[关键词] CaSin团簇;结构与稳定性;光谱性质[中图分类号] O 641 [学科代码] 150·30 [文献标识码] A 碱土金属硅团簇是近年来令人感兴趣的团簇,由于其具有特殊的物理化学性质,如:低电阻率、高温强度高、高比强度、高弹性模量、低导热系数等,可以广泛地应用于涂层材料、航空航天材料和许多其他产品[1].近几十年来,人们采用不同的理论方法研究碱土金属硅团簇.Binning等分别使用B3LYP/6-31G(2df)与G3XMP2方法研究了BenSin和Be2nSin(n=1~4)团簇的结构与相对能量[2];Fan等采用MP2(full)/6-31G(d)方法计算了MgSin±(n=2~10)团簇的结构、能量、电荷布居与偶极矩[3];Ning等采用MP2(full)/6-31G(d)方法得到了CaSin±(n=2~10)的基态结构,计算了绝热电子亲和势、离解能并与BeSin和MgSin团簇进行了比较[4].虽然目前关于碱土金属硅团簇的研究报道很多,但是关于Ca-Sin团簇稳定性与光谱性质的理论研究还未见报道.因为密度泛函理论已成功用于计算分子的结构与光谱性质[5-6],所以我们采用B3LYP方法,选取6-311G(d)基组,优化计算CaSin(n=1~10)团簇的几何结构,确定了最低能量构型,随后讨论了CaSin团簇的稳定性与光谱性质.该研究旨在为今后实验研究及混合团簇的开发提供可靠的理论依据.1 计算方法采用密度泛函理论中的B3LYP交换关联方法[7],同时选取6-311G(d)基组,使用Gaussian 09程序[8]对CaSin(n=1~10)掺杂团簇的几何结构与电子性质进行计算,计算中能量收敛精度优于10-6a.u.为了确保所选计算方法的可靠性,首先计算了二聚体Si2与Ca2的键长与振动频率,得到的结果分别为:0.216 7nm,540.37cm-1;0.442 0nm,57.87cm-1,与实验结果0.224 6nm,510.98cm-1[9];DOI:10.13563/ki.jmolsci.2014.06.0070.427 7nm,64.92cm-1[10]基本一致.通过上述对比说明选用的基组对于CaSin团簇是适用的.2 计算结果与讨论2.1 团簇的几何结构CaSin(n=1~10)团簇的构型设计如下:首先推测初始结构,随后以已报道文献Sin与XSin为基础构型,在基础构型的不同位置采用戴帽、置换、填充等方式构建初始结构,最后采用不同的自旋多重态(1,3,5),对每一个结构都进行了优化,得到了多个稳定结构,根据其能量的高低排序为na,nb与nc,na为团簇的最低能量结构.CaSin(n=1~10)团簇的结构如图1所示.图1 CaSin(1~10)团簇结构如图1所示,对于CaSi二聚体,基态为3∑,Ca—Si键长为0.305 8nm,具有C∞v对称性.当n≥2,Ca—Sin团簇的基态均为自旋单重态.CaSi2的最低能量结构2a为等腰三角形,基态为1 A1,具有C2v对称性,Si—Ca键长为0.301 4nm(腰),Si—Si键长为0.227 5nm(底边),键脚Si—Ca—Si为44.34°,可视284 分子科学学报 第30卷 为由1个Si原子戴帽于CaSi二聚体而生成的;2b是Ca原子处于中间位置的直线结构.CaSi3的最稳定构型3a为平面菱形(Cs,1 A′),它可以看做是在CaSi2的基态结构上带帽1个原子而得到;3b为畸变的三角锥构型,其能量比3a的高0.53eV.从n=4开始,CaSin团簇的最低能量构型开始由平面结构变为立体结构.当1个Ca原子戴帽于3a结构,得到2个四角锥4a和4b,它们的区别在于电子态与对称性不同,其中4b构型的能量高于4a构型0.03eV.CaSi5的基态结构5a是对称性为C1的三棱柱;5b是在3a构型基础上戴帽2个原子的立体结构,其能量比5a高0.62eV.CaSi6的基态结构6a为1个原子戴帽于5a构型,对称性为Cs;6b可以看做在四角双锥Si6的基础上戴帽1个Ca原子;6c是1个五角双锥结构;6b与6c的能量比6a分别高0.47和0.52eV.CaSi7的基态结构7a是对称性为C1的六角双锥结构;7b构型是在6c构型的基础上戴帽1个原子而得到的;7c构型类似于7a,区别在于Ca原子的位置不同;7b与7c的能量比7a分别高0.11和0.44eV.CaSi8的最低能量结构8a具有C1对称性,可以看作在六面体结构Si8的基础上带帽1个Ca原子得到;8b构型近似于8a,区别在于Ca原子的位置不同,其能量比8a的高0.41eV;8c构型是1个3-3-3的层状立体结构,其能量比8a的高0.57eV.CaSi9的最稳定构型9a是具有C1对称性的五棱柱;9b构型可以看做是在8a构型的基础上戴帽1个原子而得到的;9c构型为Ca原子位于顶点的笼状结构,9b与9c的能量比9a分别高0.40与0.71eV.CaSi10的基态结构10a为Si原子戴帽于9c构型;10b构型与10a类似,区别在于Ca原子的位置不同;10c构型为1个Ca原子戴帽于9a构型.CaSin团簇的生长行为可总结为:当n=1~3时,CaSin团簇的基态结构为平面构型,当n≥4时,CaSin团簇的基态结构则为立体构型.大多数CaSin团簇构型是在CaSin-1构型上戴帽一个原子而形成的.2.2 团簇的稳定性为了研究CaSin团簇的化学稳定性,采用相同的B3LYP方法,计算了CaSin团簇的结合能Eb、二阶能量差分Δ2E、离解能Ef与能隙Egap,并与Sin+1[11]团簇进行了比较.图2给出了CaSin团簇的结合能Eb、二阶能量差分Δ2E、离解能Ef与能隙Egap随团簇尺寸变化的趋势图.Eb,Δ2E,Ef的计算公式如下:Eb=[nE(Si)+E(Ca)-E(CaSin)]/(n+1),(1)Ef=E(CaSin-1)+E(Si)-E(CaSin),(2)Δ2E=E(CaSin+1)+E(CaSin-1)-2E(CaSin),(3)(1)-(3)式中E(Si),E(Sin),E(Ca),E(CaSin),E(CaSin-1)与E(CaSin+1)分别表示对应团簇或原子最低能量结构的总能量.图2 CaSin团簇的平均结合能Eb(a)、二阶能量差分Δ2E(b)、离解能Ef(b)、能隙Egap(c)随团簇尺寸的变化趋势图由图2(a)可以看出,随着团簇尺寸的增大,CaSin团簇的结合能Eb逐渐增加,表明团簇能够一直获取能量;但是当n≥8,Eb的增加趋于平缓,这是由于随着原子的增加,成键趋于饱和导致的.此外,掺杂团簇的Eb小于Sin+1[11]团簇的Eb,说明Ca原子的掺入可以增强团簇的化学活性.图2(b)给出了离解能Ef及能量二阶差分Δ2E随团簇尺寸变化的曲线,可以看出两条曲线的变化规律基本一致,均呈现波384第6期 白音布和,等:CaSin(n=1-10)团簇结构、稳定性与光谱性质的理论研究动变化趋势.当n=3,5,8,10时,Δ2E出现峰值,说明对应团簇相对稳定;当n=3,5,8,10时,Ef出现最大值,表明CaSi3,CaSi5,CaSi8与CaSi10团簇的稳定性相对较大.图2(c)给出了CaSin团簇与Sin+1[11]团簇的HOMO与LUMO之间的能隙Egap随团簇尺寸的变化规律,可以发现CaSin团簇的Egap随Si原子数的表现出明显的“奇-偶”振荡现象.即当n为奇数(1,3,5,7,9)时,均出现最大值,这表明CaSi,CaSi3,CaSi5,CaSi7与CaSi9团簇的稳定性相对其邻近团簇较强.与Sin+1[11]团簇的Egap对比发现掺杂团簇的能隙Egap小于纯硅团簇,反映了Ca原子的掺杂提高了团簇的化学活性,这与之前Eb的分析结果是相符的.结合图2(a-c)的分析,可以断定CaSi3与CaSi5团簇是幻数结构.2.3 团簇红外与拉曼光谱性质红外与拉曼光谱都属于分子振动光谱,偶极矩与极化率的微小变化都可以影响光谱的强烈变化[12].为了进一步研究CaSin团簇的稳定性,分析团簇基态结构的红外及拉曼光谱是十分有意义的.这里,计算了幻数结构团簇的红外及拉曼光谱并展示在图3.图3 CaSi3与CaSi5团簇的红外及拉曼光谱图3(a)和(b)分别给出了CaSi3团簇的红外及拉曼光谱.图3(a)中红外振动峰较多,有4个较强振动峰,最强的位于波数为458.63cm-1,该处的振动模式为Si2—Si3与Si2—Si4的伸缩振动.第二强峰位于波数111.43cm-1,振动模式为Si2—Ca4以Si1—Si3为轴向做上下摇摆振动.另外2个弱振动峰分别位于波数238.44与165.24cm-1处,对应的振动模式分别为Si1—Si3的伸缩振动、Si3—Ca4与Si1—Ca4做平面左右摇摆振动.图3(b)中有3个强拉曼振动峰,分别位于波数189.41,238.44与458.67cm-1处,相应的振动模式分别为Si1,Si2,Si3与Ca4之间的呼吸振动、Si1—Si3的伸缩振动与Si2—Si1和Si2—Si3之间的伸缩振动.结合图3(a)和(b)可以发现,CaSi3团簇的红外与拉曼活性在低频段均表现较好,而在高频段则拉曼活性表现较差.对于CaSi5团簇,图3(c)和(d)分别给出了其对应的红外及拉曼光谱.红外光谱图中仅有一个较强的振动峰,振动峰位于波数为217.96cm-1,振动模式为Ca6—Si3,Ca6—Si5原子间的伸缩振动.拉曼光谱中有3个强振动峰,最强的振动峰位于波数为238.82cm-1,对应振动模式同样为Ca6—Si3,Ca6—Si5原子间的伸缩振动;次强峰位于波数217.96cm-1,振动模式为Ca6—Si5,Ca6—Si1原子间的伸缩振动;第三强峰位于波数为151.73cm-1,振动模式为Si1—Si5的摇摆484 分子科学学报 第30卷 振动.观察图3(c)和(d)可以看出,CaSi5团簇的红外与拉曼活性光谱在整个频段内都表现的较好.3 结 论采用密度泛函理论(DFT)在B3LYP/6-311G(d)水平上系统研究了Ca2Sin(n=1~9)团簇的几何结构、稳定性与光谱性质.结果发现:当n≥4,CaSin团簇最低能量结构为立体构型,n=1,2,3时则为平面构型.Ca原子的掺入提高了纯Sin的化学活性.CaSi3与CaSi5团簇是幻数结构.光谱性质分析表明CaSi3团簇的红外与拉曼活性在低频段均表现较好,而在高频段则拉曼活性表现较差.CaSi5团簇的红外与拉曼活性光谱在整个频段内都表现的较好.[参 考 文 献][1] FAN H W,YANG J C,LU W,et al.[J].J Phys Chem A,2010,114:1218-1223.[2] BINNING R C,BACELO D E.[J].J Phys Chem A,2005,109:754-758.[3] FAN H W,REN Z Q,YANG J C,et al.[J].J Mol Struct:THEOCHEM,2010,958:26-32.[4] NING H M,FAN HW,YANG J C,et al.[J].Comput Theor Chem,2011,976:141-147.[5] 张秀荣,崔彦娜,洪伶俐,等.[J].分子科学学报,2009,25(3):192-199.[6] 张秀荣,崔彦娜,洪伶俐,等.[J].分子科学学报,2009,25(2):109-115.[7] LEE C,YANG W,PARR R G.[J].Phys Rev B,1988,37:785-789.[8] FRISCH M J,CLEMENTE F R,TRUCKS G W,et al.Gaussian 09[CP].Gaussian Inc,Pittsburgh PA,2009.[9] HUBER K P,HERZBERG G.Constants of Diatomic Molecules[M].New York:Van Nostrand/Reinhold,1979:212-214.[10] HAMADA H.[J].Philos Mag,1931,12:50-55.[11] QIN W,LU W C,ZHAO L Z,et al.[J].J Phys Condens Matter,2009,21:455501/1-7.[12] 张秀荣,高从花,吴礼清,等.[J].物理学报,2010,59(8):5429-5438.Theoretical study on the geometries,stabilities,and spectrum properties of CaSin(n=1~10)clustersBAIYIN Bu-he1,ZHANG Yu2,YU Qiao-nan3,Li Gen-quan*,2,ZHANG Shuai 2(1.School of Mechanical and Electric Engineering Technology,Tongliao Technical College,Tongliao 028000,China;2.School of Physics and Engineering,Nanyang Normal University,Nanyang 473061,China;3.School of Electronic and Electrical Engineering,Nanyang Institute of Technology,Nanyang 473061,China)Abstract:The geometries,stabilities and spectrum(infrared and raman)properties of CaSin(n=1~10)clusters have been systematically studied by using density functional theory with B3LYP function.The research findings indicate that the lowest-energy energy CaSingeometries are obtained by one at-om capped on the stable CaSin-1clusters;the ground state structure of CaSinclusters are favor thethree-dimensional structure for n=4~10;the impurity calcium atom in the CaSinclusters can enhancethe chemical activity of silicon clusters.CaSi3and CaSi5are the magic numbers.Moreover,the spec-trum analysis show that the CaSi3cluster presents stronger Infrared and Raman activity at low fre-quency band.Unlike CaSi3cluster,the CaSi5cluster exhibit better Infrared and Raman activity at thewhole frequency band.Keywords:CaSinclusters;geometries and stabilities;spectrum properties584第6期 白音布和,等:CaSin(n=1-10)团簇结构、稳定性与光谱性质的理论研究。