# QED の 1 ループ：真空偏極・電子の異常磁気モーメント・Ward–高橋恒等式

> 次元正則化とファインマン・パラメータの基本積分公式を示し、QED の 1 ループで真空偏極・走る結合定数とβ関数・Landau 極を導く。頂点補正から Schwinger の a=α/2π を計算し、Ward–高橋恒等式から Z₁=Z₂ を示して、赤外発散の Bloch–Nordsieck 相殺を説明する。
> https://rikai.mugen-giken.com/physics/qft/qed-one-loop

## 0. この記事の要点

電子や光子の振る舞いを精密に測ると、単純な図 1 本の計算では説明できない小さなずれが見つかります。この記事は、そのずれを「仮想的な粒子が輪になって現れる効果」（ループ補正）として計算する章です。電子の磁石としての強さが予想よりわずかに大きいこと、電荷の見え方が観測するエネルギーで変わること、そしてそれらを支えるゲージ不変性の関係式を、具体的な計算で確かめます。読み終えると、1 ループの積分をどう処理するかと、その結果が何を意味するかを説明できるようになります。

- 1 ループ積分は、**ファインマン・パラメータ**で分母を 1 つにまとめ、ループ運動量を平行移動し、Wick 回転して Euclid 空間の球対称な積分に直すと、すべて 1 つの公式 $\int\frac{d^d\ell_E}{(2\pi)^d}\frac{1}{(\ell_E^2+\Delta)^n}=\frac{1}{(4\pi)^{d/2}}\frac{\Gamma(n-d/2)}{\Gamma(n)}\Delta^{d/2-n}$ に帰着します。**次元正則化**はこの右辺を $d$ について解析接続したもので、ゲージ不変性（平行移動不変性）を保ちます。
- **Ward–高橋恒等式**（ゲージ不変性から出る、頂点補正と自己エネルギーの関係式）$q_\mu\Gamma^\mu(p+q,p)=\not{q}-\Sigma(p+q)+\Sigma(p)$ を 1 ループで証明し、そこから頂点と電子場のくりこみ定数が等しいこと $Z_1=Z_2$ を導きます。その結果、電荷のくりこみは光子の $Z_3$ だけで決まり、$e=Z_3^{1/2}e_0$ はどの荷電粒子にも共通です。
- **真空偏極**は横波 $\Pi^{\mu\nu}=(q^2g^{\mu\nu}-q^\mu q^\nu)\Pi(q^2)$ で、$\Pi_2(q^2)=-\frac{2\alpha}{\pi}\int_0^1dx\,x(1-x)\big[\frac2\varepsilon-\gamma_E+\ln4\pi-\ln\Delta\big]$ です。有効結合は $-q^2\gg m^2$ で $\alpha_{\mathrm{eff}}=\alpha/\big(1-\frac{\alpha}{3\pi}\ln\frac{-q^2}{Am^2}\big)$、$A=e^{5/3}$ となり、β関数は $\beta(e)=\frac{e^3}{12\pi^2}$ です。結合は高エネルギーで増え、形式的には **Landau 極**（結合定数が有限のエネルギーで無限大になる点）に達します。
- 頂点補正の形状因子 $F_2$ は有限で、$F_2(0)=\frac{\alpha}{2\pi}$（**Schwinger の結果**）。したがって電子の $g$ 因子は $g=2+\frac{\alpha}{\pi}+O(\alpha^2)$ です。
- 頂点補正の $F_1$ は光子質量 $\mu\to0$ で**赤外発散**（エネルギーの極めて低い光子に由来する発散）しますが、検出器で分離できない軟光子の放射を足し合わせると $\ln\mu$ が相殺します（**Bloch–Nordsieck**）。光子を 1 個も出さない排他的な断面積は、全次数で見ると $0$ です。

<Aside type="tip" title="この記事で使う事前知識">

- **ファインマン則と LSZ 簡約公式**: 場の量子論の図を積分の式に読み替える規則と、外線（入射・出射粒子）の扱い方。→ [摂動論とファインマン則](/physics/qft/feynman-rules-lsz)
- **Dirac 方程式**: スピン $\frac12$ の電子が従う相対論的な波動方程式。→ [Dirac 方程式](/physics/quantum-mechanics-2/relativistic-quantum-mechanics#def-dirac-equation)
- **γ 行列**: 反交換関係 $\{\gamma^\mu,\gamma^\nu\}=2g^{\mu\nu}$ を満たす行列で、Dirac 方程式を書くのに使います。→ [γ 行列](/physics/quantum-mechanics-2/relativistic-quantum-mechanics#def-gamma-matrices)
- **走る結合定数**: 観測するエネルギーの大きさによって値が変わる結合定数。→ [走る結合定数](/physics/qft/dimensional-regularization-rg#def-running-coupling)
- **1PI 図とくりこみ**: 1 本の線を切っても二つに分かれない図（1 粒子既約図）と、発散を裸のパラメータに吸収する考え方。

</Aside>

## 1. 動機：ツリーの先にあるもの

[前章](/physics/qft/feynman-rules-lsz)でファインマン則を導き、$e^+e^-\to\mu^+\mu^-$ の断面積をツリーレベルで計算しました。ツリーレベルの答えは古典電磁気学の延長で理解できる部分が大きく、量子論らしさが本格的に現れるのはループからです。QED の 1 ループには、歴史的に場の量子論の信頼性を決定づけた 3 つの結果が含まれています。

1. **真空偏極**：仮想的な電子・陽電子対が電荷を遮蔽し、結合定数がエネルギーとともに変化します。$\alpha^{-1}$ が低エネルギーの $137.036$ から $Z$ ボソンの質量スケールで約 $129$ に下がることは実測されています（[くりこみ理論入門](/physics/qft/renormalization)の <Ref to="physics/qft/renormalization#ex-running-couplings" />）。
2. **異常磁気モーメント**：[Dirac 方程式](/physics/quantum-mechanics-2/relativistic-quantum-mechanics#def-dirac-equation)は $g=2$ を予言しますが、1947 年の Kusch–Foley の測定は $g/2\approx1.00119$ でした。Schwinger は 1948 年にこのずれが $\alpha/2\pi$ であることを示しました。
3. **Ward 恒等式**：ループには紫外発散が現れますが、ゲージ不変性によって発散のうち独立なものが減り、電荷のくりこみが粒子の種類によらないことが保証されます。

一方で、ループ計算は 2 種類の発散に出会います。高運動量で生じる**紫外発散**はくりこみで吸収し、無質量の光子に由来する**赤外発散**は観測量の定義を見直すことで消えます。本章はこの 2 つの処理を、具体的な計算として最後まで追いかけます。符号の規約は前章と同じく Peskin–Schroeder に従い、$D_\mu=\partial_\mu+ieA_\mu$、電子の電荷を $e=-|e|$、$\alpha=e^2/4\pi$ とします（<Ref to="physics/qft/feynman-rules-lsz#rem-yukawa-qed-rules" />）。

<Figure caption="QED の 1 ループの 1PI 図。左から電子の自己エネルギー Σ、真空偏極 Π、頂点補正 δΓ。破線は内線の光子、波線は外線の光子を表す。">
<svg viewBox="0 0 600 170" width="100%" role="img" aria-label="電子の自己エネルギー、真空偏極、頂点補正の三つのファインマン図">
  <g stroke="currentColor" stroke-width="2" fill="none">
    <line x1="10" y1="110" x2="180" y2="110" />
    <path d="M 45 110 C 55 40, 135 40, 145 110" stroke="var(--sl-color-accent)" stroke-dasharray="6 4" />
    <path d="M 210 90 q 8 -10 16 0 q 8 10 16 0 q 8 -10 16 0" stroke="var(--sl-color-accent)" />
    <ellipse cx="300" cy="90" rx="42" ry="26" />
    <path d="M 342 90 q 8 -10 16 0 q 8 10 16 0 q 8 -10 16 0" stroke="var(--sl-color-accent)" />
    <line x1="440" y1="160" x2="510" y2="100" />
    <line x1="580" y1="160" x2="510" y2="100" />
    <path d="M 510 100 q -10 -8 0 -16 q 10 -8 0 -16 q -10 -8 0 -16 q 10 -8 0 -16" stroke="var(--sl-color-accent)" />
    <path d="M 463 140 L 557 140" stroke="var(--sl-color-accent)" stroke-dasharray="6 4" />
  </g>
  <g fill="currentColor" font-size="14">
    <text x="80" y="140">Σ(p)</text>
    <text x="285" y="140">Π(q)</text>
    <text x="490" y="30">γ(q)</text>
    <text x="430" y="150">p</text>
    <text x="580" y="150">p′</text>
  </g>
</svg>
</Figure>

## 2. 1 ループ積分の技法

### 2.1. ファインマン・パラメータ

ループ積分の被積分関数は、伝播関数（分母）の積です。このままでは積分しにくいので、まず分母を 1 つの式の冪にまとめます。次の補題はそのための道具です。

<Lemma id="lem-feynman-parameters" title="ファインマン・パラメータの公式">
$n\ge1$ とし、複素数 $A_1,\dots,A_n$ がすべて同じ開半平面 $\{A:\ \mathrm{Re}(e^{i\theta}A)\gt0\}$（$\theta$ は実数）に含まれるとします。このとき
$$
\frac{1}{A_1A_2\cdots A_n}=\Gamma(n)\int_0^1dx_1\cdots\int_0^1dx_n\;\frac{\delta\big(\sum_ix_i-1\big)}{\big(x_1A_1+\cdots+x_nA_n\big)^n}
$$
が成り立ちます。特に $\frac{1}{AB}=\int_0^1\frac{dx}{[xA+(1-x)B]^2}$、$\frac{1}{ABC}=2\int_0^1dx\,dy\,dz\,\frac{\delta(x+y+z-1)}{[xA+yB+zC]^3}$ です。
</Lemma>

<Proof of="lem-feynman-parameters">
両辺とも各 $A_i$ について $-1$ 次の同次式なので、すべての $A_i$ に $e^{i\theta}$ を掛けても等式は変わりません。よって $\mathrm{Re}\,A_i\gt0$ の場合に示せば十分です。このとき $\frac{1}{A_i}=\int_0^\infty dt_i\,e^{-t_iA_i}$ なので
$$
\prod_i\frac1{A_i}=\int_{[0,\infty)^n}d^nt\;e^{-\sum_it_iA_i}
$$
です。$t_i=s\,x_i$（$s\gt0$、$\sum_ix_i=1$）と変数変換すると、$d^nt=s^{n-1}ds\,d^nx\,\delta(\sum x_i-1)$ です（$s=\sum t_i$ を新しい変数に選んだ射影座標のヤコビアン）。$s$ 積分は $\int_0^\infty ds\,s^{n-1}e^{-sX}=\Gamma(n)X^{-n}$（$X=\sum x_iA_i$、$\mathrm{Re}\,X\gt0$）なので主張を得ます。
</Proof>

Minkowski 空間の伝播関数の分母 $k^2-m^2+i\varepsilon$ は虚部が正なので、$\theta=-\pi/2$ としてこの補題を使えます。

### 2.2. Wick 回転と次元正則化

ファインマン・パラメータでまとめたあと、ループ運動量を $\ell=k+(\text{外線運動量の 1 次結合})$ と平行移動すると、分母は $(\ell^2-\Delta+i\varepsilon)^n$ の形になります。$\Delta\gt0$ のとき、$\ell^0$ の極 $\pm(\sqrt{\boldsymbol\ell^2+\Delta}-i\varepsilon)$ は第 2・第 4 象限にあるので、$\ell^0$ の積分路を反時計回りに $90^\circ$ 回して虚軸に重ねられます。$\ell^0=i\ell_E^0$、$\boldsymbol\ell=\boldsymbol\ell_E$ と置くと $\ell^2=-\ell_E^2$、$d^4\ell=i\,d^4\ell_E$ です。これを **Wick 回転**と呼びます。

次の定義は、4 次元では発散する積分に、次元 $d$ を連続的に動かすことで意味を与える方法です。直感的には、積分が収束する小さい $d$ で答えを式の形で求めておき、その式を $d=4$ の近くまで延ばして読みます。式がガンマ関数で書けるため、$d$ が整数でなくても意味をもちます。

<Definition id="def-dimensional-regularization" title="次元正則化">
自然数 $d$ と $\Delta\gt0$ について、$d$ 次元 Euclid 空間の積分
$$
I_n(d,\Delta)=\int\frac{d^d\ell_E}{(2\pi)^d}\frac{1}{(\ell_E^2+\Delta)^n}
$$
は $d\lt2n$ で収束します。その値（<Ref to="thm-master-integrals" />）は $d$ の有理型関数に解析接続でき、極は $d=2n,2n+2,\dots$ にあります。**次元正則化**とは、4 次元で発散するループ積分を、Wick 回転・平行移動のあとこの解析接続で定義し、$d=4-\varepsilon$ として $\varepsilon\to0$ の極を分離する処方です。結合定数の次元を保つため、$d$ 次元では $e$ を $e\,\mu^{\varepsilon/2}$ に置き換えます（$\mu$ は任意の質量スケール）。
</Definition>

<Lemma id="lem-dimreg-shift" title="次元正則化の平行移動不変性">
$f(k)$ を、伝播関数の分母 $(k+p_i)^2-m_i^2+i\varepsilon$ の積の逆数に $k$ の多項式を掛けた 1 ループの被積分関数とし、<Ref to="def-dimensional-regularization" /> で $\int\frac{d^dk}{(2\pi)^d}f(k)$ を定める。任意の外線運動量 $q$ について
$$
\int\frac{d^dk}{(2\pi)^d}f(k+q)=\int\frac{d^dk}{(2\pi)^d}f(k).
$$
</Lemma>

<Proof of="lem-dimreg-shift">
<Ref to="def-dimensional-regularization" /> の手順では、分母を <Ref to="lem-feynman-parameters" /> でまとめると $\sum_ix_i=1$ から $\sum_ix_i\bigl((k+p_i)^2-m_i^2\bigr)=(k+P)^2-\Delta$（$P=\sum_ix_ip_i$、$\Delta=P^2-\sum_ix_ip_i^2+\sum_ix_im_i^2$ は $x_i$、$p_i$、$m_i$ だけで決まる）となり、$\ell=k+P$ と平行移動して、分子を $\ell$ の多項式として書いたうえで $\ell$ について Wick 回転と解析接続を行います。$f(k+q)$ では $p_i$ が $p_i+q$ に変わるだけなので、$P$ は $P+q$ に、$\Delta$ は $(P+q)^2-\sum_ix_i(p_i+q)^2+\sum_ix_im_i^2=P^2-\sum_ix_ip_i^2+\sum_ix_im_i^2$ で変わらず、$\ell=k+q+P$ と置くと分子の $\ell$ の多項式も $f(k)$ の場合と同じです。したがって定義に現れる $\ell$ の積分はまったく同じ式になり、値は一致します。
</Proof>

<Remark id="rem-dimreg-axioms" title="次元正則化の性質">
こうして定義した「$d$ 次元積分」は、被積分関数について線形で、**積分変数の平行移動で不変**で、スケール変換について $d$ 次元の積分と同じ性質をもちます。逆にこの 3 性質と規格化（Gauss 積分の値）から一意に定まることが知られています（Collins『Renormalization』第 4 章）。平行移動不変性は第 3 節の Ward–高橋恒等式の証明で本質的に使われます。運動量のカットオフはこの性質をもたず、ゲージ不変性を壊します（<Ref to="ex-cutoff-breaks-ward" />）。$\gamma$ 行列は $\{\gamma^\mu,\gamma^\nu\}=2g^{\mu\nu}$（$g^\mu{}_\mu=d$）と $\mathrm{tr}\,1=4$ をそのまま使います。$\gamma^5$ は 4 次元に固有なので、本章の計算には現れないようにしています。
</Remark>

<Theorem id="thm-master-integrals" title="基本積分公式">
$\Delta\gt0$、$n$ を自然数とし、$d$ を $d\lt2n$ から解析接続すると
$$
\int\frac{d^d\ell_E}{(2\pi)^d}\frac{1}{(\ell_E^2+\Delta)^n}=\frac{1}{(4\pi)^{d/2}}\frac{\Gamma(n-\frac d2)}{\Gamma(n)}\Delta^{\frac d2-n},\qquad \int\frac{d^d\ell_E}{(2\pi)^d}\frac{\ell_E^2}{(\ell_E^2+\Delta)^n}=\frac{1}{(4\pi)^{d/2}}\frac d2\frac{\Gamma(n-\frac d2-1)}{\Gamma(n)}\Delta^{\frac d2-n+1}
$$
が成り立ちます。Minkowski 空間では、Wick 回転により
$$
\int\frac{d^d\ell}{(2\pi)^d}\frac{1}{(\ell^2-\Delta)^n}=\frac{(-1)^ni}{(4\pi)^{d/2}}\frac{\Gamma(n-\frac d2)}{\Gamma(n)}\Delta^{\frac d2-n},\qquad \int\frac{d^d\ell}{(2\pi)^d}\frac{\ell^2}{(\ell^2-\Delta)^n}=\frac{(-1)^{n-1}i}{(4\pi)^{d/2}}\frac d2\frac{\Gamma(n-\frac d2-1)}{\Gamma(n)}\Delta^{\frac d2-n+1}
$$
です。
</Theorem>

<Proof of="thm-master-integrals">
$\int_{\mathbb R^d}e^{-\ell^2}d^d\ell=\pi^{d/2}$ を極座標で計算すると $|S^{d-1}|\int_0^\infty r^{d-1}e^{-r^2}dr=|S^{d-1}|\Gamma(d/2)/2$ なので、単位球面の面積は $|S^{d-1}|=2\pi^{d/2}/\Gamma(d/2)$ です。動径積分は $t=r^2/\Delta$、$u=1/(1+t)$ と置くと Beta 関数になり
$$
\int_0^\infty\frac{r^{d-1}dr}{(r^2+\Delta)^n}=\frac{\Delta^{\frac d2-n}}{2}\int_0^\infty\frac{t^{\frac d2-1}dt}{(1+t)^n}=\frac{\Delta^{\frac d2-n}}2\frac{\Gamma(\frac d2)\Gamma(n-\frac d2)}{\Gamma(n)}
$$
です。$|S^{d-1}|/(2\pi)^d$ を掛けると $\Gamma(d/2)$ が約分され、$\pi^{d/2}/(2\pi)^d=(4\pi)^{-d/2}$ から第 1 式を得ます。右辺は $d$ の有理型関数なので、これが <Ref to="def-dimensional-regularization" /> の解析接続そのものです。第 2 式は $\frac{\ell_E^2}{(\ell_E^2+\Delta)^n}=\frac{1}{(\ell_E^2+\Delta)^{n-1}}-\frac{\Delta}{(\ell_E^2+\Delta)^n}$ に第 1 式を適用し、$\frac{\Gamma(n-1-\frac d2)}{\Gamma(n-1)}-\frac{\Gamma(n-\frac d2)}{\Gamma(n)}=\frac{\Gamma(n-1-\frac d2)}{\Gamma(n)}\big[(n-1)-(n-1-\frac d2)\big]$ と整理すれば得られます。Minkowski 版は $d^d\ell=i\,d^d\ell_E$、$\ell^2-\Delta=-(\ell_E^2+\Delta)$ を代入したものです。
</Proof>

$d=4-\varepsilon$ では $\Gamma(\varepsilon/2)=\frac2\varepsilon-\gamma_E+O(\varepsilon)$、$\Delta^{-\varepsilon/2}=1-\frac\varepsilon2\ln\Delta+\cdots$ なので、$n=2$ の積分は
$$
\frac{\Gamma(2-\frac d2)}{(4\pi)^{d/2}}\Delta^{\frac d2-2}=\frac{1}{16\pi^2}\Big[\frac2\varepsilon-\gamma_E+\ln4\pi-\ln\Delta\Big]+O(\varepsilon)
$$
となります（[くりこみ理論入門](/physics/qft/renormalization)の Appendix と同じ規約です）。

<Lemma id="lem-symmetric-integration" title="対称積分と d 次元の γ 行列">
$f$ を $\ell^2$ だけの関数とすると、次元正則化された積分の中で $\ell^\mu\mapsto0$、$\ell^\mu\ell^\nu\mapsto\frac1dg^{\mu\nu}\ell^2$ と置き換えてよいです。また $d$ 次元で $\gamma^\nu\gamma_\nu=d$、$\gamma^\nu\gamma^\mu\gamma_\nu=-(d-2)\gamma^\mu$、$\gamma^\nu\not{a}\gamma^\mu\not{b}\gamma_\nu=-2\not{b}\gamma^\mu\not{a}+(4-d)\not{a}\gamma^\mu\not{b}$ が成り立ちます。
</Lemma>

<Proof of="lem-symmetric-integration">
$\ell\mapsto-\ell$ の不変性から 1 次の積分は $0$ です。$\int\ell^\mu\ell^\nu f$ は Lorentz 共変な 2 階テンソルで外部ベクトルを含まないので $Cg^{\mu\nu}$ の形であり、$g_{\mu\nu}$ で縮約すると $C=\frac1d\int\ell^2f$ です。$\gamma$ 行列の等式は反交換関係から順に出ます。$\gamma^\nu\gamma^\mu\gamma_\nu=(2g^{\nu\mu}-\gamma^\mu\gamma^\nu)\gamma_\nu=2\gamma^\mu-d\gamma^\mu$。3 つ目は $\gamma^\nu\not{a}\gamma^\mu\not{b}\gamma_\nu=\gamma^\nu\not{a}\gamma^\mu(2b_\nu-\gamma_\nu\not{b})=2\not{b}\not{a}\gamma^\mu-\gamma^\nu\not{a}\gamma^\mu\gamma_\nu\not{b}$ と書き、$\gamma^\nu\not{a}\gamma^\mu\gamma_\nu=4a^\mu-(4-d)\not{a}\gamma^\mu$（同じ方法で示せます）と $\not{a}\gamma^\mu+\gamma^\mu\not{a}=2a^\mu$ を使うと得られます。
</Proof>

## 3. くりこみ定数と Ward–高橋恒等式

この節の目的は、電荷のくりこみが電子の種類や頂点補正の詳細によらないことを示すことです。そのために、まず 3 種類の補正（電子の自己エネルギー、光子の真空偏極、頂点補正）にくりこみ定数 $Z_1,Z_2,Z_3$ を対応させ、次にゲージ不変性から出る関係式を使って $Z_1=Z_2$ を導きます。

<Definition id="def-qed-renormalization-constants" title="1PI 関数とくりこみ定数">
裸の結合 $e_0$、裸の質量 $m_0$ で書いた QED の摂動論で、外線を切り取った 1 粒子既約（1PI）図の和を次のように表します。

- 電子の自己エネルギー $\Sigma(p)$：1PI な電子 2 点図の和を $-i\Sigma(p)$ とします。完全な電子の伝播関数は $\frac{i}{\not{p}-m_0-\Sigma(p)}$ です。
- 真空偏極 $\Pi^{\mu\nu}(q)$：1PI な光子 2 点図の和を $i\Pi^{\mu\nu}(q)$ とします。
- 頂点関数 $\Gamma^\mu(p',p)$：入射電子 $p$、出射電子 $p'$、光子 $q=p'-p$ の 1PI 3 点図の和を $-ie_0\Gamma^\mu(p',p)$ とします。最低次は $\Gamma^\mu=\gamma^\mu$ です。

物理的質量 $m$ は $m-m_0-\Sigma(\not{p}=m)=0$ の解です。くりこみ定数を
$$
Z_2^{-1}=1-\frac{d\Sigma}{d\not{p}}\Big|_{\not{p}=m},\qquad Z_3^{-1}=1-\Pi(0),\qquad \bar u(p)\Gamma^\mu(p,p)u(p)=Z_1^{-1}\bar u(p)\gamma^\mu u(p)\ \ (p^2=m^2)
$$
で定めます（$\Pi(q^2)$ は <Ref to="prop-vacuum-polarization-transverse" /> の横波成分）。$Z_2,Z_3$ は完全な伝播関数の極の留数であり、<Ref to="physics/qft/feynman-rules-lsz#thm-lsz" /> の $Z$ にあたります。**物理的電荷**を $e=e_0Z_1^{-1}Z_2Z_3^{1/2}$ で定めます。$Z_1,Z_2$ は赤外発散を含むので、必要なときは光子に小さな質量 $\mu$ を与えて正則化します。
</Definition>

物理的電荷の定義は、運動量移行 $q\to0$ の散乱振幅が $e^2/q^2$ になるように選んだものです。2 つの荷電粒子が光子 1 個を交換する振幅は、各頂点の $e_0Z_1^{-1}$、LSZ 公式による外線の $\sqrt{Z_2}$ 4 個、光子の伝播関数の留数 $Z_3$ を掛けた $e_0^2Z_1^{-2}Z_2^2Z_3/q^2$ を $q\to0$ で与えるからです。1 ループでは $e_0$ と $e$ の差は高次の効果なので、図の中の結合を単に $e$ と書きます。

次の定理は、頂点補正に光子の運動量 $q$ を掛けて縮約すると、自己エネルギーの差になるという関係です。証明の方針は、$\not{q}$ を 2 つの電子の伝播関数の逆数の差として書き、頂点の図を伝播関数の差に崩すことです。

<Theorem id="thm-ward-takahashi" title="Ward–高橋恒等式（1 ループ）">
次元正則化した QED の $e^2$ の次数で、光子の伝播関数を任意の共変ゲージにとるとき、任意の（質量殻外の）$p,q$ について
$$
q_\mu\Gamma^\mu(p+q,p)=\not{q}-\Sigma(p+q)+\Sigma(p),\qquad q_\mu\Pi^{\mu\nu}(q)=0
$$
が成り立ちます。ここで $\Gamma^\mu=\gamma^\mu+\delta\Gamma^\mu$、$\Sigma$、$\Pi^{\mu\nu}$ はそれぞれ $e^2$ の次数の 1 ループ図です。
</Theorem>

<Proof of="thm-ward-takahashi">
自由電子の伝播関数 $S(k)=\frac{i}{\not{k}-m}$ について、$\not{q}=(\not{k}+\not{q}-m)-(\not{k}-m)$ と書くと
$$
S(k+q)\,(-ie\not{q})\,S(k)=-ie\cdot i\cdot i\Big[\frac{1}{\not{k}-m}-\frac{1}{\not{k}+\not{q}-m}\Big]=e\big[S(k)-S(k+q)\big]
$$
です。同様に $S(k)(-ie\not{q})S(k+q)=e[S(k)-S(k+q)]$ です。光子の運動量 $l$ を積分変数にとり、光子の伝播関数を $D_{\nu\rho}(l)$ と書くと、<Ref to="physics/qft/feynman-rules-lsz#thm-yukawa-qed-rules" /> より
$$
-i\Sigma(p)=\int\frac{d^dl}{(2\pi)^d}(-ie\gamma^\nu)S(p-l)(-ie\gamma^\rho)D_{\nu\rho}(l),\qquad -ie\,\delta\Gamma^\mu(p',p)=\int\frac{d^dl}{(2\pi)^d}(-ie\gamma^\nu)S(p'-l)(-ie\gamma^\mu)S(p-l)(-ie\gamma^\rho)D_{\nu\rho}(l)
$$
です。後者に $q_\mu$（$p'=p+q$）を掛け、上の等式を $k=p-l$ で使うと
$$
-ie\,q_\mu\delta\Gamma^\mu=e\int\frac{d^dl}{(2\pi)^d}(-ie\gamma^\nu)\big[S(p-l)-S(p+q-l)\big](-ie\gamma^\rho)D_{\nu\rho}(l)=e\big[-i\Sigma(p)+i\Sigma(p+q)\big]
$$
となり、$q_\mu\delta\Gamma^\mu=\Sigma(p)-\Sigma(p+q)$ です。$q_\mu\gamma^\mu=\not{q}$ と合わせて第 1 式を得ます。この計算は $D_{\nu\rho}$ の形を使わないので、どの共変ゲージでも成り立ちます。

真空偏極はフェルミオンの閉じたループなので、因子 $-1$ とトレースが付き
$$
i\Pi^{\mu\nu}(q)=-\int\frac{d^dk}{(2\pi)^d}\mathrm{tr}\big[(-ie\gamma^\mu)S(k)(-ie\gamma^\nu)S(k+q)\big]
$$
です。$q_\nu$ を掛けて同じ等式を使うと $q_\nu i\Pi^{\mu\nu}=ie^2\int\mathrm{tr}[\gamma^\mu(S(k)-S(k+q))]$ となり、第 2 項で $k\mapsto k-q$ と平行移動すると第 1 項と打ち消し合います。$\Pi^{\mu\nu}(q)=\Pi^{\nu\mu}(-q)$ なので $q_\mu\Pi^{\mu\nu}=0$ も同じです。最後の平行移動は <Ref to="lem-dimreg-shift" /> によって正当化されます。
</Proof>

<Remark id="rem-ward-takahashi-all-orders" title="全次数の Ward–高橋恒等式">
全次数では、完全な伝播関数 $S_{\mathrm{full}}(p)=\frac{i}{\not{p}-m_0-\Sigma(p)}$ と完全な頂点について $q_\mu\Gamma^\mu(p+q,p)=\not{q}-\Sigma(p+q)+\Sigma(p)$ が成り立ちます（Ward 1950、高橋 1957。Peskin–Schroeder 第 7 章 §7.4。証明は <Ref to="physics/qft/ward-takahashi-schwinger-dyson#cor-wt-photon-vertex" />）。根拠は電流 $\hat j^\mu=\hat{\bar\psi}\gamma^\mu\hat\psi$ の保存と等時刻交換関係 $[\hat j^0(t,\boldsymbol x),\hat\psi(t,\boldsymbol y)]=-\delta^3(\boldsymbol x-\boldsymbol y)\hat\psi(t,\boldsymbol y)$ で、これから
$$
\partial_\mu\langle\Omega|T\{\hat j^\mu(x)\hat\psi(x_1)\hat{\bar\psi}(x_2)\}|\Omega\rangle=\big[-\delta^4(x-x_1)+\delta^4(x-x_2)\big]\langle\Omega|T\{\hat\psi(x_1)\hat{\bar\psi}(x_2)\}|\Omega\rangle
$$
が出ます。時間順序積を微分すると、$\theta$ 関数の微分から交換子が現れるからです。これを Fourier 変換して外線を切り取ると上の恒等式になります。図の言葉では、<Ref to="thm-ward-takahashi" /> の証明の「運動量 $q$ の光子を内線の電子線のすべての位置に挿入した和は、端どうしの差に縮む」という操作を、すべての図で行うことにあたります。外線を質量殻に乗せると右辺の項は 1 粒子極をもたないので消え、<Ref to="physics/qft/maxwell-field-quantization#thm-ward-identity" /> の S 行列版が得られます。
</Remark>

<Corollary id="cor-z1-equals-z2" title="Z₁ = Z₂">
<Ref to="thm-ward-takahashi" /> が成り立つ次数で $Z_1=Z_2$ です。したがって物理的電荷は $e=Z_3^{1/2}e_0$ で、電子の場のくりこみにも頂点補正にもよりません。
</Corollary>

<Proof of="cor-z1-equals-z2">
<Ref to="thm-ward-takahashi" /> の両辺を $q$ で展開して 1 次の項を比べると、$\Gamma^\mu(p,p)=\gamma^\mu-\frac{\partial\Sigma}{\partial p_\mu}$ を得ます。$\Sigma$ は Lorentz 不変性とパリティから $\Sigma(p)=A(p^2)\not{p}+B(p^2)m$ の形なので
$$
\frac{\partial\Sigma}{\partial p_\mu}=A\gamma^\mu+2p^\mu\big(A'\not{p}+B'm\big)
$$
です。$p^2=m^2$ の $u(p)$ で挟むと、$\not{p}u=mu$ と Gordon 恒等式（<Ref to="lem-gordon-identity" />）の $q=0$ の場合 $\bar u\gamma^\mu u=\frac{p^\mu}{m}\bar uu$ から、$\bar u\frac{\partial\Sigma}{\partial p_\mu}u=\big[A+2m^2(A'+B')\big]\bar u\gamma^\mu u$ です。一方 $\Sigma$ を $\not{p}$ の関数 $A(\not{p}^2)\not{p}+B(\not{p}^2)m$ とみると $\frac{d\Sigma}{d\not{p}}\big|_{\not{p}=m}=A+2m^2A'+2m^2B'$ で、両者は一致します。よって $Z_1^{-1}=1-\frac{d\Sigma}{d\not{p}}\big|_{\not{p}=m}=Z_2^{-1}$ です。
</Proof>

<Remark id="rem-charge-universality" title="電荷のくりこみの普遍性">
<Ref to="cor-z1-equals-z2" /> により、電荷のくりこみは光子の自己エネルギーだけで決まります。電子と陽子のように相互作用の仕方がまったく違う粒子でも、裸の電荷が等しければ（あるいは整数比なら）物理的電荷も等しい（同じ整数比）ことが、この恒等式から保証されます。水素原子が中性であることは、強い相互作用による陽子の複雑な構造とは無関係です。
</Remark>

## 4. 真空偏極

この節では、まず光子の自己エネルギー（真空偏極）が横波になることを示し、光子が補正を受けても質量ゼロのままであることを確かめます。証明では、補正を何回挿入した項も同じ形になるので、等比級数として足し上げます。途中に出る $P^{\mu\nu}$ は、$q^\mu$ に比例する成分を取り除く射影です。そのあとで、1 ループの真空偏極を実際に計算します。

<Proposition id="prop-vacuum-polarization-transverse" title="真空偏極は横波である">
<Ref to="thm-ward-takahashi" /> の下で、$\Pi^{\mu\nu}(q)=(q^2g^{\mu\nu}-q^\mu q^\nu)\Pi(q^2)$ と書けるスカラー関数 $\Pi(q^2)$ が存在し、$q=0$ で正則です。このとき光子の完全な伝播関数は
$$
\frac{-i}{q^2\big(1-\Pi(q^2)\big)}\Big(g_{\mu\nu}-\frac{q_\mu q_\nu}{q^2}\Big)+(q_\mu q_\nu\text{ に比例する項})
$$
で、$q^2=0$ の極は動かず（光子は質量をもたない）、留数は $Z_3=(1-\Pi(0))^{-1}$ です。
</Proposition>

<Proof of="prop-vacuum-polarization-transverse">
Lorentz 共変性から $\Pi^{\mu\nu}=A(q^2)g^{\mu\nu}+B(q^2)q^\mu q^\nu$ です。$q_\mu\Pi^{\mu\nu}=(A+q^2B)q^\nu=0$ より $A=-q^2B$ で、$\Pi:=-B$ と置けば主張の形です。1 ループの $\Pi_2$ が $q=0$ で正則であることは <Ref to="thm-vacuum-polarization-one-loop" /> の明示式から分かります（$m\gt0$ なら $\Delta\gt0$）。完全な伝播関数は $D=\frac{-ig}{q^2}$（ファインマンゲージ）と $i\Pi^{\mu\nu}$ の鎖の和です。$P^{\mu\nu}=g^{\mu\nu}-q^\mu q^\nu/q^2$ は $P^2=P$ を満たす射影なので
$$
D\,(i\Pi)\,D=\frac{-i}{q^2}\,i\,q^2\Pi\,P\,\frac{-i}{q^2}=\frac{-i}{q^2}\Pi P
$$
となり、$n$ 回挿入した項は $\frac{-i}{q^2}\Pi^nP$ です（$q_\mu q_\nu$ の部分は $P$ に吸収されず別に残ります）。等比級数を足すと主張の式です。$\Pi$ が $q^2=0$ で極をもたないことが、極が動かないことの根拠です。
</Proof>

<Theorem id="thm-vacuum-polarization-one-loop" title="1 ループの真空偏極">
次元正則化（$d=4-\varepsilon$）で、1 ループの真空偏極は
$$
\Pi_2(q^2)=-\frac{8e^2}{(4\pi)^{d/2}}\int_0^1dx\,x(1-x)\frac{\Gamma(2-\frac d2)}{\Delta^{2-d/2}}=-\frac{2\alpha}{\pi}\int_0^1dx\,x(1-x)\Big[\frac2\varepsilon-\gamma_E+\ln4\pi-\ln\frac{\Delta}{\mu^2}\Big]+O(\varepsilon)
$$
です。ここで $\Delta=m^2-x(1-x)q^2$ で、$q^2\lt4m^2$ では $\Delta\gt0$ です。
</Theorem>

<Proof of="thm-vacuum-polarization-one-loop">
<Ref to="thm-ward-takahashi" /> の証明の表式で $(-ie)^2\cdot i^2=e^2$ なので
$$
i\Pi_2^{\mu\nu}(q)=-e^2\int\frac{d^dk}{(2\pi)^d}\frac{\mathrm{tr}\big[\gamma^\mu(\not{k}+m)\gamma^\nu(\not{k}+\not{q}+m)\big]}{(k^2-m^2)((k+q)^2-m^2)}
$$
です。<Ref to="physics/qft/feynman-rules-lsz#lem-gamma-traces" /> より、トレースは $4\big[k^\mu(k+q)^\nu+k^\nu(k+q)^\mu-g^{\mu\nu}(k\cdot(k+q)-m^2)\big]$ です。<Ref to="lem-feynman-parameters" /> で分母をまとめると $x((k+q)^2-m^2)+(1-x)(k^2-m^2)=(k+xq)^2-\Delta$ なので、$\ell=k+xq$ と置きます。$k=\ell-xq$ を代入し、<Ref to="lem-symmetric-integration" /> で $\ell$ の 1 次の項を落とし $\ell^\mu\ell^\nu\mapsto\ell^2g^{\mu\nu}/d$ とすると、分子は
$$
4\Big[-\Big(1-\frac2d\Big)g^{\mu\nu}\ell^2-2x(1-x)q^\mu q^\nu+g^{\mu\nu}\big(m^2+x(1-x)q^2\big)\Big]
$$
です。<Ref to="thm-master-integrals" /> の $n=2$ の 2 式を比べると、$\Gamma(2-\frac d2)=(1-\frac d2)\Gamma(1-\frac d2)$ より
$$
-\Big(1-\frac2d\Big)\int\frac{d^d\ell}{(2\pi)^d}\frac{\ell^2}{(\ell^2-\Delta)^2}=-\Delta\int\frac{d^d\ell}{(2\pi)^d}\frac{1}{(\ell^2-\Delta)^2}
$$
が成り立ちます。したがって $g^{\mu\nu}$ の係数は $-\Delta+m^2+x(1-x)q^2=2x(1-x)q^2$ となり、分子は実質的に $8x(1-x)(q^2g^{\mu\nu}-q^\mu q^\nu)$ です。$\int\frac{d^d\ell}{(2\pi)^d}\frac{1}{(\ell^2-\Delta)^2}=\frac{i\,\Gamma(2-d/2)}{(4\pi)^{d/2}\Delta^{2-d/2}}$ を掛けて $i$ で割ると第 1 の等式を得ます。$\frac{8e^2}{16\pi^2}=\frac{2\alpha}{\pi}$ と第 2 節末の展開から第 2 の等式です（$e\mapsto e\mu^{\varepsilon/2}$ により $\ln\Delta$ が $\ln(\Delta/\mu^2)$ になります）。
</Proof>

計算の途中で $g^{\mu\nu}$ の係数が $\Delta$ を含む項と打ち消し合い、$q^2$ に比例する形が残りました。これが <Ref to="thm-ward-takahashi" /> の横波性の具体的な現れです。次の例は、平行移動不変性をもたない正則化ではこの打ち消しが起こらないことを示します。

<Example id="ex-cutoff-breaks-ward" title="運動量カットオフは Ward 恒等式を破る">
$q=0$ の真空偏極を、Wick 回転したあと $|k_E|\lt\Lambda$ で切る 4 次元のカットオフで計算します。$q=0$ ではトレースは $4[2k^\mu k^\nu-g^{\mu\nu}(k^2-m^2)]$ で、4 次元の対称積分 $k^\mu k^\nu\mapsto\frac14g^{\mu\nu}k^2$ により
$$
i\Pi^{\mu\nu}(0)=-4e^2g^{\mu\nu}\int\frac{d^4k}{(2\pi)^4}\frac{-\frac12k^2+m^2}{(k^2-m^2)^2}=-4ie^2g^{\mu\nu}\int_{|k_E|\lt\Lambda}\frac{d^4k_E}{(2\pi)^4}\frac{\frac12k_E^2+m^2}{(k_E^2+m^2)^2}
$$
です。$u=k_E^2$ と置くと $\frac{d^4k_E}{(2\pi)^4}=\frac{u\,du}{16\pi^2}$、$\frac{u(\frac12u+m^2)}{(u+m^2)^2}=\frac12-\frac{m^4}{2(u+m^2)^2}$ なので
$$
\Pi^{\mu\nu}(0)=-\frac{e^2}{8\pi^2}\big(\Lambda^2-m^2\big)g^{\mu\nu}+O\Big(\frac{m^4}{\Lambda^2}\Big)
$$
となります。横波なら $\Pi^{\mu\nu}(0)=0$ のはずなので、これは <Ref to="thm-ward-takahashi" /> の破れです。<Ref to="prop-vacuum-polarization-transverse" /> の証明と同じ和をとると伝播関数の分母が $q^2-C$（$C=-\frac{e^2}{8\pi^2}(\Lambda^2-m^2)$）となり、光子の極が $q^2=0$ からずれます。次元正則化では同じ量が $0$ になります（第 2 節の $n=1,2$ の公式で $(\frac2d-1)\int\frac{1}{k^2-m^2}+\frac2dm^2\int\frac{1}{(k^2-m^2)^2}=0$ を確かめられます）。壊れたのは「積分変数の平行移動で積分が変わらない」という仮定で、カットオフは 2 次発散の積分でこれを満たしません。
</Example>

## 5. 電荷のくりこみと走る結合定数

この節では、真空偏極の結果から、観測するエネルギーによって値が変わる結合定数（[走る結合定数](/physics/qft/dimensional-regularization-rg#def-running-coupling)）を導きます。発散する部分は電荷のくりこみに吸収されるので、観測にかかるのは $q^2$ に依存する有限の差だけです。

<Ref to="cor-z1-equals-z2" /> と <Ref to="prop-vacuum-polarization-transverse" /> から、$e^2=Z_3e_0^2=\frac{e_0^2}{1-\Pi(0)}$ です。$\Pi_2(0)$ は $\frac2\varepsilon$ の極をもつので裸の電荷と物理的電荷の比は発散しますが、観測量には差 $\hat\Pi_2(q^2)=\Pi_2(q^2)-\Pi_2(0)$ だけが現れます。

<Theorem id="thm-running-alpha" title="有効結合定数">
2 つの荷電粒子が光子 1 個を交換する振幅の、真空偏極の鎖をすべて含めた光子の伝播関数の $g_{\mu\nu}$ の部分は、$\alpha^2$ の次数の誤差を除いて $\frac{-ig_{\mu\nu}}{q^2}\frac{e^2}{1-\hat\Pi_2(q^2)}$ です（$e_0^2$ を $e^2$ に吸収）。ここで
$$
\hat\Pi_2(q^2)=-\frac{2\alpha}{\pi}\int_0^1dx\,x(1-x)\ln\frac{m^2}{m^2-x(1-x)q^2}
$$
は有限です。$\alpha_{\mathrm{eff}}(q^2)=\frac{\alpha}{1-\hat\Pi_2(q^2)}$ と置くと、空間的な運動量移行 $-q^2\gg m^2$ で
$$
\alpha_{\mathrm{eff}}(q^2)=\frac{\alpha}{1-\frac{\alpha}{3\pi}\ln\frac{-q^2}{Am^2}}+O\Big(\alpha^2\frac{m^2}{-q^2}\ln\frac{-q^2}{m^2}\Big),\qquad A=e^{5/3}
$$
です。
</Theorem>

<Proof of="thm-running-alpha">
$\frac{e_0^2}{1-\Pi(q^2)}=\frac{e_0^2}{1-\Pi(0)}\cdot\frac{1-\Pi(0)}{1-\Pi(q^2)}=\frac{e^2}{1-\hat\Pi(q^2)}+O(e^6)$ です。<Ref to="thm-vacuum-polarization-one-loop" /> の差をとると極と定数が消え、$\hat\Pi_2$ の式を得ます。$Q^2=-q^2\gg m^2$ では $\ln\frac{m^2}{m^2+x(1-x)Q^2}=-\ln\frac{Q^2}{m^2}-\ln x(1-x)+O\big(\frac{m^2}{x(1-x)Q^2}\big)$ です。$\int_0^1x(1-x)dx=\frac16$、$\int_0^1x(1-x)\ln x\,dx=-\frac14+\frac19=-\frac5{36}$（対称性から $\ln(1-x)$ も同じ値）を使うと
$$
\hat\Pi_2=-\frac{2\alpha}{\pi}\Big[-\frac16\ln\frac{Q^2}{m^2}+\frac5{18}\Big]=\frac{\alpha}{3\pi}\Big(\ln\frac{Q^2}{m^2}-\frac53\Big)
$$
です。端点 $x\to0,1$ の近くでは展開が使えませんが、その幅は $m^2/Q^2$ 程度で、被積分関数は $x(1-x)\ln$ で抑えられるので、誤差は主張の大きさです。
</Proof>

$\alpha_{\mathrm{eff}}$ は $Q^2$ とともに増えます。遠くからは電荷が仮想対の分極で**遮蔽**され、近づくほど裸に近い電荷が見えるからです。この $Q$ 依存性をスケール $\mu$ の関数として表したものが β関数です（<Ref to="physics/qft/renormalization#def-beta-anomalous" />）。

<Theorem id="thm-qed-beta-one-loop" title="QED の 1 ループのβ関数">
質量 $m$ の Dirac フェルミオン 1 種類の QED を最小引き算（MS）でくりこみ、$e_0=\mu^{\varepsilon/2}Z_1Z_2^{-1}Z_3^{-1/2}e(\mu)$ でくりこまれた結合 $e(\mu)$ を定めると
$$
\beta(e)=\mu\frac{de}{d\mu}=\frac{e^3}{12\pi^2}+O(e^5),\qquad \mu\frac{d\alpha}{d\mu}=\frac{2\alpha^2}{3\pi}+O(\alpha^3)
$$
です。$n_f$ 種類のフェルミオン（電荷 $Q_ie$）があれば係数は $\sum_iQ_i^2$ 倍になります。
</Theorem>

<Proof of="thm-qed-beta-one-loop">
MS では $Z_i$ は $\frac1\varepsilon$ の極だけを含みます。<Ref to="thm-ward-takahashi" /> は次元正則化の各次数で成り立つので、極の部分にも <Ref to="cor-z1-equals-z2" /> の議論が通り $Z_1=Z_2$ です（直接の計算は演習で確かめます）。$Z_3=1+\Pi(0)+O(e^4)$ のうち極だけをとると、<Ref to="thm-vacuum-polarization-one-loop" /> から $\Pi_2$ の極は $-\frac{2\alpha}{\pi}\cdot\frac16\cdot\frac2\varepsilon=-\frac{e^2}{6\pi^2\varepsilon}$ なので、$Z_3=1-\frac{e^2}{6\pi^2\varepsilon}$ です。よって
$$
e_0=\mu^{\varepsilon/2}e\Big(1+\frac{e^2}{12\pi^2\varepsilon}+O(e^4)\Big)
$$
です。左辺は $\mu$ によらないので、$\mu\frac{d}{d\mu}$ をとって $\beta=-\frac\varepsilon2e+b\,e^3+O(e^5)$ と置くと
$$
0=\frac\varepsilon2e\Big(1+\frac{e^2}{12\pi^2\varepsilon}\Big)+\Big(-\frac\varepsilon2e+be^3\Big)\Big(1+\frac{3e^2}{12\pi^2\varepsilon}\Big)+O(e^5)
$$
です。$e$ の 1 次は打ち消し、$e^3$ の係数から $\frac{1}{24\pi^2}+b-\frac{3}{24\pi^2}=0$、すなわち $b=\frac{1}{12\pi^2}$ を得ます。$\varepsilon\to0$ で主張の式です。$\alpha$ については $\mu\frac{d\alpha}{d\mu}=\frac{2e\beta}{4\pi}=\frac{e^4}{24\pi^3}=\frac{2\alpha^2}{3\pi}$ です。フェルミオンごとに真空偏極のループが 1 つ増え、頂点に $Q_i$ が 2 回掛かるので係数は $\sum Q_i^2$ 倍です。
</Proof>

この係数は <Ref to="thm-running-alpha" /> と整合します。$\alpha_{\mathrm{eff}}^{-1}=\alpha^{-1}-\frac{1}{3\pi}\ln\frac{Q^2}{Am^2}$ を $\ln Q$ で微分すると $\frac{d\alpha_{\mathrm{eff}}}{d\ln Q}=\frac{2\alpha_{\mathrm{eff}}^2}{3\pi}$ で、MS の β関数と同じ係数です。1 ループの係数はスキームによりません。

<Example id="ex-landau-pole" title="Landau 極">
<Ref to="thm-qed-beta-one-loop" /> の方程式を解くと $\frac{1}{\alpha(\mu)}=\frac{1}{\alpha(\mu_0)}-\frac{2}{3\pi}\ln\frac{\mu}{\mu_0}$ で、右辺は
$$
\mu_{\mathrm L}=\mu_0\exp\Big(\frac{3\pi}{2\alpha(\mu_0)}\Big)
$$
で $0$ になり、$\alpha(\mu)$ が発散します。これを **Landau 極**と呼びます。電子だけを考えて $\mu_0=m_e$、$\alpha=1/137.036$ とすると指数は $\frac{3\pi\times137.036}{2}\approx645.8$ で、$\mu_{\mathrm L}\approx10^{280.5}m_e\approx10^{277}\ \mathrm{GeV}$ です。Planck 質量 $\sim10^{19}\ \mathrm{GeV}$ よりはるかに大きいので、実際の物理には影響しません。ただし原理的には、1 ループの β関数を信じる限り QED は任意に高いエネルギーまで定義できる理論ではありません。$\alpha$ が大きくなる領域では摂動論そのものが使えないので、この極が実在するかどうかは 1 ループでは決まりません。格子上の数値的研究は、QED が紫外で「自明」（くりこまれた結合が $0$ でなければ切断を外せない）であることを示唆しています。非可換ゲージ理論では β関数の符号が逆になり、この問題が起こりません（[非可換ゲージ理論の量子化](/physics/qft/nonabelian-gauge-quantization)）。
</Example>

## 6. 頂点補正と電子の異常磁気モーメント

<Lemma id="lem-gordon-identity" title="Gordon 恒等式">
$p^2=p'^2=m^2$、$q=p'-p$、$\sigma^{\mu\nu}=\frac i2[\gamma^\mu,\gamma^\nu]$ とすると
$$
\bar u(p')\gamma^\mu u(p)=\bar u(p')\Big[\frac{(p'+p)^\mu}{2m}+\frac{i\sigma^{\mu\nu}q_\nu}{2m}\Big]u(p)
$$
が成り立ちます。
</Lemma>

<Proof of="lem-gordon-identity">
$\bar u(p')\not{p}'=m\bar u(p')$、$\not{p}u(p)=mu(p)$ より $\bar u(p')\gamma^\mu u(p)=\frac1{2m}\bar u(p')(\not{p}'\gamma^\mu+\gamma^\mu\not{p})u(p)$ です。$\gamma^\nu\gamma^\mu=g^{\nu\mu}-i\sigma^{\nu\mu}$ を使うと $\not{p}'\gamma^\mu+\gamma^\mu\not{p}=(p'+p)^\mu-ip'_\nu\sigma^{\nu\mu}-ip_\nu\sigma^{\mu\nu}=(p'+p)^\mu+i\sigma^{\mu\nu}(p'-p)_\nu$ です。
</Proof>

<Definition id="def-form-factors" title="形状因子">
質量殻上の $p,p'$ について、頂点関数を
$$
\bar u(p')\Gamma^\mu(p',p)u(p)=\bar u(p')\Big[\gamma^\mu F_1(q^2)+\frac{i\sigma^{\mu\nu}q_\nu}{2m}F_2(q^2)\Big]u(p)
$$
と分解し、$F_1,F_2$ を**形状因子**と呼びます。Lorentz 共変性、[パリティ](/physics/qft/dirac-field#def-discrete-symmetries)、<Ref to="physics/qft/maxwell-field-quantization#thm-ward-identity" /> から、$\bar u(p')\Gamma^\mu u(p)$ はこの 2 項で尽くされます（$\gamma^\mu,(p+p')^\mu,q^\mu$ の 3 つの構造のうち $q^\mu$ は Ward 恒等式で禁止され、$(p+p')^\mu$ は <Ref to="lem-gordon-identity" /> で $\gamma^\mu$ と $\sigma^{\mu\nu}q_\nu$ に直せます）。最低次は $F_1=1$、$F_2=0$ です。
</Definition>

<Ref to="def-qed-renormalization-constants" /> の物理的電荷の定義により、くりこまれた頂点では $F_1(0)=1$ です。これは $q\to0$ の極限で電子が電荷 $e$ の点電荷として振る舞うことを表しています。

ここで知りたいのは、頂点補正が電子の磁石としての強さをどれだけ変えるかです。次の命題では、静的な磁場をかけたときの散乱振幅を、磁気モーメントをもつ粒子のポテンシャル散乱と見比べて、$g$ 因子を $F_1(0)$ と $F_2(0)$ で表します。

<Proposition id="prop-g-factor" title="g 因子と形状因子">
時間によらない外部ベクトルポテンシャル $\boldsymbol A_{\mathrm{cl}}(\boldsymbol x)$ による電子の散乱振幅を、非相対論的極限で $\boldsymbol q$ の 1 次まで展開すると、ポテンシャル $V(\boldsymbol x)=-\langle\boldsymbol\mu\rangle\cdot\boldsymbol B(\boldsymbol x)$ の Born 近似と一致します。ここで $\boldsymbol\mu=g\frac{e}{2m}\boldsymbol S$、$\boldsymbol S$ はスピン演算子、
$$
g=2\big[F_1(0)+F_2(0)\big]=2+2F_2(0)
$$
です。
</Proposition>

<Proof of="prop-g-factor">
静的な外場なので $A^0_{\mathrm{cl}}=0$、エネルギーは保存し $q^0=0$ です。頂点 $-ie\gamma^\mu$ の外場への結合から、振幅は $i\mathcal M=-ie\,\bar u(p')\Gamma^\mu(p',p)u(p)\,\tilde A_{\mathrm{cl},\mu}(\boldsymbol q)$ で、$\tilde A_{\mathrm{cl},\mu}=(0,-\tilde A^i_{\mathrm{cl}})$ なので $i\mathcal M=ie\,\bar u(p')\Gamma^iu(p)\,\tilde A^i_{\mathrm{cl}}(\boldsymbol q)$ です。ここで $\tilde A(\boldsymbol q)=\int d^3x\,e^{-i\boldsymbol q\cdot\boldsymbol x}A(\boldsymbol x)$、$\boldsymbol q=\boldsymbol p'-\boldsymbol p$ です。

**スピンに依存する部分。** <Ref to="def-form-factors" /> の $\gamma^iF_1$ に <Ref to="lem-gordon-identity" /> を使うと
$$
\bar u(p')\Gamma^iu(p)=\bar u(p')\Big[\frac{(p+p')^i}{2m}F_1+\frac{i\sigma^{i\nu}q_\nu}{2m}\big(F_1+F_2\big)\Big]u(p)
$$
です。第 1 項は非相対論的極限の最低次で $\bar u(p')u(p)\approx2m\,\xi'^\dagger\xi$ に比例し、スピンによりません（軌道運動と外場の結合 $\boldsymbol p\cdot\boldsymbol A$ の項です）。第 2 項で $q_0=0$、$q_j=-q^j$、$\sigma^{ij}=\epsilon^{ijk}\Sigma^k$（$\Sigma^k=\mathrm{diag}(\sigma^k,\sigma^k)$）なので $i\sigma^{i\nu}q_\nu=-i\epsilon^{ijk}q^j\Sigma^k$ です。カイラル表示で $u(p)\approx\sqrt m(\xi,\xi)^{\mathsf T}$ とすると $\bar u(p')\Sigma^ku(p)\approx2m\,\xi'^\dagger\sigma^k\xi$ なので、$\boldsymbol q$ の 1 次までで（$F_i(q^2)=F_i(0)+O(\boldsymbol q^2)$）
$$
i\mathcal M_{\mathrm{spin}}=ie\cdot\frac{F_1(0)+F_2(0)}{2m}\big(-i\epsilon^{ijk}q^j\big)\tilde A^i_{\mathrm{cl}}\cdot2m\,\xi'^\dagger\sigma^k\xi=e\big[F_1(0)+F_2(0)\big]\epsilon^{ijk}q^j\tilde A^i_{\mathrm{cl}}\,\xi'^\dagger\sigma^k\xi
$$
です。$\boldsymbol B=\nabla\times\boldsymbol A_{\mathrm{cl}}$ の Fourier 成分は $\tilde B^k=i\epsilon^{kab}q^a\tilde A^b_{\mathrm{cl}}$ で、$\epsilon^{ijk}q^j\tilde A^i=-\epsilon^{kab}q^a\tilde A^b=i\tilde B^k$ ですから、
$$
i\mathcal M_{\mathrm{spin}}=ie\big[F_1(0)+F_2(0)\big]\,\xi'^\dagger\boldsymbol\sigma\xi\cdot\tilde{\boldsymbol B}(\boldsymbol q)
$$
です。

**Born 近似との比較。** ポテンシャル $V(\boldsymbol x)$ の Born 近似の振幅は、相対論的な規格化（状態あたり $2m$）で $i\mathcal M=-i\tilde V(\boldsymbol q)\cdot2m$ です。$V=-\langle\boldsymbol\mu\rangle\cdot\boldsymbol B$、$\boldsymbol\mu=g\frac{e}{2m}\frac{\boldsymbol\sigma}{2}$ なら $-i\tilde V\cdot2m=i\frac{ge}{2}\,\xi'^\dagger\boldsymbol\sigma\xi\cdot\tilde{\boldsymbol B}$ です。上の式と比べると $\frac g2=F_1(0)+F_2(0)$、すなわち $g=2[F_1(0)+F_2(0)]$ で、くりこまれた頂点の $F_1(0)=1$ を代入すれば $g=2+2F_2(0)$ です。
</Proof>

<Remark id="rem-g-factor-proof" title="g 因子の導出の所在">
導出は Peskin–Schroeder 第 6 章 §6.2 にあります。筋道は次のとおりです。振幅は $i\mathcal M=-ie\bar u(p')\Gamma^i(p',p)u(p)\tilde A^i_{\mathrm{cl}}(\boldsymbol q)$ で、非相対論的極限 $u(p)\approx\sqrt m(\xi,\xi)^{\mathsf T}$ では $\gamma^iF_1$ の項のうちスピンに依存する部分と $\sigma^{ij}q_jF_2$ の項がともに $\xi'^\dagger\sigma^kq^j\epsilon^{ijk}\xi$ に比例し、合わせて $F_1(0)+F_2(0)$ の係数をもちます。これを磁場 $\boldsymbol B=\nabla\times\boldsymbol A_{\mathrm{cl}}$ の Fourier 成分で書き直し、ポテンシャル $-\langle\boldsymbol\mu\rangle\cdot\boldsymbol B$ の Born 振幅と比べると、$\boldsymbol\mu$ の係数として $\frac em[F_1(0)+F_2(0)]$ が読み取れます。Dirac 方程式の $g=2$ は $F_1=1$、$F_2=0$ の場合です。
</Remark>

<Theorem id="thm-schwinger" title="Schwinger の異常磁気モーメント">
QED の $\alpha$ の 1 次で、$F_2(q^2)$ は紫外・赤外ともに有限で
$$
F_2(q^2)=\frac{\alpha}{2\pi}\int_0^1dx\,dy\,dz\;\delta(x+y+z-1)\frac{2m^2z(1-z)}{m^2(1-z)^2-q^2xy}
$$
です。特に $F_2(0)=\frac{\alpha}{2\pi}$ で、電子の異常磁気モーメントは $a_e\equiv\frac{g-2}{2}=\frac{\alpha}{2\pi}+O(\alpha^2)$ です。
</Theorem>

<Proof of="thm-schwinger">
ファインマンゲージで、光子に赤外正則化の質量 $\mu$ を与えます。<Ref to="thm-ward-takahashi" /> の証明の $\delta\Gamma^\mu$ の式でループ運動量を電子側の $k=p-l$、$k'=k+q$ にとり、<Ref to="lem-symmetric-integration" /> の $\gamma^\nu\not{a}\gamma^\mu\not{b}\gamma_\nu$ の公式を使うと
$$
\bar u(p')\delta\Gamma^\mu u(p)=2ie^2\int\frac{d^4k}{(2\pi)^4}\frac{\bar u(p')\big[\not{k}\gamma^\mu\not{k}'+m^2\gamma^\mu-2m(k+k')^\mu\big]u(p)}{((k-p)^2-\mu^2)(k'^2-m^2)(k^2-m^2)}
$$
です（$(-i)(-ie)^2i^2=-ie^2$ に $\gamma$ 行列の縮約の $-2$ が掛かります）。$d$ 次元の補正 $(4-d)\not{a}\gamma^\mu\not{b}$ は、発散する $\ell^2$ の項に掛かったときだけ有限の寄与を残し、それは $\gamma^\mu$ に比例するので $F_1$ にしか効きません。以下 $F_2$ に必要な部分だけを 4 次元で計算します。

**分母。** <Ref to="lem-feynman-parameters" /> で $x(k^2-m^2)+y(k'^2-m^2)+z((k-p)^2-\mu^2)$ とまとめ、$p^2=p'^2=m^2$ から出る $p\cdot q=-q^2/2$ を使って平方完成すると、$\ell=k+yq-zp$ について分母は $(\ell^2-\Delta)^3$、
$$
\Delta=-xyq^2+(1-z)^2m^2+z\mu^2
$$
です。$\Delta$ は $x\leftrightarrow y$ について対称です。

**分子。** $k=\ell+a$、$k'=\ell+b$ で $a=(1-x)p-yp'$、$b=-xp+(1-y)p'$ です（$q=p'-p$、$x+y+z=1$ で書き直しました）。$\ell$ の 1 次の項は消え、$\ell$ の 2 次の項は $\gamma^\mu$ に比例します。残りは $\bar u(p')[\not{a}\gamma^\mu\not{b}+m^2\gamma^\mu-2m(a+b)^\mu]u(p)$ で、$\bar u(p')\not{p}'=m\bar u(p')$、$\not{p}u(p)=mu(p)$ から
$$
\bar u\not{a}\gamma^\mu\not{b}u=\bar u\big[(1-x)\not{p}-ym\big]\gamma^\mu\big[-xm+(1-y)\not{p}'\big]u
$$
です。展開した 4 項を、反交換関係 $\not{p}\gamma^\mu=2p^\mu-\gamma^\mu\not{p}$、$\gamma^\mu\not{p}'=2p'^\mu-\not{p}'\gamma^\mu$ と、$\bar u\not{p}\gamma^\mu\not{p}'u=\bar u(m-\not{q})\gamma^\mu(m+\not{q})u$ に Gordon 恒等式を適用した結果 $\bar u[(q^2-3m^2)\gamma^\mu+2m(p+p')^\mu]u$ を使って、$\gamma^\mu$、$(p+p')^\mu$、$q^\mu$ に整理します。$(p+p')^\mu\bar uu$ の係数を集めると
$$
\frac12\Big[4m(1-x)(1-y)-2m\big\{x(1-x)+y(1-y)\big\}-2m(2-2x-2y)\Big]=mz(z-1)
$$
（$(1-x)(1-y)=z+xy$、$x(1-x)+y(1-y)=z(1-z)+2xy$ を使いました）、$q^\mu\bar uu$ の係数は $m(z-2)(x-y)$ です。後者は $x\leftrightarrow y$ について反対称で $\Delta$ は対称なので、積分すると $0$ です（<Ref to="physics/qft/maxwell-field-quantization#thm-ward-identity" /> と整合します）。

**形状因子。** <Ref to="lem-gordon-identity" /> で $(p+p')^\mu=2m\gamma^\mu-i\sigma^{\mu\nu}q_\nu$（$u$ に挟んだ意味で）と置き換えると、$i\sigma^{\mu\nu}q_\nu/(2m)$ の係数は $2m^2z(1-z)$ です。<Ref to="thm-master-integrals" /> の $n=3$、$d=4$ から $\int\frac{d^4\ell}{(2\pi)^4}\frac{1}{(\ell^2-\Delta)^3}=\frac{-i}{32\pi^2\Delta}$ なので
$$
F_2(q^2)=2ie^2\cdot2\int dx\,dy\,dz\,\delta(x+y+z-1)\,\frac{-i}{32\pi^2\Delta}\,2m^2z(1-z)=\frac{\alpha}{2\pi}\int dx\,dy\,dz\,\delta(\cdots)\frac{2m^2z(1-z)}{\Delta}
$$
です。この積分は $\mu\to0$ でも有限です（$z\to1$ で分子が $1-z$ の 1 次で消え、$\Delta\approx(1-z)^2m^2$ の 1 次の零点と打ち消し合います）。$F_2$ は $\gamma^\mu$ に比例しないので、$Z_1$ による引き算の影響も受けません。$q^2=0$ では $x$ について積分（長さ $1-z$）して
$$
F_2(0)=\frac{\alpha}{2\pi}\int_0^1dz\,(1-z)\frac{2z(1-z)}{(1-z)^2}=\frac{\alpha}{2\pi}\int_0^12z\,dz=\frac{\alpha}{2\pi}
$$
を得ます。<Ref to="prop-g-factor" /> より $a_e=F_2(0)$ です。
</Proof>

<Remark id="rem-ae-experiment" title="実験との比較">
$\frac{\alpha}{2\pi}=0.0011614$ に対し、実測値は $a_e=0.00115965218\ldots$ です。差の大部分は $\alpha^2$ の項 $-0.328\frac{\alpha^2}{\pi^2}$ で説明され、現在は $\alpha^5$ までの QED の計算と測定が 10 桁以上一致しています。逆に、$a_e$ の測定は $\alpha$ の最も精密な決定法の 1 つになっています。Schwinger の 1 ループの値がすでに 0.15% の精度で合っていることは、くりこまれた摂動論が物理的な予言を与えることの最初の決定的な証拠でした。
</Remark>

## 7. 赤外発散と Bloch–Nordsieck の相殺

<Ref to="thm-schwinger" /> の計算の $\gamma^\mu$ の項、すなわち $F_1$ は事情が違います。$\mu\to0$ で $z\to1$ の角の近くから $\ln\mu$ に比例する発散が生じます。

次の命題の証明は長いですが、やることは 2 段です。まず $F_1$ の積分表示を求め、次に $\mu\to0$ で発散する部分だけを取り出して、残りが有限に収まることを評価します。

<Proposition id="prop-ir-divergence-f1" title="F₁ の赤外発散">
光子質量 $\mu$ で正則化すると、$\mu\to0$ で
$$
F_1(q^2)=1-\frac{\alpha}{2\pi}f_{\mathrm{IR}}(q^2)\ln\frac{m^2}{\mu^2}+(\mu\text{ によらない項})+O(\alpha^2),\qquad f_{\mathrm{IR}}(q^2)=\int_0^1d\xi\,\frac{m^2-q^2/2}{m^2-q^2\xi(1-\xi)}-1
$$
です。空間的な $q^2\lt0$ で $f_{\mathrm{IR}}\gt0$、$-q^2\to\infty$ で $f_{\mathrm{IR}}\approx\ln\frac{-q^2}{m^2}$ です。
</Proposition>

<Proof of="prop-ir-divergence-f1">
$q^2\lt4m^2$（対生成の閾値より下。散乱の $q^2\le0$ を含む）とし、$c(\xi)=m^2-q^2\xi(1-\xi)$ と置きます。$\xi(1-\xi)\le\frac14$ なので $c(\xi)\ge\min(m^2,m^2-\frac{q^2}4)\gt0$ です。

**$F_1$ の表示。** <Ref to="thm-schwinger" /> の証明の記号で、分子の $\gamma^\mu$ の係数を集めます。$\bar u\not{p}\gamma^\mu u=\bar u(2p^\mu-m\gamma^\mu)u$、$\bar u\gamma^\mu\not{p}'u=\bar u(2p'^\mu-m\gamma^\mu)u$ と、そこで求めた $\bar u\not{p}\gamma^\mu\not{p}'u=\bar u[(q^2-3m^2)\gamma^\mu+2m(p+p')^\mu]u$、$(p+p')^\mu$ の係数 $mz(z-1)$ を Gordon 恒等式で $\gamma^\mu$ に直した $2m^2z(z-1)$ を足すと、$\ell$ によらない分子の $\gamma^\mu$ の係数は
$$
m^2\big[1+x(1-x)+y(1-y)+xy\big]+(1-x)(1-y)(q^2-3m^2)+2m^2z(z-1)=(1-x)(1-y)q^2+(1-4z+z^2)m^2
$$
です（$x(1-x)+y(1-y)=z(1-z)+2xy$、$(1-x)(1-y)=z+xy$ を使いました）。これを $N_1(x,y,z)$ と書きます。$F_2$ のときと同じ係数の計算で、頂点補正の $\gamma^\mu$ の部分は
$$
\delta F_1(q^2)=A(q^2)+\frac{\alpha}{2\pi}\int dF\,\frac{N_1}{\Delta},\qquad\int dF=\int_0^1dx\,dy\,dz\,\delta(x+y+z-1)
$$
です。$A(q^2)$ は $\ell^2$ の項からの寄与で、次元正則化では $\frac2\varepsilon$ の極と $-\int dF\ln\Delta$ に比例する項からなります。くりこまれた頂点は $F_1(q^2)=1+\delta F_1(q^2)-\delta F_1(0)$ です（<Ref to="def-qed-renormalization-constants" /> の条件 $F_1(0)=1$）。

**$\ell^2$ の項は $\mu\to0$ で有限。** $A(q^2)-A(0)$ は極が打ち消し合い、$\int dF\,[\ln\Delta(q^2)-\ln\Delta(0)]$ に比例します。$y=w\xi$、$x=w(1-\xi)$、$z=1-w$ と変数変換すると $\int dF=\int_0^1dw\,w\int_0^1d\xi$、
$$
\Delta=w^2c(\xi)+(1-w)\mu^2
$$
です。$|\ln\Delta|\le|\ln(w^2c_{\min})|+|\ln(c_{\max}+\mu^2)|$（$\mu\le1$）は $w\,dw\,d\xi$ について可積分なので、優収束定理により $\mu\to0$ の極限が存在します。

**$N_1/\Delta$ の項。** $N_1$ は多項式で角 $x=y=0$、$z=1$ での値は $q^2-2m^2$ なので、$|N_1-(q^2-2m^2)|\le Kw$ です。差の部分は $\frac{w\cdot Kw}{w^2c_{\min}}=\frac{K}{c_{\min}}$ で抑えられるので、やはり $\mu\to0$ の極限をもちます。残りは
$$
(q^2-2m^2)\int_0^1d\xi\,I_\mu(\xi),\qquad I_\mu(\xi)=\int_0^1\frac{w\,dw}{w^2c(\xi)+(1-w)\mu^2}
$$
です。$J_\mu(\xi)=\int_0^1\frac{w\,dw}{w^2c+\mu^2}=\frac1{2c}\ln\frac{c+\mu^2}{\mu^2}$ との差は
$$
I_\mu-J_\mu=\int_0^1\frac{\mu^2w^2\,dw}{\big(w^2c+(1-w)\mu^2\big)\big(w^2c+\mu^2\big)}\ \ge0
$$
です。$\mu\le\frac12$ とすると、$w\le\mu$ の部分は分母が $\frac12\mu^2\cdot\mu^2$ 以上なので $\int_0^\mu\frac{2w^2}{\mu^2}dw=\frac{2\mu}3$ 以下、$w\ge\mu$ の部分は分母が $w^4c^2$ 以上なので $\int_\mu^1\frac{\mu^2}{c^2w^2}dw\le\frac{\mu}{c^2}$ 以下で、どちらも $0$ に収束します。よって
$$
I_\mu(\xi)=\frac{1}{2c(\xi)}\ln\frac{m^2}{\mu^2}+\frac1{2c(\xi)}\ln\frac{c(\xi)}{m^2}+o(1)
$$
です。以上をまとめると、$\mu\to0$ で
$$
\delta F_1(q^2)=-\frac{\alpha}{2\pi}\int_0^1d\xi\,\frac{m^2-q^2/2}{c(\xi)}\ln\frac{m^2}{\mu^2}+(\mu\text{ によらない項})+o(1)
$$
です。$q^2=0$ では $c=m^2$ で被積分関数は $1$ なので、$F_1(q^2)=1+\delta F_1(q^2)-\delta F_1(0)$ の $\ln\frac{m^2}{\mu^2}$ の係数は $-\frac{\alpha}{2\pi}\big[\int_0^1\frac{m^2-q^2/2}{c(\xi)}d\xi-1\big]=-\frac{\alpha}{2\pi}f_{\mathrm{IR}}(q^2)$ です。

**$f_{\mathrm{IR}}$ の性質。** $q^2=-Q^2\lt0$ なら被積分関数は $\frac{m^2+Q^2/2}{m^2+Q^2\xi(1-\xi)}$ で、$\xi(1-\xi)\le\frac14\lt\frac12$ より $1$ より大きいので $f_{\mathrm{IR}}\gt0$ です。$\epsilon=m^2/Q^2$ と置くと $f_{\mathrm{IR}}=\big(\frac12+\epsilon\big)\int_0^1\frac{d\xi}{\xi(1-\xi)+\epsilon}-1$ です。$\epsilon\to0$ では $\int_0^1\frac{d\xi}{\xi(1-\xi)+\epsilon}=2\ln\frac1\epsilon+O(1)$（主要部は $\xi=0$ と $\xi=1$ の近くの $\int_0\frac{d\xi}{\xi+\epsilon}$ の 2 つの寄与）なので、$f_{\mathrm{IR}}=\ln\frac{Q^2}{m^2}+O(1)=\ln\frac{-q^2}{m^2}+O(1)$ です。
</Proof>

<Remark id="rem-ir-divergence-proof" title="F₁ の赤外発散の導出">
導出は Peskin–Schroeder 第 6 章 §6.4 にあります。$F_1$ の被積分関数で $\mu$ に敏感なのは $z\approx1$、$x\approx y\approx0$ の角だけです。そこで分子の $x,y,z$ を角の値に置き、$y=(1-z)\xi$、$w=1-z$ と変数変換すると、$\int\frac{d(w^2)}{w^2c+\mu^2}$ の形の積分が現れて $\ln\frac{c}{\mu^2}$ を与えます。$Z_1$ による引き算（$F_1(0)=1$）の分を含めて係数を集めると $f_{\mathrm{IR}}$ になります。
</Remark>

この発散は、$\mu\to0$ で $e^-(p)\to e^-(p')$ の弾性散乱断面積が $\big(\frac{d\sigma}{d\Omega}\big)_0\big[1-\frac\alpha\pi f_{\mathrm{IR}}\ln\frac{m^2}{\mu^2}+\cdots\big]$ となって、負の無限大の補正を受けることを意味します。原因は、無質量の光子をいくらでも低いエネルギーで放出できるのに、その過程を数えていないことです。

<Proposition id="prop-soft-photon-emission" title="軟光子放射の因子化">
過程 $p\to p'$ の振幅を $\mathcal M_0$ とし、さらにエネルギー $|\boldsymbol k|$ の光子 1 個を放射する振幅を $\mathcal M$ とします。$|\boldsymbol k|\to0$ の主要項は
$$
\mathcal M(p\to p'+\gamma(k,\varepsilon))=\mathcal M_0(p\to p')\cdot e\Big(\frac{p'\cdot\varepsilon^*}{p'\cdot k}-\frac{p\cdot\varepsilon^*}{p\cdot k}\Big)
$$
で、偏極について和をとった光子 1 個の放射確率は
$$
d\sigma(p\to p'+\gamma)=d\sigma(p\to p')\cdot\frac{d^3k}{(2\pi)^32|\boldsymbol k|}\,e^2\Big(\frac{2p\cdot p'}{(p\cdot k)(p'\cdot k)}-\frac{m^2}{(p'\cdot k)^2}-\frac{m^2}{(p\cdot k)^2}\Big)
$$
です。
</Proposition>

<Proof of="prop-soft-photon-emission">
出射電子の線から光子が出る図は、$\mathcal M_0$ の出射スピノル $\bar u(p')$ を $\bar u(p')(-ie\gamma^\nu\varepsilon^*_\nu)\frac{i(\not{p}'+\not{k}+m)}{(p'+k)^2-m^2}$ に置き換えたものです。$k\to0$ では分子の $\not{k}$ を落とせ、$\bar u(p')\gamma^\nu(\not{p}'+m)=\bar u(p')[2p'^\nu-(\not{p}'-m)\gamma^\nu]=2p'^\nu\bar u(p')$、分母は $2p'\cdot k$ です（$k^2=0$）。よって因子は $(-ie)\,i\,\frac{2p'\cdot\varepsilon^*}{2p'\cdot k}=e\frac{p'\cdot\varepsilon^*}{p'\cdot k}$ です。入射電子の線からの放射は分母が $(p-k)^2-m^2=-2p\cdot k$ となるので符号が逆です。内線から出る図は $k\to0$ で特異にならないので主要項に入りません。偏極和は、$J^\mu=\frac{p'^\mu}{p'\cdot k}-\frac{p^\mu}{p\cdot k}$ が $k_\mu J^\mu=0$ を満たすので <Ref to="physics/qft/maxwell-field-quantization#cor-polarization-sum-replacement" /> により $-J\cdot J$ に置き換えられ、展開すれば主張の式です。光子の位相空間 $\frac{d^3k}{(2\pi)^32|\boldsymbol k|}$ は <Ref to="physics/qft/feynman-rules-lsz#thm-cross-section" /> の $d\Pi$ から、運動量保存の変化が $O(k)$ であることを使って分離しました。
</Proof>

括弧の中は $|\boldsymbol k|^{-2}$ に比例するので、$\frac{d^3k}{|\boldsymbol k|}$ と合わせて $\frac{d|\boldsymbol k|}{|\boldsymbol k|}$ の対数発散が生じます。光子質量 $\mu$ から検出器の分解能 $E_\ell$ までの範囲で積分すると、$\boldsymbol k$ の向きについての平均 $\int\frac{d\Omega}{4\pi}$ が $2f_{\mathrm{IR}}(q^2)$ を与え（Peskin–Schroeder 第 6 章 §6.1・§6.5）、
$$
d\sigma(p\to p'+\gamma,\ |\boldsymbol k|\lt E_\ell)=d\sigma_0\cdot\frac\alpha\pi f_{\mathrm{IR}}(q^2)\ln\frac{E_\ell^2}{\mu^2}
$$
です。

<Theorem id="thm-bloch-nordsieck" title="Bloch–Nordsieck の相殺">
電子が運動量移行 $q$ で散乱し、全エネルギーが $E_\ell$ 未満の光子をいくつ放射してもよい（検出器が区別しない）とする包含的な断面積は、$\alpha$ の 1 次で
$$
\frac{d\sigma}{d\Omega}\Big(p\to p'+(\text{軟光子 }E\lt E_\ell)\Big)=\Big(\frac{d\sigma}{d\Omega}\Big)_0\Big[1-\frac\alpha\pi f_{\mathrm{IR}}(q^2)\ln\frac{m^2}{E_\ell^2}+(\mu\text{ によらない項})+O(\alpha^2)\Big]
$$
となり、光子質量 $\mu\to0$ の極限が存在します。
</Theorem>

<Proof of="thm-bloch-nordsieck">
散乱なので $q^2\lt0$ とし、軟光子近似が使えるように $E_\ell$ は $m$ と $\sqrt{-q^2}$ に比べて小さいとします（$E_\ell$ に比例して消える補正は「$\mu$ によらない項」に含めます）。$\alpha$ の 1 次では、包含的な断面積は光子を出さない過程の 1 ループまでの断面積と、光子 1 個の放射の木の断面積の和です（光子 2 個以上の放射は $O(\alpha^2)$）。

**仮想補正。** 光子を出さない過程の振幅は、最低次の振幅の $\gamma^\mu$ を $\Gamma^\mu$ に替え、真空偏極の補正を加えたものです。$|\mathcal M|^2$ の $\alpha$ の 1 次の項のうち、$F_2$ の項は <Ref to="thm-schwinger" /> により赤外有限、真空偏極（<Ref to="thm-vacuum-polarization-one-loop" />）は電子のループで光子質量によりません。残る $F_1$ の項は $|F_1|^2=1+2\,\mathrm{Re}\,\delta F_1+O(\alpha^2)$ を与え、<Ref to="prop-ir-divergence-f1" /> から
$$
\Big(\frac{d\sigma}{d\Omega}\Big)_{\text{仮想}}=\Big(\frac{d\sigma}{d\Omega}\Big)_0\Big[1-\frac\alpha\pi f_{\mathrm{IR}}(q^2)\ln\frac{m^2}{\mu^2}+(\mu\text{ によらない項})\Big]
$$
です。

**角度平均。** 時間的な $P$（$P^0\gt|\boldsymbol P|$）について、$u=\cos\theta$ として
$$
\int\frac{d\Omega}{4\pi}\frac{1}{(P^0-|\boldsymbol P|u)^2}=\frac12\int_{-1}^1\frac{du}{(P^0-|\boldsymbol P|u)^2}=\frac{1}{2|\boldsymbol P|}\Big[\frac1{P^0-|\boldsymbol P|}-\frac1{P^0+|\boldsymbol P|}\Big]=\frac1{P^2}
$$
です。$k=|\boldsymbol k|(1,\hat{\boldsymbol n})$ と書くと $p\cdot k=|\boldsymbol k|(p^0-\boldsymbol p\cdot\hat{\boldsymbol n})$ です。<Ref to="lem-feynman-parameters" /> の $\frac1{AB}=\int_0^1\frac{d\xi}{(\xi A+(1-\xi)B)^2}$ を使い、$P_\xi=\xi p+(1-\xi)p'$（時間的なベクトルの凸結合なので時間的）と置くと
$$
|\boldsymbol k|^2\int\frac{d\Omega}{4\pi}\frac{1}{(p\cdot k)(p'\cdot k)}=\int_0^1\frac{d\xi}{P_\xi^2},\qquad P_\xi^2=m^2-\xi(1-\xi)q^2=c(\xi)
$$
です（$p\cdot p'=m^2-\frac{q^2}2$ から）。同様に $|\boldsymbol k|^2\int\frac{d\Omega}{4\pi}\frac{m^2}{(p\cdot k)^2}=\frac{m^2}{p^2}=1$ です。したがって <Ref to="prop-soft-photon-emission" /> の括弧の角度平均は
$$
\frac1{|\boldsymbol k|^2}\Big[2\Big(m^2-\frac{q^2}2\Big)\int_0^1\frac{d\xi}{c(\xi)}-2\Big]=\frac{2f_{\mathrm{IR}}(q^2)}{|\boldsymbol k|^2}
$$
です。

**実放射。** <Ref to="prop-soft-photon-emission" /> を光子の運動量について積分します。光子質量 $\mu$ の効果は、$|\boldsymbol k|\lesssim\mu$ で被積分関数の $|\boldsymbol k|^{-3}$ の特異性を均すことで、対数の係数を変えず、定数の違いだけを生みます。よって
$$
\int_{\mu\lt|\boldsymbol k|\lt E_\ell}\frac{d^3k}{(2\pi)^32|\boldsymbol k|}\,e^2\frac{2f_{\mathrm{IR}}}{|\boldsymbol k|^2}=\frac{4\pi e^2\cdot2f_{\mathrm{IR}}}{2(2\pi)^3}\int_\mu^{E_\ell}\frac{d|\boldsymbol k|}{|\boldsymbol k|}=\frac{\alpha}{\pi}f_{\mathrm{IR}}(q^2)\ln\frac{E_\ell^2}{\mu^2}
$$
（$e^2=4\pi\alpha$）で、
$$
\Big(\frac{d\sigma}{d\Omega}\Big)_{\text{実}}=\Big(\frac{d\sigma}{d\Omega}\Big)_0\Big[\frac\alpha\pi f_{\mathrm{IR}}(q^2)\ln\frac{E_\ell^2}{\mu^2}+(\mu\text{ によらない項})\Big]
$$
です。

**和。** $-\ln\frac{m^2}{\mu^2}+\ln\frac{E_\ell^2}{\mu^2}=-\ln\frac{m^2}{E_\ell^2}$ なので $\ln\mu^2$ が打ち消し合い、主張の式を得ます。残りの項は $\mu\to0$ で極限をもつので、包含的な断面積の $\mu\to0$ の極限が存在します。
</Proof>

<Remark id="rem-bloch-nordsieck-proof" title="相殺の仕組みと全次数への拡張">
$\alpha$ の 1 次では、<Ref to="prop-ir-divergence-f1" /> から来る仮想補正 $2\,\mathrm{Re}\,\delta F_1$ の $-\frac\alpha\pi f_{\mathrm{IR}}\ln\frac{m^2}{\mu^2}$ と、<Ref to="prop-soft-photon-emission" /> を積分した実放射の $+\frac\alpha\pi f_{\mathrm{IR}}\ln\frac{E_\ell^2}{\mu^2}$ を足すと $\ln\mu^2$ が打ち消し合います。$f_{\mathrm{IR}}$ が両方に共通に現れるのは偶然ではなく、どちらも外線の古典的な電流 $J^\mu$ が作る放射場で決まっているからです。全次数では、$n$ 個の軟光子の放射確率が $\frac{1}{n!}\big(\frac\alpha\pi f_{\mathrm{IR}}\ln\frac{E_\ell^2}{\mu^2}\big)^n$ となって指数関数に足し上がり、仮想光子の補正も $\exp\big(-\frac\alpha\pi f_{\mathrm{IR}}\ln\frac{m^2}{\mu^2}\big)$ に指数化して、両者の積で $\mu$ が消えます（Bloch–Nordsieck 1937、Yennie–Frautschi–Suura 1961。Peskin–Schroeder 第 6 章 §6.5）。理論的な背景は、<Ref to="physics/qft/feynman-rules-lsz#rem-lsz-derivation" /> で述べた「電子の 1 粒子極が軟光子の連続スペクトルに埋もれている」ことにあり、観測可能なのは縮退した状態をまとめた包含的な量だけです。
</Remark>

<Example id="ex-exclusive-vanishes" title="光子を出さない散乱の確率は 0 である">
<Ref to="rem-bloch-nordsieck-proof" /> の指数化を認めると、光子を 1 個も放射しない排他的な断面積は
$$
\frac{d\sigma}{d\Omega}(p\to p'\text{ のみ})=\Big(\frac{d\sigma}{d\Omega}\Big)_0\exp\Big(-\frac\alpha\pi f_{\mathrm{IR}}(q^2)\ln\frac{m^2}{\mu^2}\Big)\times(1+\cdots)
$$
で、$q^2\ne0$ なら $f_{\mathrm{IR}}\gt0$ なので $\mu\to0$ で $0$ になります。1 次の摂動論で見えた「負の無限大の補正」は、この指数関数の展開の第 1 項です。加速された電荷は必ず無限個の軟光子を放射するので、「光子を出さずに散乱する」という事象には確率がありません。包含的な断面積をとるという仮定（検出器の分解能 $E_\ell\gt0$）を外すと、観測量が定義できなくなるのです。古典電磁気学で、加速された電荷の放射スペクトル $\frac{dI}{d\omega}$ が $\omega\to0$ で有限な値に近づき、光子数 $\frac{1}{\omega}\frac{dI}{d\omega}$ が対数発散することに対応しています。
</Example>

## 8. 演習

<Exercise id="exr-z2-divergence" difficulty="標準">
ファインマンゲージで、1 ループの電子の自己エネルギーが
$$
\Sigma_2(p)=\frac{e^2}{(4\pi)^{d/2}}\int_0^1dx\,\big[-(d-2)(1-x)\not{p}+dm\big]\frac{\Gamma(2-\frac d2)}{\Delta^{2-d/2}},\qquad \Delta=xm^2-x(1-x)p^2
$$
であることを示し、その紫外発散の部分から $Z_2=1-\frac{\alpha}{4\pi}\frac2\varepsilon+(\text{有限})$ を求めてください。さらに <Ref to="thm-schwinger" /> の証明で $\gamma^\mu$ に比例する $\ell^2$ の項から $Z_1$ の極を求め、<Ref to="cor-z1-equals-z2" /> と整合することを確かめてください。

<Solution>
<Ref to="physics/qft/feynman-rules-lsz#thm-yukawa-qed-rules" /> より
$$
-i\Sigma_2(p)=\int\frac{d^dk}{(2\pi)^d}(-ie\gamma^\mu)\frac{i(\not{k}+m)}{k^2-m^2}(-ie\gamma_\mu)\frac{-i}{(p-k)^2}=-e^2\int\frac{d^dk}{(2\pi)^d}\frac{-(d-2)\not{k}+dm}{(k^2-m^2)(p-k)^2}
$$
です（<Ref to="lem-symmetric-integration" /> の $\gamma^\mu\gamma^\nu\gamma_\mu=-(d-2)\gamma^\nu$ と $\gamma^\mu\gamma_\mu=d$）。$x(k^2-m^2)+(1-x)(p-k)^2=(k-(1-x)p)^2-\Delta$ なので $\ell=k-(1-x)p$ として $\not{k}\mapsto(1-x)\not{p}$ です。<Ref to="thm-master-integrals" /> の $n=2$ を使って $i$ で割ると主張の式です。

$d=4-\varepsilon$ の極は $\frac{e^2}{16\pi^2}\frac2\varepsilon\int_0^1[-2(1-x)\not{p}+4m]dx=\frac\alpha{4\pi}\frac2\varepsilon(-\not{p}+4m)$ なので $\frac{d\Sigma_2}{d\not{p}}$ の極は $-\frac\alpha{4\pi}\frac2\varepsilon$ で、$Z_2^{-1}=1+\frac\alpha{4\pi}\frac2\varepsilon$、$Z_2=1-\frac\alpha{4\pi}\frac2\varepsilon$ です。

頂点の $\ell^2$ の項は $2ie^2\int\frac{d^d\ell}{(2\pi)^d}\int dF\frac{2}{(\ell^2-\Delta)^3}\big(-\frac12\ell^2\big)\gamma^\mu$ で、<Ref to="thm-master-integrals" /> の第 2 式（$n=3$）の極は $\frac{i}{16\pi^2}\frac2\varepsilon$、$\int dx\,dy\,dz\,\delta(x+y+z-1)=\frac12$ なので、$\delta F_1$ の極は $2ie^2\cdot2\cdot(-\frac12)\cdot\frac12\cdot\frac{i}{16\pi^2}\frac2\varepsilon=\frac{\alpha}{4\pi}\frac2\varepsilon$ です。$Z_1^{-1}=F_1(0)$（引き算前）$=1+\frac\alpha{4\pi}\frac2\varepsilon$ なので $Z_1=1-\frac\alpha{4\pi}\frac2\varepsilon$ となり、極は $Z_2$ と一致します。
</Solution>
</Exercise>

<Exercise id="exr-uehling" difficulty="標準">
$|q^2|\ll m^2$ で $\hat\Pi_2(q^2)\approx-\frac{\alpha}{15\pi}\frac{q^2}{m^2}$ を示してください。これを使い、静的な Coulomb ポテンシャル $V(r)=-\frac\alpha r$（電荷が逆符号の 2 粒子）が真空偏極によって
$$
V(\boldsymbol x)=-\frac{\alpha}{r}-\frac{4\alpha^2}{15m^2}\delta^3(\boldsymbol x)
$$
と補正されることを示してください（Uehling 項。水素原子の $2S$ 準位を下げ、Lamb シフトの一部をなします）。

<Solution>
$\ln\frac{m^2}{m^2-x(1-x)q^2}=-\ln\big(1-\frac{x(1-x)q^2}{m^2}\big)\approx\frac{x(1-x)q^2}{m^2}$ なので、$\int_0^1x^2(1-x)^2dx=B(3,3)=\frac{1}{30}$ より $\hat\Pi_2\approx-\frac{2\alpha}{\pi}\cdot\frac{1}{30}\frac{q^2}{m^2}=-\frac{\alpha}{15\pi}\frac{q^2}{m^2}$ です。

静的な散乱では $q=(0,\boldsymbol q)$、$q^2=-|\boldsymbol q|^2$ です。Born 近似で $V(r)=-\frac\alpha r$ に対応する運動量空間のポテンシャルは $\tilde V(\boldsymbol q)=-\frac{4\pi\alpha}{|\boldsymbol q|^2}$ で、<Ref to="thm-running-alpha" /> により $\alpha$ を $\alpha_{\mathrm{eff}}(-|\boldsymbol q|^2)$ に置き換えます。$\alpha_{\mathrm{eff}}\approx\alpha(1+\hat\Pi_2)=\alpha\big(1+\frac{\alpha}{15\pi}\frac{|\boldsymbol q|^2}{m^2}\big)$ なので
$$
\tilde V(\boldsymbol q)=-\frac{4\pi\alpha}{|\boldsymbol q|^2}-\frac{4\alpha^2}{15m^2}
$$
です。定数の Fourier 逆変換は $\delta^3(\boldsymbol x)$ なので主張を得ます。符号が負なので、原点で波動関数が $0$ でない $S$ 状態だけが引力的に下がります。
</Solution>
</Exercise>

<Exercise id="exr-alpha-at-mz" difficulty="標準">
<Ref to="thm-running-alpha" /> を電子・ミューオン・タウの 3 種の荷電レプトンに足し合わせ（<Ref to="thm-qed-beta-one-loop" /> の $\sum Q_i^2$ と同じ理由で寄与は加算されます）、$\alpha=1/137.036$ から空間的な運動量移行 $Q=m_Z=91.19\ \mathrm{GeV}$ での $\alpha_{\mathrm{eff}}^{-1}$ を求めてください（$m_e=0.000511$、$m_\mu=0.1057$、$m_\tau=1.777$、単位 GeV）。実測の $\alpha^{-1}(m_Z)\approx128.9$ との差は何に由来しますか。

<Solution>
1 ループで $\alpha_{\mathrm{eff}}^{-1}=\alpha^{-1}-\frac{1}{3\pi}\sum_\ell\Big(\ln\frac{Q^2}{m_\ell^2}-\frac53\Big)$ です。$\ln\frac{m_Z}{m_e}=12.09$、$\ln\frac{m_Z}{m_\mu}=6.76$、$\ln\frac{m_Z}{m_\tau}=3.94$ なので $\sum\ln\frac{Q^2}{m_\ell^2}=2\times22.79=45.58$ で、$\frac53\times3=5$ を引いて $40.58$ です。$\frac{40.58}{3\pi}=4.31$ なので $\alpha_{\mathrm{eff}}^{-1}(m_Z)\approx137.04-4.31=132.7$ です。

残りの約 $3.8$ はクォークのループ（ハドロンの真空偏極）の寄与です。軽いクォークでは質量 $m$ がよく定義されず摂動論も使えないので、$e^+e^-\to$ ハドロンの断面積の実測値から分散関係で求めます。
</Solution>
</Exercise>

<Exercise id="exr-f2-spacelike" difficulty="難">
<Ref to="thm-schwinger" /> の $F_2(q^2)$ の式で $z$ を固定して $x,y$ について積分する前に、$-q^2=Q^2\gg m^2$ での振る舞いを調べ、主要項が $F_2(q^2)\approx\frac{\alpha}{\pi}\frac{m^2}{Q^2}\ln\frac{Q^2}{m^2}$ であることを示してください。この結果は、高エネルギーでの散乱で電子の異常磁気モーメントの効果が見えにくくなることを意味します。

<Solution>
$y=(1-z)\eta$、$x=(1-z)(1-\eta)$（$0\le\eta\le1$）と置くと $dx\,dy\,\delta(\cdots)\mapsto(1-z)d\eta$、$xy=(1-z)^2\eta(1-\eta)$ で
$$
F_2=\frac\alpha{2\pi}\int_0^1dz\int_0^1d\eta\,\frac{2m^2z(1-z)^2}{(1-z)^2\big[m^2+Q^2\eta(1-\eta)\big]}=\frac\alpha{2\pi}\int_0^1dz\,2z\int_0^1\frac{m^2\,d\eta}{m^2+Q^2\eta(1-\eta)}
$$
です。$z$ 積分は $1$ です。$\eta$ 積分は、$Q^2\gg m^2$ では $\eta\approx0$ と $\eta\approx1$ の近くの幅 $m^2/Q^2$ の領域が支配し、$2\int_0^{1/2}\frac{m^2d\eta}{m^2+Q^2\eta}\approx\frac{2m^2}{Q^2}\ln\frac{Q^2}{m^2}$ です。よって $F_2\approx\frac{\alpha}{2\pi}\frac{2m^2}{Q^2}\ln\frac{Q^2}{m^2}=\frac{\alpha}{\pi}\frac{m^2}{Q^2}\ln\frac{Q^2}{m^2}$ で、$Q^2\to\infty$ で $0$ に近づきます。$q^2=0$ では $\eta$ 積分が $1$ で、<Ref to="thm-schwinger" /> の $F_2(0)=\frac{\alpha}{2\pi}$ に戻ります。
</Solution>
</Exercise>

## 参考文献

- M. E. Peskin, D. V. Schroeder『An Introduction to Quantum Field Theory』Addison-Wesley、1995 — 第 6 章（§6.2 形状因子と $g$ 因子、§6.3 頂点補正、§6.4・§6.5 赤外発散と軟光子）、第 7 章（§7.1 電子の自己エネルギー、§7.4 Ward–高橋恒等式、§7.5 真空偏極と電荷のくりこみ）、第 10 章 §10.3。
- J. Schwinger, "On quantum-electrodynamics and the magnetic moment of the electron", *Physical Review* 73 (1948), 416.
- F. Bloch, A. Nordsieck, "Note on the radiation field of the electron", *Physical Review* 52 (1937), 54.
- J. C. Ward, "An identity in quantum electrodynamics", *Physical Review* 78 (1950), 182. Y. Takahashi, "On the generalized Ward identity", *Il Nuovo Cimento* 6 (1957), 371.
- G. 't Hooft, M. Veltman, "Regularization and renormalization of gauge fields", *Nuclear Physics B* 44 (1972), 189.
- J. C. Collins『Renormalization』Cambridge University Press、1984 — 第 4 章（次元正則化）。
