Feynman QED Twelfth Lecture

\(\)

\(\textit{Twelfth Lecture}\)

第11講の結果である式 (11-6) と式 (11-7) を用いると, 低エネルギー近似 \((W-V)\ll 2mc^{2}\) が行える.それはオーダー \((v/c)^{4}\) までの項を残す近似である.近似式 \((1+x)^{-1}\approx 1-x\ \, (|x|< 1)\) を用いると,
\(
\def\ppdiff#1#2{\frac{\partial #1}{\partial #2}}
\def\pdiff#1{\frac{\partial}{\partial #1}}
\def\Bppdiff#1#2{\frac{\partial^{2} #1}{\partial #2^{2}}}
\def\Bpdiff#1{\frac{\partial^{2}}{\partial #1^{2}}}
\def\mb#1{\mathbf{#1}}
\def\mr#1{\mathrm{#1}}
\def\reverse#1{\frac{1}{#1}}
\def\ds#1{\mbox{${\displaystyle\strut #1}$}}
\def\half{\frac{1}{2}}
\newcommand\Braket[1]{\left\langle #1 \right\rangle}
\)
\begin{align*}
\frac{c^{2}}{2mc^{2}+W-V}&=\frac{c^{2}}{2mc^{2}\left\{\ds{1+\frac{W-V}{2mc^{2}}}\right\}}
=\frac{1}{2m}\left(1+\frac{W-V}{2mc^{2}}\right)^{-1}\\
&\approx \frac{1}{2m}\left\{1-\frac{W-V}{2mc^{2}}\right\}
=\frac{1}{2m}-\frac{W-V}{(2mc)^{2}}
\tag{12-1}
\end{align*}
すると, 式 (11-7) は次となる:
\begin{align*}
(W-V)\psi_a &=\mb{\sigma}\cdot\mb{\pi}\,\frac{c^{2}}{2mc^{2}+W-V}\,\mb{\sigma}\cdot\mb{\pi}\,\psi_a\\
&= (\mb{\sigma}\cdot\mb{\pi})\left\{\frac{1}{2m}-\frac{W-V}{4m^{2}c^{2}}\right\}
(\mb{\sigma}\cdot\mb{\pi})\psi_a\\
&=\frac{1}{2m}(\mb{\sigma}\cdot\mb{\pi})^{2}\psi_a
-\frac{1}{4m^{2}c^{2}}(\mb{\sigma}\cdot\mb{\pi})(W-V)(\mb{\sigma}\cdot\mb{\pi})\psi_a
\tag{12-2}
\end{align*}
他方, 規格化の必要条件 \(\displaystyle \int \bigl(\psi_a^{2}+\psi_b^{2}\bigr)\,dVol =1\) は, 式 (11-6) を

\begin{equation*}
\psi_b=\frac{c}{2mc^{2}+\cancel{W-V}}\,\mb{\sigma}\cdot\mb{\pi}\,\psi_a
\approx \frac{1}{2mc}\,\mb{\sigma}\cdot\mb{\pi}\,\psi_a =\frac{\mb{\sigma}\cdot\mb{\pi}}{2mc}\psi_a
\end{equation*}

と近似することで次となる:
\begin{align*}
\int \bigl(\psi_a^{2}+\psi_b^{2}\bigr)\,dVol
&\approx \int \left(\psi^{*}_a\psi_a
+ \psi_a^{*}\frac{\mb{\sigma}\cdot\mb{\pi}}{2mc}\cdot\frac{\mb{\sigma}\cdot\mb{\pi}}{2mc}\psi_a\right)\,dVol\\
&=\int \psi^{*}_a\left(1+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{4m^{2}c^{2}}\right)\psi_a\,dVol =1
\tag{12-3}
\end{align*}
次のような置換え
\begin{equation*}
\chi =\Omega\,\psi_a,\quad
\Omega\equiv \left(1+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\right),\quad\rightarrow\quad
\Omega^{-1} = \left(1-\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\right)
\tag{12-4}
\end{equation*}

を用いると \(v^{2}/c^{2}\) のオーダーで
\begin{align*}
\chi^{*}\chi&=\psi_a^{*}\left(1+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\right)
\left(1+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\right)\psi_a
=\psi_a^{*}\left\{1+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{4m^{2}c^{2}}
+\cancel{\frac{(\mb{\sigma}\cdot\mb{\pi})^{4}}{64m^{4}c^{4}}}\right\}\psi_a \\
&\approx \psi_a^{*}\left\{1+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{4m^{2}c^{2}}\right\}\psi_a
\end{align*}
となるので, 規格化積分の式 (12-3) は単純化することが出来て, 次のような表現となる:
\begin{equation*}
\int \chi^{*}\chi\,dVol =1
\tag{12-3′}
\end{equation*}

この置換えは, 式 (12-2) のより分かり易い解釈も可能とする.式 (12-4) から \(\psi_a=\Omega^{-1}\chi\) である.「\(\mb{\pi}\) と \(V\) とは可換ではない!」ことに注意して, まず次式を近似する:
\begin{align*}
&\Omega(W-V)\,\chi=\Omega(W-V)\Omega\,\psi_a =\left(1+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\right)(W-V)
\left(1+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\right)\psi_a \\
&=(W-V)\psi_a +(W-V)\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\psi_a
+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}(W-V)\psi_a
+\cancel{\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}(W-V)\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\psi_a}\\
&\approx (W-V)\psi_a
+(W-V)\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\psi_a
+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}(W-V)\psi_a
\end{align*}
右辺第1項に式 (12-2) を用いると,
\begin{align*}
\Omega(W-V)\,\chi
&=\frac{1}{2m}(\mb{\sigma}\cdot\mb{\pi})^{2}\psi_a
-\frac{1}{4m^{2}c^{2}}(\mb{\sigma}\cdot\mb{\pi})(W-V)(\mb{\sigma}\cdot\mb{\pi})\psi_a
+(W-V)\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\psi_a
+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}(W-V)\psi_a\\
&=\frac{1}{2m}(\mb{\sigma}\cdot\mb{\pi})^{2}\psi_a
+\frac{1}{8m^{2}c^{2}}\bigl[(\mb{\sigma}\cdot\mb{\pi})^{2}(W-V)-2(\mb{\sigma}\cdot\mb{\pi})(W-V)
(\mb{\sigma}\cdot\mb{\pi})+(W-V)(\mb{\sigma}\cdot\mb{\pi})^{2}\bigr]\psi_a
\end{align*}
この両辺に \(\Omega^{-1}\) を左から作用させる.右辺第1項では \(\Omega^{-1}\) と \((\mb{\sigma}\cdot\mb{\pi})^{2}\) とは可換であることに注意し, 右辺第2項に作用する \(\Omega^{-1}\) はオーダーの関係から \(\Omega^{-1}\approx 1\) と近似してしまうと,
\begin{align*}
&\Omega^{-1}\Omega (W-V)\,\chi=(W-V)\,\chi \\
&\quad =\frac{1}{2m}(\mb{\sigma}\cdot\mb{\pi})^{2}\Omega^{-1}\psi_a
+\Omega^{-1}\frac{1}{8m^{2}c^{2}}\bigl[(\mb{\sigma}\cdot\mb{\pi})^{2}(W-V)-2(\mb{\sigma}\cdot\mb{\pi})(W-V)
(\mb{\sigma}\cdot\mb{\pi})+(W-V)(\mb{\sigma}\cdot\mb{\pi})^{2}\bigr]\psi_a\\
&\quad =\frac{1}{2m}(\mb{\sigma}\cdot\mb{\pi})^{2}\Omega^{-1}\Omega^{-1}\chi
+\Omega^{-1}\frac{1}{8m^{2}c^{2}}\bigl[(\mb{\sigma}\cdot\mb{\pi})^{2}(W-V)-2(\mb{\sigma}\cdot\mb{\pi})(W-V)
(\mb{\sigma}\cdot\mb{\pi})+(W-V)(\mb{\sigma}\cdot\mb{\pi})^{2}\bigr]\Omega^{-1}\chi\\
&\quad \approx \frac{1}{2m}(\mb{\sigma}\cdot\mb{\pi})^{2}\Omega^{-2}\chi
+\frac{1}{8m^{2}c^{2}}\bigl[(\mb{\sigma}\cdot\mb{\pi})^{2}(W-V)-2(\mb{\sigma}\cdot\mb{\pi})(W-V)
(\mb{\sigma}\cdot\mb{\pi})+(W-V)(\mb{\sigma}\cdot\mb{\pi})^{2}\bigr]\chi
\end{align*}
さらに, この第1項目は
\begin{align*}
\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{2m}\Omega^{-2}\chi
&=\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{2m}\left(1+\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\right)^{-2}\chi
\approx \frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{2m}
\left(1-2\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{8m^{2}c^{2}}\right)\chi \\
&=\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{2m}\chi-\frac{(\mb{\sigma}\cdot\mb{\pi})^{4}}{8m^{3}c^{2}}\chi
\end{align*}
よって, 次の結果となる:
\begin{align*}
(W-V)\,\chi&=\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{2m}\chi-\frac{(\mb{\sigma}\cdot\mb{\pi})^{4}}{8m^{3}c^{2}}\chi\\
&\quad +\frac{1}{8m^{2}c^{2}}\left[(\mb{\sigma}\cdot\mb{\pi})^{2}(W-V)-2(\mb{\sigma}\cdot\mb{\pi})(W-V)(\mb{\sigma}\cdot\mb{\pi})+(W-V)(\mb{\sigma}\cdot\mb{\pi})^{2}\right]\chi
\tag{12-5}
\end{align*}
式 (12-5) をより解釈し易い形に変換するには, 演算子代数の技法が使えるであろう.特に, 次の公式を思い出すべきである:
\begin{align*}
A^{2}B -2ABA +BA^{2}&=A(AB-BA)-(AB-BA)A\\
&=A[A,B]-[A,B]A=[A,[A,B]]
\tag{F1}
\end{align*}
そして \(\displaystyle\mb{\pi}=\left(\mb{p}-\frac{e}{c}\mb{A}\right)\) であり, また \(W\) は定数であることを考慮すると,
\begin{align*}
&\bigl[\mb{\sigma}\cdot\mb{A},(W-V)\bigr]=0,\quad \mb{p}W=-i\hbar\nabla W =0,\\
&\mb{\sigma}\cdot\mb{p}\,(W-V)=\mb{\sigma}\cdot\bigl[\mb{p}\,(W-V)\bigr]+(W-V)\,\mb{\sigma}\cdot\mb{p}
=\mb{\sigma}\cdot\mb{p}\,(-V)+(W-V)\,\mb{\sigma}\cdot\mb{p}
\end{align*}
従って,
\begin{align*}
(\mb{\sigma}\cdot\mb{\pi})(W-V)-(W-V)(\mb{\sigma}\cdot\mb{\pi})
&=[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)]
=\left[\mb{\sigma}\cdot\mb{p}-\frac{e}{c}\mb{\sigma}\cdot\mb{A},\,(W-V)\right]\\
&=\left[\mb{\sigma}\cdot\mb{p},\,(W-V)\right]
=\mb{\sigma}\cdot\mb{p}\,(W-V)-(W-V)\,\mb{\sigma}\cdot\mb{p}\\
&=\mb{\sigma}\cdot\mb{p}\,(-V)+(W-V)\,\mb{\sigma}\cdot\mb{p}-(W-V)\,\mb{\sigma}\cdot\mb{p}\\
&=-\mb{\sigma}\cdot\mb{p}V=-\mb{\sigma}\cdot(-i\hbar\nabla V)=i\hbar\mb{\sigma}\cdot\nabla V
\end{align*}
ここで \(V=e\phi\) とすると, 電磁気学の公式 \(\mb{E}=-\nabla \phi\) から \(e\mb{E}=-\nabla V\) である.従って, 次式が得られる:
\begin{equation*}
(\mb{\sigma}\cdot\mb{\pi})(W-V)-(W-V)(\mb{\sigma}\cdot\mb{\pi})
=i\hbar\mb{\sigma}\cdot\nabla V =-ie\hbar\mb{\sigma}\cdot\mb{E}
\tag{1}
\end{equation*}

前述の公式 (F1) に於いて \(\mb{\sigma}\cdot\mb{\pi}=A\) そして \((W-V)=B\) とすると,
\begin{align*}
\Bigl[(\mb{\sigma}\cdot\mb{\pi}),\,\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr]\Bigr]
&=(\mb{\sigma}\cdot\mb{\pi})\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr]
-\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr](\mb{\sigma}\cdot\mb{\pi})\\
&=(\mb{\sigma}\cdot\mb{\pi})\bigl[-ie\hbar\mb{\sigma}\cdot\mb{E}\bigr]
-\bigl[-ie\hbar\mb{\sigma}\cdot\mb{E}\bigr](\mb{\sigma}\cdot\mb{\pi})\\
&=-ie\hbar\Bigl\{(\mb{\sigma}\cdot\mb{\pi})(\mb{\sigma}\cdot\mb{E})-(\mb{\sigma}\cdot\mb{E})
(\mb{\sigma}\cdot\mb{\pi})\Bigr\}
\tag{2}
\end{align*}
ここで \(\mb{\sigma}\) についての重要な公式
\begin{equation*}
(\mb{\sigma}\cdot\mb{a})(\mb{\sigma}\cdot\mb{b})=\mb{a}\cdot\mb{b}+i\mb{\sigma}(\mb{a}\times\mb{b})
\tag{F2}
\end{equation*}

を上式に用いる.演算子の後ろには \(\mb{x}\) の関数 \(\circ\) が存在することに注意すると,
\begin{align*}
(\mb{\sigma}\cdot\mb{\pi})(\mb{\sigma}\cdot\mb{E})
&=\mb{\pi}\cdot\mb{E}+i\mb{\sigma}\cdot(\mb{\pi}\times\mb{E})\\
&=\left(\mb{p}-\frac{e}{c}\mb{A}\right)\cdot\mb{E}
+i\mb{\sigma}\cdot\left\{\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times\mb{E}\right\}\\
&=\mb{p}\cdot\mb{E}\circ -\frac{e}{c}\mb{A}\cdot\mb{E}+i\mb{\sigma}\cdot(\mb{p}\times\mb{E}\circ)
-i\frac{e}{c}\mb{\sigma}\cdot(\mb{A}\times\mb{E})\\
&=(\mb{p}\cdot\mb{E})\circ + \mb{E}\cdot(\mb{p}\circ) -\frac{e}{c}\mb{A}\cdot\mb{E}
+i\mb{\sigma}\cdot\bigl\{(\mb{p}\times\mb{E})\circ -\mb{E}\times(\mb{p}\circ)\bigr\}
-i\frac{e}{c}\mb{\sigma}\cdot(\mb{A}\times\mb{E})\\
&=(\mb{p}\cdot\mb{E})\circ + \mb{E}\cdot(\mb{p}\circ) -\frac{e}{c}\mb{E}\cdot\mb{A}
+i\mb{\sigma}\cdot(\mb{p}\times\mb{E})\circ-i\mb{\sigma}\cdot\bigl\{\mb{E}\times(\mb{p}\circ)
-\frac{e}{c}\mb{E}\times\mb{A}\bigr\}\\
&=(\mb{p}\cdot\mb{E})\circ +\mb{E}\cdot(\mb{\pi}\circ)+i\mb{\sigma}\cdot(\mb{p}\times\mb{E})\circ
-i\mb{\sigma}\bigl\{\mb{E}\times(\mb{\pi}\circ)\bigr\},\\
(\mb{\sigma}\cdot\mb{E})(\mb{\sigma}\cdot\mb{\pi})
&=\mb{E}\cdot\mb{\pi}+i\mb{\sigma}\cdot(\mb{E}\times\mb{\pi})\\
&=\mb{E}\cdot(\mb{\pi}\circ)+i\mb{\sigma}\cdot\bigl\{\mb{E}\times(\mb{\pi}\circ)\bigr\},\\
\therefore\quad
(\mb{\sigma}\cdot\mb{\pi})(\mb{\sigma}\cdot\mb{E})&-(\mb{\sigma}\cdot\mb{E})(\mb{\sigma}\cdot\mb{\pi})
=(\mb{p}\cdot\mb{E})-2i\mb{\sigma}\cdot(\mb{E}\times\mb{\pi})+i\mb{\sigma}\cdot(\mb{p}\times\mb{E})
\tag{3}
\end{align*}
ここで, 磁場はほぼ定常であると仮定する.すると,
\begin{equation*}
\mathrm{rot}\,\mb{E}=-\frac{1}{c}\ppdiff{\mb{H}}{t}\approx 0\quad\rightarrow\quad
\nabla\times\mb{E}=0
\end{equation*}

従って,
\begin{equation*}
\mb{p}\times\mb{E}=-i\hbar\nabla\times\mb{E}=0
\end{equation*}

よって, 前式 (3) は,
\begin{align*}
(\mb{\sigma}\cdot\mb{\pi})(\mb{\sigma}\cdot\mb{E})-(\mb{\sigma}\cdot\mb{E})(\mb{\sigma}\cdot\mb{\pi})
&=\mb{p}\cdot\mb{E}-2i\mb{\sigma}\cdot(\mb{E}\times\mb{\pi})\\
&=-i\hbar\nabla\cdot\mb{E} -2i\mb{\sigma}\cdot(\mb{E}\times\mb{\pi})
\tag{4}
\end{align*}
そして式 (2) は,
\begin{align*}
\left[(\mb{\sigma}\cdot\mb{\pi}),\,\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr]\right]
&=-ie\hbar\Bigl\{(\mb{\sigma}\cdot\mb{\pi})(\mb{\sigma}\cdot\mb{E})-(\mb{\sigma}\cdot\mb{E})
(\mb{\sigma}\cdot\mb{\pi})\Bigr\}\\
&=-ie\hbar\Bigl\{-i\hbar\nabla\cdot\mb{E} -2i\mb{\sigma}\cdot(\mb{E}\times\mb{\pi})\Bigr\}\\
&=-e\hbar^{2}\nabla\cdot\mb{E}-2e\hbar\mb{\sigma}\cdot(\mb{E}\times\mb{\pi})
\end{align*}
以上の結果から, 式 (12-5) は,
\begin{align*}
(W-V)\,\chi&=\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{2m}\chi
-\frac{(\mb{\sigma}\cdot\mb{\pi})^{4}}{8m^{3}c^{2}}\chi\\
&\quad +\frac{1}{8m^{2}c^{2}}\left[(\mb{\sigma}\cdot\mb{\pi}),\,
\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr]\right]\chi\\
&=\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{2m}\chi
-\frac{(\mb{\sigma}\cdot\mb{\pi})^{4}}{8m^{3}c^{2}}\chi
-\frac{e}{8m^{2}c^{2}}\Bigl\{\hbar^{2}\nabla\cdot\mb{E}
+2\hbar\mb{\sigma}\cdot(\mb{E}\times\mb{\pi})\Bigr\}\chi
\tag{5}
\end{align*}
更に, 公式 (F2) を用いる.\(\displaystyle \mb{\pi}=\mb{p}-\frac{e}{c}\mb{A}\) なので, この場合の \(\mb{\pi}\times\mb{\pi}\ne0\) に注意する:
\begin{align*}
\mb{\pi}\times\mb{\pi}&=\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)
=\mb{p}\times\mb{p}-\frac{e}{c}\mb{p}\times\mb{A}
-\frac{e}{c}\mb{A}\times\mb{p}+\left(\frac{e}{c}\right)^{2}\mb{A}\times{A}\\
&=-\frac{e}{c}\bigl\{\mb{A}\times\mb{p}+\mb{p}\times\mb{A}\bigr\}
=-\frac{e}{c}\bigl\{\mb{A}\times(-i\hbar\nabla)+(-i\hbar\nabla)\times\mb{A}\bigr\}\\
&=i\frac{e\hbar}{c}\bigl\{\mb{A}\times\nabla+\nabla\times\mb{A}\bigr\}
=i\frac{e\hbar}{c}\,\nabla\times\mb{A}=i\frac{e\hbar}{c}\,\mathrm{rot}\,\mb{A}
=i\frac{e\hbar}{c}\,\mb{H}
\end{align*}
従って,
\begin{align*}
(\mb{\sigma}\cdot\mb{\pi})^{2}&=(\mb{\sigma}\cdot\mb{\pi})(\mb{\sigma}\cdot\mb{\pi})
=\mb{\pi}\cdot\mb{\pi}+i\mb{\sigma}\cdot(\mb{\pi}\times\mb{\pi})\\
&=\left(\mb{p}-\frac{e}{c}\mb{A}\right)\cdot\left(\mb{p}-\frac{e}{c}\mb{A}\right)
+i\mb{\sigma}\cdot\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times \left(\mb{p}-\frac{e}{c}\mb{A}\right)\\
&=\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}-\frac{e\hbar}{c}\mb{\sigma}\cdot\mb{H}
\end{align*}
また \((\mb{\sigma}\cdot\mb{\pi})^{4}\) についてはオーダーの関係から \(\mb{\pi}\sim \mb{p}\) として次とする:
\begin{equation*}
\bigl(\mb{\sigma}\cdot\mb{\pi}\bigr)^{4}\approx \bigl(\mb{\sigma}\cdot\mb{p}\bigr)^{4}=\mb{p}^{4}
\end{equation*}

よって, 式 (5) すなわち式 (12-5) は最終的に次の様に展開される:
\begin{align*}
W\chi &=V\chi +\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{2m}\chi
-\frac{(\mb{\sigma}\cdot\mb{\pi})^{4}}{8m^{3}c^{2}}\chi
-\frac{e}{8m^{2}c^{2}}\Bigl\{\hbar^{2}\nabla\cdot\mb{E}
+2\hbar\mb{\sigma}\cdot(\mb{E}\times\mb{\pi})\Bigr\}\chi \\
&=\Biggl[\underbrace{V}_{(1)} +\underbrace{\frac{1}{2m}\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}}_{(2)}
-\underbrace{\frac{e\hbar}{2mc}\mb{\sigma}\cdot\mb{H}}_{(3)}-\underbrace{\frac{\mb{p}^{4}}{8m^{3}c^{2}}}_{(4)}
-\underbrace{\frac{e\hbar^{2}}{8m^{2}c^{2}}\nabla\cdot\mb{E}}_{(5)}-\underbrace{\frac{e\hbar}{4m^{2}c^{2}}\mb{\sigma}\cdot
\mb{E}\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)}_{(6)}\,\Biggr]\chi
\tag{12-6}
\end{align*}
この形式では, 波動方程式は式 (12-6) の各項を別々に考えることで解釈できる.
項(1)は, 以前に登場した通常の「スカラーポテンシャルエネルギー」を与える.
項(2)は,「運動エネルギー」と解釈できる.
項(3)は,「パウリ-スピンの効果」であり, まさにパウリ方程式に出現した項である.
項(4)は,「運動エネルギーの相対論的補正」である.補正は次式に由来する:
\begin{equation*}
E=\sqrt{m^{2}c^{4}+c^{2}\mb{p}^{2}}=mc^{2}\sqrt{1+\frac{\mb{p}^{2}}{m^{2}c^{2}}}
=mc^{2}+\frac{\mb{p}^{2}}{2m}-\frac{\mb{p}^{4}}{8m^{3}c^{2}}+\dotsb
\end{equation*}

この表式の最後の項が, 項 (4) に相当している.
項(5)と項(6)は,「スピン-軌道結合」を表現している.この解釈を理解するために, 式 (6) に於ける \(\mb{\sigma}\cdot(\mb{p}\times\mb{E})\) 部分を考えてみる.逆2乗の場に於いては \(\displaystyle \mb{E}=\frac{Ze^{2}}{r^{2}}\frac{\mb{r}}{r}=\frac{Ze^{2}}{r^{3}}\mb{r}\) であるから, これは \(\mb{\sigma}\cdot(\mb{p}\times\mb{r})/r^{3}\) に比例する.因子 \(\mb{r}\times\mb{p}\) は「角運動量」\(\mb{L}\) と解釈出来るので, 従って スピン-軌道結合 \((\mb{\sigma}\cdot\mb{L})/r^{3}=(2/\hbar)\mb{S}\cdot\mb{L}/r^{3}\sim\mb{S}\cdot\mb{L}\) を得る.この項は, 電子が S 状態 \((L=0)\) に在るときには何の効果もない.他方, 項 (5) は \(\nabla\cdot\mb{E}=-\nabla^{2}\phi =4\pi Z|e|\delta(0)\) となる.これは (波動関数が \(\mb{r}=0\) でゼロでない) S 状態にのみ影響する.従って, 項 (5) と項 (6) は一緒になって, スピン-軌道結合に対する連続関数に帰着する(result in).電子の磁気モーメント \(\displaystyle \mu=\frac{e\hbar}{2mc}\) は, 項 (3) の係数として出現し, そして項 (5) の係数 \(\displaystyle \frac{\hbar}{4mc}\left(\frac{e\hbar}{2mc}\right)\) 及び項 (6) の係数 \(\displaystyle \frac{1}{2mc}\left(\frac{e\hbar}{2mc}\right)\) として再び出現する.

項(6)を解釈するのに, 古典的議論をすることが可能である.速度 \(\mb{v}\) で電場中を運動する電荷は, 有効磁場 \(\displaystyle \mb{H}=\frac{\mb{v}}{c}\times\mb{E}=\frac{1}{mc}\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times\mb{E}\) を感じる [1][ブログ註] オッペンハイマー:「電気力学」§ 11 によれば,「一様な速度 \(\mb{v}\) で動いている点電荷 \(e\) の作る電場 \(\mb{E}\) と磁場 … Continue reading.そしてこの磁場中で, 項 (6) はちょうど, エネルギー \(\displaystyle E_{(6)}=-\frac{1}{2}\mb{\mu}\cdot\mb{H}\) になっている:
\begin{align*}
E_{(6)}&=-\frac{e\hbar}{4m^{2}c^{2}}\mb{\sigma}\cdot\mb{E}\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)
=-\left(\frac{e}{2mc}\right)\cdot\frac{\hbar}{2}\mb{\sigma}\cdot\frac{1}{mc}\left(\mb{p}-\frac{e}{c}\mb{A}\right)
\times\mb{E}\\
&=-\frac{e}{2mc} \mb{S}\cdot\mb{H}=-\frac{1}{2}\,\mb{\mu}\cdot\mb{H}
\end{align*}
しかしながら, このように因子 \(2\) だけ余分に得られてしまう.ディラック方程式が開発される以前から, L.H.Thomas はこの単純な古典的議論が不完全であることを示し, 正しい項 (6) を与えていた.パウリが中性子と陽子を記述するために導入した「異常磁気モーメント」については, 状況は異なる(下の問題 (3) を見よ).パウリの変形方程式では, 項 (5) と項 (6) を掛け合わせたとき, 異常磁気モーメントは因子 \(2\) を持って現れる.


【 問題 】 (1) 式 (12-6) を水素原子に適用して, そのエネルギー準位を1次のオーダーで修正せよ.その結果を実際の実験結果と比較せよ [2] Schiff, “Quantum Mechanics,” McGraw-Hill, New York, 1949, pp. 323ff..座標原点に於ける波動関数の違いに注意せよ.しかし, その違いはスペースが限られているために重要な意味を持たない.水素様原子に対するディラック方程式の原点付近の厳密解は,

\begin{equation*}
r^{\sqrt{1-(Z\alpha)^{2}}-1}\approx r^{-\frac{1}{2}\left(\frac{Z}{137}\right)^{2}} \sim r^{-1/40,000}
\end{equation*}

に比例するが, 他方, シュレディンガー方程式では \(r\to 0\) に連れて波動関数 \(\psi\) は一定となる.

(2) \(\mb{A}\) と \(\phi\) が時間に依存すると仮定しよう.\(W=i\hbar\partial/\partial t\) として, 同じ近似のオーダーまでこの講義の手順をたどってみよ.

(3) パウリの変形方程式は, 中性子と陽子に適用することが出来る.それは, ディラック方程式に異常磁気モーメントの項を付加することで得ることが出来る [3][ブログ註] (A.P.アーヤ:基礎現代物理学 § 11.2 から抜粋) 原子スペクトルでは, 原子の磁気モーメントは「ボーア磁子」\(\mu_B\) … Continue reading.従って,

\begin{equation*}
\gamma^{\mu}\left(i\hbar\partial_{\mu} – \frac{e}{c}A_{\mu}\right)\psi\,-\, \frac{i\mu}{c}\left(\frac{1}{2}\gamma^{\mu}\gamma^{\nu}F_{\mu\nu}\right)\,\psi =mc\,\psi
\end{equation*}

この式は \(c\beta\) を掛け合わせることで, より馴染みのある「ハミルトニアン表現」に書くことが出来る:
\begin{align*}
i\hbar\pdiff{t}\psi &= H_{\mathrm{Dirac}}\,\psi
+i\mu\beta(\mb{\alpha}\cdot\mb{E}+i\,\mb{\Sigma}\cdot\mb{H})\,\psi,\quad \left(\mu=\frac{e\hbar}{2mc}\kappa\right), \\
\mathrm{where}\quad &H_{\mathrm{Dirac}}= c\mb{\alpha}\cdot\left(\mb{p}
-\frac{e}{c}\mb{A}\right) +\beta mc^{2} + e\phi,\quad
\gamma^{\mu}\gamma^{\nu}F_{\mu\nu}
=2\mb{\alpha}\cdot\mb{E}+2i\,\mb{\Sigma}\cdot\mb{H}
\end{align*}
(ただし, 朝永:「スピンはめぐる」や L.Foldy, Phys. Rev. \(\mathbf{87}\), (1952) を参照して, 原書の式中の付加項表現は修正したので注意すべし.また \(\gamma^{\mu}\gamma^{\nu}F_{\mu\nu}\) の展開は, 第9講の最後の問題を見よ).式 (12-6) を導出するのと同じ近似をすると, 今度は陽子に対しては
\begin{align*}
&\left[V+\frac{1}{2M}\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}
+\left(\mu+\frac{e}{2M}\right)\mb{\sigma}\cdot\mb{H}
+\frac{1}{8M^{3}}(\mb{p}\cdot\mb{p})^{2}\right.\\
&\left.\qquad\qquad +\frac{1}{4M^{2}}\left(2\mu+\frac{e}{2M}\right)
\left\{\nabla\cdot\mb{E}+2\mb{\sigma}\cdot\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times\mb{E}\right\}
\right]\,\psi
\tag{12-7}
\end{align*}
のような項が生成し, そして中性子の場合にも \(e=0\) ではあるが, 同様な表式が生成することを示せ.(ただし式 (12-7) は原書のままなので, 結果と合わせるには \(\hbar\) や \(c\) を補足するなどの修正が必要である).

(4) 式 (12-7) は, 原子中に於ける「電子-中性子散乱」を解釈するのに用いることが出来る.原子による中性子散乱のほとんどは, 原子核からの等方散乱(isotropic scattering)である.しかしながら, 原子内電子もやはり中性子を散乱し, そして原子核散乱と干渉する波を引き起こす.遅い中性子の場合, この効果は実験的に観測されている.それは式 (12-6) の項 (5) によって解釈される [式 (12-7) では \(e=0\) に修正される].原子核の外側では電子の電荷が存在するので, \(\nabla\cdot\mb{E}\) はゼロでない値を持つ.ボルン近似によって「中性子-電子散乱」の散乱振幅を計算するとき, 項 (5) を用いることが出来る.しかしながら, その効果が最初に発見されたとき, それはポテンシャル \(c\delta(\mb{R})\) で与えられる「電子-中性子 相互作用」(Electron-Neutron Interaction)を仮定することで説明された.ただし \(\delta\) はディラックのデルタ関数であり, そして \(\mb{R}\) は中性子-電子間の距離である.
ボルン近似によって, ポテンシャルが \(c\delta(\mb{R})\) である場合の散乱振幅を計算せよ.そして, それを項 (5) による振幅と比較せよ.次を示せ:

\begin{equation*}
c=4\pi \mu_{N}\frac{e^{2}}{4M_{N}^{2}}
\end{equation*}

\(c\delta(\mb{R})\) をポテンシャルと解釈する目的で「平均ポテンシャル \(\overline{V}\)」が定義される.それは,「古典電子半径 \(r_e=e^{2}/mc^{2}\)」の球全体に作用することで同じ効果を生成するポテンシャルである.
\(\mu_{N}=-1.9135 e\hbar/2M_{N}\) を用いると, その時の \(V\) は, 述べらている精度内で実験結果, すなわち \(4400\pm 400 \mathrm{eV}\) [4] L.Foldy, Phys. Rev., \(\mb{87}\), 693 (1952). と一致することを示せ.

(5) \(v^{2}/c^{2}\) オーダーの項を無視することで, 次式が成り立つことを示せ:

\begin{equation*}
\int \psi_{f}^{\dagger}\mb{\alpha} f(\mb{R})\,\psi_{f}\,dVol \quad\rightarrow\quad
\int \chi_{f}^{*}\left[\frac{\mb{p}f+f\mb{p}}{2m}+\frac{\mb{\sigma}}{2m}\times \nabla f\right]\,\chi_{f}\,dVol
\end{equation*}


〈 (1)の解答例 〉Landau and Lifshitz:「Quantum Electrodynamics」§34. Fine structure of levels of the hydrogen atom
の文章を抜粋することで, この解答とする.
水素原子, すなわち動かない原子核のクーロン場に於ける電子のエネルギー準位の相対論的補正を求めよう.水素原子内での電子の速度は \(v/c\sim \alpha\ll1\) である.従って必要な補正は, 摂動論を使って, 式 (12-6) に於いて電場しかない (即ち \(\mb{A}=0\) で磁場は存在せず \(\mb{H}=0\)) と仮定した時のハミルトニアン

\begin{equation*}
\hat{H}=\frac{\hat{\mb{p}}^{2}}{2m}+e\phi -\frac{\hat{\mb{p}}^{4}}{8m^{3}c^{2}}
-\frac{e\hbar^{2}}{8m^{2}c^{2}}\nabla\cdot\mb{E}
-\frac{e\hbar}{4m^{2}c^{2}}\mb{\sigma}\cdot\mb{E}\times\hat{\mb{p}}
\tag{33.12}
\end{equation*}

の相対論的項を非摂動状態(つまり非相対論的波動関数)で平均化することで計算出来る.やや一般性を持たせるために, 原子核の電荷を \(Ze\) とし \(Z\alpha\ll 1\) と仮定する.もしも電場が中心対称の場合には, 電荷 \(-e\) を持つ電子の電場中のポテンシャルエネルギーは \(V=-e\phi\) であるから,
\begin{equation*}
\mb{E}=-\frac{d\phi}{dr}\frac{\mb{r}}{r}
=\frac{1}{e}\frac{dV}{dr}\frac{\mb{r}}{r}
\end{equation*}

であり, 式 (33.12) の最後の項(これは「Thomas項」と呼ばれている)は,「スピン-軌道相互作用の演算子」として表わされる:
\begin{align*}
-\frac{e\hbar}{4m^{2}c^{2}}\mb{\sigma}\cdot\mb{E}\times\hat{\mb{p}}
&=-\frac{e\hbar}{4m^{2}c^{2}}\mb{\sigma}\cdot
\left(-\frac{1}{e}\frac{dV}{dr}\frac{\mb{r}}{r}\times\hat{\mb{p}}\right)
=+\frac{1}{2m^{2}c^{2}}\frac{1}{r}\frac{dV}{dr}\left(\frac{\hbar}{2}\mb{\sigma}\right)\cdot(\mb{r}\times\hat{\mb{p}})\\
&=\frac{1}{2m^{2}c^{2}}\frac{1}{r}\frac{dV}{dr}\mb{S}\cdot\mb{L}
\end{align*}

原子核の場(電場)は中心対称で \(\displaystyle\mb{E}=\frac{Ze}{r^{3}}\mb{r}\) であり, そのポテンシャル \(\displaystyle\phi=\frac{Ze}{r}\) は, ポアソン方程式 \(\Delta \phi=-4\pi\rho=-4\pi Ze\delta(\mb{r})\) を満たす.従って, 電子電荷は負であることを用いると式 (33.12) の第4番目の項は次となる:
\begin{align*}
&\nabla\cdot\mb{E}=\nabla\cdot(-\nabla\phi)=-\nabla^{2}\phi =-\Delta\phi=4\pi Ze\delta(\mb{r}),\\
&\rightarrow\quad -\frac{e\hbar^{2}}{8m^{2}c^{2}}\nabla\cdot\mb{E}
=-\frac{(-e)\hbar^{2}}{8m^{2}c^{2}}4\pi Ze\delta(\mb{r})
=\frac{\pi Ze^{2}\hbar^{2}}{2m^{2}c^{2}}\delta(\mb{r})
\end{align*}
この項を最初に考察したのは C.G.Darwin なので, この項は「Darwin項」と呼ばれる.また, この場合の「Thomas項」は \(\displaystyle V=-e\phi=-\frac{Ze^{2}}{r}\) として次となる:

\begin{equation*}
\frac{dV}{dr}=\frac{d}{dr}\left(\frac{-Ze^{2}}{r}\right)=\frac{Ze^{2}}{r^{2}}\quad\rightarrow\quad
\frac{1}{2m^{2}c^{2}}\frac{1}{r}\frac{dV}{dr}\mb{S}\cdot\mb{L}
=\frac{Ze^{2}}{2m^{2}c^{2}r^{3}}\mb{S}\cdot\mb{L}
\end{equation*}

従って, この電場を式 (33.12) の最後の3項に代入し, 摂動演算子は次となる:
\begin{equation*}
\hat{V}=\hat{V}_{(1)}+\hat{V}_{(2)}+\hat{V}_{(3)}
=-\frac{\hat{\mb{p}}^{4}}{8m^{3}c^{2}} + \frac{Ze^{2}}{2m^{2}c^{2}r^{3}}\hat{\mb{S}}\cdot\hat{\mb{L}}
+ \frac{\pi Ze^{2}\hbar^{2}}{2m^{2}c^{2}}\delta(\mb{r})
\tag{34.1}
\end{equation*}

非相対論的なシュレディンガー方程式 \(\displaystyle \left(\frac{\hat{\mb{p}}^{2}}{2m}+V\right)\psi_0=E_0\psi_0\) によれば,
\begin{equation*}
\hat{\mb{p}}^{2}\psi_0 =2m\left(E_0 + \frac{Ze^{2}}{r}\right)\psi_0,\quad
E_0=-\frac{Z^{2}me^{4}}{2\hbar^{2}n^{2}}
\end{equation*}

ただし \(E_0\) は非摂動エネルギー準位であり \(n\) は「主量子数」である [5][ブログ註] ランダウ=リフシッツ:「量子力学」§ 31 クーロン場の中の運動 を参照すべし..従って, 平均値は次となる:
\begin{equation*}
\overline{\mb{p}^{4}}=4m^{2}\overline{\left(E_0 +\frac{Ze^{2}}{r}\right)^{2}}
\end{equation*}

式 (34.1) の第2項, 第3項のような平均量は, 次の公式たちによって計算される (L.Pauling and E.Wilson§21c を見よ):
\begin{align*}
\overline{r^{-1}}&=\frac{Z}{a_0 n^{2}}=\frac{me^{2}Z}{\hbar^{2}n^{2}},\qquad
\overline{r^{-2}}=\frac{Z^{2}}{a_0^{2}n^{3}(l+1/2)}=\frac{(me^{2}Z)^{2}}{\ds{\hbar^{4}n^{3}\bigl(l+1/2\bigr)}},\\
\overline{r^{-3}}&=\frac{Z^{3}}{a_0^{3}n^{3}l(l+1/2)(l+1)}=\frac{(me^{2}Z)^{3}}{\ds{\hbar^{6}n^{3}l\,(l+1)\bigl(l+1/2\bigr)}},
\qquad a_0=\frac{\hbar^{2}}{me^{2}}
\tag{34.2}
\end{align*}
ただし \(a_0\) は「ボーア半径」であり, また最後の式は \(l\ne0\) のときに成り立つ.また, 固有値は \(\mb{J}^{2}=(\mb{L}+\mb{S})^{2}\) から得られる次の公式によって計算される:
\begin{equation*}
\begin{cases}
l\ne0 & \mb{L}\cdot\mb{S}=\ds{\frac{1}{2}\bigl(\mb{J}^{2}-\mb{L}^{2}-\mb{S}^{2}\bigr)}
\quad\rightarrow\quad
\mb{L}\cdot\mb{S}\,\psi=\ds{\frac{\hbar^{2}}{2}\left[j(j+1)-l(l+1)-\frac{3}{4}\right]}, \\
l=0 & \mb{L}\cdot\mb{S}=0
\end{cases}
\end{equation*}

最後に, 第3項の平均は次の公式を用いて遂行される:
\begin{equation*}
\begin{cases}
l=0 : & \ds{\psi(0)=\frac{1}{\sqrt{\pi}}\left(\frac{Z}{a_0 n}\right)^{3/2}
=\frac{1}{\sqrt{\pi}}\left(\frac{Zme^{2}}{\hbar^{2} n}\right)^{3/2}},\\
l\ne 0 : & \psi(0)=0
\end{cases}
\tag{34.3}
\end{equation*}

ただしこれは, シュレディンガー方程式の解
\begin{align*}
\psi_{1s}&=\frac{1}{\sqrt{\pi}}\left(\frac{Z}{a_0}\right)^{3/2}e^{-\sigma},\quad
\psi_{2s}=\frac{1}{4\sqrt{2\pi}}\left(\frac{Z}{a_0}\right)^{3/2}(2-\sigma)e^{-\sigma/2},\\
\psi_{3s}&=\frac{1}{81\sqrt{3\pi}}\left(\frac{Z}{a_0}\right)^{3/2}(27-18\sigma+2\sigma^{2})e^{-\sigma/3},
\quad \dotsb,\quad \mathrm{with}\quad \sigma=\frac{Z}{a_0}r
\end{align*}
によるものである.しかし J.J.Sakurai の§ 3.8 によれば,「ディラック理論による基底状態の波動関数は,
\begin{align*}
\psi_{gd}&=\frac{N}{\sqrt{\pi}}\left(\frac{Z}{a_0}\right)^{3/2}e^{-Zr/a_0}
\times\left(\frac{Zr}{a_0}\right)^{\sqrt{1-(Z\alpha)^{2}}-1}\begin{pmatrix}
\chi^{(s)} \\ \ds{\frac{i\left(1-\sqrt{1-(Z\alpha)^{2}}\right)}{Z\alpha}\frac{(\mb{\sigma}\cdot\mb{x})}{r}\chi^{(s)}}
\end{pmatrix},
\end{align*}
ただし \(\alpha\) は「微細構造定数」,そして \(\chi^{(s)}\) は「Pauliスピノル」であり \(j_{3}=1/2\) か \(j_3=-1/2\) かに依存して
\begin{equation*}
\begin{pmatrix}
1 \\ 0
\end{pmatrix}\quad\mathrm{or}\quad
\begin{pmatrix}
0 \\ 1
\end{pmatrix}
\end{equation*}

となる.また, 規格化定数 \(N\) は \(\Gamma(x)\) をガンマ関数として次である:
\begin{equation*}
N=2^{\sqrt{1-(Z\alpha)^{2}}-1}
\sqrt{\frac{1+\sqrt{1-(Z\alpha)^{2}}}{\Gamma\bigl(1+2\sqrt{1-(Z\alpha)^{2}}\bigr)}},\quad
\mathrm{where}\quad \Gamma(x)=\int_{0}^{\infty}e^{-t}t^{x-1}\,dt
\end{equation*}

この \(N\) は \(Z\alpha\to0\) のとき \(1\) に近づき, 更に \((Zr/a_0)^{\sqrt{1-(Z\alpha)^{2}}-1}\) は, 次のオーダーの距離
\begin{equation*}
r\sim \frac{137\hbar}{2mcZ}e^{-2(137)^{2}/Z^{2}}
\end{equation*}

を除いて本質的に \(1\) であることに注意すべきである.\(r\to0\) のとき \(\psi\) は穏やかな特異性を示す.しかし, それは学術的な興味に過ぎない.なぜなら, 短距離での波動関数は原子核の有限な電荷分布のために修正されるはずだからである.従って, 小さな \(Z\) の水素様原子の基底状態の場合, 上方の2成分波動関数は本質的に”シュレディンガー波動関数にPauliスピノルを掛け合わせたもの”と言える」.

式 (34.1) に於ける各摂動項の「非相対論的状態での平均値」を順番に計算して行こう.まず第1項目の平均は,
\begin{align*}
\left\langle -\frac{\hat{\mb{p}}^{4}}{8m^{3}c^{2}} \right\rangle
&=-\frac{1}{8m^{3}c^{2}}\times \overline{\mb{p}^{4}}
=-\frac{1}{8m^{3}c^{2}}\times 4m^{2}\overline{\left(E_0 +\frac{Ze^{2}}{r}\right)^{2}}\\
&=-\frac{1}{2mc^{2}}\overline{\left(-\frac{Z^{2}me^{4}}{2\hbar^{2}n^{2}} +\frac{Ze^{2}}{r}\right)^{2}}\\
&=-\frac{1}{2mc^{2}}\left(Z^{2}e^{4}\overline{r^{-2}}-\frac{mZ^{3}e^{6}}{\hbar^{2}n^{2}}\overline{r^{-1}}
+\frac{m^{2}Z^{4}e^{8}}{4\hbar^{4}n^{4}}\right)\\
&=-\frac{1}{2mc^{2}}\left(\frac{m^{2}e^{8}Z^{4}}{\hbar^{4}n^{3}\left(l+1/2\right)}
-\frac{m^{2}Z^{4}e^{8}}{\hbar^{4}n^{4}}+\frac{m^{2}Z^{4}e^{8}}{4\hbar^{4}n^{4}}\right)\\
&=-\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{l+1/2}-\frac{3}{4n}\right)
\tag{1}
\end{align*}
第2項目すなわち「Thomas項」の平均は,
\begin{align*}
\left\langle\frac{Ze^{2}}{2m^{2}c^{2}r^{3}}\hat{\mb{S}}\cdot\hat{\mb{L}}\right\rangle
&=\frac{Ze^{2}}{2m^{2}c^{2}}\,\overline{r^{-3}}\times \frac{\hbar^{2}}{2}
\left[j(j+1)-l(l+1)-\frac{3}{4}\right]\\
&=\frac{Ze^{2}}{2m^{2}c^{2}}\times\frac{(me^{2}Z)^{3}}{\hbar^{6}n^{3}l\,(l+1)(l+1/2)}
\times \frac{\hbar^{2}}{2}\left[j(j+1)-l(l+1)-\frac{3}{4}\right]\\
&=\frac{mZ^{4}e^{8}}{4\hbar^{4}c^{2}n^{3}}\times\frac{j(j+1)-l(l+1)-3/4}{l\,(l+1)(l+1/2)}
\end{align*}
ただし, この「Thomas項」は S 状態 \((l=0)\) のときにはゼロになってしまう.それ以外 \((l\ne0)\) の場合では, \(j\) の値は \(j=l+1/2\) または \(j=j-1/2\) のみを取るから,

\begin{equation*}
j(j+1)-l(l+1)-\frac{3}{4}=
\begin{cases}
l & \ :\ j=l+1/2 \\
-(l+1) & \ :\ j=l-1/2
\end{cases}
\end{equation*}

従って,
\begin{align*}
\left\langle\frac{Ze^{2}}{2m^{2}c^{2}r^{3}}\hat{\mb{S}}\cdot\hat{\mb{L}}\right\rangle
&=\begin{cases}
\ds{+\frac{mZ^{4}e^{8}}{4\hbar^{4}c^{2}n^{3}}\frac{1}{(l+1)(l+1/2)}} & \ :\ \ds{j=l+\frac{1}{2}} \\
\ds{-\frac{mZ^{4}e^{8}}{4\hbar^{4}c^{2}n^{3}}\frac{1}{l(l+1/2)}} & \ :\ \ds{j=l-\frac{1}{2}}
\end{cases}
\tag{2}
\end{align*}
第3の項, すなわち最後の項は S 状態以外 (\(l\ne0\)) ではゼロになってしまう:
\begin{align*}
\left\langle\frac{\pi Ze^{2}\hbar^{2}}{2m^{2}c^{2}}\delta(\mb{r})\right\rangle
&=\int \psi^{*}(r)\frac{\pi Ze^{2}\hbar^{2}}{2m^{2}c^{2}}\delta(\mb{r})\,\psi(r)\,dV
=\frac{\pi Ze^{2}\hbar^{2}}{2m^{2}c^{2}}|\psi(0)|^{2}\\
&=\begin{cases}
\ds{\frac{\pi Ze^{2}\hbar^{2}}{2m^{2}c^{2}}\times \frac{1}{\pi}\left(\frac{Zme^{2}}{\hbar n}\right)^{3}
=\frac{mZ^{4}e^{8}}{2\hbar^{2}c^{2}n^{3}}} & \ : l=0 \\
0 & \ : l\ne 0
\end{cases}
\tag{3}
\end{align*}
以上の結果をまとめる.まず S 状態 \((l=0)\) の場合では, 式 (1) と式 (3) を足し合わせて次となる:
\begin{equation*}
\Delta E_{(l=0)}=-\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{0+1/2}-\frac{3}{4n}\right)
+\frac{mZ^{4}e^{8}}{2\hbar^{2}c^{2}n^{3}}
=-\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{1}-\frac{3}{4n}\right)
\end{equation*}

これは \(j=1/2\) とすれば, 次のように書くことが出来ることに注意する:
\begin{equation*}
\Delta E_{(l=0)}=-\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{j+1/2}-\frac{3}{4n}\right)
\tag{4}
\end{equation*}

次に S 状態以外 \((l\ne0)\) を考えると, 式 (1) と式 (2) を足し合わせて次となる:
まず \(j=l+1/2\) の場合すなわち \(l=j-1/2\) の場合には,
\begin{align*}
\Delta E_{(j=l+1/2)} &= -\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{l+1/2}-\frac{3}{4n}\right)
+\frac{mZ^{4}e^{8}}{4\hbar^{4}c^{2}n^{3}}\frac{1}{(l+1)(l+1/2)}\\
&=-\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{l+1}-\frac{3}{4n}\right)
=-\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{j+1/2}-\frac{3}{4n}\right)
\end{align*}
また \(j=l-1/2\) の場合すなわち \(l=j+1/2\) の場合には,
\begin{align*}
\Delta E_{(j=l-1/2)} &= -\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{l+1/2}-\frac{3}{4n}\right)
-\frac{mZ^{4}e^{8}}{4\hbar^{4}c^{2}n^{3}}\frac{1}{l(l+1/2)}\\
&=-\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{l}-\frac{3}{4n}\right)
=-\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{j+1/2}-\frac{3}{4n}\right)
\end{align*}
つまり \(j=l\pm1/2\) の両方の場合が, 同じ表現として表せることが分かる:
\begin{equation*}
\Delta E_{(l\ne0)}
=-\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{j+1/2}-\frac{3}{4n}\right)
\tag{5}
\end{equation*}

以上の結果式 (4) と式 (5) とから, あらゆる \(j\) 及び \(l\) の場合で, 結果は同じ形に表せることが分かる:
\begin{align*}
\Delta E= \langle \hat{V}\rangle
&=-\frac{mZ^{4}e^{8}}{2\hbar^{4}c^{2}n^{3}}\left(\frac{1}{j+1/2}-\frac{3}{4n}\right)\\
&=\frac{E_{0}(Z\alpha)^{2}}{n}\left(\frac{1}{j+1/2}-\frac{3}{4n}\right),\quad
E_0=-\frac{Z^{2}e^{2}}{2a_0 n^{2}}=-\frac{Z^{2}me^{4}}{2\hbar^{2}n^{2}}
\tag{6}
\end{align*}
ただし \(\alpha\) は「微細構造定数」, そして \(a_0\) は「ボーア半径」である:
\begin{align*}
\alpha&=\frac{e^{2}}{\hbar c}=0.0072973525643(11),\quad \frac{1}{\alpha}=137.04\\
a_0&=\frac{\hbar^{2}}{me^{2}}=5.292\times 10^{-9}\,\mathrm{cm}
\end{align*}

〈 (2)の解答例 〉\(\mb{A}\) と \(\phi\) が時間依存するとして, 式 (12-5) 以降を修正して行く.
式 (12-5) の第3項は, 公式 (F1) に於いて \(A=\mb{\sigma}\cdot\mb{\pi}\) そして \(B=(W-V)\) とした

\begin{equation*}
\Bigl[(\mb{\sigma}\cdot\mb{\pi}),\,\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr]\Bigr]
\tag{7}
\end{equation*}

を用いることが出来る.従って, 式 (12-5) は次となる:
\begin{align*}
(W-V)\,\chi
&=\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{2m}\chi-\frac{(\mb{\sigma}\cdot\mb{\pi})^{4}}{8m^{3}c^{2}}\chi
+\frac{1}{8m^{2}c^{2}}\Bigl[(\mb{\sigma}\cdot\mb{\pi}),\,
\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr]\Bigr]\chi
\tag{12-5}
\end{align*}
今度は \(W=i\hbar\partial/\partial t\) とし \(\mb{A}\) は時間依存するのであるから, この式 (7) は修正して行かなければならない.まず, 演算子の後ろには位置 \(\mb{x}\) と時間 \(t\) の関数 \(\chi(\mb{x},t)\) が存在することを考慮すると,
\begin{align*}
\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr]
&=\left[\mb{\sigma}\cdot\left(\mb{p}-\frac{e}{c}\mb{A}\right),\,(W-V)\right]\\
&=\left[\mb{\sigma}\cdot\mb{p},\,(W-V)\right]-\frac{e}{c}\left[\mb{\sigma}\cdot\mb{A},\,(W-V)\right],\\
\mathrm{where}\quad \left[\mb{\sigma}\cdot\mb{p},\,(W-V)\right]&=\mb{\sigma}\cdot\mb{p}(W-V)-(W-V)\mb{\sigma}\cdot\mb{p}\\
&=\cancel{\mb{\sigma}\cdot\mb{p}W}-\mb{\sigma}\cdot\mb{p}V-\bcancel{V\mb{\sigma}\cdot\mb{p}}-\cancel{W\mb{\sigma}\cdot\mb{p}}
+\bcancel{V\mb{\sigma}\cdot\mb{p}}\\
&=-\mb{\sigma}\cdot(\mb{p}V),\\
\left[\mb{\sigma}\cdot\mb{A},\,(W-V)\right]&=\mb{\sigma}\cdot\mb{A}(W-V)-(W-V)\mb{\sigma}\cdot\mb{A}\\
&=\cancel{\mb{\sigma}\cdot\mb{A}W}-\bcancel{\mb{\sigma}\cdot\mb{A}V}-\mb{\sigma}\cdot(W\mb{A})
-\cancel{\mb{\sigma}\cdot\mb{A}W}+\bcancel{V\mb{\sigma}\cdot\mb{A}}\\
&=-\mb{\sigma}\cdot(W\mb{A}),\\
\therefore\quad
\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr]&=-\mb{\sigma}\cdot(\mb{p}V)+\frac{e}{c}\mb{\sigma}\cdot(W\mb{A})
\tag{8}
\end{align*}
である.すると,式 (7) は
\begin{align*}
\left[(\mb{\sigma}\cdot\mb{\pi}),\,\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr]\right]
&=\left[\mb{\sigma}\cdot\mb{\pi},\,-\mb{\sigma}\cdot(\mb{p}V)+\frac{e}{c}\mb{\sigma}\cdot(W\mb{A})\right]\\
&=-\bigl[\mb{\sigma}\cdot\mb{\pi},\,\mb{\sigma}\cdot(\mb{p}V)\bigr]
+\frac{e}{c}\bigl[\mb{\sigma}\cdot\mb{\pi},\,\mb{\sigma}\cdot(W\mb{A})\bigr]
\tag{9}
\end{align*}
このとき, 公式 (F2) から
\begin{align*}
&\bigl[\mb{\sigma}\cdot\mb{\pi},\,\mb{\sigma}\cdot(\mb{p}V)\bigr]
=(\mb{\sigma}\cdot\mb{\pi})(\mb{\sigma}\cdot\mb{p}V)-(\mb{\sigma}\cdot\mb{p}V)(\mb{\sigma}\cdot\mb{\pi})\\
&\quad =\mb{\pi}\cdot\mb{p}V+i\mb{\sigma}\cdot(\mb{\pi}\times\mb{p}V)-(\mb{p}V)\cdot\mb{\pi}
-i\mb{\sigma}\cdot\bigl\{(\mb{p}V)\times\mb{\pi}\bigr\}\\
&\quad =\left(\mb{p}-\frac{e}{c}\mb{A}\right)\cdot(\mb{p}V)-(\mb{p}V)\cdot\left(\mb{p}-\frac{e}{c}\mb{A}\right)
+i\mb{\sigma}\cdot\left\{\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times(\mb{p}V)\right\}
-i\mb{\sigma}\cdot\left\{(\mb{p}V)\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)\right\}\\
&\quad =\mb{p}\cdot(\mb{p}V)+\cancel{(\mb{p}V)\cdot\mb{p}}-\cancel{\frac{e}{c}\mb{A}\cdot(\mb{p}V)}
-\cancel{(\mb{p}V)\cdot\mb{p}}+\cancel{\frac{e}{c}(\mb{p}V)\cdot\mb{A}}
+i\mb{\sigma}\cdot\left\{\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times\mb{p}V
-(\mb{p}V)\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)\right\}\\
&\quad =\mb{p}\cdot(\mb{p}V)+i\mb{\sigma}\cdot\left\{\cancel{\mb{p}\times\mb{p}V}-(\mb{p}V)\times\mb{p}-\frac{e}{c}\mb{A}\times(\mb{p}V)
-(\mb{p}V)\times\mb{p}-\frac{e}{c}\mb{A}\times(\mb{p}V)\right\}\\
&\quad =\mb{p}\cdot(\mb{p}V)+2i\mb{\sigma}\cdot\left\{\mb{p}\times(\mb{p}V)-\frac{e}{c}\mb{A}\times(\mb{p}V)\right\}
=\mb{p}\cdot(\mb{p}V)+2i\mb{\sigma}\cdot\left\{\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times\mb{p}V\right\}
\tag{10}
\end{align*}
同様にして,
\begin{align*}
&\bigl[\mb{\sigma}\cdot\mb{\pi},\,\mb{\sigma}\cdot(W\mb{A})\bigr]
=(\mb{\sigma}\cdot\mb{\pi})(\mb{\sigma}\cdot(W\mb{A}))-(\mb{\sigma}\cdot(W\mb{A}))(\mb{\sigma}\cdot\mb{\pi})\\
&\quad =\mb{\pi}\cdot(W\mb{A})+i\mb{\sigma}\cdot\bigl\{\mb{\pi}\times(W\mb{A})\bigr\}
-(W\mb{A})\cdot\mb{\pi}-i\mb{\sigma}\cdot\bigl\{(W\mb{A})\times\mb{\pi}\bigr\}\\
&\quad =\left(\mb{p}-\frac{e}{c}\mb{A}\right)\cdot(W\mb{A})
+i\mb{\sigma}\cdot\left\{\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times(W\mb{A})\right\}
-(W\mb{A})\cdot\left(\mb{p}-\frac{e}{c}\mb{A}\right)
-i\mb{\sigma}\cdot\left\{(W\mb{A})\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)\right\}\\
&\quad = \mb{p}\cdot(W\mb{A})+\cancel{(W\mb{A})\cdot\mb{p}}-\cancel{\frac{e}{c}\mb{A}\cdot(W\mb{A})}
-\cancel{(W\mb{A})\cdot\mb{p}}+\cancel{\frac{e}{c}(W\mb{A})\cdot\mb{A}}\\
&\qquad\qquad\quad +i\mb{\sigma}\cdot\left\{\mb{p}\times(W\mb{A})-(W\mb{A})\times\mb{p}-\frac{e}{c}\mb{A}\times(W\mb{A})
-(W\mb{A})\times\mb{p}+\frac{e}{c}(W\mb{A})\times\mb{A}\right\}\\
&\quad =\mb{p}\cdot(W\mb{A})+i\mb{\sigma}\cdot\left\{\mb{p}\times(W\mb{A})
-2(W\mb{A})\times\mb{p}+2(W\mb{A})\times\frac{e}{c}\mb{A}\right\}\\
&\quad = \mb{p}\cdot(W\mb{A})+i\mb{\sigma}\cdot\left\{\mb{p}\times(W\mb{A})
-2(W\mb{A})\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)\right\}
\tag{11}
\end{align*}
式 (10) と式 (11) から, 式 (9) は,
\begin{align*}
&\Bigl[(\mb{\sigma}\cdot\mb{\pi}),\,\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr]\Bigr]
=-\bigl[\mb{\sigma}\cdot\mb{\pi},\,\mb{\sigma}\cdot(\mb{p}V)\bigr]
+\frac{e}{c}\bigl[\mb{\sigma}\cdot\mb{\pi},\,\mb{\sigma}\cdot(W\mb{A})\bigr]\\
&\quad =-\mb{p}\cdot(\mb{p}V)-2i\mb{\sigma}\cdot\left\{\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times\mb{p}V\right\}
+\frac{e}{c}\mb{p}\cdot(W\mb{A})
+i\frac{e}{c}\mb{\sigma}\cdot\left\{\mb{p}\times(W\mb{A})-2(W\mb{A})\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)\right\}\\
&\quad = -\mb{p}\cdot(\mb{p}V)+\frac{e}{c}\mb{p}\cdot(W\mb{A})
+i\mb{\sigma}\cdot\left\{\frac{e}{c}\mb{p}\times(W\mb{A})-2\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times(\mb{p}V)
+2\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times\frac{e}{c}(W\mb{A})\right\}\\
&\quad =-\mb{p}\cdot\left\{(\mb{p}V)-\frac{e}{c}(W\mb{A})\right\}
+i\mb{\sigma}\cdot\left[\frac{e}{c}\mb{p}\times(W\mb{A})-2\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times
\left\{(\mb{p}V)-\frac{e}{c}(W\mb{A})\right\}\right]
\tag{12}
\end{align*}
従って, 式 (12-5) は次となる:
\begin{align*}
(W-V)\,\chi&=\frac{(\mb{\sigma}\cdot\mb{\pi})^{2}}{2m}\chi
-\frac{(\mb{\sigma}\cdot\mb{\pi})^{4}}{8m^{3}c^{2}}\chi
-\frac{1}{8m^{2}c^{2}}\mb{p}\cdot\left\{(\mb{p}V)-\frac{e}{c}(W\mb{A})\right\}\,\chi\\
&\quad +\frac{1}{8m^{2}c^{2}}i\mb{\sigma}\cdot\left[\frac{e}{c}\mb{p}\times(W\mb{A})-2\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times
\left\{(\mb{p}V)-\frac{e}{c}(W\mb{A})\right\}\right]\chi
\tag{13}
\end{align*}
ここで \(V=e\phi\) とし, 式中の量に \(\mb{p}=-i\hbar\nabla\) と \(W=i\hbar\partial/\partial t\) を代入してみる.電場 \(\mb{E}\) はポテンシャルを用いて \(\displaystyle\mb{E}=-\frac{1}{c}\ppdiff{\mb{A}}{t}-\nabla\phi\) と表せるから,
\begin{align*}
(\mb{p}V)-\frac{e}{c}(W\mb{A})&=-i\hbar\nabla V -\frac{e}{c}\left(i\hbar\pdiff{t}\mb{A}\right)
=-i\hbar \nabla e\phi -i\hbar\frac{e}{c}\ppdiff{\mb{A}}{t}
=ie\hbar \left(-\frac{1}{c}\ppdiff{\mb{A}}{t}-\nabla \phi\right)\\
&=ie\hbar\,\mb{E}
\tag{14}
\end{align*}
同様にして,
\begin{align*}
\frac{e}{c}\mb{p}\times(W\mb{A})&=\frac{e}{c}\mb{p}\times\left(i\hbar\ppdiff{\mb{A}}{t}\right)
=-ie\hbar\mb{p}\times\left(-\frac{1}{c}\ppdiff{\mb{A}}{t}\right)
=-ie\hbar\mb{p}\times \bigl(\mb{E}+\nabla\phi\bigr)\\
&=-ie\hbar\cdot(-i\hbar\nabla)\times\mb{E}-ie\hbar\cdot(-i\hbar\nabla)\times\nabla\phi\\
&=-e\hbar^{2}\nabla\times\mb{E}-\cancel{e\hbar^{2}\nabla\times\nabla \phi}
=-e\hbar^{2}\nabla\times \mb{E}
\end{align*}
また \(\displaystyle \mathrm{rot}\,\mb{E}=\nabla\times\mb{E}=-\frac{1}{c}\ppdiff{\mb{H}}{t}\) である.従って,
\begin{equation*}
\frac{e}{c}\mb{p}\times(W\mb{A})=-e\hbar^{2}\nabla\times \mb{E}
=\frac{e\hbar^{2}}{c}\ppdiff{\mb{H}}{t}
\tag{15}
\end{equation*}

であることが分かる.よってこの場合, 本文の式 (12.6) は次のように修正される:
\begin{align*}
W\chi &=\left[V+\frac{1}{2m}\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}-\frac{e\hbar}{2mc}\mb{\sigma}\cdot\mb{H}
-\frac{\mb{p}^{4}}{8m^{3}c^{2}} -\frac{ie\hbar}{8m^{2}c^{2}}\mb{p}\cdot\mb{E}+\frac{1}{8m^{2}c^{2}}i\mb{\sigma}\cdot\left\{\frac{e\hbar^{2}}{c}\ppdiff{\mb{H}}{t}
-2ie\hbar\left(\mb{p}-\frac{e}{c}\mb{A}\right)\times\mb{E}\right\}\right]\chi \\
&=\left[V+\frac{1}{2m}\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}-\frac{e\hbar}{2mc}\mb{\sigma}\cdot\mb{H}
-\frac{\mb{p}^{4}}{8m^{3}c^{2}} -\frac{e\hbar^{2}}{8m^{2}c^{2}}\nabla\cdot\mb{E}-\frac{e\hbar}{4m^{2}c^{2}}\mb{\sigma}\cdot \mb{E}\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)
+\frac{ie\hbar^{2}}{8m^{2}c^{3}}\mb{\sigma}\cdot\ppdiff{\mb{H}}{t}\right]\chi
\tag{16}
\end{align*}
すなわち, 最後の項だけが付加することが分かる:
\begin{equation*}
\frac{ie\hbar^{2}}{8m^{2}c^{3}}\mb{\sigma}\cdot\ppdiff{\mb{H}}{t}
\tag{17}
\end{equation*}


〈 (3)の解答例 〉題意の式の場合について, 本文の式 (11-4) と式 (11-5) に相当する式を求めてみる:
\begin{align*}
i\hbar\pdiff{t}\psi =H\,\psi &= \left[c\mb{\alpha}\cdot\left(\mb{p}
-\frac{e}{c}\mb{A}\right) +\beta mc^{2} + e\phi+\frac{i\mu}{2}\beta\gamma^{\mu}\gamma^{\nu}F_{\mu\nu}\right]\,\psi\\
&= \left[c\mb{\alpha}\cdot\left(\mb{p}
-\frac{e}{c}\mb{A}\right) +\beta mc^{2} + e\phi
+ i\mu (\beta\mb{\alpha}\cdot\mb{E}+i\,\beta\mb{\Sigma}\cdot\mb{H})\right]\,\psi,\\
\therefore\quad W\psi&=\left[c\mb{\alpha}\cdot\left(\mb{p}
-\frac{e}{c}\mb{A}\right) +\beta mc^{2} + V
+i\mu (\beta\mb{\alpha}\cdot\mb{E}+i\,\beta\mb{\Sigma}\cdot\mb{H})\right]\,\psi
\tag{18}
\end{align*}
ハミルトニアン中の第4項を計算する.まず,
\begin{align*}
&\beta=\gamma^{0}=\begin{pmatrix} 1 & 0 \\ 0 & -1\end{pmatrix},\
\mb{\alpha}=\begin{pmatrix} 0 & \mb{\sigma} \\ \mb{\sigma} & 0\end{pmatrix},\
\mb{\Sigma}=\begin{pmatrix} \mb{\sigma} & 0\\ 0 & \mb{\sigma}\end{pmatrix}
\ \rightarrow\
\beta\mb{\alpha}=\begin{pmatrix} 0 & \mb{\sigma} \\ -\mb{\sigma} & 0\end{pmatrix},\quad
\beta\mb{\Sigma}=\begin{pmatrix} \mb{\sigma} & 0\\ 0 & -\mb{\sigma}\end{pmatrix}
\tag{19}
\end{align*}
従って,
\begin{align*}
i\mu (\beta\mb{\alpha}\cdot\mb{E}+i\,\beta\mb{\Sigma}\cdot\mb{H})
&=i\mu\mb{E}\cdot(\beta\mb{\alpha})-\mu\mb{H}\cdot(\beta\mb{\Sigma})
=i\mu\mb{E}\cdot\begin{pmatrix} 0 & \mb{\sigma} \\ -\mb{\sigma} & 0\end{pmatrix}
-\mu\mb{H}\cdot\begin{pmatrix} \mb{\sigma} & 0\\ 0 & -\mb{\sigma}\end{pmatrix}\\
&=-\mu \begin{pmatrix}
\mb{\sigma}\cdot\mb{H} & -i\mb{\sigma}\cdot\mb{E} \\ i\mb{\sigma}\cdot\mb{E} & -\mb{\sigma}\cdot\mb{H}
\tag{20}
\end{pmatrix}
\end{align*}
第11講の本文と同様に \(\psi\) を2成分 \(\psi_a\) 及び \(\psi_b\) で表わすと, 上式 (18) は次の2式に書くことが出来る:
\begin{align*}
W\psi_a =c\mb{\sigma}\cdot\mb{\pi}\psi_b +V\psi_a -\mu\mb{\sigma}\cdot\mb{H}\psi_a +i\mu\mb{\sigma}\cdot\mb{E}\psi_b
\tag{11-4′}\\
(2mc^{2}+W)\psi_b = c\mb{\sigma}\cdot\mb{\pi}\psi_a +V\psi_b -i\mu\mb{\sigma}\cdot\mb{E}\psi_a +\cancel{\mu\mb{\sigma}\cdot\mb{H}\psi_b}
\tag{11-5′}
\end{align*}
朝永:「スピンはめぐる」には, この問題に相当する事柄が「異常磁気能率の量子力学的取り扱い法」として解説されている.その p.247 の記述を参照して, 上式 (11-5′) 右辺の第4項目は捨てることにする.すると, 上の2式からは次が得られる:
\begin{align*}
&(W-V+\mu\mb{\sigma}\cdot\mb{H})\psi_a = (c\mb{\sigma}\cdot\mb{\pi}+i\mu\mb{\sigma}\cdot\mb{E})\psi_b,\tag{11-4”} \\
&(2mc^{2}+W-V)\psi_b \approx (c\mb{\sigma}\cdot\mb{\pi} -i\mu\mb{\sigma}\cdot\mb{E})\psi_a,
\quad \rightarrow\quad
\psi_b =\frac{c\mb{\sigma}\cdot\mb{\pi}-i\mu\mb{\sigma}\cdot\mb{E}}{2mc^{2}+W-V}\psi_a
\tag{11-6′}
\end{align*}
式 (11-6′) を式 (11-4”) に代入すると,
\begin{align*}
(W-V+\mu\mb{\sigma}\cdot\mb{H})\psi_a &=c\mb{\sigma}\cdot\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\times
\frac{\ds{c\mb{\sigma}\cdot\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)}}{2mc^{2}+W-V}\psi_a\\
&=\mb{\sigma}\cdot\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\left(\frac{c^{2}}{2mc^{2}+W-V}\right)
\mb{\sigma}\cdot\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)\psi_a
\end{align*}
第12講と同様に近似式 (12-1) を用いるならば, 式 (12-2) に相当する式として次が得られる:
\begin{align*}
(W-V)\psi_a
&=-\mu\mb{\sigma}\cdot\mb{H}\psi_a
+ \frac{1}{2m}\mb{\sigma}\cdot\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\,
\mb{\sigma}\cdot\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)\psi_a\\
&\qquad -\frac{1}{4m^{2}c^{2}}\mb{\sigma}\cdot\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)(W-V)\,
\mb{\sigma}\cdot\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)\,\psi_a
\tag{12-2′}
\end{align*}
パウリの付加項は「微小量」と考える.そこで規格化のための変換式は, 式 (12-4) を近似的に用いてしまうことにする.すると式 (12-5) に相当する式は
\begin{align*}
(W-V)\chi&=\left[-\mu\mb{\sigma}\cdot\mb{H}
+ \frac{1}{2m}\mb{\sigma}\cdot\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\,
\mb{\sigma}\cdot\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)\right.\\
&\qquad \left.-\frac{1}{8m^{3}c^{2}}\left\{\mb{\sigma}\cdot\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\right\}^{2}
\left\{\mb{\sigma}\cdot\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)\right\}^{2}
+\frac{1}{8m^{2}c^{2}}\Bigl[(\mb{\sigma}\cdot\mb{\pi}),\,\bigl[(\mb{\sigma}\cdot\mb{\pi}),\,(W-V)\bigr]\Bigr]
\right]\chi
\tag{12-5′}
\end{align*}
公式 (F1) を用いると,

\begin{equation*}
\frac{1}{2m}\mb{\sigma}\cdot\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\,
\mb{\sigma}\cdot\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)
=\frac{1}{2m}\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\cdot\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)
+\frac{i}{2m}\mb{\sigma}\cdot\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\times\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)
\tag{21}
\end{equation*}

簡単化のために「\(\mb{A}=0\)」, また「磁場は定常的である:\(\displaystyle\nabla\times\mb{E}=-\frac{1}{c}\ppdiff{\mb{H}}{t}=0\to\mb{p}\times\mb{E}=0\)」とすると,
\begin{align*}
&\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\cdot\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)
=\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}-\frac{i\mu}{c}\mb{p}\cdot\mb{E}
=\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}-\frac{\hbar\mu}{c}\nabla\cdot\mb{E},\\
\rightarrow&\quad
\frac{1}{2m}\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\cdot\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)
=\frac{1}{2m}\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}-\frac{\mu\hbar}{2mc}\nabla\cdot\mb{E},\\
&\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\times\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)
=\mb{\pi}\times\mb{\pi}-\cancel{\frac{i\mu}{c}\mb{\pi}\times\mb{E}}+\frac{2i\mu}{c}\mb{E}\times\mb{\pi}
=i\frac{e\hbar}{c}\mb{H}+\frac{2i\mu}{c}\mb{E}\times\left(\mb{p}-\frac{e}{c}\mb{A}\right),\\
\rightarrow&\quad
\frac{i}{2m}\mb{\sigma}\cdot (\mb{\pi}-2\mu\mb{E})\times(\mb{\pi}+2\mu\mb{E})
=-\frac{e\hbar}{2mc}\mb{\sigma}\cdot\mb{H}-\frac{\mu}{mc}\mb{E}\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)
\end{align*}
すると式 (21) は,
\begin{align*}
\frac{1}{2m}\mb{\sigma}\cdot\left(\mb{\pi}+\frac{i\mu}{c}\mb{E}\right)\,\mb{\sigma}
\cdot\left(\mb{\pi}-\frac{i\mu}{c}\mb{E}\right)
&=\frac{1}{2m}\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}-\frac{\mu\hbar}{2mc}\nabla\cdot\mb{E}
-\frac{e\hbar}{2mc}\mb{\sigma}\cdot\mb{H}-\frac{\mu}{mc}\mb{\sigma}\cdot\mb{E}\times
\left(\mb{p}-\frac{e}{c}\mb{A}\right)
\end{align*}
式 (12-5′) 中の第3項と第4項は式 (12-5) と同じ近似式とする.以上から, 式 (12-5′) は
\begin{align*}
(W-V)\chi&=\left[-\mu\mb{\sigma}\cdot\mb{H}+\frac{1}{2m}\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}
-\frac{\mu}{2mc}\hbar\nabla\cdot\mb{E}-\frac{e\hbar}{2mc}\mb{\sigma}\cdot\mb{H}\right.\\
&\qquad\left.-\frac{\mu}{mc}\mb{\sigma}\cdot\mb{E}\times
\left(\mb{p}-\frac{e}{c}\mb{A}\right)-\frac{\mb{p}^{4}}{8m^{3}c^{2}}
-\frac{e}{8m^{2}c^{2}}\left\{\hbar^{2}\nabla\cdot\mb{E}+2\hbar\mb{\sigma}\cdot\mb{E}\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)\right\}\right]
\end{align*}
従って \(\displaystyle\mu=\frac{e\hbar}{2mc}\kappa\) とすると, 式 (12-6) に相当する式は次となる:
\begin{align*}
W\chi &=\left[V+\frac{1}{2m}\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}
-\left(\mu + \frac{e\hbar}{2mc}\right)\mb{\sigma}\cdot\mb{H}-\frac{\mb{p}^{4}}{8m^{3}c^{2}}\right.\\
&\left.\qquad\qquad\qquad\qquad\qquad -\frac{1}{4mc}\left(2\mu +\frac{e\hbar}{2mc}\right)
\left\{\hbar\nabla\cdot\mb{E}+2\mb{\sigma}\cdot\mb{E}\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)\right\}
\right]\,\chi\\
&=\left[V+\frac{1}{2m}\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{2}
-\frac{e\hbar}{2mc}(1+\kappa)\mb{\sigma}\cdot\mb{H}-\frac{\mb{p}^{4}}{8m^{3}c^{2}}\right.\\
&\left.\qquad\qquad\qquad\qquad\qquad -\frac{1}{4mc}\frac{e\hbar}{2mc}(1+2\kappa)
\left\{\hbar\nabla\cdot\mb{E}+2\mb{\sigma}\cdot\mb{E}\times\left(\mb{p}-\frac{e}{c}\mb{A}\right)\right\}
\right]\,\chi
\tag{12-6′}
\end{align*}

以下は「J.J.Sakurai:Advanced Quantum Mechanics §3.5」からの抜粋である:
『1947 年に P.Kusch が初めて示したように, 実験で観測される電子の磁気回転比は正確に 2 ではなく

\begin{equation*}
g=2\left[1+\left(\frac{e^{2}}{4\pi\hbar c}\right)\frac{1}{2\pi}+\dotsb \right]
\end{equation*}

であり, これはミューオンでも同様である.余分な磁気モーメントの起源の納得行く説明がなされたのは 1947 年の J.Schwinger によってであった.彼は, 物理的な電子は仮想光子を放出または吸収することが出来るという事実を考慮したのだった.磁気能率が \(g=2\) によって正確に与えられない場合, 相互作用ハミルトニアンに現象論的な項を付加することが出来る:
\begin{equation*}
H_{\mathrm{anomalous}}=-\frac{e\hbar}{2mc}\kappa \left(\frac{1}{2}\gamma^{\mu}\gamma^{\nu}F_{\mu\nu}\right)
\end{equation*}

これは「異常磁気能率(またはPauli磁気能率)相互作用」と呼ばれる.すると, 全磁気能率は次となる:
\begin{equation*}
\mu = \frac{e\hbar}{2mc}\bigl(1+\kappa\bigr)
\end{equation*}

陽子の場合には \(e\) と \(m\) を \(|e|\) と陽子質量 \(m_p\) で置き換えて次となる:
\begin{equation*}
\mu_p = \frac{|e|\hbar}{2m_p c}\bigl(1+\kappa_p \bigr)=\mu_{N}\bigl(1+\kappa_p\bigr),\quad
\mu_N \equiv \frac{|e|\hbar}{2m_p c}
\end{equation*}

ただし \(\mu_N\) は「核磁子」と呼ばれる.また \(\kappa_p=1.79\) であり, これは O.Stern が最初に実験的に決定した.従って,「観測される陽子の磁気モーメントのおよそ \(60\%\) が異常である」.上式で \(e\to0\) とするとき, スピン \(1/2\) 粒子の磁気モーメントはゼロになるはずである.しかし実験的には中性子や \(\Lambda\)-ハイぺロンのような中性な粒子さえも, かなり大きな磁気モーメントを持つことが知られている.それらを現象論的に表現するには, やはり上式の形の相互作用によって行うことが出来る』.


〈 (4)の解答例 〉場 \(U(r)\) 中にある換算質量 \(m\) を持つ一粒子に向かって, 運動量 \(\mb{p}=\hbar\mb{k}\) を持って入射する粒子の波動関数を平面波 \(\psi^{(0)}=e^{i\mb{k}\cdot\mb{r}}\) とするとき, 散乱過程を記述するシュレディンガー方程式の遠距離での漸近形の解 \(\psi\) は, ボルン近似の場合では次のように書き表せる:
\begin{align*}
\psi =\psi^{(0)}+\psi^{(1)} &\approx e^{i\mb{k}\cdot\mb{r}} + \frac{f(\theta)}{r}e^{ikr}
=e^{i\mb{k}\cdot\mb{r}}-\frac{m}{2\pi\hbar^{2}}\frac{e^{ikr}}{r}
\int U(\mb{r}’)\,e^{i(\mb{k}’-\mb{k})\cdot\mb{r}’}\,dV’
\end{align*}
このとき, 散乱角 \(\theta\) の関数 \(f(\theta)\) は「散乱振幅」と呼ばれる:

\begin{equation*}
f(\theta)=-\frac{m}{2\pi\hbar^{2}}\int U(\mb{r}’)\,e^{i\mb{q}\cdot\mb{r}’}\,dV’,\qquad
\mb{q}=\mb{k}’-\mb{k},\quad q=2k\sin\frac{\theta}{2}
\end{equation*}

そこで, ポテンシャル \(U(\mb{r})\) が \(c\delta(\mb{r})\) の場合の散乱振幅 \(f_1\) を計算すると,
\begin{equation*}
f_1 =-\frac{m}{2\pi\hbar^{2}}\int c\delta(\mb{r})\,e^{i\mb{q}\cdot\mb{r}’}\,dV’
=-\frac{m}{2\pi\hbar^{2}}\times c
\tag{22}
\end{equation*}

他方, 摂動項 (5) に相当する場合の散乱振幅 \(f_2\) を考えると, 前述の問題 (3) から
\begin{equation*}
U=H^{(5)}=\frac{1}{4mc}\left(2\mu+\frac{e\hbar}{2mc}\right)\hbar\nabla\cdot\mb{E}
\end{equation*}

であるが, 中性子の場合では \(e=0\) であるからカッコ中の第2項目はゼロとして
\begin{equation*}
U=\frac{1}{4mc}(2\mu)\times \hbar\nabla\cdot\mb{E}
=\frac{1}{4m c}\left(2\kappa_n\frac{e\hbar}{2m c}\right)\times \hbar\nabla\cdot\mb{E}
\end{equation*}

である.問題の前述で示したごとく \(\nabla\cdot\mb{E}=4\pi\rho=4\pi Z|e|\delta(\mb{r})\) であったから,
\begin{equation*}
U=\frac{1}{4m c}\left(2\kappa_n\frac{e\hbar}{2mc}\right)\times
\hbar\times 4\pi\, e\, \delta(\mb{r})
=4\pi\kappa_n \times \frac{e^{2}\hbar^{2}}{4m^{2}c^{2}}\delta(\mb{r})
\end{equation*}

従って, この場合の散乱振幅 \(f_2\) は
\begin{equation*}
f_2 = -\frac{m}{2\pi\hbar^{2}}\int 4\pi\kappa_n\,\frac{e^{2}\hbar^{2}}{4m^{2}c^{2}}
\delta(\mb{r})\,e^{i\mb{q}\cdot\mb{r}’}\,dV’
=-\frac{m}{2\pi\hbar^{2}}\times 4\pi\kappa_n\,\frac{e^{2}\hbar^{2}}{4m^{2}c^{2}}
\tag{23}
\end{equation*}

式 (22) の \(f_1\) と式 (23) の \(f_2\) を等しいとして次を得る:
\begin{equation*}
c=4\pi\kappa_n \frac{e^{2}\hbar^{2}}{4m^{2}c^{2}}
=4\pi\kappa_n \left(\frac{e\hbar}{2m_pc}\right)^{2}=4\pi\kappa_n \mu_{N}^{2}
\tag{24}
\end{equation*}

平均ポテンシャル \(\overline{V}\) を求めてみる.題意から, 平均ポテンシャル \(\overline{V}\) は次式の関係にある:
\begin{equation*}
\int c\delta(\mb{r})\,dV =\frac{4\pi}{3}\left(\frac{e^{2}}{m_e c^{2}}\right)^{3}\overline{V}
\quad\rightarrow\quad c=\frac{4\pi}{3}\left(\frac{e^{2}}{m_e c^{2}}\right)^{3}\overline{V}
\end{equation*}

従って, 式 (24) を用いると
\begin{equation*}
4\pi\kappa_n \mu_{N}^{2}=\frac{4\pi}{3}\left(\frac{e^{2}}{m_e c^{2}}\right)^{3}\overline{V}
=\frac{4\pi}{3}\,r_{e}^{3}\,\overline{V}
\tag{25}
\end{equation*}

上式に \(\kappa_n=-1.9135\), 核磁子 \(\mu_N=5.05\times10^{-24}\,\mathrm{erg/gauss}\), 古典電子半径 \(r_e=2.818\times10^{-13}\,\mathrm{cm}\) を代入して \(\overline{V}\) を求めると, 次を得る:
\begin{align*}
4\pi&\times(-1.9135)\times\bigl(5.05\times 10^{-24}\bigr)^{2}
=\frac{4\pi}{3}\times \bigl(2.818\times 10^{-13}\bigr)^{3}\,\overline{V},\\
\rightarrow\quad
\overline{V}&=\frac{3\times(-1.9135)\times(5.05\times10^{-24})^{2}}{(2.818\times10^{-13})^{3}}
=-6.542\times 10^{-9}\,\mathrm{erg}
=-\frac{6.542\times 10^{-9}}{1.602\times 10^{-12}}\,\mathrm{eV}\\
&=-4084\,\mathrm{eV}
\end{align*}


〈 (5)の解答例 〉この問題は要求どうりの回答が出来なかった.式の定義がしっかり書かれていないこともあるので, 取り敢えずは, 回答とは異なるが自分なりに考えた式を示しておくことにする.後で検討を要する.

問題文の式の左辺中の被積分関数 \(\psi^{*}\mb{\alpha}f(\mb{r})\psi\) を具体的に書くと,

\begin{equation*}
\psi^{*}\mb{\alpha}f(\mb{r})\psi=\begin{pmatrix} \phi^{*} & \chi^{*} \end{pmatrix}
\begin{pmatrix} 0 & \mb{\sigma} \\ \mb{\sigma} & 0 \end{pmatrix}
\begin{pmatrix} f\phi \\ f\chi \end{pmatrix}
=\phi^{*}\mb{\sigma}f\chi + \chi^{*}\mb{\sigma}f\phi
\tag{26}
\end{equation*}

他方, § 3.7 の式 (11-6) に於いて分母を \(2mc^{2}\) だけ残すことで, 次の1次近似式を得る:
\begin{align*}
&\psi_b = \frac{c\mb{\sigma}\cdot\mb{\pi}}{2mc^{2}+W-V}\psi_a \quad\rightarrow\quad
\chi=\frac{\mb{\sigma}\cdot\mb{\pi}}{2mc}\phi,\\
&\Rightarrow\quad
\chi=\frac{1}{2mc}\mb{\sigma}\cdot\left(\mb{p}-\frac{e}{c}\mb{A}\right)\phi,\quad
\chi^{*}=\frac{1}{2mc}\left(\mb{p}-\frac{e}{c}\mb{A}\right)^{*}\phi^{*}\cdot\mb{\sigma}
\tag{27}
\end{align*}
これを式 (26) の各項に代入すると,
\begin{align*}
\phi^{*}\mb{\sigma}f\chi&=\phi^{*}\mb{\sigma}f\left(\frac{\mb{\sigma}\cdot\mb{\pi}}{2mc}\phi\right)
=\frac{\phi^{*}}{2mc}\mb{\sigma}f(\mb{\sigma}\cdot\mb{\pi})\,\phi,\\
\chi^{*}\mb{\sigma}f\phi&=\phi^{*}\frac{\mb{\pi}^{*}\cdot\mb{\sigma}}{2mc}\mb{\sigma}f\phi
=\frac{\phi^{*}}{2mc}(\mb{\sigma}\cdot\mb{\pi}^{*})\,\mb{\sigma}f\phi
\tag{28}
\end{align*}
ここで, 公式 (F2) から得られる次の公式を利用する:
\begin{align*}
&\mb{\sigma}(\mb{\sigma}\cdot\mb{a})=\mb{a}+i\mb{a}\times\mb{\sigma},\quad\rightarrow
\quad f\mb{\sigma} (\mb{\sigma}\cdot\mb{a})=\mb{\sigma}f (\mb{\sigma}\cdot\mb{a})=f\mb{a}+i f\mb{a}\times\mb{\sigma},\\
&(\mb{\sigma}\cdot\mb{a})\mb{\sigma}=\mb{a}+i\mb{\sigma}\times\mb{a},\quad\rightarrow
\quad (\mb{\sigma}\cdot\mb{a})\mb{\sigma}f=(\mb{\sigma}\cdot\mb{a}) f \mb{\sigma}=\mb{a}f +i\mb{\sigma}\times\mb{a}f
\tag{29}
\end{align*}

公式 (29) を式 (28) に用い, また問題を簡略化して \(\mb{A}=0\) の場合を考える.\(\mb{p}=-i\hbar\nabla\) とし, また, 運動量演算子は実(real)であることを用いて \(\mb{p}^{*}=\mb{p}\) とするならば,
\begin{align*}
\phi^{*}\mb{\sigma}f\chi&=\frac{\phi^{*}}{2mc}\mb{\sigma}f(\mb{\sigma}\cdot\mb{\pi})\,\phi
=\frac{\phi^{*}}{2mc}\left(f\mb{\pi}+if\mb{\pi}\times\mb{\sigma}\right)\phi
=\phi^{*}\left(\frac{f\mb{p}}{2mc} +\frac{if\mb{p}\times\mb{\sigma}}{2mc}\right)\phi\\
&=\phi^{*}\left(\frac{f\mb{p}}{2mc} +\frac{\hbar}{2mc}f\nabla\times\mb{\sigma}\right)\phi,\\
\chi^{*}\mb{\sigma}f\phi&=\frac{\phi^{*}}{2mc}(\mb{\sigma}\cdot\mb{\pi}^{*})\,\mb{\sigma}f\phi
=\frac{\phi^{*}}{2mc}\left(\mb{\pi}^{*}f+i\mb{\sigma}\times\mb{\pi}^{*}f\right)\phi
=\phi^{*}\left(\frac{\mb{p}^{*}f}{2mc} +\frac{i\mb{\sigma}\times\mb{p}^{*}f}{2mc}\right)\phi\\
&=\phi^{*}\left(\frac{\mb{p}f}{2mc}-\frac{\hbar}{2mc}\mb{\sigma}\times\nabla f\right)\phi
=\phi^{*}\left(\frac{\mb{p}f}{2mc}-\frac{\hbar}{2mc}\mb{\sigma}\times\nabla f-\frac{\hbar}{2mc}f\mb{\sigma}\times\nabla\right)\phi,\\
\therefore\quad \phi^{*}\mb{\sigma}f\chi &+ \chi^{*}\mb{\sigma}f\phi
=\phi^{*}\left(\frac{f\mb{p}+\mb{p}f}{2mc}
-\frac{\hbar}{2mc}\mb{\sigma}\times\nabla f\right)\phi
\tag{30}
\end{align*}
従って,
\begin{align*}
\int \psi^{*}\mb{\alpha}f(\mb{r})\psi\,dV
=\int \bigl(\phi^{*}\mb{\sigma}f\chi + \chi^{*}\mb{\sigma}f\phi\bigr)\,dV
=\int \phi^{*}\left(\frac{f\mb{p}+\mb{p}f}{2mc}-\frac{\hbar}{2mc}\mb{\sigma}\times\nabla f\right)\phi\,dV
\end{align*}

References

References
1 [ブログ註] オッペンハイマー:「電気力学」§ 11 によれば,「一様な速度 \(\mb{v}\) で動いている点電荷 \(e\) の作る電場 \(\mb{E}\) と磁場 \(\mb{H}\)」とは, 次の関係にある:
\begin{equation*}
\mb{H}=\frac{\mb{v}}{c}\times\mb{E}
\tag{11.7}
\end{equation*}
2 Schiff, “Quantum Mechanics,” McGraw-Hill, New York, 1949, pp. 323ff.
3 [ブログ註] (A.P.アーヤ:基礎現代物理学 § 11.2 から抜粋) 原子スペクトルでは, 原子の磁気モーメントは「ボーア磁子」\(\mu_B\) を単位として測られる.ボーア磁子とは, 角運動量 \(1\hbar\) で軌道運動する電子に付随する磁気モーメントのことである.\(\mu_B\) の大きさは,
\[\mu_B =\frac{e\hbar}{2m_e c}=0.9273\times 10^{-20}\ \mathrm{erg}/\mathrm{gauss}\]
である.ただし \(e\) は電子の電荷そして \(m_e\) は電子質量である.
全ての原子核は電荷を持つから, 原子核の全(スピン)角運動量がゼロでないなら, 原子核内部に於ける核子の運動は(原子中の電子の運動と同様に)原子核磁気モーメントを引き起こすはずである.また, 原子核内電荷分布が球対称ならば, 原子核は磁気双極子モーメントのみを生じる.核磁気モーメント \(\mu_I\) は「核磁子」の単位で測られる.「核磁子 \(\mu_N\)」は次の様に定義されている:
\[\mu_N =\frac{e\hbar}{2m_p c}=\frac{\mu_B}{1836}=5.050\times10^{-24}\ \mathrm{erg}/\mathrm{gauss}\]
ただし \(m_p\) は陽子質量であり電子質量の 1836 倍である.核磁気モーメント \(\mu_I\) の実測値は \(-3\mu_N\) から \(+10\mu_N\) の間にある.正符号は核の磁気モーメントが核スピンと同じ方向であり, 負の符号は逆向きであることを示す.特に注目すべきは, 陽子と中性子の磁気モーメントの実測値である.陽子の磁気モーメントは,
\[\mu_{\mathrm{proton}}= + 2.79276\, \mu_N\]
であり, 予想されるような 1核磁子の値とは異なる.このことを「陽子は異常磁気モーメントを持つ」と言い,これは陽子の電荷分布が一様でないことを意味している.同様に中性子は総電荷はゼロであるにもかかわらず磁気モーメントを持ち,その値は
\[\mu_{\mathrm{neutron}}= – 1.191315\, \mu_N\]
である.この事実は「中性子は一様でない電荷分布を持ち, 磁気モーメントの方向が固有のスピン角運動量ベクトルの方向と逆向きである」ことを示している.
4 L.Foldy, Phys. Rev., \(\mb{87}\), 693 (1952).
5 [ブログ註] ランダウ=リフシッツ:「量子力学」§ 31 クーロン場の中の運動 を参照すべし.