A publishing partnership

The Changing-look Quasar Mrk 590 Is Awakening

, , , , , , , , , , , and

Published 2018 October 19 © 2018. The American Astronomical Society. All rights reserved.
, , Citation S. Mathur et al 2018 ApJ 866 123 DOI 10.3847/1538-4357/aadd91

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0004-637X/866/2/123

Abstract

Mrk 590 was originally classified as a Seyfert 1 galaxy, but then it underwent dramatic changes: the nuclear luminosity dropped by over two orders of magnitude and the broad emission lines all but disappeared from the optical spectrum. Here we present follow-up observations to the original discovery and characterization of this "changing-look" active galactic nucleus (AGN). The new Chandra  and Hubble Space Telescope observations from 2014 show that Mrk 590 is awakening, changing its appearance again. While the source continues to be in a low state, its soft excess has re-emerged, though not to the previous level. The UV continuum is brighter by more than a factor of two and the broad Mg ii emission line is present, indicating that the ionizing continuum is also brightening. These observations suggest that the soft excess is not due to reprocessed hard X-ray emission. Instead, it is connected to the UV continuum through warm Comptonization. Variability of the Fe Kα emission lines suggests that the reprocessing region is within ∼10 lt-yr or 3 pc of the central source. The change in AGN type is neither due to obscuration nor due to one-way evolution from Type 1 to Type 2, as suggested in the literature, but may be related to episodic accretion events.

Export citation and abstract BibTeX RIS

1. Introduction

Active galactic nuclei (AGNs) are characterized by continuum emission across the electromagnetic spectrum and strong emission lines in the UV, optical, and near-IR. AGNs with both broad and narrow emission lines are classified as "Type 1" and those with only narrow emission lines as "Type 2," and this difference is usually attributed to our viewing angle relative to an obscuring midplane. AGNs are also known to be variable sources, showing variability in the continuum as well as in emission lines. The continuum variability is observed on both short and long timescales, but the variability amplitude is usually small, of the order of ∼20% over a year for low-luminosity AGNs (i.e., with Seyfert-like luminosity of about 1041–1044 erg s−1). The emission-line variations usually track the continuum variations, but see Goad et al. (2016).

A special class of AGN that has recently been gaining recognition has been termed "changing-look quasars."13 Their broad emission lines are observed to appear or disappear together with large changes in continuum luminosity. The changes in the emission-line spectra are such that they change type, e.g., from "Type 1" to "Type 2," and this is clearly an intrinsic change, implying that "type" is not always associated with viewing angle. This phenomenon has become of recent interest due to a combination of several serendipitous discoveries of changing-look AGNs (e.g., Shappee et al. 2014; Denney et al. 2014; LaMassa et al. 2015) and the growing availability of large spectroscopic quasar databases with long time-baselines of multi-epoch photometry and spectroscopy. These rare occurrences have the potential to contribute to our limited understanding of the central engines of active galaxies. Systematic searches of databases such as the SDSS/BOSS/TDSS have begun to uncover more changing-look quasars at larger redshifts and have led to predictions of the possible occurrence rates for quasars to change type (MacLeod et al. 2016; Runnoe et al. 2016).

Changing-look quasars are of particular interest for their potential to shed light on the physical origin of quasar variability and the details of accretion, since they have undergone such an extreme apparent change. Additionally, what is the role these objects play in the accretion history of the black hole (BH)—are we observing the beginning/end of a quasar phase? What role do these objects play in the larger context of galaxy and BH coevolution? Recent investigations of changing-look quasars have tried to use the observations before and after the change and any additional data available on the few known sources to test whether or not the changes are due to changes in accretion rate or obscuration, and to determine whether the quasar is turning on for the first time, reviving after a period of quiescence, or something else. In all cases so far, obscuration has not been identified as the major cause of changing look. Instead, a significant increase (decrease) in the mass accretion rate is the most favored explanation for the appearance (disappearance) of the broad emission lines.

The nearby, low-luminosity AGN Markarian 590 (Mrk 590 hereafter) remains an interesting case of a changing-look quasar in the local universe. Denney et al. (2014; D14 hereafter) have described the initial set of observations that catalog this object's observed history and "change," where it went from being a strong broad-line emitter three decades ago, allowing a reverberation mapping-based direct BH mass measurement to be made ((4.75 ± 0.74) × 107 M, Peterson et al. 2004), to being broad-line weak—the optical broad lines all but disappeared—in the past decade. Mrk 590 also shows the presence of ultrafast outflows in the X-ray band (Gupta et al. 2015). We have obtained additional UV and X-ray spectra since those presented by D14 that we describe in Sections 2 and 3. Implications of our results are discussed in Section 4 and we show that the "soft X-ray excess," observed in a large fraction of AGNs, is not due to the reprocessing of hard X-ray continuum, but is instead a result of thermal Comptonization of UV photons. We also argue that the changing-look phenomenon is a normal event in the duty cycle of "normal" quasars.

2. Chandra Spectroscopy

Mrk 590 was observed by Chandra  for 68.87 ks on 2014 November 15 (ObsID 16109). This is longer than the 30 ks initial Director's Discretionary Time (DDT) observation presented by D14, but the observing parameters were otherwise kept the same. We observed with ACIS-S/HETG; the grating was inserted in front of the detector in order to mitigate pile-up, if any. The data were reduced with CIAO following the normal procedure and analyzed using XSPEC. A Galactic column density of NH = 2.7 × 1020 cm−2 was included in all fits and was held fixed. All the quoted errors are for 90% confidence, unless noted otherwise. Figure 1 shows the new, 2014 spectrum, together with the 2013 spectrum presented by D14, and the previous "high state" spectrum observed in 2004. The source continues to be in a low state in 2014, so the quality of the high-resolution grating spectrum is too poor to perform useful science. Hence we present the zeroth-order spectrum here, as was done by D14.

Figure 1.

Figure 1. Chandra spectra of Mrk 590. The high-state spectrum of 2004 and the low-state spectra of 2013 and 2014 are shown. The spectra are fitted with simple power-law continuum models modified by Galactic absorption; these are shown with colored lines to guide the eye. The purpose of this figure is to show that the source continues to remain in the low state in 2014.

Standard image High-resolution image

A simple absorbed power-law model did not provide a good fit to the 2014 Chandra  spectrum. The residuals to this fit are shown in Figure 2, and clearly show a soft excess, so we added a blackbody component to the model. Clear residuals were also observed around 6 keV, so following Longinotti et al. (2007) we added two Gaussian emission lines at rest-frame 6.4 keV and 6.7 keV to represent neutral and ionized Fe Kα emission lines, respectively. The resulting fit was good (${\chi }_{\nu }^{2}=1.0$ for 104 degrees of freedom); the parameters of the fit are presented in Table 1 and the spectra are shown in Figures 3 and 4. Longinotti et al. (2007) required a third emission line at 7 keV; our data are not good enough in this energy range to require this line. They also detected a line in the soft X-ray band, at 19 Å (0.65 keV), but we see no evidence of that line.

Figure 2.

Figure 2. Residuals to the 2014 spectrum fitted with an absorbed power-law model. The soft excess is clearly visible below 1 keV. Note also the residuals between 6 and 7 keV, at the locations of Fe lines.

Standard image High-resolution image
Figure 3.

Figure 3. The best-fit spectrum and the residuals. The spectrum is fitted with an absorbed power law, a blackbody component to parameterize the soft excess, and two emission lines at 6.4 and 6.7 keV (rest frame).

Standard image High-resolution image
Figure 4.

Figure 4. The "unfolded" spectrum, i.e., the spectrum without the instrumental response folded in. The model components (power-law continuum, soft excess, and emission lines) are shown as dotted lines.

Standard image High-resolution image

Table 1.  Chandra  Best-fit Spectral Modela

Model Parameter 1 Normalizationb Parameter 3
Power law Γ = 1.6 ± 0.1 (2.4 ± 0.3) × 10−4  
Blackbody ${kT}={0.11}_{-0.02}^{+0.05}$ keV ${7}_{-4}^{+7}\times {10}^{-6}$  
Line 1 E = 6.4 keVc (5 ± 2) × 10−6 EW = 268 ± 100 eV
Line 2 E = 6.7 keVc (4 ± 2) × 10−6 EW = 176 ± 100 eV

Notes.

a ${\chi }_{\nu }^{2}$ for the joint fit is 1.0 for ν = 104 degrees of freedom. bIn units of photons s−1 cm−2 keV−1. cRest frame; parameter frozen.

Download table as:  ASCIITypeset image

The flux of the source is 2.2 × 10−12 erg s−1 cm−2 at 0.3–10 keV (1.9 × 10−12 erg s−1 cm−2 at 0.5–10 keV and 2.0 × 10−12 erg s−1 cm−2 at 0.3–8 keV). This is significantly smaller than the 2002 flux (8.4 × 10−12 erg s−1 cm−2) reported by Gallo et al. (2006), the 2004 flux (11 × 10−12 erg s−1 cm−2) reported by Longinotti et al. (2007), and the Suzaku  flux of 6.8 × 10−12 erg s−1 cm−2 at 2–10 keV reported by Rivers et al. (2012). D14 reported the low state of Mrk 590 with 0.5–10 keV flux of 1.1 × 10−12 erg s−1 cm−2, so compared to our DDT Chandra  observation in 2013, the source flux appears to have increased slightly, but as shown in Figure 1, the source continues to be in a low state. The flux in the soft excess component is 3.8 × 10−13 erg s−1 cm−2 in the 0.3–8 keV band (though it is concentrated below 2 keV), which is 19% of the total. If we restrict to the 0.3–2 keV soft band, then the total flux is 9.7 × 10−13 erg s−1 cm−2 while the soft excess flux remains the same as above, making it 39% of the total soft-band flux.

3. HST/COS NUV Spectroscopy

We obtained observations with the Cosmic Origins Spectrograph on board the Hubble Space Telescope (HST) that were coordinated with our approved Chandra  program. The goal of the COS observations was to observe the portion of the UV spectrum not covered by the initial COS G140L DDT observations presented by D14. The new observations were obtained on 2014 November 15 with the G230L grating and covered the observed-frame wavelength ranges ∼1670−2175 Å and ∼2765−3215 Å. Importantly, this coverage was selected to allow a continuum overlap region with the 2013 G140L spectrum and also to cover the Mg ii λ2800 doublet in this object of redshift z = 0.026385. The G230L spectrum was processed with the standard COS pipelines. Figure 5 shows the G230L COS spectrum, which was modeled with a power-law continuum, a Balmer continuum, and an Fe ii template (Vestergaard & Wilkes 2001) as a means to isolate the Mg ii emission, with both a narrow-line component and a broad-line component. The width and velocity shift of the narrow-line components of Mg ii and N iii] 1750 Å were tied, allowing only the strength to differ. We also corrected the spectrum for foreground reddening, assuming E(BV) = 0.0367 mag determined from the dust maps of Schlegel et al. (1998) and the extinction curve of O'Donnell (1994). As can be seen in the figure, the Mg ii λ2800 emission line clearly shows a broad component. We measured the FWHM = 10028 ± 627 km s−1 and the line dispersion σ = 4000 ± 790 km s−1 for the broad component from the emission-line models.

Figure 5.

Figure 5. HST/COS spectrum of Mrk 590 corrected for Milky Way extinction with continuum and line models. Observed flux densities are shown as a function of the rest-frame wavelengths. Top panel: the observed spectrum with the power-law continuum model (purple curve), the Balmer continuum model (blue curve), the iron emission model (red curve), and the sum of the three models (green curve) shown superimposed. Middle panel: residuals obtained by subtracting the continuum, iron, and Balmer emission models shown along with the flux errors (gray dotted curve). Indicated in black are the wavelength regions used to fit the power-law continuum. The Mg ii emission line is clearly broad. Bottom panel: emission-line spectrum (black) obtained by subtracting the continuum, iron, and Balmer emission models (shown in the top panel) from the observed spectrum with the line models superimposed. Model fits were made to N iii] 1750 Å (red curve) and the broad (magenta curve) and narrow (red curves) components of Mg ii. Red curves show the narrow-line components while brown, blue, and yellow curves represent the three Gaussian functions used to model the broad Mg ii component (magenta curve). The gray curve (around zero flux level) shows the residuals after subtracting the emission-line model. Yellow stars indicate pixels deviating more than 3σ from the rms in the spectrum, which were filtered out prior to the line modeling. The flux errors are shown by the green dotted curve.

Standard image High-resolution image

4. Discussion

4.1. Soft Excess and the Fe Lines

Over the years, X-ray observations of Mrk 590 have shown that its spectrum is well fit by a power-law continuum, a soft excess, a hard reflection component, and Fe emission lines (Section 1). Rivers et al. (2012), using Suzaku  data, first noticed that while the hard X-ray continuum of the source remained practically the same as in the historic data, its soft excess had disappeared in 2011. In the 2013 Chandra  spectrum of D14, there was no evidence of a soft excess either, and the source was in a historically low state. Our new Chandra  observations show that Mrk 590 continues to be in a low state, but a soft excess has re-emerged. The soft-excess temperature (${kT}={0.11}_{-0.02}^{+0.05}$ keV) is lower than kT = 0.18 ± 0.01 keV seen in the 2007 XMM-Newton  spectrum reported by Longinotti et al. (2007) and Rivers et al. (2012). The soft-excess luminosity of 5 × 1040 erg s−1 cm−2 is also an order of magnitude lower than the 2007 value of (8 ± 1) × 1041 erg s−1 cm−2, but higher than the upper limit derived from the 2011 Suzaku  spectra of (3 ± 1) × 1040 erg s−1 cm−2 (Rivers et al. 2012). Thus the soft excess seems to have returned, but not to the "normal" level seen in the historical high state.

In Mathur et al. (2017 and references therein) we discussed various models of the soft excess, which we briefly discuss here. The main contenders are (1) reflection of the hard X-ray source by the accretion disk (e.g., Crummy et al. 2006); (2) a warm Comptonizing medium around the accretion disk (e.g., Ross et al. 1992; Titarchuk 1994); (3) thermal emission from an accretion disk. Thermal emission from an accretion disk may appear in the soft X-ray band only for relatively low-mass BHs. The re-emergence of the soft excess in the low state, when the hard X-ray emission was weak, suggests that the soft excess is not due to reprocessed hard X-ray emission. An alternative model, in which the soft excess is generated by thermal Comptonization in the optically thick accretion disk corona, is a more likely explanation.

Using a series of observations of Mrk 509, Mehdipour et al. (2011) found that the soft excess is correlated with the UV flux as Lsoft excess = 1.261(FUVW2)2.869, where Lsoft excess is the 0.3–2 keV soft-excess luminosity in units of 1042 erg s−1 and FUVW2 is the flux in the UVW2 filter in UVOT on Swift  and OM on XMM-Newton, centered at 2030 Å, in units of 10−14 erg s−1 cm−2 Å−1. In our new observations, the soft-excess luminosity of Mrk 590 in the 0.3–2 keV band is 3.5 × 1040 erg s−1 and the continuum flux at 2030 Å is 7.7 × 10−16 erg s−1 cm−2 Å−1. The observed soft-excess luminosity is two orders of magnitude larger that that predicted by the Mehdipour relation (correcting the UV flux for the distance of NGC 5548) and shows that the Mehdipour relation between the soft excess and the UV flux is not universal.

The Fe Kα emission lines have smaller equivalent widths and smaller fluxes than the values in 2004 reported by Longinotti et al. (2007). This shows that both the neutral and the ionized lines have responded to the continuum changes between 2004 and 2014. This places the reprocessing material within 10 lt-yr or 3 pc of the continuum source, and not from the kiloparsec-scale extended emission region as suggested by Longinotti et al. (2007) for the ionized line.

4.2. The ${Mg}\,{II}$ Emission Line

Our HST  spectrum clearly shows a broad Mg ii emission line. Since we do not have a UV spectrum covering the Mg ii line in the low state when the Hβ line had vanished, the presence of the Mg ii line in the current spectrum can be interpreted in a couple of different ways.

Photoionization models and evidence from reverberation mapping of the ionization stratification of the broad-line region (BLR) support the co-spatial existence of Hβ and Mg ii emitting gas. Thus it is possible that when the broad Hβ emission line disappeared, the broad Mg ii line had also disappeared, and now that a broad Mg ii line has emerged, possibly a broad Hβ has also emerged. We should note, however, that the response of a line to continuum variations depends not only on the geometry and location of line-emitting region, but also on line responsivity. Cackett et al. (2015) have shown that the Mg ii emission line has very low responsivity compared to other high-ionization emission lines (see also Korista & Goad 2000). It is thus possible that the broad Mg ii line had never disappeared. If this is the case, then we cannot comment on the presence of a broad Hβ line in 2014. Roig et al. (2014) found a class of AGNs in the Baryon Oscillation Spectroscopic Survey showing a broad Mg ii without a broad Hβ. Wavelength-dependent scattering of the BLR is offered as an explanation by these authors as to why the broad Balmer lines could be missing but the broad Mg ii lines are strong. Our results on Mrk 590 suggest that such AGNs may have been in a transition state when Hβ disappeared, but Mg ii did not, owing to its low responsivity.

The connection between the soft excess and the emission lines is of interest. In the warm Comptonization model, the seed UV/EUV photons are Compton up-scattered into soft X-rays in the optically thick accretion disk corona. Observationally, the soft excess has been found to correlate with UV emission (e.g., Atlee & Mathur 2009). Thus the emergence of soft excess in Mrk 590 implies a corresponding emergence of the EUV ionizing continuum. This in turn may be responsible for the reappearance of the BLR, which we see with Mg ii. This is our preferred model. D14 found the continuum flux at 1450 Å to be 3.7 × 10−16 erg s−1 cm−2 Å−1, which was two orders of magnitude lower than the IUE value (D14). The 1634 Å flux in our new HST  data is 9.2 × 10−16 erg s−1 cm−2 Å−1, which is 15.1 × 10−16 erg s−1 cm−2 Å−1 when extrapolated to 1460 Å. Thus, while the source continues to be in a low state, the UV continuum is re-emerging along with the soft excess. While the presence of the broad Mg ii line can be interpreted in different ways, as noted above, it is consistent with the increasing luminosity of the ionizing continuum.

5. Conclusion

D14 presented optical, UV, and X-ray observations of the classical Seyfert 1 Mrk 590 that span the past 40+ years. This interesting object brightened by a factor of a few tens between the 1970s and 1990s and then faded by a factor of 100 or more at all continuum wavelengths between the mid-1990s and 2013. Notably, there is no evidence in the current data set that this recent, significant decline in flux is due to obscuration; in particular, the most recent X-ray observations are consistent with zero intrinsic absorption. There were similarly dramatic changes in the emission-line fluxes, the most striking of which is the complete disappearance of the broad component of the Hβ emission line, which had previously been strong (equivalent widths ∼20−60 Å).

In this paper we have presented new Chandra  and HST  observations from 2014 and we find that Mrk 590 is awakening. While the source continues to be in a low state, its soft excess has re-emerged, though not to the historical level. The broad Mg ii emission line is also present. These observations suggest that the soft excess is connected to the UV continuum through warm Comptonization and that the ionizing continuum is also on the rise, but it does not follow the relation between soft excess and UV continuum reported in Mehdipour et al. (2011).

The implications from this long time series of Mrk 590 observations are that (1) Mrk 590 is a direct challenge to the historical paradigm that AGN type is exclusively a geometrical effect, and (2) there may not be a strict, one-way evolution from Type 1 to Type 1.5−1.9 to Type 2 as recently suggested by Elitzur et al. (2014). Instead, for at least some objects, the presence of BLR emission may coincide with episodic accretion events throughout a single active phase of an AGN. If true, such behavior may be more prominent in Seyfert galaxies, where accretion is likely to be a consequence of secular processes (e.g., Mathur et al. 2012; Martin et al. 2018) and therefore likely more episodic than quasar activity, which may be governed more predominantly by major mergers.

Support for this work was provided by the National Aeronautics and Space Administration through Chandra Award Number G04-15114X to S.M. issued by the Chandra X-ray Observatory Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space Administration under contract NAS8-03060. Support for the HST  program number GO-13185 was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. K.D.D. was supported by an NSF AAPF fellowship awarded under NSF grant AST-1302093 while she worked on this project. M.V. gratefully acknowledges support from the Danish Council for Independent Research via grant number DFF 4002-00275.

Footnotes

  • 13 

    We make no distinction among the terms quasar, AGN, and Seyfert galaxy—all are objects that contain a central black hole that is actively accreting sufficient matter to shine as a point source and that, to first order, produce the same emission characteristics.

Please wait… references are loading.
10.3847/1538-4357/aadd91