第1章 波函数:粒子在观测之前在哪里?

📑 目录

"没有人真正理解量子力学。" —— 理查德·费曼

Griffiths的回应是:没关系,我们先学会计算。

这一章,我们要建立一个看似荒谬但实验证明为真的世界观:粒子在观测前,不处于确定的位置。


第1章 波函数:粒子在观测之前在哪里?

0. 前置知识:复数基础

波函数是一个复数值函数。在深入量子力学之前,我们需要确保复数运算已经熟练。

0.1 复数的基本形式

一个复数 zz 可以写成:

z=a+biz = a + bi

其中 a,ba, b 是实数,i=1i = \sqrt{-1} 是虚数单位。

模(绝对值)

z=a2+b2|z| = \sqrt{a^2 + b^2}

复共轭

z=abiz^* = a - bi

注意一个重要关系:z2=zz=a2+b2|z|^2 = z^* z = a^2 + b^2

0.2 欧拉公式

这是量子力学中最重要的公式之一:

eiθ=cosθ+isinθe^{i\theta} = \cos\theta + i\sin\theta

为什么重要? 平面波 eikxe^{ikx} 是量子力学的"基本构件",而它正是复数指数形式。任何正弦/余弦波都可以写成复指数的实部或虚部。

常见形式

  • eiπ=1e^{i\pi} = -1(欧拉恒等式,最美的数学公式之一)
  • eiπ/2=ie^{i\pi/2} = i
  • e0=1e^{0} = 1

0.3 复数的极坐标表示

z=zeiϕz = |z| e^{i\phi}

其中 z|z| 是模,ϕ\phi相位。两个复数相乘,模相乘,相位相加:

z1z2=z1z2ei(ϕ1+ϕ2)z_1 z_2 = |z_1||z_2| e^{i(\phi_1 + \phi_2)}

物理直觉:量子力学中,概率是模的平方 Ψ2|\Psi|^2,相位本身不出现在概率中。但相位在叠加干涉中至关重要——两个波函数的相位差决定了它们是相长干涉还是相消干涉。这就像声波:两个同相声波叠加会变大,两个反相声波叠加会抵消。

0.4 数值例题

例题1:化简 eiπ/4e^{i\pi/4},并求其模。

解答

eiπ/4=cos(π/4)+isin(π/4)=22+i22e^{i\pi/4} = \cos(\pi/4) + i\sin(\pi/4) = \frac{\sqrt{2}}{2} + i\frac{\sqrt{2}}{2}

模:eiπ/4=(2/2)2+(2/2)2=1/2+1/2=1|e^{i\pi/4}| = \sqrt{(\sqrt{2}/2)^2 + (\sqrt{2}/2)^2} = \sqrt{1/2 + 1/2} = 1

例题2:计算 (1+i)2(1+i)^2 的模。

解答

(1+i)2=1+2i+i2=1+2i1=2i(1+i)^2 = 1 + 2i + i^2 = 1 + 2i - 1 = 2i

模:2i=2|2i| = 2

或者利用 z2=z2=(12+12)2=2|z^2| = |z|^2 = (\sqrt{1^2+1^2})^2 = 2


1.1 故事:凌晨2点的量子实验室

凌晨2点,量子科技公司的实验室里,实习生小林盯着屏幕上那条诡异的概率分布曲线。

他们研发的"量子骰子机"出了故障。按照经典物理,骰子落地后应该有一个确定的面朝上。但机器的数据显示:在打开观测罩之前,骰子似乎同时处于1点到6点的所有状态——不是"不知道是哪个",而是真的同时存在,像一个模糊的云团。

"这不可能。"小林揉着眼睛,"骰子明明已经停下来了,怎么可能同时是1又是6?"

导师老张走过来,看了一眼数据:"这不是bug,这就是量子力学。在观测之前,粒子的状态由波函数描述——它不是一个确定的值,而是所有可能值的叠加,每种可能性带一个复数权重。"

"那打开观测罩之后呢?"

"波函数坍缩成一个确定的结果。但在那之前,它确实是所有可能性的叠加。"

小林沉默了很久,终于问出那个经典问题:"那在打开盖子之前,骰子到底是什么?"

老张笑了:"费曼说没有人真正理解量子力学。你现在知道为什么了。"

核心洞察:量子世界不是"我们不知道粒子在哪",而是"粒子在观测前确实不处于确定的位置"。


1.2 薛定谔方程:波函数的演化法则

如果粒子在观测前不是一个确定的位置,那我们如何描述它?

答案是:波函数 Ψ(x,t)\Psi(x, t)。它是一个复数值函数,给每个位置 xx 和时间 tt 分配一个复数振幅。

1.2.1 时间依赖薛定谔方程

波函数怎么随时间变化?由薛定谔方程决定:

iΨt=22m2Ψx2+VΨi\hbar \frac{\partial \Psi}{\partial t} = -\frac{\hbar^2}{2m} \frac{\partial^2 \Psi}{\partial x^2} + V\Psi

这个方程的地位,相当于经典力学中的 F=maF=ma。但有几个关键区别:

经典力学量子力学
求解 x(t)x(t)——粒子在某时刻的位置求解 Ψ(x,t)\Psi(x,t)——粒子在某时刻的"概率幅分布"
给定初始位置 x0x_0 和速度 v0v_0给定初始波函数 Ψ(x,0)\Psi(x,0)
结果是确定性的结果是概率性的
graph LR
    subgraph 经典力学
        A["初始条件
x₀, v₀"] -->|"F=ma"| B["确定轨迹
x(t)"] end subgraph 量子力学 C["初始波函数
Ψ(x,0)"] -->|薛定谔方程| D["波函数演化
Ψ(x,t)"] D -->|测量| E[概率性结果] end style B fill:#e8f5e9 style E fill:#fff3e0

历史背景:薛定谔的灵感

1925年的圣诞节,埃尔温·薛定谔带着他的情人去瑞士阿尔卑斯山度假。他在那里研究了德布罗意的物质波假说,突然产生了一个想法:如果电子真的是波,那它的波动方程应该是什么?

几周后,他发表了那篇改变物理学的论文,提出了以他名字命名的方程。有趣的是,薛定谔最初以为自己找到了"经典波"的方程,试图用它来解释原子结构——他花了好几个月才接受玻恩提出的概率诠释。

graph TD
    subgraph 薛定谔方程的诞生
        A[德布罗意物质波] --> B[薛定谔的疑问]
        B --> C["如果粒子是波
它的方程是什么?"] C --> D[1926年薛定谔方程] D --> E["薛定谔:这是经典波"] E --> F["玻恩:不,这是概率幅"] F --> G["爱因斯坦:上帝不掷骰子"] G --> H["玻尔:别告诉上帝怎么做"] end style D fill:#e8f5e9 style F fill:#fff3e0 style H fill:#ffebee

1.2.2 波函数必须满足的条件

不是所有函数都能当波函数。一个合法的波函数必须满足:

  1. 归一化:总概率为1
  2. 连续(势能有限时)
  3. 一阶导数连续(势能有限时)
  4. 平方可积Ψ2dx\int_{-\infty}^{\infty} |\Psi|^2 dx 有限

物理直觉:波函数在边界处必须"平滑连接"。如果波函数本身不连续,意味着粒子在某处有无限大的概率密度——这在物理上不可接受(除非势能无限大,如无限深势阱壁)。如果导数不连续,则意味着二阶导数存在δ函数,对应势能中有δ函数项。


1.3 玻恩诠释:波函数到底什么意思?

这是量子力学最核心的概念转变。

1.3.1 概率密度

Ψ(x,t)2|\Psi(x,t)|^2 不是粒子的"密度",而是概率密度——在位置 xx 附近找到粒子的概率。

具体来说,在区间 [a,b][a, b] 内找到粒子的概率是:

P(x[a,b])=abΨ(x,t)2dxP(x \in [a,b]) = \int_a^b |\Psi(x,t)|^2 dx

graph TD
    subgraph 波函数的概率诠释
        A["波函数 Ψ(x)"] -->|取模平方| B["概率密度 |Ψ(x)|²"]
        B -->|积分| C[在某区域找到粒子的概率]
        B -->|全空间积分| D["总概率 = 1"]
    end
    
    style D fill:#e8f5e9

例题详解:归一化常数的确定

问题:波函数 Ψ(x)=Aeax2\Psi(x) = A e^{-a x^2}a>0a > 0)是否可归一化?归一化常数 AA 是多少?

分步推导

第一步:写出归一化条件

Ψ2dx=1\int_{-\infty}^{\infty} |\Psi|^2 dx = 1

第二步:代入波函数

A2e2ax2dx=1|A|^2 \int_{-\infty}^{\infty} e^{-2a x^2} dx = 1

第三步:利用高斯积分公式

eαx2dx=πα\int_{-\infty}^{\infty} e^{-\alpha x^2} dx = \sqrt{\frac{\pi}{\alpha}}

这里 α=2a\alpha = 2a,所以:

A2π2a=1|A|^2 \sqrt{\frac{\pi}{2a}} = 1

第四步:解出 AA

A=(2aπ)1/4A = \left(\frac{2a}{\pi}\right)^{1/4}

(取实数正根,整体相位不影响物理)

物理直觉:高斯波包的宽度由 aa 控制——aa 越大,波包越窄,粒子位置越确定。但由不确定性原理,位置越确定,动量就越不确定。

graph LR
    subgraph 高斯波包的参数
        A["参数 a 大"] --> B[波包窄]
        B --> C[位置确定]
        C --> D[动量不确定]
        
        E["参数 a 小"] --> F[波包宽]
        F --> G[位置不确定]
        G --> H[动量确定]
    end
    
    style C fill:#ffebee
    style G fill:#e3f2fd

数值例题:具体波函数的归一化与概率

问题:波函数 Ψ(x)=Ae2x2\Psi(x) = A e^{-2x^2}xx 单位为 m),求归一化常数 AA,并计算在区间 [0.5,+0.5][-0.5, +0.5] m 内找到粒子的概率。

解答
第一步:写出归一化条件

A2e4x2dx=1\int_{-\infty}^{\infty} |A|^2 e^{-4x^2} dx = 1

第二步:利用高斯积分

e4x2dx=π4=π2\int_{-\infty}^{\infty} e^{-4x^2} dx = \sqrt{\frac{\pi}{4}} = \frac{\sqrt{\pi}}{2}

第三步:解出 AA

A2π2=1A=(2π)1/21.06 m1/2|A|^2 \cdot \frac{\sqrt{\pi}}{2} = 1 \Rightarrow A = \left(\frac{2}{\sqrt{\pi}}\right)^{1/2} \approx 1.06 \text{ m}^{-1/2}

第四步:计算区间概率

P=0.50.5Ψ2dx=A20.50.5e4x2dxP = \int_{-0.5}^{0.5} |\Psi|^2 dx = |A|^2 \int_{-0.5}^{0.5} e^{-4x^2} dx

换元 u=2xu = 2xdu=2dxdu = 2dx

P=A2211eu2du=1ππ2erf(1)=12erf(1)P = \frac{|A|^2}{2} \int_{-1}^{1} e^{-u^2} du = \frac{1}{\sqrt{\pi}} \cdot \frac{\sqrt{\pi}}{2} \text{erf}(1) = \frac{1}{2}\text{erf}(1)

查表得 erf(1)0.843\text{erf}(1) \approx 0.843,所以 P0.421P \approx 0.421

结论:约有42.1%的概率在 [0.5,+0.5][-0.5, +0.5] m 区间内找到粒子。


1.3.2 归一化

粒子一定在某个地方,所以:

Ψ(x,t)2dx=1\int_{-\infty}^{\infty} |\Psi(x,t)|^2 dx = 1

一个好消息:如果波函数初始时归一化,薛定谔方程会保证它永远归一化(概率守恒)。

证明思路:对薛定谔方程取复共轭,与原方程相减并积分,可以证明:

ddtΨ2dx=0\frac{d}{dt} \int_{-\infty}^{\infty} |\Psi|^2 dx = 0

这意味着概率"不会泄漏"——粒子既不会凭空消失,也不会凭空产生。薛定谔方程自动守恒概率,这正是我们想要的。

graph TD
    subgraph 概率守恒
        A[初始归一化] --> B[薛定谔方程演化]
        B --> C[始终归一化]
        D[概率守恒] --> E[粒子不会消失]
        D --> F[粒子不会凭空产生]
    end
    
    style C fill:#e8f5e9
    style E fill:#e8f5e9
    style F fill:#e8f5e9

统计诠释的完整公理

Griffiths 将量子力学的统计诠释归结为以下几条核心规则:

公理1(概率密度)Ψ(x,t)2|\Psi(x,t)|^2 是在位置 xx 处、时刻 tt 找到粒子的概率密度

公理2(区域概率):在区间 [a,b][a, b] 内找到粒子的概率为:

P(x[a,b])=abΨ(x,t)2dxP(x \in [a,b]) = \int_a^b |\Psi(x,t)|^2 dx

公理3(归一化):总概率为1:

Ψ(x,t)2dx=1\int_{-\infty}^{\infty} |\Psi(x,t)|^2 dx = 1

公理4(期望值):可观测量 Q(x,p)Q(x,p) 的期望值为:

Q=ΨQ^Ψdx\langle Q \rangle = \int_{-\infty}^{\infty} \Psi^* \hat{Q} \Psi \, dx

公理5(测量坍缩):测量后,波函数坍缩到与测量结果对应的本征态上。

这五条公理构成了第1章的完整逻辑框架。第1章的其余内容,本质上都是在展开、解释和应用这五条公理。


1.3.3 故事续集:骰子的真相

回到小林的实验室。

他们终于理解了:量子骰子机里不是一个"不知道具体数值的骰子",而是一个由波函数描述的量子系统。在观测前,波函数可能是这样的:

Ψ=16(ψ1+ψ2+ψ3+ψ4+ψ5+ψ6)\Psi = \frac{1}{\sqrt{6}}(\psi_1 + \psi_2 + \psi_3 + \psi_4 + \psi_5 + \psi_6)

每个分量 ψn\psi_n 对应"点数为 nn"的状态。系数 16\frac{1}{\sqrt{6}} 保证归一化,也让每个结果的概率都是 1/62=1/6|1/\sqrt{6}|^2 = 1/6

"所以"小林恍然大悟,"骰子不是’还没确定’,而是’同时是1到6’,观测后才变成确定的某一个?"

"对。"老张说,"而且更诡异的是:如果你把波函数改成

Ψ=12ψ1+32ψ6\Psi = \frac{1}{2}\psi_1 + \frac{\sqrt{3}}{2}\psi_6

那测出1的概率是 1/41/4,测出6的概率是 3/43/4。不是50/50,因为振幅的平方才是概率。"

关键:概率是振幅的模平方,不是振幅本身。

例题详解:概率振幅的干涉

问题:如果 Ψ=12(ψ1+ψ2)\Psi = \frac{1}{\sqrt{2}}(\psi_1 + \psi_2),和 Ψ=12(ψ1ψ2)\Psi' = \frac{1}{\sqrt{2}}(\psi_1 - \psi_2),测量得到状态1的概率分别是多少?

解答

对于 Ψ\Psi

P1=122=12P_1 = \left|\frac{1}{\sqrt{2}}\right|^2 = \frac{1}{2}

对于 Ψ\Psi'

P1=122=12P_1 = \left|\frac{1}{\sqrt{2}}\right|^2 = \frac{1}{2}

等等——看起来概率一样?但如果我们问的是"在位置空间测量",而 ψ1\psi_1ψ2\psi_2 是位置波函数,那总波函数在某些位置可能出现相消干涉!

Ψ2=12ψ1+ψ22=12(ψ12+ψ22+2Re(ψ1ψ2))|\Psi|^2 = \frac{1}{2}|\psi_1 + \psi_2|^2 = \frac{1}{2}(|\psi_1|^2 + |\psi_2|^2 + 2\text{Re}(\psi_1^* \psi_2))

第三项 2Re(ψ1ψ2)2\text{Re}(\psi_1^* \psi_2) 就是干涉项——它可以是正的(相长干涉)也可以是负的(相消干涉)。这是量子力学最核心的特征之一。

graph TD
    subgraph 概率振幅的干涉
        A["ψ₁ 概率幅"] --> B["叠加 ψ = ψ₁ + ψ₂"]
        C["ψ₂ 概率幅"] --> B
        B --> D["|ψ|² = |ψ₁|² + |ψ₂|² + 2Re(ψ₁*ψ₂)"]
        D --> E[干涉项]
        E --> F["相长: 概率增大"]
        E --> G["相消: 概率减小"]
    end
    
    style F fill:#e8f5e9
    style G fill:#ffebee

1.4 期望值与不确定性

如果每次测量结果都是概率性的,我们如何描述"典型"结果?

1.4.1 期望值

位置 xx期望值(平均值):

x=xΨ(x,t)2dx\langle x \rangle = \int_{-\infty}^{\infty} x |\Psi(x,t)|^2 dx

注意:期望值不是"最可能的结果",而是大量重复实验的平均值

graph LR
    A[制备N个相同系统] --> B[对每个系统测量x]
    B --> C[记录所有结果]
    C --> D[计算平均值]
    D --> E["≈ 期望值 ⟨x⟩"]
    
    style E fill:#e8f5e9

物理直觉:期望值就像统计力学中的系综平均。你没法对一个粒子测量1000次得到期望值——每次测量都会改变系统(波函数坍缩)。但你可以制备1000个完全相同的系统,分别测量,然后取平均。

例题详解:高斯波包的期望值

问题:对于归一化波函数 Ψ(x)=(2aπ)1/4ea(xx0)2\Psi(x) = \left(\frac{2a}{\pi}\right)^{1/4} e^{-a(x-x_0)^2},求位置的期望值。

分步推导

第一步:写出期望值积分

x=x(2aπ)1/2e2a(xx0)2dx\langle x \rangle = \int_{-\infty}^{\infty} x \left(\frac{2a}{\pi}\right)^{1/2} e^{-2a(x-x_0)^2} dx

第二步:换元 y=xx0y = x - x_0

x=(2aπ)1/2(y+x0)e2ay2dy\langle x \rangle = \left(\frac{2a}{\pi}\right)^{1/2} \int_{-\infty}^{\infty} (y + x_0) e^{-2ay^2} dy

第三步:拆分积分

=(2aπ)1/2[ye2ay2dy+x0e2ay2dy]= \left(\frac{2a}{\pi}\right)^{1/2} \left[\int_{-\infty}^{\infty} y e^{-2ay^2} dy + x_0 \int_{-\infty}^{\infty} e^{-2ay^2} dy\right]

第四步:第一项是奇函数积分,结果为0;第二项是高斯积分

=0+x01=x0= 0 + x_0 \cdot 1 = x_0

结论:高斯波包的期望值正好在波包中心 x0x_0。这符合直觉——对称的波包"平均位置"就在中心。

数值例题:三角波包的期望值

问题:粒子在一维空间中的波函数为:

\Psi(x) = \begin{cases} \sqrt{2}x, & 0 \leq x \leq 1 \\ 0, & \text{其他} \end{cases}

(单位已归一化)求位置的期望值 x\langle x \rangle

解答
第一步:验证归一化

01Ψ2dx=012x2dx=213=231\int_0^1 |\Psi|^2 dx = \int_0^1 2x^2 dx = 2 \cdot \frac{1}{3} = \frac{2}{3} \neq 1

波函数未归一化!先求归一化常数 AA

A2012x2dx=1A223=1A=32|A|^2 \int_0^1 2x^2 dx = 1 \Rightarrow |A|^2 \cdot \frac{2}{3} = 1 \Rightarrow A = \sqrt{\frac{3}{2}}

归一化波函数:Ψ(x)=3x\Psi(x) = \sqrt{3}x(在 [0,1][0,1] 内)。

第二步:计算期望值

x=01x3x2dx=301x3dx=314=34=0.75\langle x \rangle = \int_0^1 x \cdot 3x^2 dx = 3 \int_0^1 x^3 dx = 3 \cdot \frac{1}{4} = \frac{3}{4} = 0.75

结论:粒子平均位置在 x=0.75x = 0.75 处。这符合直觉——波函数在 x=1x=1 处最大,概率密度 Ψ2=3x2|\Psi|^2 = 3x^2 偏向右侧,所以平均值偏向中间偏右。


1.4.2 标准差与不确定性

单次测量结果围绕期望值的"分散程度":

σx=x2x2\sigma_x = \sqrt{\langle x^2 \rangle - \langle x \rangle^2}

这就是位置的不确定性。量子力学的一个核心发现是:不确定性不是测量技术不够精确,而是量子系统的内在属性。

例题详解:高斯波包的不确定性

问题:计算上述高斯波包的 σx\sigma_x

分步推导

第一步:计算 x2\langle x^2 \rangle

x2=x2(2aπ)1/2e2a(xx0)2dx\langle x^2 \rangle = \int_{-\infty}^{\infty} x^2 \left(\frac{2a}{\pi}\right)^{1/2} e^{-2a(x-x_0)^2} dx

第二步:换元 y=xx0y = x - x_0

x2=(2aπ)1/2(y2+2yx0+x02)e2ay2dy\langle x^2 \rangle = \left(\frac{2a}{\pi}\right)^{1/2} \int_{-\infty}^{\infty} (y^2 + 2yx_0 + x_0^2) e^{-2ay^2} dy

第三步:利用高斯积分公式

  • ye2ay2dy=0\int_{-\infty}^{\infty} y e^{-2ay^2} dy = 0(奇函数)
  • y2e2ay2dy=14aπ2a\int_{-\infty}^{\infty} y^2 e^{-2ay^2} dy = \frac{1}{4a}\sqrt{\frac{\pi}{2a}}

x2=14a+x02\langle x^2 \rangle = \frac{1}{4a} + x_0^2

第四步:计算标准差

σx=x2x2=14a+x02x02=12a\sigma_x = \sqrt{\langle x^2 \rangle - \langle x \rangle^2} = \sqrt{\frac{1}{4a} + x_0^2 - x_0^2} = \frac{1}{2\sqrt{a}}

物理意义aa 越大,波包越窄,σx\sigma_x 越小——位置越确定。


1.5 动量:波函数的另一面

位置和动量是一对"共轭变量"。在量子力学中,动量不是"质量乘速度",而是与波函数的空间频率相关。

1.5.1 动量算符

动量在量子力学中是一个算符

p^=ix\hat{p} = -i\hbar \frac{\partial}{\partial x}

为什么动量是导数?因为高频振荡的波函数意味着粒子动量大——这来自于德布罗意的物质波假设 p=h/λp = h/\lambda

动量的期望值:

p=Ψ(ix)Ψdx\langle p \rangle = \int_{-\infty}^{\infty} \Psi^* \left(-i\hbar \frac{\partial}{\partial x}\right) \Psi dx

物理直觉:为什么动量是导数?

想象一个平面波 eikxe^{ikx}。对它求导:

ixeikx=i(ik)eikx=keikx-i\hbar \frac{\partial}{\partial x} e^{ikx} = -i\hbar (ik) e^{ikx} = \hbar k \cdot e^{ikx}

而德布罗意关系告诉我们 p=kp = \hbar k。所以动量算符作用在平面波上,正好给出它的动量本征值!

graph TD
    subgraph 动量算符的物理起源
        A["平面波 e^(ikx)"] --> B["波长 λ = 2π/k"]
        B --> C["德布罗意: p = h/λ = ℏk"]
        C --> D["动量算符 p̂ = -iℏ ∂/∂x"]
        D --> E[作用在平面波上]
        E --> F["p̂ e^(ikx) = ℏk e^(ikx)"]
        F --> G["本征值 = 动量"]
    end
    
    style G fill:#e8f5e9

例题详解:高斯波包的动量期望值

问题:计算 Ψ(x)=(2aπ)1/4eax2\Psi(x) = \left(\frac{2a}{\pi}\right)^{1/4} e^{-ax^2} 的动量期望值。

分步推导

第一步:计算导数

Ψx=(2aπ)1/4(2ax)eax2\frac{\partial \Psi}{\partial x} = \left(\frac{2a}{\pi}\right)^{1/4} (-2ax) e^{-ax^2}

第二步:代入期望值公式

p=iΨΨxdx\langle p \rangle = -i\hbar \int_{-\infty}^{\infty} \Psi^* \frac{\partial \Psi}{\partial x} dx

=i(2aπ)1/2eax2(2ax)eax2dx= -i\hbar \left(\frac{2a}{\pi}\right)^{1/2} \int_{-\infty}^{\infty} e^{-ax^2} (-2ax) e^{-ax^2} dx

=2ia(2aπ)1/2xe2ax2dx= 2ia\hbar \left(\frac{2a}{\pi}\right)^{1/2} \int_{-\infty}^{\infty} x e^{-2ax^2} dx

第三步:积分是奇函数积分,结果为0

p=0\langle p \rangle = 0

物理意义:对称的高斯波包"静止"在中心——平均动量为零。如果波包中心以速度 vv 运动,则 p=mv\langle p \rangle = mv

1.5.2 海森堡不确定性原理

这是量子力学最著名的结论:

σxσp2\sigma_x \sigma_p \geq \frac{\hbar}{2}

含义:粒子位置越确定,动量就越不确定;反之亦然。这不是测量仪器的限制,而是波函数的数学性质——一个尖锐的脉冲(位置确定)必须由很多不同波长的波叠加而成(动量不确定)。

例题详解:验证高斯波包满足不确定性原理

问题:对于 Ψ(x)=(2aπ)1/4eax2\Psi(x) = \left(\frac{2a}{\pi}\right)^{1/4} e^{-ax^2},验证 σxσp2\sigma_x \sigma_p \geq \frac{\hbar}{2}

分步推导

第一步:我们已经知道 σx=12a\sigma_x = \frac{1}{2\sqrt{a}}

第二步:计算 p2\langle p^2 \rangle

p2=2Ψ2Ψx2dx\langle p^2 \rangle = -\hbar^2 \int_{-\infty}^{\infty} \Psi^* \frac{\partial^2 \Psi}{\partial x^2} dx

计算二阶导数:

2Ψx2=(2aπ)1/4(4a2x22a)eax2\frac{\partial^2 \Psi}{\partial x^2} = \left(\frac{2a}{\pi}\right)^{1/4} (4a^2 x^2 - 2a) e^{-ax^2}

第三步:代入积分

p2=2(2aπ)1/2eax2(4a2x22a)eax2dx\langle p^2 \rangle = -\hbar^2 \left(\frac{2a}{\pi}\right)^{1/2} \int_{-\infty}^{\infty} e^{-ax^2} (4a^2 x^2 - 2a) e^{-ax^2} dx

=2(2aπ)1/2[4a2x2e2ax2dx2ae2ax2dx]= -\hbar^2 \left(\frac{2a}{\pi}\right)^{1/2} \left[4a^2 \int_{-\infty}^{\infty} x^2 e^{-2ax^2} dx - 2a \int_{-\infty}^{\infty} e^{-2ax^2} dx\right]

第四步:利用高斯积分

x2e2ax2dx=14aπ2a\int_{-\infty}^{\infty} x^2 e^{-2ax^2} dx = \frac{1}{4a}\sqrt{\frac{\pi}{2a}}

e2ax2dx=π2a\int_{-\infty}^{\infty} e^{-2ax^2} dx = \sqrt{\frac{\pi}{2a}}

第五步:代入

p2=2(2aπ)1/2[4a214aπ2a2aπ2a]\langle p^2 \rangle = -\hbar^2 \left(\frac{2a}{\pi}\right)^{1/2} \left[4a^2 \cdot \frac{1}{4a}\sqrt{\frac{\pi}{2a}} - 2a \cdot \sqrt{\frac{\pi}{2a}}\right]

=2(2aπ)1/2[aπ2a2aπ2a]= -\hbar^2 \left(\frac{2a}{\pi}\right)^{1/2} \left[a \sqrt{\frac{\pi}{2a}} - 2a \sqrt{\frac{\pi}{2a}}\right]

=2(2aπ)1/2[aπ2a]= -\hbar^2 \left(\frac{2a}{\pi}\right)^{1/2} \left[-a \sqrt{\frac{\pi}{2a}}\right]

=2a= \hbar^2 a

第六步:计算 σp\sigma_p

σp=p2p2=2a0=a\sigma_p = \sqrt{\langle p^2 \rangle - \langle p \rangle^2} = \sqrt{\hbar^2 a - 0} = \hbar \sqrt{a}

第七步:验证不确定性原理

σxσp=12aa=2\sigma_x \sigma_p = \frac{1}{2\sqrt{a}} \cdot \hbar \sqrt{a} = \frac{\hbar}{2}

结论:高斯波包恰好"饱和"了不确定性原理——它是最小不确定波包!

graph TD
    subgraph 不确定性原理的本质
        A["位置精确
σₓ小"] -->|傅里叶变换| B["动量分散
σₚ大"] C["动量精确
σₚ小"] -->|傅里叶变换| D["位置分散
σₓ大"] E[高斯波包] --> F["最小不确定态
σₓσₚ = ℏ/2"] end style A fill:#ffebee style B fill:#ffebee style C fill:#e3f2fd style D fill:#e3f2fd style F fill:#e8f5e9

1.5.3 能量-时间不确定性原理

除了位置-动量不确定性,量子力学还有另一对著名的不确定性关系:

ΔEΔt2\Delta E \Delta t \geq \frac{\hbar}{2}

这里的 ΔE\Delta E 是能量的不确定度,Δt\Delta t 是时间的不确定度。

注意:能量-时间不确定性与位置-动量不确定性有本质区别。位置和动量是可观测量,有对应的算符 x^\hat{x}p^\hat{p}。但时间不是可观测量——在标准量子力学中,没有时间算符 t^\hat{t}。因此,能量-时间不确定性不能直接从对易子 [H^,t^][\hat{H}, \hat{t}] 推导。

那么它从何而来? Griffiths 提供了几种理解方式:

理解1:寿命与能级宽度
一个不稳定粒子(如激发态原子)只能存在有限时间 Δt\Delta t(寿命)。能量-时间不确定性告诉我们,它的能量不能精确确定,而是有一个"自然线宽" ΔE/(2Δt)\Delta E \sim \hbar/(2\Delta t)。寿命越短,能级越宽。这就是为什么光谱线总有一定宽度,而不是理想的δ\delta函数。

理解2:态的演化
如果一个系统处于 H^\hat{H} 的本征态,能量完全确定(ΔE=0\Delta E = 0),那么它的概率分布永不变化(定态)——相当于 Δt\Delta t \to \infty,系统"永远保持这样"。反之,如果系统处于能量本征态的叠加,概率分布会随时间演化,这种演化的时间尺度 Δt\Delta t 与能量差 ΔE\Delta E 满足不确定性关系。

理解3:测量时间
如果你测量一个系统的能量,需要花费时间 Δt\Delta t。不确定性原理说,测量时间越长,能量可以测得越精确。快速测量(Δt\Delta t 小)必然导致能量精度下降(ΔE\Delta E 大)。

数值例题:原子激发态的寿命与线宽

问题:一个原子激发态的寿命约为 10810^{-8} s。估算其发射光谱的自然线宽(能量不确定度)和对应波长的不确定度(发射波长约为 500 nm)。

解答
第一步:估算能量不确定度

ΔE2Δt=1.055×10342×1085.3×1027 J\Delta E \geq \frac{\hbar}{2\Delta t} = \frac{1.055 \times 10^{-34}}{2 \times 10^{-8}} \approx 5.3 \times 10^{-27} \text{ J}

转换为 eV:

ΔE5.3×10271.602×10193.3×108 eV33 neV\Delta E \approx \frac{5.3 \times 10^{-27}}{1.602 \times 10^{-19}} \approx 3.3 \times 10^{-8} \text{ eV} \approx 33 \text{ neV}

第二步:波长与能量的关系

E=hcλΔE=hcλ2ΔλE = \frac{hc}{\lambda} \Rightarrow \Delta E = \frac{hc}{\lambda^2} \Delta \lambda

(对 E1/λE \propto 1/\lambda 微分)

第三步:求解波长不确定度

Δλ=λ2ΔEhc=(500×109)2×5.3×1027(6.626×1034)(2.998×108)\Delta \lambda = \frac{\lambda^2 \Delta E}{hc} = \frac{(500 \times 10^{-9})^2 \times 5.3 \times 10^{-27}}{(6.626 \times 10^{-34})(2.998 \times 10^8)}

2.5×1013×5.3×10271.986×1025\approx \frac{2.5 \times 10^{-13} \times 5.3 \times 10^{-27}}{1.986 \times 10^{-25}}

1.33×10391.986×10256.7×1015 m=6.7 fm\approx \frac{1.33 \times 10^{-39}}{1.986 \times 10^{-25}} \approx 6.7 \times 10^{-15} \text{ m} = 6.7 \text{ fm}

结论:发射光谱的自然线宽约为 33 neV,对应波长不确定度约为 6.7 fm(飞米)。这是一个极小的量,但在精密光谱学中是可以测量的。

物理直觉:想象一个钟摆。如果它永远摆动(无限寿命),它的频率是精确确定的。但如果钟摆只摆动几秒钟就停了,你很难准确判断它的频率——摆动时间越短,频率越模糊。原子激发态也是这样:存在时间越短,发射光的"颜色"越模糊。


1.6 动量空间波函数:另一扇窗

我们已经知道 Ψ(x,t)\Psi(x,t) 描述粒子在空间中的概率幅。但量子力学还有一个完全等价的描述:动量空间波函数 Φ(p,t)\Phi(p,t)

1.6.1 从位置到动量的变换

两者通过傅里叶变换联系:

Φ(p,t)=12πΨ(x,t)eipx/dx\Phi(p,t) = \frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty} \Psi(x,t) e^{-ipx/\hbar} dx

逆变换:

Ψ(x,t)=12πΦ(p,t)eipx/dp\Psi(x,t) = \frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty} \Phi(p,t) e^{ipx/\hbar} dp

物理直觉:位置波函数描述"粒子在哪里",动量波函数描述"粒子以什么动量运动"。两者包含完全相同的信息,只是视角不同——就像同一物体从不同角度拍摄的照片。

graph LR
    subgraph 两种表象
        A["#quot;位置表象
Ψ(x)"] <-->|傅里叶变换| B["动量表象
Φ(p)#quot;"] A --> C["测量x的概率
|Ψ(x)|²"] B --> D["测量p的概率
|Φ(p)|²"] end style A fill:#e3f2fd style B fill:#fff3e0

1.6.2 为什么两种表象都重要?

某些问题在位置空间简单(如势阱问题),某些问题在动量空间更简单(如自由粒子)。这种"表象变换"的思想将在第3章形式主义中成为核心——量子力学的所有物理都可以在任何表象中表达。

数值例题:从位置波函数求动量分布

问题:一个粒子的位置波函数为 Ψ(x)=1aex/a\Psi(x) = \frac{1}{\sqrt{a}} e^{-|x|/a}a=1a = 1 nm),求其动量空间波函数 Φ(p)\Phi(p),并估算动量分布的半高全宽(FWHM)。

解答
第一步:写出傅里叶变换

Φ(p)=12π1aex/aeipx/dx\Phi(p) = \frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty} \frac{1}{\sqrt{a}} e^{-|x|/a} e^{-ipx/\hbar} dx

第二步:拆分积分

Φ(p)=12πa[0ex/aeipx/dx+0ex/aeipx/dx]\Phi(p) = \frac{1}{\sqrt{2\pi\hbar a}} \left[\int_0^{\infty} e^{-x/a} e^{-ipx/\hbar} dx + \int_{-\infty}^0 e^{x/a} e^{-ipx/\hbar} dx\right]

第三步:分别计算
第一个积分(x>0x > 0):

0e(1/a+ip/)xdx=11a+ip=a+iap\int_0^{\infty} e^{-(1/a + ip/\hbar)x} dx = \frac{1}{\frac{1}{a} + \frac{ip}{\hbar}} = \frac{a\hbar}{\hbar + iap}

第二个积分(x<0x < 0):

0e(1/aip/)xdx=11aip=aiap\int_{-\infty}^0 e^{(1/a - ip/\hbar)x} dx = \frac{1}{\frac{1}{a} - \frac{ip}{\hbar}} = \frac{a\hbar}{\hbar - iap}

第四步:合并

Φ(p)=12πa[a+iap+aiap]\Phi(p) = \frac{1}{\sqrt{2\pi\hbar a}} \left[\frac{a\hbar}{\hbar + iap} + \frac{a\hbar}{\hbar - iap}\right]

=a2π2()2+(ap)2= \frac{\sqrt{a\hbar}}{\sqrt{2\pi}} \cdot \frac{2\hbar}{(\hbar)^2 + (ap)^2}

=2a3π1(ap)2+2= \sqrt{\frac{2a^3\hbar}{\pi}} \cdot \frac{1}{(a p)^2 + \hbar^2}

第五步:分析动量分布

Φ(p)2=2a3π1[(ap)2+2]2|\Phi(p)|^2 = \frac{2a^3\hbar}{\pi} \cdot \frac{1}{[(ap)^2 + \hbar^2]^2}

这是一个洛伦兹型的平方分布。最大值在 p=0p = 0 处。

半高全宽(FWHM):令 Φ(p)2=12Φ(0)2|\Phi(p)|^2 = \frac{1}{2}|\Phi(0)|^2,解得:

pFWHMap_{\text{FWHM}} \approx \frac{\hbar}{a}

代入 a=1a = 1 nm:

p_{\text{FWHM}} \approx \frac{1.055 \times 10^{-34}}{10^{-9}} = 1.055 \times 10^{-25} \text{ kg·m/s}

结论:位置越局域化(aa 越小),动量越弥散(pFWHMp_{\text{FWHM}} 越大)。这正是位置-动量不确定性的直接体现——一个紧窄的位置波包必须由大范围的动量分量叠加而成。


1.7 本章总结

graph TD
    A[第1章核心概念] --> B["波函数 Ψ"]
    A --> C[薛定谔方程]
    A --> D[玻恩诠释]
    A --> E[不确定性原理]
    A --> F[动量空间]
    
    B -->|描述| G[粒子所有可能性的叠加]
    C -->|决定| H[波函数如何演化]
    D -->|说明| I["|Ψ|² = 概率密度"]
    E -->|约束| J["σₓσₚ ≥ ℏ/2"]
    F -->|等价描述| K[傅里叶变换]
    
    G --> L["测量 → 坍缩到某一本征态"]
    I --> L
    K --> B
    
    style G fill:#fff3e0
    style L fill:#e8f5e9
    style K fill:#e3f2fd

带走的三句话:

  1. 波函数不是粒子本身,而是描述粒子所有可能性的复数振幅。
  2. Ψ2|\Psi|^2 是概率密度,不是物理密度。
  3. 不确定性是量子系统的内在属性,不是测量误差。

与其他章节的联系

graph LR
    A[第1章] --> B["第2章
求解定态方程"] A --> C["第3章
形式主义与算符"] A --> D["第4章
三维波函数"] B --> E["定态解 = 能量本征态"] C --> F[动量算符是形式化的起点] D --> G[位置空间扩展到3D] style A fill:#e3f2fd style B fill:#e8f5e9 style C fill:#e8f5e9 style D fill:#e8f5e9

"量子力学要求一种全新的、彻底反直觉的世界观。但别担心——你会习惯的。" —— Griffiths


1.8 练习与思考

  1. 归一化检验:验证 Ψ(x)=Aeax\Psi(x) = A e^{-a|x|} 是否可归一化?如果是,AA 等于多少?并计算该波函数的位置期望值和不确定性 σx\sigma_x

  2. 概率计算:对于归一化的高斯波包 Ψ(x)=(2aπ)1/4eax2\Psi(x) = \left(\frac{2a}{\pi}\right)^{1/4} e^{-ax^2},计算在 [1/2a,+1/2a][-1/\sqrt{2a}, +1/\sqrt{2a}] 内找到粒子的概率。(提示:利用误差函数 erf(x)=2π0xet2dt\text{erf}(x) = \frac{2}{\sqrt{\pi}} \int_0^x e^{-t^2} dt

  3. 不确定性估算:粗略估算:如果一个电子被限制在原子核尺度(1015\sim 10^{-15} m),它的动量不确定性有多大?对应的动能是多少?(提示:这就是电子为什么不会掉进原子核的原因之一。)