Fe掺杂CdTe电磁性质的第一性原理研究
铁电体的第一性原理研究进展PPT课件
国内外研究现状
Cohen首次采用第一性原理计算BaTIO3和PbTIO3的电子密度,表明Ti的3d电子 和O的2P电子波函数有显著的交叠,而且铁电情况下的交叠比顺电情况更强, 进而得出结论:对于ABO3结构的钙钦矿型铁电体,B离子的d电子与氧离子的 2p电子之间存在比较强的轨道杂化,这种轨道杂化抑制了短程排斥力从而使 铁电性得以稳定。
4.集成铁电体的研究(铁电薄膜与半导体集成):
由于铁电存储器的诸多优点,近几年来人们对铁电薄膜与半导体集成投入了大量的研究。 铁电薄膜的极化具备两个不同的稳定状态(剩余极化强度士Pr),可分别作为信息存储的“0‘’ 和,‘l”代码。早在50年代人们就开始作了这方面的研究。当时存在的问题主要为:块材要求 电压很高,不能满足应用的要求;电滞回线的矩形度差,易发生读写错误;疲劳特性很差。80年 代以来,由于铁电薄膜制备技术的改进,新的铁电材料及电极材料的出现,铁电存储器又重新 活跃起来。
2.Gaussian98程序包。Gaussian98程序包中包含许多种计算方法,包括半经验及第一性 原理计算方法等。它是一个功能全面的计算程序包。它的主要处理对象是有机大分子体系, 计算时主要对单一大分子体系的各种性质进行计算。能给出有机分子的振动模式及反应过 程的信息。它的缺点是对含有重金属原子体系的计算几乎无法进行。
研究热点
尺寸效应和表面界面效应
金属或半导体电极间的铁电薄膜 铁电聚合物和复合材料的研究
写在最后
成功的基础在于好的学习习惯
The foundation of success lies in good habits
17
谢谢聆听
·学习就是为了达到一定目的而努力去干, 是为一个目标去 战胜各种困难的过程,这个过程会充满压力、痛苦和挫折
Fe原子薄片的磁性_第一性原理计算
047502-2
物 理 学 报
Acta Phys. Sin.
Vol. 60 ,No. 4 ( 2011 )
047502
作用的是 Fe 的 3d 电子 .
表2 s 轨道 电子占据数 / e 平面正方 平面六角 体心立方 0. 502 0. 523 0. 487 磁矩 / μ В 0. 025 0. 020 - 0. 014 三种结构在平衡晶格常数时的价电子占据数以及对磁矩的贡献 p 轨道 电子占据数 / e 0. 305 0. 367 0. 538 磁矩 / μ В - 0. 009 0. 008 - 0. 056 d 轨道 原子总的 磁矩 / μ В 2. 651 2. 536 2. 199
* 国家自然科学基金( 批准号:10774124 ) 资助的课题 . 通讯联系人 . E-mail : zzhu@ xmu. edu. cn
, 选取 了 15 × 15 × 1 的 k 网 格 . 计 算 时 使 用 了
超原胞和周期性边界条件, 单 原 子 层 薄 片 方 向 为 x, y 方向, z 方向取 为 15 厚 的 真 空 层, 以消除原子薄
( 2010 年 6 月日 13 收到;2010 年 7 月 27 日收到修改稿)
使用基于密度泛函理论的第一原理方法, 对 Fe 单层原子薄片在二维正方 、 二维 六 角 晶 格 下 的 电 子 结 构 和 磁 学 性质进行了系统研究 . 结果表明, 二维正方 、 二 维 六 角 以 及 bcc 晶 格 在 平 衡 晶 格 常 数 下 都 具 有 磁 性, 其单位原子磁 2. 54 和 2. 20 μ В . 对二维晶格在被压缩和被拉伸时的磁性计算表明, 随着 晶 格 的 被 拉 伸, 当最近邻原 矩分别为 2. 65 , 铁原子间的键合被拉断, 体系单位原子的磁矩趋于孤立 Fe 原子的磁矩 4 μ В ;随着原子键长的 子间距大于 4. 40 时, 减小, 各体系的磁矩都随着最紧邻原子间距的减 小 而 减 小 . 当 键 长 缩 短 到 一 定 的 临 界 值 时 ( 平 面 正 方 1. 80 , 平面 六角 1. 75 ) , 铁磁性都会消失 . 使用 Stoner 理论, 可以理解晶格被缩短时体系由磁性到非磁性的变化 .
Fe掺杂碳纳米管吸附Al原子的第一性原理研究-中国有色金属学报
利用第一性原理对
基金项目:国家自然科学基金资助项目(51172062,51472074);河北省引进海外高层次人才“百人计划”资助项目(E2012100005) 收稿日期:2016-05-26;修订日期:2016-11-20 通信作者:安立宝,教授,博士;电话:0315-8805440;E-mail: lan@
由于具有理想的一维结构和独特的物理、化学性 质
[1−3]
碳纳米管与 Al 接触界面进行了理论计算, 结果表明碳 纳米管与 Al 之间有很高的肖特基势垒。GAO 等[15]研 究了 Pd 和 Al 与碳纳米管接触界面的电子结构和接触 电阻, 结果表明: 提高 Al 与碳纳米管之间的杂化作用 能降低它们之间的接触电阻。在碳纳米管表面掺杂的 研究中,B、N、Si 是经常采用的掺杂元素,而掺杂过 渡金属的研究却很少, 特别是 Fe 具有与碳纳米管相近 的功函数,并且与碳纳米管的接触电阻相对较低[13]。 因此, 本文作者采用第一性原理对本征和掺杂 Fe 碳纳 米管对 Al 的吸附进行模拟,研究掺杂 Fe 原子前后碳 纳米管对 Al 原子的吸附能,以及本征和掺 Fe 碳纳米 管吸附 Al 原子前后的能带结构和电荷分布, 从而探究 掺杂 Fe 对 Al 在碳纳米管表面吸附性能的影响。
, 碳纳米管被认为是制备高性能复合材料的理想
目前, 国内外已经开展了包括多种金属基 增强体[4−5]。 碳纳米管复合材料的研究, 其中 Al 基碳纳米管复合材 料具有耐腐蚀、低密度、易加工等优良性能,具有非 常广阔的应用前景
[6−7]
。但已有研究表明,Al 金属基[9]体与碳纳米管之间因浸润性差而存在吸附强度低、界 面接触电阻高等问题 。ESPEJEL-MORALES 等 研 究发现碳纳米管表面掺杂 B 原子后形成了缺电子结 构,并减小了碳纳米管的能隙。ZHANG 等
锰酸锂掺杂Fe和Co的性质的第一性原理计算
研究为锂电池电容量研究提供理论数据的参考
关键词:LiMn2O4尖晶石氧化物;第一性原理计算;Jahn - Teller 效应;Fe/ Co 掺杂;电压
中图分类号: 文献标识码: : TM912
A DOI 2."274889:"2..../;0"-.-/".2;..;
First - principles calculations for properties of lithium manganate doped with transition metals Fe and Co
过渡金属及氧化物团簇结构、磁性、催化性质的第一性原理计算研究
magnetic moment remarkably increases from 7烛of Scl30 to 131xB of Scl30CO,whereas it reduces
from 199n of Scl3 to 51xs ofScl3CO.
Eley—Rideal㈣and The competition between the
chemical reaction.
II
东南大学博士论文
东南大学博士论文
The performance of Pt/Cu(111)surface alloy catalysts in water gas shift reaction with a few Pt atoms doping in the Cu(111)surface were explored by using spin—polarized DFT.The Pt/Cu(111) surface alloy catalysts not only remarkably decrease the adsorption energy of CO,which Can
sequence in the presence ofCO.
The electronic properties and controlled single atom catalysis of isomorphous substituted bimetallic
oxide TMV3010(TM=Sc,Ti,Cr and Co)and V4010 were studied within the framework of a gradient-corrected DFT calculations.The charged clusters V409。1∥一show different electronic property from isoelectronic subsfitution bimetallic oxide Tiv309J0/CrV309,10 clusters and the terminal TM
Fe原子薄片的磁性_第一性原理计算
关键词 : Fe ,原子薄片,磁性,从头计算
PACS : 75. 50. Bb ,75. 70. Ak ,71. 15. Nc ,73. 20. - r
研 究 了 厚 度 为 一 个 原 子 层 的 Fe 一原理平面波法,
1. 引
言
薄片在平面正方和平面六角两种晶格下的电子 结构和磁性. 计 算 了 不 同 晶 格 常 数 下 ( 即 晶 格 被 压缩和拉伸 情 况 下) 的 磁 矩 和 电 子 状 态 密 度, 试 图 通 过 Stoner 理 论 以 及 对 能 带 结 构 和 电 子 态 密 讨 论 过 渡 金 属 元 素 Fe 在 二 维 情 况 下 度的分析, 的纳米磁性.
047502-2
物 理 学 报
Acta Phys. Sin.
Vol. 60 ,No. 4 ( 2011 )
047502
作用的是 Fe 的 3d 电子 .
表2 s 轨道 电子占据数 / e 平面正方 平面六角 体心立方 0. 502 0. 523 0. 487 磁矩 / μ В 0. 025 0. 020 - 0. 014 三种结构在平衡晶格常数时的价电子占据数以及对磁矩的贡献 p 轨道 电子占据数 / e 0. 305 0. 367 0. 538 磁矩 / μ В - 0. 009 0. 008 - 0. 056 d 轨道 原子总的 磁矩 / μ В 2. 651 2. 536 2. 199
电子占据数 / e 磁矩( μ В / atom ) 电子占据数 / e 6. 263 6. 223 6. 251 2. 634 2. 508 2. 269 7. 069 7. 113 7. 277
使得判据 I × D ( E F , 费米面上的电子态密度足够大, NM ) > 1 得 到 满 足 时, 体 系 将 出 现 磁 的 不 稳 定 性, I 为 Stoner 因 子, 即 会 出 现 自 发 的 磁 化 . 这 里, D( E F , NM ) 为非 自 旋 极 化 情 况 下 得 到 的 费 米 面 上 NM ) 由大于 1 变 的电子态密度 . 反之, 当 I × D( E F , 体系将由磁化状态变为非磁性 . 图 4 以 为小于 1 时, NM ) 随 平面六角晶格 为 例, 给 出 了 参 数 I × D( E F , 平面六角结构的最 紧 邻 原 子 间 距 的 变 化 曲 线 . 可 以 看到, 相应 体 系 的 磁 性 变 化, 即键长缩短时体系由 可以通过 Stoner 理论来理解 . 磁性到非磁性的变化,
金属元素掺杂α-Fe(N)体系的电子结构及力学性能的第一性原理计算
金属元素掺杂α-Fe(N)体系的电子结构及力学性能的第一性原理计算刘香军;杨吉春;贾桂霄;杨昌桥;蔡长焜【摘要】基于第一性原理研究M(M=Ti,V,Cr,Mn,Co和Ni)掺杂α Fe(N)的结合能、电子结构及力学性能.计算结果表明,Ti和V优先占据晶胞的顶角位置,Cr和Mn优先占据晶胞的体心位置,Co和Ni与N不相邻时结构最稳定.Ti与V的掺杂加强了晶胞的稳定性,Cr,Mn与Ni的掺杂削弱了晶胞的稳定性,Co的掺杂不影响晶胞的稳定性.这些过渡金属在α-Fe晶胞中均存在金属键和离子键的共同作用,成键轨道主要来自M3d,Fe4s3p3d与N2p.与纯α-Fe体系相比,掺杂体系刚性均变强,经计算可得α-Fe(N)-V体系的弹性模量E、剪切模量G和体积模量B均为最大值,即掺杂V可显著提高材料的力学性能,V是最有效的固氮元素,与高氮钢冶炼的实验结果相吻合.【期刊名称】《材料工程》【年(卷),期】2019(047)009【总页数】6页(P72-77)【关键词】第一性原理;高氮钢;电子结构;弹性常数【作者】刘香军;杨吉春;贾桂霄;杨昌桥;蔡长焜【作者单位】内蒙古科技大学材料与冶金学院,内蒙古包头014010;内蒙古科技大学材料与冶金学院,内蒙古包头014010;内蒙古科技大学材料与冶金学院,内蒙古包头014010;内蒙古科技大学工业技术研究院,内蒙古包头014010;内蒙古科技大学材料与冶金学院,内蒙古包头014010【正文语种】中文【中图分类】TG146N对钢材性能有着显著的影响,如强度、硬度均随N含量的增加而提高,韧性却不降低[1]。
同时N还能提升高氮钢的抗摩擦、抗腐蚀能力[2-4]。
由于常压下N 在钢液中的溶解度很低[5],通常采用高压充氮[6]和添加固溶金属氮化物的方法[7-9]来提高N在钢中的溶解度。
但高压充氮对设备的要求高、成本大,因此添加固溶金属氮化物(如Cr, Mn, Mo和V等的氮化物)来提高钢中氮含量是一种经济合理的有效方法。
电子和空穴掺杂CdO热电材料的第一性原理研究
电子和空穴掺杂CdO热电材料的第一性原理研究作者:张欣然来源:《山东工业技术》2015年第09期摘要:采用密度泛函理论下的第一性原理计算方法结合波尔兹曼输运理论,讨论不同电子和空穴掺杂浓度对CdO材料的电子结构及其热电输运特性的影响。
通过能带结构以及态密度的计算,分析掺杂后CdO的导带底和价带顶在费米能级附近的变化情况;在高温区电子和空穴掺杂CdO的赛贝克系数的绝对值都是随温度的升高而升高。
空穴掺杂后的载流子浓度比电子掺杂大,且空穴掺杂对CdO赛贝克系数的提高比电子掺杂要多。
关键词:第一性原理方法;热电材料;电子结构;赛贝克系数热电材料是将热能与电能直接相互转换的一种环境友好型功能材料,在温差发电和制冷领域有着广阔的应用潜力,寻找并开发高性能的热电转换材料是当务之急。
金属氧化物陶瓷材料引起了研究人员极大的兴趣,被认为是最具潜力的热电材料之一[1]。
氧化镉(CdO)是宽带隙的金属氧化物半导体材料,存在的氧空位和镉间隙缺陷使其具有较高的载流子浓度,具有优异电学性能,而被应用于平板显示装置以及光电子装置等领域[2]。
本文中,利用波尔兹曼输运理论计算不同电子和空穴掺杂浓度的CdO赛贝克系数随温度的变化情况,这为理解CdO的热电输运性质及设计CdO基的热电材料提供了理论依据[3]。
1 计算方法CdO属于面心立方结构,其空间群为Fm-3m,晶格常数为a=4.69Å。
计算中采用基于密度泛函理下的第一性原理软件包(Wien2k)和基于波尔兹曼输运理论的基础上的BoltzTrap程序软件包,电子之间的交换关联势采用Tran和Blaha修正的Becke-Johnson(TB-MBJ)泛函,以进一步提高计算的精确度。
对第一布里渊区采用54×54×54网格K点积分求和,进行自洽电子结构计算,停止后达到收敛[4]。
2 CdO的电子结构文中利用虚晶近似方法来模拟体系的不同电子和空穴掺杂浓度。
图1(a-b)分别为采用广义梯度近似下的PBE关联泛函和TB-MBJ近似方法计算的本征CdO分波态密度图。
铁电体的第一性原理研究进展
详细描述
铁电体的第一性原理研究具有重要的意义,它可以帮 助我们深入理解铁电体的电子结构、极化机制和物理 性能。通过第一性原理计算,我们可以获得铁电体的 能带结构、态密度、电荷分布等重要信息,从而更好 地理解其铁电效应的来源和机制。此外,第一性原理 研究还可以为新型铁电材料的设计和开发提供理论支 持,促进铁电材料在电子器件、能源存储、催化等领 域的应用。
详细描述
密度泛函理论通过将多电子系统的薛定谔方程近似为单电子近似,将复杂的电 子相互作用简化为密度与能量的关系,从而能够准确地计算分子的电子结构和 物理性质。
广义梯度近似
总结词
广义梯度近似是一种用于计算固体材 料的电子结构和物理性质的方法。
详细描述
广义梯度近似考虑了电子密度的空间 变化,能够更准确地描述电子的相互 作用和能量差分关系,相对于局域密 度近似具有更高的精度和适用范围。
详细描述
第一性原理计算是一种基于量子力学原理的计算方法,它通过求解薛定谔方程来获得原子尺度的物理 性质和化学反应过程。这种方法不需要实验参数,只需要元素的种类和原子间的相互作用力即可进行 计算。第一性原理计算在材料科学、化学、物理学等领域中有着广泛的应用。
研究目的与意义
总结词
铁电体的第一性原理研究旨在深入理解其电子结构、 极化机制和物理性能,为新型铁电材料的设计和开发 提供理论支持。
02
深入研究铁电材料的微观结构和物理机制,揭示其本质规律和
性能调控机制。
开展跨学科合作,整合不同领域的资源和优势,共同推动铁电
03
材料的研究和应用。
THANKS
感谢观看
详细描述
通过第一性原理计算,可以精确地描述铁电体的电子结构,包括能级分布、电子 态密度等,从而揭示铁电体的本征性质和内在机制。
合金元素掺杂对γ-Fe磁性影响的第一性原理研究
合金元素掺杂对γ-Fe磁性影响的第一性原理研究芦博俊;陈宏涛【摘要】本文采用第一性原理计算的方法研究了合金元素掺杂对γ-Fe磁性的影响.四种磁性γ-Fe的稳定性排序为:AFMD> AFM> NM> FM.N/C在γ-Fe中的固溶会导致γ-Fe由AFMD态向FM态转变,即N/C间隙原子导致了γ-Fe晶体内局域形成铁磁性结构.置换合金元素掺杂影响了几种磁性结构γ-Fe (M)的相对能量和稳定性,较高浓度合金元素掺杂会导致γ-Fe (M)稳定的磁性结构的改变或局域内会存在亚稳定磁性结构.【期刊名称】《热处理技术与装备》【年(卷),期】2016(037)005【总页数】4页(P5-8)【关键词】γ-Fe;合金元素;磁性;第一性原理【作者】芦博俊;陈宏涛【作者单位】哈尔滨理工大学材料科学与工程学院,黑龙江哈尔滨150080;哈尔滨理工大学材料科学与工程学院,黑龙江哈尔滨150080【正文语种】中文【中图分类】TG156.9钢铁材料在高温条件下为γ相,冷却到室温时一般为α相。
磁性在钢铁材料的相稳定性上的重要性已被世人所认可,但很多理论研究只是针对特定的因瓦铁-镍合金提出的[1-2]。
C、N、Cr、Ni、Mn等合金元素的掺杂使得钢铁材料在室温条件下保持为γ相。
γ-Fe具有面心立方的晶体结构,其磁性结构主要有:非铁磁(non ferromagnetic,NM)、铁磁(ferromagnetic,FM)、反铁磁(antiferromagnetic,AFM)、双反铁磁(double antiferromagnetic,AFMD)和顺磁态(paramagnetic,PM)[3]。
奥氏体相通常为顺磁性,但热处理工艺和形变将会对其磁性态产生影响。
γ-Fe最稳定的磁性状态与合金的微观组织结构密切相关。
合金元素的掺杂对γ-Fe的磁性状态有较大的影响。
由于钢铁材料复杂的磁性状态,现有的分析测试技术很难对其进行表征。
俄罗斯的Boukhvalov等人采用第一性原理计算的方法研究了C掺杂对γ-Fe磁性态和局部形变的影响[3]。
Fe掺杂锐钛矿相TiO2的第一性原理计算
Fe掺杂锐钛矿相TiO2的第一性原理计算作者:段霄斌来源:《山东工业技术》2018年第15期摘要:本文采用基于密度泛函理论的第一性原理,对未掺杂、Fe掺杂锐钛矿相TiO2的电子能带结构、电子态密度、电荷布居、光吸收谱进行理论计算,定性分析Fe掺杂前后对TiO2电子结构和光学性质的影响。
结果表明: Fe掺杂锐钛矿相TiO2,会导致禁带宽度减小,掺杂后在可见光波长范围内的的光吸收系数增大,TiO2的吸收带产生红移,增强了TiO2的光催化活性,为光催化剂的选择与制备提供了理论依据。
关键词:密度泛函理论;第一性原理;掺杂;电子结构;光学性质DOI:10.16640/ki.37-1222/t.2018.15.166近些年来,锐钛矿相TiO2不仅在实验上而且在理论上都得到了深入广泛的研究,这是因为TiO2相对于其他半导体氧化物,它的光催化能力较强。
但是锐钛矿相TiO2的太阳能利用率偏低,其光催化活性仍需进一步提高。
上述科学问题的解决仍然是目前太阳能利用领域的研究热点[1-3]。
国际上对掺杂TiO2进行理论研究也有大量报道。
Wu等[4]利用超声辅助法制备了Fe掺杂的TiO2纳米管,发现Fe掺杂可以有效增强TiO2纳米管的可见光吸收能力,可以高效地处理降解各种废水。
Hou等[5]利用简易的溶剂加热法制备了Fe/Er共掺杂的TiO2,发现样品在可见光照射下具有良好的光催化活性,并对相关机理进行了阐释。
1997年,Luque等[6]首次提出了中间带太阳能电池的概念,即在半导体禁带之内引入一个电子半填满的中间能带。
在小于禁带宽度光子能量照射下,光子可激发电子从价带跃迁到中间带,同时光子还可激发电子从中间带到导带。
这两个光子激发电子的过程,是对半导体材料在能量大于禁带宽度Eg的光子激发电子产生本征跃迁的有效补充,这可以充分扩展半导体材料的光响应范围,提高半导体材料对太阳光子的利用率,能够提高材料的光电效率。
详细了解和深入认识材料的电子结构和光学性质对基础理论研究及新材料的研制与开发均有指导意义,因此本文将计算Fe掺杂TiO2电子结构和光学性质。
基于第一性原理的Fe(111)
第15卷第2期2024年4月有色金属科学与工程Nonferrous Metals Science and EngineeringVol.15,No.2Apr. 2024基于第一性原理的Fe (111)/Al 2O 3(0001)界面稳定性研究王浦璠, 邓道繁, 陈果, 张慧宁, 廖春发*(江西理工大学材料冶金化学学部,江西 赣州 341000)摘要:稀贵金属催化剂催化甲烷分解效果显著但成本较高,为此,探寻低成本催化剂是实现其工业化推广的重要课题。
铁基催化剂具有耐高温、价格低、性能好等优点,目前针对铁基催化剂的表面机理研究较少,催化剂界面组成及结构稳定性仍没有统一的观点。
本研究基于密度泛函理论(DFT ),利用第一性原理计算表面能、态密度、差分电荷密度,研究Fe (111)/Al 2O 3(0001)界面稳定性,并探讨不同端面的催化剂形成机制。
结果表明:O 端界面的吸附能和电荷浓度大于Al1和Al2端界面,O 原子周围电荷聚集明显;Al1和Al2端界面的Fe 和Al 之间形成微弱金属键,而O 端界面的Fe 和Al 之间形成共价键、金属键和离子键,O 端共价键强度高于Al1和Al2端界面。
因此,O 端Fe/Al 2O 3催化剂相比于其他端面结合更加紧密,同时表现出更理想的催化效果,上述研究结果为制备Fe/Al 2O 3催化剂及研究其界面的形成机制提供理论基础。
关键词:密度泛函理论;表面能;金属催化剂;态密度;差分电荷密度中图分类号:TG113.12 文献标志码:AStability of Fe (111)/Al 2O 3(0001) interface based on first principles WANG Pufan, DENG Daofan, CHEN Guo, ZHANG Huining, LIAO Chunfa *(Faculty of Materials Metallurgy and Chemistry , Jiangxi University of Science and Technology , Ganzhou 341000, Jiangxi , China )Abtsract :Rare and precious metal catalysts have remarkable effects on methane decomposition, but the cost is high, so that it is an important task to explore low-cost catalysts to realize their industrialization and popularization. Iron-based catalysts have the advantages of high-temperature resistance, low price and good performance. At present, there is little research on the surface mechanism of iron-based catalysts, and there is still no unified view of catalyst interface composition and structural stability. In this study, based on density functional theory (DFT), the stability of the Fe (111) /Al 2O 3 (0001) interface was studied by first principles, and the formation mechanism of the catalyst with different end faces was discussed. The results of surface energy, state density and differential charge density analysis showed that the adsorption energy and charge concentration at the O-terminal interface were greater than those at the interface of Al1 and Al2, and the charge accumulation around the O atom was obvious. Weak metal bonds were formed between Fe and Al at the interface of Al1 and Al2, while covalent, metallic and ionic bonds were formed between Fe and Al at the interface of the O terminal. The O terminal covalent bond strength was also higher than the Al1 and Al2 terminal interfaces. Therefore, the O-end Fe/Al 2O 3 catalyst was more tightly bound than other end faces, and would also show a better catalytic effect. The above research results provide a theoretical basis for收稿日期:2023-09-28;修回日期:2023-11-23基金项目:国家自然科学基金地区基金资助项目(52064020)通信作者:廖春发(1965— ),博士,教授,主要从事稀土、钨等稀贵金属冶金方面的研究。
基于第一性原理的Au、Cu、Sb掺杂CdTe的结构模拟和光学性能预测
第52卷第8期2023年8月人㊀工㊀晶㊀体㊀学㊀报JOURNAL OF SYNTHETIC CRYSTALS Vol.52㊀No.8August,2023基于第一性原理的Au ㊁Cu ㊁Sb 掺杂CdTe 的结构模拟和光学性能预测张平威1,林㊀龙2,张战营2(1.中建材(洛阳)新能源有限公司,洛阳㊀471000;2.河南理工大学材料科学与工程学院,焦作㊀454000)摘要:基于密度泛函理论研究了Au㊁Cu㊁Sb 掺杂CdTe 体系的电子结构和光学性能㊂Au㊁Cu㊁Sb 掺杂CdTe 体系均能稳定存在,过渡金属原子与Cd 原子轨道的杂化减小了CdTe 的带隙,提高了CdTe 对可见光的利用,同时降低了光生电子从价带跃迁到导带所需的能量,从而促进了更多的光生电子发生迁移,大大提高了其光学性能㊂三种掺杂体系中Sb /CdTe 体系在可见光范围内光吸收系数提升最显著,其光生电子和光生空穴迁移率相对于CdTe 体系分别增加5.97倍和15.54倍㊂通过计算掺杂体系的能带㊁态密度㊁电子布居㊁光吸收函数㊁载流子迁移率,从理论上揭示了Au㊁Cu㊁Sb提高CdTe 光学性能的机理㊂关键词:CdTe;第一性原理;电子性质;光学性能;载流子迁移中图分类号:O439;O469;TQ132.4+4㊀㊀文献标志码:A ㊀㊀文章编号:1000-985X (2023)08-1400-07Structural Simulation and Optical Property Prediction of Au ,Cu ,and Sb Doped CdTe Based on First-PrinciplesZHANG Pingwei 1,LIN Long 2,ZHANG Zhanying 2(1.China National Buiding Materials (Luoyang)New Energy Resources,Luoyang 471000,China;2.School of Materials Science and Engineering,Henan Polytechnic University,Jiaozuo 454000,China)Abstract :The electronic and optical properties of Au,Cu,Sb doped CdTe systems were studied based on density functional theory.Au,Cu and Sb doped CdTe systems all exist stably.The hybridization of transition metal atoms with Cd atomic orbitals reduces the band gap of CdTe and improves the utilization of visible light by CdTe.The lower energy required to jump from the valence band to the conduction band promotes the migration of more photogenerated electrons,which greatly improves the optical properties of doped CdTe.Among the three systems,Sb /CdTe system shows the most significant increase of absorption coefficient in the visible light range,with photogenerated electron and hole mobilities increasing by a factor of 5.97times and15.54times compared with CdTe system,respectively.The mechanism of the enhancement of the optical properties of Au,Cu,and Sb doped CdTe is theoretically revealed by calculating the band,density of states,electron population,optical absorption function,and carrier mobility.Key words :CdTe;first-principle;electronic property;optical property;carrier mobility ㊀㊀收稿日期:2023-02-13㊀㊀基金项目:国家自然科学基金联合基金重点项目(U2004209)㊀㊀作者简介:张平威(1968 ),男,河南省人,工程师㊂E-mail:lyblgfyxgs@㊀㊀通信作者:林㊀龙,博士,教授㊂E-mail:linlonghpu@ 0㊀引㊀㊀言近年来,CdTe 材料由于其高光吸收系数㊁窄带隙㊁良好的电荷传输等性能而受到了广泛的关注[1-3]㊂CdTe 在现代电子领域有着广泛的应用,如太阳能电池㊁辐射探测器㊁光电探测器等[4-7]㊂CdTe 良好的光学性能使其在太阳能电池的研发中占据了重要地位㊂目前,CdTe 基太阳能电池的转换效率已达到22.1%,但距离其理论效率27%~30%仍有较大改进空间[8-10]㊂这主要也是由一些固有缺陷引起,如光激发电子和空穴㊀第8期张平威等:基于第一性原理的Au㊁Cu㊁Sb掺杂CdTe的结构模拟和光学性能预测1401㊀的快速复合,大大降低了其光激发效率,带隙较宽而对可见光利用率较低,以及光化学稳定性不足引起材料的失活等[11-12]㊂因此提升CdTe的光学性能对于其在现代光电子领域的应用具有重要意义㊂大量的文献报道CdTe材料中的固有缺陷和外在缺陷可提高其性能,如空位及外来掺杂原子等可以影响材料能量结构进而提升其光学性能[13-15]㊂Sedzicki等[16]报道了P掺杂CdTe促进了材料吸收边的红移进而提升了其光学性能㊂Ding等[17]发现Ag掺杂CdTe后,其激子共振吸收增强,在带边附近的非线性折射增强了35倍以上㊂Alzaid等[18]研究了Cu掺杂含量对CdTe的影响,结果发现其消光系数随着Cu含量的增加而减小,而电导率则随着Cu含量的增加逐渐增大,进而提高了太阳能电池结构中CdTe吸收层的性能㊂Li 等[19]发现Au掺杂CdTe后,其晶体间的强偶级相互作用促进了发光峰的红移㊂Arivarasan等[20]报道了Sb 掺杂后CdTe量子点敏化的TiO2光阳极表现出比纯CdTe量子点更好的光伏响应㊂因此,利用过渡金属掺杂的方式来提高CdTe的光学性能是可行的㊂但目前为止,虽然实验上大量报道了过渡金属原子掺杂改变了CdTe光学性能,但其性能提高的机理却很少有人报道㊂本文通过密度泛函理论计算,研究Au㊁Cu㊁Sb掺杂CdTe(M/CdTe)的电子结构和光学性能,包括能带结构㊁态密度㊁光吸收函数㊁载流子迁移率等,系统地阐明其光电性能增强的机理㊂1㊀模型与方法在这项工作中,所有密度泛函理论的计算都在从头计算的模拟包VASP中进行[21-22]㊂为了描述电子和离子核之间的相互作用,采用了投影增广波势的方法㊂采用广义梯度近似(GGA)方法处理电子之间的交换相互作用[23-24]㊂所有计算均采用截止能量为400eV的平面波基进行㊂当每个原子受到的最大力小于0.02eV/Å,两个连续迭代步骤的总能量差小于10-4eV时,认为结构弛豫和自洽环实现了收敛㊂采用了3ˑ2ˑ1的以G 点为中心的网格对结构进行了几何优化㊂模拟了2ˑ2ˑ1的超胞,显示在图1中,其包含了32个原子,包括16个Cd原子,显示为米白色;16个Te原子,显示为橙色;红色位置为选择的过渡金属(Au㊁Cu㊁Sb)掺杂位置,其替代了Cd原子位置㊂为了比较不同掺杂体系的稳定性,定义形成能[25-26],公式为E f=E M/CdTe-E I-E M+E Cd(1)式中:E M/CdTe为金属原子掺杂CdTe的总能量,E I为本征CdTe的总能量,E M和E Cd分别为单个金属原子和单个Cd原子的能量㊂形成能越小表示结构越稳定㊂图1㊀2ˑ2ˑ1的过渡金属原子掺杂CdTe超胞结构Fig.1㊀2ˑ2ˑ1CdTe supercell with doped transition metal atoms图2为Au/CdTe㊁Cu/CdTe㊁Sb/CdTe体系的形成能和Bader电荷转移㊂从图中可以看出,Au掺杂CdTe(Au/CdTe)㊁Cu掺杂CdTe(Cu/CdTe)的形成能为负值,说明结构是稳定的,Sb掺杂CdTe(Sb/CdTe)的形成能为3.427eV,而在优化后的Sb/CdTe晶格结构并未发生明显变化,表明该结构同样是稳定的㊂为了进一步验证Sb/CdTe的稳定性,计算了该结构的声子谱和声子态密度,结果如图3所示㊂图3(a)显示Sb/CdTe构型声子散射曲线中不存在虚频率,说明所构建模型是合理的㊂此外,Sb/CdTe构型的总声子态密度结果绘制在图3(b)中,结果显示在声子频率处没有虚频率,这表明Sb/CdTe结构处于基态,其几何优化完全收敛㊂因此Sb/CdTe结构是稳定的㊂综上所述,三个研究的掺杂体系都是稳定可靠的㊂1402㊀研究论文人工晶体学报㊀㊀㊀㊀㊀㊀第52卷图2㊀Au /CdTe㊁Cu /CdTe㊁Sb /CdTe 体系的形成能和Bader 电荷转移Fig.2㊀Formation energy and Bader charge transfer of Au /CdTe,Cu /CdTe,Sb /CdTesystems 图3㊀Sb /CdTe 体系的声子色散曲线(a)和声子态密度曲线(b)Fig.3㊀Phonon dispersion curve (a)and phonon density of states curve (b)of Sb /CdTe system 2㊀结果与讨论2.1㊀几何结构优化首先利用GGA 的方法对本征CdTe 体系进行结构优化,计算得到的晶格常数a =4.72Å,b =8.07Å,c =7.64Å,与实验结果相似[27]㊂表1中计算的电子布居显示,在Au /CdTe 体系中,Au 原子作为受体从相邻的Te 原子中获得了0.249eV 电子,相反,Cu 和Sb 原子在掺杂体系中作为供体分别损失了0.216e 和0.426e㊂这些转移的电子促进了掺杂体系中的电子转移,也相应增强了掺杂原子与周围原子的结合而形成较强的键,表明所研究的结构是稳定的㊂2.2㊀能带理论为了定性阐明过渡金属掺杂后所引起的CdTe 中电子运动的特点,计算了CdTe㊁Au /CdTe㊁Cu /CdTe㊁Sb /CdTe 体系能带图,如图4所示㊂所计算的本征CdTe 的禁带宽度为0.613eV,小于实验值1.5eV [27],这是因为能带结构的计算方法利用GGA 近似,其并没有考虑交换相关性的不连续性,但其并不影响理论结果分析㊂而Au /CdTe㊁Cu /CdTe㊁Sb /CdTe 体系的带隙则明显减小,其值分别为0.394㊁0.504㊁0.477eV,且均在费米能级附近出现杂质能级㊂当电子被激活时,其从价带跳到导带时所需的能量会降低,而光的波长也会随着可见光响应的增强而增加,增加了CdTe 的光吸收系数㊂此外,Sb /CdTe 体系的带隙有向价带迁移,导带通过了费米能级,提高了CdTe 的半导体特性和氧化还原能力,其光催化能力得到了进一步的增强㊂2.3㊀态密度为了进一步了解Au /CdTe㊁Cu /CdTe㊁Sb /CdTe 体系光学性能增强的机理,从微观上了解其能带变化的原因,计算了CdTe㊁Au /CdTe㊁Cu /CdTe㊁Sb /CdTe 体系的态密度,包括各体系的总态密度和部分态密度,结果如图5所示㊂从图5(a)可以看出,本征CdTe 的价带主要由Te 原子贡献,而导带部分则是由Cd 原子和Te 原子共同贡献㊂图5(b)为Au /CdTe 的态密度,接近费米能级的Au /CdTe 的价带主要由Te 原子和Au 原子㊀第8期张平威等:基于第一性原理的Au㊁Cu㊁Sb 掺杂CdTe 的结构模拟和光学性能预测1403㊀贡献,两者轨道间形成了杂化,而在费米能级处出现的杂质带也减少了带隙,进而提高了其对可见光的利用㊂Cu /CdTe 体系有着和Au /CdTe 体系相似的变化㊂从图5(d)中可以看出,当Sb 原子引入到CdTe 后,在费米能级左侧价带处出现了杂质峰,进一步缩短了电子的跃迁路径,促进了载流子的迁移且增加了体系的电导率㊂此外,在此出现的杂质峰主要由Sb 原子贡献,而Cd 原子也作了一定的贡献,两者之间形成了杂化㊂这种杂化也使得该体系的导带和价带向较低的能量方向移动,而导带底部通过费米能级,也使其更容易产生光生电子,进一步增强其光学性能㊂图4㊀CdTe(a)㊁Au /CdTe(b)㊁Cu /CdTe(c)㊁Sb /CdTe(d)体系能带结构Fig.4㊀Energy band structures of CdTe (a),Au /CdTe (b),Cu /CdTe (c)and Sb /CdTe (d)system 图5㊀CdTe(a)㊁Au /CdTe(b)㊁Cu /CdTe(c)㊁Sb /CdTe(d)体系的态密度Fig.5㊀Density of states of CdTe (a),Au /CdTe (b),Cu /CdTe (c)and Sb /CdTe (d)system 2.4㊀光学性能为了更直观地评价Au㊁Cu㊁Sb 原子对CdTe 光学性能的影响,计算了各掺杂体系的光吸收系数,并将结果与本征CdTe 进行比较,结果如图6所示㊂从图中可以看出,在可见光范围内(1.63~3.10eV)[28-31],相比于本征CdTe 体系,Au /CdTe㊁Cu /CdTe㊁Sb /CdTe 体系的光吸收系数均明显增强,表明其吸收光子的能力得到了增强㊂这可能是由于这些体系能带降低,光吸收电子从价带跃迁到导带所需能量降低,进而促进了更多的1404㊀研究论文人工晶体学报㊀㊀㊀㊀㊀㊀第52卷图6㊀CdTe㊁Au /CdTe㊁Cu /CdTe㊁Sb /CdTe 体系光学吸收函数Fig.6㊀Optical absorption functions of CdTe,Au /CdTe,Cu /CdTe and Sb /CdTe systems 电子跃迁㊂其中Sb /CdTe 体系的光吸收系数最强,这是因为其产生的杂质带缩短了电子的跃迁路径,进一步促进了光生电子的跃迁,此结果也与计算的态密度和能带结构相符合㊂2.5㊀载流子迁移率材料是否具有优异的载流子迁移率在现代电子设备的应用中具有极其重要的意义㊂为了进一步了解各掺杂体系的电子性质,根据形变势理论计算了室温(T =300K)下界面的载流子迁移率(μ),其表达式为[32]μ=2e ħ3C 3k B T |m ∗|2E 2d (2)式中:e 为电子电荷,k B ㊁T 和ħ分别为玻尔兹曼常数㊁温度和还原普朗克常数,E d 定义为由应变引起的带边位移而引起的形变势常数,C 为系统的弹性模量㊂从价带顶和导带底的x 和y 方向的带结构中提取出载流子的有效质量㊂此外,在模型中分别沿x 和y 方向施加-1%到1%(每次应变为0.5%)的应变,使结构被压缩或膨胀㊂弹性模量C 通过拟合关于能量-应变关系的二次函数来计算[33]㊂之后,通过对导带底和价带顶的能量作不同应变的函数的线性拟合,得到形变势常数E d [34]㊂根据这些系数带入公式计算出μ,计算的各参数及载流子迁移率显示在表1中㊂结果显示,CdTe 的电子和空穴迁移率比较小,在x 方向为38.165和3.774cm 2㊃V -1㊃s -1,而y 方向同样如此,其电子和空穴迁移率为48.341和1.571cm 2㊃V -1㊃s -1㊂在经过过渡金属Au㊁Cu㊁Sb 掺杂CdTe 后,其电子和空穴的迁移率得到了很大的提升㊂如在y 方向上,Au /CdTe㊁Cu /CdTe㊁Sb /CdTe 体系的电子迁移率相对于CdTe 分别增长了16.38倍㊁17.22倍和5.97倍,这是由于过渡金属原子的掺入导致CdTe 发生了局部结构断裂及电子重新排列,进而使得整个材料的能量和带边位置发生了变化㊂而y 方向上,Sb /CdTe 体系的空穴迁移率更高,相对于CdTe 增长15.54倍,结合电子迁移率,Sb /CdTe 体系的电子和空穴的结合率会降低,其电子和空穴寿命也会相应增加㊂结合上述分析,Au /CdTe㊁Cu /CdTe㊁Sb /CdTe 体系均能增强CdTe 的光学性能,Sb /CdTe 体系因Sb 原子与Cd 原子在费米能级发生杂化而进一步降低载流子迁移的路径,增加了载流子迁移率,具有更好的光学性能㊂表1㊀CdTe ㊁Au /CdTe ㊁Cu /CdTe ㊁Sb /CdTe 体系的载流子有效质量(m ∗)㊁形变势常数(E d )㊁弹性模量(C )和载流子迁移率(μ)Table 1㊀Effective carrier mass (m ∗),deformation potential constant (E d ),modulus of elasticity (C )andcarrier mobility (μ)for CdTe ,Au /CdTe ,Cu /CdTe ,Sb /CdTe systemsMaterial Carrier m ∗/kg E d /eV C /(J㊃m -2)μ/(cm 2㊃V -1㊃s -1)x y x y x y x yCdTe h -0.557-0.478-3.360-6.20037.56639.150 3.7741.571CdTe e 0.0810.076-7.280-7.00037.56639.15038.16548.341Au /CdTe h -0.531-0.352-5.006-4.07432.76832.768 1.631 5.619Au /CdTe e 0.0240.017-7.372-7.07232.76832.768357.596791.996Cu /CdTe h -0.308-0.222-3.986-5.01133.42235.0577.80215.101Cu /CdTe e 0.0260.016-7.349-7.01033.42235.057322.653832.514Sb /CdTe h -0.104-0.120-6.339-6.23837.05736.23130.16524.565Sb /CdTe e 0.0290.034-7.638-6.01837.05736.231261.869288.7543㊀结㊀㊀论综上所述,通过第一性原理的方法研究了Au㊁Cu㊁Sb 掺杂CdTe 的电子结构和光学性能㊂光吸收系数结果显示,在可见光范围内,各掺杂体系的光吸收系数均得到了明显提高,且计算的载流子迁移率也表明各掺㊀第8期张平威等:基于第一性原理的Au㊁Cu㊁Sb掺杂CdTe的结构模拟和光学性能预测1405㊀杂体系光生电子的迁移率得到了明显提高㊂这主要是因为在过渡金属原子掺杂后,掺杂体系的带隙减小,降低了光生电子从价带跃迁到导带的能量,进而促进载流子迁移,增加表面电导率而增强其对光子的吸收㊂在三种掺杂体系中,Sb/CdTe体系具有最强的光学性能,这是因为Sb掺杂使得其在带隙处出现了杂质带,进一步减少了光生电子跃迁的路径,使得其转移了最多的电子(0.426e)㊂此外,Sb/CdTe体系光生空穴迁移率相对Au/CdTe㊁Cu/CdTe更高,这会降低Sb/CdTe电子和空穴的结合率,且能增加其转移电子和空穴的复合时间而进一步增强光学性能㊂因此,Sb/CdTe体系具有明显增强的光学性能,可以很好地应用于光伏材料㊂参考文献[1]㊀ARIVARASAN A,BHARATHI S,EZHIL ARASI S,et al.Investigations of rare earth doped CdTe QDs as sensitizers for quantum dots sensitizedsolar cells[J].Journal of Luminescence,2020,219:116881.[2]㊀YANG P Q,CHEN C W,WANG D F,et al.Kinetics,reaction pathways,and mechanism investigation for improved environmental remediationby0D/3D CdTe/Bi2WO6Z-scheme catalyst[J].Applied Catalysis B:Environmental,2021,285:119877.[3]㊀BAJOROWICZ B,NADOLNA J,LISOWSKI W,et al.The effects of bifunctional linker and reflux time on the surface properties andphotocatalytic activity of CdTe quantum dots decorated KTaO3composite photocatalysts[J].Applied Catalysis B:Environmental,2017,203: 452-464.[4]㊀SARACHO-GONZÁLEZ S,PÉREZ-CENTENO A,SANTANA-ARANDA M A,et al.Effect of the combination of Cu and CdTe plasmas on thestructural and optical properties of CdTe:Cu thin films deposited by laser ablation[J].Materials Science in Semiconductor Processing,2018, 87:7-12.[5]㊀ÇIRIŞA,BAŞOL B M,ATASOY Y,et al.Effect of CdS and CdSe pre-treatment on interdiffusion with CdTe in CdS/CdTe and CdSe/CdTeheterostructures[J].Materials Science in Semiconductor Processing,2021,128:105750.[6]㊀RAHMAN M F,HOSSAIN J,KUDDUS A,et al.A novel CdTe ink-assisted direct synthesis of CdTe thin films for the solution-processed CdTesolar cells[J].Journal of Materials Science,2020,55(18):7715-7730.[7]㊀ZHANG J Y,CAO H C,BAI W,et al.High-sensitivity CdTe phototransistors with the response spectrum extended to1.65μm[J].Journal ofMaterials Chemistry A,2022,10(39):20837-20846.[8]㊀MANIMOZHI T,LOGU T,ARCHANA J,et al.Enhanced photo-response of CdTe Thin film via Mo doping prepared using electron beamevaporation technique[J].Journal of Materials Science:Materials in Electronics,2020,31(23):21059-21072.[9]㊀KRASIKOV D,SANKIN I.Defect interactions and the role of complexes in the CdTe solar cell absorber[J].Journal of Materials Chemistry A,2017,5(7):3503-3513.[10]㊀BURST J M,DUENOW J N,ALBIN D S,et al.CdTe solar cells with open-circuit voltage breaking the1V barrier[J].Nature Energy,2016,1(3):1-8.[11]㊀GNATENKO Y,BUKIVSKIJ P M,GAMERNYK R V,et al.Photoluminescence of CdTe thin films doped with Yb[J].Journal ofLuminescence,2021,237:118208.[12]㊀PEDETTI S,ITHURRIA S,HEUCLIN H,et al.Type-II CdSe/CdTe core/crown semiconductor nanoplatelets[J].Journal of the AmericanChemical Society,2014,136(46):16430-16438.[13]㊀SAVCHUK A I,PARANCHYCH S Y,FRASUNYAK V M,et al.Optical and magnetooptical study of CdTe crystals doped with rare earth ions[J].Materials Science and Engineering:B,2003,105(1/2/3):161-164.[14]㊀LI D B,BISTA S S,SONG Z N,et al.Maximize CdTe solar cell performance through copper activation engineering[J].Nano Energy,2020,73:104835.[15]㊀STRONG V,URIBE-ROMO F J,BATTSON M,et al.Oriented polythiophene nanofibers grown from CdTe quantum dot surfaces[J].Small,2012,8(8):1191-1196,1125.[16]㊀SEDZICKI P,SKOWRONSKI L,SZCZESNY R,et al.Influence of phosphorus ion implantation on the optical properties of CdTe bulk crystal[J].Journal of Alloys and Compounds,2020,844:156002.[17]㊀DING S J,NAN F,LIU X N,et rgely enhanced optical nonlinear response of heavily doped AgʒCdTe nanocrystals around the excitonicband edge[J].The Journal of Physical Chemistry C,2015,119(44):24958-24964.[18]㊀ALZAID M,HADIA N M A,EL-HAGARY M,et al.Microstructural,optical,and electrical characteristics of Cu-doped CdTe nanocrystallinefilms for designing absorber layer in solar cell applications[J].Journal of Materials Science:Materials in Electronics,2021,32(11): 15095-15107.[19]㊀LI J,LI D,HONG X,et al.Unique structure and photoluminescence of Au/CdTe nanostructure materials[J].Chemical Communications(Cambridge,England),2004(8):982-983.[20]㊀ARIVARASAN A,BHARATHI S,ESSAKINAVEEN D,et al.Investigation on the Role of antimony in CdTe QDs sensitized solar cells[J].Optical Materials,2022,129:112551.1406㊀研究论文人工晶体学报㊀㊀㊀㊀㊀㊀第52卷[21]㊀LIN L,YAN L B,HE C Z,et al.A theoretical study of the ability of2D monolayer Au(111)to activate gas molecules[J].InternationalJournal of Hydrogen Energy,2021,46(21):11711-11720.[22]㊀JIN J M,CHEN J F,WANG H F,et al.Insight into room-temperature catalytic oxidation of NO by CrO2(110):a DFT study[J].ChineseChemical Letters,2019,30(3):618-623.[23]㊀LIU D D,SHI Y L,TAO L,et al.First-principles study of methanol adsorption on heteroatom-doped phosphorene[J].Chinese ChemicalLetters,2019,30(1):207-210.[24]㊀LIN L,CHEN R X,HUANG J T,et al.Adsorption of CO,H2S and CH4molecules on SnS2monolayer:a first-principles study[J].MolecularPhysics,2021,119:185429.[25]㊀CHEN R X,YAN L B,LIN L,et al.Coadsorption of CO and CH4on the Au doped SnO2(110)surface:a first principles investigation[J].Physica Scripta,2022,97(4):045403.[26]㊀LIN L,ZHU L H,ZHAO R Q,et al.First-principles study on ferromagnetism in4H-SiC codoped with Al and Mn[J].New Journal ofChemistry,2018,42(12):9393-9397.[27]㊀LANE D W.A review of the optical band gap of thin film CdS x Te1-x[J].Solar Energy Materials and Solar Cells,2006,90(9):1169-1175.[28]㊀LIN L,CHEN R X,HE C Z,et al.Magnetic and optical properties of(Mn,Co)co-doped SnO2[J].Vacuum,2020,182:109681.[29]㊀AUSTIN E,GEISLER A N,NGUYEN J,et al.Visible light.Part I:Properties and cutaneous effects of visible light[J].Journal of theAmerican Academy of Dermatology,2021,84(5):1219-1231.[30]㊀PITRE S P,YOON T P,SCAIANO J C.Titanium dioxide visible light photocatalysis:surface association enables photocatalysis with visible lightirradiation[J].Chemical Communications,2017,53(31):4335-4338.[31]㊀HEILER C,BASTIAN S,LEDERHOSE P,et al.Folding polymer chains with visible light[J].Chemical Communications,2018,54(28):3476-3479.[32]㊀YAO M L,WU T,LIU B,et al.First principle study on interfacial interaction of black phosphorus and CsBr vdW heterostructure[J].PhysicsLetters A,2020,384(25):126614.[33]㊀TANG X,LI S S,MA Y D,et al.Distorted Janus transition metal dichalcogenides:stable two-dimensional materials with sizable band gap andultrahigh carrier mobility[J].The Journal of Physical Chemistry C,2018,122(33):19153-19160.[34]㊀WANG Y J,SONG C Y,ZHANG H,et al.Study on the relationship between carrier mobility and nonlinear optical characteristics of Sb2Te3-Bi2Te3lateral heterostructure materials and its applications in fiber lasers[J].Journal of Materials Chemistry C,2022,10(33):11862-11873.。
Fe掺杂碳纳米管吸附Al原子的第一性原理研究
Fe掺杂碳纳米管吸附Al原子的第一性原理研究
代利峰;刘涛;龚亮;常春蕊;安立宝
【期刊名称】《中国有色金属学报》
【年(卷),期】2017(027)006
【摘要】采用第一性原理对Al原子在本征(5,5)型及Fe掺杂(5,5)型碳纳米管的吸附行为进行研究,研究掺杂Fe前后碳纳米管吸附Al的吸附能、能带结构和电荷分布情况.计算结果表明:掺杂Fe能够增强Al在(5,5)型碳纳米管表面的吸附能,提高Al与碳纳米管之间的电荷转移.这主要是因为掺杂Fe原子后,Fe原子与Al原子之间形成共价键结合,Al与碳纳米管之间由物理吸附转变为化学吸附.此类掺杂有望改善Al基碳纳米管复合材料界面的电接触性能.
【总页数】7页(P1182-1188)
【作者】代利峰;刘涛;龚亮;常春蕊;安立宝
【作者单位】华北理工大学机械工程学院,唐山 063009;华北理工大学理学院,唐山063009;华北理工大学机械工程学院,唐山 063009;华北理工大学理学院,唐山063009;华北理工大学机械工程学院,唐山 063009
【正文语种】中文
【中图分类】TN4
【相关文献】
1.B(N)掺杂碳纳米管吸附Cu原子的第一性原理研究 [J], 代利峰;龚亮;刘涛;安立宝
2.B(N)掺杂碳纳米管吸附Cu原子的第一性原理研究 [J], 代利峰;龚亮;刘涛;安
立宝;
3.不同浓度硼掺杂石墨烯吸附多层金原子的第一性原理研究 [J], 卢子龙;安立宝;刘扬
4.Rh掺杂单壁碳纳米管表面NO吸附反应的第一性原理研究 [J], 熊浩;周昶利;曹政钦;魏钢;刘航;邓保家
5.N掺杂碳纳米管环吸附Fe原子的第一性原理研究 [J], 杨忠华;刘贵立;曲迎东;李荣德
因版权原因,仅展示原文概要,查看原文内容请购买。
Fe304磁性纳米粒子与CdTe量子点复合物的制备和表征实验
天津工业大学综合实验报告实验名称:磁性荧光复合材料的制备及表征2014年9月21日目录第一章绪论 (3)1.1量子点 (3)1.1.5量子点的合成 (4)第二章实验部分 (8)2.1实验药品及仪器 (8)2.2实验过程 (9)2.2.1水相MPA-CdTe量子点的合成 (9)2.2.2磁球的合成过程 (9)第三章结果与讨论 (10)3.1合成的CdTe量子点的表征 (10)3.1.1紫外表征 (10)3.1.2荧光表征 (11)3.1.3红外表征 (11)第四章结论 (13)参考文献 (14)摘要本实验通过水相回流在pH=11.9下合成了巯基丙烯酸修饰CdTe量子点并通过荧光和紫外并进行了表征;通过水解法共沉淀法合成了Fe304磁性纳米粒子以及CdTe与Fe304的复合连接。
关键字:CdTe量子点;Fe304;磁性纳米粒子;复合材料;第一章绪论量子点介绍量子点是一种内部电子运动受三维限制的零维纳米粒子材料,主要由Ⅱ-Ⅵ族或Ⅲ-Ⅴ族元素组成。
因为量子点研究范围交叉多个学科,因此不同学科研究中,其名称也不同。
如,胶体化学家称其为胶体颗粒(colloid particle);材料学家称其为纳米微晶(nanocrystal);原子分子物理学家称其为大分子、团簇(cluster);而固体与理论物理学家根据其机理形象地称其为量子点(Quantum Dots)、人造原子。
目前文献中常见的量子点主要涉及的是主族II-Ⅵ(如CdSe),III-V(如InP、InAs和GaAs)、副族化合物以及Si等元素,II-Ⅵ和III-V副族化合物尤其引起人们的关注。
量子点的特殊结构导致了它具有表面效应、量子尺寸效应、介电限域效应和宏观量子隧道效应,从而派生出与宏观体系和微观体系不同的低维物性,展现出许多不同于宏观块体材料的物理化学性质和独特的发光特性。
量子点所具有的荧光特性(1)半导体纳米晶体具有较宽激发波长的范围,较窄发射波长的范围。
掺杂Fe、Cd闪锌矿电子结构的第一性原理计算
掺杂Fe、Cd闪锌矿电⼦结构的第⼀性原理计算
掺杂Fe、Cd闪锌矿电⼦结构的第⼀性原理计算
丁聪;李艳;鲁安怀
【期刊名称】《岩⽯矿物学杂志》
【年(卷),期】2015(034)003
【摘要】采⽤基于密度泛函理论的第⼀性原理计算⽅法,对掺Fe和(或)Cd的闪锌矿型ZnS的电⼦结构进⾏了计算.计算结果表明,纯闪锌矿的禁带宽度为2.85 eV;掺Fe浓度为3.125%的闪锌矿禁带宽度为2.58 eV,且Fe的3d和S的3p 轨道杂化在禁带中引⼊了两条杂质能级;随着掺Fe量的增加,杂质能级的宽度和峰⾼也随之增⼤;掺Cd闪锌矿的禁带宽度为2.68 eV,并在下价带底引⼊⼀条杂质能级;Fe/Cd共掺杂的闪锌矿禁带宽度为2.49 eV,在禁带中出现的两条施主杂质能级可提⾼闪锌矿的可见光响应及催化能⼒.计算结果为深⼊探讨天然闪锌矿的可见光催化机制提供了理论⽀持.
【总页数】5页(382-386)
【关键词】第⼀性原理计算;闪锌矿;Fe、Cd掺杂;光催化
【作者】丁聪;李艳;鲁安怀
【作者单位】造⼭带与地壳演化教育部重点实验室,北京⼤学地球与空间科学学院,北京100871;造⼭带与地壳演化教育部重点实验室,北京⼤学地球与空间科学学院,北京100871;造⼭带与地壳演化教育部重点实验室,北京⼤学地球与空间科学学院,北京100871
【正⽂语种】中⽂
【中图分类】P578.2+3;P579
【相关⽂献】。
fe(cu)掺杂ge(si)基稀磁半导体薄膜的结构与磁性分析
AbstractToday, the information processing is developing fast. The transfer and storage ofinformation are still base on the classical methods —controlling the charge and spin of electron, respectively. With the increase of information we need, the disadvantage of the classical semiconductors and ferromagnets working independently is coming forth. The diluted magnetic semiconductors (DMSs) based on both the charge and spin of electron have been one of the most interesting fields. Due to the incorporation of transition metal elements or rare earth elements, DMSs have many strange properties, compared to the classical semiconductors. DMSs have potential applications in the future.A lot of studies on Mn doped Ⅳ group DMSs have been reported on theories and experiments. In spite of the delightful results, the discussions on the origin of the ferromagnetism and on the application are requiring further studies. In this paper, we have prepared Fe and Cu doped Si/Ge semiconductors. As following works:1. We have prepared Fe-doped Ge1-x Fe x films on Si substrates, with the substrate temperature of 473 K, and then the samples were annealed at 873 K for 20 min in vacuum. XRD measurements showed that all the samples had a Ge crystal structure, without any other secondary phases. The doped Fe atoms existed as Fe2+ mostly in the sample, with only a few Fe0. A majority of Ge atoms existed as Ge0, with some Ge-O and Ge-Fe bonds. In the sample with Fe-doped concentration about 11.3%, Fe3+ ions appeared. Magnetic measurements revealed that the samples were ferromagnetic at low temperature, and the Curie temperature was about 300 K. The ferromagnetism arises from the interaction between the random magnetic Fe atoms.2. The Si(Fe,N) films were prepared by Fe ions implantation using metal vapor vacuum arc technology and N ions implantation using Kaufman technology. XRD measurements revealed that ion-implantation did not change the crystal structure of Si matrix, without any other secondary phases. XAFS suggested that Fe ions were at thesubstituted sites in the low-dose implanted sample, and formed FeSi2 complex in theVhigh-dose implanted sample. The magnetic domain was observed in MFM. Magnetic measurements showed that the maximum saturated magnetization was 0.46 µB/Fe in the lowest-dose implanted sample (2.0×1016 cm-2). In site and post annealing both reduced the magnetization. Hall effect measurements revealed that the doped samples were n-type. Furthermore, the saturated magnetization decreased with the increase of the carries concentration. The ferromagnetism of the implanted samples was not dependent on the concentration of the carries, but on the substituted magnetic Fe atoms.3. The Cu-doped Si:Cu films were prepared by Cu-ion implanted into n-Si(100) substrate using metal vapor vacuum arc technology. XRD measurements showed that the singlecrystal structure changed to polycrystal structure after Cu-ion implantation, with some defects, such as interstitial Cu, interstitial Si and Si vacancy, but without any other secondary phases. Cu ions existed as Cu+ ions in the as-implanted samples, and after rapid thermal annealing in a N2 environment, a few Cu2+ ions appeared in the sample. Hall effect measurements revealed that the Cu-ion implanted Si samples were p-type. Magnetic measurements suggested that all the samples were ferromagnetism at room-temperature. The ferromagnetism was from the interaction between the Cu+ ions and Si vacancies mediated by the hole carriers.Key Words: magnetic semiconductor ion implantation magnetic domain Hall Effect transport propertyVI学位论文原创性声明本人所提交的学位论文《Fe/Cu掺杂Ge/Si基稀磁半导体薄膜的结构及磁性研究》是在导师的指导下,独立进行研究工作所取得的原创性成果。
Fe掺杂CdTe电磁性质的第一性原理研究
文章编号:1001G9731(2016)02G02036G05F e掺杂C d T e电磁性质的第一性原理研究∗李昭宁1,2,赵新才1,王伟平1,2(1.中国工程物理研究院流体物理研究所,四川绵阳621999;2.中国工程物理研究院高能激光科学与技术重点实验室,四川绵阳621999)摘㊀要:㊀使用基于密度泛函理论的第一性原理方法,计算研究了F e掺杂C d T e的晶格结构与电磁性质变化.研究发现,掺杂体系的晶格常数与电子结构等随掺杂F e原子的位置不同而异.通过能带结构与电子态密度的分析表明,不同占据位的F e原子表现出迥异的电子能级分布㊁轨道杂化等,从而引起体系电子性质与原子磁矩的显著变化.其中,替代位与间隙位F e杂质的原子磁矩分别为3.76和3.14玻尔磁子.这一研究对深入理解掺杂C d T e类稀磁半导体的物理性质有重要意义.关键词:㊀碲化镉;铁掺杂;电子结构;第一性原理中图分类号:㊀O469文献标识码:A D O I:10.3969/j.i s s n.1001G9731.2016.02.0080㊀引㊀言近年来,稀磁半导体(d i l u t e d m a g n e t i cs e m i c o nGd u c t o r s,简称D M S)引起了人们的广泛关注与研究.这种材料同时利用了电子的电荷属性与自旋属性,在信息的存储与处理上具有更快的速度与更低的能耗,从而在自旋电子学领域展现出广阔的应用前景[1].稀磁半导体主要通过在传统半导体中进行磁性原子掺杂来获得,如分子束外延法[2]等.一些磁性材料的居里温度甚至已经达到数百开尔文[3],具备了实用基础.同时,对稀磁半导体微观电子结构的研究工作[4],也在不断加深人们对材料性质及其物理机制的认识.稀磁半导体的研究主要集中在ⅢGⅤ族半导体以及ⅡGⅥ族半导体,而掺杂原子则有M n㊁C o等[4].其中,ⅡGⅥ族半导体的重要优势是,3d磁性原子具有更高的掺杂浓度(约35%),而在ⅢGⅤ族半导体中其掺杂浓度只有约8%.目前的困难则存在于具有磁性原子特征分布的稀磁半导体的生长制备,以及材料缺陷的处理与控制上.如何获取室温以上依然具有稳定良好性能的稀磁半导体[5G6],是目前的研究热点之一.作为重要的ⅡGⅥ族半导体,C d T e中掺M n(记为(C d,M n)T e)曾经被广泛研究[7G8].但是由于应用上存在的困难,研究焦点逐渐转移到(G a,M n)A s[9]㊁(G a,M n)N[10]㊁(Z n,M n)O[11]等材料上.现在,随着C d T e材料制备与应用工艺的提升,以及理论计算方法的成熟,有必要对基于C d T e的稀磁半导体进行更深层次的研究.因此,本文使用基于密度泛函理论(D F T)的第一性原理方法对(C d,F e)T e进行了计算研究,揭示了F e原子在不同掺杂位置(替代位和间隙位)下体系的电磁性质变化.1㊀计算方法计算模型采用由2ˑ2ˑ1单胞构成的超晶胞结构,共32个原子,即C d16T e16.其中,C d T e单胞为闪锌矿立方结构,晶格常数为0.6477n m,C dGT e初始键长为0.2805n m.建立超晶胞结构后,若使用一个F e 原子替换其中的一个C d原子,即构成替代位掺杂,记为F e1C d15T e16;若将一个F e原子置入以T e原子形成的四面体间隙位,即构成间隙位掺杂,记为F e1C d16T e16.这时体系的掺杂浓度约为3.1%.计算采用基于平面波赝势密度泛函理论的C A S T E P程序[12].交换关联势使用广义梯度近似(G G A)下的P e r d e wGB u r k eGE n z e r h o f形式[13],赝势使用V a n d e r b i l t超软赝势[14].计算时的平面波截止能量设为550e V,布里渊区的M o n k h o r s tGP a c k网格[15]设为3ˑ3ˑ6.为准确描述过渡族金属F e原子,使用G G A+U方法[16],其中对F e原子的3d电子加入U=3.0e V[17]的库伦作用.G G A+U方法考虑了同一轨道上两个自旋相反电子的强相互作用,因此可对半局域近似的G G A方法予以修正.计算时,先对超晶胞结构体系进行结构优化,使原子弛豫到平衡位置,当原子间受力小于0.01e V/n m时即认为得到收敛结构.再对收敛结构进行能量计算,得到体系的基态能量与电子结构.2㊀结果与分析掺杂体系的形成能表示如下E f=E d o p e d-E p u r e+(N p u r e,C d-N d o p e,C d)μC d-μF e ㊀㊀其中,E d o p e d为掺杂体系的自旋极化总能量,E p u r e 为未掺杂体系的总能量,N p u r e,C d和N d o p e d,C d分别为掺杂前与掺杂后体系中的C d原子数,μC d和μF e分别为C d630202016年第2期(47)卷∗基金项目:中国工程物理研究院高能激光科学与技术重点实验室基金资助项目(H E L2014G11)收到初稿日期:2015G04G15收到修改稿日期:2015G09G25通讯作者:李昭宁,EGm a i l:l i z h a o n i n g0@g m a i l.c o m 作者简介:李昭宁㊀(1986-),男,陕西岐山人,助研,硕士,从事激光与物质相互作用研究.原子与F e 原子的化学势.计算发现,F e 原子在间隙位掺杂与在替代位掺杂时体系的形成能之差为ΔE =E f (F e 1C d 16T e 16)-E f (F e 1C d 15T e 16)=1.104e V .因此,F e 杂质在能量上更倾向于占据在替代位.这与过渡族金属在ⅢGⅤ族半导体中掺杂的情况相同,但与在Ⅳ族半导体中掺杂情况相反[18].F e 杂质首先使体系的晶格结构发生变化.对于F e 1C d 15T e 16体系,F e T e 键长由优化前的0.2805n m 变为优化后的0.2644n m ,即键长缩短.同时,最近邻C d T e 键长为0.2890n m ,次近邻C d T e 键长为0.2875n m ,最远C d T e 键长为0.2858n m ,表现为距离F e 原子越远,C d T e 键长越短的规律;但其值均大于纯C d T e 体系中的C d T e 键长(0.2805n m ).由表1可以看出,在结构优化过程中,最近邻T e 原子向F e原子靠近,其余原子均远离.这主要与F e 原子具有较小的原子半径以及原子间相互作用有关.对F e 1C d 16T e 16体系,F e T e 键长由优化前的0.2805n m 变为优化后的0.2847n m ,即键长变长.同时,最近邻C d T e 键长为0.2891n m ,次近邻C d T e 键长为0.2867n m ,最远C d T e 键长为0.2895n m ,表现为随着距离F e 原子变远,C d T e 键长交替增减的规律.同时,由表1可知,在结构优化过程中,所有原子均向远离F e 原子方向移动.在两种掺杂情况下,体系的晶格常数均变大,分别由掺杂前的0.6477n m 增至0.6589n m (替代掺杂)和0.6685n m (间隙掺杂).表1㊀掺杂体系中杂质F e 原子与其它原子的间距T a b l e 1D i s t a n c e s b e t w e e n t h eF e i m p u r i t y an d t h eh o s t a t o m s 掺杂体系F e 原子与第n 近邻位原子间距/n m 1234567F e 1C d 15T e 16优化前0.28050.45800.53700.64770.70580.79330.9160优化后0.26440.46190.54630.65940.71810.80740.9326F e 1C d 16T e 16优化前0.28050.32390.53700.56090.70580.72420.9716优化后0.28470.33820.55330.58620.72710.74421.0001㊀㊀其次,掺杂体系的电子结构发生了很大变化.图1给出了在不同自旋通道上的能带结构图.图1㊀F e 掺杂C d T e 体系的能带结构图F i g 1E l e c t r o n i cb a n ds t r u c t u r e s f o rF ed o p e dC d T e s ys t e m s 图1中(a )㊁(b )为替代位F e 掺杂体系,图1(c)㊁(d )为间隙位F e 掺杂体系.图1(a )㊁(c)为自旋向上电子能带结构,图1(b )㊁(d )为自旋向下电子能带结构.图中虚线为费米能级,红色曲线表示来自F e 原子的3d 电子能态,绿色为价带顶,蓝色为导带底.当F e 原子在替代位掺杂时(图1(a )㊁(b )),体系的半导体属性没有发生改变,依然为Γ点处具有最小带隙的直接带隙半导体.需要注意的是,G G A 计算会低估半导体的带隙宽度E g ,此处E g =0.59e V ,远小于未掺杂C d T e带隙的实验值(1.490e V ).但计算结果的定性描述是准确并有价值的.由图1(a ),在自旋向上通道,费米能级附近的能带结构与掺杂前相比没有太大变化,其中,价带顶能态主要由T e 原子的5p 电子态所构成.图2为电子轨道分布图.图2(a )和(b )给出体系最高分子占据轨道(HOMO )在掺杂前后的空间分布.图2(a )为未掺杂体系的HOMO 轨道,图2(b )为F e 1C d 15T e 16体系的自旋向上HOMO 轨道;图2(c )㊁(d )分别为掺杂体系F e 1C d 15T e 16中杂质F e 原子的自旋向上与自旋向下3d 电子轨道;图2(e )㊁(f )分别为掺杂体系F e 1C d 16T e 16中杂质F e 原子的自旋向上与自旋向下3d 电子轨道.很明显,掺杂前HOMO 均匀分布于T e 原子上,掺杂后则只分布于杂质F e 原子的最近邻T e 原子上,表明掺杂引起HOMO 的定域化.进一步观察能带图发现,在-4.705与-5.134e V 处出现由杂质F e 原子所引入的平直能带.图2(c )给出其中一条平直能态对应的电子轨道的空间分布,可以看到该电子轨道仅分布于F e 原子周围空间.图3(a )给出相应掺杂体系的电子态密度(D O S )图.显然,这些平直能带来自F e 原子的3d 轨道电子态.由于能量分布上具有很强的定域性,这些能级上的电子轨道没有与近邻T e 73020李昭宁等:F e 掺杂C d T e 电磁性质的第一性原理研究原子间发生轨道杂化作用.在自旋向下通道(图1(b )),首先看到约5.0e V 处平直能带的消失,而在导带底部则出现了F e 原子的杂质能带.该杂质能带对应的电子轨道分布如图2(d )所示.从D O S 图可以得到,这些新能带来自F e 原子3d 壳层未被填充的电子轨道态.其中,导带底能带在布里渊区由Γ点处的0.646e V 扩展至X 点处的1.128e V ,使得自旋向下通道带隙为0.82e V ,大于自旋向上通道的0.59e V .分析可知,F e 原子的6个3d 电子中有4个占据平直能带处的四重3d 轨道态(自旋向上),其余2个占据能量稍高的F e 3d GT e5p 杂化轨道态(一个自旋向上,一个自旋向下).因此,F e 原子磁矩为M F e =[n ʏ(x )-n ˌ(x )]μB =(5-1)μB =4.0μB .这与计算结果给出的3.76μB 相符,也与F e 原子在S i 中进行替代位掺杂时的磁矩相同[17].这时,体系的总磁矩为4.00μB .此外,由D O S 图可以得到,T e 原子的5s 电子态主要局域在-10.670e V 处;C d 原子的4d 电子态则局域在-9.044e V 处.图2㊀电子轨道分布图F i g 2Di s t r i b u t i o n s o fm o l e c u l e o r b i t a l s ㊀㊀当F e 原子在间隙位掺杂时,体系电子结构发生了很大变化.由能带图可以看到(图1(c )㊁(d )),费米能级向高能方向移动了1.416e V ,进入导带,使体系由半导体转变为导体.在自旋向上通道,没有观察到替代位F e 掺杂时约-5.0e V 处的平直能带.分析可知,这些能态向高能方向移动,主要分布于约-4.0e V 附近,并扩展分布于-6.078~-1.556e V 能量区间.图2(e)给出其中一条能态对应的电子轨道分布,可以看到,电子轨道不再像在替代位掺杂时那样局域在F e 原子周围,而是进一步扩展至近邻T e 原子空间.图3为F e 掺杂C d T e 体系的电子态密度图.其中,图3(a)为替代位F e 掺杂体系,图3(b )为间隙位F e 掺杂体系.图中还给出了F e 原子3d 轨道㊁近邻C d 原子4d 轨道㊁T e 原子5s 轨道㊁5p 轨道电子的分波态密度.由DO S 图分析可知(图3(b )),这些来自F e 原子的3d 电子态与T e 原子的5p 电子态发生轨道杂化,同时还伴随有一定的F e 4s GT e 5p 轨道杂化.这时,F e 原子同样有5个电子占据自旋向上3d 杂化轨道态.不同的是,其中4个电子来自初始的3d 轨道,1个电子来自初始的4s 轨道.剩下的1个4s 电子则以概率形式平均占据自旋向上与自旋向下的F e 4s GT e 5p 杂化态.在自旋向下通道,费米能级下0.997e V 处出现从Γ点到X 点平直延伸的能带.由D O S 图可以得出,该平直能带来自F e 原子的3d 电子态,且局域在禁带中心位置,没有与近邻原子间发生轨道杂化.图2(f )给出与该能态对应的电子轨道分布.分析可知,该能态为二重简并态,可容纳2个3d 电子占据.因此,间隙位F e 原子磁矩为M F e =[n ʏ(x )-n ˌ(x )]μB =(5-2)μB =3.0μB .这与计算结果给出的3.14μB 相符.但与S i 中间隙位F e 掺杂时2.0μB 的结论相异[17].这时,体系的总磁矩为3.16μB .此外,未占据的3d 电子态不再分布在导带底,而是位于导带中较高能量处.T e 原子的5s 态㊁C d原子的4d 态相对于价带顶的位置没有发生明显变化.由上可见,F e 原子掺杂位置不同时,体系表现出830202016年第2期(47)卷不同的电子结构,特别是F e原子的3d电子态分布,使杂质F e原子具有不同的原子磁矩.这与M n原子掺杂C d T e时体系电子结构上的稳定性是截然不同的.图4给出M n掺杂C d T e体系的D O S图.图3㊀F e掺杂C d T e体系的电子态密度图.F i g3E l e c t r o nd e n s i t y o fs t a t e sf o rF ed o p e d C d T es y s t e m s图4㊀M n掺杂C d T e体系的电子态密度图F i g4E l e c t r o nd e n s i t y o f s t a t e s f o r M nd o p e dC d T es y s t e m s其中,图4(a)为替代位M n掺杂体系,图4(b)为间隙位M n掺杂体系.图中还给出了M n原子3d轨道㊁近邻C d原子3d轨道㊁T e原子5s轨道㊁5p轨道电子的分波态密度.可以看到,替代位掺杂时,M n原子的3d电子态局域在价带深处约-4.743e V处,并且不与近邻原子轨道间发生轨道杂化.特别是,这些3d电子占据态几乎全部处于自旋向上通道,自旋向下通道的3d电子态则分布在导带中,形成未占据态.因此, M n原子的磁矩为M F e=[nʏ(x)-nˌ(x)]μB=(5-0)μB=5.0μB.这与计算结果给出的M n原子磁矩5.08μB相符.而间隙位掺杂时,除了F e r m i能级向高能方向移动进入导带外,体系电子结构并没有发生显著变化,M n原子磁矩为5.18μB.3㊀结㊀论使用基于密度泛函理论的第一性原理方法,计算研究了F e掺杂C d T e晶体的晶格结构㊁电子结构等性质变化.结果发现:(1)㊀F e掺杂体系中原子的键长㊁弛豫位置随杂质原子位置(替代位或间隙位)不同而异,但体系的晶格常数都变大.在形成能上,F e原子更倾向于占据替代位.(2)㊀替代位F e掺杂时,体系保持半导体性.F e 原子的自旋向上3d电子态主要以定域态形式分布在价带深处,自旋向下3d电子态则主要分布于导带底形成未占据态.F e原子磁矩为3.76μB.(3)㊀间隙位F e掺杂时,费米能级进入导带中,体系转变为导体.F e原子的自旋向上3d电子态以F e 3dGT e5p轨道杂化态形式存在,自旋向下3d电子态则局域在禁带中.F e原子磁矩为3.14μB.F e掺杂C d T e体系的电子结构表现出对掺杂位置的依赖性,这与M n掺杂体系明显不同.在下一步的研究中,需要对掺杂体系中F eGF e原子间的耦合作用进行研究[19],从而对(C d,F e)T e作为稀磁半导体的实用性建立更充分的理论基础.参考文献:[1]㊀O h n o H.M a k i n g n o n m a g n e t i c s e m i c o n d u c t o r s f e r r o m a gGn e t i c[J].S c i e n c e,1998,281(5379):951G956.[2]㊀O h n o H,S h e nA,M a t s u k u r aF,e t a l.(G a,M n)A s:a n e wd i l u t e d m a g n e t i cs e m i c o n d u c t o rb a s e do nG a A s[J].A p p l P h y sL e t t,1996,69(3):363G365.[3]㊀J a m e tM,B a r s k iA,D e V i l l e r sT,e t a l.H i g hGC u r i eGt e mGp e r a t u r e f e r r o m a g n e t i s m i n s e l fGo r g a n i z e d G e1Gx M n xn a n o c o l u m n s[J].N a t u r e M a t e r i a l s,2006,5(8):653G659.[4]㊀S a t o K,B e r g q v i s tL,K u d r n o v s k y J,e ta l.F i r s tGp r i n c iGp l e st h e o r y o fd i l u t e m a g n e t i cs e m i c o n d u c t o r s[J].R e vM o dP h y s,2010,82(2):1633G1690.[5]㊀B o n a n n iA.F e r r o m a g n e t i cn i t r i d eGb a s e ds e m i c o n d u c t o r sd o pe dw i t ht r a n s i t i o na n dr a r ee a r t h s[J].S e m i c o n dS c iT e c h n o l,2007,22(9):R41GR56.[6]㊀P e a r t o nSJ,H e o W H,I v i l lM,e ta l.D i l u t e m a g n e t i c93020李昭宁等:F e掺杂C d T e电磁性质的第一性原理研究s e m i c o n d u t i n g o x i d e s[J].S e m i c o n dS c iT e c h n o l,2004,19(10):R59GR74.[7]㊀F u r d y n aJK,S a m a r t hN.M a g n e t i c p r o p e r t i e so f d i l u t e d m a g n e t i cs e m i c o n d u c t o r s:ar e v i e w[J].J A p p lP h y s,1987,61(8):3526G3531.[8]㊀F u r d y n a JK.D i l u t e dm a g n e t i c s e m i c o n d u c t o r s[J].JA pGp l P h y s,1988,64(4):R29GR64.[9]㊀M i r b tS,S a n y a lB,M o h nP.M a g n e t i c p r o p e r t i e so f3di m p u r i t i e ss u b s t i t u t e di n G a A s[J].J P h y s:C o n d e n sM a t t e r,2002,14(12):3295G3302.[10]㊀L u o X,M a r t i n R M.J a h nGT e l l e rd i s t o r t i o na n df e r r oGm a g n e t i s m i n t h e d i l u t e m a g n e t i c s e m i c o n d u c t o r sG a A sʒM na n dc u b i c G a NʒM n[J].P h y s R e v B,2005,72(3):035212.[11]㊀S a t oK,K a t a y a m aGY o s h i d aH.F i r s t p r i n c i p l e sm a t e r i a l sd e s i g n f o r s e m i c o n d u c t o r s p i n t r o n i c s[J].S e m i c o n dS c iT e c h n o l,2002,17(4):367G376.[12]㊀C l a r kSJ,S e g a l lM D,P i c k a r dCJ,e t a l.F i r s t p r i n c iGp l e sm e t h o d su s i n g C A S T E P[J].Z e i t s c h r i f tFürK r iGs t a l l o g r a p h i e,2005,220(5G6):567G570.[13]㊀P e r d e wJP,B u r k eK,E r n z e r h o fM.G e n e r a l i z e d g r a d iGe n ta p p r o x i m a t i o n m a d es i m p l e[J].P h y s R e v L e t t,1996,77(18):3865G3868.[14]㊀V a n d e r b i l tD.S o f ts e l fGc o n s i s t e n t p s e u d o p o t e n t i a l s i nag e n e r a l i z e d e i g e n v a l u ef o r m a l i s m[J].P h y s R e v B,1990,41(11):R7892GR7895.[15]㊀M o n k h o r s tH J,P a c kJD.S p e c i a l p o i n t s f o rb r i l l o u i nGz o n e i n t e g r a t i o n s[J].P h y sR e vB,1976,13(12):5188G5192.[16]㊀D u d a r e vSL,B o t t o n G A,S a v r a s o vS Y,e ta l.E l e cGt r o nGe n e r g yGl o s ss p e c t r aa n dt h es t r u c t u r a ls t a b i l i t y o fn i c k l eo x i d e:a n L S D A+U s t u d y[J].P h y s R e v B,1998,57(3):1505G1509.[17]㊀Küw e nF,L e i t s m a n nR,B e c h s t e d t F.M n a n dF e d o p i n g o fb u l k S i:c o n c e n t r a t i o ni n f l u e n c e o n e l e c t r o n i ca n dm a g n e t i c p r o p e r t i e s[J].P h y s R e v B,2009,80(4):045203.[18]㊀A l Z a h r a n iAZ,S r i v a s t a v aGP,G a r g R,e t a l.A n a b i nGi t i o s t u d y o f e l e c t r o n i c a n d s t r u c t u r a l p r o p e r t i e s o fM n i na G a A se n v i r o n m e n t[J].J P h y s:C o n d e n s M a t t e r,2009,21(48):485504.[19]㊀L iZ h a o n i n g,W a n g W e i p i n g,Z h a oS h i c a o.M nGM ni nGt e r a c t i o ni n d u c e d m e t a l l i c o ri n s u l a t i n g c h a r a c t e r o fd o pe ds i l i c o n:a n a b i n i t i o s t u d y[J].C o m p u tM a t e r S c i,2015,97(1):186G192.F i r s tGp r i n c i p l e s s t u d y o f e l e c t r o n i c a n dm a g n e t i c p r o p e r t i e s o fF e d o p i n g i nC d T e c r y s t a lL I Z h a o n i n g1,2,Z H A O X i n c a i1,WA N G W e i p i n g1,2(1.I n s t i t u t e o f F l u i dP h y s i c s,C h i n aA c a d e m y o fE n g i n e e r i n g P h y s i c s,M i a n y a n g621999,C h i n a;2.T h eK e y L a b o r a t o r y o f S c i e n c e a n dT e c h n o l o g y o nH i g hE n e r g y L a s e r,C h i n aA c a d e m y o fE n g i n e e r i n g P h y s i c s,M i a n y a n g621999,C h i n a)A b s t r a c t:T h e s t r u c t u r a l,e l e c t r o n i ca n d m a g n e t i c p r o p e r t i e so fF e i m p u r i t y d o p i n g i nC d T ea r es t u d i e du s i n g f i r s t p r i n c i p l e s b a s e do n t h e d e n s i t y f u n c t i o n a l t h e o r y(D F T)i n a s u p e r c e l l a p p r o a c h.I no r d e r t om a k e a n a c c uGr a t ed e s c r i p t i o n o f t r a n s i t i o n a t o mF e,a g e n e r a l i z e d g r a d i e n t a p p r o x i m a t i o n(G G A)+U m e t h o d i s u s e d.T h e r eGs u l t s s h o wt h a t t h e r e a r e s i g n i f i c a n t c h a n g e s i n t h e b o n d l e n g t h s,m a g n e t i z a t i o n a n d e l e c t r o n i c p r o p e r t i e s.B a s e d o ne l e c t r o n i c b a n d s t r u c t u r e a n d e l e c t r o nd e n s i t y o f s t a t e s,v a r i e s i n e l e c t r o n i c s t r u c t u r e s i n d u c e db y o r b i t a l h yGb r i d i z a t i o nw i t h r e s p e c t t o d i f f e r e n t d o p i n g s i t e t y p e s(s u b s t i t u t i o n a l o r i n t e r s t i t i a l)a r e a n a l y s e d.I t i s f o u n d t h a t t h e l o c a lm a g n e t i cm o m e n t s o f s u b s t i t u t i o n a l a n d i n t e r s t i t i a lF e i m p u r i t i e s a r e3.76a n d3.14B o h rm a g n e t o n s, r e s p e c t i v e l y.T h e s e c o n c l u s i o n s p r o v i d es o m eu s e f u l r e f e r e n c e t ou n d e r s t a n dt h e p h y s i c a lm e c h a n i s m o fC d T e b a s e dd i l u t e dm a g n e t i c s e m i c o n d u c t o r s.K e y w o r d s:C d T e;F e i m p u r i t i e s;e l e c t r o n i c s t r u c t u r e;f i r s tGp r i n c i p l e s s t u d y040202016年第2期(47)卷。
Fe掺杂石墨烯表面吸附Au:第一性原理
Fe掺杂石墨烯表面吸附Au:第一性原理刘扬;安立宝;龚亮【摘要】The first principles method was used to study the adsorption properties of Au atoms on Fe-doped graphene.The adsorption energy of Au atoms on graphene before and after Fe doping,as well as the partial density of state and charge distribution of graphene werecalcutated.Results showed that doping Fe atoms could increase adsorption energy of Au atoms on graphene and enhance charge transfer between Au atoms and graphene.The improvement was mainly due to the formation of covalent bonds between Fe atoms and Au atoms after Fe doping,which converted the physical adsorption between Au and Fe atoms to a chemical one.Besides,the spin-polarized calculation results showed that spin-up and spin-down density of states were symmetrical about Fermi level for the system of Fe-doped graphene and the adsorbed Au atom,which meant that no magnetic moment existed in the system and consequently the properties of the material would not change.Therefore,doping Fe could improve the catalytic performance of graphene-supported Au nanoparticle catalysts.%采用第一性原理对Au在Fe掺杂石墨烯表面的吸附特性进行研究,计算掺杂Fe前后石墨烯对Au的吸附能以及石墨烯的局部态密度和电荷分布.结果表明:掺杂Fe增强了Au在石墨烯表面的吸附能,提高了Au与石墨烯间的电荷转移.掺杂Fe后,Fe与Au间形成了共价键,Au与石墨烯间由物理吸附转变为化学吸附.此外,自旋极化计算结果显示Fe掺杂石墨烯-Au体系中自旋向上态密度和自旋向下态密度关于费米能级对称,说明Fe掺杂不会使体系产生磁矩而影响材料性能.此类掺杂有望改善石墨烯载Au纳米颗粒催化剂的催化性能.【期刊名称】《南京工业大学学报(自然科学版)》【年(卷),期】2018(040)003【总页数】6页(P109-114)【关键词】第一性原理;石墨烯;Fe掺杂;Au吸附【作者】刘扬;安立宝;龚亮【作者单位】华北理工大学机械工程学院,河北唐山063210;华北理工大学机械工程学院,河北唐山063210;华北理工大学机械工程学院,河北唐山063210【正文语种】中文【中图分类】TB321Au纳米颗粒因较大的比表面积、优异的导热性、良好的生物相容性和独特的光电性能,使其催化机制及在催化方面的应用成为当前的研究热点[1-3]。
第一性原理的计算方法及常用软件介绍
第二章研究方法与程序介绍§2.1 全电子法和赝势法应用于铁电体的第一性原理计算方法和工具很多,根据对势函数及内层电子的处理方法不同主要分为两大类,一种是波函数中包含了高能态和内层电子,而势函数只是原子核的贡献,这称为全电子(all electron calculation)法,另一种处理方法是势函数为原子核和内层电子联合产生的势,称为离子赝势,波函数只是高能态电子的函数,这称为赝势(pseudo-potential)法。
因为内层电子对价电子的排斥作用部分地抵消了原子核对价电子的强吸引作用,所以赝势是一种比较弱和比较平坦的势。
引入赝势的要点在于,赝势对应的薛定谔方程与真实势对应的薛定谔方程有相同的能量本征值。
在这一前提下,引入赝势的方法不是唯一的。
在第一性原理计算中,用的是所谓模守恒赝势法。
这种赝势所对应的波函数有一个特点,在离开原子核一定距离的空间,它与真实势对应的波函数不但形式相同,而且幅度相等,故称模守恒。
这种方法从原子势算起,不引入任何实验参数,所以又称为从头算起(ab initio)赝势方法。
一般来说,赝势法计算量较小,但其中消去了内层电子态,相对于全电子法多引入了一个近似。
该方法的优点是较便于计算离子受到的作用力,后者等于总能量对原胞内离子位矢导数的负数,称为Hellmann-Feynman力。
赝势法用平面波展开来表示价电子态,如果晶体中的原子有2p未满壳层(如氧)或3d未满壳层(如钛),则赝势将很“硬”,为满足模守恒,需要为数很多的平面波基函数,计算量太大。
为此发展了超软赝势(ultro-soft pseudo-potential)法。
对波函数引入一个重叠算符,使赝势变软,减少了平面波基函数。
在铁电体研究中用的赝势法通常是这种方法。
全电子法表示电子态时将空间分为两部分:一是原子核附近的球形区,称为丸盒(muffin-tin)区,二是原子核间的其它区域。
在球形区,基函数、电荷密度16和势均用径向函数展开,在其它区域,这些量用平面波或球面波展开。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
文章编号:1001G9731(2016)02G02036G05F e掺杂C d T e电磁性质的第一性原理研究∗李昭宁1,2,赵新才1,王伟平1,2(1.中国工程物理研究院流体物理研究所,四川绵阳621999;2.中国工程物理研究院高能激光科学与技术重点实验室,四川绵阳621999)摘㊀要:㊀使用基于密度泛函理论的第一性原理方法,计算研究了F e掺杂C d T e的晶格结构与电磁性质变化.研究发现,掺杂体系的晶格常数与电子结构等随掺杂F e原子的位置不同而异.通过能带结构与电子态密度的分析表明,不同占据位的F e原子表现出迥异的电子能级分布㊁轨道杂化等,从而引起体系电子性质与原子磁矩的显著变化.其中,替代位与间隙位F e杂质的原子磁矩分别为3.76和3.14玻尔磁子.这一研究对深入理解掺杂C d T e类稀磁半导体的物理性质有重要意义.关键词:㊀碲化镉;铁掺杂;电子结构;第一性原理中图分类号:㊀O469文献标识码:A D O I:10.3969/j.i s s n.1001G9731.2016.02.0080㊀引㊀言近年来,稀磁半导体(d i l u t e d m a g n e t i cs e m i c o nGd u c t o r s,简称D M S)引起了人们的广泛关注与研究.这种材料同时利用了电子的电荷属性与自旋属性,在信息的存储与处理上具有更快的速度与更低的能耗,从而在自旋电子学领域展现出广阔的应用前景[1].稀磁半导体主要通过在传统半导体中进行磁性原子掺杂来获得,如分子束外延法[2]等.一些磁性材料的居里温度甚至已经达到数百开尔文[3],具备了实用基础.同时,对稀磁半导体微观电子结构的研究工作[4],也在不断加深人们对材料性质及其物理机制的认识.稀磁半导体的研究主要集中在ⅢGⅤ族半导体以及ⅡGⅥ族半导体,而掺杂原子则有M n㊁C o等[4].其中,ⅡGⅥ族半导体的重要优势是,3d磁性原子具有更高的掺杂浓度(约35%),而在ⅢGⅤ族半导体中其掺杂浓度只有约8%.目前的困难则存在于具有磁性原子特征分布的稀磁半导体的生长制备,以及材料缺陷的处理与控制上.如何获取室温以上依然具有稳定良好性能的稀磁半导体[5G6],是目前的研究热点之一.作为重要的ⅡGⅥ族半导体,C d T e中掺M n(记为(C d,M n)T e)曾经被广泛研究[7G8].但是由于应用上存在的困难,研究焦点逐渐转移到(G a,M n)A s[9]㊁(G a,M n)N[10]㊁(Z n,M n)O[11]等材料上.现在,随着C d T e材料制备与应用工艺的提升,以及理论计算方法的成熟,有必要对基于C d T e的稀磁半导体进行更深层次的研究.因此,本文使用基于密度泛函理论(D F T)的第一性原理方法对(C d,F e)T e进行了计算研究,揭示了F e原子在不同掺杂位置(替代位和间隙位)下体系的电磁性质变化.1㊀计算方法计算模型采用由2ˑ2ˑ1单胞构成的超晶胞结构,共32个原子,即C d16T e16.其中,C d T e单胞为闪锌矿立方结构,晶格常数为0.6477n m,C dGT e初始键长为0.2805n m.建立超晶胞结构后,若使用一个F e 原子替换其中的一个C d原子,即构成替代位掺杂,记为F e1C d15T e16;若将一个F e原子置入以T e原子形成的四面体间隙位,即构成间隙位掺杂,记为F e1C d16T e16.这时体系的掺杂浓度约为3.1%.计算采用基于平面波赝势密度泛函理论的C A S T E P程序[12].交换关联势使用广义梯度近似(G G A)下的P e r d e wGB u r k eGE n z e r h o f形式[13],赝势使用V a n d e r b i l t超软赝势[14].计算时的平面波截止能量设为550e V,布里渊区的M o n k h o r s tGP a c k网格[15]设为3ˑ3ˑ6.为准确描述过渡族金属F e原子,使用G G A+U方法[16],其中对F e原子的3d电子加入U=3.0e V[17]的库伦作用.G G A+U方法考虑了同一轨道上两个自旋相反电子的强相互作用,因此可对半局域近似的G G A方法予以修正.计算时,先对超晶胞结构体系进行结构优化,使原子弛豫到平衡位置,当原子间受力小于0.01e V/n m时即认为得到收敛结构.再对收敛结构进行能量计算,得到体系的基态能量与电子结构.2㊀结果与分析掺杂体系的形成能表示如下E f=E d o p e d-E p u r e+(N p u r e,C d-N d o p e,C d)μC d-μF e ㊀㊀其中,E d o p e d为掺杂体系的自旋极化总能量,E p u r e 为未掺杂体系的总能量,N p u r e,C d和N d o p e d,C d分别为掺杂前与掺杂后体系中的C d原子数,μC d和μF e分别为C d630202016年第2期(47)卷∗基金项目:中国工程物理研究院高能激光科学与技术重点实验室基金资助项目(H E L2014G11)收到初稿日期:2015G04G15收到修改稿日期:2015G09G25通讯作者:李昭宁,EGm a i l:l i z h a o n i n g0@g m a i l.c o m 作者简介:李昭宁㊀(1986-),男,陕西岐山人,助研,硕士,从事激光与物质相互作用研究.原子与F e 原子的化学势.计算发现,F e 原子在间隙位掺杂与在替代位掺杂时体系的形成能之差为ΔE =E f (F e 1C d 16T e 16)-E f (F e 1C d 15T e 16)=1.104e V .因此,F e 杂质在能量上更倾向于占据在替代位.这与过渡族金属在ⅢGⅤ族半导体中掺杂的情况相同,但与在Ⅳ族半导体中掺杂情况相反[18].F e 杂质首先使体系的晶格结构发生变化.对于F e 1C d 15T e 16体系,F e T e 键长由优化前的0.2805n m 变为优化后的0.2644n m ,即键长缩短.同时,最近邻C d T e 键长为0.2890n m ,次近邻C d T e 键长为0.2875n m ,最远C d T e 键长为0.2858n m ,表现为距离F e 原子越远,C d T e 键长越短的规律;但其值均大于纯C d T e 体系中的C d T e 键长(0.2805n m ).由表1可以看出,在结构优化过程中,最近邻T e 原子向F e原子靠近,其余原子均远离.这主要与F e 原子具有较小的原子半径以及原子间相互作用有关.对F e 1C d 16T e 16体系,F e T e 键长由优化前的0.2805n m 变为优化后的0.2847n m ,即键长变长.同时,最近邻C d T e 键长为0.2891n m ,次近邻C d T e 键长为0.2867n m ,最远C d T e 键长为0.2895n m ,表现为随着距离F e 原子变远,C d T e 键长交替增减的规律.同时,由表1可知,在结构优化过程中,所有原子均向远离F e 原子方向移动.在两种掺杂情况下,体系的晶格常数均变大,分别由掺杂前的0.6477n m 增至0.6589n m (替代掺杂)和0.6685n m (间隙掺杂).表1㊀掺杂体系中杂质F e 原子与其它原子的间距T a b l e 1D i s t a n c e s b e t w e e n t h eF e i m p u r i t y an d t h eh o s t a t o m s 掺杂体系F e 原子与第n 近邻位原子间距/n m 1234567F e 1C d 15T e 16优化前0.28050.45800.53700.64770.70580.79330.9160优化后0.26440.46190.54630.65940.71810.80740.9326F e 1C d 16T e 16优化前0.28050.32390.53700.56090.70580.72420.9716优化后0.28470.33820.55330.58620.72710.74421.0001㊀㊀其次,掺杂体系的电子结构发生了很大变化.图1给出了在不同自旋通道上的能带结构图.图1㊀F e 掺杂C d T e 体系的能带结构图F i g 1E l e c t r o n i cb a n ds t r u c t u r e s f o rF ed o p e dC d T e s ys t e m s 图1中(a )㊁(b )为替代位F e 掺杂体系,图1(c)㊁(d )为间隙位F e 掺杂体系.图1(a )㊁(c)为自旋向上电子能带结构,图1(b )㊁(d )为自旋向下电子能带结构.图中虚线为费米能级,红色曲线表示来自F e 原子的3d 电子能态,绿色为价带顶,蓝色为导带底.当F e 原子在替代位掺杂时(图1(a )㊁(b )),体系的半导体属性没有发生改变,依然为Γ点处具有最小带隙的直接带隙半导体.需要注意的是,G G A 计算会低估半导体的带隙宽度E g ,此处E g =0.59e V ,远小于未掺杂C d T e带隙的实验值(1.490e V ).但计算结果的定性描述是准确并有价值的.由图1(a ),在自旋向上通道,费米能级附近的能带结构与掺杂前相比没有太大变化,其中,价带顶能态主要由T e 原子的5p 电子态所构成.图2为电子轨道分布图.图2(a )和(b )给出体系最高分子占据轨道(HOMO )在掺杂前后的空间分布.图2(a )为未掺杂体系的HOMO 轨道,图2(b )为F e 1C d 15T e 16体系的自旋向上HOMO 轨道;图2(c )㊁(d )分别为掺杂体系F e 1C d 15T e 16中杂质F e 原子的自旋向上与自旋向下3d 电子轨道;图2(e )㊁(f )分别为掺杂体系F e 1C d 16T e 16中杂质F e 原子的自旋向上与自旋向下3d 电子轨道.很明显,掺杂前HOMO 均匀分布于T e 原子上,掺杂后则只分布于杂质F e 原子的最近邻T e 原子上,表明掺杂引起HOMO 的定域化.进一步观察能带图发现,在-4.705与-5.134e V 处出现由杂质F e 原子所引入的平直能带.图2(c )给出其中一条平直能态对应的电子轨道的空间分布,可以看到该电子轨道仅分布于F e 原子周围空间.图3(a )给出相应掺杂体系的电子态密度(D O S )图.显然,这些平直能带来自F e 原子的3d 轨道电子态.由于能量分布上具有很强的定域性,这些能级上的电子轨道没有与近邻T e 73020李昭宁等:F e 掺杂C d T e 电磁性质的第一性原理研究原子间发生轨道杂化作用.在自旋向下通道(图1(b )),首先看到约5.0e V 处平直能带的消失,而在导带底部则出现了F e 原子的杂质能带.该杂质能带对应的电子轨道分布如图2(d )所示.从D O S 图可以得到,这些新能带来自F e 原子3d 壳层未被填充的电子轨道态.其中,导带底能带在布里渊区由Γ点处的0.646e V 扩展至X 点处的1.128e V ,使得自旋向下通道带隙为0.82e V ,大于自旋向上通道的0.59e V .分析可知,F e 原子的6个3d 电子中有4个占据平直能带处的四重3d 轨道态(自旋向上),其余2个占据能量稍高的F e 3d GT e5p 杂化轨道态(一个自旋向上,一个自旋向下).因此,F e 原子磁矩为M F e =[n ʏ(x )-n ˌ(x )]μB =(5-1)μB =4.0μB .这与计算结果给出的3.76μB 相符,也与F e 原子在S i 中进行替代位掺杂时的磁矩相同[17].这时,体系的总磁矩为4.00μB .此外,由D O S 图可以得到,T e 原子的5s 电子态主要局域在-10.670e V 处;C d 原子的4d 电子态则局域在-9.044e V 处.图2㊀电子轨道分布图F i g 2Di s t r i b u t i o n s o fm o l e c u l e o r b i t a l s ㊀㊀当F e 原子在间隙位掺杂时,体系电子结构发生了很大变化.由能带图可以看到(图1(c )㊁(d )),费米能级向高能方向移动了1.416e V ,进入导带,使体系由半导体转变为导体.在自旋向上通道,没有观察到替代位F e 掺杂时约-5.0e V 处的平直能带.分析可知,这些能态向高能方向移动,主要分布于约-4.0e V 附近,并扩展分布于-6.078~-1.556e V 能量区间.图2(e)给出其中一条能态对应的电子轨道分布,可以看到,电子轨道不再像在替代位掺杂时那样局域在F e 原子周围,而是进一步扩展至近邻T e 原子空间.图3为F e 掺杂C d T e 体系的电子态密度图.其中,图3(a)为替代位F e 掺杂体系,图3(b )为间隙位F e 掺杂体系.图中还给出了F e 原子3d 轨道㊁近邻C d 原子4d 轨道㊁T e 原子5s 轨道㊁5p 轨道电子的分波态密度.由DO S 图分析可知(图3(b )),这些来自F e 原子的3d 电子态与T e 原子的5p 电子态发生轨道杂化,同时还伴随有一定的F e 4s GT e 5p 轨道杂化.这时,F e 原子同样有5个电子占据自旋向上3d 杂化轨道态.不同的是,其中4个电子来自初始的3d 轨道,1个电子来自初始的4s 轨道.剩下的1个4s 电子则以概率形式平均占据自旋向上与自旋向下的F e 4s GT e 5p 杂化态.在自旋向下通道,费米能级下0.997e V 处出现从Γ点到X 点平直延伸的能带.由D O S 图可以得出,该平直能带来自F e 原子的3d 电子态,且局域在禁带中心位置,没有与近邻原子间发生轨道杂化.图2(f )给出与该能态对应的电子轨道分布.分析可知,该能态为二重简并态,可容纳2个3d 电子占据.因此,间隙位F e 原子磁矩为M F e =[n ʏ(x )-n ˌ(x )]μB =(5-2)μB =3.0μB .这与计算结果给出的3.14μB 相符.但与S i 中间隙位F e 掺杂时2.0μB 的结论相异[17].这时,体系的总磁矩为3.16μB .此外,未占据的3d 电子态不再分布在导带底,而是位于导带中较高能量处.T e 原子的5s 态㊁C d原子的4d 态相对于价带顶的位置没有发生明显变化.由上可见,F e 原子掺杂位置不同时,体系表现出830202016年第2期(47)卷不同的电子结构,特别是F e原子的3d电子态分布,使杂质F e原子具有不同的原子磁矩.这与M n原子掺杂C d T e时体系电子结构上的稳定性是截然不同的.图4给出M n掺杂C d T e体系的D O S图.图3㊀F e掺杂C d T e体系的电子态密度图.F i g3E l e c t r o nd e n s i t y o fs t a t e sf o rF ed o p e d C d T es y s t e m s图4㊀M n掺杂C d T e体系的电子态密度图F i g4E l e c t r o nd e n s i t y o f s t a t e s f o r M nd o p e dC d T es y s t e m s其中,图4(a)为替代位M n掺杂体系,图4(b)为间隙位M n掺杂体系.图中还给出了M n原子3d轨道㊁近邻C d原子3d轨道㊁T e原子5s轨道㊁5p轨道电子的分波态密度.可以看到,替代位掺杂时,M n原子的3d电子态局域在价带深处约-4.743e V处,并且不与近邻原子轨道间发生轨道杂化.特别是,这些3d电子占据态几乎全部处于自旋向上通道,自旋向下通道的3d电子态则分布在导带中,形成未占据态.因此, M n原子的磁矩为M F e=[nʏ(x)-nˌ(x)]μB=(5-0)μB=5.0μB.这与计算结果给出的M n原子磁矩5.08μB相符.而间隙位掺杂时,除了F e r m i能级向高能方向移动进入导带外,体系电子结构并没有发生显著变化,M n原子磁矩为5.18μB.3㊀结㊀论使用基于密度泛函理论的第一性原理方法,计算研究了F e掺杂C d T e晶体的晶格结构㊁电子结构等性质变化.结果发现:(1)㊀F e掺杂体系中原子的键长㊁弛豫位置随杂质原子位置(替代位或间隙位)不同而异,但体系的晶格常数都变大.在形成能上,F e原子更倾向于占据替代位.(2)㊀替代位F e掺杂时,体系保持半导体性.F e 原子的自旋向上3d电子态主要以定域态形式分布在价带深处,自旋向下3d电子态则主要分布于导带底形成未占据态.F e原子磁矩为3.76μB.(3)㊀间隙位F e掺杂时,费米能级进入导带中,体系转变为导体.F e原子的自旋向上3d电子态以F e 3dGT e5p轨道杂化态形式存在,自旋向下3d电子态则局域在禁带中.F e原子磁矩为3.14μB.F e掺杂C d T e体系的电子结构表现出对掺杂位置的依赖性,这与M n掺杂体系明显不同.在下一步的研究中,需要对掺杂体系中F eGF e原子间的耦合作用进行研究[19],从而对(C d,F e)T e作为稀磁半导体的实用性建立更充分的理论基础.参考文献:[1]㊀O h n o H.M a k i n g n o n m a g n e t i c s e m i c o n d u c t o r s f e r r o m a gGn e t i c[J].S c i e n c e,1998,281(5379):951G956.[2]㊀O h n o H,S h e nA,M a t s u k u r aF,e t a l.(G a,M n)A s:a n e wd i l u t e d m a g n e t i cs e m i c o n d u c t o rb a s e do nG a A s[J].A p p l P h y sL e t t,1996,69(3):363G365.[3]㊀J a m e tM,B a r s k iA,D e V i l l e r sT,e t a l.H i g hGC u r i eGt e mGp e r a t u r e f e r r o m a g n e t i s m i n s e l fGo r g a n i z e d G e1Gx M n xn a n o c o l u m n s[J].N a t u r e M a t e r i a l s,2006,5(8):653G659.[4]㊀S a t o K,B e r g q v i s tL,K u d r n o v s k y J,e ta l.F i r s tGp r i n c iGp l e st h e o r y o fd i l u t e m a g n e t i cs e m i c o n d u c t o r s[J].R e vM o dP h y s,2010,82(2):1633G1690.[5]㊀B o n a n n iA.F e r r o m a g n e t i cn i t r i d eGb a s e ds e m i c o n d u c t o r sd o pe dw i t ht r a n s i t i o na n dr a r ee a r t h s[J].S e m i c o n dS c iT e c h n o l,2007,22(9):R41GR56.[6]㊀P e a r t o nSJ,H e o W H,I v i l lM,e ta l.D i l u t e m a g n e t i c93020李昭宁等:F e掺杂C d T e电磁性质的第一性原理研究s e m i c o n d u t i n g o x i d e s[J].S e m i c o n dS c iT e c h n o l,2004,19(10):R59GR74.[7]㊀F u r d y n aJK,S a m a r t hN.M a g n e t i c p r o p e r t i e so f d i l u t e d m a g n e t i cs e m i c o n d u c t o r s:ar e v i e w[J].J A p p lP h y s,1987,61(8):3526G3531.[8]㊀F u r d y n a JK.D i l u t e dm a g n e t i c s e m i c o n d u c t o r s[J].JA pGp l P h y s,1988,64(4):R29GR64.[9]㊀M i r b tS,S a n y a lB,M o h nP.M a g n e t i c p r o p e r t i e so f3di m p u r i t i e ss u b s t i t u t e di n G a A s[J].J P h y s:C o n d e n sM a t t e r,2002,14(12):3295G3302.[10]㊀L u o X,M a r t i n R M.J a h nGT e l l e rd i s t o r t i o na n df e r r oGm a g n e t i s m i n t h e d i l u t e m a g n e t i c s e m i c o n d u c t o r sG a A sʒM na n dc u b i c G a NʒM n[J].P h y s R e v B,2005,72(3):035212.[11]㊀S a t oK,K a t a y a m aGY o s h i d aH.F i r s t p r i n c i p l e sm a t e r i a l sd e s i g n f o r s e m i c o n d u c t o r s p i n t r o n i c s[J].S e m i c o n dS c iT e c h n o l,2002,17(4):367G376.[12]㊀C l a r kSJ,S e g a l lM D,P i c k a r dCJ,e t a l.F i r s t p r i n c iGp l e sm e t h o d su s i n g C A S T E P[J].Z e i t s c h r i f tFürK r iGs t a l l o g r a p h i e,2005,220(5G6):567G570.[13]㊀P e r d e wJP,B u r k eK,E r n z e r h o fM.G e n e r a l i z e d g r a d iGe n ta p p r o x i m a t i o n m a d es i m p l e[J].P h y s R e v L e t t,1996,77(18):3865G3868.[14]㊀V a n d e r b i l tD.S o f ts e l fGc o n s i s t e n t p s e u d o p o t e n t i a l s i nag e n e r a l i z e d e i g e n v a l u ef o r m a l i s m[J].P h y s R e v B,1990,41(11):R7892GR7895.[15]㊀M o n k h o r s tH J,P a c kJD.S p e c i a l p o i n t s f o rb r i l l o u i nGz o n e i n t e g r a t i o n s[J].P h y sR e vB,1976,13(12):5188G5192.[16]㊀D u d a r e vSL,B o t t o n G A,S a v r a s o vS Y,e ta l.E l e cGt r o nGe n e r g yGl o s ss p e c t r aa n dt h es t r u c t u r a ls t a b i l i t y o fn i c k l eo x i d e:a n L S D A+U s t u d y[J].P h y s R e v B,1998,57(3):1505G1509.[17]㊀Küw e nF,L e i t s m a n nR,B e c h s t e d t F.M n a n dF e d o p i n g o fb u l k S i:c o n c e n t r a t i o ni n f l u e n c e o n e l e c t r o n i ca n dm a g n e t i c p r o p e r t i e s[J].P h y s R e v B,2009,80(4):045203.[18]㊀A l Z a h r a n iAZ,S r i v a s t a v aGP,G a r g R,e t a l.A n a b i nGi t i o s t u d y o f e l e c t r o n i c a n d s t r u c t u r a l p r o p e r t i e s o fM n i na G a A se n v i r o n m e n t[J].J P h y s:C o n d e n s M a t t e r,2009,21(48):485504.[19]㊀L iZ h a o n i n g,W a n g W e i p i n g,Z h a oS h i c a o.M nGM ni nGt e r a c t i o ni n d u c e d m e t a l l i c o ri n s u l a t i n g c h a r a c t e r o fd o pe ds i l i c o n:a n a b i n i t i o s t u d y[J].C o m p u tM a t e r S c i,2015,97(1):186G192.F i r s tGp r i n c i p l e s s t u d y o f e l e c t r o n i c a n dm a g n e t i c p r o p e r t i e s o fF e d o p i n g i nC d T e c r y s t a lL I Z h a o n i n g1,2,Z H A O X i n c a i1,WA N G W e i p i n g1,2(1.I n s t i t u t e o f F l u i dP h y s i c s,C h i n aA c a d e m y o fE n g i n e e r i n g P h y s i c s,M i a n y a n g621999,C h i n a;2.T h eK e y L a b o r a t o r y o f S c i e n c e a n dT e c h n o l o g y o nH i g hE n e r g y L a s e r,C h i n aA c a d e m y o fE n g i n e e r i n g P h y s i c s,M i a n y a n g621999,C h i n a)A b s t r a c t:T h e s t r u c t u r a l,e l e c t r o n i ca n d m a g n e t i c p r o p e r t i e so fF e i m p u r i t y d o p i n g i nC d T ea r es t u d i e du s i n g f i r s t p r i n c i p l e s b a s e do n t h e d e n s i t y f u n c t i o n a l t h e o r y(D F T)i n a s u p e r c e l l a p p r o a c h.I no r d e r t om a k e a n a c c uGr a t ed e s c r i p t i o n o f t r a n s i t i o n a t o mF e,a g e n e r a l i z e d g r a d i e n t a p p r o x i m a t i o n(G G A)+U m e t h o d i s u s e d.T h e r eGs u l t s s h o wt h a t t h e r e a r e s i g n i f i c a n t c h a n g e s i n t h e b o n d l e n g t h s,m a g n e t i z a t i o n a n d e l e c t r o n i c p r o p e r t i e s.B a s e d o ne l e c t r o n i c b a n d s t r u c t u r e a n d e l e c t r o nd e n s i t y o f s t a t e s,v a r i e s i n e l e c t r o n i c s t r u c t u r e s i n d u c e db y o r b i t a l h yGb r i d i z a t i o nw i t h r e s p e c t t o d i f f e r e n t d o p i n g s i t e t y p e s(s u b s t i t u t i o n a l o r i n t e r s t i t i a l)a r e a n a l y s e d.I t i s f o u n d t h a t t h e l o c a lm a g n e t i cm o m e n t s o f s u b s t i t u t i o n a l a n d i n t e r s t i t i a lF e i m p u r i t i e s a r e3.76a n d3.14B o h rm a g n e t o n s, r e s p e c t i v e l y.T h e s e c o n c l u s i o n s p r o v i d es o m eu s e f u l r e f e r e n c e t ou n d e r s t a n dt h e p h y s i c a lm e c h a n i s m o fC d T e b a s e dd i l u t e dm a g n e t i c s e m i c o n d u c t o r s.K e y w o r d s:C d T e;F e i m p u r i t i e s;e l e c t r o n i c s t r u c t u r e;f i r s tGp r i n c i p l e s s t u d y040202016年第2期(47)卷。