logo

ゼメンコ方程式 📂電磁気学

ゼメンコ方程式

概要1

連続的な電荷分布が時間によって変化する場合の電場は次の通りである。

E(r,t)=14πϵ0[ρ(r,tr)2+ρ˙(r,tr)cJ˙(r,tr)c2]dτ \mathbf{E} (\mathbf{r},t)=\frac{1}{4\pi \epsilon_{0}} \int \left[ \frac{ \rho (\mathbf{r}^{\prime}, t_{r}) }{\cR ^2} \crH + \frac{ \dot{\rho}(\mathbf{r}^{\prime}, t_{r})}{c\cR}\crH-\frac{\dot{\mathbf{J}}(\mathbf{r}^{\prime},t_{r}) }{c^2 \cR} \right]d\tau^{\prime}

連続的な電流分布が時間によって変化する場合の磁場は次の通りである。

B(r,t)=μ04π[J(r,tr)2+J˙(r,tr)c]×dτ \mathbf{B}( \mathbf{r}, t) = \dfrac{\mu_{0}}{4\pi} \int \left[ \frac{\mathbf{J}(\mathbf{r}^{\prime},t_{r})}{\cR^2} + \dfrac{ \dot{\mathbf{J}}(\mathbf{r}^{\prime}, t_{r}) } {c\cR} \right]\times \crH d\tau^{\prime}

これら二つの式を合わせてジェフィメンコ方程式と呼ぶ。ここで、trt_{r}遅延時間で、\bcR分離ベクトルである。

導出

電場と磁場は下記の式によって計算することができる。

E=VAtB=×A \begin{align} \mathbf{E} &= -\nabla V-\frac{\partial \mathbf{A}}{\partial t} \\ \mathbf{B} &= \nabla \times \mathbf{A} \end{align}

この場合、VVA\mathbf{A}は時間によって変化する遅延ポテンシャルであり、下記の通りである。

V(r, t)=14πϵ0ρ(r, tr)dτ,A(r, t)=μ04πJ(r, tr)dτ V(\mathbf{r},\ t)=\dfrac{1}{4\pi\epsilon_{0}} \int \dfrac{ \rho (\mathbf{r}^{\prime},\ t_{r}) }{ \cR } d\tau^{\prime},\quad \mathbf{A}( \mathbf{r},\ t) = \dfrac{\mu_{0}}{4\pi} \int \dfrac{\mathbf{J}(\mathbf{r}^{\prime},\ t_{r})}{\cR}d\tau^{\prime}

\cRtrt_{r}r\mathbf{r}^{\prime}を含んでいるため、計算はそう簡単ではない。VVの勾配は、連鎖法則により、遅延時間の勾配tr=1c\nabla t_{r}=-\dfrac{1}{c} \crHで、分離ベクトルの大きさの勾配1=12\dfrac{1}{\cR} = -\dfrac{1}{\cR^{2}}\crHだから、

V=14πϵ0(ρ(r, tr))dτ=14πϵ0[1ρ(r, tr)+ρ(r, tr)(1)]dτ=14πϵ0[1ρ(r,tr)trtrρ(r,tr)2]dτ=14πϵ0[ρ˙(r,tr)cρ(r,tr)2]dτ \begin{align} \nabla V &= \dfrac{1}{4\pi\epsilon_{0}} \int \nabla \left( \dfrac{ \rho (\mathbf{r}^{\prime},\ t_{r}) }{ \cR } \right) d\tau^{\prime} \nonumber \\ &= \dfrac{1}{4\pi\epsilon_{0}} \int \left[ \dfrac{1}{\cR}\nabla \rho (\mathbf{r}^{\prime},\ t_{r}) + \rho (\mathbf{r}^{\prime},\ t_{r}) \nabla \left( \dfrac{1}{\cR} \right) \right] d\tau^{\prime} \nonumber \\ &= \frac{1}{4\pi\epsilon_{0}} \int \left[ -\dfrac{1}{\cR}\dfrac{\partial \rho (\mathbf{r}^{\prime}, t_{r})}{\partial t_{r}} \nabla t_{r} -\rho (\mathbf{r}^{\prime},t_{r}) \dfrac { \crH} {\cR ^2} \right] d\tau^{\prime} \nonumber \\ &= \frac{1}{4\pi\epsilon_{0}} \int \left[ -\dfrac{\dot{\rho}(\mathbf{r}^{\prime}, t_{r})}{c \cR} {\crH} - \dfrac{\rho (\mathbf{r}^{\prime},t_{r})}{\cR ^2} \crH \right] d\tau^{\prime} \\ \end{align}

そして、ベクトル遅延ポテンシャルの時間微分を求めると、

At=tμ04πJ(r, tr)dτ=μ04π1J(r, tr)tdτ=μ04πJ˙(r, tr)dτ \begin{align} \dfrac{\partial \mathbf{A}}{ \partial t} &= \dfrac{\partial}{\partial t} \dfrac{\mu_{0}}{4\pi} \int \dfrac{\mathbf{J}(\mathbf{r}^{\prime},\ t_{r})}{\cR}d\tau^{\prime} \nonumber \\ &= \dfrac{\mu_{0}}{4\pi} \int \dfrac{1}{\cR} \dfrac{\partial \mathbf{J}(\mathbf{r}^{\prime},\ t_{r})}{\partial t} d\tau^{\prime} \nonumber \\ &= \dfrac{\mu_{0}}{4\pi} \int \dfrac{\dot{ \mathbf{J}}(\mathbf{r}^{\prime},\ t_{r})}{\cR} d\tau^{\prime} \end{align}

時間に対する変数と空間に対する変数は互いに独立しているので、積分の順序を変えることができる。(3)(3)(4)(4)(1)(1)に代入すると、

E(r,t)=14πϵ0[ρ(r,tr)2+ρ˙(r,tr)c]dτμ04πJ˙(r, tr)dτ=14πϵ0[ρ(r,tr)2+ρ˙(r,tr)cμ0ϵ0J˙(r, tr)]dτ=14πϵ0[ρ(r,tr)2+ρ˙(r,tr)cJ˙(r, tr)c2]dτ \begin{align*} \mathbf{E} (\mathbf{r},t) &= \dfrac{1}{4\pi\epsilon_{0}} \int \left[ \dfrac { \rho (\mathbf{r}^{\prime},t_{r})} {\cR ^2}\crH +\dfrac{\dot{\rho}(\mathbf{r}^{\prime}, t_{r})}{c\cR} \crH \right]d\tau^{\prime} -\dfrac{\mu_{0}}{4\pi} \int \dfrac{\dot{ \mathbf{J}}(\mathbf{r}^{\prime},\ t_{r})}{\cR} d\tau^{\prime} \\ &= \dfrac{1}{4\pi\epsilon_{0}} \int \left[ \dfrac { \rho (\mathbf{r}^{\prime},t_{r})} {\cR ^2}\crH +\dfrac{\dot{\rho}(\mathbf{r}^{\prime}, t_{r})}{c\cR} \crH -\mu_{0}\epsilon_{0}\dfrac{\dot{ \mathbf{J}}(\mathbf{r}^{\prime},\ t_{r})}{\cR} \right]d\tau^{\prime} \\ &= \dfrac{1}{4\pi\epsilon_{0}} \int \left[ \dfrac { \rho (\mathbf{r}^{\prime},t_{r})} {\cR ^2}\crH +\dfrac{\dot{\rho}(\mathbf{r}^{\prime}, t_{r})}{c\cR} \crH -\dfrac{\dot{ \mathbf{J}}(\mathbf{r}^{\prime},\ t_{r})}{c^2\cR} \right]d\tau^{\prime} \end{align*}

最後の等号は1c2=μ0ϵ0\dfrac{1}{c^2}=\mu_{0} \epsilon_{0}によって成立する。これは、電荷密度が時間によって変化する場合のクーロンの法則である。もし電荷分布が時間に対して一定ならば、静電学で学んだクーロンの法則と同じである。

磁場B\mathbf{B}は、回転演算子が含まれており、計算がより複雑である。

×A=μ04π[×(J(r, tr)dτ)] \nabla \times \mathbf{A} = \dfrac{\mu_{0}}{4\pi} \left[ \nabla \times \left( \int \dfrac{\mathbf{J}(\mathbf{r}^{\prime},\ t_{r})}{\cR} d\tau^{\prime} \right)\right]

この場合、×\nabla \times(x,y,z)(\mathbf{x},\mathbf{y},\mathbf{z})に対する微分であり、dτ\int d\tau^{\prime}(x,y,z)(\mathbf{x}^{\prime},\mathbf{y}^{\prime},\mathbf{z}^{\prime})に対する積分であるため、互いに独立している。したがって、順序を入れ替えても良い。

×A=μ04π×J(r, tr)dτ=μ04π[1(×J(r, tr))J(r, tr)×(1)]dτ \begin{align} \nabla \times \mathbf{A} &= \dfrac{\mu_{0}}{4\pi} \int \nabla \times \dfrac{ \mathbf{J}(\mathbf{r}^{\prime},\ t_{r})}{\cR} d\tau^{\prime} \nonumber \\ &= \dfrac{\mu_{0}}{4\pi} \int \left[ \dfrac{1}{\cR}\big( \nabla \times \mathbf{J}(\mathbf{r}^{\prime},\ t_{r})\big) - \mathbf{J}(\mathbf{r}^{\prime},\ t_{r}) \times \nabla \left(\dfrac{1}{\cR}\right) \right]d\tau^{\prime} \end{align}

最後の等号は積の法則 (e) ×(fA)=f(×A)A×(f)\nabla \times (f\mathbf{A}) = f(\nabla \times \mathbf{A}) - \mathbf{A} \times (\nabla f) によって成立する。まず、×J\nabla \times \mathbf{J}の各成分を求めてみよう。

(×J)x=JzyJyz (\nabla \times \mathbf{J})_{x}=\dfrac{\partial J_{z}}{\partial y}-\dfrac{\partial J_{y}}{\partial z}

微分の連鎖法則によって、

Jzy=Jztttrtry \dfrac{\partial J_{z}}{\partial y}=\dfrac{\partial J_{z}}{\partial t} \dfrac{\partial t}{\partial t_{r}}\dfrac{\partial t_{r}}{\partial y}

ここで、tr=tct_{r}=t-\dfrac{\cR}{c}なので ttr=1\dfrac{\partial t}{\partial t_{r}}=1try=1cy\dfrac{\partial t_{r}}{\partial y}=-\dfrac{1}{c}\dfrac{\partial \cR}{\partial y}なので、

Jzy=1cJz˙y \dfrac{\partial J_{z}}{\partial y}=-\dfrac{1}{c} \dot{J_{z}}\dfrac{\partial \cR}{\partial y}

結果的に、

(×J)x=1cJz˙y1cJy˙z=1c(Jy˙zJz˙y)=1c[J˙×()]x \begin{align*} (\nabla \times \mathbf{J})_{x} &=-\dfrac{1}{c} \dot{J_{z}}\dfrac{\partial \cR}{\partial y}-\dfrac{1}{c} \dot{J_{y}}\dfrac{\partial \cR}{\partial z} \\ &= \frac{1}{c} \left( \dot{J_{y}}\dfrac{\partial \cR}{\partial z}-\dot{J_{z}}\dfrac{\partial \cR}{\partial y}\right) \\ &= \frac{1}{c} \left[ \dot{\mathbf{J}} \times (\nabla \cR) \right]_{x} \end{align*}

さらに、()=\nabla ( \acR)=\acrHなので、

×J=1c(J˙×) \begin{equation} \nabla \times \mathbf{J} = \frac{1}{c} \left( \dot {\mathbf{J}} \times \crH \right) \end{equation}

そして、(1)=12\nabla \left( \frac{1}{\acR} \right)=-\frac{1}{\acR^2}\acrHなので、(6)(6)と共に(5)(5)に代入すると、

B(r,t)=×A=μ04π[11c(J˙×)+J×12]dτ=μ04π[J(r, tr)2+J˙(r, tr)c]×dτ \begin{align*} \mathbf{B}(\mathbf{r},t)=\nabla \times \mathbf{A} &= \dfrac{\mu_{0}}{4\pi} \int \left[ \dfrac{1}{\cR} \frac{1}{c} \left( \dot {\mathbf{J}} \times \crH \right) + \mathbf{J} \times \frac{1}{\cR^2}\crH\right]d\tau^{\prime} \\ &= \dfrac{\mu_{0}}{4\pi} \int \left[\frac{ \mathbf{J} (\mathbf{r}^{\prime},\ t_{r}) }{\cR^2} + \frac{ \dot {\mathbf{J}} (\mathbf{r}^{\prime},\ t_{r}) }{c\cR} \right]\times \crH d\tau^{\prime} \end{align*}


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