Constraining the Pulsar Beaming Fraction with TeV-Selected Galactic Pulsar Wind Nebulae and unidentified TeV Sources

銀河系中目前觀測到的 4000 顆無線電脈衝星僅佔實際總量的極小部分，精算潛藏數目的關鍵在於確立波束比例。東京大學團隊利用 TeV 頻段的脈衝星風雲數據，計算出該比例落在 0.1 到 0.3 之間，這項觀測大幅修正了過往高度依賴理論模型的誤差，為恆星演化與重力波預測確立了全新的量化基準。

重估 4000 顆樣本：波束比例與重力波探測關聯

脈衝星（Pulsars）的波束比例（Beaming Fraction）定義為射電束在自轉一圈時掃過的立體角佔全天球的比例，它直接決定了隨機位置的地球觀測者能看到該星體的機率。科學界過去多半依賴觀測到的無線電與伽馬射線訊號，套用極冠模型或外隙模型來逆向推導整體族群的面貌。然而，這種做法產出的數值高度依賴理論設定，無法精準代表全銀河系脈衝星的真實幾何結構。

這項參數的誤差甚至直接影響了重力波事件的預期偵測率。過去藉由傳統波束比例推演的群體合成模型，預測雷射干涉引力波天文台（LIGO）每年能觀測到 4.6 次 雙中子星併合（DNS merger）事件。真實情況卻是在 O4 運行期間經過三年的持續運作，依然沒有偵測到任何確切的雙中子星併合重力波。

這種理論與實測的巨大落差，暗示了系統性不確定性遠比想像中高。科學家意識到，過往基於射電訊號所假設的波束比例可能存在高估，因此亟需一種不受到狹窄電磁波束方向限制的新觀測手段，來重新定義銀河系內的脈衝星總數。

跳脫模型限制：觀測各向同性的 TeV 脈衝星風雲

為了避開狹窄波束造成的觀測盲區，研究團隊將目光轉向脈衝星風雲（PWNe）。這類天體是由脈衝星自轉減慢所釋放的相對論性粒子驅動，並與超新星殘骸共同演化的巨大發光氣泡。在 TeV（兆電子伏特）能量範圍內，PWNe 的輻射機制主要來自相對論性電子在星雲內部進行的逆康普頓散射。

相較於脈衝星本體高度集中的射電束，PWNe 所在的發光區域光學厚度薄，且具備非相對論性的整體速度。這意味著其 TeV 頻段的輻射幾乎是各向同性（isotropic）的。不論脈衝星的自轉軸與磁場方向如何傾斜，只要該系統存在，地球上的高能天文台就有機會捕捉到這團均勻發光的 TeV 氣泡。

利用這項特性，PWNe 成為絕佳的間接追蹤器。團隊不再需要假設特定的發光幾何形狀，只要清點天空中所有明亮的 PWNe 數量，再對比已經確認內部存在脈衝星的 PWNe 數量，就能從宏觀統計角度直接反演出波束比例的客觀數值。

建構對比算式：將未辨識 TeV 來源視為隱藏脈衝星

研究設計了一套嚴謹的分類基礎。團隊採用了 H.E.S.S.、HAWC 與 LHAASO 三大高能伽馬射線巡天計畫的資料庫，將觀測到的 TeV 來源分為兩大類。第一類是已明確對應到脈衝星的 PWNe，代表這顆恆星的波束恰巧有掃過地球；第二類則是尚未確認發光機制的未辨識來源（Unid）。

算式的核心邏輯在於將所有位於銀河平面的 Unid 來源，一律視為「波束方向並未指向地球的脈衝星」所驅動的 PWNe。計算公式非常直觀，就是將已辨識 PWNe 的數量除以 PWNe 與 Unid 來源的總和。這樣的假設確保了運算得出的是波束比例的統計下限，因為未來若部分 Unid 來源被證實為活躍星系核（AGN）或其他天體，作為分母的數量將會減少，從而推升真實的波束比例。

在計算累積數據時，團隊針對特定物理參數設置了觀測閾值。例如自轉減慢光度（$\dot{E}$）必須大於給定極限，或是特徵年齡（$\tau_c$）低於某個標準。由於 Unid 來源缺乏詳細參數，運算時會一律假設它們皆滿足設定的極限條件，藉此還原出各個能階群體在不同波段上的表現。

H.E.S.S. 巡天計畫數據：波束比例落在 0.1 到 0.3

統計結果揭示了跨波段的一致性。在 H.E.S.S. 儀器的數據庫中，總共收錄了 25 個 PWNe 以及 35 個 Unid 來源。經過公式換算，無論是針對無線電、伽馬射線還是 X 射線波段的脈衝星，其得出的波束比例皆平穩落在 0.1 到 0.3 區間。

這項數值在自轉減慢光度大於 $10^{36} \text{ erg/s}$ 的高能量區段維持著極高的穩定性。同時，若單獨抽取出具備 X 射線脈衝訊號的樣本，與僅具備非熱輻射但無脈衝的 X 射線樣本對比，前者的推估比例甚至會縮減一半。這顯示有相當高比例的系統即使釋放了強烈的 X 射線非熱輻射，其脈衝輪廓依然無法被現有儀器解析。

然而，當研究導入 HAWC 與 LHAASO 的資料時，數據卻出現了巨大的分歧。HAWC 陣列記錄了 8 個 PWNe 與 44 個 Unid 來源，而 LHAASO 則囊括了 8 個 PWNe 與高達 60 個 Unid 來源。基於這兩套資料庫所算出的波束比例，僅有 H.E.S.S. 數據的三分之一到二分之一，這種顯著差異證明了數據背後存在強烈的儀器選擇偏誤。

儀器角分辨率與 1-10 TeV 能量區間的選擇效應

導致兩組數據落差超過兩倍的原因，在於觀測儀器的硬體規格與能量敏感度。H.E.S.S. 擁有極佳的角分辨率，在 250 GeV 以上可達 $0.048^\circ$，但其窄視場設計導致它難以捕捉擴散範圍過大的天體。它對 1-10 TeV 區段最為敏感，這恰好是年輕 PWNe 能量釋放的集中帶。

相對地，HAWC 與 LHAASO 屬於廣角空氣陣列觀測站，能夠無間斷掃描大片天空。這兩台儀器的優勢在於偵測大於 10 TeV 的高能量輻射，並且具備絕佳的擴展源捕捉力。隨著脈衝星風雲老化，同步加速器自康普頓效應減弱，其頻譜峰值會往更高能階移動；同時，氣泡體積膨脹會導致表面亮度下滑，使角直徑甚至超過 $1^\circ$，這類中年與老年演化系統幾乎只有 HAWC/LHAASO 能夠大範圍記錄。

這種先天設計使得 H.E.S.S. 的樣本庫充滿了年輕、能量密集且體積緊湊的脈衝星。反觀 HAWC 與 LHAASO 的觀測庫中，充斥著年紀老邁、自轉減慢光度較低且佔據極大視角的演化末期天體。兩者對目標天體的篩選偏好，直接造就了統計分母的巨大落差。

蒙地卡羅模型驗證：隨時間變化的波束張角假設

為驗證這兩套觀測數據的矛盾是否能被物理機制統一，團隊導入了蒙地卡羅（MC）族群合成模型。該模擬不直接套用特定的幾何發光模型，而是設定初始自轉週期與初始偶極磁場，透過力無場（force-free）方程式運算脈衝星隨時間推移的磁傾角演化與能量衰退。

模擬的核心假設在於將波束張角（opening angle）設定為隨時間變化的動態參數。這意味著脈衝星年輕時與邁入中老年時，其射電束覆蓋的立體角比例並不相同。在將 PWNe 演化軌跡與儀器觀測條件輸入系統後，模型在維持現有脈衝星族群統計特徵的前提下，順利重現了 H.E.S.S. 與 HAWC/LHAASO 之間超過兩倍的比例差距。

這套驗證不僅合理解釋了儀器選擇效應對大數據統計的影響，也證明了各向同性的 TeV 來源能夠作為評估星體幾何變化的絕佳指標，為雙中子星併合機率的校準提供了堅實的後盾。

波束比例的推算不應受限於單一幾何模型，整合儀器敏感度與脈衝星演化年齡，才能還原銀河系真實的天體基數。