3.2 薛定谔方程的求解方法
3.2 薛定谔方程的求解方法
本节摘要：定态薛定谔方程是量子力学计算的核心任务。本节系统介绍求解该方程的主要方法：分离变量法（适用于可分离势能）、级数展开法（适用于谐振子）、精确解析法（适用于氢原子）和微扰论（适用于近似求解）。每种方法的适用条件和物理直觉都将详细讨论。
定态薛定谔方程的求解策略
求解定态薛定谔方程的一般策略可以概括为三步。第一步，分析势能函数的对称性——如果势能在某些变换下不变，可以利用对应的守恒量简化方程。第二步，选择合适的坐标系和分离变量方法——利用对称性把偏微分方程分解为常微分方程。第三步，求解各常微分方程，施加边界条件，确定允许的能量本征值和对应的本征函数。
这个过程的核心思想是"利用对称性简化问题"。球对称势能允许在球坐标系中分离角度部分和径向部分；周期势能允许用布洛赫定理简化；柱对称势能允许在柱坐标系中分离。每利用一个对称性，就减少一个维度——这就是为什么氢原子（三维问题）可以简化为一维径向方程加上角动量量子数。
分离变量法：从三维到一维
分离变量法适用于势能可以写成各坐标独立函数之和的情况。最典型的例子是球对称势能：势能只依赖于到原点的距离，不依赖于角度。在这种情况下，薛定谔方程可以在球坐标系中分离为角度部分和径向部分。
角度部分的解是球谐函数——它们是角动量平方算符和角动量z分量的共同本征函数。球谐函数由两个量子数标记：角量子数和磁量子数。角量子数决定了球谐函数的总角动量，磁量子数决定了角动量在z方向的投影。球谐函数的正交归一性和完备性保证了任何球面上的函数都可以展开为球谐函数的线性组合。
径向部分的解取决于具体势能的形式。对于库仑势能（氢原子），径向方程可以精确求解，给出拉盖尔多项式形式的解。对于一般势能，径向方程通常需要数值方法或近似方法求解。无论势能的具体形式如何，径向波函数都必须满足两个边界条件：在原点处有限（对于束缚态），在无穷远处趋于零（保证归一化）。
级数展开法：量子谐振子的解
量子谐振子是量子力学中最重要的精确可解模型之一。它的势能是坐标的二次函数——一个抛物线势阱。谐振子的定态薛定谔方程可以通过级数展开法精确求解，结果是厄米多项式乘以高斯包络函数。
级数展开法的基本思路是：假设解可以展开为坐标的幂级数，代入薛定谔方程，得到递推关系。为了保证解的平方可积性，级数必须在有限项截断——这个截断条件直接给出了能量的量子化。级数截断的条件是能量必须取特定的离散值，相邻能级之间的间距为常数——这是谐振子等间距能级的数学根源。
谐振子的解有一些值得记住的重要性质。第一，基态波函数是高斯函数——这是不确定性原理最小化态的典型例子。第二，能级间距恰好等于普朗克常数乘以振子频率——与普朗克最初的量子假说一致。第三，选择定则只允许相邻能级之间的跃迁——这是电磁相互作用中光子角动量守恒的结果。
微扰论：近似求解的工具
大多数真实的量子系统无法精确求解。微扰论提供了一种系统化的近似方法：把哈密顿算符分解为一个可精确求解的部分（未扰哈密顿）加上一个小的修正项（微扰）。然后用泰勒展开的思想，把微扰的影响按阶数展开——零阶就是未扰解，一阶给出第一修正，二阶给出第二修正，以此类推。
非简并微扰论处理的是未扰能级不简并的情况。能量的一阶修正等于微扰算符在未扰态中的期望值，波函数的一阶修正涉及其他未扰态的贡献。简并微扰论处理的是未扰能级简并的情况——需要先在简并子空间中"对角化"微扰算符，才能确定能量的修正。
微扰论的一个经典应用是计算氢原子的精细结构——相对论修正和自旋-轨道耦合。这些修正比主能级间距小得多，所以可以用微扰论处理。另一个重要应用是斯塔克效应和塞曼效应——外电场和外磁场对原子能级的微扰。
图：微扰论的层次结构
变分法：能量的上限估计
变分法是另一种近似方法，特别适合估计基态能量。它的基本思想是：对于任何试探波函数，哈密顿算符在其上的期望值都不低于真实的基态能量。因此，如果你能猜到一个接近真实基态波函数的试探函数，你就能得到基态能量的一个上限——试探函数越好，上限越接近真实值。
变分法在量子化学中极其重要。分子的基态波函数通常无法精确求解，但可以通过选择合适的试探函数（如斯莱特行列式的线性组合）来获得相当精确的能量估计。现代量子化学计算（如Hartree-Fock方法和密度泛函理论）都是变分法的系统化应用。
一节小结
- 求解策略：利用对称性简化，选择合适坐标系，分离变量后求解常微分方程
- 分离变量法适用于球对称势能，角度部分给出球谐函数
- 级数展开法精确求解谐振子，截断条件给出能级量子化
- 微扰论适用于近似求解，分为非简并和简并两种情况
- 变分法给出基态能量的上限，是量子化学的基础工具
精确解法：氢原子
氢原子的薛定谔方程是量子力学中最重要的精确可解问题之一。它的势能是库仑势能，具有球对称性。在球坐标系中分离变量后，角度部分给出球谐函数，径向方程是一个具有库仑势的二阶常微分方程。
径向方程的精确解由三个量子数标记：主量子数决定能量，角量子数决定轨道角动量大小，磁量子数决定角动量z方向投影。能量只依赖于主量子数——这意味着不同角量子数的态可以是简并的（具有相同的能量）。这个简并性反映了库仑势能的特殊对称性（除了旋转对称性外还有一个额外的动力学对称性）。
径向波函数由拉盖尔多项式乘以指数衰减因子和幂函数构成。拉盖尔多项式的节点数等于主量子数减去角量子数减一——这就是为什么高角量子数的态在靠近原子核处有较小的概率密度。
氢原子解的一个重要推论是玻尔半径的概念：基态波函数的概率分布有一个特征长度尺度，叫做玻尔半径。玻尔半径大约是零点五三埃，它决定了原子的典型大小。这个长度尺度由电子质量、普朗克常数和库仑常数唯一确定——它是量子力学给出的第一个原子尺度的精确预言。
数值方法与计算量子力学
对于无法精确求解的系统，数值方法提供了有力的工具。常用的数值方法包括：有限差分法（把微分方程转化为差分方程）、有限元法（把空间离散化为网格）、变分蒙特卡洛方法（随机采样波函数空间）和密度泛函理论（用电子密度代替波函数作为基本变量）。
计算量子力学在现代科学和工程中有广泛应用。量子化学用它计算分子结构和反应能量，材料科学用它预测材料的电子结构和光学性质，凝聚态物理用它研究复杂多体系统。随着计算能力的提升，数值方法能处理的系统规模越来越大——从几个原子到几千个原子。
自洽场方法与多电子原子
多电子原子的薛定谔方程由于电子-电子相互作用的存在，通常无法精确求解。自洽场方法（如Hartree-Fock方法）提供了一种系统化的近似方案：假设每个电子在一个有效势场中运动，这个有效势场由原子核和其他电子的平均电荷分布决定。通过迭代求解，可以得到自洽的有效势场和对应的单电子轨道。
Hartree-Fock方法是量子化学的基础。它把多体问题近似为单体问题，但忽略了电子关联效应。更精确的方法（如组态相互作用、耦合簇理论）在Hartree-Fock基础上加入电子关联修正，可以给出更高的精度。密度泛函理论（DFT）则从另一个角度处理多体问题：用电子密度代替波函数作为基本变量，通过一个未知的泛函来表示总能量。Kohn-Sham方程是DFT的核心——它在形式上类似于单粒子薛定谔方程，但有效势能包含了所有多体效应。
在实际的量子力学课程中，求解薛定谔方程通常占据了大部分教学时间。这是因为虽然概念框架（公设、算符、测量）可以在几节课内讲完，但真正掌握量子力学需要大量的计算练习。每个具体系统的求解都有其独特的技巧和洞察——势阱的边界条件、谐振子的代数方法、氢原子的分离变量——这些技巧只能通过动手实践来掌握。
求解定态薛定谔方程是理解量子力学从抽象走向具体的关键一步。掌握了这些方法后，你将在下一章中看到它们如何被系统地应用于势阱、谐振子和氢原子这三个量子力学最重要的精确可解模型。
在结束本节之前，我想强调一个教学观点：量子力学的数学工具不是目的，而是手段。理解这些工具的目的只有一个——更有效地求解物理问题。如果你在理解某个数学概念时感到困难，回到物理问题来寻找直觉：这个数学工具是用来解决什么物理问题的？它是如何从物理问题中自然产生的？它在哪些简化假设下成立？
薛定谔方程求解中的每个技巧——分离变量、级数展开、微扰论、变分法——都是在特定物理条件下发展出来的。分离变量适用于有对称性的势能，级数展开适用于多项式势能，微扰论适用于弱扰，变分法适用于近似估计。当你遇到一个新的势能函数时，先分析它的对称性和渐近行为，然后选择最合适的求解方法。这种从物理到数学再回到物理的思维方式，是学习量子力学数学工具的最佳途径。