6.2 量子力学与化学

6.2 量子力学与化学

第二站进入微观世界。量子力学的中心方程——薛定谔方程——是一条偏微分方程：波函数 ψ(x, t) 是未知场，哈密顿算子替换了热传导里的拉普拉斯。表面上物理内涵天翻地覆，数学骨架却熟悉得令人安心：定态问题正是 3.1 的特征值三件套，能级就是特征值，波函数就是特征函数。本节从"波函数也是方程的解"讲起，用一维盒子清点能级并与数值对表，再沿"多体困难—密度泛函—反应扩散"的路线走到化学与生命现象的门口。

波函数也是方程的解

含时薛定谔方程写作

iħ·∂ψ/∂t = −(ħ²/2m)·∂²ψ/∂x² + V(x)·ψ。

与热传导方程并排一看：都是一阶时间、二阶空间，唯一的差别是热方程的时间项前多了个虚数单位 i。这个 i 是全部"量子怪象"的数学源头——热方程的解指数衰减（e^(−λt)），薛定谔方程的解永远振荡（e^(−iEt/ħ)）：耗散变成旋转，模长守恒，概率不灭。2.2 节性格档案里的"热核磨光"在量子世界不存在，波函数的棱角也会永远振荡下去——同一套二阶算子，因子 i 之差，两种宇宙。

定态解 ψ(x, t) = φ(x)·e^(−iEt/ħ) 代回方程，得到不依赖时间的定态薛定谔方程：

−(ħ²/2m)·φ'' + V(x)·φ = E·φ。

这就是哈密顿算子的特征值问题——与 3.1 的 X'' + λX = 0 逐字同构：边界条件选出离散的特征值（能级 E_n）与特征函数（波函数 φ_n）。分立能级不是量子力学的专利，而是"二阶算子加束缚边界"的数学共性，弦的泛音列与原子的电子壳层是同一件事的两种物理翻译。

数值实验：一维盒子的能级清点

无限深势阱（盒子内 V = 0、边界外 V 无穷，波函数边界取零）的解析能级人人会背：E_n = n²π²ħ²/(2mL²)。数值上把二阶导换成三点差分（4.2 的老朋友），哈密顿算子变成三对角矩阵，本征值交给标准算法：

import numpy as np

n, L = 200, 1.0; dx = L/(n+1)
hbar2_2m = 1.0                                  # 取原子单位且 L=1: E_n = n^2 pi^2 / 2
H = (-2*np.eye(n) + np.eye(n, k=1) + np.eye(n, k=-1)) * (-0.5/dx**2)

Eig, Vec = np.linalg.eigh(H)                     # 对称三对角，本征值升序
for k in [1, 2, 3, 4]:
    E_num = Eig[k-1]
    E_exact = k**2 * np.pi**2 / 2
    print(f"n={k}: 数值 {E_num:.4f}  解析 {E_exact:.4f}  误差 {abs(E_num-E_exact):.2e}")

psi2 = Vec[:, 0]**2                              # 基态概率密度
print("基态节点数（应为 0）:", np.sum(np.abs(psi2[1:-1]) < 1e-12))
​

四条能级的数值与解析在网格误差内重合；波函数按能量升序排列的节点数恰为 n−1（斯托姆–刘维尔定理的震荡性质）——特征函数的排序有严格理论背书，这个"节点计数"顺带成了检验特征求解器没算错的免费体检。矩阵尺寸 200 时毫秒级算完：量子能级清点在数值上就是一次对称矩阵本征分解，第四章的离散化加第七章的线性代数库，两件旧工具直接上岗。

从单粒子到化学：多体困难与绕行路线

真实化学的困难立刻降临：n 个电子的波函数是 3n 维空间里的函数，存储它需要指数资源——维数灾难是量子化学的底色。绕行路线有两条，都是 PDE 思想的延伸。密度泛函理论（DFT）证明基态能量由三维电子密度唯一决定，把 3n 维问题压回三维方程组（科恩– Sham 方程），代价是交换关联泛函只能近似——现代材料计算与分子模拟的中流砥柱，2025 年前后最先进的泛函仍在为化学精度打磨。变分蒙特卡洛与神经网络波函数：直接参数化波函数做变分优化，5.4 的 AI 思想在多体问题上的正名之战（神经网络在若干小分子上首次逼近实验精度）。

反应扩散：花纹的诞生

定态与能级之外，量子化学与生命科学共享另一类 PDE 场景：反应扩散方程——多种组分各自的扩散加彼此的非线性反应，u_t = D·Δu + f(u, v)。图灵 1952 年的惊人结论：扩散（通常抹平差异的力量）与反应（局部放大）耦合，可以把均匀状态变得不稳定，自发生成空间花纹——动物斑纹、贝壳螺旋、化学波现在都可定量模拟。1.3 的"和事佬"性格与 5.1 的非线性在一条方程里合流：扩散负责空间耦合，非线性负责失稳选择，花纹的波长由扩散系数比决定。

# 图灵花纹的最小种子：二维反应扩散（Gray-Scott 变体）确定性演化
# u_t = Du*lap(u) - u*v*v + F*(1-u);  v_t = Dv*lap(v) + u*v*v - (F+k)*v
# 这里只给框架：网格上反复执行 "扩散五点差分 + 局部反应" 两步
# 参数 F=0.035, k=0.065, Du=0.16, Dv=0.08 时出现稳定的斑点花纹
​

上述框架的完整版在四百行以内可以跑出花纹，配套检查与 4.4 一脉相承：反应项的实现必须守恒型（质量漂移会伪造花纹边界），时间步长受扩散显式格式约束（4.2 的 r 条件在二维翻倍）。有兴趣的读者把它当作本章的课后实验——参数空间的探索本身就是一次小型研究。

💡 关键直觉：量子力学给 PDE 教科书最深的回礼，是把"特征值"变成了可测量——能级、振动频率、化学键强度全是特征值的物理化名。反过来，数值特征值分析（第四章离散化加线性代数）也成了量子化学的标准器官。

本节要点

• 同构迁移：定态薛定谔方程即哈密顿算子的特征值问题，能级与泛音列是同一数学的两种翻译；
• 虚数之差：i 把耗散换成旋转，概率守恒取代热核磨光，性格档案在量子世界翻页；
• 能级清点：差分离散加对称本征分解，四条能级毫秒对表，节点计数免费体检；
• 多体困难：3n 维波函数不可存储，DFT 压回三维、神经波函数直接变分，两条绕行路线；
• 反应扩散：扩散加非线性失稳生成花纹，图灵机制的数值实现需守恒型反应项。

经典与微观之外，还有一个用方程定价的市场。下一站走进金融：从对冲组合推出 Black-Scholes，化归热方程拿闭式解，再让克兰克–尼科尔森在收益函数的棱角上接受考验。