Original Research ARTICLE
Long-Term Monitoring of the Broad-Line Region Properties in a Selected Sample of AGN
- 1Department of Astronomy, Faculty of Mathematics, University of Belgrade, Belgrade, Serbia
- 2Special Astrophysical Observatory, Russian Academy of Sciences, Nizhnii Arkhyz, Russia
- 3Astronomical Observatory, Belgrade, Serbia
- 4Instituto Nacional de Astrofísica, Óptica y Electrónica, Puebla, Mexico
- 5Institut fuer Astrophysik, Universitaet Goettingen, Göttingen, Germany
- 6Faculty of Science, University of Banjaluka, Banjaluka, Republic of Srpska, Bosnia and Herzegovina
- 7Department of Physics and Astronomy, University of Padova, Padova, Italy
We present the results of the long-term optical monitoring campaign of active galactic nuclei (AGN) coordinated by the Special Astrophysical Observatory of the Russian Academy of Science. This campaign has produced a remarkable set of optical spectra, since we have monitored for several decades different types of broad-line (type 1) AGN, from a Seyfert 1, double-peaked line, radio loud and radio quiet AGN, to a supermassive binary black hole candidate. Our analysis of the properties of the broad line region (BLR) of these objects is based on the variability of the broad emission lines. We hereby give a comparative review of the variability properties of the broad emission lines and the BLR of seven different type 1 AGNs, emphasizing some important results, such as the variability rate, the BLR geometry, and the presence of the intrinsic Baldwin effect. We are discussing the difference and similarity in the continuum and emission line variability, focusing on what is the impact of our results to the supermassive black hole mass determination from the BLR properties.
The broad line region (BLR) in active galactic nuclei (AGN) is responsible for the emission of intensive broad lines, seen in the optical spectra of type 1 AGN. We use these lines to probe the physics and kinematics of the BLR, aiming to study in more details the central engine and the supermassive black hole (SMBH) properties. This is because, on one hand the BLR is photoionized by the central continuum source, and on the other, it is gravitationally bound to the SMBH (Osterbrock and Ferland, 2006).
One way to do this is by using the variability of the broad emission lines, as they are known to be strongly variable in flux and profile. From the long-term monitoring campaign of optical broad emission lines in different type 1 AGN, it is possible to estimate the size of the BLR (reverberation analysis Lyutyi and Cherepashchuk, 1972; Blandford and McKee, 1982; Gaskell and Sparke, 1986) and the mass of the SMBH (Gaskell, 1988), and test the photoionization hypothesis by studying the line and continuum flux correlations (Osterbrock and Ferland, 2006). One particularly interesting correlation, first detected in the broad C IV line, is the Baldwin effect, i.e., the anti-correlation between the emission line equivalent width and the near-by continuum luminosity (Baldwin, 1977). There are two different types of Baldwin effect: (i) global Baldwin effect which shows the above anti-correlation on single-epoch observations of a large number of AGN; and (ii) intrinsic Baldwin effect which is the same anti-correlation detected in an individual variable AGN. It is well known that the global Baldwin effect is not detected in case of the broad Hβ line. However, when testing this correlation in case of a single object monitoring data, i.e., the intrinsic Baldwin effect, it has been seen that the anti-correlation exists even in the case of Hβ (Pogge and Peterson, 1992; Rakić et al., 2017).
Here we present some interesting findings of the long-term (more than a decade) optical monitoring campaign of a sample of seven type 1 AGN of different spectral types: Seyfert 1 galaxies (NGC 5548, NGC 4151, NGC 7469), Narrow-line Seyfert 1 galaxy - NLSy 1 (Ark 564), double-peaked line radio loud (3C 390.3) and radio quiet (Arp 102B) galaxy, and a luminous quasar (E1821+643), a supermassive binary black hole candidate. All spectral data for each objects have been first presented and analyzed separately in Shapovalova et al. (2001, 2004, 2008, 2010a, 2012, 2013, 2016, 2017). However, since the obtained spectra are observed, reduced and analyzed in the same manner, the final result is a uniform set of data, which can be compared and discussed. Here we continue the discussion given in Ilić et al. (2015), where we first tried to give a comparative review of the variability properties of type 1 AGN in our sample, only that now our sample consist of seven objects, including one high-luminosity local quasar, which makes our sample more diverse. Here we focus on discussing the possible BLR geometry and the presence of the intrinsic Baldwin effect. The paper is organized as follows: in Section 2 we briefly describe the observations, the obtained data and performed analysis; in Section 3 we present the results, and in Section 4 we give the discussion and outline our conclusions.
2. Observations and Data Analysis
The long-term monitoring campaign was coordinated by the Special Astrophysical Observatory of the Russian Academy of Science (SAO), using their 1-m and 6-m telescopes. Additionally, the campaign made use of the 2.1-m telescope of Guillermo Haro Astrophysical Observatory (Mexico), 2.1-m telescope of the Observatorio Astronomico Nacional at San Pedro Martir (Mexico), and the 3.5-m and 2.2-m telescopes of Calar Alto Observatory (Spain). The spectral data have been collected for several decades (since 1987 until 2015) and have been reduced and calibrated using the same procedures and methods, thus producing a homogeneous set of data. The information about the observations and the description of all procedures for each object are given in Shapovalova et al. (2001, 2004, 2008, 2010a, 2012, 2013, 2016, 2017). Here we outline some of the important line parameter measurements. In order to measure the broad line fluxes and full width half maximum (FWHM), we needed to subtract the underlying continuum emission, and the narrow and other near-by satellite lines, for which we performed the multi-Gaussian fittings (Popović et al., 2004). In case of Hβ line region, it is important to carefully model the Fe II multiplet emission, and this was done using a template1 described in details in Kovačević et al. (2010) and Shapovalova et al. (2012).
In Table 1 we list the spectral characteristics of the sample, together with the basic information (see caption of Table 1). The RBLR comes from the time-lags calculated using the Z-transformed Discrete Correlation Function (Alexander, 1997), which was applied on continuum and line flux light-curves. Depending on the sampling of light curves we used either observed or simulated light-curves in order to get the most reliable result (for details see references in column 9, Table 1). The listed masses are taken from the monitoring results (see column 8 in Table 1), and for details on how the mass was calculated see given references (column 9 in Table 1). In case of NGC 5548 and NGC 4151 the mass was not derived in the original monitoring campaign, so we calculated in this work their SMBH masses, using the equation , where ΔVFWHM is the orbital velocity estimated from the width of the variable part of Hβ (taken from references given in Table 1), and f is the virial factor taken to be 5.5 (Onken et al., 2004). The variability parameter Fvar gives the level of flux variability and is calculated using equation from O'Brien et al. (1998). For the luminosity at 5100 Å we calculated the luminosity distances with the online calculator of Wright (2006), taking the same cosmological parameters as in Ilić et al. (2015): Ω = 0.286, ΩΛ = 0.714, and Ωk = 0, and a Hubble constant, H0 = 69.6 km/s Mpc−1.
In this paper we present a sample of 7 type 1 AGN with different broad line profiles (Table 1, columns 3–4) and different variability rate (Table 1, column 6). Figure 1 shows the Hβ line flux light curves for all 7 objects, where the line flux and the errors-bars are shown normalized to the maximal flux for comparison. In all 7 objects, for which we have been collecting spectra for more than a decade and in some cases two decades, the line and continuum flux are varying. From Table 1 and Figure 1, we can see that: NGC 4151 and NGC 5548 are strongly variable Seyfert 1.5, NGC 7469 is also a Seyfert 1 but with lower variability in line fluxes and profiles, Arp 102B and 3C 390.3 have broad double-peaked lines, Ark 564 is a NLSy1, and E1821+643 is a quasar with strong red-asymmetry in line profiles, but almost no variability.
Figure 1. Light curves of the Hβ line flux of all object in the sample (object name given in the bottom-left), normalized to the maximal flux for better comparison.
One important prerequisite for the estimates of the radius of the BLR using the reverberation monitoring data, is the fact that the BLR is photoionized by the central continuum source (see e.g., Peterson and Horne, 2006). The simplest way to test this is by plotting the line luminosity as a function of the continuum luminosity. In Figure 2 we thus show the Hβ line luminosity for all 7 objects (shown with dots) vs. the continuum luminosity at 5100 Å. Moreover, the mean value is also shown (filled circle) with the object name denoted next to it. The important finding is the fact that the luminosity of the Hβ line is following the same trend with the continuum luminosity for all objects, from the low-luminosity Seyfert up to the high-luminosity quasar (Figure 2).
Figure 2. Hβ line luminosity for all 7 objects vs. the continuum luminosity at 5100 Å (dots). The mean value is also shown (filled circle) with the object name denoted next to it.
On the other hand, it was shown that the intrinsic Baldwin effect is detected in the Balmer lines (Hα and Hβ) in 6 objects from this sample (see the results and discussion in Rakić et al., 2017). Here we show as an example the presence of the intrinsic Baldwin effect, an anti-correlation between the Hβ line equivalent width and the continuum flux2 measured at 5100 Å) in case of Hβ of NGC 4151 (Figure 3), for which this effect is of high significance (see Rakić et al., 2017).
Figure 3. The intrinsic Baldwin effect, i.e., the Hβ line equivalent width vs. the continuum flux at 5100 Å, in case of Seyfert 1 galaxy NGC 4151.
Figure 4 shows the radius of the broad Hβ line emitting-region as a function of the continuum luminosity at 5100 Å for our sample (Table 1, columns 5 and 7, respectively). Data are fitted with the scaling relation from Bentz et al. (2006) in the form of , where α is the slope of the BLR radius-luminosity relationship and K is the scaling factor. The dashed line is plotted with the coefficients α = 0.533 and K = 1.527, given by Bentz et al. (2013), and the solid line is a simple linear best-fitting of the above equation (not considering the error-bars) through all 7 objects, obtaining the fitting parameters of α = 0.226 and K = 1.639. The linear regression fit was done using a nonlinear least-squares Marquardt-Levenberg algorithm available in gnuplot.
Figure 4. The radius of the broad Hβ line emitting-region vs. the continuum luminosity at 5100 Å. The dashed line gives the BLR radius-luminosity scaling relation from Bentz et al. (2013) with , whereas the solid line is the fit through seven objects from this study, with . We note that the fit is strongly influenced from the one point, the luminous quasar E1821+643, and without that point the relation is consistent with the Bentz et al. (2013) result.
4. Discussion and Conclusions
Here we have analyzed and compared the optical spectra for seven type 1 AGN which have different variability properties and characteristics of the broad emission line profiles. This remarkable set of spectral data has been collected for several decades (Table 1, column 1) within the SAO RAS long-term monitoring campaign. Following the comparative study given in Ilić et al. (2015), we extend here our sample to 7 objects and outline the most important results. We repeat again that all observational data sets were uniformly processed, and in that way could be compared.
If we compare the normalized light-curves shown in Figure 1, the strongest variable objects are NGC 5548, NGC 4151, and 3C 390.3 (see also in Table 1 the variability parameter Fvar which is ~40%), followed by Arp 102B and NGC 7469 (Fvar ~20%). The least variable objects are the NLSy1 Ark 564 and the highly luminous quasar E1821+643, for which the line and continuum flux are almost constant for decades (for the Hβ line, Fvar ~7%). The case of the quasar E1821+643 is specially puzzling since its broad emission lines show a peak shift (of the order ~1,000 km/s) and a strong asymmetry in the red wing, extended up to ~15,000 km/s (see Shapovalova et al., 2016, and their Figure 15). This object was suggested to be a candidate for a binary SMBH or a recoiling black hole after the collision of two SMBH (Robinson et al., 2010), however Shapovalova et al. (2016) found it difficult to explain the unchanging line-profiles on such a long time-scales with the binary SMBH hypothesis, suggesting that the broad-line emitting region may be connected with some sort of gas-rich clouds that are orbiting around the recoiling SMBH.
Interestingly, the rate of variability shows no trend with the SMBH mass (Table 1, column 8), i.e., the low variability is detected in case of AGN with small and large SMBH mass (see Table 1). This may be a result of the small number of objects in our sample, since it is well known that Seyfert galaxies, which are low luminosity type 1 AGNs have lower SMBH mass and stronger variability than quasars, which are high luminosity type 1 AGNs with much higher SMBH mass and much lower variability.
The luminosity of Hβ is following the same trend with the continuum luminosity for all objects, regardless of the luminosity. This is in support of the prediction that the photoionization by the central continuum source is mainly responsible for the heating and line production, which is one of the assumptions needed when calculating the BLR radius from reverberation mapping. However, if considering a single object, the correlation between the line and continuum emission can be violated, which is seen in the case of NGC 4151 and Ark 564 (see for details Shapovalova et al., 2008, 2012, respectively). Another correlation that is present and significant is the intrinsic Baldwin effect in the Hβ line. The physical origin of this effect is still unknown (see Rakić et al., 2017, for discussion), but some lines of evidence suggest the possibility that an additional optical continuum component, not related to the ionizing continuum, can be originated either in the BLR or in nuclear outflows, thus influencing the slope of the Hβ - continuum flux relation.
However, what is intriguing is the fact that all presented objects have different line profiles (Table 1, column 4), indicating that the geometry of the broad line-emitting region is probably different. For example, in case of 3C 390.3 it is clear that the BLR is following the accretion disk geometry (Popović et al., 2011), whereas in case of the other double-peak line emitter Arp 102B this is maybe not the case (Popović et al., 2014). Also, the Seyfert galaxies NGC 4151 and NGC 5548 could host a binary SMBH with two separate BLR (Bon et al., 2012, 2016; Li et al., 2016), which puts more uncertainty on their mass estimates from reverberation mapping. Finally, the quasar E1821+643 is much below the theoretical radius-luminosity relation predicted by the photoionization theory (), and very low-variable AGN in the optical lines and continuum, therefore this SMBH mass estimate is also more uncertain. The peculiarity of the BLR of this quasar is probably the reason why the slope of the radius-luminosity scaling relation is much smaller than Bentz et al. (2006) relation. It is worth noticing that in Figure 4, if we discard Ark 564 (with very narrow lines that are not typical for type 1 AGNs and with a very large error-bars in estimated time-delays) and E1821+643 (a quasar with an outstanding luminosity), the remaining objects will apparently closely follow the dependence of Bentz et al. (2013).
We summarize here the main conclusions of this study:
(i) the rate of variability is not connected with the geometry of the BLR and mass of the SMBH (e.g., Arp 102B and 3C390.3), however in case of the SMBH mass this may be a result of the small number of objects in our sample;
(ii) the luminosity of Hβ is correlated with the continuum luminosity for all AGN in our sample, and the intrinsic Baldwin effect is present;
(iii) the photoionization is in general the main line-production mechanism, however in some AGN (e.g., NGC 4151) there are some additional mechanism which contribute to the line (continuum) emission.
Finally we can conclude that the long-term spectroscopic monitoring campaigns are very important for the investigation of the BLR structure and the SMBH mass estimates.
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
RFBR (grants N97-02-17625 N00-02-16272, N03-02-17123, 06-02-16843, N09-02-01136, 12-02-00857a, 12-02-01237a, N15-02-02101); Ministry of Education, Science, Technology and Development, Republic of Serbia (project 176001 “Astrophysical Spectroscopy of Extragalactic Objects”); INTAS (grant N96-0328); DFG grants Ko857/32-2 and Ko857/33-1; Ministry of Science and Technology of R. Srpska (project “Investigation of supermassive binary black holes in the optical and X-ray spectra”).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer ALL declared a shared affiliation, with no collaboration, with one of the authors, VC, to the handling Editor.
2. ^Since the intrinsic Baldwin effect is detected in a single object, we can use a flux instead of luminosity.
Alexander, T. (1997). “Is AGN variability correlated with other AGN properties? ZDCF analysis of small samples of sparse light curves,” in Astronomical Time Series, Vol. 218 of Astrophysics and Space Science Library, eds D. Maoz, A. Sternberg, and E. M. Leibowitz, (Dordrecht: Springer), 163.
Bentz, M. C., Denney, K. D., Grier, C. J., Barth, A. J., Peterson, B. M., Vestergaard, M., et al. (2013). The low-luminosity end of the radius-luminosity relationship for active galactic nuclei. Astrophys. J. 767:149. doi: 10.1088/0004-637X/767/2/149
Bentz, M. C., Peterson, B. M., Pogge, R. W., Vestergaard, M., and Onken, C. A. (2006). The radius-luminosity relationship for active galactic nuclei: the effect of host-galaxy starlight on luminosity measurements. Astrophys. J. 644, 133–142. doi: 10.1086/503537
Bon, E., Jovanović, P., Marziani, P., Shapovalova, A. I., Bon, N., Borka Jovanović, V., et al. (2012). The first spectroscopically resolved sub-parsec orbit of a supermassive binary black hole. Astrophys. J. 759:118. doi: 10.1088/0004-637X/759/2/118
Bon, E., Zucker, S., Netzer, H., Marziani, P., Bon, N., Jovanović, P., et al. (2016). Evidence for periodicity in 43 year-long monitoring of NGC 5548. Astrophys. J. Suppl. Ser. 225:29. doi: 10.3847/0067-0049/225/2/29
Gaskell, C. M. (1988). Direct evidence for gravitational domination of the motion of gas Within one light-week of the central object in NGC 4151 and the determination of the mass of the probable black hole. Astrophys. J. 325, 114–118. doi: 10.1086/165986
Ilić, D., Popović, L. Č., Shapovalova, A. I., Burenkov, A. N., Chavushyan, V. H., and Kovačević, A. (2015). Line shape variability in a sample of AGN with broad lines. J. Astrophys. Astron. 36, 433–445. doi: 10.1007/s12036-015-9360-y
Jovanović, P., Popović, L. Č., Stalevski, M., and Shapovalova, A. I. (2010). Variability of the Hβ line profiles as an indicator of orbiting bright spots in accretion disks of quasars: a case study of 3C 390.3. Astrophys. J. 718, 168–176. doi: 10.1088/0004-637X/718/1/168
Kovačević, J., Popović, L. Č., and Dimitrijević, M. S. (2010). Analysis of optical Fe II emission in a sample of active galactic nucleus spectra. Astrophys. J. Suppl. 189, 15–36. doi: 10.1088/0067-0049/189/1/15
Kovačević, A., Popović, L. Č., Shapovalova, A. I., and Ilić, D. (2017). Periodicity in the continua and broad line curves of a quasar E1821+643. Astrophys. Space Sci. 362:31. doi: 10.1007/s10509-017-3009-z
Kovačević, A., Popović, L. Č., Shapovalova, A. I., Ilić, D., Burenkov, A. N., and Chavushyan, V. H. (2014). Time series analysis of active galactic nuclei: the case of Arp 102B, 3C 390.3, NGC 5548 and NGC 4051. Adv. Space Res. 54, 1414–1428. doi: 10.1016/j.asr.2014.06.025
Li, Y.-R., Wang, J.-M., Ho, L. C., Lu, K.-X., Qiu, J., Du, P., et al. (2016). Spectroscopic indication of a centi-parsec supermassive black hole binary in the galactic center of ngc 5548. Astrophys. J. 822:4. doi: 10.3847/0004-637X/822/1/4
O'Brien, P. T., Dietrich, M., Leighly, K., Alloin, D., Clavel, J., Crenshaw, D. M., et al. (1998). Steps toward determination of the size and structure of the broad-line region in active galactic nuclei. XIII. ultraviolet observations of the broad-line radio galaxy 3C 390.3. Astrophys. J. 509, 163–176. doi: 10.1086/306464
Onken, C. A., Ferrarese, L., Merritt, D., Peterson, B. M., Pogge, R. W., Vestergaard, M., et al. (2004). Supermassive black holes in active galactic nuclei. II. calibration of the black hole mass-velocity dispersion relationship for active galactic nuclei. Astrophys. J. 615, 645–651. doi: 10.1086/424655
Popović, L. Č., Mediavilla, E., Bon, E., and Ilić, D. (2004). Contribution of the Disk Emission to the Broad Emission Lines in AGNs: Two-component model. Astron. Astrophys. 423, 909–918. doi: 10.1051/0004-6361:20034431
Popović, L. Č., Shapovalova, A. I., Ilić, D., Burenkov, A. N., Chavushyan, V. H., Kollatschny, W., et al. (2014). Spectral optical monitoring of the double-peaked emission line AGN Arp 102B. II. Variability of the broad line properties. Astron. Astrophys. 572:A66. doi: 10.1051/0004-6361/201423555
Popović, L. Č., Shapovalova, A. I., Ilić, D., Kovačević, A., Kollatschny, W., Burenkov, A. N., et al. (2011). Spectral optical monitoring of 3C 390.3 in 1995-2007. II. Variability of the spectral line parameters. Astron. Astrophys. 528:A130. doi: 10.1051/0004-6361/201016317
Rakić, N., La Mura, G., Ilić, D., Shapovalova, A. I., Kollatschny, W., Rafanelli, P., et al. (2017). The intrinsic Baldwin effect in broad Balmer lines of six long-term monitored AGNs. Astron. Astrophys. 603:A49. doi: 10.1051/0004-6361/201630085
Robinson, A., Young, S., Axon, D. J., Kharb, P., and Smith, J. E. (2010). Spectropolarimetric evidence for a kicked supermassive black hole in the Quasar E1821+643. Astrophys. J. lett. 717, L122–L126. doi: 10.1088/2041-8205/717/2/L122
Shapovalova, A. I., Burenkov, A. N., Carrasco, L., Chavushyan, V. H., Doroshenko, V. T., Dumont, A. M., et al. (2001). Intermediate resolution Hβ spectroscopy and photometric monitoring of 3C 390.3. I. Further evidence of a nuclear accretion disk. Astron. Astrophys. 376, 775–792. doi: 10.1051/0004-6361:20011011
Shapovalova, A. I., Doroshenko, V. T., Bochkarev, N. G., Burenkov, A. N., Carrasco, L., Chavushyan, V. H., et al. (2004). Profile variability of the Hα and Hβ broad emission lines in NGC 5548. Astron. Astrophys. 422, 925–940. doi: 10.1051/0004-6361:20035652
Shapovalova, A. I., Popović, L. Č., Burenkov, A. N., Chavushyan, V. H., Ilić, D., Kollatschny, W., et al. (2010a). Spectral optical monitoring of 3C 390.3 in 1995-2007. I. Light curves and flux variation in the continuum and broad lines. Astron. Astrophys. 517:A42. doi: 10.1051/0004-6361/201014118
Shapovalova, A. I., Popović, L. Č., Burenkov, A. N., Chavushyan, V. H., Ilić, D., Kollatschny, W., et al. (2013). Spectral optical monitoring of a double-peaked emission line AGN Arp 102B. Variability of spectral lines and continuum. Astron. Astrophys. 559:A10. doi: 10.1051/0004-6361/201321781
Shapovalova, A. I., Popović, L. Č., Burenkov, A. N., Chavushyan, V. H., Ilić, D., Kovačević, A., et al. (2010b). Long-term variability of the optical spectra of NGC 4151. II. Evolution of the broad Hα and Hβ emission-line profiles. Astron. Astrophys. 509:A106. doi: 10.1051/0004-6361/200912311
Shapovalova, A. I., Popović, L. Č., Burenkov, A. N., Chavushyan, V. H., Ilić, D., Kovačević, A., et al. (2012). Spectral Optical Monitoring of the Narrow-line Seyfert 1 Galaxy Ark 564. Astron. Astrophys. Suppl. Ser. 202:10. doi: 10.1088/0067-0049/202/1/10
Shapovalova, A. I., Popović, L. Č., Chavushyan, V. H., Afanasiev, V. L., Ilić, D., Kovačević, A., et al. (2017). Long-term optical spectral monitoring of NGC 7469. Month. Notices R. Astronom. Soc. 466, 4759–4775. doi: 10.1093/mnras/stx025
Shapovalova, A. I., Popović, L. Č., Chavushyan, V. H., Burenkov, A. N., Ilić, D., Kollatschny, W., et al. (2016). First long-term optical spectral monitoring of a binary black hole candidate E1821+643. I. Variability of spectral lines and continuum. Astron. Astrophys. Suppl. Ser. 222:25. doi: 10.3847/0067-0049/222/2/25
Shapovalova, A. I., Popović, L. Č., Collin, S., Burenkov, A. N., Chavushyan, V. H., Bochkarev, N. G., et al. (2008). Long-term variability of the optical spectra of NGC 4151. I. Light curves and flux correlations. Astron. Astrophys. 486, 99–111. doi: 10.1051/0004-6361:20079111
Keywords: galaxies:active-galaxies, quasar:individual (Arp 102B, 3C 390.3, NGC 5548, NGC 4151, NGC 7469, Ark 564, E1821+643), quasars: supermassive black holes, quasar:emission lines, line:profiles
Citation: Ilić D, Shapovalova AI, Popović LČ, Chavushyan V, Burenkov AN, Kollatschny W, Kovačević A, Marčeta-Mandić S, Rakić N, La Mura G and Rafanelli P (2017) Long-Term Monitoring of the Broad-Line Region Properties in a Selected Sample of AGN. Front. Astron. Space Sci. 4:12. doi: 10.3389/fspas.2017.00012
Received: 10 July 2017; Accepted: 21 August 2017;
Published: 14 September 2017.
Edited by:Paola Marziani, National Institute for Astrophysics (INAF), Italy
Reviewed by:Shai Kaspi, Tel Aviv University, Israel
Anna Lia Longinotti, National Institute of Astrophysics, Optics and Electronics, Mexico
Copyright © 2017 Ilić, Shapovalova, Popović, Chavushyan, Burenkov, Kollatschny, Kovačević, Marčeta-Mandić, Rakić, La Mura and Rafanelli. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dragana Ilić, firstname.lastname@example.org