引言
量子力学是描述微观粒子行为的理论框架，诞生于 20 世纪初物理学界对黑体辐射、光电效应和原子光谱等实验现象的理论探索。1900 年普朗克提出能量量子化假说，1905 年爱因斯坦以光子概念解释光电效应，1913 年玻尔建立氢原子模型，1924 年德布罗意提出物质波假说。1926 年，奥地利物理学家薛定谔（Erwin Schrödinger）在前人工作的基础上，提出了以他名字命名的波动方程，将量子系统的时间演化纳入统一的数学框架。同年，海森堡、玻恩和约当建立了矩阵力学形式。薛定谔随后证明了波动力学与矩阵力学的等价性，完成了量子力学的数学基础构建。
薛定谔方程是量子力学的核心动力学方程，地位等同于经典力学中的牛顿第二定律。它所揭示的波粒二象性、不确定性关系、态叠加原理和量子纠缠等概念，彻底改变了人类对微观世界的认知方式。在凝聚态物理、量子化学、材料科学和量子信息等领域，薛定谔方程仍是理论计算和数值模拟的出发点。
一维无限深势阱（Particle in a 1D Box）是量子力学中最简单而又最为深刻的模型之一。尽管其物理设定看似理想化——将粒子约束在一条线段内、势垒无限高——但它能清晰地展现能量量子化、零点能和波函数驻波性质等量子力学的本质特征。该模型在研究中直接对应量子线、量子点、共轭聚合物中的电子态等诸多实际系统，是理解能带理论和量子受限效应的基础。
理论基础
含时与定态薛定谔方程
含时薛定谔方程：
$$i\hbar \frac{\partial}{\partial t} \Psi(x, t) = \left[ -\frac{\hbar^2}{2m} \frac{\partial^2}{\partial x^2} + V(x) \right] \Psi(x, t)$$
对于势能不显含时间 $V = V(x)$ 的情况，波函数可分离变量 $\Psi(x, t) = \psi(x) e^{-iEt/\hbar}$，得到定态薛定谔方程：
$$-\frac{\hbar^2}{2m} \frac{d^2 \psi}{dx^2} + V(x) \psi = E \psi$$
或写为本征值问题：
$$\hat{H} \psi_n = E_n \psi_n$$
其中 $\hat{H} = -\frac{\hbar^2}{2m} \frac{d^2}{dx^2} + V(x)$ 是哈密顿算符。
一维无限深势阱
势能函数 $V(x) = 0$（$0 < x < L$），$V(x) = \infty$（$x \le 0$ 或 $x \ge L$）。由波函数连续性条件 $\psi(0) = \psi(L) = 0$，得：
$$\psi_n(x) = \sqrt{\frac{2}{L}} \sin\left(\frac{n\pi x}{L}\right), \quad n=1,2,3,\dots$$
能量本征值：
$$E_n = \frac{\hbar^2 \pi^2 n^2}{2mL^2}$$
关键物理结论：能量是分立的（量子化），$E_n \propto n^2$。
以电子（$m = 9.11 \times 10^{-31}$ kg）限制在 $L = 1$ nm 势阱为例，前四个能级的能量为：
| 量子数 $n$ | 能量 $E_n$（eV） | 相对基态倍数 |
|---|---|---|
| 1 | 0.376 | 1 |
| 2 | 1.504 | 4 |
| 3 | 3.384 | 9 |
| 4 | 6.016 | 16 |
本征函数的性质
波函数的归一化系数 $\sqrt{2/L}$ 由 $\int_0^L |\psi_n|^2 dx = 1$ 确定。不同量子数 $n$ 的本征函数满足正交关系：
$$\int_0^L \psi_m^*(x) \psi_n(x) dx = \delta_{mn}$$
这一正交性是厄密算符本征函数的普遍性质，保证了态叠加原理的数学一致性：势阱内任意波函数可展开为本征函数的线性组合 $\Psi(x) = \sum_{n} c_n \psi_n(x)$，展开系数 $c_n = \int \psi_n^* \Psi dx$。在量子测量中，$|c_n|^2$ 给出测量能量得到 $E_n$ 的概率。
期望值与不确定性
位置和动量的期望值可解析计算。对于第 $n$ 个本征态：
$$\langle x \rangle = \int_0^L x |\psi_n|^2 dx = \frac{L}{2}$$
$$\langle x^2 \rangle = \int_0^L x^2 |\psi_n|^2 dx = L^2 \left( \frac{1}{3} - \frac{1}{2\pi^2 n^2} \right)$$
$$\sigma_x = \sqrt{\langle x^2 \rangle - \langle x \rangle^2} = \frac{L}{2\pi n} \sqrt{\frac{\pi^2 n^2}{3} - 2}$$
动量算符 $\hat{p} = -i\hbar d/dx$ 的期望值：$\langle p \rangle = 0$（驻波无净动量流），$\langle p^2 \rangle = (\hbar \pi n / L)^2$。代入海森堡不确定性关系验证：
$$\sigma_x \sigma_p \ge \frac{\hbar}{2}$$
当 $n = 1$ 时，$\sigma_x \sigma_p \approx 0.568 \hbar > \hbar/2$，满足不确定性原理。随着 $n$ 增大，不确定性乘积增大，高激发态波函数在空间上更均匀分布，动量不确定性增加。
对应原理
当量子数 $n$ 很大时，概率密度 $|\psi_n|^2 = (2/L)\sin^2(n\pi x/L)$ 在势阱内快速振荡，其空间平均趋于经典均匀分布 $1/L$。这一行为体现了玻尔对应原理：大量子数极限下量子力学过渡到经典力学。对于无限深势阱，这一过渡是严格成立的：
$$\lim_{n\to\infty} \frac{1}{x_2 - x_1} \int_{x_1}^{x_2} \frac{2}{L}\sin^2\left(\frac{n\pi x}{L}\right) dx = \frac{1}{L}$$
数值实现
import numpy as np
import matplotlib.pyplot as plt
def infinite_well_wavefunction(n, L, N=1000):
x = np.linspace(0, L, N)
psi = np.sqrt(2.0 / L) * np.sin(n * np.pi * x / L)
prob = psi**2
return x, psi, prob
# 绘制前四个能级的波函数与概率密度
L = 1.0
fig, axes = plt.subplots(4, 2, figsize=(12, 16))
for n in range(1, 5):
x, psi, prob = infinite_well_wavefunction(n, L)
axes[n-1, 0].plot(x, psi)
axes[n-1, 0].set_ylabel(f'n={n} Ψ(x)')
axes[n-1, 1].plot(x, prob)
axes[n-1, 1].set_ylabel(f'n={n} |Ψ|²')
plt.tight_layout()
plt.show()
# 能量间距分析
m_e = 9.11e-31
hbar = 1.0546e-34
eV = 1.602e-19
E = lambda n: hbar**2 * np.pi**2 * n**2 / (2 * m_e * L**2) / eV
print("能级间距 ΔEₙ = Eₙ₊₁ - Eₙ (eV):")
for n in range(1, 6):
print(f" n={n}→{n+1}: {E(n+1)-E(n):.3f} eV")能级间距 $\Delta E_n = E_{n+1} - E_n = \frac{\hbar^2\pi^2}{2mL^2}(2n+1)$ 随 $n$ 线性增大。在宏观尺度（$L$ 大）下，能级间距趋于零，量子效应消失，能量表观连续——这解释了日常宏观尺度为何观察不到量子化。
应用与扩展
量子点
半导体量子点（Quantum Dot）本质上是一个三维有限高势阱。当三个空间维度的约束尺寸均小于电子的德布罗意波长时，载流子在三个方向上的能级均分立，形成”人工原子”。CdSe 胶体量子点的尺寸一般在 2–10 nm 范围，其光致发光波长可通过控制粒径来连续调谐——粒径越小，禁带越宽，发光偏蓝。这种尺寸效应已被广泛应用于量子点电视（QLED）、生物荧光标记和单光子源等前沿领域。
有限深势阱与隧穿
实际系统中势垒总是有限高 $V_0$。对于 $E < V_0$ 的束缚态，波函数在经典禁阻区域呈指数衰减 $\psi(x) \propto e^{-\kappa|x|}$（其中 $\kappa = \sqrt{2m(V_0-E)}/\hbar$），而非严格为零。这一"隧穿尾部"使得束缚态数目有限（由 $\sqrt{2mV_0}L/\pi\hbar$ 决定）——因此有限深势阱只能容纳有限个束缚态，这与无限深势阱的无界能谱形成对比。当势垒足够薄或粒子能量足够高时，势垒两侧的波函数可通过隧穿耦合，形成共振隧穿效应。共振隧穿二极管（RTD）正是利用双势垒结构中的量子隧穿实现负微分电阻，应用于太赫兹振荡器和高速逻辑电路。
一维亚晶格与超晶格
周期性排列的势阱构成一维晶格模型。紧束缚近似下，单势阱中的 $n$ 个局域波函数通过势垒间的交叠耦合，形成宽度为 $4t$ 的能带（$t$ 为跃迁积分）。这与克勒尼希-彭尼（Kronig-Penney）模型一脉相承，是理解固体能带形成、金属-绝缘体转变和拓扑绝缘体边界态的基本出发点。半导体超晶格（如 GaAs/AlGaAs）中周期性交替生长的势阱层和势垒层，可通过分子束外延原子级精度控制 $L$ 和 $V_0$，实现能带工程的自由设计，广泛应用于量子级联激光器和红外探测器。
总结与展望
无限深势阱模型以极简的数学形式揭示了量子力学的核心特征——能量量子化、零点能、波函数驻波形式、态叠加原理以及经典极限的对应过渡。从海森堡不确定性关系的定量验证到量子点的尺寸效应，从紧束缚能带形成到共振隧穿器件，这一基础模型的影响贯穿了整个现代物理学与纳米电子学。
在数值计算层面，无限深势阱问题的解析解为验证有限差分法、有限元法和谱方法等数值方案提供了精确的基准测试（benchmark）。以一维谐振子和有限深势阱为起点，可以自然地过渡到二维量子阱、量子环和非线性薛定谔方程的数值求解。
当前研究的前沿方向包括：莫尔超晶格中的平带与关联绝缘态、拓扑量子阱中的边界态输运、以及量子点中自旋量子比特的操控与读出。无限深势阱所提供的清晰物理图像，始终是研究者拓展量子世界认知边界的第一级台阶。
发表回复