2026-10-01 · chapter ヤン–ミルズと質量ギャップ 第 1 部 数学の準備 線形代数
この本の全体 目次と読む順
第 0 部 入口 — この本の読み方
0-01 この本の読み方
0-02 一枚の絵
0-03 問題文を一語ずつ読む
0-04 数学の四次元と物理の四次元
第 1 部 数学の準備
1-01 ベクトル空間と線形写像
1-02 群とは何か
1-03 リー群とリー環
1-04 SU(2) と SU(3)
1-05 多様体と接空間
1-06 微分形式と外微分
1-07 ベクトル束と接続
1-08 確率と測度
1-09 無限次元の確率
1-10 ヒルベルト空間と自己共役作用素
1-11 フーリエ解析と分布
1-12 寄り道
第 2 部 物理の準備
2-01 ラグランジアンと作用
2-02 場という考え
2-03 電磁気学はゲージ理論である
2-04 特殊相対論と時空
2-05 量子力学の骨
2-06 調和振動子と生成消滅
2-07 経路積分の考え方
2-08 統計力学と相転移
2-09 寄り道
2-10 緩和の時間と動的指数 z
第 3 部 ヤン–ミルズ理論(古典)
3-01 ゲージ原理
3-02 非可換ゲージ場
3-03 作用と方程式
3-04 幾何としてのゲージ理論
3-05 インスタントンと位相
3-06 寄り道
3-07 標準模型の中のヤン–ミルズ
第 4 部 量子化
4-01 正準量子化とハミルトニアン
4-02 経路積分とユークリッド化
4-03 摂動論と Feynman ダイアグラム
4-04 発散と繰り込み
4-05 発散の代数
4-06 漸近自由
4-07 次元転移と Λ
4-08 ゲージ固定と Faddeev–Popov
4-09 場の量子論の公理
4-10 Osterwalder–Schrader の公理と再構成
4-11 質量ギャップの定義
4-12 寄り道
第 5 部 格子ゲージ理論
5-01 Wilson の格子
5-02 強結合展開
5-03 反射正値性と転送行列
5-04 無限体積極限とクラスター展開
5-05 弱結合と連続極限
5-06 U(1) と非可換の違い
5-07 モンテカルロ法
5-08 グルーボールと弦張力の測定
5-09 何を固定して極限を取るか
5-10 有限群の格子ゲージ理論
5-11 寄り道
第 6 部 構成的場の理論
6-01 構成的場の理論とは
6-02 二次元の可解性とヤン–ミルズ測度
6-03 スカラー場の構成
6-04 クラスター展開
6-05 繰り込み群の段の列
6-06 三次元ヤン–ミルズの紫外安定性
6-07 四次元
6-08 四次元の φ⁴ の自明性
6-09 確率量子化と正則性構造
6-10 四次元で止まる場所
6-11 発散以外の障害
6-12 寄り道
第 7 部 物理の側から
7-01 物理はどう見ているか
7-02 閉じ込めの機構
7-03 弦の絵
7-04 大 N
7-05 ひも理論と余剰次元
7-06 余剰次元が見えなくなる仕組み
7-07 ゲージ場はどこから来るか
7-08 ホログラフィー
7-09 質量ギャップが幾何になる
7-10 四次元に戻す
7-11 超対称と Seiberg–Witten
7-12 等価原理に当たる一文
7-13 物理の掘り方が数学と離れる場所
第 8 部 二つの言葉の辞書 — 物理の視点と数学の視点
8-01 辞書の読み方
8-02 辞書 A
8-03 辞書 B
8-04 辞書 C
8-05 私たちの仮定の物理側の対応
第 9 部 現在地と課題
9-01 世界はどこまで来ているか
9-02 二つの掘り方の切れ目
9-03 新しい概念の候補
9-04 課題の一覧
9-05 よくある誤解
第 10 部 質量ギャップの厳密な証明へ — この端末の検討
10-00 第 10 部の入口 — 酔歩と定規と時計
10-01 理論の構成の筋
10-02 一段の記帳
10-03 仮定 H と三つの鎖
10-04 方向の地図
10-05 方向 12〜14
10-06 方向 15・15′
10-07 Lean で閉じた言明と既存の結果の対応表
10-08 壁の一覧
10-09 ひらめき帳から
10-10 主張しないこと
10-12 つじつま合わせ
10-11 定理までの距離
10-13 小さな問い — 卒業研究の大きさで決着のつく十〜二十問
付録
A-01 記号表
A-02 用語集
A-03 文献案内
A-04 Lean と機械検査
A-05 この本の作り方
A-06 仮定の索引
ベクトル空間と線形写像 — 行列を「ベクトルを動かす規則」として読み、固有値・指数関数・エルミートとユニタリを揃える
この章で分かること — 行列を「ベクトルを動かす規則」として読むこと。固有値が「向きを変えない方向」の伸び率であること。エルミート行列(固有値が実数)とユニタリ行列(長さを保つ)が U = e i H U = e^{iH} U = e i H で結ばれること。
前提となる章 — なし。行列の積と行列式、複素数、定数係数の常微分方程式を既知とします。
先に言うこと — 有限次元に限ります。質量ギャップは無限次元の作用素の問題で、有限次元の結論がそのまま延びるとは限りません(1-10 )。定理はすべて既知の線形代数です。
01
なぜ線形代数から始めるか
ヤン–ミルズ理論では、数が入るはずの場所に行列が入ります。役は三つです。
変換 ゲージ変換
g ( x ) g(x) g ( x ) はユニタリ行列
各点で長さを保つ行列を掛ける(§07)
場 ゲージ場
A μ ( x ) A_\mu(x) A μ ( x ) の値は行列
μ \mu μ は方向(時間・空間)の番号。掛ける順で値が変わる=非可換(§03・§09)
エネルギー ハミルトニアン
H H H はエルミート
質量ギャップは最低の値と次の値の隙間に当たる(§06)。正確には、真空(最低の状態)を除いたスペクトルの下端(1-10 ・4-11 )
ゲージ変換に使う行列の大きさを N N N (N × N N \times N N × N 行列)とすると、電磁気学は N = 1 N = 1 N = 1 の場合で、ゲージ変換は絶対値 1 の数 e i α e^{i\alpha} e i α を掛けるだけです(2-03 )。数の積は順序に依らないので、物質の無いとき、ゲージ場自身の方程式(Maxwell 方程式)は線形のままです。N ≥ 2 N \ge 2 N ≥ 2 では順序が効き、物質が無くてもゲージ場どうしの積が方程式に現れます(3-02 )。この章はその土台を有限次元で揃えます。
02
ベクトル空間
ベクトル空間とは、足し算とスカラー倍ができ、結合・交換・分配の法則など数と同じ規則が成り立つ集まりです。スカラーを R \mathbb{R} R から取れば実ベクトル空間、C \mathbb{C} C から取れば複素ベクトル空間です。どのベクトルも v = c 1 v 1 + ⋯ + c n v n v = c_1v_1 + \cdots + c_nv_n v = c 1 v 1 + ⋯ + c n v n とただ一通りに書けるとき v 1 , … , v n v_1, \dots, v_n v 1 , … , v n を基底 、本数 n n n を次元 と呼びます。例は R n \mathbb{R}^n R n 、C n \mathbb{C}^n C n 、n × n n \times n n × n 行列全体(次元 n 2 n^2 n 2 )、区間上の連続関数全体(無限次元、1-10 )です。
この本で何度も戻る例は、トレース 0 の 2×2 エルミート行列の全体です。エルミートとは H † = H H^\dagger = H H † = H となることで、† \dagger † は転置して各成分の複素共役を取る操作です(§05・§06)。和も実数倍もエルミートですが、i i i 倍すると ( i H ) † = − i H (iH)^\dagger = -iH ( i H ) † = − i H となって外れます。つまり成分は複素数なのに実ベクトル空間 で、次元は 3 です。
( x y − i z y + i z − x ) = x ( 1 0 0 − 1 ) + y ( 0 1 1 0 ) + z ( 0 − i i 0 ) , x , y , z ∈ R \begin{pmatrix} x & y - iz \\ y + iz & -x \end{pmatrix} = x\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix} + y\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} + z\begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix},\qquad x, y, z \in \mathbb{R} ( x y + i z y − i z − x ) = x ( 1 0 0 − 1 ) + y ( 0 1 1 0 ) + z ( 0 i − i 0 ) , x , y , z ∈ R
右辺の三つは 1-02 §05 ・1-04 で Pauli 行列 σ 3 , σ 1 , σ 2 \sigma_3, \sigma_1, \sigma_2 σ 3 , σ 1 , σ 2 として再登場し、S U ( 2 ) \mathrm{SU}(2) SU ( 2 ) のゲージ場の値はこの空間(の i i i 倍)に入ります。
03
線形写像と行列
線形写像は和とスカラー倍を保つ写像です(L ( u + v ) = L u + L v L(u+v) = Lu + Lv L ( u + v ) = Lu + Lv 、L ( c v ) = c L v L(cv) = c\,Lv L ( c v ) = c Lv )。基底を決めると行列で書け、j j j 列目が基底 e j e_j e j の行き先 です。合成は行列の積 ( A B ) v = A ( B v ) (AB)v = A(Bv) ( A B ) v = A ( B v ) で、数と違い、順序を入れ替えると一般に変わります。
A = ( 0 1 0 0 ) , B = ( 0 0 1 0 ) : A B = ( 1 0 0 0 ) , B A = ( 0 0 0 1 ) A = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix},\ B = \begin{pmatrix} 0 & 0 \\ 1 & 0 \end{pmatrix}:\qquad AB = \begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix},\quad BA = \begin{pmatrix} 0 & 0 \\ 0 & 1 \end{pmatrix} A = ( 0 0 1 0 ) , B = ( 0 1 0 0 ) : A B = ( 1 0 0 0 ) , B A = ( 0 0 0 1 )
差 [ A , B ] = A B − B A [A, B] = AB - BA [ A , B ] = A B − B A を交換子 と呼びます。この例では [ A , B ] = diag ( 1 , − 1 ) [A, B] = \operatorname{diag}(1, -1) [ A , B ] = diag ( 1 , − 1 ) で、§02 の基底の一つ目です。交換子がリー環の骨格になる話は 1-03 で扱います。
行列式 det A \det A det A は面積(体積)の符号つきの倍率、トレース tr A \operatorname{tr} A tr A は対角成分の和です。det ( A B ) = det A det B \det(AB) = \det A \det B det ( A B ) = det A det B 、tr ( A B ) = tr ( B A ) \operatorname{tr}(AB) = \operatorname{tr}(BA) tr ( A B ) = tr ( B A ) から、どちらも基底の取り替え A ↦ P − 1 A P A \mapsto P^{-1}AP A ↦ P − 1 A P で変わらず、写像そのものの性質です。例:回転 R ( θ ) = ( cos θ − sin θ sin θ cos θ ) R(\theta) = \left(\begin{smallmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{smallmatrix}\right) R ( θ ) = ( c o s θ s i n θ − s i n θ c o s θ ) は det R = 1 \det R = 1 det R = 1 、tr R = 2 cos θ \operatorname{tr} R = 2\cos\theta tr R = 2 cos θ 。
04
固有値と固有ベクトル
写像はたいていのベクトルの向きを変えますが、向きを変えずに伸び縮みだけさせる方向 をもつことがあります。その方向のベクトルが固有ベクトル、伸び率が固有値です。
A v = λ v ( v ≠ 0 ) ⟺ det ( A − λ I ) = 0 Av = \lambda v \ \ (v \ne 0) \qquad\Longleftrightarrow\qquad \det(A - \lambda I) = 0 A v = λ v ( v = 0 ) ⟺ det ( A − λ I ) = 0
2×2 では λ 2 − ( tr A ) λ + det A = 0 \lambda^2 - (\operatorname{tr} A)\lambda + \det A = 0 λ 2 − ( tr A ) λ + det A = 0 。実の 2×2 行列では、判別式 D = ( tr A ) 2 − 4 det A D = (\operatorname{tr} A)^2 - 4\det A D = ( tr A ) 2 − 4 det A が正なら実の固有方向が二本、負なら無く(固有値は複素数の組)、0 なら重根です(重根のとき固有方向は一本だけか、スカラー行列 λ I \lambda I λ I ならすべての方向)。n n n 本の独立な固有ベクトルを並べた P P P で、P − 1 A P P^{-1}AP P − 1 A P は対角行列になります(対角化)。
v v v と A v Av A v が同じ直線に乗ることは、二つを並べた行列式 det [ v , A v ] \det[v, Av] det [ v , A v ] が 0 になることと同じです。図1 で確かめられます。
図1 つまみで単位ベクトル
v v v を回すと
A v Av A v (青)が動く。同じ直線に乗る角度が固有方向で、下の曲線
det [ v , A v ] \det[v, Av] det [ v , A v ] がそこで 0 になる。点線は単位円の像
図は JavaScript で描画します。四つの行列の固有値と実の固有方向は下の表にあります(対称 (2 1 ; 1 2) では 45° と 135° で、直交)。
行列 対称 (2 1 ; 1 2) 対称でない (1.5 1 ; 0.5 1) ずれ (1 1 ; 0 1) 60° 回転 角度 動かす
JavaScript が無効なら下の表で読めます。
計算 この図の数値はこの端末で計算しました(表は NumPy の固有値計算)。行列は行を「;」で区切っています。下の表の行を押すと、その行列が図に入ります。
行列 D D D 固有値 実の固有方向(x x x 軸から)
対称 (2 1 ; 1 2) 4 3, 1 45°, 135°(直交)
対称でない (1.5 1 ; 0.5 1) 2.25 2, 0.5 26.565°, 135°
ずれ (1 1 ; 0 1) 0 1(重根) 0° の一本だけ
60° 回転 −3 0.5 ± 0.866025 i 0.5 \pm 0.866025\,i 0.5 ± 0.866025 i なし
対称行列では固有方向が直交します(理由は §06)。ずれは固有方向が一本しかなく、対角化できません。回転には実の固有方向が無く、det [ v , A v ] \det[v, Av] det [ v , A v ] はどの向きでも sin 60 ∘ ≈ 0.866 \sin 60^\circ \approx 0.866 sin 6 0 ∘ ≈ 0.866 です。
寄り道:回転に実の固有ベクトルが無いこと
寄り道
飛ばしてよい小節です。
R ( θ ) R(\theta) R ( θ ) の固有値は e ± i θ e^{\pm i\theta} e ± i θ 、固有ベクトルは ( 1 , ∓ i ) (1, \mp i) ( 1 , ∓ i ) です。比喩で言えば、平面の回転は複素数の世界では二つの「位相の掛け算」に分かれます。正確には、C 2 \mathbb{C}^2 C 2 の基底を ( 1 , − i ) (1, -i) ( 1 , − i ) 、( 1 , i ) (1, i) ( 1 , i ) に取り替えると R ( θ ) = diag ( e i θ , e − i θ ) R(\theta) = \operatorname{diag}(e^{i\theta}, e^{-i\theta}) R ( θ ) = diag ( e i θ , e − i θ ) です。複素数では固有方程式が必ず根をもつ(代数学の基本定理)ので、どの正方行列にも固有値があります。四次元の回転の二つの角(四次元の記事 )は、固有値 e ± i α e^{\pm i\alpha} e ± i α 、e ± i β e^{\pm i\beta} e ± i β の二組です。
05
内積と随伴
長さと角度を測るのが内積です。複素数の成分でそのまま ∑ k u k v k \sum_k u_kv_k ∑ k u k v k とすると、v = ( 1 , i ) v = (1, i) v = ( 1 , i ) の「長さの二乗」が 1 + i 2 = 0 1 + i^2 = 0 1 + i 2 = 0 になってしまいます。そこで片方に複素共役を付けます(z ‾ z = ∣ z ∣ 2 ≥ 0 \overline{z}\,z = |z|^2 \ge 0 z z = ∣ z ∣ 2 ≥ 0 )。
⟨ u , v ⟩ = ∑ k u k ‾ v k , ∥ v ∥ = ⟨ v , v ⟩ , ⟨ u , v ⟩ = ⟨ v , u ⟩ ‾ \langle u, v\rangle = \sum_k \overline{u_k}\,v_k,\qquad \lVert v\rVert = \sqrt{\langle v, v\rangle},\qquad \langle u, v\rangle = \overline{\langle v, u\rangle} ⟨ u , v ⟩ = k ∑ u k v k , ∥ v ∥ = ⟨ v , v ⟩ , ⟨ u , v ⟩ = ⟨ v , u ⟩
この本は物理の流儀で、第 2 引数について線形とします(⟨ c u , v ⟩ = c ‾ ⟨ u , v ⟩ \langle cu, v\rangle = \overline{c}\langle u, v\rangle ⟨ c u , v ⟩ = c ⟨ u , v ⟩ )。数学の教科書には逆の流儀も多く(例:Axler の定義 6.2)、結論は変わりません。⟨ u , v ⟩ = 0 \langle u, v\rangle = 0 ⟨ u , v ⟩ = 0 なら直交、互いに直交する長さ 1 のベクトルの基底が正規直交基底です。
A A A の随伴 A † A^\dagger A † は、転置して成分の複素共役を取った行列です(数学では A ∗ A^* A ∗ )。内積の中で A A A を反対側へ移します。
( A † ) i j = A j i ‾ , ⟨ u , A v ⟩ = ⟨ A † u , v ⟩ , ( A B ) † = B † A † (A^\dagger)_{ij} = \overline{A_{ji}},\qquad \langle u, Av\rangle = \langle A^\dagger u, v\rangle,\qquad (AB)^\dagger = B^\dagger A^\dagger ( A † ) ij = A j i , ⟨ u , A v ⟩ = ⟨ A † u , v ⟩ , ( A B ) † = B † A †
06
エルミート行列
H † = H H^\dagger = H H † = H の行列がエルミート行列 で、実の成分なら対称行列です。量子力学では測定値(エネルギーなど)がエルミートな作用素で表されます(2-05 )。鍵は二つの性質です。
固有値は実数: H v = λ v Hv = \lambda v H v = λ v 、∥ v ∥ = 1 \lVert v\rVert = 1 ∥ v ∥ = 1 なら λ = ⟨ v , H v ⟩ = ⟨ H v , v ⟩ = λ ‾ \lambda = \langle v, Hv\rangle = \langle Hv, v\rangle = \overline{\lambda} λ = ⟨ v , H v ⟩ = ⟨ H v , v ⟩ = λ 。異なる固有値の固有ベクトルは直交: H u = μ u Hu = \mu u H u = μu なら λ ⟨ u , v ⟩ = ⟨ u , H v ⟩ = ⟨ H u , v ⟩ = μ ‾ ⟨ u , v ⟩ = μ ⟨ u , v ⟩ \lambda\langle u, v\rangle = \langle u, Hv\rangle = \langle Hu, v\rangle = \overline{\mu}\langle u, v\rangle = \mu\langle u, v\rangle λ ⟨ u , v ⟩ = ⟨ u , H v ⟩ = ⟨ H u , v ⟩ = μ ⟨ u , v ⟩ = μ ⟨ u , v ⟩ (μ \mu μ は実数)で、μ ≠ λ \mu \ne \lambda μ = λ だから ⟨ u , v ⟩ = 0 \langle u, v\rangle = 0 ⟨ u , v ⟩ = 0 。
押し進めたのがスペクトル定理 です 既知 。エルミート行列には固有ベクトルからなる正規直交基底があり、固有値 λ 1 ≤ ⋯ ≤ λ n \lambda_1 \le \cdots \le \lambda_n λ 1 ≤ ⋯ ≤ λ n は実数です。基底を並べたユニタリ行列 U U U (§07)で U † H U U^\dagger HU U † H U は実対角になります。
最小の固有値と次の固有値の差を隙間 (ギャップ)と呼びます。量子力学では最小の固有値 E 0 E_0 E 0 が基底状態のエネルギー、E 1 − E 0 E_1 - E_0 E 1 − E 0 が最初の励起に要るエネルギーです。有限次元では隙間はいつも計算できます。問題は無限次元で、極限のあとも隙間が正のまま残るかどうかです(1-10 ・4-11 )。
例:§02 の行列は tr H = 0 \operatorname{tr} H = 0 tr H = 0 、det H = − ( x 2 + y 2 + z 2 ) \det H = -(x^2 + y^2 + z^2) det H = − ( x 2 + y 2 + z 2 ) なので固有値は ± x 2 + y 2 + z 2 \pm\sqrt{x^2 + y^2 + z^2} ± x 2 + y 2 + z 2 、隙間が 0 になるのは x = y = z = 0 x = y = z = 0 x = y = z = 0 だけです。図2 は y = c y = c y = c 、z = 0 z = 0 z = 0 の場合です。
図2
H ( x ) = ( x c c − x ) H(x) = \left(\begin{smallmatrix} x & c \\ c & -x \end{smallmatrix}\right) H ( x ) = ( x c c − x ) の二つの固有値。非対角成分
c c c を上げると交差が開いて隙間が残り、
x x x を動かすと下の固有ベクトルの中身(帯、青が
e 1 e_1 e 1 の重み)が入れ替わる
図は JavaScript で描画します。c = 0.3 では二つの固有値の曲線は交わらず、x = 0 で最小の隙間 0.6 が残ります(下の表)。
非対角 c x 動かす
JavaScript が無効なら下の表で読めます。
計算 曲線はその場で ± x 2 + c 2 \pm\sqrt{x^2 + c^2} ± x 2 + c 2 から描き、表は NumPy のエルミート用固有値計算の値です(両者は表示の桁で一致)。x x x を − x -x − x にすると重みが入れ替わります。表の行を押すと、その x x x と c = 0.3 c = 0.3 c = 0.3 が図に入ります。
x x x (c = 0.3 c = 0.3 c = 0.3 )E 0 = − E 1 E_0 = -E_1 E 0 = − E 1 隙間 下の固有ベクトル e 1 : e 2 e_1 : e_2 e 1 : e 2
0.0 −0.300000 0.600000 0.5000 : 0.5000
0.2 −0.360555 0.721110 0.2226 : 0.7774
0.5 −0.583095 1.166190 0.0713 : 0.9287
1.0 −1.044031 2.088061 0.0211 : 0.9789
c = 0 c = 0 c = 0 では直線 ± x \pm x ± x が交わります。c ≠ 0 c \ne 0 c = 0 なら曲線は近づくだけで、最小の隙間 2 c 2c 2 c が残り、そこを過ぎると下の固有ベクトルの中身が入れ替わります。
寄り道:準位は避け合う
寄り道
飛ばしてよい小節です。
2×2 で固有値が重なるには、エルミートなら x = y = z = 0 x = y = z = 0 x = y = z = 0 の三条件、実対称なら二条件が要ります。一つのパラメータで動く族がそこを通るのは特別な場合だけなので、対称性の無い一般の族では、固有値の曲線はふつう交わりません 既知 。対称性があって行列が区画に分かれるときは、違う区画どうしの準位は交わりえます。一般の n × n n \times n n × n での数え上げは von Neumann と Wigner(1929)に遡るとされます。比喩では準位が「反発する」と言います。正確には、力が働くのではなく、重なる行列が条件の数だけ次元の低い薄い集合をなす、ということです。
07
ユニタリ行列
内積を保つ行列がユニタリ行列 です。長さと角度を変えない変換で、実の成分なら直交行列(回転と鏡映)です。
U † U = I ⟺ ⟨ U u , U v ⟩ = ⟨ u , v ⟩ (すべての u , v ) U^\dagger U = I \qquad\Longleftrightarrow\qquad \langle Uu, Uv\rangle = \langle u, v\rangle \ \ \text{(すべての } u, v\text{)} U † U = I ⟺ ⟨ U u , U v ⟩ = ⟨ u , v ⟩ (すべての u , v )
左から右は ⟨ U u , U v ⟩ = ⟨ u , U † U v ⟩ \langle Uu, Uv\rangle = \langle u, U^\dagger Uv\rangle ⟨ U u , U v ⟩ = ⟨ u , U † U v ⟩ から出ます。列が正規直交基底をなすことと同じで、U − 1 = U † U^{-1} = U^\dagger U − 1 = U † 。固有値は ∣ λ ∣ = 1 |\lambda| = 1 ∣ λ ∣ = 1 、∣ det U ∣ = 1 |\det U| = 1 ∣ det U ∣ = 1 です。積と逆行列もユニタリなので、全体 U ( n ) \mathrm{U}(n) U ( n ) と、det U = 1 \det U = 1 det U = 1 の部分 S U ( n ) \mathrm{SU}(n) SU ( n ) は「群」になります(1-02 )。
例:U ( 1 ) = { e i α } \mathrm{U}(1) = \{e^{i\alpha}\} U ( 1 ) = { e i α } は単位円です。∣ a ∣ 2 + ∣ b ∣ 2 = 1 |a|^2 + |b|^2 = 1 ∣ a ∣ 2 + ∣ b ∣ 2 = 1 の複素数で作った ( a − b ‾ b a ‾ ) \left(\begin{smallmatrix} a & -\overline{b} \\ b & \overline{a} \end{smallmatrix}\right) ( a b − b a ) は S U ( 2 ) \mathrm{SU}(2) SU ( 2 ) の元です(すべてがこの形であることは 1-02 §05 )。
ユニタリ行列も対角化できます。根拠は複素のスペクトル定理です 既知 。A A † = A † A AA^\dagger = A^\dagger A A A † = A † A を満たす行列(正規行列 )は、ユニタリ行列で対角化できます。ユニタリ行列は U † U = U U † = I U^\dagger U = UU^\dagger = I U † U = U U † = I なのでこれに当たり(エルミート行列も当たります)、U = V diag ( e i θ 1 , … , e i θ n ) V † U = V\operatorname{diag}(e^{i\theta_1}, \dots, e^{i\theta_n})V^\dagger U = V diag ( e i θ 1 , … , e i θ n ) V † (V V V はユニタリ、θ k \theta_k θ k は実数)と書けます。
ゲージ理論でユニタリ行列が出るのは、物質の場が C N \mathbb{C}^N C N に値をとり、長さの二乗 ψ † ψ \psi^\dagger\psi ψ † ψ が作用に現れるからです。ゲージ変換はそれを保つように選ばれます(3-01 )。格子の上では、ゲージ場そのものが各辺のユニタリ行列です(5-01 )。
08
行列の指数関数
指数関数の級数は掛け算と足し算しか使わないので、行列にもそのまま使えます。
e X = I + X + X 2 2 ! + X 3 3 ! + ⋯ = ∑ k ≥ 0 X k k ! e^X = I + X + \frac{X^2}{2!} + \frac{X^3}{3!} + \cdots = \sum_{k \ge 0} \frac{X^k}{k!} e X = I + X + 2 ! X 2 + 3 ! X 3 + ⋯ = k ≥ 0 ∑ k ! X k
収束。 作用素ノルム ∥ X ∥ = max ∥ v ∥ = 1 ∥ X v ∥ \lVert X\rVert = \max_{\lVert v\rVert = 1}\lVert Xv\rVert ∥ X ∥ = max ∥ v ∥ = 1 ∥ X v ∥ (長さ 1 のベクトルを最大何倍にするか)は ∥ X Y ∥ ≤ ∥ X ∥ ∥ Y ∥ \lVert XY\rVert \le \lVert X\rVert\lVert Y\rVert ∥ X Y ∥ ≤ ∥ X ∥ ∥ Y ∥ を満たすので、各項は ∥ X ∥ k / k ! \lVert X\rVert^k/k! ∥ X ∥ k / k ! 以下、和は e ∥ X ∥ e^{\lVert X\rVert} e ∥ X ∥ 以下です。k = N k = N k = N までの部分和 S N ( X ) S_N(X) S N ( X ) の誤差は次のとおりです。
∥ e X − S N ( X ) ∥ ≤ ∥ X ∥ N + 1 ( N + 1 ) ! e ∥ X ∥ \bigl\lVert e^X - S_N(X)\bigr\rVert \le \frac{\lVert X\rVert^{N+1}}{(N+1)!}\,e^{\lVert X\rVert} e X − S N ( X ) ≤ ( N + 1 )! ∥ X ∥ N + 1 e ∥ X ∥
性質。 (a) d d t e t X = X e t X \frac{d}{dt}e^{tX} = Xe^{tX} d t d e tX = X e tX 。よって x ( t ) = e t X x 0 x(t) = e^{tX}x_0 x ( t ) = e tX x 0 は x ′ = X x x' = Xx x ′ = X x 、x ( 0 ) = x 0 x(0) = x_0 x ( 0 ) = x 0 の解です。(b) e ( s + t ) X = e s X e t X e^{(s+t)X} = e^{sX}e^{tX} e ( s + t ) X = e s X e tX 、特に ( e X ) − 1 = e − X (e^X)^{-1} = e^{-X} ( e X ) − 1 = e − X 。(c) P e X P − 1 = e P X P − 1 Pe^XP^{-1} = e^{PXP^{-1}} P e X P − 1 = e P X P − 1 。(d) det e X = e tr X \det e^X = e^{\operatorname{tr} X} det e X = e tr X 。対角化できれば (c) から ∏ k e λ k = e ∑ k λ k \prod_k e^{\lambda_k} = e^{\sum_k \lambda_k} ∏ k e λ k = e ∑ k λ k 、一般には上三角化で示せます 既知 。どの複素正方行列 X X X も、ユニタリ行列 V V V で V † X V V^\dagger XV V † X V を上三角にでき(Schur の定理)、上三角行列の指数関数は対角成分が e λ k e^{\lambda_k} e λ k の上三角行列です。
注意。 X Y = Y X XY = YX X Y = Y X なら e X + Y = e X e Y e^{X+Y} = e^Xe^Y e X + Y = e X e Y ですが、一般には違います。§03 の A , B A, B A , B は A 2 = B 2 = 0 A^2 = B^2 = 0 A 2 = B 2 = 0 なので e A = I + A e^A = I + A e A = I + A 、e B = I + B e^B = I + B e B = I + B で、
e A e B = ( 2 1 1 1 ) , e B e A = ( 1 1 1 2 ) , e^Ae^B = \begin{pmatrix} 2 & 1 \\ 1 & 1 \end{pmatrix},\qquad e^Be^A = \begin{pmatrix} 1 & 1 \\ 1 & 2 \end{pmatrix}, e A e B = ( 2 1 1 1 ) , e B e A = ( 1 1 1 2 ) ,
e A + B = ( cosh 1 sinh 1 sinh 1 cosh 1 ) ≈ ( 1.5431 1.1752 1.1752 1.5431 ) . e^{A+B} = \begin{pmatrix} \cosh 1 & \sinh 1 \\ \sinh 1 & \cosh 1 \end{pmatrix} \approx \begin{pmatrix} 1.5431 & 1.1752 \\ 1.1752 & 1.5431 \end{pmatrix}. e A + B = ( cosh 1 sinh 1 sinh 1 cosh 1 ) ≈ ( 1.5431 1.1752 1.1752 1.5431 ) .
計算 三つとも違います。差を交換子で書く公式は 1-03 の主題です。(d) は乱数の 4×4 複素行列 1000 個で相対誤差 5.7 × 10 − 14 5.7\times 10^{-14} 5.7 × 1 0 − 14 以内でした。
例:J = ( 0 − 1 1 0 ) J = \left(\begin{smallmatrix} 0 & -1 \\ 1 & 0 \end{smallmatrix}\right) J = ( 0 1 − 1 0 ) は J 2 = − I J^2 = -I J 2 = − I なので、偶数次と奇数次に分けて e t J = cos t I + sin t J = R ( t ) e^{tJ} = \cos t\,I + \sin t\,J = R(t) e t J = cos t I + sin t J = R ( t ) 、つまり回転です。一般に X 2 = − ω 2 I X^2 = -\omega^2 I X 2 = − ω 2 I なら e t X = cos ( ω t ) I + sin ( ω t ) ω X e^{tX} = \cos(\omega t)\,I + \frac{\sin(\omega t)}{\omega}X e tX = cos ( ω t ) I + ω s i n ( ω t ) X です。
図3
X = ( a − 1 1 a ) X = \left(\begin{smallmatrix} a & -1 \\ 1 & a \end{smallmatrix}\right) X = ( a 1 − 1 a ) 、
x 0 = ( 1 , 0 ) x_0 = (1, 0) x 0 = ( 1 , 0 ) 。実線が
e t X x 0 e^{tX}x_0 e tX x 0 、青の破線が部分和
S N ( t X ) x 0 S_N(tX)x_0 S N ( tX ) x 0 。
N N N を増やすと破線が実線に重なり、
t t t を延ばすと要る
N N N が増える。
a ≠ 0 a \ne 0 a = 0 で円が渦巻きになる(点線の半径
e a t e^{at} e a t )
図は JavaScript で描画します。a = 0 では軌道は単位円です。部分和との差は下の表のとおり、N を増やすと減り、t を延ばすと増えます。
N a t 動かす
JavaScript が無効なら下の表で読めます。
計算 図の部分和はその場で級数を足しています。表は倍精度の部分和と回転行列 e t J e^{tJ} e t J (a = 0 a = 0 a = 0 )との差を作用素ノルムで測った値です。表の数値のセルを押すと、その N N N と t t t (a = 0 a = 0 a = 0 )が図に入ります。
N N N t = 1 t = 1 t = 1 t = π t = \pi t = π t = 2 π t = 2\pi t = 2 π t = 4 π t = 4\pi t = 4 π
4 8.3×10⁻³ 2.3 57 1.0×10³
8 2.7×10⁻⁶ 7.9×10⁻² 36 1.3×10⁴
16 2.8×10⁻¹⁵ 7.8×10⁻⁷ 9.9×10⁻² 1.1×10⁴
30 1.1×10⁻¹⁶ 2.5×10⁻¹⁶ 6.6×10⁻¹⁰ 1.4
誤差 < 10⁻⁶ の最小の N N N 9 16 26 43
t t t が大きいと項 t k / k ! t^k/k! t k / k ! がいったん膨らむので、要る N N N は t t t とともに増えます。t = π t = \pi t = π 、N = 16 N = 16 N = 16 では、上の誤差の上界 1.8 × 10 − 5 1.8\times 10^{-5} 1.8 × 1 0 − 5 に対し実際の誤差は 7.8 × 10 − 7 7.8\times 10^{-7} 7.8 × 1 0 − 7 です。
寄り道:級数をそのまま足すと失敗する例
寄り道
飛ばしてよい小節です。
定義として正しい式が、計算法として良いとは限りません。Moler と Van Loan の総説(1978 年。2003 年に再録と追補)の例 既知 は A = ( − 49 24 − 64 31 ) A = \left(\begin{smallmatrix} -49 & 24 \\ -64 & 31 \end{smallmatrix}\right) A = ( − 49 − 64 24 31 ) (固有値 −1 と −17)です。e A e^A e A の成分は 1 程度なのに、級数の途中の項は約 9.3 × 10 6 9.3\times 10^6 9.3 × 1 0 6 まで膨らんで打ち消し合います。この端末で値が変わらなくなるまで足すと 計算 、単精度(約 7 桁)では 61 項で右下の成分が 1.1036 のはずが −0.8392 になり、倍精度では 74 項で誤差 5.7 × 10 − 9 5.7\times 10^{-9} 5.7 × 1 0 − 9 でした。論文は、問題は級数の打ち切りではなく算術の打ち切りにある、と述べています。実用のプログラム(SciPy の expm など)は、行列を 2 の冪で割ってから近似し、二乗を繰り返して戻す scaling and squaring を使います。
09
エルミートとユニタリを結ぶ
エルミート行列とユニタリ行列は指数関数でつながります。比喩で言えば、エルミート行列は「回し方」、ユニタリ行列は「回した結果」です。正確には次のとおりです。
H † = H ⟹ U = e i H はユニタリ , det e i H = e i tr H H^\dagger = H \;\Longrightarrow\; U = e^{iH} \text{ はユニタリ},\qquad \det e^{iH} = e^{i\operatorname{tr} H} H † = H ⟹ U = e i H はユニタリ , det e i H = e i tr H
項ごとに随伴を取ると ( e i H ) † = e − i H (e^{iH})^\dagger = e^{-iH} ( e i H ) † = e − i H 。i H iH i H と − i H -iH − i H は交換するので §08 (b) から U † U = e 0 = I U^\dagger U = e^0 = I U † U = e 0 = I 、行列式は (d) から出ます。tr H = 0 \operatorname{tr} H = 0 tr H = 0 なら det U = 1 \det U = 1 det U = 1 で U ∈ S U ( n ) U \in \mathrm{SU}(n) U ∈ SU ( n ) です。逆に、どのユニタリ行列も §07 の対角化から H = V diag ( θ k ) V † H = V\operatorname{diag}(\theta_k)V^\dagger H = V diag ( θ k ) V † と置いて e i H e^{iH} e i H と書けます 既知 (θ k \theta_k θ k は 2 π 2\pi 2 π の整数倍だけ選べるので、H H H は一つに決まりません)。
例:H = ( 0 1 1 0 ) H = \left(\begin{smallmatrix} 0 & 1 \\ 1 & 0 \end{smallmatrix}\right) H = ( 0 1 1 0 ) は H 2 = I H^2 = I H 2 = I なので e i θ H = cos θ I + i sin θ H e^{i\theta H} = \cos\theta\,I + i\sin\theta\,H e i θ H = cos θ I + i sin θ H 、行列式は cos 2 θ + sin 2 θ = 1 \cos^2\theta + \sin^2\theta = 1 cos 2 θ + sin 2 θ = 1 。計算 乱数の 3×3 エルミート行列 1000 個で ∥ U † U − I ∥ ≤ 1.6 × 10 − 15 \lVert U^\dagger U - I\rVert \le 1.6\times 10^{-15} ∥ U † U − I ∥ ≤ 1.6 × 1 0 − 15 、トレース 0 なら ∣ det U − 1 ∣ ≤ 1.1 × 10 − 15 |\det U - 1| \le 1.1\times 10^{-15} ∣ det U − 1∣ ≤ 1.1 × 1 0 − 15 でした。
この対応は二か所で使われます。量子力学の時間発展 e − i t H / ℏ e^{-itH/\hbar} e − i t H /ℏ はユニタリで、確率の合計を保ちます(2-05 )。ゲージ理論では、§02 の空間を N × N N \times N N × N に広げたもの(トレース 0 の N × N N \times N N × N エルミート行列、の i i i 倍)を指数関数で送った先が S U ( N ) \mathrm{SU}(N) SU ( N ) で、前者がゲージ場の値、後者がゲージ変換の空間です(1-03 ・1-04 )。実の成分に限れば、「反対称行列 K K K の e K e^K e K は回転」という形になります。
ゲージ場の値に「i i i 倍」を付けるのは幾何の書き方(1-07 )です。物理の書き方はゲージ場をエルミートな行列のまま書き、二つは A 幾何 = − i A 物理 A_{\text{幾何}}=-iA_{\text{物理}} A 幾何 = − i A 物理 で移り合います(結合定数を分けて書く章での形は 3-02 の表)。
この章が言えている範囲
言えている 言えていない
既知 有限次元のエルミート行列は実の固有値と、固有ベクトルからなる正規直交基底をもつ。ユニタリ行列も同様に対角化でき、固有値は絶対値 1。無限次元で同じことが言えるか。固有ベクトルをもたない作用素もある(1-10 )。
既知 e X e^X e X はいつも収束し、det e X = e tr X \det e^X = e^{\operatorname{tr} X} det e X = e tr X 。H H H がエルミートなら e i H e^{iH} e i H はユニタリで、どのユニタリ行列もこの形。e X + Y e^{X+Y} e X + Y と e X e Y e^Xe^Y e X e Y の差を交換子で表す公式(1-03 )。
計算 図と表の数値、乱数行列での恒等式の確認。乱数での確認は証明ではない。
既知 有限次元では、最小固有値と次の固有値の隙間はいつも計算できる。無限次元の作用素で、極限のあとも隙間が正のまま残るか。質量ギャップの問題はこの形をしており(4-11 )、この章はそれについて何も言わない。
出典と再現
もの 種別 出典・道具
内積の流儀、スペクトル定理(実・複素、正規行列)、上三角化(Schur の定理) 既知 S. Axler, Linear Algebra Done Right , 4th ed., Springer (Undergraduate Texts in Mathematics), 2024。定義 6.2、7.29、7.31、6.38。オープンアクセス版で本文を確認。
級数による計算の失敗例、scaling and squaring 既知 C. Moler, C. Van Loan, SIAM Review 45 (2003) 3–49(Method 1・3、本文を確認)。初出 SIAM Review 20 (1978) 801–836(書誌のみ)。SciPy の expm の方式は SciPy 1.11 の説明文による(A. H. Al-Mohy, N. J. Higham, SIAM J. Matrix Anal. Appl. 31 (2009) 970–989、書誌のみ)。
固有値の交差に要る条件の数 既知 J. von Neumann, E. Wigner, Physikalische Zeitschrift 30 (1929) 467–470(書誌のみ)。
図1・図2 の表 計算 この端末の Python 3・NumPy 1.26 の eig・eigh。トレース 0 の 2×2 エルミート行列 1000 個で式 ± x 2 + y 2 + z 2 \pm\sqrt{x^2+y^2+z^2} ± x 2 + y 2 + z 2 との差 2.2 × 10 − 15 2.2\times 10^{-15} 2.2 × 1 0 − 15 以内。
図3 の表、§08 の寄り道の再計算、§08・§09 の恒等式 計算 部分和を倍精度・単精度で計算(半精度は途中の項が上限を超えて計算不能)。厳密値は回転行列と P diag ( e − 1 , e − 17 ) P − 1 P\operatorname{diag}(e^{-1}, e^{-17})P^{-1} P diag ( e − 1 , e − 17 ) P − 1 (P = ( 1 3 2 4 ) P = \left(\begin{smallmatrix} 1 & 3 \\ 2 & 4 \end{smallmatrix}\right) P = ( 1 2 3 4 ) )。SciPy 1.11 の expm と乱数行列(シード固定、各 1000 個)。
次に読む章:1-02 群とは何か — ユニタリ行列の集まりが「群」をなすことから始め、S O ( 3 ) \mathrm{SO}(3) SO ( 3 ) と S U ( 2 ) \mathrm{SU}(2) SU ( 2 ) の関係へ進みます。
改訂 2026-10-01:初版。