T

1 HW 2思考题

1.1

这个张量怎么有点眼熟,算电磁和引力辐射时貌似会用到

1.1.1

Tij=δ(r1)xixjd3x=δ(r1)xixjr2drdΩ=04πδ(r1)(xir)(xjr)r4drdΩ=0r4δ(r1)drninjdΩ=ninjdΩ

1.1.2

在旋转变换下,向量 ni 变为 ni=Riknk,所以

Tij=(Riknk)(Rjlnl)dΩ=RikRjlnknldΩ=RikRjlTkl

1.1.3

见 ppt,C=4π3

1.1.4

第一个利用 dΩ=sinθdθdϕ 即可。

第二个,根据对称性,可以直接算 z 方向计算: =zrsinθdθdϕ=rcosθrsinθdθdϕ=0

第三个,看讲义

第四个,被积函数是奇函数,结果为零

第五个,积分结果必须是 SO(3) 下的不变张量,且阶数为四,它必须由克罗内克符号 δ 组合而成。所有可能的配对方式只有三种:

ninjnknldΩ=C(δijδkl+δikδjl+δilδjk)

注意,由于 ninjnknl 对索引 i,j,k,l 的任意排列都是全对称的,所以这三项前的系数 C 必须相等。先做两次缩并,左边:

(nini)(nknk)dΩ=11dΩ=4π

右边:

(δiiδkk+δikδik+δikδik)=9+3+3=15

所以 C=4π15.

1.2

1.3

用柱坐标爆算,注意:

E=e4πε0γR[s2+γ2(zvt)2]3/2

剩下的懒得写了。

2 HW 3思考题

2.1

2.1.1

验证给定的势函数是否满足 E=φB=×A 即可

2.1.2

u=λx,则 ui=λxi

ddλFj(λx)=Fjuiduidλ=Fjuixi=1λ(xiFjxi)=1λ(x)Fj(λx)

即:ddλF(λx)=1λ(x)F(λx)

  1. 证明 B=×A
×A=01λ×(x×B(t,λx))dλ

利用 ×(A×B)=(B)A+(B)A(A)B(A)B,注意此时 B=0(x)Bλ=λddλBλ

×(x×Bλ)=x0Bλ3+BλλddλBλ=2BλλddλBλ

角标 λ 没有其他含义,只是表示 λ 作为当前函数的自变量之一。上式代入积分:

×A=01(2λBλ+λ2ddλBλ)dλ=01ddλ(λ2B(t,λx))dλ=[λ2B(t,λx)]01=B(t,x)
  1. 证明 E=φtA

利用 (AB)=A×(×B)+B×(×A)+(A)B+(B)A 计算 φ

φ=01(xEλ)dλ=01[Eλ+(x)Eλ+x×(×Eλ)]dλ

利用法拉第定律 u×E=tB,得 ×Eλ=λtBλ

φ=01[Eλ+λddλEλλx×tBλ]dλ=01ddλ(λEλ)dλ+t01λx×Bλdλφ=E(t,x)tAE=φtA

证毕。

2.2

磁矢势满足 U(1) 规范:
AA=A+Λ

将变换后的 A 同样进行分解,由于 Λ 本身是一个梯度场(无旋),它完全属于纵向分量。根据亥姆霍兹分解的唯一性:

A+A=(A+Λ)+A

横向分量 A 在规范变换下保持不变(即它是规范无关的)。

B=×A=×(A+A)=0+×A=×A

可以看到,磁场完全由 A 决定。

2.3

参考 ppt

3 HW 4思考题

3.1 磁偶极子/带电旋转球

3.1.1

球壳上的面电流密度 K(x) 是由表面电荷随球壳转动产生的,其表达式为:

K(x)=σ0v=σ0(ω×x)

其中 x 为球壳上的源点位置矢量,且 |x|=a。所以磁矢势为

A(x)=μ04πSσ0(ω×x)|xx|dσ=μ0σ04πω×Sx|xx|dσ

令积分部分为 I=Sx|xx|dσ。设 x=ar^,面元 dσ=a2dΩ。球坐标系下(参见 3-D Laplace Equation And Spherical Harmonic Function#4.1 函数 $ frac{1}{ vec{r}- vec{a} }$ 的勒让德系数):

1|xx|=l=0r<lr>l+1Pl(cosγ),γ=arccosxx|xx|, r<=min(x,x), r>=max(x,x)

I=a3l=0r<lr>l+1SPl(cosγ)r^dΩ

由于 r^ 可以看作是由 l=1 的球谐函数线性组合而成,根据正交性,只有 l=1 的项对应的积分不为零。当 l=1 时,P1(cosγ)=cosγ=r^r^,故:

I=a3r<r>2S(r^r^)r^dΩ

利用课上的 Sr^r^dΩ=4π3I,得到:

I=a3r<r>24π3r^=4πa33r<r>2r^

所以

A(x)=μ0σ04πω×(4πa33r<r>2r^)=μ0σ0a33r<r>2(ω×r^)Ain(x)=μ0σ0a33ra2(ω×xr)=μ0σ0a3(ω×x),r<aAout(x)=μ0σ0a33ar2(ω×xr)=μ0σ0a43r3(ω×x),r>a

3.1.2

磁感应强度由 B=×A 给出。利用 ×(ω×x)=ω(x)(ω)x=3ωω=2ω,可得:

Bin=×[μ0σ0a3(ω×x)]=23μ0σ0aω

对于球外场,利用 ×(fV)=f×V+f(×V),设 f(r)=r3V=ω×x

(r3)×(ω×x)=3r4r^×(ω×x)=3r5x×(ω×x)

展开三重叉积 x×(ω×x)=ω(xx)x(xω)=r2ωr2(ωr^)r^,所以

×(ω×xr3)=3r5[r2ωr2(ωr^)r^]+r3(2ω)=1r3[3(ωr^)r^ω]

代入系数,得到球外的磁场分布:

Bout=μ0σ0a43r3[3(ωr^)r^ω]

显然这是一个磁偶极子场,其磁矩为 m=4π3σ0ωa4

3.1.3

和课上带电球受力的求解类似,我们选择包围整个上半空间(z>0)的闭合曲面 S。这个闭合曲面包括无穷远处的半球面(r,z>0)和赤道平面(z=0)。

由于在无穷远处,偶极子磁场 B1/r3,所以张量分量 T1/r6。而无穷远半球面的面积元 dar2,因此当 r 时,无穷远半球面上的积分项趋近于 0

所以,我们只需要计算在赤道平面(z=0)上的积分。对于上半空间,面元向量为 da=z^da

由于对称性,赤道面上的磁力只可能有 z 分量(径向受力相互抵消):

Fz=z=0(Tda)z=z=0Tzz(da)=z=0Tzzda

麦克斯韦应力张量的纯磁场部分定义为 Tij=1μ0(BiBj12δijB2)

Tzzin=12μ0(23μ0σ0aω)2=29μ0σ02a2ω2

所以

Fzin=0a(29μ0σ02a2ω2)2πrdr=4π9μ0σ02a2ω2[12r2]0a=2π9μ0σ02a4ω2

在赤道面上(ωr^=0),外部磁场退化为 Bout=μ0σ0a4ω3r3z^

Tzzout=12μ0(μ0σ0a4ω3r3)2=μ0σ02a8ω218r6

所以

Fzout=a(μ0σ02a8ω218r6)2πrdr=π9μ0σ02a8ω2ar5dr=π36μ0σ02a4ω2

最终得到

Fz=Fzin+Fzout=π4μ0σ02a4ω2

负号代表方向朝下(南半球方向),即两半球之间相互吸引。

3.2 电子经典球壳模型

3.2.1

对于半径为 a、电量为 e 的均匀带电球壳,球内电场为 0,球外电场为 E=e4πϵ0r2σ0=e4πa2

WE=a12ϵ0E2(4πr2)dr=e28πϵ0a

利用上一题结论,内部磁场 Bin=μ0eω6πa,外部为偶极场。

WM=12AKda

allB22μ0dV=12μ0B(×A)dV

利用分部积分:

B(×A)dV=A(×B)dV+(A×B)dS

对于局域分布的源,无穷远处的表面积分为 0。代入安培环路定理 ×B=μ0J

allB22μ0dV=12μ0A(μ0J)dV=12AJdV

已知球面上的矢量势为 A(a)=μ0σ0a3(ω×x)=μ0eωasinθ12πϕ^,面电流密度为 K=σ0v=σ0asinθω=eωsinθ4πaϕ^。所以磁能:

WM=120π(μ0eωasinθ12π)(eωsinθ4πa)(2πa2sinθ)dθ=μ0e2ω2a36π

总能量

W=e28πϵ0a+μ0e2ω2a36π

3.2.2

磁矩 m=4π3σ0ωa4=13eωa2,电磁场角动量

L=r×gdV=r×(ϵ0E×B)dV

由于球内 E=0,角动量只存在于球外,计算动量部分:

E×B=(e4πϵ0r2r^)×(μ04πr3[m])=μ0e16π2ϵ0r5(m×r^)

所以,被积函数的矢量部分是:

r^×(m×r^)=m(r^r^)r^(r^m)=mr^(r^m)

代入到求解角动量的积分得到

L=μ0e16π2a1r4r2drΩ[mr^(r^m)]dΩ

拆解该积分当中的每一项有:

ar2dr=1amdΩ=4πmr^(r^m)dΩ=(Im)dΩ=4π3m

所以最终得到旋转球壳的场角动量为:

L=μ0e16π21a8π3m=μ0e6πam=μ0e2a18πω

3.2.3

通常在经典模型中,由于 ωac 的预期,磁能项 WM 远小于静电能 WE。我们先忽略 WM

We28πϵ0a=mec2a=e28πϵ0mec2 (即经典电子半径);由角动量公式 μ0e2aω18π=2ω=9πμ0e2a。以及

v=ωa=9πμ0e2=9πϵ0c2e2

利用精细结构常数 α=e24πϵ0c1137,上式可改写为:

v=94αc9×1374c308c

这个速度远远超越了光速,是没有意义的,回过头来再用 WWM=μ0e2ω2a36π=mec2 也没什么意义。另外,据说将电子视为刚体,通过转动惯量来计算角动量,也会得到类似的结论(得到赤道角速度为一百多倍光速)。

4 HW 5思考题

4.1

Δs2=gαβdxαdxβ=gρσdxρdxσ=gρσΛραΛσβdxαdxβgαβ=gρσΛραΛσβΛγβΛγα=gγμgβνΛμνΛγα=gβνgνα=δβα

4.2 SO(3)

令矩阵为 J=(0110),注意到 J2=(1001)=I2×2J3=JJ4=I2×2,所以利用泰勒级数展开:

eθJ=IθJ+θ22!J2θ33!J3+=(cosθ)I(sinθ)J=(cosθsinθsinθcosθ)

所以

Λ(θ)=(1000010000cosθsinθ00sinθcosθ)

因为 [K,K]=0

Λ(θ1)Λ(θ2)=eθ1Keθ2K=e(θ1+θ2)K+12[θ1K,θ2K]+=Λ(θ1+θ2)

证明完毕。

4.3 SO(1,3)

kx=(0110),根据 Λx=eξΩx 可以得到其 Lie 群表示, 注意到 kx2=I2×2 :

Λx(ξ)=eξkx=n=0ξnkxnn!=n=0ξ2n(2n)!kx2n=1+n=0ξ2n+1(2n+1)!kx2n+1=kx=(n=0ξ2n(2n)!)I+(n=0ξ2n+1(2n+1)!)kx=(cosh(ξ)00cosh(ξ))+(0sinh(ξ)sinh(ξ)0)=(cosh(ξ)sinh(ξ)sinh(ξ)cosh(ξ))

即 :

Λx(ϕ)=(chξshξshξchξ11)

同样的逻辑,因为生成元矩阵 Ω 与其自身对易:

Λ(ξ1)Λ(ξ2)=eξ1Ωeξ2Ω=e(ξ1+ξ2)Ω=Λ(ξ1+ξ2)

这也反映了沿同一方向的两个洛伦兹提升,其快度是直接相加的。

4.4

对于 O(1,3) 群元 Λ1,Λ2

ΛTgΛ=g ΛTgΛ=(Λ1Λ2)Tg(Λ1Λ2)=Λ2T(Λ1TgΛ1)Λ2=Λ2TgΛ2=g

所以 Λ1Λ2O(3,1)

首先确定 Λ 是可逆的。对定义式 ΛTηΛ=η 两边取行列式:

det(ΛT)det(η)det(Λ)=det(η)

由于 det(η)0(其值为 1),且 det(ΛT)=det(Λ),可得:

(detΛ)2=1detΛ=±1

行列式不为 0,说明 Λ 必定存在逆矩阵 Λ1

在等式 ΛTηΛ=η 左边乘以 (ΛT)1,右边乘以 Λ1

LHS=(ΛT)1(ΛTgΛ)Λ1=g=RHS=(ΛT)1gΛ1=(Λ1)TgΛ1

因此 Λ1O(3,1)

4.5

4.5.1

成立

  1. 封闭性(乘法):
  2. 逆元存在性:
    • 如果 detΛ=1,则 det(Λ1)=1
    • 对于性质 Λ001Λ00=g0αg0βΛαβ=g00g00Λ00=Λ001
    • 因此,逆元依然属于该子集。

4.5.2

不成立。

除了包含单位矩阵的分支 SO(1,3) 以外,其余三个子集均不构成群。比如,假设我们取子集 Σi=PSO(1,3), (det=1,Λ001),选取 Λ1,Λ2Σi,则 detΛ1=detΛ2=1,它们的乘积 Λ1Λ2 的行列式为 (1)×(1)=1,落回了单位分支 SO(1,3) 中。

3.5.3 若 Λ1,Λ2 满足 (Λ1)001(Λ2)001,则它们的乘积 Λ=Λ1Λ2 的分量 Λ00 也必然满足 Λ001

对于任何 ΛO(3,1),取 ΛTgΛ=g 的00分量得:

(Λ00)2i=13(Λ0i)2=1(Λ00)2=1+i=13(Λ0i)21

所以,可以定义对于 SO(1,3) 群元,Λ001

a=( (Λ1)01,(Λ1)02,(Λ1)03 )b=( (Λ2)10,(Λ2)20,(Λ2)30 ) 其模长满足:

|a|2=i=13((Λ1)0i)2=((Λ1)00)21,|b|2=i=13((Λ2)i0)2=((Λ2)00)21

又,乘积矩阵 Λ=Λ1Λ2 的左上角元素为:

Λ00=(Λ1)0μ(Λ2)μ0=(Λ1)00(Λ2)00+i=13(Λ1)0i(Λ2)i0=(Λ1)00(Λ2)00+ab

根据柯西-施瓦茨不等式 |ab||a||b|,有:

(Λ1)00(Λ2)00Λ00|Λ00(Λ1)00(Λ2)00|=|ab||a||b|Λ00(Λ1)00(Λ2)00|a||b|

代入刚才的模长公式:

Λ00(Λ1)00(Λ2)00((Λ1)00)21((Λ2)00)21xyx21y21,x,y1

因为

(x21y21)2=x2y2y2x2+1x2y22xy+1=(xy1)2

所以总是有 xyx21y21+1,因此 Λ001

4.5.3

当一个粒子在它的瞬时共动参考系(Momentarily Comoving Reference Frame,MCRF)中感受到的固有加速度 a 为常数时,它在实验室坐标系下的运动轨迹并不是抛物线,而是双曲线。

首先,力 F 定义为动量的变化率:

F=dpdt=ddt(mγv)

在 MCRF 中,物体感受到的力满足由牛二定律 F=ma。所以有:

ma=ddt(mγv)a=ddt(γv)

其中 dvdt 是实验系测得的加速度。设 t=0 时,v=0,对运动方程进行积分:

v1v2/c2=atv(t)=at1+(at/c)2

世界线即轨迹 x(t)。对速度 v(t)=dxdt 再次积分。设 t=0 时,x=0

x(t)=0tat1+(at/c)2dt=c2a1+(atc)2c2a

改写上述方程,可以得到 xt 的双曲线关系:

(x+c2a)2(ct)2=(c2a)2