2025-11-21T00:25:15.653809

Searching for Electromagnetic Counterpart Candidates to GW231123

He, Zhu, Liu et al.
The detection of GW231123, a gravitational-wave (GW) event with exceptionally massive and rapidly spinning black holes, suggests the possible formation within an active galactic nucleus (AGN) disk, which provides a favorable environment for potentially generating an observable electromagnetic (EM) counterpart. We conduct a search for such a counterpart by crossmatching the GW localization with a comprehensive catalog of AGN flares from the Zwicky Transient Facility. Our analysis yields six plausible optical flare candidates that are spatially and temporally coincident with GW231123 and exhibit significant deviations from their AGN baseline flux. Although these candidates represent a crucial first step, their true nature remains inconclusive. Confirming any one of these flares via future observations would provide a landmark validation of the AGN formation channel and unlock the multi-messenger potential of this extraordinary merger.
academic

Searching for Electromagnetic Counterpart Candidates to GW231123

基本信息

  • 论文ID: 2511.05144
  • 标题: Searching for Electromagnetic Counterpart Candidates to GW231123
  • 作者: Lei He, Liang-Gui Zhu, Zheng-Yan Liu, Rui Niu, Chao Wei, Bing-Zhou Gao, Ming-Shen Zhou, Run-Duo Liang, Ken Chen, Jian-Min Wang, Ning Jiang, Zhen-Yi Cai, Ji-An Jiang, Zi-Gao Dai, Ye-Fei Yuan, Jian Li, Wen Zhao
  • 机构: 中国科学技术大学天文系、国家天文台、高能物理研究所、贵州大学
  • 分类: astro-ph.HE (高能天体物理)
  • 提交时间: 2025年11月7日
  • 论文链接: https://arxiv.org/abs/2511.05144

摘要

GW231123是LIGO-Virgo-KAGRA探测到的一个极端引力波事件，其双黑洞系统具有异常大的质量和极高的自旋。这些特性表明该系统可能形成于活动星系核(AGN)盘中，这为产生可观测的电磁对应体提供了有利环境。本研究通过将引力波定位与Zwicky瞬变设施(ZTF)的AGN耀发目录进行交叉匹配，搜索潜在的电磁对应体。分析发现6个可信的光学耀发候选体，它们在空间和时间上与GW231123重合，并显示出相对于AGN基线流量的显著偏离。尽管这些候选体是重要的第一步，但其真实性质尚未确定。通过未来观测确认任何一个耀发都将为AGN形成渠道提供里程碑式的验证，并开启这次非凡合并的多信使研究潜力。

研究背景与动机

1. 要解决的核心问题

本研究旨在为GW231123这一极端引力波事件寻找可能的电磁对应体。GW231123的双黑洞系统具有以下极端特性：

  • 超大质量: 主星质量137⁺²²₋₁₇M⊙，次星质量103⁺²⁰₋₅₂M⊙，是迄今探测到的最大质量双黑洞合并
  • 极高自旋: χ₁ = 0.9⁺⁰·¹⁰₋₀.₁₉，χ₂ = 0.8⁺⁰·²⁰₋₀.₅₁，是引力波观测中报告的最高自旋值

2. 问题的重要性

科学意义：

  • 形成机制验证: 双黑洞的质量落在理论预测的对不稳定性质量间隙(pair-instability mass gap)内，挑战标准恒星演化模型，需要替代形成渠道
  • AGN环境探测: 如果确认电磁对应体，将揭示AGN盘中恒星级黑洞的起源和演化
  • 多信使天文学: 确认的对应体可作为"标准汽笛"，通过结合引力波和电磁观测独立测量哈勃常数等宇宙学参数

理论预期：

  • 层级合并(hierarchical mergers)被认为是最有说服力的解释，黑洞通过连续合并组装，自然产生能跨越质量间隙的大质量黑洞
  • AGN盘提供了组装如此大质量黑洞的理想环境：致密气体促进双星形成和演化，中心超大质量黑洞的深势阱保留合并残余，气体摩擦将倾斜的残余拖回盘中

3. 现有方法的局限性

历史搜索的挑战：

  • 之前观测运行中发现了几个候选体(如与GW190521相关的ZTF19abanrhr)，但确认确定性关联仍是重大挑战
  • 文献中存在相互矛盾的解释
  • He et al. (2025a)的重新评估发现，7个先前提出的候选体中只有2个在三年后仍与相关引力波事件统计一致

4. 研究动机

GW231123的极端特性使其成为寻找电磁对应体的理想目标：

  • 大质量黑洞产生更明亮的发射
  • 高自旋与层级合并预测一致
  • AGN盘环境预期产生可观测的电磁辐射

核心贡献

  1. 系统性搜索: 首次对GW231123进行系统的电磁对应体搜索，使用ZTF数据发布23构建的AGN耀发目录(AGNFCC)
  2. 候选体识别: 发现6个在空间和时间上与GW231123重合的可信光学耀发候选体
  3. 贝叶斯评估: 计算每个候选体的几率比(odds ratio)，大多数候选体在至少一个波形模型下显示强证据(lnO > 5)
  4. 物理约束: 通过简单的Bondi吸积模型推导AGN盘的物理参数(气体密度ρ和标高H)，与理论预测一致
  5. 方法学改进:
    • 使用强制测光(forced photometry)减少宿主星系污染
    • 扩展时间基线至2025年10月，比之前研究延长近一年
    • 应用严格的耀发检测算法，最小化AGN本征变化的污染

方法详解

任务定义

输入：

  • GW231123的天空定位和距离后验分布(来自5个不同波形模型)
  • AGN耀发目录(AGNFCC)，包含源位置、红移、耀发峰值时间
  • ZTF强制测光数据(2018年3月至2025年10月)

输出：

  • 与GW231123在空间和时间上重合的AGN耀发候选体列表
  • 每个候选体与引力波事件物理关联的几率比

约束条件：

  • 耀发峰值必须在GW触发之后
  • 位于99%置信定位区域内
  • 上升和衰减时标在3-200天之间
  • 排除已知瞬变源(超新星、潮汐瓦解事件、耀变体)

搜索流程

1. 数据准备

引力波数据：

  • 事件时间: 2023年11月23日13:54:30 UTC
  • 探测器: LIGO Hanford和Livingston
  • 信噪比: ~20
  • 90%置信区域: 400-1000 deg²
  • 光度距离: 数百至数千Mpc(依赖波形模型)

AGN耀发目录：

  • 来源: AGNFCC (He et al. 2025b)
  • 基于ZTF数据发布23
  • 包含截至2024年10月31日的AGN耀发
  • 对无红移源，补充Pan-STARRS1的测光红移

2. 交叉匹配

使用ligo.skymap Python包的crossmatch方法：

  • 有红移源: 在3D置信体积内匹配
  • 无红移源: 在2D置信区域内匹配
  • 对5个波形模型分别进行匹配
  • 要求耀发峰值时间在GW触发之后

初步结果: 找到166个候选AGN耀发

3. 光变曲线分析

数据获取：

  • 从ZTF强制测光服务检索g和r波段历史光变曲线
  • 时间跨度: 2018年3月至2025年10月
  • 优势: 减少宿主星系污染，更长时间基线

耀发检测算法 (He et al. 2025b)：

  1. 用高斯过程拟合长期AGN变化
  2. 模拟光变曲线估计耀发源自本征变化的概率
  3. 仅保留低概率事件，最小化AGN涨落污染

耀发建模： 采用高斯上升-指数衰减模型：

r_0 + A \exp\left(-\frac{(t-t_0)^2}{2t_g^2}\right), & t \leq t_0 \\ r_0 + A \exp\left(-\frac{t-t_0}{t_e}\right), & t \geq t_0 \end{cases}$$ 其中： - $r_0$: 基线流量 - $A$: 耀发振幅 - $t_0$: 峰值时间 - $t_g$, $t_e$: 上升和衰减时标 #### 4. 候选体筛选 **严格标准**： 1. 排除AGNFCC中标记的已知瞬变源 2. 耀发仅在GW触发后开始上升 3. 光变曲线中表现为单一、孤立的爆发 4. 上升和衰减时标在3-200天之间 - 排除短暂假象(考虑ZTF巡天节奏) - 排除可能归因于AGN本征活动的长期变化 ### 贝叶斯评估方法 #### 几率比计算 对每个候选体，计算关联模型vs偶然重合模型的几率比： $$O = \frac{P(\text{关联}|\text{数据})}{P(\text{偶然}|\text{数据})}$$ **输入参数**： - GW源的后验分布：赤经、赤纬、光度距离、主星质量 - 宿主AGN位置 - 99%置信体积内的AGN数量 - 耀发的显著性 **判据**： - lnO > 5 (O > 150): 非常强的关联证据 ## 实验设置 ### 数据集 **引力波数据**： - 来源: LIGO-Virgo-KAGRA合作组织 - 波形模型: 5个 - IMRPhenomTPHM - IMRPhenomXO4a - IMRPhenomXPHM-SpinTaylor - NRSur7dq4 - SEOBNRv5PHM - 模型间存在显著差异(特别是距离后验) **光学数据**： - 巡天: Zwicky Transient Facility (ZTF) - 数据发布: DR23 - 波段: g和r - 时间范围: 2018年3月至2025年10月 - 数据类型: 强制测光 **辅助光谱数据**： - SDSS DR19 - DESI DR1 - 用于3个候选体的超大质量黑洞质量测量 ### 物理模型参数 **反冲速度**： - GW231123: $v_k = 967^{+829}_{-857}$ km/s **束缚气体时标**： $$t_{\text{bound}} \sim 20 \text{ days} \left(\frac{M_{\text{BBH}}}{100M_\odot}\right)\left(\frac{v_k}{200\text{ km/s}}\right)^{-3}$$ **Bondi-Hoyle-Lyttleton吸积**： $$\dot{M}_{\text{BHL}} = \frac{4\pi G^2 M_{\text{BBH}}^2 \rho}{v_{\text{rel}}^3}$$ $$L_{\text{BHL}} = \eta \dot{M}_{\text{BHL}} c^2$$ 其中$\eta \sim 0.1$为辐射效率 **盘逃逸时标**： $$t_{\text{exit}} = \frac{H\sqrt{2\ln(\tau_{\text{mp}})}}{v_k}$$ 假设$\tau_{\text{mp}} \sim 10^{4.5}$(合并点的光学深度) ## 实验结果 ### 主要结果：6个候选体 | AGN名称 | 红移 | $t_{\text{exit}}$ (天) | $t_g$ (天) | $t_e$ (天) | $\log_{10}(E_{\text{tot}})$ (erg) | |---------|------|---------|-------|-------|------------| | J164911.49+492438.9 | 0.344 | 274.8 | 79.4 | 140.3 | 51.5 | | J131824.85+253342.0 | 0.424 | 76.9 | 42.3 | 54.5 | 51.6 | | J121530.73+080220.4 | 0.802 | 28.9 | 8.4 | 74.0 | 52.4 | | J203336.84-275303.9 | 0.149 | 194.7 | 3.2 | 23.2 | 49.4 | | J121013.51+012337.5 | 0.171 | 47.2 | 4.4 | 14.3 | 49.7 | | J140528.43+222333.2 | 0.156 | 185.8 | 33.6 | 108.7 | 51.0 | **注**: $t_{\text{exit}}$为GW触发到耀发峰值的时间延迟 ### 几率比分析 **强证据候选体**(至少一个波形模型下lnO > 5)： - **J203336.84-275303.9**: 最高几率比(IMRPhenomXPHM下lnO = 10.1) - **J131824.85+253342.0**: 多个模型下lnO = 6.8-9.9 - **J164911.49+492438.9**: lnO = 6.2-8.3 - **J121530.73+080220.4**: IMRPhenomXO4a下lnO = 8.3 - **J121013.51+012337.5**: NRSur7dq4下lnO = 8.3 **弱证据候选体**： - **J140528.43+222333.2**: 最高lnO = 4.3(可能由于相对于其本征变化，耀发不够显著) ### 物理参数推导 通过Bondi吸积模型推导的AGN盘参数： | AGN名称 | $\log_{10}(M_{\text{SMBH}})$ (M⊙) | ρ (10⁻¹⁰ g/cm³) | H (rg) | |---------|------------|------|--------| | J164911.49+492438.9 | 8.3 | 3.2 | 33 | | J131824.85+253342.0 | (8.0) | 7.6 | (9.4) | | J121530.73+080220.4 | (8.0) | 241 | (3.5) | | J203336.84-275303.9 | (8.0) | 0.63 | (24) | | J121013.51+012337.5 | 6.6 | 0.91 | 5.8 | | J140528.43+222333.2 | 7.5 | 2.4 | 23 | **注**: 括号内值使用基准SMBH质量$10^8 M_\odot$ **一致性检验**： 推导的ρ和H组合与AGN盘结构模型预测一致(Sirko & Goodman 2003; Thompson et al. 2005) ### 光变曲线特征 **颜色演化**： - **J131824.85+253342.0**和**J121530.73+080220.4**: 近似恒定颜色，与简单BHL吸积场景一致 - 其他候选体: 显示颜色变化，可能反映更复杂的物理条件 **时间延迟解释**： 对于GW231123，使用残余质量~200M⊙和中值反冲速度~1000 km/s： $$t_{\text{bound}} \sim 0.4 \text{ 天}$$ 这解释了J131824.85+253342.0、J121530.73+080220.4和J121013.51+012337.5的近即时起始。但反冲速度的不确定性很大，不能排除延迟更长的耀发。 ### ALeRCE分类器结果 - **J203336.84-275303.9**: 分类为超新星 - **J121013.51+012337.5**: 未分类 - 其余: 分类为AGN **注意**: 该分类器仅基于测光数据，不包含BBH合并类别 ## 相关工作 ### 引力波源的形成渠道 **层级合并**： - Li et al. (2025): GW231123的成分源自~6和~4个第一代黑洞的连续合并 - Delfavero et al. (2025): 最符合第四代和第三代黑洞的合并 - Li & Fan (2025a): 贝叶斯框架下2G+2G合并场景优于2G+1G或3G+2G **其他场景**： - 极大质量快速旋转氦核的孤立双星演化(Croon et al. 2025; Gottlieb et al. 2025) - 化学均匀演化(Stegmann et al. 2025; Popa & de Mink 2025) - 持续吸积(Bartos & Haiman 2025) - 强引力透镜(Shan et al. 2025) - 第三族星形成(Tanikawa et al. 2025; Liu et al. 2025; Wu et al. 2025) - 原初黑洞(Yuan et al. 2025; Luca et al. 2025; Nojiri & Odintsov 2025) ### AGN盘中的BBH合并 **环境优势** (McKernan et al. 2012; Tagawa et al. 2020)： - 致密气体促进双星形成和演化(Li et al. 2023; Dittmann et al. 2024) - 深势阱保留合并残余(Varma et al. 2022) - 气体摩擦将倾斜残余拖回盘中(Fabj et al. 2020) **电磁发射机制** (Bartos et al. 2017)： 1. **冲压剥离**: 残余束缚气体的剥离(McKernan et al. 2019) 2. **超爱丁顿吸积**: 驱动强大流出爆发(Wang et al. 2021; Kimura et al. 2021; Rodrı́guez-Ramı́rez et al. 2025) 3. **喷流发射**: 吸积残余驱动的喷流(Tagawa et al. 2024; Rodrı́guez-Ramı́rez et al. 2024; Chen & Dai 2024) ### 历史搜索 **先前候选体** (Graham et al. 2020, 2023; Cabrera et al. 2024)： - **ZTF19abanrhr**与GW190521: 最著名案例，但解释存在争议(Ashton et al. 2021; Palmese et al. 2021; Morton et al. 2023) - He et al. (2025a)重新评估: 7个候选体中仅2个在三年后仍统计一致 ### 科学价值 **多信使天文学** (Vajpeyi et al. 2022)： - 揭示AGN盘中恒星级黑洞的起源和演化 - 洞察AGN环境的物理过程 **标准汽笛** (Mukherjee et al. 2020; Chen et al. 2022; Li & Fan 2025b)： - 独立测量哈勃常数 - 宇宙学参数约束 ## 结论与讨论 ### 主要结论 1. **成功识别**: 通过系统搜索，发现6个与GW231123在空间和时间上重合的AGN耀发候选体 2. **强统计支持**: 大多数候选体在贝叶斯框架下显示强关联证据(lnO > 5)，主要由GW231123的大主星质量驱动 3. **物理一致性**: 通过简单Bondi吸积模型推导的AGN盘参数(ρ和H)与理论预测一致 4. **性质未定**: 由于缺乏同期光谱数据，这些耀发的物理起源仍不确定 5. **最多一个真实**: 6个候选体中最多一个是GW231123的真实电磁对应体，其余可能是无关背景事件 ### 局限性 **观测限制**： 1. **无同期光谱**: 无法通过不对称宽线轮廓等特征确认BBH起源 2. **单次观测**: 尚未观测到预期的重现耀发 3. **波形不确定性**: 不同模型的距离后验存在显著差异，可能源于高自旋下当前信号模型的不确定性 **方法限制**： 1. **简化模型**: 采用的Bondi吸积模型假设反冲垂直于盘，实际情况可能更复杂 2. **假阳性风险**: TDE、超新星或千新星也可能在AGN盘中产生类似耀发 3. **分类不确定**: ALeRCE分类器将J203336.84-275303.9分类为超新星，需进一步验证 **统计限制**： 1. **小样本**: 仅6个候选体，统计显著性有限 2. **背景污染**: 166个初始候选体中仅6个通过严格筛选，背景率较高 ### 未来方向 **关键后续观测**： 1. **持续光学监测**： - 使用Wide Field Survey Telescope和Vera C. Rubin天文台 - 探测预期的重现耀发 - 时标: $\sim 1.6 \text{ yr} (M_{\text{SMBH}}/10^8 M_\odot)(r/10^3 r_g)^{3/2}$ - 通过测量重现间隔约束AGN盘内的合并位置 2. **快速多波段和光谱后续**： - 耀发阶段的及时光谱观测 - 识别不对称宽线特征(BBH合并证据) - 排除超新星或TDE的典型光谱特征 3. **高能观测**： - Fermi伽马射线空间望远镜 - Einstein Probe - X射线和γ射线监测约束物理环境 - 即使非探测也提供有价值信息 **理论发展**： 1. 改进高自旋BBH的波形模型 2. 更精细的AGN盘-BBH相互作用模拟 3. 多维辐射转移模型 **长期目标**： - 建立BBH合并电磁对应体的统计样本 - 验证AGN辅助合并场景 - 开启多信使天文学新时代 ## 深度评价 ### 优点 **科学意义**： 1. **及时性强**: 针对最新、最极端的引力波事件(GW231123)进行首次系统搜索 2. **目标明确**: GW231123的超大质量和高自旋使其成为寻找电磁对应体的理想目标 3. **多信使价值**: 确认的对应体将是多信使天文学的重要里程碑 **方法学优势**： 1. **数据质量**: 使用强制测光减少宿主星系污染，扩展时间基线至2025年10月 2. **严格筛选**: 多层筛选标准(时间延迟、孤立爆发、时标约束)有效减少假阳性 3. **贝叶斯评估**: 定量评估关联可能性，提供统计支持 4. **物理建模**: 通过Bondi吸积模型连接观测和理论，推导物理参数 **分析全面**： 1. **多波形**: 考虑5个不同波形模型，处理系统不确定性 2. **光谱验证**: 利用SDSS和DESI档案光谱测量SMBH质量 3. **颜色演化**: 分析g-r颜色变化提供额外物理洞察 **写作清晰**： 1. 结构合理，逻辑清晰 2. 方法描述详细，可复现性强 3. 图表信息丰富(光变曲线、天区图、光谱) ### 不足 **观测局限**： 1. **缺乏确认**: 无同期光谱，无法确定性排除其他瞬变类型 2. **单次观测**: 未观测到重现耀发，无法验证AGN盘场景的关键预测 3. **多波段缺失**: 仅有光学数据，缺乏X射线、射电等多波段约束 **方法问题**： 1. **模型简化**: - 假设反冲垂直于盘(最快逃逸路径)，可能不现实 - 忽略盘的不均匀性和动力学复杂性 - 恒定颜色假设对某些候选体不成立 2. **参数不确定**: - 3个源使用基准SMBH质量$10^8 M_\odot$，引入系统误差 - 光学深度$\tau_{\text{mp}} \sim 10^{4.5}$的假设缺乏验证 3. **分类矛盾**: J203336.84-275303.9被ALeRCE分类为超新星，但仍作为候选体保留 **统计问题**： 1. **样本小**: 仅6个候选体，难以进行统计显著性检验 2. **背景高**: 初始166个候选体中96%被排除，表明背景污染严重 3. **多重比较**: 对5个波形模型分别计算几率比，未进行多重比较校正 **讨论不足**： 1. 未充分讨论J203336.84-275303.9的超新星分类问题 2. 对不同候选体之间的相对可信度缺乏排序 3. 未量化确认任何一个候选体所需的后续观测时间/资源 ### 影响力 **对领域的贡献**： 1. **方法学**: 建立了GW事件电磁对应体搜索的系统流程 2. **目标列表**: 为后续观测提供了明确的6个目标 3. **参数空间**: 扩展了已搜索的BBH质量和自旋参数空间 **实用价值**： 1. **观测指导**: 为未来多波段和光谱后续提供具体目标 2. **理论检验**: 提供检验AGN盘BBH合并场景的机会 3. **方法可复用**: 搜索流程可应用于未来引力波事件 **潜在影响**： - **若确认**: 将是多信使天文学的重大突破，验证AGN形成渠道 - **若排除**: 将约束AGN盘中电磁发射的物理模型 - **方法传承**: AGNFCC目录和搜索流程可持续应用 **可复现性**： 1. **数据公开**: ZTF数据公开可获取 2. **代码可用**: 使用公开软件包(ligo.skymap, PyQSOFit) 3. **方法详细**: 搜索流程和筛选标准描述清晰 4. **目录共享**: AGNFCC目录在GitHub公开 ### 适用场景 **理想应用**： 1. **极端质量BBH**: 类似GW231123的超大质量、高自旋系统 2. **近距离事件**: 光度距离<2000 Mpc，确保光学可探测性 3. **良好定位**: 90%置信区域<1000 deg²，便于交叉匹配 4. **AGN环境**: 预期形成于AGN盘的系统 **扩展应用**： 1. **中等质量BBH**: 50-150 M⊙的系统 2. **中子星-黑洞合并**: 若发生在AGN盘中 3. **历史事件**: 应用于过去探测的引力波事件 **不适用场景**： 1. **低质量系统**: <30 M⊙的BBH，电磁发射可能太弱 2. **极远距离**: >5000 Mpc，光学信号低于探测阈值 3. **差定位**: >10000 deg²，交叉匹配效率低 4. **场环境**: 非AGN环境中的合并，不预期强电磁发射 ## 参考文献 ### 关键引力波论文 - The LIGO Scientific Collaboration et al. (2025): GW231123的发现和性质 - Abbott et al. (2020): 早期BBH探测和质量间隙讨论 ### AGN盘BBH理论 - McKernan et al. (2012, 2019): AGN盘中BBH形成和电磁发射 - Bartos et al. (2017): AGN环境中电磁对应体的理论预期 - Wang et al. (2021): 超爱丁顿吸积和流出模型 ### 电磁对应体搜索 - Graham et al. (2020, 2023): ZTF19abanrhr和GW190521关联 - He et al. (2025a): 历史候选体的重新评估 - He et al. (2025b): AGNFCC目录构建 ### 层级合并研究 - Gerosa & Fishbach (2021): 层级合并黑洞的自旋分布 - Li et al. (2025): GW231123的层级合并解释 ### AGN盘结构 - Sirko & Goodman (2003): AGN盘理论模型 - Thompson et al. (2005): 盘的垂直结构和密度分布 --- **总体评价**: 这是一篇及时、系统、方法严谨的观测论文，为GW231123这一极端引力波事件提供了首次电磁对应体搜索。尽管存在观测和方法上的局限，6个候选体的发现为未来后续观测提供了明确目标。确认任何一个候选体都将是多信使天文学的重大突破，验证AGN辅助BBH合并的理论预期。论文的方法学贡献和可复现性使其对领域具有持续价值。