☰
定态薛定谔方程:从分离变量法到态叠加原理的完整解析
2026/10/3 18:34:54 网站建设 项目流程

量子力学学到"定态薛定谔方程"这一块,很多人会卡在一个微妙的地方:明明前面的章节都在讲波函数、概率幅、不确定性原理,怎么到这儿突然就开始解微分方程了?而且解出来的那些波函数既不带时间,又不体现运动,看起来跟"动力学"三个字完全不沾边。

我自己的教学经验里,十个学生至少有七八个在第一次接触定态方程时都会问同一个问题:那真实粒子的状态到底是不是定态?如果不是,我们把辛辛苦苦解出来的定态波函数拿来用,岂不是在自欺欺人?

这篇文章就是想把这几个问题讲透。我会从含时薛定谔方程出发,一步步走一遍分离变量法的完整思路,解释哈密顿量为什么是量子力学的核心对象,以及态叠加原理如何把这些看似"不动"的定态组合成真正随时间演化的物理过程。适合正在学量子力学、被定态概念绕晕的本科生,也适合想重新梳理一遍框架的自学者。

1. 分离变量法:猜一个"时间与空间各管各"的解

1.1 含时方程到底难在哪

不含时薛定谔方程的完整形式是:

[ i\hbar\frac{\partial \Psi(x,t)}{\partial t} = -\frac{\hbar^2}{2m}\frac{\partial^2 \Psi(x,t)}{\partial x^2} + V(x,t)\Psi(x,t) ]

稍微看一眼就知道,这个方程把时间的一阶导数、空间的二阶导数、以及可能同时依赖时间和空间的势能项全都搅在一起。如果势能 (V) 既随时间又随空间变化,这就是一个很难处理的偏微分方程,基本只能靠数值方法硬解。

但大部分我们在本科阶段遇到的物理系统,势能其实只跟位置有关,不随时间变化。比如无限深势阱的底部是平的,谐振子的势能是 (V(x) = \frac{1}{2}m\omega^2x^2),氢原子中的库仑势是 (V(r) = -\frac{e^2}{4\pi\varepsilon_0 r}),这些势能吃不懂时间,昨天长这样,今天还长这样。这给了我们一个极其重要的抓手。

我在实际教学中通常这样跟学生讲:当一个物理系统不随时间变化时,它的"行为模式"往往也不随时间变化,变的只是振幅或者相位。这就好比一根吉他弦,它的振动形状(驻波模式)是固定的,但拨弦之后弦振动的强弱会随时间衰减。空间上的"形状"和时间上的"演化"是两件可以分开讨论的事情。

1.2 乘积形式的假设

既然势能不依赖时间,就有理由猜测波函数可以写成一个只含空间坐标的函数和一个只含时间的函数的乘积:

[ \Psi(x,t) = \psi(x)\varphi(t) ]

这个写法叫分离变量法(separation of variables),本质上是把一个二元偏微分方程拆成两个一元常微分方程。注意这是"猜",不是推导出来的。之所以敢这么猜,是因为我们观察到系统本身具有"时间平移不变性"——时间的起点选在哪里,都不应该影响物理规律的形式。这种对称性反映在数学上,就是解可以分解成时间部分和空间部分的乘积。

把猜测的解代入含时薛定谔方程:

[ i\hbar\psi(x)\frac{d\varphi(t)}{dt} = -\frac{\hbar^2}{2m}\varphi(t)\frac{d^2\psi(x)}{dx^2} + V(x)\psi(x)\varphi(t) ]

接下来两边同时除以 (\psi(x)\varphi(t)),得到非常关键的形式:

[ i\hbar\frac{1}{\varphi(t)}\frac{d\varphi(t)}{dt} = -\frac{\hbar^2}{2m}\frac{1}{\psi(x)}\frac{d^2\psi(x)}{dx^2} + V(x) ]

等号左边是纯时间的函数,右边是纯空间的函数。左边不管空间坐标怎么变都不动,右边不管时间怎么变也不动,这两边还能永远相等,唯一的可能就是它们共同等于一个常数。这个常数记作 (E)。

1.3 两个常微分方程的拆分

于是原方程被拆成了两个独立的方程:

时间部分:

[ i\hbar\frac{d\varphi(t)}{dt} = E\varphi(t) ]

这是一个一阶线性常微分方程,解起来非常直接,通解是:

[ \varphi(t) = e^{-iEt/\hbar} ]

注意这里的 (E) 还没被赋予物理意义,它只是在分离变量过程中冒出来的一个待定常数。

空间部分:

[ -\frac{\hbar^2}{2m}\frac{d^2\psi(x)}{dx^2} + V(x)\psi(x) = E\psi(x) ]

这就是定态薛定谔方程(time-independent Schr\u00f6dinger equation)。到这里已经可以隐约看出,这个方程跟能量有关,因为两边都带着能量的量纲。但更严格的意义,需要通过算符形式来理解。

2. 哈密顿量:量子力学里"能量"的正式名分

2.1 从经典力学到量子力学的对应

经典力学里,一个不受摩擦、只受保守力作用的系统的总能量,是动能加势能:

[ H = \frac{p^2}{2m} + V(x) ]

这个 (H) 叫哈密顿量。在哈密顿力学框架中,系统的全部动力学信息都封装在 (H) 里面,给定初始位置和动量, (H) 决定系统接下来怎么演化。

量子力学建立时,有一个非常核心的规则叫"对应原理":把经典力学量换成算符,坐标乘上波函数,动量换成 (-i\hbar\frac{\partial}{\partial x}),其他的函数关系保持不变。按这个规则搬过来:

[ \hat{H} = \frac{\hat{p}^2}{2m} + V(\hat{x}) = -\frac{\hbar^2}{2m}\frac{d^2}{dx^2} + V(x) ]

这个 (\hat{H}) 就是哈密顿算符,它是动能算符和势能算符之和。这里的 (V(x)) 已经不是一个普通函数,而是"乘以势能函数"这个操作,在量子力学里它同样是一个算符。

所以定态薛定谔方程可以写成极其简洁的形式:

[ \hat{H}\psi(x) = E\psi(x) ]

这是我在教学中最强调的写法,因为它揭示了一个非常重要的本质:定态薛定谔方程是哈密顿算符的本征方程,能量 (E) 是本征值,波函数 (\psi(x)) 是本征态。量子力学里的能量不是随便什么数都能取的,必须是哈密顿算符的本征值,这就是能量量子化的数学来源。

2.2 算符本征方程的物理含义

看到 (\hat{H}\psi = E\psi) 这个形式,可以这样理解它的物理图景:如果粒子处于某个定态 (\psi(x)) 中,那么测量它的能量,测量结果一定是确定的本征值 (E),没有任何随机性。这一点跟位置和动量完全不同——位置和动量在一般状态下测量结果是不确定的。

放一个实验上的类比。在光谱学里,原子发出的谱线是分立的,这件事在量子力学里就是定态方程的必然结论:原子可以处在一些分立能量本征态上,电子在这两个态之间跃迁时,发出的光子能量就是两个能级之差 (E_j - E_i = \hbar\omega)。这个直接联系了定态方程和实验观测。

顺便提一个容易混淆的地方:定态方程里的 (E) 和含时解里那个 (e^{-iEt/\hbar}) 的 (E) 必须是同一个数。当初分离变量时假设两边共同等于常数 (E),现在空间部分解出来的本征值也同样要代回时间部分,否则得不到自洽的完整波函数。

2.3 哈密顿量为能量期望值的计算服务

除了本征方程,哈密顿算符还有另一个常见用途:计算任意状态下能量的平均值。对于任意一个已经归一化的波函数 (\Psi(x,t)),能量的期望值是:

[ \langle E \rangle = \int_{-\infty}^{\infty} \Psi^*(x,t)\hat{H}\Psi(x,t),dx ]

如果一个态不是能量本征态,这个平均值依然可以算,只是它不再等于一个确定的本征值,而是相当于对应能量可能取值的加权平均。这个细节在后面讨论态叠加原理时会起到核心作用。

3. 定态究竟定在哪里:定态波函数的三条重要性质

3.1 概率密度不随时间变化

把分离变量的解写完整:

[ \Psi(x,t) = \psi(x)e^{-iEt/\hbar} ]

现在计算概率密度:

[ |\Psi(x,t)|^2 = |\psi(x)|^2|e^{-iEt/\hbar}|^2 = |\psi(x)|^2 ]

这里用到了关键性质:(e^{-iEt/\hbar}) 的模长恒等于1,因为它是纯相位因子。所以概率密度在任意时刻都等于 (|\psi(x)|^2),完完全全不随时间改变。

这正是"定态"这个名字的由来:处于定态的粒子,它在空间各处的概率分布是固定的,不演化,不流动。这一点从直观上解释了为什么氢原子的电子轨道可以"稳定存在"——如果电子一直处在某个能量本征态上,它的电荷分布就是恒定的,不发散也不坍缩。

我给学生举过一个比较生活化的例子:把定态想象成一个驻波。吉他弦拨动之后,整根弦产生的波是往两头跑的,但如果我们只保留某一个振动模式,那根弦上每个点的振动幅度只在该点原地起伏,波峰波谷的位置固定不动。量子定态就类似一种"驻留"状态。

3.2 能量不确定性为零

如果系统处于 (\hat{H}) 的本征态,那么计算能量的标准差:

[ \Delta E = \sqrt{\langle E^2 \rangle - \langle E \rangle^2} = 0 ]

因为 (\hat{H}\psi = E\psi) 意味着 (\hat{H}^2\psi = E^2\psi),所以能量的期望值平方和能量的均方值都等于 (E^2),方差为零。

这是一个容易被忽略但非常重要的结论:定态不是"静止不动"的状态,而是能量完全确定的状态。定态里粒子仍然有动能、有动量分布,只是这些量的平均值不随时间变化,而且能量的测量结果没有任何不确定性。

3.3 所有力学量的平均值都不随时间变化

不只概率密度不随时间变,对于任何一个不显含时间的力学量算符 (\hat{A}),在定态中的期望值都是常数。证明很短:

[ \frac{d}{dt}\langle A \rangle = \frac{i}{\hbar}\langle [\hat{H},\hat{A}] \rangle = \frac{i}{\hbar}\left( E\langle A \rangle - E\langle A \rangle \right) = 0 ]

推导的每一步都用到了定态能量本征的性质。这说明定态是一个"完全静止"的统计意义上的状态:测位置的概率分布不变,测动量的概率分布不变,测任何可观测量的分布统统不变。这个情况好不好?好处是极其稳定,坏处是"无聊"——没有任何动态行为发生。

但现实世界显然不是这样的。电子会跃迁,原子会发光,粒子会运动。如果所有状态都是定态,宇宙就凝固了。这里就需要态叠加原理登场,它负责打破定态的"静止魔咒"。

4. 态叠加原理:把"静止"拼成"运动"

4.1 线性微分方程与解的叠加

定态薛定谔方程和含时薛定谔方程都是线性方程,意思是:如果 (\Psi_1) 和 (\Psi_2) 都是方程的解,那么它们的任意线性组合

[ \Psi(x,t) = c_1\Psi_1(x,t) + c_2\Psi_2(x,t) ]

仍然是方程的解,其中 (c_1, c_2) 是复常数,不需要满足任何额外条件。

这个数学性质背后是物理上的态叠加原理:如果系统可以处于状态 (\Psi_1),也可以处于状态 (\Psi_2),那么它可以处于两者的任意线性叠加态。叠加不是"混合物"的意思,不是说系统有一半概率处于状态1、一半概率处于状态2,而是说系统同时处在一种全新的、无法用经典概率描述的相干状态中。

一个经典的类比是:在双缝实验中,电子同时通过两条缝的状态,是"通过缝1的状态"和"通过缝2的状态"的量子叠加。这种叠加带来了干涉条纹——这是任何经典粒子都无法产生的现象。如果只是经典概率混合,得到的是两条单缝衍射图案的简单相加,不会出现干涉。

4.2 把几个定态叠加起来会发生什么

现在关键的一步来了。假设我们有两个不同的定态:

[ \Psi_1(x,t) = \psi_1(x)e^{-iE_1t/\hbar},\quad \Psi_2(x,t) = \psi_2(x)e^{-iE_2t/\hbar} ]

叠加后的波函数:

[ \Psi(x,t) = c_1\psi_1(x)e^{-iE_1t/\hbar} + c_2\psi_2(x)e^{-iE_2t/\hbar} ]

这个态的概率密度是:

[ |\Psi(x,t)|^2 = |c_1|^2|\psi_1(x)|^2 + |c_2|^2|\psi_2(x)|^2 + 2\text{Re}\left[c_1^c_2\psi_1^(x)\psi_2(x)e^{i(E_1-E_2)t/\hbar}\right] ]

等号右边前两项是各自的概率密度,跟时间无关;第三项是交叉项,其中包含一个随时间振荡的相位因子 (e^{i(E_1-E_2)t/\hbar})。只要 (E_1 \neq E_2),这个交叉项就是时间函数,概率密度就会随时间变化,系统的空间分布就会流动起来。

我把这件事概括为一句比较好记的话:单个定态静止不动,两个定态叠加就能"动"起来,运动的频率由两个定态之间的能量差决定,即 (\omega = (E_1 - E_2)/\hbar)。

这个关系实在重要。它告诉我们:量子力学中的所有"变化",本质上都源于不同能量本征态之间的干涉。原子发光是因为电子从高能量的定态跃迁到低能量的定态,光子的频率对应着两个态的能量差。这也解释了为什么原子光谱是分立的线状谱——因为可供跃迁的定态能级是分立的,能级差也就只取某些确定的值。

4.3 测量坍缩与叠加的关系

态叠加原理还牵出一个常被误解的问题:测量前系统处于叠加态,测量后却只得到一个确定的本征值,这个过程叫波函数坍缩。比如一个电子可以处于两个能级的叠加态,但测量它的能量,结果非此即彼,要么 (E_1) 要么 (E_2),每次测量都会"随机"给出一个结果,得到 (E_1) 的概率是 (|c_1|^2),得到 (E_2) 的概率是 (|c_2|^2)。

这里有个微妙的点:测量给出 (E_1) 之后,系统的波函数就从原来的叠加态坍缩到了 (\Psi_1),下一次测量立刻再做,结果还是 (E_1)。如果不去测量,系统会继续处于叠加态,其时间演化和干涉现象完全按前面说的公式进行。测量不是"看看系统本来是什么",而是强行改变了系统的状态。

我自己在讲解这部分时喜欢强调:叠加态本身的演化是完全确定性的,只有测量这个环节引入随机性。在量子力学中,"测量"不是一个被动观察过程,而是一个主动干预过程。这个认识对理解后续很多量子信息的内容至关重要。

5. 从一维到真实系统:分离变量法如何应对不同势能

5.1 三类最典型的势能模型

分离变量法和定态方程在具体应用中最常见的几个模型,几乎构成了量子力学的"标准课表"。

无限深方势阱。势能为零的区域里,定态方程简化为:

[ -\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2} = E\psi ]

这是一个二阶常系数线性微分方程,通解是正弦和余弦的组合。加上两端波函数必须为零的边界条件,就得到量子化的能级:

[ E_n = \frac{n^2\pi^2\hbar^2}{2mL^2},\quad n=1,2,3,\dots ]

以及本征态:

[ \psi_n(x) = \sqrt{\frac{2}{L}}\sin\left(\frac{n\pi x}{L}\right) ]

这里最让学生意外的结论是:能量不能取任意值,只能取分立的一组值,最低的 (n=1) 态能量也不是零,而是 (\frac{\pi^2\hbar^2}{2mL^2}),这称为零点能。经典粒子在势阱底部可以静止,量子粒子永远有最低的动能,因为如果动能真为零,波函数就处处为零,粒子就"不存在"了。

谐振子。势能是二次函数:

[ V(x) = \frac{1}{2}m\omega^2x^2 ]

这个模型描述的是平衡点附近的微小振动,不管是分子振动、晶格振动还是量子场论里的各种激发,本质上都归结为谐振子。它的能级是等间距的:

[ E_n = \left(n+\frac{1}{2}\right)\hbar\omega,\quad n=0,1,2,\dots ]

零点能是 (\frac{1}{2}\hbar\omega),这个结果在量子场论里演化出了真空涨落的概念。谐振子的求解常用幂级数方法或者升降算符方法,后者在量子力学学习中是一个重要的里程碑。

氢原子。库仑势是球对称的:

[ V(r) = -\frac{e^2}{4\pi\varepsilon_0 r} ]

处理氢原子需要在球坐标系下做分离变量,把径向部分和角向部分分开。角向部分对应于角动量算符的本征函数,也就是球谐函数 (Y_l^m(\theta,\phi))。径向方程解出来的能级是:

[ E_n = -\frac{13.6\text{ eV}}{n^2},\quad n=1,2,3,\dots ]

这组数值精确地解释了氢原子光谱线系。值得一提的是,分离变量法的使用在这里起到了决定性作用——先分离时间,再分离角度,再分离径向,层层递进,每一步都用了系统的对称性。

5.2 一维势垒与隧穿效应

除了束缚态,定态方程还适用于散射问题。一维势垒问题中,粒子的能量可能低于势垒高度,但定态方程的解仍然存在。求解之后,反射系数和透射系数不是零,粒子有一定概率"穿过"势垒,这就是隧穿效应。

扫描隧道显微镜(STM)其实就是利用隧穿电流来成像的,探针和样品之间隔着真空,电子靠量子隧穿"漏"过去。没有定态方程和相应的求解方法,根本没法解释这种经典物理完全不允许的现象。

所以学定态方程不是纯粹学数学技巧,它是用来解释真实实验的:光谱线的分立结构、零点能、隧穿效应、势垒反射,这些都是定态方程直接给出的预言。

6. 处理定态方程时最常见的几个坑

6.1 混淆"定态"与"基态"

"定态"是能量本征态,任何 (n=1,2,3,\dots) 的状态都叫定态;"基态"是能量最低的那个态,也就是 (n=1) 的那个态。两者是包含关系,不是等同关系。我见过不少学生在答题时把"处于定态"直接理解成"处于基态",然后漏掉了所有激发态的讨论,这个坑一定要避开。

6.2 忘记归一化

解出定态方程的本征函数之后,需要做归一化处理:通过调节常数,让波函数的模方在全空间积分为 1。这个步骤看起来简单,但实际上如果漏掉,后面的概率计算全部都会出错。比如我之前列出的无限深势阱本征态里的系数 (\sqrt{2/L}),就是归一化得到的。

归一化还有一个更深的含义:定态波函数平方在全空间的积分有限,代表粒子被束缚在势阱内。如果粒子能量高于势垒,波函数是散射态,归一化要用另一种方式(比如按 δ 函数归一化),这是很多初学introduction to QM时会忽略的地方。

6.3 把哈密顿量的符号搞混

哈密顿量算符是 (\hat{H} = -\frac{\hbar^2}{2m}\frac{d^2}{dx^2} + V(x)),其中动能部分是"负号加上二阶导数"。如果推导过程中粗心,很容易把正负号写反,或者把二阶导写成一阶导。一个比较好的自查方法就是想一下自由粒子的情形: (V=0) 时定态方程变成

[ -\frac{\hbar^2}{2m}\psi'' = E\psi ]

它的解是平面波 (e^{\pm ikx}),能量是正的。如果符号写反,能量就变成负数,物理上就不合理了。

6.4 边界条件的自觉性

定态方程的边界条件跟势能函数紧密相关。比如无限深势阱中,波函数在边界处必须为零,因为势能无限大,粒子不可能出现在边界外。有限深势阱中,波函数不会在边界处突变为零,而是会指数衰减地"漏"进势垒区域。很多人只记住"边界处为零"这一种情况,遇到有限深势阱就胡乱套用,结果整个求解逻辑全错了。

一个更隐蔽的坑是边界处波函数导数的连续性条件。势能函数有限跳跃的情况下,波函数和它的一阶导数在边界处要求连续;但如果是 δ 势这种极端情况,一阶导数会有跳跃,连续条件要修正。这些细节我会建议学生做题时先把势能画出来,再逐段分析和匹配,比直接背结论可靠得多。

6.5 叠加系数要满足归一化

叠加态 (\Psi = \sum_n c_n \psi_n e^{-iE_nt/\hbar}) 的归一化条件是:

[ \sum_n |c_n|^2 = 1 ]

这个条件同时保证了各本征态出现的总概率为 1。初学者常常只归一化单个定态波函数,或者忘记检查叠加系数,结果导致概率计算不自洽。实际上,由于定态波函数之间的正交性,这个条件可以在计算内积时自动处理。

7. 把三块拼图拼成量子力学全貌

到了这一步,整个量子力学的理论框架其实已经露出轮廓:

含时薛定谔方程给出一般的演化规则——知道初始波函数,就可以算出之后任意时刻的状态。但直接求解太难,所以我们利用势能不随时间变化这个特点,用分离变量法猜出乘积形式的解,得到一个定态薛定谔方程。定态方程是一个本征值问题,哈密顿算符作用在波函数上等于能量乘以波函数,这直接导出了能量量子化以及相对应的空间波函数。

单个定态是"静止"的,但态的叠加原理给出了最普适的描述——任意一个物理上合理的波函数,都可以写成所有定态的线性组合。在叠加态里,不同能量本征态之间的干涉使概率密度随时间演化,变化频率由能级差决定。测量能量时,系统以概率 (|c_n|^2) 坍缩到对应的本征态上。

这个结构之所以优雅,在于它把三个在形式上非常不同的概念串联成了一个逻辑闭环:

  • 哈密顿量决定定态的存在与能谱结构
  • 分离变量法提供从含时方程到定态方程的途径
  • 态叠加原理保证了量子系统能"活"起来,用定态之间的干涉来描述真实的动态行为

我还想提醒一点:这三个概念不是孤立的招式,而是一套完整的认识论工具。比如量子隧穿、能级跃迁、量子纠缠,以及以后在量子计算中看到的叠加态,本质上都在反复使用"本征态 + 叠加系数 + 测量概率"这个框架。理解了定态方程,等于拿到了一把能打开量子力学后续几乎所有内容的钥匙。

在实际学习的路上,我常跟学生说的一句话是:物理概念最怕的不是难,而是"以为自己懂了",结果做题时发现边界条件不会写、归一化系数没算。定态方程这块尤其如此,因为它的内容横跨了"物理图像"和"数学计算"两个维度。建议学完本章后,亲自把无限深势阱、谐振子、氢原子三个系统的定态方程从头到尾各推一遍,推到不需要看参考书为止。推完再回头看这篇文章讲的每一个概念环节,定会有完全不一样的感觉。

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询