tecktom.com

この章のまとめ

2026年10月15日12分で読めますAI生成教科書『虎の巻』第2部 第2章 電磁場中の荷電粒子

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第2章 §2.2.7 ★★★ 速度の1次の項を持つ LL から、Lorentz 力・p=mr˙+qA\bm{p}=m\dot{\bm{r}}+q\bm{A}・最小結合の HH・ゲージの自由度・サイクロトロン運動とドリフトが出た

この章では、荷電粒子の Lagrangian として L=12mr˙2−qϕ+qr˙⋅AL=\frac12m\dot{\bm{r}}^2-q\phi+q\dot{\bm{r}}\cdot\bm{A} を選んだ。速度に比例する力は位置だけのポテンシャルからは出ないので、速度を含む一般化ポテンシャル U=qϕ−qr˙⋅AU=q\phi-q\dot{\bm{r}}\cdot\bm{A} を使った。速度の1次の一般化ポテンシャルから出る速度に比例する力は反対称な行列の形をとり、仕事をしない。Euler–Lagrange 方程式は Lorentz 力の式 mr¨=q(E+r˙×B)m\ddot{\bm{r}}=q(\bm{E}+\dot{\bm{r}}\times\bm{B}) になり、電場と磁場は E=−∇ϕ−∂tA\bm{E}=-\nabla\phi-\partial_t\bm{A}, B=∇×A\bm{B}=\nabla\times\bm{A} の組み合わせで現れた。

正準運動量は p=mr˙+qA\bm{p}=m\dot{\bm{r}}+q\bm{A} で、Hamiltonian は H=(p−qA)2/(2m)+qϕH=(\bm{p}-q\bm{A})^2/(2m)+q\phi になる。HH の値は 12mr˙2+qϕ\frac12m\dot{\bm{r}}^2+q\phi で、磁場は値に入らない。速度の成分どうしの Poisson 括弧は {mx˙,my˙}=qBz\{m\dot x,m\dot y\}=qB_z で、円運動の中心の2つの座標は互いに共役な組になる。A→A+∇Λ\bm{A}\to\bm{A}+\nabla\Lambda, ϕ→ϕ−∂tΛ\phi\to\phi-\partial_t\Lambda のゲージ変換で LL は全微分 q dΛ/dtq\,\mathrm{d}\Lambda/\mathrm{d} t だけ変わり、運動は同じで、p\bm{p} と HH の値が変わる。この取り替えは母関数 r⋅P−qΛ\bm{r}\cdot\bm{P}-q\Lambda の正準変換で、量子力学では波動関数の位相 eiqΛ/ℏ\mathrm{e}^{\mathrm{i} q\Lambda/\hbar} になる。一様な磁場では、ゲージによって保存する正準運動量の成分が違って見えるが、保存する物理量(円の中心の座標)は同じである。Landau ゲージで Hamiltonian を1次元の調和振動子に帰着させると、サイクロトロン角振動数 qB/mqB/m と E×B\bm{E}\times\bm{B} ドリフトの速度 E×B/B2\bm{E}\times\bm{B}/B^2 が出た。

次の章(第2部第3章)では、自由粒子の部分を相対論の −mc21−v2/c2-mc^2\sqrt{1-v^2/c^2} に替え、相互作用の項 −q∫Aμ dxμ-q\int A_\mu\,\mathrm{d} x^\mu と一緒に Lorentz 変換で不変な作用として書く。

この章の道具

  • 一般化ポテンシャル U(x,x˙,t)U(x,\dot x,t)、一般化力 Qi=ddt∂U∂x˙i−∂U∂xiQ_i=\frac{\mathrm{d} }{\mathrm{d} t}\frac{\partial U}{\partial \dot x_i}-\frac{\partial U}{\partial x_i}。速度の1次の部分から出る力は仕事をしない。
  • L=12mr˙2−qϕ+qr˙⋅AL=\frac12m\dot{\bm{r}}^2-q\phi+q\dot{\bm{r}}\cdot\bm{A} ⇒\Rightarrow mr¨=q(E+r˙×B)m\ddot{\bm{r}}=q(\bm{E}+\dot{\bm{r}}\times\bm{B})、E=−∇ϕ−∂tA\bm{E}=-\nabla\phi-\partial_t\bm{A}, B=∇×A\bm{B}=\nabla\times\bm{A}。
  • p=mr˙+qA\bm{p}=m\dot{\bm{r}}+q\bm{A}、H=(p−qA)2/(2m)+qϕH=(\bm{p}-q\bm{A})^2/(2m)+q\phi、{pi−qAi,pj−qAj}=q∑kϵijkBk\{p_i-qA_i,p_j-qA_j\}=q\sum_k\epsilon_{ijk}B_k。
  • ゲージ変換 A→A+∇Λ\bm{A}\to\bm{A}+\nabla\Lambda, ϕ→ϕ−∂tΛ\phi\to\phi-\partial_t\Lambda で L→L+q dΛ/dtL\to L+q\,\mathrm{d}\Lambda/\mathrm{d} t、p→p+q∇Λ\bm{p}\to\bm{p}+q\nabla\Lambda、H→H−q∂tΛH\to H-q\partial_t\Lambda、ψ→eiqΛ/ℏψ\psi\to\mathrm{e}^{\mathrm{i} q\Lambda/\hbar}\psi。
  • 一様な磁場:ωc=qB/m\omega_{\mathrm c}=qB/m、Larmor 半径 mv⊥/(qB)mv_\perp/(qB)、ドリフト E×B/B2\bm{E}\times\bm{B}/B^2。

総合問題 磁場の中の2次元調和振動子 ★★★

オリジナル質量 mm、電荷 q>0q>0 の粒子が xyxy 平面内で、原点に向かう調和ポテンシャル V=12mω02(x2+y2)V=\frac12m\omega_0^2(x^2+y^2) と、一様な磁場 B=Bez\bm{B}=B\bm{e}_z(B>0B>0)の両方の中を動く。ωc=qB/m\omega_{\mathrm c}=qB/m、ωL≔ωc/2\omega_{\mathrm L}\coloneqq\omega_{\mathrm c}/2、ω1≔ω02+ωL2\omega_1\coloneqq\sqrt{\omega_0^2+\omega_{\mathrm L}^2} とおく。

  1. (1) (使う道具:2.2.1節・第1部第2章の極座標)対称ゲージ A=B2(−y,x,0)\bm{A}=\frac B2(-y,x,0) で LL を極座標 (r,θ)(r,\theta) で書き、prp_r, pθp_\theta を求めよ。
  2. (2) (使う道具:2.2.3節・第1部第8章)Hamiltonian が H=pr22m+pθ22mr2+12mω12r2−ωLpθH=\dfrac{p_r^2}{2m}+\dfrac{p_\theta^2}{2mr^2}+\dfrac12m\omega_1^2r^2-\omega_{\mathrm L}p_\theta になることを示し、保存量を2つ挙げよ。
  3. (3) (使う道具:第1部第8章の Hamilton 方程式・例題 2.2.2)pθp_\theta を与えたときの円軌道の半径 r0r_0 と角速度 θ˙\dot\theta を求めよ。pθ>0p_\theta>0 と pθ<0p_\theta<0 の円軌道の角速度の大きさ ω±=ω1±ωL\omega_\pm=\omega_1\pm\omega_{\mathrm L} が、複素数 x+iyx+\mathrm{i} y の運動方程式の2つの固有振動数と一致することを確かめよ。
  4. (4) (使う道具:2.2.4節・第1部第10章の準対称性)Landau ゲージ A=B(0,x,0)\bm{A}=B(0,x,0) に取り替えると、LL は回転 δx=−ϵy\delta x=-\epsilon y, δy=ϵx\delta y=\epsilon x で不変でなくなる。δL\delta L が全微分 ϵ dΛN/dt\epsilon\,\mathrm{d}\Lambda_{\mathrm N}/\mathrm{d} t になることを示し、Noether の保存量 pxKx+pyKy−ΛNp_xK_x+p_yK_y-\Lambda_{\mathrm N} が (1) の pθp_\theta と同じ量であることを示せ。
  5. (5) (使う道具:(3) の結果)ωc/(2π)=1.00 MHz\omega_{\mathrm c}/(2\pi)=1.00\,\mathrm{MHz}, ω0/(2π)=0.100 MHz\omega_0/(2\pi)=0.100\,\mathrm{MHz} のとき ω±/(2π)\omega_\pm/(2\pi) を求めよ。ω0≪ωc\omega_0\ll\omega_{\mathrm c} のとき ω+≈ωc+ω02/ωc\omega_+\approx\omega_{\mathrm c}+\omega_0^2/\omega_{\mathrm c}, ω−≈ω02/ωc\omega_-\approx\omega_0^2/\omega_{\mathrm c} となることを示し、ω0→0\omega_0\to0 で遅い方の運動が何になるかを述べよ。
解答

(1) xy˙−yx˙=r2θ˙x\dot y-y\dot x=r^2\dot\theta, x˙2+y˙2=r˙2+r2θ˙2\dot x^2+\dot y^2=\dot r^2+r^2\dot\theta^2 なので

L=12m(r˙2+r2θ˙2)+qB2r2θ˙−12mω02r2.L=\frac12m(\dot r^2+r^2\dot\theta^2)+\frac{qB}{2}r^2\dot\theta-\frac12m\omega_0^2r^2 .

pr=mr˙p_r=m\dot r、pθ=mr2θ˙+qB2r2=mr2(θ˙+ωL)p_\theta=mr^2\dot\theta+\frac{qB}2r^2=mr^2(\dot\theta+\omega_{\mathrm L})。

(2) θ˙=pθ/(mr2)−ωL\dot\theta=p_\theta/(mr^2)-\omega_{\mathrm L} を使う。pθθ˙−qB2r2θ˙−12mr2θ˙2=mr2(θ˙+ωL)θ˙−mωLr2θ˙−12mr2θ˙2=12mr2θ˙2p_\theta\dot\theta-\frac{qB}2r^2\dot\theta-\frac12mr^2\dot\theta^2=mr^2(\dot\theta+\omega_{\mathrm L})\dot\theta-m\omega_{\mathrm L}r^2\dot\theta-\frac12mr^2\dot\theta^2=\frac12mr^2\dot\theta^2 なので

H=pr22m+12mr2(pθmr2−ωL)2+12mω02r2=pr22m+pθ22mr2−ωLpθ+12m(ω02+ωL2)r2.H=\frac{p_r^2}{2m}+\frac12mr^2\Bigl(\frac{p_\theta}{mr^2}-\omega_{\mathrm L}\Bigr)^2+\frac12m\omega_0^2r^2 =\frac{p_r^2}{2m}+\frac{p_\theta^2}{2mr^2}-\omega_{\mathrm L}p_\theta+\frac12m(\omega_0^2+\omega_{\mathrm L}^2)r^2 .

HH は時刻も θ\theta も含まないので、HH と pθp_\theta が保存する。−ωLpθ-\omega_{\mathrm L}p_\theta の項は、時計回りに角速度 ωL\omega_{\mathrm L} で回る座標系から見ると磁場が消え、閉じ込めが ω1\omega_1 に強まることを表している。

(3) 動径の有効ポテンシャル pθ22mr2+12mω12r2\frac{p_\theta^2}{2mr^2}+\frac12m\omega_1^2r^2 の最小から、∂H/∂r=−pθ2mr3+mω12r=0\partial H/\partial r=-\frac{p_\theta^2}{mr^3}+m\omega_1^2r=0、

r02=∣pθ∣mω1.r_0^2=\frac{\lvert p_\theta\rvert }{m\omega_1}.

角速度は θ˙=∂H/∂pθ=pθmr02−ωL=±ω1−ωL\dot\theta=\partial H/\partial p_\theta=\frac{p_\theta}{mr_0^2}-\omega_{\mathrm L}=\pm\omega_1-\omega_{\mathrm L}(pθ≷0p_\theta\gtrless0)。pθ>0p_\theta>0 では反時計回りに ω1−ωL=ω−\omega_1-\omega_{\mathrm L}=\omega_-、pθ<0p_\theta<0 では時計回りに ω1+ωL=ω+\omega_1+\omega_{\mathrm L}=\omega_+。一方、運動方程式は mx¨=−mω02x+qBy˙m\ddot x=-m\omega_0^2x+qB\dot y, my¨=−mω02y−qBx˙m\ddot y=-m\omega_0^2y-qB\dot x で、w=x+iyw=x+\mathrm{i} y について w¨=−ω02w−iωcw˙\ddot w=-\omega_0^2w-\mathrm{i}\omega_{\mathrm c}\dot w。w=e−iωtw=\mathrm{e}^{-\mathrm{i}\omega t} を入れると ω2−ωcω−ω02=0\omega^2-\omega_{\mathrm c}\omega-\omega_0^2=0、ω=ωL±ω1\omega=\omega_{\mathrm L}\pm\omega_1。ω=ωL+ω1=ω+\omega=\omega_{\mathrm L}+\omega_1=\omega_+ は時計回り、ω=ωL−ω1=−ω−\omega=\omega_{\mathrm L}-\omega_1=-\omega_- は反時計回りの回転で、円軌道の結果と一致する。

(4) Landau ゲージでは L=12m(x˙2+y˙2)+qBxy˙−12mω02(x2+y2)L=\frac12m(\dot x^2+\dot y^2)+qBx\dot y-\frac12m\omega_0^2(x^2+y^2)。回転で運動エネルギーと VV は変わらず、

δL=qB(δx y˙+x δy˙)=qBϵ(−yy˙+xx˙)=ϵddt[qB2(x2−y2)]\delta L=qB(\delta x\,\dot y+x\,\delta\dot y)=qB\epsilon(-y\dot y+x\dot x)=\epsilon\frac{\mathrm{d} }{\mathrm{d} t}\Bigl[\frac{qB}2(x^2-y^2)\Bigr]

なので ΛN=qB2(x2−y2)\Lambda_{\mathrm N}=\frac{qB}2(x^2-y^2)。Kx=−yK_x=-y, Ky=xK_y=x, px=mx˙p_x=m\dot x, py=my˙+qBxp_y=m\dot y+qBx から

pxKx+pyKy−ΛN=m(xy˙−yx˙)+qBx2−qB2(x2−y2)=m(xy˙−yx˙)+qB2(x2+y2)=mr2θ˙+qB2r2\begin{aligned}p_xK_x+p_yK_y-\Lambda_{\mathrm N}&=m(x\dot y-y\dot x)+qBx^2-\frac{qB}2(x^2-y^2)\\ &=m(x\dot y-y\dot x)+\frac{qB}2(x^2+y^2)=mr^2\dot\theta+\frac{qB}2r^2\end{aligned}

で、(1) の pθp_\theta に等しい。ゲージを替えると回転は対称性から準対称性に変わるが、保存量は同じ量になる。

(5) 単位を 2π×1 MHz2\pi\times1\,\mathrm{MHz} とすると ωL=0.500\omega_{\mathrm L}=0.500, ω1=0.0100+0.250=0.50990\omega_1=\sqrt{0.0100+0.250}=0.50990 で、ω+/(2π)=1.0099 MHz\omega_+/(2\pi)=1.0099\,\mathrm{MHz}, ω−/(2π)=0.0099 MHz\omega_-/(2\pi)=0.0099\,\mathrm{MHz}。ω0≪ωc\omega_0\ll\omega_{\mathrm c} では ω1=ωL1+ω02/ωL2≈ωL+ω02/(2ωL)=ωL+ω02/ωc\omega_1=\omega_{\mathrm L}\sqrt{1+\omega_0^2/\omega_{\mathrm L}^2}\approx\omega_{\mathrm L}+\omega_0^2/(2\omega_{\mathrm L})=\omega_{\mathrm L}+\omega_0^2/\omega_{\mathrm c} なのでω+≈ωc+ω02/ωc\omega_+\approx\omega_{\mathrm c}+\omega_0^2/\omega_{\mathrm c}(数値では 1.0100 MHz1.0100\,\mathrm{MHz})、ω−≈ω02/ωc\omega_-\approx\omega_0^2/\omega_{\mathrm c}(0.0100 MHz0.0100\,\mathrm{MHz})。速い方は磁場によるサイクロトロン運動で、遅い方は円運動の中心が原点のまわりをゆっくり回る運動である。中心の位置 rg\bm{r}_{\mathrm g} では閉じ込めの力 −mω02rg-m\omega_0^2\bm{r}_{\mathrm g} が電場 −mω02rg/q-m\omega_0^2\bm{r}_{\mathrm g}/q の役をし、2.2.5節のドリフトの速さ E/B=mω02rg/(qB)=ω02rg/ωcE/B=m\omega_0^2r_{\mathrm g}/(qB)=\omega_0^2r_{\mathrm g}/\omega_{\mathrm c} で接線方向に動くので、角速度は ω02/ωc\omega_0^2/\omega_{\mathrm c} になる。ω0→0\omega_0\to0 で ω−→0\omega_-\to0 となり、中心は止まって例題 2.2.3 の保存する (xg,yg)(x_{\mathrm g},y_{\mathrm g}) に戻る。

木下の勉強ステータス

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

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

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