半导体物理基础：半导体能带理论_4.周期性势场与能带结构.docxVIP

PAGE1

PAGE1

4.周期性势场与能带结构

在上一节中，我们讨论了半导体的基本能级结构和电子的状态。现在，我们将进一步探讨周期性势场在半导体中的作用及其对能带结构的影响。

4.1周期性势场的基本概念

在固体物理中，晶体的原子排列是周期性的，这意味着晶体中的势场也是周期性的。周期性势场可以用以下数学表达式表示：

V

其中，R是晶体的晶格向量。这个势场的周期性对电子的运动产生了重要的影响。

4.2布洛赫定理

布洛赫定理是描述电子在周期性势场中运动的基本定理。根据布洛赫定理，电子的波函数可以表示为：

ψ

其中，k是波矢，uk

u

布洛赫定理的证明可以通过解薛定谔方程来得到。考虑电子在周期性势场中的运动，薛定谔方程可以写成：

?

由于势场Vr

V

其中，G是晶格的倒格子向量。将布洛赫波函数代入薛定谔方程中，可以得到：

?

经过一些数学推导，最终可以得到uk

4.3能带的形成

布洛赫定理告诉我们，电子在周期性势场中的波函数具有特定的形式。这使得电子的能量E与波矢k之间的关系不再是简单的线性关系，而是形成了一系列能带。能带是电子在特定波矢下的能量范围。

能带的形成可以通过求解布洛赫方程来得到。对于一维周期性势场，布洛赫方程可以简化为：

?

其中，Vx是周期性势场，u

u

这里，a是晶格常数。求解这个方程可以得到能带结构。

4.4一维周期性势场的能带结构

为了更好地理解能带的形成，我们可以通过一个具体的例子来探讨一维周期性势场的能带结构。假设势场Vx是一个周期为a

V

并且Vx

我们可以通过数值方法求解布洛赫方程。以下是一个使用Python和NumPy库的示例代码：

importnumpyasnp

importmatplotlib.pyplotasplt

fromscipy.sparseimportdiags

fromscipy.sparse.linalgimporteigsh

#定义势场

defperiodic_potential(x,a,V0):

returnV0*(xa/2)

#定义晶格常数和势垒高度

a=1.0

V0=10.0

#定义波矢范围

k_values=np.linspace(-np.pi/a,np.pi/a,100)

#定义网格和步长

N=1000

dx=a/N

x=np.linspace(0,a,N,endpoint=False)

#定义系统的哈密顿矩阵

defhamiltonian(k,a,V0,dx,N):

V=periodic_potential(x,a,V0)

kinetic=-0.5*(1/dx**2)*(2*np.eye(N)-np.eye(N,k=1)-np.eye(N,k=-1))

potential=np.diag(V)

hamiltonian_matrix=kinetic+potential+k**2*(0.5*np.eye(N))

returnhamiltonian_matrix

#计算能带

energies=[]

forkink_values:

H=hamiltonian(k,a,V0,dx,N)

eigenvalues,_=eigsh(H,k=10,which=SA)#求解前10个最小的本征值

energies.append(eigenvalues)

energies=np.array(energies)

#绘制能带图

plt.figure(figsize=(10,6))

foriinrange(10):

plt.plot(k_values,energies[:,i],label=fBand{i+1})

plt.xlabel(rWavevector$k$)

plt.ylabel(rEnergy$E$)

plt.title(rEnergyBandsfor1DPeriodicPotential)

plt.legend()

plt.grid(True)

plt.show()

这段代码首先定义了一个周期性势场Vx，然后通过数值方法构建哈密顿矩阵，并求解其本征值。最后，绘制出能带图，显示不同波矢k

4.5二维和三维周期性势场的能带结构

在二维和三维晶体中，周期性势场的能带结构更为复杂。二维晶体的能带结构可以通过类似的方法求解，但需要考虑两个方向的波矢k=kx,

对于二维周期性势场，可以使用数值方法求解布洛赫