# The Close AGN Reference Survey (CARS): Mrk 1018 halts dimming and experiences strong short-term variability

M. Krumpe, B. Husemann, G.R. Tremblay, T. Urrutia, M. Powell, T.A.

Davis, J. Scharwächter, J. Dexter, G. Busch, F. Combes, et al.

► To cite this version:

M. Krumpe, B. Husemann, G.R. Tremblay, T. Urrutia, M. Powell, et al. The Close AGN Reference Survey (CARS): Mrk 1018 halts dimming and experiences strong short-term variability. Astronomy & Astrophysics - A&A, 2017, 607, pp.L9. 10.1051/0004-6361/201731967. [hal-01669509](https://hal.science/hal-01669509v1)

HAL Id: hal-01669509

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

Submitted on 31 Oct 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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
LETTER TO THE EDITOR

# The Close AGN Reference Survey (CARS)

## Mrk 1018 halts dimming and experiences strong short-term variability*

M. Krumpe<fnref n="1" />, B. Husemann<fnref n="2" />, G. R. Tremblay<fnref n="3" />,<fnref n="4" />, T. Urrutia<fnref n="1" />, M. Powell<fnref n="4" />, T. A. Davis<fnref n="5" />, J. Scharwächter<fnref n="6" />, J. Dexter<fnref n="7" />, G. Busch<fnref n="8" />, F. Combes<fnref n="9" />, S. M. Croom<fnref n="10" />,<fnref n="11" />, A. Eckart<fnref n="8" />,<fnref n="12" />, R. E. McElroy<fnref n="10" />,<fnref n="11" />, M. Perez-Torres<fnref n="13" />, and G. Leung<fnref n="14" />

Received 18 September 2017 / Accepted 24 October 2017

## ABSTRACT

After changing optical AGN type from 1.9 to 1 in 1984, the AGN Mrk 1018 recently reverted back to its type 1.9 state. Our ongoing monitoring now reveals that the AGN has halted its dramatic dimming, reaching a minimum around October 2016. The minimum was followed by an outburst rising with ~0.25 *U*-band mag/month. The rebrightening lasted at least until February 2017, as confirmed by joint *Chandra* and *Hubble* observations. Monitoring was resumed in July 2017 after the source emerged from sunblock, at which point the AGN was found only ~0.4 mag brighter than its minimum. The intermittent outburst was accompanied by the appearance of a red wing asymmetry in broad-line shape, indicative of an inhomogeneous broad-line region. The current flickering brightness of Mrk 1018 following its rapid fading either suggests that the source has reignited, remains variable at a low level, or may continue dimming over the next few years. Distinguishing between these possibilities requires continuous multiwavelength monitoring.

**Key words.** quasars: individual: Mrk 1018 – accretion, accretion disks – galaxies: evolution

## 1. Introduction

During a galaxy's lifetime, various physical processes, many of which remain poorly understood, can either trigger or shut down a mass flow onto the central supermassive black hole (SMBH). Theoretical models (e.g., Hopkins et al. 2008) predict that an SMBH accretion reservoir can be consumed in ~$10^{4-5}$ years, such that the AGN transitions to a “normal” galaxy as the central engine shuts down. The AGN phase is likely episodic with a wide range of accretion ratios from one duty cycle to the next (e.g., Novak et al. 2011). In this scenario, AGN are better described as events than as a discrete class of objects.

The term changing-look (CL) AGN originally referred to objects converting optical spectral type from Type 1 to Type 2 or vice versa. More recently, the definition was relaxed to AGN showing significant changes in the spectral features. It is not clear whether these objects are the extremes in a distribution of

AGN variability or caused by some major discrete event in or around the central engine. A clumpy accretion disk, that is, variable gas input, would naturally explain the absence of constant or only slowly varying accretion. The corresponding timescales would be shorter than the rotation period in the accretion disk, which is on the order of ~30 years or shorter (depending on mass and distance from the SMBH). Accretion disks can also be warped, precessing, or lopsided. In the past, CL AGN have been explained either by variable obscuration along the line of sight, a tidal disruption event (TDE, Rees 1988), or simply a drastic change in the intrinsic accretion rate. Understanding CL AGN therefore offers an opportunity to improve our understanding of the AGN central engine and accretion process around SMBHs.

Only a few cases of the original definition of CL AGN are known, such as NGC 7603 (Tohline & Osterbrock 1976), NGC 7582 (Aretxaga et al. 1999), NGC 2617 (Shappee et al. 2014; Bon et al. 2016), Mrk 590 (Denney et al. 2014), HE 1135-2304 (Parker et al. 2016), and SDSS J015957.64+003310.5 (LaMassa et al. 2015; Merloni et al. 2015). Significant spectral changes are found for instance in SDSS J1011+5442 (Runnoe et al. 2016), J01264-0839 and J2336+00172 (Ruan et al. 2016), and IPTF 16bco (Gezari et al. 2017). New monitoring

<footnote n="1">¹ Leibniz-Institut für Astrophysik Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany
e-mail: mkrumpe@aip.de</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">⁴ Yale Center for Astronomy and Astrophysics, Yale University, 52 Hillhouse Ave, New Haven, CT 06511, USA</footnote>

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

<footnote n="6">⁶ Gemini Observatory, Northern Operations Center, 670 N. A'ohoku Pl., Hilo, HI, 96720, USA</footnote>

<footnote n="7">⁷ Max-Planck-Institut für extraterrestische Physik, Giessenbachstr. 1, 85748 Garching, Germany</footnote>

<footnote n="8">⁸ I. Physikalisches Institut, Universität zu Köln, Zülpicher Straße 77, 50937 Köln, Germany</footnote>

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

<footnote n="10">¹⁰ Sydney Institute for Astronomy, School of Physics, University of Sydney, NSW 2006, Australia</footnote>

<footnote n="11">¹¹ ARC Centre of Excellence for All-sky Astrophysics (CAASTRO), Australia</footnote>

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

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

<footnote n="14">¹⁴ Center for Astrophysics and Space Sciences, University of California, San Diego, 9500 Gilman Dr., La Jolla, CA 92093, USA</footnote>

<footnote>* Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme(s) 098.B-0672 and 099.B-0159. The scientific results reported in this article are based on observations made by the *Chandra* X-ray Observatory and the NASA/ESA *Hubble* Space Telescope.</footnote>
projects now actively hunt for more of these unique objects (MacLeod et al. 2016).

The rarity of the changing-look behavior in AGN raises many questions: are these events simply rare and unique, or common but short lived? Why and how does the AGN shut down, and how long does the shut-down process last? During shut-down, does the accretion disk disappear first or the X-ray corona that envelopes it? How fast can an AGN be reignited after it was shut down?

As part of the Close AGN Reference Survey (CARS<fnref n="1" />), we have recently discovered a highly unique AGN in which to investigate these questions: Mrk 1018. The host galaxy of Mrk 1018 is an ongoing major merger and is one of the first known CL AGN, which transitioned from a Seyfert 1.9 to a Seyfert 1 nucleus around 1984 (Cohen et al. 1986). The stellar mass of Mrk 1018 is $log M_*/M_\odot = 10.92$ (Koss et al. 2011) and its star formation 0.6 $M_\odot$ yr⁻<fnref n="1" /> (Shimizu et al. 2017). Bennert et al. (2011) found $log (M_{BH}/M_\odot) = 8.15$. The current $L_{bol}/L_{edd}$ is around 0.03 (McElroy et al. 2016). We note that because the broad-line region (BLR) of Mrk 1018 is currently not in an equilibrium or the virial factor has changed, recent $M_{BH}$ estimates are significantly different (e.g., current faint phase: $log (M_{BH}/M_\odot) = 7.4$; McElroy et al. 2016).

McElroy et al. (2016) and Husemann et al. (2016) reported that the nucleus significantly dimmed and changed type back to a Seyfert 1.9 between 2013–2015. Hence, Mrk 1018 is the first known CL AGN to undergo a full state change cycle, which we caught and monitored in the act of shutting down. In this letter, we present recent optical, UV, and X-ray monitoring data of Mrk 1018 that reveal that the light curve has not continued dimming. Throughout, we adopt $H_0 = 70$ km s⁻<fnref n="1" /> Mpc⁻<fnref n="1" />, $\Omega_m = 0.3$ and $\Omega_\Lambda = 0.7$ cosmological parameters, and the AB magnitude system. Uncertainties represent $1\sigma$ (68.3%) confidence intervals unless otherwise stated.

## 2. Optical monitoring

### 2.1. Dimming of Mrk 1018

Shortly after the discovery of the optical dimming by a factor of ~25 in mid-February 2016, follow-up observations with *Chandra* and HST were taken. The combined *Chandra* and archival *NuStar* spectra from February 2016 revealed that Mrk 1018 still shows no detectable absorption in the X-ray spectrum (Husemann et al. 2016), consistent with the bright Seyfert 1 phase (archival *Chandra* spectrum from 2010; ID: 12868, PI: Mushotzky). Hence, the current changing-look event cannot be caused by an obscuration event, but must instead be due to a significant drop in disk luminosity. LaMassa et al. (2017) reanalyzed the *Chandra* data from 2010 and 2016 as well as the *NuStar* data. They confirmed that Mrk 1018 shows no line-of-sight X-ray absorption and concluded that the variation in the Fe $\mathrm{K\alpha}$ equivalent width between 2010 and 2016 is due to a decreasing continuum flux and not an increase in the line flux.

The significant drop in intrinsic disk luminosity is also confirmed by the decreasing luminosity at mid-IR wavelength as observed with WISE (Sheng et al. 2017). The most recent NEOWISE data also indicate that the dimming stopped even in the mid-IR. However, no detailed quantification can be given because of the low spatial resolution and the substantial host galaxy contamination. This suggests that the circumnuclear dust, usually referred to as the torus, is exposed to a significantly lower incident radiation field as well. Since the brightness of

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

**Fig. 1.** U-band light curve of the AGN nucleus between June 2016 and October 2017 after subtraction of the host galaxy contribution. The inset shows the entire U-band light curve since 2000 for comparison as presented in (McElroy et al. 2016). The gray area in the inset highlights the coverage of the zoom-in axes. The substantial sunblock period in 2017, where it is unobservable, is indicated by the black hatched area. The date of our joint *Chandra* and HST observation is shown as the green dashed line.

Mrk 1018 appeared to be stable for several decades and the dimming does not follow an expected $t^{-5/3}$ law, a TDE is unlikely. Two possible scenarios to explain the significant luminosity variations of Mrk 1018 on a timescale of a few years were discussed in Husemann et al. (2016): 1) a binary SMBH interaction; or 2) self-regulated accretion through accretion disk winds.

### 2.2. Reaching a minimum in October 2016

In response to our discovery that Mrk 1018 has faded by an order of magnitude within just a few years, we started an optical monitoring program with the VLT-VIMOS instrument (Le Févre et al. 2003). Deep U-band images with 600 s total integration time each were obtained between June 2016 and January 2017 and from July 2017 until now with an average cadence of about two weeks. The images were reduced with the ESO VIMOS pipeline and accurate relative photometric calibration was achieved with a reference star 60'' away from the nucleus.

We performed a 2D imaging fitting with *galfit* (Peng et al. 2010) of all images to decompose the light of the AGN nucleus and the host galaxy. After masking out the tidal features of this major merger system, the main body of the host galaxy can be reasonably well described by a single Séricic profile. We fixed the parameters of the model to $n = 3.9$, $r_c = 14.37''$, $b/a = 0.85$ and $PA = -85°$. These values are guided by a best-fit model of the deep coadded SDSS Stripe 82 model that we used previously to generate the long-term light curve in McElroy et al. (2016). The resulting U-band light curve of the AGN component is shown in Fig. 1. The U-band magnitude of the host galaxy component, which is the only other free parameter left for the 2D fitting, is 16.09 mag with an rms error of 0.01 mag. We adopt a $1\sigma$ error of 0.05 mag for the AGN brightness to account for systematic uncertainties of the AGN-host galaxy deblending. A minimum brightness level of $m_{U,AGN} \sim 19.3$ mag was reached in October 2016.

### 2.3. Outburst of Mrk 1018

Unexpectedly, the nucleus began to brighten again beginning in October 2016, by ~0.25 mag/month. The rebrightening lasted at least until February 2017, when our joint *Chandra* and HST observations were executed (see Sects. 3 and 4) just before this equatorial source passed into sunblock. Monitoring was resumed in July 2017, revealing a variable brightness of ±0.1 mag around $m_{U,AGN} \sim 18.9$ mag. The current brightness level is still ~0.4 mag higher than during the minimum state in October 2016.

<footnote n="1">¹ www.cars-survey.org</footnote>
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Comparison of the X-ray spectra of Mrk 1018 at three epochs (2010, 2016, and 2017). For all data, we show the best-fit model. These observed spectra are not unfolded by the telescope response.</figcaption></figure>

<table><caption>**Table 1.** Best-fit model parameters of *Chandra* spectra.</caption><thead><tr><th>Epoch</th><th>Γ</th><th>f<sub>0.5-2 keV</sub><br>[10<sup>-12</sup> erg cm<sup>-2</sup> s<sup>-1</sup>]</th><th>f<sub>2-10 keV</sub></th></tr></thead><tbody><tr><td>2010 Nov.</td><td>1.70 ± 0.03</td><td>4.96 ± 0.15</td><td>9.26 ± 0.19</td></tr><tr><td>2016 Feb.</td><td>1.62 ± 0.03</td><td>0.58 ± 0.02</td><td>1.25 ± 0.03</td></tr><tr><td>2017 Feb.</td><td>1.63 ± 0.02</td><td>1.12 ± 0.02</td><td>2.32 ± 0.03</td></tr></tbody></table>

**Notes.** The intrinsic absorption ($N_H$) is consistent with zero at all epochs. Fluxes are corrected for Galactic absorption.

(Fig. 1). We have confirmed that the host galaxy brightness remains constant within a nominal uncertainty. The outburst was apparently very short with an unknown peak brightness.

### 3. X-ray view of Mrk 1018 during the outburst

The latest *Chandra* observation of Mrk 1018 was taken on 2017-02-17 (ID: 19560, PI: Krumpe) with ACIS-S for ~50 ks. The spectra were extracted with the CIAO package (4.9) and the latest CALDB files (4.7.3) using standard setting for point sources. Since pile-up only affects the central pixel with a likelihood of ~1%, we ignored the effects of pile-up and created a standard source and background spectrum. The spectrum was grouped with a minimum binning of 20 counts.

We fit the X-ray spectra using Xspec version 12.9.1e (Arnaud 1996) in the 0.5–8 keV energy range. We used the cosmic abundances of Wilms et al. (2000) and the photoelectric absorption cross sections provided by Verner et al. (1996). We fit a model consisting of an intrinsically absorbed power law, a narrow Gaussian line profile fixed at a rest-frame energy of 6.4 keV (width $\sigma$ = 0.1 keV), and absorption by Galactic neutral hydrogen ($N_{HI, Gal}$ = 2.43 × 10²⁰ cm⁻², Kalberla et al. 2005). The 90% upper limit on the intrinsic $N_H$ is 1 × 10²⁰ cm⁻², consistent with no neutral absorption along the line of sight in excess over the Galactic column. We tentatively detect a narrow Fe Kα line with $2.6\sigma$ (EW = 0.14±0.05 keV). The best-fit photon index (slope of the X-ray spectrum) is $\Gamma$ = 1.63 ± 0.02 ($\Delta\chi^2$/d.o.f. = 298/271, freezing the intrinsic $N_H$ = 0). This model yields Galactic absorption-corrected fluxes of $f_{0.5-2 keV}$ = 1.1 × 10⁻¹² erg cm⁻² s⁻¹ and $f_{2-10 keV}$ = 2.3 × 10⁻¹² erg cm⁻² s⁻¹, which corresponds to rest-frame luminosities of $L_{0.5-2 keV}$ = 4.6 × 10⁴² erg s⁻¹ and $L_{2-10 keV}$ = 9.5 × 10⁴² erg s⁻¹.

We reanalyzed the 2010 and 2016 *Chandra* data in the same manner as described above to minimize systematic uncertainties. Best-fit model parameters for the three *Chandra* epochs are given in Table 1. A direct comparison of all Mrk 1018 data obtained with *Chandra* is shown in Fig. 2. The photon indices in 2016 and 2017 are consistent with each other. Owning to the limited signal-to-noise ratio in all *Chandra* observations, the decrease in photon index in 2016 and 2017 compared to the 2010 observation is not significant (~2σ) when considering the uncertainties.

<caption>**Fig. 3.** Top: FUV spectra of Mrk 1018 from 1984 until 2017. Bottom: ratio of the 2017 and 2016 HST/COS spectra, with a 357-day temporal baseline. Median filtering is applied to suppress the continuum noise. The horizontal dashed line highlights a continuum scale factor of 1.5.</caption>

The data from 2017 can be equally well fit ($\Delta\chi^2$/d.o.f. = 298/272) when using the above-mentioned model tbabs*(zpowerlaw+zgauss) frozen to the 2016 best-fit model and only allowing the 2017 spectrum to vary globally in normalization, as well as allowing the Fe line flux to vary freely. The best fit reveals a scaling factor of 1.92 ± 0.02. No obvious systematic difference between fitting the data with the 2016 best-fit model or a free fit are found. The fact that the broad-band X-ray spectrum only shows a scaling in flux but not in shape is also verified when comparing the change in 0.5–2 and 2–10 keV fluxes between the two observations.

### 4. Broad emission-line variations

Far-UV (FUV) spectroscopy with HST/COS using the G140L grism was obtained quasi-simultaneously with the *Chandra* observation on 2017 February 18 (ID: 14853, PI: M. Krumpe) over two orbits. The data were processed with the standard COS pipeline, and the calibrated spectrum is shown in Fig. 3 (upper panel) compared to several archival FUV spectra. The FUV continuum flux increased by a factor of ~1.5 between 2016 and 2017 (357 days), which is similar to but slightly lower than the increase in X-ray luminosity during the outburst. Interestingly, the broad $Ly\alpha$ and C IV lines show prominent shape changes between 2016 and 2017. We highlight the changes in the lower panel of Fig. 3 by dividing the spectra. While the central $Ly\alpha$ line core dimmed, the red wing of the $Ly\alpha$ and C IV lines significantly increased in flux.

We have also undertaken spectroscopic monitoring of the broad $H\alpha$ line with VIMOS in IFU mode using the high-resolution red grism over a FoV of 27′′ × 27′′. A reference star 60′′ away was observed twice together with Mrk 1018, enabling accurate relative spectrophotometry. The VIMOS IFU data were reduced with our own Py3D IFU data reduction package (Husemann et al. 2013, 2014). In Fig. 4 we compare the spectrum observed with MUSE in 2015 with three representative spectra obtained around the minimum in 2016, the last one in January 2017 and the first spectrum after sun block in July 2017. In contrast to the FUV spectrum, which entirely scales with AGN luminosity, the optical spectra include additive components from the host galaxy and narrow emission-line regions. Thus, Fig. 3 uses the ratio between spectra and Fig. 4 the flux difference for comparison.

Despite some narrow-line residuals, the spectra reveal that the blue $H\alpha$ line asymmetry seen with MUSE in 2015 has faded away. This blue asymmetry was not present during the bright
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Top: evolution of the broad Hα line of Mrk 1018. Aperture spectra (3" diameter) from MUSE and three representative VIMOS IFU observations. Matching offset stellar continua model are shown. Bottom: difference spectra with respect to the minimum spectrum in October 2016. A constant stellar continuum is subtracted in each spectrum, but offset in flux to take the varying AGN continuum emission into account.</figcaption></figure>

Type 1 AGN phase, and has likely appeared now because the red wing of Hα disappeared amid the dimming phase, as discussed in McElroy et al. (2016). More importantly, we confirm that the prominent red wing asymmetry in Lyα and C IV is not present in Hα line only one month before the HST observations were taken. This confirms that the extreme red wing asymmetry in the FUV lines is directly associated with the outburst. It corresponds either to an asymmetry in the BLR with a distance of less than four light weeks from the SMBH, or it might correspond to a distinct accreting object with its own BLR.

## 5. Discussion and summary

Unexpectedly, the nucleus of Mrk 1018 reached a minimum brightness in October 2016 followed by a short outburst. The nucleus is currently $\Delta U \sim 0.4$ mag brighter than during the minimum a few months ago, which indicates that Mrk 1018 interrupted or ended its extraordinary dimming phase. If the global dimming really has ended, the minimum brightness of Mrk 1018 was only an order of magnitude fainter than in the bright phase. This would be a substantially different behavior than observed for CL AGN Mrk 590 (Denney et al. 2014), which dimmed by a factor of $\sim 100$ in a few years after a three-decade-long bright Type 1 AGN phase. This raises the question as to whether the same physical processes are causing the rapid disk luminosity drop in the CL AGN Mrk 1018 and Mrk 590.

The evolution in the shape of the broad emission lines encodes crucial information on the structure of the BLR around the AGN. Many AGN display significant changes in their broad-line shapes, including CL AGN, for example, NGC 5548 (Bon et al. 2016) or NGC 2617 (Oknyansky et al. 2017). In the case of Mrk 1018, the blue wing asymmetry in the broad-line shape during the global dimming phase (McElroy et al. 2016) faded away after 2015, whereas a prominent red wing asymmetry appeared less than four weeks before the brief outburst. These line-shape variations may be related to asymmetries in the BLR distribution, to accretion disk winds, or even to a close binary AGN that varies independently in accretion disk luminosity and BLR response.

Between the bright Seyfert 1 phase and 2016, the UV and X-ray flux dropped by a factor of $\sim 17$ and $\sim 8$, respectively (Husemann et al. 2016), but these observations were taken several years apart in the bright phase. The quasi-simultaneous UV and X-ray observations from 2016 and 2017 showed that the flux increases during the apparent outburst in 2017 shared similar

amplitudes (factors of 1.5 and 1.9, respectively). Thus, the X-ray corona and accretion disk are responding in a nearly identical fashion. Both regions must be within a light travel time of shorter than or around a year, and perhaps as short as a few months, given the short outburst time as seen in the U band.

The details of the underlying physical processes driving the apparently complex evolution of Mrk 1018 cannot be constrained with the currently available data. A better understanding of the central engine of Mrk 1018 now awaits reverberation mapping of the BLR over a far longer temporal baseline, combined with multiwavelength monitoring of the accretion disk emission and its interaction with the X-ray corona. That the dimming has possibly halted enables continuous monitoring of the AGN with a high signal-to-noise ratio, at least for the time being. The accretion disk of Mrk 1018 is nevertheless likely in an unstable state, and may brighten or fade with little warning. We therefore expect that ongoing radio through X-ray monitoring of Mrk 1018 will yield exciting discoveries over the coming years.

*Acknowledgements.* M.K. acknowledges support from DFG grant KR 3338/3-1 and DLR grant 50OR1802. G.R.T. acknowledges support from NASA through the 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. M.A.P.T. acknowledges support from the Spanish MINECO through grants AYA2012-38491-C02-02 and AYA2015-63939-C2-1-P. T.A.D. acknowledges support from a Science and Technology Facilities Council Ernest Rutherford Fellowship. Parts of this research were conducted by the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO), through project number CE110001020. G.L. acknowledges support provided by the National Aeronautics and Space Administration through *Chandra* Award Number G07-18111X issued by the *Chandra* X-ray Observatory Center (NAS8-03060). We thank Alison Coil for helpful discussions and suggestions.

## References

Aretxaga, I., Joguet, B., Kunth, D., Melnick, J., & Terlevich, R. J. 1999, ApJ, 519, L123

Arnaud, K. A. 1996, in Astronomical Data Analysis Software and Systems V, 101, 17

Bennert, V. N., Auger, M. W., Treu, T., Woo, J.-H., & Malkan, M. A. 2011, ApJ, 742, 107

Bon, E., Zucker, S., Netzer, H., et al. 2016, ApJS, 225, 29

Cohen, R. D., Puetter, R. C., Rudy, R. J., Ake, T. B., & Foltz, C. B. 1986, ApJ, 311, 135

Denney, K. D., De Rosa, G., Croxall, K., et al. 2014, ApJ, 796, 134

Gezari, S., Hung, T., Cenko, S. B., et al. 2017, ApJ, 835, 144

Hopkins, P. F., Hernquist, L., Cox, T. J., & Kereš, D. 2008, ApJS, 175, 356

Husemann, B., Jahnke, K., Sánchez, S. F., et al. 2013, A&A, 549, A87

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

Husemann, B., Urrutia, T., Tremblay, G. R., et al. 2016, A&A, 593, L9

Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775

Koss, M., Mushotzky, R., Veilleux, S., et al. 2011, ApJ, 739, 57

LaMassa, S. M., Cales, S., Moran, E. C., et al. 2015, ApJ, 800, 144

LaMassa, S. M., Yaqoob, T., & Kilgard, R. 2017, ApJ, 840, 11

Le Févre, O., Saisse, M., Mancini, D., et al. 2003, in Proc. SPIE, eds. M. Iye, & A. F. M. Moorwood, 4841, 1671

MacLeod, C. L., Ross, N. P., Lawrence, A., et al. 2016, MNRAS, 457, 389

McElroy, R., Husemann, B., Croom, S. M., et al. 2016, A&A, 593, L8

Merloni, A., Dwelly, T., Salvato, M., et al. 2015, MNRAS, 452, 69

Novak, G. S., Ostriker, J. P., & Ciotti, L. 2011, ApJ, 737, 26

Oknyansky, V. L., Gaskell, C. M., Huseynov, N. A., et al. 2017, MNRAS, 467, 1496

Parker, M. L., Komossa, S., Kollatschny, W., et al. 2016, MNRAS, 461, 1927

Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H. 2010, AJ, 139, 2097

Rees, M. J. 1988, Nature, 333, 523

Ruan, J. J., Anderson, S. F., Cales, S. L., et al. 2016, ApJ, 826, 188

Runnoe, J. C., Cales, S., Ruan, J. J., et al. 2016, MNRAS, 455, 1691

Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48

Sheng, Z., Wang, T., Jiang, N., et al. 2017, ApJ, 846, L7

Shimizu, T. T., Mushotzky, R. F., Meléndez, M., et al. 2017, MNRAS, 466, 3161

Tohline, J. E., & Osterbrock, D. E. 1976, ApJ, 210, L117

Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 465, 487

Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914