電子は電荷量のほかにも、角速度の次元を持ったパラメータを持っている。このパラメータのことをスピンとよぶ。このスピンを表現するために、行列式を使った固有方程式が使われる。その固有方程式中の行列のことを、スピン行列とよぶ。
この記事では、電子などの素粒子のスピン行列を導出する。
スピン量子数とスピン磁気量子数の関係
前の記事では、軌道角運動量の固有方程式における方位量子数と磁気量子数の関係が、スピン演算子の固有方程式におけるスピン量子数とスピン磁気量子数に対応していることを示した。スピン量子数は半整数になり得ることに注意すると、スピン量子数とスピン磁気量子数の関係は次のように書ける。
$$m_s=-s, -s+1,…,s-1,s$$
参考:軌道角運動量とスピン角運動量の違い
例えば、電子のスピン量子数は、実験値より
$$s=\frac{1}{2}$$
となっている。したがって、電子のスピン磁気量子数は
$$m_s= \pm \frac{1}{2}$$
の2通り存在する。\(m_s=+1/2\)は上向きスピン、\(m_s=-1/2\)は下向きスピンに対応する。
電子のスピン行列の導出
上向きスピンと下向きスピンのベクトル
ここでは、上向きスピンに対応する固有ベクトル\( | \uparrow >\)を
$$| \uparrow > \equiv \left( \begin{array}{c} 1 \\ 0 \end{array} \right)$$
下向きのスピンに対応する固有ベクトル\( | \downarrow >\)を
$$| \downarrow > \equiv \left( \begin{array}{c} 0 \\ 1 \end{array} \right)$$
と定義する。パウリのスピン行列とは、この\( | \uparrow >\)や\( | \downarrow >\)を固有ベクトルとする固有方程式に含まれる行列のことである。
参考:ブラベクトル・ケットベクトルの意味とは
スピンの昇降演算子がスピンのベクトルに与える影響
そもそも昇降演算子とは、固有ベクトルに作用することで、固有方程式の固有値を増減させるような演算子のことであった。事実、ケットベクトル\(|l,m \! >\)に軌道角運動量の昇降演算子\(\hat{L}_{\pm}\)を作用させたら、ケットベクトル中の磁気量子数\(m\)が増減した。したがって、もし昇降演算子をスピンを表すベクトル\( | \uparrow >,\)\( | \downarrow >\)に作用させたら、スピン磁気量子数が増減するだろうと予想できる。
そしてスピン磁気量子数が増減するということは、\( | \uparrow >\)から\( | \downarrow >\)への遷移やその逆の遷移が起こるということになる。\( | \uparrow >\)と\( | \downarrow >\)が対応するスピン磁気量子数はそれぞれ\(+1/2\)と\(-1/2\)だから、\( | \uparrow >\)に下降演算子\(\hat{s}_-\)を適応させれば\( | \downarrow >\)、\( | \downarrow >\)に上昇演算子\(\hat{s}_+\)を適応させれば\( | \uparrow >\)が得られると考えられる。
軌道角運動量の昇降演算子は
$$\hat{L}_{\pm}|l,m \! >=ħ \sqrt{l(l+1)-m(m \pm 1)}|l,m \pm 1 \! >$$
を満たす。そしてこの式のスピンバージョンは次のようになる。\(l\)はスピン量子数\(s\)、\(m\)はスピン磁気量子数\(m_s\)に対応していることに気を付ければ、簡単に次の式が導けるはずだ。
$$\hat{s}_+ | \downarrow >=ħ | \uparrow >$$
$$\hat{s}_- | \uparrow >=ħ | \downarrow >$$
一応前者の式の解説を載せる。この式では上昇演算子\(\hat{s}_+\)を使っているため、途中の\(m_s(m_s \pm 1)\)の\(\pm\)の符号は正になる。元の\(\hat{L}_{\pm}\)に関する固有方程式に注目すれば、\(m_s\)に代入すべきは\(+1/2\)でなく\(-1/2\)であることもわかるだろう。後者の式も同様に考えれば求まる。
\begin{eqnarray}\hat{s}_+ | \downarrow >&=&ħ \sqrt{s(s+1)-m_s(m_s+ 1)}| \uparrow >\\&=&ħ \sqrt{\frac{1}{2} \left( \frac{1}{2}+1 \right) -\left( -\frac{1}{2} \right) \left( -\frac{1}{2}+ 1 \right) }| \uparrow >\\&=&ħ | \uparrow >\end{eqnarray}
また、電子のスピンは2種類しか存在しないため、\( | \uparrow >\)に\(\hat{s}_+\)を適応させたり、\( | \downarrow >\)に\(\hat{s}_-\)を適応させたりすることには意味はない。そもそもその操作によって得られる状態が物理的に存在しないということである。したがって、次の式を得る。
$$\hat{s}_+ | \uparrow >=0$$
$$\hat{s}_- | \downarrow >=0$$
昇降演算子の具体形
スピンの昇降演算子\(\hat{s}_{\pm}\)は、具体的には次のような形になる。実際に上の式に代入して確かめてみよう。
$$\hat{s}_+ =\begin{pmatrix} 0 & ħ \\ 0 & 0 \end{pmatrix}$$
$$\hat{s}_- =\begin{pmatrix} 0 & 0 \\ ħ & 0 \end{pmatrix}$$
x,y軸方向のスピン演算子
軌道角運動量のときと同じように、昇降演算子を定義する。
$$\hat{s}_+ \equiv \hat{s}_x+i\hat{s}_y$$
$$\hat{s}_- \equiv \hat{s}_x-i\hat{s}_y$$
これを\(\hat{s}_x\)と\(\hat{s}_y\)に関する連立方程式として解くと、次の関係を得る。
$$\hat{s}_x=\frac{1}{2}(\hat{s}_+ + \hat{s}_- )$$
$$\hat{s}_y=\frac{1}{2i}(\hat{s}_+ – \hat{s}_- )$$
あとはこの式の昇降演算子\(\hat{s}_{\pm}\)に、具体形を代入すればよい。
\begin{eqnarray}\hat{s}_x&=&\frac{1}{2}(\hat{s}_+ + \hat{s}_- )\\&=&\frac{1}{2} \left[ \begin{pmatrix} 0 & ħ \\ 0 & 0 \end{pmatrix} + \begin{pmatrix} 0 & 0 \\ ħ & 0 \end{pmatrix} \right]\\&=&\frac{ħ}{2} \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} \end{eqnarray}
\begin{eqnarray}\hat{s}_y&=&\frac{1}{2i}(\hat{s}_+ – \hat{s}_- )\\&=&\frac{1}{2i} \left[ \begin{pmatrix} 0 & ħ \\ 0 & 0 \end{pmatrix} – \begin{pmatrix} 0 & 0 \\ ħ & 0 \end{pmatrix} \right]\\&=&\frac{ħ}{2} \begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix} \end{eqnarray}
z軸方向のスピン演算子
軌道角運動量演算子\(\hat{L}_z\)では次の固有方程式
$$\hat{L}_z|l,m \! >=mħ|l,m \!>$$
が成り立ったが、これのスピン演算子バージョンが
$$\hat{s}_z| \uparrow >=m_s ħ| \uparrow >$$
$$\hat{s}_z| \downarrow >=m_s ħ| \downarrow >$$
である。以降、これを満たすスピン演算子\(\hat{s}_z\)を求める。
上下スピンに対応する固有ベクトルとスピン演算子の行列
電子のスピン量子数\(s\)は
$$s=\frac{1}{2}$$
だから、電子のスピン磁気量子数\(m_s\)は
$$m_s= \pm \frac{1}{2}$$
の2通りとなる。\(+1/2\)が上向きスピン、\(-1/2\)が下向きスピンに対応する。
上の\( | \uparrow >\)と\( | \downarrow >\)の定義を考慮すると、z軸のスピン演算子を表す行列\(\hat{s}_z\)は次のように定義すれば都合が良い。
$$\hat{s}_z=\frac{ħ}{2} \begin{pmatrix} 1 & 0\\ 0 & -1 \end{pmatrix}$$
事実、上向きスピンに関しては
\begin{eqnarray}\hat{s}_z | \uparrow > &=&\frac{ħ}{2}\begin{pmatrix} 1 & 0\\ 0 & -1 \end{pmatrix} \left( \begin{array}{c} 1 \\ 0 \end{array} \right)\\&=&\frac{ħ}{2} \left( \begin{array}{c} 1 \\ 0 \end{array} \right) \\&=&\frac{ħ}{2} | \uparrow >\end{eqnarray}
となるし、下向きスピンに関しては
\begin{eqnarray}\hat{s}_z | \downarrow > &=&\frac{ħ}{2}\begin{pmatrix} 1 & 0\\ 0 & -1 \end{pmatrix} \left( \begin{array}{c} 0 \\ 1 \end{array} \right)\\&=&-\frac{ħ}{2} \left( \begin{array}{c} 0 \\ 1 \end{array} \right) \\&=&-\frac{ħ}{2} | \downarrow >\end{eqnarray}
となるため、それぞれのスピンに関して固有方程式
$$\hat{s}_z| \uparrow >=m_s ħ | \uparrow >$$
$$\hat{s}_z| \downarrow >=m_s ħ | \downarrow >$$
が成立しているといえる。
以上で、z軸方向のスピン演算子が求められた。
スピン演算子とパウリ行列
スピン演算子\(\hat{\bf s}\)は、x,y,z軸のそれぞれの成分のスピン演算子\(\hat{s}_x,\)\(\hat{s}_y,\)\(\hat{s}_z\)をまとめたものである。
$$\hat{\bf s}=(\hat{s}_x,\hat{s}_y,\hat{s}_z)$$
そして、今まで求めてきたように、\(\hat{s}_x,\)\(\hat{s}_y,\)\(\hat{s}_z\)はそれぞれ次のように定義される。
$$\hat{s}_x \equiv \frac{ħ}{2} \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}$$
$$\hat{s}_y \equiv \frac{ħ}{2} \begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}$$
$$\hat{s}_z \equiv \frac{ħ}{2} \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}$$
本によってはディラック定数\(ħ\)が抜けているものもあるため、その場合はその本や教授に従うこと。
そしてパウリのスピン行列は、スピン演算子の行列部分のことを指す。
$$σ_x=\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}$$
$$σ_y=\begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}$$
$$σ_z=\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}$$
したがって、スピン行列の記号を使ってスピン演算子を記述すると、
$$\hat{s}_x \equiv \frac{ħ}{2} σ_x$$
$$\hat{s}_y \equiv \frac{ħ}{2} σ_y$$
$$\hat{s}_z \equiv \frac{ħ}{2} σ_z$$
以上でパウリのスピン行列が求められた。
より一般的なスピン行列の求め方
ここまでは電子のような、スピン量子数が1/2のスピン演算子の求め方を紹介した。だが、もちろんすべての素粒子のスピンが1/2というわけではない。例えば光子のスピンは1であり、グラビトン(重力子)のスピンは2であるとされている。そのような粒子のスピン行列を求めようにも、上の方法では不可能である。
ただし、次の方程式を利用すれば、1/2以外のスピン量子数を持つ素粒子のスピン演算子を求められる。
$$< \! s,m_s|\hat{s}_z|s,m’_s \! >=mħ \delta _{m_s,m’_s}$$
$$< \! s,m_s|\hat{s}_{\pm}|s,m_s’ \! >=ħ \sqrt{s(s+1)-m_s(m_s \pm 1)} \delta _{m_s,m_s \mp 1}$$
$$< \! s,m_s|\hat{\bf s}^2|s,m_s’ \! >=s(s+1)ħ^2 \delta _{m_s,m’_s}$$
参考:角運動量の昇降演算子の交換関係と物理的意味について
これらの導出についてだが、ここでは一番上の式のみ紹介する。他の式も同様に求められる。単純に、元々の方程式の両辺の左側からブラベクトル\(< \! s,m’_s|\)をつければ良い。
\begin{eqnarray} < \! s,m’_s| \hat{s}_z|s,m_s \! >&=&< \! s,m’_s|mħ|s,m_s \!>\\&=&mħ< \! s,m’_s|s,m_s \!>\\&=&mħ \delta _{m_s,m’_s}\end{eqnarray}
演算子とその行列表記について
演算子\(\hat{A}\)の行列成分\(A_{ij}\)は次のように書けた。
$$A_{ij}=< \! i|{\bf A}|j \! >$$
ここからスピン行列を導出するには、上の固有方程式中のスピン磁気量子数\(m_s\)と\(m’_s\)が演算子の行列の行と列に対応していると考える。つまり、行列の行と列の数は、スピン磁気量子数の数と一致するといえる。このことから、求める行列は正方行列だとわかる。
x,y軸成分のスピン行列
1/2スピンのときと同様に、昇降演算子を経由してスピン行列を求める。
上昇演算子\(\hat{s}_+\)について、
$$< \! s,m_s|\hat{s}_+ |s,m_s’ \! >=ħ \sqrt{s(s+1)-m_s(m_s + 1)} \delta _{m_s,m_s + 1}$$
が成立する。まず、クロネッカーのデルタ\(\delta _{m_s,m’_s}\)より、列番号が行番号+1になるような要素のみ0以外の数字になり得ることがわかる。\(s\)と\(m_s\)の関係を考えると、次の関係が予想できる。
\begin{eqnarray}\hat{s}_+ &=&\begin{pmatrix} < \! s,s|\hat{s}_+ |s,s \! > & < \! s,s|\hat{s}_+ |s,s-1 \! > & \ldots & < \! s,s|\hat{s}_+ |s, -s \! > \\ < \! s,s-1|\hat{s}_+ |s,s \! > & < \! s,s-1|\hat{s}_+ |s,s-1 \! > & \ldots & < \! s,s-1|\hat{s}_+ |s, -s \! > \\ \vdots & \vdots & \ddots & \vdots \\ < \! s,-s|\hat{s}_+ |s,s \! > & < \! s,-s|\hat{s}_+ |s,s-1 \! > & \ldots & < \! s,-s|\hat{s}_+ |s, -s \! > \end{pmatrix}\\&=&\begin{pmatrix} 0 & < \! s,s|\hat{s}_+ |s,s-1 \! > & 0 & \ldots & 0 \\ 0 & 0 & < \! s,s-1|\hat{s}_+ |s,s-2 \! > & \ldots & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & 0 & \ldots & < \! s,-s+1|\hat{s}_+ |s,-s \! > \\ 0 & 0 & 0 & \ldots & 0 \end{pmatrix} \end{eqnarray}
同様にして、下降演算子\(\hat{s}_-\)も次のようになる。こちらは、行番号が列番号+1になるような要素のみ0以外の数字になり得る。
\begin{eqnarray}\hat{s}_- &=&\begin{pmatrix} < \! s,s|\hat{s}_- |s,s \! > & < \! s,s|\hat{s}_- |s,s-1 \! > & \ldots & < \! s,s|\hat{s}_- |s, -s \! > \\ < \! s,s-1|\hat{s}_- |s,s \! > & < \! s,s-1|\hat{s}_- |s,s-1 \! > & \ldots & < \! s,s-1|\hat{s}_- |s, -s \! > \\ \vdots & \vdots & \ddots & \vdots \\ < \! s,-s|\hat{s}_- |s,s \! > & < \! s,-s|\hat{s}_- |s,s-1 \! > & \ldots & < \! s,-s|\hat{s}_- |s, -s \! > \end{pmatrix}\\&=&\begin{pmatrix} 0 & 0 & \ldots & 0 & 0 \\ < \! s,s-1|\hat{s}_- |s,s \! > & 0 & \ldots & 0 & 0 \\ 0 & < \! s,s-2|\hat{s}_- |s,s-1 \! > & \ldots & 0 & 0 \\ \vdots & \vdots & \ddots & \vdots & \vdots \\ 0 & 0 & \ldots & < \! s,-s|\hat{s}_- |s, -s+1 \! > & 0 \end{pmatrix}\end{eqnarray}
例えばスピンが\(1\)の場合、上の行列中のsに1を代入すれば良い。
$$\hat{s}_+ =\begin{pmatrix} 0 & < \! 1,1|\hat{s}_+ |1,0 \! > & 0 \\ 0 & 0 & < \! 1,0|\hat{s}_- |1, -1 \! > \\ 0 & 0 & 0 \end{pmatrix}$$
$$\hat{s}_- =\begin{pmatrix} 0 & 0 & 0 \\ < \! 1,0|\hat{s}_- |1,1 \! > & 0 & 0 \\ 0 & < \! 1,-1|\hat{s}_- |1,0 \! > & 0 \end{pmatrix}$$
このように、s=1を代入すると、求める行列は3×3の正方行列となることが予想できる。したがって、この場合のスピンを表すベクトルは次の3種類が予想される。
$$\left( \begin{array}{c} 1 \\ 0 \\ 0 \end{array} \right) , \ \left( \begin{array}{c} 0 \\ 1 \\ 0 \end{array} \right) , \ \left( \begin{array}{c} 0 \\ 0 \\ 1 \end{array} \right)$$
後は、
$$\hat{s}_+ \left( \begin{array}{c} 0 \\ 0 \\ 1 \end{array} \right)=ħ \sqrt{1(1+1)-(-1)((-1) + 1)}\left( \begin{array}{c} 0 \\ 1 \\ 0 \end{array} \right)$$
$$\hat{s}_+ \left( \begin{array}{c} 0 \\ 1 \\ 0 \end{array} \right)=ħ \sqrt{1(1+1)-0(0 + 1)}\left( \begin{array}{c} 1 \\ 0 \\ 0 \end{array} \right)$$
$$\hat{s}_+ \left( \begin{array}{c} 1 \\ 0 \\ 0 \end{array} \right)=0$$
$$\hat{s}_- \left( \begin{array}{c} 1 \\ 0 \\ 0 \end{array} \right)=ħ \sqrt{1(1+1)-1(1 – 1)}\left( \begin{array}{c} 0 \\ 1 \\ 0 \end{array} \right)$$
$$\hat{s}_- \left( \begin{array}{c} 0 \\ 1 \\ 0 \end{array} \right)=ħ \sqrt{1(1+1)-0(0 – 1)}\left( \begin{array}{c} 0 \\ 0 \\ 1 \end{array} \right)$$
$$\hat{s}_- \left( \begin{array}{c} 0 \\ 0 \\ 1 \end{array} \right)=0$$
を満たすような昇降演算子を求めれば良い。
$$\hat{s}_+=ħ \begin{pmatrix} 0 & \sqrt{2} & 0 \\ 0 & 0 & \sqrt{2} \\ 0 & 0 & 0 \end{pmatrix}$$
$$\hat{s}_- =ħ \begin{pmatrix} 0 & 0 & 0 \\ \sqrt{2} & 0 & 0 \\ 0 & \sqrt{2} & 0 \end{pmatrix}$$
後は、昇降演算子とx,y成分のスピン演算子の関係
$$\hat{s}_x=\frac{1}{2}(\hat{s}_+ + \hat{s}_- )$$
$$\hat{s}_y=\frac{1}{2i}(\hat{s}_+ – \hat{s}_- )$$
を利用して、x,y成分のスピン行列を求めれば良い。
\begin{eqnarray}\hat{s}_x&=&\frac{1}{2}(\hat{s}_+ + \hat{s}_- )\\&=&\frac{ħ}{2} \left[ \begin{pmatrix} 0 & \sqrt{2} & 0 \\ 0 & 0 & \sqrt{2} \\ 0 & 0 & 0 \end{pmatrix} + \begin{pmatrix} 0 & 0 & 0 \\ \sqrt{2} & 0 & 0 \\ 0 & \sqrt{2} & 0 \end{pmatrix} \right]\\&=&\frac{ħ}{2} \begin{pmatrix} 0 & \sqrt{2} & 0 \\ \sqrt{2} & 0 & \sqrt{2} \\ 0 & \sqrt{2} & 0 \end{pmatrix}\\&=&\frac{ħ}{\sqrt{2}} \begin{pmatrix} 0 & 1 & 0 \\ 1 & 0 & 1 \\ 0 & 1 & 0 \end{pmatrix}\end{eqnarray}
\begin{eqnarray}\hat{s}_y&=&\frac{1}{2i}(\hat{s}_+ – \hat{s}_- )\\&=&\frac{ħ}{2i} \left[ \begin{pmatrix} 0 & \sqrt{2} & 0 \\ 0 & 0 & \sqrt{2} \\ 0 & 0 & 0 \end{pmatrix} – \begin{pmatrix} 0 & 0 & 0 \\ \sqrt{2} & 0 & 0 \\ 0 & \sqrt{2} & 0 \end{pmatrix} \right]\\&=&\frac{ħ}{\sqrt{2}i} \begin{pmatrix} 0 & 1 & 0 \\ -1 & 0 & 1 \\ 0 & -1 & 0 \end{pmatrix}\\&=&\frac{ħ}{\sqrt{2}} \begin{pmatrix} 0 & -i & 0 \\ i & 0 & -i \\ 0 & i & 0 \end{pmatrix}\end{eqnarray}
z軸成分のスピン行列
$$< \! s,m_s|\hat{s}_z|s,m’_s \! >=mħ \delta _{m_s,m’_s}$$
に注目する。まず、クロネッカーのデルタ\(\delta _{m_s,m’_s}\)より、求めたい行列の行番号と列番号が同じでない要素はすべて0になることがわかる。
前と同じように、\(m_s\)と\(m’_s\)がスピン行列の行と列に対応していると考える。スピン量子数\(s\)とスピン磁気量子数\(m_s\)の関係にも注意すると、求めるスピン行列\(\hat{s}_z\)は次のように書ける。
\begin{eqnarray}\hat{s}_z&=&\begin{pmatrix} < \! s,s|\hat{s}_z|s,s \! > & < \! s,s|\hat{s}_z|s,s-1 \! > & \ldots & < \! s,s|\hat{s}_z|s, -s \! > \\ < \! s,s-1|\hat{s}_z|s,s \! > & < \! s,s-1|\hat{s}_z|s,s-1 \! > & \ldots & < \! s,s-1|\hat{s}_z|s, -s \! > \\ \vdots & \vdots & \ddots & \vdots \\ < \! s,-s|\hat{s}_z|s,s \! > & < \! s,-s|\hat{s}_z|s,s-1 \! > & \ldots & < \! s,-s|\hat{s}_z|s, -s \! > \end{pmatrix}\\&=&\begin{pmatrix} < \! s,s|\hat{s}_z|s,s \! > & 0 & \ldots & 0 \\ 0 & < \! s,s-1|\hat{s}_z|s,s-1 \! > & \ldots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \ldots & < \! s,-s|\hat{s}_z|s, -s \! > \end{pmatrix}\end{eqnarray}
例えばスピンが\(1\)の場合、上の行列中の\(s\)に1を代入すれば良い。
$$\hat{s}_z=\begin{pmatrix} < \! 1,1|\hat{s}_z|1,1 \! > & 0 & 0 \\ 0 & < \! 1,0|\hat{s}_z|1,0 \! > & 0 \\ 0 & 0 & < \! 1,-1|\hat{s}_z|1, -1 \! > \end{pmatrix}$$
このように、s=1を代入すると、求める行列は3×3の正方行列となることが予想できる。したがって、この場合のスピンを表すベクトルは次の3種類が予想される。
$$\left( \begin{array}{c} 1 \\ 0 \\ 0 \end{array} \right) , \ \left( \begin{array}{c} 0 \\ 1 \\ 0 \end{array} \right) , \ \left( \begin{array}{c} 0 \\ 0 \\ 1 \end{array} \right)$$
それぞれの固有ベクトルが対応するスピン磁気量子数を、順に
$$m_s=1, \ 0, -1$$
とする。後は、
$$\hat{s}_z| s,m_s \! >=m_s ħ | s,m_s \! >$$
を満たすように行列要素を求めていけば良い。その結果、次の行列を得られる。
$$\hat{s}_z= ħ \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & -1 \end{pmatrix}$$
まとめ
・スピンとは、角運動量の次元を持つ素粒子のパラメータのことである。
・スピン行列の求め方を紹介した。
参考文献
原田勲･杉山忠男(2009)『講談社基礎物理学シリーズ6 量子力学I』,講談社.
二宮正夫･杉野文彦･杉山忠男(2010)『講談社基礎物理学シリーズ7 量子力学II』,講談社.
村上雅人(2008)『なるほど量子力学III』,海鳴社.
