スピノールの時間発展(パウリ方程式のスピン依存)

磁場中の電子に対するスピン自由度を含めた波動関数 $\Psi(\boldsymbol{r}, t)$ は、シュレーディンガー方程式から導かれるパウリの方程式と呼ばれる

\begin{align}
i\hbar \frac{\partial \Psi(\boldsymbol{r}, t)}{\partial t} = \left[ \frac{1}{2m_e} (-i\hbar\nabla + e\boldsymbol{A})^2 + V(r) + \frac{e}{m_e} \hat{\boldsymbol{S}} \cdot \boldsymbol{B} \right] \Psi(\boldsymbol{r}, t)
\end{align}

を満たすね。この方程式を満たす波動関数は空間依存部分 $\psi(\boldsymbol{r}, t)$ とスピン依存部分 $\chi(t)$ に分離して、$\Psi(\boldsymbol{r}, t)=\psi(\boldsymbol{r}, t)\chi(t)$ と表わすことで、方程式は変数分離できるね。

\begin{align}
i\hbar \frac{\partial \psi(\boldsymbol{r}, t)}{\partial t} &\, = \left[ \frac{1}{2m_e} (-i\hbar\nabla + e\boldsymbol{A})^2 + V(r) \right] \psi(\boldsymbol{r}, t)\\
i\hbar \frac{d \chi( t )}{d t} &\, = \frac{e}{m_e}\, \hat{\boldsymbol{S}} \cdot \boldsymbol{B}\, \chi(\boldsymbol{r}, t)
\end{align}

このスピン依存部分は、スピン角運動量 $\boldsymbol{S} = (S_x, S_y, S_z)$ をパウリ行列

\begin{align}
S_x = \frac{\hbar}{2} \left( \matrix{ 0 & 1 \cr 1 & 0} \right) \ , \ S_y = \frac{\hbar}{2} \left( \matrix{ 0 & -i \cr 1 & 0} \right)\ , \ S_z = \frac{\hbar}{2} \left( \matrix{ 1 & 0 \cr 0 & -1} \right)
\end{align}

で表したときのスピノールと呼ばれる縦ベクトルで表すことができるね。スピノールを

\begin{align}
\chi( t ) = \left( \matrix{\chi_{\uparrow}(t) \cr \chi_{\downarrow}(t)} \right)
\end{align}

と表しておいて、パウリ方程式に代入して、$\chi_{\uparrow}(t)$ と $\chi_{\downarrow}(t)$ の時間依存性を計算することができるね。

静磁場中のスピノールの時間依存性

まずは、静磁場を $\boldsymbol{B} = (0, 0, B_z)$ として、静磁場中のスピノールの時間依存性を計算してみよう。パウリ方程式のスピン部分は

\begin{align}
i\hbar \frac{d }{d t} \left( \matrix{\chi_{\uparrow}(t) \cr \chi_{\downarrow}(t)} \right) &\ = \frac{e\hbar B_z}{2m_e}\, \left( \matrix{ 1 & 0 \cr 0 & -1} \right) \left( \matrix{\chi_{\uparrow}(t) \cr \chi_{\downarrow}(t)} \right) \\
&\ = \frac{e\hbar B_z}{2m_e}\,\left( \matrix{ \chi_{\uparrow}(t) \cr – \chi_{\downarrow}(t)} \right)
\end{align}

となるね。これは直ちに解くことができて、

\begin{align}
\chi( t ) = \left( \matrix{\chi_{\uparrow}(0) e^{-i\omega_L t} \cr \chi_{\downarrow}(0) e^{i\omega_L t}} \right) \ , \ \omega_L = \frac{e B_z}{2m_e}
\end{align}

となるね。$ |\chi( t )|^2 = |\chi( 0 )|^2 $ を満たすから、時間が経っても初期状態から変化しないことが分かるね。

振動磁場中のスピノールの時間依存性

次は、時間依存する磁場を $\boldsymbol{B}(t) = (0, 0, B_z\cos(\omega t))$ として、振動磁場中のスピノールの時間依存性を計算してみよう。パウリ方程式のスピン部分は

\begin{align}
i\hbar \frac{d }{d t} \left( \matrix{\chi_{\uparrow}(t) \cr \chi_{\downarrow}(t)} \right) &\ = \frac{e\hbar B_z\cos(\omega t)}{2m_e}\,
\left( \matrix{ 1 & 0 \cr 0 & -1} \right) \left( \matrix{\chi_{\uparrow}(t) \cr \chi_{\downarrow}(t)} \right) \\
&\ = \frac{e\hbar B_z\cos(\omega t)}{2m_e}\,\left( \matrix{ \chi_{\uparrow}(t) \cr – \chi_{\downarrow}(t)} \right)
\end{align}

となるね。これも直ちに解くことができて、

\begin{align}
\chi( t ) = \left( \matrix{\chi_{\uparrow}(0) e^{-i\frac{\omega_L}{\omega} \sin(\omega t)} \cr \chi_{\downarrow}(0) e^{i\frac{\omega_L}{\omega} \sin(\omega t)}} \right) \ , \
\omega_L = \frac{e B_z}{2m_e}
\end{align}

となるね。振動磁場を与えても $ |\chi( t )|^2 = |\chi( 0 )|^2 $ を満たすから、時間が経っても初期状態から変化しないことが分かるね。

静磁場に垂直方向の振動磁場を加えたときのスピノールの時間依存性

今度は、時間依存する磁場を $\boldsymbol{B}(t) = (0, B_y\sin(\omega t), B_z)$ として、振動磁場中のスピノールの時間依存性を計算してみよう。パウリ方程式のスピン部分は

\begin{align}
i\hbar \frac{d }{d t} \left( \matrix{\chi_{\uparrow}(t) \cr \chi_{\downarrow}(t)} \right) &\ =\frac{e\hbar }{2m_e}\left[B_y\sin(\omega t)\left( \matrix{ 0 & -i \cr i & 0} \right) + B_z\left( \matrix{ 1 & 0 \cr 0 & -1} \right)\right]
\left( \matrix{\chi_{\uparrow}(t) \cr \chi_{\downarrow}(t)} \right) \\
&\ =\frac{e\hbar }{2m_e}\left( \matrix{ B_z& -iB_y\sin(\omega t) \cr iB_y\sin(\omega t) & -B_z} \right) \left( \matrix{\chi_{\uparrow}(t) \cr \chi_{\downarrow}(t)} \right) \\
&\ = \frac{e\hbar }{2m_e}\,\left( \matrix{ B_z\chi_{\uparrow}(t)-iB_y\sin(\omega t)\chi_{\downarrow}(t) \cr iB_y\sin(\omega t) \chi_{\uparrow}(t) -B_z\chi_{\downarrow}(t)} \right)
\end{align}

と、 $\chi_{\uparrow}(t)$ と $\chi_{\downarrow}(t)$ に関する連立微分方程式が得られるね。解析的に解けるかはわからないけれども、数値的に解くことで、スピノールの時間依存性を計算することができそうだね。異常ゼーマン効果で分かれる2つの準位間のエネルギーに相当する振動磁場(電磁波)の角振動数を与えると、共鳴が起こってスピンの向きが変わるね。これは磁気共鳴現象って呼ばれるよ。これを次回の課題にするね。

円偏光を入射したときのスピノールの時間依存性

最後に、時間依存する磁場を $\boldsymbol{B}(t) = (0, B_0\sin(\omega t), B_0\cos(\omega t))$ として、円偏光を入射したときのスピノールの時間依存性を計算してみよう。パウリ方程式のスピン部分は

\begin{align}
i\hbar \frac{d }{d t} \left( \matrix{\chi_{\uparrow}(t) \cr \chi_{\downarrow}(t)} \right) &\ =\frac{e\hbar }{2m_e}\left[B_y\left( \matrix{ 0 & -i \cr i & 0} \right) + B_z\cos(\omega t)\left( \matrix{ 1 & 0 \cr 0 & -1} \right)\right]
\left( \matrix{\chi_{\uparrow}(t) \cr \chi_{\downarrow}(t)} \right) \\
&\ =\frac{e\hbar }{2m_e}\left( \matrix{ B_z\cos(\omega t) & -iB_y \cr iB_y & -B_z\cos(\omega t)} \right) \left( \matrix{\chi_{\uparrow}(t) \cr \chi_{\downarrow}(t)} \right) \\
&\ = \frac{e\hbar }{2m_e}\,\left( \matrix{ B_z\cos(\omega t)\chi_{\uparrow}(t)-iB_y\chi_{\downarrow}(t) \cr iB_y \chi_{\uparrow}(t) -B_z\cos(\omega t)\chi_{\downarrow}(t)} \right)
\end{align}

と、 $\chi_{\uparrow}(t)$ と $\chi_{\downarrow}(t)$ に関する連立微分方程式が得られるね。解析的に解けるかはわからないけれども、数値的に解くことで、スピノールの時間依存性を計算することができそうだね。円偏光でもスピンの向きは変化しそうだね。これも次回の課題にするね。