tecktom.com

当てる道具:Euler–Lagrange 方程式から Lorentz 力へ

2026年10月15日8分で読めますAI生成教科書『虎の巻』第2部 第2章 電磁場中の荷電粒子
目次
  1. 偏微分を2つ計算する
  2. Lorentz 力の式
  3. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第2章 §2.2.2 ★★★ LL を Euler–Lagrange 方程式に入れ、ベクトルの恒等式を1つ使うと mr¨=q(E+r˙×B)m\ddot{\bm{r}}=q(\bm{E}+\dot{\bm{r}}\times\bm{B}) が出る

当てる道具は Euler–Lagrange 方程式(第1部第4章)だけである。直交座標の成分ごとに ddt∂L∂x˙i−∂L∂xi=0\frac{\mathrm{d} }{\mathrm{d} t}\frac{\partial L}{\partial \dot x_i}-\frac{\partial L}{\partial x_i}=0 を書く(x1=xx_1=x, x2=yx_2=y, x3=zx_3=z)。

偏微分を2つ計算する

式 (2.2.2) を速度で偏微分すると

∂L∂r˙=mr˙+qA(r,t).(2.2.5)\frac{\partial L}{\partial \dot{\bm{r}}}=m\dot{\bm{r}}+q\bm{A}(\bm{r},t).\tag{2.2.5}

これを粒子の軌道に沿って時間で微分する。A\bm{A} は時刻と位置の両方を通して変わるので、

dAdt=∂A∂t+(r˙⋅∇)A,成分ではdAidt=∂Ai∂t+∑jx˙j∂Ai∂xj(2.2.6)\frac{\mathrm{d} \bm{A}}{\mathrm{d} t}=\frac{\partial \bm{A}}{\partial t}+(\dot{\bm{r}}\cdot\nabla)\bm{A},\qquad\text{成分では}\quad\frac{\mathrm{d} A_i}{\mathrm{d} t}=\frac{\partial A_i}{\partial t}+\sum_j\dot x_j\frac{\partial A_i}{\partial x_j}\tag{2.2.6}

となる。右辺第2項は、時刻を止めても粒子が場所を移ることで見る A\bm{A} が変わる分である。

位置で偏微分するときは速度を止めるので、

∂L∂r=−q∇ϕ+q∇(r˙⋅A),成分では∂L∂xi=−q∂ϕ∂xi+q∑jx˙j∂Aj∂xi(2.2.7)\frac{\partial L}{\partial \bm{r}}=-q\nabla\phi+q\nabla(\dot{\bm{r}}\cdot\bm{A}),\qquad\text{成分では}\quad\frac{\partial L}{\partial x_i}=-q\frac{\partial \phi}{\partial x_i}+q\sum_j\dot x_j\frac{\partial A_j}{\partial x_i}\tag{2.2.7}

である。

Lorentz 力の式

式 (2.2.5)〜(2.2.7) を Euler–Lagrange 方程式に入れ、mr¨m\ddot{\bm{r}} だけを左辺に残すと

mr¨=−q∇ϕ−q∂A∂t+q[∇(r˙⋅A)−(r˙⋅∇)A].m\ddot{\bm{r}}=-q\nabla\phi-q\frac{\partial \bm{A}}{\partial t}+q\bigl[\nabla(\dot{\bm{r}}\cdot\bm{A})-(\dot{\bm{r}}\cdot\nabla)\bm{A}\bigr].

角括弧は、r˙\dot{\bm{r}} を位置によらない定数とみなした恒等式

∇(v⋅A)−(v⋅∇)A=v×(∇×A)(2.2.8)\nabla(\bm{v}\cdot\bm{A})-(\bm{v}\cdot\nabla)\bm{A}=\bm{v}\times(\nabla\times\bm{A})\tag{2.2.8}

で書き直すことができる(成分では ∑jvj(∂iAj−∂jAi)\sum_jv_j(\partial_iA_j-\partial_jA_i) で、i=1i=1 なら v2(∂1A2−∂2A1)−v3(∂3A1−∂1A3)=v2B3−v3B2=(v×B)1v_2(\partial_1A_2-\partial_2A_1)-v_3(\partial_3A_1-\partial_1A_3)=v_2B_3-v_3B_2=(\bm{v}\times\bm{B})_1)。よって

mr¨=q(E+r˙×B),E=−∇ϕ−∂A∂t,B=∇×A(2.2.9)m\ddot{\bm{r}}=q\bigl(\bm{E}+\dot{\bm{r}}\times\bm{B}\bigr),\qquad\bm{E}=-\nabla\phi-\frac{\partial \bm{A}}{\partial t},\quad\bm{B}=\nabla\times\bm{A}\tag{2.2.9}

が得られる。これが Lorentz 力の式である。

導出の中で A\bm{A} の役は2つに分かれた。時間微分 −∂A/∂t-\partial\bm{A}/\partial t は電場に入り、回転 ∇×A\nabla\times\bm{A} は磁場になった。磁場の項は、∂L/∂r˙\partial L/\partial\dot{\bm{r}} の中の qAq\bm{A} を軌道に沿って微分したとき出る (r˙⋅∇)A(\dot{\bm{r}}\cdot\nabla)\bm{A} と、位置の偏微分で出る ∇(r˙⋅A)\nabla(\dot{\bm{r}}\cdot\bm{A}) の差で、例題 2.2.1 の反対称な行列 bijb_{ij} と同じものである。

運動方程式 (2.2.9) には ϕ\phi と A\bm{A} そのものは現れず、組み合わせ E\bm{E} と B\bm{B} だけが現れる。この事実から、E\bm{E} と B\bm{B} を変えない (ϕ,A)(\phi,\bm{A}) の取り替えが運動を変えないことがわかる。これを2.2.4節で LL の言葉で確かめる。

ポテンシャルから自動的に出る2本の Maxwell 方程式

E\bm{E} と B\bm{B} を式 (2.2.9) の形で書くと、恒等式 ∇⋅(∇×A)=0\nabla\cdot(\nabla\times\bm{A})=0 と ∇×∇ϕ=0\nabla\times\nabla\phi=0 から

∇⋅B=0,∇×E=−∂B∂t\nabla\cdot\bm{B}=0,\qquad\nabla\times\bm{E}=-\frac{\partial \bm{B}}{\partial t}

が場の運動とは無関係に成り立つ。磁荷が無いことと Faraday の法則は、ポテンシャルで書いた時点で決まっている。残りの2本(Gauss の法則と Ampère–Maxwell の法則)は、場の Lagrangian を選んで初めて出てくる(第2部第6章)。

例題 2.2.2 一様な磁場の中の円運動 ★★★

オリジナル一様な磁場 B=Bez\bm{B}=B\bm{e}_z を A=B2(−y,x,0)\bm{A}=\frac B2(-y,x,0), ϕ=0\phi=0 で表した Lagrangian L=12m(x˙2+y˙2+z˙2)+qB2(xy˙−yx˙)L=\frac12m(\dot x^2+\dot y^2+\dot z^2)+\frac{qB}{2}(x\dot y-y\dot x) を考える(q>0q>0, B>0B>0)。

  1. (1) xx, yy, zz の Euler–Lagrange 方程式を書け。
  2. (2) ωc≔qB/m\omega_{\mathrm c}\coloneqq qB/m とおき、zz 方向の速度が 00 の解が角速度 ωc\omega_{\mathrm c} の円運動であることを示せ。+z+z の側から見た回る向きを答えよ。
  3. (3) 正準運動量 pxp_x, pyp_y を求め、mx˙m\dot x, my˙m\dot y と違うことを確かめよ。
解答

(1) ∂L/∂x˙=mx˙−qB2y\partial L/\partial\dot x=m\dot x-\frac{qB}{2}y、∂L/∂x=qB2y˙\partial L/\partial x=\frac{qB}{2}\dot y なのでmx¨−qB2y˙−qB2y˙=0m\ddot x-\frac{qB}2\dot y-\frac{qB}2\dot y=0、すなわち mx¨=qBy˙m\ddot x=qB\dot y。同様に ∂L/∂y˙=my˙+qB2x\partial L/\partial\dot y=m\dot y+\frac{qB}2x、∂L/∂y=−qB2x˙\partial L/\partial y=-\frac{qB}2\dot x から my¨=−qBx˙m\ddot y=-qB\dot x。zz については mz¨=0m\ddot z=0。これは qr˙×B=qB(y˙,−x˙,0)q\dot{\bm{r}}\times\bm{B}=qB(\dot y,-\dot x,0) と一致する。

(2) u≔x˙+iy˙u\coloneqq\dot x+\mathrm{i}\dot y とおくと、(1) から u˙=ωc(y˙−ix˙)=−iωcu\dot u=\omega_{\mathrm c}(\dot y-\mathrm{i}\dot x)=-\mathrm{i}\omega_{\mathrm c}u。解は u=u0e−iωctu=u_0\mathrm{e}^{-\mathrm{i}\omega_{\mathrm c}t} で、速さ ∣u∣=∣u0∣\lvert u\rvert =\lvert u_0\rvert は一定、速度の向きは角速度 ωc\omega_{\mathrm c} で回る。積分すると x+iy=(xg+iyg)+iu0ωce−iωctx+\mathrm{i} y=(x_{\mathrm g}+\mathrm{i} y_{\mathrm g})+\frac{\mathrm{i} u_0}{\omega_{\mathrm c}}\mathrm{e}^{-\mathrm{i}\omega_{\mathrm c}t}(xgx_{\mathrm g}, ygy_{\mathrm g} は積分定数)で、中心 (xg,yg)(x_{\mathrm g},y_{\mathrm g})、半径 ∣u0∣/ωc\lvert u_0\rvert /\omega_{\mathrm c} の円の上を動く。e−iωct\mathrm{e}^{-\mathrm{i}\omega_{\mathrm c}t} は偏角が時間とともに減るので、+z+z の側から見て時計回りである。

(3) px=mx˙−qB2yp_x=m\dot x-\frac{qB}2y, py=my˙+qB2xp_y=m\dot y+\frac{qB}2x。どちらも位置を含む項 qAxqA_x, qAyqA_y の分だけ mx˙m\dot x, my˙m\dot y と違う。円運動の上で mx˙m\dot x は周期的に符号を変えるが、pxp_x の値は yy を含むので、同じ速度でも場所によって違う値になる。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ