量子力学中的变分法

玻尔百科

核心要点

变分原理保证，由任何试探波函数计算出的能量都是量子体系真实基态能量的一个上限。
该方法涉及选择一个带有可调参数的灵活试探波函数，并通过最小化计算出的能量来找到最佳近似。
该原理是计算化学的基石，为诸如Hartree-Fock和组态相互作用（CI）等广泛使用的方法提供了理论基础。
Eckart不等式证实，随着试探能量接近真实基态能量，试探波函数本身也成为对真实波函数的更好近似。

引言

由薛定谔方程支配的量子世界提出了一个巨大的挑战：对于大多数我们感兴趣的体系，如原子和分子，找到其精确解是不可能的。这种复杂性在理论形式和实际预测之间造成了巨大的鸿沟。那么，科学家们如何确定分子的性质，以应用于从药物设计到材料科学等各个领域呢？答案在于有根据猜测的艺术，并通过强大的近似技术加以完善。在这些技术中，变分法以其独特的稳健性和直观性脱颖而出。

本文深入探讨变分法，这是现代量子力学和计算化学的基石。它提供了一种可靠的方法来估算体系的最低能量状态——基态。您不仅将学习该方法的工作原理，还将了解为什么它能保证“安全”，防止对真实能量的低估。在第一章中，我们将探讨使变分法如此有效的基本原理和数学机制。随后，我们将检验其关键应用，揭示这一优雅概念如何构成计算化学的基石，并为我们提供对物质结构的深刻见解。

原理与机制

想象你面临一个极其复杂的问题，一个谜题如此错综复杂，以至于即使是我们最强大的计算机也无法找到其精确解。你会怎么做？你可能会尝试做出一个有根据的猜测。但你如何确定你的猜测是好的呢？如果你有一条神奇的规则，一个黄金保证，告诉你任何合理的猜测都会为真实答案提供一个天花板，一个确定的上限，那会怎样？然后你就可以尝试许多不同的猜测，系统地完善它们，而给出最低天花板的那个就是你的最佳选择。在量子力学的世界里，这条神奇的规则是存在的，它被称为变分原理。它是物理学家工具库中最深刻、最实用的工具之一。

有根据猜测的艺术

变分原理的核心是找到量子体系的最低可能能量状态——其基态。除了最简单的体系外，支配量子粒子行为的薛定谔方程都极其难以精确求解。包含体系所有信息的波函数 $\psi$ 可能是一个极其复杂的函数，特别是对于像原子或分子中电子这样具有许多相互作用粒子的体系。

变分原理提供了一条优雅的出路。它指出，如果你取任何一个行为合理的试探波函数，我们称之为 $\psi_{trial}$ ，并计算能量的期望值，结果永远不会低于真实的基态能量 $E_0$ 。用数学术语表示：

E_{trial} = \frac{\langle \psi_{trial} | \hat{H} | \psi_{trial} \rangle}{\langle \psi_{trial} | \psi_{trial} \rangle} \ge E_0

这里， $\hat{H}$ 是哈密顿算符，即体系总能量的量子力学“配方”。这个不等式是一个安全网。它意味着你不可能“低估”真实的基态能量。你的猜测可能很差，给出的能量远高于 $E_0$ ，但它绝不可能通过低于 $E_0$ 而具有欺骗性的“好”。

让我们看看这在实践中意味着什么。考虑氦原子，它有一个原子核和两个活跃的电子。两个电子之间的排斥使得氦原子的薛定谔方程无法精确求解。假设我们尝试几种方法：

一个忽略电子排斥的幼稚模型给出的能量是 $-108.8$ eV。
一个使用称为一阶微扰理论的更复杂猜测给出的能量是 $-74.8$ eV。
一个使用带有灵活试探波函数的变分法的猜测给出的能量是 $-77.5$ eV。

实验测得的真实基态能量是 $-79.0$ eV。哪个结果我们可以信任？变分原理是我们的指南。它保证我们的变分结果 $E_{var}$ 必须是真实能量 $E_0$ 的一个上限。 $-77.5$ eV 是 $-79.0$ eV 的上限吗？是的！在数轴上， $-77.5$ 在 $-79.0$ 的右边，因此大于它。我们的变分计算完全符合这条基本规则。请注意，另一种常见的近似工具——微扰理论，并没有这样的保证；其结果可能高于或低于真实值。变分原理提供了一种单向的保证，这非常强大。

量子安全网：为什么猜测不会出错

为什么这个非凡的原理成立？原因深植于量子力学的结构之中。它源于一个关键的假设：一个体系哈密顿量的真实本征态——即具有确定能量的状态，如基态 ( $|\psi_0\rangle$ )、第一激发态 ( $|\psi_1\rangle$ ) 等——构成一个完备集。

把它想象成音乐。小提琴弦产生的纯音——基频及其泛音（谐波）——构成一个基。小提琴能发出的任何声音，无论多么复杂的刮擦声或尖锐声，都可以被描述为那些纯谐波音符的唯一组合，即一种叠加。

同样，你能想到的任何试探波函数 $|\psi_{trial}\rangle$ ，都可以表示为真实的、但未知的能量本征态的叠加：

|\psi_{trial}\rangle = c_0 |\psi_0\rangle + c_1 |\psi_1\rangle + c_2 |\psi_2\rangle + \dots

系数 $c_n$ 告诉我们我们的猜测中包含了“多少”每个真实的本征态。如果我们的猜测是归一化的（它应该是），那么这些系数的平方和为一： $|c_0|^2 + |c_1|^2 + |c_2|^2 + \dots = 1$ 。

当我们计算试探态的能量时，我们实际上计算的是其组成纯态能量的加权平均值：

E_{trial} = |c_0|^2 E_0 + |c_1|^2 E_1 + |c_2|^2 E_2 + \dots

根据定义，基态能量 $E_0$ 是最低的可能能量。所有其他能量都更高： $E_1 \ge E_0$ ， $E_2 \ge E_0$ ，依此类推。如果我们将上述总和中的每个 $E_n$ 都替换为可能的最低值 $E_0$ ，则总和只可能减少或保持不变：

E_{trial} = \sum_n |c_n|^2 E_n \ge \sum_n |c_n|^2 E_0 = E_0 \left(\sum_n |c_n|^2\right) = E_0(1) = E_0

就是这样。我们猜测态的计算能量是能量的加权平均值，所有这些能量都大于或等于 $E_0$ 。因此，平均值本身必须大于或等于 $E_0$ 。达到绝对最低能量 $E_0$ 的唯一方法是，如果我们的猜测从一开始就是完美的——即它不包含任何激发态的贡献，意味着 $c_1, c_2, \dots$ 都为零。在这种情况下， $|\psi_{trial}\rangle = |\psi_0\rangle$ ，我们的猜测就是真实的基态。这个优美的逻辑就是使变分法奏效的量子安全网。

方法实践：寻找最佳猜测

这个原理不仅仅是理论上的好奇心；它是一个寻找良好近似解的实用方法。策略不是进行疯狂的、随机的猜测，而是为试探波函数选择一个智能、灵活的形式，其中包含一个或多个变分参数。把这些参数想象成收音机上的调谐旋钮。我们转动旋钮，直到得到最清晰的信号——或者在我们的例子中，是最低的可能能量。

让我们来看一个例子。想象一个粒子处于一个相当奇怪的二维势场中，该势场依赖于离中心距离的对数， $V(r) = C \ln(r/a)$ 。我们很可能不知道确切的基态波函数。但是，我们可以做一个合理的猜测。由于势场是围绕原点对称的，基态可能也是对称的并且在中心处达到峰值。高斯函数 $\psi(r; \alpha) = N e^{-\alpha r^2}$ 似乎是一个合理的候选者。

这里， $\alpha$ 是我们的变分参数，我们的调谐旋钮。大的 $\alpha$ 对应一个尖锐、狭窄的波函数，而小的 $\alpha$ 对应一个宽阔、展开的波函数。对于每个 $\alpha$ 值，我们可以计算总能量 $E(\alpha)$ 。这个能量有两个部分：

动能，来自波函数的“摆动”。更窄的函数（更大的 $\alpha$ ）必须摆动得更多，所以动能随 $\alpha$ 的增加而增加。
势能，来自粒子在势阱中的平均位置。

总能量 $E(\alpha)$ 将是这两个相互竞争的项的总和。我们的任务是找到“最佳点”——使总能量最小化的 $\alpha$ 值。我们使用标准微积分来做到这一点：我们计算函数 $E(\alpha)$ ，然后找到其导数为零的点，即 $\frac{dE}{d\alpha} = 0$ 。满足此条件的 $\alpha$ 值给了我们最好的高斯近似，相应的能量 $E_{min}$ 是我们对基态能量的最佳估计。

有时，我们的猜测可能不仅仅是一个近似。对于简单的量子谐振子（一个处于抛物线势 $V(x) \propto x^2$ 中的粒子），真实的基态波函数就是一个高斯函数。如果我们使用高斯试探函数，变分法不仅仅给我们一个近似；它会引导我们找到确切的参数和确切的基态能量 $E_0 = \frac{1}{2}\hbar\omega$ 。这是一个美丽的证实：如果我们的试探函数族足够灵活，以至于包含了真实的解，那么变分原理保证我们会找到它。

了解局限：何时需谨慎使用该原理

像任何强大的工具一样，使用变分原理时必须了解其局限性。它对稳定体系非常有效，但如果其基本假设不被满足，则可能会产生误导。

如果一个体系的能级没有“底部”怎么办？物理学中描述的一些病态势允许能量骤降至负无穷大。对于这样的体系，哈密顿量没有下界。没有最低能量 $E_0$ 作为地板。变分原理的陈述 $E_{trial} \ge -\infty$ 在技术上是正确的，但对于找到一个有意义的能量值来说完全无用。幸运的是，描述大多数物理体系（如原子和分子）的哈密顿量都有下界，确保了它们的稳定性。

当一个体系没有“被捕获”或束缚态时，会出现另一个微妙之处。考虑一个处于纯排斥势中的粒子，如 $V(x) = C/x$ 且 $C > 0$ 。这个势将粒子推离原点。它不能将粒子保持在一个离散、局域的状态中。最低的可能能量是零，对应于一个静止在无穷远处、势消失的地方的粒子。所有其他状态都是具有正能量的“散射态”，代表一个穿过势场的粒子。如果在这里应用变分原理，它会正确地告诉我们任何试探态给出的能量都 $\ge 0$ 。它找到了能谱的底部，但没有离散的基态能级可以近似。该原理找到了最低可能能量的一个上限，而这个能量可能是一个连续态的起点，而不是一个单一、孤立的束缚态。

从原理到实践：现代计算化学

变分原理不仅仅是一个抽象概念；它是现代计算科学大部分内容赖以建立的基石。在量子化学中，科学家使用超级计算机来预测分子的性质，从染料的颜色到新药的疗效。这项工作的核心是求解分子电子的薛定谔方程，这项任务过于复杂，无法精确完成。

取而代之的是，化学家使用一系列近似方法，而变分原理是关键的组织概念：

Hartree-Fock (HF) 理论通常是起点。它用一个更简单的图像来近似电子复杂、相关的运动，即每个电子在所有其他电子的平均场中运动。波函数是一个单一的Slater行列式，变分原理被用来找到最好的行列式，保证Hartree-Fock能量 $E_{HF}$ 是真实能量的一个上限。
组态相互作用 (CI) 方法通过混合更复杂的电子组态（代表电子激发到更高轨道）来系统地改进Hartree-Fock方法。像CISD（包含单重和双重激发的CI）这样的方法也是变分的。因为它使用了比HF更灵活的试探波函数，其能量 $E_{CISD}$ 是一个更好（更低）的上限： $E_0 \le E_{CISD} \le E_{HF}$ 。
有趣的是，一些最准确和广泛使用的方法，如Møller-Plesset微扰理论 (MP2)和耦合簇 (CC) 理论，不是变分的。它们建立在不同的数学基础上，有时可能产生一个略低于真实值的能量。这是因为它们的数学结构虽然强大，但并不对应于计算哈密顿量的简单期望值。例如，CC理论涉及哈密顿量的非幺正变换，这导致了一个非厄米问题，超出了标准变分定理的范畴。这些方法牺牲了变分界限的绝对安全性，以换取其他理想的特性，如高精度和随分子大小的正确标度。

超越能量：猜测本身有多好？

变分原理是根据能量来陈述的。但我们通常真正关心的是波函数本身。如果我们找到一个试探波函数，它给出了一个非常好的能量估计，非常接近真实的 $E_0$ ，我们能确信我们的波函数也是真实基态波函数 $\psi_0$ 的一个好近似吗？

答案是肯定的，这要归功于另一个被称为Eckart不等式的优美结果。这个不等式在能量的准确性和波函数的质量之间提供了直接的联系。它给出了保真度 $F = |\langle \psi_0 | \psi_{trial} \rangle|^2$ 的一个下限，这个量衡量真实波函数和试探波函数之间的重叠或相似性（保真度为1意味着完美匹配）。这个不等式是：

F \ge \frac{E_1 - E_{trial}}{E_1 - E_0}

其中 $E_1$ 是第一激发态的能量。让我们来解读一下。分母 $E_1 - E_0$ 是基态和第一激发态之间的能隙，是体系的固定属性。分子 $E_1 - E_{trial}$ 追踪我们的试探能量与真实基态能量的接近程度。随着我们改进试探函数，使得 $E_{trial}$ 越来越接近 $E_0$ ，分子也越来越接近分母。这迫使保真度 $F$ 越来越接近1。

这是一个深刻的结论。寻求能量最小化的过程同时也是寻找最佳波函数的过程。变分原理不仅引导我们得到一个好的能量；它还含蓄地迫使我们的猜测越来越像真实的东西。它确保了在我们寻求最低能量的过程中，我们也在收敛到对量子态最真实的描述。这就是变分法的终极之美：它是一个简单、强大且稳健的原理，将有根据猜测的艺术转变为现代科学的基石。

应用与跨学科联系

既然我们已经掌握了变分原理的机制，我们可以退后一步，问一个最重要的问题：“这又如何？”这个原理有什么用？它仅仅是让我们解决一些人为设计的教科书问题，还是为我们打开了理解真实世界的大门？你会很高兴地发现，答案是，这个单一、优雅的思想是量子物理学家和化学家工具库中最强大、最通用的工具之一。它是现代计算科学大部分内容的根基。它是一座连接抽象量子形式主义与具体化学直觉的桥梁，也是我们从第一性原理出发预测物质行为的指路明灯。

让我们踏上一段旅程，看看这个原理是如何展开的，从最简单的“物理学家的游乐场”模型到计算化学的前沿及更远。

有根据猜测的艺术：简单的体系，深刻的见解

欣赏一个工具威力的最佳方式就是使用它。让我们从一个我们完全理解的体系开始：简谐振子。我们已经精确地解出了它，并且知道它的基态能量是 $E_0 = \frac{1}{2}\hbar\omega$ 。现在，让我们假装我们不知道答案。我们只有我们的物理直觉。基态波函数应该是什么样的？它应该以势阱底部（ $x=0$ ）为中心，并且随着我们向外移动而逐渐消失。一个对称的钟形曲线似乎是一个绝佳的猜测。高斯函数 $\psi(x) = A e^{-\alpha x^2}$ 完美地符合这个描述。

如果我们将这个猜测代入变分机制中，一个小小的奇迹发生了。在我们计算能量的期望值 $E(\alpha)$ 并相对于我们的变分参数 $\alpha$ 进行最小化之后，我们发现最小能量恰好是 $\frac{1}{2}\hbar\omega$ 。我们的猜测不仅仅是好；它是完美的！变分法不仅证实了我们的直觉，还通过找到 $\alpha$ 的最优值，给了我们确切的基态波函数。这是一个罕见但美丽的案例，我们选择猜测的函数族恰好包含了真实的答案。

但如果我们的直觉比较粗糙呢？如果我们不猜光滑的钟形曲线，而是猜一个简单的帐篷形状——一个在中心最大并线性减小到零的三角形波函数，会怎么样？这肯定不是真实的波函数；它在真实波函数光滑的地方有尖锐的扭结。然而，当我们在变分计算中运行它时，我们得到的基态能量估计值为 $\hbar\omega\sqrt{3/10} \approx 0.548\hbar\omega$ 。这惊人地接近真实值 $0.5\hbar\omega$ ！而且，正如原理所保证的，我们的结果是一个上限——我们计算出的能量确实高于真实能量。这就是该方法的真正魔力：它非常宽容。即使是一个粗略的猜测也能产生一个惊人准确的估计，而且我们还有一个内置的保证，即我们没有低估真实值。

这种稳健性使其成为解决现实世界问题的宝贵工具。考虑氢原子。虽然它可以精确求解，但它提供了一个至关重要的试验平台。如果我们使用高斯试探函数——一种数学上方便但物理上不正确的形式，因为真实的基态是指数形式的——变分原理仍然给出了一个能量上限，这个上限仅比真实值高约15%。该原理为我们提供了一个定量的衡量标准，衡量我们的猜测到底有多好（或多差），从而将直觉的艺术变成一门系统化的科学。

构建原子和分子：计算化学的基础

当我们进入多电子原子和分子的领域时，变分原理的真正威力就显现出来了，因为在这些领域，精确解是不可能的。在这里，问题不再仅仅是找到能量；而是关于理解化学结构和反应性。

想象一下，向氢原子核中加入第二个电子，形成负离子 H⁻。哈密顿量现在包含了一个描述两个电子之间排斥的项。这个项耦合了它们的运动，使得精确解变得棘手。我们如何猜测一个试探波函数？一个简单的想法是假设每个电子占据一个类氢的1s轨道。但是是哪一个呢？如果我们忽略另一个电子，每个电子都会看到核电荷 $Z=1$ 。但电子之间肯定会相互影响。在某种程度上，每个电子都会“屏蔽”另一个电子所感受到的原子核，部分抵消其正电荷。因此，每个电子感受到的有效核电荷 $Z_{eff}$ 应该小于 1。

这种物理直觉可以直接融入到变分计算中。我们使用一个由两个1s轨道构建的试探波函数，但我们将 $Z_{eff}$ 视为一个变分参数。然后，变分原理为我们完成工作：通过最小化总能量，它找到了电子-核吸引和电子-电子排斥之间的最佳平衡。结果证实了我们的直觉：当 $Z_{eff}$ 显著小于1时，能量最小化，因为这允许电子的波函数扩散开来，减少它们之间的相互排斥。变分原理刚刚为化学概念电子屏蔽提供了定量基础！

这个想法是量子化学中最重要的近似之一——Hartree-Fock (HF) 方法——的种子。Hartree-Fock方法不仅仅是像 $Z_{eff}$ 那样改变一两个参数，它是一个宏大规模的变分方法。它假设多电子波函数可以被一个单一的、正确反对称化的单电子轨道乘积（一个Slater行列式）来近似。然后，它提出问题：能够使总能量最小化的最佳可能的单电子轨道是什么？该方法改变这些轨道的整个函数形式，以找到提供最低可能能量的最优集合，但受限于一个约束条件，即电子表现为在所有其他电子的平均场中运动的独立粒子。

Hartree-Fock方法的核心是变分原理的直接而复杂的应用。它计算出的能量 $E_{\text{HF}}$ 是在平均场近似下体系可能达到的最佳能量。因为它是一种变分方法，它为真实的、精确的非相对论性能量 $E_{\text{exact}}$ 提供了一个严格的上限。这两者之间的差值， $E_{\text{corr}} = E_{\text{exact}} - E_{\text{HF}}$ ，就是化学家所说的相关能。它代表了由于电子运动相关联而带来的能量降低——即它们在瞬时基础上主动避开彼此的倾向，这一特征在平均场图像中被忽略了。由于变分原理保证 $E_{\text{HF}} \ge E_{\text{exact}}$ ，相关能总是一个负值，代表了Hartree-Fock近似中的“误差”，而更高级的方法正试图弥补这一误差。

这引导我们走向一个优美的量子化学方法“阶梯”。Hartree-Fock方法是第一级阶梯。为了得到更好的答案，我们需要一个更好的试探波函数。像组态相互作用 (CI) 这样的方法正是这样做的。它们将变分原理应用于一个更灵活的试探函数，该函数是Hartree-Fock行列式和代表激发电子组态的其他行列式的线性组合。CISD方法包括所有单重和双重激发，而Full CI包括给定基组下所有可能的激发。阶梯上的每一步都扩展了变分空间，从而能够更准确地描述电子相关。变分原理保证了能量的顺序： $E_{\text{HF}} \ge E_{\text{CISD}} \ge E_{\text{Full CI}} \ge E_{\text{exact}}$ 。这种系统性的改进是我们基本原理的直接结果。

新前沿：从波函数到密度和动力学

变分原理的影响并不仅限于基于波函数的方法。它还支撑着一种完全不同且极其成功的量子化学方法：密度泛函理论 (DFT)。构成DFT理论基础的Hohenberg-Kohn定理本身就是变分原理的一种形式。它们证明了基态能量可以通过最小化电子密度 $\rho(\mathbf{r})$ 的能量泛函来获得，而电子密度是一个比多体波函数简单得多的量。与波函数理论一样，DFT变分原理指出，对于任何试探密度，计算出的能量将大于或等于真实的基态能量。这为在传统波函数方法无法处理的极其庞大和复杂的体系上进行计算打开了大门。

此外，该原理不仅限于寻找静态基态。含时变分原理将这些思想扩展到动力学领域，通过在一组受限的波函数中找到使量子作用量泛函取极值的“路径”，来近似量子态随时间的演化。例如，这种强大的技术可以用来展示量子波包位置和动量的期望值如何根据类经典运动方程演化，从而在量子世界和经典世界之间建立了深刻的联系。

从最初作为估算简单体系能量的一种方式，变分原理已经发展成为一个统一的理论框架。它为我们的化学直觉提供了依据，为庞大的计算化学机器提供了引擎，也为我们观察量子世界动力学本身提供了一个镜头。它证明了在物理学中，有时最美丽、最强大的思想也是最简单的。游戏很简单：做出一个好的猜测。原理保证，你能把能量降得越低，你就越接近真理。