量子力学

Size: px

Start display at page:

Download "量子力学"

娇马
7 years ago
Views:

1 第章中心势场我们现在将一维量子问题推广到三维中心势场中心势场 (V ( r )=V (r)) 具有空间旋转对称性, 导致哈密顿算符 Ĥ 不含有空间方位角度信息, 从而有 [ Ĥ, ˆ L ] = 0, 即轨道角动量是守恒量因为 z 轴方位的选取并不改变物理观测量, 表现为哈密顿算符不包含 ˆL z, 此时物理体系能量存在简并, 简并度为 l + 1( l m l) 由于对称性的关系, 在球坐标中处理中心势场的问题更加便利, 定义径向动量算符 (x,y,z) (r sinθ cosφ,r sinθ sinφ,r cosθ). 可以验证 ˆp r 是厄米算符, 而且满足下面的对易关系 ( ˆp r = i h r + 1 ), r [r, ˆp r ] = i h. 因为 ˆp r = h ( r + r ), r 所以哈密顿算符可以写作如下形式 : Ĥ = ˆp r ˆ L µ + + V (r). µr

2 第章中心势场轨道角动量中心势场的一般性质球方势阱两体相互作用物理学描述各种物体之间相互作用规律的一门科学与某个物体的具体属性相比, 我们更加关心的是相互作用形式, 因为相互作用才是各种物体之间共有的性质迄今为止, 我们在自然界中已经发现存在四种基本相互作用 : 引力电磁相互作用弱相互作用和强相互作用 ; 当然这些基本相互作用在不同能量标度上表现为其他的有效形式在前面课程中我们讨论了各种势场中的一维量子问题这些势场都来源于物体之间的相互作用, 下面我们研究一下量子世界中的最简单的相互作用两体之间的相互作用首先需要了解的是如何描述两体之间的相互作用经典物理中我们通常将两体问题分解为整体运动和两个物体之间的相对运动量子力学中所有物理量都是用算符表示, 我们是否还可以像经典物理一样处理两体问题哪? 下面我们考虑氢原子的粗略结构 ( ) 所谓的粗略结构是指如下的近似 :( ) 不记及原子核或电子的自旋 ;( ) 电子在静电场中做非相对论性运动通常满足这两个近似条件的原子或离子模型都被称作粗略结构上述两个近似使得氢原子的整体哈密顿算符具有简单的形式 e Ĥ = p N + p e m N m e 4πϵ 0 x e x N, 其中 x e 和 x N 分别是电子和原子核空间坐标矢量算符, 而 p e 和 p N 则是电子和原子核的动量矢量算符在坐标表象中定态薛定谔方程是 Ĥψ( x N, x e )= [ h n h e ] e ψ( x m N m N, x e )=Eψ( x N, x e ). e 4πϵ 0 x e x N 引入个新变量 X m e x e + m N x N m e + m N, r x e x N, 其中 X 是质心系的位置矢量, 而 r 则是原子核指向电子的相对坐标位置矢量因为 = X x e x e X + r x e r = ( ) ( m e = = e x e m e + m N m e m e + m N X + r ) X + r + m e m e + m N X r.

3 氢原子同理得故而有和其中 ( ) ( N = m = N x n m e + m N 1 m e e = ) X + r m N m e + m N X r. m e (m e + m N ) X + 1 m r + e m e + m N X r, 1 m N N = m N (m e + m N ) X + 1 m N r m e + m N X r, 1 e + 1 N m e m 1 N m e + m N X + 1 µ r, µ m em N. m e + m N 此时两体问题的定态薛定谔方程化为 h Eψ = (m e + m N ) ψ h X µ r ψ Ze 4πϵ 0 r ψ } {{ }} {{ } 氢原子自由运动 Ĥ X 电子和原子核的相对运动 Ĥ r. 波函数可以分解为氢原子整体运动部分 ( X) 和电子与原子核的相对运动部分 ( r ) ψ( x e, x N )=K( X)ψ r ( r ), 代入到两体问题的定态薛定谔方程中可得如下两个薛定谔方程 h (m e + m N ) X K( x )=E kk( X), h µ r ψ r Ze 4πϵ 0 r ψ r = E r ψ r, 其中第个方程就是我们前面讨论过的中心势场问题两体问题的总能量为质心系统的动能 (E k ) 和内能 (E r ) 之和 : E total = E k + E r. 注意 : 内能依赖于两个粒子的约化质量 µ 以氢原子为例, µ = m ( em N m = m N m e + m e m N m N + m e e 另外一个例子是氢核和氘核 (): 1 m e m N + ). 氢核 : m p = 1836m e = µ pe = m e 氘核 : m d = 3670m e = µ pe = m e. 如此小的约化质量差异已经可以通过光谱线实验来观测在天文学是哪个, 人们通过观测氢原子和氘原子光谱线之间的相对强度来判断星际介质中氢原子和氘原子的相对残留丰度此信息有助于人们了解在宇宙演化早期时形成氘原子的条件

4 第章中心势场氢原子将氢原子波函数记作为则氢原子中电子的定态薛定谔方程为 ψ(r,θ,φ)= u(r) Y m r l (θ,φ), d l(l + 1) u(r) dr r u(r)+ µe Z 1 h u(r)+µe 4πϵ h u (r)=0. 引入无量纲参数 ρ = 8µE µe r, λ = h 4πϵ 0 h 定态薛定谔方程为 d dρ u l(ρ) h 8µE = 1 a 0 h µe, l(l + 1) ρ u l (ρ)+ λ ρ u l(ρ) 1 4 u l(ρ)=0. 我们注意到, 由于普朗克常数的存在, 我们可以构造具有长度量纲的常数 a 0 来标记氢原子大小, 而且能量 E 和 λ 参数之间存在一一对应的关系如果 λ 只能取特定的整数值, 那么能量就会量子化下面我们可以看到, 波函数在原点和无穷远处的自然边界条件要求 λ 只能取整数考虑在 r 0 和 r 处的波函数性质 ρ : 定态薛定谔方程为其解为 ρ 0: 定态薛定谔方程为 d d dρ u l(ρ) 1 4 u l(ρ)=0, u l (ρ) e ρ/. dρ u l(l + 1) l(ρ) ρ u l (ρ)=0, 其解为 u l (ρ) ρ l+1. 我们猜测氢原子中电子波函数为如下形式 : u l (ρ) ρ l+1 e ρ/ v l (ρ), 其中两个特殊解的乘积 ρ l e ρ/ 满足 r 0 和 r 两个边界处的薛定谔方程, 而 v l (ρ) 则描述 0 <r< 之间的波函数分布波函数收敛有界要求

5 氢原子当 ρ 0 时,v l (ρ) 趋于常数, 否则 l = 0 时波函数行为不好 ; 当 ρ 时,ρ l+1 v l (ρ) 的整体发散性要慢于 e ρ/ 将 u l (ρ) 代入到薛定谔方程中可得 ρv l + [(l + 1) ρ]v l (l + 1 λ)v l = 0 此为合流超几何函数 ( ) 级数解法下面我们仿效一维谐振子势的级数解法来求解上述方程在波函数中分离出 ρ l+1 和 e ρ/ 两种波函数的渐进行为将简化级数求解将 v(ρ) 展开为 ρ 的级数 v(ρ)= v(ρ) 的一阶和二阶导数为 a j ρ j. j=0 dv dρ = ja j ρ j 1 = j=0 (j + 1)a j+1 ρ j, j=0 d v dρ = j(j + 1)a j+1 ρ j 1, 其中我们对哑指标做了替换 j j + 1 故而, 薛定谔方程给出 j=0 j=0 j=0 j(j + 1)a j+1 ρ j + (l + 1) (j + 1)a j+1 ρ j ja j ρ j (l + 1 λ) a j ρ j = 0. 令相同幂次 ρ j 的系数相等, 则有 j=0 j=0 j(j + 1)a j+1 + (l + 1)(j + 1)a j+1 ja j (l + 1 λ)a j = 0, = [j(j + 1)+(l + 1)(j + 1)]a j+1 (j + l + 1 λ)a j = 0, j + l 1 λ = a j+1 = (j + 1)(j + l + ) a j. 注意 : ) 分离出 ρ l+1 项可以避免级数展开中出现多个零系数 ; ) 分离出 e ρ/ 项可以避免出现 a j+,a j+1,a j 的递推关系式, 否则计算非常困难 u(ρ) 在处收敛要求 lim ρ ρl+1 v(ρ) <e ρ/. 下面我们看一下 v(ρ) 级数在 ρ 处的渐进行为 a j+1 lim j a j j + l + 1 λ = lim j (j + 1)(j + l + ) 1 j.

6 因为其级数展开式 a j+1 和 a j 满足 e ρ = j=0 第章中心势场 ρ j j!, a j+1 lim j a = j 1 j! 1 (j 1)! = 1 j. 所以即 v(ρ) j e ρ, u(ρ) ρ l+1 e ρ/ e ρ ρ l+1 e ρ/ ρ 发散, 明显 v(ρ) 的无穷阶的级数展开是不满足平方可积条件的, 所以我们必须对 j 的求和截断, 从而导致能量量子化截断要求 j + l + 1 λ = 0 = λ n = j + l + 1. 因为轨道角动量量子数 l 是整数, 径向波函数级数展开的幂次 j 也是整数, 所以 λ n 必定也是整数, 其取值为 n = 1,,3,, 并且 λ = n = µe h 4πϵ 0 h 8µE = 1 h a B µe = E n = h µa. Bn 氢原子离散能级只依赖于主量子数 n, 与 n r = j 和 l 无关, 能级简并度为 n 因为角动量算符平方所代表的三维空间转动不变性, 每一个轨道角动量 l 都具有 l + 1 个简并, 所以氢原子第 n 个能级的简并度为 n 1 n 1 n 1 g n = (l + 1)= l + 1 = 1 (n 1)n + n = n. l=0 引入精细结构常数 α l=0 l=0 α = e 4πϵ 0 hc 可将氢原子能级公式简化为 λ = n = α m ec ( ) e 1 1 = E E n = 4πϵ 0 a 0 n = 1 m e (αc) 1 n. 精细结构常数刻画电磁相互作用强度的物理学参数因为 e / h 具有速度量纲, 所以在非相对论情形下, 氢原子中电子的速度 (v e ) 和玻尔半径 (a B ) 分别是 v e αc = c 137,

7 氢原子 a B = 4πϵ 0 h m e c = 1 α h m e c = h m e (αc) 氢原子能级公式中包含一个无量纲的常数 α α = e 4πϵ 1 0 hc (3), 我们无法通过改变单位来改变其数值!! 这个常数的出现令物理学家感到非常困惑任何一个物理理论都没有办法预言无量纲量的数值, 因为它是理论中的自由参数人们曾经试图通过各种常数 ( 例如 e,π, 欧拉常数 γ 等 ) 的组合来凑成这个神奇的 1/137, 例如 e π / 1 139, 1 π e / 137.3, 但明显这种做法毫无意义艾丁顿 ( ) 曾经开玩笑提出一个理论 : α 常数出现在氢原子中, 而氢原子是由质子和电子组成的在维时空中, 质子和电子都有个自由度, 这样总共有 4 4 = 16 个自由度的实对称矩阵总共有 16 (16 + 1)/ = 136 个独立矩阵元, 加上电子自旋后, 就给出 = 137 著名物理学家将艾丁顿的数论推广到理想气体从而得到绝对零度最简单的理想气体是 H, 因为每一个氢分子都有个自由度, 所以 H 具有 137 = 74 个自由度形成束缚态 H + H H 给出一个束缚条件, 所以 H 总计具有 74 1 = 73 个自由度开玩笑说, 我们需要降低到度来将理想气体的所有自由度冻结年, 提出重整化群后, 人们才意识到 α 并非是一个常数, 其数值大小依赖于实验能量标度, 这被称作为跑动效应 ( ), 例如 α(q = 100 )=1/18 氢原子波函数当主量子数 n = n r + l + 1 给定时, 径向量子数 n r ( 径向波函数节点个数 ) 和轨道角动量量子数 l 分别取 n r = 0,1,,,n r 1, l = n 1,n,n 3,,0. 径向波函数是 u nl (ρ) F( n r,l +,ρ n )ρ l+1 n e ρ n/, 其中 ρ n 是第 n 能级所对应的特征长度, 用来刻画处于第 n 能级电子相对于氢原子核的距离, ρ n = 8µEn h r = na B r.

8 第章中心势场氢原子归一化的波函数为 = ψ nlm = R nl (r)yl m (θ, φ) ( ) 3/ (n + l)! 1 na B n(n l 1)! (l + 1)! ρl ne ρn/ F( n r,l +,ρ n )Yl m (θ,φ). 基态和第一激发态的径向波函数为而完整波函数为 ( ) 3/ 1 R 10 = e r/a B, R 0 = R 1 = u 100 = u 00 = u 10 = u 11 = u 1 1 = a B ( ) 3/ ( 1 r ) a B a B ( ) 3/ 1 r e r a B a B 3 1 πa 3 B e a r B e r a B, a B, ( 1 r ) e a r B 4 πa 3 a B B 1 r e r a B cosθ 4 πa 3 a B B 1 r e r a B e iφ sinθ a 8 B 1 8 πa 3 B πa 3 B r e r a B e iφ sinθ a B 下面我们讨论氢原子的几个较低能级的径向波函数分布因为 u nl (r) =rr nl (r) 归一化简单, 4π 0 u(r) dr = 1, 我们选取 u(r) 来解电子相对于氢原子核的位置分布首先, 我们先研究 n r 的物理意义 n r = n l 1 描述径向波函数的节点数, 它标志电子的径向方向运动的激烈程度 n r 大就意味着电子在径向方向上运动活跃, n r = 0 则意味着电子在径向方向上分布变化缓慢, 但此时轨道角动量 l = n 1 达到极大值, 电子的空间方位分布变化剧烈例如, n = 1 时,n r = l = 0, 径向波函数无节点 ; n = 时,

9 氢原子 n r = 1, l = 0, 径向波函数有一个节点 ; n r = 0, l = 1, 径向波函数无节点 ; n = 3 时, n r =, l = 0, 径向波函数有两个节点 ; n r = 1, l = 1, 径向波函数有一个节点 ; n r = 0, l =, 径向波函数无节点图氢原子 n = 3 能级的径向波函数分布 ( 横坐标单位是玻尔半径 ):l =, n r = 0 ( 绿色 ),l = 1, n r = 1( 蓝色 ),l = 0, n r = ( 红色 ) 其次, 我们研究电子波函数的径向分布最概然位置因为所以电子的径向几率分布为 n r = 0 : R n,n 1 (r) r n 1 e r na B, P n,n 1 Rn,n 1 r r n e r 几率分布的最大值对应的径向位置最概然位置由下式决定 : dp n,n 1 dr na B. = 0 = r max = n a B. 这和圆周运动的经典图像相符图形 ( ) 显示氢原子 n = 10 能级的径向几率分布随电子径向半径的变化, 此几率分布在 r = 100a B 处达到极大值图氢原子 n = 10 能级的径向几率分布 ( 横坐标单位是玻尔半径 ):l = 9, n r = 0

10 第章中心势场图氢原子 n = 10 能级的径向几率分布 ( 横坐标单位是玻尔半径 ):l = 0, n r = 9 当 n 数值较大时, 根据对应性原理, 氢原子中电子的几率分布要接近经典力学情形图形 ( ) 显示 n = 10,l = 0,n r = 9 波函数的径向几率分布注意 : 最概然几率的位置 r max 比 100a B 小这可以从经典物理理解在库伦势或引力势中运动粒子的轨道形成一个封闭椭圆, 具有相同半主轴的椭圆轨道对应于相同能量因为 l = n 1 轨道是圆形, 而 l = 0 轨道是一个长细椭圆, 能量守恒要求此长细椭圆的半主轴取可能的最大值, 这就要求电子离氢原子核的距离尽可能接近圆形轨道半径的两倍, 即 100a B = 00a B 在经典物理中围绕氢原子核做椭圆运动的电子在离原子核最远处耗时较长, 因为在最远处电子的速度最小运动最慢故而, 图中电子在远离氢原子核处出现几率最大的结果完全符合经典物理图像同时我们也可以看出经典和量子的区别, 虽然 n = 10 时, 量子物理结果已经接近经典物理, 但还是和经典物理有所不同在经典物理中 r 00a B, 量子物理中电子还有一定几率出现在经典禁区 (r >00a B ) 中和电子径向位置有关物理量的平均值我们经常会遇到求解和电子径向位置有关的物理量的平均值给出一个非常有用的递推关系, 将电子径向位置的不同幂次的平均值联系起来 : k + 1 n nl r k nl (k+1)ab nl r k 1 ka [ nl + B (l + 1) k ] nl r k nl = 0. 4 首先 r 1 位力定理告诉我们, 库仑势场中 E ˆT E = E ˆV E, 所以即 E ˆV Ze E = E 1 E = E = Z e 4πϵ r 4πϵa 0 n, r 1 = Z n a 0.

11 氢原子再计算 r 3 考虑对易子 [Ĥ, ˆp r ] 的平均值因为定态中任何含有哈密顿算符的对易子的平均值为零, 0 = E [Ĥ, ˆp r ] l(l + 1) h E = E [r, ˆp µ r ] Ze E E [r 1, ˆp 4πϵ r ] E. 在利用 [r, ˆp r ]=i h 可得 E r 3 E = Z 3 a 3 0 n3 l(l + 1)(l + 1 ). 通过氢原子例子, 我们看到求解中心势场中的定态薛定谔方程就可以得到氢原子的全部信息能量量子化是要求物理解满足如下两个自然边界条件的必然结果 : r 0 时,R(r) r l ; r 时,R(r) 收敛氢原子能级的代数解法氢原子具有更高简并度的原因是氢原子具有更高的动力学对称性 SO(4)( 详细内容请参见程檀生老师教材第页附注 ) 下面我们不讲解具体的对称性, 而是采用代数解法构造升降算符来求解氢原子的能级, 并以此说明氢原子哈密顿算符中包含超出三维转动不变性的更高对称性设氢原子本征函数为 E,l, 满足本征方程我们写出哈密顿算符其中玻尔半径 a 0 为 Ĥ E,l = E E,l, ˆ L E,l = l(l + 1) h E,l. Ĥ l = ˆp r l(l + 1) h + µ µr Ze 4πϵr = ˆp r l(l + 1) h + µ µr Z ( h ), r µa 0 a 0 = 4πϵ 0 h µe. 定义算符 Â l 和 Â l 如下 : Â l a ( 0 ī h ˆp r l + 1 ) Z +, r (l + 1)a 0 Â l a ( 0 ī h ˆp r l + 1 ) Z +, r (l + 1)a 0

12 第章中心势场则氢原子的哈密顿算符 Ĥ l 可以写作为 ( Ĥ l = h Â µa l Â Z ) l 0 (l + 1). 下面我们推导这个哈密顿算符形式考虑 Â l Â l 的乘积 : A l A l = a 0 ( Z l + 1 )( i (l + 1)a 0 r Z l + 1 (l + 1)a 0 r ] ī h ˆp r + h ˆp r + = a 0 ˆp ( r h + Z l + 1 ) + i [ ˆp (l + 1)a 0 r h r, l + 1 r = a 0 ˆp r h + Z ( ) l + 1 (l + 1) a + Z (l + 1) r ra 0 0 r { = a 0 ˆp r h + (l + 1) (l + 1) r Z } Z + ra 0 (l + 1) { = a 0 ˆp r l(l + 1) + h r Z } Z + ra 0 (l + 1) { } = a 0 h ˆp r l(l + 1) h + r Z h Z + ra 0 (l + 1) { } = µa 0 ˆp r l(l + 1) h h + µ µr Z h Z + µra 0 (l + 1) = µa 0 Z h Ĥl + (l + 1). 从而给出如下的哈密顿算符 ( Ĥ l = h Â µa l Â Z ) l 0 (l + 1). 这类似于一维简谐振子势的哈密顿算符形式下面我们推导两个有用的算符关系式首先考虑 A l 和 A l 的对易子 : [A l,a l ]=a 0 因为 [ 所以 [ i h ˆp r l + 1, ī r h ˆp r l + 1 ] = ī [ r h a 0 ˆp r, l + 1 ]. r ˆp r, l + 1 ] = l + 1 r r [ ˆp r,r]=i h l + 1 r, [A l,a l ]= i a 0 h l + 1 l + 1 r (i h)=a 0 r. 此对易关系也类似于简谐振子的升降算符对易关系 [â,â ]=1 由于 ) Ĥ l+1 Ĥ l = h h [(l + 1)(l + ) l(l + 1)] = (l + 1), µr µr

13 氢原子我们有故而, l + 1 r = µ h (Ĥ l+1 Ĥ l ). [A l,a l ]=µa 0 h (Ĥ l+1 Ĥ l ). 其次考虑 Â l 和 Ĥ l 的对易子由公式 ( ) 和 ( ) 可得 [Â l,ĥ l ]= h a [Â l,â l Â l ]= h 0 µa [Â l,â l ]Â l =(Ĥ l+1 Ĥ l )Â l, 0 即有 Â l Ĥ l = Ĥ l+1 Â l. 下面我们利用上面的公式推导氢原子的能量本征值记 Ĥ l 算符的本征态为 E,l 满足如下的能量本征方程 Ĥ l E,l = E E,l. 将算符 Â l 作用在此方程两侧, Â l Ĥ l E,l = Â l E E,l = EÂ l E,l, 再利用公式 ( ) 可得所以 Â l Ĥ l E,l = Ĥ l+1 Â l E,l, Ĥ l+1 (Âl E,l ) = E ( Â l E,l ). 这说明 Â l E,l 仍然是氢原子的能量本征态, 而且它对应于相同的能量本征值 E 同时 Â l E,l 也是 Ĥ l+1 算符的本征态, 说明 Â l E,l 的轨道角动量量子数为 l + 1, 即 Â l E,l E,l + 1 或 E,l Â l E,l + 1. 位力定理告诉我们 : 氢原子能量保持不变时, 其动能也是不变的为保持动能不变, 提高轨道角动量量子数的同时,Â l 算符也必然要减少径向量子数 n r 如图 ( ) 所示以此类推, 我们将算符 Â l+1 作用在 E,l + 1 上可以得到 E,l +, 例如 E,l Â l E,l + 1 Â l+1 E,l + E,l max 1 Â l max 1 E,l max. 上面每一步操作都会增加一个轨道角动量量子数并且减少一个径向量子数, 最终会得到一个为零的径向量子数 (n r = 0) 和最大的角动量量子数 l max 此时氢原子中电子将做圆周运动, 其状态满足 Â lmax E,l max = 0.

14 第章中心势场 n r Â l Â l l 图 Â l 和 Â l 作用下氢原子角动量量子数和径向量子数变化我们可由此条件得到氢原子能量本征值的表达式因为所以 0 = Âlmax E,l max = E,lmax Â l max Â lmax E,lmax Z = E,l max a 0 h µĥ lmax + (l max + 1) E,l max = a 0 h µe + Z (l max + 1), Z E = h µa 0 (l max + 1) = Z h µa = Z e 0 n 8πϵ 0 a 0 n, 其中 n = l max + 1 是主量子数因为 Â l 作用是将 l 增加为 l + 1, 而且增大步长是 1, 所以主量子数的表达式是 n = n r + l + 1 我们看到, 氢原子能级简并度高的根源隐藏在氢原子哈密顿算符之中三维简谐振子三维各向同性简谐振子势的哈密顿量为 Ĥ = ˆp r m + ˆL mr + 1 mω r. 选取 {Ĥ, ˆL, ˆL z } 作为力学量完全集, 波函数记作为 ψ nlm 定态薛定谔方程是 ( h 1 m r r r ˆL ) h ψ r nlm + 1 mω r ψ nlm = Eψ nlm. 将波函数分解为径向和空间角度两部分 ψ nlm = R nl (r)y lm (θ,φ)= u nl(r) Y r lm (θ,φ),

15 带电粒子在磁场中运动代入到定态薛定谔方程中, 并取无量纲参数则可得径向薛定谔方程 ρ = r α, ρ u(ρ)+ λ = E hω, 其中 α = [ λ ρ h mω, ] l(l + 1) ρ u(ρ)=0. ( ) 当 ρ 时, 径向方程化为如下的近似形式其在无穷远处收敛解为 ρ u(ρ) ρ u(ρ)=0, u(ρ) e ρ. ( ) 当 ρ 0 时, 波函数的行为和氢原子在 ρ 0 时相同, 所以我们猜测波函数解具有以下形式 u(ρ) ρ l+1. u(ρ)=ρ l+1 e ρ v(ρ), 同时要求 ρ l v(ρ) 在无穷远处发散行为不超过 exp( ρ /) 将 u(ρ) 代入到径向薛定谔方程中可得 ρv (ρ)+ [ (l + 1) ρ ] v (ρ)+(λ l 3)ρv(ρ)=0. 令 y = ρ, 可将上式化作为合流超几何微分方程 [( yv (y)+ l + 3 ) ] y v (y) l + 3 λ v(y)=0. 4 在 y = 0 处有正常解要求 ( l + 3 λ v(y)=cf,l + 3 ). 4,y 为使在无穷远处 R nl (r) 0, 要求截断多项式, 即有这给出量子化能级 l + 3 λ 4 = n r λ = E hω = 4n r + l + 3. ( E = hω n r + l + 3 ) ( N + 3 ) hω, 其中 N = n r + l, 而常数 3/ 则是三维空间中谐振子势的基态能具体归一化波函数形式参见程老师书中公式 ( )

16 第章中心势场带电粒子在磁场中运动处理处于外磁场 ( B) 中的带电粒子运动问题, 我们可以将自由粒子动量替换成机械动量 ( ) p p q A, 其中右边的 p 为正则动量,q 是粒子电荷, A 是满足如下关系的电磁场矢量势 A = B, A = 1 B r. 考虑处于外电磁场中的氢原子, 氢原子中电子的哈密顿算符为 Ĥ = 1 ( p q A ) e µ r = p µ e r q ( p A + A p ) + q A µ m. 选取库伦规范 A = 0, 则有 p Aψ = i h A ψ i h A ψ = A p ψ, }{{} 即有 0 p A = A p. 代入 A = 1 B r 后得 A p = 1 ( B r ) p = 1 B ( r p ) = 1 B L. 所以当仅有外部磁场时, 系统的哈密顿算符是其中利用矢量关系式 Ĥ = Ĥ 0 q A µ p + q A µ = Ĥ 0 µ L B + q A µ, µ L = q µ L = µ B h L, µ B = q h µ ( 玻尔磁子 ) ( C D ) ( E F ) = ( C E )( D F ) ( C F )( D E ), 可得 A = 1 4 ( B r )( B r ) = 1 4 [ B r ( B r ) ]. 将之代入到哈密顿算符中得 Ĥ = p µ e r q µ B L + q 8µ [ B r ( B r ) ].

17 带电粒子在磁场中运动效应当氢原子类氢原子或碱金属原子处于较强磁场中其能级发生移动的现象称为效应 : 不考虑电子自旋效应时被称作为正常效应, 而考虑电子自旋时称作为反常效应下面我们计算正常效应因为原子空间尺寸非常小, 实验室中制备的磁场在原子尺度上可视作为均匀磁场取 B 沿着 z 轴 (B = Bê z ), 哈密顿量为 Ĥ = ˆp µ e r + e µ BˆL Z + e B 8µ (x + y ). 当磁场不强时,B 项可以忽略, 从而哈密顿算符是 Ĥ = Ĥ 0 + Bµ B h ˆL Z. 因为 [Ĥ,Ĥ 0 ]=[Ĥ, ˆL Z ]=[Ĥ 0, ˆL Z ]=0, 所以我们选取力学量完全集 {Ĥ, ˆ L, ˆL Z } 因为力学量完全集和氢原子一样, 所以波函数通解形式为相对应的本征值为 ψ n,l,m (r,θ,φ)=r nl (r)y m l (θ,φ) nlm, E nlm = nlm Ĥ nlm = nlm Ĥ0 nlm µ + B B nlm ˆL nlm h Z = E 0 + mµ B B = E 0 + m hω L, 其中 ω L = eb/µ 是拉莫 ( ) 频率与氢原子不同的是, 此时系统的能级依赖于磁量子数 m 这是因为外加磁场破坏了氢原子的三维旋转不变性原来自由氢原子的 ˆL X,Y,Z 都是等价的, 但沿 z 轴的外加磁场选取了特定的空间方位, 磁矩和磁场的相互作用解除了原来氢原子能级关于空间方位的 l + 1 维简并但此相互作用仅仅移除氢原子关于 ˆ L 子空间的简并, 但 n = n r + l + 1 简并尚未完全解除如图 ( ) 所示, 氢原子最低三个能级在外磁场下关于 m 量子数的简并解除了, 但氢原子关于 l 量子数的简并仍然存在, E nl m = E nlm (l l ). 均匀磁场中的自由电子朗道能级年, 年仅岁的朗道研究均匀磁场中的自由电子运动情况并发展了金属电子的朗道抗磁理论在文章中朗道提出了后来的量子霍尔效应理论中使用的基本概念电子在磁场中的朗道能级和能级的态密度等

18 第章中心势场 (3, 0), (3, 1), (3, 3) (3,, 1) (3, 1, 1), (3, 1) (3, 0, 0), (3, 1, 0), (3,, 0) (3,, 1), (3, 1, 1) (3,, ) E (0) 3 +µ B B E (0) 3 + µ B B E (0) 3 E (0) 3 µ B B E (0) 3 µ B B (, 0), (, 1) (1, 0) (1, 0) B =0 (, 1, 1) (, 0, 0), (, 1, 0) (, 1, 1) B 0 E (0) + µ B B E (0) E (0) µ B B E (0) 1 图正常效应的氢原子能级 : 左边括号内是量子数 (n,l), 右边括号内是量子数 (n,l,m) 朗道生不逢时, 如果他早生几年, 一定可以在量子理论创建过程中大展拳脚朗道还用数学符号对理论物理学家分类 : 拉普拉斯算符是头脑尖又坐得住者, 如爱因斯坦 ; 达朗贝尔算子代表头脑笨而坐得住者, 如他的某些同事 ; 他本人由描述, 头脑虽尖但坐不住, 最后表示头脑笨又坐不住的人, 例如他讨厌的某人描述磁场中的电子运动的薛定谔方程是 1 ( p q A ) ψ = Eψ, m 此处 B 项不可忽略在前面氢原子问题中, 因为电子所受的库仑势大于电子和磁场之间的相互作用, 因为当磁场不强时仅需考虑 B 的线性项而忽略 B 项考虑外磁场沿着 z 轴方向, 我们可以选择不同的矢势 A, A = 1 B( y,x,0) 或 A = B( y,0,0). 下面我们可以看到这两种不同的选取方案会给出相同的物理结果 A = 1 ( y,x,0) 哈密顿算符是 Ĥ = 1 m = 1 m [ ( ˆp x 1 eby) +(ˆp y + 1 ] ebx) + ˆp z ( ˆp x + ˆp y) + e B 8m (x + y )+ eb m (x ˆp y y ˆp x )+ ˆp z m.

19 带电粒子在磁场中运动明显, 粒子沿着 z 轴方向做自由运动, 所以我们可以将 ˆp z 从哈密顿算符中分离出, 从而仅需考虑电子在 x y 平面上的运动相应的哈密顿算符是 Ĥ xy = 1 ( ˆp m x + ˆp y) 1 + mω L (x + y ) +ω L ˆL z, ω L = e B m, } {{ } Ĥ 0 这里 ω L 项 ( 也即 B 的线性项 ) 描述带电粒子磁矩通外磁场之间的相互作用, 而 ω L 是 B 项, 在后边讨论中我们会看到这一项代表着反磁项注意到 Ĥ 0 形式和二维各项同性的谐振子势相同, 我们选取力学量完全集为 {Ĥ, ˆL z }, 并在极坐标系 (ρ,φ) 中写出波函数通解为将波函数带入到定态薛定谔方程中得 ( [ h m ρ + 1 ρ ρ m ρ ψ(ρ,φ)=r(ρ)e imφ, m = 0,±1,±,. ) + 1 mωl ρ ] R(ρ) = (E m hω L )R(ρ) 这正是二维谐振子势的能量本征方程选取如下变量 h ρ = αy, α =, ϵ = E, mω L hω L E R(ρ) 我们将二维谐振子势的能量本征方程转化为如下的无量纲形式 ( 注意 : 无量纲化后长度的度量单位是 α, 能量的单位是 hω L ) d R(y) dy + 1 ( ) dr(y) + ϵ y m y dy y R(y)=0. 此本征方程的求解是完全类似于我们求解氢原子本征值和本征函数, 我们首先考虑本征方程在两个自然边界条件处的行为 : 当 y 时, 上面方程简化为 ( d dy y ) R(y)=0 = R(y) e y / ; 在 y 0 时, 我们得到如下近似方程 ( d dy + 1 y ) d dy m y R(y)=0. 令 R(y) y s 代入上面方程得到 s m = 0, 故而 s = ± m 所以 R(y) y 0 y + m.

20 第章中心势场所以我们猜测径向波函数的通解形式为 R(y)=y + m e y / v(y), 其中要求 v(y) 在 y = 0 附近趋于零, 而在 y 时 y + m v(y) 的总体发散性不超过 e y / 将 R(y) 代入到径向方程中得 d ( ) v(y) m + 1 dv(y) dy + y + [ϵ ( m + 1)]v(y)=0. y dy 做变量替换 η = y, 则上面方程可以转化为熟悉的合流超几何方程 [ η d v dv m ( m + 1 η) dη dη ϵ ] v = 0. 4 在 η = 0 处有正常解要求 ( m + 1 v(η) F ϵ ). 4, m + 1,η 同时要求 η 时 R(ρ) 收敛, 我们必须截断多项式, 即要求这导致二维谐振子势的能量值量子化, m + 1 ϵ 4 = n ρ, n ρ = 0,1,, ϵ = (n ρ + m +1) = E =(N + 1) hω L, } {{ } 其中 N = n ρ + m = 0,1,, 相应的波函数为即其中 ψ nρm (ρ,φ) e imφ ρ m e ρ α F N ) ( n ρ, m + 1, ρ h h, α = = mω L eb. 从上面二维谐振子势出发, 我们可以得到在磁场中运动电子的能级为 α E = E m hω L =(n ρ + m + 1) hω L, E =(n ρ + m + m +1) hω =(N + 1) hω L. } {{ } N N = 0,,4,, n ρ = 0,1,, 当 m>0 时,N = n ρ +m, 其简并度貌似有限 N /+1 但当 m 0 时,N = n ρ 与 m 值无关, 此时简并度是无穷大所以我们发现均匀磁场中运动电子, 电子能量始终大于零且能级是离散的, 每一个能级的简并度都是无穷大

21 带电粒子在磁场中运动图外磁场中带电粒子回旋运动的经典物理图像电子能量可以看作是电子在外磁场中感应而生的磁矩和外磁场相互作用, 此相互作用可写作是 Ĥ = ˆµ z B, 因为电子能量大于零 ( 量子世界中始终存在非零的零点能 ), 所以 Ĥ>0 就要求电子磁矩 µ z = (N + 1) e h m, 始终与外磁场方向相反, 这就是电子的反磁性这种反磁性和带电粒子的电荷符号无关, 正电荷粒子也具有反磁性当 q>0 时, 正电荷粒子的能级是 E =(n ρ + m m+1) hω =(N + 1) hω L, } {{ } N 此时只有 m<0 磁量子数才贡献到离散能级, 而 m>0 的所有量子态的能量相等, 导致无穷大简并总结一下 : q>0:e =(n ρ + m m + 1) hω L, 仅有 m<0 的量子数进入到能量本征值中 ; q<0:e =(n ρ + m + m + 1) hω L, 仅有 m>0 的量子数进入到能量本征值中这两种量子行为粗看起来非常奇怪, 但如果从经典物理角度出发, 我们发现它们又是非常合理的经典物理中带电粒子在外部磁场中的运动行为是确定的 ; 当给定外磁场方向和粒子电荷时, 带电粒子所受的洛伦兹力方向也是确定的如图 ( ) 所示, 当磁场沿着 z 轴时, 正电荷粒子在 x y 平面做顺时针回旋 ( 沿 z 轴角动量为负, L z < 0), 而负电荷粒子则做逆时针回旋 ( 沿 z 轴角动量为正,L z > 0) 综上所述, 仅为 (q >0,m<0) 和 (q <0,m>0) 两种特定的量子行为才具用经典对应注意 : 在量子力学中我们求解的量子体系的空间波函数的分布, 是定态问题, 我们无法得到经典物理中的运动形式上面的讨论只是类比经典物理情况来理解下面我们看一下上面推导的外磁场中带电粒子的量子状态的经典对应因为经典物理中, 带电粒子在垂直于磁场的 x y 平面做回旋运动, 没有沿径向方向的运动, 所以我们量子体系的经典对应是 n ρ = 0 和 m 此时带电粒子在 x y 平面中出

22 第章中心势场现的几率密度是 P(ρ)= ψ0,m (ρ,φ) ρ m e ρ α, 其中 α = h 和 ω µω L = q B 0 L µ. 带电粒子在空间中出现的最概然几率发生在 dp(ρ) = dρ m ρ m 1 e ρ ( α + ρ m ρ ) α e ρ α = 0 从而得到 ρ max = m α = m h µω L m L z µω L, 其中最后一步中当 m 时体系自由度无限多, m h 对应于经典物理中轨道角动量, 物理体系的特征作用量 L z h 在经典物理中, 带电粒子在外磁场中运动状态由以下方程决定 : 从第二个公式我们可得 µv = F = qvb r 0 µ v = p q A. µ( r v ) z =( r p ) z q( r A) z. 因为 A ( B C)=( A C) B ( A B) C, 所以有我们得到由公式 ( ) 和 ( ) 可得 r 1 ( B r )= 1 ( r r ) B 1 ( 0 r B) r = 1 Br, µrv = L z 1 qbr. r = L z 经典 qb = L z. µω L 显然, 经典物理的带电粒子运动的回旋半径和量子体系的经典近似一致 A = B ( y,x,0) 代数解法上面微分方程解法复杂, 我们现在采用代数解法重新求解外磁场中带电粒子的运动问题这种解法不仅简便, 而且可以更清楚地揭示无穷大简并的物理根源所在体系在 x y 平面的哈密顿量算符为 Ĥ xy = 1 m {( ˆp x + 1 qby) +(ˆp y 1 qbx) }.

23 带电粒子在磁场中运动引入两个厄米算符 ˆΠ x = ˆp x + 1 mωy m hω, ˆΠ y = ˆp y 1 mωx m hω, 其中可以将 Ĥ xy 改写作 ω = q B m = ω L. Ĥ xy = 1 hω ( ˆΠ x + ˆΠ y). 容易验证 [ ] 1 ( ˆΠ x, ˆΠ y = [y, ˆpy ] [ ˆp h x,x] ) = i. 注意到 Ĥ xy 算符中 ˆΠ x 和 ˆΠ y 都是平方形式, 和谐振子势的无量纲的哈密顿算符形式完全相同, 所以我们定义升降算符它们满足如下对易关系 â = 1 ( ) ˆΠ x + i ˆΠ y, â = 1 ( ) ˆΠ x i ˆΠ y, [â,â ] = 1. 我们利用升降算符可以将哈密顿算符改写为 ( Ĥ xy = hω â â + 1 ) = hω L (â â + 1)= hω L (N + 1), N = 0,,4 此时, 整个物理问题简化为一维振子, 明显哈密顿算符中缺少了一个自由度这个旁观的自由度导致物理体系的能级具有无穷多的简并度 A = B( y,0,0) 我们通过代数解法将哈密顿算符转化为一维简谐振子问题, 但丢失的自由度是什么? 下面我们选取第二种矢势 A 来讨论此矢势与前面讨论过的第一种矢势仅仅相差一个规范变换 1 ( By,Bx,0)=( By,0,0)+ f (x,y,z), f = 1 Bxy. 为保证薛定谔方程满足规范不变, 新哈密顿算符的本征波函数要添加一个相位因子在新矢势 A 下体系哈密顿算符是 ψ 新 (x,y,z)=e iqf (x,y,z)/ h ψ 旧 (x,y,z). Ĥ = 1 [ ( ˆpx eby) + ˆp m y] + p z m. } {{ } Ĥ xy

24 因为 [Ĥ,pz ] = [Ĥ,px ] = [px,p z ] = 0, 第章中心势场我们可以选取 {Ĥ, ˆp x, ˆp z } 组成的力学量完全集其波函数通解为 ψ n,px,p z = 1 p x x+p z z π h ei h η n (y), 其中 p x 和 p z 分别是动量算符 ˆp x 和 ˆp z 的本征值经典的数将通解代入到二维的薛定谔方程中可得 1 [ ˆp m y +(p x eyb) ] η n (y)=ϵ n η n (y), 此为一维谐振子势本征方程, 只不过平衡点不在原点, 而是在 y 0 = p x /eb 处令可得 ( h m 能量本征值为总能量为 s = y y 0 = y p x eb, ω = e B m, d ds + 1 ) mω s η n (s)=ϵ n η n (s). ( ϵ n = n + 1 ) hω =(n + 1) hω L, ϵ N =(N + 1) hω L, N = n = 0,,4, E n,px,p z = p z m +(N + 1) hω L. 因为 E n,px,p z 与 p x 无关,p x 本征值仅仅表现在谐振子平衡点位置, 所以能量对于 p x 是简并的, 而且简并度是无穷大这说明均匀磁场中带电粒子可以出现在无穷远处 ( 即非束缚态 ), y 0 = p x eb ±, 但电子能级却是离散的通常的二维非束缚态粒子的能量是连续变化的简并根源简并度无穷大的原因在于整个物理体系对于 x y 平面没有任何限制 :x 和 y 的取值都是从到 + 这相当于假设磁场是无穷大对于无穷大的物理系统, 我们无法定义 x y 平面的原点位置整个物理体系在 x y 平面上具有平移不变性, 这导致体系具有不穷大简并当物理体系被局域在有限大小空间时, 我们可以定义物理体系的原点位置, 此时朗道能级的简并度就是有限的例如考虑一个物理体系被局限在一个边长为 L 的正方形箱子中, 并非所有的圆轨道都满足边界条件, 此时物理体系的能级简并度为 N d = qbl π h.

25 带电粒子在磁场中运动简并度正比于面积的行为可以从相空间理解关于各种受限物理体系的朗道能级的详细推导请参见所著的一书物理应用我们可以估算一下 ω L 176 ( ) ( ) B B GHz, hω 1 L = ev. 1 室温下电子热运动动能 T 0.03, 所以一般室温下电子的朗道能级不显著我们需要在低温和强磁场情况下才可以观测到朗道能级的效应中子星的磁场非常强, 约为 10 8, 是非常适合研究朗道能级的物理系统 : 当电子从一个朗道能级跃迁到另外一个朗道能级时会发出硬光子射线

幻灯片 1

幻灯片 1 第二章原子的激发态结构. 单电子 (H) 原子 ( 类 ) 氢原子的薛定谔方程 -e 3D 不含时的定态薛定谔方程其中库仑势 m + ( ) V ( ) V ψ Ze 4πε Eψ + +Ze 电子束缚在原子核的中心力场中只与电子和原子核之间的径向距离有关 Fom www.hpephsics.ph-ast.gsu.edu . 单电子 (H) 原子氢原子的薛定谔方程 Catesian coodinates