Skip to content

Feynman Diagrams for Beginners

Krešimir Kumerički

Department of Physics, Faculty of Science, University of Zagreb, Croatia

注释:这是一片有关费曼图计算的介绍,我为了入门QFT并且快速上手计算,对其进行了一次阅读批注,并且完成了其中的习题。所有文字由AI翻译成中文,并经过我的校对。

Abstract

本文对费曼图进行了简要介绍,并包含许多练习题。文本针对那些先前很少或没有接触过量子场论的学生。我们介绍了单粒子狄拉克方程、自由量子场以及使用费曼图构建费曼振幅的简明描述。作为一个例子,我们详细计算了电子和正电子湮灭成 μ 子对的截面。我们还展示了如何借助计算机进行此类计算。

Contents

1 自然单位

2 单粒子狄拉克方程

2.1 狄拉克方程

2.2 伴随狄拉克方程与狄拉克流

2.3 狄拉克方程的自由粒子解

3 自由量子场

3.1 自旋 0:标量场

3.2 自旋 1/2:狄拉克场

3.3 自旋 1:矢量场

4 衰变与散射的黄金规则

5 费曼图

6 示例:QED 中的 e+eμ+μ

6.1 对极化求和

6.2 Casimir 技巧

6.3 γ 矩阵的迹与缩并恒等式

6.4 质心系中的运动学

6.5 对二粒子相空间积分

6.6 步骤总结

6.7 Mandelstam 变量

附录:在计算机上计算费曼图

1 自然单位

为了描述某个物理系统或事件的运动学,我们可以自由选择三个基本运动学物理量的度量单位:长度 (L) 、质量 (M) 和时间 (T)。等价地,我们可以选择这些量的任意三个线性无关组合。通常选择 LTM(例如在 SI 单位制中),因为它们对于描述我们的直接经验最为方便。然而,基本粒子体验的是一个不同的世界,一个由相对论性量子力学定律支配的世界。

相对论性量子力学中的自然单位选择使得该理论的基本常数 c 都等于 1。[c]=LT1[]=ML2T1 ,为了完全确定我们的单位制,我们指定能量单位 (ML2T2)

1GeV=1.61010kgm2s2

约等于质子的质量。我们在实践中的做法是:

  • 在公式中忽略 c,只在最后(如果需要)恢复它们
  • GeVGeV1GeV2 等度量一切

示例:汤姆逊散射截面

在自然单位中,经典电磁辐射被自由电子散射(汤姆逊散射)的总截面为:

(1)σT=8πα23me2

为了恢复 c,我们将它们以一般幂次 αβ 插入上述方程,通过要求截面具有面积量纲 (L2) 来确定这些幂次:

(2)σT=8πα23me2αcβ[σ]=L2=1M2(ML2T1)α(LT1)βα=2,β=2

(3)σT=8πα23me22c2=0.6651024cm2=665mb

c 的线性无关性意味着这总是可以以唯一的方式完成。

以下转换关系通常很有用:

1fermi=5.07GeV11GeV2=0.389mb1GeV1=6.5821025s1kg=5.611026GeV1m=5.071015GeV11s=1.521024GeV1

练习 1 验证这些关系。

我们利用鲜于课上讲的关系对这些关系的正确性进行一个快速的口头估算验证。

1eV1=1μm1K=104eV1g=1023GeV

结合自然单位制的基本规律,1fermi=1015m=1023s,所以

1fermi=109μm109eV1=1GeV1

同理

1mb=1027cm2=1019μm21019eV2=10GeV2

还有

1GeV1=109eV1109μm=1015m1023s1kg=1000g1026GeV1m108s=106μm=106eV1=1015GeV1

GeV 计算要优雅得多。使用 me=0.511103GeV,我们立即得到:

(4)σT=8πα23me2=1709GeV2=665mb

练习 2 π0 粒子的衰变宽度为

(5)Γ=1τ=7.7eV

计算其寿命 τ(以秒为单位)。(顺便说一下,粒子的半衰期等于 τln2。)

同样进行估算。

τ0.1eV1=0.1μm=107m1015s

2 单粒子狄拉克方程

2.1 狄拉克方程

将相对论性能量方程

(6)E2=p2+m2

通过通常的替换

(7)pi,Eit

转化为微分方程,得到 Klein-Gordan 方程:

(8)(+m2)ψ(x)=0

将其解释为单粒子波方程时,存在有问题的负能量解。这是由于 E=±p2+m2 中的负根。也就是说,在相对论力学中可以忽略这个负根,但在量子物理中必须保留微分方程的所有完全解。

为了克服这个问题,狄拉克尝试了如下 ansatz

(9)(iβμμ+m)(iγννm)ψ(x)=0

其中 βμγν 将通过要求与 Klein-Gordan 方程的一致性来确定。这要求 γμ=βμ

(10)γμμγνν=μμ

这又意味着

(γ0)2=1,(γi)2=1{γμ,γν}γμγν+γνγμ=0forμν

这可以紧凑地写成反对易关系的形式

(11){γμ,γν}=2gμν,gμν=(1000010000100001)

这些条件显然不可能用通常的数字 γ 来满足,但我们可以通过取 γ 为(至少)4×4 矩阵来满足它们。现在,要满足 (9),只需该方程中的两个因子之一为零,按照惯例我们要求第二个因子为零。由此我们得到狄拉克方程:

(12)(iγμμm)ψ(x)=0

ψ(x) 现在有四个分量,称为狄拉克旋量。最常用的 γ 矩阵表示之一是原始的狄拉克表示

(13)γ0=(1001)γi=(0σiσi0)

其中 σi 是泡利矩阵。这种表示对于非相对论近似非常方便,因为此时主导的能量项 (iγ00m)ψ(0) 是对角化的。另外两种常用的表示是:

  • Weyl(或手征)表示——在超相对论区域(Em)方便
  • Majorana 表示——使狄拉克方程变为实数;对于反粒子等于粒子的 Majorana 费米子方便

问题:为什么我们最多只能选择一个 γ 矩阵为对角矩阵?

这是显然的,否则反对易括号就会出现非对角项非零的情况。

泡利矩阵的性质

(15)σi=σi(16)σi=(iσ2)σi(iσ2)(17)[σi,σj]=2iϵijkσk(18){σi,σj}=2δij(19)σiσj=δij+iϵijkσk

练习 3 证明对于任意三维矢量 a(σa)2=a2

σiaiσjajδij=(δij+iϵijkσk)δijaiaj=aiai

练习 4 利用泡利矩阵的性质,证明狄拉克表示中的 γ 矩阵满足 {γi,γj}=2gij=δij,符合反对易关系。(反对易关系的其他分量,(γ0)2=1{γ0,γi}=0,是平凡的证明。)

{γi,γj}=γiγj+γjγi=(0σiσi0)(0σjσj0)+(0σjσj0)(0σiσi0)=2({σi,σj}00{σi,σj})=2δij

练习 5 证明在狄拉克表示中 γ0γμγ0=γμ

γ0γμγ0=(1001)γμ(1001)

假设 γμ=(abcd),此时

(1001)γμ(1001)=(1001)(abcd)(1001)=(abcd)

带入 γi=(0σiσi0),所以

γ0γμγ0=γμ

练习 6 通过将狄拉克方程写成 iψ/t=Hψ 的形式来确定狄拉克哈密顿量。证明狄拉克哈密顿量的厄米性意味着无论表示如何,前一个练习中的关系都成立。Feynman slash 记号 /aaμγμ 经常被使用。

改写狄拉克方程

iγμμψ=mψ

iψt=(m+iγjjγ0)ψ

所以哈密顿量为

H=m+iγjγ0j

2.2 伴随狄拉克方程与狄拉克流

为了构建狄拉克流,我们需要 ψ(x) 的方程。取狄拉克方程的dagger,我们得到

ψγ0(i/+m)=0

我们定义伴随旋量 ψ¯ψγ0,得到伴随狄拉克方程

ψ¯(x)(i/+m)=0

引入 ψ¯ 不仅是为了得到美学上令人愉悦的方程,而且因为可以证明,与 ψ 不同,它在洛伦兹变换下协变地变换。

练习 7 验证流 jμ=ψ¯γμψ 是守恒的,即满足连续性关系 μjμ=0。这个相对论性四流的分量是 jμ=(ρ,j)。注意 ρ= j0=ψ¯γ0ψ=ψψ>0,即概率是正定的。

μjμ=μ(ψ¯γμψ)=μ(ψγ0γμψ)=(μψ)γ0γμψ+ψγ0γμ(μψ)=mψ¯ψ+mψ¯ψ=0

2.3 狄拉克方程的自由粒子解

由于我们正在为微扰理论计算做准备,我们只需要考虑自由粒子解。关于各种势中的解,请参阅文献。狄拉克旋量满足克莱因-戈尔登方程这一事实提示了以下 ansatz

(20)ψ(x)=u(p)eipx

将其代入狄拉克方程后,得到动量空间狄拉克方程

(21)(pm)u(p)=0.

这有两个正能量解

(22)u(p,σ)=N(χ(σ)σpE+mχ(σ)),σ=1,2,

其中

(23)χ(1)=(10),χ(2)=(01),

以及两个负能量解,它们被解释为正能量的反粒子解

(24)v(p,σ)=N(σpE+m(iσ2)χ(σ)(iσ2)χ(σ)),σ=1,2,E>0.

N 是待定的归一化常数。上面的旋量与 [1] 中的一致。反粒子解的动量空间狄拉克方程是

(25)(+m)v(p,σ)=0.

可以证明,两个解(一个 σ=1,另一个 σ=2)对应于自旋-1/2 粒子的两个自旋状态。

练习 8 确定 u¯(pp,σ)v¯(p,σ) 的动量空间狄拉克方程。

归一化

在非相对论单粒子量子力学中,波函数的归一化很简单。粒子在空间中某处的概率等于一,这转化为归一化条件 ψψdV=1。另一方面,我们最终将在多粒子量子场论中使用旋量 (22) 和 (24),因此它们的归一化不是唯一的。我们将选择归一化约定:单位体积内有 2E 个粒子

(26)unitvolumeρdV=unitvolumeψψdV=2E

这个选择是相对论协变的,因为体积元的洛伦兹收缩被能量变化所补偿。还有其他具有其他优点的归一化约定。

练习 9 根据此选择确定归一化常数 N

完备性

练习 10 使用显式表达式 (22) 和 (24) 证明

(27)σ=1,2u(p,σ)u¯(p,σ)=+m,(28)σ=1,2v(p,σ)v¯(p,σ)=m.

这些关系在计算未极化费米子的费曼图时经常需要。参见后面的章节。

宇称与双线性协变量

宇称变换:

  • P:xx,tt
  • P: ψ → γ 0 ψ

练习 11 验证流 jμ=ψ¯γμψ 在宇称下像矢量一样变换,即 j0j0jj

任何费米子流都具有 ψ¯Γψ 的形式,其中 Γ 是某个 4×4 矩阵。为了构建相互作用拉格朗日量,我们希望只使用那些具有确定洛伦兹变换性质的流。为此,我们首先定义两个新矩阵:

(29)γ5iγ0γ1γ2γ3=Diracrep.(0110),{γ5,γμ}=0(30)σμνi2[γμ,γν],σμν=σνμ

现在,如果 Γ 是下表中给出的矩阵之一,则 ψ¯Γψ 将协变地变换。下表还显示了 ψ¯Γψ 的变换性质、Γ 中不同 γ 矩阵的数量以及 Γ 的分量数。

Γ变换性质# of γ's# of components
1标量01
γμ矢量14
σμν张量26
γ5γμ轴矢量34
γ5赝标量41

这穷尽了所有可能性。分量总数为 16,这意味着集合 {1,γμ,σμν,γ5γμ,γ5} 构成任何 4×4 矩阵的完备基。这样的 ψ¯Γψ 流称为双线性协变量。


3 自由量子场

单粒子狄拉克方程(a)即使对于单粒子系统(如 H 原子)也不完全正确,并且(b)无法处理多粒子过程(如 β -衰变 npeν¯)。我们必须升级到量子场论。

任何狄拉克场都是以下完备集的某种叠加

u(p,σ)eipx,v(p,σ)eipx,σ=1,2,pR3

我们可以将其写为

(31)ψ(x)=σd3p(2π)32E[u(p,σ)a(p,σ)eipx+v(p,σ)ac(p,σ)eipx].

这里 1/(2π)32E 是归一化因子(有许多不同的约定),a(pp,σ)ac(pp,σ) 是展开系数。为了使这成为量子狄拉克场,我们通过强加反对易关系将这些系数提升为算符:

(32){a(p,σ),a(p,σ)}=δσσδ3(pp),

类似地对于 ac(pp,σ)。(对于玻色场,我们将有对易关系。)这类似于通过 [xi,pj]=iδij 对易关系将位置和动量提升为算符,这就是为什么从单粒子量子理论过渡到量子场论有时被称为二次量子化。

算符 a 作用于真空态 |0 时,产生单粒子态 |p,σ

(33)a(p,σ)|0=|p,σ,

这就是它被称为产生算符的原因。类似地,a 是湮灭算符:

(34)a(p,σ)|p,σ=|0,

acac 是反粒子态的产生和湮灭算符(ac 中的 c 代表“共轭”)。

粒子物理学中的过程大多在此类场的理论框架内计算——量子场论。该理论可以在不同严格程度上描述,但无论如何都足够复杂,超出了这些笔记的范围。

然而,量子场论关于基本粒子相互作用的预测通常可以使用相对简单的“配方”——费曼图来计算。

在我们转向描述费曼图方法之前,让我们先指定参与基本粒子物理相互作用的其他量子场。所有这些场都是自由场,相互作用被视为它们的微扰。每种粒子类型(电子、光子、希格斯玻色子……)都有自己的量子场。

3.1 自旋 0:标量场

例如,希格斯玻色子、π 介子……

(35)ϕ(x)=d3p(2π)32E[a(p)eipx+ac(p)eipx]

3.2 自旋 1/2:狄拉克场

例如,夸克、轻子

我们已经指定了狄拉克自旋-1/2 场。还有其他类型:Weyl 和 Majorana 自旋-1/2 场,但它们超出了我们的范围。

3.3 自旋 1:矢量场

可能是

• 有质量的(例如 W、Z 弱玻色子)或
• 无质量的(例如光子)

(36)Aμ(x)=λd3p(2π)32E[ϵμ(p,λ)a(p,λ)eipx+ϵμ(p,λ)a(p,λ)eipx]

ϵμ(pp,λ) 是极化矢量。对于有质量粒子,它自动满足

(37)pμϵμ(p,λ)=0

而在无质量情况下,由于规范不变性(洛伦兹规范条件),可以强加这个条件。这意味着有质量矢量粒子只有三个独立极化:λ=1,2,3λ= +,,0。在无质量情况下,可以进一步利用规范对称性消除一个极化态,只留下两个:λ=1,2λ=+,

极化矢量的归一化使得

(38)ϵ(p,λ)ϵ(p,λ)=1.

例如,对于沿 z 轴运动的有质量粒子 (p=(E,0,0,|pp|)),我们可以取

(39)ϵ(p,±)=12(01±i0),ϵ(p,0)=1m(|p|00E)

练习 12 计算

λϵμ(p,λ)ϵν(p,λ)

提示:将其写成最一般的形式 (Agμν+Bpμpν),然后确定 AB

显然,所得结果不能简单地通过极限 m0 外推到无质量情况。规范对称性使无质量极化求和稍微复杂一些,但为了进行简单的费曼图计算,可以只使用以下关系

λϵμ(p,λ)ϵν(p,λ)=gμν.

4 衰变与散射的黄金规则

粒子物理学的主要实验可观测量是

• 散射截面 σ(1+21+2++n)
• 衰变宽度 Γ(11+2++n)

另一方面,理论是用量子场的拉格朗日密度定义的,例如

L=12μϕμϕ12m2ϕ2g4!ϕ4.

如何从 L 计算 σ 和 Γ?

在非相对论量子理论中,为了计算存在相互作用势 VI 时从态 |α 到态 |β 的跃迁速率,我们有费米黄金规则

(40)(αβt r a n s i t i o n r a t e)=2π|β|VIα|2×(d e n s i t y o f f i n a lq u a n t u m s t a t e s).

这是最低阶微扰理论。对于更高阶,我们有更多相互作用势矩阵元 |VI| 的乘积项。

在量子场论中,有与这些矩阵元对应的东西——S 矩阵:

(41)β|VI|α+(h i g h e r - o r d e r t e r m s)β|S|α.

一方面,S 矩阵元可以(在知道相互作用拉格朗日量/哈密顿量的情况下)借助戴森级数进行微扰计算

(42)S=1id4x1H(x1)+(i)22!d4x1d4x2T{H(x1)H(x2)}+,

另一方面,我们有将这些矩阵元与截面和衰变宽度联系起来的“黄金规则”。

将这些黄金规则用费曼不变振幅 M 表示很方便,M 是通过从 S - 矩阵中剥离一些运动学因子得到的:

(43)β|S|α=δβαi(2π)4δ4(pβpα)Mβαi=α,β1(2π)32Ei.

现在我们有两个规则:

11+2++n 的部分衰变率

(44)dΓ=12E1|Mβα|2(2π)4δ4(p1p1pn)i=1nd3pi(2π)32Ei,

1+21+2++n 散射的微分截面

(45)dσ=1uα12E112E2|Mβα|2(2π)4δ4(p1+p2p1pn)i=1nd3pi(2π)32Ei,

其中 uα 是粒子 1 和 2 的相对速度:

(46)uα=(p1p2)2m12m22E1E2,

|M|2 是费曼不变振幅对未测量的粒子自旋的平均(参见第 6.1 节)。M 的量纲(以能量为单位)为

• 对于衰变 [M]=3n
• 对于两个粒子的散射 [M]=2n

其中 n 是产生的粒子数。

因此,计算某个可观测量包括两个阶段:

  1. 确定 |M|2。为此,我们使用下一节将介绍的费曼图方法。
  2. 对洛伦兹不变相空间积分
dL i p s=(2π)4δ4(p1+p2p1pn)i=1nd3pi(2π)32Ei.

5 费曼图

示例:ϕ4 -理论

L=12μϕμϕ12m2ϕ2g4!ϕ4

• 自由(动能)拉格朗日量(恰好包含两个场的项)决定了理论的粒子及其传播子。这里我们只有一个标量场: [图片:标量传播子]

• 相互作用拉格朗日量(包含三个或更多场的项)决定了可能的顶点。这里,同样只有一个顶点: [图片:四点顶点]

我们构建所有具有固定外粒子的可能图。例如,对于此理论中两个标量粒子的散射,我们会有

费曼图时间从左向右流动。有些人绘制费曼图时时间向上流动,这更符合我们通常在相对论物理中绘制时空的方式。

由于每个顶点对应于 (42) 中的一个相互作用拉格朗日量(哈密顿量)项,带有圈的图对应于微扰理论的更高阶。这里我们只处理最低阶,因此只使用树图。

为了实际写出费曼振幅 M,我们有一套费曼规则,将因子与费曼图的元素关联起来。具体来说,为了得到 iM,我们按以下方式构建费曼规则:

• 顶点因子就是 i 乘以(动量空间)拉格朗日量中的相互作用项,并移除所有场:

iLI=ig4!ϕ4r e m o v i n g f i e l d sig4![顶点图示]

• 传播子是 i 乘以动量空间中动能算符(由自由运动方程定义)的逆:

(48)Lf r e eE u l e r - L a r g r a n e e q .(μμ+m2)ϕ=0(K l e i n - G o r d o n e q .)

使用替换 μipμ 转到动量空间,然后取逆得到:

(49)(p2m2)ϕ=0ϕ=ip2m2

(实际上,正确的费曼传播子是 i/(p2m2+iϵ),但为了我们的目的,我们可以忽略无穷小 iϵ 项。)

• 外线由适当的极化矢量或旋量表示(在 (31)、(35)、(36) 及其共轭场中,位于相应产生或湮灭算符旁边的那个):

粒子费曼规则
入射费米子u
出射费米子u
入射反费米子v
出射反费米子v
入射光子εμ
出射光子εμ*
入射标量1
出射标量1

因此,这个 ϕ4 理论中标量-标量散射振幅的树图贡献就是

(50)iM=ig4!.

练习 13 确定电子传播子以及量子电动力学(QED)唯一顶点的费曼规则:

(51)L=ψ¯(i∂̸+eAm)ψ14FμνFμνFμν=μAννAμ.

注意也有

(52)p=iσu(p,σ)u¯(p,σ)p2m2,

即,电子传播子只是标量传播子乘以极化求和。很好的是,这可以推广到所有粒子的传播子。这非常有帮助,因为由于规范对称性的复杂性,反转动子动能算符并非易事。因此,矢量粒子的传播子是:

有质量情况: [传播子表达式] 无质量情况:

(54)p=igμνp2.

原则上,这几乎是我们需要知道的全部,以便能够计算任何给定过程的费曼振幅。注意,传播子和外线极化矢量仅由粒子类型(其自旋和质量)决定,因此上述相应规则不仅限于 ϕ4 理论和 QED,而且适用于所有标量、自旋-1 矢量玻色子和狄拉克费米子(如标准模型)的理论。我们需要的唯一附加信息是顶点因子。

上一段中的“几乎”指的是,在一般的费曼图计算中,还有几个额外的微妙之处:

• 在圈图中,一些内部动量未确定,我们必须对这些动量积分。此外,每个封闭的费米子圈有一个额外的因子 (-1)。由于我们只考虑树图,我们可以忽略这一点。
• 当相同场进入单个顶点时,存在一些组合数值因子。
• 有时图之间存在相对 (-) 号。
• 如果终态有相同粒子,则存在对称性因子。

关于这些的解释,建议读者查阅一些量子场论教科书。


6 示例:QED 中的 e+eμ+μ

只有一个贡献的树图:

[图片:s-channel 电子-正电子湮灭成μ子对的费曼图]

我们使用 QED 的费曼规则并沿着费米子线反向写下振幅。线本身的顺序无关紧要。

(55)iM=[u¯(p3,σ3)(ieγν)v(p4,σ4)]igμν(p1+p2)2[v¯(p2,σ2)(ieγμ)u(p1,σ1)],

或者,引入缩写 u1u(pp1,σ1)

(56)M=e2(p1+p2)2[u¯3γμv4][v¯2γμu1].

练习 14 为康普顿散射 γeγe 绘制费曼图并写下振幅。

6.1 对极化求和

如果我们知道所有外粒子的动量和极化,我们就可以显式计算 M。然而,实验通常使用未极化粒子进行,因此我们必须对未态粒子的极化(自旋)求和,并对初态粒子的极化(自旋)求平均:

(57)|M|2|M|2=1212σ1σ2a v g . o v e r i n i t i a l p o l .σ3σ4s u m o v e r f i n a l p o l .|M|2.

因子 1/2 是因为每个初态费米子有两个极化(自旋)态。

(问题:为什么我们对概率求和而不是对振幅求和?)

在计算 |M|2=MM 时,需要以下恒等式

(58)[u¯γμv]=[uγ0γμv]=vγμγ0u=[v¯γμu].

因此,

(59)|M|2=e44(p1+p2)4σ1,2,3,4[v¯4γμu3][u¯1γμv2][u¯3γνv4][v¯2γνu1].

6.2 Casimir 技巧

使用以下技巧可以轻松执行对极化的求和。首先,我们用显式的旋量指标 α,β,γ,δ=1,2,3,4 写出 [u¯1γμv2][v¯2γνu1]

(60)σ1σ2u¯1αγαβμv2βv¯2γγγδνu1δ.

我们现在可以将 u1δ 移到前面(u1δ 只是一个数字,是 u1 向量的元素,因此它与所有东西都对易),然后使用完备性关系 (27) 和 (28),

σ1u1δu¯1α=(p1+m1)δα,σ2v2βv¯2γ=(p2m2)βγ,

这将求和 (60) 变为

(61)(1+m1)δαγαβμ(2m2)βγγγδν=Tr[(1+m1)γμ(2m2)γν].

这意味着

(62)|M|2=e44(p1+p2)4Tr[(p1+m1)γμ(p2m2)γν]Tr[(p4m4)γμ(p3+m3)γν].

因此我们摆脱了所有旋量,只剩下 γ 矩阵的迹。这些可以使用下一节的关系来评估。

6.3 γ 矩阵的迹与缩并恒等式

所有这些都是反对易关系 {γμ,γν}=2gμν{γμ,γ5}= 0,(γ5)2=1 的结果,而不是其他任何东西!

迹恒等式

  1. 奇数个 γ 的迹为零:
Tr(γμ1γμ2γμ2n+1)=Tr(γμ1γμ2γμ2n+1γ5γ51)(m o v i n gγ5o v e r e a c hγμi)=Tr(γ5γμ1γμ2γμ2n+1γ5)(c y c l i c p r o p e r t y o f t r a c e)=Tr(γμ1γμ2γμ2n+1γ5γ5)=Tr(γμ1γμ2γμ2n+1)=0
  1. Tr 1 = 4

Trγμγν=Tr(2gμνγνγμ)=(2.)8gμνTrγνγμ=8gμνTrγμγν2Trγμγν=8gμνTrγμγν=4gμν

这也意味着:

Tr//=4ab
  1. 练习 15 计算 Tr(γμγνγργσ)。提示:使用反对易关系将 γσ 一直移到左边。然后使用 3。

作业:证明 Tr(γμ1γμ2γμ2n)(2n1)!! 项。

  1. Tr(γ5γμ1γμ2γμ2n+1)=0。这可以从 1. 以及 γ5 由偶数个 γ 组成的事实得出。
  2. Trγ5=Tr(γ0γ0γ5)=Tr(γ0γ5γ0)=Trγ5=0
  3. Tr(γ5γμγν)=0。(与上面相同的技巧,用 γαμ,ν 代替 γ0。)
  4. Tr(γ5γμγνγργσ)=4iϵμνρσ,其中 ϵ0123=1。注意:也有使用 ϵ0123=1 的约定。

缩并恒等式

γμγμ=12gμν(γμγν+γνγμ)2gμν=gμνgμν=4
γμγαγμγμγα+2gμα=4γα+2γα=2γα
  1. 练习 16 缩并 γμγαγβγμ
  2. γμγαγβγγγμ=2γγγβγα

练习 17 计算 |M|2 中的迹:

Tr[(1+m1)γμ(2m2)γν]=?Tr[(p4m4)γμ(p3+m3)γν]=?

练习 18 计算 ||M|2

6.4 质心系中的运动学

e+e 对撞机中,通常 pime,mμi=1,,4,因此我们可以取

mi0‘高能”或“极端相对论”极限

那么

(63)|M|2=8e4(p1+p2)4[(p1p3)(p2p4)+(p1p4)(p2p3)]

为了计算散射截面 σ,我们必须具体到某个特定参考系(σ 不是与参考系无关的)。对于 e+e 对撞机,最相关的是质心 (CM) 系:

[图片:质心系运动学示意图]

练习 19Eθ 表示 ||M|2

6.5 对二粒子相空间积分

现在我们可以对 1+23+4 微分散射截面使用“黄金规则” (45)

(64)dσ=1uα12E112E2|M|2dLips2

其中要积分的二粒子相空间为

(65)dLips2=(2π)4δ4(p1+p2p3p4)d3p3(2π)32E3d3p4(2π)32E4.

首先,我们对六个积分变量中的四个积分,并在一般参考系中进行。δ -函数使得对 d3p4 的积分变得平凡:

(66)dLips2=1(2π)24E3E4δ(E1+E2E3E4)d3p3p32d|p3|dΩ3

现在,我们注意到 E3E4|p3| 的函数:

E3=E3(|p3|)=p32+m32,E4=p42+m42=p32+m42,

并通过 δ -函数关系

(67)δ(E1+E2p32+m32p32+m42)=δ[f(|p3|)]=δ(|p3||p3(0)|)|f(|p3|)||p3|=|p3(0)|.

这里 |p3| 只是积分变量,|p3(0)|f(|p3|) 的零点,即第三个粒子的实际动量。在对 d|p3| 积分后,我们令 |p3(0)||p3|

由于

(68)f(|p3|)=E3+E4E3E4|p3|,

我们得到

(69)dLips2=|p3|dΩ16π2(E1+E2).

现在我们再次具体到 CM 系,并注意通量因子为

(70)4E1E2uα=4(p1p2)2m12m22=4|p1|(E1+E2),

最终给出

(71)dσCMdΩ=164π2(E1+E2)2|p3||p1||M|2.

注意,我们在每一步都保留了质量,因此这个公式对于任何 CM 散射都是普遍有效的。

对于我们特定的 ee+μμ+ 散射,这给出了微分截面的最终结果(引入精细结构常数 α=e2/(4π))

(72)dσdΩ=α216E2(1+cos2θ).

练习 20 对此积分得到总截面 σ

注意,显然 σα2,并且量纲分析要求 σ1/E2,因此只有角度依赖关系 (1+cos2θ) 检验了 QED 作为轻子和光子理论的地位。

6.6 步骤总结

概括地说,计算(未极化)散射截面(或衰变宽度)包括以下步骤:

  1. 绘制费曼图
  2. 使用费曼规则写出 iM
  3. 平方 M 并使用 Casimir 技巧得到迹
  4. 评估迹
  5. 应用所选参考系的运动学
  6. 对相空间积分

6.7 Mandelstam 变量

Mandelstam 变量 s、tu 常用于散射计算。它们(对于 1+23+4 散射)定义为

s=(p1+p2)2t=(p1p3)2u=(p1p4)2

练习 21 证明 s+t+u=m12+m22+m32+m42

这意味着只有两个 Mandelstam 变量是独立的。它们的主要优点是它们是洛伦兹不变量,这使得它们在费曼振幅计算中很方便。只有在最后,我们才能将它们换成“实验者”的变量 Eθ

练习 22 用 Mandelstam 变量表示 ee+μμ+ 散射的 |M|2


附录:在计算机上计算费曼图

有几个计算机程序可以执行费曼图计算中的部分或全部步骤。以下是与其中一个程序(FeynCalc [2],Wolfram Mathematica 的包)的一个简单会话,我们在其中计算了与正文中相同的过程,ee+μμ+。另一个仅依赖开源软件的框架是 FORM [3]。

FeynCalc 演示

这个 Mathematica 笔记本演示了使用 Feyncalc 包对 e e+μμ+ 散射进行费曼不变振幅的计算机计算。

首先,我们将 FeynCalc 加载到 Mathematica 中:

I n [ 1 ] := <<HighEnergyPhysics`fc`

FeynCalc 4.1.0.3b 评估 ?FeynCalc 获取帮助或访问 www.feyncalc.org

在使用费曼规则并应用 Casimir 技巧后,自旋平均的费曼振幅平方 |M|2

I n [ 2 ] := Msq = (e^4/(4 (p1 + p2)^4)) * Contract[ 
  Tr[(GS[p1] + me).GA[μ].(GS[p2] - me).GA[ν]] *
  Tr[(GS[p4] - mm).GA[μ].(GS[p3] + mm).GA[ν]] ]
O u t [ 2 ] = (1/(4 (p1 + p2)^4)) (e^4 (64 mm^2 me^2 + 32 p3·p4 me^2 + 32 mm^2 p1·p2 + 32 p1·p4 p2·p3 + 32 p1·p3 p2·p4))

迹被自动评估,缩并自动执行。现在我们通过替换规则引入 Mandelstam 变量:

I n [ 3 ] := prod[a_, b_] := Pair[Momentum[a], Momentum[b]];
mandelstam = {prod[p1, p2] -> (s - me^2 - me^2)/2,
              prod[p3, p4] -> (s - mm^2 - mm^2)/2,
              prod[p1, p3] -> (t - me^2 - mm^2)/2,
              prod[p2, p4] -> (t - me^2 - mm^2)/2,
              prod[p1, p4] -> (u - me^2 - mm^2)/2,
              prod[p2, p3] -> (u - me^2 - mm^2)/2,
              (p1 + p2) -> Sqrt[s]};

并将这些替换应用于我们的振幅:

I n [ 5 ] := Msq /. mandelstam
O u t [ 5 ] = (1/(4 s^2)) (e^4 (64 mm^2 me^2 + 16 (s - 2 mm^2) me^2 + 8 (- me^2 - mm^2 + t)^2 + 8 (- me^2 - mm^2 + u)^2 + 16 mm^2 (s - 2 me^2)))

通过消除一个 Mandelstam 变量可以简化这个结果:

I n [ 6 ] := Simplify[TrickMandelstam[%, s, t, u, 2 me^2 + 2 mm^2]]
O u t [ 6 ] = (2 e^4 (2 me^4 + 4 (mm^2 - u) me^2 + 2 mm^4 + s^2 + 2 u^2 - 4 mm^2 u + 2 s u))/s^2

如果我们进入超相对论极限,我们得到的结果与手动计算一致:

I n [ 7 ] := Simplify[% /. {mm -> 0, me -> 0}]
O u t [ 7 ] = (2 e^4 (t^2 + u^2))/s^2

参考文献

[1] D. Griffiths, Introduction to elementary particles, Wiley (1987)

[2] V. Shtabovenko, R. Mertig and F. Orellana, New Developments in FeynCalc 9.0, arXiv:1601.01167 [hep-ph].
[3] J. A. M. Vermaseren, New features of FORM, math-ph/0010025.

Built with VitePress