束縛電子の自己エネルギー


ファインマンの第6講の最後には「自己エネルギー」の記述があるが, そこの式 (1.6.12) の導出などは書かれていない.より詳しい議論が J.J.Sakurai の§ 2-8 にあったので, それを書いておこう.

Sakurai book cover

自己エネルギー問題(Self-energy problem)

考えるのは束縛電子または自由電子としよう.電子が存在することがまさに, その電子に付随する電磁場が存在することを意味している.古典的には, 電子の存在に依る電磁場はその電子自身に反応する可能性があり, H.ポアンケア や M.エイブラハム そして H.A.ローレンツによって詳しく議論された [1] Jackson (1962), pp. 578-597..放射の量子論では類似の効果があり, それは次の2段階のプロセスとして視覚化できる:

  1. 電子が(仮想)光子を放出する.
  2. その光子は同じ電子に吸収される.

電子自身の存在による電子と電磁場の(古典的または量子化された)相互作用エネルギーは, 電子の「自己エネルギー」と呼ばれる.現実的には電磁相互作用をオフにする方法がないため, 電子の自己エネルギーはその電子の観測可能な静止エネルギーとは切り離せない一部分になっている.この節では(観測不可能な)自己エネルギーを生み出す力学的メカニズムが, 原子状態の(観測可能な)エネルギー準位に重大な影響を及ぼすことを示そうと思う.
自己エネルギーの問題はいくつかのレベルで論じることが出来る.まず始めに, この問題は古典静電気学でも生じる.電荷分布 ρ を持つ電子モデルを考えることにする.電子が存在することで生じる静電ポテンシャルを ϕ とする.すると相互作用エネルギーは次で与えられる:

(2.222)Eint=12ρϕd3x

もし電子が原点に位置する電荷 e の点粒子であれば ρ=eδ(x) そして ϕ=e/4πr となり, 従って Eint は無限大となる.その代わりに, 一様な密度の球状の電荷分布を考えることが出来る.このようなモデルでは, 電子の静止エネルギーのかなりの部分を自己エネルギーに帰するのであれば, 電子の「半径」は古典的な電子半径 r0=(e2/4πmc2) のオーダーでなければならない.電子の自己エネルギーをディラックの相対論的電子理論を用いて量子力学的に計算すると, やはり無限大になるがその発散ははるかに小さいことが分かる.この点については §4-7 で触れることにする.

原子準位のシフト(Atomic level shifts)

ここでは束縛された原子電子の自己エネルギーについて議論するが, それはここまで議論して来た「輻射の量子論」の枠組みの中で行ない, 横方向の電磁場は量子化するが電子は非相対論的に扱かうこととする.原子状態 A を考える(それは必ずしも安定な基底状態とは仮定しない).束縛された電子の自己エネルギーは, 量子化された輻射場との相互作用による状態 A のエネルギーシフトとして出現する.先に強調したように, 輻射の量子論によれば, たとえ入射する輻射場が無くても原子状態 A は光子を放射することが出来るのだった.そこで原子が放射された光子を吸収して状態 A に戻ったとしよう.e2 のオーダーには図 2-7 に示すように2種類の過程が存在する.

fig 2-7 sakurai

図 2-7 束縛電子の自己エネルギー

この2つの図のうち, 相互作用ハミルトニアンの AA 相互作用項に起因する図 2-7(a) の結果は, 原子準位のシフトの観点からは面白くないことが分かる.図 2-7(b) の基本的な摂動の行列要素は, 次の AI+γ に対する「放射」行列要素:

(2.223)HIA(emis)=HIAeiωt

と, 次の I+γA に対する「吸収」行列要素:
(2.224)HAI(abs)=HAIeiωt

である.ただし, 時間依存しない行列要素 HIAHAI は, 前と同様に双極子近似を用いるとして次とする [2][訳註] 因子はFeynmanに合わせるために「Gauss単位系」に置き換えてある.
(2.225)HIA=(HAI)=c2πωemcI|pϵ(α)|A

時間依存する振幅 cIcA を求めるには, 次の連立微分方程式を解かなければならない:

(2.226)ic˙I=photonHIA(emis)cAei(EIEA)t/,(2.227)ic˙A=photonIHAI(abs)cIei(EAEI)t/,

量子化された場の演算子 A は, 真空の状態と「任意の」運動量と偏光を持った単一光子状態との間に非ゼロの行列要素を持つので, 記号 Σphoton の和は「全ての可能な運動量と偏光」について取り, それは光子エネルギーが EIEA に等しくない場合も含むことを意味している.
我々の関心事は状態 A のエネルギーシフトの式を得ることである.このため, 式 (2.226) と式 (2.227) を次のような前提(ansatz)で解くことを提案する:
(2.228)cA=exp[iΔEAt/]

この前提では, 完全な波動関数は次のように変化する:
(2.229)ψuA(x)exp[i(EA+ΔEA)t/]

従って波動関数の時間依存性は, エネルギー (EA+ΔEA) の量子力学的状態から期待されるものと正確に一致する.式 (2.228) を式 (2.226) に挿入して積分すると,
cI=1iphoton0tHIAexp[iωt]exp[iΔEAt]exp[i(EIEA)t]dt(2.230)=photonHIAexp[i(EIEAΔEA+ω)t/]1EA+ΔEAEIω

状態 A のエネルギーシフト ΔEA を求めるには, 式 (2.228) と式 (2.230) を式 (2.227) に代入する.すると左辺は,
ic˙A=i×(i)ΔEAexp[iΔEAt/]=ΔEAexp[iΔEAt/],

右辺は,
photonIHAIeiωtHIAexp[i(EIEAΔEA+ω)t/]1EA+ΔEAEIωei(EAEI)t/=photonI|HIA|2exp[iΔEAt/]exp[i(EAEIω)t/]EA+ΔEAEIω
従って両辺を exp[iΔEAt/] で割れば,
(2.231)ΔEA=photonI|HIA|21exp[i(EA+ΔEAEIω)t/]EA+ΔEAEIω

ΔEAe の2次で計算しているので, 式 (2.231) の右辺の ΔEA を削除するのは正当なことである.従って,
(2.232)ΔEA=photonI|HIA|21exp[i(EAEIω)t/]EAEIω

ここで, 摂動が作用する時間間隔 [0,t] を無限大にすることが出来る.しかし, そのままでは(as it stands) この式は t のとき振動することを示している.
それは, 式 (2.232) を次のような積分から求めたのだから当然である (式 (2.230) を参照):
(2.233)0teixtdt,wherex=(EIEAΔEA+ω)/

これは t のとき明らかに不適切な定義である.しかしながら, このような積分は x に小さな正の虚部を加えることで収束させることが出来る.このようにして, 一連の非常に便利な関係が得られる [3][訳註] 次が成り立つことを用いている (砂川重信:「量子力学」の第2章 §1 (3) デルタ関数 を参照すべし): \begin{equation*} -i\int_{0}^{\infty}e^{ixt}\,dt … Continue reading
limt1eixtx=limϵ0+(i0ei(x+iϵ)tdt)=limϵ0+1x+iϵ=limϵ0+[xx2+ϵ2iϵx2+ϵ2](2.234)=1xiπδ(x)

これらは t とした式 (2.232) に直に当てはめることが出来る [4][訳註] x=(EAEIω)/ の場合であるから, \begin{align*} \Delta E_A &=\sum_{\mathrm{photon}}\sum_I … Continue reading
するとエネルギーシフト ΔEA は, 実数部 Re(ΔEA) と虚数部 Im(ΔEA) の両方を持っていることが分かる:
ΔEA=Re(ΔEA)+iIm(ΔEA),(2.235)Re(ΔEA)=photonI|HIA|2EAEIω(2.236)Im(ΔEA)=πphotonI|HIA|2δ(EAEIω)

式 (2.235) の光子の和は全ての運動量と偏光の光子を対象としており, 式(2.235) の原子エネルギー準位の和 I もやはり「EI<EA を満たす必要がない」という意味で無制限的である (「実際の」放射過程とは対照的である).つまり, 式 (2.235) の光子の放射と吸収は一般にエネルギー保存を満たさないので,「実光子(actual photon)」の放射と吸収の過程 (§ 2-4 で議論した) とは異なるものである.このような光子は「仮想的(virtual)」と言われる.物理的な原子は, 一部の時間では「原子+仮想光子」という解離状態にあると想像しても良いであろう.エネルギーシフトの実部を引き起こすのは,「全ての可能な運動量と偏極を持った仮想光子のエネルギー非保存的な放射および吸収」に関連した相互作用エネルギーである.

式 (2.235) の光子和とは対照的に, 式 (2.236) の光子和はエネルギー保存 EA=EI+ω を満たす光子のみを対象としている. 言い換えれば, 式 (2.236) に現れる光子は「実光子」であって,「仮想光子」ではない.エネルギーシフト ΔEA が虚部を持つのは, エネルギー保存に違反することなく, 状態 A が自然放射によって状態 I に崩壊することが可能な場合だけである.より定量的には, 式 (2.236) の両辺に 2/ を掛け合わせた次式に注意する:

(2.237)2Im(ΔEA)=photonI2π|HIA|2δ(EAEIω)

しかしこの式の右辺は正に, エネルギー的に許容される全ての最終状態を合計し「黄金律」に従って計算した自然放射の遷移確率の式である [式(2.114)から式(2.117)を参照] [5][訳註] 「遷移率」すなわち単位時間当りの遷移確率 [w]i[n] は, 時間を含む1次の摂動論が有効なとき, \begin{equation*} w_{i\to … Continue reading.従って, それは状態 A の平均寿命 τA の逆数に等しい.よって, 次のような重要な結果が得られる:
(2.238)2Im(ΔEA)=1τA=ΓA,Im(ΔEA)=ΓA2

今や Im[ΔEA] の物理的意義(重要性)は明らかである.式 (2.229) に戻ると, 完全な波動関数は次で与えられることが分かる:
ψuA(x)exp[i(EA+ΔEA)t/]=uA(x)exp[i(EA+Re[ΔEA]+iIm[ΔEA])t/]=uA(x)exp[i(EA+Re[ΔEA])t+Im[ΔEA]t](2.239)=uA(x)exp[i(EA+Re[ΔEA])tΓAt2]

ここから, おなじみの結論が導かれる:不安定な状態 A を見出す確率は次式のように減衰する [6][訳註] 自発崩壊をする独立な系がたくさんあるとき, 系が崩壊せずに存続する時間の統計的平均値, 又は未崩壊の系が 1/e … Continue reading
(2.240)|ψ|2eΓAt/

要約すると, エネルギーシフト ΔEA の実数部(仮想光子の放射と吸収から生じる)は, 我々が通常「準位のずれ」と呼んでいるものであり, そして ΔEA の虚数部(実光子の放射と吸収から生じる)は, 問題としている不安定状態の崩壊幅, すなわち寿命の逆数を特徴づけるものである.また, 微分方程式 (2.181) の第2項が ΔEA の虚数部から生じると見做せるので, 前に §2-6 で述べた放射減衰の現象論的な取り扱いが今や正当化されることになるにも注意する.

以降は ΔEA の実数部のみに注目する.(これ以降 ΔEAという表現は, 実数部を意味するものとする).式 (2.225) を式 (2.235) に挿入すると次のようになる:
ΔEA=photonI|HIA|2EAEIω=photonI2πc2ω(emc)2|(pϵ(α))IA|21EAEIω
ここで光子和 photon は「全ての運動量 p (または波数 k) と偏極 α を持った光子の和」を表しているので次のように置き換える:

photon=pα  d3p(2π)3α=d3k(2π)3α

従って [7][訳註] k 空間に於ける体積積分は k=|k| そして dEγ=d(ω)=d(kc)=cdk および k2=(ω/c)2 … Continue reading,
ΔEA=pαI2πc2ω(emc)2|(pϵ(α))IA|21EAEIω=d3k(2π)3I2πc2ω(emc)2α|(pϵ(α))IA|21EAEIω=2πc2(emc)2Id3k(2π)31ω1EAEIωα|(pϵ(α))IA|2=2πc2(2π)3(emc)2IdEγω2c31ω1EAEIωdΩα|(pϵ(α))IA|2
ただし Eγ=ω=|k|c である.そこで ω2c31ω=ωc3=ω2c3=Eγ2c3 とすれば,
ΔEA=c2(2π)2e2m2c2IdEγEγ2c31EAEIωdΩα|(pϵ(α))IA|2(2.241)=1(2π)2(e2c)1(mc)2IdEγEγEAEIωdΩα|(pϵ(α))IA|2

偏光和と角度積分は, 以前の自然放射に関する議論と同じであり (式 (2.130) から式 (2.132) を参照),
(2.242)dΩα|(pϵ(α))IA|2=8π3|(p)IA|2

従って, 残るのはエネルギー積分だけである [8][訳註] この式(2.243)の2番目の因子は「微細構造因子α=e2cであることに注意する.
(2.243)ΔEA=23π(e2c)1(mc)2IEγ|(p)IA|2EAEIEγdEγ

このエネルギー積分は Eγ のゼロから無限大までの全ての可能な値にわたって行なう.しかし, 積分は明らかに線形に発散する.電子の非相対論的近似は Eγ>mc2 の光子の放出に対しては破綻せざるを得ないので, 非常に高いエネルギーの仮想光子からの寄与をあまり深刻に考えることは出来ないと主張することが出来る.従って, 次のように置き換える:
(2.244)0dEγ0Eγ(max)dEγ

ただし Eγ(max) は「カットオフ エネルギー」(cut-off energy) として知られている量である.残念ながら, この方法で計算されたエネルギーシフトは我々が導入した任意パラメータ Eγ(max) に非常に敏感である.

質量の繰り込み(Mass renormalization)

この段階で, 図 2-7(b) に示した過程による運動量 p が固定された「自由電子の自己エネルギー」について議論しておくと有利であろう.今度は自由電子なので, 状態 A と状態 I の波動関数は各々 ψA(x)ψp=eipx//VψI(x)ψp=eipx//V で表される.従って, 式 (2.235) 中の時間に依存しない摂動の行列要素 HIA は,「自由電子による光子の自発放射であるから, 式 (2.99) の Hint に於いて輻射場起因の因子 eiωt を省き nk,α をゼロと置いた式で更に i=1,B=I とした場合」である.よって,
(2.99)B;nk,α+1|Hint|A,nk,α=em(nk,α+1)2πωViB|eikxipiϵ(α)|Aeiωt, HIA=eiωtI|Hint|A=em2πωVψI(x)eikxpϵ(α)ψA(x)d3x=em2πωVeipx/Veikx/(pϵ(α))eipx/Vd3x=em2πωV(pϵ(α))Veipx/Vei(pk)x/Vd3x=em2πωV(pϵ(α))ψp|ψpk
ψs|ψs=δss であるから,

(2.245)HIA=2πωVempϵ(α)δp,pk

電子が非相対論的に扱えるエネルギー範囲に於いて, エネルギー分母は kp として次のように近似できる:
(2.246)EAEIω=p22m(pk)22mωp22mp22mω=ω

残りのステップは前と同じである.上式から, 式 (2.243) に於いて EAEIEγEγ とすれば,
ΔEA=23π(e2c)1(mc)2IEγ|(p)IA|2dEγEAEIEγ=23π(e2c)1(mc)2IEγ|(p)IA|2dEγEγ=23π(e2c)1(mc)2I|(p)IA|2dEγ=23π(e2c)1(mc)2(p2)AAdEγ=(e2c)23π(p2)AA(mc)2dEγ,(p2)AA=p2
よって,「自由電子の自己エネルギー」として次を得る:
(2.247)ΔEfree=(e2c)2p23π(mc)20Eγ(max)dEγ=Cp2,(2.248)whereC=(e2c)23πEγ(max)(mc)2=α23πEγ(max)m2c2

この式 (2.247) もやはり Eγ(max) のとき線形に発散することに注意する.

自由電子の場合の図 2-7(b) は, 単純に p2 に比例するエネルギーを生じることが分かる.電磁相互作用のスイッチを決して切ることが出来ない程度に, 通常の運動エネルギー項 p2/2m からこの追加エネルギーを分離することは出来ない.実際, 観測可能な運動エネルギーは, 図 2-7(b) のような相互作用が存在しない仮想世界の運動エネルギーと我々が計算した Cp2 との和である.一方, 電子を非相対論的な粒子と見做すと, 電子の質量はエネルギー運動量の関係式

(2.249)Ep2=12m

で与えられる.ただし E は自由粒子のエネルギーである.従って, 観測される電子質量 mobs は, 次式に従って C と関係している:
p22mobs=p22mbare+Cp2=p22mbare+2mbareCp22mbare=p22mbare(1+2mbareC)(2.250)=p22mbare(1+2mbareC)1p22mbare(12mbareC)

ただし mbare は, 図 2-7(b) で表される相互作用がない場合に測定される電子の質量を表わす.要するに, 図 2-7(b) の正味の効果は電子質量の定義を変えただけである:
(2.251)mbare  mobsmbare(12mC)=[1+(1137)43πEγ(max)mc2]mbare

我々が実験室で観測している質量は mbare ではなくて, 時には「繰り込み質量」(renormalized mass) と呼ばれることがある mobs である [9]固体物理学に詳しい読者であれば, 輻射の量子論に於ける「観測質量」または「繰り込み質量」の概念が, … Continue reading.式 (2.247) の ΔE(free) は負であるのに対し Δm=mobsmbare は正であることに注意する.式 (2.251) から Eγ(max)mc2 のオーダーで「カットオフ」した場合, 観測される電子質量の約 0.3 % が自己エネルギーに起因することが分かる.しかし, このような計算に大きな意味を見出すことはできない.なぜなら, ディラックの相対論的電子理論を用いると, 電子の自己エネルギーは全く異なる式になるからである.

ここで「原子状態の準位シフト」に話を戻そう.Schrödinger 程式を解いて原子のエネルギー準位を計算するときに用いる運動エネルギーは p2/2mobs であり, これにはすでに補正項 Cp2 が含まれている!.観測可能なエネルギーシフトを見積もる際には, すでに考慮に入れたエネルギー部分を必ず差し引かなければならない.なぜなら, 無摂動のエネルギー準位を計算する際には p2/2mobs ではなく p2/2mbare を用いるからである.これが H.A.Kramers によって提唱された「質量繰り込み」の考え方である.定量的には, 観測可能なエネルギーシフト ΔEA(obs) は次式で与えられる:

ΔEA(obs)=ΔEAA|(p22mobsp22mbare)|A=ΔEAA|Cp2|A(2.252)=ΔEACA|p2|A

言い換えれば, 観測可能なのは「束縛電子の自己エネルギーと自由電子の自己エネルギーの差」である.

References

References
1 Jackson (1962), pp. 578-597.
2 [訳註] 因子はFeynmanに合わせるために「Gauss単位系」に置き換えてある.
3 [訳註] 次が成り立つことを用いている (砂川重信:「量子力学」の第2章 §1 (3) デルタ関数 を参照すべし):
i0eixtdt=eixtx|0=1x(eixtx)t=limt1eixtx,δ(x)=1πlimϵ0ϵx2+ϵ2
4 [訳註] x=(EAEIω)/ の場合であるから,
ΔEA=photonI|HIA|2(1xiπδ(x))=photonI|HIA|2(EAEIωiπδ(EAEIω))=photonI|HIA|2EAEIωiphotonIπ|HIA|2δ(EAEIω)
5 [訳註] 「遷移率」すなわち単位時間当りの遷移確率 [w]i[n] は, 時間を含む1次の摂動論が有効なとき,
wi[n]=2π|Vni|2δ(EnEi)

と書ける.この式は「フェルミの黄金則」と呼ばれ, 実際上非常に重要である.
6 [訳註] 自発崩壊をする独立な系がたくさんあるとき, 系が崩壊せずに存続する時間の統計的平均値, 又は未崩壊の系が 1/e になる時間を「平均寿命」と言う.平均寿命は「崩壊定数」の逆数である.統計理論によると dt 時間中に原子核が崩壊する確率は λdt である.最初に崩壊していない原子核が N0 個あったとき dt 時間内に崩壊する数 dNdN=λdtN で与えられるから, これを積分すると時刻 t に存在する原子数 NN=N0eλt となる.このときの λ を「崩壊定数」という.
7 [訳註] k 空間に於ける体積積分は k=|k| そして dEγ=d(ω)=d(kc)=cdk および k2=(ω/c)2 として次のように書ける:
d3k=k2dkdΩ=k2dkdΩ=k2dEγcdΩ=dEγω2c3dΩ
8 [訳註] この式(2.243)の2番目の因子は「微細構造因子α=e2cであることに注意する.
9 固体物理学に詳しい読者であれば, 輻射の量子論に於ける「観測質量」または「繰り込み質量」の概念が, 固体の格子振動と相互作用する電子の「有効質量」(effective mass)の概念に似ていることにお気付きだろう.