电子运动的准经典模型

                     

贡献者: _Eden_

预备知识 近自由电子模型

   在单电子近似下,利用近自由电子模型或者紧束缚模型,可以有效地分析晶体中电子的行为。和自由电子模型相比,描述电子状态的波函数从平面波过渡到 Bloch(布洛赫)波,准连续的 $\epsilon(k)$ 也分裂为能带。电子在周期性势场中运动的本征态,对应的能量本征值 $\epsilon$,和其 Bloch 波函数的波数 $ \boldsymbol{\mathbf{k}} $($\hbar \boldsymbol{\mathbf{k}} $ 也被称为准动量)有一个对应关系,被称为色散关系,画在 $\epsilon$-${ \boldsymbol{\mathbf{k}} }$ 图上就呈现出一个曲面,由于周期性势场的影响这些曲线画在简约布里渊区中呈现出若干个能带。

   可以说,能带反映了晶体中电子的几乎一切行为,可以用来解释金属、半导体、绝缘体的导电行为。因此,近自由电子模型或紧束缚模型虽然对晶体模型作了许多简化,但在解释晶体的性质方面取得了很大成功。利用能带,我们不仅可以知道晶体中电子的行为,还可以知道对晶体加外场(例如电场或磁场)时电子的变化,从而计算得到晶体对外场的响应。

   如何分析加外场后电子的响应呢?一种方法是求解外加场的情况下薛定谔方程的能量本征解,即

\begin{equation} \left[-\frac{\hbar^2}{2m}\nabla^2+V( \boldsymbol{\mathbf{r}} )+U \right] \psi = E\psi~, \end{equation}
另一种方法是将电子的运动近似当作经典粒子来处理。下面我们介绍的就是这种分析计算方法。

1. 晶体中电子的速度

   首先我们计算 Bloch 态 $\psi_{n \boldsymbol{\mathbf{k}} }( \boldsymbol{\mathbf{r}} )=e^{i \boldsymbol{\mathbf{k}} \cdot \boldsymbol{\mathbf{r}} }u_{n \boldsymbol{\mathbf{k}} }( \boldsymbol{\mathbf{r}} )$ 的速度期待值。其中 $n$ 对应能带的标号。我们有

\begin{equation} \begin{aligned} \boldsymbol{\mathbf{v}} _n( \boldsymbol{\mathbf{k}} )&=\frac{1}{m} \left\langle \psi_{n \boldsymbol{\mathbf{k}} } \right\rvert \hat p \left\lvert \psi_{n \boldsymbol{\mathbf{k}} } \right\rangle =\frac{1}{m}\int \,\mathrm{d}^{3}{ \boldsymbol{\mathbf{x}} } \psi_{n \boldsymbol{\mathbf{k}} }^*( \boldsymbol{\mathbf{x}} )(-i\hbar\nabla)\psi_{n \boldsymbol{\mathbf{k}} }( \boldsymbol{\mathbf{x}} )\\ &=\frac{1}{m}\int \,\mathrm{d}^{3}{ \boldsymbol{\mathbf{x}} } u_{n \boldsymbol{\mathbf{k}} }^*( \boldsymbol{\mathbf{x}} )(-i\hbar\nabla+\hbar \boldsymbol{\mathbf{k}} )u_{n \boldsymbol{\mathbf{k}} }( \boldsymbol{\mathbf{x}} )~, \end{aligned} \end{equation}
而 Bloch 态的能量可以表达为
\begin{equation} \begin{aligned} \epsilon_{n}( \boldsymbol{\mathbf{k}} )&= \left\langle \psi_{n \boldsymbol{\mathbf{k}} } \right\rvert \left[-\frac{\hbar^2}{2m}\nabla^2+V \right] \left\lvert \psi_{n \boldsymbol{\mathbf{k}} } \right\rangle \\&=\frac{1}{m}\int \,\mathrm{d}^{3}{ \boldsymbol{\mathbf{x}} } u_{n \boldsymbol{\mathbf{k}} }^*( \boldsymbol{\mathbf{x}} )(\hbar^2k^2/2-\hbar^2\nabla^2/2-i\hbar \boldsymbol{\mathbf{k}} \cdot \nabla)u_{n \boldsymbol{\mathbf{k}} }( \boldsymbol{\mathbf{x}} )~. \end{aligned} \end{equation}
对比式 2 式 3 可以看到
\begin{equation} \boldsymbol{\mathbf{v}} _n( \boldsymbol{\mathbf{k}} )=\frac{1}{\hbar}\nabla_{ \boldsymbol{\mathbf{k}} }\epsilon_{n}( \boldsymbol{\mathbf{k}} )~. \end{equation}
这表明 Bloch 态电子的速度期待值就等于能带在 $ \boldsymbol{\mathbf{k}} $ 处的斜率。我们也可以通过将 $ \boldsymbol{\mathbf{k}} $ 邻域内的 Bloch 态作光滑叠加得到局域的波包,并且计算波包中心的速度得到上述结果,具体推导可以参考群速度。从上面的分析我们可以得到以下结论:能带的带底和带顶电子速度为零,这是因为它对应能带的极值点,斜率一定为 $0$。另一方面由以上结果可以知道,因为一个 Bloch 态是定态,所以晶体中的电子在无外场的情况下平均速度永远保持不变;这意味着,尽管有电子与离子实之间相互作用的周期性势场,一个理想晶体金属将具有无穷大的电导。然而事实上,任何的晶体都不具有无穷大电导,这是因为晶体结构上的不理想性(存在杂质或缺陷),而且晶格有热振动,电子总会发生散射。

2. 电子运动的准经典模型

   现在我们希望讨论的是,在两次散射之间电子的行为,尤其是在外场作用下电子的行为。为了讨论这个问题,准经典模型是非常有用的。如果对晶体中电子施加一个外力(这个外力来自于外场)$ \boldsymbol{\mathbf{F}} $,那么在 $ \,\mathrm{d}{t} $ 时间内外力对电子做功为

\begin{equation} \boldsymbol{\mathbf{F}} \cdot \boldsymbol{\mathbf{v}} _n( \boldsymbol{\mathbf{k}} ) \,\mathrm{d}{t} ~, \end{equation}
电子的能量也要有相应的变化。设 $ \boldsymbol{\mathbf{k}} $ 变化 $ \,\mathrm{d}{ \boldsymbol{\mathbf{k}} } $ 根据功能原理可以得到
\begin{equation} \boldsymbol{\mathbf{F}} \cdot \boldsymbol{\mathbf{v}} _n( \boldsymbol{\mathbf{k}} ) \,\mathrm{d}{t} = \,\mathrm{d}{ \boldsymbol{\mathbf{k}} } \cdot (\nabla_{ \boldsymbol{\mathbf{k}} }\epsilon_n( \boldsymbol{\mathbf{k}} ))~. \end{equation}
代入式 4 可以得到
\begin{equation} \left( \boldsymbol{\mathbf{F}} -\hbar \frac{\mathrm{d}{ \boldsymbol{\mathbf{k}} }}{\mathrm{d}{t}} \right) \cdot \boldsymbol{\mathbf{v}} _n( \boldsymbol{\mathbf{k}} )=0~. \end{equation}
我们称 $\hbar \boldsymbol{\mathbf{k}} $ 为准动量,则上式表明准动量守恒在平行 $ \boldsymbol{\mathbf{v}} _n( \boldsymbol{\mathbf{k}} )$ 方向上是成立的。事实上可以进一步证明垂直 $ \boldsymbol{\mathbf{v}} _n( \boldsymbol{\mathbf{k}} )$ 方向上也成立,即
\begin{equation} \frac{\mathrm{d}{(\hbar \boldsymbol{\mathbf{k}} )}}{\mathrm{d}{t}} = \boldsymbol{\mathbf{F}} ~, \end{equation}
相关证明可以参考阎守胜·《固体物理基础》 [1]。上面的推导将外场对电子作用作经典处理,这实际上需要满足一定条件:外场随空间和时间是缓慢变化的。从上述结果可以看到,如果施加恒定的外场,电子将在 $ \boldsymbol{\mathbf{k}} $ 空间中向一个方向作匀速运动。

3. 有效质量

   在准经典的非相对论模型中,讨论电子的加速度是有意义的。从式 4 式 8 出发可以得到

\begin{equation} \frac{\mathrm{d}{ \boldsymbol{\mathbf{v}} _n( \boldsymbol{\mathbf{k}} )}}{\mathrm{d}{t}} =\frac{1}{\hbar^2}\nabla_{ \boldsymbol{\mathbf{k}} } \left(\nabla_{ \boldsymbol{\mathbf{k}} }\epsilon_n( \boldsymbol{\mathbf{k}} ) \right) \cdot \hbar \frac{\mathrm{d}{ \boldsymbol{\mathbf{k}} }}{\mathrm{d}{t}} =\frac{1}{\hbar^2}\nabla_{ \boldsymbol{\mathbf{k}} }\nabla_{ \boldsymbol{\mathbf{k}} }\epsilon_n( \boldsymbol{\mathbf{k}} ) \cdot \boldsymbol{\mathbf{F}} ~, \end{equation}
上式中我们省略了指标,$\nabla_{ \boldsymbol{\mathbf{k}} }\nabla_{ \boldsymbol{\mathbf{k}} }\epsilon_n( \boldsymbol{\mathbf{k}} )$ 是一个二阶对称张量,设
\begin{equation} \left[m^{*} \right] ^{-1}{}_{ij}=\frac{1}{\hbar^2} \frac{\partial {}^2\epsilon_n( \boldsymbol{\mathbf{k}} )}{\partial k_i\partial k_j} ~, \end{equation}
那么
\begin{equation} (m^*)_{ij} \frac{\mathrm{d}{ \boldsymbol{\mathbf{v}} _n( \boldsymbol{\mathbf{k}} )}}{\mathrm{d}{t}} {}_j= \boldsymbol{\mathbf{F}} _i~. \end{equation}
由于二阶对称张量可以通过正交变换使其对角化,所以我们可以切换坐标轴为 s 三个主轴。设 $k_x,k_y,k_z$ 为主轴,那么
\begin{equation} \frac{1}{m_\alpha^*}=\frac{1}{\hbar^2} \frac{\partial^{2}{\epsilon_n( \boldsymbol{\mathbf{k}} )}}{\partial{k_\alpha}^{2}} \quad (\alpha=x,y,z)~. \end{equation}
有效质量不一定等于电子的质量,而且在不同的主轴方向上有效质量的大小也可能互不相同,这正是因为晶格中离子实所产生的周期性势场对电子的作用,使得它所 “表现” 出来的 “物理质量” 同 “裸质量” 不同。

   下面我们以紧束缚近似为例分析有效质量。三维简单立方晶格的 $s$ 轨道能带有关系式

\begin{equation} \epsilon( \boldsymbol{\mathbf{k}} )=\epsilon_s-J_0-2J_1 \left[ \cos\left(k_xa\right) + \cos\left(k_ya\right) + \cos\left(k_za\right) \right] ~. \end{equation}
那么代入式 12 可以得到
\begin{equation} m_\alpha^*=\frac{\hbar^2}{2 J_1 a^2 \cos\left(k_\alpha a\right) }~. \end{equation}
如果我们考虑能带底附近的电子,即 $k\approx 0$,那么三个方向的有效质量都为
\begin{equation} m_x^*=m_y^*=m_z^*=\frac{\hbar^2}{2J_1a^2}~. \end{equation}
能带低对应于 $\epsilon( \boldsymbol{\mathbf{k}} )$ 的极小值,因此二次微商大于 $0$,有效质量为正。相反,能带顶部电子的有效质量为负。对于能带顶部附近的电子,即 $k\approx \pm \pi/a$,有效质量为
\begin{equation} m_x^*=m_y^*=m_z^*=-\frac{\hbar^2}{2J_1a^2}~. \end{equation}
一般而言,宽的能带,能量随 $ \boldsymbol{\mathbf{k}} $ 变化更剧烈,那么二次微商更大,有效质量更小;而窄的能带有效质量更大,从紧束缚近似的角度看,此时相邻原子的电子波函数的交叠较小,因此局域性更强,所以不易受外力影响。


[1] ^ 阎守胜. 固体物理基础 北京大学出版社 第三版

                     

© 小时科技 保留一切权利