2. 量子力学基础
材料计算之所以能从原子种类和位置出发预测电子结构、总能量与力学性质，根本原因在于电子等微观粒子不服从经典牛顿力学，而服从量子力学。对材料计算而言，量子力学不是背景知识，而是所有第一性原理方法的“操作系统”。本节围绕五个核心概念展开：波函数、薛定谔方程、算符、本征态，以及把量子力学用于真实材料时必须引入的基本假设。
波函数：微观状态的数学载体
在经典力学中，一个粒子的状态由位置和动量同时确定。但在量子力学中，由于不确定性原理，位置与动量不能同时被任意精确地确定。于是，微观粒子的状态改由一个复值函数描述，这就是波函数，通常记作 Ψ(r, t)。波函数本身没有直接的物理意义，真正有物理意义的是它的模平方 |Ψ(r, t)|²，它表示在时刻 t、位置 r 附近发现该粒子的概率密度。
这意味着量子力学给出的不是“粒子一定在哪里”，而是“粒子出现在各处的概率分布”。对材料计算来说，我们最关心的是电子在原子核周围形成的电子密度分布，而电子密度正是由波函数的模平方贡献得到的。例如，一个归一化的单电子波函数满足 ∫|Ψ|²dr = 1，表示电子在整个空间中被找到的总概率为 1。归一化条件不是数学装饰，而是保证概率解释成立的前提。
波函数还必须满足一些基本要求：单值、有限、连续，且一阶导数通常连续。这些条件在求解具体问题时，会直接决定哪些能量值是允许的，哪些波函数是物理上可接受的。
薛定谔方程：量子状态如何随时间与空间演化
波函数不能任意给定，它必须满足薛定谔方程。含时薛定谔方程写作：
iħ ∂Ψ(r,t)/∂t = Ĥ Ψ(r,t)其中 ħ 是约化普朗克常数，Ĥ 是哈密顿算符，代表系统的总能量。对材料计算中最常见的定态问题，我们关心不随时间变化的稳定状态，于是得到定态薛定谔方程：
Ĥ Ψ(r) = E Ψ(r)这里 E 是系统的能量本征值。这个方程的形式非常关键：算符 Ĥ 作用在波函数 Ψ 上，得到的结果等于一个常数 E 乘以同一个波函数。满足这种关系的 Ψ 称为哈密顿算符的本征态，对应的 E 称为本征值。材料计算中求解电子结构，本质上就是在给定原子核构型下求解这个方程，得到一系列电子本征态和本征能量。
算符：物理量的数学表示
在量子力学中，每一个可观测物理量都对应一个算符。算符是一种“作用规则”，它作用在波函数上，产生新的函数。常见算符包括：
- 位置算符：r̂ = r，即乘以位置坐标；
- 动量算符：p̂ = −iħ∇；
- 动能算符：T̂ = −(ħ²/2m)∇²；
- 哈密顿算符：Ĥ = T̂ + V̂，其中 V̂ 是势能算符。
对材料计算而言，哈密顿算符的具体形式决定了问题的难度。以单电子在原子核势场中运动为例，Ĥ = −(ħ²/2m)∇² − Ze²/(4πε₀r)，第一项是电子动能，第二项是电子与原子核的库仑吸引势。若系统包含多个电子，还必须加入电子之间的库仑排斥项，问题立刻从单电子方程变成多电子耦合方程，这也是第一性原理计算必须引入近似方法的根本原因。
算符还有一个重要性质：如果两个算符对易，即 [Â, B̂] = ÂB̂ − B̂Â = 0，那么它们对应的物理量可以同时被精确确定。例如，哈密顿算符与动量算符在某些对称条件下对易，意味着能量和动量可以同时有确定值。相反，位置算符与动量算符不对易，[x̂, p̂] = iħ，这正是不确定性原理的数学根源。
本征态与本征值：材料计算最常求解的对象
本征态是量子力学中最核心的概念之一。一个算符 Â 作用在某个函数 φ 上，如果满足：
 φ = a φ那么 φ 就是  的本征态，a 是对应的本征值。本征态的特殊之处在于：在该状态下测量对应物理量，结果一定是确定值 a，而不是概率分布。对哈密顿算符而言，本征态就是能量确定的状态，本征值就是该状态的能量。
在材料计算中，我们通常把电子波函数展开成一组基函数的线性组合，把薛定谔方程转化为矩阵本征值问题：
H c = E S c其中 H 是哈密顿矩阵，S 是重叠矩阵，c 是展开系数向量。求解这个广义本征值问题，就能得到一系列本征能量 E 和对应的本征波函数。实际计算中，平面波基组、原子轨道基组或数值原子轨道基组的选择，会直接影响矩阵维数、计算成本和精度。例如，一个中等规模的半导体缺陷计算，平面波截断能常取 400–600 eV，k 点网格常取 4×4×4 或更密，这些参数本质上都是在控制本征值问题的求解精度。
材料计算依赖的基本假设
把量子力学用于真实材料时，必须明确几条基本假设，否则计算模型无法建立：
- 玻恩–奥本海默近似：原子核质量远大于电子，电子运动速度远快于原子核。因此可以把原子核视为固定不动，先求解电子在固定核构型下的电子结构，再让原子核在电子产生的势能面上运动。这是几乎所有第一性原理材料计算的第一步近似。
- 单电子近似：多电子问题无法精确求解，于是把每个电子看作在原子核和其他电子平均势场中独立运动，用单电子波函数描述。密度泛函理论中的 Kohn–Sham 方程就是这一思想的典型实现。
- 绝热近似与基态假设：多数材料计算关注基态电子结构，即能量最低的电子本征态。激发态、有限温度效应和强关联效应需要额外方法处理。
- 周期性边界条件：对晶体材料，利用平移对称性把无限大晶体约化为一个原胞加布洛赫定理，波函数写成 Ψ(r) = e^{ik·r} u(r)，其中 u(r) 具有晶格周期性。这样只需在倒空间第一布里渊区内采样 k 点，就能代表整个晶体。
把上述关系串起来看：波函数提供概率幅，薛定谔方程约束波函数的可能形式，算符把能量、动量等物理量转成数学作用，本征态给出确定测量值；而玻恩–奥本海默近似、单电子近似、基态假设和周期性边界条件，则把抽象的量子力学方程转化为材料计算中可求解的矩阵本征值问题。最终得到的电子密度、总能量和能带结构，才是我们用来判断材料稳定性、电子性质和力学响应的直接依据。