Abstract
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.
1. Introduction
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, ).
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, ; Blandford and McKee, ; Gaskell and Sparke, ) and the mass of the SMBH (Gaskell, ), and test the photoionization hypothesis by studying the line and continuum flux correlations (Osterbrock and Ferland, ). 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, ). 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, ; Rakić et al., ).
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. (, , , , , , , ). 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. (), 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. (, , , , , , , ). 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., ). 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. () and Shapovalova et al. ().
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, ), 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., ). The variability parameter Fvar gives the level of flux variability and is calculated using equation from O'Brien et al. (). For the luminosity at 5100 Å we calculated the luminosity distances with the online calculator of Wright (), taking the same cosmological parameters as in Ilić et al. (): Ω = 0.286, ΩΛ = 0.714, and Ωk = 0, and a Hubble constant, H0 = 69.6 km/s Mpc−1.
Table 1
| Object period [years] (1) | z (2) | AGN type (3) | Line Profile Shape FWHM [km/s] (4) | cτBLR Hβ [ld] (5) | Fvar Hβ (6) | λLλ (5100) [1044 erg/s] (7) | MBH [M⊙] (8) | References (9) |
|---|---|---|---|---|---|---|---|---|
| NGC 5548 | 0.0172 | Sy 1.0–1.8 | Strong shoulders | 0.33 | 0.40 ± 0.12 | 2.1 × 109 | Shapovalova et al., ; Kovačević et al., | |
| (1996–2002) | 6,300 | |||||||
| NGC 4151 | 0.0033 | Sy 1.5–1.8 | Absorption,bumps | 0.42 | 0.05 ± 0.03 | 1.6 × 108 | Shapovalova et al., , ; Rakić et al., | |
| (1996–2006) | 6,110 | |||||||
| 3C390.3 | 0.0561 | Radio loud | Double-peaked | 0.38 | 0.90 ± 0.42 | 2.1 × 109 | Shapovalova et al., , ; Jovanović et al., ; Popović et al., | |
| (1995–2007) | 11,900 | |||||||
| Ark 564 | 0.0247 | NLSy1 | Strong FeII | 0.07 | 0.36 ± 0.04 | 1.0 × 106 | Shapovalova et al., | |
| (1999–2010) | 960 | |||||||
| Arp 102B | 0.0242 | LINER | Double-peaked | 0.21 | 0.11 ± 0.01 | 1.1 × 108 | Shapovalova et al., ; Popović et al., | |
| (1987–2013) | 15,900 | |||||||
| NGC 7469 | 0.0163 | Sy 1.0 | Narrower | 0.23 | 0.52 ± 0.0.08 | 1.1 × 107 | Shapovalova et al., | |
| (1996–2015) | 2,000 | |||||||
| E1821+643 | 0.297 | Quasar | Red-asymmetry | 0.07 | 104.4 ± 19.9 | 2.6 × 109 | Shapovalova et al., ; Kovačević et al., | |
| (1990–2014) | 5,610 |
Spectral characteristics of the sample.
The columns are: (1) object name and monitoring period, (2) redshift, (3) AGN type, (4) shape of line profile and FWHM of the mean Hβ, (5) BLR radius RBLR = cτBLR of Hβ, (6) variability parameter Fvar of Hβ, (7) mean continuum luminosity at 5100Å, (8) the mass of the SMBH MBH, and (9) references.
3. Results
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
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, ). 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
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., ). 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., ).
Figure 3
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. () 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. (), 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
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. (
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.,
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.,
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.,
We summarize here the main conclusions of this study:
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;
the luminosity of Hβ is correlated with the continuum luminosity for all AGN in our sample, and the intrinsic Baldwin effect is present;
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.
Statements
Author contributions
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
Conflict of interest
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.
Footnotes
1.^For Fe II modeling and subtraction see the page of the Serbian Virtual Observatory, http://servo.aob.rs.
2.^Since the intrinsic Baldwin effect is detected in a single object, we can use a flux instead of luminosity.
References
1
AlexanderT. (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 MaozD.SternbergA.LeibowitzE. M., (Dordrecht: Springer), 163.
2
BaldwinJ. A. (1977). Luminosity indicators in the spectra of quasi-stellar objects. Astrophys. J.214, 679–684. 10.1086/155294
3
BentzM. C.DenneyK. D.GrierC. J.BarthA. J.PetersonB. M.VestergaardM.et al. (2013). The low-luminosity end of the radius-luminosity relationship for active galactic nuclei. Astrophys. J.767:149. 10.1088/0004-637X/767/2/149
4
BentzM. C.PetersonB. M.PoggeR. W.VestergaardM.OnkenC. A. (2006). The radius-luminosity relationship for active galactic nuclei: the effect of host-galaxy starlight on luminosity measurements. Astrophys. J.644, 133–142. 10.1086/503537
5
BlandfordR. D.McKeeC. F. (1982). Reverberation mapping of the emission line regions of seyfert galaxies and quasars. Astrophys. J.255, 419–439. 10.1086/159843
6
BonE.JovanovićP.MarzianiP.ShapovalovaA. I.BonN.Borka JovanovićV.et al. (2012). The first spectroscopically resolved sub-parsec orbit of a supermassive binary black hole. Astrophys. J.759:118. 10.1088/0004-637X/759/2/118
7
BonE.ZuckerS.NetzerH.MarzianiP.BonN.JovanovićP.et al. (2016). Evidence for periodicity in 43 year-long monitoring of NGC 5548. Astrophys. J. Suppl. Ser.225:29. 10.3847/0067-0049/225/2/29
8
GaskellC. 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. 10.1086/165986
9
GaskellC. M.SparkeL. S. (1986). Line variations in quasars and Seyfert galaxies. Astrophys. J.305, 175–186. 10.1086/164238
10
IlićD.PopovićL. Č.ShapovalovaA. I.BurenkovA. N.ChavushyanV. H.KovačevićA. (2015). Line shape variability in a sample of AGN with broad lines. J. Astrophys. Astron.36, 433–445. 10.1007/s12036-015-9360-y
11
JovanovićP.PopovićL. Č.StalevskiM.ShapovalovaA. 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. 10.1088/0004-637X/718/1/168
12
KovačevićJ.PopovićL. Č.DimitrijevićM. S. (2010). Analysis of optical Fe II emission in a sample of active galactic nucleus spectra. Astrophys. J. Suppl.189, 15–36. 10.1088/0067-0049/189/1/15
13
KovačevićA.PopovićL. Č.ShapovalovaA. I.IlićD. (2017). Periodicity in the continua and broad line curves of a quasar E1821+643. Astrophys. Space Sci.362:31. 10.1007/s10509-017-3009-z
14
KovačevićA.PopovićL. Č.ShapovalovaA. I.IlićD.BurenkovA. N.ChavushyanV. 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. 10.1016/j.asr.2014.06.025
15
LiY.-R.WangJ.-M.HoL. C.LuK.-X.QiuJ.DuP.et al. (2016). Spectroscopic indication of a centi-parsec supermassive black hole binary in the galactic center of ngc 5548. Astrophys. J.822:4. 10.3847/0004-637X/822/1/4
16
LyutyiV. M.CherepashchukA. M. (1972). Narrow-band photoelectric observations of Hα-line variability in the nuclei of Seyfert galaxies NGC 4151, 3516, 1068. Astronomicheskij Tsirkulyar688.
17
O'BrienP. T.DietrichM.LeighlyK.AlloinD.ClavelJ.CrenshawD. 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. 10.1086/306464
18
OnkenC. A.FerrareseL.MerrittD.PetersonB. M.PoggeR. W.VestergaardM.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. 10.1086/424655
19
OsterbrockD. E.FerlandG. J. (2006). Astrophysics of Gaseous Nebulae and Active Galactic Nuclei. Sausalito, CA: University Science Books.
20
PetersonB. M.HorneK. (2006). Reverberation mapping of active galactic nuclei, in Planets to Cosmology, eds LivioM.CasertanoS. (Cambridge, UK: Cambridge University Press), 89.
21
PoggeR. W.PetersonB. M. (1992). The intrinsic nature of the Baldwin effect. Astronom. J.103, 1084–1088. 10.1086/116127
22
PopovićL. Č.MediavillaE.BonE.IlićD. (2004). Contribution of the Disk Emission to the Broad Emission Lines in AGNs: Two-component model. Astron. Astrophys.423, 909–918. 10.1051/0004-6361:20034431
23
PopovićL. Č.ShapovalovaA. I.IlićD.BurenkovA. N.ChavushyanV. H.KollatschnyW.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. 10.1051/0004-6361/201423555
24
PopovićL. Č.ShapovalovaA. I.IlićD.KovačevićA.KollatschnyW.BurenkovA. 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. 10.1051/0004-6361/201016317
25
RakićN.La MuraG.IlićD.ShapovalovaA. I.KollatschnyW.RafanelliP.et al. (2017). The intrinsic Baldwin effect in broad Balmer lines of six long-term monitored AGNs. Astron. Astrophys.603:A49. 10.1051/0004-6361/201630085
26
RobinsonA.YoungS.AxonD. J.KharbP.SmithJ. E. (2010). Spectropolarimetric evidence for a kicked supermassive black hole in the Quasar E1821+643. Astrophys. J. lett.717, L122–L126. 10.1088/2041-8205/717/2/L122
27
ShapovalovaA. I.BurenkovA. N.CarrascoL.ChavushyanV. H.DoroshenkoV. T.DumontA. 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. 10.1051/0004-6361:20011011
28
ShapovalovaA. I.DoroshenkoV. T.BochkarevN. G.BurenkovA. N.CarrascoL.ChavushyanV. H.et al. (2004). Profile variability of the Hα and Hβ broad emission lines in NGC 5548. Astron. Astrophys.422, 925–940. 10.1051/0004-6361:20035652
29
ShapovalovaA. I.PopovićL. Č.BurenkovA. N.ChavushyanV. H.IlićD.KollatschnyW.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. 10.1051/0004-6361/201014118
30
ShapovalovaA. I.PopovićL. Č.BurenkovA. N.ChavushyanV. H.IlićD.KollatschnyW.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. 10.1051/0004-6361/201321781
31
ShapovalovaA. I.PopovićL. Č.BurenkovA. N.ChavushyanV. 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. 10.1051/0004-6361/200912311
32
ShapovalovaA. I.PopovićL. Č.BurenkovA. N.ChavushyanV. 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. 10.1088/0067-0049/202/1/10
33
ShapovalovaA. I.PopovićL. Č.ChavushyanV. H.AfanasievV. 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. 10.1093/mnras/stx025
34
ShapovalovaA. I.PopovićL. Č.ChavushyanV. H.BurenkovA. N.IlićD.KollatschnyW.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. 10.3847/0067-0049/222/2/25
35
ShapovalovaA. I.PopovićL. Č.CollinS.BurenkovA. N.ChavushyanV. H.BochkarevN. 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. 10.1051/0004-6361:20079111
36
WrightE. L. (2006). A cosmology calculator for the World Wide Web. Astronom. Soc. Pacific118, 1711–1715. 10.1086/510102
Summary
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
Volume
4 - 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
Updates

Check for updates
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ć dilic@math.rs
This article was submitted to Milky Way and Galaxies, a section of the journal Frontiers in Astronomy and Space Sciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.