星系如何被組裝：從密度漲落到宇宙網

為什麼宇宙不是一鍋均勻的湯？

入門篇告訴我們：宇宙在膨脹，星系彼此遠離，一切始於大霹靂（Big Bang）。但如果你停在那裡，會留下一個尖銳的問題——如果早期宇宙幾乎完美均勻，今天為什麼會有星系、星系團、巨大的空洞與纖維狀結構？ 均勻的宇宙不會自己長出螺旋星系，就像一杯攪拌均勻的牛奶不會自己分出奶塊。

進階篇要回答的，正是這個「從近乎無到萬有」的劇本。我們不再問「星系是什麼」，而是問「星系如何被組裝出來」：重力如何放大微小的密度漲落、暗物質（dark matter）如何先搭好骨架、氣體如何冷卻塌縮、而星系又如何透過回饋（feedback）調節自己的成長。這是一條把宇宙學與星系物理縫在一起的主線。

種子：宇宙微波背景裡的十萬分之一

一切結構的種子，藏在宇宙微波背景輻射（Cosmic Microwave Background, CMB）的溫度漲落裡。CMB 各方向的溫度幾乎一致（約 $2.725\,\mathrm{K}$），但存在極微小的起伏：

$$ \frac{\Delta T}{T} \sim 10^{-5} $$

這個 $10^{-5}$ 量級的不均勻，對應到宇宙誕生約 38 萬年（再結合時期，recombination）時，物質密度的微小起伏 $\delta \equiv \delta\rho / \bar\rho$。看似微不足道，卻是整座結構大廈的地基。

重力是一種正回饋機制：稍微密一點的區域，重力略強，會吸引更多物質，使它更密——這就是重力不穩定性（gravitational instability）。在物質主導的膨脹宇宙中，線性階段的密度對比隨宇宙標度因子（scale factor）$a$ 線性增長：

$$ \delta(a) \propto a \quad (\text{物質主導期、線性區}) $$

從 $a \approx 1/1100$（再結合時）長到 $a = 1$（今天），這個成長因子約 1000 倍。但 $10^{-5} \times 1000 = 10^{-2}$，仍遠小於 1，光靠普通物質（重子，baryon）根本來不及塌縮成星系。這正是暗物質非登場不可的理由。

暗物質：先搭骨架的隱形工人

重子在再結合之前與光子緊密耦合，光子壓力像彈簧一樣抵抗塌縮，使重子的密度漲落無法成長（只能聲學振盪，產生 CMB 上的聲學峰）。但暗物質不帶電、不與光子作用，它的漲落早在再結合之前就開始默默成長。

於是劇本是這樣：暗物質先塌縮成一個個重力位能井，稱為暗物質暈（dark matter halo）。等到再結合、光子壓力鬆綁後，重子才「掉進」這些早已挖好的位能井裡。換句話說，暗物質暈是舞台，重子（最終的恆星與氣體）是後來才上場的演員。

暗物質的另一個關鍵性質是「冷」（cold）：它在脫耦時運動速度遠低於光速。冷暗物質（Cold Dark Matter, CDM）的小尺度漲落不會被自由流動（free-streaming）抹平，因此小結構先形成、再合併成大結構——這就是 ΛCDM 模型的階層式成長（hierarchical assembly）。矮星系般大小的暈先誕生，再透過併合（merging）逐步堆疊成大星系與星系團。這與早期「先形成大結構再碎裂」的由上而下圖像（top-down）相反，而觀測（如高紅移小星系豐富、星系團晚近才組裝完成）強烈支持由下而上（bottom-up）。

從線性到非線性：球塌縮的臨界值

當 $\delta$ 還很小（$\ll 1$）時，演化是線性的、可解析。但當某區域密度漲到一定程度，重力勝過膨脹，該區域會停止隨宇宙膨脹、轉而塌縮——進入非線性階段。

最簡潔的模型是球塌縮模型（spherical collapse）。把一團物質想成一個迷你的封閉宇宙：它先隨大宇宙膨脹，到達最大半徑（turnaround）後反轉、塌縮。理論給出一個漂亮的數字：當線性外推的密度對比達到

$$ \delta_c \approx 1.686 $$

時，對應的真實（非線性）區域已經塌縮並透過弛豫（relaxation）達到准平衡，形成一個束縛的暈。維里定理（virial theorem）告訴我們，束縛後系統的動能與位能滿足 $2K + U = 0$，暈會穩定在一個平均密度約為當時宇宙背景密度 $\sim 178$ 倍（在愛因斯坦–德西特宇宙下）的狀態。

這個 $\delta_c \approx 1.686$ 是現代結構形成統計理論（Press–Schechter 形式論）的心臟：它把「線性世界裡哪些區域會塌縮」轉成一個可計算的門檻，讓我們能預測不同質量暗物質暈的數量分布。

氣體冷卻：星系真正點亮的那一步

有了暗物質暈這個位能井，重子氣體掉進去後會被重力加熱到接近維里溫度（virial temperature）：

$$ T_{\mathrm{vir}} \approx \frac{\mu m_p}{2 k_B} \frac{G M}{R} $$

其中 $\mu$ 為平均分子量、$m_p$ 為質子質量、$M$ 與 $R$ 為暈的質量與半徑。對銀河系大小的暈，$T_{\mathrm{vir}}$ 高達百萬 K。但熱氣體不會塌縮成恆星——它的壓力支撐著自己。要形成恆星，氣體必須先冷卻、把熱量輻射掉、失去壓力支撐後才能向中心沉降。

冷卻靠的是輻射過程：高溫下的軔致輻射（bremsstrahlung）、稍低溫下的譜線冷卻（line cooling，金屬與氫氦的電子躍遷）。是否能有效冷卻，取決於冷卻時間 $t_{\mathrm{cool}}$ 與重力自由落體時間 $t_{\mathrm{ff}}$ 的競賽：

• 若 $t_{\mathrm{cool}} < t_{\mathrm{ff}}$：氣體來得及冷卻、塌縮、成星——星系得以成長。
• 若 $t_{\mathrm{cool}} > t_{\mathrm{ff}}$：氣體維持熱的、瀰漫的暈，難以成星。

這個對比解釋了一個深刻的觀測事實：星系的恆星質量並非隨暗物質暈質量等比例增加。在 $\sim 10^{12}\,M_\odot$ 暈附近，恆星形成效率最高；質量更小或更大的暈，效率都明顯下降。低質量端與高質量端各有不同的「煞車」，我們接著就談這個煞車——回饋。

看一個例子：銀河系尺度的冷卻競賽

考慮一個維里溫度約 $10^6\,\mathrm{K}$、電子數密度約 $n \sim 10^{-3}\,\mathrm{cm^{-3}}$ 的熱暈氣體。在此溫區，冷卻函數約 $\Lambda(T) \sim 10^{-23}\,\mathrm{erg\,cm^3\,s^{-1}}$。冷卻時間的量級估計為

$$ t_{\mathrm{cool}} \sim \frac{3 n k_B T}{n^2 \Lambda(T)} = \frac{3 k_B T}{n\,\Lambda(T)} $$

代入（$k_B T \approx 1.38\times10^{-16}\times10^{6}\,\mathrm{erg} \approx 1.4\times10^{-10}\,\mathrm{erg}$）：

$$ t_{\mathrm{cool}} \sim \frac{3 \times 1.4\times10^{-10}}{10^{-3}\times10^{-23}}\,\mathrm{s} \approx 4\times10^{16}\,\mathrm{s} \approx 1.3\,\mathrm{Gyr} $$

約十幾億年，與星系演化的時間尺度可比。這說明在這類暈中，內側緻密氣體可以及時冷卻成星，而外側稀薄氣體則維持熱暈——這也正是現代模擬中「冷流（cold flow）」與「熱暈吸積」兩種供氣模式的分水嶺。

回饋：星系如何給自己踩煞車

如果只有冷卻與重力，模擬會做出比真實宇宙多太多的恆星——這就是著名的「過度冷卻問題」（overcooling problem）。真實星系一定有把氣體加熱、推出去、阻止成星的機制，這些統稱回饋：

• 恆星回饋（stellar feedback）：大質量恆星的紫外輻射、恆星風，以及超新星（supernova）爆炸，注入能量與動量，把冷氣體加熱、吹散，甚至驅動星系尺度的外流（galactic wind）。這主導低質量星系的調節——它們淺淺的位能井留不住被加熱的氣體。
• 活躍星系核回饋（AGN feedback）：星系中心的超大質量黑洞（supermassive black hole）吸積物質時釋放巨大能量，以噴流（jet）與輻射方式加熱周圍氣體，抑制冷卻。這主導大質量星系——尤其解釋了為何巨大的橢圓星系與星系團中心星系「又老又紅」、幾乎停止造星。

回饋的存在，使星系成長成為一個自我調節（self-regulation）的系統：氣體冷卻成星 → 恆星與黑洞釋放能量 → 加熱氣體、抑制冷卻 → 成星變慢 → 能量釋放變少 → 氣體又能冷卻……如此往復。正是這套負回饋，把恆星形成效率壓在低水準，並雕刻出我們觀測到的星系質量函數與各種尺度關係。

尺度關係：星系不是隨機長出來的

星系的物理參數之間存在驚人地緊密的相關性，這些尺度關係（scaling relations）是檢驗形成理論的試金石：

• 塔利–費雪關係（Tully–Fisher relation）：螺旋星系的光度（或重子質量）與旋轉速度 $v$ 的關係近似 $L \propto v^{4}$。旋轉速度由暗物質暈主宰，光度由恆星主宰——這條關係把可見的恆星與隱形的暈牢牢綁在一起。
• 基本面（Fundamental Plane）：橢圓星系的有效半徑、表面亮度與速度彌散度落在一個薄薄的平面上，反映維里平衡與恆星族群的規律性。
• 質量–金屬度關係：星系恆星質量愈大，氣體金屬度（metallicity，指比氦重的元素豐度）愈高。這是因為大星系位能井深、外流帶不走金屬，化學增豐（chemical enrichment）得以累積。
• $M$–$\sigma$ 關係：中心超大質量黑洞質量與星系核球的速度彌散度 $\sigma$ 強相關（約 $M_{\mathrm{BH}} \propto \sigma^{4\text{–}5}$）。黑洞質量遠小於星系，卻與整個核球「對話」——這是 AGN 回饋共同演化的化石證據。

這些關係不是巧合，而是重力、冷卻、回饋三方角力後的平衡態指紋。

宇宙網：星系住在哪裡

把視野拉到最大尺度，星系並非隨機散布，而是排列在一張宇宙網（cosmic web）上：物質聚集成纖維（filament），纖維交會處是星系團（cluster），纖維之間是幾乎空無一物的空洞（void）。這個圖像由 ΛCDM 初始漲落直接演化而來，並在大型星系巡天（如 SDSS）與數值模擬中得到漂亮驗證。

星系所處的環境會影響其命運：星系團這種高密度環境裡，星系容易被剝離氣體（如沖壓剝離，ram-pressure stripping）、被頻繁交互作用攪動，因而傾向變紅、停止造星；空洞邊緣的孤立星系則保有更多氣體、持續成長。這就是環境（nature vs. nurture）對星系演化的塑造——同樣的暈，住在不同社區，會走上不同人生。

動手算一下：暈的特徵密度與紅移

球塌縮給出一個重要結論：暈塌縮時，其平均密度凍結在「當時」宇宙背景密度的某個固定倍數。在簡單的愛因斯坦–德西特近似下，倍數約為 $\Delta \approx 178$（常近似取 200）。由於物質密度隨膨脹稀釋為 $\bar\rho_m(z) \propto (1+z)^3$，愈早（高紅移 $z$）形成的暈，特徵密度愈高：

$$ \rho_{\mathrm{halo}}(z) \approx \Delta \cdot \bar\rho_{m,0}\,(1+z)^3 $$

比較 $z=3$ 與 $z=0$ 形成的暈：

$$ \frac{\rho_{\mathrm{halo}}(z=3)}{\rho_{\mathrm{halo}}(z=0)} = (1+3)^3 = 64 $$

早期形成的暈密度高出 64 倍。這正是為什麼高紅移星系往往更緻密、更小——它們是在更密的宇宙裡組裝出來的。今天巨大鬆散的星系，許多是後來透過乾併合（dry merger，少氣體的併合）把這些緻密「種子」堆疊、膨脹而成。這個量級估計，把抽象的膨脹宇宙與具體的星系形態連了起來。

重點回顧

• 星系結構的種子是 CMB 中 $\Delta T/T \sim 10^{-5}$ 的微小漲落；重力不穩定性放大它，但線性成長 $\delta \propto a$ 不足以單靠重子塌縮成星，暗物質必須先行。
• 冷暗物質決定了階層式（由下而上）成長：小暈先形成，再併合成大結構；球塌縮給出臨界線性密度 $\delta_c \approx 1.686$ 作為塌縮門檻。
• 重子掉入暈後被加熱到維里溫度，必須先冷卻（$t_{\mathrm{cool}} < t_{\mathrm{ff}}$）才能成星；冷卻效率隨暈質量變化，造就恆星形成效率的峰值。
• 回饋（恆星／超新星主導低質量端，AGN 主導高質量端）使星系自我調節，解決過度冷卻問題，並雕刻出星系質量函數與尺度關係。
• 星系鑲嵌在宇宙網上，環境與併合史共同決定其形態、顏色與命運；高紅移暈密度依 $(1+z)^3$ 升高，解釋早期星系的緻密性。

深入探討（研究所視角）

從統計形式論到完整模擬。 Press–Schechter（1974）及其擴展（Extended Press–Schechter, EPS）把球塌縮的 $\delta_c$ 與初始功率譜結合，解析地預測暈質量函數 $\dfrac{\mathrm{d}n}{\mathrm{d}M}$。其核心是把平滑後密度場的方差 $\sigma^2(M)$ 與門檻 $\delta_c$ 比較，計算超越門檻的體積占比。Bond 等人（1991）以「excursion set（漫遊集）」與布朗運動軌跡的首次穿越問題，給出更嚴謹的推導，並自然處理「雲中雲」問題；Sheth–Tormen 則引入橢球塌縮修正，更貼合 N-body 模擬。這條線是把宇宙學初始條件轉成可觀測暈統計的橋樑。

半解析模型與流體力學模擬。 在暗物質暈的併合樹（merger tree）骨架上，半解析模型（semi-analytic models, SAM）以一組描述冷卻、成星、回饋、化學增豐的微分方程「貼」上重子物理，計算成本低、可大量探索參數。另一條路是直接的宇宙學流體力學模擬（如 IllustrisTNG、EAGLE、SIMBA），在 $\Lambda$CDM 初始條件上同時演化暗物質與氣體，內建恆星與 AGN 子網格（sub-grid）回饋模型。兩者都成功重現星系質量函數與多數尺度關係，但子網格回饋的物理仍是最大不確定性——我們知道「需要」回饋，卻尚未從第一原理完全推導出其效率。

前沿張力與待解問題。 韋伯太空望遠鏡（JWST）在 $z>10$ 發現出乎意料明亮、似乎過早成熟的星系，挑戰標準成星效率的上限，可能暗示早期成星更高效、初始質量函數（IMF）偏向大質量恆星，或塵埃／光度估計需修正。此外，「小尺度危機」——衛星星系數量（missing satellites）、暈核心密度分布（core–cusp）、過大失敗（too-big-to-fail）——持續推動人們檢視冷暗物質的細節，並激發自相互作用暗物質（self-interacting dark matter）等替代模型。把這些張力與膨脹宇宙學參數（$H_0$、$\sigma_8$）的測量放在一起看，星系形成已從「定性故事」進化為精密檢驗宇宙學與基礎物理的實驗場。對有志深入者，掌握線性微擾論、輻射轉移與冷卻物理、以及數值模擬方法，是進入這個領域的三把鑰匙。