中子星與脈衝星（進階）：物態方程、超流體與宇宙級時鐘陣列

為什麼一茶匙中子星物質「重十億噸」這句話，其實只說對了一半？

入門篇裡你大概讀過那句經典比喻：一茶匙中子星物質重達十億噸。這句話沒錯，但它把中子星講成一團均勻的「超緻密泥巴」，而真相精彩得多。中子星的內部其實是一層層性質截然不同的物質：從外殼的晶格固體，到內部的中子超流體（neutron superfluid），再到我們至今仍無法在地球實驗室複製、可能由自由夸克（free quarks）組成的核心。

更關鍵的是，這顆直徑只有約 20 公里的星體，把「核物理」「廣義相對論」「量子流體力學」「電漿磁層物理」全部塞進同一個天體裡。當你問「一茶匙到底多重」時，答案其實取決於那一茶匙取自哪一層——而那正是當代物理學最前沿的未解問題之一：中子星的物態方程（equation of state, EOS）。

這篇進階篇不再重述燈塔模型，而是帶你進入殼層結構、自轉減速的能量帳、突發的「自轉跳變」（glitch），以及把整個星系當成偵測器的脈衝星計時陣列（pulsar timing array）。

物態方程：一場關於「半徑」的世紀賭局

中子星之所以能存在，是因為中子簡併壓（neutron degeneracy pressure）——源自包立不相容原理（Pauli exclusion principle）的量子壓力——抵抗了重力的塌縮。但簡併壓只是故事的開頭。在核心密度超過原子核密度（約 $2.8\times 10^{17}\ \mathrm{kg/m^3}$）數倍的環境下，強核力（strong nuclear force）的貢獻變得舉足輕重，而我們對這個密度區間的強交互作用所知仍然有限。

物態方程描述的是壓力 $P$ 與密度 $\rho$ 的關係 $P(\rho)$。不同的核物理模型給出不同的 $P(\rho)$，而這直接決定了一顆給定質量的中子星會有多大半徑：

• 「硬」（stiff）EOS：同樣密度下壓力較大，能撐起較大半徑與較高的最大質量。
• 「軟」（soft）EOS：壓力較小，星體被壓得更緊實，半徑較小。

這不是純理論遊戲。一個典型的 $1.4\ M_\odot$（太陽質量）中子星，不同 EOS 預測的半徑落在約 $10$ 到 $14$ 公里之間。聽起來差距不大，但這幾公里的差異，反映的是核物質在極端密度下行為的根本分歧。

關鍵的判別來自最大質量。如果某個 EOS「太軟」，它撐不起重的中子星就會塌成黑洞。2010 年起陸續發現質量約 $2\ M_\odot$ 的脈衝星（如 PSR J0740+6620，約 $2.08\ M_\odot$），一舉淘汰了一大批過軟的 EOS——任何不能撐起 $2\ M_\odot$ 的模型直接出局。

動手算一下：質量上限的物理直覺

中子星存在一個質量上限，稱為 Tolman–Oppenheimer–Volkoff（TOV）極限，類似白矮星的錢德拉塞卡極限（Chandrasekhar limit），但要用廣義相對論求解。其精確值依賴 EOS，落在約 $2.2$ 到 $2.9\ M_\odot$ 之間。

我們可以用一個簡單的尺度論證感受它的存在。把中子星粗略當成簡併中子氣，重力位能量級為 $\sim GM^2/R$，而相對論性簡併氣體的內能也按 $\sim 1/R$ 縮放。當兩者都隨 $1/R$ 變化時，半徑 $R$ 就從能量平衡中消失了——這正是出現「臨界質量」的數學徵兆。粗略地說：

$$ M_{\max} \sim \left(\frac{\hbar c}{G}\right)^{3/2}\frac{1}{m_n^2} $$

代入數值會得到接近一個太陽質量的量級——與真實的 TOV 極限同數量級。這個式子的美在於：它只用了普朗克常數 $\hbar$、光速 $c$、重力常數 $G$、中子質量 $m_n$ 這四個基本量，就「猜中」了一顆星的質量上限。緻密天體的物理就是這樣，把微觀的量子常數與宏觀的重力編織在一起。

洋蔥般的內部結構

中子星不是均勻的球，而是由密度遞增的數層構成：

1. 外殼（outer crust）：富中子的原子核排成庫侖晶格（Coulomb lattice），縫隙間是簡併電子海。這是一層「固體」，剛度極高。
2. 內殼（inner crust）：密度升高到中子開始從原子核「滴出」（neutron drip），自由中子形成中子超流體，與晶格共存。著名的「核義大利麵」（nuclear pasta）相——原子核被擠壓成棒狀、片狀等奇異幾何——就出現在這一層的底部。
3. 外核（outer core）：均勻的中子液體，混入少量質子與電子。中子超流、質子可能形成超導體（superconductor）。
4. 內核（inner core）：密度可達核密度數倍。這裡可能出現超子（hyperons）、$\pi$ 或 $K$ 介子凝聚，甚至解禁的夸克物質（deconfined quark matter）。沒有人確定，因為我們造不出這種物質。

超流體與超導體不是比喻，而是真實的量子相。中子配對成庫柏對（Cooper pairs），使流體黏滯性趨近於零並以量子化的渦旋（quantized vortices）攜帶角動量。這個微觀性質，待會兒會在解釋「glitch」時派上大用場。

自轉減速：脈衝星的能量帳本

入門篇說脈衝星是「精密的時鐘」，但這時鐘其實走得越來越慢。脈衝星不斷把自轉能量輻射出去，自轉週期 $P$ 緩慢增加，週期變化率 $\dot{P}>0$。

一顆毫秒脈衝星的 $\dot{P}$ 可能小到 $10^{-20}\ \mathrm{s/s}$——也就是說，要約十億年它的週期才會增加百分之幾。正是這種驚人的穩定性，讓它成為宇宙級的時鐘。

動手算一下：自轉能量的損耗率

中子星的轉動慣量約 $I \approx 10^{38}\ \mathrm{kg\,m^2}$。自轉動能為

$$ E_{\mathrm{rot}} = \frac{1}{2}I\Omega^2,\qquad \Omega = \frac{2\pi}{P} $$

能量損耗率（即「自轉減速光度」，spin-down luminosity）為

$$ \dot{E} = -\frac{dE_{\mathrm{rot}}}{dt} = -I\Omega\dot{\Omega} = 4\pi^2 I\,\frac{\dot{P}}{P^3} $$

我們代入著名的蟹狀星雲脈衝星（Crab pulsar）：$P \approx 0.0334\ \mathrm{s}$、$\dot{P}\approx 4.2\times 10^{-13}\ \mathrm{s/s}$。

$$ \dot{E} = 4\pi^2 \times 10^{38} \times \frac{4.2\times 10^{-13}}{(0.0334)^3} $$

分母 $(0.0334)^3 \approx 3.7\times 10^{-5}$，於是

$$ \dot{E} \approx 4\pi^2 \times 10^{38} \times 1.1\times 10^{-8} \approx 4.5\times 10^{31}\ \mathrm{W} $$

也就是約 $10^{38}\ \mathrm{erg/s}$，相當於太陽光度的十萬倍。這筆能量去了哪裡？它驅動了整個蟹狀星雲（Crab Nebula）的發光——這個事實在 1960 年代給了天文學家一個漂亮的閉環驗證：脈衝星損耗的自轉能，恰好等於星雲所需的能量供給。中子星就是那台藏在星雲中央、看不見的發電機。

磁偶極輻射與「特徵年齡」

脈衝星靠什麼把自轉能輻射出去？主流的磁偶極（magnetic dipole）模型指出：一顆磁矩 $m$、與自轉軸夾角 $\alpha$、以角速度 $\Omega$ 旋轉的磁化星體，會像旋轉的天線一樣輻射電磁波，功率為

$$ \dot{E}_{\mathrm{dipole}} = \frac{\mu_0\, m^2 \Omega^4 \sin^2\alpha}{6\pi c^3} $$

把它與自轉能損耗等量齊觀，可推得 $\dot\Omega \propto -\Omega^3$，更一般地寫成制動公式（braking law）：

$$ \dot{\Omega} = -K\Omega^{\,n} $$

其中 $n$ 稱為制動指數（braking index）。純磁偶極輻射預測 $n=3$。有趣的是，實測的幾顆脈衝星制動指數多落在 $2$ 到 $3$ 之間（蟹狀星雲約 $2.5$），略低於理論值——暗示真正的能量損耗不只是純偶極輻射，還牽涉磁層中的帶電粒子風（particle wind）與磁場結構演化。這種「理論與觀測的微小落差」，往往正是新物理的入口。

由制動公式還能估出脈衝星的特徵年齡（characteristic age）：

$$ \tau_c = \frac{P}{2\dot{P}} $$

對蟹狀星雲代入，$\tau_c = 0.0334 / (2\times 4.2\times 10^{-13}) \approx 1.3\times 10^{12}\ \mathrm{s} \approx 1260$ 年。而我們從史料知道蟹狀星雲是公元 1054 年超新星的遺骸——距今約 970 年。特徵年齡 1260 年與真實 970 年同數量級，誤差來自「假設出生時自轉遠快於現在」等近似。這是少數能用歷史記載直接校驗的天文時鐘之一。

Glitch：當時鐘突然「跳快」

長期來看脈衝星在減速，但偶爾它會突然自轉變快，週期瞬間縮短，這稱為自轉跳變（glitch）。蟹狀星雲與船帆座（Vela）脈衝星都是著名的「慣犯」。相對變化量

$$ \frac{\Delta\Omega}{\Omega} \sim 10^{-9}\ \text{到}\ 10^{-6} $$

雖然數值極小，但對毫秒級精度的計時來說是劇烈事件。glitch 的成因，正好用上前面提到的中子超流體與量子化渦旋。

機制大致如下：中子星的固體外殼隨著電磁輻射穩定減速，但內部的中子超流體並不直接跟著減速——它的角動量被困在量子化渦旋裡，而這些渦旋被「釘扎」（pinning）在內殼的晶格上動彈不得。於是超流體保持較快的自轉，與外殼產生越來越大的角速度差。當這個差距累積到臨界值，大量渦旋同時「解扎」（unpinning），把儲存的角動量瞬間轉移給外殼——外殼因此被「踢」得轉快了一點，我們就觀測到一次 glitch。

換句話說，glitch 是一扇難得的窗，讓我們間接觀測中子星內部的超流體。透過分析 glitch 後的恢復行為（角速度如何緩慢鬆弛回去），物理學家能反推超流體佔星體質量的比例、渦旋與晶格的耦合強度——這些都是地球實驗室永遠做不到的核物質量測。

看一個例子：磁星（magnetar）——磁場驅動的怪物

並非所有中子星都靠自轉發電。磁星（magnetar）是一類擁有極端磁場的中子星，表面磁場可達 $10^{10}$ 到 $10^{11}$ 特斯拉（$10^{14}$–$10^{15}$ 高斯）——比一般脈衝星強一千倍，是人類最強實驗室磁場的上億倍。

這個磁場有多誇張？我們可以估一下磁場的能量密度 $u_B = B^2/(2\mu_0)$。代入 $B = 10^{11}\ \mathrm{T}$：

$$ u_B = \frac{(10^{11})^2}{2\times(4\pi\times10^{-7})} \approx 4\times 10^{27}\ \mathrm{J/m^3} $$

這個磁能密度換算成等效質量密度 $u_B/c^2$ 約為 $4\times 10^{10}\ \mathrm{kg/m^3}$——磁場本身的「重量」就已超過地球上最緻密物質的密度。在這種磁場下，原子被拉成細長的「鉛筆狀」，真空本身都會雙折射（vacuum birefringence），量子電動力學（QED）的非線性效應變得可觀。

磁星的能量來源不是自轉，而是磁場的衰變與重組。當磁應力突破外殼的剛度，會發生「星震」（starquake），釋放出短暫但極亮的 $\gamma$ 射線爆發，稱為軟 $\gamma$ 重複爆（soft gamma repeater, SGR）。2004 年一次來自 SGR 1806–20 的巨閃焰（giant flare），在零點幾秒內釋放的能量超過太陽十萬年的總輻射，即使源在約 5 萬光年外，仍擾動了地球的電離層。近年磁星也成為部分快速電波爆（fast radio burst, FRB）的有力候選來源。

重點回顧

• 物態方程（EOS）是核心未解問題：中子星半徑與最大質量直接由 $P(\rho)$ 決定；發現 $2\ M_\odot$ 脈衝星淘汰了所有「過軟」的模型。
• 內部是分層的量子流體：外殼是晶格固體，內部是中子超流體與質子超導體，核心可能藏著夸克物質。
• 脈衝星是減速的時鐘：自轉能透過磁偶極輻射損耗，$\dot{E}=4\pi^2 I\,\dot{P}/P^3$；蟹狀星雲的損耗恰好點亮整個星雲。
• 特徵年齡 $\tau_c = P/2\dot{P}$ 能用歷史超新星記載校驗，制動指數 $n$ 偏離 $3$ 透露純偶極模型的不足。
• Glitch 是內部超流體的指紋：量子化渦旋的釘扎與解扎讓我們間接「看見」中子星內部，而磁星則展示了磁場主導的另一種極端。

深入探討（研究所視角）

從 EOS 到觀測量的橋樑：多信使約束。 當代研究的主軸是用多種獨立觀測夾擊 EOS。NASA 的 NICER 任務透過 X 射線脈衝輪廓建模（pulse-profile modeling），利用熱斑（hot spot）發出的光子在強重力下的彎曲（光線追跡需解 Schwarzschild 或 Hartle–Thorne 度規），同時反推中子星的質量與半徑。另一條路來自重力波：2017 年的雙中子星合併事件 GW170817，其波形中的潮汐形變性（tidal deformability） $\Lambda$ 對 EOS 高度敏感——較軟的 EOS 使星體較難被潮汐拉扯，$\Lambda$ 較小。LIGO/Virgo 的測量給出 $\Lambda$ 上限，與 NICER 的 $M$–$R$ 量測相互一致，把容許的 EOS 區間收窄到一條相對狹窄的帶狀區。

夸克物質與相變的可能簽名。 內核是否發生從強子物質（hadronic matter）到解禁夸克物質的相變，是理論的聖杯。一階相變會在 $M$–$R$ 曲線上留下特徵性的「扭折」，甚至允許雙星支（twin stars）——兩顆質量相同但半徑不同的穩定構型。若色超導（color superconductivity）的 CFL（color–flavor–locked）相在核心實現，將對冷卻曲線與震盪模式留下印記。

漸近與震盪：g、p、r 模與連續引力波。 中子星的非徑向振盪模式（quasi-normal modes）會輻射重力波，其中 r 模（Rossby mode） 受 Chandrasekhar–Friedman–Schutz（CFS）不穩定性驅動，可能限制快速自轉中子星的最高轉速。若中子星帶有微小的質量四極矩（如「山脈」，高度或許僅幾公分），旋轉時會放出連續重力波（continuous gravitational waves）——目前尚未偵測到，其上限本身已對中子星外殼的剛度設下約束。

脈衝星計時陣列：把星系變成偵測器。 最前沿的應用之一是 脈衝星計時陣列（pulsar timing array, PTA）。把分布在天空各方向的數十顆毫秒脈衝星當成一組超穩定時鐘，長期監測它們脈衝到達時間（time of arrival）的微小相關偏差。當奈赫茲（nanohertz）重力波背景——主要源自宇宙中無數對超大質量黑洞雙星——穿過地球與脈衝星之間的時空時，會在不同脈衝星的計時殘差間留下一種特定的角度相關（Hellings–Downs 曲線）。2023 年，NANOGrav、EPTA、PPTA 等團隊同時報告了與此一致的證據。換句話說，我們正用一顆顆中子星，傾聽整個宇宙最低沉的時空嗡鳴。從一茶匙物質的重量，到把銀河系變成重力波天線，中子星不斷把「極端物理」交到我們手中。