引言
经典力学中，粒子的状态由位置 $x$ 与动量 $p$ 这一对实数描述；而在量子力学中，状态是希尔伯特空间（Hilbert space）中的矢量 $|\psi\rangle$，可观测量（observable）则是作用在该空间上的厄米算符（Hermitian operator）。从「数」到「算符」的这一跃升，是量子理论区别于经典理论的核心数学特征。测量一个物理量时，仪器只能给出该算符的某个本征值，而测量后体系塌缩到对应的本征态——因此，量子力学的计算在很大程度上归结为求解算符的本征值问题。本文系统梳理算符形式的基本公设，并以数值对角化方法演示如何在计算机上实际求解一个本征值问题。
理论基础
态矢量与可观测量公设
量子体系的状态由归一化的态矢量 $|\psi\rangle$ 描述，内积 $\langle\phi|\psi\rangle$ 给出概率幅。每个可观测量 $A$ 对应一个厄米算符 $\hat{A}$，满足 $\hat{A} = \hat{A}^\dagger$。厄米性保证其本征值全为实数，且不同本征值对应的本征矢彼此正交。位置与动量算符在坐标表象中分别取：
$$\hat{x} = x, \qquad \hat{p} = -i\hbar\frac{d}{dx}$$
由 $\hat{x}$ 与 $\hat{p}$ 可构造出其余力学量的算符，例如动能算符 $\hat{T} = \hat{p}^2/2m$、角动量算符 $\hat{L}_z = \hat{x}\hat{p}_y - \hat{y}\hat{p}_x$。
对易关系与不确定性原理
算符一般不对易。最基本的对易关系是正则对易关系：
$$[\hat{x}, \hat{p}] \equiv \hat{x}\hat{p} - \hat{p}\hat{x} = i\hbar$$
两个算符对易的充要条件是它们存在一套共同的完备本征矢，可以同时被精确测定。反之，若 $[\hat{A},\hat{B}] = i\hat{C}$，则任意态下的测量涨落满足 Robertson 不确定性关系：
$$\Delta A,\Delta B \geq \frac{1}{2}\left|\langle\hat{C}\rangle\right|$$
取 $\hat{A}=\hat{x}$、$\hat{B}=\hat{p}$ 即得熟知的 $\Delta x,\Delta p \geq \hbar/2$。这不是测量仪器的局限，而是量子态本身的内禀性质。
本征值问题与谱
求解可观测量的本征值问题：
$$\hat{A}|\phi_n\rangle = a_n|\phi_n\rangle$$
全体本征值 ${a_n}$ 构成算符的谱。束缚态体系（如无限深势阱、氢原子）的谱是分立的；散射态体系（如自由粒子）的谱是连续的。由于厄米算符的本征矢构成完备正交基，任意态可展开为 $|\psi\rangle = \sum_n c_n|\phi_n\rangle$，测量得到 $a_n$ 的概率为 $|c_n|^2$。常见力学量的算符形式归纳如下：
| 可观测量 | 经典表达式 | 算符（坐标表象） |
|---|---|---|
| 位置 | $x$ | $\hat{x} = x$ |
| 动量 | $p$ | $\hat{p} = -i\hbar,\partial_x$ |
| 动能 | $p^2/2m$ | $-\dfrac{\hbar^2}{2m}\partial_x^2$ |
| 势能 | $V(x)$ | $V(\hat{x}) = V(x)$ |
| 哈密顿量 | $H = T + V$ | $\hat{H} = -\dfrac{\hbar^2}{2m}\partial_x^2 + V(x)$ |
| 角动量 $z$ 分量 | $xp_y – yp_x$ | $-i\hbar\left(x\partial_y - y\partial_x\right)$ |
数值方法与实现
实际体系中解析可解的本征值问题屈指可数，数值对角化是通用手段。以一维束缚态为例：将区间 $[0, L]$ 离散为 $N$ 个均匀格点，间距 $\Delta x = L/(N+1)$，用三点中心差分把二阶导数离散为
$$\frac{d^2\psi}{dx^2}\bigg|_{x_i} \approx \frac{\psi_{i-1} - 2\psi_i + \psi_{i+1}}{\Delta x^2}$$
于是哈密顿算符化为 $N\times N$ 的实对称三对角矩阵，其本征值即离散化后的能级。下面的程序构造该矩阵，用 numpy.linalg.eigh 对角化，并与无限深势阱的解析解 $E_n = n^2\pi^2\hbar^2/(2mL^2)$ 对比：
import numpy as np
def solve_well_energies(L=1.0, m=1.0, hbar=1.0, N=200, n_levels=4):
"""有限差分离散 + 对角化，求一维无限深势阱的低能级。"""
dx = L / (N + 1)
diag = np.full(N, -2.0) / dx**2 # 主对角线
off = np.full(N - 1, 1.0) / dx**2 # 次对角线
H = -hbar**2 / (2 * m) * (np.diag(diag) + np.diag(off, 1) + np.diag(off, -1))
E_num = np.linalg.eigvalsh(H)[:n_levels] # 实对称矩阵专用求解器
E_exact = hbar**2 * np.pi**2 * np.arange(1, n_levels + 1)**2 / (2 * m * L**2)
return E_num, E_exact
E_num, E_exact = solve_well_energies()
for n, (en, ee) in enumerate(zip(E_num, E_exact), start=1):
print(f"n={n}: 数值 {en:.5f} 解析 {ee:.5f} 相对误差 {abs(en-ee)/ee:.2e}")取 $N=200$ 时，基态能量的相对误差约为 $10^{-6}$ 量级，且误差随格点数加密按 $O(\Delta x^2)$ 减小——这正是中心差分格式的阶数。若改用变分基组展开（如以正弦函数为基），对光滑势场可获得谱精度收敛。
工程案例分析
算符本征值方法在半导体量子器件设计中是日常工具。以 GaAs/AlGaAs 量子阱红外探测器（QWIP）为例：电子在生长方向被几十纳米厚的势垒-势阱交替结构约束，子带能量由一维定态薛定谔方程的分段常数势本征值问题给出。工程师通过调整阱宽与 Al 组分，把基态到第一激发态的跃迁能量精确对准目标红外波段（如 $8\sim10,\mu\text{m}$ 大气窗口），使探测器的光谱响应峰值落在设计值上。设计流程正是上文程序的推广：把分段常数势 $V(x)$ 写入矩阵对角元即可。类似地，量子级联激光器（QCL）的增益谱、量子点量子比特的能级劈裂，都依赖大规模厄米矩阵本征值求解——这正是量子力学公设在工程实践中落地的直接体现。
密度矩阵与开放量子系统
上述算符形式以纯态矢量 $|\psi\rangle$ 为出发点,但在实际系统中与环境的耦合不可避免。密度矩阵 $\hat{\rho} = \sum_k p_k|\psi_k\rangle\langle\psi_k|$ 统一描述了纯态和统计混合态。刘维尔-冯·诺依曼方程:
$$i\hbar\frac{d\hat{\rho}}{dt} = [\hat{H},\hat{\rho}]$$
是量子主方程的最简形式。当耦合到热库时,需加入耗散项,导出 Lindblad 形式的主方程:
$$\frac{d\hat{\rho}}{dt} = -\frac{i}{\hbar}[\hat{H},\hat{\rho}] + \sum_k \left(L_k\hat{\rho}L_k^\dagger - \frac{1}{2}{L_k^\dagger L_k, \hat{\rho}}\right)$$
Lindblad 算符 $L_k$ 描述了不同的退相干通道(自发辐射、退相位、粒子损失)。这一框架是量子计算中理解和抑制退相干的数学基础。在超导量子比特中,能量弛豫时间 $T_1$ 和相位相干时间 $T_2$ 可追溯到具体的 Lindblad 算符与谱密度。例如,Transmon 量子比特中 $T_1 \sim 50\mu\text{s}$、$T_2 \sim 100\mu\text{s}$ 的实测值,对应的主导退相干机制是二能级系统(TLS)缺陷引起的介电损耗和准粒子中毒——这两种机制在 Lindblad 框架下分别建模为纵向和横向退相干通道。
数值扩展:直接对角化实际体系的限制
上文展示的有限差分+对角化方法在单粒子一维问题中效果极佳,但面对多电子体系(维度随粒子数指数增长)时面临”指数墙” 困境。对于 $d$ 个自旋轨道的费米子系统,完整组态相互作用(Full CI)的希尔伯特空间维度为 $binom{d}{N_e}$ ——仅 30 个轨道和 15 个电子即超过 $1.5\times10^8$ 维,任何直接对角化都变得不可行。解决途径包括:
- 密度泛函理论(DFT):将多体问题映射到等效单粒子 Kohn-Sham 方程,复杂度降至 $O(N^3)$。
- 张量网络/矩阵积态(MPS):在密度矩阵重整化群(DMRG)框架下,以一维纠缠结构截断希尔伯特空间,广泛用于强关联电子系统。
- 量子变分本征求解器(VQE):在量子计算机上参数化试探波函数,经典优化器通过量子测量反馈更新参数——IBM 和 Google 已在超导量子芯片上演示了 H$_2$、LiH 分子的基态能量求解。
算符形式为这些进阶方法提供了统一的理论语言:本征值问题仍是核心,但求解方式从数值对角化扩展到变分原理、虚时演化和量子-经典混合算法。
结论与展望
算符形式将量子力学统一为希尔伯特空间上的线性代数：可观测量是厄米算符，测量结果是其本征值，态演化由酉算符生成。对易关系揭示了可同时精确测量的力学量完全集，也划定了不确定性的内禀下限。从数值角度看，绝大多数现实问题都归结为（广义）本征值问题，差分离散与基组展开是两条互补的技术路线。后续文章将看到，引入产生、湮灭算符后，许多本征值问题甚至可以绕过微分方程、纯用代数方法求解。
发表回复