水素原子に電磁波（古典・円偏光）を加える場合のハミルトニアン

直線偏光の古典電磁波を水素原子に束縛された電子のハミルトニアンは以前示したね。今度は、円偏光の場合を導出するよ。偏光を考慮したベクトルポテンシャルは

\begin{align}
\boldsymbol{A} = \left\{ \matrix{ 0 \cr A_{y0} \cos(kx-\omega t + \phi_y) \cr A_{z0} \cos(kx-\omega t+ \phi_z)} \right.
\end{align}

と与えると、電場と磁場は次のようになるね。

\begin{align}
\boldsymbol{E} = -\frac{\partial \boldsymbol{A}}{\partial t} =\left\{ \matrix{ 0 \cr -\omega A_{y0} \sin(kx-\omega t + \phi_y) \cr -\omega A_{z0} \sin(kx-\omega t+ \phi_z)} \right.
\end{align}
\begin{align}
\boldsymbol{B} = \nabla \times \boldsymbol{A} =\left\{ \matrix{ 0 \cr k A_{z0} \sin(kx-\omega t + \phi_z) \cr -k A_{y0} \sin(kx-\omega t+ \phi_y)} \right.
\end{align}

これを電磁場中の電子のハミルトニアン

\begin{align}
\hat{H} = \hat{H}_0 + \frac{e}{m_e}\,\boldsymbol{A}\cdot \boldsymbol{p} + \frac{e}{m_e} \hat{\boldsymbol{S}}\cdot \boldsymbol{B} + \frac{e^2 \boldsymbol{A}^2}{2m_e}
\end{align}

に代入すると完成だね。

電磁波による状態遷移の数値計算

時間依存性を計算するために、ハミルトニアンの時間依存部分 $\hat{V}(t)$ の空間積分の項を見てみよう。

\begin{align}
\int \varphi_{n’l’m’\sigma’}^* \hat{V}(t)\varphi_{nlm\sigma} dV &\ = \int \varphi_{n’l’m’\sigma’}^* \left[
-i\frac{e}{\hbar} [\hat{H}_0,\boldsymbol{r} ]\cdot\boldsymbol{A} + \frac{e }{m_e} \hat{\boldsymbol{S}}\cdot\boldsymbol{B} + \frac{e^2 \boldsymbol{A}^2}{2m_e} \right] \varphi_{nlm\sigma} dV \\
&\ = \int \varphi_{n’l’m’\sigma’}^* \left[
-ie \frac{ E_{n’} – E_n }{\hbar}\,\boldsymbol{r}\cdot\boldsymbol{A} + \frac{e }{m_e} \hat{\boldsymbol{S}}\cdot\boldsymbol{B} + \frac{e^2 \boldsymbol{A}^2}{2m_e} \right] \varphi_{nlm\sigma} dV
\end{align}

スピンと２次の項を無視して、さらに長波長近似（電磁波の波長が原子サイズより十分に大きいとする近似）を行うと、

\begin{align}
\int \varphi_{n’l’m’}^* \hat{V}(t)\varphi_{nlm} dV &\ \simeq -ie \int \varphi_{n’l’m’}^* \left[
\frac{ E_{n’} – E_n }{\hbar}\,(yA_y +zA_z) \right] \varphi_{nlm} dV \\
&\ = \frac{-ie}{2}\,\frac{ E_{n’} – E_n }{\hbar} \left[A_{y0}\cos(\omega t – \phi_y) \int \varphi_{n’l’m’}^* y \varphi_{nlm} dV + A_{z0}\cos(\omega t – \phi_z) \int \varphi_{n’l’m’}^* z \varphi_{nlm} dV\right]
\end{align}

となるね。入射波を円偏光とするには、 $A_{y0} = A_{z0} = A_{0}$ かつ $\phi_y = 0 , \phi_z = \pm \pi/2 $ とすれば良いので、これを代入すると、

\begin{align}
\int \varphi_{n’l’m’}^* \hat{V}(t)\varphi_{nlm} dV &\ \simeq -ie \int \varphi_{n’l’m’}^* \left[
\frac{ E_{n’} – E_n }{\hbar}\,(yA_y +zA_z) \right] \varphi_{nlm} dV \\
&\ = \frac{-ieA_0}{2}\,\frac{ E_{n’} – E_n }{\hbar} \left[\cos(\omega t) \int \varphi_{n’l’m’}^* y \varphi_{nlm} dV \pm \sin(\omega t) \int \varphi_{n’l’m’}^* z \varphi_{nlm} dV\right]
\end{align}

となるね。このハミルトニアンを元に次回は水素原子の基底状態にいる電子に円偏光電磁波を与えてみるよ。