tecktom.com

出てくるもの (3):一様な磁場と電場の中の運動

2026年10月15日9分で読めますAI生成教科書『虎の巻』第2部 第2章 電磁場中の荷電粒子
目次
  1. Hamiltonian を1次元の問題に帰着させる
  2. ドリフト
  3. 外から決めた場という近似
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第2章 §2.2.5 ★★ Landau ゲージでは保存する pxp_x を止めると1次元の調和振動子になる。円運動の角振動数 qB/mqB/m と、電場に直交する一定のドリフト速度 E/BE/B が出る

Hamiltonian を1次元の問題に帰着させる

一様な磁場 B=Bez\bm{B}=B\bm{e}_z と、それに直交する一様な電場 E=Eey\bm{E}=E\bm{e}_y があるとする(q>0q>0, B>0B>0)。ポテンシャルを ϕ=−Ey\phi=-Ey, A=(−By,0,0)\bm{A}=(-By,0,0) に選ぶ。∇×A=(0,0,−∂y(−By))=Bez\nabla\times\bm{A}=(0,0,-\partial_y(-By))=B\bm{e}_z で、−∇ϕ=Eey-\nabla\phi=E\bm{e}_y である。この A\bm{A} は前の節の AL\bm{A}_{\mathrm L} と xx, yy の役を入れ替えた Landau ゲージで、Hamiltonian (2.2.11) は

H=(px+qBy)22m+py22m+pz22m−qEy(2.2.19)H=\frac{(p_x+qBy)^2}{2m}+\frac{p_y^2}{2m}+\frac{p_z^2}{2m}-qEy\tag{2.2.19}

となる。HH は xx と zz を含まないので、pxp_x と pzp_z が保存する。pxp_x を定数として止めると、yy の運動だけが残る。ωc=qB/m\omega_{\mathrm c}=qB/m とおいて yy の部分を平方完成すると

Hy=py22m+12mωc2(y−yc)2+定数,yc=−pxqB+mEqB2(2.2.20)H_y=\frac{p_y^2}{2m}+\frac12m\omega_{\mathrm c}^2\bigl(y-y_{\mathrm c}\bigr)^2+\text{定数},\qquad y_{\mathrm c}=-\frac{p_x}{qB}+\frac{mE}{qB^2}\tag{2.2.20}

で、角振動数 ωc\omega_{\mathrm c}、中心 ycy_{\mathrm c} の調和振動子になる。磁場だけ(E=0E=0)なら、yy は角振動数 ωc\omega_{\mathrm c} で振動する。xx の速度 x˙=∂H/∂px=(px+qBy)/m=ωc(y−yc)\dot x=\partial H/\partial p_x=(p_x+qBy)/m=\omega_{\mathrm c}(y-y_{\mathrm c}) も同じ角振動数で振動し、2つを合わせると例題 2.2.2 の円運動になる。この ωc\omega_{\mathrm c} をサイクロトロン角振動数と呼ぶ。円の半径 ρ=v⊥/ωc=mv⊥/(qB)\rho=v_\perp/\omega_{\mathrm c}=mv_\perp/(qB)(v⊥v_\perp は磁場に垂直な速さ)はLarmor 半径と呼ばれる。ωc\omega_{\mathrm c} は速さによらないので、速い粒子は大きな円を、遅い粒子は小さな円を、同じ周期 2π/ωc2\pi/\omega_{\mathrm c} で回る。

E×B\bm{E}\times\bm{B} ドリフト

電場があると xx 方向の速度は

x˙=ωc(y−yc)+ωc(yc+pxqB)=ωc(y−yc)+EB(2.2.21)\dot x=\omega_{\mathrm c}(y-y_{\mathrm c})+\omega_{\mathrm c}\Bigl(y_{\mathrm c}+\frac{p_x}{qB}\Bigr)=\omega_{\mathrm c}(y-y_{\mathrm c})+\frac{E}{B}\tag{2.2.21}

となる。第1項は振動して時間平均が 00、第2項は一定である。円運動の中心は xx 方向に一定の速さ E/BE/B で動く。向きを込めると vd=E×B/B2=Eey×Bez/B2=(E/B)ex\bm{v}_{\mathrm d}=\bm{E}\times\bm{B}/B^2=E\bm{e}_y\times B\bm{e}_z/B^2=(E/B)\bm{e}_x で、電場にも磁場にも直交する。これを E×B\bm{E}\times\bm{B} ドリフトと呼ぶ。vd\bm{v}_{\mathrm d} は qq と mm を含まないので、電荷の符号や質量によらず、すべての粒子が同じ向きに同じ速さで流れる。

ωct\omega_{\mathrm c}t を時間の単位、E/(Bωc)E/(B\omega_{\mathrm c}) を長さの単位にとると、原点で静止していた粒子の軌道はx=ωct−sin⁡ωctx=\omega_{\mathrm c}t-\sin\omega_{\mathrm c}t, y=1−cos⁡ωcty=1-\cos\omega_{\mathrm c}t(サイクロイド)になる(図 2.2.2)。最初に電場で +y+y へ加速された粒子が磁場で曲げられ、電場に逆らって y=0y=0 まで戻り、そこで止まる。これを繰り返しながら xx 方向へ進む。

図

図 2.2.2 磁場 BezB\bm{e}_z(紙面の手前向き)と電場 EeyE\bm{e}_y の中の正電荷の軌道(0≤ωct≤4π0\le\omega_{\mathrm c}t\le4\pi、長さの単位 ℓ=E/(Bωc)\ell=E/(B\omega_{\mathrm c}))。点線は運動方程式の厳密な解、実線は円運動を平均して中心だけを追った運動(中心は xx 方向に速さ E/BE/B で進む)。どちらの軌道も中心は同じ速さで進み、初速度は円の半径と中心の高さだけを変える

外から決めた場という近似

この章の LL では、ϕ\phi と A\bm{A} は粒子の運動によらない与えられた関数であった。実際には加速度を持つ電荷は電磁波を放射してエネルギーを失う。非相対論的な電荷の放射の仕事率は Larmor の式 P=q2∣r¨∣2/(6πε0c3)P=q^2\lvert \ddot{\bm{r}}\rvert ^2/(6\pi\varepsilon_0c^3) で、円運動では ∣r¨∣=ωcv⊥\lvert \ddot{\bm{r}}\rvert =\omega_{\mathrm c}v_\perp なので、運動エネルギー 12mv⊥2\frac12mv_\perp^2 が減る割合は

P12mv⊥2=2τrad ωc2,τrad≔q26πε0mc3(2.2.22)\frac{P}{\frac12mv_\perp^2}=2\tau_{\mathrm{rad}}\,\omega_{\mathrm c}^2,\qquad\tau_{\mathrm{rad}}\coloneqq\frac{q^2}{6\pi\varepsilon_0mc^3}\tag{2.2.22}

になる。電子では τrad=6.27×10−24 s\tau_{\mathrm{rad}}=6.27\times10^{-24}\,\mathrm s で、B=1 TB=1\,\mathrm T では ωcτrad=1.10×10−12\omega_{\mathrm c}\tau_{\mathrm{rad}}=1.10\times10^{-12} である。エネルギーが 1/e1/\mathrm{e} に減る時間 1/(2τradωc2)1/(2\tau_{\mathrm{rad}}\omega_{\mathrm c}^2) は 2.58 s2.58\,\mathrm s で、1周の時間 3.57×10−11 s3.57\times10^{-11}\,\mathrm s の約 7×10107\times10^{10} 倍である。円運動を何周も計算する間は、放射の反作用を無視して外から決めた場の中の運動として扱ってよい。放射の反作用を入れた運動方程式は3階の微分方程式になり、小さいパラメータ ωcτrad\omega_{\mathrm c}\tau_{\mathrm{rad}} が最高階の微分に掛かる特異摂動の形をとる。

近似の裏付け

型 B(極限) | 小さいもの ωcτrad→0\omega_{\mathrm c}\tau_{\mathrm{rad}}\to0(放射の反作用の時間が運動の時間に比べて短い) | 証明 要確認(反作用を入れた方程式から外場中の運動への極限の前提を原論文で確かめていない) | 行き先 近似の数理 第7章

例題 2.2.5 1 T の磁場の中の電子と陽子 ★★

オリジナルB=1.00 TB=1.00\,\mathrm T の一様な磁場がある。電気素量 e=1.602×10−19 Ce=1.602\times10^{-19}\,\mathrm C、電子の質量 me=9.109×10−31 kgm_{\mathrm e}=9.109\times10^{-31}\,\mathrm{kg}、陽子の質量 mp=1.673×10−27 kgm_{\mathrm p}=1.673\times10^{-27}\,\mathrm{kg} とする。

  1. (1) 電子と陽子のサイクロトロン周波数 fc=ωc/(2π)f_{\mathrm c}=\omega_{\mathrm c}/(2\pi) を求めよ。
  2. (2) 磁場に垂直な運動エネルギーが 1.00 keV1.00\,\mathrm{keV} の電子の速さ v⊥v_\perp と Larmor 半径 ρ\rho を求め、v⊥/cv_\perp/c が十分小さいことを確かめよ。
  3. (3) さらに磁場に垂直な電場 E=1.00×104 V m−1E=1.00\times10^4\,\mathrm{V\,m^{-1}} を加えたときのドリフトの速さを求めよ。
解答

(1) 電子では ωc=eB/me=1.759×1011 rad s−1\omega_{\mathrm c}=eB/m_{\mathrm e}=1.759\times10^{11}\,\mathrm{rad\,s^{-1}}、fc=28.0 GHzf_{\mathrm c}=28.0\,\mathrm{GHz}。陽子では ωc=9.58×107 rad s−1\omega_{\mathrm c}=9.58\times10^7\,\mathrm{rad\,s^{-1}}、fc=15.2 MHzf_{\mathrm c}=15.2\,\mathrm{MHz}。質量の比 18361836 だけ陽子の方が遅い。

(2) 12mev⊥2=1.00×103×1.602×10−19 J\frac12m_{\mathrm e}v_\perp^2=1.00\times10^3\times1.602\times10^{-19}\,\mathrm J から v⊥=1.88×107 m s−1v_\perp=1.88\times10^7\,\mathrm{m\,s^{-1}}、v⊥/c=0.0626v_\perp/c=0.0626。相対論の補正は (v/c)2≈4×10−3(v/c)^2\approx4\times10^{-3} の大きさで、この章の非相対論の LL で扱ってよい。ρ=mev⊥/(eB)=v⊥/ωc=1.07×10−4 m=0.107 mm\rho=m_{\mathrm e}v_\perp/(eB)=v_\perp/\omega_{\mathrm c}=1.07\times10^{-4}\,\mathrm m=0.107\,\mathrm{mm}。

(3) vd=E/B=1.00×104 m s−1v_{\mathrm d}=E/B=1.00\times10^4\,\mathrm{m\,s^{-1}}。電子でも陽子でも同じで、v⊥v_\perp の約 1/19001/1900 である。ドリフトの速さが光速に近づく E∼cB=3×108 V m−1E\sim cB=3\times10^8\,\mathrm{V\,m^{-1}} では、相対論の LL(第2部第3章)が要る。

ほかの巻では

木下の勉強ステータス

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

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

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