宇宙微波背景的聲學樂譜：讀懂功率譜與聲學峰

為什麼宇宙最古老的光裡，藏著一張比針孔還小的「聲學樂譜」？

如果你已經知道宇宙微波背景（Cosmic Microwave Background, CMB）是宇宙嬰兒期遺留下來、溫度約 $2.725\,\mathrm{K}$ 的均勻輻射，那麼接下來這個問題會讓你重新審視它：為什麼科學家看著一張「幾乎完全均勻」的天空地圖，卻能讀出宇宙含有多少暗物質、空間是平的還是彎的、甚至宇宙誕生後第一秒發生了什麼？

答案藏在那張地圖上百萬分之一等級的微小溫度漣漪裡。這些漣漪不是雜訊，而是早期宇宙電漿（plasma）中一場持續了近四十萬年的「聲波」所凝固下來的快照。本篇我們不再重述 CMB 是什麼，而是要拆解它的功率譜（power spectrum）——這張「聲學樂譜」如何被測量、如何被讀懂，以及它為什麼是現代宇宙學最精密的量尺。

從溫度地圖到功率譜：把天空展開成「音階」

入門篇談的是 CMB 的整體均勻性與微小擾動。進階的關鍵一步，是學會用球諧函數（spherical harmonics）把整張球面溫度地圖分解。

天空是一個球面，任何球面上的函數都可以展開成球諧函數 $Y_{\ell m}$ 的疊加。CMB 的溫度漲落 $\Delta T / T$ 在天空方向 $\hat{n}$ 上可寫成：

$$ \frac{\Delta T}{T}(\hat{n}) = \sum_{\ell=2}^{\infty} \sum_{m=-\ell}^{\ell} a_{\ell m}\, Y_{\ell m}(\hat{n}) $$

這裡的多極矩指標 $\ell$ 對應的是天空上的角尺度：粗略地說，$\ell$ 越大代表越精細的結構。一個實用的對應關係是

$$ \theta \approx \frac{180^\circ}{\ell}. $$

所以 $\ell \approx 200$ 對應約 $1^\circ$ 的斑塊——大約是滿月寬度的兩倍。

由於早期宇宙的擾動在統計上是各向同性（isotropic）且接近高斯（Gaussian）的，所有物理資訊濃縮在每個 $\ell$ 的變異數裡，也就是角功率譜（angular power spectrum）：

$$ C_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{\ell} \left| a_{\ell m} \right|^2 . $$

天文學家習慣畫的是 $\mathcal{D}_\ell \equiv \dfrac{\ell(\ell+1)}{2\pi} C_\ell$，因為這個組合在平坦的初始條件下大致平緩，能讓那一連串「聲學峰」清楚浮現。功率譜橫軸是 $\ell$、縱軸是 $\mathcal{D}_\ell$，那一排起伏的峰，就是我們要解讀的樂譜。

聲學峰的物理：光子—重子流體的駐波

為什麼會有「峰」？這是 CMB 物理最精華的部分。

在復合（recombination）之前，宇宙是由光子、自由電子與重子（baryon，主要是質子）緊密耦合成的單一光子—重子流體（photon-baryon fluid）。暗物質則不參與電磁交互作用，只透過重力提供位能井（potential well）。

想像一個由暗物質造成的重力位能井。重子流體被重力拉進井裡壓縮，但光子的輻射壓（radiation pressure）會把流體往外推。壓縮 → 反彈 → 再壓縮，這就形成了聲學振盪（acoustic oscillation）——本質上是早期宇宙裡的聲波。

這些聲波的傳播速度是流體的聲速（sound speed） $c_s$。對純光子流體 $c_s = c/\sqrt{3}$；重子的存在會讓它略微下降：

$$ c_s = \frac{c}{\sqrt{3\,(1 + R)}}, \qquad R \equiv \frac{3\rho_b}{4\rho_\gamma}. $$

這裡 $R$ 是重子與光子的能量密度比例。當宇宙在 $z \approx 1100$ 復合、電子被質子捕獲成中性氫時，光子突然「解耦」自由飛行，振盪瞬間凍結。凍結時剛好處於最大壓縮或最大膨脹的模態，就在功率譜上留下一個峰。

• 第一峰：對應在整個宇宙年齡裡剛好完成「一次壓縮」的模態，是振幅最大的基頻（fundamental mode）。
• 第二峰、第三峰……：對應完成更多次半振盪的諧波（harmonics）。

這就是為什麼 CMB 功率譜長得像一段有基頻與泛音的樂譜——它真的是一場聲學演奏的最後一幀。

每個峰在告訴我們什麼參數

聲學峰不只是漂亮，每個特徵都精準地對應一個宇宙學參數。這是 CMB 成為「精密宇宙學」基石的原因。

第一峰的位置 → 空間幾何（曲率）。 聲波在復合前能傳播的最大距離叫聲學視界（sound horizon） $r_s$，這是一把已知物理長度的「標準尺（standard ruler）」。我們觀測到它在天空上張開的角度 $\theta_s$。已知尺長、量到角度，就能反推光走過的幾何：

$$ \theta_s = \frac{r_s}{D_A}, $$

其中 $D_A$ 是到最後散射面的角直徑距離。如果空間是平的，第一峰落在 $\ell \approx 220$；若空間正彎曲，同樣的尺會顯得更大（峰往左移），負彎曲則相反。Planck 衛星的測量結果是峰精準落在 $\ell \approx 220$，告訴我們宇宙在觀測精度內是空間平坦的。

第二峰相對第一峰的高度 → 重子密度。 重子讓流體變「重」，加強壓縮相、削弱膨脹相，造成奇數峰（壓縮）被增強、偶數峰（膨脹）被壓低的不對稱。第二峰（偶數）相對第一峰（奇數）被壓得越低，重子越多。這給出 $\Omega_b h^2 \approx 0.0224$。

第三峰的高度 → 暗物質密度。 暗物質不感受輻射壓，但它的重力位能井決定了振盪的「深度」。第三峰的相對高度對總物質密度特別敏感，給出 $\Omega_c h^2 \approx 0.120$——這是暗物質存在最乾淨的證據之一，完全獨立於星系旋轉曲線。

高 $\ell$ 的衰減尾巴 → 擴散阻尼。 復合不是瞬間發生的，光子在解耦前後會發生隨機行走（random walk），把小尺度的結構抹平。這叫席爾克阻尼（Silk damping），使功率譜在 $\ell \gtrsim 1000$ 之後指數衰減。

動手算一下：第一峰為什麼在 $\ell \approx 220$？

讓我們用標準尺的概念做個量級估算，感受一下尺度。

聲學視界在復合時的物理尺度約為 $r_s \approx 147\,\mathrm{Mpc}$（共動距離），這是一個由早期宇宙物理嚴格算出的數字。到最後散射面的共動距離約為 $D_C \approx 14{,}000\,\mathrm{Mpc}$。

第一峰對應的角尺度是聲學視界張開的角度。對共動量而言：

$$ \theta_s \approx \frac{r_s}{D_C} \approx \frac{147}{14{,}000} \approx 0.0105\ \text{弧度} \approx 0.60^\circ. $$

把它換成多極矩。第一峰實際對應的是約半個波長的振盪，習慣上對應關係約為 $\ell_1 \approx \pi / \theta_s$：

$$ \ell_1 \approx \frac{\pi}{0.0105} \approx 300. $$

這個粗估落在正確的量級。精確的數值計算（考慮投影效應、聲速隨時間變化等）會把它修正到 $\ell_1 \approx 220$。重點是：一把 $147\,\mathrm{Mpc}$ 的尺，在 $140$ 億光年外張開不到一度——而我們居然能測出它在天空上的精確位置，這正是現代宇宙學的驚人之處。

偏振：CMB 的第二層資訊

溫度只是 CMB 的一半故事。光子在最後散射時，若周圍的光子分布有四極矩各向異性（quadrupole anisotropy），湯姆森散射（Thomson scattering）會讓散射出來的光帶有線偏振（linear polarization）。

偏振模式可以分解成兩種獨立的幾何型態：

• E 模（E-mode）：偏振方向呈徑向或切向的梳狀圖樣，沒有旋性（curl-free）。E 模由純量擾動（密度漲落）產生，已被 Planck 與 ACT、SPT 等地面實驗精確測量，與溫度功率譜高度一致，是標準模型的強力交叉驗證。

• B 模（B-mode）：帶有「旋轉」的圖樣（curl-like）。純量密度擾動無法在線性階產生 B 模——這是它珍貴的原因。B 模有兩個來源： 1. 重力透鏡（gravitational lensing）：前景大尺度結構把 E 模扭曲成 B 模，已被偵測。 2. 原初重力波（primordial gravitational waves）：暴脹（inflation）若真實發生，會產生張量擾動（tensor perturbations），在大角尺度留下原初 B 模。

偵測到原初 B 模，等於直接看到暴脹時期的時空漣漪——這是當前 CMB 觀測的聖杯。它的振幅由張量對純量比（tensor-to-scalar ratio） $r$ 衡量。目前的上限約為 $r \lesssim 0.03$，BICEP/Keck、Simons Observatory 與未來的 LiteBIRD 衛星正在把這條線往下逼。

不只是峰：兩個容易被忽略的進階特徵

重子聲學振盪（BAO）的「回聲」。 CMB 裡的聲學視界 $r_s$ 不只留在天空上。同一把尺也凍結在重子的空間分布裡，在後續數十億年隨星系一起被重力放大，形成今天星系兩兩相距約 $150\,\mathrm{Mpc}$ 的統計偏好。這就是 BAO，是 CMB 標準尺在低紅移宇宙的「回聲」，讓我們能跨越宇宙歷史校準暗能量。

積分薩克斯—沃爾夫效應（Integrated Sachs-Wolfe, ISW）。 在暗能量主導的近代宇宙，光子穿越正在「變淺」的重力位能井時會淨得到能量（藍移多於紅移），在大角尺度（低 $\ell$）給功率譜一點額外貢獻。把 CMB 大尺度漲落與低紅移星系分布做交叉相關，能獨立確認暗能量的存在。

重點回顧

• CMB 的物理資訊不在那張地圖「看起來」如何，而在它的角功率譜 $C_\ell$：把溫度地圖用球諧函數展開後，每個多極矩 $\ell$ 的變異數。
• 功率譜上的聲學峰是早期光子—重子流體在暗物質位能井中振盪、於復合時凍結所留下的駐波快照。
• 峰的位置測曲率（宇宙是平的）、奇偶峰高度不對稱測重子密度、第三峰高度測暗物質密度、高 $\ell$ 衰減反映席爾克阻尼。
• 聲學視界 $r_s \approx 147\,\mathrm{Mpc}$ 是一把標準尺，在天空上張開約 $1^\circ$（第一峰 $\ell \approx 220$），也以 BAO 的形式回聲在星系分布裡。
• 偏振 E/B 模提供第二層資訊；原初 B 模若被偵測，將是暴脹存在的直接證據（以 $r$ 衡量，目前 $r \lesssim 0.03$）。

深入探討（研究所視角）

要把上述定性圖像變成 Planck 那條與資料吻合到誤差棒以內的理論曲線，需要求解線性化的愛因斯坦—玻茲曼方程組（linearized Einstein-Boltzmann equations）。核心是對每個物種（光子、微中子、重子、暗物質）寫出其相空間分布擾動的演化方程，並與度規擾動耦合。光子的部分要把分布函數對方向做多極展開 $\Theta_\ell$，形成一條無窮的玻茲曼層級（Boltzmann hierarchy）：

$$ \dot{\Theta}_\ell - k\!\left(\frac{\ell}{2\ell-1}\Theta_{\ell-1} - \frac{\ell+1}{2\ell+3}\Theta_{\ell+1}\right) = -\dot{\tau}\,\Theta_\ell + (\text{源項}), $$

其中 $\dot{\tau}$ 是湯姆森散射的微分光學深度。實務上由 CAMB 或 CLASS 這類 Boltzmann solver 數值積分，再透過線視積分法（line-of-sight integration，Seljak-Zaldarriaga 1996）把震盪源投影成今日觀測的 $C_\ell$。

兩個值得深思的精細點：

第一，緊耦合近似（tight-coupling approximation）的崩潰正是 Silk damping 的物理根源。在 $\dot{\tau}$ 很大時光子與重子近乎完美耦合，但耦合不是無限緊；展開到光子平均自由程的一階，阻尼尺度約為 $k_D^{-2} \sim \int \dfrac{1}{\dot\tau}\,(\cdots)\,d\eta$，把小尺度功率以 $e^{-(k/k_D)^2}$ 抹平。

第二，宇宙變異數（cosmic variance）的根本極限。我們只有一個宇宙、一片天空，每個 $\ell$ 只有 $2\ell+1$ 個獨立的 $a_{\ell m}$ 樣本。即使儀器完美無噪，$C_\ell$ 的估計仍有不可消除的相對不確定度：

$$ \frac{\Delta C_\ell}{C_\ell} = \sqrt{\frac{2}{2\ell+1}}. $$

這對低 $\ell$ 特別嚴重——這也是為什麼 CMB 大尺度上那幾個著名的「異常」（如低四極矩、半球不對稱）難以判定是新物理還是統計巧合：在 $\ell=2$，宇宙變異數本身就高達約 $63\%$。下一步若要超越溫度與 E 模的資訊上限，物理界正把希望寄託在原初 B 模偵測、CMB 透鏡重建出的物質分布、以及與 BAO／弱重力透鏡的聯合分析上，逼近暴脹能標與微中子質量總和 $\sum m_\nu$ 的下一位有效數字。