<figure><img src="image_1.png" /></figure>

HAL
open science

# The Close AGN Reference Survey (CARS)

M. Powell, B. Husemann, G. Tremblay, M. Krumpe, T. Urrutia, S. Baum, G.

Busch, F. Combes, S. Croom, T. Davis, et al.

<figure><img src="image_2.png" /></figure>

## To cite this version:

M. Powell, B. Husemann, G. Tremblay, M. Krumpe, T. Urrutia, et al. The Close AGN Reference Survey (CARS). Astronomy & Astrophysics - A&A, 2018, 618, pp.A27. 10.1051/0004-6361/201833424. hal-02347687

HAL Id: hal-02347687

https://hal.science/hal-02347687v1

Submitted on 6 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# The Close AGN Reference Survey (CARS)

## No evidence of galaxy-scale hot outflows in two nearby AGN*

M. C. Powell<fnref n="1" />, B. Husemann<fnref n="2" />, G. R. Tremblay<fnref n="3" />, M. Krumpe<fnref n="4" />, T. Urrutia<fnref n="4" />, S. A. Baum<fnref n="5" />,<fnref n="6" />, G. Busch<fnref n="7" />, F. Combes<fnref n="8" />,
S. M. Croom<fnref n="9" />, T. A. Davis<fnref n="10" />, A. Eckart<fnref n="7" />,<fnref n="11" />, C. P. O'Dea<fnref n="5" />,<fnref n="12" />, M. Pérez-Torres<fnref n="13" />,<fnref n="14" />, J. Scharwächter<fnref n="15" />,
I. Smirnova-Pinchukova<fnref n="2" />, and C. M. Urry<fnref n="1" />

Received 14 May 2018 / Accepted 30 June 2018

## ABSTRACT

**Aims.** We probe the radiatively-efficient, hot wind feedback mode in two nearby luminous unobscured (type 1) AGN from the Close AGN Reference Survey (CARS), which show intriguing kpc-scale arc-like features of extended [O iii]ionized gas as mapped with VLT-MUSE. We aimed to detect hot gas bubbles that would indicate the existence of powerful, galaxy-scale outflows in our targets, HE 0227–0931 and HE 0351+0240, from deep (200 ks) *Chandra* observations.

**Methods.** By measuring the spatial and spectral properties of the extended X-ray emission and comparing with the sub kpc-scale IFU data, we are able to constrain feedback scenarios and directly test if the ionized gas is due to a shocked wind.

**Results.** No extended hot gas emission on kpc-scales was detected. Unless the ambient medium density is low ($n_H \sim 1 \text{ cm}^{-3}$ at 100 pc), the inferred upper limits on the extended X-ray luminosities are well below what is expected from theoretical models at matching AGN luminosities.

**Conclusions.** We conclude that the highly-ionized gas structures on kpc scales are not inflated by a hot outflow in either target, and instead are likely caused by photoionization of pre-existing gas streams of different origins. Our nondetections suggest that extended X-ray emission from an AGN-driven wind is not universal, and may lead to conflicts with current theoretical predictions.

**Key words.** galaxies: active – galaxies: evolution

## 1. Introduction

Supermassive black holes undergo phases of extreme accretion, observed as active galactic nuclei (AGN), which give off enormous amounts energy that is deposited into the host galaxy. Tight correlations between black hole and galaxy properties suggest a coevolution between them, in which the energy couples to the surrounding environment and may influence the state of star

formation in the galaxy (Silk & Rees 1998). Such feedback is generally invoked in simulations to recover the observed local luminosity function, galaxy bimodality in color space, and the number density of massive galaxies (e.g., Di Matteo et al. 2005; Bower et al. 2006; Schaye et al. 2015; Weinberger et al. 2018).

One manifestation of AGN feedback is the radiatively-efficient hot wind mode, which has been predicted in semi-analytical (e.g., Weinberger et al. 2018) and analytical models (Faucher-Giguère & Quataert 2012; Zubovas & King 2012). This is when radiation pressure from the accretion disk pushes the surrounding gas outward, driving a hot, shocked, bipolar wind if overpressured. Such a galactic “superwind” would sweep up and shock the surrounding ambient gas as it expands, possibly causing a suppression of star formation if energetic enough

<footnote n="1">¹ Yale Center for Astronomy and Astrophysics, and Physics Department, Yale University, PO Box 2018120, New Haven, CT 06520-8120, USA
e-mail: meredith.powell@yale.edu</footnote>

<footnote n="2">² Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany</footnote>

<footnote n="3">³ Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA</footnote>

<footnote n="4">⁴ Leibniz-Institut für Astrophysik Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany</footnote>

<footnote n="5">⁵ Department of Physics & Astronomy, University of Manitoba, Winnipeg, MB R3T 2N2, Canada</footnote>

<footnote n="6">⁶ Carlson Center for Imaging Science, Rochester Institute of Technology, Rochester, NY 14623, USA</footnote>

<footnote n="7">⁷ I. Physikalisches Institut der Universität zu Köln, Zülpicher Str. 77, 50937 Köln, Germany</footnote>

<footnote n="8">⁸ LERMA, Observatoire de Paris, Collège de France, PSL Univ., CNRS, Sorbonne Univ., 75014 Paris, France</footnote>

<footnote n="9">⁹ Sydney Institute for Astronomy, School of Physics, A28, The University of Sydney, NSW 2006, Australia</footnote>

<footnote n="10">¹⁰ School of Physics & Astronomy, Cardiff University, Queens Buildings, The Parade, Cardiff CF24 3AA, UK</footnote>

<footnote n="11">¹¹ Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany</footnote>

<footnote n="12">¹² School of Physics & Astronomy, Rochester Institute of Technology, Rochester, NY 14623, USA</footnote>

<footnote n="13">¹³ Instituto de Astrofísica de Andalucía, Glorieta de la Astronomía s/n, 18008 Granada, Spain</footnote>

<footnote n="14">¹⁴ Departamento de Física Teórica, Facultad de Ciencias, Universidad de Zaragoza, 50009 Zaragoza, Spain</footnote>

<footnote n="15">¹⁵ Gemini Observatory, Northern Operations Center, 670 N. A'ohoku Pl., Hilo, Hawaii 96720, USA</footnote>

<footnote>* The scientific results reported in this article are based on observations made by the *Chandra* X-ray Observatory (PI: G. Tremblay, ID: 17700519), and on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO program(s) 094.B-0345(A) (PI: B. Husemann).</footnote>
to expel the gas from the galaxy. Observationally, this outflow would be characterized by thin ionized shells surrounding a hot X-ray-bright bubble (Heckman 1990), similar to the galaxy-scale winds driven by supernovae and/or massive stars. Analytical models of AGN-driven winds have been proposed (Faucher-Giguère & Quataert 2012) based on observations of nearby quasars, with observational signatures that can distinguish them from stellar winds. Specifically, the luminosity and morphology of spatially-extended X-ray emission from the forward shock into the surrounding interstellar medium (ISM) is predictable, given assumptions on the surrounding ambient density profile (Nims et al. 2015) and the coupling efficiency between the AGN and its environment.

There have been previous observations of extended hot X-ray emission around nearby AGN, for which radiative feedback is consistent. This includes the luminous, obscured radio-quiet quasar SDSS J135646.10+102609.0 (Greene et al. 2014), selected for its high intrinsic luminosity and strong and extended ionized gas. Additionally there is the Type 2 “Teacup AGN” (SDSS J1430+1339; Lansbury et al. 2018), also selected for its luminous, extended, as well as broad (FWHM > 700 km s⁻<fnref n="1" />) [O III] emission, along with the nearest ultraluminous infrared galaxy (ULIRG) Mrk 231 (Veilleux et al. 2014). The fad-ing obscured AGN IC 2497 (“Hanny’s Voorwerp” system; Sartori et al. 2016) also shows evidence for hot, diffuse gas from a past high accretion state. However the ubiquity of quasar-driven superbubbles is currently unclear, especially for unob-scured AGN whose X-ray emission is dominated by the central engine. Deep, spatially resolved X-ray data is currently lacking for a representative sample of such objects.

In this paper we present a study of two luminous Type 1 AGN aiming to explore the connection of ionized gas filaments on kpc-scales and a presumed expanding hot gas bubble powered by the AGN radiation. We combine 200 ks of new Chandra X-ray observations in conjunction with optical integral-field unit (IFU) observations, in order to distinguish between the hot gas outflow scenario and the case where the [O III] arcs are gas clouds illuminated by the AGN, which have been previously expelled from the galaxy. Throughout this paper we assume $H_0 = 70$ km s⁻<fnref n="1" /> Mpc⁻<fnref n="1" />, $\Omega_M = 0.3$, and $\Omega_\Lambda = 0.7$.

## 2. The Close AGN Reference Survey

The Close AGN Reference Survey (CARS; Husemann et al. 2017, and in prep.<fnref n="1" />) aims to constrain the galaxy-scale feeding vs. feedback in a larger sample of nearby luminous AGN via spatially-resolved observations across the entire X-rays to radio electromagnetic spectrum. The CARS targets are drawn from a mother sample of 99 luminous AGN detected in the Hamburg-ESO survey (HES, Wisotzki et al. 2000) up to $z < 0.06$, of which 39 galaxies were targeted with the IRAM 30m telescope to characterize the molecular gas content via ¹²CO(1–0) line observations (Bertram et al. 2007). Focusing on the Bertram et al. subsample, CARS has undertaken a full 3D (spectral + spatial) mapping campaign in the optical, primarily using the Multi-Unit Spectroscopic Explorer (MUSE, Bacon et al. 2010, 2014) at the Very Large Telescope (VLT). Those observations are being complemented with additional radio, IR, optical, UV and X-ray observations to get a complete picture of the stellar and multi-phase gas content of the host galaxies in relation to the AGN properties.

<table><caption>**Table 1.** Basic parameters for each galaxy.</caption><thead><tr><th>Galaxy</th><th>z</th><th>M<sub>⋆</sub><sup>a</sup><br>(M<sub>☉</sub>)</th><th>M<sub>BH</sub><sup>b</sup><br>(M<sub>☉</sub>)</th><th>i<sup>c</sup></th><th>kpc/''</th></tr></thead><tbody><tr><td>HE 0227–0913</td><td>0.016</td><td>1.3 × 10<sup>10</sup></td><td>1.4 × 10<sup>6</sup></td><td>67 ± 5</td><td>0.33</td></tr><tr><td>HE 0351+0240</td><td>0.036</td><td>3.2 × 10<sup>10</sup></td><td>1.3 × 10<sup>7</sup></td><td>N/A</td><td>0.72</td></tr></tbody></table>

**Notes.** (a) Stellar masses based on the QSO-subtracted $g-i$ host color and $M_i$ absolute host magnitude following Taylor et al. (2011). (b) BH masses taken from Wang et al. (2014) for HE 0227–0913 and Koss et al. (2017) for HE 0351+0240. (c) Galaxy disk inclination estimated only for HE 0227–0913 from the gas kinematic velocity.

From the CARS sample we selected two galaxies, HE 0227–0913 (also known as Mrk 1044) and HE 0351+0240, for deep pilot Chandra observations to detect extended hot diffuse gas. The high accretion rates of the AGN and detection of kpc-scale ionized gas filaments make both galaxies good candidates to test the presence of expanding X-ray bubbles as expected for the radiative AGN feedback mode. The two galaxies presented here are complementary in studying the SMBH-galaxy connection, as they have contrasting host galaxy properties and are in very different evolutionary stages. HE 0227–0913 is a narrow-line Seyfert 1 galaxy (NLS1) at $z = 0.016$ with a disk-like morphology. By reverberation mapping, its black hole mass is estimated to be $1.4 \times 10^6 M_☉$ (Wang et al. 2014). HE 0351+0240, on the other hand, is in an ongoing major merger at $z = 0.036$, with an estimated black hole mass of $1.3 \times 10^7 M_☉$ (Koss et al. 2017). The characteristics of these galaxies are summarized in Table 1.

## 3. Observations and data reduction

### 3.1. MUSE integral-field spectroscopy

We observed HE 0227–0913 and HE 0351+0240 with MUSE at the VLT the nights of 10 January 2015 and 26 December 2014, respectively. The two targets were observed for 1200 s and 1600 s integration times, each split into two independent exposures which were rotated by 90° with a small dither offset to aid in the cleaning of cosmic rays and flat-fielding artifacts. The wide-field mode (WFM) of MUSE provides a nearly rectangular field-of-view (FoV) with a size of $1' \times 1'$ and a sampling of 0.2 arcsec. The wavelength range is 4750 Å–9300 Å with a spectral resolution of ~2.5 Å that slightly varies with wavelength (Bacon et al. 2017; Guérou et al. 2017).

The MUSE data were reduced with the official ESO pipeline (Version 1.6, Weilbacher et al. 2012) performing standard calibration tasks, such as bias subtraction, slice extraction, flat-fielding, wavelength calibration using arc lines, and flux calibration based on spectrophotometric standard stars. The fully calibrated individual exposures were then combined and resampled using the “drizzle” algorithm (Fruchter & Hook 2002) to a 3D datacube with rectified wavelength range. Sky subtraction was done by creating a master sky spectrum from blank-sky regions within the MUSE FoV. We further suppressed sky-line residuals with a custom-made principal component analysis approach like ZAP (Soto et al. 2016), but optimized and simplified for the given purpose and data set (Husemann et al. in prep.).

Figure 1 shows the continuum light for HE 0227–0913 and HE 0351+0240. The B, V, and R three-color composite of HE 0227–0913 is shown, taken as part of an optical and NIR

<footnote n="1">¹ www.cars-survey.org</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Three-color composite images of the two galaxies for which we present new *Chandra* and MUSE data in this paper. **Left:** SOAR Optical Imager (SOI) *b*-, *v*-, and *r*-band composite of HE 0227-0913 (Markarian 1044). These bands have been assigned to the blue, green, and red channels, respectively. **Right:** composite made by collapsing three slices through the MUSE cube across wavelength ranges that roughly correspond to the *g*-, *r*-, and *i*-bands, again respectively assigned to the blue, red, and green channels. We have applied a logarithmic stretch to both panels to better show the low surface brightness galaxy outskirts alongside the bright type 1 nuclei. Both images have been rotated, with east left and north up.</figcaption></figure>

followup campaign for the CARS galaxies with the SOAR Optical Imager (SOI; Walker et al. 2003) on the SOAR telescope, (SOAR/NOAO program 2016A-0006 PI: G. Tremblay). For HE 0351+0240, the color image was created from the *gri* broad-band image extracted from the MUSE datacube. HE 0227-0913 appears as an inclined disk galaxy with a rather faint extended disk and spiral arms, while HE 0351+0240 appears as an ongoing major merger system where the two galaxy cores are still separated.

## 3.2. Chandra X-ray spectroscopy

We observed HE 0227-0913 with the *Chandra* Advanced CCD Imaging Spectrometer (ACIS-S; PI: G. Tremblay) on 12 October 2015 (obs. ID 18188), 10-14-2015 (obs. ID 18684), and 16 October 2015 (obs. ID 18685) for a total of 93.4 ks. HE 0351+0240 was observed on 14 December 2015 (obs. ID 18189) and 15 December 2015 (obs. ID 18726) for a total of 93.7 ks. Due to the high flux of the objects, we used a 1/8 subarray to minimize pileup. The data were reduced and analyzed with the CIAO 4.9 software (Fruscione et al. 2006). The “chandra_repro” script was used with the CALDB 4.7.3 set of calibration files to process the level 1 event files. We filtered the data to eliminate cosmic rays and bad grades, and verified no significant changes in the background with time.

# 4. Data analysis

## 4.1. Optical emission-line analysis

Since HE 0227-0913 and HE 0351+0240 are both unobscured QSOs, we first have to subtract the point-like unresolved AGN

emission before performing the emission line mapping across the two galaxies. Here, we use the AGN-host galaxy IFU deblending software QDeblend³D (Husemann et al. 2013, 2014). The tool utilizes a simple nonparametric algorithm to firstly determine the on-axis PSF through intensity mapping of the broad $H\beta$ and $H\alpha$ lines, and secondly subtract the central AGN spectrum convolved by the PSF in an iterative manner to deconvominate the AGN spectrum from the host galaxy contributions at the center of the galaxy. The original spectrum within an aperture of 3'' and the respective AGN and host galaxy contribution are shown in Fig. 2 for both galaxies. The narrow line components had a best-fit velocity dispersion of $\sigma = 131 \text{ km s}^{-1}$ for HE 0227-0913, and $65 \text{ km s}^{-1}$ for HE 0351+0240.

As a next step, we fit the stellar continuum and emission lines across the QSO-subtracted data cube using PyParadise (Husemann et al. 2016a; Weaver et al. 2018). Since the signal-to-noise ratio (S/N) for both galaxies is not enough to model the continuum for individual spaxels over the entire galaxies, we first perform a Voronoi-binning scheme (Cappellari & Copin 2003) to achieve a S/N of at least 20 per bin. The Voronoi-binned spectra are then fitted with a linear superposition of stellar spectra from the INDO-U.S. library (Valdes et al. 2004), convolved with the MUSE mean spectral resolution of 2.4 Å, and a redshifted Gaussian kernel to match the kinematics. In a second step, we model the continuum in the same way for all individual spaxels but keeping the previously inferred kinematics fixed. After subtracting the best-fit continuum model, we fit all bright emission lines such as $H\beta$, [O III] $\lambda\lambda4960,5007$, $H\alpha$, [N II] $\lambda\lambda6548,6583$ and [S II] $\lambda\lambda6716,6730$, with a Gaussian shape, again convolved with the instrumental resolution, where the redshift and intrinsic velocity dispersion are assumed to be the same for all lines. The
<figure><img src="image_5.png" /><figcaption>**Fig. 2.** Integrated central MUSE aperture spectra within a diameter of 3" for HE 0227-0913 and HE 0351+0240. The total spectrum (black lines) is decomposed into a point-like QSO and extended host-galaxy contribution using the nonparametric QSO-host galaxy decomposition algorithm implemented in QDeblend3D (Husemann et al. 2013, 2014).</figcaption></figure>

[O iii] and [N ii] doublets are fixed in their line ratio according to the theoretical expectations (Storey & Zeippen 2000). Errors on all continuum and emission-line model parameters are estimated using a Monte-Carlo approach where we re-fit the spectra about 30 times after replacing each pixel within the normal distribution around the initial values with a width guided by the pipeline errors.

The fitting results for HE 0227-0913 and HE 0351+0240 are summarized in Fig. 3 and Fig. 4, respectively, featuring the ionized gas velocity and velocity dispersion as well as the Hα and [O iii] surface brightness maps. Based on the classical emission-line diagnostic diagrams (e.g., Baldwin et al. 1981; Veilleux & Osterbrock 1987), we create an excitation map using demarcation lines in the [O iii]/Hβ vs. [N ii]/Hα to distinguish between star forming, AGN ionization, intermediate and LINER-like regions (Kewley et al. 2001; Kauffmann et al. 2003; Stasińska et al. 2008). The resulting properties of both galaxies are quite different and we describe the main features in the following.

The host galaxy of HE 0227-0913 is a disk which is clearly supported by the ionized gas velocity field. This velocity field was modeled using the KinMS software<fnref n="2" /> (Davis et al. 2013) to obtain an inclination angle of 67.4 ± 4.5°. Ongoing star formation dominates the gas excitation in a circumnuclear disk, with an inner star forming ring a few arcsec (<1 kpc) in diameter and an apparently outer star forming ring at ~20" (7 kpc) distance that is likely related to the spiral arms. The presence of regions dominated by ionization from OB-stars is consistent with star formation within the ISM of the galaxy; it has an estimated H₂ mass of 4 × 10⁸ M⊙ (Bertram et al. 2007), and an HⅠ mass of 2.6 × 10⁹ M⊙ (König et al. 2009). Another striking feature is a bright AGN-ionized arc-like feature of ionized gas about 5–15" (or 5 kpc) North of the nucleus, which is redshifted in velocity

space with respect to the underlying disk motion but nearly at zero offset from the systemic velocity. The high velocity dispersion may not entirely be intrinsic to the AGN-ionized gas, but additionally broadened due to the superposition of the emission line kinematics of the gas disk.

The interacting galaxy system of HE 0351+0240 does not show any strong signatures of ongoing star formation (there is no $^{12}$CO(1–0) detection for this galaxy; Bertram et al. 2007), as all of the ionized gas is ionized by the AGN and displays a filamentary structure with a complex kinematic pattern. The velocity dispersion is elevated close to the AGN nucleus, but not in the extended filamentary structures.

In the following we will test the hypothesis that the AGN-ionized gas filaments represent the interface of a slowly expanding hot gas bubble with the surrounding cool gas. Such hot gas bubbles are predicted by numerical simulation (Weinberger et al. 2018) and other theoretical considerations (Faucher-Giguère & Quataert 2012; Zubovas & King 2012) as a natural consequence of radiation-driven feedback scenarios.

## 4.2. X-ray analysis

We searched for extended X-ray emission by comparing our observations to the *Chandra* PSF, the accurate modeling of which requires the point-source spectrum. As Type 1 AGN, the nuclei of HE 0227-0913 and HE 0351+0240 are very bright in X-rays, with high enough fluxes for pileup to be significant in warping the spectrum. However, the shutterless ACIS detector on *Chandra* cause photons from a bright source to be detected while data are being read out, forming a linear streak through the source in the image. We can therefore use photons in this readout streak to model the point source spectrum for each source, as only the unresolved bright AGN is represented in this region without the pileup issues. We used the CIAO tool "dmextract" to extract the spectrum in each observation, and

<footnote n="2">² https://github.com/TimothyADavis/KinMS</footnote>
<figure><img src="image_6.png" /><figcaption>**Fig. 3.** Left panels: ionized gas velocity and velocity dispersion maps (upper panels) for HE 0227-0913 together with the Hα and [O III] surface brightness maps from the MUSE data after stellar continuum subtraction. Contours of the continuum emission are in logarithmic scaling are overlayed in gray. Right panels: emission-line diagnostic diagram for spaxels with a S/N > 3 in each line (upper panel). Excitation map after separating line ratios, into star forming, AGN, intermediate and LINER-like ionization.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 4.** Same as Fig. 3, for HE 0351+0240.</figcaption></figure>

the tools “mkacisrmf” and “mkarf” to make the response files at the location of the central source. We extracted the background spectra in adjacent regions next to the readout streaks, and we corrected the exposure time of the readout streak in

order to model the correct flux. Lastly, we stacked the individual observations for each galaxy together using the CIAO script “combine_spectra” and used SHERPA to model the 0.3–7 keV background-subtracted spectrum.
<figure><img src="image_8.png" /><figcaption>**Fig. 5.** Top: coadded Chandra images of HE 0227–0931 and HE 0351+0240 (with readout streaks) and the corresponding simulated PSFs. The color scale is in units of counts per second, shown by the right-hand colorbar. Bottom: spectra of the co-added, background subtracted point sources taken from the readout streaks (shown in the thick black regions above, with adjacent background regions), binned to have a minimum significance of 5 per bin. The spectra were fit to a single power-law, with the `xspec` model: (`phabs*powlaw1d`). In each model we fixed Galactic absorption to the value reported in `XSPEC`, $3.26 \times 10^{20} \text{ cm}^{-2}$ for HE 0227–0913 and $1.28 \times 10^{21} \text{ cm}^{-2}$ for HE 0351+0240. The best-fit models (solid black lines) were used to simulate the PSFs (shown above).</figcaption></figure>

AGN X-ray spectra are typically well modeled by a power-law, parametrized only by the photon index, normalization, and absorption. Fixing the absorption due to the galaxy and binning the spectrum to have a minimum S/N of 5 per bin, both galaxies were well fit by this model (“phabs*powlaw1d”). The fits are shown in Fig. 5 (lower panels) and the parameters of each model are listed in Table 2. NLS1 galaxies like HE 0227–0913 are known to have soft X-ray excess above the powerlaw continuum, and indeed, soft excess in HE 0227–0913 has been previously found and studied with higher-resolution instruments (Dewangan et al. 2007; Mallick et al. 2018). These studies require multiple components to adequately fit the full spectrum. However, the improvement of our fit to the streak spectrum by including an additional component brings the reduced $\chi^2$ below 1, and the resulting PSF does not change our conclusions. We therefore use the 1-component powerlaw fit for the AGN in both galaxies.

To simulate the PSF we used the ChaRT online ray tracing tool (Carter et al. 2003), where we input the modeled source spectrum, and then used the MARX software to project the rays onto the ACIS-S detector. We corrected for SIM drift and offset, and set the “AspectBlur” parameter to 0.19” to avoid pixelization effects in the pseudo-event file (the standard value for ACIS-S). The “marxpileup” program was used to simulate pileup, where the “FrameTime” parameter was set to 0.7 s, as it is reduced when using 1/8 of the chip. The coadded data are shown in Fig. 5, along with the simulated PSF, shown in a log color scale.

We used the CIAO tool `srcflux` to obtain the X-ray fluxes of each source using the resulting PSF models. We calculated the 0.3–10 keV flux to be $2.12 \times 10^{-11} \text{ erg cm}^{-2} \text{ s}^{-1}$ for HE 0227–0931 (consistent with Dewangan et al. 2007)

<table><caption>**Table 2.** X-ray spectral parameters for power-law model.</caption><thead><tr><th>Galaxy</th><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td rowspan="4">HE 0227-0931</td><td>Γ<sub>PL</sub></td><td>2.38 ± 0.06</td></tr><tr><td>norm<sub>PL</sub></td><td>(6.78 ± 0.21) × 10<sup>-3</sup></td></tr><tr><td>χ<sup>2</sup>/d.o.f.</td><td>0.98</td></tr><tr><td>f<sub>0.3-10 keV</sub></td><td>2.12 × 10<sup>-11</sup></td></tr><tr><td rowspan="4">HE 0351+0240</td><td>L<sub>2-10 keV</sub></td><td>9.2 × 10<sup>42</sup></td></tr><tr><td>Γ<sub>PL</sub></td><td>1.38 ± 0.08</td></tr><tr><td>norm<sub>PL</sub></td><td>(1.52 ± 0.10) × 10<sup>-3</sup></td></tr><tr><td>χ<sup>2</sup>/d.o.f.</td><td>1.06</td></tr><tr><td>f<sub>0.3-10 keV</sub></td><td>1.19 × 10<sup>-11</sup></td></tr><tr><td>L<sub>2-10 keV</sub></td><td>3.6 × 10<sup>43</sup></td></tr></tbody></table>

**Notes.** The 0.3–10 keV (absorbed) flux is in units of erg cm$^{-2}$ s$^{-1}$. The luminosities are in units of erg s$^{-1}$ and are corrected for Galactic absorption. The model fits were performed on the binned data using Sherpa.

and $1.91 \times 10^{-11} \text{ erg cm}^{-2} \text{ s}^{-1}$ for HE 0351+0240. This corresponds to luminosities of $L_{0.3-10 \text{ keV}} = 10^{43.1} \text{ erg s}^{-1}$ and $L_{0.3-10 \text{ keV}} = 10^{43.6} \text{ erg s}^{-1}$, respectively. Assuming a 2–10 keV bolometric correction of 20 (Lusso et al. 2012), we estimate the Eddington ratios ($\lambda_{\text{Edd}} \equiv L_{\text{bol}}/L_{\text{Edd}}$) to be ~1 for HE 0227–0931 and ~0.5 for HE 0351+0240.

We compared the radial surface brightness profiles of the data to the PSF images to look for evidence of extended emission. We used the CIAO tool “dmextract” to extract the events in 20 concentric annuli centered on each source from 1” to 20” (excluding the region of the readout streak). To estimate
<figure><img src="image_9.png" /><figcaption>**Fig. 6.** Radial profiles of the background-subtracted data (black circles) and of the simulated PSF (red lines) in bins of 1", shown with 95% level errors. Both 0.5–7 keV data and PSF were used, excluding the readout streak regions. *Bottom panels:* fractional residual errors. No clear evidence for extended emission is detected for either galaxy.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 7.** Difference between the data and PSF of the binned 2D images, divided by error for HE 0227–0913 (**left**) and HE 0315+0240 (**right**). No significant emission is detected beyond a 1.5σ excess in the first radial bin of HE 0227–0913 (1′′5–3′′).</figcaption></figure>

systematic errors from the incorrect pileup simulation, we estimated the fraction of counts affected by pileup via the “pileup_map” tool, which decreases as a function of radius, and added this error to the Poisson error. The surface brightness profiles of the PSF and source for both galaxies are shown in Fig. 6, along with the fractional residual errors. We do not find any significant evidence for excess emission above the PSF, except in the first bin for HE 0227–0913; we find a ~2–σ excess at R < 1 kpc, however the systematic uncertainty due to pileup in the central region makes this detection very tentative. On scales of the ionized gas structures (~3–5 kpc) no excess emission is found, even after dividing the data into soft and hard bins.

We further compared the flux in radial and azimuthal bins, to find any significant angular deviations from the PSF correlated with the ionized gas structures. Figure 7 shows the residual fluctuation maps (data minus PSF divided by the rms) in three radial bins and six azimuthal bins (the increasing size of the radial bins were to have comparable S/N in each bin). This shows there are no significant deviations from the modeled PSF when including pileup modeling errors.

# 5. Discussion

## 5.1. *Comparison to AGN wind theory*

Nims et al. (2015) investigated the observational signatures of AGN outflows, and predicted the resulting X-ray spectrum

and luminosity of shocked ambient material due to an AGN-driven wind out to 2 kpc. Assuming the kinetic outflow luminosity is 5% of the AGN luminosity, they predict 1–10 keV luminosities for various forward shock radii, which depends on the parameters of the ambient density profile (assuming a spherically symmetric profile $n_H = n_{H,0} \times (R_S/R_0)^{-1}$; Faucher-Giguère & Quataert 2012) and the bolometric luminosity of the AGN. Using fiducial values motivated by observations of nearby quasars and ULIRGs (Feruglio et al. 2010; Cicone et al. 2014), they showed that on kpc scales, the X-ray luminosity is dominated by free-free emission.

We compared their predictions to our results in the following way: we calculated the 3σ upper limits of our simulated PSFs and used PIMMS to convert the surface brightness to 1–10 keV flux, by assuming a thermal bremsstrahlung spectrum with a temperature in the analytic form for the shocked ambient medium in Nims et al. (2015), for three density values. We then scaled the Nims prediction for each density to the bolometric luminosities of the AGN, assuming a bolometric correction $\kappa_{2-10\text{keV}} = 20_{-12}^{+80}$ (Lusso et al. 2012), and extrapolated to 5 kpc. We note that this bolometric correction is conservative for Mrk 1044 because its black hole is super-Eddington; previous work looking at its SED has estimated it to have a bolometric luminosity an order of magnitude larger (~10⁴⁶ erg s⁻¹; Castelló-Mor et al. 2016), which would predict a more luminous wind. The comparison plots are shown in Fig. 8. For density values (at a radius of 100 pc) of 10–100 cm⁻<fnref n="3" />, our upper limits are well under the prediction. However for low densities ~1 cm⁻<fnref n="3" /> the prediction is much fainter, $L_{1-10\text{keV}} = 10^{39}-10^{41}$ erg s⁻¹, at or below our upper limits on the scales of the ionized gas. This suggests that either (1) the ambient medium density is low (≤1 cm⁻<fnref n="3" />) such that the extended emission is undetectable, or (2) the coupling efficiency between the AGN luminosity and the kinetic luminosity is less than predicted by Nims et al. (2015); $L_{\text{kin}} < 0.05L_{\text{AGN}}$).

For the first case, it is possible that the fiducial values for the ambient density, inferred from winds in nearby ULIRGs<fnref n="3" />

<footnote n="3">³ The fiducial density values are much lower than the measured mean molecular gas densities in the nuclear regions of ULIRGs (Downes & Solomon 1998) since winds are thought to propagate along the paths of least resistance (Faucher-Giguère & Quataert 2012).</footnote>
<figure><img src="image_11.png" /><figcaption>**Fig. 8.** Prediction for X-ray bubble luminosity as a function of shock-front radius for three values of the interstellar gas density at $R = 100$ pc ($n_{H,0}$) from the analytic model in Nims et al. (2015), assuming $n_{H} \propto R^{-1}$ (solid lines), compared to the upper limits from our observations (solid triangles). Different upper limits are shown because the predicted wind spectrum depends on density. The errors on the prediction come from bolometric correction uncertainty ($\kappa_{2-10\text{keV}} = 20_{-12}^{+80}$) when scaling to the bolometric luminosities of HE 0227-0913 and HE 0351+0240. Except for the lowest densities, our upper limits fall below the prediction.</figcaption></figure>

(Cicone et al. 2014; Feruglio et al. 2010), are very different from the host galaxies in this study. If the environmental properties of these ULIRGs turn out not to be representative of typical Type 1 AGN host galaxies, then our nondetection may suggest lower densities of at least an order of magnitude. For the second case, there have been a wide range of estimates of the coupling efficiency derived from observations, anywhere between 0.001% and 10% (Harrison et al. 2018). 5% has been a standard value in the literature, obtained first from hydrodynamic simulations as the value needed to recover the local $M_{BH}-\sigma$ relation (Di Matteo et al. 2005). This value also comes up naturally from theoretical considerations of momentum-driven nuclear outflows with speeds of ~0.1 c (King 2003; Harrison et al. 2018), and indeed, this may be appropriate for Mrk 1044, as Mallick et al. (2018) tentatively detected an ultra fast outflow (with a speed of ~0.1 c) from a deep XMM-Newton observation of its nucleus. However, several recent observations and simulations of AGN suggest much smaller coupling efficiencies (~0.5%; Husemann et al. 2016b; Veilleux et al. 2017; Richings & Faucher-Giguère 2018), which would be consistent with our results.

There have been a few previous cases of extended X-ray emission found in AGN possibly from AGN superwinds. Greene et al. (2014) found extended X-ray emission in SDSS J1356+1026, a ULIRG galaxy in an ongoing merger with a 20 kpc bipolar outflow. The non-nuclear emission was estimated to be $L_{0.3-8\text{keV}} = 7 \times 10^{41}$ erg s$^{-1}$, with densities of ~2.4 cm$^{-3}$ at the kpc scale and an AGN luminosity of $L_{\text{bol}} = 10^{46}$ erg s$^{-1}$. In another example, the teacup galaxy ($L_{\text{bol}} = 10^{45}-10^{46}$ erg s$^{-1}$) has X-ray emission ($L_{2-10\text{keV}} \sim 10^{44}$ erg s$^{-1}$), spatial with its ionized, 10 kpc bubble (Lansbury et al. 2018). The spectrum is consistent with a two-temperature ~0.4 keV and ≥2.7 keV thermal gas, with $L_{1-10\text{keV}} \sim 10^{41}$ erg s$^{-1}$. For both galaxies, we would have detected their diffuse hot gas were it present in our objects. However, the X-ray luminosities fall below the predicted values for their respective AGN luminosities. Mrk 231 ($L_{\text{bol}} = 10^{46.2}$) also has diffuse X-ray emission fainter than predicted, with $L_{2-10\text{keV}} \sim 10^{41}$ erg s$^{-1}$ within the first few kpc (Veilleux et al. 2014).

## 5.2. Other sources of extended X-ray emission

Extended X-ray emission can also come from hot halo gas, typically found in galaxy groups and clusters. The stellar masses of the two galaxies studied here are $\log(M_{\star}/M_{\odot}\text{h}^{-1}) \sim 10.1$ and 10.5, estimated based on the *g-i* color and *i* band absolute magnitude from the QSO-subtracted broadband MUSE images following the prescription of Taylor et al. (2011). This corresponds to halo masses of $\log(M_{\text{h}}/M_{\odot}\text{h}^{-1}) \sim 11.6$ and 12.0 assuming the stellar mass to halo mass relation of Behroozi et al. (2010). At these mass scales, the expected total X-ray luminosity from extended emission is $\le 10^{39}$ erg s$^{-1}$ (Anderson et al. 2015), which is well below the superbubble predictions on the largest scales, below our detection limit.

Star formation in galaxies is another source of X-rays, produced by the evolved X-ray binaries as well as supernovae. While this is a negligible contribution for HE 0351+0240 due to the lack of star forming regions, the star formation rate in the nuclear region of Mrk 1044 is ~1 $M_{\odot}$ yr$^{-1}$ within the inner 3" (inner kpc) based on the extinction-corrected $H\alpha$ luminosity following Kennicutt (1998). This SFR corresponds to an X-ray luminosity of $L_{0.5-8\text{keV}} \sim 4 \times 10^{39}$ erg s$^{-1}$ (Mineo et al. 2014). While this is below our upper limit (and the superbubble predictions for the mid and high density values), it could be contributing to the tentative detection we find in the first radial bin.

## 5.3. The origin of highly ionized gas arcs

The two AGN host galaxies were selected from CARS for our Chandra pilot observations because of the prominent arc-like structures of highly-ionized gas on kpc scales from the AGN. Based on the deep Chandra observations we argue that the arcs are most likely not related to the interplay between the hot medium of an AGN-driven outflow bubble and the cold ambient medium for these two galaxies. While hot gas AGN-driven bubbles may be important for other AGN hosts, it does not explain the existence of the prominent ionized gas arcs in HE 0227-0913 and HE 0351+0240. So what are the alternative
explanations for their appearance and shape in these two studied cases?

Photoionization of gas by the AGN on several kpc scales is usually referred to as the extended narrow-line region (ENLR) or the extended emission-line region (EELR), which appears to be a common feature around luminous radio-loud (e.g., Stockton & MacKenty 1987; Villar-Martin et al. 1997; Fu & Stockton 2009) and radio-quiet AGN (e.g., Bennert et al. 2002; Husemann et al. 2013; Liu et al. 2013; Hainline et al. 2013; Harrison et al. 2014; Sun et al. 2017). The size, geometry and shape of the ENLR potentially depends on several factors, such as the AGN luminosity (e.g., Husemann et al. 2014; Sun et al. 2017, 2018), the opening angle and orientation of the AGN ionization cones (e.g., Bennert et al. 2006), the gas density distribution intercepting the ionizing radiation field of the AGN (e.g., Dempsey & Zakamska 2018), or other intrinsic properties of the central AGN engine (Husemann et al. 2008; Keel et al. 2012). The intrinsic AGN luminosity is comparable in both galaxies studied and corresponds to expected ENLR sizes of 5–10 kpc depending on different ENLR size definitions and samples (Husemann et al. 2014; Sun et al. 2017). Since the morphologies of the host galaxies are rather different, we expect that the relation between the geometry of the AGN radiation field and the distribution of the gas plays the dominant role in shaping the appearance of the [O III] arcs, as we discuss below.

HE 0227–0931 shows a disk-dominated morphology, which is also evident in the rotation-dominated velocity field of the star-forming regions (see Fig. 3). The velocity of the [O III] arc clearly does not follow that rotational motion and is close to the systemic velocity of the galaxy. One possibility is that the gas was tidally stripped from a minor companion galaxy with a close orbit around HE 0227–0931. However, there is no signature of such a companion galaxy in the continuum light that could be spatially associated with the ionized gas arc. Unless the stellar component was tidally stripped as well, resulting in an undetectable low surface brightness companion, this scenario is unlikely. Furthermore, ionization by shocks can be excluded since the high [O III]/Hβ ~ 10 ratio can only be produced by fast shocks (>500 km s⁻¹) and their precursors (Allen et al. 2008), which are not supported by the quiescent kinematics of the gas. Instead, we noticed that the highest surface brightness regions of the [O III] arc appear at a similar radius to the star forming ring visible in the stellar continuum and Hα emission. Therefore, an alternative possibility is that the high star formation surface density in the ring has caused a gas outflow perpendicular to the galaxy disk, as predicted in the Galactic fountain picture (e.g., Shapiro & Field 1976; Fraternali & Binney 2006). Recent simulations show that the warm gas will likely fall back to the disk (e.g., Kim & Ostriker 2018). The illumination of clouds ejected from the galaxy and falling back seems a plausible scenario for HE 0227–0931, as the AGN radiation field can only significantly intercept gas above the disk at radii of a few kpc.

On the other hand, HE 0351+0240 is clearly an ongoing major merger system. During this process, the gas in the two parent galaxies have likely been significantly redistributed due to their tidal interactions. Gas shells and streams can be easily created on large scales, depending on the orbital configuration and angular momenta of interacting galaxies (e.g., Arp 1966; Hibbard & van Gorkom 1996), and predicted by numerical simulation of mergers (e.g., Toomre & Toomre 1972; Mihos & Hernquist 1996). If the AGN radiation field is intercepting those streams, they become visible as highly ionized gas as happens, for example, in the major merger NGC 7252 (Schweizer et al. 2013; Weaver et al. 2018). Since ionization by

shocks can be excluded with the same argument as for HE 0227–0913, we argue that the [O III] arcs of HE 0351+0240 are likely produced as the intersection of the large-scale distribution of tidal gas streams and the AGN radiation field. Assuming a static ionization cone, the opening angle would have to be ≥80°; however the central engine likely changes orientation with respect to the larger gas environment due to the disturbed kinematics of the merging system.

## 6. Summary and conclusions

In this study, we united MUSE IFU data with deep Chandra observations to search for quasar-driven superbubbles in two nearby AGN host galaxies from the CARS survey. Selected based on their quiescent kpc-scale emission arcs, both galaxies were promising candidates for catching the hot-wind mode of radiative AGN feedback in action. However, we found no evidence for extended, hot diffuse gas associated with a super wind. Our upper X-ray luminosity limits of ~10⁴⁰⁻⁴¹ erg s⁻¹ on kpc scales is well below the expectation (~10⁴² erg s⁻¹) based on standard theory (Nims et al. 2015). From this nondetection, we conclude the following:

- Based on theoretical predictions, either the coupling efficiency between the AGN and ambient medium must be small (<5%), or the density of the medium through which the wind travels is low ($n_H \le 1 \text{ cm}^{-3}$ at a distance of 100 pc).

- Photoionization from the AGN of clouds that were previously ejected from the galaxies is likely the origin of the ionized gas structures observed in each AGN.

- Hot winds are either not ubiquitous among moderate-luminosity AGN, or are not connected to the observed ionized kpc-scale structures that are common in such objects.

While AGN-driven winds may be present for the highest-luminosity AGN, their extent and their relation to the ionized gas component are still not well constrained. The physics of AGN outflows remains murky, and only multiphase ISM observations can reveal the full picture (Cicone et al. 2018; Husemann & Harrison 2018). For this purpose, future investigations with AGN from CARS are planned, utilizing spatially-resolved data across the full electromagnetic spectrum.

*Acknowledgements.* Support for this work was provided by the National Aeronautics and Space Administration through Chandra Award Number GO6-17093X issued by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space Administration under contract NAS8-03060. GRT acknowledges support from NASA through Chandra Award Number GO7-8128X as well as Einstein Postdoctoral Fellowship Award Number PF-150128, issued by the Chandra X-ray Observatory Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of NASA under contract NAS8-03060. MK acknowledges support by DFG grant KR 3338/3-1, and TAD acknowledges support from a Science and Technology Facilities Council Ernest Rutherford Fellowship. MPT acknowledges support by the Spanish MINECO through grants AYA2012-38491-C02-02 and AYA2015-63939-C2-1-P, co-funded with FEDER funds.

## References

- Allen, M. G., Groves, B. A., Dopita, M. A., Sutherland, R. S., & Kewley, L. J. 2008, ApJS, 178, 20

- Anderson, M. E., Gaspari, M., White, S. D. M., Wang, W., & Dai, X. 2015, MNRAS, 449, 3806

- Arp, H. 1966, ApJS, 14, 1

- Bacon, R., Accardo, M., Adjali, L., et al. 2010, SPIE Conf. Ser., 7735, 8

- Bacon, R., Vernet, J., Borisova, E., et al. 2014, The Messenger, 157, 13

- Bacon, R., Conseil, S., Mary, D., et al. 2017, A&A, 608, A1

- Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5
Behroozi, P. S., Conroy, C., & Wechsler, R. H. 2010, ApJ, 717, 379

Bennert, N., Falcke, H., Schulz, H., Wilson, A. S., & Wills, B. J. 2002, ApJ, 574, L105

Bennert, N., Jungwiert, B., Komossa, S., Haas, M., & Chini, R. 2006, New Astron. Rev., 50, 708

Bertram, T., Eckart, A., Fischer, S., et al. 2007, A&A, 470, 571

Bower, R. G., Benson, A. J., Malbon, R., et al. 2006, MNRAS, 370, 645

Cappellari, M., & Copin, Y. 2003, MNRAS, 342, 345

Carter, C., Karovska, M., Jerius, D., Glotfelty, K., & Beikman, S. 2003. in Astronomical Data Analysis Software and Systems XII, eds. H. E. Payne, R. I. Jedrzejewski, & R. N. Hook, ASP Conf. Ser., 295, 477

Castelló-Mor, N., Netzer, H., & Kaspi, S. 2016, MNRAS, 458, 1839

Cicone, C., Maiolino, R., Sturm, E., et al. 2014, A&A, 562, A21

Cicone, C., Brusa, M., Ramos Almeida, C., et al. 2018, Nat. Astron., 2, 176

Davis, T. A., Alatalo, K., Bureau, M., et al. 2013, MNRAS, 429, 534

Dempsey, R., & Zakamska, N. L. 2018, MNRAS, 477, 4615

Dewangan, G. C., Griffiths, R. E., Dasgupta, S., & Rao, A. R. 2007, ApJ, 671, 1284

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604

Downes, D., & Solomon, P. M. 1998, ApJ, 507, 615

Faucher-Giguère, C.-A., & Quataert, E. 2012, MNRAS, 425, 605

Feruglio, C., Maiolino, R., Piconcelli, E., et al. 2010, A&A, 518, L155

Fraternali, F., & Binney, J. J. 2006, MNRAS, 366, 449

Fruchter, A. S., & Hook, R. N. 2002, PASP, 114, 144

Fruscione, A., McDowell, J. C., Allen, G. E., et al. 2006, Proc. SPIE, 6270, 62701V

Fu, H., & Stockton, A. 2009, ApJ, 690, 953

Greene, J. E., Pooley, D., Zakamska, N. L., Comerford, J. M., & Sun, A.-L. 2014, ApJ, 788, 54

Guérou, A., Krajnović, D., Epinat, B., et al. 2017, A&A, 608, A5

Hainline, K. N., Hickox, R., Greene, J. E., Myers, A. D., & Zakamska, N. L. 2013, ApJ, 774, 145

Harrison, C. M., Alexander, D. M., Mullaney, J. R., & Swinbank, A. M. 2014, MNRAS, 441, 3306

Harrison, C. M., Costa, T., Tadhunter, C. N., et al. 2018, Nat. Astron., 2, 198

Heckman, T. M., Armus, L., & Miley, G. K. 1990, ApJS, 74, 833

Hibbard, J. E., & van Gorkom, J. H. 1996, AJ, 111, 655

Husemann, B., & Harrison, C. M. 2018, Nat. Astron., 2, 196

Husemann, B., Wisotzki, L., Sánchez, S. F., & Jahnke, K. 2008, A&A, 488, 145

Husemann, B., Wisotzki, L., Sánchez, S. F., & Jahnke, K. 2013, A&A, 549, A43

Husemann, B., Jahnke, K., Sánchez, S. F., et al. 2014, MNRAS, 443, 755

Husemann, B., Bennert, V. N., Scharwächter, J., Woo, J.-H., & Choudhury, O. S. 2016a, MNRAS, 455, 1905

Husemann, B., Scharwächter, J., Bennert, V. N., et al. 2016b, A&A, 594, A44

Husemann, B., Tremblay, G., Davis, T., et al. 2017, The Messenger, 169, 42

Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003, MNRAS, 346, 1055

Keel, W. C., Chojnowski, S. D., Bennert, V. N., et al. 2012, MNRAS, 420, 878

Kennicutt, R. C., Jr 1998, ApJ, 498, 541

Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & Trevena, J. 2001, ApJ, 556, 121

Kim, C.-G., & Ostriker, E. C. 2018, ApJ, 853, 173

King, A. 2003, ApJ, 596, L27

König, S., Eckart, A., García-Marín, M., & Huchtmeier, W. K. 2009, A&A, 507, 757

Koss, M., Trakhtenbrot, B., Ricci, C., et al. 2017, ApJ, 850, 74

Lansbury, G. B., Jarvis, M. E., Harrison, C. M., et al. 2018, ApJ, 856, L1

Liu, G., Zakamska, N. L., Greene, J. E., Nesvadba, N. P. H., & Liu, X. 2013, MNRAS, 430, 2327

Lusso, E., Comastri, A., Simmons, B. D., et al. 2012, MNRAS, 425, 623

Mallick, L., Alston, W. N., Parker, M. L., et al. 2018, MNRAS, 479, 615

Mihos, J. C., & Hernquist, L. 1996, ApJ, 464, 641

Mineo, S., Gilfanov, M., Lehmer, B. D., Morrison, G. E., & Sunyaev, R. 2014, MNRAS, 437, 1698

Nims, J., Quataert, E., & Faucher-Giguère, C.-A. 2015, MNRAS, 447, 3612

Richings, A. J., & Faucher-Giguère, C.-A. 2018, MNRAS, 478, 3100

Sartori, L. F., Schawinski, K., Koss, M., et al. 2016, MNRAS, 457, 3629

Schaye, J., Crain, R. A., Bower, R. G., et al. 2015, MNRAS, 446, 521

Schweizer, F., Seitzer, P., Kelson, D. D., Villanueva, E. V., & Walth, G. L. 2013, ApJ, 773, 148

Shapiro, P. R., & Field, G. B. 1976, ApJ, 205, 762

Silk, J., & Rees, M. J. 1998, A&A, 331, L1

Soto, K. T., Lilly, S. J., Bacon, R., Richard, J., & Conseil, S. 2016, MNRAS, 458, 3210

Stasińska, G., Vale Asari, N., Cid Fernandes, R., et al. 2008, MNRAS, 391, L29

Stockton, A., & MacKenty, J. W. 1987, ApJ, 316, 584

Storey, P. J., & Zeppen, C. J. 2000, MNRAS, 312, 813

Sun, A.-L., Greene, J. E., & Zakamska, N. L. 2017, ApJ, 835, 222

Sun, A.-L., Greene, J. E., Zakamska, N. L., et al. 2018, MNRAS, 480, 2302

Taylor, E. N., Hopkins, A. M., Baldry, I. K., et al. 2011, MNRAS, 418, 1587

Toomre, A., & Toomre, J. 1972, ApJ, 178, 623

Valdes, F., Gupta, R., Rose, J. A., Singh, H. P., & Bell, D. J. 2004, ApJS, 152, 251

Veilleux, S., & Osterbrock, D. E. 1987, ApJS, 63, 295

Veilleux, S., Teng, S. H., Rupke, D. S. N., Maiolino, R., & Sturm, E. 2014, ApJ, 790, 116

Veilleux, S., Bolatto, A., Tombesi, F., et al. 2017, ApJ, 843, 18

Villar-Martin, M., Tadhunter, C., & Clark, N. 1997, A&A, 323, 21

Walker, A.R., Boccas, M., Bonati, M., 2003, in Instrument Design and Performance for Optical/Infrared Ground-based Telescopes, eds. M. Iye, & A. F. M. Moorwood, Proc. SPIE, 4841, 286

Wang, J.-M., Du, P., Hu, C., et al. 2014, ApJ, 793, 108

Weaver, J. R., Husemann, B., Kuntschner, H., et al. 2018, A&A, 614, A32

Weilbacher, P. M., Streicher, O., Urrutia, T., et al. 2012, SPIE Conf. Ser., 8451, 84510B

Weinberger, R., Springel, V., Pakmor, R., et al. 2018, MNRAS, 479, 4056

Wisotzki, L., Christlieb, N., Bade, N., et al. 2000, A&A, 358, 77

Zubovas, K., & King, A. 2012, ApJ, 745, L34