A new multiwavelength census of blazars

Context. Blazars are the rarest and most powerful active galactic nuclei. They play a crucial and expanding role in the multifrequency and multimessenger astrophysics of today. Dominating the high-energy extragalactic sky, they have been recently associated with high-energy astrophysical neutrinos and they may be among the accelerators of the ultra-high energy cosmic rays. Current blazar catalogs, however, are incomplete and they are depleted, in particular, at low Galactic latitudes. Aims. Our study is aimed at augmenting the current blazar census, starting from a sample of ALMA calibrators that provides more homogeneous sky coverage, especially at low Galactic latitudes, to build a catalog of blazar candidates that can provide candidate counterparts to unassociated γ-ray sources and sources of high-energy neutrino emission or ultra-high energy cosmic rays. Methods. Based on the ALMA Calibrator Catalog, we compiled a selection of 1580 ALMA Blazar Candidates (ABCs)for which we collected multiwavelength information, including Gaia photometric, parallax, and proper motion data, as well as SDSS and LAMOST photometric and spectral data, WISE photometric data, X-ray (Swift-XRT, Chandra-ACIS and XMM-Newton-EPIC) count-rates and spectra, and Fermi-LAT fluxes and spectral slopes. We also compared our ALMA Blazar Candidates with existing blazar catalogs, such as4FGL, 3HSP, WIBRaLS2, and the KDEBLLACS. Results. The ABC catalog fills in the gaps with regard to low Galactic latitude sources in current blazar catalogs. By cross-matching this catalog with Gaia DR2, SDSS DR12, LAMOST DR5, AllWISE, and 4FGL catalogs, we obtained 805, 295, 31, 1311, and 259 matches, respectively. ALMA Blazar Candidates are significantly dimmer than known blazars in the Gaia g band, while the difference in the Gaia b - r color between the two populations is less pronounced. Also, ABC sources appear bluer in SDSS than known blazars, although with low statistical significance. Most ABCs classified as QSO and BL Lac fall into the SDSS color regions of low redshiftquasars, with some QSOs entering the regions of higher redshift quasars. We collected 110 optical spectra in SDSS DR12 and LAMOST DR5, which mostly classify the corresponding sources as QSO (98), while 12 sources turned out to be galactic objects. Regarding the WISE colors, we found that ABC sources are significantly bluer than known blazars in the w2 - w3 and w3 - w4 colors. In X-ray, we detected 173 ABCs and we were able to extract X-ray spectra for 92 of them. Our sources are, on average, similar in the X-ray band to known blazars in terms of count rates and spectral slopes, implying that our sample covers the same region of the blazar parameter space in this band. A comparison of γ-ray properties shows that the ABCs are, on average, dimmer and that their γ-ray spectra are softer than known blazars, indicating a significant contribution on the part of FSRQ sources. Making use of WISE colors, we classified 715 ABCs as candidate γ-ray blazar of different classes. Conclusions. We built a new catalog of 1580 candidate blazars with a rich multiwavelength data-set, filling in the gaps for low Galactic latitude sources in current blazar catalogs. This will be particularly important for the identification of the source population of high-energy neutrinos or ultra-high energy cosmic rays, or to verify the Gaia optical reference frame. In addition, ALMA Blazar Candidates can be investigated both through optical spectroscopic observation campaigns or through repeated photometric observations for variability studies. In this context, thedata collected by the upcoming LSST surveys will provide a key tool for investigating the possible blazar nature of these sources.