We report the serendipitous multiwavelength discovery of a candidate dual black hole system with a separation of 100 pc in the gasrich luminous infrared galaxy MCG033464 z  0.016. Hubble Space TelescopeAdvanced Camera for Surveys observations show three distinct optical centroids in the O III narrowband and F814W images. Subsequent analysis of ChandraACIS data shows two spatially resolved peaks of equal intensity in the neutral Fe K 6.26.6 keV band while highresolution radio continuum observations with the Very Large Array at 8.46 GHz 3.6 cm band show two spatially coincident radio peaks. Fast shocks as the ionizing source seem unlikely given the energies required for the production of Fe K. If confirmed the separation of 100 pc would represent the closest dual active galactic nuclei reported to date with spatially resolved multiwavelength observations.
['AGN host galaxies', 'Seyfert galaxies', 'High energy astrophysics', '2017', '1447', '739', 'Astrophysics - Astrophysics of Galaxies']
Resolving a Candidate Dual Active Galactic Nucleus with 100 pc Separation in MCG033464
{'Resolving a Candidate Dual Active Galactic Nucleus with ∼ 100pc Separation in MCG03-34-64': '<!-- image -->  \nAnna Trindade Falcão 1 , T. J. Turner 2 , S. B. Kraemer 3 , J. Reeves 3,4 , V. Braito 3,4,5 , H. R. Schmitt 6 , and L. Feuillet 3 1 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA 2 Eureka Scienti /uniFB01 c, Inc., 2452 Delmer St., Suite 100, Oakland, CA 94602, USA 3 Institute for Astrophysics and Computational Sciences, The Catholic University of America, Washington, DC 20064, USA 4 INAF - Osservatorio Astronomico di Brera, Via Bianchi 46, 23807, Merate ( LC ) , Italy 5 Dipartimento di Fisica, Università di Trento, Via Sommarive 14, Trento 38123, Italy 6  \nNaval Research Laboratory, Washington, DC 20375, USA  \nReceived 2024 May 24; revised 2024 July 20; accepted 2024 August 4; published 2024 September 9', 'Abstract': 'We report the serendipitous multiwavelength discovery of a candidate dual black hole system with a separation of ∼ 100 pc, in the gas-rich luminous infrared galaxy MCG-03-34-64 ( z = 0.016 ) . Hubble Space Telescope / Advanced Camera for Surveys observations show three distinct optical centroids in the [ OIII ] narrow-band and F814W images. Subsequent analysis of Chandra / ACIS data shows two spatially resolved peaks of equal intensity in the neutral Fe K α ( 6.2 -6.6 keV ) band, while high-resolution radio continuum observations with the Very Large Array at 8.46 GHz ( 3.6 cm band ) show two spatially coincident radio peaks. Fast shocks as the ionizing source seem unlikely, given the energies required for the production of Fe K α . If con /uniFB01 rmed, the separation of ∼ 100 pc would represent the closest dual active galactic nuclei reported to date with spatially resolved, multiwavelength observations.  \nUni /uniFB01 ed Astronomy Thesaurus concepts: AGN host galaxies ( 2017 ) ; Seyfert galaxies ( 1447 ) ; High energy astrophysics ( 739 )', '1. Introduction': "The masses of supermassive black holes ( SMBHs ) in active galactic nuclei ( AGNs ) correlate with the global properties of their host galaxies ' stellar components, such as luminosity, mass, and velocity dispersion, extending over kiloparsec scales ( e.g., Kormendy & Ho 2013 ) . This correlation highlights the need to understand the mechanisms driving SMBH growth.  \nBoth galactic evolutionary models and observations suggest that a signi /uniFB01 cant fraction of AGNs, particularly those at the center of large-scale structures, undergo major mergers ( e.g., De Lucia & Blaizot 2007; Hopkins et al. 2008; Ginol /uniFB01 et al. 2017; Castignani et al. 2020 ) . Hydrodynamical simulations further demonstrate that major mergers induce gas in /uniFB02 ows toward galactic centers, potentially triggering both star formation and accretion onto central SMBHs ( Mayer et al. 2007 ) . However, the overall impact of these events on SMBH growth throughout cosmic time remains poorly constrained.  \nSMBH pairs, often manifested as dual AGNs, provide distinctive evidence for merger-fueled SMBH growth ( e.g., Wassenhove et al. 2012 ) . Numerous dual AGN candidates have been identi /uniFB01 ed using various techniques, including optical spectroscopy with emission line ratios ( e.g., Liu et al. 2011 ) , hard X-ray emission ( e.g., Koss et al. 2011 ) , and doublepeaked narrow emission lines ( e.g., Smith et al. 2010; Koss et al. 2023 ) . Nonetheless, these methods have limitations, and multiwavelength follow-up observations have revealed a substantial number of false positives ( e.g., Fu et al. 2011b ) .  \nThe advent of gravitational-wave astronomy, with the potential for detection through pulsar timing arrays ( e.g.,  \n<!-- image -->  \nOriginal content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author ( s ) and the title of the work, journal citation and DOI.  \nVerbiest et al. 2016 ) , has heightened the importance of understanding the formation timescales of binary systems. Studying kiloparsec and subkiloparsec dual AGNs offers a unique window into the /uniFB01 nal stages of SMBH binary coalescence, a crucial process in gravitational wave astronomy.  \nDual AGNs separated by kiloparsec or subkiloparsec scales are inherently more challenging to detect and investigate than wider-separation systems ( e.g., > 3 kpc ) . This dif /uniFB01 culty arises from increased obscuration in late-stage mergers ( e.g., Koss et al. 2016; Ricci et al. 2021; De Rosa et al. 2022 ) , limitations in telescope spatial resolution ( particularly at subkiloparsec scales ) , the scarcity of detected radio-bright dual systems ( Burke-Spolaor 2011 ) , and the limitations of optical selection using double-peaked narrow emission lines ( prone to false positives; see Fu et al. 2011a ) . Existing observations of dual AGNs tentatively suggest that AGN triggering becomes more prevalent in advanced mergers with stellar bulge separations < 10 kpc ( e.g., Koss et al. 2010; Fu et al. 2018; Stemo et al. 2021 ) , aligning with simulations of SMBH accretion and evolution in such mergers ( e.g., Blecha et al. 2018 ) . Therefore, studying nearby galaxies hosting dual AGNs separated at subkiloparsec scales is crucial for advancing our understanding of the late stages of galaxy mergers, the triggering and fueling of AGN activity, and the dynamics of SMBH pairs ( Steinborn et al. 2016 ) . These close-separation systems provide a unique window into the processes leading to the eventual coalescence of SMBHs, which is a major source of gravitational waves, and plays a fundamental role in the growth of SMBHs and their host galaxies ( Dotti et al. 2012; Kharb et al. 2017 ) .  \nWhile several dual AGN candidates have been proposed at scales of hundreds of parsecs, often supported by single-wave band observations, these have frequently been challenged by subsequent studies. Notable examples include the nearby Seyfert NGC 3393 ( Fabbiano et al. 2011 ) , SDSS J101022.95 + 141300.9  \nTable 1 Multiwavelength Observations of MCG-03-34-64", 'Note.': '( Goulding et al. 2019 ) , and a third active nucleus in NGC 6240 ( Kollatschny et al. 2020 ) , later disputed in other works ( Koss et al. 2015; Veres et al. 2021; Treister et al. 2020 ) .  \nRecently, Koss et al. ( 2023 ) have used multiwavelength observations to con /uniFB01 rm the presence of two active nuclei separated by ∼ 230 pc in the center of UGC 4211, a late-stage major galaxy merger. Their discovery builds on the detection of ( i ) near-infrared ( NIR ) broad lines associated with one of the nuclei; ( ii ) two separated sources of [ OIII ] emission, unresolved at ∼ 0 1 resolution; ( iii ) unresolved and compact ( < 30 pc ) nuclear millimeter emission at the location of the two nuclei, and coincident with the unresolved emission-line sources; and ( iv ) molecular gas velocity gradient centered on the position of both nuclei.  \nIn this study, we present the serendipitous discovery of a candidate dual AGN system in MCG-03-34-64 ( IRAS 131971627 ) , a nearby early-type infrared luminous galaxy at z = 0.01654 ( ∼ 78 Mpc, from NASA / IPAC Extragalactic Database ) . 7 This galaxy is identi /uniFB01 ed as one of the hardest X-ray sources in the local Universe ( Tatum et al. 2016 ) . Earlier X-ray observations with ASCA, XMM-Newton, and BeppoSAX ( Dadina & Cappi 2004; Miniutti et al. 2007 ) revealed an extremely hard and complex source spectrum, attributed to heavy absorption from a multilayered and clumpy medium. MCG-03-34-64 also shows extended radio emission ( ∼ 300 pc ) , roughly aligned with the major axis of the host galaxy ( Schmitt et al. 2001 ) , and ∼ 2 " extent in mid-infrared aligned in the same direction as the radio structure ( Hönig et al. 2010 ) .  \nWe have obtained Hubble Space Telescope / Advanced Camera for Surveys ( ACS ) imaging of MCG-03-34-64 in 2022 June ( P.I.: Turner, proposal ID: 16847 ) , and 50 ks of Chandra / ACIS-S observations in 2023 April ( obs ids 25253, 27802, and 27803, P.I.: Turner ) . This paper presents the results of the analysis of these new data sets, combined with existing Very Large Array ( VLA ) radio, and Hubble Space Telescope ( HST ) optical imaging of the source. Throughout this paper, we adopt /uni03A9 m = 0.3, /uni03A9 Λ = 0.7, and H 0 = 70 km s -1 Mpc -1 , and a scale of 340 pc arcsec -1 , based on the redshift at the galaxy \' s distance.', '2. Multiwavelength Observations and Analysis': 'Table 1 lists all observations used in this paper, including instruments, /uniFB01 lters, observation dates, obs ids, and exposure  \ntimes. Details on the reduction of new HST / ACS and Chandra / ACIS-S observations are provided in Sections 2.1 and 2.2, respectively. For reduction and analysis of archival HST F814W, VLA 8.46 GHz imaging, and Suzaku and XMMNewton spectroscopy, see Section 2.3. All the HST data used in this paper can be found in MAST at doi:10.17909 / 53rj-fw34.', '2.1. Hubble Space Telescope Imaging': 'HST / ACS observations of MCG-03-34-64 were obtained using the linear ramp /uniFB01 lter FR505N ( narrow-band [ OIII ] ) centered at 5089.6 Å , to characterize the morphology of the emission-line gas, while a continuum medium band image was obtained using FR647M, centered at 5590 Å . These /uniFB01 lters have bandwidths of 2% and 9%, respectively. Standard HST pipeline procedures were employed for data reduction. The narrow-band and continuum images were acquired sequentially and did not require realignment. Flux calibration was performed using information available on the headers.', '2.2. Chandra Imaging and Spectroscopy.': 'MCG-03-34-64 was observed with Chandra / ACIS for a total exposure of 50 ks, split into three observations ( Table 1 ) . We use CIAO 4.15 8 for data processing and analysis, merging individual observations to create a deeper data set for imaging analysis. Each image in the 0.3 -7 keV band was /uniFB01 rst visually inspected, and the required centroid shifts were found to be less than one ACIS instrument pixel ( 0 492 ) . We created images in the 2 -7 keV band, and using obs id 27802 ( longest observation, see Table 1 ) as the astrometric reference, shifts of [ /uni0394 x , /uni0394 y ] = [ -0.7, -0.8 ] , and [ 0.8,0.3 ] pixels were applied to obs ids 25253 and 27803, respectively. All data sets were processed to enable subpixel analysis ( e.g., Wang et al. 2011a ) .  \nSubpixel imaging binning was employed to effectively oversample the Chandra point-spread function ( PSF ) and overcome the limitations of the ACIS instrumental pixel size. This method has been extensively used and validated in previous studies examining the subkiloparsec regions around nearby and obscured AGNs ( e.g., Maksym et al. 2017; Fabbiano et al. 2018a; Ma et al. 2021; Trindade Falcão 2023 ) , demonstrating excellent agreement between reconstructed ACIS-S features and those imaged with higher spatial  \nFigure 1. HST [ O III ] and F814W images of the central region of MCG-03-3464. Note the prominent diffraction spikes present in the [ O III ] image.  \n<!-- image -->  \nresolution instruments such as HST and VLA ( e.g., Wang et al. 2011b; Paggi et al. 2012; Maksym et al. 2019; Fabbiano et al. 2018b ) . The Chandra PSF was simulated using ChaRT 9 and MARX . 10 This work uses a /uniFB01 nal Chandra scale of one-eighth of the native ACIS pixel.  \nFor the spectral analysis ( Section 3.2 ) , we used CIAO specextract 11 to extract the spectra from a 2 5 circular region, while the background was subtracted from a circular region of 3 5 radius. Individual spectra and response /uniFB01 les were coadded with mathpha 12 and the standard addrmf 13 and addarf 14 functions.', '2.3. Archival Radio / Optical / X-Ray Observations': 'In addition to the new Chandra and HST data sets, we analyze archival optical, radio, and X-ray observations of MCG-03-34-64, as listed in Table 1. These data include 8.46 GHz radio imaging with VLA, optical continuum imaging with HST / ACS F814W, and X-ray spectra from Suzaku and XMM-Newton.  \nThere are four additional archival Chandra observations with MCG-03-34-64 in the /uniFB01 eld of view ( obs ids 27267, 27786, 7373, and 23690 ) . However, three of these observations are not usable due to the galaxy being located at the very edge of the /uniFB01 eld ( observations were optimized for the companion galaxy ) . The fourth available Chandra observation consists of a 7 ks snapshot ( used in Miniutti et al. 2007 ) , which has insuf /uniFB01 cient counts for meaningful imaging analysis.', '3. Results': '3.1. Imaging Analysis', '3.1.1. Hubble Space Telescope Imaging': "The [ OIII ] narrow-line region ( NLR ) in MCG-03-34-64 has a highly unusual morphology, featuring three distinct, and compact emission regions, as shown in Figure 1. The NLR extends ∼ 2.3 kpc along the NE -SW direction. In the perpendicular direction ( NW -SE ) , we observe three diffraction spikes characteristic of point sources, while one diffraction spike is observed along the NE -SW cone. These features suggest high concentrations of [ OIII ] gas within a relatively small region, a rare occurrence in the local Universe ( Fischer et al. 2018 ) .  \nTo examine the overall structure of the NLR, we model the [ OIII ] light distribution with GALFIT ( Peng et al. 2002, 2010 ) . The modeling was performed with generic Sérsic pro /uniFB01 les, and we allow all parameters to freely vary during the /uniFB01 tting process. The background level and standard deviation were determined from blank regions within the image ' s /uniFB01 eld of view.  \nThe best/uniFB01 t GALFIT model consists of four Sérsic components, one for each of the three peaks visible in the [ OIII ] surface plot in the leftmost panel of Figure 2, and an additional larger component to take into account the underlying fainter, more extended emission. In the second panel, we show a zoomed-in image of the [ OIII ] emission, and mark the position of each Sérsic component corresponding to a strong peak of [ OIII ] emission as green circles. The third panel of Figure 2 shows the best/uniFB01 t GALFIT model, while the fourth panel shows the residual images from this best/uniFB01 t model, with the different model components subtracted from the data.  \nTable 2 lists the best/uniFB01 t model parameters. The strong [ OIII ] peaks of emission are /uniFB01 t with Sérsic components with indices 0.41 /uni0084 n /uni0084 0.51, indicating that they are similar to Gaussians ( n = 0.5 ) , but in two cases show a slightly more centrally concentrated distribution. Given that these components have effective radii ∼ 2 -3 times that of an unresolved source, we opted to not include a PSF component in the /uniFB01 tting model, since it is not expected to signi /uniFB01 cantly change the results.  \nThe diffraction spikes associated with these regions ( Figure 1 ) suggest angular sizes on the order of HST / ACS ' s resolution element, 0 1 ( 2 pixels ) . Assuming the measured effective radius, Re , of each [ OIII ] component in Table 2 represents its intrinsic size, we estimate that each [ OIII ] -emitting region has an average angular diameter of ∼ 0 2, or 60 pc. These regions are separated by 79 ± 8pc ( northern to central ) and 76 ± 8pc ( central to southern ) , and are also clearly observed in the broad continuum F814W HST image ( Figure 1, right panel ) , a point further discussed in Section 4.  \nGiven the unique morphology of the optical emission observed in MCG-03-34-64 with HST, we examine other available multiwavelength observations ( Table 1 ) to obtain a more comprehensive picture of this inner region.", '3.1.2. Chandra Imaging': 'Figure 3 shows the inner ∼ 200 pc region of MCG-03-34-64, as observed with Chandra / ACIS in different energy bands. The images are binned at one-eighth of the native pixel to use the full-resolution of the instrument in the high-count inner region, and processed with 1 kernel Gaussian smoothing.  \nThe X-ray-emitting region shows a compact morphology, with the most prominent features observed within the inner ∼ 0 5 ( ∼ 170 pc ) circumnuclear region. Extended diffuse emission is detected in all bands, as shown in Figure 3; however, there are clear energy-dependent differences in the surface brightness morphology. At higher energies ( Figure 3, two rightmost panels ) , the emission from the nuclear region cannot be explained as a single point source. Instead, in the hard neutral Fe K α ( 6.2 -6.6 keV ) and 6 -7 keV bands, we observe two spatially resolved peaks of emission separated by 125 ± 21 pc.  \nWe measure nearly equal X-ray luminosities in the narrow 6.2 -6.6 keV band from the Chandra image for both Fe K α peaks, L ( 6.2 -6.6 keV ) ∼ 3.2 ± 0.6 × 10 40 erg s -1 ,for D ∼ 78 Mpc ( e.g., de Grijp et al. 1992; Table 3 ) .  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 2. First panel: surface plot of the [ O III ] /uniFB02 ux in the HST / ACS image where three distinct sources are visible. Second panel: HST / ACS [ O III ] image showing the presence of three closely separated emission peaks. We mark the position of each Sérsic component for the emission peaks as green circles. Third panel: best/uniFB01 t GALFIT model, using four Sérsic components. Fourth panel: residual image after subtraction of the best/uniFB01 t GALFIT model.  \n<!-- image -->  \nFigure 3. Chandra / ACIS-S merged data set showing the inner 200 pc region of MCG-03-34-64 in different bands ( one-eighth subpixel, and 1 kernel Gaussian smoothing ) . First panel: soft ( 0.3 -3 keV ) X-ray image. Second panel: 3 -5 keV hard continuum image. Third panel: Fe K α ( 6.2 -6.6 keV ) Chandra image. We show the location of the two Fe K α centroids as blue circles. Fourth panel: 6 -7 keV hard band Chandra image.  \n<!-- image -->  \nTable 2 Best/uniFB01 t GALFIT Parameters ( c n 2 = 5.411 )Table 3 Chandra / ACIS-S Merged Data Set -Astrometry, Fluxes, and Luminosities of the Individual Fe K α', 'Regions': 'Note. The luminosities are calculated for D = 78 Mpc.  \nTo address concerns that the dual morphology might be a spurious detection due to smoothing on scales smaller than the Chandra PSF, we examine the individual Chandra observations ( prior to merging ) listed in Table 1, in the Fe K α band ( Figure 4 ) . These images are binned at one-eighth of the native pixel and smoothed with a 1 kernel Gaussian. As shown, the dual morphology of the Fe K α band is indeed observed in all individual observations, prior to merging, con /uniFB01 rming the robustness of the detection. The differences in surface brightness between the individual Chandra exposures seen in Figure 4 for the individual nuclei are most likely due to statistical noise.', '3.1.3. VLA Imaging': 'The morphology of the 8.46 GHz radio continuum emission in MCG-03-34-64 is also analyzed, and shown in Figure 5.The positions of the two radio centroids identi /uniFB01 ed by Schmitt et al. ( 2001 ) in the inner region are shown as white circles ( see their Table 2 ) . The 8.46 GHz emission starts as a linear structure at the position of the northern radio centroid, extending ∼ 100 pc southwestward to the central radio peak, and then bending southward in the direction of the southern [ OIII ] centroid ( Figure 5 ) .Table4 lists the positions and /uniFB02 uxes of the individual  \nFigure 4. Chandra / ACIS individual exposures of MCG-03-34-64 in the Fe K α ( 6.2 -6.6 keV ) band, binned at one-eighth of the native pixel and smoothed with a 1 kernel Gaussian. As shown, the dual morphology is observed prior to merging. The distances between the two peaks in the individual observations are 0 486, 0 408, and 0 443, for obs ids 25253, 27802, and 27803, respectively. Given the 0 08 uncertainty in the Chandra data, these are all consistent within the error.  \n<!-- image -->  \nFigure 5. VLA-A 8.46 GHz ( 3.6 cm ) radio continuum image of MCG-03-3464. White circles mark the position of the two radio peaks, and the [ O III ] centroids are shown in green. The image is shown in log scale.  \n<!-- image -->  \nTable 4  \nVLA-A 3.6 cm Radio Continuum Image Decomposition -Position and Fluxes of Individual Components  \nNote. From Schmitt et al. ( 2001 ) .  \nradio components, as measured in Schmitt et al. ( 2001 ) .The separation between the two centroids is 116 ± 14 pc.', '3.1.4. Astrometry Registration': "We initially apply CIAO wavdetect 15 to the merged 0.3 -7 keV Chandra image with a > 5 σ detection threshold, detecting two faint sources near the edge of the chip array ( via comparison with the Vizier source catalog ) . However, given their faintness, these are not suitable to use as a basis for astrometry correction.  \nWe then create Chandra images in the narrow Fe K α band ( 6.2 -6.6 keV ) , known to be dominated by nuclear emission in obscured sources ( Fabbiano & Elvis 2024 ) . Assuming that the radio emission in AGNs also originates from the innermost regions around the SMBH, we correct Chandra ' s absolute astrometry by aligning the emission peaks seen in the Fe K α band with those in the 8.46 GHz radio image. This alignment method has been used in similar studies, such as in Mrk 78 ( Fornasini et al. 2022 ) .  \nWe apply a total shift of [ /uni0394 x, /uni0394 y ] = [ 0.1, 0.7 ] pixels to the Chandra / ACIS data, well within the absolute astrometry accuracy of the telescope. 16 Comparison with available HST observations con /uniFB01 rms the accuracy of VLA ' s astrometry.", '3.2. Spectroscopic Analysis': 'We also analyzed Chandra / ACIS, NuSTAR, and XMMNewton spectroscopic observations of MCG-03-34-64, spanning 14 yr ( Table 1 ) . In Figure 6 ( left panel ) , we show the combined 50 ks Chandra / ACIS spectrum of the source, which was extracted from a 2 5 circular region. This spectrum is heavily absorbed at the lower energies, below 5 keV, indicating a heavily obscured nuclear source. At higher energies, we observe a broad neutral Fe K α line, with E = 6.39 ± 0.01 keV ( the 6.2 -6.6 keV Fe K α band is shown as red lines ) , equivalent width EW = 470 ± 60 eV, and σ = 0.11 ± 0.05 keV.  \nThe overall spectral pro /uniFB01 le observed in the Chandra data is consistent across the earlier NuSTAR and XMM-Newton observations. All three spectra show an absorption feature at 6.8 keV, likely arising from Fe XXV absorption, and suggesting an out /uniFB02 ow with v ∼ 5000 km s -1 , as noted in Miniutti et al. ( 2007 ) . Below 3 keV, the soft X-ray emission is dominated by photoionized and collisionally ionized emission lines, from Ne, O, and Fe L ions ( Miniutti et al. 2007 ) .', '3.2.1. XMM-Newton and NuSTAR Spectral Fitting': "We proceed to model the X-ray spectrum of MCG-03-34-64 with MYTorus ( Murphy & Yaqoob 2009 ) , a physically motivated model built to describe the interaction of the emission from an X-ray point-source with a surrounding, and homogeneous torus of cold neutral material.  \nWe /uniFB01 t the joint XMM + NuSTAR X-ray spectrum with a source model of the form:  \nA × TBabs × [ xstar × MYTZ × zpowerlw + ( C × ( MYTS + MYTL × gsmooth ) + zpowerlaw \\_ soft + soft \\_ emiss ) ] , where TBabs describes the absorption of emission by the Galactic column density, xstar is the photoionized absorber described previously in Miniutti et al. ( 2007 ) , [ MYTZ × zpowerlw] describes the intrinsic continuum in transmission absorbed by torus, MYTS is the scattered ( re /uniFB02 ected ) toroidal component off Compton-thick matter, MYTL is the associated Fe / Ni K α / K β line emission, and gsmooth accounts for some Gaussian broadening of the MYTorus line emission, where the upper limit on the line width is σ < 65 eV. The soft X-ray components are zpowerlaw \\_ soft , which is an unabsorbed scattered power-law component, and soft \\_ emiss , which is the sum of the photo and collisionally ionized emission components described by Miniutti et al. ( 2007 ) . A is the cross normalization factor between NuSTAR and XMM ( A = 1.20 ± 0.05 ) and C is the offset between re /uniFB02 ected / line components and intrinsic continuum components  \nFigure 6. Left: Chandra / ACIS spectrum from the inner 2 5 region of MCG-03-34-64, with the y -axis giving the /uniFB02 uxes in ν F ν units, and the energy of the Fe K α line E = 6.39 ± 0.01. Note the heavy absorption of the soft photons at the lower energies, which is attributed to heavy obscuration of the source. The red dashed lines mark the 6.2 -6.6 keV Fe K α band. Center: joint XMM-Newton & NuSTAR spectrum, and MYTorus individual model components. The black curve is the total emission, blue is the transmitted continuum absorbed by Compton-thin matter along the line of sight, magenta is the scattered ( re /uniFB02 ected ) MYTorus component off Compton thick matter, cyan is the Fe K line emission, red is the ionized soft X-ray emission, and green is an unabsorbed power-law component. Right: MYTorus model /uniFB01 tto the joint XMM-Newton and NuSTAR observation. The model components have the same colors as in the center panel. The data and model components are folded through the instrumental responses.  \n<!-- image -->  \nFigure 7. Chandra / ACIS-S merged Fe K ( 6.2 -6.6 keV ) image of MCG-03-3464 ( one-eighth subpixel, and smoothed with 1 kernel Gaussian ) . Optical centroids from HST / ACS are shown in red, VLA-A 8.46 GHz centroids in white, and Chandra / ACIS Fe K centroids in blue. The circle sizes re /uniFB02 ect uncertainties in the position of the centroids.  \n<!-- image -->  \n( C is frozen at 1 ) . The results of the /uniFB01 tting are shown in Figure 6 ( center and right panels ) .  \nThis model returns a good /uniFB01 t, c = n 1.05 2 ,with Γ = 2.12 ± 0.03, =  ' N 5.4 0.3 10 H 23 eq cm -2 ( absorption column density in the equatorial plane, i.e., line-of-sight column density ) , and N H > 6.0 × 10 24 cm -2 ( out of line-of-sight column density ) .In this case, we assume that the hard intrinsic and soft scattered power-law continua have the same photon index.  \nFrom the /uniFB01 tting results, we derive an integrated Fe K α luminosity of L ( 6.2 -6.6 keV ) = 1.0 × 10 41 erg s -1 ( from MYTL ) ,a scattered ( re /uniFB02 ected ) continuum luminosity of = -L 2 10 keV scat. ( ) ' 2.1 10 41 erg s -1 ( from MYTS ) , and a scattering corrected intrinsic 2 -10 keV luminosity of = ' -L 1.5 2 10 keV int. ( ) 10 43 erg s -1 ( from MYTZ ) .  \nGiven the quality and resolution of the X-ray observations used in this work, the spectra and /uniFB01 tting models employed in our analysis account for emission within the entire inner region of MCG-03-34-64, and cannot be performed separately for the  \nindividual Fe K α peaks uncovered in the Chandra imaging data. In this case, the resulting con /uniFB01 guration suggested by MYTorus requires one where one is looking along the edge ( Compton-thin line of sight ) of a very Compton-thick absorber overall, which obscures both Fe K α regions.  \nWe also note that the difference observed between the summed Fe K α luminosities derived from the Chandra imaging, ∼ 6.4 × 10 40 erg s -1 , and the total Fe K α luminosity yielded by the spectral /uniFB01 tting with MYTorus , L ( 6.2 -6.6 keV ) = 1.0 × 10 41 erg s -1 , may be attributed to line absorption by the absorber, which in turn implies a higher intrinsic X-ray luminosity, as yielded by the results of the spectral /uniFB01 t.", '4. Discussion': 'The results presented in Section 3 reveal puzzling properties of the emission in MCG-03-34-64. Our imaging analysis identi /uniFB01 ed three [ OIII ] -emitting regions in the HST / ACS data, separated by 76 ± 8 and 79 ± 8pc ( Table 2, Figure 2 ) . In X-rays with Chandra / ACIS, two spatially resolved peaks of emission are observed in the narrow 6.2 -6.6 keV Fe K α band, separated by 125 ± 21 pc ( Table 3, Figure 3 ) . In the radio with VLA-A, Schmitt et al. ( 2001 ) previously identi /uniFB01 ed two distinct radio cores in the 8.46 GHz continuum, separated by 116 ± 14 pc ( Table 4, Figure 5 ) .  \nFigure 7 shows the Chandra / ACIS Fe K α image and the position of these multiwavelength centroids. The image is binned at one-eighth of the native ACIS-S pixel and smoothed with 1 kernel Gaussian.', '4.1. Bolometric Luminosity': "Using the [ OIII ] /uniFB02 uxes in Table 2,and D ∼ 78 Mpc ( e.g., de Grijp et al. 1992 ) , we calculate individual [ OIII ] luminosities for each emitting region and for the Sérsic components modeling the HST data, as listed in Table 6. We then calculate individual bolometric luminosities by applying a correction factor of c = 454 ( Lamastra et al. 2009 ) to L [ O III ] ( Table 6 ) . The northern [ OIII ] region has the highest [ OIII ] and bolometric luminosities, 3.5 × 10 41 erg s -1 and 1.6 × 10 44 erg s -1 , respectively, while the central region is the least luminous, 1.6 × 10 41 erg s -1 and 7.3 × 10 43 erg s -1 . We calculate an integrated bolometric luminosity of = ' L 4.8 10 bol 44 OIII [ ] erg s -1 by coadding individual values ( Table 6 ) .  \nTable 5 Joint XMM and NuSTAR Spectral Fitting Results from MYTorus  \nNotes. We use a correction factor k = 30 ( Vasudevan & Fabian 2007 ) to obtain the integrated bolometric luminosity in X-rays.  \nIn X-rays, we use the scattering-corrected intrinsic 2 -10 keV luminosity from Table 5 and an X-ray correction factor of k = 30 ( Vasudevan & Fabian 2007 ) to calculate an integrated bolometric luminosity of = ' -L 4.5 10 bol 44 2 10 keV ( ) erg s -1 , consistent with L bol OIII [ ] derived from [ OIII ] ( Table 5 ) .", '4.1.1. HST F814W Continuum Fluxes': "The three distinct regions of [ OIII ] emission observed in the continuum-subtracted image ( Figures 1 and 2 ) are also visible in the F814W wide-band image. We estimate the total ( integrated ) intrinsic /uniFB02 ux in the F814W band, FF W 814 total , assuming a typical spectral energy distribution ( SED ) as in, e.g., Revalski et al. ( 2018 ) , Trindade Falcão et al. ( 2021 ) , and the scatteringcorrected intrinsic 2 -10 keV luminosity from the MYTorus /uniFB01 t, -L 2 10 keV int. ( ) = 1.5 × 10 43 erg s -1 . Our results yield a total /uniFB02 ux in the band of ~ ' -F 1.0 10 F W 814 11 total erg s -1 cm -2 .  \nWealso calculate the individual observed /uniFB02 uxes in the F814W band for each [ OIII ] region ( f F W 814 obs. ) in the HST / ACS image, assuming the regions /uniFB01 tinthe [ OIII ] image ( see Table 2 ) are the same in the F814W band. The F814W image was /uniFB02 ux calibrated by multiplying it by 7.104 × 10 -20 erg cm -2 Å -1 e -1 and the F814W /uniFB01 lter bandwidth ( ∼ 2710 Å ) . For simplicity, we repeat the /uniFB01 tpresentedinTable2 using Gaussians for the three regions of [ OIII ] emission, which resulted in components with similar parameters and /uniFB02 ux differences of less than ∼ 5%. The F814W image was used to measure the peak of emission f F W 814 peak at the central position of each corresponding [ OIII ] component, after subtracting the contribution from the host galaxy, which is estimated based on a linear /uniFB01 t to the outer portion of the onedimensional pro /uniFB01 le cut through each of the peaks of emission. These f F W 814 peak values were converted into f F W 814 obs. using  \n= ' ' ' f b a f 1.1377 FWHM . 1 F W F W 814 2 814 obs. peak ( ) ( )  \nOur results reveal that the integrated observed /uniFB02 uxes are /uni0084 3% of the integrated intrinsic /uniFB02 ux in the band ( Table 6 ) , with the largest fraction originating from the central region. These fractions are consistent with scattered, hidden continuum ( e.g., Pier et al. 1994 ) , but could also include contributions from emission lines ( e.g., Kraemer & Crenshaw 2000 ) and recombination continuum ( Osterbrock & Robertis 1985 ) . Therefore, it is unlikely that we will detect AGN continuum emission directly in the optical, and we cannot determine which of these regions harbors the AGN, given that all three regions are consistent with scatter continua from an active nucleus ( but see below ) .  \nTable 6 [ O III ] Luminosities Calculated from the Measured [ O III ] Fluxes from Table 2, and Considering a Distance of D = 78 Mpc  \nNote We use a correction factor c = 454 ( Lamastra et al. 2009 ) to calculate the bolometric luminosity in each [ O III ] region. We also show the measured F814W /uniFB02 uxes for each emitting region.", '4.2. Multiwavelength Emission Centroids': "The high /uniFB02 uxes and luminosities found in Section 3 for individual emission regions in the optical, X-ray, and radio bands support the presence of an AGN in this system. However, pinpointing the AGN ' s location is more challenging. We discuss possible interpretations for the system ' s con /uniFB01 guration, based on our results and the limitations of the data.", '4.2.1. Single AGN + Shocked Interstellar Medium': 'One interpretation is that the active nucleus is located at the position of the northern centroids ( [ OIII ] ,FeK α , and radio; see Figure 7 ) , based on the /uniFB02 uxes of individual components ( Tables 2, 3, and 4 ) . In this single AGN scenario, the remaining emission centroids ( central Fe K α , [ OIII ] and radio centroids, and southern [ OIII ] ) may arise from the interaction between the AGN and the interstellar medium ( ISM ) . This would manifest as a mix of photoionized and collisionally ionized ( shocked ) gas from an extended NLR, similar to NGC 3393 ( e.g., Maksym et al. 2016 ) . Such an interpretation is consistent with previous spectral /uniFB01 tting results for this galaxy, which indicate a mix of photoionized and shock-ionized gas in the soft X-ray emission ( Miniutti et al. 2007 ) . Similarly, it is possible that the AGN in this system is located at the position of the central [ OIII ] ,FeK α , and radio peaks, while the remaining multiwavelength centroids may be attributed to AGN -ISM shock emission.  \nHowever, the \' single AGN + shocked ISM \' scenario in this galaxy is challenged by the high and nearly equal X-ray luminosities derived for each Fe K α region in the Chandra data. The 6.4 keV Fe K α emission line typically arises from irradiation of the neutral ( or low ionized ) iron by a hard X-ray source, and in AGNs, this line may originate from the innermost part of the accretion disk, close to the central SMBH. This emission process is known to be extremely energetically inef /uniFB01 cient, with an overall conversion factor of ∼ 1% of the intrinsic X-ray luminosity ( e.g., Murphy & Yaqoob 2009 ) . This is illustrated by the /uniFB01 tting results with MYTorus ( Section 3.2, Table 5 ) , in which the total Fe K α luminosity L ( 6.2 -6.6 keV ) is ∼ 0.7% of the total scatteringcorrected intrinsic 2 -10 keV luminosity -L 2 10 keV int. ( ) .  \nIn Section 2.2, the derived Fe K α luminosities for the individual regions are on the order of 3.2 ± 0.6 × 10 40 erg s -1 ( Table 3 ) . This suggests that two individual hard continuum sources are needed to provide the required energy input ( George & Fabian 1991; Murphy & Yaqoob 2009 ) . For reference, an extended Fe K α region in NGC 5728  \nDistances between Multiwavelength Centroids Found in This Work: HST  \nTable 7 / ACS, VLA-A, and Chandra / ACIS-S ( Fe K α )  \n( D ∼ 41 Mpc ) was reported with an integrated luminosity of L ( 6.1 -6.6 keV ) ∼ 6.7 × 10 39 erg s -1 within a 1 5 -8 " annular region ( excluding the nucleus; Trindade Falcão 2023, 2024 ) . This luminosity is 1 order of magnitude lower than those of the individual Fe K α peaks in MCG-03-34-64 ( Table 3 ) . In contrast, a luminosity of L ( 6.1 -6.6 keV ) ∼ 6.0 × 10 40 erg s -1 was found within the inner 1 5 of NGC 5728 ( Trindade Falcão 2023, 2024 ) , comparable to the luminosities of the individual Fe K α peaks in MCG-03-34-64. While the emission in NGC 5728 was attributed to shocked gas in the ISM ( Trindade Falcão 2023, 2024 ) , a similar origin in MCG-03-3464 seems unlikely due to the high luminosities of the individual regions. Therefore, it is unlikely that collisional ionization in the NLR alone powers the individual Fe K α peaks observed in MCG-03-34-64.  \nThe energetic mismatch is strengthened by the similar scattering corrected intrinsic X-ray luminosities calculated for these galaxies, = L 2 10 keV int. ( -) 1.5 × 10 43 erg s -1 for MCG-0334-64 ( Table 5 ) , versus L ( 2 -10 keV ) ∼ 1.2 × 10 43 erg s -1 for NGC 5728 ( Zhao et al. 2021 ) , suggesting that the differences in Fe K α luminosities between the two objects cannot be attributed solely to differences in the intrinsic luminosity of the nuclear source.', '4.2.2. Dual AGN + Shocked ISM': 'Following the discussion in Section 4.2.1, the high Fe K α luminosities ( Table 3 ) , and the high energies required for the production of such line emission suggest that both Fe K α regions could be powered by an active SMBH. In this scenario, the northern emission centroids would pinpoint the location of one AGN, while the central emission centroids ( radio, Fe K α , and optical ) may be associated with a second active SMBH in this system ( given the high /uniFB02 uxes found for the central optical region in the F814W continuum band; Table 6 ) . The distances measured between the different centroids attributed to each AGN are consistent across different wave bands ( Table 7 ) , supporting the dual AGN scenario.  \nIn this scenario, the southern [ OIII ] region may arise from collisionally ionized emission in the ISM. Shock emission from jet -ISM interaction at the southern optical centroid location is supported by ( 1 ) the morphology of the radio emission at the central radio peak, which is observed to bend southward ( Schmitt et al. 2001 ) , in the direction of the southern [ OIII ] peak ( see Figures 5 and 3, and Section 3.1.3 ) ; and ( 2 ) the morphology of the Chandra 0.3 -3 keV ( soft ) emission, which is also observed to bend southward in the direction of the southern [ OIII ] centroid ( see Figure 3, and Section 3.1.3 ) . The lack of a corresponding southern radio or hard X-ray  \ncounterpart is consistent with the hypothesis that the northern and central [ OIII ] centroids are powered by individual AGNs.  \nThe dual AGN scenario is strengthened by the consistent values between the estimated bolometric luminosities derived from [ OIII ] and X-rays ( Tables 5 and 6 ) , and the detection of nearly equal Fe K α emission peaks in the Chandra image, a powerful tool for identifying and con /uniFB01 rming dual AGN systems ( De Rosa et al. 2022 ) . In the 3 -5 keV Chandra image ( Figure 3 ) , the northern nucleus appears brighter than the central nucleus, although some extended emission is observed toward the central nucleus in this band. The column densities in the transmission of the two nuclei may not be the same, i.e., the AGN located at the central Fe K α region could have a higher absorbing column and appear fainter at lower energies. A higher contribution from the photoionized + thermal extended soft X-ray gas is expected at these lower energies. Given the resolution of the analyzed X-ray data, performing a separate spectral analysis of each individual Fe K α region is currently impractical.', '5. Summary and Conclusions': "We analyze new HST / ACS and Chandra / ACIS observations of the nearby Seyfert galaxy MCG-03-34-64, along with archival HST / ACS, XMM-Newton / Epic-pn, NuSTAR, and VLA-A data sets. Our analysis reveals the following:  \nIn the optical with HST: Three emission regions are uncovered in the narrow-band continuum-subtracted [ OIII ] and F814W continuum images ( Figures 1, 2 ) , separated by 79 ± 8pc ( northern to central ) ,and76 ± 8 ( central to southern ) . Diffraction spikes from the bright [ OIII ] -emitting regions suggest individual sizes of /uniF088 60 pc. We derive individual bolometric luminosities in the range ~ -' L 0.73 1.6 10 bol 44 OIII [ ] erg s -1 , and an integrated bolometric luminosity = ' L 4.8 10 bol 44 OIII [ ] erg s -1 from [ OIII ] ( Table 6 ) .  \nIn X-rays with Chandra: Two spatially resolved emission centroids are detected in the 6.2 -6.6 keV Fe K α image, separated by 125 ± 21 pc. These peaks are evident in individual exposures and the merged data set. The northern and central Fe K α regions have 36 ± 6 and 37 ± 6 counts in the narrow 6.2 -6.6 keV band, respectively, corresponding to /uni0085 6 σ detections, and nearly equal Fe K α luminosities, L ( 6.2 -6.6 keV ) ∼ 3.2 ± 0.6 × 10 40 erg s -1 .  \nIn X-rays with XMM and NuSTAR: Fitting the joint spectrum with a physically motivated model yields an integrated Fe K α luminosity of L ( 6.2 -6.6 keV ) = 1.0 × 10 41 erg s -1 , a scattered ( re /uniFB02 ected ) continuum luminosity of = ' -L 2.1 2 10 keV scat. ( ) 10 41 erg s -1 , and a total, scattering corrected intrinsic 2 -10 keV luminosity of = -L 2 10 keV int. ( ) 1.5 × 10 43 erg s -1 .  \nIn the radio with VLA: Two emission regions are observed in the 3.6 cm VLA continuum image ( Schmitt et al. 2001 ) , spatially colocated with the northern and central Fe K α and [ OIII ] regions.  \nWe propose two possible physical interpretations of our results, and discuss these in the context of our analysis:  \n- 1. The ' single AGN + shocked ISM ' scenario, which proposes the existence of a single active nucleus in the system, while the remaining multiwavelength centroids may be attributed to the interaction of the ISM with the radio jet in the NLR.  \nThis scenario is strengthened by:  \na. Previous X-ray studies on this source, which /uniFB01 nd evidence for a mix of collisionally and photoionized X-ray gas in the NLR ( Miniutti et al. 2007 ) .  \nThis scenario is challenged by:  \nb. The high Fe K α luminosities derived for individual regions and the energies required for the production of such line emission.  \n- 2. The ' dual AGN + shocked ISM ' scenario, which proposes the existence of a dual SMBH pair in this system separated by just 125 ± 21 pc.  \nThis scenario is strengthened by:  \n- a. The detection of two spatially resolved Fe K α regions in the Chandra imaging data, with high individual luminosities ( Table 3 ) .\n- b. The detection of three very bright and compact ( < 60 pc diameter ) [ OIII ] -emitting regions in the HST imaging data, and the respective individual bolometric luminosities ( Table 6 ) .\n- c. The detection of spatially coincident Fe K α , radio, and [ OIII ] centroids at the northern and central regions. This is the /uniFB01 rst time spatially resolved, multiwavelength emission centroids in X-rays, radio, and optical are detected colocated in a nearby candidate dual AGN. For comparison, the recent study of Koss et al. ( 2023 ) , which identi /uniFB01 ed the presence of a dual AGN system separated by ∼ 230 pc in UGC 4211, detected colocated optical ( HST F814W, MUSE AO [ OIII ] , and H α ) , NIR ( Keck J and K ' ) , and submillimeter ( Atacama Large Millimeter / submillimeter Array continuum at ∼ 230 GHz ) centroids at the position of the two nuclei, but with no con /uniFB01 rmation from X-rays or radio observations.  \nIn summary, although we cannot de /uniFB01 nitively con /uniFB01 rm or exclude the physical scenarios presented here, identi /uniFB01 cation of the two nuclei in a deeper Chandra exposure would help to con /uniFB01 rm a possible dual black hole system in this galaxy. Analysis of gas kinematics in the nuclear region of MCG-03-34-64 is crucial to determine the nature of the observed structures. Kinematic informationobtainedwithHST / STIS long-slit spectroscopy could reveal disturbed kinematics expected from either the individual out /uniFB02 ows of two SMBHs or the highly disturbed kinematics resulting from the merger environment. This information cannot be obtained from the archival X-Shooter data and requires the resolution of HST to probe the ∼ 100 pc region of interest.", 'ORCID iDs': 'Anna Trindade Falcão https: // orcid.org / 0000-0001-8112-3464 T. J. Turner https: // orcid.org / 0000-0003-2971-1722 S. B. Kraemer https: // orcid.org / 0000-0003-4073-8977 J. Reeves https: // orcid.org / 0000-0003-3221-6765 V. Braito https: // orcid.org / 0000-0002-2629-4989 H. R. Schmitt https: // orcid.org / 0000-0003-2450-3246 L. Feuillet https: // orcid.org / 0000-0002-5718-2402', 'References': '- Blecha, L., Snyder, G. F., Satyapal, S., & Ellison, S. L. 2018, MNRAS,\n- 478, 3056 Burke-Spolaor, S. 2011, MNRAS, 410, 2113 Castignani, G., Combes, F., & Salomé, P. 2020, A&A, 635, L10 Dadina, M., & Cappi, M. 2004, A&A, 413, 921 de Grijp, M. H. K., Keel, W. C., Miley, G. K., Goudfrooij, P., & Lub, J. 1992, A&AS, 96, 389 De Lucia, G., & Blaizot, J. 2007, MNRAS, 375, 2 De Rosa, A., Vignali, C., Severgnini, P., et al. 2022, MNRAS, 519, 5149 Dotti, M., Sesana, A., & Decarli, R. 2012, AdAst, 2012, 940568 Fabbiano, G., & Elvis, M. 2024, in The Interaction of the Active Nucleus with the Host Galaxy Interstellar Medium, ed. C. Bambi ( Singapore: Springer ) , 4399 Fabbiano, G., Paggi, A., Karovska, M., et al. 2018a, ApJ, 865, 83 Fabbiano, G., Paggi, A., Siemiginowska, A., & Elvis, M. 2018b, ApJL, 869, L36 Fabbiano, G., Wang, J., Elvis, M., & Risaliti, G. 2011, Natur, 477, 431 Fischer, T. C., Kraemer, S. B., Schmitt, H. R., et al. 2018, ApJ, 856, 102 Fornasini, F. M., Elvis, M., Maksym, W. P., et al. 2022, ApJ, 931, 65 Fu, H., Myers, A. D., Djorgovski, S. G., & Yan, L. 2011a, ApJ, 733, 103 Fu, H., Steffen, J. L., Gross, A. C., et al. 2018, ApJ, 856, 93 Fu, H., Yan, L., Myers, A. D., et al. 2011b, ApJ, 745, 67 George, I. M., & Fabian, A. C. 1991, MNRAS, 249, 352 Ginol /uniFB01 , M., Maiolino, R., Nagao, T., et al. 2017, MNRAS, 468, 3468 Goulding, A. D., Pardo, K., Greene, J. E., et al. 2019, ApJL, 879, L21 Hönig, S. F., Kishimoto, M., Gandhi, P., et al. 2010, A&A, 515, A23\n- Hopkins, A. M., McClure-Grif /uniFB01 ths, N. M., & Gaensler, B. M. 2008, ApJ, 682, L13\n- Kharb, P., Lal, D. V., & Merritt, D. 2017, NatAs, 1, 727\n- Kollatschny, W., Grupe, D., Parker, M. L., et al. 2020, A&A, 638, A91 Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511 Koss, M., Mushotzky, R., Treister, E., et al. 2011, ApJL, 735, L42 Koss, M., Mushotzky, R., Veilleux, S., & Winter, L. 2010, ApJL, 716, L125 Koss, M., Romero-Cañizales, C., Baronchelli, L., et al. 2015, ApJ, 807, 149 Koss, M. J., Assef, R., Balokovi ć , M., et al. 2016, ApJ, 825, 85 Koss, M. J., Treister, E., Kakkad, D., et al. 2023, ApJL, 942, L24 Kraemer, S. B., & Crenshaw, D. M. 2000, ApJ, 535, 763 Lamastra, A., Bianchi, S., Matt, G., et al. 2009, A&A, 504, 73 Liu, X., Shen, Y., & Strauss, M. A. 2011, ApJL, 736, L7 Ma, J., Maksym, W. P., Fabbiano, G., et al. 2021, ApJ, 908, 155 Maksym, W. P., Fabbiano, G., Elvis, M., et al. 2016, ApJ, 829, 46 Maksym, W. P., Fabbiano, G., Elvis, M., et al. 2017, ApJ, 844, 69 Maksym, W. P., Fabbiano, G., Elvis, M., et al. 2019, ApJ, 872, 94 Mayer, L., Kazantzidis, S., Madau, P., et al. 2007, Sci, 316, 1874\n- Miniutti, G., Ponti, G., Dadina, M., Cappi, M., & Malaguti, G. 2007, MNRAS, 375, 227\n- Murphy, K. D., & Yaqoob, T. 2009, MNRAS, 397, 1549\n- Osterbrock, D. E., & Robertis, M. M. D. 1985, PASP, 97, 1129\n- Paggi, A., Wang, J., Fabbiano, G., Elvis, M., & Karovska, M. 2012, ApJ, 756, 39\n- Peng, C. Y., Ho, L. C., & Impey, C. D. 2010, AJ, 139, 2097\n- Peng, C. Y., Ho, L. C., Impey, C. D., & Hans-Walter, R. 2002, AJ, 124, 266  \nPier, E. A., Antonucci, R., Hurt, T., Kriss, G., & Krolik, J. 1994, ApJ, 428, 124  \nRevalski, M., Dashtamirova, D., Crenshaw, D. M., et al. 2018, ApJ, 867, 88  \nRicci, C., Privon, G. C., Pfei  \n/uniFB02  \ne, R. W., et al. 2021, MNRAS, 506, 5935  \n- Schmitt, H. R., Ulvestad, J. S., Antonucci, R. R. J., & Kinney, A. L. 2001, ApJS, 132, 199  \nSmith, K. L., Shields, G. A., Bonning, E. W., et al. 2010, ApJ, 716, 866  \n- Steinborn, L. K., Dolag, K., Comerford, J. M., et al. 2016, MNRAS, 458, 1013\n- Stemo, A., Comerford, J. M., Barrows, R. S., et al. 2021, ApJ, 923, 36\n- Tatum, M. M., Turner, T. J., Miller, L., et al. 2016, ApJ, 818, 12\n- Treister, E., Messias, H., Privon, G. C., et al. 2020, ApJ, 890, 149\n- Trindade Falcão, A., Fabbiano, G., Elvis, M., Paggi, A., & Maksym, W. P. 2023, ApJ, 950, 143\n- Trindade Falcão, A., Fabbiano, G., Elvis, M., Paggi, A., & Maksym, W. P. 2024, ApJ, 963, 6\n- Trindade Falcão, A., Kraemer, S. B., & Fischer, T. C. 2021, MNRAS, 500, 1491\n- Vasudevan, R. V., & Fabian, A. C. 2007, MNRAS, 381, 1235\n- Verbiest, J. P. W., Lentati, L., Hobbs, G., et al. 2016, MNRAS, 458, 1267\n- Veres, P. M., Gabányi, K. É, Frey, S., et al. 2021, ApJ, 922, 99\n- Wang, J., Fabbiano, G., Elvis, M., et al. 2011a, ApJ, 736, 62\n- Wang, J., Fabbiano, G., Risaliti, G., et al. 2011b, ApJ, 729, 75\n- Wassenhove, S. V., Volonteri, M., Mayer, L., et al. 2012, ApJL, 748, L7 Zhao, X., Marchesi, S., Ajello, M., et al. 2021, A&A, 650, A57'}