孙会元教授主编的固体物理基础第五章固体的输运现象课件5.1 玻尔兹曼方程

1、5.1 玻尔兹曼方程玻尔兹曼方程本节主要内容：本节主要内容：玻尔兹曼方程和弛豫时间近似玻尔兹曼方程和弛豫时间近似5.1 玻尔兹曼方程玻尔兹曼方程 0()1( )1kbkk tfe 表示表示n电子体系电子体系在在热平衡态热平衡态(温度为温度为t)时，能量为时，能量为k的单电子本征态被一个电子占据的的单电子本征态被一个电子占据的概率，亦即该电子态的平均电子数。概率，亦即该电子态的平均电子数。0()kf 在第一章我们讨论金属自由电子气体的热性在第一章我们讨论金属自由电子气体的热性质时，已经知道在热平衡状态下，即温度均匀且质时，已经知道在热平衡状态下，即温度均匀且无外场作用时，电子系统的分布函数为费米

2、分布无外场作用时，电子系统的分布函数为费米分布函数，为了和非平衡分布函数区分，把费米分布函数，为了和非平衡分布函数区分，把费米分布函数表示为：函数表示为： 玻尔兹曼方程玻尔兹曼方程 显然，对于均匀体系，显然，对于均匀体系， 与电子所在晶体与电子所在晶体内位置内位置 r 无关。无关。0( )kf 要想讨论输运问题，必须知道微扰状态下要想讨论输运问题，必须知道微扰状态下的分布函数，如何确定非平衡状态下电子的分的分布函数，如何确定非平衡状态下电子的分布函数呢？布函数呢？ 玻尔兹曼方程是用来研究非平衡状态下电子玻尔兹曼方程是用来研究非平衡状态下电子的分布函数的方程。的分布函数的方程。0()1( )1k

3、bkk tfe 由于玻尔兹曼方程比较复杂，我们只限于讨由于玻尔兹曼方程比较复杂，我们只限于讨论论电子的等能面是球面，且电子所经历的碰撞为电子的等能面是球面，且电子所经历的碰撞为弹性散射弹性散射(散射前后能量相等，仅波矢的方向有散射前后能量相等，仅波矢的方向有所改变所改变)以及弱场的情况以及弱场的情况(亦即可将外场的作用当亦即可将外场的作用当成微扰，从而非平衡的稳态分布相对于平衡分布成微扰，从而非平衡的稳态分布相对于平衡分布偏离甚少偏离甚少)。 在上述假定下在上述假定下,非平衡的稳态分布可以表示为：非平衡的稳态分布可以表示为：011,ffff为小量。1.分布函数的变化分布函数的变化 电子分布函数

4、电子分布函数 f 是波矢是波矢 、空间坐标、空间坐标 和时和时间间t的函数的函数。亦即，在。亦即，在t时刻，在单位体积晶体时刻，在单位体积晶体内位置内位置 附近找到一个波矢为附近找到一个波矢为 的电子的几率的电子的几率是：是：krrk( , , )ff r k t 电子分布函数电子分布函数 f 与位置与位置 有关系，通常是由有关系，通常是由于于化学不均匀性化学不均匀性引起的。其原因引起的。其原因一方面一方面是由于化是由于化学成分不均匀学成分不均匀(成分不同化学势不同成分不同化学势不同)引起的化学引起的化学势势(费米能费米能)依赖于位置；依赖于位置；另一方面另一方面是由于有温度是由于有温度梯度引

5、起的化学势梯度引起的化学势(费米能费米能)依赖于位置。我们这依赖于位置。我们这里仅考虑后者里仅考虑后者,即：即：r温度梯度温度梯度f 变化变化变化变化r化学势变化化学势变化 电子分布函数电子分布函数f 与波矢与波矢 有关系，也就是与有关系，也就是与能量有关系，从费米分布函数的表达式就可以能量有关系，从费米分布函数的表达式就可以理解。理解。kd()dkfe evbt 在外电场在外电场 和磁场和磁场 中中,电子的运动规律是：电子的运动规律是：eb 电子分布函数电子分布函数f 与时间与时间t有关系，是因为外力的有关系，是因为外力的作用使得波矢依赖于时间，即：作用使得波矢依赖于时间，即： eb，f 变

6、化变化变化变化k能量变化能量变化 如果不存在碰撞如果不存在碰撞，t时刻，在相空间时刻，在相空间(以波矢以波矢 坐标坐标 为变量组成的空间为变量组成的空间) 、 处的电子必来自处的电子必来自 t - dt 时刻时刻 处，即处，即不存在碰撞时不存在碰撞时来来自于电子的自于电子的漂移过程漂移过程。所以：。所以： krrk,rrdt kkdt( , , )(,)f r k tf rrdt kkdt tdtd()dkfe evbt ( , , )(,)f r k tf rrdt kkdt tdt如果不存在碰撞如果不存在碰撞 实际上实际上,由于由于存在碰撞存在碰撞,dt 时间内时间内,从从 出发的电子未必

7、都能到达出发的电子未必都能到达 、 处处,自然，自然， 、 处处的电子由于的电子由于存在碰撞存在碰撞也并不全部来自也并不全部来自 处。所以占据几率随时间的变化包括两部分：处。所以占据几率随时间的变化包括两部分：一部分是一部分是迁移过程迁移过程(粒子的漂移和扩散粒子的漂移和扩散);另一部另一部分是与相互作用的对方分是与相互作用的对方相碰撞相碰撞而出现的而出现的.所以：所以： ,rrdt kkdtkrr,rrdt kkdt k( , , )(,)ff r k tf rrdt kkdt tdtdtt碰漂移作用引起的分布函数的变化漂移作用引起的分布函数的变化fffttt碰漂碰撞引起的分布函数的变化碰撞

8、引起的分布函数的变化所以，电子分布函数的变化可表示为：所以，电子分布函数的变化可表示为： 通常我们研究输运行为时，讨论的都是定态通常我们研究输运行为时，讨论的都是定态过程，即假定：过程，即假定：0ft 注意，这里提到的稳态情形注意，这里提到的稳态情形 和费和费米米狄拉克函数狄拉克函数 所定义的平衡态之间是有所定义的平衡态之间是有区别的。区别的。0ft0( )kf 从上面的讨论可以看出，要研究输运行为，从上面的讨论可以看出，要研究输运行为，必须对必须对漂移项漂移项(drift term)和和碰撞项碰撞项(collision term)(或或散射项散射项(scattering term)也叫也叫相

9、互作用项相互作用项)有所了解。有所了解。1).漂移项漂移项(先不考虑碰撞时先不考虑碰撞时) 按照分布函数的定义，在按照分布函数的定义，在t时刻，总电子数中时刻，总电子数中存在于相空间体积元存在于相空间体积元 内的电子数为：内的电子数为：d( )( , , )( , , ;,; )xyzf t df r k t df x y z k k k t d xyzddrdkdxdydzdk dk dkfffttt碰漂 在在 t+dt 时刻时刻,空间体积元空间体积元 含有含有这样一些这样一些电子电子,这些电子在这些电子在t 时刻的坐标是时刻的坐标是 其波矢分量为其波矢分量为 亦即在亦即在t+dt 时刻有：

10、时刻有：d,xyzx v dt y v dt z v dt,xxyyzzkk dt kk dt kk dt()f t dt d 在扩散的影响和外力的作用下，坐标空间和在扩散的影响和外力的作用下，坐标空间和动量空间中的占据几率发生变化。动量空间中的占据几率发生变化。也就是也就是这部分电子是漂移过来这部分电子是漂移过来的，所以：的，所以：(,;,; )xyzxxyyzzf x vdt y v dt z vdt kk dt kk dt kkdt t dxyzxyzxyzfffffffvvvkkktxyzkkk漂()( )ff tdtf tdtt左边(,;,; )xyzxxyyzzf x vdt y

11、v dt z vdt kkdt kk dt kkdt t右xyzxyzxyzfffffffvvvkkktxyzkkk 漂推导：推导： 利用多元函数的泰勒展开利用多元函数的泰勒展开,且只取到且只取到dt的线性项的线性项( , , ;,; ) xyzxyzf x y z k k k tv dtv dtv dtxyz ( , , ; , ; )xyzxyzxyzk dtk dtk dtf x y z k k k tkkk(,)( ,) ( ( , )f xx yyf x yxyf x yxy ， 左右相等，即可得到上式。左右相等，即可得到上式。xyzxyzxyzfffffffvvvkkktxyzkk

12、k漂 fffrktrk漂或或：漂移项漂移项=外场作用力引起的外场作用力引起的电子波矢的漂移电子波矢的漂移+速度引起的速度引起的电子位置的漂移电子位置的漂移rkrfkf 漂移描述了在两次碰撞之间的纯动力学行为，漂移描述了在两次碰撞之间的纯动力学行为，并不导致不可逆因素，为此必须考虑碰撞项。并不导致不可逆因素，为此必须考虑碰撞项。 碰撞对应于不可逆过程，它迫使系统趋于平碰撞对应于不可逆过程，它迫使系统趋于平衡分布。衡分布。 对于碰撞项的考虑，只有知道了各碰撞机制，对于碰撞项的考虑，只有知道了各碰撞机制，才能够精确地计算碰撞引起的分布函数的变化。才能够精确地计算碰撞引起的分布函数的变化。 但是，为使

13、问题简化，我们可以仿照第一章但是，为使问题简化，我们可以仿照第一章的做法，引入弛豫时间的做法，引入弛豫时间(k) 作为表征各个相互作为表征各个相互作用的唯象量。作用的唯象量。 2).相互作用项相互作用项(碰撞项或散射项碰撞项或散射项)： 第一章弛豫时间是作为碰撞的几率，或相继两第一章弛豫时间是作为碰撞的几率，或相继两次碰撞间的平均时间引进的。次碰撞间的平均时间引进的。 现在我们从分布函数的角度来看弛豫时间这现在我们从分布函数的角度来看弛豫时间这一假定的意义。一假定的意义。 设没有施加外场时设没有施加外场时,电子系统处于平衡状态电子系统处于平衡状态,在外场的作用下在外场的作用下,由于由于漂移漂移

14、,电子系统进入非平衡电子系统进入非平衡态态,外场去掉后外场去掉后,由于由于碰撞碰撞将使得系统恢复平衡。将使得系统恢复平衡。0 ffft 显然系统恢复平衡的快慢显然系统恢复平衡的快慢( ) 将与系统偏将与系统偏离平衡态的程度离平衡态的程度( )以及碰撞的频度以及碰撞的频度( )有有关关.在在对平衡态的偏离较小时对平衡态的偏离较小时,我们假定系统恢复我们假定系统恢复平衡的快慢比例于系统偏离平衡态的程度以及碰平衡的快慢比例于系统偏离平衡态的程度以及碰撞的频度撞的频度,即：即： ft0ff1一般称上述近似为一般称上述近似为线性近似线性近似，负号源于偏离将随，负号源于偏离将随时间的增加而减小。时间的增加

15、而减小。从该式出发可以进一步得到从该式出发可以进一步得到 00() fftff 即电子的分布函数偏离了平衡分布后即电子的分布函数偏离了平衡分布后,系统依系统依赖碰撞恢复平衡分布的弛豫过程随时间以指数形赖碰撞恢复平衡分布的弛豫过程随时间以指数形式变化式变化,弛豫时间弛豫时间 为这一过程的时间常数为这一过程的时间常数. 是是平衡时的费米狄拉克分布函数。平衡时的费米狄拉克分布函数。0f01fff11( )0(0)00() fttftfftff所以：00() fftff011( )(0)tfff tf te 由于上述过程是在无外场，无温度梯度下完由于上述过程是在无外场，无温度梯度下完全全由碰撞由碰撞引

16、起的，即：引起的，即：0ft漂fffttt碰漂01fffft 碰得到：得到：所以由：所以由： 由于研究输运行为时，讨论的都是定态过程由于研究输运行为时，讨论的都是定态过程,即假定：即假定：0ft01fffft 碰rkfffrkrfkftrk 漂fffttt碰漂这样我们就得到了稳态条件：这样我们就得到了稳态条件：00rkffrfkf玻尔兹曼方程玻尔兹曼方程0rkffrfkf或：2.玻尔兹曼方程玻尔兹曼方程00rkffrfkf这是所有输运现象的基本方程。这是所有输运现象的基本方程。 由于公式两边都有受微扰的占据几率由于公式两边都有受微扰的占据几率 f ,所以所以在该形势下求解很困难在该形势下求解很

17、困难.然而对于所有的实际情然而对于所有的实际情况况,外界影响给占据几率带来的变化都很小外界影响给占据几率带来的变化都很小,所以所以可近似可近似用非微扰占据几率的梯度替换方程右边用非微扰占据几率的梯度替换方程右边。0rkffrfkf或：000 rkffrkff或 ：0010rkfffrk 总之有了外场和温度梯度，系统的分布才总之有了外场和温度梯度，系统的分布才会偏离平衡，无休止的漂移；有了碰撞，就会会偏离平衡，无休止的漂移；有了碰撞，就会使漂移受到遏制，被限制在一定程度而达到稳使漂移受到遏制，被限制在一定程度而达到稳定分布。定分布。 考虑仅有外电场和温度梯度存在的情况下，考虑仅有外电场和温度梯度

18、存在的情况下，漂移项变为：漂移项变为：()ekevb 在施加外磁场时在施加外磁场时,外磁场主要是引起电子波矢外磁场主要是引起电子波矢的漂移：的漂移：000rkfffeerfkfrttk 1kr01() ()kkffekfvbkk00001kkffvffkk而：所以，在电场、磁场和温度梯度都存在时，玻尔兹曼方程为：所以，在电场、磁场和温度梯度都存在时，玻尔兹曼方程为：011() ()()kkkfffeekfvbvbkkk所以：00kvffk000rkfffeerfkfrttk 00011()0kfffffeeertvbtkk 温度梯度温度梯度电场电场磁场磁场外磁场引起的变化外磁场引起的变化 外电

19、场和温度梯度存在的情况下的变化外电场和温度梯度存在的情况下的变化 关于碰撞项的讨论，也可以从量子力学出发展开。关于碰撞项的讨论，也可以从量子力学出发展开。 由此方程可得分布函数，从而可以进一步讨论金属由此方程可得分布函数，从而可以进一步讨论金属电导率随温度的变化，以及热导率、热电效应、霍尔电导率随温度的变化，以及热导率、热电效应、霍尔效应和磁电阻效应等效应和磁电阻效应等 00011()0kfffffeeertvbtkk 由于声子或杂质的散射，粒子可以从由于声子或杂质的散射，粒子可以从 k 态跃迁到态跃迁到 k /态，也可以从态，也可以从 k /态跃迁到态跃迁到 k 态。态。 设设w(k，k /

20、)表示单位体积中在单位时间由表示单位体积中在单位时间由 k 态跃迁态跃迁到到 k / 态的散射几率，态的散射几率，w(k /，k)表示单位体积中在单表示单位体积中在单位时间由位时间由 k / 态跃迁到态跃迁到 k 态的散射几率。态的散射几率。 其中，只考虑了自旋相同态之间的跃迁，其中，只考虑了自旋相同态之间的跃迁，1-f( k /)表示表示k /态未被占据的概率。态未被占据的概率。3,1( )1()(2 )k kaf kf kwdk w(k，k /)表示由表示由 k 态跃迁到态跃迁到 k / 态的散射几率，态的散射几率， w(k /，k)表示由表示由 k / 态跃迁到态跃迁到 k 态的散射几率

21、。态的散射几率。 假定在散射过程中电子的假定在散射过程中电子的自旋不变自旋不变，那么，那么单位体积单位体积中在单位时间中在单位时间由由 k 态散射到所有自旋相同的态散射到所有自旋相同的 k /态的态的净减几率为净减几率为 同样可以得到单位体积中在单位时间由同样可以得到单位体积中在单位时间由 k /态散射到态散射到所有自旋相同的所有自旋相同的 k 态的净增几率。态的净增几率。 ( , , )f r k tbat碰 碰撞导致的分布函数的变化就是净增几率和净减几碰撞导致的分布函数的变化就是净增几率和净减几率之差，即率之差，即 3,1()1( )(2 )k kbf kf kwdk 单位体积中在单位时间由单位体积中在单位时间由 k /态散射到所有自旋相同态散射到所有自旋相同的的 k 态的净增几率为态的净增几率为 则分布函数随时间变化的玻尔兹曼方程为则分布函