108 年 國立中央大學太空科學研究所碩士班一般生《太空物理學》

本題主要考察電離層中，不同無線電波（O-wave, X-wave, Z-wave）在垂直入射時的反射頻率，以及電子密度與反射頻率的關係。

核心觀念：

1. 電離層對無線電波的反射與電子密度、波的頻率、磁場強度以及波的偏振態有關。
2. 反射頻率（或稱臨界頻率）是指當無線電波的頻率達到某個值時，電離層會開始反射該無線電波。
3. 電子密度越高，反射頻率越高。
4. 在考慮地磁場影響時，O-wave（ordinary wave）和X-wave（extraordinary wave）是兩種主要的波。Z-wave（Whistler wave）在特定條件下也會出現。

解題步驟：

(a) 計算 O-wave 的反射頻率

在垂直入射且忽略碰撞（collisionless）的情況下，無線電波的反射發生在離共振頻率（plasma frequency, $\omega_p$）或磁共振頻率（cyclotron frequency, $\omega_c$）處。對於 O-wave，其反射發生在 Plasma frequency $\omega_p$。

Plasma frequency 的定義為：
$\omega_p = \sqrt{\frac{n_e e^2}{\epsilon_0 m_e}}$
其中：
$n_e$ 是電子密度
$e$ 是基本電荷
$\epsilon_0$ 是真空介電常數
$m_e$ 是電子質量

題目給定的電子密度 $n_e = 10^6 \text{ el/cm}^3$。我們需要將其轉換為 SI 單位（m$^{-3}$）：
$n_e = 10^6 \text{ el/cm}^3 = 10^6 \times (100 \text{ cm/m})^3 \text{ el/m}^3 = 10^6 \times 10^6 \text{ el/m}^3 = 10^{12} \text{ el/m}^3$

已知常數：
$e \approx 1.602 \times 10^{-19} \text{ C}$
$m_e \approx 9.109 \times 10^{-31} \text{ kg}$
$\epsilon_0 \approx 8.854 \times 10^{-12} \text{ F/m}$

計算 $\omega_p$：
$\omega_p = \sqrt{\frac{(10^{12} \text{ m}^{-3})(1.602 \times 10^{-19} \text{ C})^2}{(8.854 \times 10^{-12} \text{ F/m})(9.109 \times 10^{-31} \text{ kg})}}$
$\omega_p = \sqrt{\frac{10^{12} \times (1.602)^2 \times 10^{-38}}{8.854 \times 9.109 \times 10^{-43}}} \text{ rad/s}$
$\omega_p = \sqrt{\frac{2.566 \times 10^{-26}}{8.064 \times 10^{-43}}} \text{ rad/s}$
$\omega_p = \sqrt{0.318 \times 10^{17}} \text{ rad/s} \approx \sqrt{3.18 \times 10^{16}} \text{ rad/s}$
$\omega_p \approx 1.78 \times 10^8 \text{ rad/s}$

頻率 $f_p$ 的單位是 Hz：
$f_p = \frac{\omega_p}{2\pi} = \frac{1.78 \times 10^8 \text{ rad/s}}{2\pi} \approx \frac{1.78 \times 10^8}{6.283} \text{ Hz} \approx 2.83 \times 10^7 \text{ Hz}$

本題考察電離層 F2 層的特性，特別是臨界頻率（critical frequency）與實際傳播頻率的關係。

核心觀念：

1. 臨界頻率 ($f_c$) 是指在垂直入射時，電離層某一層能夠反射的最高頻率。
2. 當無線電波的頻率 ($f$) 大於臨界頻率 ($f_c$) 時，電離層將無法反射該無線電波，波會穿透過去。
3. 當無線電波的頻率 ($f$) 小於或等於臨界頻率 ($f_c$) 時，電離層可以反射該無線電波。
4. 對於非垂直入射，電離層可以反射的最高頻率會更高，這個頻率稱為最大可用頻率 (Maximum Usable Frequency, MUF)。MUF 與 $f_c$ 和入射角有關，通常表示為 $MUF = f_c \sec(\theta)$，其中 $\theta$ 是波線與垂直線的夾角。

解題步驟：

題目給定 F2 層的臨界頻率 $f_{c,F2} = 10 \text{ MHz}$。
發射的 ionosonde 頻率為 $f_{transmit} = 20 \text{ MHz}$。

我們要判斷在什麼條件下，可以接收到從 F2 層反射回來的 20 MHz 的回波。

首先，我們注意到發射頻率 $f_{transmit} = 20 \text{ MHz}$ 遠大於 F2 層的臨界頻率 $f_{c,F2} = 10 \text{ MHz}$。
如果是在垂直入射的情況下，20 MHz 的無線電波將會穿透 F2 層，而無法被反射回來。

因此，為了使 20 MHz 的無線電波能夠被 F2 層反射並被接收到，必須考慮非垂直入射的情況。
在非垂直入射時，電離層能夠反射的最高頻率由 MUF 決定。

本題要求定義並描述電離層中 D, E, F1, F2 四個主要層的特徵。

核心觀念：
電離層是地球上層大氣中，高度約在 60 公里到 1000 公里之間，因受太陽紫外線和 X 射線的游離作用而含有大量自由電子和離子的區域。其電子密度隨高度變化，形成不同的層。這些層的特性（高度、電子密度、形成原因、對無線電波的影響）隨時間（日夜、季節、太陽活動週期）而變化。

解題步驟：

依據高度和電子密度特性，電離層通常被劃分為 D 層、E 層、F1 層和 F2 層。

1. D 層 (D-region)

• 高度範圍： 約 60 - 90 公里。
• 形成原因： 主要由太陽的 X 射線和極高能量的紫外線（如 Lyman-alpha）游離大氣中性分子（主要是 O2 和 N2）形成。
• 電子密度： 這一層的電子密度是所有層中最低的，特別是在夜間，幾乎消失。其電子密度隨高度的增加而增加，但總體濃度不高。
• 特徵：
  • 日夜差異大： 在白天，受太陽輻射強烈，電子密度相對較高，對中長波（AM 廣播頻率）有強烈的吸收作用，導致這些頻率的信號在日間傳播距離受限。在夜間，太陽輻射消失，游離過程停止，電子和離子迅速複合，D 層幾乎消失，因此中長波信號在夜間可以傳播更遠。
  • 對無線電波的影響： 主要表現為吸收，而不是反射。它會吸收較低頻率的無線電波。
  • 高度最低。

2. E 層 (E-region)

• 高度範圍： 約 90 - 150 公里。
• 形成原因： 主要由太陽的 X 射線（特別是波長約 10 埃的 X 射線）游離大氣中的 O2 分子形成。
• 電子密度： 具有一個相對明顯的電子密度峰值。在白天，電子密度比 D 層高，足以反射較低頻率的無線電波（如短波）。
• 特徵：
  • 日夜差異： 在白天較為明顯，夜間電子密度顯著下降，但不像 D 層那樣完全消失，仍有一部分存在。
  • 對無線電波的影響： 在白天，能夠反射中波和部分短波頻率。
  • 偶爾出現的 E_s 層 (Sporadic E)： 在某些地區和時間，E 層會出現局部、薄層但電子密度極高的區域，稱為 E_s 層。E_s 層可以反射非常高的短波頻率，甚至 VHF 頻率，導致遠距離的 VHF 通信。

3. F1 層 (F1-region)

• 高度範圍： 約 150 - 250 公里。
• 形成原因： 主要由太陽的紫外線（特別是波長約 500-1000 埃的紫外線）游離大氣中的 O 原子形成。
• 電子密度： 在白天，電子密度有一個峰值，但比 F2 層低。
• 特徵：
  • 日夜差異： 在白天存在，夜間與 F2 層合併成一個單一的 F 層。
  • 對無線電波的影響： 在白天，能夠反射短波頻率。
  • 存在時間： 主要在白天出現，通常在日出後形成，日落前消失（或與 F2 層合併）。

4. F2 層 (F2-region)

• 高度範圍： 約 250 - 1000 公里（且是電子密度最高的區域）。

本題為磁層物理的基礎概念題，考察了三個重要的磁層物理概念：三個絕熱不變量、理想 MHD 中的歐姆定律物理意義、以及磁流體動力學的「磁力線凍結」條件。

核心觀念：

1. 絕熱不變量 (Adiabatic Invariants)： 在磁場中運動的帶電粒子，如果磁場變化緩慢（相對於粒子的運動），則粒子的一些運動參數會近似保持不變。這些不變量對於描述粒子在磁場中的長期行為至關重要。
2. 理想 MHD (Ideal Magnetohydrodynamics)： 是一種描述導電流體（如電漿）與磁場相互作用的近似模型。理想 MHD 假設流體是完全導電的，沒有電阻。
3. 磁力線凍結 (Frozen-in Condition)： 是理想 MHD 的一個核心結果，它描述了在理想導電流體中，磁力線與流體之間的耦合關係。

解題步驟：

(a) 三個絕熱不變量 (Three Adiabatic Invariants)

當帶電粒子在緩慢變化的磁場中運動時，其運動的三個不同尺度（或頻率）的平均量會近似守恆。這些守恆量被稱為絕熱不變量。

1. 第一絕熱不變量：磁矩 ($\mu$)

   • 描述： 這是最著名也是最常被提及的絕熱不變量。當粒子在磁場中繞磁力線做迴旋運動（cyclotron motion）時，這個運動可以被看作是一個微小的環形電流，產生一個磁矩。如果磁場的強度在粒子迴旋的尺度上變化緩慢，則這個磁矩 $\mu$ 近似守恆。
   • 數學表示： $\mu = \frac{1}{2} m v_{\perp}^2 / B$，其中 $m$ 是粒子質量，$v_{\perp}$ 是垂直於磁場方向的速度分量，$B$ 是磁場強度。
   • 物理意義： 守恆 $\mu$ 意味著粒子在磁場強度改變時，其垂直速度分量會根據磁場強度進行調整，以保持其「磁鏡」效應。當粒子進入較強磁場時，$v_{\perp}$ 增加；進入較弱磁場時，$v_{\perp}$ 減小。這個性質是磁鏡（magnetic mirroring）現象的基礎，也是輻射帶中粒子被捕獲的關鍵。

2. 第二絕熱不變量：軌道積分 ($J$)

   • 描述： 粒子沿著磁力線方向的往返運動（漂移和磁鏡反射之間的往返）具有一個週期性。如果磁場在粒子往返運動的尺度上變化緩慢，則該軌道的某個積分量（與動量和位置相關）近似守恆。
   • 數學表示： $J = \oint v_{\parallel} dl$，其中 $v_{\parallel}$ 是平行於磁場的速度分量，$dl$ 是沿磁力線方向的微分長度，積分遍及一個往返軌道。
   • 物理意義： 守恆 $J$ 意味著粒子在沿磁力線方向的運動中，其平均動量或能量的某種形式也近似守恆。它描述了粒子在磁力線上的「拉伸」或「壓縮」行為。例如，如果磁力線被拉伸，粒子在該區域的往返時間會改變，從而影響 $J$。

3. 第三絕熱不變量：磁通量積分 ($\Phi$)

   • 描述： 當粒子在磁場中進行較大的環繞運動（例如，繞行星的磁場整體環繞，或繞地球磁層的閉合磁力線）時，如果磁場在該環繞尺度上變化緩慢，則粒子繞過的磁通量（或與之相關的某個量）近似守恆。
   • 數學表示： $\Phi = \oint \mathbf{p} \cdot d\mathbf{l}$，其中 $\mathbf{p}$ 是粒子共軛動量，$d\mathbf{l}$ 是沿著環繞軌道的微分長度。有時也用 $\oint q \mathbf{A} \cdot d\mathbf{l}$ 的形式來表示，其中 $\mathbf{A}$ 是磁向量勢。
   • 物理意義： 守恆 $\Phi$ 描述了粒子在較大尺度環繞運動中的能量或動量守恆。它與粒子在磁場中的整體漂移（如 $\mathbf{E} \times \mathbf{B}$ 漂移、梯度漂移、曲率漂移）的平均行為有關。

(b) 理想 MHD 中的歐姆定律物理意義 (Physical meaning of Ohm's law in ideal MHD equations)

理想 MHD 假設流體具有無限的導電性（電阻為零）。在這種情況下，歐姆定律的表達形式與普通導體中的不同，它包含了電場、速度和磁場的相互作用。

理想 MHD 中的歐姆定律形式為：
$\mathbf{E} + \mathbf{v} \times \mathbf{B} = 0$

其中：

• $\mathbf{E}$ 是總電場。
• $\mathbf{v}$ 是流體（電漿）的速度。
• $\mathbf{B}$ 是磁場。

物理意義：

這個方程式可以從幾個角度來理解：

1. 觀察者座標系中的電場： 對於一個靜止的觀察者（或在慣性參考系中），方程表明，在一個導電流體運動的區域，如果存在磁場 $\mathbf{B}$，則必然會產生一個感應電場 $\mathbf{E}_{ind} = -\mathbf{v} \times \mathbf{B}$。這個感應電場會抵消掉任何由其他原因（如電荷分離）引起的電場，使得在流體運動的參考系中，總的「有效」電場為零。

2. 流體內部參考系中的電場： 考慮一個跟隨流體一起運動的觀察者。在該觀察者的參考系中，流體是靜止的（$\mathbf{v} = 0$）。此時，理想歐姆定律變為 $\mathbf{E}' = 0$。這意味著在流體內部，沒有殘餘的電場。任何出現在流體內部（相對於流體是靜止的）的電場，都必須被流體的運動所「屏蔽」或「抵消」。

3. 磁力線與流體的相對運動： 方程式 $\mathbf{E} + \mathbf{v} \times \mathbf{B} = 0$ 可以改寫為 $\mathbf{E} = -\mathbf{v} \times \mathbf{B}$。這意味著，在理想導體中，電場 $\mathbf{E}$ 總是與流體的速度 $\mathbf{v}$ 和磁場 $\mathbf{B}$ 的叉乘成正比。如果我們將磁力線看作是與流體「凍結」在一起的，那麼這個關係就描述了這種「凍結」的本質：任何相對於流體的磁場變化（例如，如果流體在運動但磁場沒有，或者磁場在運動但流體沒有）都會產生一個電場，進而推動流體或磁場，使得兩者保持相對運動的「零」狀態。

4. 沒有電阻損耗： 由於電阻為零，在理想 MHD 中，電流可以在流體中無限流動而不會產生能量耗散（焦耳熱）。這個歐姆定律的形式確保了這種無損耗性。

總而言之，理想 MHD 中的歐姆定律 $\mathbf{E} + \mathbf{v} \times \mathbf{B} = 0$ 描述了在完美導體中，電場、流體速度和磁場之間的強烈耦合。它表明，在流體運動的參考系中，不存在電場，這意味著磁力線與流體高度耦合，並且沒有能量通過電阻損耗。

(c) 磁力線凍結條件 (The frozen-in condition)

「磁力線凍結」是理想 MHD 的一個核心概念，它源於理想歐姆定律 $\mathbf{E} + \mathbf{v} \times \mathbf{B} = 0$。

核心觀念

本題考查磁流體力學（MHD）中的：

1. 靜態動量平衡方程式
   $0=-\nabla p+\mathbf{J}\times\mathbf{B}$

2. 抗磁電流（diamagnetic current）與壓力梯度、磁場方向的關係
   $\mathbf{J}_{\perp}=\frac{\mathbf{B}\times\nabla p}{B^2}$

3. 磁尾（magnetotail）電流片中的跨尾電流（cross-tail current）方向，以及粒子觀點與流體觀點之間的對應。

以下採用 GSM 座標：

• $+X$：指向太陽；
• $+Y$：指向黃昏側（duskward）；
• $+Z$：指向北方。

因此，$-Y$ 方向即為由黃昏側流向黎明側（dawnward）。

---

（a）由 MHD 動量平衡推導抗磁電流

解題方法

在平衡狀態下，假設：

• 系統為靜態，$\mathbf{v}=0$；
• 忽略重力、黏滯力與其他外力；
• 壓力為各向同性標量 $p$。

MHD 動量方程式為

$$\rho\frac{d\mathbf{v}}{dt} = -\nabla p+\mathbf{J}\times\mathbf{B}.$$

平衡時 $d\mathbf{v}/dt=0$，故有

$$-\nabla p+\mathbf{J}\times\mathbf{B}=0,$$

即

$$\mathbf{J}\times\mathbf{B}=\nabla p.$$

將方程式左側與右側同時取 $\mathbf{B}$ 的叉積：

$$\mathbf{B}\times(\mathbf{J}\times\mathbf{B}) = \mathbf{B}\times\nabla p.$$

利用向量三重積公式

$$\mathbf{A}\times(\mathbf{B}\times\mathbf{C}) = \mathbf{B}(\mathbf{A}\cdot\mathbf{C}) - \mathbf{C}(\mathbf{A}\cdot\mathbf{B}),$$

得到

$$\mathbf{B}\times(\mathbf{J}\times\mathbf{B}) = \mathbf{J}B^2-\mathbf{B}(\mathbf{J}\cdot\mathbf{B}).$$

將電流分解為平行與垂直於磁場的部分：

$$\mathbf{J}=\mathbf{J}_{\parallel}+\mathbf{J}_{\perp}.$$

由於

$$\mathbf{J}_{\parallel} = \frac{\mathbf{J}\cdot\mathbf{B}}{B^2}\mathbf{B},$$

所以

$$\mathbf{J}B^2-\mathbf{B}(\mathbf{J}\cdot\mathbf{B}) = \mathbf{J}_{\perp}B^2.$$

因此

$$\mathbf{J}_{\perp}B^2 = \mathbf{B}\times\nabla p,$$

得到抗磁電流：

$$\boxed{ \mathbf{J}_{\mathrm{dia}} = \mathbf{J}_{\perp} = \frac{\mathbf{B}\times\nabla p}{B^2} }$$

其中平行電流 $\mathbf{J}_{\parallel}$ 不會產生 $\mathbf{J}\times\mathbf{B}$，因此無法由單純的靜態力平衡方程式決定。

與粒子抗磁漂移的關係

對第 $s$ 種粒子而言，壓力梯度所造成的抗磁漂移速度為

$$\mathbf{v}_{*s} = \frac{\mathbf{B}\times\nabla p_s} {q_s n_s B^2}.$$

相應的電流為

$$\mathbf{J}_{*s} = q_s n_s\mathbf{v}_{*s} = \frac{\mathbf{B}\times\nabla p_s}{B^2}.$$

將離子與電子的貢獻相加：

$$\mathbf{J}_{\mathrm{dia}} = \sum_s\mathbf{J}_{*s} = \frac{\mathbf{B}\times\nabla p}{B^2},$$

其中

$$p=\sum_s p_s.$$

這說明離子與電子的漂移速度方向相反，但由於電荷符號也相反，兩者所造成的電流方向相同。

---

（b）跨尾電流方向

流體觀點：由 $\mathbf{J}\times\mathbf{B}$ 力平衡判斷

磁尾中央為電流片，等離子體壓力在中央面附近最大，因此壓力梯度大致沿垂直方向變化：

$$\nabla p\simeq \frac{\partial p}{\partial z}\hat{\mathbf{z}}.$$

磁尾兩側的磁場方向相反。典型磁尾幾何為：

本題考察地球磁層頂的激波/不連續結構，並結合了太陽風參數計算激波速度。

核心觀念：

1. 磁層頂 (Magnetopause)： 是太陽風電漿與地球磁場之間相互作用形成的邊界。
2. 激波 (Shock Wave)： 當流體（在此為太陽風電漿）的速度超過其波速（例如阿爾芬波速 $v_A$ 或快磁聲波速 $v_{ms}$）時，會產生一種不連續的結構，能量會耗散。
3. IMF $B_z$ 方向： 太陽風攜帶的行星際磁場（IMF）的南北分量 ($B_z$) 對磁層頂的結構和耦合有重要影響。
   • 南向 $B_z$ (Southward $B_z$)： IMF 的 $B_z$ 分量為負時，IMF 與地球磁場會發生磁重聯（magnetic reconnection），導致能量和動量更容易從太陽風進入磁層。
   • 北向 $B_z$ (Northward $B_z$)： IMF 的 $B_z$ 分量為正時，磁重聯較難發生，能量耦合效率較低。
4. 連續性方程和動量方程： 用於描述流體參數（密度、速度、壓力、磁場）如何變化。
5. 激波速度計算： 根據激波前後的流體參數變化，利用守恆定律（質量、動量、能量）來計算激波的速度。

解題步驟：

(a) 磁層頂的激波/不連續結構在北向和南向 IMF $B_z$ 條件下的區別

這裡的問題似乎有點誤解。通常討論的是太陽風與磁層頂的「相互作用」或「邊界結構」，而不是磁層頂本身產生「激波」。激波（或稱弓形激波 Bow Shock）是太陽風在進入地球磁場區域之前，由於其速度超過當地波速而產生的不連續結構，它位於磁層頂的前方。

假設題目指的是太陽風與磁層頂的「邊界結構」或「不連續性」，以及太陽風本身是否可能在磁層頂附近產生激波。

• 弓形激波 (Bow Shock)： 這是太陽風在進入地球磁場區域之前產生的激波。它是一個相對穩定的結構，位於磁層頂前方數十倍地球半徑的距離。弓形激波的存在與太陽風的馬赫數（Mach number）有關。

  • IMF $B_z$ 的影響： IMF 的方向（特別是 $B_z$）主要影響的是磁層頂的耦合效率和磁層的動力學，而不是弓形激波本身的產生。弓形激波的產生主要取決於太陽風的超聲速（或超阿爾芬速）性質。

• 磁層頂的結構與 IMF $B_z$ 的關係：

  • 南向 $B_z$ (Southward IMF $B_z$)：

    • 磁重聯： 當 IMF 的 $B_z$ 為負時，IMF 與地球磁場在磁層頂會發生磁重聯。這使得太陽風的能量、動量和粒子可以高效地進入磁層。
    • 結構： 磁層頂的形狀會變得更為「動態」和「開放」。磁尾會被拉伸，並可能發生磁尾的重新連接，導致磁暴。磁層頂的邊界可能變得更為不規則。
    • 不連續性： 磁重聯本身是一個高度耗散的過程，在磁重聯點附近會出現電場、電流和粒子速度的急劇變化，這可以被視為一種局部的不連續性。

  • 北向 $B_z$ (Northward IMF $B_z$)：

    • 無磁重聯或低效耦合： 當 IMF 的 $B_z$ 為正時，IMF 與地球磁場在磁層頂的直接磁重聯較難發生（或效率很低）。能量和動量從太陽風進入磁層的效率較低。
    • 結構： 磁層頂的形狀相對更為穩定和對稱。磁尾的開放性較差。
    • 可能的不連續性： 在某些情況下，可能發生所謂的「極向吸積」（polar erosion）或「開放磁層」的形成，但這些過程與南向 $B_z$ 引起的劇烈磁重聯和磁暴相比，通常較為緩和。有時也會討論「磁層頂的噴流」（magnetopause jets），它們可能與磁場傾斜和局部重聯有關。

總結 (a)：
題目可能誤解了「激波」的產生位置。弓形激波是太陽風產生的，位於磁層頂前方。
IMF $B_z$ 的方向主要影響的是磁層頂與太陽風的「耦合」方式和磁層的「動力學」：

• 南向 $B_z$： 導致太陽風與磁層發生高效的磁重聯，過程動態，磁層頂邊界可能不規則，磁層內能量和動量注入增加，可能引發磁暴。磁重聯點是局部的不連續結構。
• 北向 $B_z$： 耦合效率低，磁層頂結構相對穩定和對稱，磁層內能量和動量注入較少。

(b) 計算激波速度

這裡的題目是明確的，要求計算一個「激波」的速度，給定了激波前後的密度和太陽風速度變化。這是一個典型的激波物理問題，可以使用激波前後的守恆定律來解決。

我們假設這是一個一維、穩態的激波，並且太陽風是一種電漿。我們需要使用激波前後的質量守恆、動量守恆和能量守恆定律（或簡化形式）。