搜索

x

留言板

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

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

不同周期结构硅锗超晶格导热性能研究

刘英光 郝将帅 任国梁 张静文

引用本文:
Citation:

不同周期结构硅锗超晶格导热性能研究

刘英光, 郝将帅, 任国梁, 张静文

Thermal conductivities of different period Si/Ge superlattices

Liu Ying-Guang, Hao Jiang-Shuai, Ren Guo-Liang, Zhang Jing-Wen
PDF
HTML
导出引用
  • 构造了均匀、梯度、随机3种不同周期分布的硅/锗(Si/Ge)超晶格结构. 采用非平衡分子动力学(NEMD)方法模拟了硅/锗超晶格在3种不同周期分布下的热导率, 并研究了样本总长度和温度对热导率的影响. 模拟结果表明: 梯度和随机周期Si/Ge超晶格的热导率明显低于均匀周期结构超晶格; 在不同的周期结构下, 声子分别以波动和粒子性质输运为主; 均匀周期超晶格热导率具有显著的尺寸效应和温度效应, 而梯度、随机周期Si/Ge超晶格的热导率对样本总长度和温度的依赖性较小.
    Thermoelectric materials, which can convert wasted heat into electricity, have attracted considerable attention because they provide a solution to energy problems. The Si/Ge superlattices have shown tremendous promise as effective thermoelectric materials. The period lengths of the Si/Ge superlattices can effectively tailor the phonon's transport behaviors and control their thermal conductivities. In this paper, three kinds of Si/Ge superlattices with different period length distributions (uniform, gradient, random) are constructed. The non-equilibrium molecular dynamics (NEMD) method is used to calculate the thermal conductivities of Si/Ge superlattices under the different period length distributions. The effect of the sample’s total length and temperature on the superlattice's thermal conductivity are studied. The simulation result shows that the thermal conductivity of gradient and random periodical Si/Ge superlattices are significantly reduced at room temperature compared with that of the uniform period Si/Ge superlattices. Phonons are transported by wave or particle properties in the different periodical superlattices. The thermal conductivity of uniform period superlattices has an obvious size effect with the increasing of the sample total length. In contrast, the thermal conductivity of gradient, random periodical Si/Ge superlattices are weakly dependent on the sample’s total length. At the same time, temperature is an important factor affecting the heat transport properties. We find that the temperature affects the thermal conductivities of the three kinds of superlattices in different ways. With the increase of the temperature, (i) the thermal conductivity of uniform periodical superlattices shows an obvious temperature effect; (ii) the thermal conductivity of the gradient and random periodical Si/Ge superlattices are nearly unchanged due to the competition between phonon localization weakness and phonon-phonon scattering enhancement. In addition, the phonon densities of states of superlattices with three different periodical length distributions are calculated. We find that in the picture of uniform periodical Si/Ge superlattices, the number of pronounced peaks quickly decreases as the period length increases, particularly at higher frequencies. This indicates that as the period length increases, fewer coherent phonons will be formed over the superlattices. Moreover, the scattering mechanisms of phonons for gradient and random periodical Si/Ge superlattices are basically the same at 100 K and 500 K. These findings provide a developmental way to further reduce the thermal conductivity of superlattices.
      通信作者: 刘英光, liuyingguang@ncepu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 52076080)、河北省自然科学基金(批准号: E2020502011)和中央高校基本科研究业务费(批准号: 2020MS105)资助的课题
      Corresponding author: Liu Ying-Guang, liuyingguang@ncepu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 52076080), the Natural Science Foundation of Hebei Province, China(Grant No. E2020502011), and the Fundamental Research Fund for the Central Universities, China (Grant No. 2020MS105)
    [1]

    Martín-González M, Caballero-Calero O, Díaz-Chao P 2013 Renew. Sust. Energ. Rev. 24 288Google Scholar

    [2]

    Feng T L, Ruan X L, Ye Z, Cao B 2015 Phys. Rev. B 91 224301Google Scholar

    [3]

    Chen Z Y, Wang R F, Wang G Y, Zhou X Y, Wang Z S, Yin C, Hu Q, Zhou B Q, Tang J, Anag R 2018 Chin. Phys. B 27 047202Google Scholar

    [4]

    Wang K X, Wang J, Li Y, Zou T, Wang X H, Li J B, Cao Z, Shi W J, Xinba Y E 2018 Chin. Phys. B 27 048401Google Scholar

    [5]

    郭敬云, 陈少平, 樊文浩, 王雅宁, 吴玉程 2020 物理学报 69 146801Google Scholar

    Guo J Y, Chen S P, Fan W H, Wang Y N, Wu Y C 2020 Acta. Phys. Sin. 69 146801Google Scholar

    [6]

    张玉, 吴立华, 曾李骄开, 刘叶烽, 张继业, 邢娟娟, 骆军 2016 物理学报 65 107201Google Scholar

    Zhang Y, Wu L H, Zengli J K, Liu Y F, Zhang J Y, Xing J J, Luo J 2016 Acta. Phys. Sin 65 107201Google Scholar

    [7]

    Lin K H, Strachan A 2013 Phys. Rev. B 87 115302Google Scholar

    [8]

    Chen Y F, Li D Y, Lukes J R, Ni Z H, Chen M H 2005 Phys. Rev. B 72 174302Google Scholar

    [9]

    Giri A, Hopkins P E, Wessel J G, Duda J C 2015 J. Appl. Phys. 118 165303Google Scholar

    [10]

    Zhou K K, Xu N, Xie G F 2018 Chin. Phys. B 27 026501Google Scholar

    [11]

    Xiong R, Yang C, Wang Q, Zhang Y, Li X 2019 Int. J. Thermophys. 40 86Google Scholar

    [12]

    Zhang C W, Zhou H, Zeng Y, Zheng L, Zhan Y L, Bi K D 2019 Int. J. Heat Mass Transf. 132 681Google Scholar

    [13]

    Samaraweera N, Larkin J M, Chan K L, Mithraratne K 2018 J. Appl. Phys. 123 244303Google Scholar

    [14]

    Juntunen T, Vanska O, Tittonen I 2019 Phys. Rev. Lett. 122 105901Google Scholar

    [15]

    Wang Y, Huang H X, Ruan X L 2014 Phys. Rev. B 90 165406Google Scholar

    [16]

    Plimpton S 1995 J. Comput. Phys. 117 1Google Scholar

    [17]

    Dickey J M, Paskin A 1969 Phys. Rev. 188 1407Google Scholar

    [18]

    Xie G F, Ding D, Zhang G 2018 Adv. Phys.-X 3 1480417Google Scholar

    [19]

    Ravichandran J, Yadav A K, Cheaito R, Rossen P B, Soukiassian A, Suresha S J, Duda J C, Foley B M, Lee C-H, Zhu Y, Lichtenberger A W, Moore J E, Muller D A, Schlom D G, Hopkins P E, Majumdar A, Ramesh R, Zurbuchen M A 2014 Nat. Mater. 13 168Google Scholar

    [20]

    Simkin M V, Mahan G D 2000 Phys. Rev. Lett. 84 927Google Scholar

    [21]

    Maldovan M 2015 Nat. Mater. 14 667Google Scholar

    [22]

    Chernatynskiy A, Grimes R W, Zurbuchen M A, Clarke D R, Phillpot S R 2009 Appl. Phys. Lett. 95 161906Google Scholar

    [23]

    Chen X K, Xie Z X, Zhou W X, Tang L M, Chen K Q 2016 Appl. Phys. Lett. 109 023101Google Scholar

    [24]

    Luckyanova M N, Mendoza J, Lu H, Song B, Huang S, Zhou J, Li M, Dong Y, Zhou H, Garlow J, Wu L, Kirby B J, Grutter A J, Puretzky A A, Zhu Y, Dresselhaus M S, Gossard A, Chen G 2018 Sci. Adv. 338 936Google Scholar

    [25]

    Schelling P K, Phillpot S R, Keblinski P 2002 Phys. Rev. B 65 144306Google Scholar

    [26]

    刘英光, 边永庆, 韩中合 2020 物理学报 69 033101Google Scholar

    Liu Y G, Bian Y Q, Han Z H 2020 Acta. Phys. Sin. 69 033101Google Scholar

    [27]

    Chen J, Zhang G, Li B W 2010 Nano Lett. 10 3978Google Scholar

    [28]

    Liang T, Zhou M, Zhang P, Yuan P, Yang D G 2020 Int. J. Heat Mass Transf. 151 119395Google Scholar

    [29]

    Zhang Z W, Chen Y P, Xie Y, Zhang S B 2016 Appl. Therm. Eng. 102 1075Google Scholar

    [30]

    Wang Y, Vallabhaneni A, Hu J, Qiu B, Chen Y P, Ruan X L 2014 Nano Lett. 14 592Google Scholar

    [31]

    Bodapati A, Schelling P K, Phillpot S R, Keblinski P 2006 Phys. Rev. B 74 245207Google Scholar

  • 图 1  不同周期长度分布的超晶格 (a)均匀周期; (b)梯度周期; (c)随机周期

    Fig. 1.  The different period length distribution of superlattices: (a) Uniform period; (c) gradient period; (c) random period.

    图 2  NEMD模拟计算热性质的示意图

    Fig. 2.  Schematic diagram of the NEMD model for calculating the thermal properties.

    图 3  热导率与周期长度的关系

    Fig. 3.  Thermal conductivity of superlattice as a function of period length.

    图 4  不同平均周期长度的均匀超晶格的声子态密度.

    Fig. 4.  The phonon density of states of uniform superlattices with different average period lengths.

    图 5  热导率与样品总长度的关系

    Fig. 5.  Thermal conductivity of superlattice as a function of sample total length.

    图 6  超晶格热导率随温度的变化

    Fig. 6.  Thermal conductivity of superlattice as a function of temperature.

    图 7  超晶格的声子参与率

    Fig. 7.  The participation ratio of superlattices.

    图 8  超晶格不同温度条件下的声子态密度: (a) 梯度周期超晶格; (b) 随机周期超晶格

    Fig. 8.  The phonon density of states with different temperatures of superlattice: (a) Gradient period superlattice; (b) random.

  • [1]

    Martín-González M, Caballero-Calero O, Díaz-Chao P 2013 Renew. Sust. Energ. Rev. 24 288Google Scholar

    [2]

    Feng T L, Ruan X L, Ye Z, Cao B 2015 Phys. Rev. B 91 224301Google Scholar

    [3]

    Chen Z Y, Wang R F, Wang G Y, Zhou X Y, Wang Z S, Yin C, Hu Q, Zhou B Q, Tang J, Anag R 2018 Chin. Phys. B 27 047202Google Scholar

    [4]

    Wang K X, Wang J, Li Y, Zou T, Wang X H, Li J B, Cao Z, Shi W J, Xinba Y E 2018 Chin. Phys. B 27 048401Google Scholar

    [5]

    郭敬云, 陈少平, 樊文浩, 王雅宁, 吴玉程 2020 物理学报 69 146801Google Scholar

    Guo J Y, Chen S P, Fan W H, Wang Y N, Wu Y C 2020 Acta. Phys. Sin. 69 146801Google Scholar

    [6]

    张玉, 吴立华, 曾李骄开, 刘叶烽, 张继业, 邢娟娟, 骆军 2016 物理学报 65 107201Google Scholar

    Zhang Y, Wu L H, Zengli J K, Liu Y F, Zhang J Y, Xing J J, Luo J 2016 Acta. Phys. Sin 65 107201Google Scholar

    [7]

    Lin K H, Strachan A 2013 Phys. Rev. B 87 115302Google Scholar

    [8]

    Chen Y F, Li D Y, Lukes J R, Ni Z H, Chen M H 2005 Phys. Rev. B 72 174302Google Scholar

    [9]

    Giri A, Hopkins P E, Wessel J G, Duda J C 2015 J. Appl. Phys. 118 165303Google Scholar

    [10]

    Zhou K K, Xu N, Xie G F 2018 Chin. Phys. B 27 026501Google Scholar

    [11]

    Xiong R, Yang C, Wang Q, Zhang Y, Li X 2019 Int. J. Thermophys. 40 86Google Scholar

    [12]

    Zhang C W, Zhou H, Zeng Y, Zheng L, Zhan Y L, Bi K D 2019 Int. J. Heat Mass Transf. 132 681Google Scholar

    [13]

    Samaraweera N, Larkin J M, Chan K L, Mithraratne K 2018 J. Appl. Phys. 123 244303Google Scholar

    [14]

    Juntunen T, Vanska O, Tittonen I 2019 Phys. Rev. Lett. 122 105901Google Scholar

    [15]

    Wang Y, Huang H X, Ruan X L 2014 Phys. Rev. B 90 165406Google Scholar

    [16]

    Plimpton S 1995 J. Comput. Phys. 117 1Google Scholar

    [17]

    Dickey J M, Paskin A 1969 Phys. Rev. 188 1407Google Scholar

    [18]

    Xie G F, Ding D, Zhang G 2018 Adv. Phys.-X 3 1480417Google Scholar

    [19]

    Ravichandran J, Yadav A K, Cheaito R, Rossen P B, Soukiassian A, Suresha S J, Duda J C, Foley B M, Lee C-H, Zhu Y, Lichtenberger A W, Moore J E, Muller D A, Schlom D G, Hopkins P E, Majumdar A, Ramesh R, Zurbuchen M A 2014 Nat. Mater. 13 168Google Scholar

    [20]

    Simkin M V, Mahan G D 2000 Phys. Rev. Lett. 84 927Google Scholar

    [21]

    Maldovan M 2015 Nat. Mater. 14 667Google Scholar

    [22]

    Chernatynskiy A, Grimes R W, Zurbuchen M A, Clarke D R, Phillpot S R 2009 Appl. Phys. Lett. 95 161906Google Scholar

    [23]

    Chen X K, Xie Z X, Zhou W X, Tang L M, Chen K Q 2016 Appl. Phys. Lett. 109 023101Google Scholar

    [24]

    Luckyanova M N, Mendoza J, Lu H, Song B, Huang S, Zhou J, Li M, Dong Y, Zhou H, Garlow J, Wu L, Kirby B J, Grutter A J, Puretzky A A, Zhu Y, Dresselhaus M S, Gossard A, Chen G 2018 Sci. Adv. 338 936Google Scholar

    [25]

    Schelling P K, Phillpot S R, Keblinski P 2002 Phys. Rev. B 65 144306Google Scholar

    [26]

    刘英光, 边永庆, 韩中合 2020 物理学报 69 033101Google Scholar

    Liu Y G, Bian Y Q, Han Z H 2020 Acta. Phys. Sin. 69 033101Google Scholar

    [27]

    Chen J, Zhang G, Li B W 2010 Nano Lett. 10 3978Google Scholar

    [28]

    Liang T, Zhou M, Zhang P, Yuan P, Yang D G 2020 Int. J. Heat Mass Transf. 151 119395Google Scholar

    [29]

    Zhang Z W, Chen Y P, Xie Y, Zhang S B 2016 Appl. Therm. Eng. 102 1075Google Scholar

    [30]

    Wang Y, Vallabhaneni A, Hu J, Qiu B, Chen Y P, Ruan X L 2014 Nano Lett. 14 592Google Scholar

    [31]

    Bodapati A, Schelling P K, Phillpot S R, Keblinski P 2006 Phys. Rev. B 74 245207Google Scholar

  • [1] 刘英光, 薛新强, 张静文, 任国梁. 基于界面原子混合的材料导热性能. 物理学报, 2022, 71(9): 093102. doi: 10.7498/aps.71.20211451
    [2] 刘英光, 任国梁, 郝将帅, 张静文, 薛新强. 含有倾斜界面硅/锗超晶格的导热性能. 物理学报, 2021, 70(11): 113101. doi: 10.7498/aps.70.20201807
    [3] 黄诗浩, 谢文明, 汪涵聪, 林光杨, 王佳琪, 黄巍, 李成. 双能谷效应对N型掺杂Si基Ge材料载流子晶格散射的影响. 物理学报, 2018, 67(4): 040501. doi: 10.7498/aps.67.20171413
    [4] 李柱松, 朱泰山. 超晶格和层状结构传热特性的连续模型及其在能源材料设计中的应用. 物理学报, 2016, 65(11): 116802. doi: 10.7498/aps.65.116802
    [5] 黄丛亮, 冯妍卉, 张欣欣, 李静, 王戈, 侴爱辉. 金属纳米颗粒的热导率. 物理学报, 2013, 62(2): 026501. doi: 10.7498/aps.62.026501
    [6] 周青春, 狄尊燕. 声子对隧穿量子点分子辐射场系统量子相位的影响. 物理学报, 2013, 62(13): 134206. doi: 10.7498/aps.62.134206
    [7] 罗晓华, 何为, 吴木营, 罗诗裕. 准周期激励与应变超晶格的动力学稳定性. 物理学报, 2013, 62(24): 247301. doi: 10.7498/aps.62.247301
    [8] 鲍华. 固体氩的晶格热导率的非简谐晶格动力学计算. 物理学报, 2013, 62(18): 186302. doi: 10.7498/aps.62.186302
    [9] 冯现徉, 逯瑶, 蒋雷, 张国莲, 张昌文, 王培吉. In掺杂ZnO超晶格光学性质的研究. 物理学报, 2012, 61(5): 057101. doi: 10.7498/aps.61.057101
    [10] 蒋雷, 王培吉, 张昌文, 冯现徉, 逯瑶, 张国莲. 超晶格SnO2掺Cr的电子结构和光学性质的研究. 物理学报, 2011, 60(9): 093101. doi: 10.7498/aps.60.093101
    [11] 邓艳平, 吕彬彬, 田强. 非对称方势阱中的激子及其与声子的相互作用. 物理学报, 2010, 59(7): 4961-4966. doi: 10.7498/aps.59.4961
    [12] 高当丽, 张翔宇, 张正龙, 徐良敏, 雷瑜, 郑海荣. 调控声子提高Tm3+掺杂体系的频率上转换荧光. 物理学报, 2009, 58(9): 6108-6112. doi: 10.7498/aps.58.6108
    [13] 丁凌云, 龚中良, 黄平. 声子摩擦能量耗散机理研究. 物理学报, 2009, 58(12): 8522-8528. doi: 10.7498/aps.58.8522
    [14] 李志华, 王文新, 刘林生, 蒋中伟, 高汉超, 周均铭. As保护下的生长中断时间对AlSb/InAs超晶格界面粗糙度的影响. 物理学报, 2007, 56(3): 1785-1789. doi: 10.7498/aps.56.1785
    [15] 邓成良, 邵明珠, 罗诗裕. 带电粒子同超晶格的相互作用与系统的混沌行为. 物理学报, 2006, 55(5): 2422-2426. doi: 10.7498/aps.55.2422
    [16] 夏志林, 范正修, 邵建达. 激光作用下薄膜中的电子-声子散射速率. 物理学报, 2006, 55(6): 3007-3012. doi: 10.7498/aps.55.3007
    [17] 姚 鸣, 朱卡的, 袁晓忠, 蒋逸文, 吴卓杰. 声子辅助的电磁感应透明和超慢光效应的研究. 物理学报, 2006, 55(4): 1769-1773. doi: 10.7498/aps.55.1769
    [18] 顾培夫, 陈海星, 秦小芸, 刘 旭. 基于薄膜光子晶体超晶格理论的偏振带通滤波器. 物理学报, 2005, 54(2): 773-776. doi: 10.7498/aps.54.773
    [19] 徐 权, 田 强. 一维分子链中激子与声子的相互作用和呼吸子解 . 物理学报, 2004, 53(9): 2811-2815. doi: 10.7498/aps.53.2811
    [20] 张启义, 田强. 超晶格中电场单极畴与偶极畴的形成和输运. 物理学报, 2002, 51(8): 1804-1807. doi: 10.7498/aps.51.1804
计量
  • 文章访问数:  1803
  • PDF下载量:  46
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-10-28
  • 修回日期:  2020-11-16
  • 上网日期:  2021-03-29
  • 刊出日期:  2021-04-05

不同周期结构硅锗超晶格导热性能研究

    基金项目: 国家自然科学基金(批准号: 52076080)、河北省自然科学基金(批准号: E2020502011)和中央高校基本科研究业务费(批准号: 2020MS105)资助的课题

摘要: 构造了均匀、梯度、随机3种不同周期分布的硅/锗(Si/Ge)超晶格结构. 采用非平衡分子动力学(NEMD)方法模拟了硅/锗超晶格在3种不同周期分布下的热导率, 并研究了样本总长度和温度对热导率的影响. 模拟结果表明: 梯度和随机周期Si/Ge超晶格的热导率明显低于均匀周期结构超晶格; 在不同的周期结构下, 声子分别以波动和粒子性质输运为主; 均匀周期超晶格热导率具有显著的尺寸效应和温度效应, 而梯度、随机周期Si/Ge超晶格的热导率对样本总长度和温度的依赖性较小.

English Abstract

    • 目前, 人类对能源的利用效率较低, 很大一部分能源以废热的形式被浪费掉. 热电技术能够将废热直接转换成电能, 促进废热的二次利用, 是提高能源利用效率的有效手段[1-5]. 热电材料的热电效率由无量纲优值${\rm{ZT}} = {{\left( {{S^2}\sigma T} \right)} / \kappa }$来评估, 其中, S是塞贝克系数, $\sigma $是电导率, $\kappa $是热导率, T是绝对温度. 为了提高材料的热电优值, 需要在保持电导率和塞贝克系数的不变的情况下, 降低热导率的值[6]. 在过去几十年中, 基于增加声子散射, 减小声子平均自由程的材料纳米化是增加材料热电优值的有效途径. 超晶格(superlattice)材料是由相同或不同的材料周期性交替排列构建的纳米结构材料, 由于其超低热导率, 被认为是热电领域优秀的候选材料. 更重要的是, 由于半导体超晶格材料可以控制其周期长度进而控制声子的散射机理得到广大研究学者的青睐[7].

      通常人们通过降低超晶格的周期长度(增加其界面密度)来提高声子的界面散射率, 从而降低材料的导热率, 并且研究的周期长度都是均匀分布[8-10]. 最近研究结果表明, 在超晶格界面密度固定的条件下, 改变超晶格的周期长度分布可以进一步降低其热导率. 如Xiong等[11]通过格林函数法研究了非均匀周期分布超晶格结构的声子散射. 结果显示, 由于不同周期长度之间产生的一系列堆叠声子带隙, 除了非常低频率的声子外, 大部分声子都被有效的散射了, 从而使材料的热导率相比均匀周期降低了一个数量级. Zhang等[12]采用分子动力学方法研究了Si/Ge超晶格纳米线的非均匀周期性排列方式对热输运的影响. 他们的研究也表明在室温下, 非周期性Si/Ge超晶格纳米线的热导率会大幅降低, 与周期性Si/Ge超晶格纳米线和纯硅纳米线相比, 分别降低47.4%和4.4%. 同时他们指出周期性分布的改变会带来声子态密度的移位, 从而达到散射不同波长的声子的效果.

      Samaraweera等[13]采用非平衡分子动力学结合晶格动力学, 研究了周期长度随机排列的纳米线的热导率. 研究发现在任何不同几何和物理参数下, 随机周期排列的纳米线的绝热性能都优于均匀排列纳米线; 同时其研究指出, 虽然不同周期排列下声子的运动模式都一样, 但由于不同周期排列下声子的局域化程度不同, 导致了材料的导热性能的差异. Juntunen等[14]使用分子动力学计算了平均周期长度为2.8 nm的非周期超晶格的热导率. 研究表明, 由于声子的安德森局域化效应, 随机性超晶格的热导率比周期性超晶格的热导率低98%. Wang等[15]采用非平衡动力学方法模拟了随机周期长度超晶格结构材料的导热性能. 通过对不同周期分布下声子参与率的计算, 研究指出在均匀周期超晶格中声子的相干输运和非相干输运同时对材料的导热起作用, 但在非均匀周期条件下, 相干性声子会局域化, 仅非相干声子对导热起作用, 从而引起其热导率显著降低.

      从先前的研究可以看出, 非周期性超晶格与周期性超晶格相比, 具有急剧降低晶格热导率的作用, 这对于热电应用而言至关重要. 为了进一步了解不同周期分布对于超晶格热导率的影响, 本文建立了均匀周期、梯度周期、随机周期的Si/Ge超晶格结构, 采用非平衡分子动力学模拟方法研究了3种周期长度分布对材料导热性能的影响. 此外, 我们还研究了样品总长度和温度对超晶格导热性能的影响.

    • 图1, 我们建立了均匀周期(周期长度相同)、梯度周期(周期长度渐变排列)和随机周期(周期长度随机排列)3种Si/Ge超晶格结构. 为了避免因横截面积过小带来热导率误差, 在体系中垂直热流方向的横截面积大小定为4 × 4个单位晶胞. 其中, 均匀周期超晶格的周期长度定义为一对连续的Si层和Ge层的总厚度, 而对于梯度周期超晶格和随机周期超晶格, 周期计算为整个结构中Si层和Ge层的平均长度之和. 本文采用的纳米尺度传热试样的原理模型如图2所示, 超晶格被放置在两边的两个热浴区域之间, 模型的最外层设置了绝热壁, 它的作用是减少导热层内的粒子蒸发, 防止与外界产生热量交换, 并且设定该区域粒子的速率为0. 各个方向均采用周期边界条件, 其中, 热流传输方向为Z方向.

      图  1  不同周期长度分布的超晶格 (a)均匀周期; (b)梯度周期; (c)随机周期

      Figure 1.  The different period length distribution of superlattices: (a) Uniform period; (c) gradient period; (c) random period.

      图  2  NEMD模拟计算热性质的示意图

      Figure 2.  Schematic diagram of the NEMD model for calculating the thermal properties.

      所有的模拟过程都是基于LAMMPS软件进行的[16]. 在模拟过程中, Tersoff势函数被用于Si-Si, Ge-Ge和Si-Ge原子之间的相互作用. 我们将时间步长设置为0.001 ps. 首先, 在零温的条件下对体系进行能量子最小化, 优化原子位置; 然后, 将体系放置在NVT系综下进行1000 ps, 使系统在期望温度下的能量分布状态达到均匀; 最后, 使用朗之万恒温器在体系中产生一个温度梯度, 放置在微正则系综(NVE)下获得热流和计算热导率. 热通量(${J_z}$)的计算方法如下:

      ${J_z} = {{\left( {\frac{{\Delta Q}}{{\Delta t}}} \right)} \bigg/ A},$

      式中, A为横截面积, $\Delta Q$为模拟过程的能量差, $\Delta t$为模拟过程的时间间隔.

      温度梯度$\left(\dfrac{{\partial T}}{{\partial {\rm{z}}}}\right)$采用拟合温度分布的线性部分来计算. 温度分布通过沿热传输方向截断原子的速度来得到, 如(2)式:

      $T = {{\left\langle {\sum\limits_{i = 1}^N {{m_i}{\nu _i} \cdot {\nu _i}} } \right\rangle } \Bigg/ {(3N{k_{\rm{B}}})}},$

      式中, N为体系划分的每一部分原子数.

      最后, 基于傅里叶定律的热导率表达式为

      $\kappa = - {{{J_z}} \bigg/ {\left( {\frac{{\partial T}}{{\partial z}}} \right)}}.$

      周期长度的分布对声子的传输有很大的影响. 为了说明声子在3种Si/Ge超晶格中的传输机制, 通过对速度自相关函数进行快速傅里叶变换计算出了声子态密度(PDOS)[17]:

      $D\left( \omega \right) = \int {{{\rm{e}}^{ - {\rm{i}}\omega t}}\frac{1}{N}} \sum\limits_{i = 1}^N {\left\langle {\nu \left( t \right)\nu \left( 0 \right)} \right\rangle } {\rm{d}}t,$

      式中, $D(\omega )$表示频率为$\omega $时的PDOS, $\nu \left( t \right)$为原子在时间t时的速度, $\nu \left( 0 \right)$为原子初速度.

    • 通常, 声子在纳米结构材料中以两种不同方式进行输运, 第一种是粒子性质的非相干输运, 声子在散射中丢失掉原来的相位信息, 如倒逆过程(Umklapp过程)的声子-声子散射、声子-杂质散射、声子-界面散射等, 这些过程被称为扩散散射或者非相干散射. 另一种是声子波动性质的相干输运, 指声子在输运散射过程中保留其相位信息[18]. 因此, 声子在超晶格中的热传递特性可以用相干输运模式或非相干输运模式来解释.

      我们首先研究了周期长度和热导率的关系. 周期长度从0.5 nm增加至27.15 nm, 样本总长度固定为34.76 nm. 如图3, 均匀周期Si/Ge超晶格随着周期长度的增加呈现非单调趋势. 周期长度从0.5 nm增加至4.34 nm过程中, 热导率急剧减小, 达到一个极小值. 这主要是因为在这个周期长度范围内, 声子的平均自由程(MFP)与周期长度相当, 声子的相干输运占据主导地位. 由于能带折叠, 声子的色散会形成系列微带, 引起声子的群速度降低. 当周期长度增加时, 微带数量相应增加, 声子的群速度降低, 从而引起热导率的降低. 当周期长度大于4.34 nm时, 此时周期长度大于声子的平均自由程, 声子进行非相干性的扩散输运, 即声子主要与周期界面发生扩散散射, 声子的平均自由程与周期长度成正比, 所以热导率随周期长度增加(界面密度降低)而增加. 热导率极小值的存在提供了声子输运由波动性的相干输运向粒子性的非相干输运转变的证据[19-22]. 这一模拟结果与Ravichandran等[19]实验观察到的热导率对周期长度的依赖曲线一致. 为了进一步说明这一现象, 我们计算了不同周期长度的均匀超晶格的PDOS. 从图4可以看出, 随着周期长度的增加, 显著峰的数量迅速减少, 特别是在较高频率下峰的数量减少尤为明显, 明显的峰在PDOS表明相应的振动模式具有相干的声子性质. 这表明, 随着周期长度的增加, 在超晶格上形成的相干声子会减少, 这与图3所示的结果一致[23].

      图  3  热导率与周期长度的关系

      Figure 3.  Thermal conductivity of superlattice as a function of period length.

      图  4  不同平均周期长度的均匀超晶格的声子态密度.

      Figure 4.  The phonon density of states of uniform superlattices with different average period lengths.

      然而, 对于梯度、随机周期Si/Ge超晶格与均匀周期Si/Ge超晶格热导率趋势是不同的. 他们的热导率值始终低于均匀周期Si/Ge超晶格, 同时, 随周期长度的增加呈现单调递增的趋势. 这主要归因于两个方面: 其一, 由于周期性的破坏, 相干效应被抑制, 界面扩散散射在声子热输运中占主导地位. 另一方面, 由于梯度和随机周期Si/Ge超晶格具有不同的周期长度, 涉及的声子波长更多, 导致界面散射剧烈, 进一步是热导率降低[12].

    • 图5显示了不同样品总长度下, 均匀、梯度和随机周期Si/Ge超晶格在300 K温度下热导率的变化, 其中平均周期长度为4.34 nm, 即它对应超晶格的最小热导率值. 从图中我们可以观察到, 随着样本总长度的增加, 均匀周期Si/Ge超晶格的热导率先增加, 最终趋于饱和. 这种情况的产生是由于样品总长度较短时, 长波长的声子不能被这个界面有效地散射. 它们在保持相的同时, 相干地穿过较短的超晶格, 导致导热系数的增加. 然而, 原先的长波长声子在足够长的样品长度下可以看作进行非相干输运, 这些相干的长波长声子最终会被非弹性声子-声子过程散射[24]. 这是热导率分子动力学模拟中典型的尺寸效应[25].

      图  5  热导率与样品总长度的关系

      Figure 5.  Thermal conductivity of superlattice as a function of sample total length.

      对于梯度和随机Si/Ge超晶格, 热导率对样本长度的依赖性有很大的不同. 热导率在20—30 nm之间稍微增加之后, 随体系尺寸的变大, 其值几乎保持不变, 这一结果与先前的研究一致[15]. 值得注意的是, 最短长度样品(20 nm)的导热系数在所有3种超晶格类型之间几乎是相同的. 这表明, 在所有3个系统中的有效MFP在这一点上是由系统大小控制的. 然而, 随着系统规模的增加, 可以明显地看到, 均匀周期Si/Ge超晶格的声子MFP更大, 而梯度和随机周期Si/Ge超晶格的声子MFP基本不变.

    • 图6中显示的是温度对热导率的影响. 体系的周期长度和样本总长度分别被固定为4.34 nm和34.76 nm. 从图中可以看出温度对3种超晶格的影响是不同的. 均匀周期Si/Ge超晶格的热导率值随着温度的升高而降低, 表现出晶体的典型行为. 随着温度的增加, 晶格振动加剧会激发更多高频短波声子, 声子的扩散输运逐渐对导热起主导作用, 同时三声子散射随温度增加而增加, 使声子的平均自由程更小, 从而导致热导率进一步降低. 这也意味着, 界面散射不是整体MFP的重要贡献者, 与系统大小的分析一致[26].

      图  6  超晶格热导率随温度的变化

      Figure 6.  Thermal conductivity of superlattice as a function of temperature.

      与均匀周期Si/Ge超晶格不同的是, 梯度、随机周期Si/Ge超晶格对温度表现出很小的依赖. 这表明, 梯度、随机周期Si/Ge超晶格对所有波长的声子的散射效率更高. 它们的热导率对温度的依赖性与非晶固体一致, 热导率几乎与温度无关. 此外, 热导率对温度的低敏感性也可能与声子局域化有关, 我们通过声子态密度进一步计算了声子参与率(PPR), 如图7所示. 声子参与率是深入了解声子活动的有效方法, 同时在定量描述声子局域化效应得到很好的应用[27-30], 表达式为[31]

      图  7  超晶格的声子参与率

      Figure 7.  The participation ratio of superlattices.

      ${\rm{PPR}} = \frac{1}{N}\frac{{\left( {\displaystyle\sum\nolimits_i {{\rm{PDO}}{{\rm{S}}_i}\left( \omega \right)_{}^2} } \right)_{}^2}}{{\displaystyle\sum\nolimits_i {{\rm{PDO}}{{\rm{S}}_i}\left( \omega \right)^4} }},$

      式中N为所有原子的总数.

      图7中, 我们认为PPR < 0.4时, 意味着声子在梯度、随机周期超晶格中被局域化, 随着温度的升高, 声子的相干性的恶化将削弱局域化机制, 从而导致热导率的增加; 相反, 超晶格声子输运中非相干声子的比重与温度成反比, 这将导致超晶格的热导率下降. 因此, 在两种机制的竞争下, 梯度、随机周期超晶格几乎保持不变. 同时, 图8显示了100 K和500 K时均匀和随机周期Si/Ge超晶格的PDOS. 值得注意的是, PDOS曲线在不同温度下基本保持着一致的趋势, 这与图6中所示的热导率趋势一致.

      图  8  超晶格不同温度条件下的声子态密度: (a) 梯度周期超晶格; (b) 随机周期超晶格

      Figure 8.  The phonon density of states with different temperatures of superlattice: (a) Gradient period superlattice; (b) random.

    • 本文采用NEMD方法研究了均匀、梯度、随机周期Si/Ge超晶格的热输运性质. 系统地分析了周期分布、样本长度、温度对热导率的影响. 模拟结果表明:

      1)改变周期长度分布可以有效降低热导率. 均匀周期Si/Ge超晶格热导率呈现非单调趋势, 最小值的存在表明了声子输运由波动性的相干输运向粒子性非相干输运转变, 声子态密度图证实了这一结论. 梯度、随机周期Si/Ge超晶格热导率随周期长度增加而增加, 其归因于体系周期性被破坏, 涉及声子波长增加, 相干效应被抑制, 界面散射变剧烈, 界面扩散散射在声子输运中占据主导位置.

      2)随着样本总长度的增加, 均匀周期Si/Ge超晶格热导率存在明显的尺寸效应. 而梯度、随机周期Si/Ge超晶格热导率在20—30 nm之间稍微增加之后, 其值几乎保持不变.

      3)温度对3种周期分布的Si/Ge超晶格热导率有不同的影响方式. 均匀周期Si/Ge超晶格中随温度的增加, 声子-声子本征散射增强, 热导率值降低. 而梯度和随机周期的Si/Ge超晶格热导率对温度的低敏感性归因于局域化和界面散射两种机制的竞争.

参考文献 (31)

目录

    /

    返回文章
    返回