本文系统地介绍量子力学的基本原理、数学框架和核心概念,涵盖波粒二象性、薛定谔方程、算符理论、测量问题以及量子态的数学描述。通过具体的数值例子和实验数据,帮助读者建立对微观世界运作规律的直观理解。
量子力学(Quantum Mechanics)是描述微观粒子(原子、电子、光子等)运动规律的理论体系,与相对论并列为现代物理学的两大支柱。它彻底改变了人类对物质世界的认知——在微观尺度下,粒子表现出与宏观物体截然不同的行为:它们可以同时处于多个状态的叠加中,位置和动量不能同时精确确定,甚至可以通过"隧穿"穿越经典物理认为不可逾越的势垒。
量子力学的诞生不是一蹴而就的,而是经历了约四分之一个世纪的积累:
时间
里程碑
贡献者
意义
1900年
普朗克黑体辐射公式
M. Planck
提出能量量子化假说 E = h ν E = h\nuE = h ν ,开创量子时代
1905年
光电效应解释
A. Einstein
提出光量子概念,证实光的粒子性
1913年
玻尔原子模型
N. Bohr
用量子化轨道解释氢原子光谱
1923年
康普顿散射
A. Compton
实验证实光子具有动量
1924年
物质波假说
L. de Broglie
提出粒子也具有波动性
1925年
矩阵力学
W. Heisenberg
第一个完整的量子力学形式体系
1926年
波动力学
E. Schrödinger
建立薛定谔方程,提出波函数解释
1927年
不确定性原理
W. Heisenberg
揭示量子测量的根本限制
1927年
哥本哈根诠释
Bohr等
提出波函数的统计解释和互补原理
量子力学描述的对象涵盖了从亚原子粒子到宏观量子现象的广阔尺度:
原子核(10⁻¹⁵m) → 原子(10⁻¹⁰m) → 分子(10⁻⁹m) → 纳米结构(10⁻⁹~10⁻⁶m) → 超导宏观量子态
↓ ↓ ↓ ↓ ↓
核物理 原子物理 量子化学 纳米技术 超导/超流
波粒二象性(Wave-Particle Duality)是量子力学的第一个核心概念:微观粒子既表现出粒子的局域性,又表现出波的干涉和衍射特性。在经典物理中,波动和粒子是截然不同的两种存在形式,但在量子世界中,它们是同一实体的不同表现。
19世纪末,物理学家试图用经典理论解释黑体辐射的能谱分布,但遭遇了"紫外灾难"。
瑞利-金斯公式(经典理论) :
u ( ν , T ) = 8 π ν 2 c 3 k B T u(\nu, T) = \frac{8\pi\nu^2}{c^3} k_B T
u ( ν , T ) = c 3 8 π ν 2 k B T
这个公式在低频与实验吻合,但在高频(紫外区域)发散至无穷大,显然是荒谬的。
普朗克公式(量子假说) :
u ( ν , T ) = 8 π ν 2 c 3 ⋅ h ν e h ν / k B T − 1 u(\nu, T) = \frac{8\pi\nu^2}{c^3} \cdot \frac{h\nu}{e^{h\nu/k_B T} - 1}
u ( ν , T ) = c 3 8 π ν 2 ⋅ e h ν / k B T − 1 h ν
数值对比 :假设温度为 T = 5000 K T = 5000\,\text{K}T = 5 0 0 0 K (大约相当于太阳表面温度),计算 ν = 1 0 15 Hz \nu = 10^{15}\,\text{Hz}ν = 1 0 1 5 Hz (可见光频段)处的能量密度:
公式
计算结果
与实验对比
瑞利-金斯(经典)
u = 5.77 × 1 0 − 16 J/m 3 u = 5.77 \times 10^{-16}\,\text{J/m}^3u = 5 . 7 7 × 1 0 − 1 6 J/m 3
严重高估(发散趋势)
普朗克(量子)
u = 7.23 × 1 0 − 17 J/m 3 u = 7.23 \times 10^{-17}\,\text{J/m}^3u = 7 . 2 3 × 1 0 − 1 7 J/m 3
与实验精确吻合
维恩近似
u = 6.88 × 1 0 − 17 J/m 3 u = 6.88 \times 10^{-17}\,\text{J/m}^3u = 6 . 8 8 × 1 0 − 1 7 J/m 3
在高温段略有偏差
普朗克大胆假设:辐射能量不是连续的,而是一份一份的,每份能量为:
E = h ν E = h\nu
E = h ν
其中 h = 6.626 × 1 0 − 34 J ⋅ s h = 6.626 \times 10^{-34}\,\text{J}\cdot\text{s}h = 6 . 6 2 6 × 1 0 − 3 4 J ⋅ s 为普朗克常数。这一假说标志着量子概念的诞生。
爱因斯坦进一步提出,光的能量量子(即光子)不仅存在于辐射过程中,而且本身就是光的本质形态。光电效应实验证实了这一点:
当紫外光照射到金属表面时,会立即释放出电子(光电子)。经典波动理论无法解释以下三个关键实验事实:
阈值频率的存在 :对每种金属,只有当入射光频率 ν > ν 0 \nu > \nu_0ν > ν 0 时才会产生光电子。低于阈值频率时,无论光强多大都不会产生光电子。
瞬时性 :只要 ν > ν 0 \nu > \nu_0ν > ν 0 ,光电子的产生是瞬时的(< 1 0 − 9 s 10^{-9}\,\text{s}1 0 − 9 s ),无需能量积累时间。
光电子动能与光强无关 :光电子的最大动能只取决于频率,与光强无关。
爱因斯坦的光量子解释:
E k = h ν − W E_k = h\nu - W
E k = h ν − W
其中 W WW 是金属的逸出功。对于钠金属(W = 2.28 eV W = 2.28\,\text{eV}W = 2 . 2 8 eV ):
入射光频率 (× 1 0 14 Hz \times 10^{14}\,\text{Hz}× 1 0 1 4 Hz )
波长 (nm)
光子能量 (eV)
光电子最大动能 (eV)
4.0
750
1.66
0(低于阈值)
5.5
545
2.28
0(恰好阈值)
7.5
400
3.11
0.83
10.0
300
4.14
1.86
15.0
200
6.21
3.93
阈值频率 ν 0 = W / h = 2.28 × 1.602 × 1 0 − 19 / 6.626 × 1 0 − 34 ≈ 5.5 × 1 0 14 Hz \nu_0 = W/h = 2.28 \times 1.602 \times 10^{-19} / 6.626 \times 10^{-34} \approx 5.5 \times 10^{14}\,\text{Hz}ν 0 = W / h = 2 . 2 8 × 1 . 6 0 2 × 1 0 − 1 9 / 6 . 6 2 6 × 1 0 − 3 4 ≈ 5 . 5 × 1 0 1 4 Hz ,对应波长 545nm(绿光)。
1924年,法国物理学家德布罗意(Louis de Broglie)在他的博士论文中提出了一个革命性的观点:正如光具有波粒二象性一样,实物粒子(电子、质子等)也具有波动性。他给出了物质波的波长与动量之间的关系:
λ = h p = h m v \lambda = \frac{h}{p} = \frac{h}{mv}
λ = p h = m v h
其中 λ \lambdaλ 是德布罗意波长,p pp 是粒子的动量,m mm 是质量,v vv 是速度。
具体数值示例 :
粒子
质量 (kg)
速度 (m/s)
动量 (kg·m/s)
德布罗意波长 (m)
电子(1V加速)
9.11 × 1 0 − 31 9.11 \times 10^{-31}9 . 1 1 × 1 0 − 3 1
5.93 × 1 0 5 5.93 \times 10^55 . 9 3 × 1 0 5
5.4 × 1 0 − 25 5.4 \times 10^{-25}5 . 4 × 1 0 − 2 5
1.23 × 1 0 − 9 1.23 \times 10^{-9}1 . 2 3 × 1 0 − 9 (约1.23nm)
电子(100V加速)
9.11 × 1 0 − 31 9.11 \times 10^{-31}9 . 1 1 × 1 0 − 3 1
5.93 × 1 0 6 5.93 \times 10^65 . 9 3 × 1 0 6
5.4 × 1 0 − 24 5.4 \times 10^{-24}5 . 4 × 1 0 − 2 4
1.23 × 1 0 − 10 1.23 \times 10^{-10}1 . 2 3 × 1 0 − 1 0 (约0.123nm)
质子(1MeV)
1.67 × 1 0 − 27 1.67 \times 10^{-27}1 . 6 7 × 1 0 − 2 7
1.39 × 1 0 7 1.39 \times 10^71 . 3 9 × 1 0 7
2.32 × 1 0 − 20 2.32 \times 10^{-20}2 . 3 2 × 1 0 − 2 0
2.86 × 1 0 − 14 2.86 \times 10^{-14}2 . 8 6 × 1 0 − 1 4
网球(100km/h)
5.8 × 1 0 − 2 5.8 \times 10^{-2}5 . 8 × 1 0 − 2
27.8
1.61
4.1 × 1 0 − 34 4.1 \times 10^{-34}4 . 1 × 1 0 − 3 4 (不可观测)
地球绕太阳
5.97 × 1 0 24 5.97 \times 10^{24}5 . 9 7 × 1 0 2 4
2.98 × 1 0 4 2.98 \times 10^42 . 9 8 × 1 0 4
1.78 × 1 0 29 1.78 \times 10^{29}1 . 7 8 × 1 0 2 9
3.72 × 1 0 − 63 3.72 \times 10^{-63}3 . 7 2 × 1 0 − 6 3 (完全不可观测)
从上表可见,宏观物体的德布罗意波长极小,无法被探测到——这就是为什么我们感受不到桌子的波动性。而电子的波长与X射线相当,可以被晶体衍射实验证实。
1927年,戴维森(Davisson)和革末(Germer)在贝尔实验室进行了关键的电子衍射实验。他们将一束已知能量的电子射向镍晶体表面,测量不同角度上的散射电子强度。
实验参数 :
电子能量:54 eV
对应的德布罗意波长:λ = h 2 m E = 6.626 × 1 0 − 34 2 × 9.11 × 1 0 − 31 × 54 × 1.602 × 1 0 − 19 ≈ 1.67 × 1 0 − 10 m = 0.167 nm \lambda = \frac{h}{\sqrt{2mE}} = \frac{6.626\times 10^{-34}}{\sqrt{2\times 9.11\times 10^{-31} \times 54\times 1.602\times 10^{-19}}} \approx 1.67\times 10^{-10}\,\text{m} = 0.167\,\text{nm}λ = 2 m E h = 2 × 9 . 1 1 × 1 0 − 3 1 × 5 4 × 1 . 6 0 2 × 1 0 − 1 9 6 . 6 2 6 × 1 0 − 3 4 ≈ 1 . 6 7 × 1 0 − 1 0 m = 0 . 1 6 7 nm
镍晶体的晶格间距:d = 0.215 nm d = 0.215\,\text{nm}d = 0 . 2 1 5 nm (根据X射线衍射确定)
实验结果 :在散射角 ϕ = 5 0 ∘ \phi = 50^\circϕ = 5 0 ∘ 处观察到明显的衍射峰。根据布拉格衍射定律:
d sin ϕ = n λ d\sin\phi = n\lambda
d sin ϕ = n λ
代入 n = 1 n=1n = 1 ,d = 0.215 nm d=0.215\,\text{nm}d = 0 . 2 1 5 nm ,ϕ = 5 0 ∘ \phi=50^\circϕ = 5 0 ∘ :
λ = 0.215 × 1 0 − 9 × sin 5 0 ∘ = 0.215 × 0.766 = 0.165 nm \lambda = 0.215\times 10^{-9}\times \sin 50^\circ = 0.215\times 0.766 = 0.165\,\text{nm}
λ = 0 . 2 1 5 × 1 0 − 9 × sin 5 0 ∘ = 0 . 2 1 5 × 0 . 7 6 6 = 0 . 1 6 5 nm
与德布罗意波长的理论值 0.167 nm 0.167\,\text{nm}0 . 1 6 7 nm 相差不到 2%。这一完美吻合是物质波假说最直接的实验证据。
双缝实验是理解量子力学最"离奇"性质的关键实验。想象一束电子射向一块有两个狭缝的挡板,后面的探测器屏幕会记录电子到达的位置。
当只开一个缝时 :
只开缝1:屏幕上的分布像粒子穿过一条狭缝的衍射图案
只开缝2:类似,但从第二个缝的位置
结果类似于宏观粒子的弹道分布
当两个缝同时打开时 :
令人惊讶的是,屏幕上并不是两个单缝图案的简单叠加,而是出现了明暗相间的干涉条纹 !
不同粒子的双缝实验数据 :
粒子类型
典型波长
缝间距
屏幕距离
条纹间距
所需探测器
可见光
500nm
0.1mm
1m
5mm
肉眼/CCD
电子(50eV)
0.17nm
1μm
1m
1.7μm
电子倍增器
中子(0.025eV)
0.18nm
100nm
10m
18μm
中子计数器
富勒烯(C₆₀)
2.5pm
100nm
1m
25nm
激光电离探测器
即使每次只发射一个电子(逐个发射),经过足够多的积累后,屏幕上出现的仍然是干涉条纹。这意味着单个电子同时通过了两个缝 并与自身发生了干涉——这是经典物理完全无法理解的行为。
用波函数来描述:当两个缝都打开时,电子的波函数是:
Ψ = Ψ 1 + Ψ 2 \Psi = \Psi_1 + \Psi_2
Ψ = Ψ 1 + Ψ 2
屏幕上的概率密度为:
∣ Ψ ∣ 2 = ∣ Ψ 1 ∣ 2 + ∣ Ψ 2 ∣ 2 + Ψ 1 ∗ Ψ 2 + Ψ 2 ∗ Ψ 1 |\Psi|^2 = |\Psi_1|^2 + |\Psi_2|^2 + \Psi_1^*\Psi_2 + \Psi_2^*\Psi_1
∣ Ψ ∣ 2 = ∣ Ψ 1 ∣ 2 + ∣ Ψ 2 ∣ 2 + Ψ 1 ∗ Ψ 2 + Ψ 2 ∗ Ψ 1
其中前两项是各自单独打开时的概率分布,后两项是干涉项——正是它们产生了干涉条纹。如果我们在缝口放置探测器试图"看到"电子从哪个缝通过,干涉条纹就会消失,波函数坍缩为其中一个缝:
Ψ → { Ψ 1 (如果探测到电子从缝1通过) Ψ 2 (如果探测到电子从缝2通过) \Psi \to \begin{cases} \Psi_1 & \text{(如果探测到电子从缝1通过)}\\ \Psi_2 & \text{(如果探测到电子从缝2通过)} \end{cases}
Ψ → { Ψ 1 Ψ 2 (如果探测到电子从缝 1 通过) (如果探测到电子从缝 2 通过)
波函数 Ψ ( r , t ) \Psi(\mathbf{r}, t)Ψ ( r , t ) 是整个量子力学的核心对象。它是一个复值函数,随空间和时间的改变而改变。对波函数的统计诠释 (由玻恩提出)是量子力学的基石之一:
P ( r , t ) = ∣ Ψ ( r , t ) ∣ 2 d 3 r P(\mathbf{r}, t) = |\Psi(\mathbf{r}, t)|^2 d^3r
P ( r , t ) = ∣ Ψ ( r , t ) ∣ 2 d 3 r
表示在 t tt 时刻、r \mathbf{r}r 附近的体积元 d 3 r d^3rd 3 r 内发现粒子的概率。
波函数的归一化条件 :
∫ − ∞ ∞ ∣ Ψ ( r , t ) ∣ 2 d 3 r = 1 \int_{-\infty}^{\infty} |\Psi(\mathbf{r}, t)|^2 d^3r = 1
∫ − ∞ ∞ ∣ Ψ ( r , t ) ∣ 2 d 3 r = 1
这意味着粒子必然存在于空间的某处。对于满足归一化的波函数,我们可以定义各种物理量的期望值。
举例:高斯波包
一个自由粒子的高斯波包(最常见的基础波函数之一):
Ψ ( x , t ) = ( 1 2 π σ 2 ) 1 / 4 1 1 + i ℏ t / ( 2 m σ 2 ) exp [ − ( x − p 0 t / m ) 2 4 σ 2 ( 1 + i ℏ t / ( 2 m σ 2 ) ) + i p 0 ℏ ( x − p 0 t 2 m ) ] \Psi(x, t) = \left(\frac{1}{2\pi\sigma^2}\right)^{1/4} \frac{1}{\sqrt{1 + i\hbar t/(2m\sigma^2)}} \exp\left[ -\frac{(x - p_0 t/m)^2}{4\sigma^2(1 + i\hbar t/(2m\sigma^2))} + \frac{ip_0}{\hbar}\left(x - \frac{p_0 t}{2m}\right) \right]
Ψ ( x , t ) = ( 2 π σ 2 1 ) 1 / 4 1 + i ℏ t / ( 2 m σ 2 ) 1 exp [ − 4 σ 2 ( 1 + i ℏ t / ( 2 m σ 2 ) ) ( x − p 0 t / m ) 2 + ℏ i p 0 ( x − 2 m p 0 t ) ]
其中 σ \sigmaσ 是初始位置的展宽(标准差),p 0 p_0p 0 是平均动量。
具体数值模拟 :假设一个电子(m = 9.11 × 1 0 − 31 kg m = 9.11 \times 10^{-31}\,\text{kg}m = 9 . 1 1 × 1 0 − 3 1 kg ),初始位置标准差 σ = 1.0 × 1 0 − 10 m \sigma = 1.0 \times 10^{-10}\,\text{m}σ = 1 . 0 × 1 0 − 1 0 m (约一个原子大小),平均速度 v 0 = 1.0 × 1 0 6 m/s v_0 = 1.0 \times 10^6\,\text{m/s}v 0 = 1 . 0 × 1 0 6 m/s 。
时间 (s)
波包中心位置 (m)
波包宽度 (m)
宽度的变化
t = 0 t = 0t = 0
0
1.0 × 1 0 − 10 1.0\times 10^{-10}1 . 0 × 1 0 − 1 0
—
t = 1 0 − 16 t = 10^{-16}t = 1 0 − 1 6
1.0 × 1 0 − 10 1.0\times 10^{-10}1 . 0 × 1 0 − 1 0
1.5 × 1 0 − 10 1.5\times 10^{-10}1 . 5 × 1 0 − 1 0
展宽了约50%
t = 1 0 − 15 t = 10^{-15}t = 1 0 − 1 5
1.0 × 1 0 − 9 1.0\times 10^{-9}1 . 0 × 1 0 − 9
4.4 × 1 0 − 10 4.4\times 10^{-10}4 . 4 × 1 0 − 1 0
展宽了4倍多
t = 1 0 − 14 t = 10^{-14}t = 1 0 − 1 4
1.0 × 1 0 − 8 1.0\times 10^{-8}1 . 0 × 1 0 − 8
4.3 × 1 0 − 9 4.3\times 10^{-9}4 . 3 × 1 0 − 9
展宽了40多倍
可见,量子粒子的波包随时间的推移会迅速展宽——这就是自由粒子"扩散"的量子效应。宏观物体由于质量巨大而不会出现可观测的扩散(同样的公式,一个1kg的物体的波包需要比宇宙年龄还长的时间才会显著展宽)。
在经典力学中,粒子运动由牛顿方程 F = m a F = maF = m a 描述。在量子力学中,波函数的演化由薛定谔方程控制。
含时薛定谔方程 (TDSE):
i ℏ ∂ ∂ t Ψ ( r , t ) = [ − ℏ 2 2 m ∇ 2 + V ( r , t ) ] Ψ ( r , t ) i\hbar \frac{\partial}{\partial t} \Psi(\mathbf{r}, t) = \left[ -\frac{\hbar^2}{2m} \nabla^2 + V(\mathbf{r}, t) \right] \Psi(\mathbf{r}, t)
i ℏ ∂ t ∂ Ψ ( r , t ) = [ − 2 m ℏ 2 ∇ 2 + V ( r , t ) ] Ψ ( r , t )
其中 ℏ = h / 2 π ≈ 1.055 × 1 0 − 34 J ⋅ s \hbar = h/2\pi \approx 1.055 \times 10^{-34}\,\text{J}\cdot\text{s}ℏ = h / 2 π ≈ 1 . 0 5 5 × 1 0 − 3 4 J ⋅ s ,∇ 2 \nabla^2∇ 2 是拉普拉斯算子,V VV 是势能函数。
定态薛定谔方程 (TISE):当势能 V VV 不显含时间时,波函数可以分离变量为:
Ψ ( r , t ) = ψ ( r ) e − i E t / ℏ \Psi(\mathbf{r}, t) = \psi(\mathbf{r}) e^{-iEt/\hbar}
Ψ ( r , t ) = ψ ( r ) e − i E t / ℏ
代入含时方程得到:
− ℏ 2 2 m ∇ 2 ψ ( r ) + V ( r ) ψ ( r ) = E ψ ( r ) -\frac{\hbar^2}{2m} \nabla^2\psi(\mathbf{r}) + V(\mathbf{r})\psi(\mathbf{r}) = E\psi(\mathbf{r})
− 2 m ℏ 2 ∇ 2 ψ ( r ) + V ( r ) ψ ( r ) = E ψ ( r )
或者写成哈密顿算符形式:
H ^ ψ = E ψ \hat{H}\psi = E\psi
H ^ ψ = E ψ
其中哈密顿算符 H ^ = − ℏ 2 2 m ∇ 2 + V ( r ) \hat{H} = -\frac{\hbar^2}{2m}\nabla^2 + V(\mathbf{r})H ^ = − 2 m ℏ 2 ∇ 2 + V ( r ) ,E EE 是能量本征值。
这是量子力学中最简单的模型之一,可以严格求解,展示量子化的本质。
设势阱宽度为 a aa ,势能为:
V ( x ) = { 0 , 0 < x < a ∞ , 其他 V(x) = \begin{cases} 0, & 0 < x < a \\ \infty, & \text{其他} \end{cases}
V ( x ) = { 0 , ∞ , 0 < x < a 其他
求解定态薛定谔方程得到:
ψ n ( x ) = 2 a sin ( n π x a ) , n = 1 , 2 , 3 , … \psi_n(x) = \sqrt{\frac{2}{a}} \sin\left(\frac{n\pi x}{a}\right), \quad n=1,2,3,\ldots
ψ n ( x ) = a 2 sin ( a n π x ) , n = 1 , 2 , 3 , …
E n = n 2 π 2 ℏ 2 2 m a 2 = n 2 E 1 E_n = \frac{n^2\pi^2\hbar^2}{2ma^2} = n^2 E_1
E n = 2 m a 2 n 2 π 2 ℏ 2 = n 2 E 1
具体数值 :考虑一个电子(m e = 9.11 × 1 0 − 31 kg m_e = 9.11 \times 10^{-31}\,\text{kg}m e = 9 . 1 1 × 1 0 − 3 1 kg )在宽度 a = 0.1 nm a = 0.1\,\text{nm}a = 0 . 1 nm (约一个原子大小)的势阱中:
量子数 n nn
能量 E n E_nE n
能量(eV)
与经典对比
1(基态)
6.02 × 1 0 − 19 J 6.02 \times 10^{-19}\,\text{J}6 . 0 2 × 1 0 − 1 9 J
3.76 eV
经典允许任意连续能量
2(第一激发态)
2.41 × 1 0 − 18 J 2.41 \times 10^{-18}\,\text{J}2 . 4 1 × 1 0 − 1 8 J
15.0 eV
—
3
5.42 × 1 0 − 18 J 5.42 \times 10^{-18}\,\text{J}5 . 4 2 × 1 0 − 1 8 J
33.8 eV
—
4
9.63 × 1 0 − 18 J 9.63 \times 10^{-18}\,\text{J}9 . 6 3 × 1 0 − 1 8 J
60.1 eV
—
10
6.02 × 1 0 − 17 J 6.02 \times 10^{-17}\,\text{J}6 . 0 2 × 1 0 − 1 7 J
376 eV
—
关键发现:
能量是量子化的 ,只能取离散值 E n = n 2 E 1 E_n = n^2 E_1E n = n 2 E 1 ——这与经典物理中能量连续的观念完全不同
零点能 :即使 n = 1 n=1n = 1 的最低能态,能量也不为零(E 1 > 0 E_1 > 0E 1 > 0 ),这是海森堡不确定性原理的必然结果
波函数有 n − 1 n-1n − 1 个节点(过零点),n nn 越大振荡越快
作为对比,如果是一个1g的宏观物体在1cm宽的"势阱"中,相邻能级差为 Δ E ≈ 1 0 − 63 eV \Delta E \approx 10^{-63}\,\text{eV}Δ E ≈ 1 0 − 6 3 eV ,远远超出任何实验的探测能力——量子效应完全被"抹平"了。
谐振子是量子力学中极为重要的模型,因为它可以近似任何势能在极小值附近的行为(泰勒展开到二阶)。
势能:V ( x ) = 1 2 m ω 2 x 2 V(x) = \frac{1}{2}m\omega^2 x^2V ( x ) = 2 1 m ω 2 x 2 ,其中 ω \omegaω 是经典角频率。
薛定谔方程:
− ℏ 2 2 m d 2 ψ d x 2 + 1 2 m ω 2 x 2 ψ = E ψ -\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2} + \frac{1}{2}m\omega^2 x^2\psi = E\psi
− 2 m ℏ 2 d x 2 d 2 ψ + 2 1 m ω 2 x 2 ψ = E ψ
能量本征值 :
E n = ( n + 1 2 ) ℏ ω , n = 0 , 1 , 2 , … E_n = \left(n + \frac{1}{2}\right)\hbar\omega, \quad n=0,1,2,\ldots
E n = ( n + 2 1 ) ℏ ω , n = 0 , 1 , 2 , …
基态波函数 (n = 0 n=0n = 0 ):
ψ 0 ( x ) = ( m ω π ℏ ) 1 / 4 e − m ω x 2 / ( 2 ℏ ) \psi_0(x) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4} e^{-m\omega x^2/(2\hbar)}
ψ 0 ( x ) = ( π ℏ m ω ) 1 / 4 e − m ω x 2 / ( 2 ℏ )
数值示例 :一个质量为电子的谐振子,ω = 1 0 15 rad/s \omega = 10^{15}\,\text{rad/s}ω = 1 0 1 5 rad/s (约可见光频率):
n nn
E n E_nE n (eV)
波函数节点数
经典反转点 x c x_cx c (m)
0
0.329
0
± 3.04 × 1 0 − 11 \pm 3.04 \times 10^{-11}± 3 . 0 4 × 1 0 − 1 1
1
0.987
1
± 5.26 × 1 0 − 11 \pm 5.26 \times 10^{-11}± 5 . 2 6 × 1 0 − 1 1
2
1.645
2
± 6.79 × 1 0 − 11 \pm 6.79 \times 10^{-11}± 6 . 7 9 × 1 0 − 1 1
3
2.303
3
± 8.03 × 1 0 − 11 \pm 8.03 \times 10^{-11}± 8 . 0 3 × 1 0 − 1 1
零能点能量(n = 0 n=0n = 0 ) :E 0 = 1 2 ℏ ω ≠ 0 E_0 = \frac{1}{2}\hbar\omega \neq 0E 0 = 2 1 ℏ ω = 0 。这一"零点能"被实验证实——例如在绝对零度下,液氦仍然保持液态而不凝固,就是因为零点能抵消了原子间吸引力。即使在 T = 0 K T=0\,\text{K}T = 0 K ,量子粒子仍然在"抖动",无法完全静止。
在经典物理中,如果一个粒子的能量 E EE 低于势垒高度 V 0 V_0V 0 ,它将永远无法穿越势垒。但在量子力学中,由于波函数的波动性质,粒子有一定概率"穿"过势垒——这就是量子隧穿效应 。
考虑一个宽度为 a aa 、高度为 V 0 V_0V 0 的矩形势垒,粒子从左侧入射,能量 E < V 0 E < V_0E < V 0 。
透射系数 (近似公式,对于高势垒或宽势垒):
T ≈ 16 E ( V 0 − E ) V 0 2 e − 2 κ a , κ = 2 m ( V 0 − E ) ℏ T \approx \frac{16E(V_0 - E)}{V_0^2} e^{-2\kappa a}, \quad \kappa = \frac{\sqrt{2m(V_0 - E)}}{\hbar}
T ≈ V 0 2 1 6 E ( V 0 − E ) e − 2 κ a , κ = ℏ 2 m ( V 0 − E )
透射概率随势垒宽度的增加呈指数衰减。
数值例子 :电子(m = 9.11 × 1 0 − 31 kg m = 9.11 \times 10^{-31}\,\text{kg}m = 9 . 1 1 × 1 0 − 3 1 kg )穿过一个 V 0 = 5 eV V_0 = 5\,\text{eV}V 0 = 5 eV 的势垒,E = 3 eV E = 3\,\text{eV}E = 3 eV (低于势垒):
势垒宽度 a aa (nm)
κ \kappaκ (m− 1 ^{-1}− 1 )
e − 2 κ a e^{-2\kappa a}e − 2 κ a
透射概率 T TT
0.1
7.24 × 1 0 9 7.24 \times 10^97 . 2 4 × 1 0 9
2.34 × 1 0 − 1 2.34 \times 10^{-1}2 . 3 4 × 1 0 − 1
≈ 0.150 \approx 0.150≈ 0 . 1 5 0 (15%)
0.2
7.24 × 1 0 9 7.24 \times 10^97 . 2 4 × 1 0 9
5.49 × 1 0 − 2 5.49 \times 10^{-2}5 . 4 9 × 1 0 − 2
≈ 0.035 \approx 0.035≈ 0 . 0 3 5 (3.5%)
0.5
7.24 × 1 0 9 7.24 \times 10^97 . 2 4 × 1 0 9
7.10 × 1 0 − 4 7.10 \times 10^{-4}7 . 1 0 × 1 0 − 4
≈ 4.5 × 1 0 − 4 \approx 4.5 \times 10^{-4}≈ 4 . 5 × 1 0 − 4 (0.045%)
1.0
7.24 × 1 0 9 7.24 \times 10^97 . 2 4 × 1 0 9
5.04 × 1 0 − 7 5.04 \times 10^{-7}5 . 0 4 × 1 0 − 7
≈ 3.2 × 1 0 − 7 \approx 3.2 \times 10^{-7}≈ 3 . 2 × 1 0 − 7 (0.000032%)
随着势垒宽度增加到1nm,透射概率已经降到了亿分之几。但对于原子尺度的势垒(~0.1nm),隧穿概率可高达15%,这是扫描隧道显微镜(STM)工作的物理基础。
α衰变中的隧穿 :一个α粒子要在原子核中克服库仑势垒才能逃逸。典型的势垒参数:V 0 ≈ 25 MeV V_0 \approx 25\,\text{MeV}V 0 ≈ 2 5 MeV ,E ≈ 4 MeV E \approx 4\,\text{MeV}E ≈ 4 MeV ,α粒子的质量约为 6.64 × 1 0 − 27 kg 6.64 \times 10^{-27}\,\text{kg}6 . 6 4 × 1 0 − 2 7 kg 。计算得到的隧穿概率极低,但正是这个较小的概率决定了α衰变的半衰期在纳秒到数十亿年之间不等。
在量子力学中,每个物理可观测量(位置、动量、能量、角动量等)都对应一个线性厄米算符 。算符作用在波函数上,给出该物理量的可能取值和概率分布。
物理量
经典表达式
量子算符
作用
位置
r \mathbf{r}r
r ^ = r \hat{\mathbf{r}} = \mathbf{r}r ^ = r
乘以坐标:x ^ ψ = x ψ \hat{x}\psi = x\psix ^ ψ = x ψ
动量
p \mathbf{p}p
p ^ = − i ℏ ∇ \hat{\mathbf{p}} = -i\hbar\nablap ^ = − i ℏ ∇
微分:p ^ x ψ = − i ℏ ∂ ψ ∂ x \hat{p}_x\psi = -i\hbar\frac{\partial\psi}{\partial x}p ^ x ψ = − i ℏ ∂ x ∂ ψ
动能
p 2 / ( 2 m ) p^2/(2m)p 2 / ( 2 m )
T ^ = − ℏ 2 2 m ∇ 2 \hat{T} = -\frac{\hbar^2}{2m}\nabla^2T ^ = − 2 m ℏ 2 ∇ 2
二阶微分
势能
V ( r ) V(\mathbf{r})V ( r )
V ^ = V ( r ) \hat{V} = V(\mathbf{r})V ^ = V ( r )
乘以势能函数
能量(哈密顿)
p 2 / ( 2 m ) + V p^2/(2m) + Vp 2 / ( 2 m ) + V
H ^ = − ℏ 2 2 m ∇ 2 + V \hat{H} = -\frac{\hbar^2}{2m}\nabla^2 + VH ^ = − 2 m ℏ 2 ∇ 2 + V
—
角动量
r × p \mathbf{r} \times \mathbf{p}r × p
L ^ = − i ℏ ( r × ∇ ) \hat{\mathbf{L}} = -i\hbar(\mathbf{r} \times \nabla)L ^ = − i ℏ ( r × ∇ )
微分算符
电子波函数 ψ ( x ) = A e i p x / ℏ \psi(x) = Ae^{ipx/\hbar}ψ ( x ) = A e i p x / ℏ (德布罗意平面波):
p ^ ψ = − i ℏ d d x A e i p x / ℏ = − i ℏ ⋅ i p ℏ A e i p x / ℏ = p ψ \hat{p}\psi = -i\hbar\frac{d}{dx} Ae^{ipx/\hbar} = -i\hbar \cdot \frac{ip}{\hbar} Ae^{ipx/\hbar} = p\psi
p ^ ψ = − i ℏ d x d A e i p x / ℏ = − i ℏ ⋅ ℏ i p A e i p x / ℏ = p ψ
这说明德布罗意平面波正是动量算符的本征函数,本征值就是该粒子的动量 p pp 。
算符的乘法一般不满足交换律。两个算符 A ^ \hat{A}A ^ 和 B ^ \hat{B}B ^ 的对易子 定义为:
[ A ^ , B ^ ] = A ^ B ^ − B ^ A ^ [\hat{A}, \hat{B}] = \hat{A}\hat{B} - \hat{B}\hat{A}
[ A ^ , B ^ ] = A ^ B ^ − B ^ A ^
位置和动量算符之间的对易关系:
[ x ^ , p ^ x ] = i ℏ [\hat{x}, \hat{p}_x] = i\hbar
[ x ^ , p ^ x ] = i ℏ
[ x ^ , p ^ y ] = 0 [\hat{x}, \hat{p}_y] = 0
[ x ^ , p ^ y ] = 0
[ p ^ x , p ^ y ] = 0 [\hat{p}_x, \hat{p}_y] = 0
[ p ^ x , p ^ y ] = 0
验证 :将 [ x ^ , p ^ x ] [\hat{x}, \hat{p}_x][ x ^ , p ^ x ] 作用在一个测试波函数 ψ ( x ) \psi(x)ψ ( x ) 上:
[ x ^ , p ^ x ] ψ = − i ℏ [ x d ψ d x − d d x ( x ψ ) ] [\hat{x}, \hat{p}_x]\psi = -i\hbar\left[ x\frac{d\psi}{dx} - \frac{d}{dx}(x\psi) \right]
[ x ^ , p ^ x ] ψ = − i ℏ [ x d x d ψ − d x d ( x ψ ) ]
= − i ℏ [ x d ψ d x − ψ − x d ψ d x ] = i ℏ ψ = -i\hbar\left[ x\frac{d\psi}{dx} - \psi - x\frac{d\psi}{dx} \right] = i\hbar\psi
= − i ℏ [ x d x d ψ − ψ − x d x d ψ ] = i ℏ ψ
因此 [ x ^ , p ^ x ] = i ℏ [\hat{x}, \hat{p}_x] = i\hbar[ x ^ , p ^ x ] = i ℏ 成立。
这个非零的对易子就是海森堡不确定性原理的数学根源。
量子力学中的本征值方程为:
A ^ ψ n = a n ψ n \hat{A}\psi_n = a_n\psi_n
A ^ ψ n = a n ψ n
其中 a n a_na n 是本征值(可能的测量结果),ψ n \psi_nψ n 是本征函数。对于厄米算符:
本征值是实数(对应可观测量的测量值)
不同本征值对应的本征函数正交
所有本征函数构成一个完备集(任何波函数都可以展开为它们的线性组合)
举例:动量本征函数
p ^ ψ p ( x ) = p ψ p ( x ) \hat{p}\psi_p(x) = p\psi_p(x)
p ^ ψ p ( x ) = p ψ p ( x )
解为:ψ p ( x ) = 1 2 π ℏ e i p x / ℏ \psi_p(x) = \frac{1}{\sqrt{2\pi\hbar}}e^{ipx/\hbar}ψ p ( x ) = 2 π ℏ 1 e i p x / ℏ (平面波)
虽然平面波不是平方可积的(不能严格归一化),但它们是数学上的便利工具。任何实际的波函数都可以通过傅里叶变换表示为平面波的叠加:
ψ ( x ) = 1 2 π ℏ ∫ − ∞ ∞ ϕ ( p ) e i p x / ℏ d p \psi(x) = \frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty} \phi(p) e^{ipx/\hbar} dp
ψ ( x ) = 2 π ℏ 1 ∫ − ∞ ∞ ϕ ( p ) e i p x / ℏ d p
其中 ϕ ( p ) \phi(p)ϕ ( p ) 是动量空间的波函数。波函数在坐标空间和动量空间之间存在对偶关系,这与不确定性原理密切相关。
对于一个给定波函数 ψ \psiψ ,算符 A ^ \hat{A}A ^ 的期望值(统计平均)为:
⟨ A ⟩ = ∫ ψ ∗ A ^ ψ d 3 r \langle A \rangle = \int \psi^* \hat{A} \psi \, d^3r
⟨ A ⟩ = ∫ ψ ∗ A ^ ψ d 3 r
具体算例 :计算一维无限深势阱基态的位置期望值和动量期望值。
基态波函数:ψ 1 ( x ) = 2 / a sin ( π x / a ) \psi_1(x) = \sqrt{2/a}\sin(\pi x/a)ψ 1 ( x ) = 2 / a sin ( π x / a )
位置期望值:
⟨ x ⟩ = ∫ 0 a x ∣ ψ 1 ( x ) ∣ 2 d x = 2 a ∫ 0 a x 2 ( π x a ) d x = a 2 \langle x \rangle = \int_0^a x |\psi_1(x)|^2 dx = \frac{2}{a} \int_0^a x\sin^2\left(\frac{\pi x}{a}\right) dx = \frac{a}{2}
⟨ x ⟩ = ∫ 0 a x ∣ ψ 1 ( x ) ∣ 2 d x = a 2 ∫ 0 a x sin 2 ( a π x ) d x = 2 a
这正是势阱的中心位置,符合直觉。
位置平方期望值:
⟨ x 2 ⟩ = 2 a ∫ 0 a x 2 2 ( π x a ) d x = a 2 ( 1 3 − 1 2 π 2 ) \langle x^2 \rangle = \frac{2}{a} \int_0^a x^2\sin^2\left(\frac{\pi x}{a}\right) dx = a^2\left(\frac{1}{3} - \frac{1}{2\pi^2}\right)
⟨ x 2 ⟩ = a 2 ∫ 0 a x 2 sin 2 ( a π x ) d x = a 2 ( 3 1 − 2 π 2 1 )
不确定性:
Δ x = ⟨ x 2 ⟩ − ⟨ x ⟩ 2 = a 1 12 − 1 2 π 2 ≈ 0.181 a \Delta x = \sqrt{\langle x^2 \rangle - \langle x \rangle^2} = a\sqrt{\frac{1}{12} - \frac{1}{2\pi^2}} \approx 0.181a
Δ x = ⟨ x 2 ⟩ − ⟨ x ⟩ 2 = a 1 2 1 − 2 π 2 1 ≈ 0 . 1 8 1 a
对于 a = 0.1 nm a = 0.1\,\text{nm}a = 0 . 1 nm ,Δ x ≈ 0.0181 nm \Delta x \approx 0.0181\,\text{nm}Δ x ≈ 0 . 0 1 8 1 nm 。
动量期望值:⟨ p ⟩ = 0 \langle p \rangle = 0⟨ p ⟩ = 0 (由于对称性,粒子向左和向右的概率相等)。
动量平方期望值:
⟨ p 2 ⟩ = ∫ 0 a ψ 1 ∗ ( x ) ( − ℏ 2 d 2 d x 2 ) ψ 1 ( x ) d x = ℏ 2 π 2 a 2 \langle p^2 \rangle = \int_0^a \psi_1^*(x)\left(-\hbar^2\frac{d^2}{dx^2}\right)\psi_1(x) dx = \frac{\hbar^2\pi^2}{a^2}
⟨ p 2 ⟩ = ∫ 0 a ψ 1 ∗ ( x ) ( − ℏ 2 d x 2 d 2 ) ψ 1 ( x ) d x = a 2 ℏ 2 π 2
\Delta p = \sqrt\langle p^2 \rangle - \langle p \rangle^2 = \frac\hbar\pia \approx \frac1.05\times 10^-34 \times 3.140.1\times 10^-9 \approx 3.30 \times 10^-24\,\textkg·m/s
验证不确定性关系:
Δ x ⋅ Δ p ≈ 0.181 a × ℏ π a = 0.181 π ℏ ≈ 0.569 ℏ > ℏ 2 \Delta x \cdot \Delta p \approx 0.181a \times \frac{\hbar\pi}{a} = 0.181\pi\hbar \approx 0.569\hbar > \frac{\hbar}{2}
Δ x ⋅ Δ p ≈ 0 . 1 8 1 a × a ℏ π = 0 . 1 8 1 π ℏ ≈ 0 . 5 6 9 ℏ > 2 ℏ
满足海森堡不确定性原理 Δ x ⋅ Δ p ≥ ℏ / 2 \Delta x \cdot \Delta p \geq \hbar/2Δ x ⋅ Δ p ≥ ℏ / 2 。
不确定性原理是量子力学最著名的结果之一,由海森堡在1927年提出。它揭示了量子世界中一个根本性的限制:某些成对的物理量(如位置和动量、能量和时间)不能同时被精确测量。
对于任何两个不对易的厄米算符 A ^ \hat{A}A ^ 和 B ^ \hat{B}B ^ ,有:
Δ A ⋅ Δ B ≥ 1 2 ∣ ⟨ [ A ^ , B ^ ] ⟩ ∣ \Delta A \cdot \Delta B \geq \frac{1}{2} \left| \langle [\hat{A}, \hat{B}] \rangle \right|
Δ A ⋅ Δ B ≥ 2 1 ∣ ∣ ∣ ∣ ⟨ [ A ^ , B ^ ] ⟩ ∣ ∣ ∣ ∣
其中 Δ A = ⟨ A ^ 2 ⟩ − ⟨ A ^ ⟩ 2 \Delta A = \sqrt{\langle \hat{A}^2 \rangle - \langle \hat{A} \rangle^2}Δ A = ⟨ A ^ 2 ⟩ − ⟨ A ^ ⟩ 2 是标准差。
对于位置和动量:
Δ x ⋅ Δ p x ≥ ℏ 2 \Delta x \cdot \Delta p_x \geq \frac{\hbar}{2}
Δ x ⋅ Δ p x ≥ 2 ℏ
不确定性原理不是 说测量仪器不够精确——它是自然界的根本极限,与测量技术无关。即使我们有完美的测量仪器,也无法同时精确知道一个粒子的位置和动量。
思想实验 :要"看到"一个电子,你必须用光子照射它。光子能量越高(波长越短),定位越精确——但高能光子的动量 p = h / λ p = h/\lambdap = h / λ 越大,会显著扰动电子的动量。
期望分辨率
所需波长
光子动量
对电子的动量扰动
原子尺度(1 0 − 10 m 10^{-10}\,\text{m}1 0 − 1 0 m )
X射线(1 0 − 10 m 10^{-10}\,\text{m}1 0 − 1 0 m )
6.6 × 1 0 − 24 6.6 \times 10^{-24}6 . 6 × 1 0 − 2 4
与电子动量相当
核尺度(1 0 − 14 m 10^{-14}\,\text{m}1 0 − 1 4 m )
伽马射线(1 0 − 14 m 10^{-14}\,\text{m}1 0 − 1 4 m )
6.6 × 1 0 − 20 6.6 \times 10^{-20}6 . 6 × 1 0 − 2 0
远大于电子动量
原子核内部(1 0 − 15 m 10^{-15}\,\text{m}1 0 − 1 5 m )
高能伽马
6.6 × 1 0 − 19 6.6 \times 10^{-19}6 . 6 × 1 0 − 1 9
完全改变电子状态
典型数值应用 :
原子中电子的基态能量估算(无需解薛定谔方程) :
设电子被束缚在大小为 a aa 的区域中,则 Δ x ≈ a \Delta x \approx aΔ x ≈ a ,因此 Δ p ≈ ℏ / a \Delta p \approx \hbar/aΔ p ≈ ℏ / a 。总能量:
E ≈ ( Δ p ) 2 2 m − e 2 4 π ε 0 a ≈ ℏ 2 2 m a 2 − e 2 4 π ε 0 a E \approx \frac{(\Delta p)^2}{2m} - \frac{e^2}{4\pi\varepsilon_0 a} \approx \frac{\hbar^2}{2ma^2} - \frac{e^2}{4\pi\varepsilon_0 a}
E ≈ 2 m ( Δ p ) 2 − 4 π ε 0 a e 2 ≈ 2 m a 2 ℏ 2 − 4 π ε 0 a e 2
求 d E / d a = 0 dE/da = 0d E / d a = 0 得到最优 a aa :
ℏ 2 m a 3 = e 2 4 π ε 0 a 2 ⇒ a = 4 π ε 0 ℏ 2 m e 2 ≈ 0.053 nm \frac{\hbar^2}{ma^3} = \frac{e^2}{4\pi\varepsilon_0 a^2} \Rightarrow a = \frac{4\pi\varepsilon_0\hbar^2}{me^2} \approx 0.053\,\text{nm}
m a 3 ℏ 2 = 4 π ε 0 a 2 e 2 ⇒ a = m e 2 4 π ε 0 ℏ 2 ≈ 0 . 0 5 3 nm
这正是玻尔半径!代入得最低能量约为 − 13.6 eV -13.6\,\text{eV}− 1 3 . 6 eV ——完全正确。不确定性原理直接给出了氢原子基态的能量尺度。
原子核中不可能存在电子 :核力作用范围约 1 0 − 15 m 10^{-15}\,\text{m}1 0 − 1 5 m 。如果电子被限制在这个区域内:
\Delta p \approx \frac\hbar2\Delta x \approx \frac1.05\times 10^-342\times 10^-15 \approx 5.3 \times 10^-20\,\textkg·m/s
对应的动能 E ≈ ( Δ p ) 2 / ( 2 m e ) ≈ 9.6 × 1 0 − 10 J ≈ 6 GeV E \approx (\Delta p)^2/(2m_e) \approx 9.6\times 10^{-10}\,\text{J} \approx 6\,\text{GeV}E ≈ ( Δ p ) 2 / ( 2 m e ) ≈ 9 . 6 × 1 0 − 1 0 J ≈ 6 GeV 。而核反应放出的能量约在MeV量级——电子如果被禁闭在原子核内,其零点能就高达GeV量级,远超核反应的能量尺度。这就是为什么核物理中考虑的是质子、中子,而不是电子。
Δ E ⋅ Δ t ≥ ℏ 2 \Delta E \cdot \Delta t \geq \frac{\hbar}{2}
Δ E ⋅ Δ t ≥ 2 ℏ
这里 Δ t \Delta tΔ t 的解释较为微妙。它通常理解为系统从一个状态演化到可分辨的另一状态所需的时间。这导致:
虚粒子与相互作用 :允许在短时间内 Δ t ∼ ℏ / ( 2 Δ E ) \Delta t \sim \hbar/(2\Delta E)Δ t ∼ ℏ / ( 2 Δ E ) 内违反能量守恒(如费曼图中的虚粒子)
能级宽度 :激发态寿命 τ \tauτ 与能级宽度 Γ \GammaΓ 的关系:Γ τ ∼ ℏ \Gamma\tau \sim \hbar
Γ τ ∼ ℏ
激发态
寿命 τ \tauτ
能级宽度 Γ \GammaΓ
观测
氢原子2p态
1.6 × 1 0 − 9 s 1.6 \times 10^{-9}\,\text{s}1 . 6 × 1 0 − 9 s
4.1 × 1 0 − 7 eV 4.1 \times 10^{-7}\,\text{eV}4 . 1 × 1 0 − 7 eV
极窄的谱线
钠D线
1.6 × 1 0 − 8 s 1.6 \times 10^{-8}\,\text{s}1 . 6 × 1 0 − 8 s
4.1 × 1 0 − 8 eV 4.1 \times 10^{-8}\,\text{eV}4 . 1 × 1 0 − 8 eV
很窄
Δ \DeltaΔ 共振态(粒子物理)
5.6 × 1 0 − 24 s 5.6 \times 10^{-24}\,\text{s}5 . 6 × 1 0 − 2 4 s
117 MeV 117\,\text{MeV}1 1 7 MeV
很宽的共振峰
态叠加原理是量子力学的核心原则之一:如果 ψ 1 \psi_1ψ 1 和 ψ 2 \psi_2ψ 2 是系统的两个可能状态,那么它们的任意线性组合 Ψ = c 1 ψ 1 + c 2 ψ 2 \Psi = c_1\psi_1 + c_2\psi_2Ψ = c 1 ψ 1 + c 2 ψ 2 也是一个可能状态(其中 ∣ c 1 ∣ 2 + ∣ c 2 ∣ 2 = 1 |c_1|^2 + |c_2|^2 = 1∣ c 1 ∣ 2 + ∣ c 2 ∣ 2 = 1 )。
这与经典叠加有本质区别。在经典物理中,叠加意味着"混合"——比如50%黑墨水和50%红墨水,每个墨水分子要么是黑的要么是红的。但在量子叠加中,一个粒子可以同时 处于两个状态之中,直到测量时"坍缩"到其中一个。
薛定谔的猫 是一个著名的思想实验,用来展示量子叠加对宏观世界的含义:盒子里有一只猫、一瓶毒药和放射性物质。如果原子衰变(概率50%),毒药释放,猫死;如果没有衰变,猫活。在打开盒子之前,猫处于"死"和"活"的量子叠加态。这个思想实验并不是说猫真的既死又活,而是揭示了量子力学在宏观尺度上诠释的困境。
如果两个或更多粒子以某种方式相互作用,它们的量子态可能变得"纠缠"(Entanglement)——即不能独立地描述每个粒子的状态,只能描述它们整体的状态。
贝尔态 (最大纠缠态):
∣ Φ + ⟩ = 1 2 ( ∣ 00 ⟩ + ∣ 11 ⟩ ) |\Phi^+\rangle = \frac{1}{\sqrt{2}}(|00\rangle + |11\rangle)
∣ Φ + ⟩ = 2 1 ( ∣ 0 0 ⟩ + ∣ 1 1 ⟩ )
∣ Φ − ⟩ = 1 2 ( ∣ 00 ⟩ − ∣ 11 ⟩ ) |\Phi^-\rangle = \frac{1}{\sqrt{2}}(|00\rangle - |11\rangle)
∣ Φ − ⟩ = 2 1 ( ∣ 0 0 ⟩ − ∣ 1 1 ⟩ )
∣ Ψ + ⟩ = 1 2 ( ∣ 01 ⟩ + ∣ 10 ⟩ ) |\Psi^+\rangle = \frac{1}{\sqrt{2}}(|01\rangle + |10\rangle)
∣ Ψ + ⟩ = 2 1 ( ∣ 0 1 ⟩ + ∣ 1 0 ⟩ )
∣ Ψ − ⟩ = 1 2 ( ∣ 01 ⟩ − ∣ 10 ⟩ ) |\Psi^-\rangle = \frac{1}{\sqrt{2}}(|01\rangle - |10\rangle)
∣ Ψ − ⟩ = 2 1 ( ∣ 0 1 ⟩ − ∣ 1 0 ⟩ )
在 ∣ Φ + ⟩ |\Phi^+\rangle∣ Φ + ⟩ 态中,如果测量第一个粒子的自旋为自旋向上,则第二个粒子的自旋必然也是向上;如果第一个被测量为向下,则第二个也是向下。这种关联是即时的(但不能 用于超光速通信,因为测量结果是随机的)。
纠缠的量化对比 :
性质
经典关联
量子纠缠
关联来源
共同历史或因果关系
不可分离的量子态
测量结果
联合概率 P ( a , b ) P(a,b)P ( a , b ) 可分解为 P ( a ) P ( b ∣ a ) P(a)P(b|a)P ( a ) P ( b ∣ a )
联合概率不可分解
远程影响
需要交换信息
无需信息传递(但无因果信号)
空间分离后的行为
互不影响
测量一个立即影响另一个
数学描述
概率论
张量积空间中的矢量
爱因斯坦认为量子纠缠意味着量子力学是不完备的,他称之为"鬼魅般的超距作用"。1935年,他与波多尔斯基、罗森共同提出了EPR悖论,认为存在"隐变量"可以解释量子关联。
1964年,约翰·贝尔推导出贝尔不等式 :如果存在隐变量理论,某些测量的相关性必须满足以下不等式之一:
CHSH不等式 :∣ E ( a , b ) + E ( a , b ′ ) + E ( a ′ , b ) − E ( a ′ , b ′ ) ∣ ≤ 2 |E(a,b) + E(a,b') + E(a',b) - E(a',b')| \leq 2∣ E ( a , b ) + E ( a , b ′ ) + E ( a ′ , b ) − E ( a ′ , b ′ ) ∣ ≤ 2
其中 E ( a , b ) E(a,b)E ( a , b ) 是在方向 a aa 和 b bb 上测量自旋的期望值。
实验结果 :从1980年代到2015年,多个独立实验进行了贝尔不等式的检验:
实验
年份
测得的CHSH值
量子力学预测
不等式上限
结论
Aspect等
1982
2.70 ± \pm± 0.05
2 2 ≈ 2.828 2\sqrt{2} \approx 2.8282 2 ≈ 2 . 8 2 8
2
违背贝尔不等式
Weihs等
1998
2.75 ± \pm± 0.09
2.828
2
明确违背
Zeilinger(Hensen)
2015
2.42 ± \pm± 0.20
2.828
2
无漏洞检验
Giustina
2015
2.38 ± \pm± 0.10
2.828
2
光子无漏洞
所有实验结果都违背了贝尔不等式,与量子力学的预测一致——量子纠缠确实存在,不是隐变量可以解释的。量子力学对物理实在的描述比经典的局部实在论更准确。
量子测量的标准描述由冯·诺依曼提出,包含两个截然不同的过程:
幺正演化(U过程) :由薛定谔方程描述,波函数按确定性的、可逆的方式演化
波函数坍缩(R过程) :当进行测量时,波函数不可预测地、"非确定性"地跳转到一个本征态
测量过程图示 :
初始态 Ψ = c₁ψ₁ + c₂ψ₂
↓ (进行测量)
┌──────┴──────┐
↓ ↓
概率 |c₁|² 概率 |c₂|²
↓ ↓
本征态 ψ₁ 本征态 ψ₂
测得值 a₁ 测得值 a₂
关键的诠释问题 :
什么是"测量"?在哪个点上坍缩发生?
测量仪器也是由量子粒子组成的,为什么它不会处于叠加态?
是否存在"量子-经典边界"?
这些问题至今没有统一的答案,而是形成了不同的诠释流派。
诠释
坍缩机制
波函数的本体论地位
关键争议点
哥本哈根诠释
测量导致坍缩,不追问其机制
工具主义——波函数是计算工具
没有定义"测量"的边界
多世界诠释
没有坍缩,每次测量分支为平行世界
波函数是唯一的物理实在
世界数目的"费用"问题
退相干理论
与环境纠缠导致"表现"坍缩
波函数是物理的
需要补充选择规则
德布罗意-玻姆导航波
隐藏的确定性轨迹
波函数是真实场 + 粒子
非局域性与相对论兼容性
客观坍缩模型
自发随机坍缩(如GRW模型)
波函数是物理的
需要修改薛定谔方程
退相干——"表现"坍缩的物理机制 :
退相干是最被广泛接受的解释:一个量子系统几乎不可能完全与环境隔离。系统与环境粒子相互作用后,它们的量子态纠缠在一起。当你测量系统时,实际上也需要测量环境——这导致不同分支之间的相干性被"稀释"到海量的环境自由度中,使得干涉效应无法被观测到。
退相干时间尺度举例:
系统
环境耦合
退相干时间
悬浮冷原子
极低(磁光阱)
> 1 s > 1\,\text{s}> 1 s
离子阱中的离子
极低
> 100 ms > 100\,\text{ms}> 1 0 0 ms
单个分子(室温液体)
中等(溶剂碰撞)
∼ 1 0 − 12 s \sim 10^{-12}\,\text{s}∼ 1 0 − 1 2 s
宏观物体(1μm灰尘)
强(空气分子)
∼ 1 0 − 20 s \sim 10^{-20}\,\text{s}∼ 1 0 − 2 0 s
宏观物体的退相干时间极短——这就是为什么我们看不到桌子的波动性,也看不到猫的死活叠加态。
当系统处于纯态 Ψ = ∑ n c n ψ n \Psi = \sum_n c_n \psi_nΨ = ∑ n c n ψ n 时,密度算符为:
ρ ^ = ∣ Ψ ⟩ ⟨ Ψ ∣ \hat{\rho} = |\Psi\rangle\langle\Psi|
ρ ^ = ∣ Ψ ⟩ ⟨ Ψ ∣
对于具体的纯态 Ψ = 1 2 ( ∣ 0 ⟩ + ∣ 1 ⟩ ) \Psi = \frac{1}{\sqrt{2}}(|0\rangle + |1\rangle)Ψ = 2 1 ( ∣ 0 ⟩ + ∣ 1 ⟩ ) :
ρ 纯 = 1 2 ( 1 1 1 1 ) \rho_{\text{纯}} = \frac{1}{2}\begin{pmatrix}1 & 1\\ 1 & 1\end{pmatrix}
ρ 纯 = 2 1 ( 1 1 1 1 )
非对角元 (ρ 01 = ρ 10 = 1 / 2 \rho_{01} = \rho_{10} = 1/2ρ 0 1 = ρ 1 0 = 1 / 2 ) 代表量子相干性。
当系统与环境纠缠导致退相干后,密度矩阵变为混合态——非对角元消失:
ρ 混合 = 1 2 ( 1 0 0 1 ) \rho_{\text{混合}} = \frac{1}{2}\begin{pmatrix}1 & 0\\ 0 & 1\end{pmatrix}
ρ 混合 = 2 1 ( 1 0 0 1 )
这代表50%概率在 ∣ 0 ⟩ |0\rangle∣ 0 ⟩ 、50%概率在 ∣ 1 ⟩ |1\rangle∣ 1 ⟩ 的经典统计混合——不再有量子相干性。
量子力学中,角动量是一个厄米算符,满足对易关系:
[ L ^ x , L ^ y ] = i ℏ L ^ z [\hat{L}_x, \hat{L}_y] = i\hbar\hat{L}_z
[ L ^ x , L ^ y ] = i ℏ L ^ z
(以及循环置换)
角动量平方算符 L ^ 2 = L ^ x 2 + L ^ y 2 + L ^ z 2 \hat{L}^2 = \hat{L}_x^2 + \hat{L}_y^2 + \hat{L}_z^2L ^ 2 = L ^ x 2 + L ^ y 2 + L ^ z 2 与任何分量对易:[ L ^ 2 , L ^ z ] = 0 [\hat{L}^2, \hat{L}_z] = 0[ L ^ 2 , L ^ z ] = 0 。
这意味着 L ^ 2 \hat{L}^2L ^ 2 和 L ^ z \hat{L}_zL ^ z 可以有共同的本征态。本征值方程:
L ^ 2 ∣ l , m ⟩ = ℏ 2 l ( l + 1 ) ∣ l , m ⟩ , l = 0 , 1 , 2 , … \hat{L}^2 |l, m\rangle = \hbar^2 l(l+1) |l, m\rangle, \quad l = 0, 1, 2, \ldots
L ^ 2 ∣ l , m ⟩ = ℏ 2 l ( l + 1 ) ∣ l , m ⟩ , l = 0 , 1 , 2 , …
L ^ z ∣ l , m ⟩ = ℏ m ∣ l , m ⟩ , m = − l , − l + 1 , … , l − 1 , l \hat{L}_z |l, m\rangle = \hbar m |l, m\rangle, \quad m = -l, -l+1, \ldots, l-1, l
L ^ z ∣ l , m ⟩ = ℏ m ∣ l , m ⟩ , m = − l , − l + 1 , … , l − 1 , l
球谐函数 是轨道角动量的本征函数:
Y l m ( θ , ϕ ) = ( − 1 ) m 2 l + 1 4 π ( l − m ) ! ( l + m ) ! P l m ( cos θ ) e i m ϕ Y_l^m(\theta, \phi) = (-1)^m \sqrt{\frac{2l+1}{4\pi}\frac{(l-m)!}{(l+m)!}} P_l^m(\cos\theta) e^{im\phi}
Y l m ( θ , ϕ ) = ( − 1 ) m 4 π 2 l + 1 ( l + m ) ! ( l − m ) ! P l m ( cos θ ) e i m ϕ
轨道角动量的可视化 :
量子数 l ll
通称
光谱记号
磁量子数范围
简并度
0
s轨道
l = 0 l=0l = 0
m = 0 m=0m = 0
1
1
p轨道
l = 1 l=1l = 1
m = − 1 , 0 , 1 m=-1,0,1m = − 1 , 0 , 1
3
2
d轨道
l = 2 l=2l = 2
m = − 2 , − 1 , 0 , 1 , 2 m=-2,-1,0,1,2m = − 2 , − 1 , 0 , 1 , 2
5
3
f轨道
l = 3 l=3l = 3
m = − 3 , − 2 , − 1 , 0 , 1 , 2 , 3 m=-3,-2,-1,0,1,2,3m = − 3 , − 2 , − 1 , 0 , 1 , 2 , 3
7
角动量本征态的几何图像 :角动量方向不是任意的,而是量子化的——L z L_zL z 只能取 m ℏ m\hbarm ℏ 的离散值。对于一个给定的 l ll ,角动量矢量 L \mathbf{L}L 在以z轴为中心的圆锥上旋转,使得它的z分量是确定的,但x和y分量完全不确定。
自旋是粒子的一种内禀角动量,与空间运动无关。它没有经典对应物。
自旋算符满足与轨道角动量相同的对易关系 :
[ S ^ x , S ^ y ] = i ℏ S ^ z [\hat{S}_x, \hat{S}_y] = i\hbar\hat{S}_z
[ S ^ x , S ^ y ] = i ℏ S ^ z
但对于自旋 s ss 的粒子,s ss 可以是半整数。
粒子
自旋 s ss
S z S_zS z 可能取值
统计性质
电子
1/2
± ℏ / 2 \pm\hbar/2± ℏ / 2
费米子(费米-狄拉克统计)
质子
1/2
± ℏ / 2 \pm\hbar/2± ℏ / 2
费米子
光子
1
− ℏ , 0 , + ℏ -\hbar, 0, +\hbar− ℏ , 0 , + ℏ
玻色子(玻色-爱因斯坦统计)
π介子
0
0
玻色子
引力子(理论)
2
− 2 ℏ , − ℏ , 0 , + ℏ , + 2 ℏ -2\hbar, -\hbar, 0, +\hbar, +2\hbar− 2 ℏ , − ℏ , 0 , + ℏ , + 2 ℏ
玻色子
电子自旋可以用泡利矩阵完全描述。令 S ^ i = ℏ 2 σ i \hat{S}_i = \frac{\hbar}{2}\sigma_iS ^ i = 2 ℏ σ i :
σ x = ( 0 1 1 0 ) , σ y = ( 0 − i i 0 ) , σ z = ( 1 0 0 − 1 ) \sigma_x = \begin{pmatrix}0 & 1\\ 1 & 0\end{pmatrix}, \quad \sigma_y = \begin{pmatrix}0 & -i\\ i & 0\end{pmatrix}, \quad \sigma_z = \begin{pmatrix}1 & 0\\ 0 & -1\end{pmatrix}
σ x = ( 0 1 1 0 ) , σ y = ( 0 i − i 0 ) , σ z = ( 1 0 0 − 1 )
自旋上态 ∣ ↑ ⟩ = ( 1 0 ) |\uparrow\rangle = \begin{pmatrix}1\\0\end{pmatrix}∣ ↑ ⟩ = ( 1 0 ) 和自旋下态 ∣ ↓ ⟩ = ( 0 1 ) |\downarrow\rangle = \begin{pmatrix}0\\1\end{pmatrix}∣ ↓ ⟩ = ( 0 1 ) 是 S ^ z \hat{S}_zS ^ z 的本征态。
施特恩-格拉赫实验 是自旋存在的最直接证据。当银原子束通过非均匀磁场时,经典物理预期原子束会连续散开(因为磁矩方向在空间中连续变化),但实验结果表明原子束只分裂为两束——对应自旋向上和自旋向下。
经典计算机的信息基本单位是比特(bit),取值0或1。量子计算机的信息基本单位是量子比特(qubit),它可以处于 ∣ 0 ⟩ |0\rangle∣ 0 ⟩ 和 ∣ 1 ⟩ |1\rangle∣ 1 ⟩ 的叠加态:
∣ ψ ⟩ = α ∣ 0 ⟩ + β ∣ 1 ⟩ , ∣ α ∣ 2 + ∣ β ∣ 2 = 1 |\psi\rangle = \alpha|0\rangle + \beta|1\rangle, \quad |\alpha|^2 + |\beta|^2 = 1
∣ ψ ⟩ = α ∣ 0 ⟩ + β ∣ 1 ⟩ , ∣ α ∣ 2 + ∣ β ∣ 2 = 1
量子比特与经典比特的对比 :
性质
经典比特
量子比特
状态数
1个状态(0或1)
无限多个叠加态
读取
可直接读
测量后坍缩到0或1
复制
可以任意复制
不可克隆定理——不能精确复制
存储
独立存在
可纠缠(非局域关联)
物理实现
晶体管电压
光子偏振、电子自旋、超导电路
量子门操作是幺正的(保持概率守恒),对应矩阵运算:
门
符号
矩阵表示
作用
Hadamard
H HH
1 2 ( 1 1 1 − 1 ) \frac{1}{\sqrt{2}}\begin{pmatrix}1&1\\1&-1\end{pmatrix}2 1 ( 1 1 1 − 1 )
∣ 0 ⟩ → ∣ 0 ⟩ + ∣ 1 ⟩ 2 |0\rangle \to \frac{|0\rangle+|1\rangle}{\sqrt{2}}∣ 0 ⟩ → 2 ∣ 0 ⟩ + ∣ 1 ⟩ 创建叠加态
Pauli-X
X XX
( 0 1 1 0 ) \begin{pmatrix}0&1\\1&0\end{pmatrix}( 0 1 1 0 )
量子版的NOT门
CNOT
C X CXC X
( 1 0 0 0 0 1 0 0 0 0 0 1 0 0 1 0 ) \begin{pmatrix}1&0&0&0\\0&1&0&0\\0&0&0&1\\0&0&1&0\end{pmatrix}⎝ ⎜ ⎜ ⎜ ⎛ 1 0 0 0 0 1 0 0 0 0 0 1 0 0 1 0 ⎠ ⎟ ⎟ ⎟ ⎞
受控非门,纠缠量子比特
Toffoli
C C X CCXC C X
8x8矩阵
通用量子门
Deutsch-Jozsa算法 展示了量子计算的优势。问题:给定一个未知函数 f : { 0 , 1 } n → { 0 , 1 } f: \{0,1\}^n \to \{0,1\}f : { 0 , 1 } n → { 0 , 1 } ,已知它要么是常函数(输出恒定),要么是平衡函数(一半输出0一半输出1)。经典算法在最坏情况下需要 2 n − 1 + 1 2^{n-1}+12 n − 1 + 1 次查询,而量子算法只需1次查询。
方案
物理系统
量子比特数(2024年)
退相干时间
主要挑战
超导电子学
LC谐振器、约瑟夫森结
1000+(IBM、Google)
10-200μs
需要极低温(~15mK)
离子阱
囚禁在电磁场中的离子
50-100(IonQ、Quantinuum)
秒量级
门操作较慢
拓扑量子计算
任意子(Majorana粒子)
实验验证阶段
理论上很长
实现难度极大
光量子
光子路径编码
几百(USTC、PsiQuantum)
ms-s
光子损耗
硅自旋量子比特
硅中的电子自旋
10-50
数百μs
制造精度要求高
截至2024年,量子计算机还处于含噪中等规模量子(NISQ)阶段——量子比特数在几十到一千多之间,但错误率还不足以运行纠错码和大规模算法。Google在2023年实现了"超越经典"的随机电路采样,但实际有实用价值的量子计算机(如Shor算法破解RSA密码系统)还需要 10⁵-10⁶ 个逻辑量子比特。
量子力学是半导体物理学的基础,也是现代电子产业的根基。
PN结的量子解释 :
电子和空穴的能带结构决定了半导体行为
量子隧穿效应在隧道二极管中是最基本的工作机制
量子阱和超晶格结构用于LED和激光器中
摩尔定律的量子限制 :
当前最先进芯片制程为3nm(约含22个硅原子宽度)
当栅极长度接近量子隧穿极限(~2nm)时,电子会不可控地隧穿过栅极——这是摩尔定律放缓的根本物理原因
量子效应在纳米尺度器件中不再是忽略不计的修正,而是主导物理
量子力学在化学计算中的应用 :
H Ψ = E Ψ (需求解多电子系统的薛定谔方程) H\Psi = E\Psi \text{(需求解多电子系统的薛定谔方程)}
H Ψ = E Ψ (需求解多电子系统的薛定谔方程)
精确解仅对氢原子可得。对于多电子系统,使用近似方法:
方法
缩写
计算复杂度
适用范围
精度
哈特里-福克
HF
O ( N 4 ) O(N^4)O ( N 4 )
小分子(<100原子)
~1eV误差
密度泛函理论
DFT
O ( N 3 ) O(N^3)O ( N 3 )
大分子、固体(~1000原子)
~0.3eV误差
耦合簇方法
CCSD(T)
O ( N 7 ) O(N^7)O ( N 7 )
极小分子(<20原子)
化学精度
量子蒙特卡洛
QMC
O ( N 3 N sample ) O(N^3N_\text{sample})O ( N 3 N sample )
中等体系
高精度,可标度好
量子力学为新材料设计提供了理论预测能力 :
高温超导体的微观机制仍在探索中,但量子计算模拟可能提供关键突破
光伏材料的能带工程完全基于量子力学的能带理论
量子点用于生物成像和量子计算
利用海森堡不确定性原理和量子不可克隆定理,实现了理论上的无条件安全通信。
BB84协议流程 :
Alice随机选择基({ + , × } \{+, \times\}{ + , × } )和比特值(0或1),制备并发送光子
Bob随机选择基测量接收到的光子
Alice和Bob公开比较基(不公开比特值),只保留基匹配的位
随机抽取部分比对错误率,检查是否有窃听
剩余比特作为安全密钥
安全性保证 :如果存在窃听者Eve,她必须测量光子(因为不可克隆),这会导致量子态扰动——根据不确定性原理,她无法在不留下痕迹的情况下获取信息。检测到的错误率 > 11 % > 11\%> 1 1 % 即表明存在窃听。
实际QKD系统性能(2024年) :
系统
传输距离
密钥率
信道类型
光纤QKD
50km
~1 Mbps
光纤
光纤QKD
100km
~100 kbps
光纤
光纤QKD
420km(中微子通信)
~1 bps
超低损耗光纤
卫星QKD(墨子号)
1200km
~10 kbps
自由空间
STM利用量子隧穿效应直接观察和操纵原子表面。当探针针尖与样品表面距离在1nm以内时,电子可以通过隧穿形成电流。隧穿电流对距离极其敏感——距离每增加0.1nm,电流下降约一个数量级。
I ∝ V ⋅ e − 2 κ d , κ = 2 m ϕ ℏ I \propto V \cdot e^{-2\kappa d}, \quad \kappa = \frac{\sqrt{2m\phi}}{\hbar}
I ∝ V ⋅ e − 2 κ d , κ = ℏ 2 m ϕ
其中 ϕ \phiϕ 是功函数,d dd 是探针-样品距离。
典型参数:ϕ ≈ 4 eV \phi \approx 4\,\text{eV}ϕ ≈ 4 eV ,d ≈ 0.5 nm d \approx 0.5\,\text{nm}d ≈ 0 . 5 nm ,V = 0.1 V V = 0.1\,\text{V}V = 0 . 1 V 时:
I ≈ 1 nA × e − 2 × 1.02 × 1 0 10 × 0.5 × 1 0 − 9 ≈ 0.036 nA I \approx 1\,\text{nA} \times e^{-2 \times 1.02\times 10^{10} \times 0.5\times 10^{-9}} \approx 0.036\,\text{nA}
I ≈ 1 nA × e − 2 × 1 . 0 2 × 1 0 1 0 × 0 . 5 × 1 0 − 9 ≈ 0 . 0 3 6 nA
通过反馈控制保持电流恒定,可以精确追踪原子表面的形貌,达到 0.01nm 的垂直分辨率。STM不仅用于成像,还可以精确操作单个原子——1990年IBM用35个氙原子拼出了"IBM"字样,标志着原子尺度操纵的开端。
Griffiths, D. J. (2018). Introduction to Quantum Mechanics (3rd ed.). Cambridge University Press. —— 最受本科生欢迎的量子力学教材
Sakurai, J. J., & Napolitano, J. (2017). Modern Quantum Mechanics (2nd ed.). Cambridge University Press. —— 研究生级别标准教材
Feynman, R. P. (2011). The Feynman Lectures on Physics, Vol. III: Quantum Mechanics . Basic Books. —— 费曼的直觉式讲解,经典之作
Nielsen, M. A., & Chuang, I. L. (2010). Quantum Computation and Quantum Information (10th ed.). Cambridge University Press. —— 量子信息领域圣经
Dirac, P. A. M. (1981). The Principles of Quantum Mechanics (4th ed.). Oxford University Press. —— 量子力学奠基人的经典论述
Aspect, A., Dalibard, J., & Roger, G. (1982). Experimental test of Bell's inequalities using time-varying analyzers. Physical Review Letters , 49(25), 1804.
Hensen, B., et al. (2015). Loophole-free Bell inequality violation using electron spins separated by 1.3 km. Nature , 526, 682-686.
Arute, F., et al. (2019). Quantum supremacy using a programmable superconducting processor. Nature , 574, 505-510.