背後的物理
帶電粒子在磁場中所受的磁力稱為勞侖茲力,方向由右手定則決定(負電荷相反):
\[ \vec F = q\,\vec v \times \vec B \qquad F = |q|\,vB\sin\theta \]
本模擬的磁場垂直紙面。速度在紙面內(\(\theta = 90^\circ\))時,磁力的大小固定為 \(|q|vB\),方向永遠垂直於速度。磁力只改變速度的方向、不改變速率,因此磁力不做功,動能不變。這個「大小固定、永遠垂直於速度」的力,正好提供等速圓周運動的向心力:
\[ |q|vB = \frac{mv^{2}}{r} \quad\Longrightarrow\quad r = \frac{mv}{|q|B} \]
\[ T = \frac{2\pi r}{v} = \frac{2\pi m}{|q|B} \qquad f = \frac{1}{T} = \frac{|q|B}{2\pi m} \]
週期 \(T\) 與速率 \(v\) 無關:跑得快的粒子圓也比較大,繞一圈所需時間恰好相同(迴旋加速器就是利用這個性質)。
螺旋運動:速度與磁場不垂直
如果速度與 \(\vec B\) 的夾角為 \(\theta\),把速度分解成垂直於磁場與平行於磁場的兩個分量:
\[ v_{\perp} = v\sin\theta \qquad v_{\parallel} = v\cos\theta \]
磁力 \(q\vec v\times\vec B\) 只與 \(v_{\perp}\) 有關,而且垂直於 \(\vec B\),所以不會改變 \(v_{\parallel}\)。沿著磁場方向看,粒子以 \(v_{\perp}\) 做等速圓周運動;沿著磁場方向則以 \(v_{\parallel}\) 等速前進,合起來就是螺旋線。每轉一圈沿 \(\vec B\) 前進的距離稱為螺距:
\[ r = \frac{m v_{\perp}}{|q|B} \qquad T = \frac{2\pi m}{|q|B} \qquad p = v_{\parallel}T = \frac{2\pi m\, v_{\parallel}}{|q|B} \]
週期與圓周運動相同。\(\theta = 90^\circ\) 時螺距為零,螺旋線變成圓;\(\theta = 0^\circ\) 或 \(180^\circ\) 時速度與磁場平行,完全不受磁力,粒子沿 \(\vec B\) 直線前進。2D 畫面看到的是投影到紙面(垂直於 \(\vec B\))上的運動,所以螺旋線看起來是半徑 \(r\) 的圓;3D 檢視可以看到完整的螺旋。
- 極光:太陽風帶來的帶電粒子沿地球磁力線螺旋前進。磁力線在磁極附近伸入大氣,粒子在高緯度地區撞擊空氣分子使其發光。
- 磁瓶與范艾倫輻射帶:磁力線變密的地方磁場較強,螺旋前進的粒子會被推回磁場較弱的一側,像被「反射」一樣(磁鏡)。地球磁場因此困住在南北半球之間來回彈跳的帶電粒子,形成范艾倫輻射帶;核融合實驗也用類似的磁瓶來約束高溫電漿。
速度選擇器
再加上一個與 \(\vec v\)、\(\vec B\) 都垂直的電場,讓電力 \(q\vec E\) 與磁力 \(q\vec v\times\vec B\) 方向相反。兩力大小相等時合力為零,粒子直線前進:
\[ |q|E = |q|vB \quad\Longrightarrow\quad v = \frac{E}{B} \]
條件中沒有 \(q\) 和 \(m\):只要速率剛好等於 \(E/B\),不論正負電荷、質量多少都能直線通過;比較快的粒子磁力較大、比較慢的粒子電力較大,都會偏離。(若 \(\theta \ne 90^\circ\),只有 \(v_{\perp}\) 受磁力,條件變成 \(v_{\perp} = E/B\)。)
延伸:E×B 漂移
速率不等於 \(E/B\) 時,粒子一邊繞圈、一邊以固定的速度沿 \(\vec E\times\vec B\) 方向前進(由靜止釋放時軌跡是擺線):
\[ \vec v_{d} = \frac{\vec E \times \vec B}{B^{2}} \qquad |\vec v_{d}| = \frac{E}{B} \]
在以 \(\vec v_{d}\) 移動的座標系中看,電場的效果剛好被抵消,粒子只做半徑 \(m|\vec v - \vec v_{d}|/(|q|B)\) 的圓周運動。畫面上的虛線圓與圓心標記就是這個「導引中心」。
觀察重點
- 旋轉方向:磁場射出紙面時,正電荷順時針繞圈、負電荷逆時針。先用右手定則預測,再按「發射」驗證;把 B 改成射入紙面又如何?
- 半徑:把 \(v\) 或 \(m\) 加倍、或把 \(|q|\) 或 \(B\) 加倍,半徑各變成幾倍?畫面上的 1 cm 方格可以直接量。
- 週期與速率無關:改變發射速率,比較「經過時間」與「速度方向轉過的圈數」,算出週期是否改變。
- 動能不變:純磁場中「目前速率」和「動能」始終不變;加上電場後動能就會改變,因為電力會做功。
- 螺旋線:選「螺旋運動」並切換到 3D 檢視。改變 \(\theta\):半徑是否隨 \(\sin\theta\)、螺距是否隨 \(\cos\theta\) 變化?週期有沒有改變?\(\theta = 0^\circ\)、\(90^\circ\)、\(180^\circ\) 時又如何?\(\theta > 90^\circ\) 時粒子沿 \(-z\) 前進,繞圈方向有沒有改變?
- 速率與 \(v_{\parallel}\) 守恆:螺旋運動時「目前速率」與「動能」不變,而且 \(z\) 每經過一個週期增加相同的量,正好是一個螺距 \(p\)。
- 射入磁場區:垂直邊界射入的粒子走半圓,在距離入射點 \(2r\) 處反向射出,在場區內的時間是 \(T/2\)。質譜儀就是用 \(2r = 2mv/(|q|B)\) 來分辨不同質量的離子。
- 速度選擇器:在「速度選擇器」情境中,只改 \(q\) 的正負或 \(m\),粒子仍直線前進;只改 \(v\),粒子就會偏轉。往哪一邊偏?
模擬方法與單位
- 使用真實的 SI 單位:電荷以基本電荷 \(e = 1.602\times 10^{-19}\ \mathrm{C}\) 為單位,質量以原子質量單位 \(1\ \mathrm{u} = 1.661\times 10^{-27}\ \mathrm{kg}\) 為單位,速率以 \(10^{5}\ \mathrm{m/s}\) 為單位,磁場以 \(\mathrm{mT}\)(\(10^{-3}\ \mathrm{T}\))為單位。例如 α 粒子以 \(3\times 10^{5}\ \mathrm{m/s}\) 進入 \(0.1\ \mathrm{T}\) 的磁場,半徑約 \(6.2\ \mathrm{cm}\)、週期約 \(1.3\ \mu\mathrm{s}\)。
- 真實的運動只需百萬分之一秒,所以用慢動作播放:「慢動作」設定動畫中 1 秒對應的真實時間。
- 運動方程式 \(m\,\dfrac{d\vec v}{dt} = q(\vec E + \vec v\times\vec B)\) 以 Boris 法在三維空間中數值積分:每一小步先半個電場加速、再讓速度繞 \(\vec B\) 旋轉、最後再半個電場加速。純磁場時這一步是純旋轉,所以速率精確守恆、\(v_{\parallel}\) 不變,軌跡點投影到紙面也正好落在半徑 \(mv_{\perp}/(|q|B)\) 的圓上。
- 自適應步長:每一步速度方向最多轉 0.01 rad、位移不超過 0.01 cm;粒子跨越場區邊界時再把步長細分,讓進出場區的位置準確。粒子離開 24 cm × 16 cm 的畫面,或沿 \(\vec B\) 移動超過 48 cm 時停止。
- 忽略重力(比磁力小好幾百萬倍以上)與粒子加速時的電磁輻射,也不考慮相對論效應(速率遠小於光速)。
- 3D 檢視使用 three.js 函式庫(MIT 授權),只有切換到 3D 時才會下載。