搜索

x

留言板

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

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

缺陷铁纳米环体系的磁特性研究

叶晴莹 王文静 邓楚楚 陈水源 张鑫源 王雅婧 黄秋怡 黄志高

引用本文:
Citation:

缺陷铁纳米环体系的磁特性研究

叶晴莹, 王文静, 邓楚楚, 陈水源, 张鑫源, 王雅婧, 黄秋怡, 黄志高

Magnetic dynamic properties of defective iron nanorings

Ye Qing-Ying, Wang Wen-Jing, Deng Chu-Chu, Chen Shui-Yuan, Zhang Xin-Yuan, Wang Ya-Jing, Huang Qiu-Yi, Huang Zhi-Gao
PDF
HTML
导出引用
  • 采用Monte Carlo方法与快速傅里叶变换微磁学方法相结合的方式, 模拟含不同缺陷的铁纳米环的磁滞回线、组态、剩磁等磁特性. 研究发现: 缺陷的大小与位置明显影响系统的磁化过程. 当缺陷较小时, 系统存在双稳态特征, 此性质与无缺陷系统类似; 当缺陷增大时, 系统过渡状态增加, 双稳态特征不再明显. 进一步的研究发现, 缺陷系统的剩磁随缺陷半径D的增大而增大. 上述结果与非对称纳米环系统的磁特性类似, 并可以通过零场状态下的系统自旋组态的变化加以解释. 当系统圆心与缺陷中心的间距Y增加时, 剩磁与Y的关系是非线性的: 剩磁先随Y的增大而增大, 后随Y的增大而减小. 模拟结果可用零场状态下不同Y值的组态变化进行详细解释. 上述研究结果表明, 缺陷可以明显影响铁纳米环的磁特性.
    Magnetic nanorings can be high-density integrated because their stray field is low in vortex states. In this paper, the magnetic dynamic properties of the defective Fe nanorings are studied. For convenience, we assume the defect to be round in shape, whose coordinate is (0, Y). Based on the Monte Carlo method and fast Fourier transformation micromagnetism method, the magnetic properties of the defective Fe nanorings, such as hysteresis loops, spin configurations, remanence, etc., are studied. The simulation results indicate that the magnetization process of the system can be affected by the sizes and locations of the defects. When the defects are small, the system has a bistable state, which is similar to the system without defects. The transition state of the system increases as the defects are enlarged, and the bistable state will be no longer so visible. The system becomes open when the defects are big enough. Meanwhile, its hysteresis loop presents a rectangular shape which is similar to cluster’s or quantum dot’s. The remanence increases with the radius of defect increasing. These results are in accord with the magnetic properties of asymmetric magnetic nanoring. In order to explain the above results, the spin configurations of the system are shown. The spins of defective nanorings are divided into two parts, i.e., upper half part and lower half part, which are represented as blue and black spins respectively. When the system does not have any defects, the number of blue spins is equal to black spins’. Therefore the remanence is zero when the system is in a vortex state. It is found that the number of blue spins decreases as the radius of defect increases. This situation results in the total magnetic moment increasing, which leads the remanence to increase. However, the relationship between remanence and Y (the distance between center of nanoring and center of defect) is nonlinear. The remanence first increases and then decreases with Y increasing. The simulation results can be explained by changing the spin configuration. By analyzing the spins of the upper and lower part, the magnetic moment of the system is analyzed. It is found that the number of the spins and the local vortexes can affect the remanence significantly. The results show that the magnetic properties of Fe nanorings can be affected by the defect.
      通信作者: 叶晴莹, qyye@fjnu.edu.cn ; 黄志高, zghuang@fjnu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 61574037)和福建省自然科学基金(批准号: 2017J01553, 2016J01007)资助的课题.
      Corresponding author: Ye Qing-Ying, qyye@fjnu.edu.cn ; Huang Zhi-Gao, zghuang@fjnu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 61574037) and the Natural Science Foundation of Fujian Province, China (Grant Nos. 2017J01553, 2016J01007).
    [1]

    Wang Y G, Zhou K J, Huang G, Hensley C, Huang X N, Ma X P, Zhao T, Baran D S, Ralph J D, Gao J M 2014 Nat. Mater. 13 204Google Scholar

    [2]

    Li H, Cao Z M, Lin J Y, Zhao H, Jiang Q R, Jiang Z Y, Liao H G, Qin K, Xie Z X 2018 Nanoscale 10 1930Google Scholar

    [3]

    Kim D, Lee D R, Choi Y, Metlushko V 2012 Appl. Phys. Lett. 101 192404Google Scholar

    [4]

    王早, 张国峰, 李斌, 陈瑞云, 秦成兵, 肖连团, 贾锁堂 2015 物理学报 64 247803Google Scholar

    Wang Z, Zhang G F, Li B, Chen R Y, Qin C B, Xiao L T, Jia S T 2015 Acta Phys. Sin. 64 247803Google Scholar

    [5]

    Fatemi M, Mollania N, Momeni-Moghaddam M, Sadeghifar F 2018 J. Biotechnol. 270 1Google Scholar

    [6]

    李建昌, 邵思佳 2017 物理学报 66 017101Google Scholar

    Li J C, Shao S J 2017 Acta Phys. Sin. 66 017101Google Scholar

    [7]

    何学敏, 钟伟, 都有为 2018 物理学报 67 227501Google Scholar

    He X M, Zhong W, Du Y W 2018 Acta Phys. Sin. 67 227501Google Scholar

    [8]

    林枝钦 2009 硕士学位论文(福州: 福建师范大学)

    Lin Z Q 2009 M. S. Thesis (Fzhou: Fujian Normal University) (in Chinese)

    [9]

    Yoon S, Lee S H, Kwak, Nam C, Kim W B, Cho B K 2014 J. Appl. Phys. 115 17B511Google Scholar

    [10]

    Liang Y Z, Li L M, Lu M D, Yuan H Z 2018 Nanoscale 10 548Google Scholar

    [11]

    Parkinson P, Kamonsutthipaijit N, Anderson H L, Herz L M 2016 ACS Nano 10 5933Google Scholar

    [12]

    Yannouleas C, Romanovsky I, Landman U 2015 J. Phys. Chem. C 119 11131Google Scholar

    [13]

    张中月, 孙中华, 王红艳, 张志东 2011 物理学报 60 047808Google Scholar

    Zhang Z Y, Sun Z H, Wang Y H, Zhang Z D 2011 Acta Phys. Sin. 60 047808Google Scholar

    [14]

    王同标, 刘念华, 于天宝, 徐旭明, 廖清华 2014 物理学报 63 017301Google Scholar

    Wang T B, Liu N H, Yu T B, Xu X M, Liao Q H 2014 Acta Phys. Sin. 63 017301Google Scholar

    [15]

    吕江涛, 王凤文, 马振鹤, 司光远 2013 物理学报 62 057804Google Scholar

    Lü J T, Wang F W, Ma Z H, Si G Y 2013 Acta Phys. Sin. 62 057804Google Scholar

    [16]

    Chen X, Qin J, Han X F, Liu Y 2018 Appl. Phys. Lett. 113 142406Google Scholar

    [17]

    Liu H, Wei H, Han X F, Yu G, Zhan W, Gall S, Lu Y, Hehn M, Mangin S, Sun M, Liu Y H, Cheng H 2018 Phys. Rev. Appl. 10 054013Google Scholar

    [18]

    Singh N, Goolaup S, Tan W, Adeyeye A O, Balasubramaniam N 2007 Phys. Rev. B 75 104407Google Scholar

    [19]

    Palma J L, Morales-Concha C, Leighton B, Escrig D J, Altbir J 2012 J. Magn. Magn. Mater. 324 637Google Scholar

    [20]

    Avila J I, Tumelero M A, Pasa A A, Viegas A D C 2015 J. Appl. Phys. 117 103901Google Scholar

    [21]

    Zhu F Q, Chern G W, Tchernyshyov O, Zhu X C, Zhu J G, Chien C L 2006 Phys. Rev. Lett. 96 027205Google Scholar

    [22]

    钟克华, 冯倩, 翁臻臻, 黄志高 2005 计算物理 22 534Google Scholar

    Zhong K H, Feng Q, Weng Z Z, Huang Z G 2005 Chin. J. Comput. Phys. 22 534Google Scholar

    [23]

    Huang Z G, Chen Z G, Peng K, Wang D H, Zhang W Y, Zhang F M, Du Y W 2004 Phys. Rev. B 69 094420Google Scholar

    [24]

    Huang Z, Chen Z, Zhang F, Du Y 2004 Eur. Phys. J. B 37 177

    [25]

    Huang Z, Chen Z, Li S, Feng Q, Zhang F, Du Y 2006 Eur. Phys. J. B 51 65Google Scholar

    [26]

    Ye Q, Feng Q, Chen S, Zhang J, Huang Z 2009 J. Nanosci. Nanotechnol. 9 1635Google Scholar

  • 图 1  缺陷铁纳米环模型

    Fig. 1.  Sketch map of defective Fe nanoring.

    图 2  不同D值的纳米环磁滞回线(Y = 30 nm, R = 100 nm, r = 40 nm)

    Fig. 2.  Hysteresis loops of defective Fe nanorings with different D (Y = 30 nm, R = 100 nm, r = 40 nm).

    图 3  不同Y值条件下系统剩磁随D值的变化(R = 100 nm, r = 40 nm)

    Fig. 3.  The relation between the remanence and D with different Y (R = 100 nm, r = 40 nm).

    图 4  不同D值的铁纳米环零场下的自旋组态图(Y = 30 nm, R = 100 nm, r = 40 nm)

    Fig. 4.  The spin configurations of Fe nanorings for different D with zero field (Y = 30 nm, R = 100 nm, r = 40 nm)

    图 5  剩磁随缺陷Y值变化曲线(D = 30 nm, R = 100 nm, r = 40 nm)

    Fig. 5.  The relation between the remanence and Y (D = 30 nm, R = 100 nm, r = 40 nm).

    图 6  不同Y值的铁纳米环的自旋组态图(D = 30 nm, R = 100 nm, r = 40 nm)

    Fig. 6.  The spin configurations of Fe nanorings for different Y (D = 30 nm, R = 100 nm, r = 40 nm).

  • [1]

    Wang Y G, Zhou K J, Huang G, Hensley C, Huang X N, Ma X P, Zhao T, Baran D S, Ralph J D, Gao J M 2014 Nat. Mater. 13 204Google Scholar

    [2]

    Li H, Cao Z M, Lin J Y, Zhao H, Jiang Q R, Jiang Z Y, Liao H G, Qin K, Xie Z X 2018 Nanoscale 10 1930Google Scholar

    [3]

    Kim D, Lee D R, Choi Y, Metlushko V 2012 Appl. Phys. Lett. 101 192404Google Scholar

    [4]

    王早, 张国峰, 李斌, 陈瑞云, 秦成兵, 肖连团, 贾锁堂 2015 物理学报 64 247803Google Scholar

    Wang Z, Zhang G F, Li B, Chen R Y, Qin C B, Xiao L T, Jia S T 2015 Acta Phys. Sin. 64 247803Google Scholar

    [5]

    Fatemi M, Mollania N, Momeni-Moghaddam M, Sadeghifar F 2018 J. Biotechnol. 270 1Google Scholar

    [6]

    李建昌, 邵思佳 2017 物理学报 66 017101Google Scholar

    Li J C, Shao S J 2017 Acta Phys. Sin. 66 017101Google Scholar

    [7]

    何学敏, 钟伟, 都有为 2018 物理学报 67 227501Google Scholar

    He X M, Zhong W, Du Y W 2018 Acta Phys. Sin. 67 227501Google Scholar

    [8]

    林枝钦 2009 硕士学位论文(福州: 福建师范大学)

    Lin Z Q 2009 M. S. Thesis (Fzhou: Fujian Normal University) (in Chinese)

    [9]

    Yoon S, Lee S H, Kwak, Nam C, Kim W B, Cho B K 2014 J. Appl. Phys. 115 17B511Google Scholar

    [10]

    Liang Y Z, Li L M, Lu M D, Yuan H Z 2018 Nanoscale 10 548Google Scholar

    [11]

    Parkinson P, Kamonsutthipaijit N, Anderson H L, Herz L M 2016 ACS Nano 10 5933Google Scholar

    [12]

    Yannouleas C, Romanovsky I, Landman U 2015 J. Phys. Chem. C 119 11131Google Scholar

    [13]

    张中月, 孙中华, 王红艳, 张志东 2011 物理学报 60 047808Google Scholar

    Zhang Z Y, Sun Z H, Wang Y H, Zhang Z D 2011 Acta Phys. Sin. 60 047808Google Scholar

    [14]

    王同标, 刘念华, 于天宝, 徐旭明, 廖清华 2014 物理学报 63 017301Google Scholar

    Wang T B, Liu N H, Yu T B, Xu X M, Liao Q H 2014 Acta Phys. Sin. 63 017301Google Scholar

    [15]

    吕江涛, 王凤文, 马振鹤, 司光远 2013 物理学报 62 057804Google Scholar

    Lü J T, Wang F W, Ma Z H, Si G Y 2013 Acta Phys. Sin. 62 057804Google Scholar

    [16]

    Chen X, Qin J, Han X F, Liu Y 2018 Appl. Phys. Lett. 113 142406Google Scholar

    [17]

    Liu H, Wei H, Han X F, Yu G, Zhan W, Gall S, Lu Y, Hehn M, Mangin S, Sun M, Liu Y H, Cheng H 2018 Phys. Rev. Appl. 10 054013Google Scholar

    [18]

    Singh N, Goolaup S, Tan W, Adeyeye A O, Balasubramaniam N 2007 Phys. Rev. B 75 104407Google Scholar

    [19]

    Palma J L, Morales-Concha C, Leighton B, Escrig D J, Altbir J 2012 J. Magn. Magn. Mater. 324 637Google Scholar

    [20]

    Avila J I, Tumelero M A, Pasa A A, Viegas A D C 2015 J. Appl. Phys. 117 103901Google Scholar

    [21]

    Zhu F Q, Chern G W, Tchernyshyov O, Zhu X C, Zhu J G, Chien C L 2006 Phys. Rev. Lett. 96 027205Google Scholar

    [22]

    钟克华, 冯倩, 翁臻臻, 黄志高 2005 计算物理 22 534Google Scholar

    Zhong K H, Feng Q, Weng Z Z, Huang Z G 2005 Chin. J. Comput. Phys. 22 534Google Scholar

    [23]

    Huang Z G, Chen Z G, Peng K, Wang D H, Zhang W Y, Zhang F M, Du Y W 2004 Phys. Rev. B 69 094420Google Scholar

    [24]

    Huang Z, Chen Z, Zhang F, Du Y 2004 Eur. Phys. J. B 37 177

    [25]

    Huang Z, Chen Z, Li S, Feng Q, Zhang F, Du Y 2006 Eur. Phys. J. B 51 65Google Scholar

    [26]

    Ye Q, Feng Q, Chen S, Zhang J, Huang Z 2009 J. Nanosci. Nanotechnol. 9 1635Google Scholar

  • [1] 胡志良, 周斌, 曾智蓉, 梁天骄. 粒子(E45 MeV)核内级联Monte Carlo模拟程序研究. 物理学报, 2016, 65(23): 232501. doi: 10.7498/aps.65.232501
    [2] 李少波, 殷春浩, 徐振坤, 李佩欣, 吴彩平, 冯铭扬. 基于电子顺磁共振的锶铁氧体磁特性研究. 物理学报, 2015, 64(10): 107502. doi: 10.7498/aps.64.107502
    [3] 刘凤金, 陈水源, 黄志高. Ba掺杂及工艺对BiFeO3体系结构和磁特性的影响. 物理学报, 2014, 63(8): 085101. doi: 10.7498/aps.63.085101
    [4] 何永周, 周巧根. 上海光源低温波荡器永磁铁在低温下的磁特性研究. 物理学报, 2013, 62(4): 044106. doi: 10.7498/aps.62.044106
    [5] 周飞, 丁天怀. 散射介质中层间杂质检测效率的影响因素及分析. 物理学报, 2010, 59(12): 8451-8458. doi: 10.7498/aps.59.8451
    [6] 余波, 应阳君, 许海波. 惯性约束聚变的中子半影成像诊断系统的优化研究. 物理学报, 2010, 59(6): 4100-4109. doi: 10.7498/aps.59.4100
    [7] 李天富, 陈东风, 王洪立, 孙凯, 刘蕴韬. 超薄Fe(4?)膜磁特性极化中子反射研究. 物理学报, 2009, 58(11): 7993-7997. doi: 10.7498/aps.58.7993
    [8] 宫 野, 张建红, 王晓东, 吴 迪, 刘金远, 刘 悦, 王晓钢, 马腾才. 强流脉冲离子束辐照双层靶能量沉积的数值模拟. 物理学报, 2008, 57(8): 5095-5099. doi: 10.7498/aps.57.5095
    [9] 邱东江, 王 俊, 丁扣宝, 施红军, 郏 寅. 退火对Mn和N共掺杂的Zn0.88Mn0.12O:N薄膜特性的影响. 物理学报, 2008, 57(8): 5249-5255. doi: 10.7498/aps.57.5249
    [10] 高 湉, 曹世勋, 李文娟, 康保娟, 袁淑娟, 张金仓. Cu掺杂LaMn1-xCuxO3体系的磁转变和导电行为研究. 物理学报, 2006, 55(7): 3692-3697. doi: 10.7498/aps.55.3692
    [11] 关治强, 薛岩频, 林 海, 何贵丽, 吴晨旭. 钠离子浓度对核小体纤维结构影响的Monte Carlo模拟. 物理学报, 2006, 55(1): 460-464. doi: 10.7498/aps.55.460
    [12] 康保娟, 曹世勋, 王新燕, 李领伟, 黎文峰, 刘 芬, 曹桂新, 郁黎明, 敬 超, 张金仓. 混合场中 (Pr1-yNdy)2/3Sr1/3MnO3体系磁转变行为研究. 物理学报, 2005, 54(2): 902-906. doi: 10.7498/aps.54.902
    [13] 郑 宏, 王绍青, 成会明. 微孔对单壁纳米碳管储氢性能的影响. 物理学报, 2005, 54(10): 4852-4856. doi: 10.7498/aps.54.4852
    [14] 邵元智, 钟伟荣, 林光明. 三维X-Y模型的滞后标度和动态相变行为. 物理学报, 2003, 52(9): 2309-2313. doi: 10.7498/aps.52.2309
    [15] 郭宝增. 用全带Monte Carlo方法模拟纤锌矿相GaN和ZnO材料的电子输运特性. 物理学报, 2002, 51(10): 2344-2348. doi: 10.7498/aps.51.2344
    [16] 陈敏, 魏合林, 刘祖黎, 姚凯伦. 沉积粒子能量对薄膜早期生长过程的影响. 物理学报, 2001, 50(12): 2446-2451. doi: 10.7498/aps.50.2446
    [17] 杨 宁, 陈光华, 张 阳, 公维宾, 朱鹤孙. 薄膜生长的理论模型与Monte Carlo模拟. 物理学报, 2000, 49(11): 2225-2229. doi: 10.7498/aps.49.2225
    [18] 张国民, 杨传章. 铁磁键稀疏Blume-Capel模型相图的Monte Carlo研究. 物理学报, 1993, 42(1): 128-133. doi: 10.7498/aps.42.128
    [19] 应和平, 季达人. 一种有效的量子Monte Carlo模拟方法及其关于二维反铁磁Heisenberg模型的研究. 物理学报, 1993, 42(11): 1845-1850. doi: 10.7498/aps.42.1845
    [20] 何延才, 曹立群. Monte Carlo方法计算微颗粒的X射线强度. 物理学报, 1984, 33(2): 241-249. doi: 10.7498/aps.33.241
计量
  • 文章访问数:  3255
  • PDF下载量:  14
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-12-26
  • 修回日期:  2019-03-25
  • 上网日期:  2019-05-01
  • 刊出日期:  2019-05-20

缺陷铁纳米环体系的磁特性研究

  • 1. 福建师范大学物理与能源学院, 福建省量子调控与新能源材料重点实验室, 福州 350117
  • 2. 福建省半导体光电材料及其高效转换器件协同创新中心, 厦门 361005
  • 通信作者: 叶晴莹, qyye@fjnu.edu.cn ; 黄志高, zghuang@fjnu.edu.cn
    基金项目: 国家自然科学基金(批准号: 61574037)和福建省自然科学基金(批准号: 2017J01553, 2016J01007)资助的课题.

摘要: 采用Monte Carlo方法与快速傅里叶变换微磁学方法相结合的方式, 模拟含不同缺陷的铁纳米环的磁滞回线、组态、剩磁等磁特性. 研究发现: 缺陷的大小与位置明显影响系统的磁化过程. 当缺陷较小时, 系统存在双稳态特征, 此性质与无缺陷系统类似; 当缺陷增大时, 系统过渡状态增加, 双稳态特征不再明显. 进一步的研究发现, 缺陷系统的剩磁随缺陷半径D的增大而增大. 上述结果与非对称纳米环系统的磁特性类似, 并可以通过零场状态下的系统自旋组态的变化加以解释. 当系统圆心与缺陷中心的间距Y增加时, 剩磁与Y的关系是非线性的: 剩磁先随Y的增大而增大, 后随Y的增大而减小. 模拟结果可用零场状态下不同Y值的组态变化进行详细解释. 上述研究结果表明, 缺陷可以明显影响铁纳米环的磁特性.

English Abstract

    • 磁性纳米材料因具有超顺磁性和磁热效应等特性吸引了许多科研工作者的注意, 在生物医学、磁记录等领域发挥着重要作用, 可广泛应用于磁存储、核磁共振成像、磁热疗、纳米药物载体等方面[1-7]. 在种类繁多的纳米磁性材料中, 纳米环以其独特的物理、化学性质引起了科研人员的关注[8-15]. 研究发现: 磁性纳米环在“涡旋状态”时几乎没有杂散场, 因此它们能够被密集地排列而不互相影响, 从而使得一定空间内可存储的信息量明显增加[8]; 在复合纳米环材料中, 若将复合Co/Pt 双分子层纳米环放置于垂直的碳纳米管内, 系统的磁性质比单纯的Co/Pt薄膜具有更大的矫顽力, 其磁相互作用产生台阶状磁滞回线[9]; 在阵列系统中, 金属纳米环阵列的表面敏感性与体积和表面敏感度之间存在着明显的关系, 可为各种尺寸的生物分子提供一种高性能传感平台[10]; 在纳米环的能量传输研究中发现仿生纳米环复合物的能量传送率独立于纳米环尺度[11]; 在拓扑方法的研究中还得到石墨烯纳米环的磁特性存在着典型的量子化特征的结论[12]; 2018年, 韩秀峰课题组通过微磁学模拟研究了纳米环磁性隧道结的磁特性, 发现纳米环磁性隧道结可作为高存储密度、低功耗的存储单元, 模拟结果与实验结果吻合较好[16,17].

      在研究纳米环磁化动力学的过程中, 科研人员发现纳米环的几何形状对纳米环组态有重要影响[18]. 而在纳米环大规模量产的过程中, 难以避免地出现部分含缺陷的纳米环. 这些缺陷纳米环将使系统的磁特性发生相应的变化[19]. 因此, 有必要对含缺陷的纳米环系统的磁特性进行研究.

      在研究缺陷纳米环的过程中, 大部分文献采用的是实验方法或微磁学方法[19-21], 而少见基于Monte Carlo (MC)方法的研究. 其原因主要在于, 传统的MC程序在模拟偶极能的过程中需要耗费大量的时间, 这一情况直接影响了程序的计算效率, 也将使大体系的计算难以进行. 为了克服这一问题, 本研究在传统的MC方法的基础上, 发展出一种MC方法与快速傅里叶变换微磁学(FFTM)方法相结合的计算方式, 这种方法可大幅提高计算效率, 将计算时间缩减为传统MC方法的10%—20%[22]. 本文在对称铁纳米环的基础上构建缺陷, 采用数值模拟的方式研究缺陷铁纳米环磁化动力学特征.

    • 图1为缺陷铁纳米环模型: 纳米环外径为R, 内径为r, 缺陷为一个半径为D的圆形. 若以纳米环圆心为原点, 缺陷圆心坐标可用坐标表示为(0, Y). 此缺陷纳米环模型中R = 100 nm, r = 40 nm, 厚度t = 10 nm, 设Y值变化范围为0—140 nm, D值变化范围为0—50 nm.

      图  1  缺陷铁纳米环模型

      Figure 1.  Sketch map of defective Fe nanoring.

      模拟中利用Heisenberg (海森伯)模型[8], 将体系的哈密顿量表示为

      $\begin{split}E =\,&\! -\!\sum\limits_{ij} J {{{S}}_i} \!\cdot \! {{{S}}_j}\!+\! D\!\!\left[\!\sum\limits_{ij} {\frac{{{{{S}}_i} \cdot {{{S}}_j}}}{{{{r}}_{ij}^3}}\!-\!3\frac{{({{{S}}_i}\!\cdot\!{{{r}}_{ij}})({{{S}}_j}\!\cdot \!{{{r}}_{ij}})}}{{{{r}}_{ij}^5}}}\!\right]\\& - K{v_0}\sum\limits_i {{{({{{S}}_i} \cdot {{{u}}_i})}^2} - {M_s}{v_0}H\sum\limits_i {{{{S}}_i} \cdot {{h}}} },\end{split}$

      其中, ${{{S}}_i}$是纳米环点内第$i$个原子自旋, $J{ { = }}2A{d_1}{\text{π}}/6$(${d_1}$是小球体直径), $A$是球间的交换相互作用常数; $ D = {({M_{\rm{s}}}v_0)^2}/2$是偶极相互作用常数(其中${M_{\rm{s}}}$是饱和磁化强度); ${{{r}}_{ij}}$是第$i$个自旋与第$j$个自旋之间的位移矢量; $K$是磁各向异性常数, ${v_0} = {\text{π}}d_1^3/6$是小球的体积; ${{{u}}_i}$是第$i$个自旋的磁各向异性易轴的单位矢量; $H$是外加磁场. 对于铁纳米环模拟计算, 采用的参数为$A{{ = }}1.51 \times {10^{ - 11}}{\rm{ J/m}}$, ${M_{\rm{s}}} = 1.7 \times {10^6}{\rm{ A/m}}$, $K = {{4}}{{.8}} \times {{1}}{{{0}}^{{4}}}{\rm{ J/}}{{\rm{m}}^{{3}}}$[8].

      令磁场水平向右, 范围从–1500—1500 Oe, 温度取300 K, 采用MC方法与FFTM方法相结合的方式对系统进行模拟[22,23]. 其中MC方法用于模拟系统的交换能、各向异性能、塞曼能, FFTM方法用于模拟系统的偶极能[22].

    • 图2R = 100 nm, r = 40 nm, Y = 30 nm, D分别为0—70 nm的纳米环的磁滞回线. 值得注意的是, 由于Y = 30 nm, 因此当D = 10 nm(图2(a))时, 系统并无缺陷, 即此时纳米环为对称体系, 其磁滞回线与对称铁磁纳米环磁滞回线一致, 存在典型的双稳态[18]. 从图2中还可以发现, 当D = 20—40 nm, 即Dr时(如图2(b)-(d)), 系统的磁滞回线均呈现出较为明显的台阶, 说明此时系统虽然存在缺陷, 但依然保持着磁性纳米环系统的双稳态特征[16,18]; 随着D的增大, 系统的过渡态逐渐增多; 当D = 50—60 nm, 即D > r时(如图2(e)(f)), 由于存在众多的过渡态, 系统的台阶不再明显, 即双稳态特征已逐渐模糊; 当D = 70 nm时, 系统由于缺陷太大而从上部断开, 此时系统不再是一个圆环状态, 磁滞回线的双稳态特征消失, 其形状接近矩形, 此时系统类似于原子团或量子点的磁化行为[24-26]. 从图2中可见, 随着D值的增加, 系统的剩磁也逐渐增大, 这一结论与不对称磁性纳米环的特征类似[18,19].

      图  2  不同D值的纳米环磁滞回线(Y = 30 nm, R = 100 nm, r = 40 nm)

      Figure 2.  Hysteresis loops of defective Fe nanorings with different D (Y = 30 nm, R = 100 nm, r = 40 nm).

      为了研究系统剩磁与缺陷半径D的关系, 我们模拟了不同Y值条件下系统的剩磁随D值变化的曲线, 如图3所示. 从图3中可以发现, 对于无缺陷系统(Y = 30 nm, D = 0—10 nm), 系统的剩磁接近于零. 而当系统的D值增加时, 不同Y值的系统剩磁都随缺陷半径的增大而明显增大.

      图  3  不同Y值条件下系统剩磁随D值的变化(R = 100 nm, r = 40 nm)

      Figure 3.  The relation between the remanence and D with different Y (R = 100 nm, r = 40 nm).

      为了分析图3中剩磁变化的原因, 我们模拟了Y = 30 nm时, 零场状态下系统的组态, 如图4所示. 为了便于分析, 我们将系统自旋分为上下两个部分, 上半部分自旋用蓝色表示, 下半部分自旋用黑色表示. 从图4可以看出, 当外场为零时, 不同D值的系统均为“涡旋态”. 当D = 0—10 nm时(如图4(a)), 系统处于无缺陷状态, 这时自旋组态的上半部分原子个数与下半部分原子个数相同, 但自旋方向相反(蓝色部分偏向左侧, 黑色部分偏向右侧). 设蓝色部分合磁矩矢量为$\sum {{{{m}}_1}} $, 黑色部分合磁矩矢量为$\sum {{{{m}}_2}} $. 由于$\sum {{{{m}}_1}} $$\sum {{{{m}}_2}} $的方向相反, 则系统合磁矩$\sum {{m}} = {\sum {{m}} _1} + {\sum {{m}} _1} \approx 0$, 因此系统磁化强度接近于零; 随着D值从20 nm逐渐增大至60 nm(如图4(b)(f)), 缺陷面积逐渐增大, 系统蓝色部分原子个数逐渐减少(即$\sum {{{{m}}_1}} $减小), 而黑色部分原子个数不变($\sum {{{{m}}_2}} $基本不变), 因此$\sum {{m}} $增大, 导致系统剩磁增大; 当D = 70 nm时, 系统断开(如图4(g)), 不再构成一个完整的环, 此时由于缺陷面积很大, 蓝色部分原子数远小于黑色部分, 因此$\sum {{m}} $继续增加, 从而出现了图3中剩磁随D值的增大而增大的现象.

      图  4  不同D值的铁纳米环零场下的自旋组态图(Y = 30 nm, R = 100 nm, r = 40 nm)

      Figure 4.  The spin configurations of Fe nanorings for different D with zero field (Y = 30 nm, R = 100 nm, r = 40 nm)

      为了研究系统剩磁与Y值的关系, 我们模拟了D = 30 nm 时, 剩磁随Y值变化曲线, 如图5所示. 从图5中可以看出, 剩磁随着Y值的增加而明显的变化: 当Y为0—10 nm时, 系统剩磁保持不变; 当Y为20—90 nm时, 剩磁随Y的增加而增加, 并在Y = 90 nm时达到最大值; 当Y为90—140 nm时, 剩磁随Y的增加而减小.

      图  5  剩磁随缺陷Y值变化曲线(D = 30 nm, R = 100 nm, r = 40 nm)

      Figure 5.  The relation between the remanence and Y (D = 30 nm, R = 100 nm, r = 40 nm).

      为了解释Y值的增加导致剩磁先增大后减小的现象, 我们模拟了零场时与图5相对应的组态(如图6). 从图6中可以看到, 当Y为0—10 nm时, 系统处于无缺陷状态. 此时纳米环为“涡旋态”. 根据上述分析可知, 此时系统的磁化强度接近于0. 随着Y值增加, 缺陷开始上移(如图6(b)), 此时系统蓝色部分原子个数随着Y的增加而减少(当Y < 70 nm时), 即$\sum {{{{m}}_1}} $减小, 而$\sum {{{{m}}_2}} $几乎不变, 因此系统合磁矩$\sum {{m}} $增加, 从而导致系统剩磁增加.

      图  6  不同Y值的铁纳米环的自旋组态图(D = 30 nm, R = 100 nm, r = 40 nm)

      Figure 6.  The spin configurations of Fe nanorings for different Y (D = 30 nm, R = 100 nm, r = 40 nm).

      按上述理论分析, 系统剩磁的最大值应出现于蓝色部分原子个数最少(缺陷面积最大的位置), 即Y = 70 nm处. 但从图5中可以看到, 系统的剩磁最大值出现于Y = 90 nm处, 为什么?

      分析图6(d)图6(e)可以发现, 当Y为80 nm和90 nm时, 系统的蓝色部分出现了“局部涡旋态”[8]. 根据前面的分析, 处于“涡旋”状态的原子磁矩矢量和接近于零, 因此“局部涡旋态”处原子磁矩矢量和接近于零. 由上述分析可知, 这一情况将导致偏向左侧的蓝色自旋个数减少, 即$\sum {{{{m}}_1}} $减小, 而$\sum {{{{m}}_2}} $并没有明显变化, 因此系统合磁矩$\sum {{m}} $Y = 80 nm处继续增大, 即剩磁继续增加.

      而在Y = 90 nm处, 从图6(e)中可以看出, 局部涡旋的面积增大, 因此$\sum {{{{m}}_1}} $继续减小, 合磁矩$\sum {{m}} $进一步增大, 并在此处达到剩磁最大值.

      随着Y的继续增加, 缺陷继续向上移动, 并逐渐从上方移出纳米环. 在这一过程中, 随着蓝色原子数的增加, $\sum {{{{m}}_1}} $也逐渐增加, 因此$\sum {{m}} $逐渐减少, 系统剩磁亦逐渐减少, 从而产生了图5中系统剩磁随Y值的变化曲线.

    • 通过MC方法与FFTM方法相结合的方式模拟了缺陷铁纳米环的磁特性. 研究发现: 当Dr时, 系统存在双稳态特征, 这与对称纳米环的磁特性类似; 当D > r时, 系统的双稳态特征不再明显, 当系统由于缺陷的增大而出现开口时, 系统的磁滞回线呈矩形; 当D从10—70 nm的变化过程中, 系统的剩磁逐渐增加; 研究系统剩磁与Y值的关系可以发现, 当D = 30 nm, Y = 0—10 nm时, 系统剩磁接近于0, 并保持不变; 当Y = 20—90 nm时, 剩磁随Y的增加而增加; 当Y = 90—140 nm时, 剩磁随Y的增加而减小. 以上各现象均可用系统的自旋组态进行解释.

参考文献 (26)

目录

    /

    返回文章
    返回