搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

双磁性中心内嵌富勒烯Y2C2@C82-C2(1)中的超快自旋动力学行为

黄瑞 李春 金蔚 GeorgiosLefkidis WolfgangHübner

双磁性中心内嵌富勒烯Y2C2@C82-C2(1)中的超快自旋动力学行为

黄瑞, 李春, 金蔚, GeorgiosLefkidis, WolfgangHübner
PDF
HTML
导出引用
导出核心图
  • 自旋翻转和自旋转移是实现基于内嵌富勒体系自旋逻辑功能器件设计的先决条件. 本文以双磁性中心内嵌富勒烯Y2C2@C82-C2(1)体系为例, 采用第一性原理计算方法, 结合Λ进程理论模型和自编的遗传算法程序, 在该内嵌富勒烯体系中分别实现了亚皮秒时间尺度内的自旋翻转和自旋转移过程. 计算结果表明, 优化后的内嵌Y2C2团簇结构和实验得到的各项数据基本吻合, 并且会对外部的C82-C2(1)笼结构产生一定的排斥力, 但由于富勒烯笼状结构具有很强的稳定性, 所以整个体系仍然保持碳笼结构的完整性. 通过对自旋密度分布与激光脉冲作用下自旋期望值演化的具体分析, 经由Λ进程的自旋翻转是基于两个Y元素的整体自旋翻转; 自旋转移则源自两个磁性中心以及碳笼之间在激光脉冲作用下的自旋密度重新分布. 本文结果揭示了Y2C2@C82-C2(1)体系中的超快自旋动力学机理, 可望为基于实际内嵌富勒烯分子的自旋逻辑功能器件设计提供理论依据.
      通信作者: 李春, lichun@nwpu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11572251,11872309,11504223)、陕西省自然科学基础研究计划(批准号: 2017JM1033)、中央高校基本科研业务费(批准号: 3102017JC1003, 3102017JC11001, GK2018-01009)、西北工业大学研究生创意创新种子基金(批准号: ZZ2018011)和德国研究基金(DFG): SFB/TRR 173 “Spin+X”资助的课题.
    [1]

    Li J L, Yang G W 2009 Appl. Phys. Lett. 95 085411

    [2]

    Li J L, Yang G W 2009 J. Phys. Chem. C 113 18292

    [3]

    Wang J, Ma L, Liang Y, Gao M, Wang G 2014 J. Theor. Comput. Chem. 13 162

    [4]

    Cox B J, Thamwattana N, Hill J M 2007 Proc. Math. Phys. Eng. Sci. 463 461

    [5]

    Li C, Liu J, Zhang S, Lefkidis G, Hübner W 2015 Carbon 87 153

    [6]

    Jin P, Hou Q, Tang C, Chen Z 2015 Theor. Chem. Acc. 134 1

    [7]

    Xiao Y, Zhu S E, Liu D J, Suzuki M, Lu X, Wang G W 2014 Angew. Chem. Int. Edit. 126 3050

    [8]

    Zhang N, Zhang Y, Yang M Q, Tang Z R, Xu Y J 2013 J. Catal. 299 210

    [9]

    Ren J M, Subbiah J, Zhang B, Ishitake K, Satoh K, Kamigaito M, Qiao G G, Wong E H, Wong W W 2016 Chem. Commun. 52 3356

    [10]

    Johnston H J, Hutchison G R, Christensen F M, Aschberger K, Stone V 2010 Toxicol. Sci. 114 162

    [11]

    Shu C, Corwin F D, Zhang J, Chen Z, Reid J E, Sun M, Xu W, Sim J H, Wang C, Fatouros P P 2009 Bioconjugate. Chem. 20 1186

    [12]

    Chai Y, Guo T, Jin C, Haufler R E, Chibante L P F, Fure J, Wang L, Alford J M, Smalley R E 1991 J. Phys. Chem. 95 557

    [13]

    Wang C R, Kai T, Dr T T, Yoshida T, Dr Y K, Dr E N, Dr M T, Sakata M, Dr H S 2001 Angew. Chem. Int. Ed. Engl. 40 397

    [14]

    Chen N, Chaur M N, Moore C, Pinzón J R, Valencia R, Rodríguezfortea A, Poblet J M, Echegoyen L 2010 Chem. Commun. 46 4818

    [15]

    Li F F, Chen N, Muletgas M, Triana V, Murillo J, Rodríguezfortea A, Poblet J M, Echegoyen L 2013 Chem. Sci. 4 3404

    [16]

    Jin P, Tang C, Chen Z 2014 Coordin. Chem. Rev. 270-271 89

    [17]

    Dunsch L, Yang S, Zhang L, Svitova A, Oswald S, Popov A A 2010 J. Am. Chem. Soc. 132 5413

    [18]

    Harneit W 2002 Phys. Rev. A 65 184

    [19]

    Benjamin S C, Ardavan A, Briggs G A D, Britz D A, Gunlycke D, Jefferson J, Jones M A G, Leigh D F, Lovett B W, Khlobystov A N 2005 J. Phys-Condens. Matt. 18 1599

    [20]

    Ju C, Suter D, Du J 2011 Phys. Lett. A 375 1441

    [21]

    Beaurepaire E, Merle J, Daunois A, Bigot J 1996 Phys. Rev. Lett. 76 4250

    [22]

    Koopmans B, Ruigrok J J, Longa F D, de Jonge W J 2005 Phys. Rev. Lett. 95 267207

    [23]

    Bigot J Y, Vomir M, Beaurepaire E 2009 Nat. Phys. 5 515

    [24]

    Battiato M, Carva K, Oppeneer P M 2010 Phys. Rev. Lett. 105 027203

    [25]

    Stöhr J, Siegmann H C 2006 Magnetism-From Fundamentals to Nanoscale Dynamics (Berlin Heidelberg: Springer-Verlag) p753

    [26]

    Li C, Liu J, Zhang S, Lefkidis G, Hübner W 2015 IEEE. T. Magn. 51 11

    [27]

    Li C, Liu J, Lefkidis G, Hübner W 2017 Phys. Chem. Chem. Phys. 19 673

    [28]

    Jin F, Tamm N B, Troyanov S I, Yang S 2018 J. Am. Chem. Soc. 140 3496

    [29]

    李春, 杨帆, Georgios Lefkidis, Wolfgang Hübner 2011 物理学报 60 017802

    Li C, Yang F, Lefkidis G, Hübner W 2011 Acta Phys. Sin. 60 017802

    [30]

    Li C, Jin W, Xiang H, Lefkidis G, Hübner W 2011 Phys. Rev. B 84 2250

    [31]

    李春, 张少斌, 金蔚, Georgios Lefkidis, Wolfgang Hübner 2012 物理学报 61 177502

    Li C, Zhang S B, Jin W, Lefkidis G, Hubner W 2012 Acta Phys. Sin. 61 177502

    [32]

    Zhang N, Du H, Chang J, Jin W, Li C, Lefkidis G, Hubner W 2018 Phys. Rev. B 98 104431

    [33]

    Frisch M J, Trucks G W, Schlegel H B, Scuseria G E, Robb M A, Cheeseman J R 2009 Gaussian 09, Revision A.1. (Wallingford: Gaussian Inc.)

    [34]

    Nakatsuji H 1979 Chem. Phys. Lett. 67 329

    [35]

    Hay P J 1985 J. Chem. Phys. 82 299

  • 图 1  双磁性中心内嵌富勒烯Y2C2@C82-C2(1)体系 (a) 内嵌富勒烯Y2C2@C82C2(1); (b) 内嵌的双金属碳化物团簇(Y2C2)

    Fig. 1.  Double-magnetic-center endohedral fullerene Y2C2@C82-C2(1) system: (a) Endohedral fullerene Y2C2@C82C2(1); (b) endohedral bimetallic carbide clusters (Y2C2).

    图 2  激光脉冲作用下, Y2C2@C82-C2(1)体系中的超快自旋翻转过程 (a) 自旋翻转过程中所涉及的初始态(黑色虚线)、最终态(红色实线)和中间态(彩色实线)随时间变化的占据情况; (b) 自旋角动量分量的期望值随时间的变化情况; (c) 实现自旋翻转所用的激光脉冲包络线, 插图为结构优化后的内嵌Y2C2团簇, 碳笼用虚线表示; (d)SAC-CI计算得到的能级图(考虑SOC)

    Fig. 2.  Ultrafast spin-switching scenario achieved in Y2C2@C82-C2(1) system under laser pulses: (a) Time evolution of the occupation of the initial state (dashed black line), final state (solid red line), and the intermediate states (colorized solid lines) involved in spin-switching scenario; (b) variation of the expectation values of the spin angular momentum components along the x, y and z axis with time; (c) laser pulse envelope of the spin-switching scenario, inset represents structurally optimized endohedral Y2C2 cluster, the dashed circle represents the carbon cage; (d)energy levels of Y2C2@C82-C2(1) calculated by SAC-CI (including SOC).

    图 3  激光脉冲作用下, Y2C2@C82-C2(1)体系中的超快自旋转移过程 (a) 自旋转移过程中所涉及的初始态(黑色虚线)、最终态(红色实线)和中间态(彩色实线)随时间变化的占据情况; (b) 自旋角动量分量的期望值随时间的变化情况; (c) 实现自旋转移所用的激光脉冲包络线, 插图为结构优化后的内嵌Y2C2团簇, 碳笼用虚线表示; (d)SAC-CI计算得到的能级图(考虑SOC)

    Fig. 3.  Ultrafast spin-transfer scenario achieved in Y2C2@C82-C2(1) system under laser pulses: (a) Time evolution of the occupation of the initial state (dashed black line), final state (solid red line), and the intermediate states (colorized solid lines) involved in spin-transfer scenario; (b) variation of the expectation values of the spin angular momentum components along the x, y and z axis with time; (c) laser pulse envelope of the spin-transfer scenario, inset represents structurally optimized endohedral Y2C2 cluster, the dashed circle represents the carbon cage; (d) energy levels of Y2C2@C82-C2(1) calculated by SAC-CI (including SOC).

    图 4  Y2C2@C82-C2(1)体系自旋翻转和自旋转移过程中初始态$\left| 2 \right\rangle $和最终态$\left| 11 \right\rangle $的电子吸收光谱图

    Fig. 4.  Electronic absorption spectra of the initial (state $\left| 2 \right\rangle $) and final (state $\left| 11 \right\rangle $) state of the spin-switching and spin-transfer scenarios in Y2C2@C82-C2(1).

    表 1  Y2C2@C82-C2(1)体系中能量最低的五个三重态的自旋密度分布情况(未考虑SOC)

    Table 1.  Spin density distributions of the five lowest triplet states in Y2C2@C82-C2(1) system (without SOC).

    AtomState 1State 2State 3State 4State 5
    Y10.4480.4790.1250.3680.198
    Y20.6060.8310.1480.8010.173
    C (max)0.0580.0460.0590.0420.066
    下载: 导出CSV

    表 2  自旋翻转过程涉及能态的局域化位置、能量和自旋期望值(考虑SOC)

    Table 2.  Localization positions, energies and spin expectation values of the involved states in spin-switching scenario (including SOC).

    StateSpin localizationEnergySxSySzS
    eV
    312.4440.0010.0010.9180.918
    302.443−0.003−0.0030.0000.004
    292.4430.0020.001−0.9180.918
    4Y1 = Y20.863−0.001−0.0430.9620.962
    30.8620.0020.0830.0000.083
    2Y1 = Y20.862−0.001−0.0400.9620.963
    下载: 导出CSV

    表 3  Y2C2@C82-C2(1)体系中自旋密度分布情况(考虑SOC)

    Table 3.  Spin density distribution in Y2C2@C82-C2(1) system (including SOC).

    Energy/eVY1Y2C (max)Total spin density (C atom)
    State $\left| 2 \right\rangle $0.862−0.504−0.680−0.016−1.062
    State $\left| 11 \right\rangle $1.586−0.414−0.900−0.032−0.933
    Difference value0.7240.090−0.2200.129
    下载: 导出CSV
  • [1]

    Li J L, Yang G W 2009 Appl. Phys. Lett. 95 085411

    [2]

    Li J L, Yang G W 2009 J. Phys. Chem. C 113 18292

    [3]

    Wang J, Ma L, Liang Y, Gao M, Wang G 2014 J. Theor. Comput. Chem. 13 162

    [4]

    Cox B J, Thamwattana N, Hill J M 2007 Proc. Math. Phys. Eng. Sci. 463 461

    [5]

    Li C, Liu J, Zhang S, Lefkidis G, Hübner W 2015 Carbon 87 153

    [6]

    Jin P, Hou Q, Tang C, Chen Z 2015 Theor. Chem. Acc. 134 1

    [7]

    Xiao Y, Zhu S E, Liu D J, Suzuki M, Lu X, Wang G W 2014 Angew. Chem. Int. Edit. 126 3050

    [8]

    Zhang N, Zhang Y, Yang M Q, Tang Z R, Xu Y J 2013 J. Catal. 299 210

    [9]

    Ren J M, Subbiah J, Zhang B, Ishitake K, Satoh K, Kamigaito M, Qiao G G, Wong E H, Wong W W 2016 Chem. Commun. 52 3356

    [10]

    Johnston H J, Hutchison G R, Christensen F M, Aschberger K, Stone V 2010 Toxicol. Sci. 114 162

    [11]

    Shu C, Corwin F D, Zhang J, Chen Z, Reid J E, Sun M, Xu W, Sim J H, Wang C, Fatouros P P 2009 Bioconjugate. Chem. 20 1186

    [12]

    Chai Y, Guo T, Jin C, Haufler R E, Chibante L P F, Fure J, Wang L, Alford J M, Smalley R E 1991 J. Phys. Chem. 95 557

    [13]

    Wang C R, Kai T, Dr T T, Yoshida T, Dr Y K, Dr E N, Dr M T, Sakata M, Dr H S 2001 Angew. Chem. Int. Ed. Engl. 40 397

    [14]

    Chen N, Chaur M N, Moore C, Pinzón J R, Valencia R, Rodríguezfortea A, Poblet J M, Echegoyen L 2010 Chem. Commun. 46 4818

    [15]

    Li F F, Chen N, Muletgas M, Triana V, Murillo J, Rodríguezfortea A, Poblet J M, Echegoyen L 2013 Chem. Sci. 4 3404

    [16]

    Jin P, Tang C, Chen Z 2014 Coordin. Chem. Rev. 270-271 89

    [17]

    Dunsch L, Yang S, Zhang L, Svitova A, Oswald S, Popov A A 2010 J. Am. Chem. Soc. 132 5413

    [18]

    Harneit W 2002 Phys. Rev. A 65 184

    [19]

    Benjamin S C, Ardavan A, Briggs G A D, Britz D A, Gunlycke D, Jefferson J, Jones M A G, Leigh D F, Lovett B W, Khlobystov A N 2005 J. Phys-Condens. Matt. 18 1599

    [20]

    Ju C, Suter D, Du J 2011 Phys. Lett. A 375 1441

    [21]

    Beaurepaire E, Merle J, Daunois A, Bigot J 1996 Phys. Rev. Lett. 76 4250

    [22]

    Koopmans B, Ruigrok J J, Longa F D, de Jonge W J 2005 Phys. Rev. Lett. 95 267207

    [23]

    Bigot J Y, Vomir M, Beaurepaire E 2009 Nat. Phys. 5 515

    [24]

    Battiato M, Carva K, Oppeneer P M 2010 Phys. Rev. Lett. 105 027203

    [25]

    Stöhr J, Siegmann H C 2006 Magnetism-From Fundamentals to Nanoscale Dynamics (Berlin Heidelberg: Springer-Verlag) p753

    [26]

    Li C, Liu J, Zhang S, Lefkidis G, Hübner W 2015 IEEE. T. Magn. 51 11

    [27]

    Li C, Liu J, Lefkidis G, Hübner W 2017 Phys. Chem. Chem. Phys. 19 673

    [28]

    Jin F, Tamm N B, Troyanov S I, Yang S 2018 J. Am. Chem. Soc. 140 3496

    [29]

    李春, 杨帆, Georgios Lefkidis, Wolfgang Hübner 2011 物理学报 60 017802

    Li C, Yang F, Lefkidis G, Hübner W 2011 Acta Phys. Sin. 60 017802

    [30]

    Li C, Jin W, Xiang H, Lefkidis G, Hübner W 2011 Phys. Rev. B 84 2250

    [31]

    李春, 张少斌, 金蔚, Georgios Lefkidis, Wolfgang Hübner 2012 物理学报 61 177502

    Li C, Zhang S B, Jin W, Lefkidis G, Hubner W 2012 Acta Phys. Sin. 61 177502

    [32]

    Zhang N, Du H, Chang J, Jin W, Li C, Lefkidis G, Hubner W 2018 Phys. Rev. B 98 104431

    [33]

    Frisch M J, Trucks G W, Schlegel H B, Scuseria G E, Robb M A, Cheeseman J R 2009 Gaussian 09, Revision A.1. (Wallingford: Gaussian Inc.)

    [34]

    Nakatsuji H 1979 Chem. Phys. Lett. 67 329

    [35]

    Hay P J 1985 J. Chem. Phys. 82 299

  • [1] 李春, 张少斌, 金蔚, Georgios Lefkidis, Wolfgang Hübner. 线性磁性分子离子中由激光诱导的超快自旋转移. 物理学报, 2012, 61(17): 177502. doi: 10.7498/aps.61.177502
    [2] 李君, 刘立胜, 徐爽, 张金咏. 单轴压缩下Ti3B4的力学、电学性能及变形机制的第一性原理研究. 物理学报, 2020, 69(4): 043102. doi: 10.7498/aps.69.20191194
    [3] 邓小清, 孙琳, 李春先. 界面铁掺杂锯齿形石墨烯纳米带的自旋输运性能. 物理学报, 2016, 65(6): 068503. doi: 10.7498/aps.65.068503
    [4] Georgios Lefkidis, Wolfgang Hübner, 杨帆, 李春. 磁性纳米结构中由激光引起的超快自旋动力学研究. 物理学报, 2011, 60(1): 017802. doi: 10.7498/aps.60.017802
    [5] 向天, 程亮, 齐静波. 拓扑绝缘体中的超快电荷自旋动力学. 物理学报, 2019, 68(22): 227202. doi: 10.7498/aps.68.20191433
    [6] 田文, 袁鹏飞, 禹卓良, 陶斌凯, 侯森耀, 叶聪, 张振华. 掺杂六角形石墨烯电子输运特性的研究. 物理学报, 2015, 64(4): 046102. doi: 10.7498/aps.64.046102
    [7] 金硕, 孙璐. 带有碳杂质的钨中氢稳定性的第一性原理研究 . 物理学报, 2012, 61(4): 046104. doi: 10.7498/aps.61.046104
    [8] 王欣欣, 张颖, 周洪波, 王金龙. 铌对钨中氦行为影响的第一性原理研究. 物理学报, 2014, 63(4): 046103. doi: 10.7498/aps.63.046103
    [9] 郭志超, 张桐耀, 张靖. 微米气室铯原子自旋噪声谱. 物理学报, 2020, 69(3): 037201. doi: 10.7498/aps.69.20191623
    [10] 刘绍军, 原鹏飞, 徐济安, 祝文军, 经福谦. BeO高压相变和声子谱的第一性原理计算. 物理学报, 2010, 59(12): 8755-8761. doi: 10.7498/aps.59.8755
    [11] 苟秉聪, 顾娟, 王山鹰. Au和3d过渡金属元素混合团簇结构、电子结构和磁性的研究. 物理学报, 2009, 58(5): 3338-3351. doi: 10.7498/aps.58.3338
    [12] 李金, 桂贵, 孙立忠, 钟建新. 单轴大应变下二维六角氮化硼的结构变化. 物理学报, 2010, 59(12): 8820-8828. doi: 10.7498/aps.59.8820
    [13] 徐素芬, 冯世平, 秦吉红. 准一维强关联Zigzag型材料的自旋动力学. 物理学报, 2006, 55(10): 5511-5515. doi: 10.7498/aps.55.5511
    [14] 张淑亭, 孙志, 赵磊. 石墨烯纳米片大自旋特性第一性原理研究. 物理学报, 2018, 67(18): 187102. doi: 10.7498/aps.67.20180867
    [15] 张竹霞, 赵彦亮, 闫新, 韩培德, 许并社, 刘旭光, 郝玉英. 富勒烯衍生物苯基C71-丁酸甲酯的结构和电学性质第一性原理研究. 物理学报, 2009, 58(13): 204-S209. doi: 10.7498/aps.58.204
    [16] 龚新高. 高温及高压下液体镓的结构——第一性原理分子动力学方法研究. 物理学报, 1995, 44(6): 885-896. doi: 10.7498/aps.44.885
    [17] 陈军, 王晓中, 林理彬, 何捷. 第一性原理方法研究He掺杂Al晶界力学性质. 物理学报, 2011, 60(7): 077104. doi: 10.7498/aps.60.077104
    [18] 崔红保, 历长云, 李旭升, 王狂飞, 王海燕. AlAs相变及热动力学性质的第一性原理研究. 物理学报, 2009, 58(8): 5598-5603. doi: 10.7498/aps.58.5598
    [19] 忻晓桂, 陈香, 周晶晶, 施思齐. LiFePO4 晶格动力学性质的第一性原理研究. 物理学报, 2011, 60(2): 028201. doi: 10.7498/aps.60.028201
    [20] 龚新高, 郑庆祺. 镓原子集团结构的第一性原理分子动力学研究. 物理学报, 1993, 42(2): 244-251. doi: 10.7498/aps.42.244
  • 引用本文:
    Citation:
计量
  • 文章访问数:  583
  • PDF下载量:  2
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-10-23
  • 修回日期:  2018-12-07
  • 上网日期:  2019-08-27
  • 刊出日期:  2019-02-01

双磁性中心内嵌富勒烯Y2C2@C82-C2(1)中的超快自旋动力学行为

  • 1. 西北工业大学力学与土木建筑学院, 西安 710072
  • 2. 陕西师范大学物理学与信息技术学院, 西安 710119
  • 3. 德国凯泽斯劳滕工业大学物理系, OPTIMAS研究中心, 凯泽斯劳滕 67653, 德国
  • 通信作者: 李春, lichun@nwpu.edu.cn
    基金项目: 国家自然科学基金(批准号: 11572251,11872309,11504223)、陕西省自然科学基础研究计划(批准号: 2017JM1033)、中央高校基本科研业务费(批准号: 3102017JC1003, 3102017JC11001, GK2018-01009)、西北工业大学研究生创意创新种子基金(批准号: ZZ2018011)和德国研究基金(DFG): SFB/TRR 173 “Spin+X”资助的课题.

摘要: 自旋翻转和自旋转移是实现基于内嵌富勒体系自旋逻辑功能器件设计的先决条件. 本文以双磁性中心内嵌富勒烯Y2C2@C82-C2(1)体系为例, 采用第一性原理计算方法, 结合Λ进程理论模型和自编的遗传算法程序, 在该内嵌富勒烯体系中分别实现了亚皮秒时间尺度内的自旋翻转和自旋转移过程. 计算结果表明, 优化后的内嵌Y2C2团簇结构和实验得到的各项数据基本吻合, 并且会对外部的C82-C2(1)笼结构产生一定的排斥力, 但由于富勒烯笼状结构具有很强的稳定性, 所以整个体系仍然保持碳笼结构的完整性. 通过对自旋密度分布与激光脉冲作用下自旋期望值演化的具体分析, 经由Λ进程的自旋翻转是基于两个Y元素的整体自旋翻转; 自旋转移则源自两个磁性中心以及碳笼之间在激光脉冲作用下的自旋密度重新分布. 本文结果揭示了Y2C2@C82-C2(1)体系中的超快自旋动力学机理, 可望为基于实际内嵌富勒烯分子的自旋逻辑功能器件设计提供理论依据.

English Abstract

    • 内嵌富勒烯由于其独特的空心笼状结构可以内嵌入不同的原子、分子、离子或团簇, 使其不但具有富勒烯碳笼的物理化学性质, 而且还兼具了内嵌团簇的各种优良特性, 在磁学[13]、力学[4,5]、化学[6,7]、半导体材料[8,9]和生物医学[10,11]等领域拥有巨大的潜在应用价值, 被认为是未来纳尺度功能器件中重要的备选材料之一. 1991年, 通过激光蒸发氧化镧/石墨复合棒, 第一个宏观量的内嵌富勒烯La@C82在实验中被合成出来[12], 进而得到了内嵌富勒烯配合物的详尽波谱. 随着科研人员对富勒烯配合物研究的不断深入, 越来越多新内嵌富勒烯结构被发现和实验分离出来[1317], 其内嵌化合物也从最开始的单个原子发展到分子、离子甚至分子团簇. 内嵌富勒烯也由于其内嵌团簇的特殊电磁特性, 在纳米级逻辑功能器件中具有广泛的潜在应用前景. 人们发现内嵌富勒烯可以作为量子信息的载体, 通过对其电荷和自旋的操控, 设计出可扩展量子计算机体系[1820], 从而有望大幅度提高计算设备的运算速度和存储密度.

      目前, 基于电荷自由度的传统电子学仍在遵循着摩尔定律发展. 然而随着集成电路上线路密度的增加, 其差错率和复杂性也将呈指数增长. 由于人们对于高性能计算条件的需求不断提高, 必须考虑采用新技术和新材料来满足不断增长的需求, 因此人们将目光聚集到通过操控自旋自由度来提高计算器件的运算速度和读写速度等关键性能上. 1996年, 法国科学家Beaurepaire等[21]首次发现了铁磁金属中由激光诱导的超快磁化动力学行为, 即通过实验实现了激光诱导的Ni薄膜中的超快退磁过程, 这一现象的发现马上引起了国际上研究人员对于各种光驱动的磁化操控机理的研究热情. 2005年, Koopmans团队[22]提出了声子调制Elliott-Yafet机制, 利用磁光抽运探测技术研究了铁磁镍中电荷和自旋的弛豫过程. 2009年, Bigot研究组[23]通过研究激光脉冲与铁磁金属中的电子自旋耦合作用, 提出了相对论量子电动力学理论. 2010年, Battiato等[24]提出了由飞秒激光诱导退磁导致的自旋极化激发电子的超扩散机制. Stöhr和他的同事[25]通过大量的理论研究, 提出了Λ进程理论模型, 用来描述激光与磁性体系的相互作用, 以及激光诱导的超快光学跃迁过程. 通过对这一理论进行修正与完善, Li等[5,26,27]应用第一性原理方法实现了内嵌磁性元素富勒烯中亚皮秒时间尺度的超快自旋翻转和自旋转移过程. 然而, 之前的研究内容主要涉及内嵌富勒烯的理论预测体系, 并未将现有理论模型应用于实验中已合成的实际内嵌富勒烯分子.

      本文选取双磁性中心内嵌富勒烯Y2C2@C82-C2(1)体系[28]为研究对象, 采用基于Λ进程理论模型的第一性原理计算方法和自编程序, 对该分子体系在磁光耦合作用下的超快自旋动力学行为进行了系统研究, 并实现了在同一种体系中的超快自旋翻转和自旋转移. 本文结果可望为未来纳尺度磁性存储单元和自旋逻辑器件设计提供理论依据和新的思路.

    • 本文实现的超快自旋翻转和自旋转移是基于Λ进程理论模型所得到的. 该理论模型, 即由相干激光场驱动的发生在磁性材料中经由中间态参与的初末态之间的间接光学跃迁机理, 已经在我们前期研究工作中进行了详细的描述[2932]. 其中, 当选取从同一个三重态通过Zeeman作用分裂为具有相反自旋方向的两个简并态分别作为初始态和最终态时, 在特定激光脉冲作用后, 若体系从初始态激发至自旋混合的中间态, 再由中间态返回至最终态, 即认为实现了体系自旋方向的翻转过程, 其翻转速度可以达到亚皮秒级别. 另外, 当初始态和最终态的自旋密度高度局域化在不同的磁性原子上时, 则Λ进程也可以实现体系的自旋转移过程. 具体理论计算方法介绍如下.

      首先, 在无任何外场作用时, 可以利用量子化学计算软件Gaussian 09[33]软件包对体系的非相对论哈密顿量进行求解, 体系的哈密顿量为

      $\begin{split} \!\!\!\!\!{\hat H^{(0)}} =& - \frac{1}{2}\sum\limits_{i = 1}^{{N_{{\rm{el}}}}} {\nabla _i^2} - \sum\limits_{i = 1}^{{N_{{\rm{el}}}}} {\sum\limits_{a = 1}^{{N_{{\rm{at}}}}} {\frac{{{Z_a}}}{{\left| {{{{R}}_{{a}}} - {{{r}}_i}} \right|}}} } \\ & + \sum\limits_{i = 1}^{{N_{{\rm{el}}}}} {\sum\limits_{j = 1}^{{N_{{\rm{el}}}}} {\frac{1}{{\left| {{{{r}}_{{i}}} - {{{r}}_{{j}}}} \right|}}} } + \sum\limits_{a = 1}^{{N_{{\rm{at}}}}} {\sum\limits_{b = 1}^{{N_{{\rm{at}}}}} {\frac{{{Z_a}{Z_b}}}{{\left| {{{{R}}_{{a}}} - {{{R}}_{{b}}}} \right|}}} } , \end{split}$

      (1)式中的NatNel分别表示原子和电子的数目; ${{{R}}_{{a}}}$${{{r}}_{{i}}}$分别为原子和电子的坐标矢量; Za为核电荷数. 由于磁性体系较强的电子相互作用, 要得到较为精确的多体电子基态和激发态, 需要采用高精度的post-Hartree-Fock方法, 本文计算时采用了对称性匹配簇-组态相互作用方法(symmetry adapted cluster-configuration interaction, SAC-CI)[34]. 在Λ进程理论模型中, 考虑到体系自旋转换过程中初始态、中间态和最终态的选取, 需要对体系施加一个合适的外加静磁场和加入自旋轨道耦合效应(spin-orbit coupling, SOC)来实现Zeeman分裂, 进而产生自旋混合态, 相应的哈密顿量为

      ${\hat H^{(1)}} = \sum\limits_{i = 1}^{{N_{{\rm{el}}}}} {\frac{{Z_a^{{\rm{eff}}}}}{{2{c^2}R_i^3}}{\hat{ L}} \cdot {\hat{ S}} + } \sum\limits_{i = 1}^{{N_{{\rm{el}}}}} {{\mu _L}{\hat{ L}} \cdot {{B}}} + \sum\limits_{i = 1}^{{N_{{\rm{el}}}}} {{\mu _S}{\hat{ S}} \cdot {{B}}} ,$

      (2)式中的$\displaystyle\sum\limits_{i = 1}^{{N_{{\rm{el}}}}} {\dfrac{{Z_a^{{\rm{eff}}}}}{{2{c^2}R_i^3}}{\hat{ L}} \cdot {\hat{ S}}} $项为SOC效应项; ${{B}}$表示外加静态磁场; c为光速; $Z_a^{{\rm{eff}}}$表示有效核电荷数; ${\hat{ S}}$${\hat{ L}}$表示自旋角动量和轨道角动量算子; ${\mu _S}$${\mu _L}$分别为自旋角动量和轨道角动量的旋磁比. 然后, 对体系施加激光脉冲作用, 驱动体系各能态随时间演化的自旋动力学. 当施加激光脉冲作用后, 需要引入一个含时的激光矢量势${{{A}}_{{\rm{laser}}}}(t)$, 则此时的哈密顿量为

      ${\hat H^{(2)}}(t) = {\hat{ D}} \cdot {{{A}}_{{\rm{laser}}}}(t) + {\mu _S}{\hat{ S}} \cdot {{{B}}_{{\rm{laser}}}}(t) + {\mu _L}{\hat{ L}} \cdot {{{B}}_{{\rm{laser}}}}(t),$

      其中 ${{{B}}_{{\rm{laser}}}}(t)$为激光脉冲作用下的磁场; ${\hat{ D}}$为电偶极矩算子. 然后再采用含时的微扰理论求解体系的薛定谔方程:

      ${\rm{i}}\hbar \frac{\partial }{{\partial t}}{c_n}(t) = \sum\limits_k {\left\langle {{\psi _n}} \right|{{\hat H}^{(2)}}(t)\left| {{\psi _k}} \right\rangle {{\rm{e}}^{ - \frac{{\rm{i}}}{h}({E_k} - {E_n})t}}{c_k}(t)} , $

      (4)式中, EnEk分别为态$\left| {{{\rm{\psi }}_n}} \right\rangle $$\left| {{{\rm{\psi }}_k}} \right\rangle $的本征能量. 最后, 利用自编的遗传算法程序优化激光参数, 进而更有效地获得期望的自旋功能.

    • 目前, 对于内嵌富勒烯中的超快自旋动力学行为研究大部分仅限于理论预测的单磁性中心内嵌富勒烯体系. 理论上来说, 自旋转移和自旋翻转是实现内嵌富勒体系自旋逻辑功能器件设计必不可少的先决条件. 单磁性中心内嵌富勒烯体系可以实现其自旋翻转操控, 由于这种体系只具有一个磁性中心, 所以单磁性中心体系是无法实现自旋转移的. 而自旋逻辑器件设计中很重要的一个方面是可以将自旋翻转、自旋转移以及各种输入输出参数相结合, 以实现相对复杂的逻辑运算功能. 本文选取双磁性中心内嵌富勒烯Y2C2@C82-C2(1)作为研究对象, 基于Λ进程理论模型和第一性原理计算, 并在前期研究工作的基础上, 集中探讨该分子体系在磁光耦合作用下的自旋动力学行为.

      图  1  双磁性中心内嵌富勒烯Y2C2@C82-C2(1)体系 (a) 内嵌富勒烯Y2C2@C82C2(1); (b) 内嵌的双金属碳化物团簇(Y2C2)

      Figure 1.  Double-magnetic-center endohedral fullerene Y2C2@C82-C2(1) system: (a) Endohedral fullerene Y2C2@C82C2(1); (b) endohedral bimetallic carbide clusters (Y2C2).

      在体系结构选择上, 首先对C82碳笼结构的9种异构体C2(1), Cs(2), C2(3), Cs(4), C2(5), Cs(6), C3v(7), C3v(8), C2v(9)进行结构优化计算并得到相应的稳定结构, 再将蝶形双金属碳化物团簇(Y2C2)内嵌到优化好的C82结构中, 然后再次进行结构优化计算得到了7种稳定的Y2C2@C82异构体C2(1), Cs(2), C2(3), Cs(4), Cs(6), C3v(8), C2v(9), 但是在进行自旋密度计算时只有C2(1)构型中自旋密度是高度局域化在两个Y原子上, 即最适合本文采用的Λ进程理论模型进行分析求解, 因此本文选取的研究对象为Y2C2@C82-C2(1)异构体. 在计算过程中, 为了保证计算效率和计算精度, 在选取基组时, 对于电子结构相对简单的C原子, 选取STO-3G基组, 对于电子结构相对较为复杂的Y原子, 选取带有有效核势(ECP)的LANL2DZ基组[35]. 本文计算中所采用的收敛标准为: 计算得到的每个原子上的力和力的均方根分别小于4.5 × 10−4和3 × 10−4 Hartree/Bohr, 其相应的位移和位移的均方根则分别小于1.8 × 10−3和1.2 × 10−3 Bohr.

      对于优化后的内嵌富勒烯Y2C2@C82-C2(1)结构, 计算得到其外部C82-C2(1)笼结构中C—C共价键平均键长为0.1439 nm, 而单独的C82-C2(1)富勒烯结构中C—C共价键平均键长为0.1432 nm, 说明在Y2C2@C82-C2(1)结构中, 内嵌的Y2C2团簇会对外部的C82-C2(1)笼结构产生一定的排斥力, 但由于富勒烯笼状结构具有很强的稳定性, 所以碳笼的整体结构不会发生很大变化, 只会发生轻微的变形, 并且仍然保持碳笼结构的完整性. 对于内嵌的Y2C2团簇, Y—Y键长为0.3203 nm, C—C键长为0.1260 nm, Y—C平均键长为0.2478 nm, 这些计算结果与前人的理论研究结果(Y—Y键长为0.313 nm, C—C键长为0.127 nm, Y—C平均键长为0.205 nm)[16]吻合, 同时也与实验中测得的结果非常接近(Y—Y键长为0.3620 nm, C—C键长为0.1.24 nm, Y—C平均键长为0.239—0.246 nm)[28].

    • 内嵌金属富勒烯作为自旋逻辑功能器件潜在材料的最大优势是其外部的富勒烯笼状结构可以为内嵌磁性团簇提供一个较为稳定的化学环境, 从而保证内嵌磁性团簇的自旋状态不受外部环境因素的干扰. 但是在实际应用中, 内嵌磁性团簇有可能会与外部富勒烯笼结构发生相互作用, 使自旋密度从磁性团簇分散到外部富勒烯笼结构上, 使得内嵌磁性团簇的自旋密度局域化程度降低, 导致对其自旋操控变得更为困难. 因此, 需要首先对内嵌富勒烯Y2C2@C82-C2(1)进行自旋密度分布计算. 本文采用SAC-CI作用方法, 计算得到了该体系中各个原子上的自旋密度分布情况, 如表1所列. 通过表1可以看出, 在Y2C2@C82-C2(1)结构基态三重态中, 整个体系中的自旋密度极少分布在C原子上, 而三重态$\left| 1 \right\rangle $, $\left| 2 \right\rangle $$\left| 4 \right\rangle $上的自旋密度高度局域化在两个Y原子上, 这样在选取合适的初始态、中间态和最终态的前提下, 可以通过Λ进程, 实现在激光诱导作用下该体系中Y原子上的自旋翻转过程.

      AtomState 1State 2State 3State 4State 5
      Y10.4480.4790.1250.3680.198
      Y20.6060.8310.1480.8010.173
      C (max)0.0580.0460.0590.0420.066

      表 1  Y2C2@C82-C2(1)体系中能量最低的五个三重态的自旋密度分布情况(未考虑SOC)

      Table 1.  Spin density distributions of the five lowest triplet states in Y2C2@C82-C2(1) system (without SOC).

      依据Λ进程理论, 实现整个体系的自旋翻转过程最简单有效的是选择两个加入SOC作用之前简并的、且自旋方向相反的自旋态作为初始态和最终态. 施加SOC和外磁场作用后, 在特定激光脉冲的驱动下, 体系若可以经由中间态完成从初始态到最终态的跃迁, 则实现了自旋翻转过程. 本文通过对Y2C2@C82-C2(1)体系施加一个磁场强度为2.35 T的外加静磁场, 使得该体系中所有的三重态发生Zeeman分裂, 第一个三重态分裂得到两个自旋方向相反的态$\left| 2 \right\rangle $$\left| 4 \right\rangle $, 作为自旋翻转过程的初末态. 在合适激光脉冲作用下, 得到了如图2所示的自旋翻转动力学过程, 其中自旋态$\left| 2 \right\rangle $ (如图2(d)中黑色虚线)自旋方向向下, 其自旋平均值为−0.962, 自旋态$\left| 4 \right\rangle $ (如图2(d)中红色实线)自旋方向向上, 其自旋平均值为0.963. 整个自旋翻转过程所涉及到的中间态包括$\left| 3 \right\rangle $, $\left| 29 \right\rangle $, $\left| 30 \right\rangle $$\left| 31 \right\rangle $, 其中主要参与自旋翻转过程的中间态为能态$\left| 29 \right\rangle $, $\left| 30 \right\rangle $$\left| 31 \right\rangle $(如图2(d)中绿色虚线). 图2(a)展示了该体系中各能态随时间变化的占据情况. 可以看出, 在外加激光脉冲作用下时, 初始态$\left| 2 \right\rangle $的占据率随时间逐渐由1降低为0, 而最终态$\left| 4 \right\rangle $的占据率则由最初的0上升至最终的0.990, 实现了该体系中Y原子上的自旋状态由自旋向上转变为自旋向下的超快自旋翻转过程, 整个翻转过程耗时约1000 fs, 其保真度高达97.8%. 图2(b)展示了x, y, z 3个坐标轴方向, 体系自旋平均值分量随时间的变化情况(蓝色实线为z轴分量, 红色实线为y轴分量, 黑色实线为x轴分量), 其中分量$\langle$Sz$\rangle$主导着整个Λ进程中体系自旋方向的变化状态. 在此进一步对比体系总自旋平均值的变化情况, 在激光脉冲作用之前, 体系总自旋平均值为0.962, 其方向向量为(−0.001, −0.043, 0.962), 激光作用后, 总自旋自旋平均值变为0.963, 其方向向量为(−0.001, −0.040, 0.962), 二者夹角大小为176.80°, 非常接近180°, 进一步地证实了该体系在激光作用下实现了自旋翻转过程. 另外, 如图2(c)所示, 实现超快自旋翻转过程所采用的激光为线偏振光, 激光能量Elaser为1.59 eV, 半高宽FWHM为300.00 fs, 激光方向在球坐标系下的入射角$\theta $$\varphi $分别为142.87°和247.39°, 激光偏振方向与光学平面间的夹角$\gamma $为254.08°.

      StateSpin localizationEnergySxSySzS
      eV
      312.4440.0010.0010.9180.918
      302.443−0.003−0.0030.0000.004
      292.4430.0020.001−0.9180.918
      4Y1 = Y20.863−0.001−0.0430.9620.962
      30.8620.0020.0830.0000.083
      2Y1 = Y20.862−0.001−0.0400.9620.963

      表 2  自旋翻转过程涉及能态的局域化位置、能量和自旋期望值(考虑SOC)

      Table 2.  Localization positions, energies and spin expectation values of the involved states in spin-switching scenario (including SOC).

      图  2  激光脉冲作用下, Y2C2@C82-C2(1)体系中的超快自旋翻转过程 (a) 自旋翻转过程中所涉及的初始态(黑色虚线)、最终态(红色实线)和中间态(彩色实线)随时间变化的占据情况; (b) 自旋角动量分量的期望值随时间的变化情况; (c) 实现自旋翻转所用的激光脉冲包络线, 插图为结构优化后的内嵌Y2C2团簇, 碳笼用虚线表示; (d)SAC-CI计算得到的能级图(考虑SOC)

      Figure 2.  Ultrafast spin-switching scenario achieved in Y2C2@C82-C2(1) system under laser pulses: (a) Time evolution of the occupation of the initial state (dashed black line), final state (solid red line), and the intermediate states (colorized solid lines) involved in spin-switching scenario; (b) variation of the expectation values of the spin angular momentum components along the x, y and z axis with time; (c) laser pulse envelope of the spin-switching scenario, inset represents structurally optimized endohedral Y2C2 cluster, the dashed circle represents the carbon cage; (d)energy levels of Y2C2@C82-C2(1) calculated by SAC-CI (including SOC).

    • 实现自旋转移是在内嵌富勒烯结构中实现自旋逻辑功能的先决条件之一, 然而最近的理论研究表明, 经由Λ进程自旋动力学理论机理, 自旋转移相比于自旋翻转更难实现. 双磁性中心内嵌富勒烯Y2C2@C82-C2(1)体系拥有两个磁性中心Y1和Y2原子, 考虑SOC作用后, 得到了该体系中各个原子上自旋密度重新排布情况(如表3所示). 能态$\left| 2 \right\rangle $上Y1和Y2原子的自旋密度分别为0.504和0.680 (自旋方向向下), 可以看作对等分布. 单个C原子上自旋密度最大值仅为0.016, 84个C原子上总的自旋密度为1.062, 说明自旋密度高度局域化在Y原子上. 而能态$\left| 11 \right\rangle $上Y1和Y2原子的自旋密度分别为0.414和0.900 (自旋向下), 单个C原子上自旋密度最大值仅为0.032, 84个C原子上总的自旋密度为0.933, 可以认为自旋密度局域化在Y2原子上. 根据Λ进程理论, 选取能态$\left| 2 \right\rangle $作为初始态, 能态$\left| 11 \right\rangle $作为最终态.

      Energy/eVY1Y2C (max)Total spin density (C atom)
      State $\left| 2 \right\rangle $0.862−0.504−0.680−0.016−1.062
      State $\left| 11 \right\rangle $1.586−0.414−0.900−0.032−0.933
      Difference value0.7240.090−0.2200.129

      表 3  Y2C2@C82-C2(1)体系中自旋密度分布情况(考虑SOC)

      Table 3.  Spin density distribution in Y2C2@C82-C2(1) system (including SOC).

      图3通过Λ进程实现了初始态$\left| 2 \right\rangle $到最终态$\left| 11 \right\rangle $上的自旋转移过程. 在外加激光脉冲作用下, 能态$\left| 2 \right\rangle $的占据率随时间迅速从1降低为0, 而能态$\left| 11 \right\rangle $的占据率则由0上升至0.963, 整个自旋转移过程耗时约200 fs, 其保真度高达95.1%. 结合表3中自旋密度分布情况, 在整个自旋转移过程中, Y1原子和碳笼上的自旋密度值均有所减小, 而Y2原子上的自旋密度值明显增大, 因此, 可以看作是实现了从Y1原子到Y2原子上的自旋转移过程. 值得注意的是, 整个自旋转移过程涉及到的中间态较多(态$\left| 1 \right\rangle $, $\left| 3 \right\rangle $, $\left| 4 \right\rangle $, $\left| 12 \right\rangle $, $\left| 13 \right\rangle $, $\left| 14 \right\rangle $, $\left| 15 \right\rangle $, $\left| 16 \right\rangle $, $\left| 29 \right\rangle $, $\left| 30 \right\rangle $$\left| 31 \right\rangle $, 如图3(d)中绿色虚线所示), 但实际主要参与的中间态仅有态$\left| 14 \right\rangle $$\left| 29 \right\rangle $. 另外, 如图3(c)所示, 实现超快自旋转移过程所采用的激光为线偏振光, 激光能量Elaser为0.83 eV, 半高宽FWHM为50.75 fs, 激光方向在球坐标系下的入射角$\theta $$\varphi $分别为253.72°和305.81°, 激光偏振方向与光学平面间的夹角$\gamma $为20.41°.

      图  3  激光脉冲作用下, Y2C2@C82-C2(1)体系中的超快自旋转移过程 (a) 自旋转移过程中所涉及的初始态(黑色虚线)、最终态(红色实线)和中间态(彩色实线)随时间变化的占据情况; (b) 自旋角动量分量的期望值随时间的变化情况; (c) 实现自旋转移所用的激光脉冲包络线, 插图为结构优化后的内嵌Y2C2团簇, 碳笼用虚线表示; (d)SAC-CI计算得到的能级图(考虑SOC)

      Figure 3.  Ultrafast spin-transfer scenario achieved in Y2C2@C82-C2(1) system under laser pulses: (a) Time evolution of the occupation of the initial state (dashed black line), final state (solid red line), and the intermediate states (colorized solid lines) involved in spin-transfer scenario; (b) variation of the expectation values of the spin angular momentum components along the x, y and z axis with time; (c) laser pulse envelope of the spin-transfer scenario, inset represents structurally optimized endohedral Y2C2 cluster, the dashed circle represents the carbon cage; (d) energy levels of Y2C2@C82-C2(1) calculated by SAC-CI (including SOC).

      图4展示了Y2C2@C82-C2(1)体系初始态$\left| 2 \right\rangle $和最终态$\left| 11 \right\rangle $的电子吸收光谱图, 横轴表示体系各能态与选取的初始态/最终态之间的能量差值, 纵轴表示振子强度. 电子吸收光谱峰意味着相应态与初始态$\left| 2 \right\rangle $/最终态$\left| 11 \right\rangle $的跃迁概率不为0, 多峰值情况表示有多个可与初末态发生跃迁的电子态存在. 但是, 考虑到实现Λ进程对两分支(初态—中间态跃迁与中间态—末态跃迁)对称性的要求以及激光失谐等因素, 并非所有的峰值所对应的能态都能作为自旋翻转过程的中间态. 对于自旋翻转过程, 图4(a)给出了初始态$\left| 2 \right\rangle $的电子吸收光谱(最终态$\left| 4 \right\rangle $由于与初始态$\left| 2 \right\rangle $从同一个三重态分裂出来, 其吸收光谱基本相同, 故没有显示). 其中蓝色虚线对应的峰值所对应的能态为自旋翻转过程中涉及到的主要中间态$\left| 29 \right\rangle $, $\left| 30 \right\rangle $$\left| 31 \right\rangle $, 与初始态$\left| 2 \right\rangle $之间的能量差值为1.58 eV (其能级位置也可由图2(d)看出). 该能量差值与实现该自旋翻转过程所需激光脉冲的能量 (1.59 eV)基本相同, 表明激光失谐较小. 对于自旋转移过程, 图4(b)给出了最终态$\left| 11 \right\rangle $的电子吸收光谱. 由于选取的初始态$\left| 2 \right\rangle $和最终态$\left| 11 \right\rangle $来自两个不同的三重态, 其光学性质不尽相同. 经分析发现参与该过程的中间态只要有态$\left| 14 \right\rangle $和态$\left| 29 \right\rangle $, 其位置分别由红色虚线(图4(a)能量差0.76 eV处及图4(b)能量差0.03 eV处)和蓝色虚线(图4(a)能量差1.58 eV处及图4(b)能量差0.857 eV处)标注. 其中, 该自旋转移过程所需外加激光脉冲的能量为0.85 eV, 与E$\left| 14 \right\rangle $$\left| 29 \right\rangle $能量差基本吻合, 因此易于实现.

      图  4  Y2C2@C82-C2(1)体系自旋翻转和自旋转移过程中初始态$\left| 2 \right\rangle $和最终态$\left| 11 \right\rangle $的电子吸收光谱图

      Figure 4.  Electronic absorption spectra of the initial (state $\left| 2 \right\rangle $) and final (state $\left| 11 \right\rangle $) state of the spin-switching and spin-transfer scenarios in Y2C2@C82-C2(1).

      通过体系电子吸收光谱可以发现, 自旋翻转和自旋转移过程中涉及到的主要中间态与初末态均有跃迁. 因此, 通过分析初始态对应的电子吸收光谱, 可以预判出自旋翻转和自旋转移所采用的激光脉冲能量及其实现自旋动力学过程的难易程度.

    • 本文以双磁性中心内嵌富勒烯Y2C2@C82-C2(1)体系为例, 采用第一性原理计算方法, 结合Λ进程理论模型和自编的遗传算法程序, 系统地研究了该分子体系在磁光耦合作用下的超快自旋动力学行为. 计算结果表明, 优化后的内嵌Y2C2团簇结构和实验得到的各项数据基本吻合, 并且会对外部的C82-C2(1)笼结构产生一定的排斥力, 但由于富勒烯笼状结构具有很强的稳定性, 所以整个体系仍然保持碳笼结构的完整性. 在Y2C2@C82-C2(1)体系中, 通过计算发现自旋密度高度局域化在两个Y原子上, 极少地分布在碳笼上. 经由Λ进程理论模型, 在同一内嵌富勒烯体系中分别实现自旋翻转和自旋转移过程. 在自旋翻转过程中, 本文实现了该体系中Y原子上的自旋状态由自旋向上转变为自旋向下的超快自旋翻转过程, 整个翻转过程耗时约1000 fs, 其保真度高达97.8%; 同时本文也实现了该体系下的超快自旋转移过程, 整个自旋转移过程耗时仅200 fs, 其保真度高达95.1%. 另外, 在进行电子吸收光谱计算过程中发现, 当跃迁所需能量差值与外加激光脉冲能量的失谐程度越小, 实现自旋翻转和自旋转移过程的概率越大. 本文的理论研究结果可为基于多磁性中心内嵌富勒烯体系的自旋逻辑功能器件设计奠定理论基础并提供新的设计思路.

参考文献 (35)

目录

    /

    返回文章
    返回