- 1Osservatorio Astronomico di Padova (INAF), Padua, Italy
- 2Istituto Nazionale di Fisica Nucleare, Sezione di Padova, Padua, Italy
- 3Osservatorio Astronomico di Brera (INAF), Merate, Italy
- 4Università degli Studi dell'Insubria, Varese, Italy
- 5Instituto de Astrofísica de Canarias, Santa Cruz de Tenerife, Spain
- 6Departamento de Astrofsica, Universidad de La Laguna, San Cristóbal de La Laguna, Spain
The Fermi Gamma-ray Observatory discovered about a thousand extragalactic sources emitting energy from 100 MeV to 100 GeV. The majority of these sources belong to the class of blazars characterized by a quasi-featureless optical spectrum (BL Lac Objects). This hampers the determination of their redshift and therefore hinders the characterization of this class of objects. To investigate the nature of these sources and to determine their redshift, we are carrying out an extensive campaign using the 10 m Gran Telescopio Canarias to obtain high S/N ratio optical spectra. These observations allow us to confirm the blazar nature of the targets, to find new redshifts or to set stringent limits on the redshift based on the minimum equivalent width of specific absorption features that can be measured in the spectrum and are expected from their host galaxy, assuming it is a massive elliptical galaxy. These results are of importance for the multi-frequencies emission models of the blazars, to test their extreme physics, to shed light on their cosmic evolution and abundance in the far Universe. These gamma emitters are also of great importance for the characterization of the extragalactic background light through the absorption by the IR-optical background photons.
1. Introduction
A blazar is a jetted active galactic nucleus (AGN) with the relativistic jet that points along the line of sight of the observer. These kind of objects are bright emitters at all frequencies (from radio to TeV regime), are characterized by high variability at all frequencies and large polarization, and are often dominated by the γ-ray emission especially during the outbursts.
Their spectral energy distribution (SED) shows the typical doubled-humped structure with two broad peaks: the first bump is located at low energies, typically in the infrared to X-ray region, and is interpreted as due to synchrotron emission produced by electrons of the jet spiraling along the lines of force of the magnetic field, instead the second peak is placed at higher frequency, between the X-ray and the MeV-TeV energies, and as suggested in most leptonic models, can be due to Compton scattering of the same electrons (e.g., Maraschi et al., 1992; Dermer and Schlickeiser, 1993; Ghisellini and Tavecchio, 2009).
Blazars are commonly classified in two categories, BL Lac Objects (BLLs) and the Flat Spectrum Radio Quasars (FSRQs), and this classification depends on the strength of their broad emission lines respect to the continuum. A more physical distinction is based on to the comparison between the luminosity of the broad line region (BLR) and the Eddington luminosity (e.g., Ghisellini et al., 2017): FSRQs have radially efficient accretion disk, while the BL Lac objects are not able to photoionise gas of the BLR, causing the absence of features in the majority of their optical spectra. This classification needs to know the mass of the accreting black hole and the redshift of the source, which for broad emission line AGNs can be determined by spectroscopy, while for the BLLs is arduous due to the faintness or absence of the emission/absorption lines, showing a completely featureless optical spectra (Falomo et al., 2014).
The advent of the Fermi γ-ray observatory, starting observations in 2008 (Atwood et al., 2009), with its systematic scanning of the entire sky every 3 h at the high energy band (HE; > 20 MeV), has revolutionized the blazar studies, previously performed with radio and X-ray surveys. The extragalactic γ-ray sky is dominated by blazars (Acero et al., 2015) and in the third AGN Fermi/LAT catalogue (Ackermann et al., 2015, 3LAC), 1738 blazars are reported, compared with the 3033 γ-ray detected sources:662 are classified as BLL, 491 as FSRQ, while the remaining blazars are as of uncertain type. It is worth to note that a growing sub-sample of the GeV blazars are also emitters at the TeV energies (VHE; E > 100 GeV), detected by the Cherenkov telescopes as MAGIC, VERITAS and the HESS arrays, that can sample energies down to 30–50 GeV. The majority of them are BLLs (in the TeVcat1 there are 60 BLLs against 6 FSRQs), implying that the BLL class dominates the extragalactic TeV sky.
Although BLLs are the most numerous extragalactic class in the HE and VHE bands, for a large fraction of them the redshift is still unknown or highly uncertain, because contrarily to most AGNs, the BLLs are characterized by featureless (or quasi) featureless spectra. On the basis of the recent statistics, it was proposed that on average BLLs have lower redshift and smaller high energy γ-ray luminosity than FSRQs (Ghisellini et al., 2017). However, this could be due to a bias since the number of robustly detected high redshift BLLs is significantly limited due to the difficulty to measure their distance. Hence the determination of the blazar redshift is crucial to calculate their luminosity, to build and characterize realistic emission models and to allow us a sound comparison of the multi-frequency SEDs between the two blazar classes (see the blazar sequence, Fossati et al., 1998; Ghisellini et al., 2017).
The estimation of the BLL redshifts is also essential to determine the properties of the extragalactic background light (EBL). The BLL γ-rays of high energy can interact with the EBL infrared-optical photons and produce pairs e−/e+, resulting in an absorption in the GeV-TeV BLL spectrum starting at frequencies and with optical depth that depend on the redshift of the γ-ray source and is more pronounced in the 0.5 < z < 2 interval (Franceschini et al., 2008). At higher z the absorption due to pair production moves to Fermi energies, completely extinguishing the source in the VHE regime. Although a significant number of FSRQ detections, up to z > 4 already exist (Ackermann et al., 2017), at the TeV energies, due to their Compton inverse peak position, only a small fraction of them are detected. Therefore the identification of high redshift BLLs at these energies is particularly challenging in order to study the earliest EBL components due to the first-light sources (Population III stars, galaxies or quasars) in the universe (e.g., Franceschini and Rodighiero, 2017).
These considerations motivated us to carry out an extensive campaign at the 10.4 m Gran Telescopio Canarias aimed to obtain high signal-to-noise ratio (SNR) optical spectra and to estimate the redshift of BLLs. The results of this project are shown in Landoni et al. (2015) for S20954+65, Paiano et al. (2016) for the TeV BLL S20109+22, Falomo et al. (2017) for the blazar B0218+357, Paiano et al. (2017b) for a sample of 15 TeV BLL and 7 TeV candidates with unknown or uncertain redshift and in Paiano et al. (2017a) focused on 10 BLLs detected by Fermi satellite suggested to be at relatively high redshift by previous works. Moreover, till now, we observed 40 unassociated γ-ray objects detected by Fermi and candidate to be blazars (Paiano et al., submitted), 20 BLLs optically selected among the SDSS blazars (Landoni et al., submitted) and 10 hard Fermi sources.
2. Sample, Data Reduction and Analysis
This spectroscopic program involved a conspicuous sample of BLLs, for a total of ~100 objects and for which the selection followed different criteria (see details in Paiano et al., 2017a,b). All of these objects were also selected to have unknown or uncertain redshift for which conflicting estimates are published in the literature, mainly due to spectra with low S/N.
The observations were gathered at the GTC using the medium resolution spectrograph OSIRIS (Cepa et al., 2003), configured with the grisms R1000B and R1000R and covering the spectral range 4100-10000 Å. Details about the observational strategy, the data reduction and analysis procedure are reported in Paiano et al. (2017a,b), where all spectra, corrected for atmospheric extinction and flux calibrated, are reported. They are also reported in the website http://www.oapd.inaf.it/zbllac/. On average, the S/N ranges from 100 to 200 at 4500 and 8000 Å, to a peak of 300 at ~6000 Å.
3. Results
Our spectra allow us to confirm the blazar nature of the observed targets and, on the basis on the spectroscopic properties, they exhibit four characteristic features: (1) weak emission lines typical of low density gas; (2) spectral lines of stars from the host galaxy; (3) intervening absorption systems due to cold gas; and (4) a pure featureless continuum. While the first three types can co-exist and from them a redshift can be derived when emission and/or absorption lines are identified, in the latter case, thank to high SNR achieved, we can set a stringent redshift lower limit on the basis of the minimum Equivalent Width (EW) of absorption lines expected from the starlight emission of the blazar host galaxy.
• Spectra with emission lines characteristic of low-density gas
Although the BLL optical spectra are mainly characterized by a featureless continuum, very weak emission lines can be seen, with an intrinsic luminosity lower than those observed in quasars. Owing to the weakness of these lines, their detectability depends on the state of the source, with faint state of the source that favors the detection of intrinsic emission lines, and on the quality of the observations, especially in terms of S/N ratio.
In particular from our spectra, we detect weak emission lines due to [O III] (5007 Å) in the spectrum of 1ES1215+303, W-Comae, MS1221.8+2452 and PKS1424+240, and the [O II] (3727 Å) emission in 1ES0033+595, 1ES1215+303, and PKS1424+240. For the first time, we are able to determine the spectroscopic redshift of z = 0.467 and z = 0.6047 for 1ES0033+595 and PKS1424+240 respectively (in Figure 1 the spectrum of PKS1424+240, one of the farthest TeV BLLs).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/309154/fspas-04-00045-HTML/image_m/fspas-04-00045-g001.jpg)
Figure 1. Spectrum of the TeV source PKS 1424+240 (z = 0.6047) obtained at GTC. Top: Flux-calibrated and dereddered spectrum, and normalized spectrum. The main telluric bands are indicated by ⊕, the absorption features from interstellar medium of our galaxies are labeled as IS (Inter-Stellar). Bottom: Close-up of the normalized spectrum around the found emission spectral lines, marked by line identification.
• Spectra with absorption lines due to the stellar population of the host galaxy
BLLs are located in the nuclei of giant/massive elliptical galaxies with a prominent spheroidal component and with average luminosity in the R band of MR ~ 23 (Falomo et al., 2014, and references therein). Their stellar population is composed by old stars and therefore the main observable absorption lines are Ca II (3934, 3967 Å), G-band (4304 Å), Hβ (4861 Å), Mg I (5175 Å) and Na I (5875 Å). These absorption features can be detected over the non-thermal component and this depends on the signal-to-noise ratio and the spectral resolution of the spectrum, and finally on the relative power of the non-thermal nucleus and the host galaxy. In very few cases, high quality spectroscopy with adequate high spatial resolution can probe star formation of the host galaxy and allow us to detect narrow emission lines from it.
In our sample, we were able to detect host galaxy absorption lines for W Comae (z = 0.102), MS 1221.8+2452 (z = 0.218), 3FGL J0505.5+0416 (z = 0.423), and 3FGL J0814.5+2943 (z = 0.703, see Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/309154/fspas-04-00045-HTML/image_m/fspas-04-00045-g002.jpg)
Figure 2. Spectrum with S/N ~ 160 of the 3FGL BLL J0814.5+2943. In the close-up image, the absorption doublet at 6699–6759 Å identified as Ca II (3934, 3967 Å) at z = 0.703 from the the host galaxy.
• Spectra with intervening absorption lines
As for the high redshift quasars, gas along the BLL direction can produce systems of absorption at redshifts lower than the target redshift. Given that BLLs show quasi featureless continuous spectrum, they are ideal sources for studying absorptions from these intervening systems (Landoni et al., 2012, and references therein). The redshift of the intervening systems yields a lower limit to the BLL redshift, while an upper limit can be set from the expected distribution of the absorbers (Zhu and Ménard, 2013).
For six of our targets, 3FGL J0008.0+4713, BZB J1243+3627, and 3FGL J1450.9+5200 with an uncertain redshift in literature, and BZB J1540+8155, 3FGL J1107.5+0222, and BZB J2323+4210 with unknown redshift, we detected intervening absorption systems that allow us to set spectroscopic redshift lower limits. In Figure 3 two examples of spectra of our sample that exhibit Lyα forest, C IV (1548 Å) and Mg II (2800 Å) absorption systems.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/309154/fspas-04-00045-HTML/image_m/fspas-04-00045-g003.jpg)
Figure 3. Spectrum of the 3FGL BLLs J0008.0+4713, that shows an intervening doublet absorption system at 7440 Å identified as Mg II (2800 Å) at z = 1.659, and the spectrum of the 3FGL J1450.9+5200 that shows the absorption line at 5372 Å due to C IV (1548 Å) intervening gas and a strong Lyα (1216 Å) absorption at the same redshift (z = 2.470). A second C IV intervening system at z = 2.312 is associated to the absorption feature detected at 5127 Å.
• Featureless spectra
In spite of the high quality of the optical spectra for 18 observed targets the spectrum appears very featureless. This occurs when the emission from the underlying nebulosity of the host galaxy is over-shined by the non-thermal continuum.
As example, Figure 4 shows the featureless spectra of two TeV BLLs RGB J0136+391 and 3C 66A. For the latter, a previous redshift of 0.444 was proposed in literature (Miller et al., 1978) and reported in the NED. On the basis of our high S/N ratio and featureless spectrum, we cannot support this value and thus the redshift of 3C 66A is still unknown and the emission models of the source assuming the wrong redshift should be revised.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/309154/fspas-04-00045-HTML/image_m/fspas-04-00045-g004.jpg)
Figure 4. Featureless spectra of the TeV BLL RGB J0136+391 and 3C 66A. No emission/absorption lines are identified, hence for these sources the redshift is still unknown.
In the case of no detection of spectral lines, following the scheme outlined in Appendix A of Paiano et al. (2017b), we can set redshift lower limit using the minimum Equivalent Width method, based on the assumption that the BL Lac objects are hosted in massive elliptical galaxies and that hence the detection of their stellar absorption features depends on the SNR of their optical spectra, on the spectral resolution, and on the flux ratio between the nucleus and the host galaxy. The redshift lower limit found for 18 objects with featureless spectra are summarized in Table 3 of Paiano et al. (2017b) and Paiano et al. (2017a).
4. Conclusion
We secured high quality and S/N spectra for 32 γ-ray BL Lac objects, selected for being TeV sources, TeV candidates or high redshift Fermi objects, with unknown redshift or with conflicting redshift values in literature. We determined the redshift for 8 objects and spectroscopic redshift lower limit for 6 sources, including three of the farthest BLLs detected in the GeV and TeV regime, but for the remaining sources, in spite of the very high S/N and of the improvement of spectral resolution, their spectrum is featureless.
Although these observations represent the state of the art of the observing facilities (telescopes with large aperture and modern instrumentation) and the capabilities to perform the spectroscopy of BLLs, the redshift determination of this class of objects remains still rather arduous. For this reason, the detection of the spectral features in their spectra requires observations with very high SNR and an adequate resolution that seems feasible only with the next generation of extremely large (30–50 m class) telescopes, such as E-ELT (the European Extremely Large Telescope). An idea about the performances of these new class of telescopes is given in Landoni et al. (2014) where simulations of BLL spectra are made using an X-Shooter-like instrument (with an spectral resolution ~3000) coupled to the E-ELT and including a moderate adaptive optics module. With these new facilities, it should be possible to measure the redshift of BLLs having extreme nucleus-to-host galaxy ratios between ~300 to ~2,500, depending on the redshift and the target brightness.
Author Contributions
SP is the corresponding author of this contribution and presented the corresponding talk at the conference “QSO at all epochs 2017.” All authors contributed to data analysis, text revision and improvement.
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 handling Editor declared a shared affiliation, though no other collaboration, with the authors SP and RF.
Acknowledgments
The authors of this work thank the organizers of the “Quasars at All Cosmic Epoch” conference.
Footnotes
References
Acero, F., Ackermann, M., Ajello, M., Albert, A., Atwood, W. B., Axelsson, M., et al. (2015). Fermi large area telescope third source catalog - Fermi-LAT collaboration. Astrophys. J. Suppl. 218:23. doi: 10.1088/0067-0049/218/2/23
Ackermann, M., Ajello, M., Atwood, W. B., Baldini, L., Ballet, J., Barbiellini, G., et al. (2015). The third catalog of active galactic nuclei detected by the fermi large area telescope. Astrophys. J. 810:4. doi: 10.1088/0004-637X/810/1/14
Ackermann, M., Ajello, M., Baldini, L., Ballet, J., and Barbiellini, G. (2017). Gamma-ray blazars within the first 2 billion years. Astrophys. J. Lett. 837:L5. doi: 10.3847/2041-8213/aa5fff
Atwood, W. B., Abdo, A. A., Ackermann, M., Althouse, W., Anderson, B., Axelsson, M., et al. (2009). The large area telescope on the fermi gamma-ray space telescope mission. Astrophys. J. 697, 1071–1102. doi: 10.1088/0004-637X/697/2/1071
Cepa, J., Aguiar-Gonzalez, M., Bland-Hawthorn, J., Castaneda, H., Cobos, F. J., Correa, S., et al. (2003). OSIRIS tunable imager and spectrograph for the GTC. Instrument status. Proc. SPIE 4841, 1739–1749. doi: 10.1117/12.460913
Dermer, C. D., and Schlickeiser, R. (1993). Model for the high-energy emission from blazars. Astrophys. J. 416:458. doi: 10.1086/173251
Falomo, R., Pian, E., and Treves, A. (2014). An optical view of BL Lacertae objects. Astron. Astrophys. Rev. 22:73. doi: 10.1007/s00159-014-0073-z
Falomo, R., Treves, A., Scarpa, R., Paiano, S., and Landoni, M. (2017). On the lensed blazar B0218+357. Month. Notices R. Astron. Soc. 470, 2814–2821. doi: 10.1093/mnras/stx1411
Fossati, G., Maraschi, L., Celotti, A., Comastri, A., and Ghisellini, G. (1998). A unifying view of the spectral energy distributions of blazars. Month. Notices R. Astron. Soc. 299, 433–448. doi: 10.1046/j.1365-8711.1998.01828.x
Franceschini, A., and Rodighiero, G. (2017). The extragalactic background light revisited and the cosmic photon-photon opacity. Astron. Astrophys. 603:A34. doi: 10.1051/0004-6361/201629684
Franceschini, A., Rodighiero, G., and Vaccari, M. (2008). Extragalactic optical-infrared background radiation, its time evolution and the cosmic photon-photon opacity. Astron. Astrophys. 487, 837–852. doi: 10.1051/0004-6361:200809691
Ghisellini, G., Righi, C., Costamante, L., and Tavecchio, F. (2017). The Fermi blazar sequence. Month. Notices R. Astron. Soc. 469, 255–266. doi: 10.1093/mnras/stx806
Ghisellini, G., and Tavecchio, F. (2009). Canonical high-power blazars. Month. Notices R. Astron. Soc. 397, 985–1002. doi: 10.1111/j.1365-2966.2009.15007.x
Landoni, M., Falomo, R., Treves, A., and Sbarufatti, B. (2014). Spectroscopy of BL Lacertae objects of extraordinary luminosity. Astron. Astrophys. 570:A126. doi: 10.1051/0004-6361/201424232
Landoni, M., Falomo, R., Treves, A., Sbarufatti, B., Decarli, R., Tavecchio, F., et al. (2012). On the redshift of the bright BL Lacertae object PKS 0048-097. Astron. Astrophys. 543:A116. doi: 10.1051/0004-6361/201219114
Landoni, M., Falomo, R., Treves, A., Scarpa, R., and Reverte Payá, D. (2015). What Is the Redshift of the Gamma-ray BL Lac Source S4 0954+65? Astron. J. 150:181. doi: 10.1088/0004-6256/150/6/181
Maraschi, L., Ghisellini, G., and Celotti, A. (1992). A jet model for the gamma-ray emitting blazar 3C 279. Astrophys. J. Lett. 397, L5–L9. doi: 10.1086/186531
Miller, J. S., French, H. B., and Hawley, S. A. (1978). “Optical spectra of BL Lacertae objects,” in BL Lac Objects, ed A. M. Wolfe (Pittsburgh, PA: University of Pittsburgh), 176–187.
Paiano, S., Landoni, M., Falomo, R., Scarpa, R., and Treves, A. (2016). On the redshift of the very high-energy gamma-ray BL Lac object S2 0109+22. Month. Notices R. Astron. Soc. 458, 2836–2839. doi: 10.1093/mnras/stw472
Paiano, S., Landoni, M., Falomo, R., Treves, A., and Scarpa, R. (2017a). Spectroscopy of 10 γ-ray BL Lac objects at high redshift. Astrophys. J. 844:120. doi: 10.3847/1538-4357/aa7aac
Paiano, S., Landoni, M., Falomo, R., Treves, A., Scarpa, R., and Righi, C. (2017b). On the Redshift of TeV BL Lac Objects. Astrophys. J. 837:144. doi: 10.3847/1538-4357/837/2/144
Keywords: blazars, BL Lac objects, optical spectroscopy, redshift, γ-ray astronomy
Citation: Paiano S, Falomo R, Landoni M, Treves A and Scarpa R (2017) An Optical View of Extragalactic γ-Ray Emitters. Front. Astron. Space Sci. 4:45. doi: 10.3389/fspas.2017.00045
Received: 06 September 2017; Accepted: 31 October 2017;
Published: 23 November 2017.
Edited by:
Paola Marziani, Osservatorio Astronomico di Padova (INAF), ItalyReviewed by:
Brian Punsly, International Center for Relativistic Astrophysics, ItalyDawei Xu, National Astronomical Observatories (CAS), China
Copyright © 2017 Paiano, Falomo, Landoni, Treves and Scarpa. 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: Simona Paiano, simona.paiano@oapd.inaf.it