磁铁吸引金属时发生了什么?

当你目睹磁铁吸起一块铁时,你实际上正在见证经典物理学中最深刻的“漏洞(Glitch)”之一。为了理解那块金属为什么会动,我们必须拆毁经典直觉,并利用狄拉克(Dirac)、泡利(Pauli)和海森堡(Heisenberg)的工具来重建物理大厦。

1. 悖论:为什么洛伦兹力是“懒惰”的

在经典电磁学中,一个电荷 q 以速度 v 在磁场 B 中运动时,所受的力由洛伦兹力公式给出:

Fmag=q(v×B)

根据定义,这个力永远垂直于速度方向(Fv=0)。由于功率是力做功的速率(P=Fv),因此磁场力不做功

然而,当磁铁吸引金属时,金属明明获得了动能。如果磁场力不能做功,到底是什么在拉动金属?答案在于:磁性并非一种纯粹的经典力,而是一种相对论性的量子修正

2. 死胡同:玻尔-范雷文定理 (Bohr-van Leeuwen Theorem)

在寻找答案之前,我们必须先证明经典物理学根本无法产生磁性。这就是著名的玻尔-范雷文定理

考虑一个由 N 个电子组成的系统。磁矩 μ 由电荷的轨道运动定义:

μ=q2r×v=q2m(r×mv)=q2mL

关键在于,z 方向的总磁化强度是广义速度 q˙i 的线性函数:

Mz=i=13Nai(q1,,q3N)q˙i

在电磁场中,系统的哈密顿量(总能量)H 为:

H=i=1N(piqA)22m+qV(q1,,qN)

根据哈密顿方程,速度与动量的关系为 q˙i=Hpi。磁化强度的热力学平均值 Mz 由配分函数给出:

Mz=dq1dq3Ndp1dp3NMzeH/kBTdq1dq3Ndp1dp3NeH/kBT

让我们聚焦于单项 aiq˙i 的动量积分。由于 ai 仅依赖于坐标,我们可以先对动量 pi 进行积分:

+dpiq˙ieH/kBT=+dpiHpieH/kBT

这个积分本质上是对一个全微分的积分。注意到 HpieβH=1βpi(eβH)(其中 β=1/kBT),我们可以将其写为:

H(p=)H(p=+)dHeH/kBT=[kBTeH/kBT]pi=pi=+

由于动能项 (pqA)22m 当动量 p± 时趋于无穷大,玻尔兹曼因子 eH/kBT 在边界处趋于零。因此:

[kBTeH/kBT]+=00=0

结论: M=0。物理图像是无情的:虽然磁矢势 A 平移了动量,但由于我们在整个动量空间(从 +)上进行积分,这种平移对总积分没有任何贡献。经典物理预言,磁铁是不可能存在的。

3. 超越薛定谔方程

薛定谔方程建立在非相对论能量-动量关系 E=p2/2m 之上。如果我们试图通过替换量子算符 Eitpi 来直接使用相对论关系 E2=c2p2+m2c4,我们将得到克莱因-戈登方程 (Klein-Gordon Equation)

1c22ψt22ψ+(mc)2ψ=0

这个方程可以用达朗贝尔算符 (d'Alembertian) =1c22t22 更紧凑地写出,这是拉普拉斯算符在闵可夫斯基时空中的推广:

(+(mc)2)ψ=0

然而,当克莱因-戈登方程最初被提出时,它几乎被抛弃了,因为如果像薛定谔方程那样将 ψ 解释为单粒子波函数,会出现两个严重的概念问题。

一方面,相对论能量关系是能量的二次方程,导致两个解:E=±(pc)2+(mc2)2负能量态的存在是个大麻烦。它暗示粒子可以通过跃迁到更低的负能级来持续辐射能量,这意味着物质不存在稳定的基态。

另一方面,量子力学中的波函数用于构建守恒的概率流 jμ=(ρ,j)。对于薛定谔方程,概率密度 ρ=|ψ|2 总是非负的。而对于克莱因-戈登方程,导出的守恒密度为:

ρ=i2mc2(ψψtψψt)

由于方程对时间是二阶的,ψψt 的初始值可以独立选择。这允许 ρ 在某些条件下为负值,这对于概率密度来说是荒谬的。粒子不可能以“负概率”出现在某处。

可以看出,这些问题源于时间导数的二阶性质。狄拉克(Dirac)的目标是找到一个新方程,它既像薛定谔方程那样对时间是一阶的,又是相对论协变的。自然的方法是构造一个动量的一阶形式,系数待定:

H^=c(αxp^x+αyp^y+αzp^z)+βmc2

其中 p^k=ixk 是动量算符。对应的波动方程为:

iψt=H^ψ=(ick=13αkxk+βmc2)ψ=(icα+βmc2)ψ

为了使该方程在相对论上自洽,它的任何自由粒子解也必须满足相对论能量关系(即也是克莱因-戈登方程的解)。这要求当我们将狄拉克算符作用两次(“平方”该方程)时,能还原出克莱因-戈登算符:

(cαipi+βmc2)2=c2pi2+m2c4

展开左式(注意 pi 彼此对易,但系数矩阵 αi,β 不一定对易):

c2iαi2pi2+c2i<j{αi,αj}pipj+mc3i{αi,β}pi+β2m2c4

通过对比左右两边的系数,我们得到了对 αiβ 的代数要求,称为克利福德代数 (Clifford Algebra)

  1. αi2=I
  2. β2=I
  3. {αi,αj}=0(ij)
  4. {αi,β}=0

为了让方程在洛伦兹变换下表现得更显式协变,我们定义 Gamma 矩阵 γμ

代入后得到简洁的协变形式:

(iγμμmc)ψ=0,其中 μ=(1ct,)

利用前面 αiβ 的性质,我们可以推导出 γμ 的核心性质:

综合以上结果,Gamma 矩阵满足著名的反对易关系 (Anti-commutation Relation)

{γμ,γν}=2ημνI

其中 ημν=diag(1,1,1,1) 是闵可夫斯基度规。

满足此代数要求的最小不可约表示是 4 阶矩阵,我们可以通过排除法来推导:

在四阶矩阵中,最常用的是狄拉克表示(Dirac Representation),它使用了 2×2 的单位矩阵 I 和泡利矩阵 σi

γ0=(I00I),γi=(0σiσi0),γi=βαi=(I00I)(0σiσi0)=(0σiσi0)σx=(0110),σy=(0i+i0),σz=(1001)

这是一个深刻的结果:为了得到线性的、相对论性的波动方程,波函数 ψ 被迫必须具有四个分量。我们现在知道,这分别描述了电子(自旋上/下)和它的反粒子——正电子(自旋上/下)。狄拉克试图解决负概率问题,却意外地打开了通往自旋的大门。这说明磁性不是物质的附加属性,而是时空对称性的必然结果。

4. 耦合电磁场

现在我们考虑存在电磁势 Aμ=(Φ/c,A) 的情况。我们需要在狄拉克方程中进行局部规范变换。动量替换: 将动量 p 替换为动力学动量(Kinetic Momentum):ppqA,其中 A 是磁矢势。能量替换: 将能量算符 it 替换为:ititqΦ,其中 Φ 是电标势。最终得到的方程为:

iψt=[cα(pqA)+βmc2+qΦ]ψ

为了更直观地处理四维时空,我们通常将其写成协变形式。引入 Gamma 矩阵定义:γ0=β,γi=βαi。结合四维动量算符 pμ=iμ=(i1ct,i) 和四维势 Aμ=(Φ/c,A),方程整理为:

γ0(i1ctqΦc)ψ+γ(i+qA)ψmcψ=0

再利用协变导数 Dμ=μ+iqAμ,可以写成极简形式:

(iγμDμmc)ψ=0γμ(pμqAμ)ψ=mcψ

我们现在希望求解这个四分量旋量方程。Gamma 矩阵的简单分块形式使得我们可以把四分量两两分组。同时,因为在非相对论极限下,静止能量 mc2 是主导项,我们将波函数写成“大分量” ϕ 和“小分量” χ 的组合(均为 2 分量旋量),并显式提取出静能随时间的演化项:

ψ(r,t)=(ϕ(r,t)χ(r,t))eimc2t/

对左边的时间导数项应用乘法法则:

iψt=(it(ϕχ)+mc2(ϕχ))eimc2t/

代入狄拉克方程,并使用狄拉克标准表示下的矩阵形式,方程两边消去指数项,得到:

it(ϕχ)+mc2(ϕχ)=c(0σσ0)π(ϕχ)+(mc200mc2)(ϕχ)+qΦ(ϕχ)

其中 π=pqA 是动力学动量。我们将上式分解成两行耦合方程:

iϕt+mc2ϕ=c(σπ)χ+mc2ϕ+qΦϕiχt+mc2χ=c(σπ)ϕmc2χ+qΦχ

在非相对论极限下,电子的动能和势能远小于静能,即 |itχ||2mc2χ||qΦχ||2mc2χ|。于是小分量方程近似为:

2mc2χc(σπ)ϕχσπ2mcϕ

这说明了为什么 χ 被称为“小分量”:它的大小约为大分量的 v/c 倍。将 χ 的表达式代回大分量方程:

iϕt=[(σπ)22m+qΦ]ϕ

这是自旋出现的关键点。利用泡利矩阵恒等式 (σA)(σB)=AB+iσ(A×B)

(σπ)2=ππ+iσ(π×π)

计算算符叉乘 π×π(作用于测试函数 f):

(π×π)f=(pqA)×(pqA)f=(p×pqp×AqA×p+q2A×A)f

其中 p×p=0A×A=0。剩余项为:

q(p×A+A×p)f=q[(i)×(Af)+A×(if)]=iq[(×A)fA×(f)+A×(f)]=iq(×A)f=iqBf

因此:

(σπ)2=π2+iσ(iqB)=π2qσB

最终我们完整推导出了包含塞曼项(Zeeman term)U=μB泡利方程 (Pauli Equation)

iϕt=[(pqA)22m+qΦq2mσB]ϕ

显然其哈密顿量中多了最后一项:

H=(pqA)22m+qΦq2mσB

这就是自旋的由来。如果我们没有引入矩阵 σ(即自旋为 0 的克莱因-戈登情况),那么算符项就是简单的 π2,这一项不会产生与磁场 B 的线性耦合,会退化到正常的无自旋时的电磁耦合哈密顿量 HSpin0=(pqA)22m+qΦ。正是因为电子波函数的相对论协变性要求它必须有矩阵结构,才导致了 π 的不相容性被转化为磁能项。简单来说:引入 σ 矩阵是为了在数学上给粒子增加一个“内部转动轴”;如果不引入,粒子就是一个在空间中彻底“各向同性”的质点。它定义了算符的不对易性,在自旋为 0 的世界里,动量的三个分量 πx,πy,πz 虽然在磁场下不对易,但它们的作用是“外部”的(改变粒子的轨道),引入 σ 矩阵后,它把这种动量的不对易性(磁场强度 B)直接耦合到了波函数的内部维度上。

5. SU2 and SO3 之间的万有覆盖

以上我们已经从动力学(Dynamics)层面看到了自旋的“涌现”。但是自旋作为一种内秉性质,按理讲应该是无需借助外界电磁场来体现的,本质上它应该与对称性与变换有所关联所以才会存在。为什么波函数必须由 σ 矩阵来操作?为什么它是二分量的旋量?这需要我们深入到运动学(Kinematics)的底层——即时空对称性的群论基础。即使不考虑电磁场,电子作为一个物理对象,当我们把实验室的坐标系旋转一个角度 θ 时,它的波函数 ψ 必然会发生改变。正是这种“旋转下的变换规则”定义了自旋。在进入物理之前,我们先简要介绍一下最重要的 SU2 与 SO3 群,他们属于李群,在物理中占有极其重要的地位。

SU2 的基本定义是:

SU(2){U|UGL(2,C),UU=12×2,|U|=1}{[abba]|a,bC,|a|2+|b|2=1}{U(n,ω)=eiω2nσ|ω[0,π],n 为全体三维实单位矢量 }

若采用实参数 xiR,令 a=x4+ix3,b=x2+ix1 来描述:

U=[abba]=[x4+ix3x2+ix1x2+ix1x4ix3]

则约束条件变为 x12+x22+x32+x42=1 说明 SU2 作为流形是 S3,一个三维球面,属于超球面,其 T2-fibration 描述为:

{x1=sinθcosφx2=sinθsinφ,{x3=cosθcosχx4=cosθsinχ, 其中 θ[0,π/2];φ,χ[0,2π]

这可以用两个反向圆锥面 χ,φ 与一根轴描述 θ, 在 θ=0θ=π/2 时会有一个参数失效。也可以用一个中间会随着 θ 放缩于端点处使某参数失效的甜甜圈来描述,θ=0 时甜甜圈就是一个无宽度圆,只有 χ 沿着圆有效,φ 因为圆没有宽度而失效,而 θ=π/2 时甜甜圈变成了一个球,不再有空洞,所以 χ 失效,只有 φ 生效。

球坐标系描述为:

 对 ω[0,2π],θ[0,π],φ[0,2π] 有 {x1=sinω2sinθcosφx2=sinω2sinθsinφx3=sinω2cosθx4=cosω2

SU2 还可以由 Pauli 矩阵表达:

U=[x4+ix3x2+ix1x2+ix1x4ix3]=x412×2+ix1σ1+ix2σ2+ix3σ3.

进一步结合球坐标有:

U(n,ω)=eiω2nσ=12×2cosω2+inaσasinω2n=(sinθcosφ,sinθsinφ,cosθ);ω[0,2π],θ[0,π],φ[0,2π].

SO3 的基本定义是:

SO(3){R|RGL(3,R),RTR=13×3,|R|=1}{R(ω)ω=ωn,n=(cosφsinθ,sinφsinθ,cosθ)ω[0,π],θ[0,π],φ[0,2π]}.

SO3 很直观,就是我们看得见摸得着的旋转操作。作为流形,SO3 可以看作是 ω 端点构成的半径为 π 且对径认同的实心球,其抽象就在于对径认同,虽然其来源与绕固定轴逆时针转 180 度和顺时针转 180 度是一样的这样一个简单事实,但是对径认同导致了这个实心球虽然是一个连通流形但不是一个单连通流形,即空间内任意一条闭合曲线(环路)都能连续地收缩到一个点。对径认同的实心球有它的名字, 叫三维实射影空间 (real projective space), 记作 RP3

数学上,群 G 在向量空间 V 上的表示,是一个从群 G 到一般线性群 GL(V)(即 V 上所有可逆变换构成的群)的同态映射 g1,g2G,D(g1g2)=D(g1)D(g2),对于李群来讲还要求映射连续。而射影表示则是因为在量子力学中,物理状态由希尔伯特空间中的射线(Ray)描述,即 |ψeiα|ψ 代表同一物理状态。因此,群的乘法法则只需要在“相差一个相位因子”的意义下成立:

D(g1)D(g2)=ω(g1,g2)D(g1g2)

其中 ω(g1,g2) 是一个模为 1 的复数(相位因子)称为群指数(group exponent)。Bargmann定理(1954)为射影表示提供了严格的数学框架:对于满足 H2(g,R)=0 的李群 G(包括 SO3 和 Lorentz 群),其所有连续射影酉表示均可"提升"为某个中心扩张群 G~ 的普通酉表示。对于 SO(3) 这种通过简单的“相位消除”无法变成普通表示的情况,我们可以通过寻找它的万有覆盖群 SU(2),将 SO(3) 的射影表示转化为 SU(2) 的普通表示。

这里引入的万有覆盖群的概念点出了 SO3 和 SU2 深刻的联系。在拓扑学中,一个空间 X 的万有覆盖空间 X~ 就像是它的“升级版”,具有两个核心特征:单连通性: X~ 里面所有的圈都能缩成一个点(没有拓扑洞)。局部一致性: 在局部区域,X~ 看起来和 X 一模一样,但从整体上看,X~ 往往比 X 更“大”,并以 n:1 的方式覆盖在 X 上。

为什么 SU(2) 的矩阵能产生 SO(3) 的旋转?这里给出一个经典的构造方法。我们把三维空间中的向量 x=(x,y,z) 映射为一个二阶无迹厄米矩阵 X

X=xσ1+yσ2+zσ3=(zxiyx+iyz)

(其中 σi 是泡利矩阵)。注意,det(X)=(x2+y2+z2)=x2。让 SU(2) 中的矩阵 U 通过以下方式作用在 X 上:X=UXU。由于 U 是幺正的且行列式为 1,这个变换保持了 X 的无迹性和厄米性,并且保持了行列式不变det(X)=det(UXU)=det(X)。这意味着 x2=x2,即这个变换保持了向量的长度,因此它描述的是一个三维旋转。观察变换公式 X=UXU。如果你把 U 换成 U(U)X(U)=(1)2UXU=UXU 会发现 UU 产生了完全相同的旋转效果。这就是 2:1 覆盖的代数根源,每一个 SO(3) 中的旋转,在 SU(2) 中都有两个点对应。这解释了为什么旋转 2π 回不到原点(在 SU(2) 中只跑了半圈),而旋转 4π 才能回到原点(在 SU(2) 中跑了一整圈)。。也就是说,SO(3) 本身是有“洞”的(它的基本群是 Z2),而SU(2)(即 S3)是单连通的(基本群为 0),没有任何拓扑洞。

而万有覆盖中的第二个条件局部一致则引出了 SU2 和 SO3 的另一大特点,在单位元附近是局部同构的,意思就是:如果你只看“无穷小”的旋转,或者只看旋转一点点的情况,这两个群是完全一模一样的,没有任何区别。只有当你转了很大一圈(比如转了 2π)去探索群的“全貌”时,你才会发现它们其实不一样(一个回到了起点,另一个到了 I)。从数学上严格来讲,是因为它们拥有完全相同的李代数(Lie Algebra),即单位元的切空间同构:su(2)so(3)。在求解薛定谔方程的过程中我们知道 SO3 的 3 个生成元是 Jx,Jy,Jz(对应绕 x, y, z 轴的微小旋转),它们的对易关系(Commutation Relation)是:[Ji,Jj]=iϵijkJk,这是定义三维旋转本质的公式,所谓生成元也就是李代数的一个基底。

而对于 SU2 群,我们现在来求解其生成元。假设一个微小的变换:

U(ϵ)=IiϵS

还需要属于 SU2,那么有幺正性约束: (I+iϵS)(IiϵS)=IIiϵ(SS)=IS=SS 必须是厄米矩阵,物理上可观测。还有特殊性约束:利用公式 det(eA)=eTr(A)det(U)=det(eiϵS)=eiϵTr(S)=1Tr(S)=0S 必须是无迹矩阵。符合条件的正好就是 Pauli matrices,σx,σy,σz 构成了 SU(2) 李代数 su(2) 的一组完备基底。所以,生成元 S 必然正比于 σ,而 [σi,σj]=2iϵijkσk,我们会发现这与 [Li,Lj]=iϵijkLk 就差了一个系数 12,这就已经说明 SU2 和 SO3 的李代数同构了,而如果如果取 S=12σ,则可以变成标准的 [Si,Sj]=iϵijkSkJ 的对易关系完全相同,这样我们会获得一个自恰的角动量理论,物理世界里存在总角动量 L=J+S。如果自旋 S 想要有资格被称为“角动量”,并且能和 J 加在一起构成一个守恒量,那么 S 必须遵循和 J 一模一样的代数规则,而轨道角动量 J=r×p 是由空间坐标定义的,它的对易关系是没得选的(由 xp 的对易关系推导出来),所以只能让 S=12σ。除了理论上的自洽,现实世界的实验结果也是这样,当我们去做施特恩-格拉赫实验(Stern-Gerlach Experiment)时,我们测量电子在磁场中的偏转,测量到的物理值是 ±12。这就直接证明了,代表物理观测量的算符 S,它的特征值必须是 ±1/2。 而在数学上,只有矩阵 12σz 的特征值才是 ±1/2(因为 σz 的特征值是 ±1)。

回到表示,对于李群来说,生成元的一个很好的性质是,任意有限变换 D(θ) 都可以由生成元 J 通过指数映射得到。如果 J 是李代数 g 的元素,那么群中的元素可以写成:

D(θ)=exp(iθnJ)

这里的 J 就是我们前面讨论的角动量算符(矩阵)。注意这个式子涉及到了物理学符号中一种常见的“混用”或者说“简写”。严格数学上,指数映射 exp 将抽象的李代数元素映射为抽象的李群元素。 但是,物理公式 D(θ)=exp(iθnJ) 实际上是在表示空间(矩阵空间)内部进行的操作。这里的 J 已经是生成元的矩阵表示,所以算出来的结果 D(θ) 自然就是群元素的矩阵表示。

不管我们是在处理 SO(3) 还是 SU(2),它们的李代数是同构的。这意味着它们共享同一套生成元对易关系:[Ji,Jj]=iϵijkJk(这里取 =1)。我们要寻找这套代数规则允许的所有有限维不可约表示。定义阶梯算符:

J±=Jx±iJy

引入 Jz 的本征态 |j,m,使得:

Jz|j,m=m|j,m,J2|j,m=λ|j,m

计算对易子:

[Jz,J±]=[Jz,Jx]±i[Jz,Jy]=iJy±i(iJx)=±(Jx±iJy)=±J±

这意味着 J± 是本征值的“升降梯”:

Jz(J±|j,m)=(J±Jz+[Jz,J±])|j,m=(m±1)(J±|j,m)

如果 m 是特征值,那么 m±1 也是特征值。但是因为我们寻找的是有限维表示,特征值谱必须有上界 mmax 和下界 mmin

J+|j,mmax=0,J|j,mmin=0

利用算符恒等式 JJ+=J2Jz2Jz,作用在最高权态 |j,mmax 上:

0=(λmmax2mmax)|mmaxλ=mmax(mmax+1)

为了方便,我们将最大权重标记为 j,即 mmaxj。所以 Casimir 算符的特征值为 j(j+1)。同理,利用 J+J=J2Jz2+Jz 作用在最低权态 |j,mmin 上:

0=(j(j+1)mmin2+mmin)|mmin

解方程 mmin2mminj(j+1)=0,得到两个解:

mmin=jmmin=j+1

因为 mminmmax=j,所以只能取 mmin=j。从 mmin=j 爬升到 mmax=j,每次加 1,必须经过整数步 k 才能到达:

mmaxmmin=j(j)=2j=k(kZ)j=k2

我们获得了结论,仅凭李代数结构,允许的 j 值为 0,1/2,1,3/2,2。但是李代数结构只是局部性质,现在我们要把这些结果放回全局群结构中去检验。检验的核心标准是单值性: 如果我们把群元素沿着一条闭合路径转一圈回到起点(单位元),它的表示矩阵必须也回到单位矩阵(对于普通表示)。SO(3) 是三维空间中的旋转群,绕任意轴(设为 z 轴)旋转 2π360),物理空间完全复原:Rz(2π)=Rz(0)=1 这是群单位元。而对于 SO(3) 的普通表示 D,必须满足:D(Rz(2π))=D(1)=I 这是单位矩阵。代入李代数导出的公式,在 z 轴基底下,Jz 是对角的,对角元为 mD(2π)=exp(i2πJz)=diag(ei2πm,),要让这个矩阵等于单位阵 I,必须每一个对角元都等于 1:ei2πm=1mZ,如果 j 是整数(0,1,),则 m 是整数,条件满足,而如果 j 是半整数(1/2,3/2,),则 m 是半整数,ei2πm=11。 因此,SO(3) 的普通表示中,严格禁止半整数自旋。

但是 SU(2) 的几何结构与 SO(3) 不同。它是 SO(3) 的万有覆盖群(2:1 覆盖)。 在 SU(2) 中,参数 θ=2π 对应的群元素不是单位元,而是 U(2π)=II,只有旋转 4π 才是单位元。我们会发现 D(2π) 的行为完美匹配了 SU(2) 群本身的行为,我们实际上获得了 Bargmann 定理在 SO3 和 SU2 上的实例,SO3 非单连通李群的射影表示等价于其万有覆盖群 SU2 的普通表示,最终的映射关系表为:

自旋 j 在李代数中 在 SU(2) 中 在 SO(3) 中 物理粒子
整数 (0,1,) 存在 普通表示 (但不忠实,无法区分 ±I) 普通表示 玻色子 (光子等)
半整数 (1/2,) 存在 普通表示 (忠实表示) 射影表示 (多值) 费米子 (电子等)

“自旋”之所以能从这套抽象的数学结构中“涌现”出来,是因为量子力学对“物理状态”的定义比经典力学更宽容,从而释放了被经典物理屏蔽掉的拓扑自由度。

现在我们具体计算一下不同的 j 的表示是什么。我们要利用指数映射

D(j)(n^,θ)=k=0(iθ)kk!(n^J(j))k

其中 n^ 是旋转轴单位向量,θ 是旋转角。

最简单的,j=0 标量表示,维数只有 d=2(0)+1=1 维,基底只有一个态 |0,0,生成元因为 m 只能取 0,所以 Jz=[0]。 升降算符作用在最高/最低权态上都是 0,所以 J+=[0],J=[0],所以 Jx=0,Jy=0,Jz=0。指数映射得到表示 D(θ)=eiθn0=1,也称为平凡表示。这就是标量,不管怎么旋转,数值永远乘以 1,保持不变。

j=12,维数为 2j+1=2。我们已经知道生成元为 J=12σn^J=12(n^σ),为了计算 (n^J) 的高次幂。 回顾泡利矩阵性质:(n^σ)2=I(单位阵)。 因此生成元的幂次规律为:

(n^J)2=(12n^σ)2=14(n^σ)2=14I,(n^J)3=(n^J)2(n^J)=14(n^J)

其通项公式为:

(n^J)2k=(14)kI=(12)2kI(n^J)2k+1=(12)2k(n^J)=(12)2k+1(n^σ)

于是做级数求和,将指数级数拆分为偶数部分和奇数部分:

D(1/2)=k=0(iθ)kk!(n^J)k=m=0(iθ)2m(2m)!(12)2mI偶数项+m=0(iθ)2m+1(2m+1)!(12)2m+1(n^σ)奇数项

偶数项系数:(1)m(2m)!(θ2)2m=cos(θ2)。奇数项系数:i(1)m(2m+1)!(θ2)2m+1=isin(θ2)。最终得到:

D(1/2)(n^,θ)=cos(θ2)Iisin(θ2)(n^σ)

这就是 j=1/2 的表示。它把旋转操作映射到了 2×2 复矩阵检验 2π 代入 θ=2πcos(π)=1,sin(π)=0。结果是 I

j=1,维数为 2j+1=3。我们需要 3×3 的矩阵。在物理学定义的角动量基底(笛卡尔基底)下,生成元满足 (Jk)ab=iϵkab。例如绕 z 轴旋转的生成元 Jz

Jz=(0i0i00000)

对于任意轴 n^,令矩阵 K=n^J。直接计算 Jz 的幂(其他方向同理):

Jz2=(0i0i00000)(0i0i00000)=(100010000)(注意:这不是 I)Jz3=Jz2Jz=(100010000)(0i0i00000)=(0i0i00000)=Jz

发现规律,对于 j=1 的生成元,满足特征方程 (n^J)3=(n^J)。这意味着:奇数项 (k=1,3,5):(n^J)k=(n^J),偶数项 (k=2,4,6):(n^J)k=(n^J)2k=0 项:I (单位阵)。同样展开泰勒级数,但这次要把 I 单独拿出来,因为 J2I

D(1)=I+odd k(iθ)kk!(n^J)+even k2(iθ)kk!(n^J)2

奇数项系数:i(θθ33!+)=isinθ,偶数项系数:(θ22!+θ44!)=cosθ1。结果:

D(1)(n^,θ)=Iisinθ(n^J)+(cosθ1)(n^J)2

这就是 j=1 的表示(物理学形式的罗德里格斯旋转公式)。它把旋转操作映射到了 3×3 实矩阵(虽然 J 包含 i,但 iJ 是实矩阵)。检验 2π 代入 θ=2πsin(2π)=0,cos(2π)=1D(1)=I0+(11)()=I

一般地,我们可以建立所有 j 的生成元和表示。所有 j 的生成元构建都依赖于量子力学中角动量算符的三个核心矩阵元公式。只要有这三个公式,我们就能写出 j=0,3/2,2 甚至 j=100 的矩阵。首先 Jz 是对角矩阵:

j,m|Jz|j,m=mδmm

然后 J+ (升算符) 是次对角矩阵:

j,m+1|J+|j,m=j(j+1)m(m+1)

J (降算符) 是下次对角矩阵: 它是 J+ 的转置(实数情况下)。而 JxJyJ± 组合而成:

Jx=12(J++J),Jy=12i(J+J)

比如 j=3/2:自旋 3/2 表示 (Spin-3/2) 这属于费米子,和电子类似,但它有 4 个分量。常见于 Δ 重子或超引力中的引力微子。维数:d=2(3/2)+1=4 维。构建生成元 Jz (对角线):

Jz=(3/200001/200001/200003/2)

J+ (升算符系数) 则需要计算 j(j+1)m(m+1),其中 j=3/2,也就是 3.75m(m+1)m=1/23/2: 3.750.75=3m=1/21/2: 3.75(0.25)=4=2m=3/21/2: 3.750.75=3。所以:

J+=(0300002000030000)

利用 Jx=12(J++J+)

Jx=12(0300302002030030)

表示:这是一个 4×4 的幺正矩阵。旋转 2π 时,因为对角元是半整数,它会变成 I4×4。所以 j=3/2 也是 SU(2) 的忠实表示,SO(3) 的射影表示。

至此,我们从对称性与群论的视角完整讲述了自旋的本质,不过似乎并不像上一章一样需要相对论修正,但是不像上一章我们从 Dirac 方程出发,直接推出了 4 分量波函数。我们这里的从从空间旋转(SU(2)/SO(3))出发,只推出了 j=1/2 的 2 分量旋量(Pauli 旋量),另两个分量去哪里了?其原因同样在于相对论效应,目前的对称性分析只考虑了空间旋转,没有考虑洛伦兹变换(Boost)。只有引入洛伦兹群(Lorentz Group),才能解释为什么电子必须是“左右手”两个 SU(2) 表示的直和(2+2=4),从而完美闭环回到 Dirac 方程的结构。

6. 洛伦兹群

为了获得完整的自旋,我们必须开始考虑真正的对称性群,是包含了旋转(Rotation)和伪转动(Boost)的洛伦兹群 SO(1,3)。当我们试图寻找洛伦兹群的“基本表示”时,会发生一件极其奇妙的事情:代数结构裂开了。

洛伦兹群的定义为:

O(1,3){ΛΛGL(4,R),gμνΛμρΛνσ=gρσ}dimO(1,3)=6,g=diag(1,1,1,1)

本质上是闵氏时空的保度规线性变换群,从保度规条件 ΛTgΛ=ggμνΛμρΛνσ=gρσ 可以得到一条对分量 Λ00 的限制:

1=gμνΛ0μΛ0ν=(Λ00)2i(Λ0i)2(Λ00)2=1+i(Λ0i)21

这就是说洛伦兹变换必有 Λ1Λ1,所以已经不连通了,由此可以将洛伦兹群氛围 O+(1,3)O(1,3) 两个流形,前者称为时向 (orthochronous) 洛伦兹群,后者不存在恒等元所以不构成群,称为洛伦兹群的反时向 (antichronous) 分支。

从保度规条件还可以确定行列式的取值:

|ΛTgΛ|=|g||Λ|2|g|=|g||Λ|2=1, 即有 |Λ|=±1

其中 |Λ|=1 的记作 SO(1,3) 称为固有 (proper) 洛伦兹群,而 |Λ|=1 就称之为洛伦兹群的非固有 (improper) 分支。结合上述两点考虑我们可以将洛伦兹群 O(1,3) 分成四块儿连通流形,但实际上我们只需要研究固有保时性分支 SO+(1,3),这是因为另外三块儿都可以通过两个确定的洛伦兹变换作用在 SO+(1,3) 上来得到,更重要的是,现实世界的参考系变换绝对是保时向且固有的,而这两个确定的洛伦兹变换分别为:时间反演变换 T=T1=diag(1,1,1,1) 和宇称变换 P=P1=diag(1,1,1,1)

我们关注固有保时洛伦兹群 SO+(1,3)。在这个连通分量中,任何一个变换都可以写成单位元的指数映射。如同 SO(3) 有 3 个转动生成元,SO+(1,3) 共有 6 个自由度(3 个转动 + 3 个伪转动),对应 6 个生成元。考虑无穷小变换 ΛIiϵX,类似上面的过程,我们可以写出两组生成元:转动生成元 J=(J1,J2,J3):对应空间转动,它们就是我们熟悉的角动量算符;伪转动(Boost)生成元 K=(K1,K2,K3):对应沿 x,y,z 轴的速度变换。这 6 个生成元满足的洛伦兹李代数 so(1,3) 如下:

[Ji,Jj]=iϵijkJk [Ji,Kj]=iϵijkKk [Ki,Kj]=iϵijkJk

此时,这个代数结构看起来还是耦合在一起的(JK 互相缠绕)。为了找到不可约表示,我们引入一个非幺正的基底变换(Complexification)。定义两组新的算符 N+N

N+=12(J+iK),N=12(JiK)

让我们计算这两组新算符的对易关系。 首先看 N+ 内部:

=14[Ji+iKi,Jj+iKj]=14([Ji,Jj]+i[Ji,Kj]+i[Ki,Jj][Ki,Kj])=14(iϵijkJk+i(iϵijkKk)+i(iϵijkKk)(iϵijkJk))=14(2iϵijkJk2ϵijkKk)=iϵijk12(Jk+iKk)=iϵijkNk+

同样地,我们可以验证 [Ni,Nj]=iϵijkNk。最令人震惊的结果在于 N+N 之间:

[Ni+,Nj]=14[Ji+iKi,JjiKj]==0

这意味着:洛伦兹群的李代数在复化后,分裂成了两个相互独立的 su(2) 代数的直和。

so(1,3)Csu(2)Lsu(2)R

这在群论上是一个巨大的简化。既然我们已经对 su(2) 的表示了如指掌(由自旋 j 标记),那么洛伦兹群的不可约表示就可以由一对半整数或者整数 (jL,jR) 来唯一标记。根据这种分解,最基本的旋量表示不再是唯一的,而是有了两种最基础的选择(基本表示),它们分别对应其中一个 su(2)j=1/2,另一个取 j=0。这就引出了手性的概念。

左手外尔旋量 (Left-handed Weyl Spinor) 对应标记 (1/2,0)。它对 N 表现为自旋 1/2,对 N+ 表现为标量。这是一个 2 分量复向量,记作 ψL。在此表示下,N=12σN+=0。由此反解出物理生成元:

J=N++N=12σ,K=i(N+N)=i12σ

右手外尔旋量 (Right-handed Weyl Spinor) 对应标记 (0,1/2)。它对 N 表现为标量,对 N+ 表现为自旋 1/2。这也是一个 2 分量复向量,记作 ψR。在此表示下,N=0N+=12σ。物理生成元为:

J=12σ,K=i12σ

注意这里 K 的符号差异!这说明虽然 ψLψR 在空间旋转 (J) 下行为一样(都是自旋 1/2),但在 Lorentz Boost (K) 下,它们的变换性质截然相反。

既然 ψLψR 都是 2 分量的,为什么我们需要 4 分量?原因在于宇称 (Parity)。宇称变换 P 将空间坐标反演 xxJ 是轴矢量(r×p),在 P 下不变:JJK 是极矢量(v),在 P 下变号:KK。将此代入 N± 的定义,我们发现宇称变换互换了这两个代数:

P:N+N

这意味着宇称将左手表示 (1/2,0) 变成了右手表示 (0,1/2)。如果我们要描述一个像电子这样既有自旋又有质量、且遵循宇称守恒(在电磁力下)的粒子,我们不能只取其中一个。我们必须把它们“直和”在一起。于是,Dirac 旋量 Ψ 作为 SO+(1,3) 扩张了宇称算符后的表示,正是这两个基本表示的直和:

Ψ=(ψLψR)(12,0)(0,12)

这就解释了为什么是 4 分量:两个分量来自左手扇区,两个分量来自右手扇区,它们通过质量项和宇称变换紧密耦合在一起。 我们之前用 SU(2) 看到的 Pauli 旋量,仅仅是这个相对论性对象在静止系(或非相对论极限)下的一个侧影。

顺便提一下,这里的 Weyl 旋量和刚才提到的 gamma matrices 中的 γ5 关联深刻,实际上,厄米算符 γ5iγ0γ1γ2γ3 就是那个用来“识别”和“定义” Weyl 旋量的算符。如果没有 γ5,我们就无法在数学上区分什么是“左手(Left-handed)”,什么是“右手(Right-handed)”。在 Dirac 理论中,γ5 被称为手征算符(Chirality Operator)。它有一个至关重要的代数性质:它与所有的 γμ 反对易 {γ5,γμ}=0。但它与洛伦兹变换的生成元 Sμν=i4[γμ,γν] 对易[γ5,Sμν]=0,这意味着 γ5 是洛伦兹群表示的一个守恒量(对于无质量粒子),我们可以用它的本征值来给旋量分类。右手 Weyl 旋量是 γ5 本征值为 +1 的态 γ5ψR=+ψR,左手 Weyl 旋量是 γ5 本征值为 -1 的态 γ5ψL=ψL,所以,物理上说的“左手性”和“右手性”,在数学上指的就是 γ5 的特征值是 -1 还是 +1。既然 Dirac 旋量 Ψ 是左手和右手的混合体(直和),我们如何从一个混杂的 Ψ 中把左手部分和右手部分单独“筛”出来呢?这就需要用到基于 γ5投影算符

PL=1γ52,PR=1+γ52

这两个算符具有投影算符的标准性质(P2=P,PLPR=0,PL+PR=1)。它们的作用就是杀掉一手分量,只留下另一手分量:

PLΨ=1γ52(ψL+ψR)=1(1)2ψL+112ψR=ψL,PRΨ=ψR

在粒子物理(特别是弱相互作用)的计算中,你经常会看到 1γ52 这样的项,这就是在告诉你:“这个相互作用只和左手旋量玩,右手旋量请走开。”(这正是宇称不守恒的数学表达)。为了让这种关系一目了然,我们可以选择一组特殊的 Gamma 矩阵形式,称为Weyl 表象(或手征表象),这区别于之前提到的 Dirac 表示,在 Weyl 表象下,Gamma 矩阵是对角分块的,这使得 Dirac 旋量显式地分裂为上下两个 Weyl 旋量 Ψ=(ψL,ψR)T,这是高能物理中最常用的视角;而在低能凝聚态物理中,我们常用 Dirac(标准)表象,那里 ψLψR 深度混合,更能体现‘大分量’与‘小分量’的非相对论近似。

7. 磁梯度力

我们已经完成了一段漫长的数学旅程。从狄拉克方程的推导到洛伦兹群的表示,我们已经确立了电子必须具有自旋,并且是一个 4 分量的相对论性客体。现在,让我们回到最初的那个悖论:如果不做功的洛伦兹力不能解释磁铁吸铁,那到底是谁在做功?为了回答这个问题,我们需要将微观的自旋与宏观的力联系起来。

通过对狄拉克方程进行非相对论近似,我们自然地得到了泡利方程。让我们重新审视那个凭空出现的“塞曼项”(Zeeman term):

HZeeman=e2m(σB)

在经典物理中,我们定义磁矩 μ 与角动量 L 的关系为回磁比。对于轨道角动量,这一关系为:

μL=e2mL

如果我们试图用同样的逻辑来定义“自旋磁矩”,我们需要将自旋算符 S 与磁矩联系起来。回顾 Sec 5 中自旋算符的定义:

S=2σσ=2S

σ 代回塞曼项 HZeeman,我们得到:

HZeeman=e2m(2S)B=2e2mSB

我们将磁势能的一般形式写作 U=μSB。对比上式,我们可以读出电子的自旋磁矩 μS

μS=2e2mS

如果将其写成通用的朗德因子(Landé g-factor)形式 μ=ge2mS,我们可以直接得到结论:

g=2

这个 g=2 不是实验凑出来的参数,它是狄拉克方程时空对称性的直接数学结果。它意味着:电子自旋产生磁矩的效率是经典轨道运动的两倍。

有了磁矩 μ,我们终于可以解释“谁在做功”了。经典洛伦兹力 FLorentz=q(v×B) 确实不做功。但是,对于一个具有内秉磁矩的物体,其动力学由势能 U 决定。根据哈密顿力学,力是势能的负梯度:

F=U=(μB)=(μB)

由于自旋 μ 是内秉属性,在空间微分中保持常数,我们得到:

FGradient=(μ)B

这就是那个做功的力,这个力被称为梯度力(Gradient Force)。如果磁场是均匀的 (B=0),则受力为 0,只有力矩。而现实中的磁铁产生的磁场是非均匀的(磁感线发散),B0。这是一个保守力,它将磁场与磁矩耦合的势能转化为物体的动能。所以,磁铁吸铁,本质上是量子化的自旋磁矩在非均匀磁场中受到梯度力的作用。洛伦兹力负责偏转,梯度力负责做功。

所以真的有两种磁力,而从对称性的角度,这对应着两种完全不同的对称性机制,这已经触及了量子场论的边缘。在量子场论中,相互作用是由对称性规定的:

既然我们已经站在了量子场论的视角上,必须指出 Dirac 的 g=2 虽然辉煌,但并非终极真理。在 Dirac 方程中,我们将电子看作一个经典的场在与电磁场耦合。但在完全的量子电动力学(QED)中,真空不是空的。电子在传播过程中,会不断地发射和吸收虚光子(Virtual Photons),甚至产生正负电子对。这意味着电子与磁场的相互作用顶点(Vertex),不再仅仅是一个简单的点(Tree level),而是包含了无穷多的圈图修正(Loop corrections)。朱利安·施温格(Julian Schwinger)在 1948 年计算了第一阶修正(上图所示的单圈图),给出了著名的公式:

g=2(1+α2π+O(α2))

其中 α1/137 是精细结构常数。 这使得 g 的理论值变为 g2.002319...,这与实验测量的结果惊人地一致(精确到小数点后 12 位)。这个微小的偏差(Anomalous magnetic moment)不仅证实了 g=2 的相对论起源,更揭示了磁性背后的深层物理:当我们感受到磁铁的吸力时,我们不仅见证了时空的几何属性(自旋),甚至触摸到了真空中虚粒子沸腾的海洋。

8. Heisenberg Model

我们现在知道每个电子都是一个小磁针(g=2)。但是,如果你把一堆小磁针放在一起,室温下的热骚动(Thermal Agitation)足以将它们的方向完全打乱,使得宏观磁矩为零(顺磁性)。要形成铁磁性,自旋之间必须存在一种极强的“耦合力”,强制它们整齐排列。经典磁偶极相互作用(Dipole-Dipole Interaction)太弱了,大约只有热能的万分之一。真正的力量来自于我们之前提到的全同粒子统计性(Spin-Statistics)库仑相互作用的结合。这被称为交换相互作用(Exchange Interaction)

为了展示本质,我们考虑最简单的模型:双电子系统(例如 He 原子或者是两个相邻的 Fe 原子上的电子)。假设有两个电子 1 和 2,以及两个空间轨道 ψa(r)ψb(r)ψa(r) 局域在原子 A 附近,ψb(r) 局域在原子 B 附近,并假设这两个轨道是正交归一的:ψa|ψb=0,注意这里的正交是我们接下来将要得到的 Heisenberg model 的前提假设,即使不正交归一我们也可以通过基底变换造出两个正交归一的新轨道,认为二者的重叠积分很小,这会使得双电子系统的计算大大简化,但是这并不意味着两个电子之间就没有交互了,因为交互作用中会伴随着交互作用项积分,一般不为零。

这个系统的总哈密顿量为:H=H0+Hint,其中 H0 是单电子部分(动能+核势能),Hint 是两个电子之间的库仑相互作用Hint=e2|r1r2|,注意:这里完全没有磁相互作用项,只有纯粹的静电斥力。

根据量子力学全同粒子费米统计假设,电子的总波函数 Ψ(1,2) 必须在粒子交换算符 P12 作用下变号(反对称):

P12Ψ(1,2)=Ψ(1,2)

因为总波函数由空间部分 ϕ(r1,r2)自旋部分 χ(s1,s2) 组成:Ψ=ϕχ。我们首先推导自旋部分,然后根据总体的反对称性得到空间部分。我们有两个自旋为 1/2 的粒子(比如两个电子),总共有 2×2=4 种可能的直积状态(未耦合基底),我们要把它们的自旋加在一起,看看总自旋 Stot=S1+S2 会是什么样子,找到总自旋算符 S^2 和总磁量子数 S^z 的共同本征态 |S,M。根据角动量加法规则,两个 1/2 自旋合成的总自旋 S 可以取:S=1/2+1/2=1 (三重态,有 3 个 M 值:+1,0,1)和 S=1/21/2=0 (单态,有 1 个 M 值:0)。

三重态对应三个分量。总磁量子数 M 是两个粒子磁量子数之和:M=m1+m2。 要让 M=1,唯一的可能是两个电子都向上:1/2+1/2=1。 所以,三重态的第一个成员是确定的:|1,1=|↑↑,然后使用降算符 (S) 得到中间态 (M=0) 从 |1,1 推导出 |1,0。 利用总降算符 S^=S^1+S^2 作用在态 |j,m 上:J|j,m=j(j+1)m(m1)|j,m1,对左边的耦合态作用得到:

S^|1,1=1(1+1)1(11)|1,0=2|1,0

对右边的直击态作用得到:

(S^1+S^2)|↑↑=(S^1|1)|2+|1(S^2|2) =|1|2+|1|2 =|↓↑+|↑↓

于是得到:

|1,0=12(|↑↓+|↓↑)

最低权态 (M=1) 也很简单,只有两个都向下才能得到 1|1,1=|↓↓。于是三重态自旋部分是对称的,交换都不变号,那么空间部分必须是反对称的。

单态的量子数是 S=0,M=0。 它必须是 M=0 的直积态 |↑↓|↓↑ 的某种线性组合:|0,0=a|↑↓+b|↓↑,因为不同量子数的本征态必须是正交的,单态 |0,0 必须与三重态中的 |1,0 正交,可以解出:

|0,0=12(|↑↓|↓↑)

单态(Singlet)的自旋部分是反对称的,对应的空间部分必须是对称的。总结为:

我们的目标是:寻找一个只包含自旋算符 SiSj 的数学表达式 H^eff,使得它作用在单态和三重态上时,能够自动给出对应的能量 ESET。要构建这个哈密顿量,最自然的积木块是两个自旋的点积 SiSj。我们需要先算出这个算符在单态和三重态下的特征值分别是多少。定义双电子系统的总自旋算符Stot=Si+Sj,对总自旋算符取平方反解出点积项:SiSj=12(Stot2Si2Sj2)。利用量子力学中角动量平方算符的特征值公式 S^2|s=s(s+1)|s(此处为了简洁暂时省略 2,最后结果不受影响,或者视为无量纲自旋):

我们假设等效哈密顿量 H^ij 具有如下线性形式(这是最普遍的旋转对称形式):H^ij=C0+C1(SiSj)C0 只是一个与自旋构型无关的常数能量平移(Energy Shift),在研究相变和自旋动力学时,我们可以将其丢弃(或者重新定义能量零点),然后将单态和三重态的数据代入,有 ES=C1(34),ET=C1(14),定义常数 JESET,就得到了双粒子哈密顿量:

H^ij=J(SiSj)

现在我们将这个双粒子相互作用推广到整个晶格。假设每个电子 i 只和它的最近邻(Neighbors)发生相互作用。我们需要对晶格中所有的原子求和。为了修正这种重复计数:

Hexchange=J2i,j neighborSiSj

最后,我们还要考虑每个电子自旋与外部均匀磁场 B 的相互作用。这是单体相互作用,不需要考虑邻居。回顾我们之前从 Dirac 方程导出的结论,电子具有自旋磁矩:

μS=ge2mS

为了让公式更加简洁和通用,物理学家定义了一个自然常数组合,称为玻尔磁子(Bohr Magneton)。玻尔磁子是原子物理中磁矩的自然单位。它的定义为:

μBe2me

这个物理量包含了三个基本常数:基本电荷 e、普朗克常数 和电子质量 me。它代表了一个经典电子在氢原子基态轨道上运动所产生的轨道磁矩大小,这是国际标准公制,而在高斯单位制下,还包含光速 μBe2mec。如果我们将自旋算符 S 视为无量纲算符(即特征值为 1/2 而不是 /2),那么真实的物理角动量是 S。将这个 提取出来与前面的常数结合:

μS=ge2me(Sdimensionless)=g(e2me)S=gμBS

当电子处于外部磁场 B 中时,其势能(Zeeman Energy)由经典电磁学公式 U=μB 给出。将上述磁矩表达式代入:

UZeeman=(μS)B(注意符号:通常势能定义为 μB)=(gμBS)B=gμBSB

关于符号的物理约定: 在凝聚态物理中,我们通常希望哈密顿量反映能量极小值。电子带负电,磁矩 μ 与自旋 S 反向。能量最低态是磁矩 μ 与磁场 B 平行。这意味着自旋 S 与磁场 B 反平行。为了避免处理繁琐的负号,或者为了让自旋看起来“顺着”磁场(定义 S 指向磁矩方向而非角动量方向),有时文献会调整定义。但最标准的推导(保留电子电荷负号)给出的塞曼项通常写为 +gμBSBμB。然而,在海森堡模型的习惯写法中,为了数学上的对称性及方便讨论(例如假设 g 为负或重新定义自旋方向),通常将塞曼项写为负号形式,表示自旋倾向于沿着场排列(这是一种唯象的处理):

HZeeman=gμBiSiB

(这意味着能量最低时,SiB 同向。这蕴含了我们重新定义了自旋方向,或者 g 取负值。在唯象模型中,我们只关心:磁场倾向于把自旋拉向哪个方向。)

现在,我们将两部分拼图合二为一:内部相互作用(由泡利原理和库仑力产生的交换能)与外部相互作用(由相对论量子效应产生的磁矩与外场耦合),将它们相加,我们终于得到了描述固体磁性最核心的哈密顿量——海森堡模型(Heisenberg Model)

H=J2i,jSiSjgμBiSiB

这个公式是现代磁学的基石。第一项(J)解释了为什么磁铁会有磁性(自发磁化,自旋整齐排列)。第二项(B)解释了磁铁如何被外界控制(磁化过程,磁滞回线)。μBg 将微观的量子常数(,e,me)与宏观的可观测磁场联系在了一起。

9. Ising Model

我们已经完成了微观机制的构建(Dirac Spin Exchange Interaction Heisenberg Hamiltonian)。现在,我们要从微观走向宏观,为了做到这一点,我们需要处理 Heisenberg 模型。但在二维或三维中精确求解 Heisenberg 模型极度困难(因为它包含非对易算符)。因此,我们需要引入Ising 模型作为近似,并使用平均场论(Mean-Field Theory)来展示对称性是如何破缺的。

我们对 Heisenberg model 进行无量纲化处理,并分析其结构。假设外磁场沿着 z 方向 B=(0,0,B)。我们将自旋算符 S 的点积展开为纵向(z)和横向(x,y)分量:

SiSj=SizSjz+(SixSjx+SiySjy)

为了看清横向部分的物理意义,我们引入升降算符(Ladder Operators):

Si+=Six+iSiy,Si=SixiSiy

于是,海森堡模型可以重写为两部分:

H=[J2i,jSizSjzhiSiz]Ising Part+[J4i,j(Si+Sj+SiSj+)]Flip Part

(其中 h=gμBB)。这两部分有着截然不同的物理含义:

在许多真实的磁性材料中,由于晶格结构的对称性,存在磁各向异性(Magnetic Anisotropy)。这意味着自旋在某些方向(如 z 轴,易轴)上的能量比在 x,y 平面内更低。如果各向异性足够强,或者我们只关心经典极限下的相变行为,我们可以忽略 Flip Part(量子涨落项),只保留纵向项。这便是著名的 Ising Model

此时,我们将算符 Siz 替换为经典变量 σi=±1(吸收系数到 J 中):

HIsing=J2i,jσiσjhiσi

这是一个巨大的简化:我们将一个非对易的量子矩阵问题,变成了一个经典的统计排列组合问题。2D Heisenberg model 至今是没有办法精确解的,我们可以计算方法求解。不过其中的 Ising model 相对简单一些,1D Ising model 的解很简单,在 1D 没有相变(Ising, 1924),但是 2D 也是很难解,直到有一个奇人 Lars Onsager 发了个很简单的论文说这个他解了,但是没有给出怎么解的,但是他给出了那个 critical temperature,且指出低温时铁磁,高温是顺磁,这是第一个可以精确解出的铁磁相变的 model,他一直没有公布解答过程,直到杨振宁看到这篇论文,给出了解答,不过一样很困难。为了直观理解物理图像,我们这里采用平均场近似

我们面临的问题是多体耦合:σi 的状态取决于 σj,而 σj 又取决于 σk... 这种连锁反应使得配分函数难以计算。平均场思想:当我们关注第 i 个原子时,我们不关心邻居 j 到底是在 +1 还是 1 之间跳变,我们只关心邻居的平均影响。我们将 σj 写为平均值加涨落:σj=σ+δσj。 忽略二阶涨落项 δσiδσj0,哈密顿量可以线性化为单体形式:

HMFA=iσi(Jjneighσ+h)heff

这就定义了有效分子场(Molecular Field) heff

heff=Jzσ+h

其中 z 是配位数(每个原子的邻居数量)。现在,问题变成了一个自旋在“外场” heff 中的统计分布。根据玻尔兹曼分布,该自旋向上的概率正比于 eβheff,向下的概率正比于 eβheffβ=1/kBT)。 于是,该自旋的热力学平均值 σi 为:

σi=(+1)eβheff+(1)eβheffeβheff+eβheff=tanh(βheff)

因为晶格是均匀的,σi 必须等于场源本身的平均值 m=σ。代入 heff,我们得到了著名的自洽方程(Self-consistent Equation)

m=tanh(Jzm+hkBT)

让我们考虑最关键的情况:没有外部磁场 (h=0)。方程简化为:

m=tanh(TcTm)

其中我们将常数打包定义为 居里温度(Curie Temperature) Tc=Jz/kB。这是一个超越方程,我们可以通过作图法(寻找 y=my=tanh(TcTm) 的交点)来分析解的行为:

这就是自发对称性破缺(Spontaneous Symmetry Breaking): 哈密顿量在 h=0 时具有 σσ 的翻转对称性。但是,当温度低于 Tc 时,大自然被迫在“全体向上”和“全体向下”中选择一个。这种选择不是由外力强加的,而是系统为了降低能量(由交换相互作用 J 驱动)而自发做出的集体决策。这正是磁铁吸铁这一宏观现象背后的统计力学本质。

10. Magnetic Domains

我们根据 Ising 模型和平均场论得出了结论:当温度低于居里温度 Tc 时,电子自旋会自发排列,产生巨大的宏观磁化强度 M。但是,这立刻引出了一个新的悖论:如果你去五金店买一颗铁钉(室温显然远低于铁的居里温度 1043K),它并没有磁性。它不会去吸别的东西。这是因为我们忽略了最后一项能量竞争。

我们之前的哈密顿量只考虑了交换能(Exchange Energy)和塞曼能(Zeeman Energy)。但在宏观尺度上,还有一个经典的静磁能(Magnetostatic Energy)。如果一块铁中的所有 1023 个原子都整齐朝上排列,这块磁铁会在周围空间建立一个巨大的磁场。磁场包含了能量密度 B2/2μ0。把所有磁力线撑开需要耗费巨大的能量。为了降低这个静磁能,材料会自发地分裂成许多微小的区域,称为磁畴(Magnetic Domains)。虽然在磁畴内部,交换相互作用让自旋整齐排列(满足微观铁磁性);但在整体上,各个磁畴的磁矩矢量和为零(Mi=0),外部没有磁感线,从而极大地降低了静磁能。

磁畴之间的边界被称为畴壁(Domain Wall)。在畴壁内,自旋不是突然反转的,而是渐进旋转的。这又是一次能量的博弈:交换能希望自旋平行,不想让它们转弯(希望畴壁越宽越好);磁各向异性(Anisotropy)希望自旋沿着易轴,不想让它们指向中间方向(希望畴壁越窄越好)。 这两者的平衡决定了畴壁的厚度(通常是几百个原子层)。磁畴的形成并不是像自旋那样来自“第一性原理”的推导,而是来自连续介质场论下的能量最小化(Energy Minimization)属于微磁学(Micromagnetics)的范畴,这里就不展开了。

现在,我们终于可以完整描述“磁铁吸铁”的宏观过程了:

Conclusion: The Deep Symmetry of the Universe

当你手中把玩两块磁铁,感受它们之间的排斥与吸引时,你感受到的不仅仅是力。你是在用手触摸量子力学的本质宇宙演化的秘密。让我们回顾这段旅程,看看我们是如何重建物理直觉的:

  1. 经典的崩塌:我们发现洛伦兹力不做功,且经典统计物理禁止磁性(Bohr-van Leeuwen 定理)。
  2. 相对论的修正:Dirac 方程揭示了电子必须是 4 分量的旋量,并携带 g=2 的内秉磁矩。磁性,是相对论效应在低速世界的残留。
  3. 量子统计的强力:泡利不相容原理结合库仑斥力,产生了一种等效的“交换相互作用”,迫使自旋平行排列。
  4. 对称性的破缺:Ising 模型告诉我们,当温度降低时,系统为了生存(降低能量),不得不破坏旋转对称性,选择一个方向。

最后,值得一提的是,我们在 Ising 模型中看到的自发对称性破缺(SSB),其意义远超固体物理。它是现代物理学理解宇宙的核心范式。在宇宙大爆炸之初(极高温),物理规律具有极高的对称性。所有的基本粒子都没有质量,就像高温下的铁块没有磁性一样(顺磁相)。 随着宇宙冷却,当温度跌破某个临界值时,充满宇宙的**希格斯场(Higgs Field)**发生了相变。就像电子自旋突然选择朝向一个方向一样,希格斯场在真空中获得了一个非零的期望值(Vacuum Expectation Value)。

所以,下一次当你看到磁铁吸起一枚回形针时,请意识到:你正在目睹一个微缩版的宇宙创生时刻。那个让铁钉获得磁性的机制,正是让你身体里的夸克和电子获得质量、让这个宇宙得以存在的同一机制。

磁力不做功,是时空的几何在做功;磁铁的吸力,是量子的幽灵在宏观世界的舞蹈。