Abstract
ASASSN-14ko is a recently discovered periodically flaring transient at the center of the active galactic nucleus (AGN) ESO 253−G003 with a slowly decreasing period. Here, we show that the flares originate from the northern, brighter nucleus in this dual-AGN, post-merger system. The light curves for the two flares that occurred in 2020 May and September are nearly identical over all wavelengths. For both events, Swift observations showed that the UV and optical wavelengths brightened in unison. The effective temperature of the UV/optical emission rises and falls with the increase and subsequent decline in the luminosity. The X-ray flux, by contrast, first rapidly drops over ∼2.6 days, rises for ∼5.8 days, drops again over ∼4.3 days, and then recovers. The X-ray spectral evolution of the two flares differ, however. During the 2020 May peak the spectrum softened with increases in the X-ray luminosity, while we observed the reverse for the 2020 September peak. We found a small change in the period derivative, which seems to indicate that the system does not have a static period derivative and there is some stochasticity in its evolution.
1. Introduction
Active galactic nuclei (AGNs) normally vary stochastically across the electromagnetic spectrum, and the variability can be reasonably well modeled as a damped random walk (e.g., Kelly et al. 2008; Kozłowski et al. 2010; MacLeod et al. 2010; Zu et al. 2013). Occasionally, AGNs can undergo flares, or outbursts, in which their brightness increases dramatically over a finite period of time. These intense brightening events have been attributed to changes in inflow or accretion disk instabilities surrounding the central supermassive black hole (e.g., Kawaguchi et al. 1998), or to tidal disruption events (TDEs), in which a star is ripped apart as it passes within the tidal radius of its host SMBH (Hills 1975; Rees 1988; Evans & Kochanek 1989; Phinney 1989). These nuclear transients show diverse behavior and potentially new classes are still being discovered (e.g., Trakhtenbrot et al. 2019a; Frederick et al. 2021; Neustadt et al. 2020; van Velzen et al. 2021).
In Payne et al. (2021) we reported the discovery of ASASSN-14ko, which has shown flares with a mean period of P0 = 114.2 ±0.4 days and a period derivative of (Payne et al. 2021). These periodic flares have been observed since 2014 by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014; Kochanek et al. 2017). No prior available data appear to constrain when these flares first began, so it is possible that they have been occurring undetected for decades, or longer.
MUSE data presented in Tucker et al. (2021) revealed that the host, ESO 253−G003, is a complex merger remnant involving two AGNs and a larger tidal arm. Both nuclei are classified as an AGN by several diagnostics, but the brighter northeastern nucleus exhibits asymmetric broad-line emission while the fainter southwestern nucleus shows narrow-line emission.
As discussed in Payne et al. (2021), a likely explanation of the available data is that ASASSN-14ko is a repeating partial TDE as opposed to an SMBH binary or SMBH+perturbing star binary. Suková et al. (2021) also found some similarities between ASASSN-14ko's flares and the accretion rate simulations they found in GRMHD simulations of an orbiting star interacting with an accretion disk. Many aspects of the data, such as the similarities in the optical flares since 2014 and the modest mass-loss rates, seem most consistent with the repeating partial TDE explanation, although interpreting each flare as one pericenter passage means each flare cannot be truly identical because precession would slowly change the orbital geometry, and over tens of encounters the total mass loss ceases to be modest. Making ASASSN-14ko an SMBH binary system seems unlikely due to the lack of emission line velocity shifts, but this could also be explained if the system is face-on or it is a result of poor temporal sampling.
Payne et al. (2021) analyzed data from the six years of the ASAS-SN survey and multiwavelength data from the 2020 May flare. In this paper, we present analysis of the subsequent flare, which occurred as predicted in 2020 September. This new flare presented the opportunity to compare the X-ray, UV, and optical emission between individual flares for the first time, and to further examine how the flares fit into the model as a repeating partial TDE. In Section 2, we discuss the photometric and spectroscopic data used in this paper. In Section 3, we show that the northern, brighter nucleus in ESO 253−G003 is the origin of the flares. In Section 4, we analyze the 2020 September light curves and spectra, and we compare them to the previous flare in 2020 May. Finally, in Sections 5 and 6, we discuss the two flares in the context of TDEs and summarize the results of this paper. Throughout this paper, we assume a flat Ωm = 0.3 universe, leading to a luminosity distance of ≈188 Mpc and a projected scale of . We also assume a Galactic extinction of AV = 0.118 mag (Schlafly & Finkbeiner 2011).
2. Observations
Using our mean period P0 = 114.2 ± 0.4 days and period derivative from Payne et al. (2021), we predicted the subsequent flare would peak in the g band on UT 2020-09-7.4 ± 1.1, and the optical peak time measured in Section 4.1 was consistent with the prediction. Anticipating that the multiwavelength evolution would follow a similar trend as the 2020 May event, we scheduled X-ray, UV, and optical observations around the predicted optical peak.
2.1. ASAS-SN Photometry
ASAS-SN is an ongoing all-sky survey to discover supernovae and other transient phenomena. The 20 robotic telescopes at five sites in both the Northern and Southern Hemispheres are hosted by the Las Cumbres Observatory Global Telescope (LCOGT; Brown et al. 2013). Each of the telescopes consists of four 14 cm aperture Nikon telephoto lenses with 80 pixels and a 45 × 45 field of view. All telescopes obtain data in the g band.
The data were reduced using a fully automated pipeline based on the ISIS image subtraction package (Alard & Lupton 1998; Alard 2000). Under most circumstances, a single photometric epoch represents three combined dithered 90 s image exposures with a 45 × 45 field of view that are then subtracted from a reference image. We then used the IRAF package apphot to perform aperture photometry with a 2 pixel, or approximately 160, radius aperture on each subtracted image, resulting in a differential light curve. The AAVSO Photometric All-Sky Survey (APASS; Henden et al. 2015) was used to calibrate the photometry. All low-quality ASAS-SN images were inspected by eye to remove data affected by clouds or other systematic problems.
2.2. Swift XRT and UVOT Photometry
Using our prediction of the optical peak time and the expected multiwavelength trend based on the 2020 May flare, we requested Swift UltraViolet/Optical Telescope (UVOT; Roming et al. 2005) target-of-opportunity (ToO) observations (ToO ID: 14389, 14488, PI: Payne). The UVOT data were obtained in six filters (Poole et al. 2008): V (5425 Å), B (4350 Å), U (3467 Å), UVW1 (2581 Å), UVM2 (2246 Å), and UVW2 (2055 Å). We used the HEAsoft ( HEASARC 2014) software task uvotsource to extract the source counts using a 160 radius aperture and used a sky region of ∼400 radius to estimate and subtract the sky background. This aperture size was used to match the ASAS-SN photometry. All fluxes were aperture corrected and converted into magnitudes and fluxes using the most recent UVOT calibration (Poole et al. 2008; Breeveld et al. 2010), and were corrected for Galactic extinction. To properly isolate and measure the transient flux in each epoch, the quiescent host fluxes were measured in the same aperture and subtracted. We converted the Swift UVOT B and V magnitudes to Johnson B and V magnitudes using the standard color corrections. 16
2.3. Las Cumbres Observatory Global Telescope Photometry and Spectroscopy
We obtained photometric observations from the LCOGT (Brown et al. 2013) using the 1 m telescope at Siding Spring Observatory in New South Wales, Australia. The data consisted of B, V, , and observations. Aperture magnitudes were obtained using a 160 radius aperture using the IRAF apphot package using an annulus to estimate and subtract background counts. The data were calibrated using stars with APASS DR 10 magnitudes, and the aperture magnitudes were corrected for Galactic extinction. The quiescent host fluxes were measured in the same aperture and subtracted.
2.4. Pan-STARRS Photometry
Located atop Haleakalā on Maui, the Pan-STARRS1 telescope has a 1.8 m diameter primary mirror with a wide-field 1.4 gigapixel camera consisting of sixty Orthogonal Transfer Array devices, which has carried out a set of synoptic imaging sky surveys with resulting catalog data products (Chambers et al. 2016; Flewelling et al. 2020). Pan-STARRS1 data are obtained using grizyP1 filters, similar to the Sloan Digital Sky Survey filters (Abazajian et al. 2009). The Pan-STARRS1 photometric system is discussed in detail in Tonry et al. (2012).
The Image Processing Pipeline (IPP; see details in Magnier et al. 2020) processes the Pan-STARRS1 data. The processing steps include device "de-trending," a flux-conserving warping to a sky-based image plane, masking and artifact location involving bias and dark correction, flat-fielding, and illumination correction through rastering sources across the field of view (Waters et al. 2020). There were 28 epochs of ASASSN-14ko Pan-STARRS1 data between 2014 and 2020. These data covered times both during quiescence and during separate flares that peaked in the optical on MJD , , and (Payne et al. 2021). We used the difference imaging technique created by the IPP to obtain stacked Pan-STARRS1 3π data reference images using data taken in quiescence. Those reference images were then subtracted from the images taken during the prior outbursts to isolate the transient (e.g., Huber et al. 2015).
2.5. SALT Spectroscopy
We used the 10 m Southern African Large Telescope (SALT; Buckley et al. 2006) with the Robert Stobie Spectrograph (RSS; Burgh et al. 2003; Kobulnicky et al. 2003) to obtain optical spectra during the flare. The data were obtained with the 15 slit on UT 2020-09-04 and UT 2020-09-16, which corresponded to ∼1 day prior and ∼11 days after the measured optical g-band peak, respectively. The slit position was oriented with a position angle of −1° for the first spectrum and 30° for the second spectrum. The data were then reduced using standard procedures, including bias subtraction, flat-fielding, wavelength calibration, sky subtraction, and flux calibration.
2.6. Amateur Astronomer Photometry
Data were collected using a 41 cm telescope at Savannah Skies Observatory from Queensland, Australia. The 180 s exposures were obtained using the B, V, and RC filters, and the data were bias and dark subtracted. Data were also taken at Moondyne Observatory located east of Perth, Australia, using a 0.4 m telescope with AOX adaptive optics. Guided BC -, V-, GS -, and RS -band images were reduced and calibrated with bias and dark subtraction and flat-field normalization. Data were taken twice nightly over the duration of the event with 120 and 600 s exposures, and then aligned using background stars.
3. The Brighter Nucleus Is the Origin of the Flares
Using the spatially resolved MUSE IFU spectra, Tucker et al. (2021) showed that the nuclear region of ESO 253−G003 houses two AGNs separated by 17 as part of a post-merger system. While ASAS-SN's astrometry is generally accurate to 10, it is not ideal for conclusively determining which nucleus is the origin of the outbursts.
Figure 1 shows the Pan-STARRS1 reference image of the nuclear region of ESO 253−G003. The ASAS-SN position for ASASSN-14ko (Holoien et al. 2014a) is marked and consistent with the northern nucleus. The Pan-STARRS1 data also show that the location of the flares corresponds to the northern nucleus. The average and standard deviation of the Pan-STARRS1 difference image source positions is shown by the ellipse in Figure 1. Its coordinates in R.A. and decl. are 813255 ± 00002 and −460056 ± 00002.
4. X-Ray, UV, and Optical Evolution
The 2020 September flare was the first opportunity to directly compare the X-ray and UV evolution between events. In addition, this was the first outburst where the UV rise was observed. These observations allow us to characterize the high-energy properties of ASASSN-14ko's flares in a wider context, instead of as a singular episode.
4.1. UV/Optical Evolution
Figure 2 shows the host-subtracted data from 2020 September and also given in Table 1. The host magnitudes were determined from quiescent 2020 May data derived in Payne et al. (2021). We used the September flare's peak timing to update the parameters in the timing model. Following the same procedure as in Payne et al. (2021), we fit the ASAS-SN g-band light curve with a fifth-order polynomial and found the peak time of MJD . The errors on the peak time were measured by bootstrap resampling the light curve. Combining this time with the previous seventeen reported in Payne et al. (2021) and fitting them using the peak timing equation from Payne et al. (2021), the time of the flare peak n is
where t0 = 2,456,854.9 ± 2.6 is the reference time, P0 = 114.4 ±0.4 days is the mean period, and is the period derivative. The original detection is defined as n = 0 and the current flare is n = 20. The errors on the peak times were expanded in quadrature to get a χ2 per degree of freedom of unity. These updated parameters predict the next flare to peak in the g band on MJD 59,208.76 ± 1.1.
Download figure:
Standard image High-resolution imageTable 1. Light-curve Data of ASASSN-14ko Used in This Analysis
JD | Band | X-Ray Energy | LX or λ Lλ |
---|---|---|---|
Central λ | (1043 erg s−1) | ||
2,459,089.823 | X-ray | 0.3–10.10 keV | 0.174 ± 0.057 |
2,459,090.746 | W2 | 2055 Å | 2.196 ± 0.540 |
2,459,090.752 | M2 | 2246 Å | 1.375 ± 0.413 |
2,459,089.823 | W1 | 2581 Å | 0.584 ± 0.362 |
2,459,089.824 | U | 3467 Å | 0.374 ± 0.308 |
2,459,091.806 | B | 4350 Å | 1.428 ± 0.501 |
2,458,982.732 | g | 4640 Å | 1.029 ± 0.128 |
2,459,087.100 | V | 5425 Å | 1.121 ± 0.440 |
2,459,087.771 | R/r | 6166 Å | 0.942 ± 0.181 |
Note. Only one observation in each band is shown here to demonstrate its form and content. R/r is the combined , RC , and RS data from the different observatories. Table to be published in its entirety in machine-readable form in the online journal.
A machine-readable version of the table is available.
Download table as: DataTypeset image
Figure 2 also shows the photometry from the 2020 May flare, where we have aligned them using the predicted peak times from the updated timing model. The similarities of the flares across the UV/optical is remarkable. We then fit each epoch using a blackbody model to track the change in blackbody luminosity, temperature, and radius for the duration of the flare. As shown in Figure 3, the blackbody temperature and luminosity evolved rapidly during the rise, peaking within two days. Considering this rapid early evolution, subdaily or a 6 hr cadence will be a crucial requirement for future observing campaigns to provide better constraints on any potential short-timescale changes and to look for differences between the peak times of different wavelengths.
Download figure:
Standard image High-resolution imageFigure 4 compares the SALT spectra and the MUSE spectrum obtained during quiescence in 2015 that was analyzed in Tucker et al. (2021). We extracted a spectrum from the MUSE data cube using a rectangular region of width 15 matching the geometry of the SALT slit on UT 2020-09-04. The blue continuum is similar but shallower than what was observed for the original discovery spectrum taken seven days post peak (Holoien et al. 2014a). There was also a small degree of broadening in Hβ during the flare. However, we caution that the nuclear region is especially complex (Tucker et al. 2021) and the extracted spectra depend on seeing, slit width, slit angle, and the extraction region as illustrated by the right panel of Figure 4. Future work to robustly measure changes in the line profiles during an outburst will require integral field unit spectroscopy or careful consideration of the slit placement.
Download figure:
Standard image High-resolution image4.2. X-Ray Light-curve Evolution
As with the UV/optical emission, the X-ray flux evolved in a consistent manner between the 2020 May and September flares, as shown in Figure 2. However, the X-ray evolution contrasts sharply with the UV/optical behavior. The X-ray flux first decreased during the UV/optical rise, and then increased as the UV/optical began to decline. The initial drop in the X-ray flux was very rapid, by a factor of ∼9 in ∼2.6 days. The second dip in this "double-dip" behavior was similarly rapid, with a decrease in the flux by a factor of ∼4 in ∼4.3 days. Comparing the 2020 May and September events indicates that the drop in X-ray flux during the UV/optical rise is a characteristic feature of both flares.
Similarly to the evolution of the total X-ray flux between the two flares, the hardness ratios evolve similarly, as shown in Figure 5. In 2020 May, the peak X-ray luminosities corresponded to the lowest hardness ratios implying relatively softer X-ray emission when the source is X-ray brightest, and this trend also occurred in 2020 September. In addition, the X-ray flux minimum during the UV/optical peak coincided with an increased hardness ratio in 2020 May and 2020 September. The hardness ratio at the optical peak in 2020 May was −0.24 ± 0.08, whereas the hardness ratio at the optical peak in 2020 September was 0.14 ± 0.04, as shown by the larger points in Figure 5. The correlation is, however, rather noisy and more observations are needed to fully confirm it.
Download figure:
Standard image High-resolution image4.3. X-Ray Spectra
Given the numbers of counts, we can only study the average X-ray spectrum in detail. Figure 6 compares the stacked Swift spectra of the 2020 May and September flares. The best fit to the 2020 September X-ray data was the power law plus blackbody xspec model:
with a photon index of 0.98 ± 0.20, a temperature of 0.13 ± 0.03 keV, and a reduced χ2 of 1.06. For fitting the spectra, the absorbing column density was set to the Galactic column density along the line of sight, 3.5 × 1020 cm−2, because allowing it to be a free parameter did not improve the fits. However, the data were equally well fit by two power laws with photon indices of 2.17 ± 0.12 and 0.10 ± 0.15 with a reduced χ2 of 1.20. Assuming a single power law was not a very good fit, with a reduced χ2 of 1.70. All of the X-ray spectral model parameters are summarized in Table 2.
Download figure:
Standard image High-resolution imageTable 2. Best-fit Parameters for the 2020 September Swift X-Ray Spectra Compared with the 2020 May Swift Spectra and XMM-Newton Spectra Reported in Payne et al. (2021)
Instrument | kT (keV) | Blackbody Normalization (K) | Photon Index Γ | χ2 per dof | Absorbed Flux [0.3–10.0 keV] (erg cm− 2 s−1) |
---|---|---|---|---|---|
Swift 2020 Sep | 0.13 ± 0.03 | 0.98 ± 0.20 | 1.06 | (1.10 ± 0.13) × 10−12 | |
Swift 2020 May | 0.13 ± 0.03 | 1.07 | |||
XMM+NuSTAR | 0.15 ± 0.01 | 0.87 ± 0.04 | 1.64 |
Note. All models used a redshift of 0.0425 and a neutral hydrogen column density NH of 3.5 × 1020 cm−2 fixed to the Galactic column density (HI4PI Collaboration et al. 2016). The power-law normalization for the Swift 2020 September spectrum is (6.70 ± 1.40) × 10−5. The Swift 2020 September spectrum fit had 76 degrees of freedom, the Swift 2020 May spectrum fit had 56 degrees of freedom, and the XMM-Newton+NuSTAR spectrum fit had 470 degrees of freedom.
Download table as: ASCIITypeset image
As reported in Payne et al. (2021), the best-fit model to the merged Swift data from 2020 May was a single power law with photon index and the fit was not significantly improved by adding an additional blackbody with temperature 0.13 ± 0.03 keV. However, the data from 2020 May also had fewer total counts due to the lower total integration time (22.5 ks versus 56.0 ks). Overall, even though the blackbody plus power law and the two power laws fit the spectrum well, a power law with a photon index of 0.10 plus blackbody is more aligned with the expectation of AGNs that have blackbody temperatures of >0.1 keV.
Archival XMM-Newton data of ESO 253−G003 taken during ASASSN-14ko's quiescence in 2015 show a prominent Fe Kα line near 6.4 keV, as shown in Figure 6. There is also a small increase in flux at this energy in the Swift XRT spectra; however, there are too few counts in the Swift observations to properly constrain the presence of this line. Due to the spatial resolution of Swift and XMM-Newton, we are unable to resolve the two distinct AGNs in ESO 253−G003, making it difficult to disentangle the X-ray emission between the two sources. Observing changes in the Fe Kα emission with XMM-Newton would indicate that the source is compact rather than extended, but higher resolution Chandra images are needed to differentiate between the two nuclei.
We also separately considered the average X-ray spectra before and after the UV/optical peak. Before the peak, we found that an absorbed power law plus blackbody model is a better fit than two power laws or a single power-law model. This model consists of a photon index of 0.91 ± 0.1 and a temperature of 0.16 ± 0.03 keV with a reduced χ2 of 1.5. The best post-UV/optical peak model is also an absorbed power law plus blackbody but with a photon index of 0.62 ± 0.09 and a temperature of 0.18 ± 0.02 keV with a reduced χ2 of 1.1. Comparing the spectra between these two epochs shows that after UV/optical peak, the photon index decreased and became harder. This is consistent with the hardness ratio evolution shown in Figure 2.
5. Discussion
The outbursts of ASASSN-14ko continue to follow a periodic trend. However, the updated timing model parameters t0, P0, and changed compared to those found in Payne et al. (2021), although the g-band optical peak occurred within the error bars on the prediction in Payne et al. (2021). Compared to 2020 May, changed by −0.0003. Seeing a change in between events seems to indicate that the system does not have a global and the orbit of a partial TDE should show some stochasticity in its evolution. Regardless of this variation, the timing model still describes the system well and was able to correctly predict the next flare in 2020 December (A. V. Payne et al. 2022, in preparation).
We now see in two flares that the rapid rise and fall observed in the optical over the last 6 yr is also typical of the UV emission. This ∼20 day timescale over which the flares' blackbody luminosity decays is much shorter than for all previously studied TDEs. In the Hinkle et al. (2021, 2020) studies of the correlation between peak luminosity and decline rates in TDEs, all of the examples decline on timescales of over 100 days.
While the optical, UV, and X-ray luminosity evolution of the individual flares are nearly identical, the X-ray spectra of the May and September flares evolved differently. In May, the spectra softened as the X-ray luminosity increased, while in September they hardened. This is the first property found that distinguishes individual flares. Auchettl et al. (2017) systematically studied a sample of TDEs in order to uniformly characterize the X-ray emission of TDEs as a whole. Auchettl et al. (2017) found that all TDE candidates had soft X-ray emission with HRs ranging from −1.0 to +0.3 and the majority of objects had HRs between −1.0 and 0 while highly variable AGNs had hardness ratios ranging from −1 to +1. Thus, both flares are consistent with previously detected TDE candidates even though the HR evolutions are distinct.
ASASSN-14ko's X-ray behavior is similar to another nuclear transient discovered by the ASAS-SN survey, ASASSN-18el (Nicholls et al. 2018). ASASSN-18el displayed strong Balmer emission lines that broadened after the galaxy's nucleus brightened and was subsequently described as a changing-look AGN (Trakhtenbrot et al. 2019b). Further X-ray monitoring by Ricci et al. (2020) revealed strong variability on longer timescales. The X-ray luminosity decreased rapidly ∼160 days after the optical event and then slowly increased to eventually exceed its pre-flare level. This led to an X-ray light curve characterized by an asymmetric trough sharply contrasting with the optical brightening that triggered the event. While the overall pattern is similar, the timing of the X-ray changes relative to the UV/optical is different in ASASSN-14ko. ASASSN-18el's X-ray trough occurred much later after the UV/optical brightening, whereas ASASSN-14ko's X-ray flare decreased during the UV/optical rise and peak, with a secondary dip later in the UV/optical decline. In addition, ASASSN-18el occurred over hundreds of days as opposed to ASASSN-14ko's much faster evolution where the X-ray flux decrease occurred over only ∼14 days.
Ricci et al. (2020) also found that the X-ray flux variability of ASASSN-18el and spectral variability were tightly connected, becoming harder at times of increased X-ray flux. ASASSN-14ko did not show this trend, and instead showed the opposite. ASASSN-18el's spectrum was dominated by a blackbody-like component after the optical brightening, while the previously observed power-law component faded during the decreased X-ray flux (Ricci et al. 2020). When the X-ray flux increased again, the power-law component returned. Ricci et al. (2020) argued that the disappearance and reappearance of the power-law component was evidence that the X-ray corona was destroyed by a TDE and then gradually reformed. This is another facet that contrasts with ASASSN-14ko, whose X-ray spectrum was consistently modeled as a power law plus blackbody for the duration of the 2020 September event. That ASASSN-14ko's power-law component remains present throughout the event may suggest that the innermost regions of the accretion flow remains intact as opposed to being destroyed by the TDE.
Changes in spectral lines may be clear discriminators between TDE or AGN activity, or key features defining new classes of objects (e.g., Trakhtenbrot et al. 2019a; Frederick et al. 2021; Neustadt et al. 2020; van Velzen et al. 2021). Reverberation mapping (Blandford & McKee 1982; Peterson 1993; Peterson et al. 2004) is a technique that has been used for decades as a powerful method to probe the innermost regions of AGNs. The timings of ASASSN-14ko's future flares can be predicted, allowing high-cadence spectroscopic campaigns to be organized around future events that would be able to monitor changes in the spectral lines over the flare and how they correlate with changes in the continuum. AGN broad lines typically become narrower as they become more luminous due to the lower Keplerian velocity at greater distances. TDEs, however, show the opposite behavior. Previously studied TDEs have shown that optical emission lines strengthen and broaden with increased luminosity. For example, both a blue continuum and prominent, broad emission lines He ii λ4686, Hβ, and Hα formed in the TDEs ASASSN-14ae (Holoien et al. 2014b) and ASASSN-14li (Jose et al. 2014; Holoien et al. 2016a). For ASASSN-14ae they narrowed and then faded completely after ∼750 days (Brown et al. 2016), while they were still prominent at ∼500 days for ASASSN-14li (Brown et al. 2018). This behavior was also apparent in "helium-rich" TDE ASASSN-15oi (Holoien et al. 2016b, 2018), and Bowen TDE ASASSN-18pg (Leloudas et al. 2019; Holoien et al. 2020).
Future spectroscopic campaigns need to carefully consider the complex host galaxy background at the center of ESO 253−G003. A spectroscopic observing campaign using long-slit spectrographs would require careful consideration of the slit orientations to either avoid contamination from the southern nucleus or try to always have the same level of contamination. Ideally, the observations would use an IFU spectrograph where the extracted aperture and background can be selected.
Overall, the new data presented here support the partial TDE hypothesis. In Payne et al. (2021), an SMBH binary system and an SMBH+perturbing star were discussed as possible alternate explanations. The new data presented here do not provide evidence to support these other scenarios. There was a change in , but this updated value is still over an order of magnitude larger than predicted for an SMBH binary system with ASASSN-14ko's black hole mass. In the case of a star in orbit around a single SMBH that perturbs its accretion disk, the star alternates between downward plunges into the disk and upward escapes out of the disk as viewed by the observer. This difference ought to imprint differences in the flares between even and odd events, such as changes in their amplitudes or profiles. However, there were no discernible differences detected in the long-term ASAS-SN g-band data as reported in Payne et al. (2021), and here, the X-ray/UV light curves between the two individual flares are also indistinguishable. Continued observations and monitoring are crucial to conclusively determine if there are any long-term differences or trends in the X-ray/UV flare structure.
6. Conclusions
ASASSN-14ko continues to be a periodic nuclear transient at the center of AGN ESO 253−G003 whose flares can be predicted with a timing model consisting of a mean period and a nonzero, negative period derivative. The 2020 September flare event was the first opportunity to directly compare individual outbursts at X-ray and UV wavelengths. In this work, we show that:
- 1.The brighter, northeastern nucleus of ESO 253−G003 is the source of the flares.
- 2.ASASSN-14ko's updated timing model parameters are t0 = 2,456,854.9 ± 2.6, P0 = 114.4 ± 0.4 days, and after including the latest peak timing from 2020 September, which occurred on MJD . The model predicted a next peak for 2020 December on MJD 59,208.76 ± 1.1 that occurred on MJD 59,205.65 ± 0.19 (A.V. Payne et al., 2022, in preparation).
- 3.When comparing the 2020 May and September X-ray and UV light curves, the two flares show a rapid UV/optical brightening coincident with an equally rapid decrease in X-ray flux.
- 4.The X-ray hardness ratio evolution was consistent between the two flares. The flares are softer at peak X-ray luminosity.
Future flares will continue to be important to determine which characteristics change and which remain the same, especially in the high-energy regime. One aspect in particular is determining if there will be any continued difference in the spectral properties. In addition, future reverberation mapping campaigns centered around the predicted flare timings could provide an accurate black hole mass measurement, especially when using IFU instruments to avoid contamination from the southwestern nucleus of the host galaxy.
We thank the anonymous referee for their constructive comments. We thank Jennifer van Saders for helpful discussion. We thank the Las Cumbres Observatory and its staff for its continuing support of the ASAS-SN project. ASAS-SN is supported by the Gordon and Betty Moore Foundation through grant GBMF5490 to the Ohio State University, and NSF grants AST-1515927 and AST-1908570. Development of ASAS-SN has been supported by NSF grant AST-0908816, the Mount Cuba Astronomical Foundation, the Center for Cosmology and AstroParticle Physics at the Ohio State University, the Chinese Academy of Sciences South America Center for Astronomy (CASSACA), and the Villum Foundation.
A.V.P. acknowledges support from the NASA Graduate Fellowship through grant 80NSSC19K1679. B.J.S., C.S.K., and K.Z.S. are supported by NSF grant AST-1907570. B.J.S. is also supported by NASA grant 80NSSC19K1717 and NSF grants AST-1920392 and AST-1911074. C.S.K. and K.Z.S. are supported by NSF grant AST-181440. J.T.H. is supported by NASA award 80NSSC21K0136. Support for T.W.-S.H. was provided by NASA through the NASA Hubble Fellowship grant HST-HF2-51458.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. Parts of this research were supported by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. T.A.T. is supported in part by Scialog Scholar grant 24215 from the Research Corporation. M.A.T. acknowledges support from the DOE CSGF through grant DE-SC0019323. The SALT observations presented here were made through Rutgers University program 2020-1-MLT-007 and supported by NSF grant AST-1615455 to S.W.J.
Operation of the Pan-STARRS telescopes is primarily supported by the National Aeronautics and Space Administration under grant No. NNX12AR65G and grant No. NNX14AM74G issued through the SSO Near-Earth Object Observations Program. The PanSTARRS1 Surveys (PS1) and the PS1 public science archive have been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, the Queens University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation grant No. AST1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory, and the Gordon and Betty Moore Foundation.
Software: astropy (Astropy Collaboration et al. 2018), ftools (Blackburn 1995), HEAsoft (HEASARC 2014), IRAF (Tody 1986, 1993), numpy (Harris et al. 2020), matplotlib (Hunter 2007).