太初核合成進階：把宇宙最初三分鐘變成一台精密儀器

如果中子壽命再長 10 秒，宇宙會是什麼模樣？

你已經知道太初核合成(Big Bang Nucleosynthesis, BBN)在大霹靂後約三分鐘鑄造了氫與氦，也知道那個漂亮的 $Y_p \approx 0.25$。但這裡有一個值得認真追問的問題：這整套計算到底有多「精密」？我們常聽說 BBN 是「精密宇宙學(precision cosmology)」的支柱，可是一個發生在 138 億年前、持續不到二十分鐘、溫度高達十億度的事件，憑什麼能算到小數點後第三位、還能拿來檢驗粒子物理？

答案藏在三件入門篇沒細談的事情裡：核反應網路的動力學細節、宇宙膨脹與溫度演化的精確耦合，以及那些「輸入參數」的不確定性如何傳播到最終豐度。這篇進階文章不再重述瓶頸與凍結的故事，而是把 BBN 當成一台精密儀器拆開來看——看它的齒輪如何咬合，看它對哪些參數敏感、對哪些麻木，以及為什麼一個叫「中子壽命」的實驗室數字，會直接寫進宇宙的氦含量。

把 BBN 寫成一組微分方程

入門篇把核合成描述成「氘瓶頸一突破，反應就瞬間連鎖完成」。這個圖像對，但太粗。真實的 BBN 是一組剛性常微分方程(stiff ODE system)：每一種核種 $i$ 的數量豐度 $Y_i \equiv n_i / n_b$（相對於重子總數）隨時間的變化，由所有產生與消滅它的反應決定。

對一個典型的雙體反應 $i + j \to k + l$，其貢獻形如

$$ \frac{dY_k}{dt} = n_b\,\langle \sigma v \rangle_{ij}\, Y_i Y_j - (\text{逆反應與其他通道}) $$

其中 $\langle \sigma v \rangle$ 是熱平均反應率(thermally averaged reaction rate)——把核反應截面 $\sigma$ 對馬克士威-波茲曼速度分布做平均的結果。整個 BBN 標準計算（始於 1967 年 Wagoner 的程式，今天的後裔是 PArthENoPE、AlterBBN、PRIMAT 等公開程式碼）追蹤大約十幾種核種、上百條反應通道，把它們和宇宙膨脹方程一起數值積分。

這裡的關鍵洞見是：BBN 的不確定性，主要不是來自宇宙學，而是來自核物理實驗室。$\langle \sigma v \rangle$ 是從加速器上量到的核反應截面外推來的。有些反應（例如 $p(n,\gamma)d$，質子捕獲中子生成氘）在 BBN 相關能量下測得很準，有些（例如生成鋰-7 的 $^3\text{He}(\alpha,\gamma)^7\text{Be}$）則仍有可觀誤差。現代 BBN 的「精密」，本質上是核反應率資料庫的精密。

溫度與時間：宇宙的計時器如何運轉

要解上面那組方程，你必須知道每個時刻的溫度 $T$、重子數密度 $n_b$ 與膨脹率 $H$。這三者不是獨立的，而是被弗里德曼方程(Friedmann equation)鎖在一起。

在輻射主導(radiation-dominated)的早期宇宙，能量密度由相對論性粒子主宰：

$$ \rho \;=\; \frac{\pi^2}{30}\, g_*(T)\, \frac{(k_B T)^4}{(\hbar c)^3} $$

代入弗里德曼方程 $H^2 = \frac{8\pi G}{3}\rho$，可以推出一個極為實用的近似式——溫度與宇宙年齡的對應關係：

$$ t \;\approx\; \frac{2.4}{\sqrt{g_*}}\;\left(\frac{1\ \text{MeV}}{k_B T}\right)^2 \ \text{秒} $$

這條式子是 BBN 的「碼錶」。讓我們用它感受一下尺度感：當 $k_B T = 1\ \text{MeV}$（約 $1.16 \times 10^{10}$ K）時，宇宙才約 1 秒大；當溫度降到 $k_B T \approx 0.07\ \text{MeV}$（氘瓶頸突破點）時，

$$ t \approx \frac{2.4}{\sqrt{g_*}}\left(\frac{1}{0.07}\right)^2 \approx \frac{2.4 \times 204}{\sqrt{g_*}} \ \text{秒} $$

取核合成期的 $g_* \approx 3.4$，得 $t \approx 265$ 秒，正是那個著名的「約三分鐘」。注意整個 BBN 視窗 $T$ 從 $\sim 0.1$ 掉到 $\sim 0.03\ \text{MeV}$，對應時間從約 3 分鐘到約 20 分鐘——溫度只變了三倍，時間卻拉長了一個數量級，因為 $t \propto T^{-2}$。

膨脹同時稀釋了重子密度。原初核反應之所以「來不及」越過質量數 8 的鴻溝，根本原因就是 $n_b \propto a^{-3} \propto T^3$ 隨膨脹急速下降，三體碰撞率正比於 $n_b^2$，在這種密度下完全可以忽略。宇宙不是「不想」造碳，而是膨脹得太快、稀釋得太徹底。

中子壽命：一個實驗室數字如何寫進宇宙

現在來回答開篇的問題。氦豐度 $Y_p$ 對好幾個輸入量敏感，其中最直接的兩個是：重子-光子比 $\eta$（入門篇談過）與中子平均壽命 $\tau_n$。後者特別有趣，因為它是一個能在地球實驗室裡量、卻會直接改變宇宙化學組成的數字。

回憶一下邏輯鏈：中子質子比在凍結後並非定值，而是隨自由中子的 $\beta$ 衰變持續下降。從凍結（$t \approx 1$ 秒）到核合成啟動（$t \approx 3$ 分鐘），中子有大約 $\Delta t \approx 200$ 秒在衰變。存活比例是

$$ \frac{n_n(t_{\text{nuc}})}{n_n(t_{\text{freeze}})} \;=\; \exp\!\left(-\frac{\Delta t}{\tau_n}\right) $$

中子壽命越長，越多中子撐到核合成、越多被鎖進氦-4，$Y_p$ 越高。一條被廣泛引用的近似關係是

$$ \frac{\Delta Y_p}{Y_p} \;\approx\; 0.72\,\frac{\Delta \tau_n}{\tau_n} $$

也就是中子壽命變動 1%，氦豐度跟著變動約 0.72%。

動手算一下：中子壽命之爭對 $Y_p$ 的影響

這不是純學術問題。中子壽命的兩種主流測量法——「瓶法(bottle method)」（把超冷中子關在容器裡數剩幾個）與「束法(beam method)」（讓中子束飛過、數產生了幾個質子）——長期給出不一致的結果：瓶法約 $878.4\ \text{s}$，束法約 $887.7\ \text{s}$，差距約 $9\ \text{s}$，遠超各自的誤差棒。這個「中子壽命謎題(neutron lifetime puzzle)」至今懸而未決。

我們來估算這 $9\ \text{s}$ 的差異對宇宙氦含量的影響。以 $\tau_n \approx 880\ \text{s}$ 為基準：

$$ \frac{\Delta \tau_n}{\tau_n} = \frac{9}{880} \approx 1.02\% $$

$$ \frac{\Delta Y_p}{Y_p} \approx 0.72 \times 1.02\% \approx 0.74\% $$

$$ \Delta Y_p \approx 0.74\% \times 0.247 \approx 0.0018 $$

換句話說，兩種中子壽命測量法的差異，會讓預測的原初氦質量豐度相差約 $0.0018$，也就是 $Y_p$ 在 $0.246$ 與 $0.248$ 之間擺動。這個量級已經逼近天文觀測的精度（最佳的氦豐度測量誤差約 $\pm 0.003$）。一場關於中子在罐子裡能活多久的桌上實驗，最終會影響我們對整個宇宙氦含量的理論預測——這正是 BBN 把粒子物理與宇宙學焊在一起的具體例證。

為什麼氘是「敏感的」、氦是「飽和的」

入門篇提過氘是好用的密度計，氦則穩定地落在 25%。現在我們用反應動力學的語言，把這個差異講透，因為它揭示了 BBN 為何能同時又精密又穩健。

關鍵在於每種核種對 $\eta$ 的反應「徹底程度」不同：

• 氦-4 是「燃燒終點」：它束縛能最高，是反應網路的吸引子(attractor)。只要核合成有發生，幾乎所有中子最終都進到氦-4，所以 $Y_p$ 幾乎只取決於「有多少中子可用」（由 $\tau_n$ 與凍結決定），而對 $\eta$ 只有微弱的對數依賴。$Y_p$ 因此飽和、穩健，是檢驗「中子數與膨脹率」的探針。
• 氘是「殘留的中間產物」：它是通往氦的中繼站，本質上是「沒被燒完的灰燼」。$\eta$ 越大、密度越高，氘被進一步燒成氦越徹底，殘留越少。氘豐度因此對 $\eta$ 呈陡峭的冪次依賴，近似

$$ \frac{\text{D}}{\text{H}} \;\propto\; \eta^{-1.6} $$

這條冪律就是氘成為靈敏密度計的數學根據：$\eta$ 變 10%，D/H 變約 16%。

• 鋰-7 最詭異：它有兩條生成路徑，在低 $\eta$ 區主要由 $^3\text{H}$ 直接合成，在高 $\eta$ 區則主要經由 $^7\text{Be}$ 後來電子捕獲衰變而成。兩條路徑對 $\eta$ 的依賴相反，疊加起來形成一條非單調的「U 形谷」曲線。真實宇宙的 $\eta$ 不幸正好落在能讓兩條路徑都貢獻的區間，使鋰-7 的預測既敏感又困難——這也是「鋰問題」格外棘手的技術原因之一。

理解這層差異，你就懂了 BBN 觀測策略的分工：用氘量 $\eta$（密度），用氦量 $N_{\text{eff}}$（膨脹率與粒子種類），用鋰找麻煩（新物理線索）。三種輕元素各司其職，互不替代。

重點回顧

• BBN 在運算上是一組剛性微分方程，追蹤十幾種核種與上百條反應通道；它的「精密」主要受限於實驗室測得的核反應率 $\langle \sigma v \rangle$，而非宇宙學參數。
• 溫度與時間由弗里德曼方程鎖死，$t \propto g_*^{-1/2}\,T^{-2}$ 是 BBN 的碼錶；整個合成視窗溫度僅變約三倍，時間卻拉長近一個數量級。
• 中子壽命 $\tau_n$ 直接決定有多少中子撐到核合成，$\Delta Y_p/Y_p \approx 0.72\,\Delta\tau_n/\tau_n$；瓶法與束法約 9 秒的爭議，會讓 $Y_p$ 預測擺動約 $0.0018$。
• 不同核種對 $\eta$ 的依賴截然不同：氦-4 飽和而穩健、氘呈 $\eta^{-1.6}$ 陡峭冪律、鋰-7 是非單調 U 形谷；這決定了各自的觀測用途。
• 三種輕元素分工互補——氘測重子密度、氦測膨脹率與粒子種類、鋰提供超越標準模型的線索。

深入探討（研究所視角）

要把 BBN 推到當代研究前沿，需要把上面的圖像放進更嚴謹的動力學與宇宙學框架。以下三條線索，是研究所與當前文獻真正在處理的問題。

凍結的本質：玻茲曼方程與離開平衡

入門篇用 $\Gamma_{\text{weak}} \sim H$ 的尺度比較解釋凍結，這是啟發式的。嚴格處理要解中子質子比 $X_n \equiv n_n/(n_n+n_p)$ 的玻茲曼方程(Boltzmann equation)：

$$ \frac{dX_n}{dt} \;=\; -\lambda_{n\to p}\,X_n \;+\; \lambda_{p\to n}\,(1 - X_n) $$

其中 $\lambda_{n\to p}$、$\lambda_{p\to n}$ 是把弱交互作用各通道（$n+\nu_e \to p+e^-$、$n+e^+ \to p+\bar\nu_e$、中子衰變等）的相空間積分做完後的總轉換率，正比於 $G_F^2$ 並對溫度有複雜依賴。在高溫時兩率龐大、$X_n$ 緊貼平衡解 $X_n^{\text{eq}} = [1+e^{\Delta m c^2/k_BT}]^{-1}$；當率掉到 $H$ 以下，$X_n$ 偏離平衡並「凍結」。值得強調的是凍結不是瞬間的——它是一段漸進偏離平衡的過程，必須數值求解才能得到精確的凍結豐度。現代計算還要加上一階弱作用輻射修正、有限溫度 QED 修正、以及微中子不完全退耦(incomplete neutrino decoupling)的效應，這些一起把 $Y_p$ 的理論誤差壓到 $\sim 0.0001$ 量級。

$N_{\text{eff}}$ 作為暗輻射的探針

入門篇提到 $Y_p$ 對有效相對論性自由度 $g_*$ 敏感，能限制微中子種類。研究層級把這量化為有效微中子數 $N_{\text{eff}}$，它衡量早期宇宙中除光子外所有相對論性「輻射」的總量：

$$ \rho_{\text{rad}} = \rho_\gamma\left[1 + \frac{7}{8}\left(\frac{4}{11}\right)^{4/3} N_{\text{eff}}\right] $$

標準模型因微中子退耦與電子-正子湮滅的細節，精確預測 $N_{\text{eff}}^{\text{SM}} = 3.044$（不是整數 3）。任何額外的輕自由度——惰性微中子(sterile neutrino)、暗光子、極輕的新粒子——都會抬高 $N_{\text{eff}}$，增大 $H$，使中子凍結提前、$Y_p$ 上升。BBN 對 $Y_p$ 與 D/H 的聯合分析給出 $N_{\text{eff}} = 2.9 \pm 0.3$ 量級的約束，與 CMB 的獨立測量一致。這讓 BBN 成為時間上最早的暗輻射偵測器：它在大霹靂後一秒就「秤」過宇宙的輻射成分，對任何在 MeV 能標下仍為相對論性的新粒子設下界限。值得注意的是，BBN 與 CMB 探測的是不同紅移（$z \sim 10^9$ 對 $z \sim 1100$），兩者一致本身就是對「$N_{\text{eff}}$ 不隨時間演化」的檢驗，能約束如衰變暗物質等情境。

鋰問題的可能出路與新物理視窗

原初鋰問題（理論預測高出觀測約 2–3 倍）至今是標準 BBN 唯一的明顯張力，當前研究沿三條路徑進攻。其一是核物理：重新檢視 $^7\text{Be}$ 的生成與摧毀反應率，特別是共振態(resonance)貢獻是否被低估；但近年實驗大致排除了用核反應率單獨解決問題的可能。其二是天文物理：「Spite 平台」恆星可能透過大氣對流、原子擴散或自轉混合，把表面鋰沉降耗損掉，使觀測值低於真正的原初值——這條路徑在貧金屬恆星模型中漸獲支持。其三是新物理：若早期宇宙存在衰變粒子或時變基本常數，可在不破壞氘與氦一致性的前提下選擇性地壓低鋰-7，這也是約束超對稱長壽命粒子等模型的窗口。BBN 的精密恰恰使這 2–3 倍的偏差無法被輕易掩蓋——一個成功到能秤量宇宙、能數粒子種類的理論，仍在它最輕的一個元素上，為下一代研究者留著一道未解的縫。從一組微分方程、一個實驗室裡的中子壽命、到一條約束暗輻射的界限，太初核合成展示了精密宇宙學的真諦：把最古老的事件，變成最嚴格的實驗。