The hydrodynamical response of cold circumgalactic clouds to quasar radiation

Nicolas Ledos, Sebastiano Cantalupo, Titouan Lazeyras, Gabriele Pezzulli, Kentaro Nagamine, et al.

View Original ↗
AI 導讀 technology general 重要性 4/5

當冷氣體雲遭遇特定強度的類星體極紫外光,會觸發流體動力學的「火箭效應」,使萊曼阿爾法光度暴增達 10 倍。

  • 斯特龍根數(St)是決定氣體雲命運的核心指標,準確劃分光學薄、火箭效應與輻射屏蔽三種物理區間。
  • 在火箭效應區間,氣體雲內部產生極端雙峰密度分佈,中性氫聚束因子可飆升至 3000。
  • 高達 50-60% 的萊曼阿爾法光度來自氫複合機制,推翻過去單一靜態模型對黯淡類星體周圍氣體的假設。

天文學家在觀測環星系介質(CGM)時,往往難以釐清劇烈的類星體輻射如何改變冷氣體的物理型態。最新高解析度輻射流體動力學模擬指出,當溫度約 10^4 K 的冷氣體雲暴露在特定強度的類星體極紫外光下時,會觸發獨特的流體動力學響應,使其萊曼阿爾法(Lyα)輻射光度暴增高達 10 倍。這項研究不但量化了氣體游離狀態的演化條件,更表明在未達自屏蔽的氣體雲中,高達 50-60% 的發射光譜實際上來自氫複合機制。

環星系介質與類星體輻射的物理角力

探討星系演化的宇宙學放大模擬(zoom-in simulations)已經證實,星系周圍廣達數十萬秒差距的環星系介質中,存在著結構極度複雜的冷氣體網絡。過去十幾年來,透過 MUSE(多單元光譜探測儀)等大型觀測設備,研究人員能夠藉由類星體的強烈照射,直接觀測這些氣體在紅移大於 2 的宇宙深處所發出的萊曼阿爾法射線。然而,從觀測到的輻射反推氣體的真實密度與動力學狀態卻極度困難,因為發射機制受到氣體熱力學狀態的嚴重制約。

現有文獻針對恆星輻射如何蒸發星際雲已有諸多探討,但這些研究多半聚焦於密度極高(>10^3 cm^-3)且溫度極低(<100 K)的氣體,並不適用於類星體周圍的環星系介質。為填補這項空白,研究團隊建立了一套嚴謹的分析框架,並搭配三維輻射流體動力學模擬,專門針對暴露在極紫外光(EUV)下的 CGM 冷氣體雲進行推演。這套模型考慮了單向照射的游離光子通量如何與氫原子的重組機制達到動態平衡。

建構這套框架的核心在於定義無因次的斯特龍根數(Strömgren number,縮寫為 St)。St 被定義為氣體雲直徑與斯特龍根長度(游離光子被複合機制完全消耗的距離)的比值。研究團隊同時引入了一個與 St 正交的全新無因次參數 Υ,用來獨立追蹤游離輻射場的絕對強度。透過這兩個參數的組合,氣體雲在面臨輻射時的命運得以被精確劃分。

斯特龍根數決定的三大演化光譜區間

依照公式推導與模擬驗證,氣體雲受到輻射照射後會走向三種截然不同的演化途徑。第一種是光學薄區間(optically thin regime,條件為 St < 1)。在這種狀態下,強烈的輻射會在極短時間內完全穿透並游離整團氣體雲。氣體被均勻加熱後會向外膨脹,直到與周圍的熱星系介質達到壓力平衡。模擬顯示,這類氣體雲的最終溫度會穩定在約 7×10^4 K,密度則降至初始值的 6% 左右。

當氣體雲密度夠高或輻射相對微弱時,則會進入輻射屏蔽區間(radiation-shielded regime)。研究團隊計算出電離氣體與中性氣體的密度反差平方值(Δ^2,在給定初始溫度下約為 37)作為臨界閾值。當 St > 37 時,輻射只能穿透氣體雲表層極淺的區域,游離鋒面完全停滯。氣體雲內部幾乎不受影響,繼續維持約 9×10^3 K 的相對低溫,密度也保留在初始狀態的 80% 以上。

介於兩者之間的 1 < St < 37 則是動力學反應最劇烈的火箭效應區間(rocket-effect regime)。光子會在氣體雲迎光面建立一層光致游離區(PIR)。這個區域的高溫膨脹會對後方未游離的中性氣體產生強烈擠壓。當游離鋒面持續推進,被壓縮的殘餘冷氣體(稱為彗星狀球體)會因為前端膨脹氣體的反作用力而被向後加速,整個過程猶如火箭噴射推進器一般,徹底改變了氣體的結構與密度分佈。

50秒差距氣體雲的高解析度數值模擬

執行模擬時,團隊採用了 RAMSES-RT 輻射流體動力學程式碼,並設定氣體雲初始半徑為 50 pc(秒差距),這是宇宙學模擬中常見的緻密 CGM 團塊尺寸。為涵蓋完整的演化光譜,入射光子通量被設定為小於 10^9 cm^-2 s^-1,這大約相當於距離微弱類星體 10 千秒差距(pkpc)的輻射環境。周圍熱介質的溫度被設定為冷氣體雲的 400 到 1000 倍,以確保在 2000 萬年的模擬時間內背景氣體不會因為輻射冷卻而崩塌。

為求精準捕捉游離鋒面的動態,網格自適應加密(AMR)的最高解析度達到了氣體雲直徑的 1/80。輻射傳輸模型涵蓋了氫與氦的 6 個光子能量群組,並使用 M1 閉包近似法來處理輻射壓力。同時,為滿足降速光速近似法(reduced speed-of-light approximation)的物理條件,光速被設定為真實光速的 1%,這足以涵蓋光子穿越時間極短且游離鋒面推進速度不超過降速光速的物理情境。

分析氫原子數量密度的機率分佈函數(PDF)後,研究人員證實了火箭效應區間會產生極端的雙峰分佈。在這種狀態下,密度分佈會在 0.1 倍初始密度(對應游離區)與大於 1 倍初始密度(對應被壓縮的中性氣體)出現兩個極值。這種劇烈的密度波動直接反映在氫氣聚束效應因子(clumping factor)上,中性氫的聚束因子在火箭效應區間內可飆升至 2 到 3000 不等,而其他兩個演化區間的聚束因子則幾乎維持在 1 附近。

萊曼阿爾法輻射光度暴增的觀測意涵

從光度計算的結果來看,流體動力學響應對觀測數據的詮釋具有決定性影響。在光學薄區間,由於氣體完全游離且大幅膨脹,儘管發光體積增加,但密度平方依賴性導致複合輻射效率大跌,其標準化光度僅剩最大螢光光度的 4.6%。然而,一旦氣體落入火箭效應區間,強烈的內部壓縮與部分游離狀態會大幅拉高局部電子與氫離子密度。

數據顯示,在 St 接近 6 的中段火箭效應區間,冷氣體的萊曼阿爾法光度會比光學薄區間高出整整一到兩個數量級,峰值可達理想最大螢光光度的 70%。更關鍵的是,團隊將輻射拆分為碰撞激發與複合兩種機制後發現,在未達自屏蔽的氣體雲中,重組複合機制至少貢獻了 50-60% 的發射總量,打破了過去認為碰撞激發主導低密度氣體發光機制的刻板印象。

將這套理論套用到真實觀測環境中,團隊得出結論:在明亮類星體(特定光度達 10^31.6 erg/s/Hz)周圍的冷 CGM 氣體,幾乎無可避免會被完全游離並進入光學薄狀態;相反地,黯淡類星體(特定光度約 10^28.6 erg/s/Hz)的輻射強度恰好落在觸發火箭效應的甜蜜點。這意味著在分析不同亮度星體周圍的介質時,不能套用單一的靜態光度轉換公式。

觀測黯淡類星體周圍的環星系介質時,必須將輻射驅動的火箭效應納入考量,否則劇烈波動的氣體密度將導致嚴重的質量與運動狀態誤判。

Abstract

Recent simulations increasingly resolve the small-scale structure of the circumgalactic medium (CGM), but the dynamical impact of ionising radiation on its cold $10^4$ K component remains poorly understood. We investigate the evolution of cold gas structures exposed to quasars' EUV radiation. We develop an analytical framework to describe the evolution of such clouds, introducing a new threshold that defines when a cloud becomes radiation-shielded. The framework is validated using radiation-hydrodynamic simulations of single static clouds. It predicts three evolutionary paths: (i) an optically thin regime, in which radiation uniformly ionises the cloud; (ii) a radiation-shielded regime, where the cloud remains largely unaffected; and (iii) a rocket-effect regime, in which the propagation of the ionisation front ionises the illuminated side while compressing the opposite side, later accelerating the surviving cold clump. In the latter regime, the cloud's Ly$α$ luminosity can be up to one order of magnitude higher than the optically thin case. Such luminosities are as high as $70\%$ of the values obtained from a fluorescent regime without considering hydrodynamical response. Unless the cloud is shielded, at least $\sim 50$-$60\,\%$ of Ly$α$ emission arises from recombination. Applying this framework to both a ray crossing a population of clouds, and a ray propagating inside a cold stream, we find that the cold CGM around bright quasars ($L_{\mathrm{ν,LL}} \sim 10^{31.6} \, \mathrm{erg\, s^{-1}\, Hz^{-1}}$) is likely fully ionised, whereas the one around faint quasars ($L_{\mathrm{ν,LL}} \sim 10^{28.6} \, \mathrm{erg\, s^{-1}\, Hz^{-1}}$) predominantly experiences a rocket-effect regime. These results imply that the hydrodynamical response of cold CGM structures to quasar radiation must be considered when deriving their physical properties, particularly for faint quasars.