logo

磁気双極子が生成する磁場 📂電磁気学

磁気双極子が生成する磁場

説明

596C5A923.jpg

磁気双極子 m\mathbf{m} によるベクトルポテンシャルは、磁気双極子モーメントで次のように与えられる。

Adip(r)=μ04πm×r^r2=μ04πmsin⁡θr2ϕ^ \mathbf{A}_{\text{dip}}(\mathbf{r}) = \dfrac{\mu_{0}}{4 \pi} \dfrac{\mathbf{m} \times \hat{\mathbf{r}}}{r^2} = \dfrac{\mu_{0}}{4 \pi} \dfrac{m\sin\theta}{r^{2}} \hat{\boldsymbol{\phi}}

今、m\mathbf{m}が原点にあり、そしてzz軸と並行だとすると、上の図のようになる。磁場はベクトルポテンシャルのカールであるから、球面座標系で次のようになる。

B=∇×Adip=1rsin⁡θ[∂(sin⁡θAϕ)∂θ−∂Aθ∂ϕ]r^+1r[1sin⁡θ∂Ar∂ϕ−∂(rAϕ)∂r]θ^+1r[∂(rAθ)∂r−∂Ar∂θ]ϕ^ \begin{align*} \mathbf{B} = \nabla \times \mathbf{A}_{\text{dip}} &= \frac{1}{r\sin\theta} \left[\frac{\partial (\sin\theta A_\phi)}{\partial \theta} - \frac{\partial A_\theta}{\partial \phi} \right]\mathbf{\hat{\mathbf{r}}} + \frac{1}{r}\left[\frac{1}{\sin\theta} \frac{\partial A_{r}}{\partial \phi} - \frac{\partial (rA_\phi)}{\partial r} \right] \boldsymbol{\hat{\boldsymbol{\theta}}} \\[1em] &\quad+ \frac{1}{r} \left[\frac{\partial (rA_\theta)}{\partial r}-\frac{\partial A_{r}}{\partial \theta} \right]\boldsymbol{\hat \phi} \end{align*}

各成分を計算すると以下のようになる。Adip\mathbf{A}_{\text{dip}}の成分はAϕA_{\phi}のみ存在するので、

Br=1rsin⁡θ∂(sin⁡θAϕ)∂θ=1rsin⁡θ∂∂θ(μ04πmsin⁡2θr2)=μ04π2mcos⁡θr3Bθ=−1r∂(rAϕ)∂r=−1r∂∂r(μ04πmsin⁡θr)=μ04πmsin⁡θr3Bϕ=0 \begin{align*} B_{r} &= \frac{1}{r\sin\theta} \frac{\partial (\sin\theta A_\phi)}{\partial \theta} = \frac{1}{r\sin\theta}\frac{\partial }{\partial \theta}\left( \dfrac{\mu_{0}}{4 \pi} \dfrac{m\sin^{2}\theta}{r^{2}} \right) = \dfrac{\mu_{0}}{4 \pi} \dfrac{2m\cos\theta}{r^{3}} \\[1em] B_{\theta} &= - \frac{1}{r} \frac{\partial (rA_\phi)}{\partial r} = - \frac{1}{r} \frac{\partial }{\partial r} \left( \dfrac{\mu_{0}}{4 \pi} \dfrac{m\sin\theta}{r} \right) = \dfrac{\mu_{0}}{4 \pi} \dfrac{m\sin\theta}{r^{3}} \\[1em] B_{\phi} &= 0 \end{align*}

したがって、磁気双極子が作る磁場は次のようになる。

Bdip(r,θ)=μ04πmr3(2cos⁡θr^+sin⁡θθ^) \begin{equation} \mathbf{B}_{\text{dip}} (r,\theta)=\frac{\mu_{0} }{4 \pi }\frac{m}{r^3}(2\cos\theta \hat{\mathbf{r}} + \sin\theta \hat{\boldsymbol{\theta}}) \end{equation}

面白いことに、これは電気双極子が作る電場と正確に同じ式である。

公式

上の式を座標系に依存しない形に変えると、次のようになる。

Bdip(r)=μ04π1r3[3(m⋅r^)r^−m] \mathbf{B}_{\text{dip}}(\mathbf{r}) = \frac{\mu_{0}}{4 \pi}\frac{1}{r^3}[3 (\mathbf{m} \cdot \hat{\mathbf{r}}) \hat{\mathbf{r}} - \mathbf{m}]

導出

まず、球面座標系の単位ベクトルを直交座標系の単位ベクトルで表すと、次のようになる。

r^= cos⁡ϕsin⁡θx^+sin⁡ϕsin⁡θy^+cos⁡θz^θ^= cos⁡ϕcos⁡θx^+sin⁡ϕcos⁡θy^−sin⁡θz^ \begin{align*} \hat{\mathbf{r}} =&\ \cos\phi \sin\theta \hat{\mathbf{x}} + \sin\phi \sin\theta\hat{\mathbf{y}} + \cos\theta\hat{\mathbf{z}} \\ \hat{\boldsymbol{\theta}} =&\ \cos\phi \cos\theta \hat{\mathbf{x}} + \sin\phi \cos\theta\hat{\mathbf{y}} - \sin\theta\hat{\mathbf{z}} \end{align*}

したがって、(1)(1)の括弧内の式を計算すると、次のようになる。

2cos⁡θr^+sin⁡θθ^= 2cos⁡ϕsin⁡θcos⁡θx^+2sin⁡ϕsin⁡θcos⁡θy^+2cos⁡2θz^+cos⁡ϕsin⁡θcos⁡θx^+sin⁡ϕcos⁡θsin⁡θy^−sin⁡2θz^= 3cos⁡ϕsin⁡θcos⁡θx^+3sin⁡ϕsin⁡θcos⁡θy^+3cos⁡2θz^−(sin⁡2θ+cos⁡2θ)z^= 3cos⁡θ(cos⁡ϕsin⁡θx^+sin⁡ϕsin⁡θy^+cos⁡θz^)−z^= 3(m^⋅r^)r^−z^ \begin{align*} & 2\cos\theta \hat{\mathbf{r}} + \sin \theta \hat{\boldsymbol{\theta}} \\ =&\ 2 \cos\phi \sin\theta \cos\theta \hat{\mathbf{x}} + 2 \sin\phi \sin\theta \cos\theta \hat{\mathbf{y}} + 2 \cos^2 \theta \hat{\mathbf{z}} \\ & + \cos\phi \sin\theta \cos\theta \hat{\mathbf{x}} + \sin\phi \cos\theta \sin\theta \hat{\mathbf{y}} -\sin^2\theta \hat{\mathbf{z}} \\ =&\ 3\cos\phi \sin\theta \cos\theta \hat{\mathbf{x}} + 3 \sin\phi \sin\theta \cos\theta \hat{\mathbf{y}} + 3 \cos^2 \theta \hat{\mathbf{z}} -(\sin^2\theta + \cos^2\theta)\hat{\mathbf{z}} \\ =&\ 3 \cos\theta (\cos\phi \sin\theta \hat{\mathbf{x}} + \sin\phi \sin\theta\hat{\mathbf{y}} + \cos\theta\hat{\mathbf{z}}) - \hat{\mathbf{z}} \\ =&\ 3 (\hat{\mathbf{m}} \cdot \hat{\mathbf{r}}) \hat{\mathbf{r}} - \hat{\mathbf{z}} \end{align*}

最後の等式はcos⁡θ=m^⋅r^\cos\theta = \hat{\mathbf{m}} \cdot \hat{\mathbf{r}}のため成り立つ。これで、次の結果を得る。

Bdip(r,θ)=μ04πmr3(2cos⁡θr^+sin⁡θθ^)=μ04πmr3[3(m^⋅r^)r^−z^]=μ04π1r3[3(m⋅r^)r^−mz^]=μ04π1r3[3(m⋅r^)r^−m]=Bdip(r)) \begin{align*} \mathbf{B}_{\text{dip}}(r,\theta) &=\frac{\mu_{0} }{4 \pi }\frac{m}{r^3}(2\cos\theta \hat{\mathbf{r}} + \sin\theta \hat{\boldsymbol{\theta}}) \\[1em] &= \frac{\mu_{0} }{4 \pi }\frac{m}{r^3}[3 (\hat{\mathbf{m}}\cdot \hat{\mathbf{r}}) \hat{\mathbf{r}} - \hat{\mathbf{z}}] \\[1em] &= \frac{\mu_{0} }{4 \pi }\frac{1}{r^3}[3 (\mathbf{m} \cdot \hat{\mathbf{r}}) \hat{\mathbf{r}} - m \hat{\mathbf{z}}] \\[1em] &= \frac{\mu_{0}}{4 \pi }\frac{1}{r^3}[3 (\mathbf{m} \cdot \hat{\mathbf{r}}) \hat{\mathbf{r}} - \mathbf{m}] \\[1em] &= \mathbf{B}_{\text{dip}}(\mathbf{r})) \end{align*}

■