Issue |
A&A
Volume 638, June 2020
|
|
---|---|---|
Article Number | A13 | |
Number of page(s) | 7 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/202037532 | |
Published online | 29 May 2020 |
A flare in the optical spotted in the changing-look Seyfert NGC 3516⋆
1
Department of Astronomy, Faculty of Mathematics, University of Belgrade, Studentski trg 16, 11000 Belgrade, Serbia
e-mail: dilic@matf.bg.ac.rs
2
Sternberg Astronomical Institute, Moscow M.V. Lomonosov State University, Universitetskij Pr., 13, Moscow 119234, Russia
3
Astronomical Observatory Belgrade, Volgina 7, 11000 Belgrade, Serbia
4
Department of Physics and Astronomy, University of Turku, 20014 Turku, Finland
5
Space Research Institute of the Russian Academy of Sciences, Profsoyuznaya Str. 84/32, Moscow 117997, Russia
6
Faculty of Natural Sciences and Mathematics, University of Banjaluka, Mladena Stojanovića 2, 78000 Banjaluka, Republic of Srpska, Bosnia and Herzegovina
7
Faculty of Physics, Moscow M.V. Lomonosov State University, Leninskie Gory 1, Moscow 119991, Russia
Received:
20
January
2020
Accepted:
2
April
2020
Context. We present observations from the short-term intensive optical campaign (from September 2019 to January 2020) of the changing-look Seyfert NGC 3516. This active galactic nucleus is known to have strong optical variability and has changed its type in the past. It has been in the low-activity state in the optical since 2013, with some rebrightening from the end of 2015 to the beginning of 2016, after which it remained dormant.
Aims. We aim to study the photometric and spectral variability of NGC 3516 from the new observations in U- and B-bands and examine the profiles of the optical broad emission lines in order to demonstrate that this object may be entering a new state of activity.
Methods. NGC 3516 has been monitored intensively for the past 4 months with an automated telescope in U and B filters, enabling accurate photometry of 0.01 precision. Spectral observations were triggered when an increase in brightness was spotted. We support our analysis of past-episodes of violent variability with the UV and X-ray long-term light curves constructed from the archival Swift/UVOT and Swift/XRT data.
Results. An increase of the photometric magnitude is seen in both U and B filters to a maximum amplitude of 0.25 mag and 0.11 mag, respectively. During the flare, we observe stronger forbidden high-ionization iron lines ([Fe VII] and [Fe X]) than reported before, as well as the complex broad Hα and Hβ lines. This is especially seen in Hα, which appears to be double-peaked. It seems that a very broad component of ∼10 000 km s−1 in width in the Balmer lines is appearing. The trends in the optical, UV, and X-ray light curves are similar, with the amplitudes of variability being significantly larger in the case of UV and X-ray bands.
Conclusions. The increase of the continuum emission, the variability of the coronal lines, and the very broad component in the Balmer lines may indicate that the AGN of NGC 3516 is finally leaving the low-activity state in which it has been for the last ∼3 years.
Key words: line: profiles / galaxies: active / galaxies: Seyfert / quasars: individual: NGC 3516 / quasars: emission lines
Photometric measurements are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/638/A13
© ESO 2020
1. Introduction
The open question concerning the mechanism by which activity is triggered and evolves in active galactic nuclei (AGNs) is still pertinent, especially because we know that growth of supermassive black holes in the centers of galaxies happens during the AGN stage. We hope to get closer to answering the above question by studying the extremely variable AGNs, that is, the class of so-called changing-look AGNs (CL AGNs). These objects show extreme changes of emission line intensities, with sometimes almost complete disappearance and reappearance of the broad component (see, e.g., Lyuty et al. 1984; Kollatschny & Fricke 1985; Denney et al. 2014; Oknyansky et al. 2019a,b, etc.). Recently, the CL AGNs have come into immediate focus, with large monitoring and spectroscopic surveys (e.g., Runco et al. 2016; MacLeod et al. 2016). There are even speculations that each strongly variable AGN could almost be a CL AGN if constantly observed (see discussion in Oknyansky et al. 2017).
There are several physical processes that could cause such dramatic changes. The main ones are: changes of accretion disk structure (e.g., Stern et al. 2018), broad-line region (BLR) occultation by obscuring material (e.g., Elitzur 2012), tidal disruption events (TDEs), supernova explosions, tidal stripping of stars (see, e.g., Wang et al. 2012; Campana et al. 2015; Komossa et al. 2017), or a combination of the above. So far, it seems that the significant change of the accretion rate is responsible for the change (Noda & Done 2018; Śniegowska & Czerny 2019), whereas obscuration is less favored as a major cause (Oknyansky et al. 2017; Mathur et al. 2018). However monitoring and multiwavelength observations of CL AGNs are needed to discriminate among the different scenarios (e.g., MacLeod et al. 2019).
NGC 3516 is one of the original six Seyfert galaxies (Seyfert 1943). It is a face-on spiral galaxy at a distance of between D ≈ 38 Mpc and 66 Mpc1. NGC 3516 has previously been the subject of several optical monitoring campaigns (e.g., Lyutyi & Doroshenko 1993; Wanders et al. 1993; Wanders & Horne 1994; Maoz et al. 2002; Noda et al. 2016; De Rosa et al. 2018; Shapovalova et al. 2019, etc.), which showed that the object is variable on short and long timescales. Reverberation mapping analysis showed that the Balmer lines respond with a delay of ∼10–15 days to the continuum variations (Wanders et al. 1993; De Rosa et al. 2018; Shapovalova et al. 2019) with even shorter delays in the line wings (De Rosa et al. 2018). NGC 3516 was recently confirmed to be an optical CL AGN as a result of a long-term (∼22 years, 1996–2018) optical monitoring campaign (Shapovalova et al. 2019). The campaign showed that the broad Balmer emission lines almost completely disappeared in 2014, and that at the end of the campaign in 2018 a weak, blueshifted and asymmetric broad component started to reappear (Shapovalova et al. 2019). This AGN is known to show dramatic optical variability and has changed its type in the past several times (Andrillat & Souffrin 1968; Shapovalova et al. 2019). However, it has been in a low-activity state in the optical since 2013, with some rebrightening occurring at the end of 2015 until the beginning of 2016, after which it remained in a dormant state. For NGC 3516, it is well known that the broad emission lines show complex, multicomponent profiles, indicating complex kinematics of the BLR, with evidence of outflowing, infalling, and a disk-like emitting region that mostly contributes to the line wings (e.g., Popović et al. 2002; Denney et al. 2010; Storchi-Bergmann et al. 2017; De Rosa et al. 2018).
NGC 3516 shows even more interesting properties in the high-energy band (e.g., Edelson et al. 2000; Netzer et al. 2002; Maoz et al. 2002; Turner et al. 2005; Mehdipour et al. 2010; Liu et al. 2014; Huerta et al. 2014; Noda et al. 2016). Edelson et al. (2000) showed that there are uncorrelated trends on longer timescales between the X-ray and optical bands, which was supported with the later findings of Maoz et al. (2002). This is not unusual for Seyfert galaxies, as discussed by Edelson et al. (2000). Recently, Noda et al. (2016) found from a simultaneous X-ray and optical monitoring program during the faintest phase of NGC 3516 (2013–2014) that the X-ray flux and B-band variations were comparable and significantly correlated, with the B-band delayed by ∼2 days relative to the X-rays. Both results indicate that the standard X-ray reprocessing model in which the X-ray source heats a stratified accretion disk, which then re-emits in the optical and ultraviolet, may not be applicable in the case of NGC 3516 and that these two emitting regions are powered primarily by different processes (Edelson et al. 2000; Noda et al. 2016). Also, NGC 3516 is known for the presence of complex absorption, showing the presence of eight kinematic components in the UV (Kraemer et al. 2002), supported by at least three distinct outflowing components in the X-rays (Mehdipour et al. 2010; Turner et al. 2011; Huerta et al. 2014). Interestingly, Mehdipour et al. (2010) concluded, similarly to Netzer et al. (2002), that the X-ray variability is better understood as a consequence of changes in the source continuum emission than in the warm absorber.
Here we report on the very recent flare discovered during the latest intensive photometric monitoring in U and B filters that may be an indicator that NGC 3516 is transitioning from a low to a high state of activity. We also acquired an optical spectrum during this period, covering Hβ and Hα lines, in order to find differences between their profiles (intensities) with respect to the latest published ones (Shapovalova et al. 2019). In order to analyze the previous episodes of violent variability, our optical observations are supported with the archival UV data of Swift/Ultraviolet/Optical Telescope (UVOT) and X-ray Swift/XRT data of NGC 3516. This paper is organized as follows: in Sect. 2, we describe observations and data reduction; in Sect. 3, we list the main results; in Sect. 4, we discuss our findings; and in Sect. 5, we outline our conclusions and future work.
2. Observations and data reduction
2.1. Photometry
The intensive photometric monitoring was performed in U and B filters from 2019 Sept 16 to 2020 Jan 04 with the RC600, a new automated telescope of the Caucasus Mountain Observatory at the Sternberg Astronomical Institute (SAI CMO) of Moscow State University. The RC600 has a 600 mm main mirror with a focal length of 4200 mm and is equipped with an Andor iKon-L BV camera (installed filters U, B, V, Rc, Ic, g′, r′, i′). The field of view is 22′ × 22′, with scaling of 0.67″ per pixel (for more details see Berdnikov et al. 2020). During these 4 months, data were obtained during good photometric nights with several exposures (2−6) in each band, with typical exposure times of 300 s in U and 20−30 s in B filter. The photometry was done with a small aperture of 6.7″ in diameter, and using comparison stars from Lyutyi & Doroshenko (1993).
To compare the new data with the long-term trend in NGC 3516, we use: (i) data from the UBV photoelectric monitoring initiated in the 1970s (see details in Lyutyi 1972), with the aperture of 14.3″ in diameter using a photoelectrical photometer with EMI 9789 photomultiplier attached to a a 60 cm telescope in the Crimean astronomical station SAI MSU (the same as was used by Lyutyi & Doroshenko 1993), and the same comparison stars as in this campaign, and (ii) data from the long-term campaign with the aperture size of 10″, which are described in Shapovalova et al. (2019). To match these two long-term light curves, we used the overlapping points to reduce all data to the same aperture of 10″. The long-term photometric light curve in B filter is presented in Fig. 1.
![]() |
Fig. 1. Multiwavelength light curves of NGC 3516. Top panel: long-term photometric optical and UV light curves in B and UVW2 filters: circles denote data from Shapovalova et al. (2019), squares denote the result from this paper (SAI CMO), diamonds stand for the data from SAI Crimea campaign shown for continuity, and asterisks denote UVW2 UVOT Swift photometry. All magnitudes are reduced to the same aperture of 10″. A dashed horizontal line at magnitude 14.1 is given at the bottom of the figure to guide the eye, and the arrow indicates epochs when spectra were acquired. Bottom panel: Swift/XRT 0.5−10 keV X-ray flux light curve in erg cm−2 s−1. |
2.2. Spectroscopy
The new spectra were obtained on 2019 Dec 07 (JD 2458824.57) and 2019 Dec 09 (JD 2458826.59), with the 2.5 m telescope of the SAI CMO equipped with the medium-resolution optical double-beam spectrograph, the Andor Newton 940P cameras with CCDs E2V CCD42-10, and the volume phase holographic grating. The slit width was 1″ and 1.5″ for 2019 Dec 07 and 2019 Dec 09, with exposure times of 600 s and 1000 s, respectively, giving a spectral resolution of ∼4 Å and signal-to-noise ratio (S/N) of 30−50 near Hβ, and ∼3 Å and S/N of 40−90 near Hα. The spectrophotometric data reduction was carried out using our python3 scripts, including bias, flat-field, and dark corrections, cosmic ray removal, two-dimensional wavelength linearization, background subtraction, and relative flux calibration based on spectrophotometric standard stars.
2.3. Archival UV and X-ray data
The Neil Gehrels Swift Observatory (Gehrels et al. 2004) provided simultaneous XRT and UVOT monitoring of NGC 3516 from 2006 to 2018, with the most intensive campaign taking place in 2012–2014. Some of the Swift data used in our study were previously published by Buisson et al. (2017). In order to trace the behavior of NGC 3516 on a longer timescale, data from both the XRT and UVOT telescopes were re-reduced consistently from the raw data available in the archive (for details see Buisson et al. 2017; Oknyansky et al. 2017). The X-ray spectra were obtained through the standard online tools provided by the UK Swift Science Data Centre2 (Evans et al. 2009). To account for low count statistics, the XRT spectra in the 0.5−10 keV range were binned to assure at least one count per energy bin and fitted using the W-statistic (Wachter et al. 1979). To get the source flux in physical units, spectra were fitted with a simple absorbed power-law model leaving the photon index free to vary and freezing the equivalent absorption column NH at 3 × 1020 cm−2, the Galactic column density along the line of sight to the galaxy (Kalberla et al. 2005).
The image analysis of the Swift/UVOT data in different bands (V, B, U, UVW1, UVW2, UVM2) was done again following the procedure provided by the UK Swift Science Data Centre. Photometry was performed with the UVOTSOURCE tool with aperture diameter of 10″ and 40″ for the source and background, respectively. The archive provides the largest number of epochs for the UVW2 magnitude. The background was chosen with the center about 2′ away from the galaxy for all filters. The resulting light curves of the XRT X-ray flux in the 0.5−10 keV band and the UVOT UVW2 magnitude are shown in Fig. 1.
3. Analysis and results
From the long-term light curve in Fig. 1, it is clear that NGC 3516 experienced several low and high-activity states. The object was the brightest in 2007, followed with the minimum in 2014 (Shapovalova et al. 2019). Several weaker changes in the optical flux happened afterwards, with the strongest subsequent flare occurring at the end of 2015 when the broad component of Hα line started to increase (see Fig. 12 in Shapovalova et al. 2019), after which the object entered the low-activity state. The latest brightening probably started from the end of 2018, followed by the local minimum at the end of July 2019, after which there was again brightening with the maximum in January 2020. This is more obvious when we plot only the latest data, now both in U and B filters (Fig. 2). The accuracy of photometric measurements extracted from the small-size aperture allows for the short-term variability to be seen. The increase in the U magnitude is more prominent, and is 0.25 mag, whereas in the B filter, it is 0.11 mag (Fig. 2). Figure 2 shows the trend of increase, and compared to the minimum seen in the first half of 2019 (Fig. 1) the change in the B magnitude (scaled to 10″ aperture) is ∼0.15. The variability amplitude appears lower in the larger-size aperture, therefore the gradient in the smaller-size aperture is probably larger.
![]() |
Fig. 2. Short-term light curves in U and B bands (aperture size 6.7″) from the high-cadence monitoring campaign from 2019 Sept 16 to 2020 Jan 4. The arrow indicates epochs when spectra were acquired. |
The U and B magnitudes are closely correlated. However, the strong increase in the U magnitude, which could be interpreted as the continuum emission mostly coming from the accretion disk with the contribution of the Balmer continuum coming from the BLR, does not have the same amplitude in the B magnitude, which also contains the Hβ line. This is most likely due to the host-galaxy contribution, which is much stronger in the B band. Looking at the light curves it seems there is little or no time delay between the U and B band. A more detailed analysis of the light-curve properties (e.g., time delay and oscillations; see Kovačević et al. 2018) will be given in a separate publication.
The most striking change in the flux is seen when the UV and X-ray data are considered (Fig. 1). In the 2012−2014 intensive monitoring campaign, the UVW2 magnitude changed by more than 2 mag, while B changed by only ∼0.3 mag. The largest amplitude in the available data was ∼3 mag from 2006 to 2015. The X-ray data show similar behavior; that is, the difference between the minimum in 2014 and a rise in 2016 in the optical band is very small, on the order of ∼0.2 mag, while the X-ray flux changed significantly by about a factor of five. Most importantly, the optical variations are similar and closely follow the trends seen in the UV and X-ray, therefore we may expect that in the UV/X-ray bands the object is experiencing a much more abrupt transition.
We obtained the optical spectra of NGC 3516 in order to follow the flare spotted in the photometric data. To extract the broad line profile, we subtracted the host-galaxy spectrum using the principal component analysis (PCA), a statistical method that was described in Vanden Berk et al. (2006), and was used in the previous campaign (for details see Shapovalova et al. 2019). Figure 3 shows the result of the host galaxy fitting for the spectrum obtained on Dec 09, 2019. The positions of prominent lines are also marked in Fig. 3.
![]() |
Fig. 3. Host galaxy fitting using PCA in the case of the spectrum observed on 2019 Dec 07. The position of all prominent emission lines is marked. |
After subtracting the host galaxy, several forbidden high-ionization lines are clearly seen (see Fig. 3), which are also visible in the case of spectra with different spectral resolution (see Fig. 4). The strongest are [Fe VII] λ4893, λ5721, λ6087, [Fe X] λ6374, and [Ne II] λ3869. These so-called coronal lines are probably originating from the inner narrow-line region or from the inner edge of the torus (e.g., Peterson 1988; Rose et al. 2015), and for their creation most likely a hot X-ray wind is needed (Pier & Voit 1995) or some special geometry, like polar conical regions (Oknyanskii et al. 1991). The variability of coronal lines has been detected before (e.g., Landt et al. 2015a,b), and their strengthening has been shown to be a signature of transition from a low to a high-activity state in AGNs (Oknyansky et al. 2019b). Several other examples of CL AGNs that have displayed variable coronal line emission are known, such as for example J000904.54−103428.7, J013203 (MacLeod et al. 2019), NGC 5548 (Landt et al. 2015b; Fausnaugh et al. 2016), Mkn 110 (Kollatschny & Bischoff 2002), NGC 4151 (Oknyanskii et al. 1991; Landt et al. 2015a), and ZTF18aajupnt (Frederick et al. 2019).
![]() |
Fig. 4. Two spectra observed in 2019, the Hβ (left) and Hα (right) spectral regions are shown. We note that the spectra are with the original different spectral resolution and are scaled for comparison. |
Furthermore, we modeled the broad line profiles with multiGaussian components (following the procedures and parameter criteria described in e.g., Popović et al. 2004, 2014; Ilić et al. 2006; Dimitrijević et al. 2007; Kovačević et al. 2010; Shapovalova et al. 2012, and references therein), using the newly developed python-based AGN line-fitting code (see Rakić, in prep.). This code simultaneously fits the underlying continuum and all emission lines, and is based on the Sherpa python package (Burke et al. 2019). In Fig. 5 we plot the result of the multicomponent fitting of the Hα (up) and Hβ (bottom). Due to its complex profile, the Hα broad line is modeled with four components (blue and magenta lines in Fig. 5), whereas Hβ was reproduced with only two components due to its low S/N. Both lines show asymmetric profiles with a blueshifted peak and very extended wings that are better seen in Hα line because of the higher S/N.
![]() |
Fig. 5. Multicomponent fitting of the Hα (top) and Hβ (bottom) lines. Below the observed (blue) and modeled (orange) spectrum, the residual (gray) and all Gaussian components are shown (blue, magenta – broad components of Hα (up) and Hβ; green – narrow lines; red – forbidden [Fe X] λ6374 and [Fe VII] λ4893 lines). The thick dashed line indicates the modeled broad line; this is the sum of all broad Gaussian components. |
The blended profile of [O I] λ6363 and [Fe X] λ6374 was also fitted (see Fig. 5, upper panel). From the fit, the estimated ratio of [O I]+[Fe X] to [O I] λ6300 is ∼1.2, similar that observed in NGC 1566 during the brightening phase (Oknyansky et al. 2019b), but we note that it is strongly dependent on the estimated level of the underlying continuum and Hα extended wing. From the clear single profile of [Fe VII] λ6087 (see Fig. 3), we estimated that the line fluxes have increased by the factor of approximately two compared to the last published data of Shapovalova et al. (2019). The presence of strong coronal [Fe VII] and [Fe X] lines could be connected with the TDEs (e.g., Yang et al. 2013) but more intense [Fe X] would be expected Wang et al. (2012). On the other hand, their strengthening can be a tracer of an awaking phase, as recently shown in a “turn-on” of broad-line AGNs detected in a normal galaxy (Yan et al. 2019).
Figure 6 shows the broad-line profiles of Hα and Hβ lines; these have been normalized for comparison. The lines show the same double-peaked profile that is better seen in the Hα line because of higher S/N in the spectra. They have a similar full width half maximum (FWHM) of ∼3400 km s−1 (Fig. 6). Compared to the results from our previous long-term campaign (broad Hα profile is also shown in Fig. 6), the FWHM is on the order of the average FWHM during the whole 22 year campaign (Shapovalova et al. 2019). Both broad components are still significantly shifted to the blue, as was noticed in 2017 when the broad component started to appear. The lines were then blueshifted for around 1000 km s−1, and it may be that the blue peak has changed its position, which requires verification with newer continuous spectral observations. In the Hα profile, signatures of the blue and red bumps are noticeable, whereas in both lines the very broad component of ∼10 000 km s−1 in width seems to appear. These features indicate that the BLR of NGC 3516 is complex and that broad line emission is present both from the outflowing material and the disk, as suggested by Popović et al. (2002).
![]() |
Fig. 6. Comparison of the normalized Hα and Hβ broad-line profiles from 2019. The wavelength has been converted to the velocity scale. The best-quality broad-line profile of Hα from the Shapovalova et al. (2019) campaign obtained in February 2017 is also plotted. |
4. Discussion
During the high-cadence fourth-month optical monitoring campaign of the prototypical changing-look Seyfert galaxy NGC 3516, we detected an increase in the optical photometric light curves and spectral features including the brightening of coronal Fe lines and a proposed emergent broad Balmer feature. We strengthened our analysis with previous optical photometric data (1999–2018), as well as with Swift X-ray and UV archival monitoring data of NGC 3516 from 2006–2018.
We find that the change in the U and B band is ∼ 0.25 mag and 0.11 mag, respectively. This could be considered as ordinary broad-line AGN variability (e.g., Vanden Berk et al. 2004; Sesar et al. 2007) and one could expect an increase of around ∼1 mag or even more (e.g., Graham et al. 2017, 2020; MacLeod et al. 2019; Rumbaugh et al. 2018) in the case of a peculiar flaring (or an extreme brightening). However, a high increase of ∼1 mag is more likely to happen in bright quasars, and this more strict limit is used as a flaring criterion in large surveys, such as for example the Sloan Digital Sky Survey in which the magnitude uncertainties can be about 0.2 mag or more. However, here we performed high-precision photometry toward NGC 3516 achieving ∼0.01 mag uncertainties, and therefore we are able to detect real magnification in the object brightness. Moreover, one would expect the changes in the U and B of the CL AGN to be significantly smaller than that in high-energy bands, i.e. the dramatic change in the X-ray and far UV band is probably followed by a much lower change in the optical, since the host galaxy has a very strong contribution in the optical spectral band. This is clearly seen in the case of NGC 1566, which is the closest CL AGN (Oknyansky et al. 2019b). At the distance of NGC 3516, using the same aperture as for NGC 1566, there is a significantly higher contribution of the host galaxy to the NGC 3516 optical flux. Comparing the NGC 1566 dramatic variations of ∼1 mag observed in 2018 (see Oknyansky et al. 2019b) with the NGC 3516 ones reported in this paper, if NGC 1566 were at the same distance as NGC 3516, then for the same aperture the NGC 1566 B-band brightness would change by much less than 0.2 mag. Furthermore, we would not be able to detect these large variations of ∼1 mag, but only a smaller change of 0.2–0.3 mag. This suggests that real small-amplitude optical variations in NGC 3516 could be supported with much stronger variations of a few magnitudes in the X-ray and far UV band, and could indicate that the object is going through a transition. Therefore, further multiwavelength observations are required to distinguish that this optical variability is beyond ordinary AGN behavior and that it could be linked to a change in the state of AGN activity.
One of the signatures of a transition from low to high activity is the presence of strong coronal lines, of which [Fe X] is of particular interest. Intense X-ray flux is needed for this coronal line to be strong, because we need 0.2 keV energy photons to ionize Fe8+ to Fe9+. We detected this line earlier, in our long-term spectroscopic campaign (see, e.g. Figs. 2 and 3 in Shapovalova et al. 2019), although it was not mentioned in the publication. It appears that in the highest state of activity in 2007 (see Fig. 12 in Shapovalova et al. 2019) the line was not particularly prominent. This could be due to either the poorer quality of the spectra, and the fact that Hα is much stronger, making it difficult to detect [Fe X], or some physical phenomenon, which future monitoring of coronal-line variability will reveal.
5. Conclusions
The changing-look AGN NGC 3516 was monitored intensively in U and B filters over a period of four months, and these observations were supported with additional spectral observations once the increase in brightness was spotted. We performed accurate photometric measurements and performed a preliminary analysis of the line profiles. We support our analysis of the past behavior of NGC 3516 with archival Swift data in UV and X-ray bands. Our main findings can be summarized as follows.
-
We detected a general increase in brightness compared to the last three years of photometric data. The latest flare, at the end of 2019, is seen in both U and B magnitude to a maximum amplitude of 0.25 mag and 0.11 mag, respectively. The variations in the U and B magnitude are well correlated, but the amplitude of variability in U is significantly greater than in B, which is partially connected with the stronger contribution of the host galaxy in the B band.
-
Several coronal lines, such as the forbidden high-ionization iron lines [Fe VII] and [Fe X], are clearly detected in the NGC 3516 spectrum observed in 2019, and increased in flux compared to the end of Shapovalova et al. (2019) campaign. The strengthening of these lines may indicate that the AGN is entering a more active phase (see, e.g., Yan et al. 2019).
-
The broad Balmer lines show a complex double-peaked structure, with an asymmetric profile and prominent blue peak (similar to in Shapovalova et al. 2019). However, it seems that a very broad component of ∼10 000 km s−1 width is starting to appear, which may also indicate that the AGN is experiencing a transition to a higher activity phase.
-
The optical variations are similar and closely follow the trends seen in UV and X-rays, and therefore we may expect that in UV and X-rays the object is experiencing a more abrupt transition. However, it is not sufficiently empirically demonstrated that a flare in the X-rays is coupled to a transition in the presence or absence of emission lines (LaMassa et al. 2017), and therefore a final conclusion may only be possible with future observations and additional analysis.
The photometric and spectroscopic observations presented in this paper indicate that the AGN of NGC 3516 may be in a transition phase, changing from a low-activity state to one of high activity. Further intensive multiwavelength (optical, UV, and X-ray) monitoring in photometry and spectroscopy is needed and may help to elucidate the processes behind the changing-look phenomenon in AGNs.
Acknowledgments
This work made use of data supplied by the UK Swift Science Data Centre at the University of Leicester. Part of this work was supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia through the project number 176001, and from the M.V. Lomonosov Moscow State University Program of Development (RC600 & TDS). BAA, EOM, SAP, NIS, AVD acknowledge the support by the Program of development of M.V. Lomonosov Moscow State University (Leading Scientific School “Physics of stars, relativistic objects and galaxies”). EOM, NIS, AVD, AMT, SGJ are also supported by the RSF grant 17-12-01241. ST acknowledges financial support from the Russian Foundation for Basic Research project 17-52-80139 BRICS-a. We also express our thanks to the anonymous referee for helpful comments and suggestions that greatly improved this work.
References
- Andrillat, Y., & Souffrin, S. 1968, Astrophys. Lett., 1, 111 [NASA ADS] [Google Scholar]
- Berdnikov, L. N., Belinski, A. A., Shatsky, N. I., et al. 2020, Astron. Rep., 64, 310 [NASA ADS] [CrossRef] [Google Scholar]
- Buisson, D. J. K., Lohfink, A. M., Alston, W. N., & Fabian, A. C. 2017, MNRAS, 464, 3194 [NASA ADS] [CrossRef] [Google Scholar]
- Burke, D., Laurino, O., dtnguyen2, et al. 2019, https://doi.org/10.5281/zenodo.3358134 [Google Scholar]
- Campana, S., Mainetti, D., Colpi, M., et al. 2015, A&A, 581, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Denney, K. D., Peterson, B. M., Pogge, R. W., et al. 2010, ApJ, 721, 715 [NASA ADS] [CrossRef] [Google Scholar]
- Denney, K. D., De Rosa, G., Croxall, K., et al. 2014, ApJ, 796, 134 [NASA ADS] [CrossRef] [Google Scholar]
- De Rosa, G., Fausnaugh, M. M., Grier, C. J., et al. 2018, ApJ, 866, 133 [NASA ADS] [CrossRef] [Google Scholar]
- Dimitrijević, M. S., Popović, L. Č., Kovačević, J., Dačić, M., & Ilić, D. 2007, MNRAS, 374, 1181 [NASA ADS] [CrossRef] [Google Scholar]
- Edelson, R., Koratkar, A., Nandra, K., et al. 2000, ApJ, 534, 180 [NASA ADS] [CrossRef] [Google Scholar]
- Elitzur, M. 2012, ApJ, 747, L33 [NASA ADS] [CrossRef] [Google Scholar]
- Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2009, MNRAS, 397, 1177 [NASA ADS] [CrossRef] [Google Scholar]
- Fausnaugh, M. M., Denney, K. D., Barth, A. J., et al. 2016, ApJ, 821, 56 [NASA ADS] [CrossRef] [Google Scholar]
- Frederick, S., Gezari, S., Graham, M. J., et al. 2019, ApJ, 883, 31 [NASA ADS] [CrossRef] [Google Scholar]
- Gehrels, N., Chincarini, G., Giommi, P., et al. 2004, ApJ, 611, 1005 [NASA ADS] [CrossRef] [Google Scholar]
- Graham, M. J., Djorgovski, S. G., Drake, A. J., et al. 2017, MNRAS, 470, 4112 [NASA ADS] [CrossRef] [Google Scholar]
- Graham, M. J., Ross, N. P., Stern, D., et al. 2020, MNRAS, 491, 4925 [NASA ADS] [Google Scholar]
- Huerta, E. M., Krongold, Y., Nicastro, F., et al. 2014, ApJ, 793, 61 [NASA ADS] [CrossRef] [Google Scholar]
- Ilić, D., Popović, L. Č., Bon, E., Mediavilla, E. G., & Chavushyan, V. H. 2006, MNRAS, 371, 1610 [NASA ADS] [CrossRef] [Google Scholar]
- Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kollatschny, W., & Bischoff, K. 2002, A&A, 386, L19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kollatschny, W., & Fricke, K. J. 1985, A&A, 146, L11 [NASA ADS] [Google Scholar]
- Komossa, S., Grupe, D., Schartel, N., et al. 2017, in New Frontiers in Black Hole Astrophysics, ed. A. Gomboc, IAU Symp., 324, 168 [NASA ADS] [Google Scholar]
- Kovačević, J., Popović, L. Č., & Dimitrijević, M. S. 2010, ApJS, 189, 15 [NASA ADS] [CrossRef] [Google Scholar]
- Kovačević, A. B., Pérez-Hernández, E., Popović, L. Č., et al. 2018, MNRAS, 475, 2051 [NASA ADS] [CrossRef] [Google Scholar]
- Kraemer, S. B., Crenshaw, D. M., George, I. M., et al. 2002, ApJ, 577, 98 [NASA ADS] [CrossRef] [Google Scholar]
- LaMassa, S. M., Yaqoob, T., & Kilgard, R. 2017, ApJ, 840, 11 [NASA ADS] [CrossRef] [Google Scholar]
- Landt, H., Ward, M. J., Steenbrugge, K. C., & Ferland, G. J. 2015a, MNRAS, 449, 3795 [NASA ADS] [CrossRef] [Google Scholar]
- Landt, H., Ward, M. J., Steenbrugge, K. C., & Ferland, G. J. 2015b, MNRAS, 454, 3688 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, T., Wang, J.-X., Yang, H., Zhu, F.-F., & Zhou, Y.-Y. 2014, ApJ, 783, 106 [NASA ADS] [CrossRef] [Google Scholar]
- Lyutyi, V. M. 1972, Astron. Zh, 49, 930 [Sov. Astron. 16, 763 (1972)] [NASA ADS] [Google Scholar]
- Lyutyi, V. M., & Doroshenko, V. T. 1993, Astron. Lett., 19, 405 [NASA ADS] [Google Scholar]
- Lyuty, V. M., Oknyansky, V. L., & Chuvaev, K. K. 1984, Pisma v Astron. Zh, 10, 803 [NASA ADS] [Google Scholar]
- MacLeod, C. L., Ross, N. P., Lawrence, A., et al. 2016, MNRAS, 457, 389 [NASA ADS] [CrossRef] [Google Scholar]
- MacLeod, C. L., Green, P. J., Anderson, S. F., et al. 2019, ApJ, 874, 8 [NASA ADS] [CrossRef] [Google Scholar]
- Maoz, D., Markowitz, A., Edelson, R., & Nandra, K. 2002, AJ, 124, 1988 [NASA ADS] [CrossRef] [Google Scholar]
- Mathur, S., Denney, K. D., Gupta, A., et al. 2018, ApJ, 866, 123 [NASA ADS] [CrossRef] [Google Scholar]
- Mehdipour, M., Branduardi-Raymont, G., & Page, M. J. 2010, A&A, 514, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Netzer, H., Chelouche, D., George, I. M., et al. 2002, ApJ, 571, 256 [NASA ADS] [CrossRef] [Google Scholar]
- Noda, H., & Done, C. 2018, MNRAS, 480, 3898 [NASA ADS] [CrossRef] [Google Scholar]
- Noda, H., Minezaki, T., Watanabe, M., et al. 2016, ApJ, 828, 78 [NASA ADS] [CrossRef] [Google Scholar]
- Oknyanskii, V. L., Lyutyi, V. M., & Chuvaev, K. K. 1991, Sov. Astron. Lett., 17, 100 [NASA ADS] [Google Scholar]
- Oknyansky, V. L., Gaskell, C. M., Huseynov, N. A., et al. 2017, MNRAS, 467, 1496 [NASA ADS] [Google Scholar]
- Oknyansky, V. L., Winkler, H., Tsygankov, S. S., et al. 2019a, Odessa Astron. Publ., 32, 75 [NASA ADS] [CrossRef] [Google Scholar]
- Oknyansky, V. L., Winkler, H., Tsygankov, S. S., et al. 2019b, MNRAS, 483, 558 [NASA ADS] [CrossRef] [Google Scholar]
- Peterson, B. M. 1988, PASP, 100, 18 [NASA ADS] [CrossRef] [Google Scholar]
- Pier, E. A., & Voit, G. M. 1995, in Oxford Torus Workshop, ed. M. J. Ward, 93 [Google Scholar]
- Popović, L. Č., Mediavilla, E. G., Kubičela, A., & Jovanović, P. 2002, A&A, 390, 473 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Popović, L. Č., Mediavilla, E., Bon, E., & Ilić, D. 2004, A&A, 423, 909 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Popović, L. Č., Shapovalova, A. I., Ilić, D., et al. 2014, A&A, 572, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rose, M., Elvis, M., & Tadhunter, C. N. 2015, MNRAS, 448, 2900 [NASA ADS] [CrossRef] [Google Scholar]
- Rumbaugh, N., Shen, Y., Morganson, E., et al. 2018, ApJ, 854, 160 [NASA ADS] [CrossRef] [Google Scholar]
- Runco, J. N., Cosens, M., Bennert, V. N., et al. 2016, ApJ, 821, 33 [NASA ADS] [CrossRef] [Google Scholar]
- Sesar, B., Ivezić, Ž., Lupton, R. H., et al. 2007, AJ, 134, 2236 [NASA ADS] [CrossRef] [Google Scholar]
- Seyfert, C. K. 1943, ApJ, 97, 28 [NASA ADS] [CrossRef] [Google Scholar]
- Shapovalova, A. I., Popović, L. Č., Burenkov, A. N., et al. 2012, ApJS, 202, 10 [NASA ADS] [CrossRef] [Google Scholar]
- Shapovalova, A. I., Popović, L. Č., Afanasiev, V. L., et al. 2019, MNRAS, 485, 4790 [NASA ADS] [CrossRef] [Google Scholar]
- \’{S}niegowska, M., & Czerny, B. 2019, ApJL, submitted [arXiv:1904.06767] [Google Scholar]
- Stern, D., McKernan, B., Graham, M. J., et al. 2018, ApJ, 864, 27 [NASA ADS] [CrossRef] [Google Scholar]
- Storchi-Bergmann, T., Schimoia, J. S., Peterson, B. M., et al. 2017, ApJ, 835, 236 [NASA ADS] [CrossRef] [Google Scholar]
- Turner, T. J., Kraemer, S. B., George, I. M., Reeves, J. N., & Bottorff, M. C. 2005, ApJ, 618, 155 [NASA ADS] [CrossRef] [Google Scholar]
- Turner, T. J., Miller, L., Kraemer, S. B., & Reeves, J. N. 2011, ApJ, 733, 48 [NASA ADS] [CrossRef] [Google Scholar]
- Vanden Berk, D. E., Wilhite, B. C., Kron, R. G., et al. 2004, ApJ, 601, 692 [NASA ADS] [CrossRef] [Google Scholar]
- Vanden Berk, D. E., Shen, J., Yip, C.-W., et al. 2006, AJ, 131, 84 [NASA ADS] [CrossRef] [Google Scholar]
- Wachter, K., Leach, R., & Kellogg, E. 1979, ApJ, 230, 274 [Google Scholar]
- Wanders, I., & Horne, K. 1994, A&A, 289, 76 [NASA ADS] [Google Scholar]
- Wanders, I., van Groningen, E., Alloin, D., et al. 1993, A&A, 269, 39 [NASA ADS] [Google Scholar]
- Wang, T.-G., Zhou, H.-Y., Komossa, S., et al. 2012, ApJ, 749, 115 [NASA ADS] [CrossRef] [Google Scholar]
- Yan, L., Wang, T., Jiang, N., et al. 2019, ApJ, 874, 44 [NASA ADS] [CrossRef] [Google Scholar]
- Yang, C.-W., Wang, T.-G., Ferland, G., et al. 2013, ApJ, 774, 46 [CrossRef] [Google Scholar]
All Figures
![]() |
Fig. 1. Multiwavelength light curves of NGC 3516. Top panel: long-term photometric optical and UV light curves in B and UVW2 filters: circles denote data from Shapovalova et al. (2019), squares denote the result from this paper (SAI CMO), diamonds stand for the data from SAI Crimea campaign shown for continuity, and asterisks denote UVW2 UVOT Swift photometry. All magnitudes are reduced to the same aperture of 10″. A dashed horizontal line at magnitude 14.1 is given at the bottom of the figure to guide the eye, and the arrow indicates epochs when spectra were acquired. Bottom panel: Swift/XRT 0.5−10 keV X-ray flux light curve in erg cm−2 s−1. |
In the text |
![]() |
Fig. 2. Short-term light curves in U and B bands (aperture size 6.7″) from the high-cadence monitoring campaign from 2019 Sept 16 to 2020 Jan 4. The arrow indicates epochs when spectra were acquired. |
In the text |
![]() |
Fig. 3. Host galaxy fitting using PCA in the case of the spectrum observed on 2019 Dec 07. The position of all prominent emission lines is marked. |
In the text |
![]() |
Fig. 4. Two spectra observed in 2019, the Hβ (left) and Hα (right) spectral regions are shown. We note that the spectra are with the original different spectral resolution and are scaled for comparison. |
In the text |
![]() |
Fig. 5. Multicomponent fitting of the Hα (top) and Hβ (bottom) lines. Below the observed (blue) and modeled (orange) spectrum, the residual (gray) and all Gaussian components are shown (blue, magenta – broad components of Hα (up) and Hβ; green – narrow lines; red – forbidden [Fe X] λ6374 and [Fe VII] λ4893 lines). The thick dashed line indicates the modeled broad line; this is the sum of all broad Gaussian components. |
In the text |
![]() |
Fig. 6. Comparison of the normalized Hα and Hβ broad-line profiles from 2019. The wavelength has been converted to the velocity scale. The best-quality broad-line profile of Hα from the Shapovalova et al. (2019) campaign obtained in February 2017 is also plotted. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.