原内容请参考哈尔滨工业大学何飞教授:https
文章目录
材料物理基础知识
k空间
可用波矢
(
k
x
,
k
y
,
k
z
)
(k_x,k_y,k_z)
(kx,ky,kz)表示电子的状态,称为k空间或波矢空间。
由量子数n与波矢k的关系
k
=
2
π
n
L
k=\frac{2\pi n}{L}
k=L2πn可知,自由电子的能量为:
E
=
ℏ
2
2
m
(
k
x
2
+
k
y
2
+
k
z
2
)
=
h
2
2
m
L
2
(
n
x
2
+
n
y
2
+
n
z
2
)
E=\frac{\hbar^2}{2m}(k_x^2+k_y^2+k_z^2)=\frac{h^2}{2mL^2}(n_x^2+n_y^2+n_z^2)
E=2mℏ2(kx2+ky2+kz2)=2mL2h2(nx2+ny2+nz2)
- k空间中的每个代表点表示一种能量状态
- k空间中相邻代表点间距为 2 π L \frac{2\pi}{L} L2π
- 每个代表点均匀分布在k空间中
k空间的几何关系
将k空间内的任意代表点看作边长为 2 π L \frac{2\pi}{L} L2π的立方体的中心,则:
- 每个代表点在k空间内的体积为: ( 2 π L ) 3 (\frac{2\pi}{L})^3 (L2π)3
- 将单位体积视为1,则k空间单位体积内状态的数目为: 1 / ( 2 π L ) 3 = ( L 2 π ) 3 = V ( 2 π ) 3 1/(\frac{2\pi}{L})^3=(\frac{L}{2\pi})^3=\frac{V}{(2\pi)^3} 1/(L2π)3=(2πL)3=(2π)3V
其中
V
=
L
3
V=L^3
V=L3
k空间上的等能面
在k空间内,通过能量E和波矢k的数学关系可以发现,具有相同大小的波矢,不论其方向如何,都具有相同的能量。则在k空间内从原点到k空间某点之间,如果距离相同,则这些距离相同的所有代表点都具有相同的能量
E
=
ℏ
2
k
2
2
m
E=\frac{\hbar^2k^2}{2m}
E=2mℏ2k2。此时,以原点为球心,半径为
k
=
2
m
E
ℏ
k=\frac{\sqrt{2mE}}{\hbar}
k=ℏ2mE的球面便是能量E在k空间上的等能面。等能面上有多少个代表点,就说明存在多少个简并。
能级密度(状态密度)
球壳层中可容纳的电子状态总数为:
d
N
=
2
⋅
V
(
2
π
)
3
⋅
4
π
k
2
d
k
\text{d}N=2\cdot\frac{V}{(2\pi)^3}\cdot4\pi k^2\text{d}k
dN=2⋅(2π)3V⋅4πk2dk
其中,
C
=
4
π
V
(
2
m
h
2
)
3
2
C=4\pi V(\frac{2m}{h^2})^{\frac{3}{2}}
C=4πV(h22m)23对于确定的研究对象而言应为常数。
能级密度的物理意义
单位能量间隔范围内所能容纳的电子状态数。
能级密度随能量变化的曲线
注意,上述讨论均处于自由电子情形下,在晶体中,变化情况会变得非常复杂。
费米-狄拉克分布函数
在热平衡条件下,电子处在能量为E状态的概率用F-D分布函数描述: f ( E ) = 1 e x p ( E − E F k B T ) + 1 f(E)=\frac{1}{exp({\frac{E-E_F}{k_BT}})+1} f(E)=exp(kBTE−EF)+11
- f ( E ) f(E) f(E):表征具有能量E的状态被电子占据的概率;
- k B = 1.38 × 1 0 − 23 J / K k_B=1.38\times 10^{-23}\;J/K kB=1.38×10−23J/K;
-
E
F
E_F
EF:费米能,是固体的固有属性。
F-D分布函数受温度的影响
T=0K时的费米-狄拉克分布函数
T=0K时,对应的费米能记为 E F 0 E_F^0 EF0
- 当 E > E F 0 E>E_F^0 E>EF0时, e x p ( E − E F 0 k B T ) → ∞ ⇒ f ( E ) = 0 exp(\frac{E-E_F^0}{k_BT})\rightarrow\infty\Rightarrow f(E)=0 exp(kBTE−EF0)→∞⇒f(E)=0
- 当 E < E F 0 E<E_F^0 E<EF0时, e x p ( E − E F 0 k B T ) → 0 ⇒ f ( E ) = 1 exp(\frac{E-E_F^0}{k_BT})\to0\Rightarrow f(E)=1 exp(kBTE−EF0)→0⇒f(E)=1
即,当温度为绝对零度时,费米能级以下的能量一定被电子占据(存在电子)。
T>0K时的费米-狄拉克分布函数
- 当 E = E F E=E_F E=EF时, f ( E ) = 1 2 f(E)=\frac{1}{2} f(E)=21——能量等于费米能级的电子占有和不占有的概率相等
- 当 E > E F E>E_F E>EF时, f ( E ) < 1 2 f(E)<\frac{1}{2} f(E)<21——能量高于费米能级的状态少部分被电子所占据
- 当
E
<
E
F
E<E_F
E<EF时,
f
(
E
)
>
1
2
f(E)>\frac{1}{2}
f(E)>21——能量低于费米能级的状态大部分被电子所占据
费米分布函数随温度变化是电子受热激发的结果。当T>0开始,自由电子受到热激发,将具有高于
E
F
E_F
EF的能量,但只有在费米能级
E
F
E_F
EF小范围附近内的电子才能够吸收能量,从低于
E
F
E_F
EF的能级跃迁至高于
E
F
E_F
EF的能级。即只有少量高于
E
F
E_F
EF能量的电子存在于这个
E
F
E_F
EF附近的较小范围内。通过计算可得,这个范围从数值上等于
k
B
T
k_BT
kBT。因此,随着温度升高,将有更多电子具有高于
E
F
E_F
EF的能量。
这一结果表明,金属在熔点以下,虽然自由电子都受到热激发,但只有能量在费米能级附近范围内的一小部分电子受到温度的影响,获得更高的能量,而大部分电子仍旧保持较低的能量水平。
费米能级的温度变化
严格来说,随着温度的升高,费米能级会减小,但是由于费米能级随温度变化的数值改变很小,因此在上述费米分布函数图中忽略掉费米能的变化(左移)。
费米能的意义
费米能
E
F
0
E_F^0
EF0是热力学温度0K时电子填充的最高能级。在此温度下,全部电子均占据费米能级以下的能级。对于不同材料而言具有不同的费米能,费米能级是材料自身的属性。
T=0K和T>0K时能级被电子占据的情况示意图
泡利不相容原理
在费米子组成的系统中,不能有两个或两个以上的粒子处在完全相同的状态。即每个能级最多容纳自旋方向相反的两个电子。并且从能量的最低能级逐渐向高能级填充,直到填满所有能级。因此低能级上的电子其实不能够获得热能而激发,而对于稍低于费米能级的电子来说,则有可能获得能量跃迁至费米能级之上的空能级。这说明金属中虽然存在大量的自由电子,但只有在费米能级附近(大约
k
B
T
k_BT
kBT)范围内的少数电子才能够受热激发而跃迁至更高的能级。
费米能级的计算原理
自由电子系统中能量在 E E E和 E + d E E+\text{d}E E+dE之间的电子状态数
d
N
=
Z
(
E
)
d
E
=
C
E
1
/
2
d
E
\text{d}N=Z(E)\text{d}E=CE^{1/2}\text{d}E
dN=Z(E)dE=CE1/2dE
能量为 E E E的状态被电子占有的概率为 f ( E ) f(E) f(E),则 E E E和 E + d E E+\text{d}E E+dE之间的电子数
能够占据能量为E的状态的电子数=存在电子状态数乘以电子在该能级上出现的概率:
d
N
=
Z
(
E
)
f
(
E
)
d
E
=
C
E
1
/
2
d
E
e
x
p
[
(
E
−
E
F
)
/
k
B
T
]
+
1
\text{d}N=Z(E)f(E)\text{d}E=\frac{CE^{1/2}\text{d}E}{exp[(E-E_F)/k_BT]+1}
dN=Z(E)f(E)dE=exp[(E−EF)/kBT]+1CE1/2dE
推导费米能级
该系统中总电子数应当是从零到
E
F
E_F
EF能级之间的电子状态数
d
N
\text{d}N
dN的积分:
当T=0K时,f(E)=1,则
N
=
∫
0
E
F
C
E
d
E
=
2
3
C
(
E
F
)
3
2
=
8
π
V
3
h
3
(
2
m
E
F
)
3
2
N=\int_{0}^{E_F}C\sqrt{E}\text{d}E=\frac{2}{3}C(E_F)^{\frac{3}{2}}=\frac{8\pi V}{3h^3}(2mE_F)^{\frac{3}{2}}
N=∫0EFCEdE=32C(EF)23=3h38πV(2mEF)23
以
n
=
N
V
n=\frac{N}{V}
n=VN表示单位体积中的自由电子数(即电子浓度)。则费米能级的表达式为:
E
F
=
h
2
2
m
(
3
n
8
π
)
2
3
E_F=\frac{h^2}{2m}(\frac{3n}{8\pi})^{\frac{2}{3}}
EF=2mh2(8π3n)32
- 金属的费米能级为几eV到几十eV
T=0 K时自由电子系统内每个电子的平均能量
E
‾
0
=
∫
0
E
F
0
E
d
N
N
=
∫
0
E
F
0
C
E
3
/
2
d
E
N
=
3
5
E
F
0
\overline{E}_0=\frac{\int_0^{E_F^0}{E\text{d}N}}{N}=\frac{\int_{0}^{E_F^0}{CE^{3/2}\text{d}E}}{N}=\frac{3}{5}E_F^0
E0=N∫0EF0EdN=N∫0EF0CE3/2dE=53EF0
这说明0 K时自由电子仍具有相当大的动能,这是电子满足泡利不相容原理的缘故。
T>0 K时费米能级的计算
能量高于
E
F
E_F
EF的能级上可能被电子占有,能量低于
E
F
E_F
EF的能级上可能未被电子占满,故
f
(
E
F
)
f(E_F)
f(EF)不是常数。
结论:
当T>0K时,
E
F
=
E
F
0
[
1
−
π
2
12
(
k
B
T
E
F
0
)
2
]
E_F=E_F^0[1-\frac{\pi^2}{12}(\frac{k_BT}{E_F^0})^2]
EF=EF0[1−12π2(EF0kBT)2]
由此可见,温度升高,
E
F
E_F
EF略微降低。由于费米能级的量级为
1
0
−
5
10^{-5}
10−5,因此温度造成费米能级下降的影响很小。
T>0 K时电子的平均能量
E
‾
=
3
5
E
F
0
[
1
+
5
π
2
12
(
k
B
T
E
F
0
)
2
]
\overline E=\frac{3}{5}E_F^0[1+\frac{5\pi^2}{12}(\frac{k_BT}{E_F^0})^2]
E=53EF0[1+125π2(EF0kBT)2]
费米面与费米半径
原子结合成晶体时的电子共有化运动
- 在构成晶体之前,孤立原子中的电子只受原子核的束缚绕核运动,该电子称为原子束缚态电子。
- 孤立原子靠近构成晶体时,一个原子的电子会受到相邻原子核的作用转移到相邻原子中,从而可以在整个晶体中运动(即自由电子,不再受某个原子核的束缚),称为电子共有化运动。
- 共有化电子的波函数有一定的交叠(因为电子壳层有一定的交叠),进而引起电子的共有化运动。
- 相邻原子外层电子的波函数交叠较大,因此相邻原子最外层交叠最多,内壳层交叠较少,最外层电子共有化最显著。
晶体中电子波的传播
一维晶体中势能的周期性方程
U
(
x
)
=
U
(
x
+
n
a
)
U(x)=U(x+na)
U(x)=U(x+na)
晶体中电子运动特点
- 从一个原子过渡到另一个原子的共有化运动——自由电子
- 在一个原子附近运动——原子束缚态电子
通常将在晶体内运动的电子称为准自由电子:同时具备自由电子和原子束缚态电子的运动特点。
布洛赫定理
一维晶体中准自由电子的薛定谔方程
−
ℏ
2
2
m
d
2
φ
(
x
)
d
x
2
+
U
(
x
)
φ
(
x
)
=
E
φ
(
x
)
-\frac{\hbar^2}{2m}\frac{\text{d}^2\varphi(x)}{\text{d}x^2}+U(x)\varphi(x)=E\varphi(x)
−2mℏ2dx2d2φ(x)+U(x)φ(x)=Eφ(x)
其中
U
(
x
)
=
U
(
x
+
n
a
)
U(x)=U(x+na)
U(x)=U(x+na)
布洛赫定理的解(布洛赫波)
周期性势场的薛定谔方程的解必定有如下的形式:
φ
k
(
x
)
=
u
k
(
x
)
e
x
p
(
i
k
x
)
\varphi_k(x)=u_k(x)exp(ikx)
φk(x)=uk(x)exp(ikx)
这个解就是在周期性势场中描述电子状态的波函数。该函数是由平面波
e
x
p
(
i
k
x
)
exp(ikx)
exp(ikx)和周期性函数
u
k
(
x
)
u_k(x)
uk(x)相乘的调幅平面波,称为布洛赫波,该函数称为布洛赫函数。
其中,
u
k
(
x
)
=
u
k
(
x
+
n
a
)
u_k(x)=u_k(x+na)
uk(x)=uk(x+na)
自由电子的平面波
若
u
k
(
x
)
u_k(x)
uk(x)为常数,则反映的是自由电子的平面波形式。
晶体中准自由电子出现的概率的周期性变化
∣ φ k ( x ) ∣ 2 = ∣ u k ( x ) ∣ 2 = ∣ u k ( x + n a ) ∣ 2 = ∣ φ k ( x + n a ) ∣ 2 |\varphi_k(x)|^2=|u_k(x)|^2=|u_k(x+na)|^2=|\varphi_k(x+na)|^2 ∣φk(x)∣2=∣uk(x)∣2=∣uk(x+na)∣2=∣φk(x+na)∣2
- 因此自由电子在空间中等概率出现,即作自由运动。
晶体势场近似
周期性势场中电子波函数的具体表达式很难求解,因此需要对势场进行近似:
- 近自由电子近似
- 紧束缚近似
近自由电子近似
该近似理论认为,晶格内的势场起伏小,晶体内的准自由电子近似在一个恒定的势场内运动,仅受到离子实势场的微扰。
- 适用范围:价电子比较自由的金属,比如碱金属、金、银、铜等最外层电子。
周期性势场 U ( x ) U(x) U(x)的傅里叶级数展开分成的两个部分
U ( x ) = U 0 + U ′ U(x)=U_0+U' U(x)=U0+U′
- U 0 U_0 U0:势能的平均值——对应自由电子的运动
-
U
′
U'
U′:势能随坐标周期性变化的部分(微扰)
准自由电子的微扰理论
{
φ
k
=
φ
k
(
0
)
+
Δ
φ
k
E
k
=
E
k
(
0
)
+
Δ
E
k
\begin{cases} \varphi_k=\varphi_k^{(0)}+\Delta\varphi_k\\ E_k=E_k^{(0)}+\Delta E_k \end{cases}
{φk=φk(0)+ΔφkEk=Ek(0)+ΔEk
φ
k
(
0
)
\varphi_k^{(0)}
φk(0)和
E
k
(
0
)
E_k^{(0)}
Ek(0)是零级近似的波函数和能量。
Δ
φ
k
\Delta \varphi_k
Δφk和
Δ
E
k
\Delta E_k
ΔEk是考虑微扰后波函数和能量的修正值。
- 零级近似:周期性势场中,仅考虑恒定势场 U 0 U_0 U0部分作用的结果。
(零级近似实际忽略了晶体的周期性势场,即量子自由电子理论中自由电子的波函数和能量)
- 微扰后的修正值使准自由电子的能量E与波矢k的取值对应关系发生了变化。
布洛赫函数对应的能量受周期性势场的影响
E
k
=
E
k
(
0
)
+
Δ
E
=
{
ℏ
2
2
m
k
2
k
≠
±
n
π
a
⇒
Δ
E
=
0
(
自由电子
)
ℏ
2
2
m
k
2
±
Δ
E
k
=
±
n
π
a
⇒
Δ
E
≠
0
E_k=E_k^{(0)}+\Delta E= \begin{cases} \frac{\hbar^2}{2m}k^2\hspace{5mm}k\neq\pm\frac{n\pi}{a}\Rightarrow\Delta E=0\;(自由电子)\\ \frac{\hbar^2}{2m}k^2\pm\Delta E\hspace{5mm}k=\pm\frac{n\pi}{a}\Rightarrow\Delta E\neq0 \end{cases}
Ek=Ek(0)+ΔE={2mℏ2k2k=±anπ⇒ΔE=0(自由电子)2mℏ2k2±ΔEk=±anπ⇒ΔE=0
一维k空间的布里渊区及其特点
- 当波矢 k = ± n π a k=\pm\frac{n\pi}{a} k=±anπ时,能量出现不连续
- 能量不连续的点把k空间分成许多区域,称为布里渊区。这些能量不连续的点称为布里渊区的边界。
一维k空间的布里渊区的特点
L
L
L:晶体长度
a
a
a:点阵常数
N
N
N:晶体内原胞数量
- 每个波矢在k空间占有的线度为 2 π L = 2 π N a \frac{2\pi}{L}=\frac{2\pi}{Na} L2π=Na2π
- 每个布里渊区含有N个状态:由于每个布里渊区的宽度为 2 π a \frac{2\pi}{a} a2π,故每个布里渊区含有的状态个数为 ( 2 π a ) / ( 2 π N ) = N (\frac{2\pi}{a})/(\frac{2\pi}{N})=N (a2π)/(N2π)=N
- 由2可知,每个布里渊区含有 N N N个能级
- 如果考虑电子自旋,每个能带则能够容纳
2
N
2N
2N个电子。
二维k空间的布里渊区倒格子绘制
每个布里渊区的面积相等。
二维k空间能带交叠现象
二维k空间有可能发生能带之间的交叠,因此不一定存在禁带。
三维k空间的布里渊区特点
- 布里渊区内能量准连续分布,边界处发生突变
- 同一晶格内,尽管各布里渊区形状不同,但体积一定相同
- 每个布里渊区对应一个能带,每个波矢k对应一个能级
- 三维晶体布里渊区同样会出现能带交叠