过渡金属双原子掺杂硅团簇的理论研究

论文摘要

团簇作为连接微观原子分子和宏观凝聚态物质间的桥梁,具有许多奇特的性质,其物理、化学性质随尺寸发生显著变化。人们希望通过对团簇的研究,实现以团簇作为基元组装出各种纳米功能材料和器件的目标。硅基半导体材料在微电子工业上的重要地位,使得硅团簇倍受重视。然而由于硅原子具有sp3杂化性质,不适合作为稳定的组装单元,为了稳定硅团簇,人们提出通过掺杂一些其它原子的方法来饱和硅团簇表面的悬键,而过渡金属原子以其独特的电磁特性,成为被研究最多的掺杂原子。由于计算量与结构复杂程度等因素的制约,以往的研究以掺杂单个过渡金属原子为主,对于多个金属原子共同掺杂的情况研究较少。本文通过基于密度泛函的第一性原理方法,对3d过渡金属双原子掺杂双棱柱硅笼MM’Si18 （M,M’=Ti,V,Cr,Mn,Fe,Co,Ni, Cu,Zn）以及过渡金属钨原子掺杂硅团簇WmSin（m=1,2;n=1-18）的几何结构、电磁性质、稳定性与电子数规则等进行了系统的研究。采用第一性原理计算方法,我们对同种3d金属双原子掺杂双棱柱硅团簇M2Si18（M =Ti,V,Cr,Mn,Fe,Co,Ni,Cu,Zn）进行了系统地几何优化、电荷分析和能量计算。从几何结构来看,只有当掺杂金属原子的3d电子数在半满以下时.掺杂后团簇能够形成类双六棱柱的结构。自旋多重度和对称性与两过渡金属原子间的距离之间相互影响。除MMn外的3d过渡金属掺杂的M2Si18其金属原子间的距离都大于金属二聚物键长。结构相似的异构体中,自旋多重度较高的构型其过渡金属原子间距离较大。对称（?）较高的,特别是类双六棱柱的结构,其过渡金属原子间距离小于对称性较低的异构体。除Zn原子外,电荷总是从硅原子向过渡金属原子转移。除Co2Si18和Cr2Si18为自旋三重态外,其余M2Si18的磁矩均湮灭了,这与掺杂后体系的几何对称性下降有关。过渡金属原子闭壳层18电子数规则不能准确判定MM’Si18的稳定性,而基于电子气闭壳层模型的电子数幻数规则结合自旋守恒约束条件则成功预言了Fe2Si18的特殊稳定性。为了丰富磁性与电子数样本,我们又进一步对异种3d金属双原子混合掺杂MM’Si18进行了系统研究。按照几何结构特点,MM’Si18团簇可以分为七类。除了某些Ti和V原子与Cr、Mn、Fe、Co等磁性原子混合掺杂外,其余情况下团簇中过渡金属间的距离均大于相应金属二聚物中的键长值。一些磁性较高的二聚物如VMn和CrMn:惨杂硅团簇MM’Si18后,体系的磁矩都下降到6μB以内,说明硅笼的存在导致磁性下降。在CrZuSi18以及所有不包含Zn原子的团簇中,电荷总是从硅原子向过渡金属原子转移。磁性方面,过渡金属原子的组合方式与团簇的磁性间存在明显的周期规律性,异种原子混合掺杂的团簇中只有3种磁性完全湮灭,其余团簇都具有一定的磁性,同时硅原子对整体磁性的贡献有所增加,特别是出现了硅原子对磁矩贡献大于金属原子的情况,这可能主要是由于体系对称性下降造成的。含有Ti原子的MM’Si18团簇往往比较稳定,而含有Zn的团簇则稳定性较差。过渡金属闭壳层的18电子数规则不能正确判定MM’Si18的稳定性,而基于电子气闭壳层的电子数规则结合自旋守恒条件约束,虽然在同种掺杂时成功预言了满足幻数34价电子要求的Fe2Si18的特殊稳定性,但是进一步研究发现在异种金属混合情况下幻数价电子的团簇稳定性明显低于Fe2Si18,因此这一电子数规则的确定性也令人质疑。除电子数规则外,团簇的对称性、掺杂原子的种类、过渡金属间距离、平均配位数、电子转移、磁矩分布都与MM’Si18团簇稳定性具有密切关系。随后我们对单个钨原子掺杂的’WSin（n=1-12）进行了系统研究。随着硅原子数从1增加到12,团簇中钨原子的平衡位置从团簇的顶点、表而逐渐向团簇内部转移,到硅原子数目n=9时钨原子完全落入由硅原子构成的笼形结构中心位置,构成了包裂钨原子的硅笼结构。在这一转变临界尺度前后,WSin团簇的几何结构参数、电荷密度分布、稳定性等都发生了明显的变化。电荷总是从硅原子向钨原子转移,随着钨原了内嵌,钨原子上的电子数目明显增加。磁性方而,除子WSi和WSi2之外,其余基态结构均为自旋单重态无磁矩。掺杂钨原子后的硅团簇稳定性高于相应尺寸的纯硅团簇,最稳定的是六棱柱结构的WSi12,这与实验和之前的计算结果符合。最后我们还对两个钨原子掺杂的W2Sin（n=1-18）进行了系统的计算。随着硅原子数目从7增加到18,两个钨原子依次嵌入硅笼,当n=10时一个钨原子完全嵌入硅笼中而另一个仍位于笼外,当n=14时,另一个的钨原子也开始嵌入硅笼,到n>16后两个钨原子完全嵌入硅笼。在钨原子依次嵌入的两个临界点前后,体系的几何结构、电荷分布和稳定性等许多物理量都出现了明显变化。随着尺寸的增加,团簇中两个钨原子的距离总体上呈现增加的趋势,体系表现出以钨原子为中心的笼形局域化分布趋势。虽然双钨原子掺杂后体系对称性较低,但存在着一些棱柱或层状结构随硅原子增加而连续演化的路径。电荷总是从硅原子向钨原子转移。内嵌的钨原子上的电荷数目大于笼外的钨原子,这一差异随着笼外钨原子的逐步内嵌而减小。在n=7-18的尺度范围之内,W2Sin的基态结构均为自旋单重态,净磁矩为零。结构最稳定的是W2Si18。在两个内嵌点前后,体系的平均结合能与嵌入能出现明显增加,说明内嵌有利于团簇稳定性的提高。同时片段能在内嵌点出现了明显的峰值,能量出现了明显地突变,这可以作为判定团簇尺度结构相变的重要依据。

论文目录

• 摘要
• Abstract
• 第一章 引言
• 1.1 团簇与团簇物理学
• 1.1.1 团簇的定义与性质
• 1.1.2 团簇物理学研究的主要内容与意义
• 1.1.3 团簇物理学的研究历史
• 1.2 硅团簇的研究发展现状
• 1.3 过渡金属原子掺杂硅团簇实验方面的进展
• 1.4 过渡金属原子掺杂硅团簇理论方面的进展
• 1.4.1 寻找金属内嵌的临界点
• 1.4.2 寻找影响硅团簇稳定性的电子数规则
• 1.4.3 多金属掺杂的硅团簇、硅纳米管等
• 1.5 本文中的主要内容
• 第二章 理论基础和计算方法
• 2.1 第一性原理与密度泛函方法
• 2.1.1 Hartree-Fock从头计算方法
• 2.1.2 密度泛函理论
• 2.1.3 交换-关联泛函的显式表达式
• 2.1.4 相对论效应、有效芯势
• 2.1.5 量子化学与材料科学中的计算软件包
• 2.2 结构的优化与相关性质计算
• 2.2.1 几何结构的优化与结构相关计算
• 2.2.2 能量的相关计算
• 2.2.3 电荷分布的相关计算
• 第三章 3d过渡金属双原子掺杂双棱柱硅团簇的理论研究
• 3.1 3d过渡金属双原子掺杂双棱柱硅团簇的研究进展
• 3.2 计算方法
• 2Si₁₈（M=Ti,...,Zn）的理论研究'>3.3 同种3d过渡金属双原子掺杂双棱柱硅团簇M₂Si₁₈（M=Ti,...,Zn）的理论研究
• 3.3.1 几何结构
• 3.3.2 电荷分布
• 3.3.3 稳定性
• 18（M,M’=Ti,...,Zn）理论研究'>3.4 异种3d过渡金属双原了混合掺杂双棱柱硅团簇MM’Si₁₈（M,M’=Ti,...,Zn）理论研究
• 3.4.1 几何结构
• 3.4.2 电荷分布
• 3.4.3 稳定性
• 3.5 小结
• mSi_n（m=1-2,n,=1-18）的理论研究'>第四章 钨掺杂硅团簇W_mSi_n（m=1-2,n,=1-18）的理论研究
• 4.1 钨掺杂硅团簇的研究进展
• 4.2 计算方法
• n（n=1-12）的计算结果'>4.3 单个钨原子掺杂硅团簇WSi_n（n=1-12）的计算结果
• 4.3.1 几何结构
• 4.3.2 电荷分布
• 4.3.3 稳定性
• 2Si_n（n=7-18）的计算结果'>4.4 两个钨原子掺杂硅团簇W₂Si_n（n=7-18）的计算结果
• 4.4.1 几何结构
• 4.4.2 电荷分布
• 4.4.3 稳定性
• 4.5 小结
• 第五章 结论与展望
• 参考文献
• 攻读博士学位期间的研究成果及项目
• 致谢