水素原子スペクトルの微細構造


ファインマンは, 第 9 講 (Ninth Lecture) の最初の問題の直前で次のように書いていた:

次に E=mc2+W と置く.ただし Wmc2 である.そして V=Ze2/r を代入すると,

(9-2”)(WV)χ+22mχ=(WV)22mχ

左辺第1項との比較から右辺の項を無視すると, 通常のシュレディンガー方程式となる.(WV)2/2mc2 を摂動ポテンシャル ΔV とすることで, 学生諸君は水素原子スペクトルの微細構造を得て, それを正確な数値と比較すべきである.

そこで, シュポルスキー:「原子物理学」§ 204 微細構造の公式 の文章や, シッフ:「量子力学」の文章の抜粋を示すことで, この課題の解答としようと思う.

微細構造の公式

ケプラー問題をディラック方程式を用いて解くと, バルマーのエネルギー準位以外にその補正項が得られる.この補正項は「微細構造の公式」と呼ばれている.この補正項をディラック方程式に依らず, 初等的なやり方で導出してみる.
シュレディンガー方程式によるケプラー問題を解く際に, さらに次の2つの重要な因子を考慮する必要がある:

  1. 質量の速度に対する相対論的依存性.
  2. スピンの存在.

これら2つの因子の各々に対して, この問題の古典論的ハミルトニアン関数に於ける付加項が対応する.

  1. 質量の速度に対する依存性の考慮に関連した付加項の起源は,「電子の相対論的運動は,ケプラー楕円運動 H0 に歳差運動 Hr が付随したものとなる」ことによる:
    H=H0+Hr,H0=12μ0(pr2+1r2pϕ2)Ze2r

    ただし μ0 は静止質量であり, Hr は相対論的歳差運動に対応した付加的なエネルギーである.
  2. スピンと関連した磁気モーメントの存在のために, 電子はその「軌道運動」によって生じている磁場中に在る「磁気双極子」として振る舞う.即ち磁気的相互作用の補足的なエネルギーを持っている.この付加項を「スピン-軌道相互作用」のエネルギー Hls という.

従って, 全ハミルトニアン関数は次の3項から構成される:

H=H0+Hr+Hls

その際に H0 はエネルギーの主要部分即ち「バルマー項」を与える.そして HrHls は小さい「摂動」と考えることが出来, これがスペクトル項 T0 に補正 ΔT を与える[1][ブログ註] 「スペクトル項」 T は,「波数」(波長の逆数)を表わす量である.振動数を ν, 光速度を c とするとき, 光の波長 λ は … Continue reading
T0=T(n)Enhc,T0+ΔT,ΔT=ΔTr+ΔTls

まず, 相対論的な補正 ΔTr を摂動論を用いて計算すると次の公式が得られる:
(204.1)ΔTr=Rα2Z4n3(1l+1234n)

ただし R は「リュードベリ定数」,そしてαは「微細構造定数」である:
R=me44π3c=109,737.303cm1,α=e2c=1137

さて, 今度は第2の補正, すなわち「スピンと軌道角運動量との間の相互作用による補正」の考察に向かおう.その補正を半古典論的な考察によって求めよう.電子の軌道運動によって作られる磁場中に於ける磁石としての電子が持つ付加的なエネルギーを,「磁場中に置かれた磁気双極子の位置エネルギー V=MH」として計算しよう.ここで M は電子の固有磁気モーメントであり, H は電子の軌道運動によって生じる磁場の強さである.磁場の強さの値 H は次のようにして計算できる.核と結び付いた「動かない」座標系から, 電子と一緒に動く座標系に移ろう.この座標系では電子の重心は静止しており, 核 Ze は電子の速度 v と絶対値は等しいが方向は反対の速度 v で運動している.この運動は強さ Zev の電流を生じる (図 3.1参照).「ビオ-サバールの法則」から,
H(x)=1cj(x)×(xx)|xx|3d3x,orH=1cj(x)×RR3dVH=1cj×RR3dV=(jRdV)×RcR2,jRdV=Ze(v)RH=(Zev)R×RcR2=ZecR3v×R
よって電子が存在する場所 r=R に於いて, 電流 (Zev) が生起する磁場は次である:

(204.3)H=Zecr3v×r=Zecr3r×v

スピン-軌道相互作用

図 3-1. 電子から見たとき, 核 Ze は電子に対し周回しているので, 磁場 H を電子の位置に生起する.

ベクトル積 r×v は, 容易に角運動量 L を使って表される.つまり, 質量 μ の電子の角運動量は

(204.4)L=μr×v

であることを考慮に入れると, 次が得られる:
(204.5)H=Zecr3μμr×v=Zeμcr3L

H 中の磁気モーメント M の双極子の持つエネルギーは MH であるから, 付加的エネルギーは次で与えられる:
(204.6)MH=M(Zeμcr3L)=Zeμcr3ML

この式は, 同時に電子のラーモア歳差運動の運動エネルギーにも等しくなっている.しかし, 式 (204.6) で表されているエネルギーは「電子の重心」が静止している様な座標系に於けるラーモアの歳差運動に対応していることを念頭に置く必要がある.「核」が静止している様な座標系, あるいはもっと正確に言えば「原子全体の重心」が静止している様な座標系に戻るためには, もう一度ローレンツ変換をする必要がある.最終結果として得られる磁気的相互作用の付加エネルギーは, J.I.フレンケル が示したところによれば, 式 (204.6) の半分に等しい.すなわち,
(204.7)Hls=Ze2μcr3ML

この付加エネルギーを, ハミルトニアン関数の主要部分 H0 に対する小さな「摂動」と考えることが出来る.摂動理論によれば, 摂動によって生じる付加エネルギーはハミルトン関数の摂動項を無摂動状態で平均した「平均値」に等しい.すなわち,
ΔEls=Hls

式 (204.7) に含まれている量は r を除けば全て定数であるから, 平均はただ 1/r3 についてだけ取ればよい:
(204.8)ΔEls=Ze2μc(1r3)ML

平均は無摂動状態についてとられるのであるから, その計算の実行のためには, ケプラー問題の固有関数を用いる必要がある[2][ブログ註] ケプラー問題は, シュポルスキーの § 184 に解説されている.
(1r3)=ψnl1r3ψnldτ

その計算の結果, 次式が得られる:
(204.9)(1r3)=Z3a13n3l(l+12)(l+1)

ただし a1 は「第1ボーア半径」であって, 次に等しい:
(204.10)a1=2μe2=0.528×108 cm

さて今度は, 式 (204.8) に含まれているスカラー積 ML の計算に向かおう.スピンに対応する磁気モーメント M とスピン角運動量 S との比は 2×e2μc=eμc に等しい.M の方向は, 電子の電荷がマイナスであるから S とは反対向きである.よって,
(204.11)M=eμcS

スピン角運動量および軌道角運動量の数値は次に等しい:
(204.12,13)|S|=s(s+1),|L|=l(l+1)

記法を簡潔にするために, 今後途中の計算では s(s+1) 及び l(l+1) をそれぞれ s 及び l と書くことにしよう.式 (204.9) 及び式 (204.11) を式 (204.8) に代入し, 式 (204.10) を a1 に用い, 次の「リュードベり定数R と「微細構造定数α を持ち込む [3][ブログ註] 「ゾンマーフェルトの微細構造定数」と呼ばれているこの定数 α は, 極めて注目すべきものである.3つの基本的な普遍定数 e, … Continue reading
R=μe44π3c,α=e2c

すると簡単な運算の後, 次を得るであろう:
(204.14)ΔEls=Rα2cZ42πn3(l+12)(l+1)lscosθ

今や最後の段階, すなわち軌道角運動量とスピン角運動量との間の角 θ の 余弦 cosθ の計算が残るだけとなった.この計算のためには角運動量 LS との合成を考慮する必要がある.

図 3-2. 軌道角運動量ベクトル L とスピン角運動量ベクトル S は (a) 平行, または (b) 反平行に加えることができる.全角運動量ベクトル J=L+S の大きさは J=j(j+1) であり,
(a) j=l+1/2 および (b) j=l1/2 である.

角運動量 L 及び S は決して平行にもならず, それらの間の角度は 0 とも π とも異なることが結論される.次の様な半古典論的な考察方法を用いることが出来る.LS とは普通の平行四辺形の規則に従って合成され, その結果, 全角運動量 J が得られる.ところがベクトル LS とはそれらに対応している磁場を媒介として結び付いているから, それらは J に対して, ちょうど弾性的な糸で結ばれている2つのジャイロスコープが, その不変な全角運動量の方向の周りに歳差運動をするのと同じように歳差運動をする.この純然たる「古典論的な」考察の中に,「量子力学」にとって特徴的な次の様な補正が持ち込まれる.第一に LS の「方向量子化」の要求によって制限される.すなわちベクトル LJ に対して, ただ J の方向の L の成分が ml に等しくなる様な方向しか取ることが出来ず, ベクトル S は, 場の方向の成分が ms に等しくなる様な方向しか取ることが出来ない.これによって LS のなす角度 θ は飛び飛びの選ばれた値の一組に制限される.第二に, 全角運動量ベクトル J の大きさは j に等しくはなくて j(j+1) に等しい.ただし j=l±1/2 である.つまり j は全角運動量の量子数なのである.このような考察方法は、明らかに論理的に不完全であるにも関わらず, 驚くべき正確さで正しい結果をもたらす.

図 3-3. 角運動量ベクトル LS は, 全角運動量ベクトル J の周りを歳差運動する.
右図は, 全角運動量べクトル J の空間量子化の様子.


図 3-4. 種々の j の値に対する全角運動量ベクトル J の可能な方位.

以上のことを考慮に入れつつ, 今や我々は容易に cosθ を計算出来る.三角法の定理によって,

j2=l2+s22lscos(πθ)=l2+s2+2lscosθ

よって,
(204.15)lscosθ=j2l2s22=12{j(j+1)l(l+1)s(s+1)}

この結果を式 (204.14) に代入すると,
ΔEls=Rα2cZ42πn3l(l+12)(l+1)j(j+1)l(l+1)s(s+1)2

が得られる.従って「スピン-軌道相互作用」によるスペクトル項の補正は次となる:
(204.16)ΔTls=ΔEls2πc=Rα2Z4n3l(l+12)(l+1)j(j+1)l(l+1)s(s+1)2

j の値は j=l+1/2 及び j=l1/2 のみをとるから,
j(j+1)l(l+1)s(s+1)2={l2:j=l+12l+12:j=l12

従って, 式 (204.16) は次のような2つの値を持つ:
(204.17)ΔTls={Rα2Z4n32(l+12)(l+1):j=l+12+Rα2Z4n32l(l+12):j=l12

スペクトル項の全補正を求めるには, 式 (204.1) による相対論的補正 ΔTr と式 (204.17) による補正 ΔTls との和を見出すべきである.その結果は次となる:

  1. j=l+1/2 に対して( l=j1/2),
    ΔT=ΔTr+ΔTls=Rα2Z4n3(1l+1234n)Rα2Z42n3(l+12)(l+1)=Rα2Z4n3(1l+134n)=Rα2Z4n3(1j12+134n)
  2. j=l1/2 に対して (l=j+1/2),
    ΔT=ΔTr+ΔTls=Rα2Z4n3(1l+1234n)+Rα2Z4n32l(l+12)=Rα2Z4n3(1l34n)=Rα2Z4n3(1j+1234n)

これらは j=l±12 であることに留意すると一つにまとめることが出来る.そのとき両方のカッコ中の第1項の分母は共に j+12 となるので, 次の結果となる:

(204.18)ΔT=Rα2Z4n3(1j+1234n)

この式から j がただ一つの値を持っているだけの S 項 (n=1 の場合) に対しては, 条件 [nj12,j12] から, 補正はただ一つの値 (j=1/2) しか持たないことが分かる:
ΔT=Rα2Z4n3(112+1234n)=Rα2Z4n3(134n)

従って補正は S 項をずらせるだけで, それを分岐させないことが分かる.その他 (n=2,3,4,) の項 (P,D,F,) に対しては j は2つずつの値 l±1/2 を持ち, 従ってこれらの項に対応する各エネルギー準位は2つの副準位に分かれる.更に, 主量子数が n=2 の際には2つの準位 2S 及び 2P があり, 分岐は3つの副準位 (1つのS 項と2つの P項) を与える.しかしこのとき, 条件 [nj12,j12]から, jS 項 (l=0) に対しては値 1/2 を持ち, P 項 (l=1) に対しては 1/2 及び 3/2 を持つ.このとき, S 項の副準位 2S1/2(j=0+1/2=1/2)P 項の副準位 2P1/2(j=11/2=1/2) とは互いに一致するが, 他方 2P1/2 準位と 2P3/2 準位 (j=1+1/2=3/2) とは分岐してしまう.

式 (204.18) から 2P3/2 準位と 2P1/2 準位を計算してみよう.「リュードベり定数R と「微細構造定数α の値を次とする:

R=μe44π3c=109,737.303cm1,α=e2c=1137,α2=5.32×105

すると各準位の補正量 ΔT は,「波数」として次となる:

  1. 2P3/2の場合(n=2,l=1,j=1+12=32)ΔT3/2=Rα28(1234×2),
  2. 2P1/2の場合(n=2,l=1,j=112=12)ΔT1/2=Rα28(134×2),

従って, 2つの分岐量 (補正量の差) は次となる:
Δ=ΔT1/2ΔT3/2=Rα28×(112)=1.09737×105×5.32×1058×12=0.364cm1


[ 参考 ] 上式 (204.1) を示すために, シッフ:「量子力学」§ 51 の文章の抜粋を示す.
電磁ポテンシャル A(r,t), ϕ(r,t) が存在する場合, これを波動方程式の中に取り入れるためには ϕ(1/c)A とが E 及び p と同一のLorentz変換性を持っていることを利用してやればよい.非相対論的な式 (24.29)

(24.29)H=12m(pecA)2+eϕ

の形から類推して, 電荷 e の粒子に対しては式 (51.2):E2=c2p2+m2c4 を,
(51.9)(Eeϕ)2=(cpeA)2+m2c4

で置き換える.このときの p は一般化運動量 L/v の方であり, 通常の3次元運動量mvではない[4][ブログ註] … Continue reading.そこで, 置き換え Eit, pi をすると
(51.10)(22t22ieϕtieϕt+e2ϕ2)ψ=[2c22+2iecA+iec(A)+e2A2+m2c4]ψ

ポテンシャル A,ϕ が時間に関係しない場合には, 波動関数 (51.10) は rt とに関して分離することが出来る.即ち, その時には
ψ(r,t)=u(r)eiEt/

と置いて式 (51.10) に代入する.左辺だけが変わり itu(r)eiEt/=Eu(r)eiEt/ となるので,
(51.12)(Eeϕ)2u=[2c22+2iecA+iec(A)+e2A2+m2c4]u

さて次は「A=0 で, かつ ϕ(r) が球対称, つまり ϕ(r)=ϕ(r) という特別な場合」を考えてみる.すると式 (51.12) は
(51.12′)(2c22+m2c4)u(r)=[Eeϕ(r)]2u(r)

となる.これは球座標を使って変数分離が可能であって次となる(シッフ:「量子力学」§14を見よ):
u(r,θ,ϕ)=R(r)Ylm(θ,φ),(51.14)[1r2ddr(r2ddr)+l(l+1)r2]R=(Eeϕ)2m2c42c2R,l=0,1,2,

ここで E=mc2+E と置き, さらに Eeϕ とは mc2 に比べて無視できるものと仮定すると, この式は動径 r に対する非相対論的な方程式に帰着する.そのとき式 (51.14) の右辺の分数部分は当然期待されるように 2m2(Eeϕ) になる[5][ブログ註] 式 (51.12′) を変形すると, \begin{equation*} -\nabla u(\mb{r}) = \frac{(E-e\phi)^{2}-m^{2}c^{4}}{\hbar^{2}c^{2}}u(\mb{r}),\quad\rightarrow\quad -\frac{\hbar^{2}}{2m}\nabla … Continue reading

さらに水素様原子の場合すなわちクーロン場でのエネルギー準位を考えてみる.ポテンシャルを eϕ=Ze2/r としたときの式 (51.14) の正確な解は, §16 の結果を使って容易に求められる[6][ブログ註] シュレディンガー方程式を球座標そして相対座標の表現にすると, 角運動量の量子数 l … Continue reading.式 (51.10) で表される粒子はスピンを持たないから, もちろん電子では有り得ないが, このことを度外視すればこの状況は水素様原子を表していると考えて良い.
ρ=αrと置くと, 式 (51.14) は次のようになる:

1ρ2ddρ(ρ2dRdρ)+[λρ14l(l+1)γ2ρ2]R=0,(51.15)whereγZe2c,α24(m2c4E2)2c2,λ2Eγcα

l(l+1)l(l+1)γ2 で置き換えられている点以外は, これはまさしく水素原子に対する式 (16.7) と一致している.パラメータ λρ= での R の境界条件で決まり, 式 (51.15) の最後の2式から α を消去すると, Eλ によって次のように表される:
(51.16)E=mc2(1+γ2λ2)1/2

式 (16.7) の解き方を調べれば分かることだが, ρ=0ρ= とで式 (51.15) の有限な解が存在するのは
(51.17)λ=n+s+1

のときに限られる.ただし n はゼロか正の整数で, s は2次方程式
(51.18)s(s+1)=l(l+1)γ2

の解のうち負でないものである.この方程式は2つの解
(51.19)s=12±12(2l+1)24γ2

を持ち, l>0 に対してそのうち一つは正, 他方は負である.

途中を少し省略する

そこで l の全ての価に対して + の符号を取ると, 式 (51.17) から次となる:

(51.20)λ=n+12+(l+12)2γ2

式 (51.16) と式 (51.20) とは式 (16.15) の非相対論的エネルギー準位に対する「微細構造」を与える.これを調べるには, エネルギー準位を表す式を γ2 のベキで展開すればよい.この結果を γ4 まで取った結果は次のようになる:
(51.21)E=mc2[1γ22n2γ42n4(nl+1234)](51.21′)W=Emc2=mc2γ22n2mc2γ42n4(nl+1234)

ただし, この γ は「微細構造定数」と呼ばれ, その大きさは γ1/137 である.また n=n+l+1 は全量子数で, 正の整数値を取ることが出来る.式 (51.21) の右辺の第1項は静止エネルギーである.第2項は水素原子に対するエネルギー準位の式 En と一致する:
En0=mc2γ22n2=mZ2e422n2

第3項は相対論的な補正 ΔEr であり, n が同じで l が異なる状態の間に有った縮退が除かれる.式 (51.21) をスペクトル項の公式に変形すると, 上式 (204.1) に一致する:
ΔTr=ΔErhc=mc4πn4(Ze2c)4(nl+1234)=me44π3c(e2c)2Z4n4(nl+1234)(204.1)=Rα2Z4n4(nl+1234)=Rα2Z4n3(1l+1234n)


References

References
1 [ブログ註] 「スペクトル項T は,「波数」(波長の逆数)を表わす量である.振動数を ν, 光速度を c とするとき, 光の波長 λc=νλ の関係にある.光子のエネルギーを E=hν とすると,
1λ=νc=hνhc=Ehc,T(n)Enhc cm1

原子中のエネルギー準位は負の値となるので, スペクトル項は上の波数 1/λ の式中のエネルギーにマイナス符号を付けて定義しているのである.
2 [ブログ註] ケプラー問題は, シュポルスキーの § 184 に解説されている.
3 [ブログ註] 「ゾンマーフェルトの微細構造定数」と呼ばれているこの定数 α は, 極めて注目すべきものである.3つの基本的な普遍定数 e, , 及び c の組み合わせであるこの定数が「ディメンションのない数」であることは注目に値する.
実際, e2 はディメンション [erg×cm] を持っている.これは次のことから分かる.e2/r はクーロンの位置エネルギーであり, 従ってエルグを単位として表される.それ故, [e2]=[エネルギー×長さ], すなわち [erg×cm]である.c のディメンションは [エネルギー×時間×速度] である.すなわち
[c]=[erg×sec×cmsec]=[erg×cm].

結局 ce2 と同じディメンションを持っている.すなわち α はディメンションの無い量である.α の数値は e, , c の数値から求められ, 次のものに等しい:
α=7.0×102,α2=5.32×105,1α=137.03,
4 [ブログ註] このような置き換えをして,電磁場との相互作用を導入することを「minimalな電磁相互作用を導入する」という.minimal結合を導入する相対論的な処方は, 例えば Claude Itzykson, Jean-Bernard Zuber, Quantum Field Theory の §2-2-3 を参照すべし.
5 [ブログ註] 式 (51.12′) を変形すると,
(1)u(r)=(Eeϕ)2m2c42c2u(r),22mu=(Eeϕ)2m2c42mc2u

ここで E=mc2+E と置くと右辺の分子は,
(Eeϕ)2m2c4=(E+mc2eϕ)2m2c4=m2c4+2mc2(Eeϕ)+(Eeϕ)2m2c4=2mc2(Eeϕ)+(Eeϕ)2
従って, 上式 (1) は
22mu=12mc2{2mc2(Eeϕ)+(Eeϕ)2}u=(Eeϕ)u+(Eeϕ)22mc2u[22m+eϕ(Eeϕ)22mc2]u=Eu
ここで uχ, EW と書き換え, そして V=eϕ と置くならば,
[22m+V(WV)22mc2]χ=Wχ,(WV)χ+22mχ=(WV)22mc2χ

となり, これは本文の式 (9-2”) に一致することが分かる.
6 [ブログ註] シュレディンガー方程式を球座標そして相対座標の表現にすると, 角運動量の量子数 l に対する動径方向の方程式は次となる:
(16-6)22μ1r2ddr(r2dRdr)Ze2rR+22μl(l+1)r2R=ER

さらに, 無次元の独立変数 ρ=αr となるように次のような αλ を導入する:
α2=8μ|E|2α=8μ|E|2,λ=2μZe2α2=Ze2μ2|E|

すると方程式 (16-6) はディメンションの無い形に書き直される:
(16-7)1ρ2ddρ(ρ2dRdρ)+(λρ14l(l+1)ρ2)R=0