1、第 38 卷哈尔滨师范大学自然科学学报Vol 38，No 6 2022第 6 期NATUAL SCIENCE JOUNAL OF HABIN NOMAL UNIVESITY单卟啉铁分子自旋输运性质第一性原理研究*吴慧洁，尹海涛(哈尔滨师范大学，光电带隙材料教育部重点实验)【摘要】分子器件是指构建在单个分子或有限个分子上的具有特定功能的器件，它是电子器件向小体积化发展的极限，而它也被证明将会是未来电子元器件的替代者 利用密度泛函理论与非平衡态格林函数相结合的方法，研究了单卟啉铁分子与金纳米线电极构成的分子器件的自旋输运的性质 研究发现:该分子器件中存在自旋相干的量子现象;在非平衡态下，两种自旋态
2、电流都随偏压的增大而非线性增大，自旋过滤效应在偏压的调控实现正负的转化【关键词】自旋输运;自旋过滤;量子相干【中图分类号】O469【文献标识码】A【文章编号】1000 5617(2022)06 0062 05收稿日期:2022 09 13*基金项目:黑龙江自然基金资助项目(LH2020A013)通讯作者0引言随着器件尺度进入到亚 3 nm 尺度以下，量子效应越发突显，传统的硅基器件即将到达极限 寻找新的材料来代替传统硅基材料，发展处理速度更快、更强大的新一代微电子器件是必然趋势 基于分子的器件由于其电子传输的特性，被作为取代传统电子器件的最佳候选器件，而分子自旋电子器件因其自旋输运特性，具有更
3、低的能耗、更快的信息处理速度以及更高的集成度等传统电子器件所不具备的一系列优点而受到广泛关注1 除此之外，分子器件具有不同于传统电子器件的诸多依赖电子自旋的独特性质，例如自旋开关2、自旋过滤3、自旋阀4、自旋整流5、自旋 ND 效应6 等等，基于此许多具有各种自旋输运性质的分子器件，如自旋存储器件7、自旋过滤器件等也已被理论模拟和实验制备出来分子器件中，处于器件中间区域的分子是器件的关键组成部分，它在很大程度上决定着分子器件整体的电学性质 寻找合适的中心分子构建分子器件，探究分子电子器件的电子输运性质，已经成为分子电子器件的研究热点之一 卟啉分子因具有环型共轭平面结构，而且内环中心可与金属配位
4、生成金属卟啉而广受研究者青睐 目前，对于卟啉类分子器件的自旋输运性质研究的工作已经进行了很多 例如，利用非平衡格林函数结合密度泛函理论，曾晶、陈克求等学者研究了由两个锰卟啉分子与对苯 乙烯基连接而成的锰卟啉基自旋电子器件的自旋输运性质，在该装置中可以观察到有趣的自旋滤波和磁电阻效应8;赵朋、刘德胜等学者研究了夹在两个扶椅式单壁碳纳米管电极之间的 Fe P2TA 的自旋极化输运性质，该器件通过借助磁场调制，呈现高效磁电阻、自旋滤波和低偏压负微分电阻效应9;Bhattacharyya 等学者研究了一种基于卟啉纳米带的自旋滤波器件，通过对其输运性质的计算，自旋极化率可达到 100%等10 由于卟啉分
5、第 6 期单卟啉铁分子自旋输运性质第一性原理研究子器件还处于发展阶段，具有很好的发展前景，该文将采用密度泛函理论(DFT)和非平衡格林函数(NEGF)相结合的方法，以单卟啉铁分子器件为例，研究磁性分子器件中自旋输运性质1模型及其计算图 1 为单卟啉铁分子器件模型图，模型分成3 个部分:左电极，中心散射区，右电极 其中，左、右电极由金纳米线构成，中心区包含了卟啉铁分子和一部分金电极 为了增强电极和分子之间的耦合强度，用硫原子代替原来的氢原子，以使电极和分子形成强的共价键，有利于电子输运 在对模型优化后，采用基于密度泛函理论(DFT)和非平衡态格林函数(NEGF)相结合方法的 Nanodcal 第
6、一性原理软件研究自旋相关的输运性质11 计算中采用双 极化的原子轨道为基组，交 换 关 联 泛 函 采 用 的 广 义 梯 度 近 似(GGA)12，在自洽过程中，实空间格点取的是 166，190，330，由于器件没有周期性，k 点采样为 1，1，1 器件的温度设置为 100 K图 1单卟啉铁分子器件模型示意图通过分子器件的电流可以采用 Landauer Bttiker 公式计算13:I(V)=ehLdET(E，V)fL(E，L)f(E，)(1)式中，e是电子电荷，h是普朗克常数，fL/是左右电极的费米 狄拉克分布函数，L和 为左右极的化学势，其与施加在器件左右电极上的偏压V 满足关系:L=e
7、V(2)T(E，V)为透射系数，其表示为T Tr Im(rL)GrIm(r)Ga(3)rL(r)为左(右)电极的延迟自能矩阵，表示中心区与左、右电极的作用强度，其中 Gr(Ga)为延迟(超前)格林函数，定义为:Gr=E+S H rL=r1(4)其中 E+=lima0(E+ia)，Ga与 Gr满足关系 Ga=Gr+S 是中心散射区的重叠矩阵，H 是哈密顿量自旋过滤效应(SFE)定义为:I I/(I+I)，当 I大于 I，SFE 为正值，反之为负 由于在零偏压时没有电流通过器件，可以使用费米能级处的透射系数值代替电流计算此状态下的自旋过滤效应2结果与讨论首先，研究了卟啉分子器件在平衡态下的性质 图
8、2(a)和2(b)分别为器件在平衡态下的透射谱(Transmission)以及投影到卟啉铁分子的态密度(PDOS)在平衡态下，自旋向上的 PDOS峰值分别位于 0 9、0 67、0 05 eV 附近，在态密度峰值对应的能量处，电子透射谱也呈现峰值;自旋 向 下 的 PDOS 峰 值 分 别 位 于 0 89、0 65、0 16、0、0 11、0 29 eV 附近，在此附近电子透射谱也同样呈现峰值，这是由于这些状态为电子隧穿提供了相应的通道 有趣的是自旋向36哈尔滨师范大学自然科学学报2022 年 第 38 卷下量子态中出现了 Fano 共振现象 该共振通常来源于连续态(或扩展态)和局域态之间的
9、相干 在研究的结构中，可看作强耦合量子态(相当于扩展态)和弱耦合量子态(相当于局域态)之间的量子干涉形成的 强耦合量子态被视为参考通道，而弱耦合量子态被看做为离散的能级，当离散能级嵌入参考通道内，它们之间的量子干涉就会产生 Fano 效应 主要原因是参考通道的相位变化很小，而离散能级的相位变化了14 为了更直观的说明，图3 给出了图2(a)中标注 A 和 B峰处能量的透射本征态，即散射态 对于 A 峰对应能量的散射态，分子和电极部分散射态值最大，表明这些量子态和电极具有强耦合特性;而对于 B 峰能量的散射态，其最大值在分子上，电极和分子之间呈现的是弱耦合态 这里需要注意的是，对应 A 和 B
10、峰能量处，散射态是简并的，共有 8 个不同的散射态 由于第 8 散射态对输运起主要贡献，其他态对输运贡献不大，这里给出了第5 个 散射态和第 8 个散射态，其他散射态和第5 个散射态无太大差别图 2平衡态下的自旋相关性质(a)平衡态下的透射谱;(b)投影态密度图 3来自左电极的入射电子的散射态(a)、(c)分别为透射谱中能量位置 A 的第 5 散射态和第 8 散射态;(b)、(d)分别为透射谱中能量位置 B 的第 5 散射态和第 8 散射态在分析平衡态自旋相关的性质后，进一步研究器件在非平衡态下性质 图 4(a)和(b)分别为器件在0 1 V偏压下的自旋相关电流以及自旋过滤效应 计算结果表明两
11、种自旋态的电流均随偏压的增大而非线性增加，在 0 0 2 V 偏压范围内，自旋向下电流大于自旋向上电流，在0 2 1 V偏压范围内自旋向上电流大于自旋向下电流 在 0 0 2 V 偏压范围内由于自旋向下电流大于自旋向上电流，故其 SFE 为负值，且在0 0 1 V 自旋向下电流增大的趋势一直大于自旋向上电流的增大趋势，导致在这个偏压范围内的 SFE 一直在减小;在0 1 0 2 V 偏压范围内46第 6 期单卟啉铁分子自旋输运性质第一性原理研究自旋向下电流增大趋势减小，自旋向上电流增大趋势增大，二者电流差值逐渐减小，导致在此范围内SFE有所增大;在0 2 1 V偏压范围内，自旋向上电流大于自旋
12、向下电流，故在此范围内的SFE 为正值，且自旋向上电流整体增大趋势大于自旋向下电流增大趋势，这也导致在这个偏压范围内的 SFE 整体呈增大趋势 上述结果表明偏压对器件的自旋过滤效应起到了一定的调控效果，对于器件的自旋过滤效应，不仅可以调控其大小，还可以调其正负图 4分子器件的自旋输运性质(a)电流 电压曲线;(b)自旋过滤效应图 5不同偏压下的自旋相关透射谱(a)自旋向上态电子透射谱;(b)自旋向下态电子透射谱为了解释以上特性，计算了不同偏压下的自旋相关透射系数，其中黑色虚线围成的部分代表偏压窗口，红色实线代表进入偏压窗口的透射谱，如图 5 所示 对于器件的两种自旋态，随着偏压的增大，进入偏压
13、窗口的透射峰均在逐渐的增大，从而导致了两种自旋态的电流在偏压范围内逐渐的增大 从图中还可以看出，在 0 1 V 偏压下，进入偏压窗口的透射峰自旋向下大于自旋向上，在 0 2 V 偏压以后，进入偏压窗口的透射峰自旋向上超过自旋向下，这导致了在 0 0 2 V偏压下，自旋向下电流大于自旋向上电流，在0 2 1 V 偏压下，自旋向上电流大于自旋向下电流3结束语该文利用 DFT NEGF 方法研究了单卟啉铁分子器件的自旋输运性质 研究发现:该分子器件中存在自旋相关的 Fano 共振现象;器件的两种自旋态电流均随着偏压的增大而非线性增大;自旋过滤效应在偏压的调控下可以实现正负56哈尔滨师范大学自然科学学
14、报2022 年 第 38 卷的转换，在 0 1V 偏压下自旋过滤效应达到最大值为 10%，在 0 2V 偏压下自旋过滤效应最弱为0 9%;通过不同偏压下的透射谱分析得知，由于器件透射峰随着偏压而移动，进入偏压窗口内的透射峰发生变化从而影响了器件电流的变化参考文献 1韩伟，吴镝，罗锋，等 走向自旋的未来信息时代 J 中国科学:物理学 力学 天文学，2022，52(06):5 15 2Zhang Z，Wang Y，Wang H，et al Controllable spin switc-hing in a single molecule magnetic tunneling junction J
15、Nanoscale esearch Letters，2021，16(1):1 8 3陈伟 单分子功能器件设计的第一性原理研究 D 南京邮电大学，2017 4谷现荣，郭立丹，秦阳，等 基于分子自旋阀的自旋电子学J 科学通报，2018，63(35):3689 3696 5Seeger L，Garcia W J S，Dugato D A，et al Wireless pow-er transfer exploring spin rectification and inverse spin HalleffectsJ Journal of Physics D:Applied Physics，2018，5
16、1(16):165002 6Wu X，Xiong L，Feng Y，et al The half metallicity and thespin filtering，ND and spin Seebeck effects in 2D Ag doped SnSe2monolayer J The Journal of chemical phys-ics，2019，150(6):064701 7袁野 基于新型磁性材料的存储器件输运特性研究D 南京邮电大学，2021 8Zeng J，Chen K Q Spin filtering，magnetic and electronicswitching behaviors in manganese porphyrin based spintronicdevices J Journal of Materials Chemistry C，2013，1(25):4014 4019 9Zhao P，Liu D S，Chen G Spin polarized transport proper-ties of Fe oligoporphyrin dimer bas