磁石を近づけると、原子が出す光の色がわずかに変わります。磁石が電子に仕事をしないはずなのに、なぜエネルギーが変わるのでしょうか。この記事は、その問いを量子力学の側から考えます。磁場の中の電子がどんなエネルギーをとるか、磁場がゼロの場所でも磁束の影響が残る不思議な現象、そして水素原子のスペクトルに現れる細かい分裂までを、順に計算します。読み終えると、これらの量を自分で導けるようになります。
電磁場の中の荷電粒子のハミルトニアン は、運動量を p → π = p − q A \boldsymbol{p}\to\boldsymbol{\pi}=\boldsymbol{p}-q\boldsymbol{A} p → π = p − q A に置き換える 最小結合 (電磁場との相互作用を、運動量の置き換えだけで書く方法)で与えられます。ゲージ変換は波動関数の位相 e i q χ / ℏ e^{iq\chi/\hbar} e i q χ /ℏ の掛け算で吸収され、物理量は π \boldsymbol{\pi} π を通してゲージ不変に書けます。
π \boldsymbol{\pi} π の成分どうしは [ π x , π y ] = i ℏ q B [\pi_x,\pi_y]=i\hbar qB [ π x , π y ] = i ℏ q B という交換関係を満たすため、一様磁場中の平面運動は調和振動子と同型になり、エネルギーは ℏ ω c ( n + 1 2 ) \hbar\omega_c(n+\tfrac12) ℏ ω c ( n + 2 1 ) (Landau 準位)、各準位の縮退度(同じエネルギーをもつ状態の数)は面積あたり ∣ q B ∣ / h |qB|/h ∣ q B ∣/ h です。
磁場がゼロの領域でも、そこに囲まれた磁束 があると干渉が変わります(Aharonov–Bohm 効果)。リング上の粒子では、エネルギーが E n = ℏ 2 2 m R 2 ( n − Φ / Φ 0 ) 2 E_n=\frac{\hbar^2}{2mR^2}(n-\Phi/\Phi_0)^2 E n = 2 m R 2 ℏ 2 ( n − Φ/ Φ 0 ) 2 と磁束に周期的に依存します。ベクトルポテンシャル は単なる計算の道具ではありません。
Wigner–Eckart の定理は、テンソル演算子の行列要素を Clebsch–Gordan 係数 と「換算行列要素」(m m m によらない定数の部分)の積に分けます。これで選択則と、Zeeman 効果の Landé の g g g 因子 g J = 1 + j ( j + 1 ) − l ( l + 1 ) + s ( s + 1 ) 2 j ( j + 1 ) g_J=1+\frac{j(j+1)-l(l+1)+s(s+1)}{2j(j+1)} g J = 1 + 2 j ( j + 1 ) j ( j + 1 ) − l ( l + 1 ) + s ( s + 1 ) が出ます。
水素の微細構造は、相対論的補正・スピン軌道相互作用・Darwin 項の和で、E n j ( 1 ) = − m c 2 α 4 2 n 4 ( n j + 1 / 2 − 3 4 ) E_{nj}^{(1)}=-\frac{mc^2\alpha^4}{2n^4}\left(\frac{n}{j+1/2}-\frac34\right) E nj ( 1 ) = − 2 n 4 m c 2 α 4 ( j + 1/2 n − 4 3 ) と j j j だけに依存します。
この記事で使う事前知識
ゲージ変換 : 電場と磁場を変えずに、ポテンシャル A , ϕ \boldsymbol{A},\phi A , ϕ を取り替える操作。→ ゲージ変換
電磁ポテンシャル : 電場と磁場を生み出すもとになるスカラー関数 ϕ \phi ϕ とベクトル関数 A \boldsymbol{A} A の組。→ 電磁ポテンシャル
全角運動量 : 軌道角運動量とスピンの和 J \boldsymbol{J} J 。→ 全角運動量
Clebsch–Gordan 係数 : 2 つの角運動量を合成したときの、状態の組み合わせ方を決める係数。→ 角運動量とスピン
縮退のある摂動論 : 同じエネルギーの状態が複数あるときに、摂動で準位がどう分かれるかを調べる方法。→ 摂動論
古典力学では、磁場がかかった荷電粒子に働くローレンツ力 q v × B q\boldsymbol{v}\times\boldsymbol{B} q v × B は速度に垂直で、仕事をしません。それにもかかわらず、量子力学では磁場がエネルギー準位を動かします。原子を磁場に入れるとスペクトル線が分かれる Zeeman 効果が、その最も古い例です。そして Aharonov–Bohm 効果では、磁場がゼロの場所にいる粒子の干渉縞が、遠く離れた磁束に反応します。
この謎の鍵は、量子力学ではハミルトニアンが電場・磁場ではなく ポテンシャル A , ϕ \boldsymbol{A},\phi A , ϕ で書かれることにあります。ハミルトン力学の章 で見たように、荷電粒子の正準運動量 p \boldsymbol{p} p は m v m\boldsymbol{v} m v とは異なり、ゲージによって値が変わります(例 6.3[ハミルトン形式の力学] )。この章では、まず量子論の側でゲージ不変性を確かめ、その帰結を Landau 準位と Aharonov–Bohm 効果で見ます。後半では、原子の準位を調べる道具として Wigner–Eckart の定理を証明し、Zeeman 効果、Stark 効果、水素の微細構造を計算します。
質量 m m m 、電荷 q q q の粒子の、電磁ポテンシャル ( ϕ , A ) (\phi,\boldsymbol{A}) ( ϕ , A ) の中でのハミルトニアンを
H = 1 2 m ( p − q A ( r , t ) ) 2 + q ϕ ( r , t ) , p = − i ℏ ∇ H=\frac{1}{2m}\left(\boldsymbol{p}-q\boldsymbol{A}(\boldsymbol{r},t)\right)^2+q\phi(\boldsymbol{r},t),\qquad \boldsymbol{p}=-i\hbar\nabla H = 2 m 1 ( p − q A ( r , t ) ) 2 + q ϕ ( r , t ) , p = − i ℏ∇
とします(スピンは §5 で扱う)。電磁ポテンシャルの定義と、ゲージ変換 で E , B \boldsymbol{E},\boldsymbol{B} E , B が変わらないことは、電磁気学の章のとおりです(定義 4.1[電磁ポテンシャルとゲージ変換] )。
定義 2.1 (運動量(kinetic momentum) ) ↗
π = p − q A \boldsymbol{\pi}=\boldsymbol{p}-q\boldsymbol{A} π = p − q A を 運動量 (古典の m v m\boldsymbol{v} m v に対応する)と呼ぶ。ハミルトニアンは H = π 2 / 2 m + q ϕ H=\boldsymbol{\pi}^2/2m+q\phi H = π 2 /2 m + q ϕ である。
定理 2.2 (量子力学のゲージ共変性 ) ✓
スカラー関数 χ ( r , t ) \chi(\boldsymbol{r},t) χ ( r , t ) によるゲージ変換 A ′ = A + ∇ χ \boldsymbol{A}'=\boldsymbol{A}+\nabla\chi A ′ = A + ∇ χ 、ϕ ′ = ϕ − ∂ t χ \phi'=\phi-\partial_t\chi ϕ ′ = ϕ − ∂ t χ に対して、ψ \psi ψ が i ℏ ∂ t ψ = H ψ i\hbar\partial_t\psi=H\psi i ℏ ∂ t ψ = H ψ を満たすことと、ψ ′ = e i q χ / ℏ ψ \psi'=e^{iq\chi/\hbar}\psi ψ ′ = e i q χ /ℏ ψ が A ′ , ϕ ′ \boldsymbol{A}',\phi' A ′ , ϕ ′ のハミルトニアン H ′ H' H ′ について i ℏ ∂ t ψ ′ = H ′ ψ ′ i\hbar\partial_t\psi'=H'\psi' i ℏ ∂ t ψ ′ = H ′ ψ ′ を満たすことは同値である。さらに π ′ ψ ′ = e i q χ / ℏ π ψ \boldsymbol{\pi}'\psi'=e^{iq\chi/\hbar}\boldsymbol{\pi}\psi π ′ ψ ′ = e i q χ /ℏ π ψ であり、確率密度 ∣ ψ ∣ 2 |\psi|^2 ∣ ψ ∣ 2 と速度の期待値 ⟨ π ⟩ / m \langle\boldsymbol{\pi}\rangle/m ⟨ π ⟩ / m はゲージによらない。
証明(定理 2.2) 積の微分で p ( e i q χ / ℏ ψ ) = e i q χ / ℏ ( p + q ∇ χ ) ψ \boldsymbol{p}\left(e^{iq\chi/\hbar}\psi\right)=e^{iq\chi/\hbar}\left(\boldsymbol{p}+q\nabla\chi\right)\psi p ( e i q χ /ℏ ψ ) = e i q χ /ℏ ( p + q ∇ χ ) ψ なので、
π ′ ψ ′ = ( p − q A − q ∇ χ ) e i q χ / ℏ ψ = e i q χ / ℏ ( p + q ∇ χ − q A − q ∇ χ ) ψ = e i q χ / ℏ π ψ \boldsymbol{\pi}'\psi'=(\boldsymbol{p}-q\boldsymbol{A}-q\nabla\chi)\,e^{iq\chi/\hbar}\psi=e^{iq\chi/\hbar}(\boldsymbol{p}+q\nabla\chi-q\boldsymbol{A}-q\nabla\chi)\psi=e^{iq\chi/\hbar}\boldsymbol{\pi}\psi π ′ ψ ′ = ( p − q A − q ∇ χ ) e i q χ /ℏ ψ = e i q χ /ℏ ( p + q ∇ χ − q A − q ∇ χ ) ψ = e i q χ /ℏ π ψ です。これを 2 回使うと π ′ 2 ψ ′ = e i q χ / ℏ π 2 ψ \boldsymbol{\pi}'^2\psi'=e^{iq\chi/\hbar}\boldsymbol{\pi}^2\psi π ′2 ψ ′ = e i q χ /ℏ π 2 ψ です。時間微分は i ℏ ∂ t ψ ′ = e i q χ / ℏ ( i ℏ ∂ t ψ − q ∂ t χ ψ ) i\hbar\partial_t\psi'=e^{iq\chi/\hbar}\left(i\hbar\partial_t\psi-q\,\partial_t\chi\,\psi\right) i ℏ ∂ t ψ ′ = e i q χ /ℏ ( i ℏ ∂ t ψ − q ∂ t χ ψ ) で、q ϕ ′ ψ ′ = e i q χ / ℏ ( q ϕ − q ∂ t χ ) ψ q\phi'\psi'=e^{iq\chi/\hbar}(q\phi-q\partial_t\chi)\psi q ϕ ′ ψ ′ = e i q χ /ℏ ( q ϕ − q ∂ t χ ) ψ なので、i ℏ ∂ t ψ ′ − H ′ ψ ′ = e i q χ / ℏ ( i ℏ ∂ t ψ − H ψ ) i\hbar\partial_t\psi'-H'\psi'=e^{iq\chi/\hbar}\left(i\hbar\partial_t\psi-H\psi\right) i ℏ ∂ t ψ ′ − H ′ ψ ′ = e i q χ /ℏ ( i ℏ ∂ t ψ − H ψ ) です。両辺が同時に 0 0 0 になるので同値です。残りの主張は、位相因子が ∣ ψ ∣ 2 |\psi|^2 ∣ ψ ∣ 2 を変えないことと、π ′ ψ ′ \boldsymbol{\pi}'\psi' π ′ ψ ′ の変換則から従います。
∎
ゲージ変換は粒子ごとの座標依存の位相回転です。確率分布 ∣ ψ ∣ 2 |\psi|^2 ∣ ψ ∣ 2 、運動量 π \boldsymbol{\pi} π の期待値は変わらないので、物理量はゲージによりません。一方、p \boldsymbol{p} p の期待値は A \boldsymbol{A} A の選び方に依存します。
命題 2.3 (運動量の交換関係 ) ✓
[ π i , π j ] = i ℏ q ( ∂ i A j − ∂ j A i ) = i ℏ q ε i j k B k [\pi_i,\pi_j]=i\hbar q\left(\partial_iA_j-\partial_jA_i\right)=i\hbar q\,\varepsilon_{ijk}B_k [ π i , π j ] = i ℏ q ( ∂ i A j − ∂ j A i ) = i ℏ q ε ij k B k 。
証明(命題 2.3) [ p i , p j ] = 0 [p_i,p_j]=0 [ p i , p j ] = 0 、[ A i , A j ] = 0 [A_i,A_j]=0 [ A i , A j ] = 0 より [ π i , π j ] = − q [ p i , A j ] − q [ A i , p j ] = − q [ p i , A j ] + q [ p j , A i ] [\pi_i,\pi_j]=-q[p_i,A_j]-q[A_i,p_j]=-q[p_i,A_j]+q[p_j,A_i] [ π i , π j ] = − q [ p i , A j ] − q [ A i , p j ] = − q [ p i , A j ] + q [ p j , A i ] 。[ p i , A j ] = − i ℏ ∂ i A j [p_i,A_j]=-i\hbar\partial_iA_j [ p i , A j ] = − i ℏ ∂ i A j なので [ π i , π j ] = i ℏ q ( ∂ i A j − ∂ j A i ) [\pi_i,\pi_j]=i\hbar q(\partial_iA_j-\partial_jA_i) [ π i , π j ] = i ℏ q ( ∂ i A j − ∂ j A i ) です。最後は B k = ε k i j ∂ i A j B_k=\varepsilon_{kij}\partial_iA_j B k = ε k ij ∂ i A j の定義です。
∎
運動量どうしが非可換になる理由は 磁場がある ことです。この単純な事実から、次の節の Landau 準位がすべて従います。
一様な磁場 B = B z ^ \boldsymbol{B}=B\hat{\boldsymbol{z}} B = B z ^ (B ≠ 0 B\ne0 B = 0 ) の中の粒子で、電場はないとします。z z z 方向は自由運動 ℏ 2 k z 2 / 2 m \hbar^2k_z^2/2m ℏ 2 k z 2 /2 m なので、x y xy x y 平面内の運動を考えます。
定理 3.1 (Landau 準位 ) ✓
x y xy x y 平面内の運動のハミルトニアン H ⊥ = π x 2 + π y 2 2 m H_\perp=\dfrac{\pi_x^2+\pi_y^2}{2m} H ⊥ = 2 m π x 2 + π y 2 のスペクトルは、サイクロトロン周波数 ω c = ∣ q B ∣ / m \omega_c=|qB|/m ω c = ∣ q B ∣/ m を使って
E n = ℏ ω c ( n + 1 2 ) , n = 0 , 1 , 2 , … E_n=\hbar\omega_c\left(n+\tfrac12\right),\qquad n=0,1,2,\ldots E n = ℏ ω c ( n + 2 1 ) , n = 0 , 1 , 2 , … である。各準位は(平面全体では)無限に縮退し、辺の長さ L x × L y L_x\times L_y L x × L y の長方形(y y y 方向に周期境界)では、縮退度は N ϕ = ∣ q B ∣ L x L y h N_\phi=\dfrac{|qB|L_xL_y}{h} N ϕ = h ∣ q B ∣ L x L y 、すなわち面積あたり ∣ q B ∣ / h |qB|/h ∣ q B ∣/ h である。
次の証明は 2 段に分かれます。前半は、π x , π y \pi_x,\pi_y π x , π y の交換関係だけから調和振動子と同じ形が出ることを示し、エネルギーの形を決めます。後半は、具体的なゲージを選んで実際に状態を作り、それが何個あるかを数えます。
証明(定理 3.1) 代数的な構造。 s = s g n ( q B ) s=\mathrm{sgn}(qB) s = sgn ( q B ) と置き、a = π x + i s π y 2 ℏ ∣ q B ∣ a=\dfrac{\pi_x+is\pi_y}{\sqrt{2\hbar|qB|}} a = 2ℏ∣ q B ∣ π x + i s π y を定めます。命題 2.3 より [ π x , π y ] = i ℏ q B [\pi_x,\pi_y]=i\hbar qB [ π x , π y ] = i ℏ q B なので
[ a , a † ] = 1 2 ℏ ∣ q B ∣ [ π x + i s π y , π x − i s π y ] = − 2 i s [ π x , π y ] 2 ℏ ∣ q B ∣ = 2 s q B 2 ∣ q B ∣ = 1 [a,a^\dagger]=\frac{1}{2\hbar|qB|}[\pi_x+is\pi_y,\ \pi_x-is\pi_y]=\frac{-2is\,[\pi_x,\pi_y]}{2\hbar|qB|}=\frac{2s\,qB}{2|qB|}=1 [ a , a † ] = 2ℏ∣ q B ∣ 1 [ π x + i s π y , π x − i s π y ] = 2ℏ∣ q B ∣ − 2 i s [ π x , π y ] = 2∣ q B ∣ 2 s q B = 1 です。また a † a = ( π x − i s π y ) ( π x + i s π y ) 2 ℏ ∣ q B ∣ a^\dagger a=\dfrac{(\pi_x-is\pi_y)(\pi_x+is\pi_y)}{2\hbar|qB|} a † a = 2ℏ∣ q B ∣ ( π x − i s π y ) ( π x + i s π y ) で、分子は π x 2 + π y 2 + i s [ π x , π y ] = π x 2 + π y 2 − s ℏ q B = π x 2 + π y 2 − ℏ ∣ q B ∣ \pi_x^2+\pi_y^2+is[\pi_x,\pi_y]=\pi_x^2+\pi_y^2-s\hbar qB=\pi_x^2+\pi_y^2-\hbar|qB| π x 2 + π y 2 + i s [ π x , π y ] = π x 2 + π y 2 − s ℏ q B = π x 2 + π y 2 − ℏ∣ q B ∣ なので、
H ⊥ = ℏ ω c ( a † a + 1 2 ) , ℏ ω c = ℏ ∣ q B ∣ m H_\perp=\hbar\omega_c\left(a^\dagger a+\tfrac12\right),\qquad \hbar\omega_c=\frac{\hbar|qB|}{m} H ⊥ = ℏ ω c ( a † a + 2 1 ) , ℏ ω c = m ℏ∣ q B ∣ です。[ a , a † ] = 1 [a,a^\dagger]=1 [ a , a † ] = 1 だけから a † a a^\dagger a a † a の固有値が 0 , 1 , 2 , … 0,1,2,\ldots 0 , 1 , 2 , … に限られることは、調和振動子と同じ代数の議論です(定理 5.3[1次元の簡単な系] )。これでスペクトルが E n E_n E n の形に限られることが、ゲージによらずに示せました。
存在と縮退。 各 n n n が実際に現れること、縮退度を調べるため、ランダウゲージ A = ( 0 , B x , 0 ) \boldsymbol{A}=(0,Bx,0) A = ( 0 , B x , 0 ) を取ります。H ⊥ = p x 2 2 m + ( p y − q B x ) 2 2 m H_\perp=\dfrac{p_x^2}{2m}+\dfrac{(p_y-qBx)^2}{2m} H ⊥ = 2 m p x 2 + 2 m ( p y − q B x ) 2 は y y y に依存しないので、p y p_y p y の固有値 ℏ k \hbar k ℏ k で ψ = e i k y φ ( x ) \psi=e^{iky}\varphi(x) ψ = e ik y φ ( x ) と分離できます。
[ p x 2 2 m + 1 2 m ω c 2 ( x − x k ) 2 ] φ = E φ , x k = ℏ k q B \left[\frac{p_x^2}{2m}+\frac12m\omega_c^2\left(x-x_k\right)^2\right]\varphi=E\varphi,\qquad x_k=\frac{\hbar k}{qB} [ 2 m p x 2 + 2 1 m ω c 2 ( x − x k ) 2 ] φ = E φ , x k = q B ℏ k で、これは中心 x k x_k x k の 1 次元調和振動子です。固有値は ℏ ω c ( n + 1 2 ) \hbar\omega_c(n+\frac12) ℏ ω c ( n + 2 1 ) で、すべての n ≥ 0 n\ge0 n ≥ 0 が現れます。固有関数は k k k で指定された中心 x k x_k x k にある波束で、k k k は連続(または周期境界 y ∼ y + L y y\sim y+L_y y ∼ y + L y では k = 2 π j / L y k=2\pi j/L_y k = 2 π j / L y 、j ∈ Z j\in\mathbb{Z} j ∈ Z )の値をとり、準位は k k k について縮退します。中心 x k x_k x k が箱の内部 0 < x k < L x 0<x_k<L_x 0 < x k < L x にある k k k の個数は、ℏ k = q B x k \hbar k=qBx_k ℏ k = q B x k なので k k k の幅 ∣ q B ∣ L x / ℏ |qB|L_x/\hbar ∣ q B ∣ L x /ℏ を k k k の間隔 2 π / L y 2\pi/L_y 2 π / L y で割って ∣ q B ∣ L x L y / 2 π ℏ = ∣ q B ∣ L x L y / h |qB|L_xL_y/2\pi\hbar=|qB|L_xL_y/h ∣ q B ∣ L x L y /2 π ℏ = ∣ q B ∣ L x L y / h です。
∎
縮退度は 磁束量子 Φ 0 = h / ∣ q ∣ \Phi_0=h/|q| Φ 0 = h /∣ q ∣ を使うと N ϕ = B L x L y / Φ 0 N_\phi=B L_xL_y/\Phi_0 N ϕ = B L x L y / Φ 0 、つまり面積を貫く磁束の中に含まれる磁束量子の数に等しく、1 つの Landau 準位が磁束量子 1 本あたり 1 状態を収容します。
例 3.2 (1 テスラの中の電子 )
B = 1 T B=1\ \mathrm{T} B = 1 T 、電子では ω c = e B / m e = 1.76 × 10 11 r a d / s \omega_c=eB/m_e=1.76\times10^{11}\ \mathrm{rad/s} ω c = e B / m e = 1.76 × 1 0 11 rad/s 、ℏ ω c = 0.116 m e V \hbar\omega_c=0.116\ \mathrm{meV} ℏ ω c = 0.116 meV (温度では約 1.3 K 1.3\ \mathrm{K} 1.3 K )、縮退度は e B / h = 2.42 × 10 14 m − 2 eB/h=2.42\times10^{14}\ \mathrm{m^{-2}} e B / h = 2.42 × 1 0 14 m − 2 です。波動関数の広がり(磁気長)は ℓ B = ℏ / ∣ q B ∣ = 25.7 n m \ell_B=\sqrt{\hbar/|qB|}=25.7\ \mathrm{nm} ℓ B = ℏ/∣ q B ∣ = 25.7 nm で、これは 1 個の電子が磁束量子 Φ 0 \Phi_0 Φ 0 の面積 2 π ℓ B 2 2\pi\ell_B^2 2 π ℓ B 2 を占有していることに対応します。半導体界面の 2 次元電子系を強磁場・低温に置くと Landau 準位が実際に観測され、整数・分数量子 Hall 効果の舞台になります。
磁場がゼロの領域に、ベクトルポテンシャルだけが残る場合を考えます。z z z 軸に沿った無限に長く細いソレノイドを考え、内部に磁束 Φ \Phi Φ を閉じ込めたとします。外側では B = 0 \boldsymbol{B}=0 B = 0 ですが、A = Φ 2 π r φ ^ \boldsymbol{A}=\dfrac{\Phi}{2\pi r}\hat{\boldsymbol{\varphi}} A = 2 π r Φ φ ^ は 0 0 0 になりません(例 3.5[電磁ポテンシャルとゲージ変換] )。
次の定理は、ソレノイドの外側では A \boldsymbol{A} A が単なるゲージ変換の形になっていて、波動関数が自由粒子の解に位相因子を掛けたものになる、ということを述べます。単連結な領域とは、穴がなく、領域内の輪を一点まで縮められる領域のことです。軸を含まない領域を選ぶのは、軸のまわりを一周する輪を含めないためです。
定理 4.1 (Aharonov–Bohm 位相 ) ✓
ソレノイドの外側の単連結な領域 D D D (軸を含まない)では、A = ∇ χ \boldsymbol{A}=\nabla\chi A = ∇ χ 、χ = Φ 2 π φ \chi=\dfrac{\Phi}{2\pi}\varphi χ = 2 π Φ φ (φ \varphi φ は方位角の D D D 上の 1 価な分枝)である。したがって D D D 上の定常解は ψ = e i q χ / ℏ ψ 0 \psi=e^{iq\chi/\hbar}\psi_0 ψ = e i q χ /ℏ ψ 0 と書け、ψ 0 \psi_0 ψ 0 は 磁場のない自由粒子の 方程式の解である。ソレノイドを時計回りと反時計回りに迂回する 2 つの経路の波が、1 周分の差で受ける位相差は
Δ θ = q Φ ℏ = 2 π Φ Φ 0 , Φ 0 = h q \Delta\theta=\frac{q\Phi}{\hbar}=2\pi\frac{\Phi}{\Phi_0},\qquad \Phi_0=\frac{h}{q} Δ θ = ℏ q Φ = 2 π Φ 0 Φ , Φ 0 = q h である。
証明(定理 4.1) ∇ φ = φ ^ / r \nabla\varphi=\hat{\boldsymbol{\varphi}}/r ∇ φ = φ ^ / r なので、∇ χ = Φ φ ^ / 2 π r = A \nabla\chi=\Phi\hat{\boldsymbol{\varphi}}/2\pi r=\boldsymbol{A} ∇ χ = Φ φ ^ /2 π r = A です。A = ∇ χ \boldsymbol{A}=\nabla\chi A = ∇ χ なら B = 0 \boldsymbol{B}=0 B = 0 で、定理 2.2 により A = ∇ χ \boldsymbol{A}=\nabla\chi A = ∇ χ のハミルトニアンの解は、χ \chi χ のゲージ変換で A = 0 \boldsymbol{A}=0 A = 0 の解 ψ 0 \psi_0 ψ 0 から ψ = e i q χ / ℏ ψ 0 \psi=e^{iq\chi/\hbar}\psi_0 ψ = e i q χ /ℏ ψ 0 と得られます。経路を D D D の中で 1 周分ずらすと φ \varphi φ が 2 π 2\pi 2 π だけ変わるので、2 つの経路上の因子 e i q χ / ℏ e^{iq\chi/\hbar} e i q χ /ℏ の比は e i q Φ / ℏ e^{iq\Phi/\hbar} e i q Φ/ℏ です。
∎
物理的に重要なのは、位相差 q Φ / ℏ q\Phi/\hbar q Φ/ℏ が 干渉縞のずれ として観測できる点です。実験では、電子線を 2 つに分けて、極細の磁石(トロイダル磁石)を挟んで干渉させ、磁束が電子に触れていないにもかかわらず縞が動くことが、1986 年に外村彰らの実験で確かめられました。
リング状の導体なら、この効果はエネルギー準位として厳密に現れます。
命題 4.2 (磁束を貫くリング ) ✓
半径 R R R のリング上を動く質量 m m m 、電荷 q q q の粒子で、リングの内側を磁束 Φ \Phi Φ が貫いている(リング上で B = 0 \boldsymbol{B}=0 B = 0 )とする。ハミルトニアン H = 1 2 m R 2 ( − i ℏ ∂ φ − q Φ 2 π ) 2 H=\dfrac{1}{2mR^2}\left(-i\hbar\partial_\varphi-\dfrac{q\Phi}{2\pi}\right)^2 H = 2 m R 2 1 ( − i ℏ ∂ φ − 2 π q Φ ) 2 (周期境界条件 ψ ( φ + 2 π ) = ψ ( φ ) \psi(\varphi+2\pi)=\psi(\varphi) ψ ( φ + 2 π ) = ψ ( φ ) )の固有値は
E n = ℏ 2 2 m R 2 ( n − Φ Φ 0 ) 2 , n ∈ Z E_n=\frac{\hbar^2}{2mR^2}\left(n-\frac{\Phi}{\Phi_0}\right)^2,\qquad n\in\mathbb{Z} E n = 2 m R 2 ℏ 2 ( n − Φ 0 Φ ) 2 , n ∈ Z で、固有関数は e i n φ / 2 π e^{in\varphi}/\sqrt{2\pi} e in φ / 2 π である。スペクトル全体が Φ → Φ + Φ 0 \Phi\to\Phi+\Phi_0 Φ → Φ + Φ 0 で不変である。
証明(命題 4.2) リング上で A φ = Φ / 2 π R A_\varphi=\Phi/2\pi R A φ = Φ/2 π R 、方位方向の運動量は p φ = − i ℏ R − 1 ∂ φ p_\varphi=-i\hbar R^{-1}\partial_\varphi p φ = − i ℏ R − 1 ∂ φ なので π φ = − i ℏ R − 1 ∂ φ − q Φ / 2 π R \pi_\varphi=-i\hbar R^{-1}\partial_\varphi-q\Phi/2\pi R π φ = − i ℏ R − 1 ∂ φ − q Φ/2 π R です。e i n φ e^{in\varphi} e in φ に作用させると π φ = ℏ R ( n − q Φ 2 π ℏ ) \pi_\varphi=\dfrac{\hbar}{R}(n-\dfrac{q\Phi}{2\pi\hbar}) π φ = R ℏ ( n − 2 π ℏ q Φ ) で、q Φ / 2 π ℏ = Φ / Φ 0 q\Phi/2\pi\hbar=\Phi/\Phi_0 q Φ/2 π ℏ = Φ/ Φ 0 なので固有値が決まります。周期条件から n n n は整数です。Φ → Φ + Φ 0 \Phi\to\Phi+\Phi_0 Φ → Φ + Φ 0 では n → n + 1 n\to n+1 n → n + 1 の付け替えで全体が元に戻ります。
∎
磁場がリング上でゼロでも、エネルギーが磁束で変わります。エネルギーを磁束で微分すると 持続電流 I = − ∂ E / ∂ Φ I=-\partial E/\partial\Phi I = − ∂ E / ∂ Φ が流れます。量子力学のポテンシャル A \boldsymbol{A} A は、閉じた経路に沿った積分 ∮ A ⋅ d l \oint\boldsymbol{A}\cdot d\boldsymbol{l} ∮ A ⋅ d l (ホロノミー)という形で、物理的な意味をもちます。ゲージ変換で変わらず、Stokes の定理で囲まれた磁束に等しいからです。
原子に弱い磁場や電場をかけたときの準位のずれは、摂動の行列要素 ⟨ α ′ j ′ m ′ ∣ V ∣ α j m ⟩ \langle\alpha' j'm'|V|\alpha jm\rangle ⟨ α ′ j ′ m ′ ∣ V ∣ α j m ⟩ で決まります。摂動は一様な場のとき、角運動量の言葉で自然に分類できる形(ベクトルなど)を持っています。この分類を使うと、m m m に関する行列要素の依存性が、角運動量の合成則だけで決まります。次の定義は、ベクトルが回転で成分を混ぜ合わせる決まり方を、より高い階数のものにも拡げたものです。回転を生成する J \boldsymbol{J} J との交換関係で、その混ざり方を述べています。
定義 5.1 (球テンソル演算子 ) ↗
全角運動量 J \boldsymbol{J} J に対し、2 k + 1 2k+1 2 k + 1 個の演算子の組 T q k T^k_q T q k (q = − k , … , k q=-k,\ldots,k q = − k , … , k )が 階数 k k k の球テンソル演算子 であるとは、
[ J z , T q k ] = ℏ q T q k , [ J ± , T q k ] = ℏ k ( k + 1 ) − q ( q ± 1 ) T q ± 1 k [J_z,T^k_q]=\hbar q\,T^k_q,\qquad [J_\pm,T^k_q]=\hbar\sqrt{k(k+1)-q(q\pm1)}\;T^k_{q\pm1} [ J z , T q k ] = ℏ q T q k , [ J ± , T q k ] = ℏ k ( k + 1 ) − q ( q ± 1 ) T q ± 1 k を満たすことをいう(J ± = J x ± i J y J_\pm=J_x\pm iJ_y J ± = J x ± i J y )。T k + 1 k = T − k − 1 k = 0 T^k_{k+1}=T^k_{-k-1}=0 T k + 1 k = T − k − 1 k = 0 とする。
たとえばベクトル演算子 [ J i , V j ] = i ℏ ε i j k V k [J_i,V_j]=i\hbar\varepsilon_{ijk}V_k [ J i , V j ] = i ℏ ε ij k V k の球面成分 V 0 = V z V_0=V_z V 0 = V z 、V ± 1 = ∓ ( V x ± i V y ) / 2 V_{\pm1}=\mp(V_x\pm iV_y)/\sqrt2 V ± 1 = ∓ ( V x ± i V y ) / 2 は階数 1 1 1 です(一例として [ J + , V z ] = − ℏ ( V x + i V y ) = ℏ 2 V + 1 [J_+,V_z]=-\hbar(V_x+iV_y)=\hbar\sqrt2\,V_{+1} [ J + , V z ] = − ℏ ( V x + i V y ) = ℏ 2 V + 1 が確かめられる)。位置 r \boldsymbol{r} r 、運動量、角運動量 J \boldsymbol{J} J 自身がそうです。定義の関係式は、定理 3.6[角運動量とスピン] で ∣ j , m ⟩ |j,m\rangle ∣ j , m ⟩ に対する J ± J_\pm J ± の作用を決めたものと同じ形です。つまり、T q k T^k_q T q k は ∣ k q ⟩ |kq\rangle ∣ k q ⟩ と同じように J \boldsymbol{J} J のもとで変換します 。
定理 5.2 (Wigner–Eckart の定理 ) ✓
{ ∣ α j m ⟩ } \{|\alpha jm\rangle\} { ∣ α j m ⟩} を、各 α , j \alpha,j α , j について標準位相(J ± J_\pm J ± の行列要素が 定理 3.6[角運動量とスピン] の形)で選んだ正規直交基底とする(α \alpha α は j , m j,m j , m 以外の量子数)。T q k T^k_q T q k を階数 k k k の球テンソル演算子とすると、
⟨ α ′ j ′ m ′ ∣ T q k ∣ α j m ⟩ = ⟨ k q ; j m ∣ j ′ m ′ ⟩ ⟨ α ′ j ′ ∥ T k ∥ α j ⟩ \langle\alpha'j'm'|T^k_q|\alpha jm\rangle=\langle kq;\,jm\,|\,j'm'\rangle\;\langle\alpha'j'\|T^k\|\alpha j\rangle ⟨ α ′ j ′ m ′ ∣ T q k ∣ α j m ⟩ = ⟨ k q ; j m ∣ j ′ m ′ ⟩ ⟨ α ′ j ′ ∥ T k ∥ α j ⟩ が成り立つ。⟨ k q ; j m ∣ j ′ m ′ ⟩ \langle kq;jm|j'm'\rangle ⟨ k q ; j m ∣ j ′ m ′ ⟩ は Clebsch–Gordan 係数、換算行列要素 ⟨ α ′ j ′ ∥ T k ∥ α j ⟩ \langle\alpha'j'\|T^k\|\alpha j\rangle ⟨ α ′ j ′ ∥ T k ∥ α j ⟩ は m , m ′ , q m,m',q m , m ′ , q によらない定数である。
証明の方針は、T q k T^k_q T q k を角運動量 k k k の状態 ∣ k q ⟩ |kq\rangle ∣ k q ⟩ と同じように扱うことです。∣ k q ⟩ |kq\rangle ∣ k q ⟩ と ∣ α j m ⟩ |\alpha jm\rangle ∣ α j m ⟩ の積に T q k T^k_q T q k を対応させる写像を考え、(i) それが角運動量の作用と両立すること、(ii) 合成則で積空間を分解すること、(iii) 行列要素を読み取ること、の 3 段で進みます。
証明(定理 5.2) V k = s p a n { ∣ k q ⟩ } V_k=\mathrm{span}\{|kq\rangle\} V k = span { ∣ k q ⟩} 、M = s p a n { ∣ α j m ⟩ } m \mathcal{M}=\mathrm{span}\{|\alpha jm\rangle\}_m M = span { ∣ α j m ⟩ } m とし、線形写像 Ψ : V k ⊗ M → H \Psi:V_k\otimes\mathcal{M}\to\mathcal{H} Ψ : V k ⊗ M → H を Ψ ( ∣ k q ⟩ ⊗ ∣ α j m ⟩ ) = T q k ∣ α j m ⟩ \Psi(|kq\rangle\otimes|\alpha jm\rangle)=T^k_q|\alpha jm\rangle Ψ ( ∣ k q ⟩ ⊗ ∣ α j m ⟩) = T q k ∣ α j m ⟩ で定めます。
(i) Ψ \Psi Ψ は角運動量の作用と可換。 V k ⊗ M V_k\otimes\mathcal{M} V k ⊗ M には J t o t = J ( k ) ⊗ 1 + 1 ⊗ J \boldsymbol{J}^{\rm tot}=\boldsymbol{J}^{(k)}\otimes1+1\otimes\boldsymbol{J} J tot = J ( k ) ⊗ 1 + 1 ⊗ J が作用します。J z J_z J z については、J z T q ∣ j m ⟩ = ( T q J z + ℏ q T q ) ∣ j m ⟩ = ℏ ( m + q ) T q ∣ j m ⟩ J_zT_q|jm\rangle=(T_qJ_z+\hbar qT_q)|jm\rangle=\hbar(m+q)T_q|jm\rangle J z T q ∣ j m ⟩ = ( T q J z + ℏ q T q ) ∣ j m ⟩ = ℏ ( m + q ) T q ∣ j m ⟩ 。J ± J_\pm J ± については、定義の関係式と J ± ∣ j m ⟩ = ℏ j ( j + 1 ) − m ( m ± 1 ) ∣ j m ± 1 ⟩ J_\pm|jm\rangle=\hbar\sqrt{j(j+1)-m(m\pm1)}|jm\pm1\rangle J ± ∣ j m ⟩ = ℏ j ( j + 1 ) − m ( m ± 1 ) ∣ j m ± 1 ⟩ から
J ± T q ∣ j m ⟩ = [ J ± , T q ] ∣ j m ⟩ + T q J ± ∣ j m ⟩ = ℏ k ( k + 1 ) − q ( q ± 1 ) T q ± 1 ∣ j m ⟩ + ℏ j ( j + 1 ) − m ( m ± 1 ) T q ∣ j m ± 1 ⟩ J_\pm T_q|jm\rangle=[J_\pm,T_q]|jm\rangle+T_qJ_\pm|jm\rangle=\hbar\sqrt{k(k+1)-q(q\pm1)}\,T_{q\pm1}|jm\rangle+\hbar\sqrt{j(j+1)-m(m\pm1)}\,T_q|jm\pm1\rangle J ± T q ∣ j m ⟩ = [ J ± , T q ] ∣ j m ⟩ + T q J ± ∣ j m ⟩ = ℏ k ( k + 1 ) − q ( q ± 1 ) T q ± 1 ∣ j m ⟩ + ℏ j ( j + 1 ) − m ( m ± 1 ) T q ∣ j m ± 1 ⟩ で、これは Ψ ( J ± t o t ∣ k q ⟩ ∣ j m ⟩ ) \Psi\bigl(J^{\rm tot}_\pm|kq\rangle|jm\rangle\bigr) Ψ ( J ± tot ∣ k q ⟩ ∣ j m ⟩ ) にほかなりません。
(ii) 既約分解。 定理 6.2[角運動量とスピン] により V k ⊗ M = ⨁ J = ∣ k − j ∣ k + j W J V_k\otimes\mathcal{M}=\bigoplus_{J=|k-j|}^{k+j}W_J V k ⊗ M = ⨁ J = ∣ k − j ∣ k + j W J で、W J W_J W J は標準位相の基底 ∣ J M ⟩ = ∑ q , m ⟨ k q ; j m ∣ J M ⟩ ∣ k q ⟩ ∣ j m ⟩ |JM\rangle=\sum_{q,m}\langle kq;jm|JM\rangle|kq\rangle|jm\rangle ∣ J M ⟩ = ∑ q , m ⟨ k q ; j m ∣ J M ⟩ ∣ k q ⟩ ∣ j m ⟩ を持ちます。ϕ M ( J ) = Ψ ∣ J M ⟩ \phi_M^{(J)}=\Psi|JM\rangle ϕ M ( J ) = Ψ∣ J M ⟩ と置くと、(i) により J z ϕ M = ℏ M ϕ M J_z\phi_M=\hbar M\phi_M J z ϕ M = ℏ M ϕ M 、J ± ϕ M = ℏ J ( J + 1 ) − M ( M ± 1 ) ϕ M ± 1 J_\pm\phi_M=\hbar\sqrt{J(J+1)-M(M\pm1)}\phi_{M\pm1} J ± ϕ M = ℏ J ( J + 1 ) − M ( M ± 1 ) ϕ M ± 1 です。したがって J 2 ϕ M = ℏ 2 J ( J + 1 ) ϕ M \boldsymbol{J}^2\phi_M=\hbar^2J(J+1)\phi_M J 2 ϕ M = ℏ 2 J ( J + 1 ) ϕ M (J 2 \boldsymbol{J}^2 J 2 は J z , J ± J_z,J_\pm J z , J ± で書ける)で、ϕ M ( J ) \phi_M^{(J)} ϕ M ( J ) は 0 であるか、標準位相の角運動量 J J J の多重項をなします。
(iii) 行列要素。 ⟨ α ′ j ′ m ′ ∣ ϕ M ( J ) ⟩ \langle\alpha'j'm'|\phi_M^{(J)}\rangle ⟨ α ′ j ′ m ′ ∣ ϕ M ( J ) ⟩ は、J 2 \boldsymbol{J}^2 J 2 の固有値が異なる(J ≠ j ′ J\ne j' J = j ′ )と 0 0 0 です。J = j ′ J=j' J = j ′ のときは、ϕ j ′ ( j ′ ) \phi_{j'}^{(j')} ϕ j ′ ( j ′ ) の ⟨ α ′ j ′ j ′ ∣ \langle\alpha'j'j'| ⟨ α ′ j ′ j ′ ∣ 成分を c α ′ c_{\alpha'} c α ′ と置くと、両辺に同じ降下演算子 J − J_- J − を作用させて同じ係数が出るので、すべての M M M について ⟨ α ′ j ′ M ∣ ϕ M ( j ′ ) ⟩ = c α ′ \langle\alpha'j'M|\phi_M^{(j')}\rangle=c_{\alpha'} ⟨ α ′ j ′ M ∣ ϕ M ( j ′ ) ⟩ = c α ′ 、M ≠ m ′ M\ne m' M = m ′ では 0 0 0 です。最後に ∣ k q ⟩ ∣ j m ⟩ = ∑ J , M ⟨ J M ∣ k q ; j m ⟩ ∣ J M ⟩ |kq\rangle|jm\rangle=\sum_{J,M}\langle JM|kq;jm\rangle|JM\rangle ∣ k q ⟩ ∣ j m ⟩ = ∑ J , M ⟨ J M ∣ k q ; j m ⟩ ∣ J M ⟩ なので
⟨ α ′ j ′ m ′ ∣ T q ∣ α j m ⟩ = ⟨ α ′ j ′ m ′ ∣ Ψ ( ∣ k q ⟩ ∣ j m ⟩ ) ⟩ = ⟨ k q ; j m ∣ j ′ m ′ ⟩ c α ′ \langle\alpha'j'm'|T_q|\alpha jm\rangle=\langle\alpha'j'm'|\Psi\bigl(|kq\rangle|jm\rangle\bigr)\rangle=\langle kq;jm|j'm'\rangle\,c_{\alpha'} ⟨ α ′ j ′ m ′ ∣ T q ∣ α j m ⟩ = ⟨ α ′ j ′ m ′ ∣Ψ ( ∣ k q ⟩ ∣ j m ⟩ ) ⟩ = ⟨ k q ; j m ∣ j ′ m ′ ⟩ c α ′ (Clebsch–Gordan 係数は実数)で、c α ′ c_{\alpha'} c α ′ が換算行列要素です。
∎
系 5.3 (選択則 ) ✓
⟨ α ′ j ′ m ′ ∣ T q k ∣ α j m ⟩ ≠ 0 \langle\alpha'j'm'|T^k_q|\alpha jm\rangle\ne0 ⟨ α ′ j ′ m ′ ∣ T q k ∣ α j m ⟩ = 0 ならば m ′ = m + q m'=m+q m ′ = m + q かつ ∣ j − k ∣ ≤ j ′ ≤ j + k |j-k|\le j'\le j+k ∣ j − k ∣ ≤ j ′ ≤ j + k である。
証明(系 5.3) Clebsch–Gordan 係数 ⟨ k q ; j m ∣ j ′ m ′ ⟩ \langle kq;jm|j'm'\rangle ⟨ k q ; j m ∣ j ′ m ′ ⟩ は m ′ = m + q m'=m+q m ′ = m + q のときだけ 0 0 0 でなく、j ′ j' j ′ は 定理 6.2[角運動量とスピン] の範囲 ∣ j − k ∣ ≤ j ′ ≤ j + k |j-k|\le j'\le j+k ∣ j − k ∣ ≤ j ′ ≤ j + k に入るときだけ現れます。
∎
これは、たとえば電気双極子遷移(位置 r \boldsymbol{r} r が階数 1 1 1 のテンソル)の Δ j = 0 , ± 1 \Delta j=0,\pm1 Δ j = 0 , ± 1 、Δ m = 0 , ± 1 \Delta m=0,\pm1 Δ m = 0 , ± 1 という選択則を一般的に示します。パリティの議論を加えると Δ l = ± 1 \Delta l=\pm1 Δ l = ± 1 も出ます。
命題 5.4 (射影定理 ) ✓
ベクトル演算子 V \boldsymbol{V} V と、固定された α , j \alpha,j α , j の多重項 M \mathcal{M} M への射影 P P P について、P V P = λ J P P\boldsymbol{V}P=\lambda\boldsymbol{J}P P V P = λ J P で、係数は λ = ⟨ J ⋅ V ⟩ ℏ 2 j ( j + 1 ) \lambda=\dfrac{\langle\boldsymbol{J}\cdot\boldsymbol{V}\rangle}{\hbar^2j(j+1)} λ = ℏ 2 j ( j + 1 ) ⟨ J ⋅ V ⟩ である。ここで ⟨ J ⋅ V ⟩ \langle\boldsymbol{J}\cdot\boldsymbol{V}\rangle ⟨ J ⋅ V ⟩ は M \mathcal{M} M の任意の状態での期待値(m m m によらない)。
証明(命題 5.4) V \boldsymbol{V} V も J \boldsymbol{J} J も階数 1 1 1 なので、定理 5.2 により M \mathcal{M} M 内の行列要素はどちらも ⟨ 1 q ; j m ∣ j m ′ ⟩ \langle 1q;jm|jm'\rangle ⟨ 1 q ; j m ∣ j m ′ ⟩ に比例し、比例定数は m m m に依りません。したがって P V q P = λ J q P PV_qP=\lambda J_qP P V q P = λ J q P (λ \lambda λ は定数)です。J \boldsymbol{J} J は P P P と交換するので P ( J ⋅ V ) P = ∑ i J i P V i P = λ ∑ i J i J i P = λ ℏ 2 j ( j + 1 ) P P(\boldsymbol{J}\cdot\boldsymbol{V})P=\sum_iJ_iPV_iP=\lambda\sum_iJ_iJ_iP=\lambda\hbar^2j(j+1)P P ( J ⋅ V ) P = ∑ i J i P V i P = λ ∑ i J i J i P = λ ℏ 2 j ( j + 1 ) P で、λ \lambda λ が決まります。
∎
一様磁場 B \boldsymbol{B} B 中の電子(q = − e q=-e q = − e )を、A = 1 2 B × r \boldsymbol{A}=\tfrac12\boldsymbol{B}\times\boldsymbol{r} A = 2 1 B × r (クーロンゲージ)で扱います。A ⋅ p = 1 2 ( B × r ) ⋅ p = 1 2 B ⋅ L \boldsymbol{A}\cdot\boldsymbol{p}=\tfrac12(\boldsymbol{B}\times\boldsymbol{r})\cdot\boldsymbol{p}=\tfrac12\boldsymbol{B}\cdot\boldsymbol{L} A ⋅ p = 2 1 ( B × r ) ⋅ p = 2 1 B ⋅ L なので、H = p 2 2 m + V − q 2 m B ⋅ L + q 2 A 2 2 m H=\dfrac{p^2}{2m}+V-\dfrac{q}{2m}\boldsymbol{B}\cdot\boldsymbol{L}+\dfrac{q^2A^2}{2m} H = 2 m p 2 + V − 2 m q B ⋅ L + 2 m q 2 A 2 です。第 4 項(反磁性項)は B 2 B^2 B 2 に比例し、実験室の磁場 B ≲ 10 T B\lesssim10\ \mathrm{T} B ≲ 10 T では第 3 項に比べて無視できます。スピンの磁気モーメント(注意 5.5[角運動量とスピン] 、g e ≃ 2 g_e\simeq2 g e ≃ 2 )を加えて、Zeeman ハミルトニアン
H Z = μ B ℏ B ⋅ ( L + 2 S ) , μ B = e ℏ 2 m H_Z=\frac{\mu_B}{\hbar}\,\boldsymbol{B}\cdot(\boldsymbol{L}+2\boldsymbol{S}),\qquad \mu_B=\frac{e\hbar}{2m} H Z = ℏ μ B B ⋅ ( L + 2 S ) , μ B = 2 m e ℏ
を得ます。スピンの g = 2 g=2 g = 2 は第 6 章の Dirac 方程式から出ます。
磁場が微細構造より弱いとき 、各 j j j の多重項の内部で摂動を考えれば十分です(縮退摂動論、定理 5.2[摂動論] )。
系 5.5 (Landé の g 因子 ) ✓
L + 2 S \boldsymbol{L}+2\boldsymbol{S} L + 2 S は多重項 ( l , s , j ) (l,s,j) ( l , s , j ) 内で g J J g_J\,\boldsymbol{J} g J J と同じ行列要素を持ち、
g J = 1 + j ( j + 1 ) − l ( l + 1 ) + s ( s + 1 ) 2 j ( j + 1 ) g_J=1+\frac{j(j+1)-l(l+1)+s(s+1)}{2j(j+1)} g J = 1 + 2 j ( j + 1 ) j ( j + 1 ) − l ( l + 1 ) + s ( s + 1 ) である。したがって B = B z ^ \boldsymbol{B}=B\hat{\boldsymbol{z}} B = B z ^ のとき準位は E n l j + g J μ B B m j E_{nlj}+g_J\mu_BB\,m_j E n l j + g J μ B B m j へ等間隔に分裂する。
証明(系 5.5) V = L + 2 S = J + S \boldsymbol{V}=\boldsymbol{L}+2\boldsymbol{S}=\boldsymbol{J}+\boldsymbol{S} V = L + 2 S = J + S に 命題 5.4 を使います。J ⋅ V = J 2 + J ⋅ S \boldsymbol{J}\cdot\boldsymbol{V}=\boldsymbol{J}^2+\boldsymbol{J}\cdot\boldsymbol{S} J ⋅ V = J 2 + J ⋅ S で、L = J − S \boldsymbol{L}=\boldsymbol{J}-\boldsymbol{S} L = J − S の 2 乗から J ⋅ S = 1 2 ( J 2 + S 2 − L 2 ) \boldsymbol{J}\cdot\boldsymbol{S}=\tfrac12(\boldsymbol{J}^2+\boldsymbol{S}^2-\boldsymbol{L}^2) J ⋅ S = 2 1 ( J 2 + S 2 − L 2 ) です。固有値を入れると ⟨ J ⋅ V ⟩ = ℏ 2 [ j ( j + 1 ) + 1 2 ( j ( j + 1 ) + s ( s + 1 ) − l ( l + 1 ) ) ] \langle\boldsymbol{J}\cdot\boldsymbol{V}\rangle=\hbar^2\left[j(j+1)+\tfrac12\bigl(j(j+1)+s(s+1)-l(l+1)\bigr)\right] ⟨ J ⋅ V ⟩ = ℏ 2 [ j ( j + 1 ) + 2 1 ( j ( j + 1 ) + s ( s + 1 ) − l ( l + 1 ) ) ] となり、ℏ 2 j ( j + 1 ) \hbar^2j(j+1) ℏ 2 j ( j + 1 ) で割れば g J g_J g J です。H Z = ( μ B / ℏ ) g J B J z H_Z=(\mu_B/\hbar)g_JBJ_z H Z = ( μ B /ℏ ) g J B J z の固有値は g J μ B B m j g_J\mu_BBm_j g J μ B B m j です。
∎
例 5.6 (ナトリウムの D 線 )
3 2 S 1 / 2 3^2S_{1/2} 3 2 S 1/2 (l = 0 , j = 1 2 l=0,j=\frac12 l = 0 , j = 2 1 )は g = 2 g=2 g = 2 、3 2 P 1 / 2 3^2P_{1/2} 3 2 P 1/2 は g = 2 3 g=\frac23 g = 3 2 、3 2 P 3 / 2 3^2P_{3/2} 3 2 P 3/2 は g = 4 3 g=\frac43 g = 3 4 です(公式にそれぞれ代入)。したがって D 1 D_1 D 1 線(P 1 / 2 → S 1 / 2 P_{1/2}\to S_{1/2} P 1/2 → S 1/2 )は Δ m = 0 , ± 1 \Delta m=0,\pm1 Δ m = 0 , ± 1 の組合せで 4 本、D 2 D_2 D 2 線(P 3 / 2 → S 1 / 2 P_{3/2}\to S_{1/2} P 3/2 → S 1/2 )は 6 本に分かれます。スピンのない原子では g J = 1 g_J=1 g J = 1 (s = 0 s=0 s = 0 、j = l j=l j = l )で、すべての線が同じ間隔 μ B B \mu_BB μ B B で分かれる 正常 Zeeman 効果 になります。実際の原子で観測される複雑な分裂は、g J g_J g J が準位ごとに異なるために起こる 異常 Zeeman 効果 です。
一様電場 F = F z ^ \boldsymbol{F}=F\hat{\boldsymbol{z}} F = F z ^ は摂動 H S = − q F z H_S=-qFz H S = − q F z を与えます。z z z は階数 1 1 1 、q = 0 q=0 q = 0 成分のテンソルなので、系 5.3 により ⟨ n ′ l ′ m ′ ∣ z ∣ n l m ⟩ ≠ 0 \langle n'l'm'|z|nlm\rangle\ne0 ⟨ n ′ l ′ m ′ ∣ z ∣ n l m ⟩ = 0 には m ′ = m m'=m m ′ = m と ∣ l ′ − l ∣ ≤ 1 |l'-l|\le1 ∣ l ′ − l ∣ ≤ 1 が必要で、さらにパリティ(z z z は奇)から l ′ = l ± 1 l'=l\pm1 l ′ = l ± 1 です。したがって、縮退のない状態の 1 次補正 ⟨ z ⟩ \langle z\rangle ⟨ z ⟩ は消え、効果は電場の 2 乗から始まります(水素の n ≥ 2 n\ge2 n ≥ 2 は l l l の違う状態が縮退しているので別で、例 5.5[摂動論] に 1 次の効果を示しました)。
例 5.7 (水素の基底状態の分極率 )
原子単位(ℏ = m = e = a 0 = 1 \hbar=m=e=a_0=1 ℏ = m = e = a 0 = 1 )で H 0 = − 1 2 ∇ 2 − 1 r H_0=-\tfrac12\nabla^2-\tfrac1r H 0 = − 2 1 ∇ 2 − r 1 、摂動 H S = F z H_S=Fz H S = F z 、ψ 0 = e − r / π \psi_0=e^{-r}/\sqrt\pi ψ 0 = e − r / π 、E 0 = − 1 2 E_0=-\tfrac12 E 0 = − 2 1 とします。定理 3.1[摂動論] の 2 次補正は無限和ですが、1 次の波動関数 ψ 1 \psi_1 ψ 1 が ( H 0 − E 0 ) ψ 1 = − F z ψ 0 (H_0-E_0)\psi_1=-Fz\,\psi_0 ( H 0 − E 0 ) ψ 1 = − F z ψ 0 を満たすことを使えば、E ( 2 ) = ⟨ ψ 0 ∣ F z ∣ ψ 1 ⟩ E^{(2)}=\langle\psi_0|Fz|\psi_1\rangle E ( 2 ) = ⟨ ψ 0 ∣ F z ∣ ψ 1 ⟩ です。ψ 1 = g ( r ) cos θ e − r / π \psi_1=g(r)\cos\theta\,e^{-r}/\sqrt\pi ψ 1 = g ( r ) cos θ e − r / π と置くと、( H 0 − E 0 ) (H_0-E_0) ( H 0 − E 0 ) の作用は
( − 1 2 g ′ ′ + g ′ − g ′ r + g r 2 ) cos θ e − r π \left(-\tfrac12g''+g'-\tfrac{g'}{r}+\tfrac{g}{r^2}\right)\cos\theta\;\frac{e^{-r}}{\sqrt\pi} ( − 2 1 g ′′ + g ′ − r g ′ + r 2 g ) cos θ π e − r に書けます(∇ 2 ( χ cos θ ) = ( χ ′ ′ + 2 χ ′ / r − 2 χ / r 2 ) cos θ \nabla^2(\chi\cos\theta)=(\chi''+2\chi'/r-2\chi/r^2)\cos\theta ∇ 2 ( χ cos θ ) = ( χ ′′ + 2 χ ′ / r − 2 χ / r 2 ) cos θ を使う)。これを − F r cos θ e − r / π -Fr\cos\theta\,e^{-r}/\sqrt\pi − F r cos θ e − r / π に等しいと置き、g = a r + b r 2 g=ar+br^2 g = a r + b r 2 を試すと左辺は a − 2 b + 2 b r a-2b+2br a − 2 b + 2 b r になり、b = − F / 2 b=-F/2 b = − F /2 、a = − F a=-F a = − F で一致します。つまり g = − F ( r + r 2 / 2 ) g=-F(r+r^2/2) g = − F ( r + r 2 /2 ) が厳密な解です(ψ 1 ⊥ ψ 0 \psi_1\perp\psi_0 ψ 1 ⊥ ψ 0 は cos θ \cos\theta cos θ の角度積分から)。2 次補正は
E ( 2 ) = F π ⋅ 4 π 3 ∫ 0 ∞ r 3 g ( r ) e − 2 r d r = − 4 F 2 3 ∫ 0 ∞ ( r 4 + r 5 2 ) e − 2 r d r = − 4 F 2 3 ⋅ 27 16 = − 9 4 F 2 E^{(2)}=\frac{F}{\pi}\cdot\frac{4\pi}{3}\int_0^\infty r^3g(r)e^{-2r}dr=-\frac{4F^2}{3}\int_0^\infty\left(r^4+\tfrac{r^5}{2}\right)e^{-2r}dr=-\frac{4F^2}{3}\cdot\frac{27}{16}=-\frac94F^2 E ( 2 ) = π F ⋅ 3 4 π ∫ 0 ∞ r 3 g ( r ) e − 2 r d r = − 3 4 F 2 ∫ 0 ∞ ( r 4 + 2 r 5 ) e − 2 r d r = − 3 4 F 2 ⋅ 16 27 = − 4 9 F 2 です。E ( 2 ) = − 1 2 α F 2 E^{(2)}=-\tfrac12\alpha F^2 E ( 2 ) = − 2 1 α F 2 と書くと分極率は α = 9 2 a 0 3 ( × 4 π ε 0 ) ≈ 0.667 A ˚ 3 \alpha=\tfrac92a_0^3\ (\times4\pi\varepsilon_0)\approx0.667\ \mathrm{\AA}^3 α = 2 9 a 0 3 ( × 4 π ε 0 ) ≈ 0.667 A ˚ 3 です。符号は 系 3.3[摂動論] の通りで、基底状態は電場で必ず下がります。電場は原子に双極子モーメント α F \alpha F α F を誘起し、その双極子と電場との相互作用が − 1 2 α F 2 -\tfrac12\alpha F^2 − 2 1 α F 2 です。
水素原子 のエネルギーは主量子数 n n n だけで決まりましたが、実際のスペクトルは細かく分裂しています。これは相対論的な補正のためで、α = e 2 / 4 π ε 0 ℏ c ≃ 1 / 137.036 \alpha=e^2/4\pi\varepsilon_0\hbar c\simeq1/137.036 α = e 2 /4 π ε 0 ℏ c ≃ 1/137.036 の 2 乗の大きさ(Bohr のエネルギー ∼ m c 2 α 2 \sim mc^2\alpha^2 ∼ m c 2 α 2 に対し m c 2 α 4 mc^2\alpha^4 m c 2 α 4 )です。補正は
H r e l = − p 4 8 m 3 c 2 , H S O = 1 2 m 2 c 2 1 r d V d r L ⋅ S , H D = ℏ 2 8 m 2 c 2 ∇ 2 V H_{\rm rel}=-\frac{p^4}{8m^3c^2},\qquad H_{\rm SO}=\frac{1}{2m^2c^2}\frac1r\frac{dV}{dr}\,\boldsymbol{L}\cdot\boldsymbol{S},\qquad H_{\rm D}=\frac{\hbar^2}{8m^2c^2}\nabla^2V H rel = − 8 m 3 c 2 p 4 , H SO = 2 m 2 c 2 1 r 1 d r d V L ⋅ S , H D = 8 m 2 c 2 ℏ 2 ∇ 2 V
の 3 つです。H r e l H_{\rm rel} H rel は p 2 c 2 + m 2 c 4 \sqrt{p^2c^2+m^2c^4} p 2 c 2 + m 2 c 4 を展開した第 3 項、H S O H_{\rm SO} H SO は電子が核の電場の中で動くと感じる磁場とスピンの磁気モーメントの相互作用(Thomas 歳差の因子 1 2 \tfrac12 2 1 を含む)、H D H_{\rm D} H D は Darwin 項です。この 3 項が Dirac 方程式の非相対論的展開から出ることは第 6 章で触れ、ここでは これを与えられた摂動として 、クーロンポテンシャル V = − ℏ c α / r V=-\hbar c\alpha/r V = − ℏ c α / r について 1 次の摂動論を実行します。この結果は第 6 章で Dirac 方程式の厳密解から再現されるので、結論の正しさは別に確かめられます。
必要な期待値をまず求めます。以下では原子単位(a 0 = 1 a_0=1 a 0 = 1 )で考え、∣ n l m ⟩ |nlm\rangle ∣ n l m ⟩ の動径部分を u ( r ) = r R ( r ) u(r)=rR(r) u ( r ) = r R ( r ) とします。
補題 6.1 (Kramers の漸化式 ) ✓
水素の束縛状態 ∣ n l m ⟩ |nlm\rangle ∣ n l m ⟩ について、s s s が s > − 2 l − 1 s>-2l-1 s > − 2 l − 1 を満たす範囲で
s + 1 n 2 ⟨ r s ⟩ − ( 2 s + 1 ) ⟨ r s − 1 ⟩ + s 4 [ ( 2 l + 1 ) 2 − s 2 ] ⟨ r s − 2 ⟩ = 0 \frac{s+1}{n^2}\langle r^s\rangle-(2s+1)\langle r^{s-1}\rangle+\frac s4\left[(2l+1)^2-s^2\right]\langle r^{s-2}\rangle=0 n 2 s + 1 ⟨ r s ⟩ − ( 2 s + 1 ) ⟨ r s − 1 ⟩ + 4 s [ ( 2 l + 1 ) 2 − s 2 ] ⟨ r s − 2 ⟩ = 0 が成り立つ。とくに ⟨ r − 1 ⟩ = 1 n 2 \langle r^{-1}\rangle=\dfrac1{n^2} ⟨ r − 1 ⟩ = n 2 1 、l ≥ 1 l\ge1 l ≥ 1 のとき ⟨ r − 3 ⟩ = ⟨ r − 2 ⟩ l ( l + 1 ) \langle r^{-3}\rangle=\dfrac{\langle r^{-2}\rangle}{l(l+1)} ⟨ r − 3 ⟩ = l ( l + 1 ) ⟨ r − 2 ⟩ である。
この補題は、r r r のべき乗の期待値どうしを結ぶ関係式です。証明では、動径方程式を、掛ける関数 f f f を自由に選んで 2 通りに積分し、両者を比べます。f f f の取り方で、好きなべき乗の間の関係が一度に得られます。
証明(補題 6.1) u ′ ′ = Q u u''=Qu u ′′ = Q u 、Q = l ( l + 1 ) r 2 − 2 r + 1 n 2 Q=\dfrac{l(l+1)}{r^2}-\dfrac2r+\dfrac1{n^2} Q = r 2 l ( l + 1 ) − r 2 + n 2 1 です。任意の関数 f ( r ) f(r) f ( r ) について、u ′ ′ = Q u u''=Qu u ′′ = Q u に f u ′ fu' f u ′ を掛けて積分し、右辺を f Q u u ′ = 1 2 f Q ( u 2 ) ′ fQuu'=\tfrac12fQ(u^2)' f Q u u ′ = 2 1 f Q ( u 2 ) ′ として部分積分すると ∫ f ′ u ′ 2 = ∫ ( f Q ) ′ u 2 \int f'u'^2=\int(fQ)'u^2 ∫ f ′ u ′2 = ∫ ( f Q ) ′ u 2 です。次に f ′ u f'u f ′ u を掛けて積分し、− ∫ ( f ′ u ) ′ u ′ = − ∫ f ′ ′ u u ′ − ∫ f ′ u ′ 2 = 1 2 ∫ f ′ ′ ′ u 2 − ∫ f ′ u ′ 2 -\int(f'u)'u'=-\int f''uu'-\int f'u'^2=\tfrac12\int f'''u^2-\int f'u'^2 − ∫ ( f ′ u ) ′ u ′ = − ∫ f ′′ u u ′ − ∫ f ′ u ′2 = 2 1 ∫ f ′′′ u 2 − ∫ f ′ u ′2 を使うと ∫ f ′ u ′ 2 = 1 2 ∫ f ′ ′ ′ u 2 − ∫ f ′ Q u 2 \int f'u'^2=\tfrac12\int f'''u^2-\int f'Qu^2 ∫ f ′ u ′2 = 2 1 ∫ f ′′′ u 2 − ∫ f ′ Q u 2 です(境界項は、f = r s + 1 f=r^{s+1} f = r s + 1 で s s s が上の範囲なら原点で u ∼ r l + 1 u\sim r^{l+1} u ∼ r l + 1 のため消える)。二式を等しいと置くと
∫ [ f Q ′ + 2 f ′ Q − 1 2 f ′ ′ ′ ] u 2 d r = 0 \int\left[fQ'+2f'Q-\tfrac12f'''\right]u^2\,dr=0 ∫ [ f Q ′ + 2 f ′ Q − 2 1 f ′′′ ] u 2 d r = 0 です。f = r s + 1 f=r^{s+1} f = r s + 1 を代入して r s − 2 , r s − 1 , r s r^{s-2},r^{s-1},r^s r s − 2 , r s − 1 , r s の係数を整理すると、s 2 [ ( 2 l + 1 ) 2 − s 2 ] \tfrac s2[(2l+1)^2-s^2] 2 s [( 2 l + 1 ) 2 − s 2 ] 、− 2 ( 2 s + 1 ) -2(2s+1) − 2 ( 2 s + 1 ) 、2 ( s + 1 ) / n 2 2(s+1)/n^2 2 ( s + 1 ) / n 2 になり、2 で割れば主張です。s = 0 s=0 s = 0 で ⟨ r − 1 ⟩ = n − 2 \langle r^{-1}\rangle=n^{-2} ⟨ r − 1 ⟩ = n − 2 、s = − 1 s=-1 s = − 1 で ⟨ r − 2 ⟩ − l ( l + 1 ) ⟨ r − 3 ⟩ = 0 \langle r^{-2}\rangle-l(l+1)\langle r^{-3}\rangle=0 ⟨ r − 2 ⟩ − l ( l + 1 ) ⟨ r − 3 ⟩ = 0 です。
∎
補題 6.2 (⟨r⁻²⟩ と s 状態の原点での値 ) ✓
⟨ r − 2 ⟩ = 1 n 3 ( l + 1 2 ) \langle r^{-2}\rangle=\dfrac{1}{n^3(l+\frac12)} ⟨ r − 2 ⟩ = n 3 ( l + 2 1 ) 1 である。また l = 0 l=0 l = 0 のとき ∣ ψ n 00 ( 0 ) ∣ 2 = 1 π n 3 |\psi_{n00}(0)|^2=\dfrac{1}{\pi n^3} ∣ ψ n 00 ( 0 ) ∣ 2 = π n 3 1 である。
証明(補題 6.2) 動径方程式 [ − 1 2 d 2 d r 2 + l ( l + 1 ) 2 r 2 − 1 r ] u = E u \left[-\tfrac12\tfrac{d^2}{dr^2}+\tfrac{l(l+1)}{2r^2}-\tfrac1r\right]u=Eu [ − 2 1 d r 2 d 2 + 2 r 2 l ( l + 1 ) − r 1 ] u = E u の固有値は E = − 1 2 ( n r + l + 1 ) 2 E=-\dfrac1{2(n_r+l+1)^2} E = − 2 ( n r + l + 1 ) 2 1 です(n r n_r n r は節の数)。この導出は、べき級数の終端(補題 5.1[水素原子] )が l l l を r l + 1 r^{l+1} r l + 1 の指数と漸化式にしか使わないため、l l l を連続な実数に拡げても成り立ちます。ヘルマン–ファインマン(命題 3.5[摂動論] )を λ = l ( l + 1 ) \lambda=l(l+1) λ = l ( l + 1 ) について使うと ∂ E / ∂ λ = 1 2 ⟨ r − 2 ⟩ \partial E/\partial\lambda=\tfrac12\langle r^{-2}\rangle ∂ E / ∂ λ = 2 1 ⟨ r − 2 ⟩ で、一方 ∂ E / ∂ l = n − 3 \partial E/\partial l=n^{-3} ∂ E / ∂ l = n − 3 、d λ / d l = 2 l + 1 d\lambda/dl=2l+1 d λ / d l = 2 l + 1 なので ⟨ r − 2 ⟩ = 2 / [ ( 2 l + 1 ) n 3 ] \langle r^{-2}\rangle=2/[(2l+1)n^3] ⟨ r − 2 ⟩ = 2/ [( 2 l + 1 ) n 3 ] です。
l = 0 l=0 l = 0 のとき、u ′ ′ = 2 ( V − E ) u u''=2(V-E)u u ′′ = 2 ( V − E ) u に 2 u ′ 2u' 2 u ′ を掛けて ( u ′ 2 ) ′ = 2 ( V − E ) ( u 2 ) ′ (u'^2)'=2(V-E)(u^2)' ( u ′2 ) ′ = 2 ( V − E ) ( u 2 ) ′ とし、r = 0 r=0 r = 0 から ∞ \infty ∞ まで積分します。右辺を部分積分すると境界項は r = 0 r=0 r = 0 (u ( 0 ) = 0 u(0)=0 u ( 0 ) = 0 と V u 2 → 0 Vu^2\to0 V u 2 → 0 )でも ∞ \infty ∞ でも消えて − u ′ ( 0 ) 2 = − 2 ∫ V ′ u 2 d r -u'(0)^2=-2\int V'u^2dr − u ′ ( 0 ) 2 = − 2 ∫ V ′ u 2 d r 、つまり u ′ ( 0 ) 2 = 2 ⟨ V ′ ⟩ u'(0)^2=2\langle V'\rangle u ′ ( 0 ) 2 = 2 ⟨ V ′ ⟩ です。ψ = R / 4 π = u / ( r 4 π ) \psi=R/\sqrt{4\pi}=u/(r\sqrt{4\pi}) ψ = R / 4 π = u / ( r 4 π ) なので ∣ ψ ( 0 ) ∣ 2 = u ′ ( 0 ) 2 / 4 π = ⟨ V ′ ⟩ / 2 π |\psi(0)|^2=u'(0)^2/4\pi=\langle V'\rangle/2\pi ∣ ψ ( 0 ) ∣ 2 = u ′ ( 0 ) 2 /4 π = ⟨ V ′ ⟩ /2 π です。V = − 1 / r V=-1/r V = − 1/ r なら V ′ = r − 2 V'=r^{-2} V ′ = r − 2 で、l = 0 l=0 l = 0 の ⟨ r − 2 ⟩ = 2 / n 3 \langle r^{-2}\rangle=2/n^3 ⟨ r − 2 ⟩ = 2/ n 3 から ∣ ψ ( 0 ) ∣ 2 = 1 / ( π n 3 ) |\psi(0)|^2=1/(\pi n^3) ∣ ψ ( 0 ) ∣ 2 = 1/ ( π n 3 ) を得ます。
∎
定理 6.3 (水素の微細構造 ) ✓
H r e l + H S O + H D H_{\rm rel}+H_{\rm SO}+H_{\rm D} H rel + H SO + H D の 1 次の摂動により、水素の準位 E n = − m c 2 α 2 2 n 2 E_n=-\dfrac{mc^2\alpha^2}{2n^2} E n = − 2 n 2 m c 2 α 2 は
E n j = E n − m c 2 α 4 2 n 4 ( n j + 1 2 − 3 4 ) E_{nj}=E_n-\frac{mc^2\alpha^4}{2n^4}\left(\frac{n}{j+\frac12}-\frac34\right) E nj = E n − 2 n 4 m c 2 α 4 ( j + 2 1 n − 4 3 ) となる。補正は l l l によらず n n n と全角運動量 j j j だけで決まる。
証明(定理 6.3) 縮退のある準位ですが、3 項はいずれも J 2 , J z , L 2 \boldsymbol{J}^2,J_z,\boldsymbol{L}^2 J 2 , J z , L 2 と交換する(球対称)ので、基底 ∣ n l j m j ⟩ |nljm_j\rangle ∣ n l j m j ⟩ では摂動行列が対角になります(( l , j , m j ) (l,j,m_j) ( l , j , m j ) が異なる 2 つの基底ベクトルは L 2 , J 2 , J z \boldsymbol{L}^2,\boldsymbol{J}^2,J_z L 2 , J 2 , J z のどれかの固有値が異なるので、補題 5.3[摂動論] から行列要素は 0 0 0 )。期待値を順に計算します。以下 ε 0 = m c 2 α 4 \varepsilon_0=mc^2\alpha^4 ε 0 = m c 2 α 4 と書き、長さは a 0 = ℏ / ( m c α ) a_0=\hbar/(mc\alpha) a 0 = ℏ/ ( m c α ) で測ります。
相対論補正。 p 4 = 4 m 2 ( H 0 − V ) 2 p^4=4m^2(H_0-V)^2 p 4 = 4 m 2 ( H 0 − V ) 2 で、H 0 H_0 H 0 は固有状態に作用して E n E_n E n になるので ⟨ H r e l ⟩ = − 1 2 m c 2 ⟨ ( E n − V ) 2 ⟩ = − 1 2 m c 2 [ E n 2 − 2 E n ⟨ V ⟩ + ⟨ V 2 ⟩ ] \langle H_{\rm rel}\rangle=-\dfrac{1}{2mc^2}\langle(E_n-V)^2\rangle=-\dfrac{1}{2mc^2}\left[E_n^2-2E_n\langle V\rangle+\langle V^2\rangle\right] ⟨ H rel ⟩ = − 2 m c 2 1 ⟨( E n − V ) 2 ⟩ = − 2 m c 2 1 [ E n 2 − 2 E n ⟨ V ⟩ + ⟨ V 2 ⟩ ] です。V = − ℏ c α / r V=-\hbar c\alpha/r V = − ℏ c α / r 、補題 6.1 の ⟨ r − 1 ⟩ \langle r^{-1}\rangle ⟨ r − 1 ⟩ から ⟨ V ⟩ = − m c 2 α 2 / n 2 \langle V\rangle=-mc^2\alpha^2/n^2 ⟨ V ⟩ = − m c 2 α 2 / n 2 、補題 6.2 から ⟨ V 2 ⟩ = ( m c 2 α 2 ) 2 / [ n 3 ( l + 1 2 ) ] \langle V^2\rangle=(mc^2\alpha^2)^2/[n^3(l+\frac12)] ⟨ V 2 ⟩ = ( m c 2 α 2 ) 2 / [ n 3 ( l + 2 1 )] です。代入すると
⟨ H r e l ⟩ = − ε 0 2 n 3 [ 1 l + 1 2 − 3 4 n ] \langle H_{\rm rel}\rangle=-\frac{\varepsilon_0}{2n^3}\left[\frac{1}{l+\frac12}-\frac{3}{4n}\right] ⟨ H rel ⟩ = − 2 n 3 ε 0 [ l + 2 1 1 − 4 n 3 ] です。
スピン軌道。 1 r d V d r = ℏ c α / r 3 \frac1r\frac{dV}{dr}=\hbar c\alpha/r^3 r 1 d r d V = ℏ c α / r 3 で、⟨ L ⋅ S ⟩ = ℏ 2 2 [ j ( j + 1 ) − l ( l + 1 ) − 3 4 ] \langle\boldsymbol{L}\cdot\boldsymbol{S}\rangle=\frac{\hbar^2}{2}\left[j(j+1)-l(l+1)-\frac34\right] ⟨ L ⋅ S ⟩ = 2 ℏ 2 [ j ( j + 1 ) − l ( l + 1 ) − 4 3 ] (系 6.5[角運動量とスピン] )です。l ≥ 1 l\ge1 l ≥ 1 では 補題 6.1 と 補題 6.2 から ⟨ r − 3 ⟩ = 1 / [ n 3 l ( l + 1 2 ) ( l + 1 ) ] \langle r^{-3}\rangle=1/[n^3l(l+\frac12)(l+1)] ⟨ r − 3 ⟩ = 1/ [ n 3 l ( l + 2 1 ) ( l + 1 )] なので
⟨ H S O ⟩ = ε 0 4 n 3 j ( j + 1 ) − l ( l + 1 ) − 3 4 l ( l + 1 2 ) ( l + 1 ) \langle H_{\rm SO}\rangle=\frac{\varepsilon_0}{4n^3}\,\frac{j(j+1)-l(l+1)-\frac34}{l(l+\frac12)(l+1)} ⟨ H SO ⟩ = 4 n 3 ε 0 l ( l + 2 1 ) ( l + 1 ) j ( j + 1 ) − l ( l + 1 ) − 4 3 です。l = 0 l=0 l = 0 では L = 0 \boldsymbol{L}=0 L = 0 で ⟨ H S O ⟩ = 0 \langle H_{\rm SO}\rangle=0 ⟨ H SO ⟩ = 0 とします。
Darwin 項。 ∇ 2 V = 4 π ℏ c α δ 3 ( r ) \nabla^2V=4\pi\hbar c\alpha\,\delta^3(\boldsymbol{r}) ∇ 2 V = 4 π ℏ c α δ 3 ( r ) なので、⟨ H D ⟩ = π ℏ 3 α 2 m 2 c ∣ ψ ( 0 ) ∣ 2 \langle H_{\rm D}\rangle=\dfrac{\pi\hbar^3\alpha}{2m^2c}|\psi(0)|^2 ⟨ H D ⟩ = 2 m 2 c π ℏ 3 α ∣ ψ ( 0 ) ∣ 2 で、これは l = 0 l=0 l = 0 だけに寄与します。補題 6.2 の ∣ ψ ( 0 ) ∣ 2 = 1 / ( π n 3 a 0 3 ) |\psi(0)|^2=1/(\pi n^3a_0^3) ∣ ψ ( 0 ) ∣ 2 = 1/ ( π n 3 a 0 3 ) を入れて ⟨ H D ⟩ = ε 0 2 n 3 δ l 0 \langle H_{\rm D}\rangle=\dfrac{\varepsilon_0}{2n^3}\delta_{l0} ⟨ H D ⟩ = 2 n 3 ε 0 δ l 0 です。
足し合わせ。 (a) l = 0 l=0 l = 0 、j = 1 2 j=\frac12 j = 2 1 :⟨ H r e l ⟩ + ⟨ H D ⟩ = − ε 0 2 n 3 [ 2 − 3 4 n − 1 ] = − ε 0 2 n 3 [ 1 − 3 4 n ] \langle H_{\rm rel}\rangle+\langle H_{\rm D}\rangle=-\dfrac{\varepsilon_0}{2n^3}\left[2-\dfrac3{4n}-1\right]=-\dfrac{\varepsilon_0}{2n^3}\left[1-\dfrac{3}{4n}\right] ⟨ H rel ⟩ + ⟨ H D ⟩ = − 2 n 3 ε 0 [ 2 − 4 n 3 − 1 ] = − 2 n 3 ε 0 [ 1 − 4 n 3 ] 。(b) l ≥ 1 l\ge1 l ≥ 1 、j = l + 1 2 j=l+\frac12 j = l + 2 1 :スピン軌道の分子は l l l なので ⟨ H S O ⟩ = ε 0 4 n 3 1 ( l + 1 2 ) ( l + 1 ) \langle H_{\rm SO}\rangle=\dfrac{\varepsilon_0}{4n^3}\dfrac{1}{(l+\frac12)(l+1)} ⟨ H SO ⟩ = 4 n 3 ε 0 ( l + 2 1 ) ( l + 1 ) 1 で、1 ( l + 1 2 ) ( l + 1 ) ⋅ 1 2 − 1 l + 1 2 = − 1 l + 1 \dfrac1{(l+\frac12)(l+1)}\cdot\dfrac12-\dfrac1{l+\frac12}=-\dfrac{1}{l+1} ( l + 2 1 ) ( l + 1 ) 1 ⋅ 2 1 − l + 2 1 1 = − l + 1 1 より ∑ = − ε 0 2 n 3 [ 1 l + 1 − 3 4 n ] \sum=-\dfrac{\varepsilon_0}{2n^3}\left[\dfrac1{l+1}-\dfrac3{4n}\right] ∑ = − 2 n 3 ε 0 [ l + 1 1 − 4 n 3 ] 。(c) j = l − 1 2 j=l-\frac12 j = l − 2 1 :分子は − ( l + 1 ) -(l+1) − ( l + 1 ) で ⟨ H S O ⟩ = − ε 0 4 n 3 1 l ( l + 1 2 ) \langle H_{\rm SO}\rangle=-\dfrac{\varepsilon_0}{4n^3}\dfrac1{l(l+\frac12)} ⟨ H SO ⟩ = − 4 n 3 ε 0 l ( l + 2 1 ) 1 、− 1 2 l ( l + 1 2 ) − 1 l + 1 2 = − 1 l -\dfrac{1}{2l(l+\frac12)}-\dfrac1{l+\frac12}=-\dfrac1l − 2 l ( l + 2 1 ) 1 − l + 2 1 1 = − l 1 より ∑ = − ε 0 2 n 3 [ 1 l − 3 4 n ] \sum=-\dfrac{\varepsilon_0}{2n^3}\left[\dfrac1l-\dfrac3{4n}\right] ∑ = − 2 n 3 ε 0 [ l 1 − 4 n 3 ] 。いずれの場合も、l = 0 l=0 l = 0 の (a) で j + 1 2 = 1 j+\frac12=1 j + 2 1 = 1 、(b) で j + 1 2 = l + 1 j+\frac12=l+1 j + 2 1 = l + 1 、(c) で j + 1 2 = l j+\frac12=l j + 2 1 = l なので、結果は − ε 0 2 n 3 [ 1 j + 1 2 − 3 4 n ] -\dfrac{\varepsilon_0}{2n^3}\left[\dfrac{1}{j+\frac12}-\dfrac{3}{4n}\right] − 2 n 3 ε 0 [ j + 2 1 1 − 4 n 3 ] に統一されます。これは主張の式と一致します。
∎
例 6.4 (n = 2 の微細構造 )
n = 2 n=2 n = 2 では、2 S 1 / 2 2S_{1/2} 2 S 1/2 と 2 P 1 / 2 2P_{1/2} 2 P 1/2 (どちらも j = 1 2 j=\frac12 j = 2 1 )が縮退したまま、2 P 3 / 2 2P_{3/2} 2 P 3/2 (j = 3 2 j=\frac32 j = 2 3 )が上に分かれます。分裂幅は
E 2 , 3 / 2 − E 2 , 1 / 2 = m c 2 α 4 32 = 1.449 × 10 − 3 e V 32 = 4.53 × 10 − 5 e V ( 10.9 G H z ) E_{2,3/2}-E_{2,1/2}=\frac{mc^2\alpha^4}{32}=\frac{1.449\times10^{-3}\ \mathrm{eV}}{32}=4.53\times10^{-5}\ \mathrm{eV}\ \ (10.9\ \mathrm{GHz}) E 2 , 3/2 − E 2 , 1/2 = 32 m c 2 α 4 = 32 1.449 × 1 0 − 3 eV = 4.53 × 1 0 − 5 eV ( 10.9 GHz ) で、基本の準位間隔(10.2 e V 10.2\ \mathrm{eV} 10.2 eV )の 4 × 10 − 6 4\times10^{-6} 4 × 1 0 − 6 倍にすぎません。これが「微細」の意味です。2 S 1 / 2 2S_{1/2} 2 S 1/2 と 2 P 1 / 2 2P_{1/2} 2 P 1/2 の縮退は Dirac 方程式でも厳密に残りますが、量子電磁力学の補正(Lamb シフト、約 1.06 G H z 1.06\ \mathrm{GHz} 1.06 GHz )で解けます。
演習 7.1 易
面積 1 c m 2 1\ \mathrm{cm}^2 1 cm 2 の 2 次元電子系に B = 5 T B=5\ \mathrm{T} B = 5 T をかける。(1) Landau 準位 1 本に入れる軌道状態の数(スピンを除く)を求めよ。(2) 下から 2 本の Landau 準位をスピンも含めて満たすために必要な面密度(ν = 4 \nu=4 ν = 4 )はいくらか。
解答 (1) 定理 3.1 より N ϕ = e B A / h = 2.418 × 10 14 m − 2 T − 1 × 5 T × 10 − 4 m 2 = 1.21 × 10 11 N_\phi=eBA/h=2.418\times10^{14}\ \mathrm{m^{-2}\,T^{-1}}\times5\ \mathrm{T}\times10^{-4}\ \mathrm{m^2}=1.21\times10^{11} N ϕ = e B A / h = 2.418 × 1 0 14 m − 2 T − 1 × 5 T × 1 0 − 4 m 2 = 1.21 × 1 0 11 。(2) スピンを含めると各準位は 2 N ϕ 2N_\phi 2 N ϕ 個の電子を収容でき、2 本で 4 N ϕ = 4.8 × 10 11 4N_\phi=4.8\times10^{11} 4 N ϕ = 4.8 × 1 0 11 個、面密度では 4.8 × 10 11 c m − 2 4.8\times10^{11}\ \mathrm{cm^{-2}} 4.8 × 1 0 11 c m − 2 です。
演習 7.2 標準
(1) 3 P 2 ^3P_2 3 P 2 (S = 1 , L = 1 , J = 2 S=1,L=1,J=2 S = 1 , L = 1 , J = 2 )の g J g_J g J を求めよ。(2) 2 P 3 / 2 ^2P_{3/2} 2 P 3/2 の準位に磁場 B B B をかけたとき、m j = − 3 2 , … , 3 2 m_j=-\frac32,\ldots,\frac32 m j = − 2 3 , … , 2 3 の 4 つの準位のエネルギーを μ B B \mu_BB μ B B を単位にして書け。
解答 (1) 系 5.5 に j = 2 , l = 1 , s = 1 j=2,l=1,s=1 j = 2 , l = 1 , s = 1 を代入して g J = 1 + 6 − 2 + 2 12 = 3 2 g_J=1+\dfrac{6-2+2}{12}=\dfrac32 g J = 1 + 12 6 − 2 + 2 = 2 3 です。(2) g J = 4 3 g_J=\dfrac43 g J = 3 4 なので、エネルギーは E 0 + 4 3 m j μ B B E_0+\dfrac43m_j\mu_BB E 0 + 3 4 m j μ B B 、すなわち − 2 , − 2 3 , + 2 3 , + 2 -2,-\frac23,+\frac23,+2 − 2 , − 3 2 , + 3 2 , + 2 の順に μ B B \mu_BB μ B B をかけた値です。隣り合う準位の間隔は等間隔の 4 3 μ B B \frac43\mu_BB 3 4 μ B B です。
演習 7.3 標準
命題 4.2 のリングで、0 ≤ Φ ≤ Φ 0 0\le\Phi\le\Phi_0 0 ≤ Φ ≤ Φ 0 の範囲で基底状態を与える n n n を求め、持続電流 I = − ∂ E / ∂ Φ I=-\partial E/\partial\Phi I = − ∂ E / ∂ Φ を Φ \Phi Φ の関数として書け。Φ = Φ 0 / 2 \Phi=\Phi_0/2 Φ = Φ 0 /2 で何が起きるか。
解答 E n E_n E n が最小になる n n n は ∣ n − Φ / Φ 0 ∣ |n-\Phi/\Phi_0| ∣ n − Φ/ Φ 0 ∣ が最小の整数で、0 ≤ Φ < Φ 0 / 2 0\le\Phi<\Phi_0/2 0 ≤ Φ < Φ 0 /2 で n = 0 n=0 n = 0 、Φ 0 / 2 < Φ ≤ Φ 0 \Phi_0/2<\Phi\le\Phi_0 Φ 0 /2 < Φ ≤ Φ 0 で n = 1 n=1 n = 1 です。I = − ∂ E n / ∂ Φ = ℏ 2 m R 2 Φ 0 ( n − Φ Φ 0 ) I=-\partial E_n/\partial\Phi=\dfrac{\hbar^2}{mR^2\Phi_0}\left(n-\dfrac{\Phi}{\Phi_0}\right) I = − ∂ E n / ∂ Φ = m R 2 Φ 0 ℏ 2 ( n − Φ 0 Φ ) なので、0 ≤ Φ < Φ 0 / 2 0\le\Phi<\Phi_0/2 0 ≤ Φ < Φ 0 /2 では I = − ℏ 2 m R 2 Φ 0 Φ Φ 0 I=-\dfrac{\hbar^2}{mR^2\Phi_0}\dfrac{\Phi}{\Phi_0} I = − m R 2 Φ 0 ℏ 2 Φ 0 Φ で磁束に比例して逆向きに流れ、Φ = Φ 0 / 2 \Phi=\Phi_0/2 Φ = Φ 0 /2 で n = 0 n=0 n = 0 と n = 1 n=1 n = 1 の準位が交差して電流の向きが反転します(のこぎり波状の電流)。電流は Φ \Phi Φ に対して周期 Φ 0 \Phi_0 Φ 0 です。
演習 7.4 標準
定理 6.3 を使って、n = 2 n=2 n = 2 の 2 P 1 / 2 2P_{1/2} 2 P 1/2 と 2 P 3 / 2 2P_{3/2} 2 P 3/2 のエネルギー補正を m c 2 α 4 mc^2\alpha^4 m c 2 α 4 を単位にそれぞれ求め、差が m c 2 α 4 / 32 mc^2\alpha^4/32 m c 2 α 4 /32 になることを確かめよ。
解答 n = 2 n=2 n = 2 で E ( 1 ) = − m c 2 α 4 32 ( 2 j + 1 2 − 3 4 ) E^{(1)}=-\dfrac{mc^2\alpha^4}{32}\left(\dfrac{2}{j+\frac12}-\dfrac34\right) E ( 1 ) = − 32 m c 2 α 4 ( j + 2 1 2 − 4 3 ) です。j = 1 2 j=\frac12 j = 2 1 :2 1 − 3 4 = 5 4 \dfrac2{1}-\dfrac34=\dfrac54 1 2 − 4 3 = 4 5 で E ( 1 ) = − 5 128 m c 2 α 4 E^{(1)}=-\dfrac{5}{128}mc^2\alpha^4 E ( 1 ) = − 128 5 m c 2 α 4 。j = 3 2 j=\frac32 j = 2 3 :2 2 − 3 4 = 1 4 \dfrac22-\dfrac34=\dfrac14 2 2 − 4 3 = 4 1 で E ( 1 ) = − 1 128 m c 2 α 4 E^{(1)}=-\dfrac{1}{128}mc^2\alpha^4 E ( 1 ) = − 128 1 m c 2 α 4 。差は 4 128 = 1 32 \dfrac{4}{128}=\dfrac1{32} 128 4 = 32 1 (単位 m c 2 α 4 mc^2\alpha^4 m c 2 α 4 )です。
J. J. Sakurai, J. Napolitano, Modern Quantum Mechanics , 3rd ed., Cambridge University Press, 2020 — テンソル演算子と Wigner–Eckart の定理、微細構造、Aharonov–Bohm 効果、Landau 準位。
D. J. Griffiths, D. F. Schroeter, Introduction to Quantum Mechanics , 3rd ed., Cambridge University Press, 2018 — 水素の微細構造、Zeeman 効果、Stark 効果、Aharonov–Bohm 効果の章。
L. D. Landau, E. M. Lifshitz, Quantum Mechanics: Non-relativistic Theory , 3rd ed., Pergamon, 1977 — 磁場中の運動(Landau 準位)、原子の準位(微細構造)。
Y. Aharonov, D. Bohm, “Significance of Electromagnetic Potentials in the Quantum Theory”, Phys. Rev. 115 (1959), 485.
A. Tonomura et al., “Evidence for Aharonov-Bohm Effect with Magnetic Field Completely Shielded from Electron Wave”, Phys. Rev. Lett. 56 (1986), 792.
H. A. Kramers, Quantentheorie des Elektrons und der Strahlung , 1938. Kramers の漸化式の出典。
執筆: Claude Code(Sonnet 5.5)、Claude Code(Opus 5.5)・監査 3 回(Claude Code)・最終 2026-10-12