搜索

x

留言板

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

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

表面低配位原子对声子的散射机制

霍龙桦 谢国锋

表面低配位原子对声子的散射机制

霍龙桦, 谢国锋
PDF
HTML
导出引用
  • 由于纳米结构具有极高的表体比, 声子-表面散射机制对声子的热输运性质起到关键作用. 提出了表面低配位原子对声子的散射机制, 并且结合量子微扰理论与键序理论推导出该机制的散射率. 由于散射率正比于材料的表体比, 这种散射机制对声子输运的重要性随着纳米结构尺寸的减小而增大. 散射率正比于声子频率的4次方, 所以这种散射机制对高频声子的作用远远强于对低频声子的作用. 基于声子玻尔兹曼输运方程, 计算了硅纳米薄膜和硅纳米线的热导率, 发现本文模型比传统的声子-边界散射模型更接近实验值. 此发现不仅有助于理解声子-表面散射的物理机制, 也有助于应用声子表面工程调控纳米结构的热输运性质.
      通信作者: 谢国锋, xieguofeng@hnust.cn
    • 基金项目: 国家自然科学基金(批准号: 11874145)资助的课题.
    [1]

    Wang J, Xie F, Cao X H, An S C, Zhou W X, Tang L M, Chen K Q 2017 Sci. Rep. 7 41418

    [2]

    Ding Z D, An M, Mo S Q, Yu X X, Jin Z L, Liao Y X, Esfarjani K, Lü J T, Shiomi J, Yang N 2019 J. Mater. Chem. A 7 2114

    [3]

    Tang L P, Tang L M, Geng H, Yi Y P, Wei Z M, Chen K Q, Deng H X 2018 Appl. Phys. Lett. 112 012101

    [4]

    Ouyang T, Jiang E L, Tang C, Li J, He C Y, Zhong J X 2018 J. Mater. Chem. A 6 21532

    [5]

    Xie G F, Ju Z F, Zhou K K, Wei X L, Guo Z X, Cai Y Q, Zhang G 2018 npj Comput. Mater. 4 21

    [6]

    Liu Y Y, Zeng Y J, Jia P Z, Cao X H, Jiang X W, Chen K Q 2018 J. Phys.: Condens. Matter 30 275701

    [7]

    Hu S Q, Zhang Z W, Wang Z T, Zeng K Y, Cheng Y, Chen J, Zhang G 2018 ES Energy Environ 1 74

    [8]

    Liu C Q, Chen M, Yu W, He Y 2018 ES Energy Environ 2 31

    [9]

    Wang H, Hu S, Takahashi K, Zhang X, Takamatsu H, Chen J 2017 Nat. Commun. 8 15843

    [10]

    Chen X K, Xie Z X, Zhou W X, Tang L M, Chen K Q 2016 Carbon 100 492

    [11]

    Yang N, Xu X, Zhang G, Li B W 2012 AIP Adv. 2 041410

    [12]

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

    [13]

    Ziman J M 1962 Electrons and Phonons: The Theory of Transport Phenomena in Solids (Oxford: Clarendon) p168

    [14]

    Hochbaum A I, Chen R, Delgado R D, Liang W J, Garnett E C, Najarian M, Majumdar A, Yang P D 2008 Nature 451 163

    [15]

    Maldovan M 2015 Nat. Mater. 14 667

    [16]

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

    [17]

    Yu J K, Mitrovic S, Tham D, Varghese J, Heath J R 2010 Nat. Nanotech. 5 718

    [18]

    Hopkins P E, Reinke C M, Su M F, Olsson III R H, Shaner E A, Leseman Z C, Serrano J R, Phinney L M, El-Kady I 2010 Nano Lett. 11 107

    [19]

    Alaie S, Goettler D F, Su M, Leseman Z C, Reinke C M, El-Kady I 2015 Nat. Commun. 6 7228

    [20]

    Lee J, Lee W, Wehmeyer G, Dhuey S, Olynick D L, Cabrini S, Dames C, Urban J J, Yang P D 2017 Nat. Commun. 8 14054

    [21]

    Maire J, Anufriev R, Yanagisawa R, Ramiere A, Volz S, Nomura M 2017 Sci. Adv. 3 e1700027

    [22]

    Wagner M R, Graczykowski B, Reparaz J S, El Sachat A, Sledzinska M, Alzina F, Sotomayor Torres C M 2016 Nano Lett. 16 5661

    [23]

    Dechaumphai E, Chen R 2012 J. Appl. Phys. 111 073508

    [24]

    Ravichandran N K, Minnich A J 2014 Phys. Rev. B 89 205432

    [25]

    Tesanovic Z, Jaric M V, Maekawa S 1986 Phys. Rev. Lett. 57 2760

    [26]

    Liu X J, Zhou Z F, Yang L W, Li J W, Xie G F, Fu S Y, Sun C Q 2011 J. Appl. Phys. 109 074319

    [27]

    Klemens P G 1955 Proc. Phys. Soc. A 68 1113

    [28]

    Pauling L 1947 J. Am. Chem. Soc. 69 542

    [29]

    Bahn S R, Jacobsen K W 2001 Phys. Rev. Lett. 87 266101

    [30]

    Sun C Q 2007 Prog. Solid State Chem. 35 1

    [31]

    Pan L K, Sun C Q, Li C M 2004 J. Phys. Chem. B 108 3404

    [32]

    McGaughey A J H, Landry E S, Sellan D P, Amon C H 2011 Appl. Phys. Lett. 99 131904

    [33]

    Xie G F, Guo Y, Wei X L, Zhang K W, Sun L Z, Zhong J X, Zhang G, Zhang Y W 2014 Appl. Phys. Lett. 104 233901

    [34]

    Xie G F, Guo Y, Li B H, Yang L W, Zhang K W, Tang M H, Zhang G 2013 Phys. Chem. Chem. Phys. 15 14647

    [35]

    Song I H, Peter Y A, Meunier M 2007 J. Micromech. Microeng. 17 1593

    [36]

    Li D Y, Wu Y Y, Kim P, Shi L, Yang P D, Majumdar A 2003 Appl. Phys. Lett. 83 2934

    [37]

    Ju Y S, Goodson K E 1999 Appl. Phys. Lett. 74 3005

    [38]

    Liu W, Asheghi M 2006 J. Heat Transfer 128 75

    [39]

    Cuffe J, Eliason J K, Maznev A A, Collins K C, Johnson J A, Shchepetov A, Prunnila M, Ahopelto J, Torres C M S, Chen G, Nelson K A 2015 Phys. Rev. B 91 245423

  • 图 1  薄膜和纳米线示意图

    Fig. 1.  Schematic illustration of thin film and nanowire.

    图 2  硅纳米结构热导率计算值与实验值的比较 (a)纳米线; (b)薄膜

    Fig. 2.  Comparison of thermal conductivities between models and experimental data for silicon nanostructures: (a) Nanowire; (b) thin film.

    图 3  室温下声子各种散射率与角频率的关系

    Fig. 3.  Various scattering rates of phonons as functions of phonon angular frequency at room temperature.

    图 4  R与温度的关系

    Fig. 4.  R as a function of temperature.

  • [1]

    Wang J, Xie F, Cao X H, An S C, Zhou W X, Tang L M, Chen K Q 2017 Sci. Rep. 7 41418

    [2]

    Ding Z D, An M, Mo S Q, Yu X X, Jin Z L, Liao Y X, Esfarjani K, Lü J T, Shiomi J, Yang N 2019 J. Mater. Chem. A 7 2114

    [3]

    Tang L P, Tang L M, Geng H, Yi Y P, Wei Z M, Chen K Q, Deng H X 2018 Appl. Phys. Lett. 112 012101

    [4]

    Ouyang T, Jiang E L, Tang C, Li J, He C Y, Zhong J X 2018 J. Mater. Chem. A 6 21532

    [5]

    Xie G F, Ju Z F, Zhou K K, Wei X L, Guo Z X, Cai Y Q, Zhang G 2018 npj Comput. Mater. 4 21

    [6]

    Liu Y Y, Zeng Y J, Jia P Z, Cao X H, Jiang X W, Chen K Q 2018 J. Phys.: Condens. Matter 30 275701

    [7]

    Hu S Q, Zhang Z W, Wang Z T, Zeng K Y, Cheng Y, Chen J, Zhang G 2018 ES Energy Environ 1 74

    [8]

    Liu C Q, Chen M, Yu W, He Y 2018 ES Energy Environ 2 31

    [9]

    Wang H, Hu S, Takahashi K, Zhang X, Takamatsu H, Chen J 2017 Nat. Commun. 8 15843

    [10]

    Chen X K, Xie Z X, Zhou W X, Tang L M, Chen K Q 2016 Carbon 100 492

    [11]

    Yang N, Xu X, Zhang G, Li B W 2012 AIP Adv. 2 041410

    [12]

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

    [13]

    Ziman J M 1962 Electrons and Phonons: The Theory of Transport Phenomena in Solids (Oxford: Clarendon) p168

    [14]

    Hochbaum A I, Chen R, Delgado R D, Liang W J, Garnett E C, Najarian M, Majumdar A, Yang P D 2008 Nature 451 163

    [15]

    Maldovan M 2015 Nat. Mater. 14 667

    [16]

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

    [17]

    Yu J K, Mitrovic S, Tham D, Varghese J, Heath J R 2010 Nat. Nanotech. 5 718

    [18]

    Hopkins P E, Reinke C M, Su M F, Olsson III R H, Shaner E A, Leseman Z C, Serrano J R, Phinney L M, El-Kady I 2010 Nano Lett. 11 107

    [19]

    Alaie S, Goettler D F, Su M, Leseman Z C, Reinke C M, El-Kady I 2015 Nat. Commun. 6 7228

    [20]

    Lee J, Lee W, Wehmeyer G, Dhuey S, Olynick D L, Cabrini S, Dames C, Urban J J, Yang P D 2017 Nat. Commun. 8 14054

    [21]

    Maire J, Anufriev R, Yanagisawa R, Ramiere A, Volz S, Nomura M 2017 Sci. Adv. 3 e1700027

    [22]

    Wagner M R, Graczykowski B, Reparaz J S, El Sachat A, Sledzinska M, Alzina F, Sotomayor Torres C M 2016 Nano Lett. 16 5661

    [23]

    Dechaumphai E, Chen R 2012 J. Appl. Phys. 111 073508

    [24]

    Ravichandran N K, Minnich A J 2014 Phys. Rev. B 89 205432

    [25]

    Tesanovic Z, Jaric M V, Maekawa S 1986 Phys. Rev. Lett. 57 2760

    [26]

    Liu X J, Zhou Z F, Yang L W, Li J W, Xie G F, Fu S Y, Sun C Q 2011 J. Appl. Phys. 109 074319

    [27]

    Klemens P G 1955 Proc. Phys. Soc. A 68 1113

    [28]

    Pauling L 1947 J. Am. Chem. Soc. 69 542

    [29]

    Bahn S R, Jacobsen K W 2001 Phys. Rev. Lett. 87 266101

    [30]

    Sun C Q 2007 Prog. Solid State Chem. 35 1

    [31]

    Pan L K, Sun C Q, Li C M 2004 J. Phys. Chem. B 108 3404

    [32]

    McGaughey A J H, Landry E S, Sellan D P, Amon C H 2011 Appl. Phys. Lett. 99 131904

    [33]

    Xie G F, Guo Y, Wei X L, Zhang K W, Sun L Z, Zhong J X, Zhang G, Zhang Y W 2014 Appl. Phys. Lett. 104 233901

    [34]

    Xie G F, Guo Y, Li B H, Yang L W, Zhang K W, Tang M H, Zhang G 2013 Phys. Chem. Chem. Phys. 15 14647

    [35]

    Song I H, Peter Y A, Meunier M 2007 J. Micromech. Microeng. 17 1593

    [36]

    Li D Y, Wu Y Y, Kim P, Shi L, Yang P D, Majumdar A 2003 Appl. Phys. Lett. 83 2934

    [37]

    Ju Y S, Goodson K E 1999 Appl. Phys. Lett. 74 3005

    [38]

    Liu W, Asheghi M 2006 J. Heat Transfer 128 75

    [39]

    Cuffe J, Eliason J K, Maznev A A, Collins K C, Johnson J A, Shchepetov A, Prunnila M, Ahopelto J, Torres C M S, Chen G, Nelson K A 2015 Phys. Rev. B 91 245423

  • [1] 华钰超, 曹炳阳. 多约束纳米结构的声子热导率模型研究. 物理学报, 2015, 64(14): 146501. doi: 10.7498/aps.64.146501
    [2] 李静, 冯妍卉, 张欣欣, 黄丛亮, 杨穆. 考虑界面散射的金属纳米线热导率修正. 物理学报, 2013, 62(18): 186501. doi: 10.7498/aps.62.186501
    [3] 黄丛亮, 冯妍卉, 张欣欣, 李静, 王戈, 侴爱辉. 金属纳米颗粒的热导率. 物理学报, 2013, 62(2): 026501. doi: 10.7498/aps.62.026501
    [4] 侯泉文, 曹炳阳, 过增元. 碳纳米管的热导率:从弹道到扩散输运. 物理学报, 2009, 58(11): 7809-7814. doi: 10.7498/aps.58.7809
    [5] 张程宾, 程启坤, 陈永平. 分形结构纳米复合材料热导率的分子动力学模拟研究. 物理学报, 2014, 63(23): 236601. doi: 10.7498/aps.63.236601
    [6] 程笃庆, 关庆丰, 朱健, 邱东华, 程秀围, 王雪涛. 强流脉冲电子束诱发纯镍表层纳米结构的形成机制. 物理学报, 2009, 58(10): 7300-7306. doi: 10.7498/aps.58.7300
    [7] 王丹, 贺永宁, 叶鸣, 崔万照. 金纳米结构表面二次电子发射特性. 物理学报, 2018, 67(8): 087902. doi: 10.7498/aps.67.20180079
    [8] 曹炳阳, 董若宇, 孔杰, 陈恒, 徐雁, 容启亮, 蔡岸. 纳米孔模板浸润法制备聚乙烯纳米线阵列的热导率的实验研究. 物理学报, 2012, 61(4): 046501. doi: 10.7498/aps.61.046501
    [9] 王建立, 顾明, 张兴, 熊国平, 梁吉. 多壁碳纳米管/聚丙烯复合材料热导率研究. 物理学报, 2009, 58(7): 4536-4541. doi: 10.7498/aps.58.4536
    [10] 杨平, 吴勇胜, 许海锋, 许鲜欣, 张立强, 李培. TiO2/ZnO纳米薄膜界面热导率的分子动力学模拟. 物理学报, 2011, 60(6): 066601. doi: 10.7498/aps.60.066601
    [11] 杨平, 王晓亮, 李培, 王欢, 张立强, 谢方伟. 氮掺杂和空位对石墨烯纳米带热导率影响的分子动力学模拟. 物理学报, 2012, 61(7): 076501. doi: 10.7498/aps.61.076501
    [12] 李威, 冯妍卉, 陈阳, 张欣欣. 碳纳米管中点缺陷对热导率影响的正交试验模拟分析. 物理学报, 2012, 61(13): 136102. doi: 10.7498/aps.61.136102
    [13] 郑伯昱, 董慧龙, 陈非凡. 基于量子修正的石墨烯纳米带热导率分子动力学表征方法. 物理学报, 2014, 63(7): 076501. doi: 10.7498/aps.63.076501
    [14] 贺慧芳, 陈志权. 用正电子湮没研究纳米碲化铋的缺陷及其对热导率的影响. 物理学报, 2015, 64(20): 207804. doi: 10.7498/aps.64.207804
    [15] 许路加, 胡明, 杨海波, 杨孟琳, 张洁. 基于微结构参数建模的多孔硅绝热层热导率研究. 物理学报, 2010, 59(12): 8794-8800. doi: 10.7498/aps.59.8794
    [16] 刘天庆, 孙玮, 李香琴, 孙相彧, 艾宏儒. 纳米结构表面上部分润湿液滴合并诱导弹跳的理论研究. 物理学报, 2014, 63(8): 086801. doi: 10.7498/aps.63.086801
    [17] 胡梦丹, 张庆宇, 孙东科, 朱鸣芳. 纳米结构超疏水表面冷凝现象的三维格子玻尔兹曼方法模拟. 物理学报, 2019, 68(3): 030501. doi: 10.7498/aps.68.20181665
    [18] 黄金华, 张 琨, 潘 楠, 高志伟, 王晓平. 表面修饰ZnO纳米线紫外光响应的增强效应. 物理学报, 2008, 57(12): 7855-7859. doi: 10.7498/aps.57.7855
    [19] 李威, 冯妍卉, 唐晶晶, 张欣欣. 碳纳米管Y形分子结的热导率与热整流现象. 物理学报, 2013, 62(7): 076107. doi: 10.7498/aps.62.076107
    [20] 李玲红, 李培刚, 符秀丽, 陈雷明, 唐为华, 张海英. PED沉积La-Sr-Cu-O薄膜表面的有序纳米结构. 物理学报, 2005, 54(6): 2843-2846. doi: 10.7498/aps.54.2843
  • 引用本文:
    Citation:
计量
  • 文章访问数:  1086
  • PDF下载量:  7
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-02-12
  • 修回日期:  2019-03-06
  • 上网日期:  2019-04-01
  • 刊出日期:  2019-04-20

表面低配位原子对声子的散射机制

  • 1. 湘潭大学物理与光电工程学院, 湘潭 411105
  • 2. 湖南科技大学, 新能源储存与转换先进材料湖南省重点实验室, 湘潭 411201
  • 通信作者: 谢国锋, xieguofeng@hnust.cn
    基金项目: 国家自然科学基金(批准号: 11874145)资助的课题.

摘要: 由于纳米结构具有极高的表体比, 声子-表面散射机制对声子的热输运性质起到关键作用. 提出了表面低配位原子对声子的散射机制, 并且结合量子微扰理论与键序理论推导出该机制的散射率. 由于散射率正比于材料的表体比, 这种散射机制对声子输运的重要性随着纳米结构尺寸的减小而增大. 散射率正比于声子频率的4次方, 所以这种散射机制对高频声子的作用远远强于对低频声子的作用. 基于声子玻尔兹曼输运方程, 计算了硅纳米薄膜和硅纳米线的热导率, 发现本文模型比传统的声子-边界散射模型更接近实验值. 此发现不仅有助于理解声子-表面散射的物理机制, 也有助于应用声子表面工程调控纳米结构的热输运性质.

English Abstract

    • 对纳米结构热输运的理解和调控对于热电器件的开发[1-6]以及纳米电子器件的热管理[7-10]起到非常重要的作用. 参与导热的声子主要集中在太赫兹频段, 由于纳米结构的特征尺寸可以小于太赫兹频段声子倒逆散射的自由程, 声子-表面散射机制占支配地位, 从而对纳米结构热输运产生极为重要的影响, 所以当材料尺寸进入纳米尺度, 热输运呈现出诸多不同于块体材料的新特征, 例如反常热导率、量子限制效应和傅里叶定律失效等[11]. 声子-表面散射机制影响和决定纳米结构热输运的诸多特征, 因此利用表面效应成为声子工程调控纳米结构热导率的有效手段. 然而目前人们对于声子-表面散射的物理机制不是十分清楚. 声子与表面作用有两种物理图像, 即基于粒子属性的非相干散射和基于波动属性的相干散射[12]. 传统的声子-表面非相干散射模型将声子被表面的散射视为粒子与容器边界的碰撞, 散射率公式为$ \tau _{\rm{B}}^{ - 1} = \dfrac{{1 - P}}{{1 + P}}\dfrac{v}{\varLambda }$, $ v$为声子速度, $ \varLambda $为纳米结构的特征尺寸, P为镜面因子. 早期的文献中将镜面因子作为一个常数来处理, 后来人们用Ziman公式加以完善[13], Ziman公式表明镜面因子随声子波长$ \lambda $的减小或者表面粗糙度$ \delta $的增加而急剧降低. 虽然传统的声子-表面散射模型较好地解释了一些纳米结构热导率的实验现象和规律, 但是基于这种声子-表面散射模型的热导率计算结果往往比实验数据高, 而且纳米结构的尺寸越小, 计算值与实验值相对差别越大. 特别是当纳米结构的表面非常粗糙时, 例如电化学腐蚀方法制备的粗糙表面纳米线, 实验测量表明其热导率远远低于相同直径的平整表面纳米线[14]. 这种由于表面粗糙度增加而使纳米结构热导率远远低于Casimir极限(即$ P = 0$)的现象, 无法用传统的声子-表面非相干散射模型来解释. 根据传统的声子-表面非相干散射模型, 纳米结构的热导率只能无限接近但不能低于Casimir极限. 声子不但具有粒子性, 还具有波动性, 声子被周期性分布的表面或界面反射可以产生干涉效应, 声子的相干性导致声子色散关系改变, 从而使声子群速度和态密度都发生变化, 并且相干效应还能产生声子带隙, 因此可以利用周期边界对声子的相干作用来调控纳米结构的热输运性质, 并取得经典调控方法无法达到的效果[15, 16]. 对于多孔硅薄膜声子晶体, 在$ \lambda < p < {\varLambda _{\rm{U}}}$的尺寸范围(对于室温下工作的硅材料而言, 这正是目前纳米制备技术到达的尺寸范围, 其中$ p$为声子晶体的周期间距尺寸, $ \lambda $为占支配地位声子的波长, $ {\varLambda _{\rm{U}}}$为占支配地位声子的倒逆散射自由程), 近年来科学家们就声子相干输运是否占据支配地位这一问题展开了激烈的争论. Yu等[17]和Hopkins等[18]认为即使在室温下声子-表面相干散射也是导致多孔硅薄膜声子晶体热导率极度降低的原因, 但是基于完全相干模型得到的热导率计算值明显低于实验值[19]. 最近Lee等[20]通过实验发现, 当温度超过14 K, 多孔硅薄膜中的孔是否周期分布对热导率没有影响, 因为非周期分布的孔破坏了声子的干涉效应, 所以他们认为多孔硅薄膜声子晶体中的声子-表面相干散射机制在温度14 K以上时是不重要的. Maire等[21]和Wagner等[22]给出了类似的结论, 即在$ \lambda < p < {\varLambda _{\rm{U}}}$的尺寸范围, 声子-表面非相干散射仍然占据主导地位. 然而在理论研究方面, 玻尔兹曼输运方程方法[23]和蒙特卡罗模拟[24]都表明应用传统的声子-表面非相干散射模型计算多孔硅薄膜声子晶体的热导率, 计算值明显高于实验值. 以上这些争议和矛盾都表明人们对于声子-表面散射的物理机制还不是很清楚, 需要对其做进一步研究.

    • 材料的表面并不等同于数学意义上的边界, 材料表面有非常丰富的物理和化学效应. 例如表面原子配位数降低导致表面原子键长收缩以及键能增强, 使得表面势阱变深, 从而导致电子在表面局域化, 表面散射降低了电导率[25]. 我们认为现有声子-表面散射机制研究的不足之处在于只考虑了表面的数学属性, 即表面的几何限制对粒子或者波的运动产生突变, 忽视了表面附近微观物理属性不同于内部, 以及由此产生对声子的作用. 表面原子键长收缩键能增强[26], 对系统的哈密顿量产生微扰, 同理应该对声子输运性质产生重要影响, 因此我们提出表面低配位原子对声子的散射机制. 下面结合量子微扰理论与键序理论推导这种声子-表面散射机制的散射率.

      完整晶格的哈密顿量为

      $ {{{H}}_0} = \frac{1}{2}M\mathop \sum \limits_k \omega _k^2\left( {{{a}}_k^ \dagger {{{a}}_k} + {{{a}}_k}{{a}}_k^ \dagger } \right), $

      式中M为原子质量; $ {\omega _k}$为简振模式$ k$的角频率; $ {{{a}}^ \dagger }$, $ {{a}}$分别为声子产生和湮灭算符, 在量子力学的矩阵表达中非零元素如下:

      $ {\rm{}}{\left( {{{a}}_k^ \dagger } \right)_{N,N - 1}} = {\left( {\frac{\hbar }{{M{\omega _k}}}N} \right)^{1/2}}, $

      $ {\left( {{{{a}}_k}} \right)_{N,N + 1}} = {\left[ {\frac{\hbar }{{M{\omega _k}}}\left( {N + 1} \right)} \right]^{1/2}}. $

      N为矩阵的行号. 如果哈密顿量的微扰是$ {{{u}}_x}$或者$ {{\dot {{u}}}_x}$的二次函数($ {{{u}}_x}$指第x个原子偏离平衡位置的位移, $ {{\dot {{u}}}_x}$则是相对应的速度), 那么哈密顿量的微扰算符可以表达为

      $ {{H}}' = \mathop \sum \limits_{k,k'} {c_{k,k'}}{{a}}_k^ \dagger {{{a}}_{k'}},\;\;\;{c_{k,k'}} = c_{k',k}^*. $

      上面求和表达式中每一项意味着声子模式$ k'$的湮灭以及声子模式$ k$的产生, 即声子从$ k'$$ k$的散射. Klemens根据量子力学微扰理论, 推导出声子散射率表达式为[27]

      $ \begin{split} {\tau ^{ - 1}} = & - \frac{1}{n}\cdot\frac{{{\rm{d}}N}}{{{\rm{d}}t}} = \frac{{3V}}{{{\text{π}}{v^3}{M^2}}}\int {c_{k,k'}^2\left( {1 - \frac{{n'}}{n}} \right)} \\ & \times {\text{δ}} \left( {\omega ' - \omega } \right){\rm{d}}{\omega ^{'}} ,\end{split} $

      其中$ V$为晶体体积, $ \omega $为角频率, $ \dfrac{{n'}}{n} $等于散射角余弦的平均值, 对于球对称散射其值为0. 只要找到相应的微扰源, 并且求出(4)式中$ {c_{k,k'}}$的表达式, 就可以通过(5)式求出该物理机制的声子散射率. 如果材料中有一个原子的力常数不同于其他原子, 根据Klemens理论[27]:

      $ c_{k,k'}^2 = \frac{2}{3}\cdot\frac{{{M^2}}}{{{G^2}}}{\omega ^2}{{\omega '}^2}{\left( {\frac{{{\text{δ}}v}}{v}} \right)^2}, $

      其中$ G$为晶体中原子的总数. 因为$ v \propto \sqrt f $($ f$是原子间的力常数), 所以$ {\left( {\dfrac{{{\text{δ}}v}}{v}} \right)^2} = \dfrac{1}{4}{\left( {\dfrac{{{\text{δ}}f}}{f}} \right)^2}$, $ {\text{δ}}f$是力常数的变化量. 结合(5)式和(6)式, 可以得到单个力常数缺陷对声子的散射率为

      $ \tau {\left( \omega \right)^{ - 1}} = \frac{V}{{2{\text{π}}{G^2}{v^3}}}{\left( {\frac{{{\text{δ}}f}}{f}} \right)^2}{\omega ^4}. $

      如果晶体内部随机分布着多个力常数缺陷, 它们对角频率为$ \omega $的声子的总散射率正比于缺陷的数目, 如(8)式所示:

      $ \tau {\left( \omega \right)^{ - 1}} = \frac{{V{N_{\rm{I}}}}}{{2{\text{π}}{G^2}{v^3}}}{\left( {\frac{{{\text{δ}}f}}{f}} \right)^2}{\omega ^4} = x\frac{V}{{2{\text{π}}G{v^3}}}{\left( {\frac{{{\text{δ}}f}}{f}} \right)^2}{\omega ^4}, $

      其中$ {N_{\rm{I}}}$为力常数缺陷的数目; $ x$为力常数缺陷的浓度.

      Pauling[28]指出原子的配位数降低导致原子半径收缩, 低配位原子的化学键变短变强[29], 即低配位原子之间的力常数增加. Sun[30]提出了键弛豫理论, 给出了原子键收缩系数以及单键能计算式, 即

      $ \left\{ \begin{aligned} &{C_z} = {d_z}/{d_0} = 2/\left\{ {1 + {\rm{exp}}\left[ {\left( {12 - z} \right)/\left( {8z} \right)} \right]} \right\},\\ &{E_z} = C_z^{ - m}{E_{\rm{b}}}, \end{aligned} \right. $

      式中$ z$为原子的配位数; $ d_z$是配位数为$ z$的原子键长; $ E_z$是配位数为$ z$的原子的单键能; $ C_z$为键收缩系数; $ {E_{\rm{b}}}$$ {d_0}$分别为体内原子的单键能和键长; m是一个与材料性质有关的因子, 对于硅, m = 4.88. 键弛豫理论将材料表面最外层的两层原子视为表面层, 表面层原子的配位数小于体内原子, 由量纲分析可得, 原子间力常数$ {f_z} = {\left. {{{\rm{d}}^2}u\left( r \right)/{\rm{d}}{r^2}} \right|_{r = {d_z}}}$$ \propto {E_z}/d_z^2 $, $ {u\left( r \right)}$为原子间作用势, $ r$为原子间距离. 根据(9)式可知表面层低配位原子间力常数大于体内原子间力常数, 从而对晶格振动系统的哈密顿量产生微扰, 所以表面低配位原子对声子产生散射作用. 下面根据(8)式和键弛豫理论来推导这种机制的声子散射率.

      键弛豫理论将材料表面最外层的两层原子视为表面层, 表面第一层原子配位数为4, 第二层原子配位数为6, 我们做一个近似处理, 将表面层原子的平均配位数设为(4 + 6)/2 = 5, 表面层原子数占总原子数的比例为$ x = 2{C_5}{d_0} \cdot $SVR, SVR是材料的表体比, 即表面积除以体积. 根据(8)式, 表面低配位原子对声子的散射率为

      $ \tau _{{\rm{Bond}}}^{ - 1} = \frac{1}{2}x{\left( {\frac{{{f_5} - {f_{\rm{b}}}}}{{{f_{\rm{b}}}}}} \right)^2}{\left( {\frac{5}{6}} \right)^2}\frac{V}{{2{\text{π}}G{v^3}}}{\omega ^4}, $

      其中$ {{f_{\rm{b}}}}$为体内原子间的力常数, $ {{f_5}}$为表层配位数为5的原子间力常数. 式中的因子1/2, 是因为表面层内的每个原子键连接着两个低配位原子, 所以只有一半的低配位原子对系统的哈密顿量产生微扰. 另一个因子$ {\left( {\dfrac{5}{6}} \right)^2}$则是因为Klemens在推导(7)式过程中假设每个原子有6个最近邻原子, 而在本文模型中表层原子配位数为5, 即等效为每个原子有5个最近邻原子, 原子力常数改变引起的哈密顿微扰量矩阵元的平方应该乘以因子$ {\left( {\dfrac{5}{6}} \right)^2}$. 由量纲分析可得, 原子间力常数${k_z} = $$ {\left. {{{\rm{d}}^2}u\left( r \right)/{\rm{d}}{r^2}} \right|_{r = {d_z}}} \propto {E_z}/d_z^2$, 基于(9)式, (10)式可写为

      $ \tau _{{\rm{Bond}}}^{ - 1} = {C_5}{d_0} \cdot {\rm{SVR}} \cdot {\left( {C_5^{ - m - 2} - 1} \right)^2}{\left( {\frac{5}{6}} \right)^2}\frac{V}{{2{\text{π}}G{v^3}}}{\omega ^4}. $

      (11)式给出了材料表面低配位原子对声子的散射率, 它表明该机制的声子散射率正比于材料的表体比. 对于宏观材料, 由于表体比极低, 这种声子散射机制的作用非常微小, 但是对于纳米材料, 表体比非常高, 而且随着尺寸的减小, 表体比越来越大, 该声子散射机制的作用越来越重要. 另外(11)式还表明, 表面低配位原子对声子的散射率正比于声子频率的4次方, 是强烈频率相关的, 对高频声子的散射作用远远强于对低频声子的散射作用. 对于硅, $ m$ = 4.88, $ {d_0} = 0.263$ nm[31], $ v = $$ 6733$ m/s[32], $ V/G = {a^3}/8$, 晶格常数$ a = 0.543$ nm, 所以硅纳米结构表面低配位原子对声子的散射率为

      $ \tau _{{\rm{Bond}}}^{ - 1} = 1.33 \times {10^{ - 51}} \cdot {\rm{SVR}} \cdot {\omega ^4}, $

      其中表体比SVR的单位为m–1.

    • 下面将这种表面低配位原子对声子的散射机制引入到声子玻尔兹曼输运方程, 来研究典型硅纳米结构, 如硅纳米薄膜平面内方向的热导率和硅纳米线的热导率. 图1为薄膜和纳米线的示意图. 基于近弛豫时间近似以及德拜模型, 根据声子玻尔兹曼输运方程, 纳米结构在i方向的热导率可以由积分公式(13)式来计算[33]:

      图  1  薄膜和纳米线示意图

      Figure 1.  Schematic illustration of thin film and nanowire.

      $ \begin{split} \kappa =\; & \dfrac{{3{k_{\rm{B}}}}}{{8{\rm{}}{{\rm{}}^3}}}\int_0^{{\omega _{\rm{D}}}} {\int_0^{{\rm{}}} {\int_0^{2{\rm{}}} {\dfrac{{{{\left( {\hbar \omega /\left( {{k_{\rm{B}}}T} \right)} \right)}^2}{{\rm{e}}^{\hbar \omega /\left( {{k_{\rm{B}}}T} \right)}}}}{{{{\left( {{{\rm{e}}^{\hbar \omega /\left( {{k_{\rm{B}}}T} \right)}} - 1} \right)}^2}}}} } } \\ & \times \dfrac{{v_i^2}}{{{v^3}}}{\omega ^2}\tau {\rm{sin}}\theta {\rm{d}}\omega {\rm{d}}\theta {\rm{d}}\psi , \end{split} $

      式中$ {{\omega _{\rm{D}}}}$为德拜频率; $ {{k_{\rm{B}}}}$是玻尔兹曼常数; $ v_i$为群速度vi方向的分量, 对于图1所示薄膜平面内方向$ {v_i} = {v_X} = v{\rm{sin}}\theta {\rm{cos}}\psi $, 而对于图1所示纳米线$ {v_i} = {v_Z} = v{\rm{cos}}\theta $, $ \theta $为声子波失与Z轴的夹角, $ \psi $为声子波失在XY平面的投影矢量与X轴的夹角, $ \tau $为声子平均弛豫时间. 根据Matthiessen’s规则, $ {\tau ^{ - 1}} = \tau _{\rm{U}}^{ - 1} + \tau _{\rm{B}}^{ - 1} + \tau _{{\rm{Bond}}}^{ - 1}$, $ \tau _{\rm{U}}^{ - 1}$, $ \tau _{\rm{B}}^{ - 1}$$ \tau _{{\rm{Bond}}}^{ - 1}$分别为声子-声子U过程散射率, 传统的声子-边界散射率, 以及(12)式所示的表面低配位原子对声子的散射率.

      声子-声子U过程散射率为

      $ \tau _{\rm{U}}^{ - 1} = A\left( T \right){\omega ^2}, $

      式中$ A\left( T \right)$是一个与温度有关的参数, 是由体材料硅的热导率来确定, 假设温度为T时体硅的热导率为$ {\kappa _{{\rm{bulk}}}}$, $ A\left( T \right)$的计算式为[34]

      $ A\left( T \right) = \frac{{{k_{\rm{B}}}}}{{2{{\text{π}}^2}v{\kappa _{{\rm{bulk}}}}}}\int_0^{{\omega _{\rm{D}}}} {\frac{{{{\left( {\hbar \omega /\left( {{k_B}T} \right)} \right)}^2}{e^{\hbar \omega /\left( {{k_{\rm{B}}}T} \right)}}}}{{{{\left( {{{\rm{e}}^{\hbar \omega /\left( {{k_{\rm{B}}}T} \right)}} - 1} \right)}^2}}}{\rm{d}}\omega } . $

      传统的声子-边界碰撞自由程应用McGaughey等[32]的公式:

      $ \varLambda = \left\{ {\begin{align} & {\dfrac{L}{{2\left| {\cos \theta } \right|}},\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\text{薄膜}},}\\ & {\dfrac{{4D}}{{3{\text{π}}\left| {\sin \theta } \right|}},\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\text{纳米线}},} \end{align}} \right. $

      式中$ L$为薄膜的厚度, $ D$为纳米线的直径. 传统的声子-边界散射率为

      $ \tau _{\rm{B}}^{ - 1} = \frac{{1 - P}}{{1 + P}}\frac{v}{\varLambda }, $

      鉴于实验制备的纳米线或薄膜表面粗糙度通常不低于1 nm[35], 所以在我们的计算中镜面因子P = 0, 即完全漫散射.

      由(12)式可知, 纳米结构表面低配位原子对声子的散射率正比于结构的表体比. 假设纳米线的长度远远大于直径, 薄膜的长和宽都远远大于厚度, 那么纳米线和薄膜的表体比分别为

      $ {\rm{SVR}} = \left\{ {\begin{align} & {\dfrac{2}{L},\;\;\;\;\;\;\;\;\;{\text{薄膜}},}\\ & {\dfrac{4}{D},\;\;\;\;\;\;\;\;{\text{纳米线}}.} \end{align}} \right. $

      应用上面的方法计算了室温下(T = 300 K)硅纳米线和硅薄膜的热导率, 分别如图2(a), (b)所示. 图中黑色曲线是考虑表面低配位原子对声子散射的计算结果, 而红色曲线是不考虑这种表面散射机制的传统声子-边界模型的计算结果. 其中图2(a)的实验数据来自文献[36], 图2(b)的实验数据来自文献[3739]. 本文模型的计算结果与实验测量吻合较好, 而传统的传统声子-边界模型的计算结果则高估了纳米线和薄膜的热导率, 这种高估是由于传统模型没有考虑表面低配位原子对声子的散射机制. 而且从图2可以看出, 当纳米线的直径(或者薄膜的厚度)越小, 传统声子-边界模型的计算结果与实验值的相对偏差越大, 这是因为纳米结构的特征尺寸越小, 其表体比越大, 表面低配位原子对声子的散射率越高, 如果在计算中不考虑这种散射机制, 必然导致计算值的偏差越大.

      图  2  硅纳米结构热导率计算值与实验值的比较 (a)纳米线; (b)薄膜

      Figure 2.  Comparison of thermal conductivities between models and experimental data for silicon nanostructures: (a) Nanowire; (b) thin film.

      为了进一步说明表面低配位原子对声子的散射机制对于纳米结构声子输运的重要性, 在图3中给出了室温下直径为100 nm的硅纳米线中声子-声子U过程散射率$ \tau _{\rm{U}}^{ - 1}$、传统的声子-边界散射率$ \tau _{\rm{B}}^{ - 1}$, 以及表面低配位原子对声子的散射率$ \tau _{{\rm{Bond}}}^{ - 1}$随声子角频率的关系. 将$ \tau _{\rm{B}}^{ - 1}$$ \tau _{{\rm{Bond}}}^{ - 1}$两条曲线的交点横坐标定义为特征频率$ {\omega _{\rm{C}}}$, 从图3可知$ {\omega _{\rm{C}}}$ = 33.7 Trad/s, 对于角频率小于$ {\omega _{\rm{C}}}$的声子, 在上述三种散射机制中传统的声子-边界散射占据支配地位, 而对于角频率大于$ {\omega _{\rm{C}}}$的声子, 表面低配位原子对声子的散射机制最为重要. 根据德拜线性色散假设以及声子服从Bose-Einstein分布, 角频率大于$ {\omega _{\rm{C}}}$的声子占所有声子的比例为

      图  3  室温下声子各种散射率与角频率的关系

      Figure 3.  Various scattering rates of phonons as functions of phonon angular frequency at room temperature.

      $ R\left( T \right) = \frac{{\displaystyle\int_{{\omega _{\rm C}}}^{{\omega _{\rm D}}} {{\left( {{{\rm{e}}^{\hbar \omega /\left( {{k_{\rm{B}}}T} \right)}} - 1} \right)}^{ - 1}}{\omega ^2}{\rm{d}}\omega }}{{\displaystyle\int_0^{{\omega _{\rm D}}} {{\left( {{{\rm{e}}^{\hbar \omega /\left( {{k_{\rm{B}}}T} \right)}} - 1} \right)}^{ - 1}}{\omega ^2}{\rm{d}}\omega }}, $

      对于$ {\omega _{\rm C}}$ = 33.7 Trad/s, T = 300 K, 从(19)式可以得到R = 71.3%, 这表明在室温下对于硅纳米线中的大部分声子, 表面低配位原子对它们的散射率高于其他两种机制, 而传统的声子输运模型却忽视了这种极为重要的散射机制.

      从(12)式和(17)式可知, $ {\rm{\tau}} _{{\rm{Bond}}}^{ - 1}$$ \tau _{\rm{B}}^{ - 1}$这两个散射率都与温度无关, 所以特征频率$ {\omega _{\rm{C}}}$也与温度无关, 只取决于纳米结构的几何形状. 在$ {\omega _{\rm{C}}}$给定的情况下, (19)式所示的R与温度有关. 图4给出了横截面为圆形的硅纳米线R与温度T的关系. 由于声子服从Bose-Einstein分布, 在低温的情况下, 大部分的声子布居在低频态. 从图3可知, 对于低频声子, 表面低配位原子对声子的散射率不但低于传统的声子-边界散射率, 也低于声子U散射率, 只有较少比例的高频声子, 表面低配位原子对其的散射率高于其他两种机制, 例如: 当T = 80 K时, R = 30.5%, 随着温度的进一步降低, R还会快速下降. 随着温度的升高, 更多的声子布居在高频态, 表面低配位原子对声子的散射机制对于声子输运的影响增大, 当T = 120 K时, 已经有50%的声子, $ \tau _{{\rm{Bond}}}^{ - 1}$大于其他两种散射率.

      图  4  R与温度的关系

      Figure 4.  R as a function of temperature.

    • 表面对声子的散射机制对于纳米结构中的声子输运行为起到关键作用. 纳米结构的表面有丰富而复杂的物理和化学效应, 对声子运动的影响不仅仅是简单的几何限制. 由于表面原子配位数降低, 原子键长收缩键能增强, 对晶格振动系统的哈密顿量产生微扰, 从而对声子产生散射. 基于键弛豫理论和量子微扰理论, 推导出了表面低配位原子对声子的散射率. 这种声子-表面散射机制的散射率正比于材料的表体比, 所以当纳米结构的特征尺寸越小表体比越大, 这种机制对声子输运的影响越强. 表面低配位原子对声子的散射率正比于声子频率的4次方, 所以这种散射机制对高频声子的作用特别强烈. 由于声子服从Bose-Einstein分布, 在极低的温度下, 高频声子所占的比例很小, 表面低配位原子散射机制对于声子输运的影响很小, 但是随着温度升高, 高频态声子的比例增加, 这种散射机制对声子输运的影响增强, 在室温下, 这种散射机制对硅纳米结构中的声子输运起到支配性作用. 应用声子玻尔兹曼输运方程, 计算了硅纳米薄膜和硅纳米线的热导率, 考虑了表面低配位原子对声子的散射机制, 热导率的计算值与实验值符合得很好. 纳米结构的热导率随着尺寸的减小而降低, 不仅是因为声子-边界散射更加频繁, 也因为表体比增加导致表面低配位原子对声子特别是对高频声子的散射更加强烈. 表面低配位原子对声子的散射不仅有助于理解声子-表面散射的物理机制, 也有助于应用声子表面工程调控纳米结构的热输运性质.

参考文献 (39)

目录

    /

    返回文章
    返回