After the AllSky Automated Survey for SuperNovae discovered a significant brightening of the inner region of NGC 2617 we began a 70 day photometric and spectroscopic monitoring campaign from the Xray through nearinfrared NIR wavelengths. We report that NGC 2617 went through a dramatic outburst during which its Xray flux increased by over an order of magnitude followed by an increase of its opticalultraviolet UV continuum flux by almost an order of magnitude. NGC 2617 classified as a Seyfert 1.8 galaxy in 2003 is now a Seyfert 1 due to the appearance of broad optical emission lines and a continuum blue bump. Such changing look active galactic nuclei AGNs are rare and provide us with important insights about AGN physics. Based on the H line width and the radiusluminosity relation we estimate the mass of central black hole BH to be 4  1  10SUP7SUP M SUBSUB. When we crosscorrelate the light curves we find that the disk emission lags the Xrays with the lag becoming longer as we move from the UV 23 days to the NIR 69 days. Also the NIR is more heavily temporally smoothed than the UV. This can largely be explained by a simple model of a thermally emitting thin disk around a BH of the estimated mass that is illuminated by the observed variable Xray fluxes.
['galaxies: active', 'galaxies: nuclei', 'galaxies: Seyfert', 'line: formation', 'line: profiles', 'Astrophysics - High Energy Astrophysical Phenomena', 'Astrophysics - Cosmology and Nongalactic Astrophysics']
The Man behind the Curtain XRays Drive the UV through NIR Variability in the 2013 Active Galactic Nucleus Outburst in NGC 2617
{'The Man Behind the Curtain: X-rays Drive the UV through NIR Variability in the 2013 AGN Outburst in NGC 2617': 'B. J. Shappee 1 , 2 , J. L. Prieto 3 , 4 , D. Grupe 5 , C. S. Kochanek 1 , 6 , K. Z. Stanek 1 , 6 , G. De Rosa 1 , S. Mathur 1 , 6 , Y. Zu 1 , B. M. Peterson 1 , 6 , R. W. Pogge 1 , 6 , S. Komossa 7 , M. Im 8 , J. Jencson 1 , T.W-S. Holoien 1 , U. Basu 1 , 9 , J. F. Beacom 1 , 6 , 10 , D. M. Szczygieł 11 , J. Brimacombe 12 , S. Adams 1 , A. Campillay 13 , C. Choi 8 , C. Contreras 13 , M. Dietrich 14 , M. Dubberley 15 , M. Elphick 15 , S. Foale 15 , M. Giustini 16 , C. Gonzalez 13 , E. Hawkins 15 , D. A. Howell 15 , 17 , E. Y. Hsiao 13 , M. Koss 18 , K. M. Leighly 19 , N. Morrell 13 , D. Mudd 1 , D. Mullins 15 , J. M. Nugent 19 , J. Parrent 15 , M. M. Phillips 13 , G. Pojmanski 11 , W. Rosing 15 , R. Ross 15 , D. Sand 20 , D. M. Terndrup 1 , S. Valenti 15 , 17 , Z. Walker 15 , and Y. Yoon 8  \nshappee@astronomy.ohio-state.edu  \nAugust 20, 2018', 'ABSTRACT': "After the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered a significant brightening of the inner region of NGC 2617, we began a ∼ 70 day photometric and spectroscopic monitoring campaign from the X-ray through near-infrared (NIR) wavelengths. We report that NGC 2617 went through a dramatic outburst, during which its X-ray flux increased by over an order of magnitude followed by an increase of its optical/ultraviolet (UV) continuum flux by almost an order of magnitude. NGC 2617, classified as a Seyfert 1.8 galaxy in 2003, is now a Seyfert 1 due to the appearance of broad optical emission lines and a continuum blue bump. Such 'changing look Active Galactic Nuclei (AGN)' are rare and provide us with important insights about AGN physics. Based on the H β line width and the radius-luminosity relation, we estimate the mass of central black hole to be (4 ± 1) × 10 7 M /circledot . When we crosscorrelate the light curves, we find that the disk emission lags the X-rays, with the lag becoming longer as we move from the UV (2 -3 days) to the NIR (6 -9 days). Also, the NIR is more heavily temporally smoothed than the UV. This can largely be explained by a simple model of a thermally emitting thin disk around a black hole of the estimated mass that is illuminated by the observed, variable X-ray fluxes.  \nSubject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert - line: formation line: profiles", '1. Introduction': "NGC 2617 is a face-on Sc galaxy at a redshift of z = 0 . 0142 (Paturel et al. 2003), corresponding to a Hubble flow corrected distance of 64 Mpc. Moran et al. (1996) classified it as a Seyfert 1.8 galaxy. The galaxy displayed only narrow emission lines except for a weak, broad base to the H α line. It lacked both the broad emission lines that are characteristic of Seyfert 1 galaxies, and a bright power-law continuum. Here we report on the dramatic transformation of NGC 2617 into a Seyfert 1 and the monitoring of the (presumably) associated X-ray-NIR outburst.  \nIt is well known that both the broad line and continuum components of active galactic nuclei (AGN) are time variable at all wavelengths. This enables, for example, the reverberation mapping (RM) technique for studying black hole (BH) masses (e.g., Peterson et al. 2004; Grier et al. 2012b,a, 2013b) or measuring time delays in gravitational lenses (e.g. Tewes et al. 2013). However, there are only a few cases where the changes are so large that the broad-line components of H α and H β appear in a Seyfert 2 or vanish from a Seyfert 1. Vanishing broad-line components have been reported for NGC 7603 (Tohline & Osterbrock 1976), NGC 4151 (Penston & Perez 1984), Mrk 372 (Gregory et al. 1991), and 3C390.3 (Penston & Perez 1984; Veilleux & Zheng 1991), and broad-line components have appeared in Mrk 6 (Khachikian & Weedman 1971), Mrk 1018 (Cohen et al. 1986), Mrk 993 (Tran et al. 1992), NGC 1097 (Storchi-Bergmann et al. 1993), and NGC 7582 (Aretxaga et al. 1999). 'Changing-look' AGN have been observed in X-rays as well  \n(e.g., Matt et al. 2003). The X-ray absorption in IRAS 09104 + 4109, known as a Seyfert 2 AGN, changed from being Compton-thick to Compton-thin, making it essentially a Seyfert 1 AGN (e.g., Lusso et al. 2013). Similar variability was subsequently observed in a several other AGN (Marchese et al. 2012). These strong variations must be telling us something fundamental about the structure of the broad-line region or the obscuring material.  \nAdditionally, some AGN show large changes in their (X-ray) continuum emission. An example is IC 3599, which flared by approximately a factor of a hundred in soft X-rays, accompanied by the appearance and then disappearance of broad and narrow optical emission lines (Brandt et al. 1995; Grupe et al. 1995b; Komossa & Bade 1999). Other AGN show high-amplitude X-ray variability without significant emission line variability. Among the most extreme cases are WPVS007 (e.g., Grupe et al. 1995a, 2013), E 1615+061 (Piro et al. 1997), Mkn 335 (Grupe et al. 2007, 2008, 2012), and PHL 1092 (Miniutti et al. 2009). There are examples of still more extreme changes in X-ray emission, up to factors of several thousand, which have been detected in optically quiescent galaxies. These events have been interpreted as stellar tidal disruption flares, which are thought to occur when stars are disrupted and subsequently accreted by SMBHs in galaxy cores (see Komossa 2012 for a recent review).  \nBroad emission lines are one of the defining characteristics of Seyfert 1 AGN. In the unification model of AGN, it is believed that Seyfert 2 AGN are Seyfert 1 AGN where the broad emission line region (BLR) is being blocked from view by an obscuring torus of gas and dust. However, the intermediate Seyfert types (1.8 and 1.9) do not easily fit into this picture. Goodrich (1989) suggested that the difference between Seyfert 1 AGN and these intermediate Seyfert types is simply that the intermediate Seyferts are viewed through a screen of dust and that the variability of these lines is due to the changing optical depth of this material. Beyond viewing angles and obscuration, the fundamental parameters determining the presence of a BLR are still uncertain, but are likely related to the BH mass and accretion rate. For example, Chakravorty et al. (2014) have proposed that the presence of a BLR depends on the BH mass, with objects below a critical mass showing no BLR. Nicastro (2000) offers a different model, arguing that the accretion rate relative to Eddington is the fundamental parameter. In this model, the inner radius of the BLR is the transition radius at which the accretion disk goes from being radiation-to-gas pressure dominated and below a critical mass accretion rate no BLR should form. Czerny & Hryniewicz (2011) offer a third model. They note that the accretion disk temperature at the radius of the BLR is about 1000 K in all AGN. This is close to the dust sublimation temperature, which may allow the formation of a radiation-driven wind which then forms the BLR.  \nOn the other hand, the origins of the optical/UV AGN continuum are known to be thermal emission from an accretion disk around a supermassive BH (Shields 1978). The exact structure of the disk is uncertain, and disk size measurements using gravitational microlensing generally find sizes larger than expected from simple thin disk theory (e.g., Morgan et al. 2010).  \nThe nature of the X-ray continuum and its relation to the UV continuum, however, is still a matter of debate. The X-rays are generally believed to arise from Compton up-scattering of disk photons in a hot corona (e.g., Reynolds & Nowak 2003). Alternatively, optical/UV continuum may be produced from down-  \nscattering of X-ray photons (e.g., Cameron et al. 2012 and references therein). Theoretical models have postulated a broad range of physical sizes for this hot corona. However, recent X-ray gravitational microlensing results consistently find sizes comparable to the inner edge of the accretion disk (e.g., Morgan et al. 2008; Mosquera et al. 2013), thus ruling out models with very extended coronae.  \nThe causal ordering of X-ray and disk variability is, however, unclear, since attempts to measure time lags between these emission components have generally yielded only null results or tentative detections (e.g. Breedt et al. 2009; Cameron et al. 2012). These lags, however, could be a powerful tool to distinguish between the two physical processes most likely responsible for strong X-ray-NIR variability observed in NGC 2617 which are 1) locally-generated (e.g. viscous) perturbations propagating inward through the accretion disk or 2) an increase in the X-ray flux coming from the smaller, central hot corona which heats the inside of the disk first and moves outward. Thus, changes in accretion rate might be expected to produce 'outside-in' variations, moving from red to blue to X-rays, while changes in coronal emission might be expected to produce 'inside-out' variations with the increased X-ray flux irradiating the disk and driving increased blue and then red emission.  \nDelays between wavelengths might also be measurable for the UV/optical emission. The best case to date is that of NGC 7469 (Wanders et al. 1997; Collier et al. 1998; Peterson et al. 1998; Kriss et al. 2000), where longer wavelengths increasingly lag behind those in the 1315 Å continuum. The lag detections are significant only at about the 2-sigma level, however. Surprisingly, the X-ray variations (Nandra et al. 1998) did not correlate particularly well with those in the UV/optical, and in retrospect this is probably attributable to variable absorption of the X-ray spectrum (e.g., Risaliti 2011). Lags between shorter- and longer-wavelength continuum variations have been detected in several sources (e.g., Sergeev et al. 2005; Cackett et al. 2007) at marginal significance, and always in the sense that the short-wavelength variations lead those at longer wavelengths.  \nThe dramatic state change and outburst we observe in NGC 2617 provides a laboratory for exploring the relations between these AGN components. In §2 we describe our NIR -X-ray observations during the outburst of NGC 2617. In §3 we search for variations in the broad H β and He II λ 4686 lines, estimate the central BH mass, and cross-correlate the X-ray to UV-NIR light curves to measure the time delays between the variability in different bands. We present a simple physical model that reproduces the observed UV-IR variability in §4 and summarize our results in §5. Throughout the paper all dates are given in UT.", '2. Observations': "Wehave been conducting the All-Sky Automated Survey for SuperNovae (ASAS-SN 21 or 'Assassin'), which scans the extragalactic sky visible from Hawaii roughly once every five nights in the V -band. On 2013 Apr. 10.27, a transient source alert was triggered by the brightening of the central region of NGC 2617. Follow-up observations confirmed that the nucleus was the source of the outburst and a follow-up spec-  \ntrum showed that the source now has the strong broad emission lines and continuum shape characteristic of a Seyfert 1. Given the strong outburst and dramatic change in the spectral properties, we requested and obtained Swift (Gehrels et al. 2004) target of opportunity (ToO) observations to monitor NGC 2617 in the optical, UV, and X-rays, which found that NGC 2617 continued to brighten. During this outburst we were able to obtain an extensive time series of optical spectroscopy and NIR/optical/UV/X-ray photometry that are presented in this paper. The INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL; Winkler et al. 2003) also detected the outburst in hard X-rays (17 -60 keV), as has the European VLBI Network real-time Interferometer (e-VLBI) and the Very Long Baseline Array (VLBA) in the radio at 1.6 GHz and 5 and 1.7 GHz, respectively. Finally, XMM-Newton observations were obtained for NGC 2617 during the outburst, and these observations will be discussed in a companion paper (Giustini et al. in prep.). In the following subsections we describe the observations of NGC 2617 in detail.", '2.1. ASAS-SN Discovery of NGC 2617 Outburst': "ASAS-SNis a long-term project to monitor the whole extragalactic sky to find nearby supernovae (SNe) and other transient sources. We began running our real-time search for variable sources in April 2013 with our first unit, Brutus. Brutus presently consists of two telescopes on a common mount which is hosted by Las Cumbres Observatory Global Telescope Network (LCOGT; Brown et al. 2013) in the Faulkes Telescope North enclosure on Mount Haleakala, Hawaii. Each telescope consists of a 14-cm aperture Nikon telephoto lens and an FLI ProLine CCD camera with a Fairchild Imaging 2k × 2k thinned CCD, giving a 4 . 47 × 4 . 47 degree field-of-view and a 7 '' . 8 pixel scale. The cameras are mounted to have moderately overlapping fieldsof-view. On a typical clear night, Brutus can survey more than 5000 deg 2 , and with the effects of weather we can survey the extragalactic sky visible from Hawaii roughly once every 5 days. The data are reduced in real-time, and we can search for transient candidates about an hour after the data are taken using an automated difference imaging pipeline. We are now meeting, and frequently exceeding, our current depth goal of V ≈ 16 mag, corresponding to the apparent brightness at maximum light of core-collapse SNe within 30 Mpc and SNe Ia out to 100 Mpc.  \nOur transient source detection pipeline was triggered by an image taken at 2013 Apr. 10.27, detecting a ∼ 10% flux increase from the inner region of NGC 2617 (Figure 1), equivalent to detecting a new point source source of V ∼ 16 . 7 mag superimposed on the image of the galaxy (Shappee et al. 2013b). It seemed likely to be AGN variability, but the AGN and the host are not separately resolved given ASAS-SN's 16 '' resolution. We obtained a follow-up image on Apr. 24.14 with the Ohio State Multi-Object Spectrograph (OSMOS; Martini et al. 2011) at the 2.4-m telescope at MDM Observatory, which showed that the central region of NGC 2617 continued to brighten and was a factor of ∼ 3 . 3 brighter than the archival Sloan Digital Sky Survey (SDSS; York et al. 2000) Data Release 9 (Ahn et al. 2012) g -band image, as shown in Figure 1. This increase in flux was within 0 '' . 03 of the center of the galaxy, corresponding to a projected physical distance of less than ∼ 9 pc at NGC 2617.  \nWe then requested director's discretionary (DD) time on the Apache Point Observatory (APO) 3.5-m  \nFig. 1.- Comparison between our ASAS-SN discovery image (top-left), the archival SDSS image (topright), our follow-up MDM image (bottom-left), and a difference image showing the brightening of the central region between the SDSS and MDM images (bottom-right). All images are on the same angular scale. This increase in flux is associated with the center of the galaxy to within 0 '' . 03 or ∼ 9 pc at NGC 2617.  \n<!-- image -->  \ntelescope and obtained an optical spectrum (range 3500 -9600 Å) using the Dual Imaging Spectrograph (DIS) on Apr. 25.1. Figure 2 shows a comparison between this spectrum, a spectrum from the 6dF Galaxy Survey (6dFGS) (Jones et al. 2004, 2009) taken on 2003 Apr. 25, and the spectrum from Moran et al. (1996) taken in 1994. This figure illustrates the dramatic spectral changes of NGC 2617. Along with the increase in the continuum flux, NGC 2617 changed spectral type from a Seyfert 1.8 to a Seyfert 1, with the appearance of strong broad-line emission in H α and H β . The spectrum also shows the appearance of a strong 'blue bump' shortward of 4000 Å. It is possible that t he spectral changes observed in NGC 2617 are not correlated with the large AGN outburst discovered by ASAS-SN. However, given the rarity of both 'changing-look' AGN and powerful X-ray flares, this seems a remote possibility even if the available data cannot prove an association.", '2.2. Swift ToO Optical, Ultraviolet, and X-ray Observations': "The first 2 ks Swift X-ray Telescope (XRT; Burrows et al. 2005) observation obtained on 2013 Apr. 30.93 detected NGC 2617 with a flux of 3 . 1 × 10 -11 ergs cm -2 s -1 (0.3-10 keV; Mathur et al. 2013). This X-ray flux is significantly higher than the fluxes measured for NGC 2617 by both the XMM-Newton slew survey in 2007 (Saxton et al. 2008) of 9 . 0 × 10 -12 ergs cm -2 s -1 (0.212 keV) and the Röntgen Satellite (ROSAT; Truemper 1982) of (6 . 9 ± 1 . 0) × 10 -12 ergs cm -2 s -1 (0.1-2.4 keV; Boller et al. 1992). Based on this initial result we requested and were granted further Swift ToO observations. By 2013 May 8.04, continued monitoring with Swift XRT showed an increase of another factor of four in the X-ray flux, so the AGN was now an order of magnitude brighter than during the XMM-Newton Slew Survey observation in 2007. Parallel optical and UV observations obtained with Swift 's UltraViolet/Optical Telescope (UVOT; Roming et al. 2005) also showed an increase in flux, but on a longer timescale (Shappee et al. 2013a). In total, we obtained Swift XRT and UVOT observations spanning ∼ 50 days with an almost daily cadence from 2013 May 01 through June 20 when it became unobservable due to Swift 's Sun-constraint.", '2.2.1. Swift UVOT Observations': "The Swift UVOT observations of NGC 2617 were obtained in 6 filters (Poole et al. 2008): V (5469 Å), B (4392 Å), U (3465 Å), UVW 1 (2600 Å), UVM 2 (2246 Å), and UVW 2 (1928 Å). We used the UVOT software task uvotsource to extract the source and background counts from a 5 '' . 0 radius region and a 5 '' . 0 to 20 '' . 0 annulus, respectively. The UVOT count rates were converted into magnitudes and fluxes based on the most recent UVOT calibration (Poole et al. 2008; Breeveld et al. 2010). The UVOT UV/optical photometry is reported in Table 1 and shown in Figure 3.  \nFig. 2.- Comparison between the APO spectrum from 2013 April 25 (top, red), the 6dFGS spectrum from 2003 Dec. 30 (middle, black), and the Moran et al. (1996) spectrum from 1994 (bottom, blue). Note that the absolute flux calibration is arbitrary. NGC 2617 underwent a dramatic spectral evolution changing from a Seyfert 1.8 to a Seyfert 1 sometime between 2003 and 2013.  \n<!-- image -->  \nTable 1. Photometric Observations  \nNote. Swift , U , B , V , R , I , Z , Y , J , H , and K -band photometry are calibrated in the Vega magnitude system. u , g , r , and i -band photometry are calibrated in the AB magnitude system. X-Ray fluxes (0.3 -10 keV) are in units of 10 -11 ergs/cm 2 /s. Only the first observation in each band is shown here to demonstrate its form and content. Table to be published in its entirety in machine-readable form.  \nFig. 3.- NIR-UV photometric observations of NGC 2617 spanning ∼ 60 days from 2013 Apr. 24 through June 20. The filled circles are our new observations while the large stars to the left are archival magnitudes from SDSS and 2MASS. Error bars are shown, but they are usually smaller than the points.  \n<!-- image -->", '2.2.2. Swift XRT Observations': "The XRT was usually operating in Windowed Timing (WT) mode (Hill et al. 2004) but it was in photon counting (PC) mode for some observations at the beginning and middle of the monitoring campaign. The data were reduced by the task xrtpipeline version 0.12.6, which is included in the HEASOFT package 6.12. For the PC mode data the source counts were collected in a circle with a radius of 70 '' . 7. Background counts were estimated from a source-free region 212 '' in radius using the task xselect (version 2.4b). The WT source and background data were selected in a box with a length of 40 pixels. Auxiliary response files were created using the XRT task xrtmkarf . The spectra were rebinned with 20 counts per bin using the task grppha . We applied the response files swxpc0to12s6\\_20010101v013.rmf and swxwt0to2s6\\_20070901v012.rmf for the PC and WT data, respectively. The re-binned 0 . 3 -10 keV spectra were modeled with XSPEC v.12.7 using a single power law and Galactic absorption corresponding to a hydrogen column density of N H = 3.64 × 10 20 cm -2 (Kalberla et al. 2005). The Swift XRT X-ray fluxes and XRT X-ray photon spectral indices ( Γ ) are also reported in Table 1. Figure 4 shows both the X-ray flux and UV magnitudes on the same relative scale as well as the X-ray spectral index. The X-ray spectral index becomes softer as the X-ray flux increases, as has been seen in other AGN (e.g., Romano et al. 2002; Grupe et al. 2012).", '2.3. Additional Photometric Follow-up': "We obtained optical and NIR photometry from a number of ground-based observatories. In an attempt to homogenize the photometric data, we performed photometry using a 5 '' radius aperture on all our data to minimize the differences between the Swift UVOT observations, which have a large pixel scale, and our ground-based observations, which have a much smaller pixel scale. We measured the background within a annulus of 30 '' . 0 to 60 '' . 0.  \nThe Z , Y , J , H , and K -band observations were obtained with the 3.8-m United Kingdom Infrared Telescope (UKIRT). The data are calibrated using Two Micron All-Sky Survey (2MASS; Skrutskie et al. 2006) stars within 7 ' from NGC 2617 and then converted to the WFCAM system following Hodgkin et al. (2009).  \nB , V , R , and I -band observations were taken at the LCOGT 1-m facilities in Cerro Tololo, Chile and Sutherland, South Africa (Brown et al. 2013). All images were reduced following standard procedures and cosmic rays were removed using LACOSMIC (van Dokkum 2001). The photometric measurements include a small aperture-corrections ( /lessorsimilar 0 . 02 mag). The data were calibrated using SDSS stars within 1 ' . 5 of NGC 2617 and transformed onto the Johnson-Cousins magnitude system using procedures described by Lupton (2005). 22  \nSDSS g , r , and i -band observations were obtained with the Retractable Optical Camera (RETROCAM;  \nFig. 4.Lower Panel : The Swift XRT 0 . 3 -10 keV X-ray (black circles), Swift UVOT UV (magenta, violet, and purple circles), and INTEGRAL IBIS 17 -60 keV X-ray (apricot rectangle) observation(s) of NGC 2617. Upper Panel : The Swift XRT X-ray photon spectral index ( Γ ). The apricot rectangle represents the INTEGRAL IBIS hard X-ray flux measurement, where the width and height represent the span of the observation and the measurement uncertainty, respectively. The large star and square to the left are archival XMM-Newton measurement (0.2-12 keV) from 2007 and archival ROSAT measurement (0.1-2.4 keV) from 1990-1991, respectively. The apricot upper limit shown to the left is an archival flux limit from the 70-month Swift /BAT survey (14-195 keV; see §2.5). As NGC 2617 brightened the X-ray spectrum softened. Both the UV and X-ray light curves show a similar double humped feature during the outburst, but with a 2-3 day delay (see §3.2).  \n<!-- image -->  \nMorgan et al. 2005) and SDSS g and r -band observations were obtained with OSMOS both on the MDM 2.4-m telescope. Small aperture-corrections ( /lessorsimilar 0 . 02 mag) were applied to the aperture magnitudes, which were then calibrated using SDSS stars within 1 ' . 5 of NGC 2617.  \nB , V , u , g , r , and i -band observations were obtained with the Direct CCD Camera (CCD) at the Swope 1m telescope at Las Campanas Observatory. Small aperture-corrections ( /lessorsimilar 0 . 02 mag) were applied to the aperture magnitudes, which were then calibrated using SDSS stars within 1 ' . 5 of NGC 2617 and transformed onto the Johnson magnitude system as described by Lupton (2005).  \nFigures 3 and 4 show that the amplitude of the flare was largest in the X-rays, followed by the UV, and then the optical/NIR. We note, however, that there is significant host galaxy contamination at longer wavelengths which is discussed in §3.3. For comparison, we also show the archival SDSS and 2MASS magnitudes of NGC 2617 from 2006 Nov 18 and 1998 Nov 27, respectively. The 2MASS magnitudes were translated onto the WFCAM system following Hodgkin et al. (2009).", '2.4. Spectroscopic Follow-up': "We obtained several low-resolution optical spectra of NGC 2617 between 2013 April and June using different telescopes, instruments, and spectroscopic setups. Table 2 shows a summary of all the spectra, including date (UT+JD), telescope/instrument, wavelength range, spectral resolution, slit width, seeing estimate, slit position angle, airmass at the beginning of the observation, and exposure time.  \nThe single slit spectra from DIS at the APO 3.5-m, Boller & Chivens CCD Spectrograph (CCDS) at the MDM 2.4-m and 1.3-m, OSMOS at the MDM 2.4-m, and the Inamori-Magellan Areal Camera & Spectrograph (IMACS; Dressler et al. 2011) at the Magellan 6.5-m were all reduced with standard routines in the IRAF twodspec and onedspec packages. The reductions included bias subtraction, flat-fielding, 1D spectral extraction, wavelength calibration using an arc-lamp, and flux calibration using a spectroscopic standard usually taken the same night.  \nThe integral field unit spectrum obtained with the Supernova Integral Field Spectrograph (SNIFS; Lantz et al. 2004) at the University of Hawaii (UH) 2.2-m was reduced and extracted using the SNIFS Data Reduction Package (Aldering et al. 2006). A 3 '' . 2 aperture was used to extract the final 1D spectrum.  \nThe spectrum obtained with the FLOYDS Spectrograph (Sand et al., in preparation) at the robotic 2-m Faulkes Telescope South (FTS; Brown et al. 2013) was reduced by an automated reduction pipeline that performs order rectification, flat-fielding, wavelength and flux calibration, and object extraction. Direct checks of the wavelength solution are also performed using telluric absorption features and night sky lines.  \nThe spectra were taken under various conditions, many of which were non-photometric. Additionally, many of the spectra were taken at high airmass and/or during astronomical twilight because the object was close to setting for the season. These conditions caused the normal flux calibration of the spectra to be poor. We perform an additional step and calibrated each spectrum onto an absolute flux scale under the assumption that the flux of the [O III] λ 5007 narrow emission line is constant as discussed in §3.1. We chose  \nthe MDM spectrum taken on 2013 May 12.13 as a reference spectrum for the absolute flux calibration, since it was obtained under clear conditions. In this spectrum the line flux is F ([O III] λ 5007) = 5 . 5 × 10 -14 erg cm -2 s -1 . These calibrated spectra are shown in Figure 5 and our search for temporal variations in the H β and He II λ 4686 lines is described in §3.1. Note, the spectra have not been corrected for telluric absorption, but the spectral region near H β is not significantly affected.", '2.5. Hard X-Rays': 'Between 2013 Apr. 22 and 2013 Jun. 2, INTEGRAL discovered a new source, IGR J08357 -0404, coincident with NGC 2617 during their Galactic Latitude Scan observations of the Puppis region (Tsygankov et al. 2013). They detected this source at a 9-sigma significance with a flux of 3 . 3 ± 0 . 4 mCrab (4 . 7 × 10 -11 erg cm -2 s -1 ) in the 17 -60 keV energy range. Correcting for the distance to NGC 2617, this flux corresponds to a luminosity of L = 1 . 86 × 10 43 erg s -1 . This observation is shown in Figure 4 by an apricot rectangle where the width represents the span of the observation and the height shows its uncertainty.  \nTsygankov et al. (2013) also report that this source was undetected, with a 1-sigma upper limit of 1.6 mCrab (2 . 3 × 10 -11 erg cm -2 s -1 ), in the INTEGRAL/Imager on-Board the INTEGRAL Satellite (IBIS; Ubertini et al. 2003) nine-year Galactic plane survey (Krivonos et al. 2012). Correcting for the distance to NGC 2617, this flux upper limit corresponds to a upper limit on the luminosity of L < 9 . 0 × 10 42 erg s -1 .  \nFinally, this source is also not present in the 70-month Swift /Burst Alert Telescope (BAT; Barthelmy et al. 2005) survey (Baumgartner et al. 2013). Conservatively scaling from Figure 1 in Baumgartner et al. (2013) for the position of NGC 2617, we see that from Dec. 2004 - Sept. 2010 the exposure time by the 70-month Swift /BAT survey was > 10 megaseconds. Then from Figure 11 in Baumgartner et al. (2013) we see that this exposure time leads to a 5-sigma upper limit on the hard X-ray flux (14-195 keV) of < 0 . 6 mCrab (8 . 5 × 10 -12 erg cm -2 s -1 ). This upper limit is shown on Figure 4. Correcting for the distance to NGC 2617, this flux upper limit corresponds to a upper limit on the luminosity of L < 3 . 4 × 10 42 erg s -1 .', '2.6. Radio': 'On 2013 Jun. 7, Yang et al. (2013) obtained a 7-hour European VLBI Network ToO observation of NGC 2617 at 1.6 GHz and detected a coincident compact source with a flux density of 1 . 5 ± 0 . 3 mJy, corresponding to a radio luminosity of L 1 . 6 = 4 . 8 × 10 37 erg s -1 . From these observations, Yang et al. (2013) constrain the source size to be less than 4 milliarcseconds, corresponding to a physical size of < 1 . 2 pc at NGC 2617.  \nOn 2013 Jun. 29, Jencson et al. (2013) conducted VLBA observations of the nucleus of NGC 2617 at 5 and 1.7 GHz. We detected radio emission with flux densities of 1 . 5 ± 0 . 1 and 1 . 6 ± 0 . 3 mJy at 5 and 1.7 GHz, respectively. These measurements indicate little to no variability since the earlier e-EVN observations (Yang et al. 2013), and are consistent with a flat spectrum for the nucleus of NGC 2617. Jencson et al.  \nFig. 5.- Spectroscopic observations of NGC 2617 spanning ∼ 70 days from 2013 Apr. 25 - Jun. 30, 2013. Details of the observations are giving in Table 2. Individual spectra are calibrated by the flux in their narrow [O III] λ 5007 lines. Spectra marked with an asterisk are those employed in the spectral analysis presented in §3.1. The UH spectrum, excluded from our analysis, has some calibration problems leading to distorted line ratios.  \n<!-- image -->  \n(2013) may have resolved the source at 5 GHz, finding major and minor axes of 2 . 6 ± 0 . 6 and 1 . 5 ± 0 . 2 mas, respectively. These correspond to physical sizes of 0 . 81 ± 0 . 19 and 0 . 47 ± 0 . 06 pc at NGC 2617.', '3. Analysis': 'We use our time-series spectra to search for variations in the broad H β and He II λ 4686 lines and to estimate the central BH mass in §3.1. We cross-correlate our photometric light-curves in §3.2. Finally, we construct a spectral energy distribution (SED) of NGC 2617 from the hard X-rays through the radio in §3.3.', '3.1. Emission-Line Variability and Black Hole Mass Estimate': "Accurate flux calibration is necessary to analyze intrinsic variations in AGN spectra. We therefore can only search for variability using the six well-calibrated, high signal-to-noise ratio (SNR) spectra that were collected under similar observing conditions (e.g., seeing, airmass; Peterson et al. 1995). These six spectra are marked with asterisks in Table 2 and Figure 5. We restricted our analysis to the H β emission-line region, since this spectral region is covered by all six spectra. The region near H α is either outside the wavelength range or is located near the wavelength limits where the flux calibration is less accurate. Although the spectra in this subset were taken under similar observing conditions, there are still differences in the apertures and seeing that will lead to spectral differences. However, under the assumption that the [O III] narrow-line flux is constant over the time scale of interest (e.g., Peterson et al. 2013) and that the [O III] emitting region is not resolved, we can calibrate the spectra onto the same absolute scale using the observed [O III] λ 5007 flux. The same is not true for the AGN continuum, which is contaminated by the emission of the host galaxy and, therefore, will still be subject to aperture effects.  \nWe degraded the six spectra to a common resolution of 9 Å (FWHM ). We adopted the spectrum taken on 2013 May 12.13 (under clear conditions) as the reference spectrum and inter-calibrated the spectra using the χ 2 minimization algorithm of van Groningen & Wanders (1992), under the assumption that the [O III] λ 5007 flux is constant. From the inter-calibrated spectra, we then computed the mean and the rootmean-square (RMS) residual spectra shown in Figure 6. The SNR of the RMS spectrum is very low due to the small number of spectra in the sample and the relatively small flux variations. Although a detailed analysis of the line profiles is precluded by the low SNR of the spectrum, both H β and He II changed significantly during the period of interest.  \nFrom the RMS spectrum, we can estimate the width of the variable part of the emission lines. Using a Monte Carlo approach, we generated 1000 sets of simulated spectra with randomized fluxes (see Peterson et al. 2004). We then created an RMS spectrum for each simulated data set, and measured the FWHMand line dispersion ( σ line ) following Peterson et al. (2004). Based on these procedures we find that the FWHM of the H β and He II λ 4686 lines are 5490 ± 1140kms -1 and 10050 ± 2590kms -1 and their velocity dispersions are 2200 ± 230kms -1 and 4180 ± 450kms -1 , respectively. The uncertainties are relatively large due to the low SNR of the RMS spectrum.  \nTable 2. Spectroscopic Observations  \nNote. - Observational and derived properties of our time series spectra of NGC 2617. Seeing gives the FWHM of the spatial profile in the 2d spectra. The airmass is reported for the beginning of the exposures. Spectra with UT Date marked with a * are used in the spectral analysis presented in §3.1.  \nFig. 6.- Mean (red solid, right flux scale) and RMS (blue dashed, left flux scale) spectra of NGC 2617 near the H β and He II λ 4686 broad-lines. Note that the RMS spectrum flux scale is 10 times lower than for the mean spectrum. The RMS spectrum H β and He II λ 4686 FWHM velocities are 5490 ± 1140kms -1 and 10050 ± 2590kms -1 , respectively, where the velocity dispersions are σ = 2200 ± 230kms -1 and σ = 4180 ± 450kms -1 , respectively.  \n<!-- image -->  \nThe width of the H β emission line and the mean continuum luminosity can then be used to estimate the mass of the black hole ( M BH ) powering the AGN. The width of the broad emission lines should primarily be due to the orbital motion of the emitting gas in the gravitational potential of the BH. Under this 'virial' assumption, M BH can be estimated by  \nM BH = f G -1 R BLR ∆ v 2 , (1)  \nwhere ∆ v is the velocity dispersion of the emitting gas estimated from the width of the emission line, R BLR is the distance of the broad line emitting gas from the central black hole, G is the gravitational constant, and f is a dimensionless scale factor that depends on inclination, geometry and kinematics of the BLR (e.g., Onken et al. 2004). Reverberation mapping studies of the H β emission line in local AGN indicate that there is a tight correlation between the AGN continuum luminosity at 5100 Å and R BLR (Kaspi et al. 2005; Bentz et al. 2009, 2013; Zu et al. 2011). Although we do not have enough spectral epochs to perform a reverberation mapping measurement, we can still estimate R BLR through the radius-luminosity relation. Measuring λ L λ (5100Å) from the mean spectrum and using the R BLR -λ L λ (5100Å) relation from Bentz et al. (2013), we obtain R BLR =9 ± 1 light days. Contamination by the host galaxy to the continuum luminosity can lead to an overestimate of the BLR radius. However, we estimate in §4 that the V -band contamination from the host galaxy in a 5 '' radius aperture is /similarequal 30% of the flux. Since the spectroscopic apertures are smaller and M BH ∝ L 1 / 2 our R BLR estimate is little affected by host galaxy contamination. Assuming f =4 . 31 ± 1 . 05 (Grier et al. 2013a) and using the H β σ line estimate from the RMS spectrum, we obtain M BH = (4 ± 1) × 10 7 M /circledot . For this estimate of the black hole mass, the Eddington luminosity is L Edd = 10 45 . 7 erg s -1 .", '3.2. Wavelength-Dependent Variability Lags and Temporal Smoothing': "Visual inspection of Figures 3 and 4 clearly show a lag between the X-ray and UV-NIR variability. It also seems reasonably clear that the longer-wavelength (e.g., NIR) light curves are also smoother than the shorter-wavelength light curves (e.g., UV). In this subsection, we quantify these lags and smoothing time scales. A large number of our observations were obtained in filters that differ only slightly in wavelength. As a first step we combined the u and U swift light curves; the g , B , and B swift light curves; the V and V swift light curves; the R and r light curves; and the I and i light curves to improve the overall time sampling for each wavelength range. To do this, we have added a constant to each band to match it to the partner band with the largest number of epochs. Each of the resulting light curves are shown in Figure 7 in arbitrary flux units and they are scaled to emphasize the light-curve structure. Figure 7 further emphasizes the temporal ordering of the variability.  \nTo quantitatively estimate the time lags and the temporal smoothing, we use the reverberation mapping and light curve analysis package JAVELIN 23 (formerly known as SPEAR; Zu et al. 2011, 2013b,a). As with a normal cross-correlation analysis, JAVELIN assumes all emission line light curves are scaled and shifted versions of the continuum light curve. However, JAVELIN differs from simple cross-correlation methods  \nFig. 7.- Scaled NIR-UV photometric observations of NGC 2617 spanning ∼ 60 days from 2013 Apr. 24 through June 20. Similar filters have been scaled and combined into a common light curve. Error bars are drawn but they are usually smaller than the points. The coral vertical lines from 6415-6445 indicate the epochs where SEDs are shown in Figure 10 and aid in seeing the lag from shorter to longer wavelengths. Black lines show light curves of a simple X-ray irradiation model described in §4. While not a perfect match, the general amplitudes, lags, and smoothness of the light curves as a function of wavelength are successfully reproduced. Finally, note how the relative strengths of the two peaks in the X-ray light curve change in the redward filters. This behavior is not captured in our our simple model with instantaneous reradiation.  \n<!-- image -->  \nin two respects. First, JAVELIN explicitly builds a model of the light curve and transfer function and fits it to both the continuum and the line data, maximizing the likelihood of the model and then computing uncertainties using the (Bayesian) Markov Chain Monte Carlo (MCMC) method. Second, as part of this process it models the continuum light curve using a damped random walk (DRW) model, which has been demonstrated to be a good statistical model of optical quasar variability using large ( ∼ 10 4 ) samples of quasar light curves (e.g., Kelly et al. 2009; Kozłowski et al. 2010; MacLeod et al. 2010; Zu et al. 2013b). The parameters of the DRW model are included in the fits and their uncertainties, as is a simple top-hat model of the transfer function and the light curve means (or trends if desired).  \nFollowing Zu et al. (2011), we apply the JAVELIN method to each of the light curves shown in Figure 7, using either the Swift XRTX-ray or the Swift UVOT UVW 2 band light curve as the reference ('continuum') light curve and the other light curves at longer wavelengths are the responding ('line') light curve. We fit fluxes rather than the magnitudes, which makes the results independent of any contamination by the host galaxy since the mean flux is subtracted from the light curves as part of the analysis. Also, to avoid a spurious peaks at ∼ 15 days, the spacing between the distinct two peaks visible in light curves, we forced JAVELIN to only return lags of < 10 days. Finally, since X-ray variability typically has a different power spectrum than the UV/optical bands (Cameron et al. 2012), modeling the X-ray light curve as a DRW process may underestimate the uncertainties on the time lag of the paired light curves. However, the estimate of the time lag is determined mainly by the variability features manifestly present in the light curves, so the result should be insensitive to the statistics assumed for the underlying variability models.  \nThe measured time lags are presented in Table 3. The upper panel (lower panel) of Figure 8 shows the measured time lags between the Swift 's UVOT uvw 2-band ( Swift 's XRT X-ray) flux and the redward photometric bands' flux as a function of wavelength. There is a clear trend that redder bands have larger lags, with the UV and NIR lagging 2 -3 and 6 -9 days behind the X-rays, respectively. We then fit a power law and a power law with a zero-point offset to the X-ray lags. The best-fits are  \nt lag = (0 . 38 ± 0 . 19) × ( λ nm ) 0 . 37 ± 0 . 08 days (2)  \nand  \nt lag = ( (2 . 43 ± 0 . 46) + (2 . 5 ± 0 . 7) × ( λ µ m ) 1 . 18 ± 0 . 33 ) days (3)  \nwhere t lag is the time lag as a function of wavelength ( λ ), respectively. The chi-square per degree of freedom for the fits are 1 . 48 and 0 . 97, respectively. The F -test confidence level for the addition of the temporal offset parameter is 94 . 2%. The zero-point of the second fit (2 . 43 ± 0 . 46 days) is in agreement with the non-geometric delay required by the simple quantitative model presented in §4 to reproduce the observed UV-NIR light curves. Furthermore, depending on whether dissipation or irradiation dominates, we would expect t lag ∝ λ α with 4 / 3 ≤ α ≤ 2 (see Equation 4). The results for Equation 3 are consistent with the shallower slopes where dissipation remains the dominant heating mechanism.  \nThe JAVELIN models also smooth the reference light curves with a top hat filter, and the wavelengthdependent top hat widths are presented in Table 3 and shown in Figure 9. The upper panel (lower panel) of  \nFig. 8.Upper panel: ( Lower panel: ) Measured time lags between variations in the Swift 's UVOT uvw 2band ( Swift 's XRT X-ray) flux and the redder photometric bands as a function of wavelength. There is a clear trend that redder bands have larger lags, with the UV, optical, and NIR lagging 2 -3, 3 -6, and 6 -9 days behind the X-rays, respectively. The black dashed and solid lines show the best-fit power law and best-fit power law with a zero-point offset, respectively.  \n<!-- image -->  \nFigure 9 shows the width of the transfer function between the Swift 's UVOT uvw 2-band ( Swift 's XRT X-ray) flux and the redward photometric bands' flux as a function of wavelength. Even though the widths become increasingly uncertain at longer wavelengths, there is still a clear trend that redder bands have larger widths, and thus smoother light curves. If the X-ray to UV lag and transfer function width measurements are robust, they have interesting implications for the origin of variability in AGN which are discussed in §4.", '3.3. Spectral Energy Distribution': 'Using the combined fluxes from §3.2, the Integral 17 -60 keV observation discussed in §2.5, and the EVN and VLBA radio observations discussed in §2.6 we produce SEDs of NGC 2617. To put the X-ray observations on this SED we assigned an effective wavelength for the Swift X-ray band, using the measured energy spectral index, and for the Integral X-ray band, assuming a flat energy spectral index. We then multiplied the measured X-ray fluxes by this effective wavelength and divided by the wavelength range integrated in each measurement. Figure 10 shows the X-ray through NIR SED of NGC 2617 at a number of epochs in the left panel and a hard X-ray through radio SED in the right panel. For the right panel of Figure 10, we interpolated the X-ray through NIR light curves on the JD= 2 , 456 , 420 . 0, we show the Integral 17 -60 keV X-ray observation which was obtained during this epoch, and we show the EVN and VLBA observations taken 30 and 52 days later, respectively. During the beginning of the flare there was a marked increase in the X-ray and UV flux while the increase at longer wavelengths is less pronounced.  \nIn the right panel of Figure 10 we also show the composite SED for radio-quiet quasars of Shang et al. (2011) scale to match in the Swift UV fluxes. We normalized the SED to the UV luminosity to minimize the effects of host-galaxy contamination which, given the shape of the SED in the optical-NIR, is large. In §4 we include a model of this contamination and find it is roughly consistent with the SED of an early type spiral. Finally, we see that NGC 2617 was extremely X-ray bright during the flare, being 10-25 times brighter than the composite SED.  \nIn Figure 10, we also see that the SED roughly has ν L ν ∼ 0 . 01 L Edd , although there is significant host contamination as we move to the redder optical and NIR bands (see §4 and Figure 11). Assuming the existence of an unobserved hard UV accretion peak, this suggests that the AGN flare peaked at ∼ 10% of L Edd . This estimate is consistent estimates ( ∼ 0.03-0.15) obtained using Equation 1 from Grupe et al. (2010) and the measured X-ray spectral indices reported in Table 1 and shown in Figure 4.', '4. A Simple Quantitative Model': "Since the sequence of variability seems to lead with the X-rays, followed by the ultraviolet and eventually the NIR, the natural explanation of the observations seems to be variability driven by X-ray irradiation of the BH accretion disk. We know from gravitational microlensing of lensed X-ray emission from quasars (e.g., Mosquera et al. 2013) that the X-ray emission has a half-light radius of order ten times the gravitational radius ( r g ). This is small enough, 10 r g / c /similarequal 0 . 5 hours, such that we can treat the X-ray emission as a point  \nFig. 9.Upper panel: ( Lower panel: ) Measured widths for the transfer function between the Swift 's UVOT uvw 2-band ( Swift 's XRT X-ray) flux and the redder photometric bands as a function of wavelength. There is a reasonably clear trend that redder bands have wider transfer functions and thus smoother light curves.  \n<!-- image -->  \nTable 3. Measured Time Lags and Smoothing Scales  \nFig. 10.- SED of NGC 2617. Left panel: SED as a function of time as determined from our combined light curves shown in Figure 7. Right panel: SED stretching from hard X-rays to Radio. X-ray through NIR observations were interpolated on the JD date 2,456,420.0. The Integral 17 -60 keV X-ray observation was obtained during this epoch. Radio observations shown were obtained 30-52 days later. The black line shows the composite SED for radio-quiet quasars from Shang et al. (2011) scaled to match in the UV. Fractions of L Edd assuming log M BH / M /circledot = 7 . 6 are shown for reference. The left panel shows that the X-ray and UV fluxes rapidly strengthen at the beginning of the flare. The right panel shows that NGC 2617 was extremely X-ray bright during the flare, NGC 2617 is radio quiet, and our SED contains a significant amount of contamination from the host galaxy in the optical and NIR wavelength ranges (see Figure 11).  \n<!-- image -->  \n<!-- image -->  \nsource. Therefore, we consider a simple 'lamp post' model (e.g., Berkley et al. 2000, Kazanas & Nayakshin 2001) for the NGC 2617 outburst which consists of a flat, non-relativistic thin disk with an accretion luminosity L A and an inner disk edge at R in = 3 r g illuminated by an X-ray point source a distance h X above the disk center.  \nIn this model we assume that all X-rays incident on the disk are absorbed and instantaneously heat the disk surface. In this case the local disk temperature at observed time t is  \nσ T ( R , t ) 4 = 3 L A R in 2 π R 3 [ 1 -( R in R ) 1 / 2 ] + L X ( t -τ ( R )) 4 π ( R 2 + h 2 X ) h X ( R 2 + h 2 X ) 1 / 2 (4)  \nwhere L X ( t ) is the X-ray luminosity and τ ( R ) = τ 0 + h X / c + ( R 2 + h 2 X ) 1 / 2 / c is (for simplicity) the time delay between light from radius R and the direct X-ray emission. To synthesize emission at wavelength λ we simply assume black body emission at each radius and integrate over the disk,  \nν F ν = 4 π h ν 4 c 2 D 2 ∫ ∞ Rin RdR exp( h ν/ kT ( R , t )) -1 . (5)  \nWe experimented with adding a disk opening angle and a bolometric correction for the X-ray fluxes, but these had little effect.  \nAs seen in §3.3, our fairly large aperture photometry contains significant host contamination. We experimented with several approaches to remove this contamination short of fully modeling the images. For example, the contamination is reasonably well modeled by the Sbc SED template from Assef et al. (2010). However, there are systematic offsets that still dominate the χ 2 of the resulting fits, which could either be due to the host or our use of a simple, non-relativistic thin disk model. Since such questions are peripheral to the temporal changes in the fluxes, we simply included an independent, constant 'host luminosity' for each band found by doing a χ 2 fit for fixed AGN parameters. We then experimented with the parameters for the accretion disk luminosity L A , the 'quiescent' X-ray luminosity L X 0 used where we lack Swift data, the X-ray source height h X and the non-geometric delay τ 0 keeping the black hole mass fixed at log( M BH / M /circledot ) = 7 . 6.  \nGiven the simplicity of the model, we can obtain surprisingly good (but not perfect) fits. First, Figure 11 shows the SED of the 'host galaxy' contamination model as compared to the SED templates of Assef et al. (2010). Although we did not impose any prior that these estimates resemble the SED of a galaxy, we see that the estimates are broadly consistent with the Sbc template except near the H α line. Because the broad lines should respond less and more slowly than the continuum, it is not surprising that they are absorbed into the host-galaxy parameters. The match to the galaxy templates (or combinations of them) is not perfect due to both observational issues (e.g. seeing variations between observations will lead to variations in the aperture flux of the host) and any shortcomings in our very simple disk model. There is also a tendency to over-estimate the host flux at the longest wavelengths.  \nFigure 7 shows the model UV, optical and NIR light curves compared with the observed light curves. The model naturally reproduces the patterns seen in the data, namely that the shorter wavelengths have shorter lags and are less smooth than longer wavelengths, simply because of geometry. Shorter wavelength  \nemission comes from smaller radii with shorter light crossing times to the X-ray source, leading to the correlation of the lags with wavelength. Furthermore, the size of the region emitting the shorter wavelengths is smaller, so there is less temporal smoothing of the X-ray emission. This is somewhat exaggerated at the longer wavelengths because the host galaxy contribution also reduces the fractional variability, which is why the analysis in §3.2 was carried out using fluxes rather than magnitudes. For example, the tendency to overshoot the NIR host flux leads to the less than observed NIR variability seen in Figure 7.  \nThe best fit model has log L A =42 . 6 and log L X 0 =42 . 6 (in ergs s -1 ) which is 10 -3 . 1 L Edd for the estimated black hole mass. The X-ray source can be at a very low altitude, log( h x / cm) = 13 . 5 which is 4 . 8 r g where r g = GM BH / c 2 is the gravitational radius of the black hole, consistent with the results from microlensing (e.g. Mosquera et al. 2013). We do not think an attempt to assign uncertainties to the model parameters is presently warranted.  \nThe three main weaknesses of our model's fit to the data are: 1) The model NIR light curves provide a poor fit to the observed data. This could be due to a number of reasons. First, as mentioned above, accounting for the host galaxy flux is especially important in the NIR where the galaxy contribution is large and our over estimation of the host galaxy's contribution at the longest wavelength reduces the expected variability. Second, our model does not include any dust reprocessing. If there is significant reprocessing of UV and optical photons on scales larger than the BLR, then the NIR would respond to the increase in X-ray flux on two distinct timescales, whereas only one is currently included in our models. 2) We do not reproduce the second bump seen in the X-ray light curve around JD=2,456,451 in the UV-NIR light curves. The reason for this is unclear. 3) We were unable to find successful solutions without the non-geometric lag τ 0 /similarequal 2 days. Trying to make it part of the geometric lag by putting the X-ray source at h x /similarequal τ 0 / c above the disk did not work and would be an odd geometry in any case, with the X-ray emission arising from ∼ 10 3 r g above the disk. The basic problem is that the UV emission arises from such small radii that there is no geometric lag relative to a central X-ray source, so that portion of the delayed response has to be supplied by τ 0 . We do not have a physical interpretation for this non-geometric lag.", '5. Summary': "On 2013 Apr. 10.27, ASAS-SN was triggered by a ∼ 10% relative flux increase from the inner regions of NGC 2617. Follow-up spectroscopy found that NGC 2617, previously classified as a Seyfert 1.8 galaxy, had also changed its spectral type to become a Seyfert 1 with broad optical emission lines and a continuum blue bump. We then began a ∼ 70 day photometric and spectroscopic monitoring campaign from the NIR X-rays, observing NGC 2617 through a dramatic outburst in which its X-ray flux increased by over an order of magnitude followed by an increase in its optical/ultraviolet (UV) continuum flux by almost an order of magnitude. Although these spectral changes could have occurred at any time in the intervening decade, it seems most likely that the change in Seyfert type was associated with the strong outburst given the rarity of both 'changing-look' AGN and powerful X-ray flares.  \nWhen analyzing these observations we found significant variations in the broad emission lines H β and  \nFig. 11.- The SED for the 'host galaxy' contamination in the best fit model as compared to the elliptical (solid), Sbc (dashed) and irregular (dotted) galaxy templates from Assef et al. (2010). The host is well matched to an SED to the Sbc template, with deviations near H α .  \n<!-- image -->  \nHeII λ 4686. When cross-correlating our X-ray light curve and the other photometric light-curves we found that the disk emission lags the X-rays, with the lag becoming longer as we move from the UV (2 -3 days) to the NIR (6 -9 days). We also see that the NIR is more heavily temporally smoothed than the UV. Since the sequence of variability seems to lead with the X-rays and is followed by redward bands, the natural explanation of the observations seems to be variability driven by X-ray irradiation of the disk, possibly due to coronal activity. A simple physical model using X-ray illumination to heat a thin disk roughly reproduces the UV-NIR light curves. That such a simple model works at all strongly suggests that a more detailed and realistic modeling study is warranted.  \nWe believe this is the first case where there is a clear and convincing determination of whether X-ray variability drives UV → NIR variability or the reverse. In this outburst, the X-ray variability unambiguously drives the UV → NIR variability, almost certainly by irradiating and heating the disk. While quasar and AGN variability is ubiquitous, it has proved difficult to measure time delays between thermal disk and nonthermal X-ray emission or between different wavelengths of the thermal emission, leaving the origin of the variability unclear. We can now quantify quasar variability well, using the DRW model (e.g. Kelly et al. 2009, Kozłowski et al. 2010, MacLeod et al. 2010, MacLeod et al. 2012, Zu et al. 2013b) and the parameters of the DRW model are correlated with properties of the AGN. Determining the physical origins of these parameters and their scalings has been more difficult. For example, Kelly et al. (2009) argue that the DRW time scale could be related to the thermal time scale of the disk, while MacLeod et al. (2010) note that the scalings of both the DRW time scale and amplitude with wavelength, black hole mass, and luminosity all seem too weak compared to the natural scalings one would expect for a simple disk model. A program systematically triggering high cadence Swift monitoring campaigns after optically detecting a significant luminosity change can likely provide similar results to this case for larger numbers of systems to test whether the pattern we saw here is unique or ubiquitous, and the time ordering of the variability should lead to a clear physical model for understanding quasar variability more generally. Ideally, one should also try to have far better spectroscopic sampling than we achieved here so that one simultaneously carries out a reverberation mapping campaign to directly determine the size of the broad line region and black hole mass.  \nWe thank the APO director, S. Hawley, for granting us DD Time for this observation, and the observatory crew. Additionally, we would like to thank N. Gehrels for approving our Swift ToO observation requests and the Swift ToO team for promptly scheduling and executing these observations. We also thank Jennifer van Saders, Todd Thompson, Michael Fausnaugh, and Joe Antognini for discussions and encouragement. We also thank the anonymous referee for his/her helpful comments and suggestions improving this manuscript. B.J.S. was supported by a Graduate Research Fellowship from the National Science Foundation. Development of ASAS-SN has been supported by NSF grant AST-0908816 and the Center for Cosmology and AstroParticle Physics at The Ohio State University. Swift at PSU is supported by NASA contract NAS5-00136. C.S.K. and B.J.S are supported by NSF grant AST-1009756. G.D.R, B.M.P, and R.W.P. are supported by NSF grant AST-1008882. S.K. would like to thank the Kavli Institute for Theoretical Physics (KITP) for their hospitality and support during the program on 'A Universe of Black Holes'. S.K. is supported by NSF grant NSF PHY11-25915. M.I., C.C., and Y.Y. acknowledge the support from the Creative Initiative program No. 2008-0060544 of the National Research Foundation of Korea (NRFK)  \nfunded by the Korea government (MSIP). J.F.B. is supported by NSF grant PHY-1101216. This paper uses the data taken by the United Kingdom Infrared Telescope which is operated by the Joint Astronomy Centre on behalf of the Science and Technology Facilities Council of the U.K. This paper includes data gathered with the 6.5 meter Magellan Telescopes located at Las Campanas Observatory, Chile. The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research has made use of NASA's Astrophysics Data System Bibliographic Services. This paper makes use of data from the AAVSO Photometric All Sky Survey, whose funding has been provided by the Robert Martin Ayers Sciences Fund.", 'REFERENCES': "Ahn C. P., Alexandroff R., Allende Prieto C., Anderson S. F., Anderton T., Andrews B. H., Aubourg É., Bailey S., Balbinot E., Barnes R., et al., 2012, ApJS, 203, 21  \nAldering G., Antilogus P., Bailey S., Baltay C., Bauer A., Blanc N., Bongard S., Copin Y., Gangler E., Gilles S., Kessler R., Kocevski D., Lee B. C., Loken S., Nugent P., Pain R., Pécontal E., Pereira R., Perlmutter S., Rabinowitz D., Rigaudier G., Scalzo R., Smadja G., Thomas R. C., Wang L., Weaver B. A., Nearby Supernova Factory, 2006, ApJ, 650, 510  \nAretxaga I., Joguet B., Kunth D., Melnick J., Terlevich R. J., 1999, ApJ, 519, L123  \nAssef R. J., Kochanek C. S., Brodwin M., Cool R., Forman W., Gonzalez A. H., Hickox R. C., Jones C., Le Floc'h E., Moustakas J., Murray S. S., Stern D., 2010, ApJ, 713, 970  \nBarthelmy S. D., Barbier L. M., Cummings J. R., Fenimore E. E., Gehrels N., Hullinger D., Krimm H. A., Markwardt C. B., Palmer D. M., Parsons A., Sato G., Suzuki M., Takahashi T., Tashiro M., Tueller J., 2005, Space Sci. Rev., 120, 143  \nBaumgartner W. H., Tueller J., Markwardt C. B., Skinner G. K., Barthelmy S., Mushotzky R. F., Evans P. A., Gehrels N., 2013, ApJS, 207, 19  \nBentz M. C., Denney K. D., Grier C. J., Barth A. J., Peterson B. M., Vestergaard M., Bennert V. N., Canalizo G., De Rosa G., Filippenko A. V., Gates E. L., Greene J. E., Li W., Malkan M. A., Pogge R. W., Stern D., Treu T., Woo J.-H., 2013, ApJ, 767, 149  \nBentz M. C., Peterson B. M., Netzer H., Pogge R. W., Vestergaard M., 2009, ApJ, 697, 160  \nBerkley A. J., Kazanas D., Ozik J., 2000, ApJ, 535, 712  \nBoller T., Meurs E. J. A., Brinkmann W., Fink H., Zimmermann U., Adorf H.-M., 1992, A&A, 261, 57  \nBrandt W. N., Pounds K. A., Fink H., 1995, MNRAS, 273, L47  \nBreedt E., Arévalo P., McHardy I. M., Uttley P., Sergeev S. G., Minezaki T., Yoshii Y., Gaskell C. M., Cackett E. M., Horne K., Koshida S., 2009, MNRAS, 394, 427  \nBreeveld A. A., Curran P. A., Hoversten E. A., Koch S., Landsman W., Marshall F. E., Page M. J., Poole T. S., Roming P., Smith P. J., Still M., Yershov V., Blustin A. J., Brown P. J., Gronwall C., Holland S. T., Kuin N. P. M., McGowan K., Rosen S., Boyd P., Broos P., Carter M., Chester M. M., Hancock B., Huckle H., Immler S., Ivanushkina M., Kennedy T., Mason K. O., Morgan A. N., Oates S., de Pasquale M., Schady P., Siegel M., vanden Berk D., 2010, MNRAS, 406, 1687  \nBrown T. M., Baliber N., Bianco F. B., Bowman M., Burleson B., Conway P., Crellin M., Depagne É., De Vera J., Dilday B., Dragomir D., Dubberley M., Eastman J. D., Elphick M., Falarski M., Foale S., Ford M., Fulton B. J., Garza J., Gomez E. L., Graham M., Greene R., Haldeman B., Hawkins E., Haworth B., Haynes R., Hidas M., Hjelstrom A. E., Howell D. A., Hygelund J., Lister T. A., Lobdill R., Martinez J., Mullins D. S., Norbury M., Parrent J., Paulson R., Petry D. L., Pickles A., Posner V., Rosing W. E., Ross R., Sand D. J., Saunders E. S., Shobbrook J., Shporer A., Street R. A., Thomas D., Tsapras Y., Tufts J. R., Valenti S., Vander Horst K., Walker Z., White G., Willis M., 2013, PASP, 125, 1031  \nBurrows D. N., Hill J. E., Nousek J. A., Kennea J. A., Wells A., Osborne J. P., Abbey A. F., Beardmore A., Mukerjee K., Short A. D. T., Chincarini G., Campana S., Citterio O., Moretti A., Pagani C., Tagliaferri G., Giommi P., Capalbi M., Tamburelli F., Angelini L., Cusumano G., Bräuninger H. W., Burkert W., Hartner G. D., 2005, Space Sci. Rev., 120, 165  \nCackett E. M., Horne K., Winkler H., 2007, MNRAS, 380, 669  \nCameron D. T., McHardy I., Dwelly T., Breedt E., Uttley P., Lira P., Arevalo P., 2012, MNRAS, 422, 902  \nChakravorty S., Elvis M., Ferland G., 2014, MNRAS, 437, 740  \nCohen R. D., Puetter R. C., Rudy R. J., Ake T. B., Foltz C. B., 1986, ApJ, 311, 135  \nCollier S. J., Horne K., Kaspi S., Netzer H., Peterson B. M., Wanders I., Alexander T., Bertram R., Comastri A., Gaskell C. M., Malkov Y. F., Maoz D., Mignoli M., Pogge R. W., Pronik V. I., Sergeev S. G., Snedden S., Stirpe G. M., Bochkarev N. G., Burenkov A. N., Shapovalova A. I., Wagner R. M., 1998, ApJ, 500, 162  \nCzerny B., Hryniewicz K., 2011, A&A, 525, L8  \nDressler A., Bigelow B., Hare T., Sutin B., Thompson I., Burley G., Epps H., Oemler A., Bagish A., Birk C., Clardy K., Gunnels S., Kelson D., Shectman S., Osip D., 2011, PASP, 123, 288  \nGehrels N., Chincarini G., Giommi P., Mason K. O., Nousek J. A., Wells A. A., White N. E., Barthelmy S. D., Burrows D. N., Cominsky L. R., Hurley K. C., Marshall F. E., Mészáros P., Roming P. W. A., Angelini L., Barbier L. M., Belloni T., Campana S., Caraveo P. A., Chester M. M., Citterio O., Cline T. L., Cropper M. S., Cummings J. R., Dean A. J., Feigelson E. D., Fenimore E. E., Frail D. A., Fruchter A. S., Garmire G. P., Gendreau K., Ghisellini G., Greiner J., Hill J. E., Hunsberger S. D., Krimm H. A., Kulkarni S. R., Kumar P., Lebrun F., Lloyd-Ronning N. M., Markwardt C. B., Mattson B. J., Mushotzky R. F., Norris J. P., Osborne J., Paczynski B., Palmer D. M., Park H.-S., Parsons A. M., Paul J., Rees M. J., Reynolds C. S., Rhoads J. E., Sasseen T. P., Schaefer B. E., Short A. T., Smale A. P., Smith I. A., Stella L., Tagliaferri G., Takahashi T., Tashiro M., Townsley L. K., Tueller J., Turner M. J. L., Vietri M., Voges W., Ward M. J., Willingale R., Zerbi F. M., Zhang W. W., 2004, ApJ, 611, 1005  \nGoodrich R. W., 1989, ApJ, 340, 190  \nGregory S. A., Tifft W. G., Cocke W. J., 1991, AJ, 102, 1977  \nGrier C. J., Martini P., Watson L. C., Peterson B. M., Bentz M. C., Dasyra K. M., Dietrich M., Ferrarese L., Pogge R. W., Zu Y., 2013a, ApJ, 773, 90  \nGrier C. J., Peterson B. M., Horne K., Bentz M. C., Pogge R. W., Denney K. D., De Rosa G., Martini P., Kochanek C. S., Zu Y., Shappee B., Siverd R., Beatty T. G., Sergeev S. G., Kaspi S., Araya Salvo C., Bird J. C., Bord D. J., Borman G. A., Che X., Chen C., Cohen S. A., Dietrich M., Doroshenko V. T., Efimov Y. S., Free N., Ginsburg I., Henderson C. B., King A. L., Mogren K., Molina M., Mosquera A. M., Nazarov S. V., Okhmat D. N., Pejcha O., Rafter S., Shields J. C., Skowron J., Szczygiel D. M., Valluri M., van Saders J. L., 2013b, ApJ, 764, 47  \nGrier C. J., Peterson B. M., Pogge R. W., Denney K. D., Bentz M. C., Martini P., Sergeev S. G., Kaspi S., Minezaki T., Zu Y., Kochanek C. S., Siverd R., Shappee B., Stanek K. Z., Araya Salvo C., Beatty T. G., Bird J. C., Bord D. J., Borman G. A., Che X., Chen C., Cohen S. A., Dietrich M., Doroshenko V. T., Drake T., Efimov Y. S., Free N., Ginsburg I., Henderson C. B., King A. L., Koshida S., Mogren K., Molina M., Mosquera A. M., Nazarov S. V., Okhmat D. N., Pejcha O., Rafter S., Shields J. C., Skowron J., Szczygiel D. M., Valluri M., van Saders J. L., 2012a, ApJ, 755, 60  \nGrier C. J., Peterson B. M., Pogge R. W., Denney K. D., Bentz M. C., Martini P., Sergeev S. G., Kaspi S., Zu Y., Kochanek C. S., Shappee B. J., Stanek K. Z., Araya Salvo C., Beatty T. G., Bird J. C., Bord D. J., Borman G. A., Che X., Chen C., Cohen S. A., Dietrich M., Doroshenko V. T., Efimov Y. S., Free N., Ginsburg I., Henderson C. B., Horne K., King A. L., Mogren K., Molina M., Mosquera A. M., Nazarov S. V., Okhmat D. N., Pejcha O., Rafter S., Shields J. C., Skowron J., Szczygiel D. M., Valluri M., van Saders J. L., 2012b, ApJ, 744, L4  \nGrupe D., Beuerman K., Mannheim K., Thomas H.-C., Fink H. H., de Martino D., 1995a, A&A, 300, L21  \nGrupe D., Beuermann K., Mannheim K., Bade N., Thomas H.-C., de Martino D., Schwope A., 1995b, A&A, 299, L5  \n- Grupe D., Komossa S., Gallo L. C., 2007, ApJ, 668, L111\n- Grupe D., Komossa S., Gallo L. C., Fabian A. C., Larsson J., Pradhan A. K., Xu D., Miniutti G., 2008, ApJ, 681, 982\n- Grupe D., Komossa S., Gallo L. C., Lia Longinotti A., Fabian A. C., Pradhan A. K., Gruberbauer M., Xu D., 2012, ApJS, 199, 28\n- Grupe D., Komossa S., Leighly K. M., Page K. L., 2010, ApJS, 187, 64\n- Grupe D., Komossa S., Scharwächter J., Dietrich M., Leighly K. M., Lucy A., Barlow B. N., 2013, AJ, 146, 78\n- Hill J. E., Burrows D. N., Nousek J. A., Abbey A. F., Ambrosi R. M., Bräuninger H. W., Burkert W., Campana S., Cheruvu C., Cusumano G., Freyberg M. J., Hartner G. D., Klar R., Mangels C., Moretti A., Mori K., Morris D. C., Short A. D. T., Tagliaferri G., Watson D. J., Wood P., Wells A. A., 2004, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 5165, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Flanagan K. A., Siegmund O. H. W., eds., pp. 217-231  \nHodgkin S. T., Irwin M. J., Hewett P. C., Warren S. J., 2009, MNRAS, 394, 675  \n- Jencson J., Kundert K., Mioduszewski A., Lucy A., Kadowaki J., Mellon S., 2013, The Astronomer's Telegram, 5347, 1\n- Jones D. H., Read M. A., Saunders W., Colless M., Jarrett T., Parker Q. A., Fairall A. P., Mauch T., Sadler E. M., Watson F. G., Burton D., Campbell L. A., Cass P., Croom S. M., Dawe J., Fiegert K., Frankcombe L., Hartley M., Huchra J., James D., Kirby E., Lahav O., Lucey J., Mamon G. A., Moore L., Peterson B. A., Prior S., Proust D., Russell K., Safouris V., Wakamatsu K.-I., Westra E., Williams M., 2009, MNRAS, 399, 683\n- Jones D. H., Saunders W., Colless M., Read M. A., Parker Q. A., Watson F. G., Campbell L. A., Burkey D., Mauch T., Moore L., Hartley M., Cass P., James D., Russell K., Fiegert K., Dawe J., Huchra J., Jarrett T., Lahav O., Lucey J., Mamon G. A., Proust D., Sadler E. M., Wakamatsu K.-i., 2004, MNRAS, 355, 747  \nKalberla P. M. W., Burton W. B., Hartmann D., Arnal E. M., Bajaja E., Morras R., Pöppel W. G. L., 2005, A&A, 440, 775  \n- Kaspi S., Maoz D., Netzer H., Peterson B. M., Vestergaard M., Jannuzi B. T., 2005, ApJ, 629, 61  \nKazanas D., Nayakshin S., 2001, ApJ, 550, 655  \nKelly B. C., Bechtold J., Siemiginowska A., 2009, ApJ, 698, 895  \nKhachikian E. Y., Weedman D. W., 1971, ApJ, 164, L109  \nMosquera A. M., Kochanek C. S., Chen B., Dai X., Blackburne J. A., Chartas G., 2013, ApJ, 769, 53  \nNandra K., Clavel J., Edelson R. A., George I. M., Malkan M. A., Mushotzky R. F., Peterson B. M., Turner T. J., 1998, ApJ, 505, 594  \nNicastro F., 2000, ApJ, 530, L65  \n- Onken C. A., Ferrarese L., Merritt D., Peterson B. M., Pogge R. W., Vestergaard M., Wandel A., 2004, ApJ, 615, 645  \nPaturel G., Theureau G., Bottinelli L., Gouguenheim L., Coudreau-Durand N., Hallet N., Petit C., 2003, A&A, 412, 57  \nPenston M. V., Perez E., 1984, MNRAS, 211, 33P  \nPeterson B. M., Denney K. D., De Rosa G., Grier C. J., Pogge R. W., Bentz M. C., Kochanek C. S., Vestergaard M., Kilerci-Eser E., Dalla Bontà E., Ciroi S., 2013, ApJ, 779, 109  \nPeterson B. M., Ferrarese L., Gilbert K. M., Kaspi S., Malkan M. A., Maoz D., Merritt D., Netzer H., Onken C. A., Pogge R. W., Vestergaard M., Wandel A., 2004, ApJ, 613, 682  \nPeterson B. M., Pogge R. W., Wanders I., Smith S. M., Romanishin W., 1995, PASP, 107, 579  \nPeterson B. M., Wanders I., Horne K., Collier S., Alexander T., Kaspi S., Maoz D., 1998, PASP, 110, 660  \nPiro L., Balucinska-Church M., Fink H., Fiore F., Matsuoka M., Perola G. C., Soffitta P., 1997, A&A, 319, 74  \nPoole T. S., Breeveld A. A., Page M. J., Landsman W., Holland S. T., Roming P., Kuin N. P. M., Brown P. J., Gronwall C., Hunsberger S., Koch S., Mason K. O., Schady P., vanden Berk D., Blustin A. J., Boyd P., Broos P., Carter M., Chester M. M., Cucchiara A., Hancock B., Huckle H., Immler S., Ivanushkina M., Kennedy T., Marshall F., Morgan A., Pandey S. B., de Pasquale M., Smith P. J., Still M., 2008, MNRAS, 383, 627  \nReynolds C. S., Nowak M. A., 2003, Phys. Rep., 377, 389  \nRisaliti G., 2011, Memorie della Societa Astronomica Italiana Supplementi, 17, 118  \nRomano P., Turner T. J., Mathur S., George I. M., 2002, ApJ, 564, 162  \nRoming P. W. A., Kennedy T. E., Mason K. O., Nousek J. A., Ahr L., Bingham R. E., Broos P. S., Carter M. J., Hancock B. K., Huckle H. E., Hunsberger S. D., Kawakami H., Killough R., Koch T. S., McLelland M. K., Smith K., Smith P. J., Soto J. C., Boyd P. T., Breeveld A. A., Holland S. T., Ivanushkina M., Pryzby M. S., Still M. D., Stock J., 2005, Space Sci. Rev., 120, 95  \nSaxton R. D., Read A. M., Esquej P., Freyberg M. J., Altieri B., Bermejo D., 2008, A&A, 480, 611  \nSergeev S. G., Doroshenko V. T., Golubinskiy Y. V., Merkulova N. I., Sergeeva E. A., 2005, ApJ, 622, 129  \nShang Z., Brotherton M. S., Wills B. J., Wills D., Cales S. L., Dale D. A., Green R. F., Runnoe J. C., Nemmen R. S., Gallagher S. C., Ganguly R., Hines D. C., Kelly B. J., Kriss G. A., Li J., Tang B., Xie Y., 2011, ApJS, 196, 2  \nShappee B. J., Grupe D., Mathur S., Prieto J. L., Kochanek C. S., Stanek K. Z., Jencson J., Talabere A., Beacom J. F., Szczygiel D., Pojmanski G., Dubberley M., Elphick M., Foale S., Hawkins E., Mullens D., Rosing W., Ross R., Walker Z., 2013a, The Astronomer's Telegram, 5059, 1  \nShappee B. J., Prieto J. L., Nugent J., Kochanek C. S., Stanek K. Z., Jencson J., Talabere A., Beacom J. F., Szczygiel D., Pojmanski G., Dubberley M., Elphick M., Foale S., Hawkins E., Mullens D., Rosing W., Ross R., Walker Z., 2013b, The Astronomer's Telegram, 5010, 1  \nShields G. A., 1978, Nature, 272, 706  \nSkrutskie M. F., Cutri R. M., Stiening R., Weinberg M. D., Schneider S., Carpenter J. M., Beichman C., Capps R., Chester T., Elias J., Huchra J., Liebert J., Lonsdale C., Monet D. G., Price S., Seitzer P., Jarrett T., Kirkpatrick J. D., Gizis J. E., Howard E., Evans T., Fowler J., Fullmer L., Hurt R., Light R., Kopan E. L., Marsh K. A., McCallon H. L., Tam R., Van Dyk S., Wheelock S., 2006, AJ, 131, 1163  \nStorchi-Bergmann T., Baldwin J. A., Wilson A. S., 1993, ApJ, 410, L11  \nTewes M., Courbin F., Meylan G., Kochanek C. S., Eulaers E., Cantale N., Mosquera A. M., Magain P., Van Winckel H., Sluse D., Cataldi G., Vörös D., Dye S., 2013, A&A, 556, A22  \nTohline J. E., Osterbrock D. E., 1976, ApJ, 210, L117  \nTran H. D., Osterbrock D. E., Martel A., 1992, AJ, 104, 2072  \nTruemper J., 1982, Advances in Space Research, 2, 241  \nTsygankov S., Krivonos R., Sazonov S., Lutovinov A., Sunyaev R., 2013, The Astronomer's Telegram, 5103, 1  \nUbertini P., Lebrun F., Di Cocco G., Bazzano A., Bird A. J., Broenstad K., Goldwurm A., La Rosa G., Labanti C., Laurent P., Mirabel I. F., Quadrini E. M., Ramsey B., Reglero V., Sabau L., Sacco B., Staubert R., Vigroux L., Weisskopf M. C., Zdziarski A. A., 2003, A&A, 411, L131  \nvan Dokkum P. G., 2001, PASP, 113, 1420  \nvan Groningen E., Wanders I., 1992, PASP, 104, 700  \nVeilleux S., Zheng W., 1991, ApJ, 377, 89  \nWanders I., Peterson B. M., Alloin D., Ayres T. R., Clavel J., Crenshaw D. M., Horne K., Kriss G. A., Krolik J. H., Malkan M. A., Netzer H., O'Brien P. T., Reichert G. A., Rodriguez-Pascual P. M., Wamsteker W., Alexander T., Anderson K. S. J., Benitez E., Bochkarev N. G., Burenkov A. N.,  \nCheng F.-Z., Collier S. J., Comastri A., Dietrich M., Dultzin-Hacyan D., Espey B. R., Filippenko A. V., Gaskell C. M., George I. M., Goad M. R., Ho L. C., Kaspi S., Kollatschny W., Korista K. T., Laor A., MacAlpine G. M., Mignoli M., Morris S. L., Nandra K., Penton S., Pogge R. W., Ptak R. L., Rodriguez-Espinoza J. M., Santos-Lleo M., Shapovalova A. I., Shull J. M., Snedden S. A., Sparke L. S., Stirpe G. M., Sun W.-H., Turner T. J., Ulrich M.-H., Wang T.-G., Wei C., Welsh W. F., Xue S.-J., Zou Z.-L., 1997, ApJS, 113, 69  \nWinkler C., Courvoisier T. J.-L., Di Cocco G., Gehrels N., Giménez A., Grebenev S., Hermsen W., MasHesse J. M., Lebrun F., Lund N., Palumbo G. G. C., Paul J., Roques J.-P., Schnopper H., Schönfelder V., Sunyaev R., Teegarden B., Ubertini P., Vedrenne G., Dean A. J., 2003, A&A, 411, L1  \nYang J., Paragi Z., Komossa S., van Bemmel I., Oonk R., 2013, The Astronomer's Telegram, 5125, 1  \nYork D. G., Adelman J., Anderson Jr. J. E., Anderson S. F., Annis J., Bahcall N. A., Bakken J. A., Barkhouser R., Bastian S., Berman E., Boroski W. N., Bracker S., Briegel C., Briggs J. W., Brinkmann J., Brunner R., Burles S., Carey L., Carr M. A., Castander F. J., Chen B., Colestock P. L., Connolly A. J., Crocker J. H., Csabai I., Czarapata P. C., Davis J. E., Doi M., Dombeck T., Eisenstein D., Ellman N., Elms B. R., Evans M. L., Fan X., Federwitz G. R., Fiscelli L., Friedman S., Frieman J. A., Fukugita M., Gillespie B., Gunn J. E., Gurbani V. K., de Haas E., Haldeman M., Harris F. H., Hayes J., Heckman T. M., Hennessy G. S., Hindsley R. B., Holm S., Holmgren D. J., Huang C.-h., Hull C., Husby D., Ichikawa S.-I., Ichikawa T., Ivezi'c Ž., Kent S., Kim R. S. J., Kinney E., Klaene M., Kleinman A. N., Kleinman S., Knapp G. R., Korienek J., Kron R. G., Kunszt P. Z., Lamb D. Q., Lee B., Leger R. F., Limmongkol S., Lindenmeyer C., Long D. C., Loomis C., Loveday J., Lucinio R., Lupton R. H., MacKinnon B., Mannery E. J., Mantsch P. M., Margon B., McGehee P., McKay T. A., Meiksin A., Merelli A., Monet D. G., Munn J. A., Narayanan V. K., Nash T., Neilsen E., Neswold R., Newberg H. J., Nichol R. C., Nicinski T., Nonino M., Okada N., Okamura S., Ostriker J. P., Owen R., Pauls A. G., Peoples J., Peterson R. L., Petravick D., Pier J. R., Pope A., Pordes R., Prosapio A., Rechenmacher R., Quinn T. R., Richards G. T., Richmond M. W., Rivetta C. H., Rockosi C. M., Ruthmansdorfer K., Sandford D., Schlegel D. J., Schneider D. P., Sekiguchi M., Sergey G., Shimasaku K., Siegmund W. A., Smee S., Smith J. A., Snedden S., Stone R., Stoughton C., Strauss M. A., Stubbs C., SubbaRao M., Szalay A. S., Szapudi I., Szokoly G. P., Thakar A. R., Tremonti C., Tucker D. L., Uomoto A., Vanden Berk D., Vogeley M. S., Waddell P., Wang S.-i., Watanabe M., Weinberg D. H., Yanny B., Yasuda N., SDSS Collaboration, 2000, AJ, 120, 1579  \nZu Y., Kochanek C. S., Kozłowski S., Peterson B. M., 2013a, ArXiv e-prints  \nZu Y., Kochanek C. S., Kozłowski S., Udalski A., 2013b, ApJ, 765, 106  \nZu Y., Kochanek C. S., Peterson B. M., 2011, ApJ, 735, 80"}