📂電磁気学

動く点電荷が作る電場

概要 ¹

動く点電荷が作る電磁場は次の通りだ。

$$ \begin{align*} \mathbf{E}(\mathbf{r}, t) &= \frac{q}{4\pi\epsilon_{0}} \frac{\cR} {( \bcR\cdot \mathbf{u} )^3 } \left[(c^2-v^2)\mathbf{u} +\bcR\times (\mathbf{u} \times \mathbf{a} ) \right] \\ \mathbf{B} (\mathbf{ r}, t) &=\frac{1}{c} \crH\times \mathbf{ E } (\mathbf{ r}, t) \end{align*} $$

説明

電場についての誘導過程を紹介する。

誘導

動く点電荷が作る電場、磁場はリエナール-ビーケルト電位を使って求めることができる。

$$ V(\mathbf{r}, t)= \frac{1}{4\pi \epsilon_{0}} \frac{qc}{ (\cR c -\bcR\cdot \mathbf{v})} ,\quad \mathbf{A}(\mathbf{r}, t) = \frac{ \mathbf{v} } {c^2} V(\mathbf{r}, t) $$

さらに電磁場は下の式で求めることができる。

$$ \mathbf{ E} = -\nabla V -\frac{\partial \mathbf{ A} }{\partial t},\quad \mathbf{B}=\nabla \times \mathbf{A} $$

$\nabla V$と$\dfrac{\partial \mathbf{A}}{\partial t}$を順に求めてみよう。まず$\nabla V$を計算してみると

$$ \nabla V=\frac{qc}{4 \pi \epsilon_{0}} \nabla \dfrac{1}{(\cR c - \bcR \cdot \mathbf{v}) } $$

この時$\dfrac{ d}{dx} \left( \dfrac{1}{f(x)} \right)=\dfrac{-1}{ [f(x)]^2} f^{\prime}(x)$なので

$$ \begin{equation} \nabla V = \frac{qc}{4\pi\epsilon_{0}} \frac{-1}{ (\cR c -\bcR \cdot \mathbf{v} )^2}\nabla (\cR c -\bcR \cdot \mathbf{v} ) \end{equation} $$

従って$\nabla \cR$と$\nabla (\bcR \cdot \mathbf{v})$を計算すればいい。まず$\nabla ( \bcR \cdot \mathbf{v} )$から計算してみると乗法規則2によって

$$ \begin{equation} \nabla (\bcR \cdot \mathbf{v}) = (\bcR \cdot \nabla ) \mathbf{v}+(\mathbf{v} \cdot \nabla ) \bcR + \bcR\times (\nabla \times \mathbf{v}) +\mathbf{v}\times (\nabla \times \bcR) \end{equation} $$

第一項

$$ \begin{align*} (\bcR \cdot \nabla) \mathbf{v} &= \left( \cR_{x} \frac{\partial}{\partial x} + \cR _{y} \frac{\partial}{\partial y}+\cR_{z}\frac{\partial }{\partial z} \right)\mathbf{v} \\[1em] &= \cR_{x}\frac{d \mathbf{v} }{dt_{r}}\frac{\partial t_{r}}{\partial x} +\cR_{y}\frac{d \mathbf{v} }{dt_{r}}\frac{\partial t_{r}}{\partial y} + \cR_{z}\frac{d \mathbf{v} }{dt_{r}}\frac{\partial t_{r}}{\partial z} \\[1em] &= \frac{d \mathbf{v} }{ dt_{r}} \left( \cR_{x}\frac{\partial t_{r}}{\partial x} + \cR_{y}\frac{\partial t_{r}}{\partial y}+\cR_{z}\frac{\partial t_{r}}{\partial z}\right) \\[1em] &= \mathbf{a} ( \bcR \cdot \nabla t_{r}) \end{align*} $$

ここで$\mathbf{a}$は遅れ時間での粒子（点電荷）の加速度だ。

第二項

$\bcR=\mathbf{r}-\mathbf{w}$なので

$$ \begin{equation} (\mathbf{v} \cdot \nabla ) \bcR=(\mathbf{v}\cdot \nabla)\mathbf{r} -(\mathbf{v}\cdot \nabla ) \mathbf{w} \end{equation} $$

任意のベクトル$\mathbf{A}=(A_{x}, A_{y}, A_{z})$に対して

$$ \begin{align*} (\mathbf{A} \cdot \nabla ) \mathbf{r} &= \left( A_{x}\frac{\partial}{\partial x}+A_{y}\frac{\partial}{\partial y}+A_{z}\frac{\partial}{\partial z} \right) ( x \hat{\mathbf{x}} + y \hat{\mathbf{y}} +z\hat{ \mathbf{z} }) \\ &= A_{x}\hat{\mathbf{x}}+A_{y} \hat{\mathbf{y}} + A_{z}\hat{\mathbf{z}} \\ &= \mathbf{A} \end{align*} $$

そして

$$ \begin{align*} (\mathbf{A} \cdot \nabla ) \mathbf{w}(t_{r}) &= \left( A_{x}\frac{\partial}{\partial x}+A_{y}\frac{\partial}{\partial y}+A_{z}\frac{\partial}{\partial z} \right) \mathbf{w} \\[1em] &= A_{x}\frac{\partial \mathbf{w} } {\partial x}+A_{y}\frac{\partial \mathbf{w} }{\partial y}+A_{z}\frac{\partial \mathbf{w} }{\partial z} \\[1em] &= A_{x}\frac{\partial \mathbf{w} }{\partial t_{r}}\frac{\partial t_{r}} {\partial x}+A_{y}\frac{\partial \mathbf{w} }{\partial t_{r}}\frac{\partial t_{r}}{\partial y}+A_{z}\frac{\partial \mathbf{w} }{\partial t_{r}}\frac{\partial t_{r}}{\partial z} \\[1em] &= A_{x} \mathbf{v}\frac{\partial t_{r}} {\partial x}+A_{y}\mathbf{v}\frac{\partial t_{r}}{\partial y}+A_{z}\mathbf{v}\frac{\partial t_{r}}{\partial z} \\[1em] &= \mathbf{v} \left( A_{x} \frac{\partial t_{r}} {\partial x}+A_{y}\frac{\partial t_{r}}{\partial y}+A_{z}\frac{\partial t_{r}}{\partial z}\right) \\[1em] &= \mathbf{v}(\mathbf{A} \cdot \nabla t_{r}) \end{align*} $$ なので$(3)$は

$$ (\mathbf{v} \cdot \nabla ) \bcR=\mathbf{v} -\mathbf{v}(\mathbf{v}\cdot \nabla t_{r} ) $$

第三項

まず次が成り立つ。

$$ \begin{align*} \nabla \times \mathbf{v} &= \left( \frac{\partial v_{z}}{\partial y} - \frac{\partial v_{y}}{\partial z} \right)\hat{\mathbf{x}} + \left( \frac{\partial v_{x}}{\partial z} - \frac{\partial v_{z}}{\partial x} \right)\hat{\mathbf{y}} + \left( \frac{\partial v_{y}}{\partial x} - \frac{\partial v_{x}}{\partial y} \right)\hat{\mathbf{z}} \\[1em] &= \left( \frac{\partial v_{z}}{\partial t_{r}} \frac{t_{r}}{\partial y} - \frac{\partial v_{y}}{\partial t_{r}} \frac{t_{r}}{\partial z} \right)\hat{\mathbf{x}} + \left( \frac{\partial v_{x}}{\partial t_{r}} \frac{t_{r}}{\partial z} - \frac{\partial v_{z}}{\partial t_{r}} \frac{t_{r}}{\partial x} \right)\hat{\mathbf{y}} + \left( \frac{\partial v_{y}}{\partial t_{r}} \frac{t_{r}}{\partial x} - \frac{\partial v_{x}}{\partial t_{r}} \frac{t_{r}}{\partial y} \right)\hat{\mathbf{z}} \\[1em] &= \Big( a_{z} (\nabla t_{r})_{y} - a_{y} (\nabla t_{r})_{z} \Big)\hat{\mathbf{x}} + \Big( a_{x} (\nabla t_{r})_{z} - a_{z} (\nabla t_{r})_{x} \Big)\hat{\mathbf{y}} + \Big( a_{y} (\nabla t_{r})_{x} - a_{x} (\nabla t_{r})_{y} \Big)\hat{\mathbf{z}} \\[1em] &= -\mathbf{a} \times \nabla t_{r} \end{align*} $$

従って

$$ \bcR \times ( \nabla \times \mathbf{v} )= -\bcR\times (\mathbf{a} \times \nabla t_{r}) $$

第四項

$\bcR=\mathbf{r}-\mathbf{w}$であり、$\nabla \times \mathbf{r}=0$なので

$$ \nabla \times \bcR = \nabla \times \mathbf{r} -\nabla \times \mathbf{w}=-\nabla \times \mathbf{w} $$

上で計算した$\nabla \times \mathbf{v} = -\mathbf{a} \times \nabla t_{r}$の結果を利用すると$\nabla \times \mathbf{w} = -\mathbf{v} \times \nabla t_{r}$であることが分かる。従って

$$ \mathbf{v}\times (\nabla \times \bcR)=\mathbf{v} \times (- \nabla \times \mathbf{w})=\mathbf{v} \times (\mathbf{v} \times \nabla t_{r}) $$

これで上の計算結果を$(2)$に代入すると

$$ \begin{align} \nabla (\bcR \cdot \mathbf{v} ) &= \mathbf{a} ( \bcR \cdot \nabla t_{r}) + \mathbf{v} -\mathbf{v}(\mathbf{v} \cdot \nabla t_{r}) -\color{blue}{\abcR\times(\mathbf{a} \times \nabla t_{r} )} +\color{green}{\mathbf{v} \times (\mathbf{v} \times \nabla t_{r})} \nonumber \\[1em] &= \mathbf{a} ( \bcR \cdot \nabla t_{r}) + \mathbf{v} -\mathbf{v}(\mathbf{v} \cdot \nabla t_{r}) -\color{blue}{\left[\mathbf{a} ( \abcR \cdot \nabla t_{r}) - \nabla t_{r}(\abcR \cdot \mathbf{a} ) \right]}+ \color{green}{\left[ \mathbf{v}(\mathbf{v} \cdot \nabla t_{r} ) -\nabla t_{r} (\mathbf{v} \cdot \mathbf{v}) \right] } \nonumber \\[1em] &=\mathbf{v} + \nabla t_{r}(\bcR \cdot \mathbf{a}) -\nabla t_{r}( \mathbf{v} \cdot \mathbf{v}) \nonumber \\[1em] &= \mathbf{v}+( \bcR \cdot \mathbf{a} -v^2) \nabla t_{r} \end{align} $$

2番目の等号はBAC-CAB公式により成立する。

結論

$\cR=c(t-t_{r})$なので$\nabla \cR = -c\nabla t_{r}$である。これと$(4)$を$(1)$に代入すると

$$ \begin{align} \nabla V &= \frac{qc}{4\pi\epsilon_{0}} \frac{-1}{ (\cR c -\bcR \cdot \mathbf{v} )^2}\nabla (\cR c -\bcR \cdot \mathbf{v} ) \nonumber \\[1em] &= \frac{qc}{4\pi\epsilon_{0}} \frac{-1}{ (\cR c -\bcR \cdot \mathbf{v} )^2} \Big[ c\nabla \cR -\nabla (\bcR \cdot \mathbf{v} ) \Big] \nonumber \\[1em] &=\frac{qc}{4\pi\epsilon_{0}} \frac{-1}{ (\cR c -\bcR \cdot \mathbf{v} )^2} \Big[ -c^2\nabla t_{r} -\mathbf{v}-(\bcR \cdot \mathbf{a} -v^2)\nabla t_{r} \Big] \nonumber \\[1em] &=\frac{qc}{4\pi\epsilon_{0}} \frac{1}{ (\cR c -\bcR \cdot \mathbf{v} )^2} \Big[ \mathbf{v} + (c^2 -v^2+\bcR \cdot \mathbf{a} )\nabla t_{r}\Big] \end{align} $$

そして遅れ時間の勾配は下の通りだ。

$$ \nabla t_{r} = \frac {-\bcR }{\cR c -\bcR \cdot \mathbf{v} } $$

これを$(5)$に代入すると

$$ \begin{align*} \nabla V &= \frac{qc}{4\pi\epsilon_{0}} \frac{1}{ (\cR c -\bcR \cdot \mathbf{v} )^2} \left[ \mathbf{v} + (c^2 -v^2+\bcR \cdot \mathbf{a} )\nabla t_{r}\right] \\ &= \frac{qc}{4\pi\epsilon_{0}} \frac{1}{ (\cR c -\bcR \cdot \mathbf{v} )^2} \left[ \mathbf{v} + (c^2 -v^2+\bcR \cdot \mathbf{a} ) \frac {-\bcR }{\cR c -\bcR \cdot \mathbf{v} } \right] \\ &= \frac{qc}{4\pi\epsilon_{0}} \frac{1}{ (\cR c -\bcR \cdot \mathbf{v} )^3} \Big[ (\cR c -\bcR \cdot \mathbf{v})\mathbf{v} - (c^2 -v^2+\bcR \cdot \mathbf{a} ) \bcR \Big] \end{align*} $$

最後に$\dfrac{\partial \mathbf{A}}{\partial t}$を求めれば終わりだ。計算してみると結果は下の通りだ。

$$ \frac{ \partial \mathbf{A}}{ \partial t }=\frac{qc}{4\pi \epsilon_{0}} \frac{1}{(\cR c -\bcR \cdot \mathbf{v})^{3}}\left[ (\cR c +\bcR\cdot \mathbf{v})(\cR\mathbf{a}/c-\mathbf{v})+ \frac{\cR}{c}\mathbf{v}\left( c^{2} -v^{2}+\bcR\cdot \mathbf{a}\right) \right] $$

$\mathbf{u}\equiv c \crH-\mathbf{v}$として整理すると

$$ \begin{align*} \mathbf{E}(\mathbf{r},t) &= -\nabla V - \frac{ \partial \mathbf{A}}{ \partial t} \\ &= \frac{q}{4\pi\epsilon_{0}}\frac{\cR}{(\bcR\cdot \mathbf{u})^{3}}\left[ (c^{2}-v^{2})\mathbf{u} + \mathbf{u}(\bcR\cdot\mathbf{a})-\mathbf{a}(\bcR\cdot\mathbf{u}) \right] \\ &= \frac{q}{4\pi\epsilon_{0}}\frac{\cR}{(\bcR\cdot \mathbf{u})^{3}}\left[ (c^{2}-v^{2})\mathbf{u} + \bcR \times (\mathbf{u}\times \mathbf{a}) \right] \end{align*} $$

■