Issue |
A&A
Volume 638, June 2020
|
|
---|---|---|
Article Number | A128 | |
Number of page(s) | 13 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/201936928 | |
Published online | 24 June 2020 |
The γ-ray sky seen at X-ray energies
II. The Swift hunt of Fermi BL Lac objects among unidentified γ-ray sources
1
Dipartimento di Fisica, Università degli Studi di Torino, Via Pietro Giuria 1, 10125 Turin, Italy
e-mail: ejmarchesini@gmail.com
2
Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, Paseo del Bosque, B1900FWA La Plata, Argentina
3
INFN – Istituto Nazionale di Fisica Nucleare, Sezione di Torino, Via Pietro Giuria 1, 10125 Turin, Italy
4
Instituto de Astrofísica de La Plata, CONICET–UNLP, CCT La Plata, Paseo del Bosque, B1900FWA La Plata, Argentina
5
INAF – Osservatorio di Astrofisica e Scienza dello Spazio, Via Gobetti 93/3, 40129 Bologna, Italy
6
Departamento de Ciencias Fisicas, Universidad Andres Bello, Fernandez Concha 700, Las Condes, Santiago, Chile
7
Center for Astrophysics | Harvard & Smithsonian, 60 Garden St, Cambridge, MA 02138, USA
Received:
15
October
2019
Accepted:
9
April
2020
Context. Nearly 50% of all sources detected by the Fermi Large Area Telescope are classified as blazars or blazar candidates, one of the most elusive classes of active galaxies. Additional blazars can also be hidden within the sample of unidentified or unassociated γ-ray sources (UGSs) that constitute about one-third of all gamma-ray sources detected to date. We recently confirmed that the large majority of Fermi blazars of the BL Lac subclass have an X-ray counterpart.
Aims. Using the X-ray properties of a BL Lac training set and combining these with archival multifrequency information, we aim to search for UGSs that could have a BL Lac source within their γ-ray positional uncertainty regions.
Methods. We reduced and analyzed the Swift X-ray observations of a selected sample of 327 UGSs. We then compared the X-ray fluxes and hardness ratios of all sources detected in the pointed fields with those of known Fermi BL Lacs.
Results. We find at least one X-ray source, lying within the γ-ray positional uncertainty at 95% confidence level, for 223 UGSs and a total of 464 X-ray sources in all fields analyzed. The X-ray properties of a large fraction of them, eventually combined with radio, infrared, and optical information, exhibit BL Lac multi-frequency behavior, thus allowing us to select high-confidence BL Lac candidates; some of them were recently observed during our optical spectroscopic campaign which confirmed their nature.
Conclusions. We find that out of 50 X-ray sources that were confirmed as BL Lacs through optical spectroscopy, 12 do not show canonical mid-infrared or radio BL Lac properties. This indicates that the selection of X-ray BL Lac candidates is a strong method to find new counterparts within Fermi UGSs. Finally, we pinpoint a sample of 32 Swift/XRT candidate counterparts to Fermi UGSs that are most likely BL Lac objects.
Key words: galaxies: active / galaxies: nuclei / galaxies: jets / BL Lacertae objects: general / X-rays: galaxies / gamma rays: galaxies
© ESO 2020
1. Introduction
Since its launch in 2008, the Large Area Telescope (LAT, Atwood et al. 2009) onboard the Fermi satellite, thanks to its superior angular resolution and larger collecting area with respect to previous γ-ray telescopes, has catalogued ∼3000 sources detected to date above ∼4σ significance and listed in the Third Fermi Large Area Telescope Source Catalog (3FGL, Acero et al. 2015).
Recently, it has been announced that this number will increase to more than 5000 in the upcoming release of the Fourth Fermi Large Area Telescope Source Catalog (4FGL, Abdollahi et al. 2020). Nearly one-third of all sources associated with a low-energy counterpart in the 3FGL belong to the type of active galaxies known as blazars, as was the case with the previous Fermi catalogs: 1FGL (Abdo et al. 2010) and 2FGL (Nolan et al. 2012).
Blazars are a peculiar class of active galactic nuclei characterized by emission arising from a relativistic jet oriented at small angles (less than a few degrees; see e.g., Pushkarev et al. 2009; Lister et al. 2013) with respect to the line of sight, which in most cases overwhelms the radiation of its host galaxy (Blandford & Rees 1978). They are divided into two main classes: flat-spectrum radio quasars with broad emission lines present in their optical spectra, and BL Lac sources, which show featureless optical spectra (Stickel et al. 1991; Landoni et al. 2014).
Two different subclasses for BL Lacs were defined on the basis of the position of the spectral energy distribution (SED) peak for their first component: low-frequency and high-frequency peaked BL Lacs (i.e., HBLs and LBLs; Padovani & Giommi 1995). This criterion is also equivalent to a certain threshold in the ratio between X-ray and radio flux (see also Maselli et al. 2010, and references therein).
Another ∼20% of the sources listed in the 3FGL catalog belong to the Blazar Candidate of Uncertain Type (BCU) class, including γ-ray celestial objects having counterparts that show blazar-like characteristics at other wavelengths but still need confirmation of their nature, mainly via optical spectroscopic observations (Álvarez Crespo et al. 2016a; Peña-Herazo et al. 2019; D’Abrusco et al. 2019). Furthermore, an additional ∼10% of the sources listed in 3FGL have Galactic origin, mainly pulsars, pulsar wind nebulae, and supernova remnants (see e.g., Massaro et al. 2015a,b, and references therein).
The remaining one-third of the 3FGL sources are still unassociated/unidentified. These sources, which also amount to nearly one-fourth in the preliminary release of the 4FGL, are labeled as Unidentified/Unassociated Gamma-ray Sources (UGSs), lacking an assigned low-energy counterpart (Massaro et al. 2013a). Searching for low-energy sources potentially associated with the remaining UGSs is still a challenging task that requires extensive multifrequency follow-up campaigns (e.g., Acero et al. 2013; Marchesini et al. 2016; Lico et al. 2016) alongside statistical analyses (Ackermann et al. 2012a; Hassan et al. 2013; D’Abrusco et al. 2014, 2019).
Several methods have been developed to search for UGS counterparts. These are based, for example, on optical polarization (Blinov et al. 2018), optical spectroscopy (Sandrinelli et al. 2013; Paiano et al. 2017a), γ-ray spectral and variability characteristics (Ackermann et al. 2012b), or on machine learning algorithms (Doert & Errando 2014; Salvetti et al. 2017). There are also methods that rely on radio (Healey et al. 2007; Ghirlanda et al. 2010; Hovatta et al. 2012; Nori et al. 2014; Massaro et al. 2013a) or infrared (D’Abrusco et al. 2012, 2013) observations. In particular, the latter two methods are based on well-known connections between the blazar emission in these last two bands, and their γ-ray properties (see, e.g., Taylor et al. 2007; Mahony et al. 2010; Ackermann et al. 2011; Massaro & D’Abrusco 2016, and references therein).
Nevertheless, optical spectroscopy is still the only method that can determine the real nature of the associated counterpart, and spectroscopic campaigns confirm that most of the UGSs are indeed blazars, in particular of the BL Lac kind (Paggi et al. 2014; Landoni et al. 2015, 2018; Ricci et al. 2015; Falomo et al. 2017; Kaur et al. 2019; Franceschini et al. 2019).
X-ray follow up observations have also been widely used as a tool to search for X-ray counterparts that could be blazars for the UGSs (e.g., Stephen et al. 2010; Takahashi et al. 2012; Maselli et al. 2013; Paggi et al. 2013; Masetti et al. 2013; Massaro et al. 2014; Landi et al. 2015; Marchesini et al. 2016; Paiano et al. 2017a). These X-ray observations were carried out even if a connection between X-ray and γ-ray emission in blazars was not established as in other energy ranges.
We recently analyzed a uniform sample of Fermi BZBs defined by Roma-BZCAT (Massaro et al. 2015c) as BL Lac objects that have been confirmed as such through optical spectroscopy. We proved that, above the γ-ray flux (i.e., Fγ) threshold of the order of 10−11 erg cm−2 s−1, all known γ-ray BZBs have an X-ray counterpart detected by Swift/XRT (Burrows et al. 2005), with ∼5 ks exposure time and with signal to noise ratios (S/Ns) larger than 3 (Marchesini et al. 2019a, hereinafter Paper I). This X-ray-γ-ray connection, even if different from those observed at radio and mid-infrared frequencies, supports ongoing and future X-ray follow up observations of UGSs, aiming to find BZB-like counterparts1 (Stroh & Falcone 2013; Falcone et al. 2014). It is worth noting that follow-up spectroscopic observations of the X-ray source lying within the positional uncertainty region of each UGS are always necessary to confirm its nature and to potentially obtain its redshift (Massaro et al. 2015d, 2016; Paiano et al. 2017b, 2019).
In this work, we analyze soft X-ray (i.e., in the 0.5−10 keV energy range) observations of Fermi UGSs carried out thanks to the ongoing Swift/XRT follow-up campaign (Stroh & Falcone 2013; Falcone et al. 2014), searching for X-ray sources that show X-ray spectral behaviour similar to that of known and previously analyzed Fermi BZBs. We also aim to provide a catalog of X-ray sources that could be targeted by follow up spectroscopic observations to search for UGS counterparts.
This paper is organized as follows. In Sect. 2 we describe the sample-selection criteria, while in Sect. 3 we describe the data-reduction procedure. Section 4 is devoted to report our results. Finally, in Sect. 5 we summarize our main results and report our conclusions.
Throughout the paper, we adopt cgs units and a flat cosmology with H0 = 72 km s−1 Mpc−1, ΩΛ = 0.74, and Ωm = 0.26 (Dunkley et al. 2009). Spectral indices α were defined so that flux density Sν ∝ ν−α, considering α < 0.5 as flat spectra, especially at radio frequencies around 1.4 GHz. The Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010) magnitudes from the AllWISE catalog, in the 3.4 μm, 4.6 μm, and 12 μm nominal filters, are in the Vega system, and are not corrected for Galactic extinction since this correction is negligible for Galactic latitudes above and below 10° (D’Abrusco et al. 2013).
2. Sample selection
We selected all Fermi UGSs listed in the 3FGL (Acero et al. 2015), amounting to a total of 1010 sources. We then chose those that had at least one Swift/XRT observation performed in photon counting (PC) mode, lying within a circular region of 6 arcmin radius around UGSs positions, for a total of 706 Fermi sources. A radius of 6 arcmin was chosen because this is the average semi-major axis of the elliptical positional uncertainty region of UGSs at 95% level of confidence in the 3FGL, as performed in our previous analysis of Fermi BZBs reported in Paper I.
We then selected those with total exposure time in the range between 1 and 10 ks, in agreement with Paper I, in which we avoided both under- and over-exposed observations, because the latter are not snapshots. In any case, the aforementioned Swift X-ray campaign of UGSs is performed with a nominal 5 ks exposure time (Stroh & Falcone 2013), which is well sampled within our exposure range. Thus, the number of UGSs decreased to 636.
Finally, we discarded all UGSs lying within the Galactic plane (|b| < 10°). This reduces the contamination from Galactic sources. Our final sample consists of 327 Fermi UGSs, with an average exposure time of 4.2 ks. The flow chart shown in Fig. 1 summarizes all selection steps reported above.
![]() |
Fig. 1. Flow chart to highlight all steps followed to build our final UGS sample. |
3. Swift/XRT data reduction
3.1. Data processing
For all Swift/XRT observations analyzed here, we adopted the same data-reduction procedure as in Paper I, following Massaro et al. (2008, 2012a), Paggi et al. (2013), and references therein. Here we report only basic details.
All Swift/XRT observations of our sample of Fermi UGSs were reduced with standard procedures, obtaining clean event files with the use of the XRTPIPELINE task version 0.13.4, which is part of the Swift X-ray Telescope Data Analysis Software (XRTDAS, Capalbi et al. 2005). We used the High Energy Astrophysics Science Archive Research Center (HEASARC) calibration database (CALDB) version 1.0.2.
We excluded all time intervals with count rates exceeding 40 counts per second using the XSELECT task, and all time intervals during which the CCD temperature exceeds −50°C in regions located at the CCD edge, following recommendations from D’Elia et al. (2013). Our data-reduction procedure is similar to the one adopted by the Swift XRT Point Source catalog (1SXPS, Evans et al. 2014), with our results showing differences of only a few percent (see Paper I for details).
After extracting all clean event files, we merged exposures of the same source using XSELECT. The same was done with the corresponding exposure maps with the use of the XIMAGE software.
3.2. Source detection and photometry
We performed a first detection run over merged event files with the DET algorithm in XIMAGE. This way we obtained pixel positions for every detection with S/N larger than 3.
We then used the SOSTA task, available within the XIMAGE package, with the pixel positions obtained from DET. Since SOSTA takes into account the local background and foreground of each X-ray source to evaluate the source photometry, it achieves more precise results than the DET algorithm. This way we obtained counts and count rates for each source in the full 0.5−10 keV band, and also in the soft (0.5−2 keV) and hard (2−10 keV) bands. We chose a S/N minimum threshold of 3 to claim a source detection, as in Paper I.
We detected 464 X-ray sources within 223 UGS fields. No X-ray sources were found for the remaining 104 UGS, and were consequently excluded from the present analysis. We obtained positions, photon counts, and count rates corrected by the exposure map for every X-ray source. No X-ray sources detected in our sample show pile-up. We then derived the hardness ratio (HRX) for each X-ray source, and full band X-ray fluxes (FX). The hardness ratio was computed as (H−S)/(H+S), where H are counts in the hard band (2 to 10 keV) and S those in the soft band (0.5 to 2 keV), respectively. We remark as in our previous work that the exposure map does not change when splitting the data into bands, meaning that using counts or count rates to obtain HRX is equivalent because the exposure map does not vary between the two energy ranges.
To derive fluxes, we used PIMMS (Mukai 1993), assuming a power-law model with a photon index of 2.02 Galactic column density values were taken from the Kalberla et al. (2005). Applying the same criteria reported in Paper I, we excluded from our sample extended sources and spurious detections due to artifacts and bad pixels, leading to a final sample of 397 XRT sources detected lying within the fields of 223 UGSs. For each X-ray source with equatorial celestial right ascension HH:MM:SS.s and declination ±DD:MM:SS, we adopt a designation of the form SWXRT JHHMMSS.s±DDMMSS, following Paggi et al. (2013) as in D’Elia et al. (2013). In Fig. 2 we show the XRT merged image of the Fermi UGS 3FGL J0121.8−3917, for which we find one XRT counterpart lying within its positional uncertainty ellipse. In Table 1, we present results obtained from our X-ray data.
![]() |
Fig. 2. Full-band (0.5−10 keV) merged XRT image corresponding to the Fermi UGS 3FGL J0121.8−3917. The yellow dashed line indicates the 3FGL positional uncertainty ellipse at 95% level of confidence, while the white circles show the XRT sources detected with S/N greater than 3. There is only one XRT counterpart within the positional uncertainty ellipse of 3FGL J0121.8−3917. The image was smoothed with a Gaussian kernel radius of 5 pixels. |
First ten rows of the table summarizing the results of our UGS X-ray analysis.
4. Results
4.1. Crossmatches
Searching for BZB that could be potentially associated with selected UGSs, we first restricted the sample considering X-ray sources lying within the Fermi positional uncertainty ellipses drawn at 95% confidence level. This crossmatch was done with an XRT positional uncertainty circle with a radius of 5.6 arcsec, which is the maximum positional uncertainty reported for 97% of all 160250 sources in the 1SXPS catalog that lay outside the Galactic plane and were detected with S/N greater than 3 (Evans et al. 2014). This resulted in 197 X-ray sources within the positional uncertainty regions of 154 UGSs.
We found a single X-ray source lying within the positional uncertainty ellipse for 121 UGSs, while for the remaining 33 UGSs there were multiple X-ray counterparts. In particular, 27 out of these 33 UGSs had two X-ray potential counterparts, while 6 UGSs had 3 to 6 X-ray source counterparts within their γ-ray positional uncertainty.
In the present analysis, UGSs with more than two X-ray sources lying within the Fermi positional uncertainty regions were not investigated, and will be analyzed in a forthcoming paper together with those lying at Galactic latitudes (i.e., |b| < 10°). As an example, in Fig. 3 we show the XRT merged image of the Fermi UGS 3FGL J0809.3−0941, for which we find two X-ray counterparts lying within its positional uncertainty ellipse at 95% confidence level.
![]() |
Fig. 3. Full-band (0.5−10 keV) merged XRT image corresponding to the Fermi UGS 3FGL J0809.3−0941. The yellow dashed line indicates the 3FGL positional uncertainty ellipse at 95% confidence level, while the white circles show all XRT sources detected with S/N greater than 3. We find two X-ray sources lying within the Fermi ellipse, which we label both as XBCs; the remaining one is discarded from this analysis. The image was smoothed with a Gaussian kernel of radius 5 pixels. |
Our sample lists 148 Fermi UGSs having between one and two X-ray sources within their positional uncertainty ellipse at 95% confidence level, for a total of 175 X-ray sources. This sample is labeled as X-ray blazar candidates (XBCs).
On the other hand, there are also 69 Fermi UGSs with at least one X-ray source detected in the XRT field that does not lie within their positional uncertainty regions. This second sample, labeled as outliers (hereinafter marked as OUTs), lists 105 X-ray sources. We investigate this sample separately from the XBCs, since it is expected that ∼5% of the potential X-ray counterparts do not lie within the 3FGL positional uncertainty ellipse (Acero et al. 2015; Massaro et al. 2015a). The remaining OUTs are probably not associated to the Fermi sources, given the fact that we chose to use the positional uncertainty of the 3FGL at 95% confidence level.
In Fig. 4, we show the XRT merged image corresponding to the field of the Fermi UGS 3FGL J0216.0+0300, in which we find three OUTs. In Fig. 5, we show the distribution of the angular separation in arcminutes between the XRT and the 3FGL positions for both subsamples.
![]() |
Fig. 4. Full-band (0.5−10 keV) merged XRT image corresponding to the Fermi UGS 3FGL J0216.0+0300. The yellow dashed line indicates the 3FGL positional uncertainty ellipse at 95% confidence level, while the white circles show all XRT sources detected with S/N greater than 3. We find three OUTs in this field, since no X-ray source is detected within the Fermi ellipse. The image was smoothed with a Gaussian kernel of radius 5 pixels. |
![]() |
Fig. 5. Angular separation between Swift/XRT and 3FGL positions for XBCs (in purple) and OUTs (in black) in arcminutes. There is an overlap between XBCs and OUTs, implying the latter could be associated in future Fermi catalogs. |
It is worth noting that there are X-ray sources in the OUT sample lying at an angular separation comparable to that of XBCs. These could eventually become associated in future releases of the Fermi catalogs. Some BZBs could lie outside of the Fermi positional uncertainty regions in the current version of the catalog, but a subsequent refined analysis with larger exposure or improved source detection accuracy (Arsioli & Chang 2017) could eventually improve source localization (e.g., due to the improvements in the diffuse model of the gamma-ray background). This could provide future associations, as has occurred in previous versions of the Fermi catalogs (Abdo et al. 2010; Nolan et al. 2012; Acero et al. 2015).
For both subsamples, XBCs and OUTs, we looked for mid-infrared magnitudes in the AllWISE catalog (Wright et al. 2010), and for optical magnitudes in the Sloan Digital Sky Survey (SDSS, Ahn et al. 2012). We also crossmatched both subsamples with the National Radio Astronomy Observatory Very Large Array Sky Survey catalog (NVSS, Condon et al. 1998) and the Sydney University Molonglo Sky Survey (SUMSS, Mauch et al. 2003).
These crossmatches were done by combining 5.6 arcsec (which is the average positional uncertainty of all sources listed in the 1SXPS catalog Evans et al. 2014), with the positional uncertainty of the corresponding catalog. In the case of radio catalogs, we used 10.3 and 7.4 arcsec for the NVSS and SUMSS catalogs, respectively, as in D’Abrusco et al. (2014, 2019). For WISE we considered an uncertainty of 3.3 arcsec, following D’Abrusco et al. (2013), while for SDSS we considered an uncertainty of 1.8 arcsec, following Massaro et al. (2014).
In a total of nine cases between XBCs and OUTs, the crossmatch with WISE and SDSS resulted in two counterparts associated to the same XRT source. We chose to keep only the closest one. In total, we obtained 161, 80, and 49 WISE, radio, and SDSS unique counterparts for XBCs, and 100, 22, and 34 WISE, radio, and SDSS for OUTs, respectively. In particular, 76 and 21 XBCs and OUTs, respectively, have both a radio and a mid-infrared counterpart; and 25 and 5 have optical, mid-infrared, and radio counterparts, respectively. The results for all crossmatches are shown in Table 2, and in Fig. 6 we summarize the selection process for the two samples in a flow chart.
![]() |
Fig. 6. Flow chart to highlight all steps followed to build our final XBC and OUT samples. |
Results from the positional crossmatch of XBCs and OUTs with radio, mid-infrared, and optical catalogs.
4.2. X-ray selection of candidate γ-ray BL Lacs within the XBCs
We compared the X-ray properties of XBCs with that of BZBs analyzed in Paper I. Although the 3FGL average positional uncertainty of UGSs is of 6.2 ± 2.1 arcmin, our sample of XBCs show an average angular separation with the 3FGL position of 2.8 ± 1.3 arcmin, as can be seen in Fig. 5. In particular, 93 XBCs lie within less than 3 arcmin of the 3FGL positional centroid (28% of all UGSs in our sample). Half of all associated identified BZBs in the 3FGL catalog also lie within 3 arcmin of the 3FGL positional centroid, suggesting that a fraction of our XBCs could be associated in future releases of the Fermi catalogs as explained in Sect. 4.1.
In Fig. 7 we show HRX versus FX for all XBCs (in purple filled circles). We also plot these values for Fermi BZBs that we analyzed in Paper I as a comparison; these are also reported as background for comparison (i.e., marked as grey filled squares) in all figures shown in the present analysis.
![]() |
Fig. 7. HRX versus FX in the 0.5−10 keV band for all sources in the XBC sample. XBCs are plotted in empty purple circles, while Fermi BZBs are plotted in grey filled squares. When two XBCs lie within the positional uncertainty area of the same UGS, the most distant one is marked with filled black circles. |
X-ray Blazar Candidates are on average one order of magnitude fainter in X-rays (average flux FX = (3.3 ± 1.8)×10−13 erg cm−2 s−1) than BZBs (FX = (1.2 ± 0.9)×10−12 erg cm−2 s−1). This is expected, given the characteristics of each corresponding parent sample: UGSs are in general fainter than identified associated objects in the 3FGL. This is shown and discussed in more detail in Fig. 8.
![]() |
Fig. 8. γ-ray energy flux in the 100 MeV−100 GeV band vs. FX from the Fermi 3FGL catalog for all XBCs in our sample. The symbol coding is the same as in Fig. 7. The Fγ thresholds above which 100%, 98%, and 96% of the BZB sample would be detected in X-rays when their exposure times have been scaled down to 5 ks are plotted in solid, dashed, and dotted black lines, respectively. |
The sample of XBCs shows an average FX compatible at 1σ level with the FX of what we defined in Paper I as “background/foreground objects”, that is XRT counterparts not associated with BZBs. These show an average flux of FX = (1.0 ± 0.4)×10−13 erg cm−2 s−1.
On the other hand, XBCs show an average HRX of −0.56 ± 0.18, while BZBs show an average HRX = −0.56 ± 0.13, and background/foreground objects show HRX = −0.41 ± 0.22.
We applied a Kolmogorov-Smirnov test to the FX and HRX distributions of the XBCs, BZBs, and background/foreground objects in order to test the null hypotesis that the two compared distributions are randomly sampled from a common parent distribution. We found that XBCs and background/foreground objects do not share the same FX distribution, with a negligible p-chance of p < 10−13, while when comparing XBCs to BZBs we obtain a p-chance of p < 10−7. For background/foreground objects and BZBs, the resulting p-chance is also negligible, being less than 1 × 10−15.
Regarding the distributions of HRX, we find a p-chance of p < 0.002 when comparing XBCs to background/foreground objects, and p < 10−7 when comparing BZBs to background/foreground objects. However, we cannot discard that XBCs and BZBs share a common parent distribution, since when comparing their HRX distributions we find a p-chance of p ∼ 0.6.
The subsample of XBCs with a radio counterpart has an average HRX = −0.60 ± 0.14, compatible with that of XBCs. This is an indication that XBCs could belong to the HBL class, since in Paper I we found they follow the same pattern, with an average HRX = −0.63 ± 0.09. When compared with a Kolmogorov-Smirnov test, we obtain a p-chance of p = 0.13 for these two distributions, which is consistent with them sharing a common parent distribution.
The discovery of new HBL sources among samples of X-ray and radio-selected objects has already been indicated in recent studies: the HSP catalogs (Arsioli et al. 2015; Chang et al. 2017, 2019) list HBL sources and candidate sources based on their X-ray-to-radio flux ratios and infrared properties, many of which (∼25%) are also γ-ray emitters (Arsioli & Chang 2017; Arsioli et al. 2018, 2020). Indeed, a subsample of our XBCs was listed in the latest version of the HSP catalog, 3HSP. We discuss this sample in more detail in Sect. 4.4.1.
Moreover, for the 27 cases of Fermi UGSs with two XBCs lying within the same 3FGL uncertainty ellipse, we plot both sources but we indicate with filled black circles those with the largest angular separation between the X-ray and the UGS γ-ray position.
For those XBCs with two X-ray sources within the Fermi positional uncertainty regions, the most distant X-ray counterpart tends to show FX/HRX values that differ from the majority of the XBCs, while the opposite is true for the closest counterparts.
In Fig. 8 we show the Fγ in the 100 MeV−100 GeV band versus FX for all XBCs. We also compare them to our sample of Fermi BZBs. We also report the Fγ thresholds above which 100%, 98%, and 96% of the BZB sample would be detected in X-rays when their exposure times have been scaled down to 5 ks, as shown in Paper I. These thresholds are plotted in solid, dashed, and dotted black lines, respectively. They serve as a comparison with BZBs since our XBC sample has an average exposure time of 4.2 ks.
If the counterpart of an UGS is indeed a BL Lac, above this Fγ threshold we expect to see its X-ray counterpart. Those UGSs emitting above this threshold and lacking X-ray counterparts could probably not be BZBs.
XBCs are fainter in γ-ray flux than BZBs, as shown in Fig. 8. As Fermi BZBs analyzed in Paper I were the brightest in γ-rays, while in this work we are sampling the faint tail of the γ-ray flux distribution, we expect putative BL Lacs in our sample of UGSs to be less bright in X-rays.
In particular, 69% of the XBCs lie above the 96% BZB Fγ threshold line. The average Fγ for XBCs is Fγ = (4.2 ± 1.1)×10−12 erg cm−2 s−1, while for the Fermi BZBs analyzed in Paper I is Fγ = (8.1 ± 3.9)×10−12 erg cm−2 s−1. Since the largest fraction of the XBC sample is fainter in γ-rays on average than BZBs, this also translates into greater γ-ray positional uncertainties, which could be the main reason why it is challenging to find lower energy counterparts. There is also no visible trend for those with two X-ray sources within their γ-ray positional uncertainty region.
In Fig. 9 we show the γ-ray spectral index versus HRX. There is no clear trend for the general XBC sample or for fields with two XBCs.
![]() |
Fig. 9. γ-ray spectral index from the 3FGL vs, HRX for all XBCs. The symbol coding is the same as in Fig. 7. |
4.3. Searching for additional candidate γ-ray BL Lacs within the OUT sample
In this section we present the results on Fermi UGSs that do not show any X-ray counterpart within their 95% confidence positional uncertainty ellipse, but show at least one X-ray source out of it.
From our sample of 69 OUTs, 8 (∼12%) present analysis flags in the 3FGL, all of them indicating uncertainties or artifacts in the γ-ray analysis due to the adopted γ-ray diffuse model. However, this is in agreement with XBCs, which also present ∼12% of 3FGL analysis flags. We note that 31 out of 69 OUTs (∼45%) were listed in the 2FGL catalog, including all 13 (∼19%) that were listed in the 1FGL catalog.
In Fig. 10, we show HRX versus FX for all OUTs, and we also show BZBs from Paper I. Only a handful of sources lie in the area in which HBLs lie, which corresponds to FX ≥ 10−12. Their average flux is FX = (1.9 ± 0.7)×10−13 erg cm−2 s−1, and their average hardness ratio is HRX = −0.49 ± 0.26. These are compatible with the XBCs behavior in X-rays, albeit with larger uncertainties.
![]() |
Fig. 10. HRX vs. FX in the 0.5−10 keV band for all the sample of OUTs. OUTs are plotted as empty black circles, while Fermi BZBs are plotted as grey filled squares. |
In Fig. 11, we show Fγ in the 100 MeV−100 GeV band against the full-band FX, again with the thresholds above which 100%, 98%, and 96% of the BZBs should be detected when observed for 5 ks. Of the 105 OUTs, 79 (75%) display Fγ values above the 96% threshold. In particular, there are 8 sources for which the Fγ is relatively high, lying above the 100% threshold line, but their FX is among the lowest values in the sample. These sources are most probably contaminants of the sample and are likely not associated to UGSs.
![]() |
Fig. 11. Fγ in the 100 MeV−100 GeV band vs. FX in the 0.5−10 keV band, from the Fermi 3FGL catalog, for all OUTs in our sample. The symbol coding is the same as in Fig. 10. The Fγ thresholds above which 100%, 98%, and 96% of the BZB sample would be detected in X-rays when their exposure times have been scaled down to 5 ks, are plotted in solid, dashed and dotted black lines, respectively. |
4.4. Multifrequency analysis of candidate γ-ray BL Lacs
As previously stated, we crossmatched both the XBC and the OUT subsamples with the AllWISE, NVSS and SUMSS catalogs: 161 and 100 of XBCs and OUTs, respectively, have a WISE mid-infrared counterpart, while 80 and 22 have a radio counterpart. Of the latter, 76 of the XBCs and 21 of the OUTs have both infrared and radio counterparts. Optical, mid-infrared, and radio counterparts are present for 25 XBCs and 5 OUTs (see Table 2).
4.4.1. X-ray blazar candidates
In Fig. 12 we show the [3.4]−[4.6] μm mid-infrared color versus FX. In Paper I, we found that BZBs are bluer in mid-infrared when brighter in X-rays. This is also the case for XBCs, although no strict correlation is found in this case.
![]() |
Fig. 12. [3.4]−[4.6] μm mid-infrared color vs. FX in the 0.5−10 keV band for all XBCs of the sample with a WISE counterpart. The symbol coding is the same as in Fig. 7. Sources with μm colors close to zero are likely due to contamination from normal elliptical galaxies. |
More than 80% of XBCs lie within the BZB area in Fig. 12, with only a handful of sources showing [3.4]−[4.6] μm colors close to zero, which are likely due to contamination from normal elliptical galaxies (Massaro et al. 2012b).
We looked for WISE counterparts of XBCs that showed BZB characteristics. D’Abrusco et al. (2019) built two catalogs of likely mid-infrared counterparts of BZB sources, WIBRALS2, and KDEBLLACS, depending on their infrared and radio properties. We selected WISE counterparts of XBCs that were compatible only with the infrared method used to build WIBRALS2 and KDEBLLACS. In particular, WIBRALS2 depends on all four WISE bands, including the 22 μm band, while KDEBLLACS relies only on the 3.4 μm, 4.6 μm, and 12 μm bands. We note that out of the 161 XBCs with a WISE counterpart, 74 are detected in all four WISE bands while the remaining 87 are not. Out of these, 33 and 49, respectively, are compatible with the infrared color model used to find likely BZBs. Indeed, 23 and 13 are even listed in the WIBRALS2 and KDEBLLACS catalogs, meaning they also comply with the radio selection criteria applied by D’Abrusco et al. (2019). Thus, in total, 82 XBCs show BZB-like infrared colors.
We also checked the optical colors of the XBCs. In Fig. 13 we show the (u − r) optical color distribution for all 49 XBCs for which we found an SDSS counterpart. There are two sources that show a (u − r) color index greater than 5, not shown in Fig. 13. Following Massaro et al. (2012b), (u − r) colors lower than 1.4 are a signature of BL Lac sources. Of this subsample, 34 sources (69%) indeed show BL Lac colors, the average being (u − r) = 0.9 ± 0.3. Of these, 33 have a WISE counterpart, 18 (55%) of which are compatible with BZB-like WISE colors.
![]() |
Fig. 13. Distribution of the (u − r) optical color taken from SDSS for all XBCs with an optical counterpart. The vertical solid line marks the 1.4 limit below which sources are expected to be of the BL Lac type following Massaro et al. (2012b). |
Of the 18 XBCs with optical and mid-infrared BZB-like counterparts, 11 also have a radio counterpart. These 11 sources with radio counterparts, (u − r) colors of BL Lacs, and a WISE counterpart with BZB-like mid-infrared colors are very likely to be classified as BZBs.
Moreover, 50 of the XBCs were further classified: 3 as quasars (QSO), 2 as BZQs (flat spectrum radio quasars, as defined by Roma-BZCAT), and 45 as BZBs, all confirmed through optical spectroscopy; the distinction between BZQs and QSOs was made with the use of radio data. These numbers are consistent with our expectations of finding BZBs.
Of these 50 XBCs, 27 were pointed during our ongoing optical spectroscopical campaign aimed at associating counterparts to Fermi UGSs (Paggi et al. 2014; Massaro et al. 2015e; Landoni et al. 2015; Ricci et al. 2015; Álvarez Crespo et al. 2016b,c; Peña-Herazo et al. 2017, 2019; Marchesini et al. 2019b). In these works, we classified sources as BZBs when they showed featureless optical spectra with a dominant blue continuum (i.e. with emission lines with restframe equivalent widths of less than 5 Å; Laurent-Muehleisen et al. 1998). We also distinguished BL Lac sources from those with a strong host galaxy contribution, by measuring their relative flux depression towards the CaII line break following Stocke et al. (1991). Thus, we could ensure the pure BL Lac nature of our sources, all lying within the positional uncertainty ellipse of Fermi UGSs.
The remaining 23 XBCs were classified through optical spectroscopy in the literature (Paggi et al. 2014; Marchesini et al. 2016; Paiano et al. 2017a, 2019; Marchesi et al. 2018; Desai et al. 2019).
We find that for this subsample of XBCs already associated to BZBs, (u − r) = 1.0 ± 0.2, FX = (7.1 ± 5.1)×10−13 erg cm−2 s−1, and HRX = −0.60 ± 0.12, on average. These values are indeed all consistent with the BZB population, as shown in Paper I, and with the whole XBC sample analyzed in this work. This can be seen in Fig. 14, where we show again HRX versus FX for the XBC sample, but marking the sources that have been classified as BZBs with black circles, the BZQs with triangles, and the QSOs with diamonds.
![]() |
Fig. 14. HRX for all the XBC sample vs. FX in the 0.5−10 keV band. The symbol coding is the same as in Fig. 7. XBCs that have been classified as BZBs are plotted as empty black circles, while BZQs and QSOs as black triangles and diamonds, respectively. The solid black line represents a threshold in FX above which 68% of all XBCs have already been identified and only as BZBs. |
All 45 XBCs classified as BZBs through optical spectroscopy have a mid-infrared counterpart. However, 12 of these BZBs do not display typical mid-infrared BZB colors. Moreover, they are not listed in either the second WISE Blazar-like Radio-Loud Sources (WIBRaLS2) catalog or in the KDEBLLACS catalog (D’Abrusco et al. 2019), which are based on a generalization of the mid-infrared and radio properties of the γ-ray BL Lac population. We show these 12 BZBs in Table 3. The majority of these sources, 10 out of 12, are listed in the latest catalog of high-synchrotron peaked blazars, 3HSP (Chang et al. 2019). This catalog, an updated version of the AllWISE based 1WHSP (Arsioli et al. 2015) and 2WHSP (Chang et al. 2017) catalogs, includes sources with typical HBL X-ray-to-radio flux ratios. Given their selection process, it is probable that these 10 candidates are not only blazars but also of the extreme HBL kind, which would explain their peculiar SED characteristics. This further strengthens the selection of XBCs as a method to find γ-ray BZBs within Fermi UGSs, as ∼27% of the XBCs classified as BZBs do not show canonical mid-infrared or radio properties.
Fermi, Swift, and WISE designations for XBCs that were confirmed as BZBs through optical spectroscopy and have a mid-infrared counterpart that is not compatible with the WIBRALS2 or KDEBLLACS color models for typical BZBs.
Moreover, 68% of XBCs with fluxes FX ≥ 8 × 10−13 erg cm−2 s−1 were associated with BZBs. We plot this threshold as a solid black line in Fig. 14. We expect that the remaining 32% of XBCs to be likely BZBs, which could be confirmed through optical follow-up observations.
In Fig. 15 we show the distribution of Fγ in the 100 MeV−100 GeV band for the XBCs already classified through optical spectroscopy. We mark again the detection thresholds as in Fig. 8. We note that 73% of this subsample is above the 96% threshold. This is agreement with the behaviour of the parent XBC sample shown in Fig. 8.
![]() |
Fig. 15. Distribution of Fγ in the 100 MeV−100 GeV band for the XBCs already classified through optical spectroscopy. The thresholds above which 100%, 98%, and 96% of the BZBs are expected to be detected in X-ray when observing for 5 ks or more are plotted as solid, dashed, and dot-dashed black lines, respectively. |
Finally, in Fig. 16 we show the (u − r) color distribution for the 29 classified XBCs that have a counterpart in SDSS, as in Fig. 15. There are six sources above the 1.4 limit, two of which are BZQs and four are BZBs. The remaining sources lie below the 1.4 threshold. Following Massaro et al. (2012b), the ∼21% of the sample showing (u − r) color indices above the 1.4 limit can be either high-redshift objects (with z > 0.5) or galaxy-dominated BZBs, meaning their jets are undergoing a low-activity phase and therefore the host galaxy becomes the dominant feature in its optical spectrum.
![]() |
Fig. 16. Distribution of the (u − r) optical color taken from SDSS for all XBCs with an optical counterpart associated through optical spectroscopy and classified as either BZBs, BZQs, or QSOs. The vertical solid line marks the limit below which sources are expected to be of the BL Lac type following Massaro et al. (2012b). |
The fact that XBCs closely follow the BZB trends in mid-infrared colors, X-ray hardness ratio values, and optical color indices indicates that they share a similar spectral shape in all these bands with BZBs. The main difference between XBCs and BZBs lies in their fluxes, XBCs being fainter than BZBs in both γ-rays and X-rays.
We selected a subsample of XBCs that are likely BZB sources with a high confidence degree. We selected all XBCs that satisfied at least two of the following criteria:
-
XBCs that have a radio counterpart,
-
XBCs that have an optical counterpart with color index (u − r) < 1.4,
-
XBCs that have a mid-infrared counterpart which is compatible with either the WIBRALS2 or the KDEBLLACS mid-infrared color model for BZBs.
Thus, we obtained a subsample of 54 XBCs. Finally, we discarded 35 sources that had already been associated in the literature or during our follow-up campaigns. The remaining 19 form part of our catalog of likely BZBs that will be part of our future optical spectroscopic follow-up campaigns.
We highlight the fact that with this selection process, all XBCs correspond to a single UGS field. This means that, in cases in which we found two XBCs within the same UGS uncertainty ellipse, at least one of them did not satisfy two of the three criteria stated above. We list the whole sample of 30 X-ray sources that are candidate γ-ray BL Lacs in Table 4. In Col. 1 we report the 3FGL name, in Col. 2 the Swift/XRT source designation, in Col. 3 the angular separation between the 3FGL and XRT positions, and in Col. 4 the class we assigned to it. This table includes sources selected from the OUT sample (see Sect. 4.4.2). We note that nine of the sources included in Table 4 also belong to the 3HSP catalog (Chang et al. 2019).
Refined sample of XBCs and outliers that are with high confidence likely BL Lac sources.
4.4.2. Outliers
The sample of OUTs show greater variation in their properties. In Fig. 17, we show the [3.4]−[4.6] μm mid-infrared color versus FX, again comparing with BZBs and their trends as explained in Fig. 12.
![]() |
Fig. 17. [3.4]−[4.6] μm mid-infrared color FX in the 0.5−10 keV band, for all the sample of OUTs with a WISE counterpart. The symbol coding is the same as in Fig. 10. |
There are 18 sources that show [3.4]−[4.6] μm ≈ 0, which is the typical color of elliptical galaxies, or of BL Lacs for which the contribution of the host galaxy is dominant. The remaining sources show values similar to BZBs, with all but nine of the OUTs lying towards the fainter side in X-rays, and being redder in mid-infrared (FX < 1.0 × 10−12, [3.4]−[4.6]< 1.5). Again, sources with a radio counterpart tend to group together in the area covered by BZBs.
When checking their mid-infrared colors, 23 OUTs are compatible with the WIBRALS2 color model, while another 13 sources with no detection at 22 μm are compatible with the KDEBLLACS infrared color model (D’Abrusco et al. 2019). This means that a total of 36 OUTs display infrared colors of BZBs, with 8 of these having also a radio counterpart; these latter are therefore listed in both the WIBRALS2 and KDEBLLACS catalogs of infrared BZB-like sources.
The majority of the OUT sample, that is, 64 sources, are not compatible with either infrared color model. However, we stress that there is still a non-negligible (at least 30%) number of potential BZBs within the OUTs sample.
Regarding the optical colors of OUTs, in Fig. 18 we show their (u − r) distribution. There are 34 OUTs with a SDSS counterpart; 14 of these show (u − r) < 1.4 mag, which is typical of BZBs. All 34 show large scatter, their average color index being (u − r) = 1.7 ± 1.2 mag.
![]() |
Fig. 18. Distribution of the (u − r) optical color taken from SDSS, for all OUTs with an optical counterpart. The vertical solid line marks the limit below which sources are expected to be of the BL Lac type, following Massaro et al. (2012b). |
When crossmatching the OUTs sample with sources identified or associated in our campaigns or in the literature, we found a classification for only five of them: two BZQs, two QSOs and one BZB. We show them plotted in Fig. 19; they all indeed lie closer to the BZB region of the plot than the majority of the OUTs. We note that there are 13 OUTs lying within the BZB area and above the 68% FX threshold reported in Fig. 14. None of these have been associated through optical spectroscopy, and are likely BZB candidates for future optical follow-up campaigns.
As for XBCs, we selected a subsample of OUTs that are likely BZBs based on the same criteria stated in Sect. 4.4.1 and discarding any sources that had already been associated. We obtained a subsample of 11 OUTs, which we listed in Table 4.
5. Summary and conclusions
In this work we analyze the X-ray properties of Fermi UGSs, in particular to search for BZB candidates through the recently proposed X-ray–γ-ray connection, and comparing it with the already known connections for BZBs between γ-rays and radio, infrared, and optical wavelengths. To this aim, we built a sample of 327 UGSs listed by Fermi that were observed by the Swift/XRT telescope, collected up to December 2018.
There are 223 out of 327 Fermi UGSs observed by Swift/XRT that present at least one X-ray detection in the XRT field. Out of these 223 UGSs, 148 UGSs have at least one XRT source within their positional uncertainty ellipse: We call these X-ray blazar candidates, or XBCs. They amount to 175 sources. A further 69 UGSs do not have any X-ray source within their positional uncertainty ellipse, but have at least one X-ray source in the field. This second sample, labeled as OUTs, lists 105 X-ray sources.
XBCs present X-ray fluxes, hardness ratio values, and mid-infrared and optical colors similar to BZBs, strongly suggesting a BZB nature. OUTs show a combination of BZB and nonBZB characteristics, indicating there are a number of BZBs within the sample but with a high degree of contamination.
There are 45 XBCs which were already classified as BZBs through optical spectroscopy. Of these, 12 (∼27%) do not show canonical mid-infrared or radio properties. This indicates that the selection of XBCs is a strong method to find γ-ray BZB candidates, that are not found with other multiwavelength methods, within Fermi UGSs.
Based on a combination of their X-ray and multiwavelength properties, we selected a sample of X-ray sources that are very likely BL Lac candidates. These constitute a list of 32 sources that merit follow-up optical spectroscopic observations to confirm their nature.
The estimate of the X-ray flux only slightly depends on the choice of the photon index for the values typically observed in BZBs as in other radio-loud active galaxies, and it can be considered within the uncertainty of the X-ray flux itself (see e.g., Massaro et al. 2010, 2013b, 2015b, and references therein).
Acknowledgments
E. J. Marchesini would like to thank the anonymous referee on behalf of all the authors for the insights given to improve this manuscript. E. J. Marchesini would also like to thank Dr. Rocío I. Páez and Dr. M. Victoria Reynaldi for useful discussions on this work. This work is supported by the “Departments of Excellence 2018 – 2022” Grant awarded by the Italian Ministry of Education, University and Research (MIUR) (L. 232/2016). This research has made use of resources provided by the Compagnia di San Paolo for the grant awarded on the BLENV project (S1618_L1_MASF_01) and by the Ministry of Education, Universities and Research for the grant MASF_FFABR_17_01. F. M. acknowledges financial contribution from the agreement ASI-INAF n.2017-14-H.0 A. P. acknowledges financial support from the Consorzio Interuniversitario per la fisica Spaziale (CIFS) under the agreement related to the grant MASF_CONTR_FIN_18_02. This research has made use of data obtained from the high- energy Astrophysics Science Archive Research Center (HEASARC) provided by NASA’s Goddard Space Flight Center. Part of this work is based on the NVSS (NRAO VLA Sky Survey). The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under contract with the National Science Foundation. The Molonglo Observatory site manager, Duncan Campbell- Wilson, and the staff, Jeff Webb, Michael White, and John Barry, are responsible for the smooth operation of the Molonglo Observatory Synthesis Telescope (MOST) and the day-to-day observing program of SUMSS. SUMSS is dedicated to Michael Large, whose expertise and vision made the project possible. The MOST is operated by the School of Physics with the support of the Australian Research Council and the Science Foundation for Physics within the University of Sydney. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. TOPCAT (http://www.star.bris.ac.uk/~mbt/topcat/) (Taylor 2005) and STILTS (Taylor 2006) were used for the preparation and manipulation of the images and the tabular data.
References
- Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, ApJS, 188, 405 [NASA ADS] [CrossRef] [Google Scholar]
- Abdollahi, S., Acero, F., Ackermann, M., et al. 2020, ApJS, 247, 33 [NASA ADS] [CrossRef] [Google Scholar]
- Acero, F., Donato, D., Ojha, R., et al. 2013, ApJ, 779, 133 [NASA ADS] [CrossRef] [Google Scholar]
- Acero, F., Ackermann, M., Ajello, M., et al. 2015, ApJS, 218, 23 [Google Scholar]
- Ackermann, M., Ajello, M., Allafort, A., et al. 2011, ApJ, 741, 30 [NASA ADS] [CrossRef] [Google Scholar]
- Ackermann, M., Ajello, M., Ballet, J., et al. 2012a, ApJ, 751, 159 [NASA ADS] [CrossRef] [Google Scholar]
- Ackermann, M., Ajello, M., Allafort, A., et al. 2012b, ApJ, 753, 83 [NASA ADS] [CrossRef] [Google Scholar]
- Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al. 2012, ApJS, 203, 21 [NASA ADS] [CrossRef] [Google Scholar]
- Álvarez Crespo, N., Massaro, F., D’Abrusco, R., et al. 2016a, Ap&SS, 361, 316 [NASA ADS] [CrossRef] [Google Scholar]
- Álvarez Crespo, N., Masetti, N., Ricci, F., et al. 2016b, AJ, 151, 32 [NASA ADS] [CrossRef] [Google Scholar]
- Álvarez Crespo, N., Massaro, F., Milisavljevic, D., et al. 2016c, AJ, 151, 95 [NASA ADS] [CrossRef] [Google Scholar]
- Arsioli, B., & Chang, Y. L. 2017, A&A, 598, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Arsioli, B., Fraga, B., Giommi, P., Padovani, P., & Marrese, P. M. 2015, A&A, 579, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Arsioli, B., Barres de Almeida, U., Prandini, E., Fraga, B., & Foffano, L. 2018, MNRAS, 480, 2165 [NASA ADS] [CrossRef] [Google Scholar]
- Arsioli, B., Chang, Y. L., & Musiimenta, B. 2020, MNRAS, 493, 2438 [CrossRef] [Google Scholar]
- Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 2009, ApJ, 697, 1071 [NASA ADS] [CrossRef] [Google Scholar]
- Blandford, R. D., & Rees, M. J. 1978, in BL Lac Objects, ed. A. M. Wolfe, 328 [Google Scholar]
- Blinov, D., Pavlidou, V., Papadakis, I., et al. 2018, MNRAS, 474, 1296 [NASA ADS] [CrossRef] [Google Scholar]
- Burrows, D. N., Hill, J. E., Nousek, J. A., et al. 2005, Space Sci. Rev., 120, 165 [NASA ADS] [CrossRef] [Google Scholar]
- Capalbi, M., Perri, M., Saija, B., & Tamburelli, F. 2005, ASI Science Data Center, 1 [Google Scholar]
- Chang, Y. L., Arsioli, B., Giommi, P., & Padovani, P. 2017, A&A, 598, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Chang, Y. L., Arsioli, B., Giommi, P., Padovani, P., & Brandt, C. H. 2019, A&A, 632, A77 [CrossRef] [EDP Sciences] [Google Scholar]
- Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 115, 1693 [NASA ADS] [CrossRef] [Google Scholar]
- D’Abrusco, R., Massaro, F., Ajello, M., et al. 2012, ApJ, 748, 68 [NASA ADS] [CrossRef] [Google Scholar]
- D’Abrusco, R., Massaro, F., Paggi, A., et al. 2013, ApJS, 206, 12 [NASA ADS] [CrossRef] [Google Scholar]
- D’Abrusco, R., Massaro, F., Paggi, A., et al. 2014, ApJS, 215, 14 [NASA ADS] [CrossRef] [Google Scholar]
- D’Abrusco, R., Álvarez Crespo, N., Massaro, F., et al. 2019, ApJS, 242, 4 [NASA ADS] [CrossRef] [Google Scholar]
- D’Elia, V., Perri, M., Puccetti, S., et al. 2013, A&A, 551, A142 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Desai, A., Marchesi, S., Rajagopal, M., & Ajello, M. 2019, ApJS, 241, 5 [CrossRef] [Google Scholar]
- Doert, M., & Errando, M. 2014, ApJ, 782, 41 [NASA ADS] [CrossRef] [Google Scholar]
- Dunkley, J., Komatsu, E., Nolta, M. R., et al. 2009, ApJS, 180, 306 [NASA ADS] [CrossRef] [Google Scholar]
- Evans, P. A., Osborne, J. P., Beardmore, A. P., et al. 2014, ApJS, 210, 8 [NASA ADS] [CrossRef] [Google Scholar]
- Falcone, A., Stroh, M., & Pryal, M. 2014, Am. Astron. Soc. Meet. Abstr., 223, 301.05 [NASA ADS] [Google Scholar]
- Falomo, R., Treves, A., Scarpa, R., Paiano, S., & Landoni, M. 2017, MNRAS, 470, 2814 [CrossRef] [Google Scholar]
- Franceschini, A., Foffano, L., Prandini, E., & Tavecchio, F. 2019, A&A, 629, A2 [CrossRef] [EDP Sciences] [Google Scholar]
- Ghirlanda, G., Ghisellini, G., Tavecchio, F., & Foschini, L. 2010, MNRAS, 407, 791 [NASA ADS] [CrossRef] [Google Scholar]
- Hassan, T., Mirabal, N., Contreras, J. L., & Oya, I. 2013, MNRAS, 428, 220 [NASA ADS] [CrossRef] [Google Scholar]
- Healey, S. E., Romani, R. W., Taylor, G. B., et al. 2007, ApJS, 171, 61 [NASA ADS] [CrossRef] [Google Scholar]
- Hovatta, T., Lister, M. L., Aller, M. F., et al. 2012, AJ, 144, 105 [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]
- Kaur, A., Ajello, M., Marchesi, S., & Omodei, N. 2019, ApJ, 871, 94 [NASA ADS] [CrossRef] [Google Scholar]
- Landi, R., Bassani, L., Stephen, J. B., et al. 2015, A&A, 581, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Landoni, M., Falomo, R., Treves, A., & Sbarufatti, B. 2014, A&A, 570, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Landoni, M., Massaro, F., Paggi, A., et al. 2015, AJ, 149, 163 [NASA ADS] [CrossRef] [Google Scholar]
- Landoni, M., Paiano, S., Falomo, R., Scarpa, R., & Treves, A. 2018, ApJ, 861, 130 [CrossRef] [Google Scholar]
- Laurent-Muehleisen, S. A., Kollgaard, R. I., Ciardullo, R., et al. 1998, ApJS, 118, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Lico, R., Giroletti, M., Orienti, M., & D’Ammando, F. 2016, A&A, 594, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lister, M. L., Aller, M. F., Aller, H. D., et al. 2013, AJ, 146, 120 [NASA ADS] [CrossRef] [Google Scholar]
- Mahony, E. K., Sadler, E. M., Murphy, T., et al. 2010, ApJ, 718, 587 [NASA ADS] [CrossRef] [Google Scholar]
- Marchesi, S., Kaur, A., & Ajello, M. 2018, AJ, 156, 212 [CrossRef] [Google Scholar]
- Marchesini, E. J., Masetti, N., Chavushyan, V., et al. 2016, A&A, 596, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Marchesini, E. J., Paggi, A., Massaro, F., et al. 2019a, A&A, 631, A150 (Paper I) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Marchesini, E. J., Peña-Herazo, H. A., Álvarez Crespo, N., et al. 2019b, Ap&SS, 364, 5 [NASA ADS] [CrossRef] [Google Scholar]
- Maselli, A., Massaro, E., Nesci, R., et al. 2010, A&A, 512, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Maselli, A., Massaro, F., Cusumano, G., et al. 2013, ApJS, 206, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Masetti, N., Sbarufatti, B., Parisi, P., et al. 2013, A&A, 559, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Massaro, F., & D’Abrusco, R. 2016, ApJ, 827, 67 [NASA ADS] [CrossRef] [Google Scholar]
- Massaro, F., Giommi, P., Tosti, G., et al. 2008, A&A, 489, 1047 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Massaro, F., Harris, D. E., Tremblay, G. R., et al. 2010, ApJ, 714, 589 [NASA ADS] [CrossRef] [Google Scholar]
- Massaro, F., D’Abrusco, R., Tosti, G., et al. 2012a, ApJ, 752, 61 [NASA ADS] [CrossRef] [Google Scholar]
- Massaro, E., Nesci, R., & Piranomonte, S. 2012b, MNRAS, 422, 2322 [NASA ADS] [CrossRef] [Google Scholar]
- Massaro, F., D’Abrusco, R., Giroletti, M., et al. 2013a, ApJS, 207, 4 [NASA ADS] [CrossRef] [Google Scholar]
- Massaro, F., Harris, D. E., Tremblay, G. R., et al. 2013b, ApJS, 206, 7 [NASA ADS] [CrossRef] [Google Scholar]
- Massaro, F., Masetti, N., D’Abrusco, R., Paggi, A., & Funk, S. 2014, AJ, 148, 66 [NASA ADS] [CrossRef] [Google Scholar]
- Massaro, F., D’Abrusco, R., Landoni, M., et al. 2015a, ApJS, 217, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Massaro, F., Harris, D. E., Liuzzo, E., et al. 2015b, ApJS, 220, 5 [NASA ADS] [CrossRef] [Google Scholar]
- Massaro, E., Maselli, A., Leto, C., et al. 2015c, Ap&SS, 357, 75 [NASA ADS] [CrossRef] [Google Scholar]
- Massaro, F., Thompson, D. J., & Ferrara, E. C. 2015d, A&ARv, 24, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Massaro, F., Landoni, M., D’Abrusco, R., et al. 2015e, A&A, 575, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Massaro, F., Álvarez Crespo, N., D’Abrusco, R., et al. 2016, Ap&SS, 361, 337 [NASA ADS] [CrossRef] [Google Scholar]
- Mauch, T., Murphy, T., Buttery, H. J., et al. 2003, MNRAS, 342, 1117 [NASA ADS] [CrossRef] [Google Scholar]
- Mukai, K. 1993, Legacy, 3, 21 [Google Scholar]
- Nolan, P. L., Abdo, A. A., Ackermann, M., et al. 2012, ApJS, 199, 31 [NASA ADS] [CrossRef] [Google Scholar]
- Nori, M., Giroletti, M., Massaro, F., et al. 2014, ApJS, 212, 3 [NASA ADS] [CrossRef] [Google Scholar]
- Padovani, P., & Giommi, P. 1995, ApJ, 444, 567 [NASA ADS] [CrossRef] [Google Scholar]
- Paggi, A., Massaro, F., D’Abrusco, R., et al. 2013, ApJS, 209, 9 [NASA ADS] [CrossRef] [Google Scholar]
- Paggi, A., Milisavljevic, D., Masetti, N., et al. 2014, AJ, 147, 112 [NASA ADS] [CrossRef] [Google Scholar]
- Paiano, S., Falomo, R., Franceschini, A., Treves, A., & Scarpa, R. 2017a, ApJ, 851, 135 [CrossRef] [Google Scholar]
- Paiano, S., Franceschini, A., & Stamerra, A. 2017b, MNRAS, 468, 4902 [NASA ADS] [CrossRef] [Google Scholar]
- Paiano, S., Falomo, R., Treves, A., Franceschini, A., & Scarpa, R. 2019, ApJ, 871, 162 [CrossRef] [Google Scholar]
- Peña-Herazo, H. A., Marchesini, E. J., Álvarez Crespo, N., et al. 2017, Ap&SS, 362, 228 [NASA ADS] [CrossRef] [Google Scholar]
- Peña-Herazo, H. A., Massaro, F., Chavushyan, V., et al. 2019, Ap&SS, 364, 85 [NASA ADS] [CrossRef] [Google Scholar]
- Pushkarev, A. B., Kovalev, Y. Y., Lister, M. L., & Savolainen, T. 2009, A&A, 507, L33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ricci, F., Massaro, F., Landoni, M., et al. 2015, AJ, 149, 160 [NASA ADS] [CrossRef] [Google Scholar]
- Salvetti, D., Chiaro, G., La Mura, G., & Thompson, D. J. 2017, MNRAS, 470, 1291 [NASA ADS] [CrossRef] [Google Scholar]
- Sandrinelli, A., Treves, A., Falomo, R., et al. 2013, AJ, 146, 163 [CrossRef] [Google Scholar]
- Stephen, J. B., Bassani, L., Landi, R., et al. 2010, MNRAS, 408, 422 [NASA ADS] [CrossRef] [Google Scholar]
- Stickel, M., Padovani, P., Urry, C. M., Fried, J. W., & Kuehr, H. 1991, ApJ, 374, 431 [NASA ADS] [CrossRef] [Google Scholar]
- Stocke, J. T., Morris, S. L., Gioia, I. M., et al. 1991, ApJS, 76, 813 [NASA ADS] [CrossRef] [Google Scholar]
- Stroh, M. C., & Falcone, A. D. 2013, ApJS, 207, 28 [NASA ADS] [CrossRef] [Google Scholar]
- Takahashi, Y., Kataoka, J., Nakamori, T., et al. 2012, ApJ, 747, 64 [CrossRef] [Google Scholar]
- Taylor, M. B. 2005, in Astronomical Data Analysis Software and Systems XIV, eds. P. Shopbell, M. Britton, & R. Ebert, ASP Conf. Ser., 347, 29 [Google Scholar]
- Taylor, M. B. 2006, in Astronomical Data Analysis Software and Systems XV, eds. C. Gabriel, C. Arviset, D. Ponz, & S. Enrique, ASP Conf. Ser., 351, 666 [Google Scholar]
- Taylor, G. B., Healey, S. E., Helmboldt, J. F., et al. 2007, ApJ, 671, 1355 [NASA ADS] [CrossRef] [Google Scholar]
- Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Results from the positional crossmatch of XBCs and OUTs with radio, mid-infrared, and optical catalogs.
Fermi, Swift, and WISE designations for XBCs that were confirmed as BZBs through optical spectroscopy and have a mid-infrared counterpart that is not compatible with the WIBRALS2 or KDEBLLACS color models for typical BZBs.
Refined sample of XBCs and outliers that are with high confidence likely BL Lac sources.
All Figures
![]() |
Fig. 1. Flow chart to highlight all steps followed to build our final UGS sample. |
In the text |
![]() |
Fig. 2. Full-band (0.5−10 keV) merged XRT image corresponding to the Fermi UGS 3FGL J0121.8−3917. The yellow dashed line indicates the 3FGL positional uncertainty ellipse at 95% level of confidence, while the white circles show the XRT sources detected with S/N greater than 3. There is only one XRT counterpart within the positional uncertainty ellipse of 3FGL J0121.8−3917. The image was smoothed with a Gaussian kernel radius of 5 pixels. |
In the text |
![]() |
Fig. 3. Full-band (0.5−10 keV) merged XRT image corresponding to the Fermi UGS 3FGL J0809.3−0941. The yellow dashed line indicates the 3FGL positional uncertainty ellipse at 95% confidence level, while the white circles show all XRT sources detected with S/N greater than 3. We find two X-ray sources lying within the Fermi ellipse, which we label both as XBCs; the remaining one is discarded from this analysis. The image was smoothed with a Gaussian kernel of radius 5 pixels. |
In the text |
![]() |
Fig. 4. Full-band (0.5−10 keV) merged XRT image corresponding to the Fermi UGS 3FGL J0216.0+0300. The yellow dashed line indicates the 3FGL positional uncertainty ellipse at 95% confidence level, while the white circles show all XRT sources detected with S/N greater than 3. We find three OUTs in this field, since no X-ray source is detected within the Fermi ellipse. The image was smoothed with a Gaussian kernel of radius 5 pixels. |
In the text |
![]() |
Fig. 5. Angular separation between Swift/XRT and 3FGL positions for XBCs (in purple) and OUTs (in black) in arcminutes. There is an overlap between XBCs and OUTs, implying the latter could be associated in future Fermi catalogs. |
In the text |
![]() |
Fig. 6. Flow chart to highlight all steps followed to build our final XBC and OUT samples. |
In the text |
![]() |
Fig. 7. HRX versus FX in the 0.5−10 keV band for all sources in the XBC sample. XBCs are plotted in empty purple circles, while Fermi BZBs are plotted in grey filled squares. When two XBCs lie within the positional uncertainty area of the same UGS, the most distant one is marked with filled black circles. |
In the text |
![]() |
Fig. 8. γ-ray energy flux in the 100 MeV−100 GeV band vs. FX from the Fermi 3FGL catalog for all XBCs in our sample. The symbol coding is the same as in Fig. 7. The Fγ thresholds above which 100%, 98%, and 96% of the BZB sample would be detected in X-rays when their exposure times have been scaled down to 5 ks are plotted in solid, dashed, and dotted black lines, respectively. |
In the text |
![]() |
Fig. 9. γ-ray spectral index from the 3FGL vs, HRX for all XBCs. The symbol coding is the same as in Fig. 7. |
In the text |
![]() |
Fig. 10. HRX vs. FX in the 0.5−10 keV band for all the sample of OUTs. OUTs are plotted as empty black circles, while Fermi BZBs are plotted as grey filled squares. |
In the text |
![]() |
Fig. 11. Fγ in the 100 MeV−100 GeV band vs. FX in the 0.5−10 keV band, from the Fermi 3FGL catalog, for all OUTs in our sample. The symbol coding is the same as in Fig. 10. The Fγ thresholds above which 100%, 98%, and 96% of the BZB sample would be detected in X-rays when their exposure times have been scaled down to 5 ks, are plotted in solid, dashed and dotted black lines, respectively. |
In the text |
![]() |
Fig. 12. [3.4]−[4.6] μm mid-infrared color vs. FX in the 0.5−10 keV band for all XBCs of the sample with a WISE counterpart. The symbol coding is the same as in Fig. 7. Sources with μm colors close to zero are likely due to contamination from normal elliptical galaxies. |
In the text |
![]() |
Fig. 13. Distribution of the (u − r) optical color taken from SDSS for all XBCs with an optical counterpart. The vertical solid line marks the 1.4 limit below which sources are expected to be of the BL Lac type following Massaro et al. (2012b). |
In the text |
![]() |
Fig. 14. HRX for all the XBC sample vs. FX in the 0.5−10 keV band. The symbol coding is the same as in Fig. 7. XBCs that have been classified as BZBs are plotted as empty black circles, while BZQs and QSOs as black triangles and diamonds, respectively. The solid black line represents a threshold in FX above which 68% of all XBCs have already been identified and only as BZBs. |
In the text |
![]() |
Fig. 15. Distribution of Fγ in the 100 MeV−100 GeV band for the XBCs already classified through optical spectroscopy. The thresholds above which 100%, 98%, and 96% of the BZBs are expected to be detected in X-ray when observing for 5 ks or more are plotted as solid, dashed, and dot-dashed black lines, respectively. |
In the text |
![]() |
Fig. 16. Distribution of the (u − r) optical color taken from SDSS for all XBCs with an optical counterpart associated through optical spectroscopy and classified as either BZBs, BZQs, or QSOs. The vertical solid line marks the limit below which sources are expected to be of the BL Lac type following Massaro et al. (2012b). |
In the text |
![]() |
Fig. 17. [3.4]−[4.6] μm mid-infrared color FX in the 0.5−10 keV band, for all the sample of OUTs with a WISE counterpart. The symbol coding is the same as in Fig. 10. |
In the text |
![]() |
Fig. 18. Distribution of the (u − r) optical color taken from SDSS, for all OUTs with an optical counterpart. The vertical solid line marks the limit below which sources are expected to be of the BL Lac type, following Massaro et al. (2012b). |
In the text |
![]() |
Fig. 19. Same as Fig. 10, but including the associated sources. |
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.