#rst2hooktail_source ============================================================ ファラデー効果 ============================================================ 掲示板で要望があったので書いてみました。 ファラデー効果とは、「分極を起こす物質に磁場をかけると、その磁場に平行な方向に進む 直線偏光光は、その偏光方向が回転していく。」という現象です。 それを古典振動子モデルで説明します。 物質中のマクスウェル方程式 ============================= これから、物質中のイオンに束縛される電子の運動方程式を考えたいのですが、その前に物質中の物理量の 満たす関係式を復習しておきます。 電場 $\bm{E}$ と電束密度 $\bm{D}$ 、分極 $\bm{P}$ の間には、 <tex> \bm{D} = \varepsilon _0 \bm{E} + \bm{P} \tag{##} </tex> という関係が、 磁場 $\bm{H}$ と磁束密度 $\bm{B}$ 、磁化 $\bm{M}$ の間には、 <tex> \bm{H} = \frac{1}{\mu_0}\bm{B} - \bm{M} </tex> これは物質の透磁率 $\mu$ を用いて、 <tex> \bm{B} = \mu_0 \bm{H} + \mu_0 \bm{M} = \mu \bm{H} \tag{##} </tex> という関係が成り立ちます。 一方、マクスウェル方程式は、原子に束縛されていない自由な電荷を添え字の $f$ で表わすと、 <tex> \nabla \cdot \bm{D} = \rho_f </tex> <tex> \nabla \times \bm{E} = - \frac{\partial\bm{B}}{\partial t} </tex> <tex> \nabla \cdot \bm{B} = 0 </tex> <tex> \nabla \times \bm{H} = \bm{J}_f + \frac{\partial\bm{D}}{\partial t} </tex> 今回は、物質中の自由な電荷はなく、かつ、磁化は弱いと考え無視します。 よって、扱うマクスウェル方程式を書きなおすと、 <tex> \nabla \cdot \bm{D} = 0 \tag{##} </tex> <tex> \nabla \times \bm{E} = - \mu \frac{\partial\bm{H}}{\partial t} \tag{##} </tex> <tex> \nabla \cdot \bm{B} = 0 \tag{##} </tex> <tex> \nabla \times \bm{H} = \frac{\partial\bm{D}}{\partial t} \tag{##} </tex> 式 $( 4 )$ に $\nabla \times$ を施して、式 $( 6 )$ を利用すると、 <tex> \nabla \times (\nabla \times \bm{E}) &= \mu \frac{\partial}{\partial t} (\nabla \times \bm{H} ) \\ &= - \mu \frac{\partial^2 \bm{D}}{\partial t^2} \tag{##} </tex> この式の左辺は、 <tex> \nabla \times (\nabla \times \bm{E}) = \nabla (\nabla \cdot \bm{E}) - \triangle \bm{E} </tex> のように変形出来るので、結局、式 $( 7 )$ は、 <tex> \nabla (\nabla \cdot \bm{E}) - \triangle \bm{E} = - \mu \frac{\partial^2 \bm{D}}{\partial t^2} \tag{##} </tex> となります。 この式 $( 8 )$ は、左辺第一項がなければ波動方程式になります。 実はこの節を書いたのは、式 $( 8 )$ をあとで使うことになるからです。 この式は重要ですので、こういう式があったな程度には覚えておいてください。 振動子と光のの古典的モデル ============================= 少々長いので、何をしたいかを簡潔に話しておくと、 $\bm{E}$ と $\bm{D}$ の間の関係を求めて 式 $( 8 )$ に代入したいのです。 さて、下の図のように透明な物質にコイルを巻き磁場Bを作り、そこに直線偏光の光を 磁場と平行に入射することを考えます。 .. image :: chromel-faradayEffect-01.png ここで、物質中の正イオン(原子核)に束縛された最外殻電子について考えます。 古典的な運動方程式は、下の様なモデルで考えることができます。 座標は、磁場の方向をz軸、上方向をy軸、画面手前方向をx軸とする右手系とします。 <tex> m(\ddot{\bm{r}} + \gamma \dot{{r}} + \omega_0^2 \bm{r} ) = \bm{F} </tex> $\gamma$ は正の値を持つ摩擦による減衰項、 $\omega_0$ はイオンとの引力が作る復元力の角振動数 です。 $\bm{F}$ は外力であり、 今回の場合、 $\bm{F}=-e(\bm{E}+\bm{v} \times \bm{B})$ となります。 成分ごとに式を書いてみます。 <tex> \ddot{x} + \gamma \dot{x} + \omega_0^2 x = -\frac{e}{m}(E_x + B \dot{y}) \tag{##} </tex> <tex> \ddot{y} + \gamma \dot{y} + \omega_0^2 y = -\frac{e}{m}(E_y - B \dot{x}) \tag{##} </tex> <tex> \ddot{z} + \gamma \dot{z} + \omega_0^2 z = 0 </tex> z方向には特別な動きがないことがわかります。議論の中心となるのは、xy平面での運動です。 この方程式を扱うちょっとしたテクニックがあります。 現在、二変数が二つの式にちらばっていますが、ふたつの式の線形結合をとることによって、 独立な一変数の式、ふたつに変えるのです。 それには線形代数の対角化を利用します。 演算子を $D$ と書き、あたかも係数のように書くことができます。 <tex> f(D)=D^2 + \gamma D + \omega_0^2 </tex> とすると、 式 $( 9 )$ 、 $( 10 )$ は、行列を用いて、 <tex> \begin{pmatrix} f(D) & \frac{eB}{m} \\ -\frac{eB}{m} & f(D) \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} -\frac{eE_x}{m} \\ -\frac{eE_y}{m} \end{pmatrix} \tag{##} </tex> これを対角化すれば、うまく変数をまとめることができます。 上の式の二次正方行列をAとすると、 <tex> det(A-\lambda I)=(\lambda-f(D))^2 + (\frac{eB}{m})^2 = 0 </tex> より、 <tex> \lambda = f(D) \pm i \frac{eB}{m} </tex> ここで、式 $( 11 )$ に左から二次の行ベクトルをかけることは、式 $( 9 )$ 、 $( 10 )$ の 線形結合をとることに対応します。よって、あまり聞きなれませんが、 今回求めるのは、Aの固有行ベクトルを求めればいいことになります。 <tex> \begin{pmatrix} 1 & X \end{pmatrix} (A-\lambda I) = \begin{pmatrix} 1 & X \end{pmatrix} \begin{pmatrix} \mp i \frac{eB}{m} & \frac{eB}{m} \\ - \frac{eB}{m} & \mp i \frac{eB}{m} \end{pmatrix} = \begin{pmatrix} 0 & 0 \end{pmatrix} </tex> をみたすのは、 <tex> \begin{pmatrix} 1 & \mp i \end{pmatrix} </tex> であることが分りました。これから $(9) \ \pm i(10)$ を作ると、確かに、 <tex> (f(D)-i\frac{eB}{m}D)(x+iy)=-\frac{e(E_x+iE_y)}{m} \tag{##} </tex> <tex> (f(D)+i\frac{eB}{m}D)(x-iy)=-\frac{e(E_x-iE_y)}{m} \tag{##} </tex> となって、きれいにまとめることができました。 特解と一般解 ============================= これをまともに解くのは大変なので、定性的な説明をしておきます。 この微分方程式は、式 $( 12 )$ 、 $( 13 )$ の右辺をゼロとした式を斉次方程式といい、 その解を一般解といいます。一般解の線形結合は、やはり一般解です。 そして、右辺がゼロでない場合、その方程式は非斉次方程式といわれ、 その解の一つを特解といいます。この方程式のすべての解は、一般解の線形結合+特解 で尽くされます。 今回の場合、式 $( 12 )$ 、 $( 13 )$ の右辺が振動項で、 右回り円偏光 [*]_ $E_x+iE_y=E_R e^{i(kz-\omega t)}$ の時や、 左回り円偏光 $E_x-iE_y=E_L e^{i(kz-\omega t)}$ の時を考えます。 .. [*] なぜ円偏光と言えるのか考えてみましょう。右回りの場合、電場が実数であることを考えれば、オイラーの公式 $e^{i\theta}=\cos \theta + i \sin \theta$ より、 $E_x=E_R \cos (kz -\omega t)$ かつ $E_y=E_R \sin (kz -\omega t)$ となり、z=一定の場所で変化を追うと、進行方向に向かって右回りであることがわかります。 非斉次項が振動項の場合、この方程式の解はどんな初期条件であっても、 時間がたつと振動項に対応する定常応答に落ち着きます。 一方、一般解は大きさが指数関数的に減少する項であり減衰項と呼ばれます。 これは初期条件から定常応答に移行する際のつじつま合わせの役割をします。 以降、定常応答だけを考えます。 右回り円偏光の入射(rとEの関係) ================================= まず、右回り偏光 $E_x+iE_y=E_R e^{i(kz-\omega t)}$ の時だけ考えます。 このとき、複素共役を $*$ で表わすと、 $E_x-iE_y=(E_x+iE_y)^*= E_R e^{-i(kz-\omega t)}$ となりますから、式 $( 12 )$ については、 $x+iy = \alpha E_R e^{i(kz - \omega t)} =E_x + iE_y$ が定常応答になり、 その結果、微分演算子 $D$ を複素数 $-i \omega$ で置き換えることができます。 一方、式 $( 13 )$ については、 $x-iy = \alpha^\prime E_R e^{-i(kz - \omega t)}$ が定常応答になり、 こっちの式では、微分演算子 $D$ を複素数 $i \omega$ で置き換えることができます。 $\alpha$ も $\alpha^\prime$ も複素数であることに注意してください。 書きなおすと式 $( 12 )$ は、 <tex> (f(-i \omega)+(-i\frac{eB}{m})(-i \omega)) \alpha e^{i(kz-\omega t)}=-\frac{eE_R}{m}e^{i(kz-\omega t)} \tag{##} </tex> <tex> (f(i \omega)+(i\frac{eB}{m})(i \omega)) \alpha^\prime e^{-i(kz-\omega t)}=-\frac{eE_R}{m}e^{-i(kz-\omega t)} \tag{##} </tex> $f(D)=D^2 + \gamma D + \omega_0^2$ を代入してやると、 <tex> (- \omega^2-\frac{eB}{m}\omega + \omega_0^2 - i \gamma \omega) \alpha e^{i(kz-\omega t)}=-\frac{eE_R}{m}e^{i(kz-\omega t)} \tag{##} </tex> <tex> (- \omega^2-\frac{eB}{m}\omega + \omega_0^2 + i \gamma \omega) \alpha^\prime e^{-i(kz-\omega t)}=-\frac{eE_R}{m}e^{-i(kz-\omega t)} \tag{##} </tex> となり、 $\alpha_1 \ , \ \alpha_2$ を実数として、 <tex> \alpha &= \frac{e}{m(- \omega^2-\frac{eB}{m}\omega + \omega_0^2 - i \gamma \omega)} \\ &= \frac{e(- \omega^2-\frac{eB}{m}\omega + \omega_0^2 + i \gamma \omega)}{m((- \omega^2-\frac{eB}{m}\omega + \omega_0^2)^2 + (\gamma \omega)^2)} \\ &= \alpha_1 + i \alpha_2 = \alpha^{\prime *} </tex> ひとつ注意しておくと、 $\gamma$ がゼロ(抵抗力が働かない)時は、 $\alpha_2$ がゼロとなり、 $\alpha$ が実数になります。 ですから、 $x \ ,\ y\ ,\ E_i$ はすべて実数であることを考慮すると、 式 $( 16 )$ 、 $( 17 )$ は、どちらか一つだけで十分な意味を持ち、 以下のように、 <tex> x+iy &= \alpha (E_x + iE_y) \\ &= (\alpha_1 + i \alpha_2)(E_x + i E_y) \\ &= (\alpha_1 E_x -\alpha_2 E_y)+i(\alpha_1 E_y + \alpha_2 E_x) </tex> を実部、虚部に分離できて、 <tex> x = \alpha_1 E_x -\alpha_2 E_y </tex> <tex> y = \alpha_1 E_y + \alpha_2 E_x </tex> 左回り円偏光の入射(rとEの関係) ================================= 左回り円偏光についても同様に、こんどは、 $E_x + i E_y = E_L e^{ -i ( kz-\omega t)}$ の時を考えます。 式 $( 14 )$ に対応する式は、 <tex> (f(i \omega)+(-i\frac{eB}{m})(i \omega)) \beta e^{-i(kz-\omega t)}=-\frac{eE_L}{m}e^{-i(kz-\omega t)} </tex> ですから、 <tex> (- \omega^2 + \frac{eB}{m}\omega + \omega_0^2 + i \gamma \omega) (x+iy)=-\frac{e}{m}(E_x+iE_y) \tag{##} </tex> となり、 $\beta$ は以下のようになり、 <tex> \beta &= \frac{e}{m(- \omega^2 + \frac{eB}{m}\omega + \omega_0^2 + i \gamma \omega)} \\ &= \frac{e(- \omega^2+\frac{eB}{m}\omega + \omega_0^2 - i \gamma \omega)}{m((- \omega^2+\frac{eB}{m}\omega + \omega_0^2)^2 + (\gamma \omega)^2)} \\ &= \beta_1 + i \beta_2 = \beta^{\prime *} </tex> 結局、位置 $\bm{r}$ と電場 $\bm{E}$ の関係は、 <tex> x = \beta_1 E_x - \beta_2 E_y </tex> <tex> y = \beta_1 E_y + \beta_2 E_x </tex> となります。 電束密度D、分極P、電場Eの関係 ================================ 分極 $\bm{P}=(P_x,P_y,P_z)$ は、 $x\ ,\ y$ に $-ne$ を掛けたものであり、 <tex> \bm{P} &=-ne\bm{r} \\ &=-ne \begin{pmatrix} x \\ y \\ 0 \end{pmatrix} </tex> となりますから、まず右回り偏光の時は、 <tex> \bm{D} &= \varepsilon_0 \bm{E} + \bm{P} \\ &=\varepsilon_0 \begin{pmatrix} E_x \\ E_y \\ 0 \end{pmatrix} -ne \begin{pmatrix} x \\ y \\ 0 \end{pmatrix} \\ &=\varepsilon_0 \begin{pmatrix} E_x \\ E_y \\ 0 \end{pmatrix} -ne \begin{pmatrix} \alpha_1 E_x -\alpha_2 E_y \\ \alpha_1 E_y + \alpha_2 E_x \\ 0 \end{pmatrix} \\ &= \begin{pmatrix} \varepsilon_0 -ne \alpha_1 & ne \alpha_2 & 0 \\ -ne \alpha_2 & \varepsilon_0 -ne \alpha_1 & 0 \\ 0&0&\varepsilon_0 \end{pmatrix} \begin{pmatrix} E_x \\ E_y \\ 0 \end{pmatrix} &= \begin{pmatrix} \varepsilon_{R1} & -\varepsilon_{R2} & 0 \\ \varepsilon_{R2} & \varepsilon_{R1} & 0 \\ 0&0&\varepsilon_0 \end{pmatrix} \begin{pmatrix} E_x \\ E_y \\ 0 \end{pmatrix} </tex> これで、やっと式 $( 8 )$ が利用できます。 電場 $\bm{E}=(E_x,E_y,0)$ と進行方向 $\bm{k}=(0,0,k)$ は垂直なので、 <tex> \nabla (\nabla \cdot \bm{E}) &= \bm{k}(\bm{k}\cdot\bm{E}) \\ &= 0 </tex> ですから、式 $( 8 )$ は、以下のように波動方程式になります。 <tex> - \triangle E_x &= - \mu \frac{\partial^2}{\partial t^2} D_x \\ &= - \mu \frac{\partial^2}{\partial t^2} (\varepsilon_{R1}E_x -\varepsilon_{R2} E_y) \tag{##} </tex> <tex> - \triangle E_y &= - \mu \frac{\partial^2}{\partial t^2} D_y \\ &= - \mu \frac{\partial^2}{\partial t^2} (\varepsilon_{R2}E_x + \varepsilon_{R1} E_y) \tag{##} </tex> さっき行った対角化の方法によると、 式について、 $(19) \pm i (20)$ をとれば、二つの共役な式が出てきて、一方だけを書くと、 <tex> - \triangle (E_x+iE_y) &= - \mu (\varepsilon_{R1} + i\varepsilon_{R2})\frac{\partial^2}{\partial t^2} (E_x + i E_y) \tag{##} </tex> 右回り偏光なので、 $E_x+iE_y = E_R e^{i(kz -\omega t)}$ を代入すると、 分散関係が出てきます。 <tex> |\bm{k}|^2 = k_z^2 = \mu (\varepsilon_{R1} + i\varepsilon_{R2}) \omega^2 </tex> <tex> E_x+iE_y &= E_R e^{i(kz -\omega t) } \\ &= E_R e^{-i \omega (t-\sqrt{\mu (\varepsilon_{R1} + i \varepsilon_{R2})} z )} </tex> この式は、 $\varepsilon_{R2}$ がノンゼロの時減衰することが解ります。 つまり、それは抵抗項 $\gamma$ がノンゼロの時を表しますから、 電子に抵抗があると光は減衰しながら伝わっていくことになり、 物理的感覚に合っています。 左回り円偏光の場合は、前は $\alpha$ に対応して $\varepsilon$ が決定されたのに対し、 $\beta$ に対応して $\varepsilon$ が決定されます。その分散関係は、 <tex> |\bm{k}|^2 = k_z^2 = \mu (\varepsilon_{L1} + i \varepsilon_{L2}) \omega^2 </tex> であり、 <tex> E_x+iE_y &= E_L e^{-i(kz -\omega t)} \\ &= E_L e^{i \omega (t-\sqrt{\mu (\varepsilon_{L1} + i \varepsilon_{L2})}z)} </tex> ここまで来て注目するべきところは、 抵抗項 $\gamma$ がゼロの時を考えると、 光の振動数 $\omega$ が同じでも、右回り円偏光と左回り円偏光で、 波数kが異なってくることです。 .. image:: chromel-faradayEffect-02.png 直線偏光は、右回りと左回りの円偏光の重ね合わせと考えられますから、 その二つの偏光の波数が異なるということは、位置と共に波の位相がずれていく、 つまり、直線偏光の方向が、異なってくるということになります。 余談 ======= 今回の話は、磁場によるローレンツ力という、速度に垂直に働く力が原因です。 そういえば、フーコーの振り子もコリオリ力という速度に垂直に働く力が原因で、 振り子の回転方向が回転していく現象なので、方程式としては似ているなと思いました。 @@author:クロメル@@ @@accept:2007-09-26@@ @@accept:2009-05-04@@ @@category:電磁気学@@ @@id:faradayEffect@@