logo

動く点電荷が作る電場 📂電磁気学

動く点電荷が作る電場

概要 1

5D6FB6AA3.png

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

E(r,t)=q4πϵ0(⋅u)3[(c2−v2)u+×(u×a)]B(r,t)=1c×E(r,t) \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(r,t)=14πϵ0qc(c−⋅v),A(r,t)=vc2V(r,t) 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)

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

E=−∇V−∂A∂t,B=∇×A \mathbf{ E} = -\nabla V -\frac{\partial \mathbf{ A} }{\partial t},\quad \mathbf{B}=\nabla \times \mathbf{A}

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

∇V=qc4πϵ0∇1(c−⋅v) \nabla V=\frac{qc}{4 \pi \epsilon_{0}} \nabla \dfrac{1}{(\cR c - \bcR \cdot \mathbf{v}) }

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

∇V=qc4πϵ0−1(c−⋅v)2∇(c−⋅v) \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と∇(⋅v)\nabla (\bcR \cdot \mathbf{v})を計算すればいい。まず∇(⋅v)\nabla ( \bcR \cdot \mathbf{v} )から計算してみると乗法規則2によって

∇(⋅v)=(⋅∇)v+(v⋅∇)+×(∇×v)+v×(∇×) \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}

第一項

(⋅∇)v=(x∂∂x+y∂∂y+z∂∂z)v=xdvdtr∂tr∂x+ydvdtr∂tr∂y+zdvdtr∂tr∂z=dvdtr(x∂tr∂x+y∂tr∂y+z∂tr∂z)=a(⋅∇tr) \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*}

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

第二項

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

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

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

(A⋅∇)r=(Ax∂∂x+Ay∂∂y+Az∂∂z)(xx^+yy^+zz^)=Axx^+Ayy^+Azz^=A \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*}

そして

(A⋅∇)w(tr)=(Ax∂∂x+Ay∂∂y+Az∂∂z)w=Ax∂w∂x+Ay∂w∂y+Az∂w∂z=Ax∂w∂tr∂tr∂x+Ay∂w∂tr∂tr∂y+Az∂w∂tr∂tr∂z=Axv∂tr∂x+Ayv∂tr∂y+Azv∂tr∂z=v(Ax∂tr∂x+Ay∂tr∂y+Az∂tr∂z)=v(A⋅∇tr) \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)(3)は

(v⋅∇)=v−v(v⋅∇tr) (\mathbf{v} \cdot \nabla ) \bcR=\mathbf{v} -\mathbf{v}(\mathbf{v}\cdot \nabla t_{r} )

第三項

まず次が成り立つ。

∇×v=(∂vz∂y−∂vy∂z)x^+(∂vx∂z−∂vz∂x)y^+(∂vy∂x−∂vx∂y)z^=(∂vz∂trtr∂y−∂vy∂trtr∂z)x^+(∂vx∂trtr∂z−∂vz∂trtr∂x)y^+(∂vy∂trtr∂x−∂vx∂trtr∂y)z^=(az(∇tr)y−ay(∇tr)z)x^+(ax(∇tr)z−az(∇tr)x)y^+(ay(∇tr)x−ax(∇tr)y)z^=−a×∇tr \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*}

従って

×(∇×v)=−×(a×∇tr) \bcR \times ( \nabla \times \mathbf{v} )= -\bcR\times (\mathbf{a} \times \nabla t_{r})

第四項

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

∇×=∇×r−∇×w=−∇×w \nabla \times \bcR = \nabla \times \mathbf{r} -\nabla \times \mathbf{w}=-\nabla \times \mathbf{w}

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

v×(∇×)=v×(−∇×w)=v×(v×∇tr) \mathbf{v}\times (\nabla \times \bcR)=\mathbf{v} \times (- \nabla \times \mathbf{w})=\mathbf{v} \times (\mathbf{v} \times \nabla t_{r})

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

∇(⋅v)=a(⋅∇tr)+v−v(v⋅∇tr)−×(a×∇tr)+v×(v×∇tr)=a(⋅∇tr)+v−v(v⋅∇tr)−[a(⋅∇tr)−∇tr(⋅a)]+[v(v⋅∇tr)−∇tr(v⋅v)]=v+∇tr(⋅a)−∇tr(v⋅v)=v+(⋅a−v2)∇tr \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公式により成立する。

結論

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

∇V=qc4πϵ0−1(c−⋅v)2∇(c−⋅v)=qc4πϵ0−1(c−⋅v)2[c∇−∇(⋅v)]=qc4πϵ0−1(c−⋅v)2[−c2∇tr−v−(⋅a−v2)∇tr]=qc4πϵ01(c−⋅v)2[v+(c2−v2+⋅a)∇tr] \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}

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

∇tr=−c−⋅v \nabla t_{r} = \frac {-\bcR }{\cR c -\bcR \cdot \mathbf{v} }

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

∇V=qc4πϵ01(c−⋅v)2[v+(c2−v2+⋅a)∇tr]=qc4πϵ01(c−⋅v)2[v+(c2−v2+⋅a)−c−⋅v]=qc4πϵ01(c−⋅v)3[(c−⋅v)v−(c2−v2+⋅a)] \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*}

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

∂A∂t=qc4πϵ01(c−⋅v)3[(c+⋅v)(a/c−v)+cv(c2−v2+⋅a)] \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]

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

E(r,t)=−∇V−∂A∂t=q4πϵ0(⋅u)3[(c2−v2)u+u(⋅a)−a(⋅u)]=q4πϵ0(⋅u)3[(c2−v2)u+×(u×a)] \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*}

■

  1. David J. Griffiths, 기초전자기학(Introduction to Electrodynamics, 김진승 역) (4th Edition1 2014), p494-498 ↩︎