2019年晶格热导率的推导.docVIP

晶格热导率的推导 　　记录下推导晶格热导率的过程，备忘！ 　　------范登栋-武汉大学物理科学与技术学院 　　1064540094@ 　　2015.4.28 　　　　晶格热导率kl的定义，单位时间dt通过单位面积dS单位温度梯度?T的热量。在下面的推导过程中，为简单化处理，仅考虑一个方向的热传导。 　　　　首先，我们定义单位时间dt通过单位面积dS的热量为热流J。这样晶格热导率可表达为 　　J??kl?T 将每一种声子模式?看做一个携带能量为??粒子，其速度为v?，则有： 　　J??N???v? ? 　　式中N?为单位体积内，声子模式为?的声子数。N?可表达为： 　　N??nf 这里f表示在一个能级(状态)上声子的占据数，很接近玻色-爱因斯坦分布fBE： 　　fBE? 　　1 　　? 　　ekBT?1 　　(在计算热导率的过程中，由于温度梯度的存在，分布函数f会偏离平衡分布fBE，这里将留到波尔兹曼输运理论来解决。)n为倒空间中体积元dq内所包含的声子状态数。单位体积倒空间中声子的状态数密度为： 　　?? 　　1 　　?2?? 　　3 　　因此， 　　n??dq? 　　dq 　　?2?? 　　3 　　最终，将 (2) 式写为积分形式，对整个布里渊区积分得到： 　　J? 　　?dq 　　?v?f 　　q?? 　　BZ 　　?2?? 　　3 　　结合式(1)和(7)可知，要求晶格热导率kl，关键在于求出声子的分布函数f。 　　(1) 　　(2) 　　(3) 　　(4) 　　(5) 　　(6) 　　(7) 　　　　接下来，我们将通过Boltzmann输运理论来解决这个问题。 　　Boltzmann输运理论是解决输运问题的重要工具，无论是对于声子还是电子或者其它粒子的输运问题都可借助Boltzmann理论来分析，得到定量的表达。 　　Boltzmann输运理论需要用到体系的分布函数f，对于电子就是Fermi-Dirac分布，而声子则是Bose-Einstein分布。对于一个由温度驱动的状态，当体系达到稳定的非平衡态，即所谓动态平衡时，有： 　　?f? 　　?f?0 (8) 　　?t 　　Temp drive 　　?t 　　diffusion 　　第一项表示由于温度的驱动所引起的平衡态的偏离，第二项表示声子的散射项，由于声子的散射使得非平衡态能够达到动态平衡。第一项我们可以进一步写为： 　　?f?T?x?t 　　? 　　?f?T?x?t?v?T?f 　　?T 　　Temp drive 　　则我们的分布函数f可写为： 　　?f?f 　　?t 　　??v?T 　　?T 　　diffusion 　　由于分布函数f只是对平衡分布函数fBE的微小偏移，引入声子弛豫时间??概念，表示系统有非平衡态回复到平衡态的平均时间。(10)式的左边可表达为： 　　?f?t 　　? 　　f?fBE 　　diffusion 　　? ? 　　结合(10)式和(11)式，可以得到： 　　f?f?f 　　BE?v??T 　　?T 　　从式(12)中可以看到，温度梯度会引起平衡分布函数f?f 　　BE的偏移，偏移量为?v?T?T 　　, 　　由于分布函数f只是对平衡分布函数fBE的微小偏移，将(12)式中的分布函数f用平衡分布函数fBE替代，得到： 　　f?ffBE 　　BE?v??T 　　??T 　　式代入(7)式，并考虑(13)式中的平衡分布函数fBE对热输运没有贡献，我们得到： 　　J? 　　v?(?v??f??T 　　?BEq?? 　　??BZ 　　?T)dq 　　?2??3 　　结合式(14)和式(1)，可得到热导率表达式为： 　　(9) 　　(10) (11) 　　(12) (13) (14) 　　kl? 　　最后将 　　q?BZ 　　? 　　??v?v?? 　　?fBEdq 　　(15) 3 　　?T?2?? 　　?fBE 　　的表达式代入，得到： ?T 　　kl?kB 　　将(16)式写为离散形式， 　　????dq 　　vv?f(f?1) (16) ?????BEBE3?kT?2???q?BZ?B 　　2 　　???k 　　kl?B????v?v???fBE(fBE?1) (17) 　　VN??kBT? 　　由积分形式写成求和形式需要注意： 　　2 　　? 　　q 　　f(q) 　　dq 　　?2?? 　　3 　　? 　　? 　　?2??