Feynman QED Eighth Lecture

\(\)

\(\textit{Eighth Lecture}\)

自由空間に於けるマックスウェル方程式の解

自由空間(すなわち \(\rho=0,\,\mathbf{j}=0\) の真空)に於ける波動方程式, つまり式 (2.7.21′) で \(\rho=0,\,\mathbf{j}=0\) とした式, または \(j^{\,\mu}=0\) と置いた式 (2.7.23) のマクスウェル方程式
\begin{align*}
\def\ppdiff#1#2{\frac{\partial #1}{\partial #2}}
\def\pdiff#1{\frac{\partial}{\partial #1}}
\def\Bppdiff#1#2{\frac{\partial^{2} #1}{\partial #2^{2}}}
\def\Bpdiff#1{\frac{\partial^{2}}{\partial #1^{2}}}
\def\mb#1{\mathbf{#1}}
\def\mr#1{\mathrm{#1}}
\def\reverse#1{\frac{1}{#1}}
\def\ds#1{\mbox{${\displaystyle\strut #1}$}}
\nabla^{2} \phi -\frac{1}{c^{2}}\Bppdiff{\phi}{t}=0,\quad\mathrm{and}\quad
\nabla^{2}\mb{A}-\frac{1}{c^{2}}\Bppdiff{\mb{A}}{t}=0,
\quad\rightarrow\quad \square^{2}A^{\,\mu}=0
\tag{2.8.1}
\end{align*}
の平面波解は [1][訳註] (ランダウ:「力学・場の力学」より) ポテンシャルの持つ任意性のため, … Continue reading,
\begin{align*}
A^{\,\mu}=\epsilon^{\,\mu}e^{-ik\cdot X}
\tag{2.8.2}
\end{align*}
ただし\(\epsilon^{\,\mu}\)と\(k^{\,\mu}\)は定数ベクトルであり, そして \(k^{\,\mu}\) は次の条件に従う:
\begin{align*}
k_{\mu}k^{\,\mu}=k\cdot k = 0
\tag{2.8.3}
\end{align*}
これは \(\partial_{\nu}\) が \(e^{-ik\cdot X}\) に作用すると \(-ik_{\nu}\) を掛け合わせる効果を持つことから理解できるであろう (直交座標系なので \(\partial_{\nu}\) は \(\epsilon^{\,\mu}\) には作用しない).従って,
\begin{align*}
\square^{2} A^{\,\mu} &= \partial^{\,\nu}(\partial_{\nu}\,A^{\,\mu})
=\partial^{\,\nu}\bigl(\partial_{\nu}\,\epsilon^{\,\mu}e^{-ik^{\,\nu}X_{\nu}}\bigr)\\
&=\partial^{\,\nu}\bigl(\epsilon^{\,\mu}(-ik_{\nu})\,e^{-ik_{\nu}X^{\,\nu}}\bigr)
=\epsilon^{\,\mu}\,(-ik_{\nu})\,\bigl(\partial^{\,\nu}\,e^{-ik^{\,\nu}X_{\nu}}\bigr)\\
&=(-ik_{\nu})(-ik^{\,\nu})\epsilon^{\,\mu}\,e^{-ik\cdot X}
=-(k_{\nu}k^{\,\nu})\,A^{\,\mu}=0
\end{align*}
この演算に於いて\(\partial_{\nu}A^{\,\mu}\)は, 実際は「2階テンソル」, そして \(\partial^{\,\nu}(\partial_{\nu}A^{\,\mu})\) は「3階テンソル」を形成しており, 従って添字 \(\nu\) の縮約により「1階テンソル」即ち「ベクトル」となることに注意する.

\(k^{\,\mu}\) は「伝搬ベクトル」(propagation vector) であって, 次の成分を持つ:
\begin{align*}
k^{\,\mu}=\left(\frac{\omega}{c},K_x,K_y,K_z\right)=\left(\frac{\omega}{c},\mb{K}\right)
\tag{2.8.4}
\end{align*}
従って, 通常の表記では [2][訳註] 前講から4元座標は \(\mb{X}=(ct,\,\mb{X})\) であるから, \begin{align*} k\cdot X =k^{\,\mu}X_{\mu}=\left(\frac{\omega}{c},\,\mb{K}\right)\cdot \left(ct,\,-\mb{X}\right)=\omega t … Continue reading
\begin{align*}
\exp(-ik\cdot X)=\exp\bigl\{-i\bigl(\omega t-\mb{K}\cdot\mb{X}\bigr)\bigr\}
\tag{2.8.5}
\end{align*}
そして, 条件 \(k\cdot k=0\) は次を意味する:
\begin{align*}
k\cdot k = k^{\,\mu}k_{\mu}=\left(\frac{\omega}{c},\mb{K}\right)\cdot\left(\frac{\omega}{c},-\mb{K}\right)
=\frac{\omega^{2}}{c^{2}}-\mb{K}\cdot\mb{K}=0
\tag{2.8.6}
\end{align*}


【 問題 】 ローレンツ条件: \(\partial^{\,\mu}A_{\mu}=0\) は, \(k\cdot\epsilon=0\) を含意していることを示せ.

(解答例) ローレンツ条件式 \(\partial^{\,\mu}A_{\mu}=0\) に \(A_{\mu}=\epsilon_{\mu}\,e^{-ik\cdot X}\) を代入してみると,

\begin{equation*}
\partial^{\,\mu}A_{\mu}=\partial^{\,\mu}\bigl(\epsilon_{\mu}\,e^{-ik\cdot X}\bigr)
=\epsilon_{\mu}\bigl(\partial^{\,\mu}\,e^{-ik^{\,\mu}\,X_{\mu}}\bigr)
=\epsilon_{\mu}(-ik^{\,\mu})\,e^{-ik\cdot X}=-i(\epsilon_{\,\mu}k^{\,\mu})e^{-ik\cdot X}=0
\end{equation*}

これが満たされるためには \(\epsilon_{\,\mu}k^{\,\mu}=0\) すなわち \(\epsilon\cdot k=0\) が成り立っていれば良い.


3次元で作業する場合, 偏極ベクトル \(\mb{\epsilon}\) は \(\mb{K}\cdot\mb{\epsilon}=0\) となるように取り, そして「スカラーポテンシャルは \(\phi=0\) とする」 のが通例である.しかしそれが唯一の条件ではない.つまり, これは相対論的に不変ではなく, ある一つの座標系に於いてのみ成り立つ.これは \(\mb{K}\cdot\mb{e}=0\) の系に唯一性を与えるパラドックスであり, 相対性理論とは相容れない状況のように思われる.しかしこの「パラドックス」は, いわゆる「ゲージ変換は常に行うことが出来る」という事実によって解決される.その変換では, 場 \(F_{\mu\nu}\) は変化しないが \(\mb{\epsilon}\) は「変化する」.従って, 特定の系で \(\mb{K}\cdot\mb{\epsilon}=0\) を選択するのは「特定のゲージを選択すること」に等しい.(それを以下で示す).

ゲージ変換の式 (2.7.19) は次の2式に書ける:
\begin{align*}
&A’_{\mu}=A_{\mu}-\partial_{\mu}\,\chi,\quad\mathrm{where}\quad A_{\mu}=\bigl(\phi,-\mb{A}\bigr),\quad
\partial_{\mu}=\left(\frac{1}{c}\pdiff{t},\nabla \right),\\
&\Rightarrow\quad \phi’ =\phi -\frac{1}{c}\ppdiff{\chi}{t},\quad \mb{A}’=\mb{A}+\nabla \chi
\tag{2.8.7}
\end{align*}
ただし \(\chi\) はスカラーである.式 (2.7.20) のローレンツ条件 \(\partial^{\,\mu}A_{\mu}= 0\) の両辺に \(\partial^{\,\mu}\) を作用させると,
\begin{align*}
\partial^{\,\mu} A’_{\mu}=\partial^{\,\mu}\bigl(A_{\mu}-\partial_{\mu}\chi\bigr)
=\partial^{\,\mu}A_{\mu}-\partial^{\,\mu}\partial_{\mu}\chi
=\partial^{\,\mu}A_{\mu}-\square^{2}\chi =0
\tag{2.8.8}
\end{align*}
従って, もし
\begin{align*}
\square^{2}\chi=0
\tag{2.8.9}
\end{align*}
が成り立つならば,

\begin{equation*}
\partial^{\,\mu} A’_{\mu}=\partial^{\,\mu}A_{\mu}=0
\end{equation*}

となるから,「ゲージ変換をしてもローレンツ条件の式 (2.7.20) は依然として成立する」であろう.方程式 (2.8.9) は, 解として \(\chi=\alpha\,e^{-ik\cdot X}\) を持つ.ただし \(\alpha\) は任意定数である.従って,
\begin{align*}
A’_{\mu}=\epsilon’_{\mu}\,e^{-ik\cdot X}
&=A_{\mu}-\partial_{\mu}\,\chi=A_{\mu}-\partial_{\mu}(\alpha\,e^{-ik_{\mu}X^{\,\mu}})
=\epsilon_{\mu}\,e^{-ik\cdot X}-\alpha (-ik_{\mu})\,e^{-ik\cdot X}\\
&=(\epsilon_{\mu}+i\alpha\,k_{\mu})\,e^{-ik\cdot X}
\end{align*}
従って,
\begin{align*}
\epsilon’_\mu = \epsilon_{\mu}+ i\alpha k_{\mu}
\tag{2.8.9}
\end{align*}
は, ゲージ変換によって得られる新たな偏極ベクトルである.通常の表記では, これは次の2式となる:
\begin{align*}
\mb{\epsilon}’ = \mb{\epsilon}+i\alpha \mb{K},\quad \epsilon’_{0}=\epsilon_0 +i\alpha \frac{\omega}{c}
\tag{2.8.10}
\end{align*}
この両辺に \(\mb{K}\) を掛け合わせ, 式 (2.8.6) の \(\displaystyle \frac{\omega^{2}}{c^{2}}-\mb{K}\cdot\mb{K}=0\) を用いると,
\begin{align*}
\mb{K}\cdot\mb{\epsilon}’ = \mb{K}\cdot\mb{\epsilon}+i\alpha\,\mb{K}\cdot\mb{K}
=\mb{K}\cdot\mb{\epsilon}+i\alpha\frac{\omega^{2}}{c^{2}}
\tag{2.8.11}
\end{align*}
従って, どんな座標系を用いている場合でも, 定数 \(\alpha\) を適切に選択することで
\begin{equation*}
\mb{K}\cdot\mb{\epsilon’}=0
\end{equation*}

とすることが出来る.

場は明らかにゲージ変換を受けても変わらない.なぜなら, 微分の順序は重要ではないからである:
\begin{align*}
F’_{\mu\nu}&=\partial_{\mu} A’_{\nu}-\partial_{\nu} A’_{\mu}
=\partial_{\mu}\bigl(A_{\nu}-\partial_{\nu}\,\chi\bigr)-\partial_{\nu}
\bigl(A_{\mu}-\partial_{\mu}\,\chi\bigr)\\
&=\partial_{\mu}A_{\nu}-\partial_{\mu}\partial_{\nu}\,\chi
-\partial_{\nu}A_{\mu}+\partial_{\nu}\partial_{\mu}\,\chi
=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}-\partial_{\mu}\partial_{\nu}\,\chi
+\partial_{\nu}\partial_{\mu}\,\chi\\
&=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}-\partial_{\mu}\partial_{\nu}\,\chi
+\partial_{\mu}\partial_{\nu}\,\chi\\
&=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}= F_{\mu\nu}
\end{align*}

相対論的素粒子論(Relativistic Particle Mechanics)

通常の(3次元)速度」の成分は 4-ベクトルの成分のようには変換しない.しかし, 別の量
\begin{align*}
u^{\,\mu}\equiv \frac{dz^{\,\mu}}{d\tau}
=\left(\frac{cdt}{d\tau},\frac{dx}{d\tau},\frac{dy}{d\tau},\frac{dz}{d\tau}\right),
\ \mathrm{where}\ \ dz^{\,\mu}=(cdt,dx,dy,dz)
\tag{2.8.13}
\end{align*}
は, 粒子の経路要素である (この \(u^{\,\mu}\) が「4元速度」である).そして \(d\tau\) は次式で定義される「固有時間」(proper time) である:
\begin{align*}
ds^2=dt^{2} -dx^{2} -dy^{2} -dz^{2}\equiv c^2 d\tau^{2}
\tag{2.8.14}
\end{align*}
通常の速度 \(\mb{v}=(v_x,v_y,v_z)\) との関係は次式のようになる: [3][訳註] (ランダウ:「力学・場の理論」§35より) ある慣性基準系 \(K\) から速さ \(v\) で運動している時計(その慣性系を \(K’\) とする) … Continue reading
\begin{align*}
\frac{dx}{d\tau}=\frac{dx}{dt}\frac{dt}{d\tau}&=\frac{v_x}{\sqrt{1-(v/c)^{2}}},\qquad
\frac{dy}{d\tau}=\frac{v_y}{\sqrt{1-(v/c)^{2}}},\\
\frac{dz}{d\tau}&=\frac{v_z}{\sqrt{1-(v/c)^{2}}},\qquad
\frac{dt}{d\tau}=\frac{1}{\sqrt{1-(v/c)^{2}}}
\tag{2.8.15}
\end{align*}
従って \(u^{\,\mu}=(u^{\,0},u^{\,1},u^{\,2},u^{\,3})\) は,
\begin{align*}
u^{\,0}=c\frac{dt}{d\tau}=\frac{c}{\sqrt{1-(v/c)^{2}}},\quad
u^{\,k}=\frac{v^{\,k}}{\sqrt{1-(v/c)^{2}}}\quad (k=1,2,3)
\tag{2.8.16}
\end{align*}
\(u^{2}= u^{\,\mu}u_{\mu}=c^{2}\) であることは明らかである.なぜなら,
\begin{align*}
u^{\,\mu}u_{\mu}&=\left(\frac{cdt}{d\tau},\frac{dx}{d\tau},\frac{dy}{d\tau},\frac{dz}{d\tau}\right)
\cdot\left(\frac{cdt}{d\tau},-\frac{dx}{d\tau},-\frac{dy}{d\tau},-\frac{dz}{d\tau}\right)\\
&=c^{2}\left(\frac{dt}{d\tau}\right)^2-\left(\frac{dx}{d\tau}\right)^2
-\left(\frac{dy}{d\tau}\right)^2-\left(\frac{dz}{d\tau}\right)^2\notag\\
&=\frac{c^2}{1-(v/c)^{2}}-\frac{v_x^{2}}{1-(v/c)^{2}}-\frac{v_y^{2}}{1-(v/c)^{2}}
-\frac{v_z^{2}}{1-(v/c)^{2}}\\
&=c^2 \frac{1-(v/c)^{2}}{1-(v/c)^{2}}=c^2
\tag{2.8.17}
\end{align*}
4元運動量」 \(p^{\,\mu}\) は次式で定義される:
\begin{align*}
p^{\,\mu}&=mu^{\,\mu}=\bigl(p^0,\,p^1,\,p^2,\,p^3\bigr)
=\left(\frac{E}{c},\,p_x,\,p_y,\,p_z\right)\\
&=\left(\frac{mc}{\sqrt{1-(v/c)^{2}}},\frac{mv_x}{\sqrt{1-(v/c)^{2}}},\frac{mv_y}{\sqrt{1-(v/c)^{2}}},
\frac{mv_z}{\sqrt{1-(v/c)^{2}}}\right)
\tag{2.8.18}
\end{align*}
\(cp^0\) は「全エネルギー」\(E\) であることに注意する:
\begin{align*}
p^0=\frac{E}{c}=\frac{mc}{\sqrt{1-(v/c)^{2}}},\quad\rightarrow\quad
E=cp^0 = \frac{mc^2}{\sqrt{1-(v/c)^{2}}}=\gamma mc^2
\tag{2.8.19}
\end{align*}
従って, 通常の表記での「3次元運動量」\(\mb{p}\) は次で与えられる:
\begin{align*}
\frac{E}{c^2}=\frac{1}{\sqrt{1-\beta^2}}=\gamma m,\quad\rightarrow\quad
\mb{p}=\frac{m\mb{v}}{\sqrt{1-\beta^2}}=\gamma m\mb{v}=\frac{E}{c^2}\mb{v}
\tag{2.8.20}
\end{align*}
ここで \(\mb{v}\) は「通常の3次元速度」である.

速度と同様に, 通常の力の成分は \(d\mb{p}/dt\) で定義されるが, これは4元ベクトルの成分を形成することは出来ない.しかし次の量は4元ベクトルを成す:
\begin{align*}
f^{\,\mu}=\frac{d p^{\,\mu}}{d\tau}
\tag{2.8.21}
\end{align*}
この成分は次となる:
\begin{align*}
f^{\,k}&=\frac{d p^{\,k}}{d\tau}=\frac{d p^{\,k}}{dt}\cdot\frac{dt}{d\tau}
=F^{\,k}\frac{1}{\sqrt{1-\beta^2}}
=\frac{dt}{d\tau}\cdot\frac{d}{dt}\,p^{\,k}\\
&=\frac{1}{\sqrt{1-\beta^2}}\frac{d}{dt}\left(\frac{mv^{\,k}}{\sqrt{1-\beta^2}}\right)
\quad (k =1,2,3),\\
f^{\,0}&=\frac{d p^{\,0}}{d\tau}=\frac{d p^{\,0}}{dt}\frac{dt}{d\tau}
=\frac{d (E/c)}{dt}\frac{dt}{d\tau}=\frac{1}{c}\frac{dE}{dt}\frac{dt}{d\tau}\\
&=\frac{(\mathrm{Power})}{c}\frac{1}{\sqrt{1-\beta^2}}
=\frac{(\mathrm{rate\ of\ change\ of\ Energy})}{c}\cdot\gamma\\
&=\frac{dt}{d\tau}\cdot \frac{d}{dt}\,p^{\,0}
=\gamma\cdot \frac{d}{dt}\left(\frac{mc^2}{\sqrt{1-\beta^2}}\right)
\tag{2.8.22}
\end{align*}
ただし \(F^{\,k}\) は通常の力である:
\begin{align*}
F^{\,k}=\frac{d}{dt}\,p^{\,k}=\frac{d}{dt}\left(\frac{mv^{\,k}}{\sqrt{1-\beta^2}}\right)
\end{align*}
以上のことは, \(\displaystyle \frac{mc^2}{\sqrt{1-\beta^2}}\) が全エネルギー \(E\) であることから, または, 次の通常の恒等式からも分かることである:
\begin{align*}
(\mathrm{Power})=\mb{F}\cdot\mb{v}&=\left(\frac{d}{dt}\,\mb{p}\right)\cdot\mb{v}
=\mb{v}\cdot\left[\frac{d}{dt}\left(\frac{m\mb{v}}{\sqrt{1-\beta^2}}\right)\right]\\
&=\frac{\mb{v}}{\sqrt{1-\beta^2}}\cdot\frac{d }{dt}\bigl(m\mb{v}\bigr)
+m\,\mb{v}\cdot\mb{v}\,\left(\frac{d}{dt}\frac{1}{\sqrt{1-\beta^2}}\right)\\
&=\frac{m}{\sqrt{1-\beta^2}}\left(\mb{v}\cdot\frac{d\mb{v}}{dt}\right)
+mv^2\frac{1}{c^2(1-\beta^2)^{3/2}}v\frac{dv}{dt}
\tag{2.8.23}
\end{align*}
ここで, \(\mb{v}\cdot\mb{v}=v^2\) の両辺を時間微分すると
\begin{align*}
\frac{d}{dt}\bigl(\mb{v}\cdot\mb{v}\bigr)&=\frac{d\mb{v}}{dt}\cdot\mb{v}
+\mb{v}\cdot\frac{d\mb{v}}{dt}=2\mb{v}\cdot\frac{d\mb{v}}{dt}
=\frac{d}{dt}u^2=2u\frac{du}{dt},\\
\rightarrow\quad
&\mb{v}\cdot\frac{d\mb{v}}{dt}=v\frac{dv}{dt}
\tag{2.8.24}
\end{align*}
この関係を前式 (2.8.23) に代入すると,
\begin{align*}
\mathrm{(Power)} &= \frac{m}{\sqrt{1-\beta^2}}v\frac{dv}{dt}
+\frac{m\beta^{2}}{(1-\beta^2)\sqrt{1-\beta^2}}v\frac{dv}{dt}\\
&=\frac{m}{\sqrt{1-\beta^2}}\left(1-\frac{\beta^2}{1-\beta^2}\right)v\frac{dv}{dt}\\
&=\frac{mv}{\sqrt{1-\beta^2}}\frac{dv}{dt}
=\frac{d}{dt}\left(\frac{mc^2}{\sqrt{1-\beta^2}}\right)
\tag{2.8.25}
\end{align*}
従って, ニュートン方程式に類似した相対論的な式は次となる:[4][訳註] (C.メラー §27 より) 粒子の運動量が変化していれば, 粒子は単位時間当たりの運動量の変化に等しい力 \(\mb{F}\) … Continue reading
\begin{align*}
\frac{d}{d\tau}p^{\,\mu}=f^{\,\mu}=m\frac{d^2 z^{\,\mu}}{d\tau^{2}}
=m\frac{d u^{\,\mu}}{d\tau}
\tag{8-1}
\end{align*}
通常のローレンツ力は,
\begin{align*}
\mb{F}=e\mb{E}+\frac{e}{c}\mb{v}\times\mb{B}
\tag{2.8.26}
\end{align*}
そして, エネルギーの変化率(Power) は,
\begin{align*}
\mathrm{(Power)} = \mb{F}\cdot\mb{v}=e\mb{E}\cdot\mb{v}
\tag{2.8.27}
\end{align*}
すると, 式 (2.8.22) の \(f^{\,\mu}\) すなわち「4元力」(four-force) は,
\begin{align*}
\mb{f}&=\frac{\mb{F}}{\sqrt{1-\beta^2}}=\frac{e}{\sqrt{1-\beta^2}}
\bigl(\mb{E}+\frac{1}{c}\mb{v}\times\mb{B}\bigr),\\
f^{0}&=\frac{(\mathrm{Power})}{c}\frac{1}{\sqrt{1-\beta^2}}
=\frac{e\mb{E}\cdot\mb{v}}{c}\frac{1}{\sqrt{1-\beta^2}}
=\frac{e}{c\sqrt{1-\beta^2}}\mb{E}\cdot\mb{v}
\tag{2.8.28}
\end{align*}
となる.


【 問題 】 \(\mb{f}\) と \(f^{\,0}\) で与えられる表現式に等価なのは次式であることを示せ:
\begin{align*}
f_{\mu} = \frac{e}{c}F_{\mu\nu}\,u^{\,\nu}
\tag{2.8.29}
\end{align*}
従って, ニュートンの運動方程式に相当する相対論的な式は次となる:
\begin{align*}
m\frac{d^2 z^{\,\mu}}{d\tau^2}=\frac{e}{c}\frac{dz^{\,\nu}}{d\tau}\,F_{\mu\nu}
\tag{8-3}
\end{align*}
この式は, 次式が成り立つことを必要(前提)とする(imply):[5][訳註] 式(2.8.30)は, 原書では次となっている: \begin{align*} \frac{d}{d\tau}\left(\frac{dz^{\,\mu}}{d\tau}\right)^{2} =0 \end{align*} しかし, … Continue reading
\begin{align*}
\frac{d}{d\tau}\,u^2 = \frac{d}{d\tau}\,\bigl(u^{\,\nu}u_{\nu}\bigr)=0
\tag{2.8.30}
\end{align*}


〈解答例〉 ある慣性系 \(S\) に於いて速度 \(\mb{u}\) で与えられた電磁場内を運動する荷電粒子 \(e\) に作用する力を決めることにする.通常の3元ベクトルでは式 (2.8.27) のローレンツ力になる:
\begin{align*}
\mb{f}=e\mb{E}+\frac{e}{c}\mb{v}\times\mb{B}
\tag{1}
\end{align*}
これを踏まえて, この粒子が瞬間的に静止しているような慣性系 \(S_0\) を導入する.\(S_0\) 系では, 荷電粒子に作用する通常の力 \(\mb{f}’\) は, 上式 (1) に\(\mb{v}=0\) を代入して次となる:
\begin{align*}
\mb{f}’=e\mb{E}
\tag{2}
\end{align*}
他方, \(S_0\) 系での4元力 \(f_{\mu}\) は式 (2.8.29) に \(\mb{v}=0\) と \(\beta=v/c=0\) を代入して次となる:
\begin{align*}
f_{\mu}=\left(\frac{e\mb{E}\cdot\mb{v}}{c\sqrt{1-\beta^2}},\,
\frac{-e}{\sqrt{1-\beta^2}}\bigl(\mb{E}+\frac{1}{c}\mb{v}\times\mb{B}\bigr)\right)
=\bigl(0,\,-e\mb{E}\bigr)
\tag{3}
\end{align*}
以上のことを踏まえて, 問題文の式 (2.8.29) の4元量を考える.静止系 \(S’\) に於ける粒子の4元速度は, 式 (2.8.16) で \(\mb{v}=0\) 及び \(\beta=0\) として次である:
\begin{align*}
u^{\,\mu}=\left(\frac{c}{\sqrt{1-\beta^2}},\,\frac{\mb{v}}{\sqrt{1-\beta^2}}\right),
\quad\rightarrow\quad
u’^{\,\mu}=\bigl(u^{\,0},\,\mb{u}\bigr)=\bigl(c,\,0 \bigr)
\tag{4}
\end{align*}
よって \(u^{0}=c\), \(u^{\,k}=0\) であるから, これを式 (2.8.29) の4元量に代入してみると,
\begin{align*}
f’_{\mu} &= \frac{e}{c}F_{\mu\nu}\,u^{\,\nu}
=\frac{e}{c}\left(F_{\mu\,0}\,u^{\,0}+F_{\mu\,k}\,u^{\,k}\right)
=\frac{e}{c}u^{\,0}\,F_{\mu\,0}=\frac{e}{c}\,c\,F_{\mu\,0}= eF_{\mu\,0}
\tag{5}
\end{align*}
前節の結果式 (2.7.16) から, \(F_{\mu\,0}\) は
\begin{align*}
F_{0\,0}=0,\quad F_{0\,1}=-E_x \quad F_{0\,2}=-E_y \quad F_{0\,3}=-E_z
\tag{6}
\end{align*}
であったから, これらを式 (5) に用いれば
\begin{align*}
f’_{0}=eF_{0\,0}=0,\quad f’_1=eF_{1\,0}=-eE_x,\quad f’_2=eF_{2\,0}=-eE_y,
\quad f’_3=eF_{3\,0}=-eE_z
\tag{7}
\end{align*}
すなわち, 静止系 \(S’\) に於ける粒子の4元力は
\begin{align*}
f’_{\mu}=\bigl(f’_0,\,\mb{f}’\bigr)=\bigl(0,\,-e\mb{E}\bigr)
\tag{8}
\end{align*}
以上の結果から, 式 (3) の4元力 \(f_{\mu}\) と式 (8) の静止系の4元力 \(f’_{\mu}\) とは相等しいことが分かる.そして「2つの4元ベクトルがある座標系で成分が相等しいならば, それらは任意の座標系でも相等しい」はずである.よって, 任意の座標系で式 (2.8.30) が成り立つと言える:
\begin{align*}
f’_{\mu}=f_{\mu}\quad\rightarrow\quad \frac{e}{c}F_{\mu\nu}\,u^{\,\nu}=f_{\mu}
\tag{9}
\end{align*}
以上は C.メラー:「相対性理論」 § 58 を参照した.
前述の式 (2.8.17) から \(u^{2}=u^{\,\nu}u_{\nu}=c^2\) である.この両辺を固有時間 \(\tau\) で微分すると, 右辺 \(c^2\) は定数なのでゼロとなる.よって,
\begin{align*}
\frac{d}{d\tau}\,\bigl(u^{\,\nu}u_{\nu}\bigr)
= \frac{d u^{\,\nu}}{d\tau}u_{\nu}
+ u^{\,\nu}\frac{d u_{\nu}}{d\tau}
=u^{\,\nu}\frac{d u_{\nu}}{d\tau}+u^{\,\nu}\frac{d u_{\nu}}{d\tau}
=2u^{\,\nu}\frac{d u_{\nu}}{d\tau}=0
\tag{10}
\end{align*}
この式(10)から
\begin{align*}
u^{\,\nu}\frac{d u_{\nu}}{d\tau}=0
\tag{11}
\end{align*}
である.この式 (11) と式 (8-1) から, 4元速度 \(u^{\,\nu}\) と 4元力 \(f^{\,\nu}\) とは互いに直交することが分かる:
\begin{align*}
f\cdot u =f^{\,\nu}\,u_{\nu}=u_{\nu}\,\left(m \frac{d u^{\,\nu}}{d\tau}\right)
=m\left(u^{\,\nu}\frac{d u_{\nu}}{d\tau}\right)=0
\end{align*}


通常の3次元的な力の表式は,
\begin{align*}
\mb{F}=\frac{d}{dt}\mb{p}
=\frac{d}{dt}\left(\frac{m\mb{v}}{\sqrt{1-\beta^2}}\right)
\end{align*}
であった.これが式 (2.8.26) のローレンツ力である場合, 運動方程式は次となる:
\begin{align*}
\frac{d}{dt}\left(\frac{m\mb{v}}{\sqrt{1-\beta^2}}\right)
=e\mb{E}+\frac{e}{c}\mb{v}\times\mb{B}
\tag{8-4}
\end{align*}
ラグランジュ方程式
\begin{align*}
\frac{d}{dt}\left(\ppdiff{L}{u^{\,\mu}}\right)-\ppdiff{L}{x^{\,\mu}}=0
\tag{2.8.31}
\end{align*}
を直接適用して運動方程式 (8-4) を導くには, ラグランジアンを
\begin{align*}
L = -mc^2 \sqrt{1-\beta^2}-e\phi+\frac{e}{c}\mb{A}\cdot\mb{v},\qquad
\beta=\frac{v}{c}
\tag{8-5}
\end{align*}
とすればよいことを示すのは容易である.さらに, \(\mb{x}\) に「共役な運動量」\(\mb{P}\) は \(\partial L/\partial \mb{v}\) であり, それは次で与えられる:
\begin{align*}
\mb{P}=\ppdiff{L}{\mb{v}}=\frac{m\mb{v}}{\sqrt{1-\beta^2}} + \frac{e}{c}\mb{A}
\tag{2.8.32}
\end{align*}
これらに相当するハミルトニアンは
\begin{align*}
H=e\phi + \left[c^2\left(\mb{P}-\frac{e}{c}\mb{A}\right)^2+m^2c^4\right]^{1/2}
\tag{8-6}
\end{align*}
これは次式を満たす:
\begin{align*}
\left(\frac{\mathscr{H}-e\phi}{c}\right)^2 – \left(\mb{P}-\frac{e}{c}\mb{A}\right)^2 = m^2c^2
\tag{2.8.33}
\end{align*}
ハミルトニアンの考え方を「共変量」または「4次元の定式化」に変換するのは難しい.しかし,「作用 \(S\) は最小になるべきである」ことを述べた「最小作用の原理」:

\begin{equation*}
S=\int L\,dt \quad\rightarrow\quad \delta S = \delta\int L\,dt =0
\end{equation*}

は運動方程式の相対論的な形を導出する.それは, 作用 \(S\) を次とするときである: [6][訳註] 原書では式 (2.8.34)\(\sim\)(2.8.36) で \(\alpha\) がなぜか時間 \(t\) の代わりに用いられている.よく分からないが, ここでは時間 \(t\) … Continue reading
\begin{align*}
S=\int L\,dt&=-mc\int ds -\frac{e}{c}\int A_{\mu}dz^{\,\mu}\\
&= -\int \left\{mc\left(\frac{dz}{dt}\cdot\frac{dz}{dt}\right)^{1/2}
+\frac{e}{c}A_{\mu}\frac{dz^{\,\mu}}{dt}\right\}\,dt
\tag{2.8.34}
\end{align*}
このとき, 定義から次であることに注意する:
\begin{align*}
ds^2=dz^{\,\mu}dz_{\mu}\quad\rightarrow\quad
\left(\frac{ds}{dt}\right)^2
=\frac{dz^{\,\mu}}{dt}\cdot\frac{dz_{\mu}}{dt}
\tag{2.8.35}
\end{align*}
面白いことに, 次式で定義される別の「作用」 \(S’\) も前述の \(S\) と同じ結果を導出する:
\begin{align*}
S’=\frac{m}{2}\int \left(\frac{dz^{\,\mu}}{dt}\right)^2\,dt
+\frac{e}{c}\int A_{\,\mu}(z)\frac{dz^{\,\mu}}{dt}\,dt
\tag{2.8.36}
\end{align*}


【 問題 】 (1) 式 (8-5) のラグランジアンは式 (8-4) の運動方程式を導出し, 相当するハミルトニアンは式 (8-6) であることを示せ.また, \(\mb{P}\) に対する表式を見出せ.
(2) \(S\) がこの与えられた作用とするとき, \(S\) の変分 \(\delta S=0\) は同じ方程式を導くことを示せ.


〈解答例〉 (1) ランダウの §2 と §16 を要約することで答えに代える.
与えられた電磁場中で運動する粒子に対する作用は2つの部分から成る.すなわち, 自由粒子の作用 \(S_m\) と, 粒子の場との相互作用を記述する項 \(S’\) とから成っている.後者の \(S’\) は, 粒子を特徴付ける量と場を特徴付ける量との両方を含むはずである.電磁場との相互作用に関する粒子の性質は, 粒子の持つ電荷 \(e\) で規定される.場の性質は 4元ポテンシャル \(A^{\,\mu}\) によって特徴付けられる.これらの量は次の形の項として作用の中に現れる:
\begin{align*}
S_m &= -mc\int_a^{b}\,ds,\quad ds=cd\tau=cdt\sqrt{1-\beta^2},\\
S’&=-\frac{e}{c}\int_a^{b} A_{\mu}(x)\,dx^{\,\mu},\quad dx^{\,\mu}=(cdt,d\mb{x})
\ A^{\,\mu}=\bigl(\phi,\mb{A}\bigr)
\end{align*}
ここで \(A_{\mu}(x)\,dx^{\,\mu}=c\phi\,dt – \mb{A}\cdot d\mb{x}\) である.
従って, 電磁場中の荷電粒子に対する作用関数の形は次のように成る:
\begin{align*}
S= S_m + S’ &= -mc\int_a^{b} ds\,
-\frac{e}{c}\int_a^{b} A_{\mu}(x)\,dx^{\,\mu} \tag{16.1}\\
&= -mc^2dt\int_a^{b} \sqrt{1-\beta^2}\ -\frac{e}{c}\int_a^{b}
\bigl(\phi cdt-\mb{A}\cdot d\mb{x}\bigr)\\
&=\int_a^{b}\left\{-mc^2\sqrt{1-\beta^2}
-e\phi+\frac{e}{c}\mb{A}\cdot\frac{d\mb{x}}{dt}\right\}\,dt
=\int_a^{b} L\,dt \tag{16.3}
\end{align*}
粒子の速度は \(\mb{v}=d\mb{x}/dt\) である.従って,「電磁場中の電荷のラグランジアン」は次である:
\begin{align*}
L=-mc^2\sqrt{1-\beta^2} -e\phi + \frac{e}{c}\mb{A}\cdot\mb{v},
\qquad \beta=\frac{v}{c}
\tag{16.4}
\end{align*}
ラグランジュ方程式」は次である:
\begin{align*}
\frac{d}{dt}\left(\ppdiff{L}{\dot{\mb{x}}}\right)-\ppdiff{L}{\mb{x}}=0,\quad
\rightarrow\quad
\frac{d}{dt}\left(\ppdiff{L}{\mb{v}}\right)-\ppdiff{L}{\mb{x}}=0
\tag{17.1}
\end{align*}
このとき,「一般化運動量」は,
\begin{align*}
\mb{P}=\ppdiff{L}{\dot{\mb{x}}}=\ppdiff{L}{\mb{v}}
=\frac{m\mb{v}}{\sqrt{1-\beta^2}} +\frac{e}{c}\mb{A}=\mb{p}+\frac{e}{c}\mb{A}
\tag{16.5}
\end{align*}
ただし \(\mb{p}\) は「通常の(3次元の)運動量」である.また,

\begin{equation*}
\ppdiff{L}{\mb{x}}=\nabla L =\frac{e}{c}\nabla\,\mb{A}\cdot\mb{v}
-e\nabla \phi
\end{equation*}

ここで, よく知られた次のベクトル解析の公式を用いる:
\begin{equation*}
\nabla (\mb{a}\cdot\mb{b})=(\mb{a}\cdot\nabla)\mb{b}+(\mb{b}\cdot\nabla)\mb{a}
+\mb{a}\times\mathrm{rot}\,\mb{b} + \mb{b}\times\mathrm{rot}\,\mb{a}
\end{equation*}

ただし \(\mb{a},\mb{b}\) は任意の2つのベクトルである.この公式を \(\mb{A}\cdot\mb{v}\) に適用し,「\(\mb{x}\) についての微分は \(\mb{v}\) を一定にして行われる」ことに注意すれば,
\begin{equation*}
\ppdiff{L}{\mb{x}}=\frac{e}{c}(\mb{v}\cdot\nabla)\mb{A}
+\frac{e}{c}\mb{v}\times \mathrm{rot}\,\mb{A}-e\nabla\phi
\end{equation*}

となる.従ってラグランジュ方程式は次となる:
\begin{align*}
\frac{d}{dt}\left(\mb{p}+\frac{e}{c}\mb{A}\right)
=\frac{d\mb{p}}{dt}+\frac{e}{c}\frac{d\mb{A}}{dt}
=\frac{e}{c}(\mb{v}\cdot\nabla)\mb{A}+\frac{e}{c}\mb{v}\times\mathrm{rot}\,\mb{A}
-e\nabla\phi
\tag{17.1-a}
\end{align*}
ところで, 微分 \(d\mb{A}\) は時間 \(t\) と位置 \(\mb{x}\) の関数と考えられるから, 次のように書ける:
\begin{equation*}
d\mb{A}=\ppdiff{\mb{A}}{t}dt + \ppdiff{\mb{A}}{\mb{x}}\cdot d\mb{x}
=\ppdiff{\mb{A}}{t}dt + (\nabla\mb{A})\cdot d\mb{x}
\end{equation*}

従って,
\begin{equation*}
\frac{d\mb{A}}{dt}=\ppdiff{\mb{A}}{t}+(\nabla\mb{A})\cdot\frac{d\mb{x}}{dt}
=\ppdiff{\mb{A}}{t}+\mb{v}\cdot(\nabla\mb{A})
=\ppdiff{\mb{A}}{t}+(\mb{v}\cdot\nabla)\mb{A}
\end{equation*}

これを上式 (17.1-a) に用いると
\begin{align*}
\frac{d\mb{p}}{dt}+\frac{e}{c}\frac{d\mb{A}}{dt}
&=\frac{d\mb{p}}{dt}+\frac{e}{c}\left\{\ppdiff{\mb{A}}{t}+(\mb{v}\cdot\nabla)\mb{A}\right\}\\
&=\frac{d\mb{p}}{dt}+\frac{e}{c}\ppdiff{\mb{A}}{t}+\frac{e}{c}(\mb{v}\cdot\nabla)\mb{A}
=\frac{e}{c}(\mb{v}\cdot\nabla)\mb{A}+\frac{e}{c}\mb{v}\times\mathrm{rot}\,\mb{A}-e\nabla\phi
\end{align*}
従って, 次式が得られる:
\begin{align*}
\frac{d\mb{p}}{dt}=-\frac{e}{c}\ppdiff{\mb{A}}{t}-e\nabla\phi +\frac{e}{c}\mb{v}\times
\mathrm{rot}\,\mb{A}
\tag{17.2}
\end{align*}
さらに, 電場 \(\mb{E}\) と磁場 \(\mb{B}\) は次であった:
\begin{align*}
\mb{E}=-\frac{1}{c}\ppdiff{\mb{A}}{t} -\nabla\phi,\qquad \mb{B}=\mathrm{rot}\,\mb{A}
\tag{17.4}
\end{align*}
これを上式 (17.2) に代入すれば,「電磁場中の電荷の運動方程式」として次が得られる:
\begin{align*}
\frac{d\mb{p}}{dt}=e\mb{E} + \frac{e}{c}\mb{v}\times\mb{B}
\tag{17.5}
\end{align*}
このときの右辺の表現は「ローレンツ力」と呼ばれる.
一般公式
\begin{equation*}
\mathscr{H}=\mb{v}\cdot\ppdiff{L}{\mb{v}}-L
\end{equation*}

によって, ラグランジアンから「場中の粒子のハミルトニアン」を見出すことが出来る.式 (16.4) を代入すると, まず
\begin{equation*}
\ppdiff{L}{\mb{v}}=-mc^2\pdiff{\mb{v}}\left(1-\frac{\mb{v}^2}{c^2}\right)^{1/2}+\frac{e}{c}\mb{A}
=\frac{m\mb{v}}{\sqrt{1-\beta^2}}+\frac{e}{c}\mb{A}
\end{equation*}

従って,
\begin{align*}
\mathscr{H}&=\frac{m\mb{v}\cdot\mb{v}}{\sqrt{1-\beta^2}}+\frac{e}{c}\mb{A}\cdot\mb{v}+mc^{2}\sqrt{1-\beta^2}
-\frac{e}{c}\mb{A}\cdot\mb{v}+e\phi\\
&= \frac{mv^2}{\sqrt{1-\beta^2}}+mc^{2}\sqrt{1-\beta^2}+e\phi\\
&=\frac{mc^2}{\sqrt{1-\beta^2}}+e\phi
\tag{16.6}
\end{align*}
ところが, ハミルトニアンは「一般化運動量」に表さなければならない.式 (16.5) と (16.6) から, 場が無い (\(\phi=0,\mb{A}=0\)) とき
\begin{equation*}
\mb{P}=\mb{p},\qquad \mathscr{H}=\frac{mc^2}{\sqrt{1-\beta^2}}=E
\end{equation*}

となるから,
\(E\,(=\mathscr{H}-e\phi)\) と \(\displaystyle \mb{p}\,(=\mb{P}-\frac{e}{c}\mb{A})\) の間の関係すなわち \(\displaystyle \frac{E^{2}}{c^2}=\mb{p}+m^2c^2\) は, \(\mathscr{H}\) と \(\mb{p}\) との間の関係と同じであることが分かる.すなわち,
\begin{align*}
\frac{E^{2}}{c^2}=\mb{p}+m^2c^2\quad\rightarrow\quad
\left(\frac{\mathscr{H}-e\phi}{c}\right)^2 = \left(\mb{P}-\frac{e}{c}\mb{A}\right)^2 +m^2c^2
\tag{16.7}
\end{align*}
これを書き直すと次となる:
\begin{align*}
\mathscr{H}=e\phi + \sqrt{m^2c^4+c^2\left(\mb{P}-\frac{e}{c}\mb{A}\right)^2}
\tag{16.8}
\end{align*}

(2) (ランダウの §2 から) 「最小作用の原理」は次の形に書ける:
\begin{align*}
\delta S = \delta \int_{t_1}^{t_2} L(q,\dot{q},t)\,dt =0
\tag{2.4}
\end{align*}
あるいは「変分」を実行して,

\begin{equation*}
\int_{t_1}^{t_2} \left(\ppdiff{L}{q}\delta q + \ppdiff{L}{\dot{q}}\delta \dot{q}\right)\,dt=0
\end{equation*}

\(\displaystyle \delta \dot{q}=d\left(\frac{\delta q}{dt}\right)=\frac{d}{dt}\delta q\) を用い, 第2項を部分積分すると,
\begin{equation*}
\int_{t_1}^{t_2} \ppdiff{L}{\dot{q}}\delta \dot{q}\,dt=\int_{t_1}^{t_2} \ppdiff{L}{\dot{q}}\dot{(\delta q)}\,dt
=\left.\ppdiff{L}{\dot{q}}\right|_{t_1}^{t_2} -\int_{t_1}^{t_2} \frac{d}{dt}\left(\ppdiff{L}{\dot{q}}\right)\,\delta q\,dt
\end{equation*}

従って, 次式を得る:
\begin{align*}
\delta S=\left[\ppdiff{L}{\dot{q}}\delta q\right]_{t_1}^{t_2}
+\int_{t_1}^{t_2} \left(\ppdiff{L}{q}-\frac{d}{dt}\ppdiff{L}{\dot{q}}\right)\,\delta q
\,dt = 0
\tag{2.5}
\end{align*}
端点では変分はゼロであるという条件 \(\delta q(t_1)=\delta q(t_2)=0\) から第1項は消える.残った積分は, 任意に選ばれた \(\delta q\) に対してゼロでなけらばならない.このことは被積分関数が常にゼロである時にだけ可能である.このようにしてラグランジュ方程式が得られる:
\begin{equation*}
\frac{d}{dt}\ppdiff{L}{\dot{q}}-\ppdiff{L}{q}=0
\end{equation*}

このラグランジュ方程式を解くと 運動方程式 (8-4) を導くことが出来るのであるから, 結局式 (2.4) の \(\delta S=0\) も同じ運動方程式を導くのは明らかである.


References

References
1 [訳註] (ランダウ:「力学・場の力学」より) ポテンシャルの持つ任意性のため, 我々は常に付加条件をポテンシャルに課することが出来る.この理由のために, ここでは電磁波のポテンシャルをスカラー・ポテンシャルが方程式\(\phi=0\)を満たすような具合に選ぶことにする.
2 [訳註] 前講から4元座標は \(\mb{X}=(ct,\,\mb{X})\) であるから,
\begin{align*}
k\cdot X =k^{\,\mu}X_{\mu}=\left(\frac{\omega}{c},\,\mb{K}\right)\cdot
\left(ct,\,-\mb{X}\right)=\omega t -\mb{K}\cdot\mb{X}
\end{align*}
3 [訳註] (ランダウ:「力学・場の理論」§35より) ある慣性基準系 \(K\) から速さ \(v\) で運動している時計(その慣性系を \(K’\) とする) を観察する.微小時間 \(dt\) の間に運動している時計は距離 \(dr=\sqrt{dx^{2}+dy^{2}+dz^{2}}\) だけ進む.従って \(v=dr/dt\) である.このとき運動している時計が刻む時間間隔 \(d\tau\) を求めよう.座標系 \(K’\) では時計は静止して見える.従って \(dx’=dy’=dz’=0\) である.世界間隔 \(ds\) の不変性から次が言える:
\begin{align*}
&ds^{2}=c^{2}dt^{2}-dx^{2}-dy^{2}-dz^{2}=ds’=c^{2}d\tau^{2},\quad \rightarrow\quad
d\tau=dt\sqrt{1-\frac{dx^{2}+dy^{2}+dz^{2}}{c^{2}dt^{2}}},\\
&\mathrm{and}\quad v^{2}=\frac{dx^{2}+dy^{2}+dz^{2}}{dt^{2}}=v_x^2+v_y^2+v_z^2,\quad
v_x=\frac{dx}{dt},\ v_y=\frac{dy}{dt},\ v_z=\frac{dz}{dt}\\
&\therefore\quad
d\tau=\frac{ds}{c}=dt\sqrt{1-\frac{v^{2}}{c^{2}}}\quad \rightarrow\quad
dt=\frac{d\tau}{\sqrt{1-\beta^{2}}}=\gamma\,d\tau,\quad \beta\equiv \frac{v}{c}
\end{align*}
4 [訳註] (C.メラー §27 より) 粒子の運動量が変化していれば, 粒子は単位時間当たりの運動量の変化に等しい力 \(\mb{F}\) の作用を受けているとする.即ち,
\begin{align*}
\mb{F}=\frac{d\mb{p}}{dt}
\tag{24}
\end{align*}
方程式 (24) は, 粒子の速度が小さいときには Newton の運動の第2法則に等しくなるが, この「方程式 (24) を相対論的領域に於ける力の定義を与える式」と見做すことにする.なぜなら, 粒子の運動量の変化を生じる原因である系の物理的状態に対する力 \(\mb{F}\) の依存関係が与えられておれば, それを運動方程式と見做して差し支えないからである.
Newton 力学に於けると同様に, 単位時間当たりの力のする仕事 \(A\) は,
\begin{align*}
A=\mb{F}\cdot\mb{u}
\tag{25}
\end{align*}
によって定義する.ここで \(\mb{u}\) は粒子の速度である.更に, 粒子の運動エネルギー \(E\) は,
\begin{align*}
\frac{dE}{dt}=A=\mb{F}\cdot\mb{u}
\tag{26}
\end{align*}
によって定義される.この式は単位時間当たりの運動エネルギーの変化(Power)が, 仕事 \(A\) に等しいことを表している.4元力(Minkowskiの力) \(\mb{f}\) を導入すれば, これらの式は次となる:
\begin{align*}
\frac{d\mb{p}}{d\tau}=\mb{f},\qquad \frac{dE}{d\tau}=\mb{f}\cdot\mb{u}
\tag{42}
\end{align*}
従って, 4元運動量 \(p^{\,\mu}=(E/c,\mb{p})=mu^{\,\mu}\) を固有時間 \(\tau\) で微分したものは, 質量 \(m\) を定数とすると次となる:
\begin{align*}
&\frac{d p^{\,\mu}}{d\tau}=m\frac{d\,u^{\,\mu}}{d\tau}
=\left(\frac{1}{c}\frac{d E}{d\tau},\,\frac{d \mb{p}}{d\tau}\right)
=\left(\frac{1}{c}\mb{f}\cdot\mb{u},\,\mb{f}\right)\equiv f^{\,\mu}
=(f^{\,0},\,\mb{f}) \tag{55}\\
&\rightarrow\quad \frac{1}{c}\frac{d E}{d\tau}=f^{\,0}
=\frac{1}{c}\frac{\mb{F}\cdot\mb{u}}{\sqrt{1-\beta^2}},\quad \mathrm{and}\quad
\frac{d\mb{p}}{d\tau}=\mb{f}=\frac{\mb{F}}{\sqrt{1-\beta^2}}
\end{align*}