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

Basic Information

  • Paper ID: 2511.05144
  • Title: Searching for Electromagnetic Counterpart Candidates to GW231123
  • Authors: 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
  • Institutions: Department of Astronomy, University of Science and Technology of China; National Astronomical Observatories; Institute of High Energy Physics; Guizhou University
  • Classification: astro-ph.HE (High-Energy Astrophysics)
  • Submission Date: November 7, 2025
  • Paper Link: https://arxiv.org/abs/2511.05144

Abstract

GW231123 is an extreme gravitational wave event detected by LIGO-Virgo-KAGRA, whose binary black hole system exhibits anomalously large masses and extremely high spins. These characteristics suggest the system may have formed within an active galactic nucleus (AGN) disk, providing a favorable environment for producing observable electromagnetic counterparts. This study searches for potential electromagnetic counterparts by cross-matching the gravitational wave localization with the AGN flare catalog from the Zwicky Transient Facility (ZTF). The analysis identifies six credible optical flare candidates that spatially and temporally coincide with GW231123, showing significant deviations from the AGN baseline flux. Although these candidates represent an important first step, their true nature remains undetermined. Future observational confirmation of any flare would provide milestone verification for the AGN formation channel and unlock the multi-messenger research potential of this extraordinary merger.

Research Background and Motivation

1. Core Problem to be Addressed

This study aims to identify possible electromagnetic counterparts to GW231123, an extreme gravitational wave event. The binary black hole system of GW231123 exhibits the following extreme characteristics:

  • Supermassive masses: Primary mass 137⁺²²₋₁₇M⊙, secondary mass 103⁺²⁰₋₅₂M⊙, representing the most massive binary black hole merger detected to date
  • Extremely high spins: χ₁ = 0.9⁺⁰·¹⁰₋₀.₁₉, χ₂ = 0.8⁺⁰·²⁰₋₀.₅₁, the highest spin values reported in gravitational wave observations

2. Importance of the Problem

Scientific Significance:

  • Formation mechanism verification: The binary black hole masses fall within the theoretically predicted pair-instability mass gap, challenging standard stellar evolution models and requiring alternative formation channels
  • AGN environment detection: If electromagnetic counterparts are confirmed, it will reveal the origin and evolution of stellar-mass black holes in AGN disks
  • Multi-messenger astronomy: Confirmed counterparts can serve as "standard sirens," enabling independent measurement of cosmological parameters such as the Hubble constant through combined gravitational wave and electromagnetic observations

Theoretical Expectations:

  • Hierarchical mergers are considered the most compelling explanation, where black holes assemble through successive mergers, naturally producing massive black holes capable of spanning the mass gap
  • AGN disks provide an ideal environment for assembling such massive black holes: dense gas promotes binary formation and evolution, the deep potential well of the central supermassive black hole retains merger remnants, and gas friction drags tilted remnants back into the disk

3. Limitations of Existing Methods

Historical Search Challenges:

  • Previous observing runs identified several candidates (e.g., ZTF19abanrhr associated with GW190521), but confirming definitive associations remains a major challenge
  • Contradictory interpretations exist in the literature
  • He et al. (2025a)'s reassessment found that only 2 of 7 previously proposed candidates remained statistically consistent with the relevant gravitational wave events three years later

4. Research Motivation

The extreme characteristics of GW231123 make it an ideal target for electromagnetic counterpart searches:

  • More massive black holes produce brighter emissions
  • High spin is consistent with hierarchical merger predictions
  • AGN disk environment is expected to produce observable electromagnetic radiation

Core Contributions

  1. Systematic Search: First systematic electromagnetic counterpart search for GW231123 using the AGN Flare Candidate Catalog (AGNFCC) constructed from ZTF data release 23
  2. Candidate Identification: Discovery of six credible optical flare candidates that spatially and temporally coincide with GW231123
  3. Bayesian Assessment: Calculation of odds ratios for each candidate, with most candidates showing strong evidence (lnO > 5) under at least one waveform model
  4. Physical Constraints: Derivation of AGN disk physical parameters (gas density ρ and scale height H) through simple Bondi accretion models, consistent with theoretical predictions
  5. Methodological Improvements:
    • Use of forced photometry to reduce host galaxy contamination
    • Extended time baseline to October 2025, approximately one year longer than previous studies
    • Application of rigorous flare detection algorithms to minimize contamination from AGN intrinsic variability

Detailed Methodology

Task Definition

Inputs:

  • Sky localization and distance posterior distribution of GW231123 (from 5 different waveform models)
  • AGN Flare Candidate Catalog (AGNFCC), containing source positions, redshifts, and flare peak times
  • ZTF forced photometry data (March 2018 to October 2025)

Outputs:

  • List of AGN flare candidates coinciding spatially and temporally with GW231123
  • Odds ratios for physical association of each candidate with the gravitational wave event

Constraints:

  • Flare peak must occur after GW trigger
  • Located within 99% confidence localization region
  • Rise and decay timescales between 3-200 days
  • Exclusion of known transient sources (supernovae, tidal disruption events, blazars)

Search Workflow

1. Data Preparation

Gravitational Wave Data:

  • Event time: November 23, 2023, 13:54:30 UTC
  • Detectors: LIGO Hanford and Livingston
  • Signal-to-noise ratio: ~20
  • 90% confidence region: 400-1000 deg²
  • Luminosity distance: hundreds to thousands of Mpc (dependent on waveform model)

AGN Flare Catalog:

  • Source: AGNFCC (He et al. 2025b)
  • Based on ZTF data release 23
  • Contains AGN flares through October 31, 2024
  • For sources without redshift, supplemented with Pan-STARRS1 photometric redshifts

2. Cross-Matching

Using the ligo.skymap Python package crossmatch method:

  • Sources with redshift: Matching within 3D confidence volume
  • Sources without redshift: Matching within 2D confidence region
  • Separate matching for 5 waveform models
  • Requirement that flare peak time occurs after GW trigger

Initial Result: 166 candidate AGN flares identified

3. Light Curve Analysis

Data Retrieval:

  • Historical light curves in g and r bands retrieved from ZTF forced photometry service
  • Time span: March 2018 to October 2025
  • Advantage: Reduced host galaxy contamination, longer time baseline

Flare Detection Algorithm (He et al. 2025b):

  1. Fit long-term AGN variability with Gaussian process
  2. Simulate light curves to estimate probability that flare originates from intrinsic variability
  3. Retain only low-probability events to minimize AGN fluctuation contamination

Flare Modeling: Gaussian rise-exponential decay model:

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}$$ Where: - $r_0$: baseline flux - $A$: flare amplitude - $t_0$: peak time - $t_g$, $t_e$: rise and decay timescales #### 4. Candidate Selection **Strict Criteria**: 1. Exclude known transient sources marked in AGNFCC 2. Flare rise begins only after GW trigger 3. Light curve exhibits single, isolated outburst 4. Rise and decay timescales between 3-200 days - Excludes brief artifacts (accounting for ZTF survey cadence) - Excludes long-term variations potentially attributable to AGN intrinsic activity ### Bayesian Assessment Method #### Odds Ratio Calculation For each candidate, calculate the odds ratio of association model versus chance coincidence model: $$O = \frac{P(\text{association}|\text{data})}{P(\text{chance}|\text{data})}$$ **Input Parameters**: - Posterior distribution of GW source: right ascension, declination, luminosity distance, primary mass - Host AGN position - Number of AGN within 99% confidence volume - Flare significance **Criterion**: - lnO > 5 (O > 150): Very strong association evidence ## Experimental Setup ### Datasets **Gravitational Wave Data**: - Source: LIGO-Virgo-KAGRA Collaboration - Waveform models: 5 - IMRPhenomTPHM - IMRPhenomXO4a - IMRPhenomXPHM-SpinTaylor - NRSur7dq4 - SEOBNRv5PHM - Significant differences exist between models (particularly in distance posteriors) **Optical Data**: - Survey: Zwicky Transient Facility (ZTF) - Data Release: DR23 - Bands: g and r - Time Range: March 2018 to October 2025 - Data Type: Forced photometry **Auxiliary Spectroscopic Data**: - SDSS DR19 - DESI DR1 - Used for supermassive black hole mass measurements for 3 candidates ### Physical Model Parameters **Recoil Velocity**: - GW231123: $v_k = 967^{+829}_{-857}$ km/s **Bound Gas Timescale**: $$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 Accretion**: $$\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$$ Where $\eta \sim 0.1$ is the radiative efficiency **Disk Escape Timescale**: $$t_{\text{exit}} = \frac{H\sqrt{2\ln(\tau_{\text{mp}})}}{v_k}$$ Assuming $\tau_{\text{mp}} \sim 10^{4.5}$ (optical depth at merger point) ## Experimental Results ### Main Results: Six Candidates | AGN Name | Redshift | $t_{\text{exit}}$ (days) | $t_g$ (days) | $t_e$ (days) | $\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 | **Note**: $t_{\text{exit}}$ is the time delay from GW trigger to flare peak ### Odds Ratio Analysis **Strong Evidence Candidates** (lnO > 5 under at least one waveform model): - **J203336.84-275303.9**: Highest odds ratio (lnO = 10.1 under IMRPhenomXPHM) - **J131824.85+253342.0**: lnO = 6.8-9.9 across multiple models - **J164911.49+492438.9**: lnO = 6.2-8.3 - **J121530.73+080220.4**: lnO = 8.3 under IMRPhenomXO4a - **J121013.51+012337.5**: lnO = 8.3 under NRSur7dq4 **Weak Evidence Candidate**: - **J140528.43+222333.2**: Maximum lnO = 4.3 (possibly due to flare being insufficiently prominent relative to its intrinsic variability) ### Derived Physical Parameters AGN disk parameters derived through Bondi accretion model: | AGN Name | $\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 | **Note**: Values in parentheses use fiducial SMBH mass $10^8 M_\odot$ **Consistency Check**: Derived combinations of ρ and H are consistent with AGN disk structure model predictions (Sirko & Goodman 2003; Thompson et al. 2005) ### Light Curve Characteristics **Color Evolution**: - **J131824.85+253342.0** and **J121530.73+080220.4**: Approximately constant color, consistent with simple BHL accretion scenario - Other candidates: Show color variations, possibly reflecting more complex physical conditions **Time Delay Interpretation**: For GW231123, using remnant mass ~200M⊙ and median recoil velocity ~1000 km/s: $$t_{\text{bound}} \sim 0.4 \text{ days}$$ This explains the near-immediate onset for J131824.85+253342.0, J121530.73+080220.4, and J121013.51+012337.5. However, the large uncertainty in recoil velocity cannot exclude flares with longer delays. ### ALeRCE Classifier Results - **J203336.84-275303.9**: Classified as supernova - **J121013.51+012337.5**: Unclassified - Others: Classified as AGN **Note**: This classifier is based solely on photometric data and does not include a BBH merger category ## Related Work ### Formation Channels for Gravitational Wave Sources **Hierarchical Mergers**: - Li et al. (2025): Component sources of GW231123 originate from successive mergers of ~6 and ~4 first-generation black holes - Delfavero et al. (2025): Best fit to merger of fourth-generation and third-generation black holes - Li & Fan (2025a): 2G+2G merger scenario preferred over 2G+1G or 3G+2G in Bayesian framework **Alternative Scenarios**: - Isolated binary evolution of extremely massive, rapidly rotating helium cores (Croon et al. 2025; Gottlieb et al. 2025) - Chemically homogeneous evolution (Stegmann et al. 2025; Popa & de Mink 2025) - Sustained accretion (Bartos & Haiman 2025) - Strong gravitational lensing (Shan et al. 2025) - Third-generation star formation (Tanikawa et al. 2025; Liu et al. 2025; Wu et al. 2025) - Primordial black holes (Yuan et al. 2025; Luca et al. 2025; Nojiri & Odintsov 2025) ### BBH Mergers in AGN Disks **Environmental Advantages** (McKernan et al. 2012; Tagawa et al. 2020): - Dense gas promotes binary formation and evolution (Li et al. 2023; Dittmann et al. 2024) - Deep potential well retains merger remnants (Varma et al. 2022) - Gas friction drags tilted remnants back into disk (Fabj et al. 2020) **Electromagnetic Emission Mechanisms** (Bartos et al. 2017): 1. **Ram Pressure Stripping**: Stripping of bound gas remnants (McKernan et al. 2019) 2. **Super-Eddington Accretion**: Driving powerful outburst flows (Wang et al. 2021; Kimura et al. 2021; Rodrı́guez-Ramı́rez et al. 2025) 3. **Jet Emission**: Jets driven by accreting remnants (Tagawa et al. 2024; Rodrı́guez-Ramı́rez et al. 2024; Chen & Dai 2024) ### Historical Searches **Previous Candidates** (Graham et al. 2020, 2023; Cabrera et al. 2024): - **ZTF19abanrhr** associated with GW190521: Most famous case, but interpretations remain controversial (Ashton et al. 2021; Palmese et al. 2021; Morton et al. 2023) - He et al. (2025a) reassessment: Only 2 of 7 candidates remained statistically consistent three years later ### Scientific Value **Multi-Messenger Astronomy** (Vajpeyi et al. 2022): - Reveals origin and evolution of stellar-mass black holes in AGN disks - Provides insights into physical processes in AGN environments **Standard Sirens** (Mukherjee et al. 2020; Chen et al. 2022; Li & Fan 2025b): - Independent measurement of Hubble constant - Cosmological parameter constraints ## Conclusions and Discussion ### Main Conclusions 1. **Successful Identification**: Systematic search identifies six AGN flare candidates spatially and temporally coinciding with GW231123 2. **Strong Statistical Support**: Most candidates show strong association evidence (lnO > 5) in Bayesian framework, primarily driven by GW231123's large primary mass 3. **Physical Consistency**: AGN disk parameters (ρ and H) derived through simple Bondi accretion model are consistent with theoretical predictions 4. **Undetermined Nature**: Due to lack of contemporaneous spectroscopic data, the physical origin of these flares remains uncertain 5. **At Most One Real**: At most one of the six candidates is a genuine electromagnetic counterpart to GW231123; the others may be unrelated background events ### Limitations **Observational Constraints**: 1. **No Contemporaneous Spectroscopy**: Cannot confirm BBH origin through features like asymmetric broad line profiles 2. **Single Observation**: No detection of expected recurrent flares 3. **Waveform Uncertainty**: Significant differences in distance posteriors between models, possibly stemming from uncertainties in current signal models at high spin **Methodological Limitations**: 1. **Simplified Model**: Bondi accretion model assumes recoil perpendicular to disk; actual situation may be more complex 2. **False Positive Risk**: TDEs, supernovae, or kilonovae could produce similar flares in AGN disks 3. **Classification Uncertainty**: ALeRCE classifier categorizes J203336.84-275303.9 as supernova, requiring further verification **Statistical Limitations**: 1. **Small Sample**: Only 6 candidates, limited statistical significance 2. **Background Contamination**: Only 6 of 166 initial candidates pass strict filtering, indicating high background rate ### Future Directions **Key Follow-up Observations**: 1. **Continuous Optical Monitoring**: - Using Wide Field Survey Telescope and Vera C. Rubin Observatory - Detection of expected recurrent flares - Timescale: $\sim 1.6 \text{ yr} (M_{\text{SMBH}}/10^8 M_\odot)(r/10^3 r_g)^{3/2}$ - Recurrence interval measurements constrain merger location within AGN disk 2. **Rapid Multi-Wavelength and Spectroscopic Follow-up**: - Timely spectroscopic observations during flare phase - Identification of asymmetric broad line features (BBH merger evidence) - Exclusion of typical spectroscopic features of supernovae or TDEs 3. **High-Energy Observations**: - Fermi Gamma-ray Space Telescope - Einstein Probe - X-ray and γ-ray monitoring constraining physical environment - Non-detections also provide valuable information **Theoretical Development**: 1. Improved waveform models for high-spin BBH 2. More refined simulations of AGN disk-BBH interactions 3. Multi-dimensional radiative transfer models **Long-term Goals**: - Establish statistical sample of electromagnetic counterparts to BBH mergers - Verify AGN-assisted merger scenario - Usher in new era of multi-messenger astronomy ## In-Depth Evaluation ### Strengths **Scientific Significance**: 1. **High Timeliness**: First systematic search for electromagnetic counterparts to the latest and most extreme gravitational wave event (GW231123) 2. **Clear Objectives**: GW231123's supermassive masses and high spin make it an ideal target for electromagnetic counterpart searches 3. **Multi-Messenger Value**: Confirmed counterparts would be important milestones in multi-messenger astronomy **Methodological Advantages**: 1. **Data Quality**: Forced photometry reduces host galaxy contamination; time baseline extended to October 2025 2. **Rigorous Filtering**: Multi-layer filtering criteria (time delay, isolated outburst, timescale constraints) effectively reduce false positives 3. **Bayesian Assessment**: Quantitative evaluation of association probability provides statistical support 4. **Physical Modeling**: Bondi accretion model connects observations to theory, deriving physical parameters **Comprehensive Analysis**: 1. **Multiple Waveforms**: Considers 5 different waveform models, addressing systematic uncertainties 2. **Spectroscopic Verification**: Utilizes SDSS and DESI archival spectra for SMBH mass measurements 3. **Color Evolution**: Analyzes g-r color variations providing additional physical insights **Clear Writing**: 1. Logical structure with clear reasoning 2. Detailed method descriptions ensuring reproducibility 3. Information-rich figures (light curves, sky maps, spectra) ### Weaknesses **Observational Limitations**: 1. **Lack of Confirmation**: No contemporaneous spectroscopy; cannot definitively exclude other transient types 2. **Single Observation**: No detection of recurrent flares; cannot verify key predictions of AGN disk scenario 3. **Missing Multi-Wavelength**: Only optical data available; lacks X-ray, radio, and other multi-wavelength constraints **Methodological Issues**: 1. **Model Simplification**: - Assumes recoil perpendicular to disk (fastest escape path), potentially unrealistic - Ignores disk inhomogeneity and dynamical complexity - Constant color assumption fails for some candidates 2. **Parameter Uncertainty**: - 3 sources use fiducial SMBH mass $10^8 M_\odot$, introducing systematic errors - Optical depth assumption $\tau_{\text{mp}} \sim 10^{4.5}$ lacks verification 3. **Classification Contradiction**: J203336.84-275303.9 classified as supernova by ALeRCE but retained as candidate **Statistical Issues**: 1. **Small Sample**: Only 6 candidates, limiting statistical significance testing 2. **High Background**: 96% of initial 166 candidates excluded, indicating severe background contamination 3. **Multiple Comparisons**: Odds ratios calculated for 5 waveform models without multiple comparison correction **Insufficient Discussion**: 1. Insufficient discussion of J203336.84-275303.9's supernova classification 2. Lack of relative credibility ranking among different candidates 3. No quantification of observing time/resources needed to confirm any candidate ### Impact **Contributions to the Field**: 1. **Methodology**: Establishes systematic workflow for electromagnetic counterpart searches for GW events 2. **Target List**: Provides six specific targets for follow-up observations 3. **Parameter Space**: Extends searched BBH mass and spin parameter space **Practical Value**: 1. **Observational Guidance**: Provides specific targets for future multi-wavelength and spectroscopic follow-up 2. **Theory Testing**: Offers opportunity to test AGN disk BBH merger scenario 3. **Method Reusability**: Search workflow applicable to future gravitational wave events **Potential Impact**: - **If Confirmed**: Major breakthrough in multi-messenger astronomy, verifying AGN formation channel - **If Excluded**: Constrains physical models of electromagnetic emission in AGN disks - **Method Legacy**: AGNFCC catalog and search workflow enable continued application **Reproducibility**: 1. **Open Data**: ZTF data publicly accessible 2. **Available Code**: Uses public software packages (ligo.skymap, PyQSOFit) 3. **Detailed Methods**: Search workflow and filtering criteria clearly described 4. **Shared Catalog**: AGNFCC catalog publicly available on GitHub ### Applicable Scenarios **Ideal Applications**: 1. **Extreme-Mass BBH**: Systems similar to GW231123 with supermassive masses and high spins 2. **Nearby Events**: Luminosity distance <2000 Mpc ensuring optical detectability 3. **Good Localization**: 90% confidence region <1000 deg², facilitating cross-matching 4. **AGN Environment**: Systems expected to form in AGN disks **Extended Applications**: 1. **Intermediate-Mass BBH**: 50-150 M⊙ systems 2. **Neutron Star-Black Hole Mergers**: If occurring in AGN disks 3. **Historical Events**: Application to previously detected gravitational wave events **Inapplicable Scenarios**: 1. **Low-Mass Systems**: <30 M⊙ BBH with potentially undetectable electromagnetic emission 2. **Extreme Distances**: >5000 Mpc with optical signals below detection threshold 3. **Poor Localization**: >10000 deg² with low cross-matching efficiency 4. **Field Environment**: Mergers in non-AGN environments with no expected strong electromagnetic emission ## References ### Key Gravitational Wave Papers - The LIGO Scientific Collaboration et al. (2025): Discovery and properties of GW231123 - Abbott et al. (2020): Early BBH detections and mass gap discussion ### AGN Disk BBH Theory - McKernan et al. (2012, 2019): BBH formation in AGN disks and electromagnetic emission - Bartos et al. (2017): Theoretical expectations for electromagnetic counterparts in AGN environments - Wang et al. (2021): Super-Eddington accretion and outflow models ### Electromagnetic Counterpart Searches - Graham et al. (2020, 2023): ZTF19abanrhr and GW190521 association - He et al. (2025a): Reassessment of historical candidates - He et al. (2025b): AGNFCC catalog construction ### Hierarchical Merger Studies - Gerosa & Fishbach (2021): Spin distribution of hierarchical merger black holes - Li et al. (2025): Hierarchical merger interpretation of GW231123 ### AGN Disk Structure - Sirko & Goodman (2003): AGN disk theoretical models - Thompson et al. (2005): Disk vertical structure and density distribution --- **Overall Assessment**: This is a timely, systematic, and methodologically rigorous observational paper providing the first electromagnetic counterpart search for GW231123, an extreme gravitational wave event. Despite observational and methodological limitations, the discovery of six candidates provides clear targets for future follow-up observations. Confirmation of any candidate would represent a major breakthrough in multi-messenger astronomy, verifying theoretical expectations for AGN-assisted BBH mergers. The paper's methodological contributions and reproducibility ensure continued value to the field.