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 中的
6.1 对极化求和
6.2 Casimir 技巧
6.3
6.4 质心系中的运动学
6.5 对二粒子相空间积分
6.6 步骤总结
6.7 Mandelstam 变量
附录:在计算机上计算费曼图
1 自然单位
为了描述某个物理系统或事件的运动学,我们可以自由选择三个基本运动学物理量的度量单位:长度
相对论性量子力学中的自然单位选择使得该理论的基本常数
约等于质子的质量。我们在实践中的做法是:
- 在公式中忽略
和 ,只在最后(如果需要)恢复它们 - 用
、 、 等度量一切
示例:汤姆逊散射截面
在自然单位中,经典电磁辐射被自由电子散射(汤姆逊散射)的总截面为:
为了恢复
即
以下转换关系通常很有用:
练习 1 验证这些关系。
我们利用鲜于课上讲的关系对这些关系的正确性进行一个快速的口头估算验证。
结合自然单位制的基本规律,
,所以 同理
还有
用
练习 2
计算其寿命
同样进行估算。
2 单粒子狄拉克方程
2.1 狄拉克方程
将相对论性能量方程
通过通常的替换
转化为微分方程,得到 Klein-Gordan 方程:
将其解释为单粒子波方程时,存在有问题的负能量解。这是由于
为了克服这个问题,狄拉克尝试了如下 ansatz
其中
这又意味着
这可以紧凑地写成反对易关系的形式
这些条件显然不可能用通常的数字
其中
- Weyl(或手征)表示——在超相对论区域(
)方便 - Majorana 表示——使狄拉克方程变为实数;对于反粒子等于粒子的 Majorana 费米子方便
问题:为什么我们最多只能选择一个
这是显然的,否则反对易括号就会出现非对角项非零的情况。
泡利矩阵的性质
练习 3 证明对于任意三维矢量
练习 4 利用泡利矩阵的性质,证明狄拉克表示中的
练习 5 证明在狄拉克表示中
假设
,此时 带入
,所以
练习 6 通过将狄拉克方程写成
改写狄拉克方程
即
所以哈密顿量为
2.2 伴随狄拉克方程与狄拉克流
为了构建狄拉克流,我们需要
我们定义伴随旋量
引入
练习 7 验证流
2.3 狄拉克方程的自由粒子解
由于我们正在为微扰理论计算做准备,我们只需要考虑自由粒子解。关于各种势中的解,请参阅文献。狄拉克旋量满足克莱因-戈尔登方程这一事实提示了以下 ansatz
将其代入狄拉克方程后,得到动量空间狄拉克方程
这有两个正能量解
其中
以及两个负能量解,它们被解释为正能量的反粒子解
可以证明,两个解(一个
练习 8 确定
归一化
在非相对论单粒子量子力学中,波函数的归一化很简单。粒子在空间中某处的概率等于一,这转化为归一化条件
这个选择是相对论协变的,因为体积元的洛伦兹收缩被能量变化所补偿。还有其他具有其他优点的归一化约定。
练习 9 根据此选择确定归一化常数
完备性
练习 10 使用显式表达式 (22) 和 (24) 证明
这些关系在计算未极化费米子的费曼图时经常需要。参见后面的章节。
宇称与双线性协变量
宇称变换:
- P:
- P: ψ → γ 0 ψ
练习 11 验证流
任何费米子流都具有
现在,如果
| Γ | 变换性质 | # of γ's | # of components |
|---|---|---|---|
| 1 | 标量 | 0 | 1 |
| γμ | 矢量 | 1 | 4 |
| σμν | 张量 | 2 | 6 |
| γ5γμ | 轴矢量 | 3 | 4 |
| γ5 | 赝标量 | 4 | 1 |
这穷尽了所有可能性。分量总数为 16,这意味着集合
3 自由量子场
单粒子狄拉克方程(a)即使对于单粒子系统(如 H 原子)也不完全正确,并且(b)无法处理多粒子过程(如
任何狄拉克场都是以下完备集的某种叠加
我们可以将其写为
这里
类似地对于
算符
这就是它被称为产生算符的原因。类似地,
而
粒子物理学中的过程大多在此类场的理论框架内计算——量子场论。该理论可以在不同严格程度上描述,但无论如何都足够复杂,超出了这些笔记的范围。
然而,量子场论关于基本粒子相互作用的预测通常可以使用相对简单的“配方”——费曼图来计算。
在我们转向描述费曼图方法之前,让我们先指定参与基本粒子物理相互作用的其他量子场。所有这些场都是自由场,相互作用被视为它们的微扰。每种粒子类型(电子、光子、希格斯玻色子……)都有自己的量子场。
3.1 自旋 0:标量场
例如,希格斯玻色子、π 介子……
3.2 自旋 1/2:狄拉克场
例如,夸克、轻子
我们已经指定了狄拉克自旋-1/2 场。还有其他类型:Weyl 和 Majorana 自旋-1/2 场,但它们超出了我们的范围。
3.3 自旋 1:矢量场
可能是
• 有质量的(例如 W、Z 弱玻色子)或
• 无质量的(例如光子)
而在无质量情况下,由于规范不变性(洛伦兹规范条件),可以强加这个条件。这意味着有质量矢量粒子只有三个独立极化:
极化矢量的归一化使得
例如,对于沿
练习 12 计算
提示:将其写成最一般的形式
显然,所得结果不能简单地通过极限
4 衰变与散射的黄金规则
粒子物理学的主要实验可观测量是
• 散射截面
• 衰变宽度
另一方面,理论是用量子场的拉格朗日密度定义的,例如
如何从
在非相对论量子理论中,为了计算存在相互作用势
这是最低阶微扰理论。对于更高阶,我们有更多相互作用势矩阵元
在量子场论中,有与这些矩阵元对应的东西——S 矩阵:
一方面,S 矩阵元可以(在知道相互作用拉格朗日量/哈密顿量的情况下)借助戴森级数进行微扰计算
另一方面,我们有将这些矩阵元与截面和衰变宽度联系起来的“黄金规则”。
将这些黄金规则用费曼不变振幅
现在我们有两个规则:
•
•
其中
而
• 对于衰变
• 对于两个粒子的散射
其中
因此,计算某个可观测量包括两个阶段:
- 确定
。为此,我们使用下一节将介绍的费曼图方法。 - 对洛伦兹不变相空间积分
5 费曼图
示例:
• 自由(动能)拉格朗日量(恰好包含两个场的项)决定了理论的粒子及其传播子。这里我们只有一个标量场: [图片:标量传播子]
• 相互作用拉格朗日量(包含三个或更多场的项)决定了可能的顶点。这里,同样只有一个顶点: [图片:四点顶点]
我们构建所有具有固定外粒子的可能图。例如,对于此理论中两个标量粒子的散射,我们会有
费曼图时间从左向右流动。有些人绘制费曼图时时间向上流动,这更符合我们通常在相对论物理中绘制时空的方式。
由于每个顶点对应于 (42) 中的一个相互作用拉格朗日量(哈密顿量)项,带有圈的图对应于微扰理论的更高阶。这里我们只处理最低阶,因此只使用树图。
为了实际写出费曼振幅
• 顶点因子就是
• 传播子是
使用替换
(实际上,正确的费曼传播子是
• 外线由适当的极化矢量或旋量表示(在 (31)、(35)、(36) 及其共轭场中,位于相应产生或湮灭算符旁边的那个):
| 粒子 | 费曼规则 |
|---|---|
| 入射费米子 | u |
| 出射费米子 | u |
| 入射反费米子 | v |
| 出射反费米子 | v |
| 入射光子 | εμ |
| 出射光子 | εμ* |
| 入射标量 | 1 |
| 出射标量 | 1 |
因此,这个
练习 13 确定电子传播子以及量子电动力学(QED)唯一顶点的费曼规则:
注意也有
即,电子传播子只是标量传播子乘以极化求和。很好的是,这可以推广到所有粒子的传播子。这非常有帮助,因为由于规范对称性的复杂性,反转动子动能算符并非易事。因此,矢量粒子的传播子是:
有质量情况: [传播子表达式] 无质量情况:
原则上,这几乎是我们需要知道的全部,以便能够计算任何给定过程的费曼振幅。注意,传播子和外线极化矢量仅由粒子类型(其自旋和质量)决定,因此上述相应规则不仅限于
上一段中的“几乎”指的是,在一般的费曼图计算中,还有几个额外的微妙之处:
• 在圈图中,一些内部动量未确定,我们必须对这些动量积分。此外,每个封闭的费米子圈有一个额外的因子 (-1)。由于我们只考虑树图,我们可以忽略这一点。
• 当相同场进入单个顶点时,存在一些组合数值因子。
• 有时图之间存在相对 (-) 号。
• 如果终态有相同粒子,则存在对称性因子。
关于这些的解释,建议读者查阅一些量子场论教科书。
6 示例:QED 中的
只有一个贡献的树图:
[图片:s-channel 电子-正电子湮灭成μ子对的费曼图]
我们使用 QED 的费曼规则并沿着费米子线反向写下振幅。线本身的顺序无关紧要。
或者,引入缩写
练习 14 为康普顿散射
6.1 对极化求和
如果我们知道所有外粒子的动量和极化,我们就可以显式计算
因子
(问题:为什么我们对概率求和而不是对振幅求和?)
在计算
因此,
6.2 Casimir 技巧
使用以下技巧可以轻松执行对极化的求和。首先,我们用显式的旋量指标
我们现在可以将
这将求和 (60) 变为
这意味着
因此我们摆脱了所有旋量,只剩下
6.3 矩阵的迹与缩并恒等式
所有这些都是反对易关系
迹恒等式
- 奇数个
的迹为零:
Tr 1 = 4
这也意味着:
- 练习 15 计算
。提示:使用反对易关系将 一直移到左边。然后使用 3。
作业:证明
。这可以从 1. 以及 由偶数个 组成的事实得出。 。(与上面相同的技巧,用 代替 。) ,其中 。注意:也有使用 的约定。
缩并恒等式
- 练习 16 缩并
练习 17 计算
练习 18 计算
6.4 质心系中的运动学
在
那么
为了计算散射截面
[图片:质心系运动学示意图]
练习 19 用
6.5 对二粒子相空间积分
现在我们可以对
其中要积分的二粒子相空间为
首先,我们对六个积分变量中的四个积分,并在一般参考系中进行。
现在,我们注意到
并通过
这里
由于
我们得到
现在我们再次具体到 CM 系,并注意通量因子为
最终给出
注意,我们在每一步都保留了质量,因此这个公式对于任何 CM 散射都是普遍有效的。
对于我们特定的
练习 20 对此积分得到总截面
注意,显然
6.6 步骤总结
概括地说,计算(未极化)散射截面(或衰变宽度)包括以下步骤:
- 绘制费曼图
- 使用费曼规则写出
- 平方
并使用 Casimir 技巧得到迹 - 评估迹
- 应用所选参考系的运动学
- 对相空间积分
6.7 Mandelstam 变量
Mandelstam 变量 s、
练习 21 证明
这意味着只有两个 Mandelstam 变量是独立的。它们的主要优点是它们是洛伦兹不变量,这使得它们在费曼振幅计算中很方便。只有在最后,我们才能将它们换成“实验者”的变量
练习 22 用 Mandelstam 变量表示
附录:在计算机上计算费曼图
有几个计算机程序可以执行费曼图计算中的部分或全部步骤。以下是与其中一个程序(FeynCalc [2],Wolfram Mathematica 的包)的一个简单会话,我们在其中计算了与正文中相同的过程,
FeynCalc 演示
这个 Mathematica 笔记本演示了使用 Feyncalc 包对
首先,我们将 FeynCalc 加载到 Mathematica 中:
I n [ 1 ] := <<HighEnergyPhysics`fc`FeynCalc 4.1.0.3b 评估 ?FeynCalc 获取帮助或访问 www.feyncalc.org
在使用费曼规则并应用 Casimir 技巧后,自旋平均的费曼振幅平方
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 /. mandelstamO 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.
