logo

リエナール-ヴィーヘルトポテンシャルの時間微分 📂電磁気学

リエナール-ヴィーヘルトポテンシャルの時間微分

概要

5E561AEA0.png

リエナール-ビケールトポテンシャルの時間に対する導関数は以下の通りである。

Vt=qc4πϵ01(cv)2(c2c2ttrv2+a)trtAt=qc4πϵ01(cv)3[(c+v)(a/cv)+cv(c2v2+a)] \begin{align*} \frac{ \partial V}{ \partial t} &= \frac{qc}{4\pi \epsilon_{0}} \frac{1}{(\cR c -\bcR \cdot \mathbf{v})^{2}} \left( c^{2} -c^{2}\frac{ \partial t}{ \partial t_{r}}-v^{2}+\bcR\cdot \mathbf{a} \right)\frac{ \partial t_{r}}{ \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] \end{align*}

補題

遅延時間の時間微分は以下の通りである。

trt=ccv=cu \frac{ \partial t_{r}}{ \partial t}=\frac{\cR c}{\cR c-\bcR\cdot \mathbf{v}}=\frac{ \cR c}{\bcR \cdot \mathbf{u}}

この時、u=cv\mathbf{u}=c\crH-\mathbf{v}である。

証明

遅延時間の定義によると、

c(ttr)= c(t-t_{r})=\cR

従って、両辺を二乗すると、

c2(ttr)2=2= c^{2}(t-t_{r})^{2}=\cR ^{2}=\bcR \cdot \bcR

両辺をttで微分すると、

2c2(ttr)(1trt)=2t(1) 2c^{2}(t-t_{r}) \left( 1-\frac{ \partial t_{r}}{ \partial t } \right)=2\frac{ \partial \bcR}{ \partial t }\cdot \bcR \tag{1}

この時、=rw(tr)\bcR=\mathbf{r}-\mathbf{w}(t_{r})であるので、

t=wt=wtrtrt=vtrt \frac{ \partial \bcR}{ \partial t}=-\frac{ \partial \mathbf{w}}{ \partial t}=-\frac{ \partial \mathbf{w}}{ \partial t_{r}}\frac{ \partial t_{r}}{ \partial t}=-\mathbf{v}\frac{ \partial t_{r} }{ \partial t}

である。これを(1)(1)に代入して整理すると、

c(1trt)=vtrt    cctrt=vtrt    trt=ccv=cu \begin{align*} && c\cR \left( 1-\frac{ \partial t_{r}}{ \partial t } \right) &= -\bcR\cdot \mathbf{v}\frac{ \partial t_{r}}{ \partial t } \\[1em] \implies && c\cR -c\cR \frac{ \partial t_{r}}{ \partial t} &= -\bcR\cdot\mathbf{v}\frac{ \partial t_{r} }{ \partial t } \\[1em] \implies && \frac{ \partial t_{r}}{ \partial t } &= \frac{c\cR}{c\cR -\bcR\cdot \mathbf{v}} \\[1em] && &= \frac{c\cR}{\bcR \cdot \mathbf{u}} \end{align*}

証明

リエナール-ビケールトポテンシャル

遅延時間 trt_{r}で速度v\mathbf{v}で移動する点電荷に対するポテンシャルは、以下の通りである。

V(r,t)=14πϵ0qc(cv)A(r,t)=μ04πqcv(cv)=vc2V(r,t) \begin{align*} V(\mathbf{r}, t) &= \frac{1}{4\pi \epsilon_{0}} \frac{qc}{ (\cR c -\bcR\cdot \mathbf{v})} \\ \mathbf{A}(\mathbf{r}, t) &= \frac{\mu_{0}}{4 \pi}\frac{qc \mathbf{v} }{(\cR c - \bcR\cdot \mathbf{v} )}=\frac{\mathbf{v}}{c^2}V(\mathbf{r}, t) \end{align*}

この時、=rw(tr)\bcR=\mathbf{r} -\mathbf{w}(t_{r})は遅延位置から観測点までのベクトルであり、w(tr)\mathbf{w}(t_{r})は遅延時間での点電荷の位置、つまり遅延位置である。

  • パート1. VV

スカラーポテンシャル VVttで微分すると、

Vt=qc4πϵ0t(1cv)=qc4πϵ0tr(1cv)trt \begin{align*} \frac{ \partial V}{ \partial t } &= \frac{qc}{4\pi \epsilon_{0}}\frac{\partial}{\partial t} \left( \frac{1}{\cR c -\bcR\cdot \mathbf{v}} \right) \\ &= \frac{qc}{4\pi \epsilon_{0}}\frac{\partial}{\partial t_{r}} \left(\frac{1}{\cR c -\bcR\cdot \mathbf{v}}\right) \frac{ \partial t_{r} }{ \partial t } \end{align*}

この時、ddx(1f(x))=ddf(x)(1f(x))df(x)dx=1(f(x))2f(x)\dfrac{ d }{ dx }\left( \dfrac{1}{f(x)} \right)=\dfrac{ d }{ d f(x) }\left(\dfrac{1}{f(x)} \right)\dfrac{ d f(x)}{ d x }=\dfrac{-1}{\left(f(x)\right)^{2}}f^{\prime}(x)であるので、

Vt=qc4πϵ01(cv)2(cv)trtrt=qc4πϵ01(cv)2(ctrtrvvtr)trt \begin{align*} \frac{ \partial V}{ \partial t }&=\frac{qc}{4\pi \epsilon_{0}} \frac{-1}{(\cR c -\bcR \cdot \mathbf{v})^{2}}\frac{ \partial (\cR c -\bcR \cdot \mathbf{v})}{ \partial t_{r} }\frac{ \partial t_{r}}{ \partial t } \\ &= \frac{qc}{4\pi \epsilon_{0}} \frac{-1}{(\cR c -\bcR \cdot \mathbf{v})^{2}} \left( c\frac{ \partial \cR}{ \partial t_{r}} - \frac{ \partial \bcR}{ \partial t_{r} }\cdot \mathbf{v}-\bcR\cdot \frac{ \partial \mathbf{v}}{ \partial t_{r} }\right)\frac{ \partial t_{r}}{ \partial t } \end{align*}

=rw(tr)\bcR=\mathbf{r}-\mathbf{w}(t_{r})であるので、

tr=w(tr)tr=v(tr) \frac{ \partial \bcR}{ \partial t_{r}}=-\frac{ \partial \mathbf{w}(t_{r})}{ \partial t_{r} }=-\mathbf{v}(t_{r})

==c(ttr)| \bcR|=\cR=c(t-t_{r})であるので、

tr=cttrc \frac{ \partial \cR}{ \partial t_{r}}=c\frac{ \partial t}{ \partial t_{r} }-c

従って、

Vt=qc4πϵ01(cv)2(c2+c2ttr+v2a)trt=qc4πϵ01(cv)2(c2c2ttrv2+a)trt \begin{align*} \frac{ \partial V}{ \partial t} &=\frac{qc}{4\pi \epsilon_{0}} \frac{-1}{(\cR c -\bcR \cdot \mathbf{v})^{2}} \left( -c^{2} +c^{2}\frac{ \partial t}{ \partial t_{r}}+v^{2}-\bcR\cdot \mathbf{a} \right)\frac{ \partial t_{r}}{ \partial t } \\ &= \frac{qc}{4\pi \epsilon_{0}} \frac{1}{(\cR c -\bcR \cdot \mathbf{v})^{2}} \left( c^{2} -c^{2}\frac{ \partial t}{ \partial t_{r}}-v^{2}+\bcR\cdot \mathbf{a} \right)\frac{ \partial t_{r}}{ \partial t } \end{align*}

  • パート2. A\mathbf{A}

補題と(a)(a)を使ってうまく整理すると、容易に得られる。

At=t(vc2V(r,t))=1c2(vtV+vVt)=1c2[vtrtrtV+vqc4πϵ01(cv)2(c2c2ttrv2+a)trt]=1c2trt[aqc4πϵ01cv+vqc4πϵ01(cv)2(c2c2ttrv2+a)]=1c2trtqc4πϵ01(cv)2[a(c+v)+v(c2c2ttrv2+a)] \begin{align*} \frac{ \partial \mathbf{A}}{ \partial t } &= \frac{ \partial }{ \partial t }\left( \frac{ \mathbf{v}}{c^{2}}V(\mathbf{r},t) \right) \\[1em] &= \frac{1}{c^{2}} \left( \frac{ \partial \mathbf{v}}{ \partial t }V +\mathbf{v}\frac{ \partial V }{ \partial t} \right) \\[1em] &= \frac{1}{c^{2}} \left[ \frac{ \partial \mathbf{v}}{ \partial t_{r} } \frac{ \partial t_{r}}{ \partial t }V + \mathbf{v}\frac{qc}{4\pi \epsilon_{0}} \frac{1}{(\cR c -\bcR \cdot \mathbf{v})^{2}} \left( c^{2} -c^{2}\frac{ \partial t}{ \partial t_{r}}-v^{2}+\bcR\cdot \mathbf{a} \right)\frac{ \partial t_{r}}{ \partial t }\right] \\[1em] &= \frac{1}{c^{2}}\frac{\partial t_{r}}{\partial t}\left[ \mathbf{a}\frac{qc}{4\pi \epsilon_{0}} \frac{1}{\cR c- \bcR \cdot \mathbf{v}}+ \mathbf{v}\frac{qc}{4\pi \epsilon_{0}} \frac{1}{(\cR c -\bcR \cdot \mathbf{v})^{2}} \left( c^{2} -c^{2}\frac{ \partial t}{ \partial t_{r}}-v^{2}+\bcR\cdot \mathbf{a} \right)\right] \\[1em] &= \frac{1}{c^{2}}\frac{\partial t_{r}}{\partial t}\frac{qc}{4\pi \epsilon_{0}} \frac{1}{(\cR c -\bcR \cdot \mathbf{v})^{2}}\left[ \mathbf{a}(\cR c +\bcR\cdot \mathbf{v})+ \mathbf{v} \left( c^{2} -c^{2}\frac{ \partial t}{ \partial t_{r}}-v^{2}+\bcR\cdot \mathbf{a} \right)\right] \end{align*}

補題によれば、

At=1c2ccvqc4πϵ01(cv)2[a(c+v)+v(c2c2cvcv2+a)]=cqc4πϵ01(cv)3[a(c+v)+c2vc2v(cvc)v2v+(a)v]=cqc4πϵ01(cv)3[(c+v)(ac2vc)+v(c2v2+a)]=qc4πϵ01(cv)3[(c+v)(a/cv)+cv(c2v2+a)] \begin{align*} \frac{ \partial \mathbf{A}}{ \partial t }&= \frac{1}{c^{2}}\frac{\cR c}{\cR c-\bcR\cdot \mathbf{v}}\frac{qc}{4\pi \epsilon_{0}} \frac{1}{(\cR c -\bcR \cdot \mathbf{v})^{2}}\left[ \mathbf{a}(\cR c +\bcR\cdot \mathbf{v})+ \mathbf{v} \left( c^{2} -c^{2}\frac{\cR c-\bcR\cdot \mathbf{v}}{\cR c}-v^{2}+\bcR\cdot \mathbf{a} \right)\right] \\[1em] &=\frac{\cR}{c} \frac{qc}{4\pi \epsilon_{0}} \frac{1}{(\cR c -\bcR \cdot \mathbf{v})^{3}}\left[ \mathbf{a}(\cR c +\bcR\cdot \mathbf{v})+ c^{2}\mathbf{v} -c^{2}\mathbf{v}\left( \frac{\cR c-\bcR\cdot \mathbf{v}}{\cR c} \right)-v^{2 }\mathbf{v}+(\bcR\cdot \mathbf{a}) \mathbf{v}\right] \\[1em] &=\frac{\cR}{c} \frac{qc}{4\pi \epsilon_{0}} \frac{1}{(\cR c -\bcR \cdot \mathbf{v})^{3}}\left[ (\cR c +\bcR\cdot \mathbf{v})(\mathbf{a}-\frac{c^{2}\mathbf{v}}{\cR c})+ \mathbf{v}\left( c^{2} -v^{2}+\bcR\cdot \mathbf{a}\right) \right] \\[1em] &= \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] \end{align*}