非相対論的量子力学は、粒子数 N N N を固定した Hilbert 空間 L 2 ( R 3 N ) L^2(\mathbb{R}^{3N}) L 2 ( R 3 N ) の上に作られています。粒子の生成・消滅を書く場所が枠組みに存在しません。
相対論的な分散関係 E 2 = p 2 c 2 + m 2 c 4 E^2 = \boldsymbol{p}^2c^2 + m^2c^4 E 2 = p 2 c 2 + m 2 c 4 を素朴に演算子へ置き換えるとクライン・ゴルドン方程式が得られますが、負エネルギー解が必ず現れ、保存カレントの時間成分は負にもなります。確率密度として読めません。
ディラック方程式は確率密度の正値性を回復しますが、負エネルギー解は消えません。「ディラックの海」は実質的に無限個の粒子を持ち込む多体理論であり、ボソンには使えません。
さらに深刻なのは因果律です。H = − ∇ 2 + m 2 H = \sqrt{-\nabla^2 + m^2} H = − ∇ 2 + m 2 による一粒子の時間発展はコンパクト台を瞬時に壊し、空間的に隔たった点への伝播振幅もゼロになりません。
出口はひとつ、場 ϕ ( x ) \phi(x) ϕ ( x ) を量子化することです。自由場は各モードが独立な調和振動子になり、その励起量子が分散関係 E = p 2 + m 2 E = \sqrt{\boldsymbol{p}^2 + m^2} E = p 2 + m 2 をもつ粒子として現れます。
因果律は「振幅がゼロ」ではなく「空間的隔たりで交換子がゼロ」として実現されます。この打ち消しには反粒子が要り、反粒子の存在は相対論と量子力学の帰結です。
素朴な疑問から始めます。励起した原子が基底状態に落ちるとき、光子が 1 個放出されます。放出の前には光子は存在せず、後には存在します。ところが シュレーディンガー方程式と波動関数 の枠組み(状態と波動関数の定義(定義 2.1)[シュレーディンガー方程式と波動関数] )では、状態は N N N 個の粒子の座標の関数 ψ ( x 1 , … , x N , t ) \psi(\boldsymbol{x}_1,\ldots,\boldsymbol{x}_N,t) ψ ( x 1 , … , x N , t ) であり、N N N は方程式を書いた瞬間に固定されます。粒子が 1 個増える過程を書く場所がどこにもありません。
同じことは電子と陽電子の対消滅 e − e + → γ γ e^- e^+ \to \gamma\gamma e − e + → γ γ でも起こります。始状態は質量をもつ粒子 2 個、終状態は質量ゼロの粒子 2 個です。粒子の個数だけでなく種類まで変わります。
歴史的にも、この困難は量子力学の完成とほぼ同時に自覚されていました。シュレーディンガーが 1925 年末に最初に書いた波動方程式は相対論的なもの(現在クライン・ゴルドン方程式と呼ばれるもの)でしたが、水素原子の微細構造が実験と合わないため彼はこれを捨て、非相対論的な方程式を発表します。1926 年にクライン、ゴルドン、フォックらが独立に再発見した同じ方程式には、確率密度が負になるという病がついていました。
一方、ディラックは 1927 年の論文で電磁場そのものを量子化し、原子による光の放出・吸収を第一原理から計算してみせます。これが最初の場の量子論です。電磁場は古典的にも「場」であり、その量子化から光子という粒子が出てきました。自然な問いが立ちます。
電子もまた、何らかの場の量子なのではないか。
答えは「そうです」。ただしその結論に到達するには、一粒子の相対論的量子力学がどのように 破綻するかを正確に知る必要があります。以下では破綻を三つの層(負エネルギー、確率解釈、因果律)に分けて示し、場の量子化がそれらを同時に解くことを見ます。
flowchart TD
A["非相対論的量子力学<br/>粒子数は保存量"] --> B["相対論的分散関係を代入"]
B --> C["クライン・ゴルドン方程式"]
B --> D["ディラック方程式"]
C --> E["負エネルギー解<br/>確率密度が負になる"]
D --> F["負エネルギー解<br/>ディラックの海はボソンに使えない"]
E --> G["空間的隔たりでも伝播振幅がゼロにならない"]
F --> G
G --> H["場そのものを量子化する"]
H --> I["粒子 = 場の励起量子<br/>反粒子・生成消滅・ミクロ因果律"] 一粒子の相対論的量子力学が破綻し、場の量子化へ至る道筋
以下、特に断らない限り自然単位系 ℏ = c = 1 \hbar = c = 1 ℏ = c = 1 を使い、計量は η μ ν = d i a g ( + 1 , − 1 , − 1 , − 1 ) \eta_{\mu\nu} = \mathrm{diag}(+1,-1,-1,-1) η μν = diag ( + 1 , − 1 , − 1 , − 1 ) と取ります。時空点は x = ( t , x ) x = (t,\boldsymbol{x}) x = ( t , x ) 、内積は p ⋅ x = p 0 t − p ⋅ x p\cdot x = p^0 t - \boldsymbol{p}\cdot\boldsymbol{x} p ⋅ x = p 0 t − p ⋅ x 、質量殻上のエネルギーは ω p = p 2 + m 2 \omega_{\boldsymbol{p}} = \sqrt{\boldsymbol{p}^2 + m^2} ω p = p 2 + m 2 と書きます。□ = ∂ μ ∂ μ = ∂ t 2 − ∇ 2 \Box = \partial_\mu\partial^\mu = \partial_t^2 - \nabla^2 □ = ∂ μ ∂ μ = ∂ t 2 − ∇ 2 です。
非相対論的量子力学は、次の二つを暗黙に仮定しています。破綻はどちらもここから来ます。
(A1) 粒子数の固定。 状態空間は H N = L 2 ( R 3 N ) \mathcal{H}_N = L^2(\mathbb{R}^{3N}) H N = L 2 ( R 3 N ) (同種粒子なら対称化または反対称化した部分空間)です。ハミルトニアンは H N \mathcal{H}_N H N を H N \mathcal{H}_N H N に写すので、粒子数 N N N は定義により 保存します。保存則を証明したのではなく、保存しない状況を書けない空間を選んだだけです。
(A2) 時間と空間の非対称性。 位置 x ^ \hat{\boldsymbol{x}} x ^ は演算子ですが、時刻 t t t は演算子ではなく単なるパラメータです。正準交換関係は [ x ^ i , p ^ j ] = i δ i j [\hat{x}^i, \hat{p}^j] = i\delta^{ij} [ x ^ i , p ^ j ] = i δ ij であり、t t t はここに現れません。ところが ローレンツ変換 (ブーストの変換則(定理 6.1)[ローレンツ変換] )は t t t と x \boldsymbol{x} x を混ぜます。「演算子である量」と「パラメータである量」が変換で混ざってしまうので、この区別はローレンツ不変ではありません。
非相対論では E = p 2 / 2 m E = \boldsymbol{p}^2/2m E = p 2 /2 m に E → i ∂ t E \to i\partial_t E → i ∂ t 、p → − i ∇ \boldsymbol{p}\to -i\nabla p → − i ∇ を代入してシュレーディンガー方程式を得ました。相対論的な関係式 E 2 = p 2 + m 2 E^2 = \boldsymbol{p}^2 + m^2 E 2 = p 2 + m 2 (相対論的力学 の エネルギー・運動量関係式(定理 4.4)[相対論的力学] )に同じ置き換えを施します。
定義 3.1 (クライン・ゴルドン方程式 )
質量 m > 0 m > 0 m > 0 の複素スカラー場(あるいは一粒子波動関数の候補)ϕ ( t , x ) \phi(t,\boldsymbol{x}) ϕ ( t , x ) に対する方程式
( □ + m 2 ) ϕ = 0 , □ = ∂ t 2 − ∇ 2 \left(\Box + m^2\right)\phi = 0,
\qquad \Box = \partial_t^2 - \nabla^2 ( □ + m 2 ) ϕ = 0 , □ = ∂ t 2 − ∇ 2 をクライン・ゴルドン方程式と呼びます。
平面波 ϕ = N e − i ( E t − p ⋅ x ) \phi = N e^{-i(Et - \boldsymbol{p}\cdot\boldsymbol{x})} ϕ = N e − i ( E t − p ⋅ x ) を代入すると − E 2 + p 2 + m 2 = 0 -E^2 + \boldsymbol{p}^2 + m^2 = 0 − E 2 + p 2 + m 2 = 0 、すなわち
E = ± ω p = ± p 2 + m 2 E = \pm\,\omega_{\boldsymbol{p}} = \pm\sqrt{\boldsymbol{p}^2 + m^2} E = ± ω p = ± p 2 + m 2
の二つの枝が出ます。負の枝は「エネルギーが下に有界でない」ことを意味し、系を摂動すればいくらでも低い状態へ落ち続けられます。
「負エネルギー解は非物理的だから捨てればよい」と言いたくなりますが、捨てられません。理由は初期値問題の構造にあります。クライン・ゴルドン方程式は時間について 2 階なので、Cauchy データは ϕ ( 0 , x ) \phi(0,\boldsymbol{x}) ϕ ( 0 , x ) と ϕ ˙ ( 0 , x ) \dot\phi(0,\boldsymbol{x}) ϕ ˙ ( 0 , x ) の二つ の任意関数です。空間 Fourier 変換すると各 p \boldsymbol{p} p について ϕ ~ ¨ = − ω p 2 ϕ ~ \ddot{\tilde\phi} = -\omega_{\boldsymbol p}^2\tilde\phi ϕ ~ ¨ = − ω p 2 ϕ ~ となり、一般解は
ϕ ~ ( t , p ) = a ( p ) e − i ω p t + b ( p ) e + i ω p t \tilde\phi(t,\boldsymbol{p}) = a(\boldsymbol{p})e^{-i\omega_{\boldsymbol p}t} + b(\boldsymbol{p})e^{+i\omega_{\boldsymbol p}t} ϕ ~ ( t , p ) = a ( p ) e − i ω p t + b ( p ) e + i ω p t
です。正振動数だけに制限する(b ≡ 0 b \equiv 0 b ≡ 0 とする)ことは、ϕ ~ ˙ ( 0 , p ) = − i ω p ϕ ~ ( 0 , p ) \dot{\tilde\phi}(0,\boldsymbol{p}) = -i\omega_{\boldsymbol{p}}\tilde\phi(0,\boldsymbol{p}) ϕ ~ ˙ ( 0 , p ) = − i ω p ϕ ~ ( 0 , p ) という拘束条件 を課すことに他なりません。位置空間では ϕ ˙ ( 0 , ⋅ ) = − i − ∇ 2 + m 2 ϕ ( 0 , ⋅ ) \dot\phi(0,\cdot) = -i\sqrt{-\nabla^2+m^2}\,\phi(0,\cdot) ϕ ˙ ( 0 , ⋅ ) = − i − ∇ 2 + m 2 ϕ ( 0 , ⋅ ) という非局所的な関係で、初期条件を局所的に自由に指定する権利を放棄したことになります。この非局所性が何を壊すかは 命題 5.1 で見ます。
シュレーディンガー方程式では ρ = ∣ ψ ∣ 2 ≥ 0 \rho = \lvert\psi\rvert^2 \ge 0 ρ = ∣ ψ ∣ 2 ≥ 0 が連続の方程式を満たし(確率の連続の方程式(定理 4.2)[シュレーディンガー方程式と波動関数] )、確率密度として読めました。クライン・ゴルドン方程式でも保存カレントは作れますが、正値性が失われます。
命題 3.2 (クライン・ゴルドンのカレントとその非正値性 )
ϕ \phi ϕ が 定義 3.1 の解であるとき、
j μ = i 2 m ( ϕ ∗ ∂ μ ϕ − ( ∂ μ ϕ ∗ ) ϕ ) j^\mu = \frac{i}{2m}\left(\phi^{*}\partial^\mu\phi - (\partial^\mu\phi^{*})\phi\right) j μ = 2 m i ( ϕ ∗ ∂ μ ϕ − ( ∂ μ ϕ ∗ ) ϕ ) は ∂ μ j μ = 0 \partial_\mu j^\mu = 0 ∂ μ j μ = 0 を満たす。さらに平面波解 ϕ = N e − i ( E t − p ⋅ x ) \phi = Ne^{-i(Et-\boldsymbol{p}\cdot\boldsymbol{x})} ϕ = N e − i ( E t − p ⋅ x ) (E = ± ω p E = \pm\omega_{\boldsymbol p} E = ± ω p )に対して
j 0 = E m ∣ N ∣ 2 j^0 = \frac{E}{m}\lvert N\rvert^2 j 0 = m E ∣ N ∣ 2 となる。したがって負エネルギー解では j 0 < 0 j^0 < 0 j 0 < 0 であり、j 0 j^0 j 0 は確率密度として解釈できない。
証明(命題 3.2) まず保存を確かめます。積の微分から
∂ μ j μ = i 2 m ( ∂ μ ϕ ∗ ∂ μ ϕ + ϕ ∗ □ ϕ − ( □ ϕ ∗ ) ϕ − ∂ μ ϕ ∗ ∂ μ ϕ ) = i 2 m ( ϕ ∗ □ ϕ − ( □ ϕ ∗ ) ϕ ) . \partial_\mu j^\mu = \frac{i}{2m}\left(\partial_\mu\phi^{*}\partial^\mu\phi + \phi^{*}\Box\phi - (\Box\phi^{*})\phi - \partial_\mu\phi^{*}\partial^\mu\phi\right)
= \frac{i}{2m}\left(\phi^{*}\Box\phi - (\Box\phi^{*})\phi\right). ∂ μ j μ = 2 m i ( ∂ μ ϕ ∗ ∂ μ ϕ + ϕ ∗ □ ϕ − ( □ ϕ ∗ ) ϕ − ∂ μ ϕ ∗ ∂ μ ϕ ) = 2 m i ( ϕ ∗ □ ϕ − ( □ ϕ ∗ ) ϕ ) . 第 1 項と第 4 項が相殺しました。ここで 定義 3.1 より □ ϕ = − m 2 ϕ \Box\phi = -m^2\phi □ ϕ = − m 2 ϕ 、複素共役を取って □ ϕ ∗ = − m 2 ϕ ∗ \Box\phi^{*} = -m^2\phi^{*} □ ϕ ∗ = − m 2 ϕ ∗ ですから
∂ μ j μ = i 2 m ( − m 2 ϕ ∗ ϕ + m 2 ϕ ∗ ϕ ) = 0 \partial_\mu j^\mu = \frac{i}{2m}\left(-m^2\phi^{*}\phi + m^2\phi^{*}\phi\right) = 0 ∂ μ j μ = 2 m i ( − m 2 ϕ ∗ ϕ + m 2 ϕ ∗ ϕ ) = 0 です。次に j 0 j^0 j 0 を計算します。計量の取り方から ∂ 0 = η 00 ∂ 0 = ∂ t \partial^0 = \eta^{00}\partial_0 = \partial_t ∂ 0 = η 00 ∂ 0 = ∂ t なので j 0 = i 2 m ( ϕ ∗ ϕ ˙ − ϕ ˙ ∗ ϕ ) j^0 = \frac{i}{2m}(\phi^{*}\dot\phi - \dot\phi^{*}\phi) j 0 = 2 m i ( ϕ ∗ ϕ ˙ − ϕ ˙ ∗ ϕ ) です。平面波では ϕ ˙ = − i E ϕ \dot\phi = -iE\phi ϕ ˙ = − i E ϕ 、ϕ ˙ ∗ = + i E ϕ ∗ \dot\phi^{*} = +iE\phi^{*} ϕ ˙ ∗ = + i E ϕ ∗ ですから
j 0 = i 2 m ( ∣ N ∣ 2 ( − i E ) − ( i E ) ∣ N ∣ 2 ) = i 2 m ( − 2 i E ∣ N ∣ 2 ) = E m ∣ N ∣ 2 . j^0 = \frac{i}{2m}\left(\lvert N\rvert^2(-iE) - (iE)\lvert N\rvert^2\right) = \frac{i}{2m}\left(-2iE\lvert N\rvert^2\right) = \frac{E}{m}\lvert N\rvert^2 . j 0 = 2 m i ( ∣ N ∣ 2 ( − i E ) − ( i E ) ∣ N ∣ 2 ) = 2 m i ( − 2 i E ∣ N ∣ 2 ) = m E ∣ N ∣ 2 . E = − ω p < 0 E = -\omega_{\boldsymbol p} < 0 E = − ω p < 0 の枝では右辺は負です。密度が負になる量は確率密度ではありません。
∎
例 3.3 (非相対論極限では正しい確率密度に戻る )
j 0 j^0 j 0 が完全に無意味というわけではありません。静止エネルギーを取り出して ϕ ( t , x ) = e − i m t ψ ( t , x ) \phi(t,\boldsymbol{x}) = e^{-imt}\psi(t,\boldsymbol{x}) ϕ ( t , x ) = e − im t ψ ( t , x ) と書き、非相対論極限 ∣ ψ ˙ ∣ ≪ m ∣ ψ ∣ \lvert\dot\psi\rvert \ll m\lvert\psi\rvert ∣ ψ ˙ ∣ ≪ m ∣ ψ ∣ を仮定します。すると
ϕ ˙ = e − i m t ( ψ ˙ − i m ψ ) ≃ − i m e − i m t ψ \dot\phi = e^{-imt}\left(\dot\psi - im\psi\right) \simeq -im\,e^{-imt}\psi ϕ ˙ = e − im t ( ψ ˙ − im ψ ) ≃ − im e − im t ψ なので、命題 3.2 の j 0 j^0 j 0 の表式に代入して
j 0 ≃ i 2 m ( ψ ∗ ( − i m ψ ) − ( i m ψ ∗ ) ψ ) = i 2 m ( − 2 i m ∣ ψ ∣ 2 ) = ∣ ψ ∣ 2 j^0 \simeq \frac{i}{2m}\left(\psi^{*}(-im\psi) - (im\psi^{*})\,\psi\right) = \frac{i}{2m}\left(-2im\lvert\psi\rvert^2\right) = \lvert\psi\rvert^2 j 0 ≃ 2 m i ( ψ ∗ ( − im ψ ) − ( im ψ ∗ ) ψ ) = 2 m i ( − 2 im ∣ ψ ∣ 2 ) = ∣ ψ ∣ 2 が得られます。シュレーディンガー方程式の確率密度に一致します。つまり j 0 j^0 j 0 は非相対論極限で確率密度に「見える」量であり、相対論的な領域で符号を変えるのです。
ディラックは 1928 年、確率密度の非正値性の原因を「時間について 2 階であること」に求めました。1 階の方程式なら、シュレーディンガー方程式と同様に ρ = ψ † ψ ≥ 0 \rho = \psi^\dagger\psi \ge 0 ρ = ψ † ψ ≥ 0 が保存量になるはずです。
定義 4.1 (ディラック方程式 )
N N N 成分の波動関数 ψ ( t , x ) ∈ C N \psi(t,\boldsymbol{x}) \in \mathbb{C}^N ψ ( t , x ) ∈ C N に対する 1 階の方程式
i ∂ t ψ = H ^ ψ , H ^ = α ⋅ p ^ + β m i\,\partial_t \psi = \hat H\psi,
\qquad
\hat H = \boldsymbol{\alpha}\cdot\hat{\boldsymbol{p}} + \beta m i ∂ t ψ = H ^ ψ , H ^ = α ⋅ p ^ + β m をディラック方程式と呼びます。ここで α 1 , α 2 , α 3 , β \alpha^1,\alpha^2,\alpha^3,\beta α 1 , α 2 , α 3 , β は N × N N\times N N × N のエルミート行列で、H ^ 2 = p ^ 2 + m 2 \hat H^2 = \hat{\boldsymbol{p}}^2 + m^2 H ^ 2 = p ^ 2 + m 2 が成り立つよう
{ α i , α j } = 2 δ i j 1 , { α i , β } = 0 , β 2 = 1 \{\alpha^i,\alpha^j\} = 2\delta^{ij}\mathbb{1},
\qquad
\{\alpha^i,\beta\} = 0,
\qquad
\beta^2 = \mathbb{1} { α i , α j } = 2 δ ij 1 , { α i , β } = 0 , β 2 = 1 を満たすものとします。
この関係式が H ^ 2 = p ^ 2 + m 2 \hat H^2 = \hat{\boldsymbol{p}}^2 + m^2 H ^ 2 = p ^ 2 + m 2 と同値であることは、展開すれば見えます。
H ^ 2 = α i α j p ^ i p ^ j + m ( α i β + β α i ) p ^ i + β 2 m 2 = 1 2 { α i , α j } p ^ i p ^ j + m { α i , β } p ^ i + β 2 m 2 . \hat H^2 = \alpha^i\alpha^j \hat p_i\hat p_j + m\left(\alpha^i\beta + \beta\alpha^i\right)\hat p_i + \beta^2 m^2
= \tfrac{1}{2}\{\alpha^i,\alpha^j\}\hat p_i\hat p_j + m\{\alpha^i,\beta\}\hat p_i + \beta^2m^2 . H ^ 2 = α i α j p ^ i p ^ j + m ( α i β + β α i ) p ^ i + β 2 m 2 = 2 1 { α i , α j } p ^ i p ^ j + m { α i , β } p ^ i + β 2 m 2 .
第 1 項で p ^ i p ^ j \hat p_i \hat p_j p ^ i p ^ j が i ↔ j i \leftrightarrow j i ↔ j について対称なので α i α j \alpha^i\alpha^j α i α j の対称部分だけが効き、それが反交換子の半分です。よって上の三条件のもとで H ^ 2 = δ i j p ^ i p ^ j + m 2 = p ^ 2 + m 2 \hat H^2 = \delta^{ij}\hat p_i\hat p_j + m^2 = \hat{\boldsymbol{p}}^2 + m^2 H ^ 2 = δ ij p ^ i p ^ j + m 2 = p ^ 2 + m 2 となります。
H ^ \hat H H ^ がエルミートなので、ρ = ψ † ψ \rho = \psi^\dagger\psi ρ = ψ † ψ と j = ψ † α ψ \boldsymbol{j} = \psi^\dagger\boldsymbol{\alpha}\psi j = ψ † α ψ が連続の方程式 ∂ t ρ + ∇ ⋅ j = 0 \partial_t\rho + \nabla\cdot\boldsymbol{j} = 0 ∂ t ρ + ∇ ⋅ j = 0 を満たし、しかも ρ ≥ 0 \rho \ge 0 ρ ≥ 0 です。命題 3.2 の病その 2 は解決しました。ところが病その 1 は残ります。
定理 4.2 (ディラック・ハミルトニアンは下に有界でない )
定義 4.1 の条件を満たす N × N N\times N N × N エルミート行列 α i , β \alpha^i,\beta α i , β が存在するとする。運動量固有値 p \boldsymbol{p} p に対する行列 H ( p ) = α ⋅ p + β m H(\boldsymbol{p}) = \boldsymbol{\alpha}\cdot\boldsymbol{p} + \beta m H ( p ) = α ⋅ p + β m について、次が成り立つ。
tr α i = 0 \operatorname{tr}\alpha^i = 0 tr α i = 0 (i = 1 , 2 , 3 i=1,2,3 i = 1 , 2 , 3 )かつ tr β = 0 \operatorname{tr}\beta = 0 tr β = 0 。したがって tr H ( p ) = 0 \operatorname{tr}H(\boldsymbol{p}) = 0 tr H ( p ) = 0 。
H ( p ) H(\boldsymbol{p}) H ( p ) の固有値は + ω p +\omega_{\boldsymbol p} + ω p と − ω p -\omega_{\boldsymbol p} − ω p のみで、両者の重複度は等しく N / 2 N/2 N /2 である。特に N N N は偶数で、H ^ \hat H H ^ のスペクトルは ( − ∞ , − m ] ∪ [ m , ∞ ) (-\infty,-m]\cup[m,\infty) ( − ∞ , − m ] ∪ [ m , ∞ ) となり下に有界でない。
証明(定理 4.2) 1 の証明。 { α i , β } = 0 \{\alpha^i,\beta\}=0 { α i , β } = 0 より β α i = − α i β \beta\alpha^i = -\alpha^i\beta β α i = − α i β ですから、両辺に右から β \beta β を掛けて β 2 = 1 \beta^2=\mathbb{1} β 2 = 1 を使うと
β α i β = − α i β 2 = − α i . \beta\alpha^i\beta = -\alpha^i\beta^2 = -\alpha^i . β α i β = − α i β 2 = − α i . 両辺のトレースを取ります。左辺はトレースの巡回性から tr ( β α i β ) = tr ( α i β 2 ) = tr α i \operatorname{tr}(\beta\alpha^i\beta) = \operatorname{tr}(\alpha^i\beta^2) = \operatorname{tr}\alpha^i tr ( β α i β ) = tr ( α i β 2 ) = tr α i です。よって tr α i = − tr α i \operatorname{tr}\alpha^i = -\operatorname{tr}\alpha^i tr α i = − tr α i 、すなわち tr α i = 0 \operatorname{tr}\alpha^i = 0 tr α i = 0 。同様に ( α 1 ) 2 = 1 (\alpha^1)^2 = \mathbb{1} ( α 1 ) 2 = 1 と { α 1 , β } = 0 \{\alpha^1,\beta\}=0 { α 1 , β } = 0 から α 1 β α 1 = − β ( α 1 ) 2 = − β \alpha^1\beta\alpha^1 = -\beta(\alpha^1)^2 = -\beta α 1 β α 1 = − β ( α 1 ) 2 = − β となり、トレースを取って tr β = − tr β = 0 \operatorname{tr}\beta = -\operatorname{tr}\beta = 0 tr β = − tr β = 0 を得ます。線形性より tr H ( p ) = p i tr α i + m tr β = 0 \operatorname{tr}H(\boldsymbol{p}) = p_i\operatorname{tr}\alpha^i + m\operatorname{tr}\beta = 0 tr H ( p ) = p i tr α i + m tr β = 0 です。
2 の証明。 定義 4.1 の直後で確かめたとおり H ( p ) 2 = ( p 2 + m 2 ) 1 = ω p 2 1 H(\boldsymbol{p})^2 = (\boldsymbol{p}^2+m^2)\mathbb{1} = \omega_{\boldsymbol p}^2\mathbb{1} H ( p ) 2 = ( p 2 + m 2 ) 1 = ω p 2 1 です。H ( p ) H(\boldsymbol p) H ( p ) はエルミートなので対角化でき、固有値 λ \lambda λ は λ 2 = ω p 2 \lambda^2 = \omega_{\boldsymbol p}^2 λ 2 = ω p 2 、すなわち λ = ± ω p \lambda = \pm\omega_{\boldsymbol p} λ = ± ω p に限られます。重複度を n + , n − n_+,n_- n + , n − とすると n + + n − = N n_++n_-=N n + + n − = N であり、1 より
tr H ( p ) = n + ω p − n − ω p = 0. \operatorname{tr}H(\boldsymbol p) = n_+\omega_{\boldsymbol p} - n_-\omega_{\boldsymbol p} = 0 . tr H ( p ) = n + ω p − n − ω p = 0. m > 0 m > 0 m > 0 より ω p ≥ m > 0 \omega_{\boldsymbol p} \ge m > 0 ω p ≥ m > 0 ですから n + = n − = N / 2 n_+ = n_- = N/2 n + = n − = N /2 です。N ≥ 1 N \ge 1 N ≥ 1 なので n − ≥ 1 n_- \ge 1 n − ≥ 1 、つまり負の固有値は必ず存在します。p \boldsymbol{p} p を全空間で動かせば − ω p -\omega_{\boldsymbol p} − ω p は − m -m − m から − ∞ -\infty − ∞ まで連続的に値を取り、スペクトルは下に有界でありません。
∎
つまり方程式を 1 階にしても負エネルギー解は消えません。定理 4.2 は、消えないどころかちょうど半分 が負エネルギーであることを言っています。
例 4.3 (ディラックの海とその代償 )
ディラックの処方はこうです。電子はフェルミオンなのでパウリ排他律が働く。そこで真空を「すべての負エネルギー準位が電子で埋まった状態」と定義すると、正エネルギー電子は排他律により負エネルギー準位へ落ちられません。海に空いた穴(ホール)は電荷 + e +e + e 、正エネルギーの粒子のようにふるまいます。これが 1932 年にアンダーソンが霧箱で発見した陽電子で、予言としては見事な成功でした。
しかし代償は大きいものです。
真空が無限個の電子を含み、電荷密度もエネルギー密度も発散します。観測量は「真空からの差」として定義し直さねばなりません。
排他律に依存するので、ボソンには使えません 。スピン 0 のパイ中間子にも反粒子(π + \pi^+ π + と π − \pi^- π − )があることは実験事実です。
決定的なのは、これがもはや一粒子理論ではないことです。無限個の粒子を持ち込んだ時点で記述は多体理論になっており、この処方自身が「一粒子の相対論的量子力学は成り立たない」ことを認めています。
負エネルギー解の問題には、まだ「正振動数部分だけを取り出せばよい」という逃げ道が残っているように見えます。すなわち Hilbert 空間を正エネルギー解に制限し、ハミルトニアンとして
H = − ∇ 2 + m 2 H = \sqrt{-\nabla^2 + m^2} H = − ∇ 2 + m 2
を採用するのです。これは L 2 ( R 3 ) L^2(\mathbb{R}^3) L 2 ( R 3 ) 上の自己共役作用素でスペクトルは [ m , ∞ ) [m,\infty) [ m , ∞ ) 、確かに下に有界で、定理 4.2 の困難は形式的に回避されます。しかしここでより深い病が現れます。相対論の核心である因果律 が破れるのです。
命題 5.1 (正エネルギー粒子は局在を保てない )
m > 0 m > 0 m > 0 、H = − ∇ 2 + m 2 H = \sqrt{-\nabla^2+m^2} H = − ∇ 2 + m 2 、ψ t = e − i t H ψ 0 \psi_t = e^{-itH}\psi_0 ψ t = e − i t H ψ 0 とする。ψ 0 ∈ L 2 ( R 3 ) \psi_0 \in L^2(\mathbb{R}^3) ψ 0 ∈ L 2 ( R 3 ) がコンパクトな台をもち ψ 0 ≠ 0 \psi_0 \ne 0 ψ 0 = 0 ならば、t ≠ 0 t \ne 0 t = 0 であるすべての時刻で ψ t \psi_t ψ t はコンパクトな台をもたない。
証明(命題 5.1) 空間 Fourier 変換すると ψ ^ t ( p ) = e − i t ω p ψ ^ 0 ( p ) \hat\psi_t(\boldsymbol{p}) = e^{-it\omega_{\boldsymbol p}}\hat\psi_0(\boldsymbol{p}) ψ ^ t ( p ) = e − i t ω p ψ ^ 0 ( p ) 、ω p = p 2 + m 2 \omega_{\boldsymbol p} = \sqrt{\boldsymbol{p}^2+m^2} ω p = p 2 + m 2 です。
背理法で示します。ある t ≠ 0 t\ne 0 t = 0 で ψ t \psi_t ψ t もコンパクトな台をもったとします。Paley–Wiener–Schwartz の定理により、コンパクトな台をもつ L 2 L^2 L 2 関数の Fourier 変換は C 3 \mathbb{C}^3 C 3 上の整関数(指数型)に解析接続されます。よって ψ ^ 0 \hat\psi_0 ψ ^ 0 と ψ ^ t \hat\psi_t ψ ^ t はともに整関数 F 0 , F t F_0, F_t F 0 , F t に延びます。
ψ 0 ≠ 0 \psi_0 \ne 0 ψ 0 = 0 より F 0 F_0 F 0 は恒等的にゼロではなく、F 0 F_0 F 0 は実点上で連続なので F 0 ( p 0 ) ≠ 0 F_0(\boldsymbol{p}_0)\ne 0 F 0 ( p 0 ) = 0 となる実ベクトル p 0 \boldsymbol{p}_0 p 0 が取れます。実単位ベクトル e \boldsymbol{e} e を固定し、複素直線 p = p 0 + z e \boldsymbol{p} = \boldsymbol{p}_0 + z\boldsymbol{e} p = p 0 + z e (z ∈ C z\in\mathbb{C} z ∈ C )に制限して f 0 ( z ) = F 0 ( p 0 + z e ) f_0(z) = F_0(\boldsymbol{p}_0+z\boldsymbol{e}) f 0 ( z ) = F 0 ( p 0 + z e ) 、f t ( z ) = F t ( p 0 + z e ) f_t(z) = F_t(\boldsymbol{p}_0+z\boldsymbol{e}) f t ( z ) = F t ( p 0 + z e ) とおきます。どちらも z z z の整関数で、f 0 ( 0 ) ≠ 0 f_0(0)\ne 0 f 0 ( 0 ) = 0 です。
この直線上で
Q ( z ) : = ( p 0 + z e ) 2 + m 2 = z 2 + 2 ( p 0 ⋅ e ) z + p 0 2 + m 2 Q(z) := (\boldsymbol{p}_0+z\boldsymbol{e})^2 + m^2 = z^2 + 2(\boldsymbol{p}_0\cdot\boldsymbol{e})z + \boldsymbol{p}_0^2 + m^2 Q ( z ) := ( p 0 + z e ) 2 + m 2 = z 2 + 2 ( p 0 ⋅ e ) z + p 0 2 + m 2 です。判別式は 4 ( p 0 ⋅ e ) 2 − 4 ( p 0 2 + m 2 ) 4(\boldsymbol{p}_0\cdot\boldsymbol{e})^2 - 4(\boldsymbol{p}_0^2+m^2) 4 ( p 0 ⋅ e ) 2 − 4 ( p 0 2 + m 2 ) で、Cauchy–Schwarz より ( p 0 ⋅ e ) 2 ≤ p 0 2 < p 0 2 + m 2 (\boldsymbol{p}_0\cdot\boldsymbol{e})^2 \le \boldsymbol{p}_0^2 < \boldsymbol{p}_0^2+m^2 ( p 0 ⋅ e ) 2 ≤ p 0 2 < p 0 2 + m 2 (ここで m > 0 m>0 m > 0 を使いました)なので負です。よって Q Q Q は相異なる二つの非実根 z ± z_\pm z ± をもち、Q ( z ) = ( z − z + ) ( z − z − ) \sqrt{Q(z)} = \sqrt{(z-z_+)(z-z_-)} Q ( z ) = ( z − z + ) ( z − z − ) は z ± z_\pm z ± に真の分岐点をもちます。
実軸上では f t ( z ) = g ( z ) f 0 ( z ) f_t(z) = g(z)f_0(z) f t ( z ) = g ( z ) f 0 ( z ) 、g ( z ) = e − i t Q ( z ) g(z) = e^{-it\sqrt{Q(z)}} g ( z ) = e − i t Q ( z ) が成り立ちます。f 0 , f t f_0, f_t f 0 , f t はともに整関数で f 0 ≢ 0 f_0\not\equiv 0 f 0 ≡ 0 ですから、商 f t / f 0 f_t/f_0 f t / f 0 は C \mathbb{C} C 上の一価 の有理型関数です。一方 g g g を z + z_+ z + のまわりで 1 周解析接続すると Q \sqrt{Q} Q の符号が反転し、g → e + i t Q g \to e^{+it\sqrt{Q}} g → e + i t Q に移ります。t ≠ 0 t\ne 0 t = 0 のとき e − 2 i t Q ( z ) = 1 e^{-2it\sqrt{Q(z)}} = 1 e − 2 i t Q ( z ) = 1 となるのは Q ( z ) ∈ ( π / t ) Z \sqrt{Q(z)} \in (\pi/t)\mathbb{Z} Q ( z ) ∈ ( π / t ) Z を満たす孤立点だけなので、g g g は z + z_+ z + の近傍で一価ではありません。一価な有理型関数と一価でない関数が実軸上で一致することはできず、矛盾します。
∎
台がコンパクトでなくなるということは、初期時刻に有限領域に閉じ込められていた粒子が、任意に小さい t > 0 t > 0 t > 0 の後には任意に遠い場所に非ゼロの振幅をもつということです。光速の制限が効いていません。次にその振幅の大きさを見積もります。
例 5.2 (光円錐の外への伝播振幅 )
時刻 0 0 0 に点 x \boldsymbol{x} x にあった粒子が、時刻 t t t に点 y \boldsymbol{y} y に見出される振幅は
U ( t , r ) = ⟨ y ∣ e − i H t ∣ x ⟩ = ∫ d 3 p ( 2 π ) 3 e i p ⋅ ( y − x ) e − i ω p t , r = ∣ y − x ∣ U(t,r) = \langle \boldsymbol{y}\,\lvert\, e^{-iHt}\,\rvert\,\boldsymbol{x}\rangle
= \int\frac{d^3p}{(2\pi)^3}\;e^{i\boldsymbol{p}\cdot(\boldsymbol{y}-\boldsymbol{x})}\,e^{-i\omega_{\boldsymbol p}t},
\qquad r = \lvert\boldsymbol{y}-\boldsymbol{x}\rvert U ( t , r ) = ⟨ y ∣ e − i H t ∣ x ⟩ = ∫ ( 2 π ) 3 d 3 p e i p ⋅ ( y − x ) e − i ω p t , r = ∣ y − x ∣ です。角度積分を実行すると(∫ d Ω e i p r cos θ = 4 π sin ( p r ) / ( p r ) \int d\Omega\, e^{ipr\cos\theta} = 4\pi\sin(pr)/(pr) ∫ d Ω e i p r c o s θ = 4 π sin ( p r ) / ( p r ) )
U ( t , r ) = 1 2 π 2 r ∫ 0 ∞ d p p sin ( p r ) e − i ω p t U(t,r) = \frac{1}{2\pi^2 r}\int_0^\infty dp\;p\,\sin(pr)\,e^{-i\omega_p t} U ( t , r ) = 2 π 2 r 1 ∫ 0 ∞ d p p sin ( p r ) e − i ω p t となります。p 2 + m 2 \sqrt{p^2+m^2} p 2 + m 2 の分岐点 p = ± i m p = \pm im p = ± im のまわりに積分路を変形すると(計算は Appendix)、空間的隔たり r > t > 0 r > t > 0 r > t > 0 に対して
U ( t , r ) = i 2 π 2 r ∫ m ∞ d ρ ρ e − ρ r sinh ( t ρ 2 − m 2 ) U(t,r) = \frac{i}{2\pi^2 r}\int_m^\infty d\rho\;\rho\,e^{-\rho r}\,\sinh\!\left(t\sqrt{\rho^2-m^2}\right) U ( t , r ) = 2 π 2 r i ∫ m ∞ d ρ ρ e − ρ r sinh ( t ρ 2 − m 2 ) という表式が得られます。積分区間の内部で被積分関数は厳密に正 です(ρ > m \rho > m ρ > m で sinh ( t ρ 2 − m 2 ) > 0 \sinh(t\sqrt{\rho^2-m^2}) > 0 sinh ( t ρ 2 − m 2 ) > 0 )。したがって U ( t , r ) ≠ 0 U(t,r) \ne 0 U ( t , r ) = 0 です。光円錐の外でも振幅は消えません。
大きさは、m r ≫ 1 mr \gg 1 m r ≫ 1 かつ r ≫ t r \gg t r ≫ t のとき積分の下端 ρ ≃ m \rho\simeq m ρ ≃ m が支配することから読めます。ρ = m + u \rho = m+u ρ = m + u とおいて ρ 2 − m 2 ≃ 2 m u \sqrt{\rho^2-m^2}\simeq\sqrt{2mu} ρ 2 − m 2 ≃ 2 m u 、sinh x ≃ x \sinh x\simeq x sinh x ≃ x とし、∫ 0 ∞ d u u e − u r = π / ( 2 r 3 / 2 ) \int_0^\infty du\,\sqrt{u}\,e^{-ur} = \sqrt{\pi}/(2r^{3/2}) ∫ 0 ∞ d u u e − u r = π / ( 2 r 3/2 ) を使うと
∣ U ∣ ∼ 2 π m 3 / 2 t 4 π 2 r 5 / 2 e − m r \lvert U\rvert \sim \frac{\sqrt{2\pi}\;m^{3/2}\,t}{4\pi^2\,r^{5/2}}\,e^{-mr} ∣ U ∣ ∼ 4 π 2 r 5/2 2 π m 3/2 t e − m r です。より丁寧な鞍点評価では指数がローレンツ不変な形 e − m r 2 − t 2 e^{-m\sqrt{r^2-t^2}} e − m r 2 − t 2 になります。振幅はコンプトン波長 1 / m 1/m 1/ m を超える距離で指数関数的に小さいものの、決してゼロではありません。
原点の事象から空間的に隔たった点 P へも、一粒子の伝播振幅は指数関数的に小さいだけでゼロにならない。灰色は空間的領域。
例 5.2 は具体的な計算でしたが、実はこれは H ≥ 0 H \ge 0 H ≥ 0 という条件だけから従う一般的な事実です。次の定理はヘーゲルフェルト(1974 年)によるものです。
定理 5.3 (ヘーゲルフェルトの定理 )
H \mathcal{H} H を Hilbert 空間、H H H を H \mathcal{H} H 上の自己共役作用素で H ≥ 0 H \ge 0 H ≥ 0 (スペクトルが [ 0 , ∞ ) [0,\infty) [ 0 , ∞ ) に含まれる)とする。ψ 0 ∈ H \psi_0\in\mathcal{H} ψ 0 ∈ H 、ψ t = e − i H t ψ 0 \psi_t = e^{-iHt}\psi_0 ψ t = e − i H t ψ 0 とし、A A A を有界な正作用素(A ≥ 0 A \ge 0 A ≥ 0 )とする。このとき、ある開区間 I ⊂ R I\subset\mathbb{R} I ⊂ R 上で
⟨ ψ t , A ψ t ⟩ = 0 ( ∀ t ∈ I ) \langle \psi_t, A\psi_t\rangle = 0 \quad (\forall t\in I) ⟨ ψ t , A ψ t ⟩ = 0 ( ∀ t ∈ I ) が成り立つならば、すべての t ∈ R t\in\mathbb{R} t ∈ R で ⟨ ψ t , A ψ t ⟩ = 0 \langle\psi_t,A\psi_t\rangle = 0 ⟨ ψ t , A ψ t ⟩ = 0 である。
証明(定理 5.3) A ≥ 0 A\ge 0 A ≥ 0 かつ有界なので、有界な平方根 A 1 / 2 ≥ 0 A^{1/2}\ge 0 A 1/2 ≥ 0 が存在し ⟨ ψ t , A ψ t ⟩ = ∥ A 1 / 2 ψ t ∥ 2 \langle\psi_t,A\psi_t\rangle = \lVert A^{1/2}\psi_t\rVert^2 ⟨ ψ t , A ψ t ⟩ = ∥ A 1/2 ψ t ∥ 2 です。したがって仮定は「t ∈ I t\in I t ∈ I に対し A 1 / 2 ψ t = 0 A^{1/2}\psi_t = 0 A 1/2 ψ t = 0 」と同値です。
任意に χ ∈ H \chi\in\mathcal{H} χ ∈ H を固定し、下半平面 C − = { z : Im z < 0 } \mathbb{C}_- = \{z : \operatorname{Im}z < 0\} C − = { z : Im z < 0 } 上で
g ( z ) = ⟨ χ , A 1 / 2 e − i z H ψ 0 ⟩ g(z) = \left\langle \chi,\; A^{1/2}e^{-izH}\psi_0\right\rangle g ( z ) = ⟨ χ , A 1/2 e − i z H ψ 0 ⟩ を考えます。H H H のスペクトル測度を E ( ⋅ ) E(\cdot) E ( ⋅ ) とし、複素測度 μ ( ⋅ ) = ⟨ χ , A 1 / 2 E ( ⋅ ) ψ 0 ⟩ \mu(\cdot) = \langle \chi, A^{1/2}E(\cdot)\psi_0\rangle μ ( ⋅ ) = ⟨ χ , A 1/2 E ( ⋅ ) ψ 0 ⟩ (全変動は ∥ χ ∥ ∥ A 1 / 2 ∥ ∥ ψ 0 ∥ \lVert\chi\rVert\,\lVert A^{1/2}\rVert\,\lVert\psi_0\rVert ∥ χ ∥ ∥ A 1/2 ∥ ∥ ψ 0 ∥ 以下の有限値)を使うと
g ( z ) = ∫ [ 0 , ∞ ) e − i z λ d μ ( λ ) g(z) = \int_{[0,\infty)} e^{-iz\lambda}\,d\mu(\lambda) g ( z ) = ∫ [ 0 , ∞ ) e − i z λ d μ ( λ ) と書けます。H ≥ 0 H\ge 0 H ≥ 0 よりスペクトルは λ ≥ 0 \lambda\ge 0 λ ≥ 0 に限られ、z = t − i s z = t-is z = t − i s (s > 0 s>0 s > 0 )に対して ∣ e − i z λ ∣ = e − s λ ≤ 1 \lvert e^{-iz\lambda}\rvert = e^{-s\lambda}\le 1 ∣ e − i z λ ∣ = e − s λ ≤ 1 です。よって g g g は C − \mathbb{C}_- C − 上で有界であり、被積分関数が正則かつ一様可積分な優関数をもつことから、Fubini と Morera の定理により正則です。優収束定理により、境界値 g ( t ) = ⟨ χ , A 1 / 2 ψ t ⟩ g(t) = \langle\chi, A^{1/2}\psi_t\rangle g ( t ) = ⟨ χ , A 1/2 ψ t ⟩ をもちます。
仮定より境界値は I I I 上で消えます。I I I は正の Lebesgue 測度をもつ集合です。半平面上の有界正則関数(Hardy 空間 H ∞ H^\infty H ∞ の元)で、境界値が正測度の集合上で消えるものは恒等的に 0 0 0 である、という境界一意性定理(F. and M. Riesz–Privalov の定理)により g ≡ 0 g\equiv 0 g ≡ 0 を得ます。
χ \chi χ は任意でしたから、すべての実 t t t で A 1 / 2 ψ t = 0 A^{1/2}\psi_t = 0 A 1/2 ψ t = 0 、すなわち ⟨ ψ t , A ψ t ⟩ = ∥ A 1 / 2 ψ t ∥ 2 = 0 \langle\psi_t,A\psi_t\rangle = \lVert A^{1/2}\psi_t\rVert^2 = 0 ⟨ ψ t , A ψ t ⟩ = ∥ A 1/2 ψ t ∥ 2 = 0 です。
∎
系 5.4 (厳密な局在と因果的伝播は両立しない )
H ≥ 0 H\ge 0 H ≥ 0 とし、有界開集合 W ⊂ R 3 W\subset\mathbb{R}^3 W ⊂ R 3 に対して「粒子が W W W の中に見出される」ことを表す射影 P W P_W P W (W ⊂ W ′ W\subset W' W ⊂ W ′ ならば P W ≤ P W ′ P_W \le P_{W'} P W ≤ P W ′ )が与えられているとする。状態 ψ 0 \psi_0 ψ 0 が有界領域 V V V に厳密に局在し、かつ伝播が因果的である、すなわち
dist ( W , V ) > ∣ t ∣ ⟹ ⟨ ψ t , P W ψ t ⟩ = 0 \operatorname{dist}(W,V) > \lvert t\rvert \;\Longrightarrow\; \langle\psi_t, P_W\psi_t\rangle = 0 dist ( W , V ) > ∣ t ∣ ⟹ ⟨ ψ t , P W ψ t ⟩ = 0 が成り立つとする(c = 1 c=1 c = 1 )。このとき、V V V から正の距離にある任意の有界開集合 W W W について、すべての 時刻 t ∈ R t\in\mathbb{R} t ∈ R で ⟨ ψ t , P W ψ t ⟩ = 0 \langle\psi_t,P_W\psi_t\rangle = 0 ⟨ ψ t , P W ψ t ⟩ = 0 である。
証明(系 5.4) d = dist ( W , V ) > 0 d = \operatorname{dist}(W,V) > 0 d = dist ( W , V ) > 0 とおきます。因果的伝播の仮定より、開区間 I = ( − d , d ) I = (-d,d) I = ( − d , d ) 上で ⟨ ψ t , P W ψ t ⟩ = 0 \langle\psi_t,P_W\psi_t\rangle = 0 ⟨ ψ t , P W ψ t ⟩ = 0 です。射影 P W P_W P W は有界かつ P W = P W 2 = P W † ≥ 0 P_W = P_W^2 = P_W^\dagger \ge 0 P W = P W 2 = P W † ≥ 0 なので、定理 5.3 を A = P W A = P_W A = P W に適用でき、結論はすべての t t t に拡張されます。
∎
つまり粒子は永遠に V ˉ \bar V V ˉ の外側では見出されないことになります。ところが H = − ∇ 2 + m 2 H=\sqrt{-\nabla^2+m^2} H = − ∇ 2 + m 2 のスペクトルは絶対連続なので、RAGE の定理により任意の有界領域 K K K で ∥ P K ψ t ∥ → 0 \lVert P_K\psi_t\rVert \to 0 ∥ P K ψ t ∥ → 0 (t → ∞ t\to\infty t → ∞ )、粒子は必ず広がります。矛盾です。
結論はこうです。「H ≥ 0 H\ge 0 H ≥ 0 」「厳密な局在の概念」「因果的伝播」の三つは同時には成り立ちません。 H ≥ 0 H\ge 0 H ≥ 0 を捨てれば 定理 4.2 の困難に戻り、因果律を捨てれば相対論を捨てることになり、局在を捨てれば「粒子がどこにあるか」という問い自体が意味を失います。
ここまでは論理的な破綻でした。物理的にも、一粒子理論が成り立たない理由ははっきりしています。相対論では質量とエネルギーが交換可能であり(相対論的力学 )、十分なエネルギーがあれば粒子は作られてしまうのです。
定義 6.1 (換算コンプトン波長 )
質量 m m m の粒子に対して
λ ˉ C = ℏ m c \bar\lambda_C = \frac{\hbar}{mc} λ ˉ C = m c ℏ を換算コンプトン波長と呼びます(2 π 2\pi 2 π を掛けたものがコンプトン波長 λ C = h / m c \lambda_C = h/mc λ C = h / m c です)。
このスケールの意味は、ハイゼンベルクの不確定性原理(系 5.4)[演算子と物理量] から読めます。粒子を大きさ Δ x \Delta x Δ x の領域に閉じ込めると Δ p ≳ ℏ / Δ x \Delta p \gtrsim \hbar/\Delta x Δ p ≳ ℏ/Δ x の運動量の不定性が生じ、相対論的にはこれが
Δ E ≳ c Δ p ≳ ℏ c Δ x \Delta E \gtrsim c\,\Delta p \gtrsim \frac{\hbar c}{\Delta x} Δ E ≳ c Δ p ≳ Δ x ℏ c
のエネルギーの不定性を意味します。Δ x ≲ λ ˉ C \Delta x \lesssim \bar\lambda_C Δ x ≲ λ ˉ C になると Δ E ≳ m c 2 \Delta E \gtrsim mc^2 Δ E ≳ m c 2 となり、粒子・反粒子対を作るのに十分なエネルギーが「借りられて」しまいます。
例 6.2 (電子の場合の数値 )
ℏ = 1.055 × 10 − 34 J s \hbar = 1.055\times10^{-34}\ \mathrm{J\,s} ℏ = 1.055 × 1 0 − 34 J s 、m e = 9.109 × 10 − 31 k g m_e = 9.109\times10^{-31}\ \mathrm{kg} m e = 9.109 × 1 0 − 31 kg 、c = 2.998 × 10 8 m / s c = 2.998\times10^{8}\ \mathrm{m/s} c = 2.998 × 1 0 8 m/s を使います。m e c = 2.731 × 10 − 22 k g m / s m_e c = 2.731\times10^{-22}\ \mathrm{kg\,m/s} m e c = 2.731 × 1 0 − 22 kg m/s なので
λ ˉ C = 1.055 × 10 − 34 2.731 × 10 − 22 m = 3.86 × 10 − 13 m = 386 f m \bar\lambda_C = \frac{1.055\times10^{-34}}{2.731\times10^{-22}}\ \mathrm{m} = 3.86\times10^{-13}\ \mathrm{m} = 386\ \mathrm{fm} λ ˉ C = 2.731 × 1 0 − 22 1.055 × 1 0 − 34 m = 3.86 × 1 0 − 13 m = 386 fm です。一方、対生成のしきい値は 2 m e c 2 = 1.022 M e V 2m_ec^2 = 1.022\ \mathrm{MeV} 2 m e c 2 = 1.022 MeV です。つまり電子を 400 f m 400\ \mathrm{fm} 400 fm 程度より狭い領域に押し込もうとすると、そのために注ぎ込むエネルギーで新しい電子・陽電子対ができてしまいます。「1 個の電子の波動関数」を λ ˉ C \bar\lambda_C λ ˉ C より細かい分解能で語ることに意味がなくなるのです。
これは 例 5.2 で見た振幅の減衰長 1 / m 1/m 1/ m (自然単位系)と同じスケールです。因果律を破る振幅が効く距離と、対生成が始まる距離が一致するのは偶然ではありません。どちらも「一粒子描像が壊れる」という同じ事実の別の顔です。
ここまでの三つの破綻(負エネルギー、確率密度、因果律)と粒子数非保存を、一挙に解く処方があります。古典場を力学変数として正準量子化する ことです。
出発点は ラグランジュ形式の力学 を場に拡張したものです(詳しくは 古典場の理論とラグランジアン で扱います)。実スカラー場 ϕ ( t , x ) \phi(t,\boldsymbol{x}) ϕ ( t , x ) の最も簡単な Lagrangian は
L = ∫ d 3 x L , L = 1 2 ϕ ˙ 2 − 1 2 ( ∇ ϕ ) 2 − 1 2 m 2 ϕ 2 L = \int d^3x\;\mathcal{L},
\qquad
\mathcal{L} = \tfrac12\dot\phi^2 - \tfrac12(\nabla\phi)^2 - \tfrac12 m^2\phi^2 L = ∫ d 3 x L , L = 2 1 ϕ ˙ 2 − 2 1 ( ∇ ϕ ) 2 − 2 1 m 2 ϕ 2
で、Euler–Lagrange 方程式はまさに 定義 3.1 のクライン・ゴルドン方程式です(古典場としての導出と分散関係(例 5.1)[古典場の理論とラグランジアン] )。ここが決定的な視点の転換点です。クライン・ゴルドン方程式を「一粒子の波動関数の方程式」ではなく「古典場の運動方程式 」として読み直すのです。
定理 7.1 (自由スカラー場のモード分解と量子化 )
一辺 L L L の立方体(体積 V = L 3 V = L^3 V = L 3 、周期境界条件)の中の実スカラー場 ϕ \phi ϕ が上の Lagrangian に従うとする。Fourier 展開
ϕ ( t , x ) = 1 V ∑ k q k ( t ) e i k ⋅ x , k ∈ 2 π L Z 3 , q − k = q k ‾ \phi(t,\boldsymbol{x}) = \frac{1}{\sqrt{V}}\sum_{\boldsymbol{k}} q_{\boldsymbol{k}}(t)\,e^{i\boldsymbol{k}\cdot\boldsymbol{x}},
\qquad
\boldsymbol{k}\in\frac{2\pi}{L}\mathbb{Z}^3,
\qquad
q_{-\boldsymbol{k}} = \overline{q_{\boldsymbol{k}}} ϕ ( t , x ) = V 1 k ∑ q k ( t ) e i k ⋅ x , k ∈ L 2 π Z 3 , q − k = q k を代入すると
L = ∑ k 1 2 ( ∣ q ˙ k ∣ 2 − ω k 2 ∣ q k ∣ 2 ) , ω k = k 2 + m 2 L = \sum_{\boldsymbol{k}}\frac12\left(\lvert\dot q_{\boldsymbol{k}}\rvert^2 - \omega_{\boldsymbol{k}}^2\,\lvert q_{\boldsymbol{k}}\rvert^2\right),
\qquad
\omega_{\boldsymbol{k}} = \sqrt{\boldsymbol{k}^2+m^2} L = k ∑ 2 1 ( ∣ q ˙ k ∣ 2 − ω k 2 ∣ q k ∣ 2 ) , ω k = k 2 + m 2 となる。すなわち自由場は、モードごとに角振動数 ω k \omega_{\boldsymbol{k}} ω k をもつ独立な調和振動子の集まりと等価である。これを正準量子化すると、[ a ^ k , a ^ k ′ † ] = δ k k ′ [\hat a_{\boldsymbol k},\hat a^\dagger_{\boldsymbol k'}] = \delta_{\boldsymbol{k}\boldsymbol{k}'} [ a ^ k , a ^ k ′ † ] = δ k k ′ 、他の交換子はゼロ、として
H ^ = ∑ k ω k ( a ^ k † a ^ k + 1 2 ) , P ^ = ∑ k k a ^ k † a ^ k \hat H = \sum_{\boldsymbol{k}}\omega_{\boldsymbol{k}}\left(\hat a^\dagger_{\boldsymbol{k}}\hat a_{\boldsymbol{k}} + \tfrac12\right),
\qquad
\hat{\boldsymbol{P}} = \sum_{\boldsymbol{k}}\boldsymbol{k}\,\hat a^\dagger_{\boldsymbol{k}}\hat a_{\boldsymbol{k}} H ^ = k ∑ ω k ( a ^ k † a ^ k + 2 1 ) , P ^ = k ∑ k a ^ k † a ^ k が得られ、エネルギー固有値は(零点エネルギーを除いて)∑ k n k ω k \sum_{\boldsymbol{k}} n_{\boldsymbol{k}}\,\omega_{\boldsymbol{k}} ∑ k n k ω k 、n k ∈ { 0 , 1 , 2 , … } n_{\boldsymbol{k}}\in\{0,1,2,\ldots\} n k ∈ { 0 , 1 , 2 , … } である。
証明(定理 7.1) 周期境界条件のもとで 1 V ∫ V d 3 x e i ( k + k ′ ) ⋅ x = δ k , − k ′ \frac{1}{V}\int_V d^3x\;e^{i(\boldsymbol{k}+\boldsymbol{k}')\cdot\boldsymbol{x}} = \delta_{\boldsymbol{k},-\boldsymbol{k}'} V 1 ∫ V d 3 x e i ( k + k ′ ) ⋅ x = δ k , − k ′ が成り立つことを繰り返し使います。まず
∫ d 3 x ϕ ˙ 2 = ∑ k , k ′ q ˙ k q ˙ k ′ δ k , − k ′ = ∑ k q ˙ k q ˙ − k = ∑ k ∣ q ˙ k ∣ 2 \int d^3x\;\dot\phi^2 = \sum_{\boldsymbol{k},\boldsymbol{k}'}\dot q_{\boldsymbol{k}}\dot q_{\boldsymbol{k}'}\,\delta_{\boldsymbol{k},-\boldsymbol{k}'} = \sum_{\boldsymbol{k}}\dot q_{\boldsymbol{k}}\dot q_{-\boldsymbol{k}} = \sum_{\boldsymbol{k}}\lvert\dot q_{\boldsymbol{k}}\rvert^2 ∫ d 3 x ϕ ˙ 2 = k , k ′ ∑ q ˙ k q ˙ k ′ δ k , − k ′ = k ∑ q ˙ k q ˙ − k = k ∑ ∣ q ˙ k ∣ 2 です(最後で実数条件 q − k = q k ‾ q_{-\boldsymbol{k}}=\overline{q_{\boldsymbol{k}}} q − k = q k を使いました)。同様に ∇ e i k ⋅ x = i k e i k ⋅ x \nabla e^{i\boldsymbol{k}\cdot\boldsymbol{x}} = i\boldsymbol{k}\,e^{i\boldsymbol{k}\cdot\boldsymbol{x}} ∇ e i k ⋅ x = i k e i k ⋅ x より
∫ d 3 x ( ∇ ϕ ) 2 = ∑ k ( i k q k ) ⋅ ( − i k q − k ) = ∑ k k 2 ∣ q k ∣ 2 , ∫ d 3 x ϕ 2 = ∑ k ∣ q k ∣ 2 . \int d^3x\;(\nabla\phi)^2 = \sum_{\boldsymbol{k}}(i\boldsymbol{k}q_{\boldsymbol{k}})\cdot(-i\boldsymbol{k}q_{-\boldsymbol{k}}) = \sum_{\boldsymbol{k}}\boldsymbol{k}^2\lvert q_{\boldsymbol{k}}\rvert^2,
\qquad
\int d^3x\;\phi^2 = \sum_{\boldsymbol{k}}\lvert q_{\boldsymbol{k}}\rvert^2 . ∫ d 3 x ( ∇ ϕ ) 2 = k ∑ ( i k q k ) ⋅ ( − i k q − k ) = k ∑ k 2 ∣ q k ∣ 2 , ∫ d 3 x ϕ 2 = k ∑ ∣ q k ∣ 2 . 三つを合わせると主張の L L L の表式になります。k 2 + m 2 = ω k 2 \boldsymbol{k}^2 + m^2 = \omega_{\boldsymbol{k}}^2 k 2 + m 2 = ω k 2 です。
次に、複素変数 q k q_{\boldsymbol{k}} q k は独立ではないので実の自由度に直します。k \boldsymbol{k} k と − k -\boldsymbol{k} − k の一方だけを含む半分の集合 K + K_+ K + を取り、k ∈ K + \boldsymbol{k}\in K_+ k ∈ K + に対して q k = ( x k + i y k ) / 2 q_{\boldsymbol{k}} = (x_{\boldsymbol{k}} + i y_{\boldsymbol{k}})/\sqrt2 q k = ( x k + i y k ) / 2 と書きます。∣ q − k ∣ = ∣ q k ∣ \lvert q_{-\boldsymbol{k}}\rvert = \lvert q_{\boldsymbol{k}}\rvert ∣ q − k ∣ = ∣ q k ∣ なので ± k \pm\boldsymbol{k} ± k の寄与は等しく、
L = ∑ k ∈ K + 1 2 [ ( x ˙ k 2 + y ˙ k 2 ) − ω k 2 ( x k 2 + y k 2 ) ] + 1 2 ( q ˙ 0 2 − m 2 q 0 2 ) L = \sum_{\boldsymbol{k}\in K_+}\frac12\left[(\dot x_{\boldsymbol{k}}^2 + \dot y_{\boldsymbol{k}}^2) - \omega_{\boldsymbol{k}}^2(x_{\boldsymbol{k}}^2+y_{\boldsymbol{k}}^2)\right] + \frac12\left(\dot q_{\boldsymbol{0}}^2 - m^2 q_{\boldsymbol{0}}^2\right) L = k ∈ K + ∑ 2 1 [ ( x ˙ k 2 + y ˙ k 2 ) − ω k 2 ( x k 2 + y k 2 ) ] + 2 1 ( q ˙ 0 2 − m 2 q 0 2 ) となります。これは質量 1 1 1 、角振動数 ω k \omega_{\boldsymbol{k}} ω k の実の調和振動子が k \boldsymbol{k} k 1 個につき 1 個(K + K_+ K + の各元が 2 個ずつ、k = 0 \boldsymbol{k}=\boldsymbol{0} k = 0 が 1 個)並んだものです。
各振動子を ハミルトン形式の力学 に移して [ q ^ , p ^ ] = i [\hat q,\hat p]=i [ q ^ , p ^ ] = i を課し、標準的な生成消滅演算子を導入すると H ^ osc = ω ( a ^ † a ^ + 1 2 ) \hat H_{\text{osc}} = \omega(\hat a^\dagger\hat a + \tfrac12) H ^ osc = ω ( a ^ † a ^ + 2 1 ) を得ます(調和振動子のスペクトル(定理 5.3)[1次元の簡単な系] )。x k , y k x_{\boldsymbol{k}},y_{\boldsymbol{k}} x k , y k の組を複素の組み合わせに戻すと、a ^ k \hat a_{\boldsymbol{k}} a ^ k と a ^ − k \hat a_{-\boldsymbol{k}} a ^ − k という運動量 ± k \pm\boldsymbol{k} ± k のラベルをもつ演算子になり、全ハミルトニアンは主張の形にまとまります。運動量については、Noether の定理(対称性と保存則 )から得られる P = − ∫ d 3 x ϕ ˙ ∇ ϕ \boldsymbol{P} = -\int d^3x\,\dot\phi\,\nabla\phi P = − ∫ d 3 x ϕ ˙ ∇ ϕ を同じ手順で書き直すと ∑ k k a ^ k † a ^ k \sum_{\boldsymbol{k}}\boldsymbol{k}\,\hat a^\dagger_{\boldsymbol{k}}\hat a_{\boldsymbol{k}} ∑ k k a ^ k † a ^ k になります(零点項は k → − k \boldsymbol{k}\to-\boldsymbol{k} k → − k の対称性で相殺します)。
固有状態は各モードの数状態の積で、固有値は ∑ k ω k ( n k + 1 2 ) \sum_{\boldsymbol{k}}\omega_{\boldsymbol{k}}(n_{\boldsymbol{k}}+\tfrac12) ∑ k ω k ( n k + 2 1 ) です。
∎
この結果の読み方が、場の量子論そのものです。a ^ k † \hat a^\dagger_{\boldsymbol{k}} a ^ k † を真空 ∣ 0 ⟩ \lvert 0\rangle ∣ 0 ⟩ に 1 回作用させた状態 a ^ k † ∣ 0 ⟩ \hat a^\dagger_{\boldsymbol{k}}\lvert 0\rangle a ^ k † ∣ 0 ⟩ は、定理 7.1 より
E = ω k = k 2 + m 2 , P = k E = \omega_{\boldsymbol{k}} = \sqrt{\boldsymbol{k}^2+m^2},
\qquad
\boldsymbol{P} = \boldsymbol{k} E = ω k = k 2 + m 2 , P = k
をもちます。これは質量 m m m の相対論的な粒子 1 個のエネルギー・運動量そのものです。n n n 回作用させれば n n n 個の粒子です。粒子とは場の励起量子である という描像が、計算の帰結として出てきました。
定義 7.2 (Fock 空間 )
真空 ∣ 0 ⟩ \lvert 0\rangle ∣ 0 ⟩ (すべての k \boldsymbol{k} k で a ^ k ∣ 0 ⟩ = 0 \hat a_{\boldsymbol{k}}\lvert 0\rangle = 0 a ^ k ∣ 0 ⟩ = 0 )に生成演算子を作用させて得られる状態全体が張る空間
F = ⨁ n = 0 ∞ H ( n ) , H ( n ) = span { a ^ k 1 † ⋯ a ^ k n † ∣ 0 ⟩ } \mathcal{F} = \bigoplus_{n=0}^{\infty}\mathcal{H}^{(n)},
\qquad
\mathcal{H}^{(n)} = \operatorname{span}\left\{\hat a^\dagger_{\boldsymbol{k}_1}\cdots \hat a^\dagger_{\boldsymbol{k}_n}\lvert 0\rangle\right\} F = n = 0 ⨁ ∞ H ( n ) , H ( n ) = span { a ^ k 1 † ⋯ a ^ k n † ∣ 0 ⟩ } を Fock 空間と呼びます。粒子数演算子は N ^ = ∑ k a ^ k † a ^ k \hat N = \sum_{\boldsymbol{k}}\hat a^\dagger_{\boldsymbol{k}}\hat a_{\boldsymbol{k}} N ^ = ∑ k a ^ k † a ^ k です。
Fock 空間は粒子数の異なる部分空間の直和なので、粒子数を変える演算子が書けます 。(A1) の制約が外れました。しかも状態は占有数 { n k } \{n_{\boldsymbol{k}}\} { n k } だけで決まるので、粒子の入れ替えに対する対称性は自動的に成り立ちます。同種粒子の対称化を手で課す必要がありません。
例 7.3 (相互作用は粒子数を保存しない )
Lagrangian に − λ 4 ! ϕ 4 -\frac{\lambda}{4!}\phi^4 − 4 ! λ ϕ 4 を加えます。対応する相互作用ハミルトニアンは H ^ int = λ 4 ! ∫ d 3 x ϕ ^ 4 \hat H_{\text{int}} = \frac{\lambda}{4!}\int d^3x\;\hat\phi^4 H ^ int = 4 ! λ ∫ d 3 x ϕ ^ 4 です。場の演算子は生成部分と消滅部分の和 ϕ ^ ∼ a ^ + a ^ † \hat\phi \sim \hat a + \hat a^\dagger ϕ ^ ∼ a ^ + a ^ † ですから、ϕ ^ 4 \hat\phi^4 ϕ ^ 4 を展開すると 2 4 = 16 2^4 = 16 2 4 = 16 個の項が出て、含まれる生成演算子の個数に応じて分類できます。
ϕ ^ 4 ∼ a ^ † a ^ † a ^ † a ^ † ⏟ Δ N = + 4 + a ^ † a ^ † a ^ † a ^ ⏟ Δ N = + 2 + a ^ † a ^ † a ^ a ^ ⏟ Δ N = 0 + a ^ † a ^ a ^ a ^ ⏟ Δ N = − 2 + a ^ a ^ a ^ a ^ ⏟ Δ N = − 4 \hat\phi^4 \sim \underbrace{\hat a^\dagger\hat a^\dagger\hat a^\dagger\hat a^\dagger}_{\Delta N = +4} + \underbrace{\hat a^\dagger\hat a^\dagger\hat a^\dagger\hat a}_{\Delta N=+2} + \underbrace{\hat a^\dagger\hat a^\dagger\hat a\hat a}_{\Delta N = 0} + \underbrace{\hat a^\dagger\hat a\hat a\hat a}_{\Delta N=-2} + \underbrace{\hat a\hat a\hat a\hat a}_{\Delta N=-4} ϕ ^ 4 ∼ Δ N = + 4 a ^ † a ^ † a ^ † a ^ † + Δ N = + 2 a ^ † a ^ † a ^ † a ^ + Δ N = 0 a ^ † a ^ † a ^ a ^ + Δ N = − 2 a ^ † a ^ a ^ a ^ + Δ N = − 4 a ^ a ^ a ^ a ^ (各括弧は運動量の並び替えを含むすべての組み合わせを表します。)Δ N = ± 4 , ± 2 \Delta N = \pm 4, \pm 2 Δ N = ± 4 , ± 2 の項があるので [ H ^ , N ^ ] ≠ 0 [\hat H,\hat N]\ne 0 [ H ^ , N ^ ] = 0 、粒子数は保存しません。真空から 4 個の粒子が生まれる過程も、2 個が 2 個に散乱する過程も、同じ 1 個の項から出てきます。定義 7.2 の枠組みで初めて、こうした過程が記述できるようになりました。
定理 7.1 の H ^ \hat H H ^ から零点エネルギー ∑ k ω k / 2 \sum_{\boldsymbol{k}}\omega_{\boldsymbol{k}}/2 ∑ k ω k /2 (定数)を差し引くと H ^ = ∑ k ω k a ^ k † a ^ k ≥ 0 \hat H = \sum_{\boldsymbol{k}}\omega_{\boldsymbol{k}}\hat a^\dagger_{\boldsymbol{k}}\hat a_{\boldsymbol{k}} \ge 0 H ^ = ∑ k ω k a ^ k † a ^ k ≥ 0 です。負エネルギーはどこへ行ったのでしょうか。答えは「解の負振動数部分 e + i ω t e^{+i\omega t} e + iω t は、負エネルギー状態ではなく生成演算子 の係数になった」です。モード展開では
ϕ ^ ( x ) ∼ a ^ k e − i ω k t ⏟ 正振動数 + a ^ k † e + i ω k t ⏟ 負振動数 \hat\phi(x) \;\sim\; \underbrace{\hat a_{\boldsymbol{k}}\,e^{-i\omega_{\boldsymbol{k}}t}}_{\text{正振動数}} \;+\; \underbrace{\hat a^\dagger_{\boldsymbol{k}}\,e^{+i\omega_{\boldsymbol{k}}t}}_{\text{負振動数}} ϕ ^ ( x ) ∼ 正振動数 a ^ k e − i ω k t + 負振動数 a ^ k † e + i ω k t
という形で両方が現れますが、a ^ k † \hat a^\dagger_{\boldsymbol{k}} a ^ k † はエネルギーを + ω k +\omega_{\boldsymbol{k}} + ω k だけ増やす 演算子です。同じ数式が、解釈を変えるだけで下に有界なスペクトルを与えます。
複素スカラー場では 命題 3.2 の j μ j^\mu j μ が電荷カレントになり、保存電荷は
Q ^ = ∑ k ( a ^ k † a ^ k − b ^ k † b ^ k ) \hat Q = \sum_{\boldsymbol{k}}\left(\hat a^\dagger_{\boldsymbol{k}}\hat a_{\boldsymbol{k}} - \hat b^\dagger_{\boldsymbol{k}}\hat b_{\boldsymbol{k}}\right) Q ^ = k ∑ ( a ^ k † a ^ k − b ^ k † b ^ k )
の形になります。a ^ † \hat a^\dagger a ^ † が作る粒子と b ^ † \hat b^\dagger b ^ † が作る粒子は、質量が同じで電荷の符号が逆です。後者が反粒子 です。j 0 < 0 j^0 < 0 j 0 < 0 は「負の確率」ではなく「反粒子の存在」だったのです(注意 3.4 で予告したとおりです)。
最後に残った因果律を見ます。場の量子論では、因果律は「振幅がゼロになること」ではなく「空間的に隔たった二点での場の演算子が可換であること 」として実現されます。可換なら、一方の点での測定がもう一方の点での測定結果に影響しません。
定義 7.4 (自由実スカラー場の演算子 )
連続体極限(V → ∞ V\to\infty V → ∞ )では
ϕ ^ ( x ) = ∫ d 3 p ( 2 π ) 3 1 2 ω p ( a ^ p e − i p ⋅ x + a ^ p † e + i p ⋅ x ) , p 0 = ω p \hat\phi(x) = \int\frac{d^3p}{(2\pi)^3}\;\frac{1}{\sqrt{2\omega_{\boldsymbol{p}}}}\left(\hat a_{\boldsymbol{p}}\,e^{-ip\cdot x} + \hat a^\dagger_{\boldsymbol{p}}\,e^{+ip\cdot x}\right),
\qquad p^0 = \omega_{\boldsymbol{p}} ϕ ^ ( x ) = ∫ ( 2 π ) 3 d 3 p 2 ω p 1 ( a ^ p e − i p ⋅ x + a ^ p † e + i p ⋅ x ) , p 0 = ω p と書きます。交換関係は [ a ^ p , a ^ q † ] = ( 2 π ) 3 δ 3 ( p − q ) [\hat a_{\boldsymbol{p}},\hat a^\dagger_{\boldsymbol{q}}] = (2\pi)^3\delta^3(\boldsymbol{p}-\boldsymbol{q}) [ a ^ p , a ^ q † ] = ( 2 π ) 3 δ 3 ( p − q ) 、他はゼロです。
定理 7.5 (ミクロ因果律 )
定義 7.4 の場について、D ( x ) : = ∫ d 3 p ( 2 π ) 3 2 ω p e − i p ⋅ x D(x) := \displaystyle\int\frac{d^3p}{(2\pi)^3\,2\omega_{\boldsymbol{p}}}e^{-ip\cdot x} D ( x ) := ∫ ( 2 π ) 3 2 ω p d 3 p e − i p ⋅ x とおくと
[ ϕ ^ ( x ) , ϕ ^ ( y ) ] = D ( x − y ) − D ( y − x ) [\hat\phi(x),\hat\phi(y)] = D(x-y) - D(y-x) [ ϕ ^ ( x ) , ϕ ^ ( y )] = D ( x − y ) − D ( y − x ) であり、( x − y ) 2 < 0 (x-y)^2 < 0 ( x − y ) 2 < 0 (空間的隔たり)ならばこれは 0 0 0 である。
証明(定理 7.5) 第 1 段:交換子の表式。 定義 7.4 を代入すると、生き残るのは [ a ^ p , a ^ q † ] [\hat a_{\boldsymbol{p}},\hat a^\dagger_{\boldsymbol{q}}] [ a ^ p , a ^ q † ] と [ a ^ p † , a ^ q ] = − ( 2 π ) 3 δ 3 ( p − q ) [\hat a^\dagger_{\boldsymbol{p}},\hat a_{\boldsymbol{q}}] = -(2\pi)^3\delta^3(\boldsymbol{p}-\boldsymbol{q}) [ a ^ p † , a ^ q ] = − ( 2 π ) 3 δ 3 ( p − q ) を含む項だけです。
[ ϕ ^ ( x ) , ϕ ^ ( y ) ] = ∫ d 3 p d 3 q ( 2 π ) 6 4 ω p ω q ( [ a ^ p , a ^ q † ] e − i p ⋅ x + i q ⋅ y + [ a ^ p † , a ^ q ] e i p ⋅ x − i q ⋅ y ) [\hat\phi(x),\hat\phi(y)] = \int\frac{d^3p\,d^3q}{(2\pi)^6\sqrt{4\omega_{\boldsymbol{p}}\omega_{\boldsymbol{q}}}}
\left([\hat a_{\boldsymbol{p}},\hat a^\dagger_{\boldsymbol{q}}]e^{-ip\cdot x + iq\cdot y} + [\hat a^\dagger_{\boldsymbol{p}},\hat a_{\boldsymbol{q}}]e^{ip\cdot x - iq\cdot y}\right) [ ϕ ^ ( x ) , ϕ ^ ( y )] = ∫ ( 2 π ) 6 4 ω p ω q d 3 p d 3 q ( [ a ^ p , a ^ q † ] e − i p ⋅ x + i q ⋅ y + [ a ^ p † , a ^ q ] e i p ⋅ x − i q ⋅ y ) = ∫ d 3 p ( 2 π ) 3 2 ω p ( e − i p ⋅ ( x − y ) − e i p ⋅ ( x − y ) ) = D ( x − y ) − D ( y − x ) . = \int\frac{d^3p}{(2\pi)^3\,2\omega_{\boldsymbol{p}}}\left(e^{-ip\cdot(x-y)} - e^{ip\cdot(x-y)}\right) = D(x-y) - D(y-x). = ∫ ( 2 π ) 3 2 ω p d 3 p ( e − i p ⋅ ( x − y ) − e i p ⋅ ( x − y ) ) = D ( x − y ) − D ( y − x ) . 第 2 段:測度のローレンツ不変性。 δ ( p 2 − m 2 ) = δ ( ( p 0 ) 2 − ω p 2 ) \delta(p^2-m^2) = \delta\!\left((p^0)^2 - \omega_{\boldsymbol{p}}^2\right) δ ( p 2 − m 2 ) = δ ( ( p 0 ) 2 − ω p 2 ) の零点は p 0 = ± ω p p^0 = \pm\omega_{\boldsymbol{p}} p 0 = ± ω p で、そこでの微分の絶対値は 2 ω p 2\omega_{\boldsymbol{p}} 2 ω p ですから
δ ( p 2 − m 2 ) = 1 2 ω p [ δ ( p 0 − ω p ) + δ ( p 0 + ω p ) ] . \delta(p^2-m^2) = \frac{1}{2\omega_{\boldsymbol{p}}}\left[\delta(p^0-\omega_{\boldsymbol{p}}) + \delta(p^0+\omega_{\boldsymbol{p}})\right]. δ ( p 2 − m 2 ) = 2 ω p 1 [ δ ( p 0 − ω p ) + δ ( p 0 + ω p ) ] . したがって任意の関数 f f f について
∫ d 4 p ( 2 π ) 3 δ ( p 2 − m 2 ) θ ( p 0 ) f ( p ) = ∫ d 3 p ( 2 π ) 3 2 ω p f ( ω p , p ) \int\frac{d^4p}{(2\pi)^3}\,\delta(p^2-m^2)\,\theta(p^0)\,f(p) = \int\frac{d^3p}{(2\pi)^3\,2\omega_{\boldsymbol{p}}}f(\omega_{\boldsymbol{p}},\boldsymbol{p}) ∫ ( 2 π ) 3 d 4 p δ ( p 2 − m 2 ) θ ( p 0 ) f ( p ) = ∫ ( 2 π ) 3 2 ω p d 3 p f ( ω p , p ) が成り立ちます。左辺の d 4 p d^4p d 4 p 、δ ( p 2 − m 2 ) \delta(p^2-m^2) δ ( p 2 − m 2 ) はローレンツ不変であり、θ ( p 0 ) \theta(p^0) θ ( p 0 ) も正質量殻上では固有直交ローレンツ変換 Λ ∈ S O + ( 1 , 3 ) \Lambda\in SO^+(1,3) Λ ∈ S O + ( 1 , 3 ) で不変です(Λ \Lambda Λ は p 0 p^0 p 0 の符号を保ちます)。よって D ( Λ x ) = D ( x ) D(\Lambda x) = D(x) D ( Λ x ) = D ( x ) が任意の Λ ∈ S O + ( 1 , 3 ) \Lambda\in SO^+(1,3) Λ ∈ S O + ( 1 , 3 ) について成り立ちます。
第 3 段:空間的ベクトルの反転。 z = x − y z = x-y z = x − y が空間的、すなわち z 2 = ( z 0 ) 2 − ∣ z ∣ 2 < 0 z^2 = (z^0)^2 - \lvert\boldsymbol{z}\rvert^2 < 0 z 2 = ( z 0 ) 2 − ∣ z ∣ 2 < 0 とします。すると ∣ z ∣ > ∣ z 0 ∣ \lvert\boldsymbol{z}\rvert > \lvert z^0\rvert ∣ z ∣ > ∣ z 0 ∣ なので、v = z 0 / ∣ z ∣ v = z^0/\lvert\boldsymbol{z}\rvert v = z 0 / ∣ z ∣ は ∣ v ∣ < 1 \lvert v\rvert < 1 ∣ v ∣ < 1 を満たします。z ^ \hat{\boldsymbol{z}} z ^ 方向への速度 v v v のブーストを施すと
z ′ 0 = γ ( z 0 − v ∣ z ∣ ) = 0 z'^0 = \gamma\left(z^0 - v\lvert\boldsymbol{z}\rvert\right) = 0 z ′0 = γ ( z 0 − v ∣ z ∣ ) = 0 となり、z ′ = ( 0 , z ′ ) z' = (0,\boldsymbol{z}') z ′ = ( 0 , z ′ ) という同時刻ベクトルになります。次に z ′ \boldsymbol{z}' z ′ に直交する軸のまわりの角度 π \pi π の回転を行うと z ′ → − z ′ \boldsymbol{z}'\to-\boldsymbol{z}' z ′ → − z ′ 、すなわち z ′ → − z ′ z'\to -z' z ′ → − z ′ です。ブーストも回転も S O + ( 1 , 3 ) SO^+(1,3) S O + ( 1 , 3 ) の元なので、これらの合成 Λ \Lambda Λ は Λ z = − z \Lambda z = -z Λ z = − z を満たします。
第 2 段より D ( z ) = D ( Λ z ) = D ( − z ) D(z) = D(\Lambda z) = D(-z) D ( z ) = D ( Λ z ) = D ( − z ) です。第 1 段に戻すと [ ϕ ^ ( x ) , ϕ ^ ( y ) ] = D ( z ) − D ( − z ) = 0 [\hat\phi(x),\hat\phi(y)] = D(z)-D(-z) = 0 [ ϕ ^ ( x ) , ϕ ^ ( y )] = D ( z ) − D ( − z ) = 0 を得ます。
∎
注意 7.6 (因果律が反粒子を要求する )
定理 7.5 の証明で本質的だったのは、D ( x − y ) D(x-y) D ( x − y ) と D ( y − x ) D(y-x) D ( y − x ) という二つの項が引き算で現れた ことです。実際 D D D そのものは空間的隔たりでゼロになりません。同時刻 x − y = ( 0 , r ) x-y = (0,\boldsymbol{r}) x − y = ( 0 , r ) では
D ( 0 , r ) = 1 4 π 2 r ∫ 0 ∞ d p p sin ( p r ) p 2 + m 2 = m 4 π 2 r K 1 ( m r ) > 0 D(0,\boldsymbol{r}) = \frac{1}{4\pi^2 r}\int_0^\infty dp\;\frac{p\sin(pr)}{\sqrt{p^2+m^2}} = \frac{m}{4\pi^2 r}K_1(mr) > 0 D ( 0 , r ) = 4 π 2 r 1 ∫ 0 ∞ d p p 2 + m 2 p sin ( p r ) = 4 π 2 r m K 1 ( m r ) > 0 (K 1 K_1 K 1 は第 2 種変形 Bessel 関数で、K 1 ( u ) ∼ π / 2 u e − u K_1(u) \sim \sqrt{\pi/2u}\,e^{-u} K 1 ( u ) ∼ π /2 u e − u )と、例 5.2 と同じ e − m r e^{-mr} e − m r 型の裾をもちます。ゼロになるのは差だけです。
複素スカラー場ではこの差の意味がはっきりします。D ( x − y ) D(x-y) D ( x − y ) は「y y y で生成された粒子 が x x x で消滅する振幅」、D ( y − x ) D(y-x) D ( y − x ) は「x x x で生成された反粒子 が y y y で消滅する振幅」です。空間的隔たりでは二つの事象に不変な時間順序が付かないため両者は同じ大きさになり、ちょうど打ち消し合います。反粒子が無ければ打ち消しは起こらず因果律が破れます。反粒子の存在は、相対論・量子力学・因果律から強制される帰結です (演習 9.4 )。
フェルミオンでは同じ打ち消しが反交換子で起こります。整数スピンは交換子、半整数スピンは反交換子で量子化しなければ因果律かエネルギーの正値性が壊れる、というのがスピン統計定理です。ディラックの海が必要としたパウリ排他律も、ここでは定理になります。
問題 一粒子の相対論的量子力学 場の量子論 負エネルギー解 ハミルトニアンが下に有界でない(定理 4.2 ) 負振動数の係数が生成演算子になり H ^ ≥ 0 \hat H\ge 0 H ^ ≥ 0 密度の符号 j 0 j^0 j 0 が負になる(命題 3.2 )j 0 j^0 j 0 は電荷密度。両符号は粒子と反粒子因果律 空間的隔たりでも振幅が非ゼロ(例 5.2 ) 空間的隔たりで交換子がゼロ(定理 7.5 ) 粒子数 変える演算子を書けない Fock 空間で自然に扱える(例 7.3 ) 同種粒子 対称化・反対称化を手で課す 同じ場の量子なので自動的に同一 スピンと統計 独立な仮定 スピン統計定理として導かれる
以降の章で、この記事で開いた扉の先を歩きます。
古典場の理論とラグランジアン :場の Lagrangian 形式と Noether の定理。定理 7.1 の出発点を厳密に整えます。
スカラー場の正準量子化 :等時刻交換関係 [ ϕ ^ ( t , x ) , π ^ ( t , y ) ] = i δ 3 ( x − y ) [\hat\phi(t,\boldsymbol{x}),\hat\pi(t,\boldsymbol{y})] = i\delta^3(\boldsymbol{x}-\boldsymbol{y}) [ ϕ ^ ( t , x ) , π ^ ( t , y )] = i δ 3 ( x − y ) から Fock 空間と伝播関数へ。
経路積分量子化 :演算子を使わない定式化。ゲージ理論で威力を発揮します。
くりこみ理論入門 :相互作用する場に現れる発散が、理論の有効なスケールを示す仕組みであること。
ゲージ理論と自発的対称性の破れ :標準模型の骨格。
摂動論 の計算技法はそのまま基礎になります。例 7.3 の相互作用項を摂動として扱うのが Feynman 図の出発点です。
演習 9.1 標準
自然単位系 ℏ = c = 1 \hbar=c=1 ℏ = c = 1 、計量 η μ ν = d i a g ( + 1 , − 1 , − 1 , − 1 ) \eta_{\mu\nu} = \mathrm{diag}(+1,-1,-1,-1) η μν = diag ( + 1 , − 1 , − 1 , − 1 ) とし、ϕ \phi ϕ をクライン・ゴルドン方程式 ( □ + m 2 ) ϕ = 0 (\Box+m^2)\phi=0 ( □ + m 2 ) ϕ = 0 の複素解とします。
カレント j μ = i 2 m ( ϕ ∗ ∂ μ ϕ − ( ∂ μ ϕ ∗ ) ϕ ) j^\mu = \frac{i}{2m}\left(\phi^{*}\partial^\mu\phi - (\partial^\mu\phi^{*})\phi\right) j μ = 2 m i ( ϕ ∗ ∂ μ ϕ − ( ∂ μ ϕ ∗ ) ϕ ) の空間成分 j \boldsymbol{j} j を ϕ \phi ϕ と ∇ ϕ \nabla\phi ∇ ϕ で書き下し、非相対論的な確率流 1 2 m i ( ψ ∗ ∇ ψ − ψ ∇ ψ ∗ ) \frac{1}{2mi}(\psi^{*}\nabla\psi - \psi\nabla\psi^{*}) 2 mi 1 ( ψ ∗ ∇ ψ − ψ ∇ ψ ∗ ) と同じ形であることを確かめてください。
重ね合わせ ϕ = N + e − i ω t + N − e + i ω t \phi = N_+e^{-i\omega t} + N_-e^{+i\omega t} ϕ = N + e − iω t + N − e + iω t (p = 0 \boldsymbol{p}=\boldsymbol{0} p = 0 、ω = m \omega=m ω = m )に対して j 0 j^0 j 0 を計算し、∣ N + ∣ = ∣ N − ∣ \lvert N_+\rvert = \lvert N_-\rvert ∣ N + ∣ = ∣ N − ∣ のとき j 0 = 0 j^0 = 0 j 0 = 0 となることを示してください。これは何を意味しますか。
解答 1. ∂ i = η i i ∂ i = − ∂ / ∂ x i \partial^i = \eta^{ii}\partial_i = -\partial/\partial x^i ∂ i = η ii ∂ i = − ∂ / ∂ x i ですから(i i i について和を取らない)、空間成分は
j = i 2 m ( ϕ ∗ ( − ∇ ϕ ) − ( − ∇ ϕ ∗ ) ϕ ) = − i 2 m ( ϕ ∗ ∇ ϕ − ( ∇ ϕ ∗ ) ϕ ) \boldsymbol{j} = \frac{i}{2m}\left(\phi^{*}(-\nabla\phi) - (-\nabla\phi^{*})\phi\right) = -\frac{i}{2m}\left(\phi^{*}\nabla\phi - (\nabla\phi^{*})\phi\right) j = 2 m i ( ϕ ∗ ( − ∇ ϕ ) − ( − ∇ ϕ ∗ ) ϕ ) = − 2 m i ( ϕ ∗ ∇ ϕ − ( ∇ ϕ ∗ ) ϕ ) です。− i / ( 2 m ) = 1 / ( 2 m i ) -i/(2m) = 1/(2mi) − i / ( 2 m ) = 1/ ( 2 mi ) なので
j = 1 2 m i ( ϕ ∗ ∇ ϕ − ϕ ∇ ϕ ∗ ) \boldsymbol{j} = \frac{1}{2mi}\left(\phi^{*}\nabla\phi - \phi\nabla\phi^{*}\right) j = 2 mi 1 ( ϕ ∗ ∇ ϕ − ϕ ∇ ϕ ∗ ) となり、非相対論的な確率流と同じ形です。相対論化で形が変わったのは時間成分だけだったことがわかります。
2. ω = m \omega = m ω = m とします。ϕ ˙ = − i m ( N + e − i m t − N − e + i m t ) \dot\phi = -im\left(N_+e^{-imt} - N_-e^{+imt}\right) ϕ ˙ = − im ( N + e − im t − N − e + im t ) と ϕ ∗ = N + ‾ e + i m t + N − ‾ e − i m t \phi^{*} = \overline{N_+}e^{+imt} + \overline{N_-}e^{-imt} ϕ ∗ = N + e + im t + N − e − im t を掛け、w : = N − ‾ N + e − 2 i m t w := \overline{N_-}N_+e^{-2imt} w := N − N + e − 2 im t (このとき N + ‾ N − e 2 i m t = w ˉ \overline{N_+}N_-e^{2imt} = \bar w N + N − e 2 im t = w ˉ )とおくと
ϕ ∗ ϕ ˙ = − i m [ ( ∣ N + ∣ 2 − ∣ N − ∣ 2 ) + ( w − w ˉ ) ] = − i m ( ∣ N + ∣ 2 − ∣ N − ∣ 2 ) + 2 m Im w \phi^{*}\dot\phi = -im\left[\left(\lvert N_+\rvert^2-\lvert N_-\rvert^2\right) + (w-\bar w)\right]
= -im\left(\lvert N_+\rvert^2-\lvert N_-\rvert^2\right) + 2m\operatorname{Im}w ϕ ∗ ϕ ˙ = − im [ ( ∣ N + ∣ 2 − ∣ N − ∣ 2 ) + ( w − w ˉ ) ] = − im ( ∣ N + ∣ 2 − ∣ N − ∣ 2 ) + 2 m Im w です(w − w ˉ = 2 i Im w w-\bar w = 2i\operatorname{Im}w w − w ˉ = 2 i Im w を使いました)。第 2 項は実数なので、ϕ ˙ ∗ ϕ = ϕ ∗ ϕ ˙ ‾ \dot\phi^{*}\phi = \overline{\phi^{*}\dot\phi} ϕ ˙ ∗ ϕ = ϕ ∗ ϕ ˙ との差を取ると相殺し、
j 0 = i 2 m ( ϕ ∗ ϕ ˙ − ϕ ˙ ∗ ϕ ) = i 2 m ⋅ ( − 2 i m ) ( ∣ N + ∣ 2 − ∣ N − ∣ 2 ) = ∣ N + ∣ 2 − ∣ N − ∣ 2 j^0 = \frac{i}{2m}\left(\phi^{*}\dot\phi - \dot\phi^{*}\phi\right) = \frac{i}{2m}\cdot\left(-2im\right)\left(\lvert N_+\rvert^2-\lvert N_-\rvert^2\right) = \lvert N_+\rvert^2 - \lvert N_-\rvert^2 j 0 = 2 m i ( ϕ ∗ ϕ ˙ − ϕ ˙ ∗ ϕ ) = 2 m i ⋅ ( − 2 im ) ( ∣ N + ∣ 2 − ∣ N − ∣ 2 ) = ∣ N + ∣ 2 − ∣ N − ∣ 2 が残ります。よって ∣ N + ∣ = ∣ N − ∣ \lvert N_+\rvert=\lvert N_-\rvert ∣ N + ∣ = ∣ N − ∣ のとき j 0 = 0 j^0 = 0 j 0 = 0 です。
意味:正エネルギー成分と負エネルギー成分が等量ある状態では、j 0 j^0 j 0 を「粒子がそこにいる確率密度」と読むと恒等的にゼロという不合理が生じます。注意 3.4 のとおり、j 0 j^0 j 0 は電荷密度であり、正電荷と負電荷が等量あって中性になっている、と読むのが正しい解釈です。
演習 9.2 標準
α 1 , α 2 , α 3 , β \alpha^1,\alpha^2,\alpha^3,\beta α 1 , α 2 , α 3 , β をエルミートな N × N N\times N N × N 行列で { α i , α j } = 2 δ i j 1 \{\alpha^i,\alpha^j\}=2\delta^{ij}\mathbb{1} { α i , α j } = 2 δ ij 1 、{ α i , β } = 0 \{\alpha^i,\beta\}=0 { α i , β } = 0 、β 2 = 1 \beta^2=\mathbb{1} β 2 = 1 を満たすものとします。
N N N が偶数であることを示してください。
N = 2 N=2 N = 2 が不可能であることを示し、最小次元が 4 4 4 であることを結論してください。
解答 1. 定理 4.2 で示したとおり tr β = 0 \operatorname{tr}\beta = 0 tr β = 0 です。一方 β 2 = 1 \beta^2=\mathbb{1} β 2 = 1 かつ β \beta β はエルミートなので対角化でき、固有値は ± 1 \pm1 ± 1 のみです。重複度を n ± n_\pm n ± とすると n + + n − = N n_++n_-=N n + + n − = N かつ tr β = n + − n − = 0 \operatorname{tr}\beta = n_+-n_-=0 tr β = n + − n − = 0 、よって n + = n − = N / 2 n_+=n_-=N/2 n + = n − = N /2 です。n ± n_\pm n ± は整数なので N N N は偶数です。
2. N = 2 N=2 N = 2 と仮定します。2 × 2 2\times2 2 × 2 のエルミート行列でトレースが 0 0 0 のものは、実係数で σ 1 , σ 2 , σ 3 \sigma^1,\sigma^2,\sigma^3 σ 1 , σ 2 , σ 3 の線形結合として一意に書けます(この空間は実 3 次元)。定理 4.2 より α i \alpha^i α i も β \beta β もトレースが 0 0 0 なので、実ベクトル a 1 , a 2 , a 3 , b ∈ R 3 \boldsymbol{a}_1,\boldsymbol{a}_2,\boldsymbol{a}_3,\boldsymbol{b}\in\mathbb{R}^3 a 1 , a 2 , a 3 , b ∈ R 3 を使って
α i = a i ⋅ σ , β = b ⋅ σ \alpha^i = \boldsymbol{a}_i\cdot\boldsymbol{\sigma},\qquad \beta = \boldsymbol{b}\cdot\boldsymbol{\sigma} α i = a i ⋅ σ , β = b ⋅ σ と書けます。パウリ行列の恒等式 ( u ⋅ σ ) ( v ⋅ σ ) = ( u ⋅ v ) 1 + i ( u × v ) ⋅ σ (\boldsymbol{u}\cdot\boldsymbol{\sigma})(\boldsymbol{v}\cdot\boldsymbol{\sigma}) = (\boldsymbol{u}\cdot\boldsymbol{v})\mathbb{1} + i(\boldsymbol{u}\times\boldsymbol{v})\cdot\boldsymbol{\sigma} ( u ⋅ σ ) ( v ⋅ σ ) = ( u ⋅ v ) 1 + i ( u × v ) ⋅ σ を使うと、u × v \boldsymbol{u}\times\boldsymbol{v} u × v が反対称なので
{ u ⋅ σ , v ⋅ σ } = 2 ( u ⋅ v ) 1 \{\boldsymbol{u}\cdot\boldsymbol{\sigma},\,\boldsymbol{v}\cdot\boldsymbol{\sigma}\} = 2(\boldsymbol{u}\cdot\boldsymbol{v})\,\mathbb{1} { u ⋅ σ , v ⋅ σ } = 2 ( u ⋅ v ) 1 が成り立ちます。したがって仮定の三条件は、a 1 , a 2 , a 3 , b \boldsymbol{a}_1,\boldsymbol{a}_2,\boldsymbol{a}_3,\boldsymbol{b} a 1 , a 2 , a 3 , b が
a i ⋅ a j = δ i j , a i ⋅ b = 0 , b ⋅ b = 1 \boldsymbol{a}_i\cdot\boldsymbol{a}_j = \delta_{ij},\qquad \boldsymbol{a}_i\cdot\boldsymbol{b} = 0,\qquad \boldsymbol{b}\cdot\boldsymbol{b} = 1 a i ⋅ a j = δ ij , a i ⋅ b = 0 , b ⋅ b = 1 を満たすこと、すなわち R 3 \mathbb{R}^3 R 3 の中に互いに直交する単位ベクトルが 4 本あることと同値です。R 3 \mathbb{R}^3 R 3 の次元は 3 3 3 なので、直交する非ゼロベクトルは高々 3 本しか取れません。矛盾します。よって N ≠ 2 N\ne 2 N = 2 。
1 より N N N は偶数、いま N = 2 N=2 N = 2 が排除され、N = 0 N=0 N = 0 は行列として無意味ですから N ≥ 4 N\ge 4 N ≥ 4 です。そして N = 4 N=4 N = 4 は実際に実現します(ディラック表現 α i = ( 0 σ i σ i 0 ) \alpha^i = \begin{pmatrix}0&\sigma^i\\\sigma^i&0\end{pmatrix} α i = ( 0 σ i σ i 0 ) 、β = ( 1 0 0 − 1 ) \beta=\begin{pmatrix}\mathbb{1}&0\\0&-\mathbb{1}\end{pmatrix} β = ( 1 0 0 − 1 ) が三条件を満たすことは、{ σ i , σ j } = 2 δ i j 1 \{\sigma^i,\sigma^j\}=2\delta^{ij}\mathbb{1} { σ i , σ j } = 2 δ ij 1 を使ってブロックごとに確かめられます)。よって最小次元は 4 4 4 です。
演習 9.3 易
電子の換算コンプトン波長 λ ˉ C = ℏ / ( m e c ) \bar\lambda_C = \hbar/(m_ec) λ ˉ C = ℏ/ ( m e c ) と Bohr 半径 a 0 = ℏ / ( m e c α ) a_0 = \hbar/(m_ec\alpha) a 0 = ℏ/ ( m e c α ) の比を求め、数値を出してください(α − 1 = 137.04 \alpha^{-1}=137.04 α − 1 = 137.04 )。
水素原子中の電子の典型的な速度と運動エネルギーを α \alpha α と m e c 2 = 511 k e V m_ec^2 = 511\ \mathrm{keV} m e c 2 = 511 keV で表し、非相対論的量子力学が良い近似である理由を述べてください。
電子を λ ˉ C \bar\lambda_C λ ˉ C 程度の領域に閉じ込めようとすると何が起こりますか。
解答 1. 定義から直ちに a 0 / λ ˉ C = 1 / α = 137.04 a_0/\bar\lambda_C = 1/\alpha = 137.04 a 0 / λ ˉ C = 1/ α = 137.04 です。例 6.2 の λ ˉ C = 3.86 × 10 − 13 m \bar\lambda_C = 3.86\times10^{-13}\ \mathrm{m} λ ˉ C = 3.86 × 1 0 − 13 m を使うと
a 0 = 137.04 × 3.86 × 10 − 13 m = 5.29 × 10 − 11 m a_0 = 137.04\times 3.86\times10^{-13}\ \mathrm{m} = 5.29\times10^{-11}\ \mathrm{m} a 0 = 137.04 × 3.86 × 1 0 − 13 m = 5.29 × 1 0 − 11 m で、既知の Bohr 半径に一致します。
2. 電子が a 0 a_0 a 0 程度の広がりに閉じ込められているので、不確定性関係から典型的な運動量は p ∼ ℏ / a 0 = α m e c p\sim\hbar/a_0 = \alpha\, m_ec p ∼ ℏ/ a 0 = α m e c です。よって
v c ∼ α ≃ 1 137 , E kin ∼ p 2 2 m e = 1 2 α 2 m e c 2 = 5.325 × 10 − 5 2 × 511 k e V = 13.6 e V . \frac{v}{c}\sim\alpha\simeq\frac{1}{137},
\qquad
E_{\text{kin}}\sim\frac{p^2}{2m_e} = \frac{1}{2}\alpha^2 m_ec^2 = \frac{5.325\times10^{-5}}{2}\times 511\ \mathrm{keV} = 13.6\ \mathrm{eV}. c v ∼ α ≃ 137 1 , E kin ∼ 2 m e p 2 = 2 1 α 2 m e c 2 = 2 5.325 × 1 0 − 5 × 511 keV = 13.6 eV . これは水素原子の Rydberg エネルギーそのものです。運動エネルギーは静止エネルギーの 2.7 × 10 − 5 2.7\times10^{-5} 2.7 × 1 0 − 5 倍にすぎず、相対論補正は ( v / c ) 2 ∼ α 2 (v/c)^2\sim\alpha^2 ( v / c ) 2 ∼ α 2 の相対精度、すなわち 10 − 5 10^{-5} 1 0 − 5 程度の効果です(これが微細構造の大きさです)。したがって非相対論的量子力学が良い近似になります。原子物理が量子力学だけでほぼ足りるのは、微細構造定数が小さいおかげです。
3. Δ x ∼ λ ˉ C \Delta x\sim\bar\lambda_C Δ x ∼ λ ˉ C とすると Δ p ≳ ℏ / λ ˉ C = m e c \Delta p\gtrsim\hbar/\bar\lambda_C = m_ec Δ p ≳ ℏ/ λ ˉ C = m e c 、よって Δ E ≳ m e c 2 = 511 k e V \Delta E\gtrsim m_ec^2 = 511\ \mathrm{keV} Δ E ≳ m e c 2 = 511 keV です。対生成のしきい値 2 m e c 2 = 1.022 M e V 2m_ec^2 = 1.022\ \mathrm{MeV} 2 m e c 2 = 1.022 MeV と同じ桁なので、閉じ込めに使ったエネルギーで電子・陽電子対が生成されます。「もとの 1 個の電子」と「新しく生まれた電子」は区別できないため、一粒子の波動関数という記述そのものが破綻します。定義 7.2 のように粒子数の異なる状態を含む空間が必要になります。
演習 9.4 難
複素スカラー場の演算子を
ϕ ^ ( x ) = ∫ d 3 p ( 2 π ) 3 1 2 ω p ( a ^ p e − i p ⋅ x + b ^ p † e + i p ⋅ x ) \hat\phi(x) = \int\frac{d^3p}{(2\pi)^3}\frac{1}{\sqrt{2\omega_{\boldsymbol{p}}}}\left(\hat a_{\boldsymbol{p}}e^{-ip\cdot x} + \hat b^\dagger_{\boldsymbol{p}}e^{+ip\cdot x}\right) ϕ ^ ( x ) = ∫ ( 2 π ) 3 d 3 p 2 ω p 1 ( a ^ p e − i p ⋅ x + b ^ p † e + i p ⋅ x ) とします([ a ^ p , a ^ q † ] = [ b ^ p , b ^ q † ] = ( 2 π ) 3 δ 3 ( p − q ) [\hat a_{\boldsymbol{p}},\hat a^\dagger_{\boldsymbol{q}}] = [\hat b_{\boldsymbol{p}},\hat b^\dagger_{\boldsymbol{q}}] = (2\pi)^3\delta^3(\boldsymbol{p}-\boldsymbol{q}) [ a ^ p , a ^ q † ] = [ b ^ p , b ^ q † ] = ( 2 π ) 3 δ 3 ( p − q ) 、a ^ \hat a a ^ 型と b ^ \hat b b ^ 型の間の交換子はすべてゼロ)。
[ ϕ ^ ( x ) , ϕ ^ † ( y ) ] = D ( x − y ) − D ( y − x ) [\hat\phi(x),\hat\phi^\dagger(y)] = D(x-y) - D(y-x) [ ϕ ^ ( x ) , ϕ ^ † ( y )] = D ( x − y ) − D ( y − x ) を示し、空間的隔たりでゼロになることを結論してください(D D D は 定理 7.5 のもの)。
D ( x − y ) D(x-y) D ( x − y ) と D ( y − x ) D(y-x) D ( y − x ) がそれぞれ物理的に何の振幅かを述べ、打ち消しの意味を説明してください。
反粒子が存在しない(b ^ † \hat b^\dagger b ^ † の項が無い)と仮定すると何が壊れますか。
解答 1. ϕ ^ † ( y ) = ∫ d 3 q ( 2 π ) 3 1 2 ω q ( a ^ q † e + i q ⋅ y + b ^ q e − i q ⋅ y ) \hat\phi^\dagger(y) = \int\frac{d^3q}{(2\pi)^3}\frac{1}{\sqrt{2\omega_{\boldsymbol{q}}}}\left(\hat a^\dagger_{\boldsymbol{q}}e^{+iq\cdot y} + \hat b_{\boldsymbol{q}}e^{-iq\cdot y}\right) ϕ ^ † ( y ) = ∫ ( 2 π ) 3 d 3 q 2 ω q 1 ( a ^ q † e + i q ⋅ y + b ^ q e − i q ⋅ y ) です。交換子を展開すると、a ^ \hat a a ^ 型と b ^ \hat b b ^ 型が交わる項はゼロなので
[ ϕ ^ ( x ) , ϕ ^ † ( y ) ] = ∫ d 3 p d 3 q ( 2 π ) 6 4 ω p ω q ( [ a ^ p , a ^ q † ] e − i p ⋅ x + i q ⋅ y + [ b ^ p † , b ^ q ] e + i p ⋅ x − i q ⋅ y ) [\hat\phi(x),\hat\phi^\dagger(y)] = \int\frac{d^3p\,d^3q}{(2\pi)^6\sqrt{4\omega_{\boldsymbol{p}}\omega_{\boldsymbol{q}}}}\left([\hat a_{\boldsymbol{p}},\hat a^\dagger_{\boldsymbol{q}}]e^{-ip\cdot x+iq\cdot y} + [\hat b^\dagger_{\boldsymbol{p}},\hat b_{\boldsymbol{q}}]e^{+ip\cdot x - iq\cdot y}\right) [ ϕ ^ ( x ) , ϕ ^ † ( y )] = ∫ ( 2 π ) 6 4 ω p ω q d 3 p d 3 q ( [ a ^ p , a ^ q † ] e − i p ⋅ x + i q ⋅ y + [ b ^ p † , b ^ q ] e + i p ⋅ x − i q ⋅ y ) = ∫ d 3 p ( 2 π ) 3 2 ω p ( e − i p ⋅ ( x − y ) − e + i p ⋅ ( x − y ) ) = D ( x − y ) − D ( y − x ) . = \int\frac{d^3p}{(2\pi)^3 2\omega_{\boldsymbol{p}}}\left(e^{-ip\cdot(x-y)} - e^{+ip\cdot(x-y)}\right) = D(x-y)-D(y-x). = ∫ ( 2 π ) 3 2 ω p d 3 p ( e − i p ⋅ ( x − y ) − e + i p ⋅ ( x − y ) ) = D ( x − y ) − D ( y − x ) . 第 2 項の符号は [ b ^ p † , b ^ q ] = − ( 2 π ) 3 δ 3 ( p − q ) [\hat b^\dagger_{\boldsymbol{p}},\hat b_{\boldsymbol{q}}] = -(2\pi)^3\delta^3(\boldsymbol{p}-\boldsymbol{q}) [ b ^ p † , b ^ q ] = − ( 2 π ) 3 δ 3 ( p − q ) から来ています。あとは 定理 7.5 の第 2 段・第 3 段がそのまま使えます。( x − y ) 2 < 0 (x-y)^2<0 ( x − y ) 2 < 0 ならば Λ ( x − y ) = − ( x − y ) \Lambda(x-y) = -(x-y) Λ ( x − y ) = − ( x − y ) となる Λ ∈ S O + ( 1 , 3 ) \Lambda\in SO^+(1,3) Λ ∈ S O + ( 1 , 3 ) が存在し、D D D の不変性から D ( x − y ) = D ( y − x ) D(x-y)=D(y-x) D ( x − y ) = D ( y − x ) 、よって交換子はゼロです。
2. 真空期待値を取ると ⟨ 0 ∣ ϕ ^ ( x ) ϕ ^ † ( y ) ∣ 0 ⟩ = D ( x − y ) \langle 0\rvert\hat\phi(x)\hat\phi^\dagger(y)\lvert 0\rangle = D(x-y) ⟨ 0 ∣ ϕ ^ ( x ) ϕ ^ † ( y ) ∣ 0 ⟩ = D ( x − y ) です。ϕ ^ † ( y ) \hat\phi^\dagger(y) ϕ ^ † ( y ) は y y y で粒子(a ^ \hat a a ^ 型)を 1 個生成し、ϕ ^ ( x ) \hat\phi(x) ϕ ^ ( x ) は x x x でそれを消滅させますから、D ( x − y ) D(x-y) D ( x − y ) は「y y y で生まれた粒子が x x x で消える振幅」です。一方 ⟨ 0 ∣ ϕ ^ † ( y ) ϕ ^ ( x ) ∣ 0 ⟩ = D ( y − x ) \langle 0\rvert\hat\phi^\dagger(y)\hat\phi(x)\lvert 0\rangle = D(y-x) ⟨ 0 ∣ ϕ ^ † ( y ) ϕ ^ ( x ) ∣ 0 ⟩ = D ( y − x ) で、ϕ ^ ( x ) \hat\phi(x) ϕ ^ ( x ) は x x x で反粒子(b ^ \hat b b ^ 型)を生成し、ϕ ^ † ( y ) \hat\phi^\dagger(y) ϕ ^ † ( y ) が y y y でそれを消滅させます。すなわち「x x x で生まれた反粒子が y y y で消える振幅」です。
空間的に隔たった 2 点にはローレンツ不変な時間順序が付きません。「y y y が先」と見る観測者には粒子が y → x y\to x y → x と伝わったように見え、「x x x が先」と見る観測者には反粒子が x → y x\to y x → y と伝わったように見えます。二つの振幅は等しく、交換子では引き算されて消えます。観測可能な因果的影響が残らない仕組みです。
3. b ^ † \hat b^\dagger b ^ † の項が無ければ ϕ ^ ( x ) = ∫ d 3 p ( 2 π ) 3 2 ω p a ^ p e − i p ⋅ x \hat\phi(x) = \int\frac{d^3p}{(2\pi)^3\sqrt{2\omega_{\boldsymbol p}}}\hat a_{\boldsymbol{p}}e^{-ip\cdot x} ϕ ^ ( x ) = ∫ ( 2 π ) 3 2 ω p d 3 p a ^ p e − i p ⋅ x となり、[ ϕ ^ ( x ) , ϕ ^ † ( y ) ] = D ( x − y ) [\hat\phi(x),\hat\phi^\dagger(y)] = D(x-y) [ ϕ ^ ( x ) , ϕ ^ † ( y )] = D ( x − y ) です。注意 7.6 のとおり D D D は空間的隔たりで m 4 π 2 r K 1 ( m r ) > 0 \frac{m}{4\pi^2r}K_1(mr)>0 4 π 2 r m K 1 ( m r ) > 0 となりゼロになりません。よってミクロ因果律が破れ、空間的に隔たった二つの観測が互いに影響し合うことになります。
つまり相対論的な局所場を作るには、消滅演算子と生成演算子の両方を同じ ϕ ^ ( x ) \hat\phi(x) ϕ ^ ( x ) に含めねばなりません。場が電荷をもつ(ϕ ^ ≠ ϕ ^ † \hat\phi\ne\hat\phi^\dagger ϕ ^ = ϕ ^ † )とき、後者は必然的に反対符号の電荷をもつ粒子、すなわち反粒子の生成演算子になります。中性の実スカラー場(定義 7.4 )では b ^ = a ^ \hat b = \hat a b ^ = a ^ で、粒子が自分自身の反粒子です。
M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory , Addison-Wesley, 1995 — 第 1 章および §2.1(一粒子相対論的量子力学の破綻と光円錐外の伝播振幅)。
S. Weinberg, The Quantum Theory of Fields, Volume I: Foundations , Cambridge University Press, 1995 — 第 1 章(歴史的序説)、第 5 章(場の必然性とスピン統計)。
坂本眞人『場の量子論 — 不変性と自由場を中心にして』裳華房、2014 — 第 1 章〜第 3 章(一粒子理論の困難、正準量子化)。
P. A. M. Dirac, “The Quantum Theory of the Electron”, Proceedings of the Royal Society A 117 (1928), 610–624. DOI: 10.1098/rspa.1928.0023
C. D. Anderson, “The Positive Electron”, Physical Review 43 (1933), 491–494. DOI: 10.1103/PhysRev.43.491
G. C. Hegerfeldt, “Remark on causality and particle localization”, Physical Review D 10 (1974), 3320–3321. DOI: 10.1103/PhysRevD.10.3320
例 5.2 で使った表式を導きます。出発点は
U ( t , r ) = 1 2 π 2 r ∫ 0 ∞ d p p sin ( p r ) e − i ω p t , ω p = p 2 + m 2 U(t,r) = \frac{1}{2\pi^2 r}\int_0^\infty dp\;p\,\sin(pr)\,e^{-i\omega_p t},
\qquad \omega_p = \sqrt{p^2+m^2} U ( t , r ) = 2 π 2 r 1 ∫ 0 ∞ d p p sin ( p r ) e − i ω p t , ω p = p 2 + m 2
です。被積分関数のうち p sin ( p r ) p\sin(pr) p sin ( p r ) と ω p \omega_p ω p はどちらも p p p の偶関数なので、積分区間を R \mathbb{R} R 全体に広げて 1 2 \frac12 2 1 を掛けられます。さらに sin ( p r ) = ( e i p r − e − i p r ) / ( 2 i ) \sin(pr) = (e^{ipr}-e^{-ipr})/(2i) sin ( p r ) = ( e i p r − e − i p r ) / ( 2 i ) とし、e − i p r e^{-ipr} e − i p r の項で p → − p p\to-p p → − p と置換すると e i p r e^{ipr} e i p r の項と同じものになるので、
U ( t , r ) = 1 4 π 2 i r ∫ − ∞ ∞ d p p e i p r e − i p 2 + m 2 t U(t,r) = \frac{1}{4\pi^2 i\,r}\int_{-\infty}^{\infty}dp\;p\,e^{ipr}\,e^{-i\sqrt{p^2+m^2}\,t} U ( t , r ) = 4 π 2 i r 1 ∫ − ∞ ∞ d p p e i p r e − i p 2 + m 2 t
を得ます。ここで p 2 + m 2 \sqrt{p^2+m^2} p 2 + m 2 を複素 p p p 平面へ解析接続します。分岐点は p = ± i m p=\pm im p = ± im で、切断を [ i m , i ∞ ) [im, i\infty) [ im , i ∞ ) と ( − i ∞ , − i m ] (-i\infty,-im] ( − i ∞ , − im ] に取り、実軸上で p 2 + m 2 > 0 \sqrt{p^2+m^2}>0 p 2 + m 2 > 0 となる枝を選びます。
r > 0 r>0 r > 0 なので e i p r e^{ipr} e i p r は上半平面で減衰します。大きな半径 R R R の上半円弧では ∣ p ∣ → ∞ \lvert p\rvert\to\infty ∣ p ∣ → ∞ で p 2 + m 2 ≃ ± p \sqrt{p^2+m^2}\simeq \pm p p 2 + m 2 ≃ ± p となり、被積分関数の指数部は ∣ e i p ( r − t ) ∣ = e − ( r − t ) Im p \lvert e^{ip(r-t)}\rvert = e^{-(r-t)\operatorname{Im}p} ∣ e i p ( r − t ) ∣ = e − ( r − t ) Im p のようにふるまいます。空間的隔たり r > t > 0 r>t>0 r > t > 0 のときに限り これは減衰し、円弧の寄与が消えます。ここが決定的で、時間的隔たり r < t r<t r < t では以下の変形はできません。
積分路を上へ持ち上げると切断 [ i m , i ∞ ) [im,i\infty) [ im , i ∞ ) に引っかかります。結果として、切断の左側を i ∞ i\infty i ∞ から i m im im まで下り、i m im im を回って右側を i ∞ i\infty i ∞ まで上る積分路が残ります。p = i ρ p = i\rho p = i ρ (ρ > m \rho>m ρ > m )とおくと d p = i d ρ dp = i\,d\rho d p = i d ρ 、e i p r = e − ρ r e^{ipr} = e^{-\rho r} e i p r = e − ρ r です。p = ϵ + i ρ p=\epsilon+i\rho p = ϵ + i ρ として p 2 + m 2 ≃ ( m 2 − ρ 2 ) + 2 i ϵ ρ p^2+m^2 \simeq (m^2-\rho^2) + 2i\epsilon\rho p 2 + m 2 ≃ ( m 2 − ρ 2 ) + 2 i ϵ ρ を見ると、ρ > m \rho>m ρ > m で実部は負、虚部の符号は ϵ \epsilon ϵ の符号と同じです。よって s : = ρ 2 − m 2 s:=\sqrt{\rho^2-m^2} s := ρ 2 − m 2 とおくと、切断の右側(ϵ > 0 \epsilon>0 ϵ > 0 )では p 2 + m 2 = + i s \sqrt{p^2+m^2} = +is p 2 + m 2 = + i s 、左側(ϵ < 0 \epsilon<0 ϵ < 0 )では − i s -is − i s となります。したがって
e − i p 2 + m 2 t = { e + t s ( 右側 ) e − t s ( 左側 ) e^{-i\sqrt{p^2+m^2}\,t} = \begin{cases} e^{+ts} & (\text{右側}) \\ e^{-ts} & (\text{左側}) \end{cases} e − i p 2 + m 2 t = { e + t s e − t s ( 右側 ) ( 左側 )
です。左側(ρ : ∞ → m \rho:\infty\to m ρ : ∞ → m )と右側(ρ : m → ∞ \rho: m\to\infty ρ : m → ∞ )の寄与を足すと
∫ − ∞ ∞ d p p e i p r e − i ω p t = ∫ ∞ m ( i ρ ) e − ρ r e − t s i d ρ + ∫ m ∞ ( i ρ ) e − ρ r e + t s i d ρ \int_{-\infty}^{\infty}dp\;p\,e^{ipr}e^{-i\omega_pt}
= \int_{\infty}^{m}(i\rho)e^{-\rho r}e^{-ts}\,i\,d\rho + \int_{m}^{\infty}(i\rho)e^{-\rho r}e^{+ts}\,i\,d\rho ∫ − ∞ ∞ d p p e i p r e − i ω p t = ∫ ∞ m ( i ρ ) e − ρ r e − t s i d ρ + ∫ m ∞ ( i ρ ) e − ρ r e + t s i d ρ
= ∫ m ∞ d ρ ρ e − ρ r ( e − t s − e + t s ) = − 2 ∫ m ∞ d ρ ρ e − ρ r sinh ( t s ) . = \int_m^\infty d\rho\;\rho\,e^{-\rho r}\left(e^{-ts} - e^{+ts}\right)
= -2\int_m^\infty d\rho\;\rho\,e^{-\rho r}\sinh(ts). = ∫ m ∞ d ρ ρ e − ρ r ( e − t s − e + t s ) = − 2 ∫ m ∞ d ρ ρ e − ρ r sinh ( t s ) .
これを U U U の表式に戻すと、− 2 / ( 4 i ) = i / 2 -2/(4i) = i/2 − 2/ ( 4 i ) = i /2 より
U ( t , r ) = i 2 π 2 r ∫ m ∞ d ρ ρ e − ρ r sinh ( t ρ 2 − m 2 ) U(t,r) = \frac{i}{2\pi^2 r}\int_m^\infty d\rho\;\rho\,e^{-\rho r}\,\sinh\!\left(t\sqrt{\rho^2-m^2}\right) U ( t , r ) = 2 π 2 r i ∫ m ∞ d ρ ρ e − ρ r sinh ( t ρ 2 − m 2 )
が得られます。ρ → ∞ \rho\to\infty ρ → ∞ で被積分関数は ρ e − ( r − t ) ρ / 2 \rho\,e^{-(r-t)\rho}/2 ρ e − ( r − t ) ρ /2 のようにふるまうので、r > t r>t r > t で収束します。t > 0 t>0 t > 0 のとき被積分関数は区間内で厳密に正なので、積分は正、すなわち U ≠ 0 U\ne 0 U = 0 です。t → 0 t\to 0 t → 0 では sinh → 0 \sinh\to0 sinh → 0 となり U → 0 U\to0 U → 0 で、r > 0 r>0 r > 0 における δ 3 ( r ) \delta^3(\boldsymbol{r}) δ 3 ( r ) の値と整合します。