本地极亮红外星系中的高中红外光变发生率

《Publications of the Astronomical Society of Japan》：A high incidence of mid-infrared variability in local ultraluminous infrared galaxies

　　

研究人员利用广域红外巡天探测器（Wide-Field Infrared Survey Explorer，WISE）的NEOWISE多历元星表约11年的测光数据，探索了本地极亮红外星系（ultraluminous infrared galaxies，ULIRGs；红外光度L_IR>10¹²L_⊙）中的中红外（mid-infrared，MIR）光变。研究人员识别出30个具有统计显著MIR光变的ULIRGs。该光变在数年时标上被观测到，意味着可变MIR辐射区域是致密的（?几秒差距）。30个可变ULIRGs的最大与最小W2（4.6 μm）波段光度之差（ΔL_W2）范围为ΔL_W2 = 7×10⁴²到5×10⁴⁴ erg s^?1。30个可变ULIRGs中有25个的ΔL_W2大于ΔL_W2 = 10⁴³ erg s^?1，超过了已知超新星（supernovae，SNe；L_3.6μm和L_4.5μm < 10^42.3 erg s^?1）所观测到的MIR光度范围。因此，这25个ULIRGs中的MIR光变极有可能由潮汐瓦解事件（tidal disruption events，TDEs）或活动星系核（active galactic nuclei，AGN）活动驱动。该样本包括硬X射线探测到的AGN（如UGC 05101）和先前报道的TDE候选体（IRAS F01004?2237、IRAS 05189?2524）。这25个源除了MIR光变外，均表现出至少一种AGN特征，表明即使MIR光变起源于TDEs，所涉及的黑洞很可能也是AGN。研究结果表明，MIR光变是探测掩埋型AGN的有效工具，并凸显了ULIRGs中强烈的核活动。

**论文解读：利用中红外光变探测本地极亮红外星系中的掩埋型活动星系核**

**1. 研究背景与科学动机**

几乎所有大质量星系中心都存在一个超大质量黑洞（supermassive black hole，SMBH），且SMBH质量与核球恒星质量存在相关性。极亮红外星系（ultraluminous infrared galaxies，ULIRGs）被认为是富含气体的星系合并产物，是核球和SMBH增长的重要场所。然而，在ULIRGs中，增长中的SMBH往往被尘埃深度包裹，这些尘埃阻挡了电离光子到达窄线区，使得诸如Baldwin–Phillips–Terlevich（BPT）图等光学诊断方法无法探测到这些掩埋型活动星系核（active galactic nuclei，AGN）。因此，需要在受尘埃消光影响较小的波段进行观测。

中红外（mid-infrared，MIR）观测受消光影响远小于光学波段，为研究ULIRGs中的AGN提供了有效工具。Stern等人（2012）的研究表明，MIR颜色W1?W2（[3.4]?[4.6]）是区分AGN与正常星系的有效指标，因为AGN加热的热尘埃辐射比星系恒星辐射表现出系统性更红的MIR颜色。然而，在某些情况下，恒星形成区的尘埃辐射也可能模拟出类似的红色，这一局限性促使研究人员将MIR光变作为额外诊断手段。光变能够清晰区分AGN活动与恒星形成：AGN在全电磁波谱上均表现出光变，其红外波段的可变辐射源于亚秒差距到秒差距尺度上的热尘埃对中心吸积盘紫外/光学辐射的再处理；而恒星形成通常延伸至千秒差距尺度，除超新星（supernovae，SNe）外，预期不会在可观测时标上产生光变。即使在MIR探测到的最亮SNe，其Spitzer 4.5μm波段光度也仅达到10^42.3 erg s^?1，因此应用光度阈值可以有效消除亮SNe的污染。此外，虽然邻近ULIRGs中观测到非常致密的（约100 pc）核恒星暴，但其光穿越时标仍超过几×10²年，比WISE探测的时标大一个数量级。

AGN的光变通常归因于吸积盘的变化，可能源于吸积盘不稳定性、非均匀温度涨落，或吸积率的极端增强（常见于变脸AGN，changing-look AGNs，CLAGNs）。吸积率增加可能来源于潮汐瓦解事件（tidal disruption events，TDEs）——当恒星进入SMBH的潮汐半径内被瓦解并随后吸积碎片时触发。在ULIRGs中已发现TDE候选体，并在MIR波段探测到红外回响。

此前Reynolds等人（2022）对ULIRGs进行了部分系统性MIR光变分析，但主要聚焦于在亮红外星系（luminous infrared galaxies，LIRGs）和ULIRGs中识别TDE候选体，而非探测源自AGN的光变，且其ULIRG样本局限于Sanders等人（2003）的星系列表。相比之下，Nardini等人（2010）后期发表的星表包含了更多ULIRGs。本研究将原本为矮星系开发的统计光变探测技术应用于Nardini等人（2010）更大、更完整的ULIRG星表，系统性地搜索MIR光变，以探测掩埋型AGN活动。

**2. 样本与关键技术方法**

- **样本来源**：选取Nardini等人（2010）中的全部164个ULIRGs，该样本基于红外天文卫星（Infrared Astronomical Satellite，IRAS）PSCz巡天和60 μm 1 Jy ULIRG样本，覆盖完整本地ULIRG光度范围，对AGN活动无强烈偏置。所有源均具有Spitzer光谱观测数据。
- **测光数据与光变曲线构建**：利用WISE/NEOWISE任务约11年（2013–2024）的W1（3.4 μm）和W2（4.6 μm）单次曝光测光数据。通过AllWISE星表进行位置交叉匹配，提取NEOWISE-R单次曝光（L1b）源表数据，并按约半年一个历元进行加权平均分箱，同时采用三倍中位数绝对偏差进行异常值剔除。
- **光变判据**：应用W2<14.5 mag的星等截断以消除爱丁顿偏差，计算每颗ULIRG的W1与W2波段无偏Pearson相关系数r，以及两波段的方差σ_W1和σ_W2。光变判定标准为r>0.75、σ_W1>0.010且σ_W2>0.034，后者对应NEOWISE系统误差阈值。
- **光度阈值比较**：计算可变源的最大与最小W2波段光度差ΔL_W2，并与已知最亮MIR超新星的光度（10^42.3 erg s^?1）进行比较，以排除SNe污染。

**3. 研究结果**

**3.1 光变检测与光学分类相关性**

对155个具有足够数据的ULIRGs计算统计量后，r值分布表现出明显的光学类型依赖性：H ii星系和LINERs大多位于较低r值（r?0.65），而赛弗特星系趋向于较高r值。赛弗特2型星系构成r>0.9源的大多数，11个赛弗特1型星系中有9个也超过该阈值。根据光变判定标准，共识别出30个可变ULIRGs。

**3.2 光变幅度与超新星比较**

30个可变ULIRGs的ΔL_W2范围从7×10⁴²到5×10⁴⁴ erg s^?1。其中25个源的ΔL_W2>10⁴³ erg s^?1，超出已知最亮MIR超新星光度的约5倍，表明这些光变无法由恒星活动解释，而更可能源于AGN固有光度变化或TDEs。样本中包含先前报道的TDE候选体IRAS F01004?2237和IRAS 05189?2524，以及硬X射线探测到的AGN（如UGC 05101）。

**3.3 AGN特征的多波段验证**

所有25个ΔL_W2>10⁴³ erg s^?1的可变ULIRGs均满足至少一种AGN判据：14个被光学分类为赛弗特1型或2型星系；3个非赛弗特星系通过X射线观测被确认为AGN（IRAS 20551?4250、IRAS 12514+1027和IRAS 09320+6134/UGC 05101）；除IRAS 01388?4618外，其余24个源满足Stern等人（2012）的MIR颜色判据W1?W2>0.8，而IRAS 01388?4618的Spitzer光谱证实了MIR光谱中的AGN贡献。这表明即使MIR光变起源于TDEs，所涉及的黑洞也很可能是活跃的AGN。

**3.4 辐射区域尺度约束**

对于在数年时标上呈现明显上升和下降光变曲线的ULIRGs（如IRAS 01572+0009、IRAS 07246+6125、IRAS 08559+1053），光变时标隐含R?几秒差距的致密辐射区域，与AGN环状尘埃内区的典型热尘埃分布尺度一致。仅呈现单调上升或下降的源无法设定上限。

**3.5 光变发生率**

155个ULIRGs中识别出30个可变源，光变发生率约为20%，比矮星系中的光变发生率（?1%）高一个数量级。另有22个ULIRGs满足r>0.75但未达到方差阈值，若将其计入，发生率可能高达约30%。这表明保守的系统误差剔除可能遗漏了一些实际可变源。

**4. 讨论与结论**

**4.1 光变起源的物理机制**

MIR光变可源于热辐射（尘埃对UV/光学的再处理）或非热辐射（相对论性喷流的同步辐射）。少数源（如IRAS 11223?1244）表现出剧烈的上下波动光变，可能来自喷流的非热辐射；大多数源的光变曲线较为平滑，与热再处理过程更一致。

对于平滑的MIR光变，由于这些系统的大MIR光度和空间尺度，由遮挡引起的光变可能性较低，变化更可能源于吸积功率的内在变化。相关物理机制包括AGN吸积流的随机涨落、TDEs以及CLAGNs中讨论的吸积盘状态转变。随机涨落表现为剧烈的上下波动，而TDEs和CLAGNs通常表现为单一耀发事件。

**4.2 个别源的特殊性质**

IRAS 05189?2524在初始爆发后的最后几个历元还表现出光度增加，可能对应另一个瞬变事件。IRAS 00456?2904和IRAS 17028+5817缺乏光学、X射线、MIR颜色或MIR光谱AGN特征，但前者通过ALMA观测发现升高的HCN/HCO⁺通量比，暗示存在光学上难以探测的AGN；后者表现出两个明显的瞬变事件，光变曲线表面类似于重复性部分TDEs。IRAS 01003?2238和IRAS 01572+0009（Mrk 1014）表现出数年时标上的重复增亮和减弱，两者均被Efstathiou等人（2022）基于SED拟合提出可能拥有双AGN。

**4.3 尘埃环回响模型**

对于IRAS 01003?2238、IRAS 01572+0009和IRAS 12071?0444这三个具有缓慢上升或双峰MIR光变曲线的源，研究人员构建了考虑几何延迟的尘埃环回响模型。中心源模型能够再现三个耀发的总体光变形状，但无法完全再现IRAS 12071?0444的双峰结构；而偏离中心的源模型能更好地再现该双峰发射特征。该模型说明，即使单一的UV/光学耀发，当照射源偏离环中心时也能产生双峰MIR特征。利用尘埃升华半径关系估算的UV光度（L_UV,sub）与总IR光度之比L_IR/L_UV,sub为0.1–0.4，表明推断的UV光度在能量上足以驱动观测到的IR光度。

**5. 研究结论**

研究人员分析了Nardini等人（2010）中164个ULIRGs的11年NEOWISE测光数据，以MIR光变作为探测掩埋型AGN的探针。研究在30个ULIRGs中探测到显著的光变，其中25个源的光变幅度超过ΔL_W2=10⁴³ erg s^?1，比已知MIR中最亮的SNe更亮。这些源包括先前报道的TDE候选体和硬X射线探测到的AGN。由于光变在数年时标上被观测到，对于光度同时出现增加和减少的ULIRGs，MIR辐射区域被约束为?几秒差距的致密结构。大光变幅度无法由SNe解释，极可能源于AGN固有光度变化或发生在宿主AGN系统中的TDEs。所有ΔL_W2>10⁴³ erg s^?1的可变ULIRGs均表现出至少一种独立的AGN特征，支持了这一解释。这项工作证明了MIR光变可以开辟揭示ULIRGs中掩埋核活动的新领域。