塞曼效應與星際磁場的光譜學測量原理
塞曼效應 — 原子或分子在外磁場中譜線分裂為多個偏振分量的現象 — 是目前唯一能夠直接測量星際磁場強度的方法。當一個具有總角動量量子數 J 的原子態置於磁場 B 中時，其能級分裂為 2J+1 個子能級，能級位移為 ΔE = g_J·μ_B·m_J·B，其中 g_J 為朗德 g 因子，μ_B 為玻爾磁子，m_J 為磁量子數。譜線的塞曼分裂圖樣取決於上下能級的 g 因子和選擇規則，最常見的是正常塞曼效應（分裂為三條線：π 線和 σ± 線）和反常塞曼效應（更複雜的分裂圖樣）。
關鍵在於這些塞曼分量的偏振特性：π 分量的電場振盪平行於磁場方向，σ± 分量的電場振盪垂直於磁場方向並以相反方向圓偏振。通過測量譜線的圓偏振斯托克斯 V 參數，可以直接確定視線方向的磁場分量 B_los，其關係為 V ∝ dI/dν·B_los·cos θ，其中 θ 為磁場與視線的夾角。
斯托克斯偏振光譜學
塞曼測量的核心實驗技術是斯托克斯偏振光譜學。四個斯托克斯參數 I、Q、U、V 完整描述了電磁輻射的偏振狀態：I 為總強度，Q/U 為線偏振，V 為圓偏振。對於塞曼測量，主要關注 V 譜（圓偏振）與 I 譜的頻率導數之間的相關性：V(ν) ∝ B_los·cos θ · dI/dν。這種關係的物理直觀是：磁場引入的頻移在譜線翼部產生了過量的圓偏振。
在實際觀測中，V 譜的信號通常極其微弱（ΔV/I ~ 10⁻³ 至 10⁻⁵），需要長時間積分和高頻譜解析度。現代射電望遠鏡（如 ALMA、VLA、Arecibo）將圓偏振接收機與數位相關器結合，頻道寬度可達亞 kHz 級別，使 GHz 頻率的 OH、CN 和 H₂O 分子塞曼測量成為可能。
從恆星形成到星系動力學
星際磁場在星際介質的幾乎所有物理過程中都扮演著關鍵角色。在恆星形成的初始階段，磁場提供了抵抗重力塌縮的額外壓力（磁臨界質量概念）：如果分子雲核的質量-磁通量比超過臨界值 (M/Φ)_crit = 1/(2π√G)，塌縮才能發生。塞曼測量直接約束這一比值，揭示磁場在恆星形成效率中的調節作用。
在更大的星系尺度上，磁場通過與宇宙射線和湍流的交互作用影響星際介質的能量平衡和多相結構。磁場的壓力張量各向異性在螺旋星系中驅動了從旋臂向旋臂間的氣體流動，而磁場的再連接則是非熱粒子加速的可能機制。塞曼測量與偏振塵埃輻射和法拉第旋轉測量互補，共同繪製了從分子雲核心到星系暈的多尺度磁場結構。
塞曼光譜數據的擬合模擬
import numpy as np from scipy.optimize import curve_fit from scipy.special import voigt_profile class ZeemanSpectrumFitter: # Stokes V spectrum fitting for magnetic field measurement def __init__(self, freq_rest, g_lande=1.0): self.nu0 = freq_rest self.g = g_lande self.z = 1.4 # Zeeman splitting factor [Hz/μG] def _stokes_v_model(self, velocity, amp, v_lsr, sigma, B_los): # Derivative-of-Gaussian model for Stokes V dv = velocity - v_lsr I_model = amp * np.exp(-dv**2 / (2 * sigma**2)) dI_dv = -amp * dv / sigma**2 * np.exp(-dv**2 / (2 * sigma**2)) delta_nu = self.z * self.g * B_los delta_v = -delta_nu / self.nu0 * 2.9979e5 return delta_v * dI_dv, I_model def fit(self, velocity, I_obs, V_obs): def forward_model(v, amp, v_lsr, sigma, B_los): V_model, I_model = self._stokes_v_model(v, amp, v_lsr, sigma, B_los) return V_model amp0 = np.max(I_obs); v0 = velocity[np.argmax(I_obs)] sigma0 = np.std(velocity[:10]) * 3; B0 = 50 # initial B guess [μG] popt, _ = curve_fit(forward_model, velocity, V_obs, p0=[amp0, v0, sigma0, B0]) return popt[3], popt # B_los best-fit value # Simulated OH 1667 MHz Zeeman observation vel = np.linspace(-10, 10, 256) I_syn = 5 * np.exp(-vel**2 / (2 * 1.5**2)) + 0.1 * np.random.randn(256) V_syn = 0.01 * np.gradient(I_syn, vel[1]-vel[0]) + 0.005 * np.random.randn(256) fitter = ZeemanSpectrumFitter(freq_rest=1667.35903e6, g_lande=0.28) B_est, params = fitter.fit(vel, I_syn, V_syn) print(f"Estimated B_los = {B_est:.1f} μG")
磁場宇宙學的未來
星際和星系間磁場的起源是現代天體物理學最深刻的未解問題之一。銀河系的微高斯級磁場是否來源於早期宇宙中的原初種子場（通過 Battery 機制或暴脹產生），還是完全由天體物理發電機過程（星系旋轉+湍流+宇宙射線）後期放大？新一代射電望遠鏡 — 包括平方公里陣列及其先導項目 — 將通過大規模偏振巡天進行數十萬個源的塞曼和法拉第旋轉測量，從統計上區分這兩種場景。此外，對高紅移星系的塞曼觀測將首次揭示宇宙正午時期（z ~ 2）的星系磁場演化歷史，為理解磁場在星系形成與演化中的作用提供決定性約束。
本文內容僅供學術研究參考。塞曼光譜學測量的靈敏度受限於望遠鏡口徑和觀測時間，文中引用的磁場數值為典型值，具體天體的磁場強度可能因環境而異。