The study of machine learning (ML) techniques for the autonomous classification of astrophysical sources is of great interest, and we explore its applications in the context of a multifrequency data-frame. We test the use of supervised ML to classify blazars according to its synchrotron peak energy, either lower or higher than 10$^{15}$Hz. We select a blazar sample with 4178 objects labelled as 1279 high synchrotron peak sources (HSP: $ \nu$-peak > 10$^{15}$Hz) and 2899 low synchrotron peak sources (LSP: $ \nu$-peak < 10$^{15}$Hz). A set of multifrequency input features were defined to represent each source, that includes spectral slopes ($\alpha_{\nu_1, \nu_2}$) between the radio, infra-red, optical, and X-ray bands, also considering IR colours. We describe the optimisation of five ML algorithms: Random Forests (RF), Support Vector Machine (SVM), K-Nearest Neighbours (KNN), Gaussian Naive Bayes (GNB) and the Ludwig auto-ML framework, and report on accurate models for blazar classification into LSP and HSP. In our particular case, the SVM algorithm had the best performance, reaching 93% of balanced-accuracy. A joint-feature permutation test revealed that the spectral slopes alpha-radio-IR and alpha-radio-optical are the most relevant for the ML modelling, followed by the IR colours. This work shows that ML algorithms can distinguish multifrequency spectral characteristics and handle the classification task of blazars into LSPs and HSPs, and is a hint for the potential use of ML algorithms for the autonomous determination of broadband spectral parameters (as the synchrotron $\nu$-peak), or even to search for new blazars in all-sky databases.

We present a systematic study of gamma-ray blazar candidates based on a sample of 40 objects taken from the WIBR catalogue. By using a likelihood analysis, 26 of the 40 sources showed significant gamma-ray signatures $\geq3\sigma$. Using high-energy test statistics (TS) maps, we confirm 8 sources, which are completely new, and show another 15 promising $\gamma$-ray candidates. The results from this analysis show that a multi-frequency approach can help to improve the current description of the gamma-ray sky.

This paper presents the results of a $\rm \gamma$-ray likelihood analysis over all the extreme and high synchrotron peak blazars (EHSP & HSP) from the 3HSP catalogue. We investigate 2013 multifrequency positions under the eyes of Fermi Large Area Telescope, considering 11 years of observations in the energy range between 500 MeV to 500 GeV, which results in 1160 $\rm \gamma$-ray signatures detected down to the TS = 9 threshold. The detections include 235 additional sources concerning the Fermi Large Area Telescope Fourth Source Catalog (4FGL), all confirmed via high-energy TS maps, and represent an improvement of ~25% for the number of EHSP & HSP currently described in $\rm \gamma$-rays. We build the $\rm \gamma$-ray spectral energy distribution for all the 1160 2BIGB sources, plot the corresponding $\rm \gamma$-ray logN-logS, and measure their total contribution to the extragalactic gamma-ray background, which reaches up to ~33% at 100 GeV. Also, we show that the $\rm \gamma$-ray detectability improves according to the synchrotron peak flux as represented by the Figure of Merit (FOM) parameter, and note that the search for TeV peaked blazars may benefit from considering HSP and EHSP as a whole, instead of EHSPs only. The 2BIGB acronym stands for `Second Brazil-ICRANet Gamma-ray Blazars' catalogue, and all the broadband models and spectral energy distribution data-points will be available on public data repositories (OpenUniverse, GitHub, and Brazilian Science Data Center-BSDC)

High synchrotron peaked blazars (HSPs or HBLs) play a central role in very high energy (VHE) gamma-ray astronomy, and likely in neutrino astronomy. Currently, the largest compilation of HSP blazars, the 2WHSP sample, includes 1691 sources, but it is not complete neither in the radio nor in the X-ray band. In order to provide a more sizable and more accurate set of HSP blazars that is useful for future statistical studies and to plan for VHE/TeV observations, we present the largest sample of extreme and high synchrotron peaked (EHSP;HSP) blazars and blazar candidates: the 3HSP catalogue. The 3HSP catalogue includes 2013 sources, 88% of which with a redshift estimation, a much larger percentage than in any other list of HSP blazars. All new {\gamma}-ray detections are described in the First and Second Brazil ICRANet gamma-ray blazar catalogues (1BIGB & 2BIGB) also taking into account the 4FGL list of {\gamma}-ray sources published by the Fermi large area telescope (Fermi-LAT) team. Moreover, the cross-matching between the 2WHSP, 2FHL HSPs, and IceCube neutrinos position has suggested that HSPs are likely counterparts of neutrino events, which implies the 3HSP catalogue is useful also in that respect. The 3HSP catalogue features improved completeness compared to its predecessors, the 1WHSP and 2WHSP catalogues, and follows the track of their increasing relevance for VHE astronomy.

We present the dissection in space, time, and energy of the region around the IceCube-170922A neutrino alert. This study is motivated by: (1) the first association between a neutrino alert and a blazar in a flaring state, TXS 0506+056; (2) the evidence of a neutrino flaring activity during 2014 - 2015 from the same direction; (3) the lack of an accompanying simultaneous $\gamma$-ray enhancement from the same counterpart; (4) the contrasting flaring activity of a neighbouring bright $\gamma$-ray source, the blazar PKS 0502+049, during 2014 - 2015. Our study makes use of multi-wavelength archival data accessed through Open Universe tools and includes a new analysis of Fermi-LAT data. We find that PKS 0502+049 contaminates the $\gamma$-ray emission region at low energies but TXS 0506+056 dominates the sky above a few GeV. TXS 0506+056, which is a very strong (top percent) radio and $\gamma$-ray source, is in a high $\gamma$-ray state during the neutrino alert but in a low though hard $\gamma$-ray state in coincidence with the neutrino flare. Both states can be reconciled with the energy associated with the neutrino emission and, in particular during the low/hard state, there is evidence that TXS 0506+056 has undergone a hadronic flare with very important implications for blazar modelling. All multi-messenger diagnostics reported here support a single coherent picture in which TXS 0506+056, a very high energy $\gamma$-ray blazar, is the only counterpart of all the neutrino emissions in the region and therefore the most plausible first non-stellar neutrino and, hence, cosmic ray source.

From the early days in gamma-ray astronomy, locating the origin of GeV emission within the core of an active galactic nucleus (AGN) persisted as an open question; the problem is to discern between near- and far-site scenarios with respect to the distance from the super massive central engine. We investigate this question under the light of a complete sample of low synchrotron peak (LSP) blazars which is fully characterized along many decades in the electromagnetic spectrum, from radio up to tens of GeV. We consider the high-energy emission from bright radio blazars and test for synchrotron self-Compton (SSC) and external Compton (EC) scenarios in the framework of localizing the $\gamma$-ray emission sites. Given that the inverse Compton (IC) process under the EC regime is driven by the abundance of external seed photons, these photons could be mainly ultraviolet (UV) to X-rays coming from the accretion disk region and the broad-line region (BLR), therefore close to the jet launch base; or infrared (IR) seed photons from the dust torus and molecular cloud spine-sheath, therefore far from jet launch base. We use enhanced SED information from recent works (including new $\gamma$-ray detections) to refine the study of Syn to IC peak correlations, which points to a particular $\gamma$-ray emission site. We show that SSC alone is not enough to account for the observed SEDs. Our analysis favors an EC scenario under the Thomson scattering regime, with a dominant IR external photon field. Therefore, the far-site (i.e., far from the jet launch) is probably the most reasonable scenario to account for the population properties of bright LSP blazars in cases modeled with a pure leptonic component. We calculate the photon energy density associated with the external field at the jet comoving frame to be $\rm U'_{ext}=1.69 \times 10^{-2}$ erg/cm$^3$, finding good agreement to other correlated works.

We present the 1-100 GeV spectral energy distribution for a population of 148 high-synchrotron-peaked blazars (HSPs) recently detected with Fermi-LAT as part of the First Brazil-ICRANet Gamma-ray Blazar catalogue (1BIGB). Most of the 1BIGB sources do not appear in previous Fermi-LAT catalogues and their gamma-ray spectral properties are presented here for the first time, representing a significant new extension of the gamma-ray blazar population. Since our sample was originally selected from an excess signal in the 0.3-500 GeV band, the sources stand out as promising TeV blazar candidates, potentially in reach of the forthcoming very-high-energy (VHE) gamma-ray observatory, CTA. The flux estimates presented here are derived considering PASS8 data, integrating over more than 9 years of Fermi-LAT observations. We also review the full broadband fit between 0.3-500 GeV presented in the original 1BIGB paper for all sources, updating the power-law parameters with currently available Fermi-LAT dataset. The importance of these sources in the context of VHE population studies with both current instruments and the future CTA is evaluated. To do so, we select a subsample of 1BIGB sources and extrapolate their gamma-ray SEDs to the highest energies, properly accounting for absorption due to the extragalactic background light. We compare those extrapolations to the published CTA sensitivity curves and estimate their detectability by CTA. Two notable sources from our sample, namely 1BIGB J224910.6-130002 and 1BIGB J194356.2+211821, are discussed in greater detail. All gamma-ray SEDs, which are shown here for the first time, are made publicly available via the Brazilian Science Data Center (BSDC) service, maintained at CBPF, in Rio de Janeiro.

We study the $\gamma$-ray and broadband spectral energy distribution (SED) properties of a complete sample of 104 bright, radio-selected low-synchrotron peaked (LSP) blazars. Most of the sources have already been detected in the $\gamma$-ray band by Fermi-LAT, however almost 20% of these blazars have no counterpart in any of the Fermi catalogs so far. Using the Fermi Science Tools, we look for $\gamma$-ray emission for those objects not yet reported in any Fermi-LAT catalog. We performed a binned likelihood analysis in the 0.3-500 GeV energy band with Fermi-LAT Pass 8 data, integrating over 7.5 years of observations. We studied $\gamma$-ray light curves and test statistic (TS) maps to validate new detections and associations. We fit the synchrotron (Syn) and inverse Compton (IC) components for all sources using all available historical data, enhancing the study of Syn to IC peak-power correlations. We derive the distribution of the Compton dominance (CD) along with population properties such as Syn and IC peak power, and frequency distributions. We deliver a unique characterization in $\gamma$-rays for a complete sample of LSP blazars. We show that three previously unidentified 3FGL sources can be associated with blazars when using improved $\gamma$-ray positions obtained from TS maps. Six previously unreported $\gamma$-ray sources are detected at TS>20 level, while another three show TS values between 10-20. We evaluate two cases in which source confusion is likely present. In four cases there is no significant $\gamma$-ray signature, however short-lived flares have been detected in these sources. Finally, the log(CD) has a Gaussian-like distribution with median log(CD)=0.1, implying that on average the peak-power for the Syn and IC components in LSP blazars are similar.

A direct search of $\gamma$-ray emission centered on multi-frequency selected candidates is a valuable complementary approach to the standard search adopted in current $\gamma$-ray Fermi-LAT catalogs. Our candidates are part of the 2WHSP sample, that was assembled with the aim of providing targets for Imaging Atmospheric Cherenkov Telescopes (IACT), and is currently the largest set of high synchrotron peaked (HSP) blazars. We perform a likelihood analysis with the Fermi Science Tools using positions from 400 2WHSP blazars as seeds of tentative $\gamma$-ray sources. This enabled us to detect 150 $\gamma$-ray excess signals that have not yet been reported in previous $\gamma$-ray catalogs (1FGL, 2FGL, and 3FGL). By identifying new sources, we solve a fraction of the extragalactic isotropic $\gamma$-ray background (IGRB) composition, improving the description of the $\gamma$-ray sky. Our analysis considers the 0.3-500 GeV energy band, integrating over 7.2 yrs of Fermi-LAT observation and making use of the Pass 8 data release. Within the 150 2WHSPs that showed excess $\gamma$-ray signature: 85 are high-significance detections with test statistic (TS)$>$25, and 65 are lower-significance detections with TS between 10 to 25. We study the $\gamma$-ray photon spectral index distribution, the likelihood of detection according to the synchrotron peak brightness, and plot the measured $\gamma$-ray LogN-LogS of HSP blazars, also discussing the portion of the IGRB that has been resolved by this work. We report on four cases where we could resolve source confusion and find counterparts for unassociated 3FGL sources with the help of high-energy TS maps together with multi-frequency data. The 150 new $\gamma$-ray sources are named with the acronym 1BIGB for the first version of the Brazil ICRANet Gamma-ray Blazar catalog, in reference to the cooperation agreement supporting this work.

Aims. High Synchrotron Peaked blazars (HSPs) dominate the -ray sky at energies larger than a few GeV; however, only a few hundred blazars of this type have been catalogued so far. In this paper we present the 2WHSP sample, the largest and most complete list of HSP blazars available to date, which is an expansion of the 1WHSP catalog of gamma-ray source candidates off the Galactic plane. Methods. We cross-matched a number of multi-wavelength surveys (in the radio, infrared and X-ray bands) and applied selection criteria based on the radio to IR and IR to X-ray spectral slopes. To ensure the selection of genuine HSPs we examined the SED of each candidate and estimated the peak frequency of its synchrotron emission ($\nu_{\rm peak}$) using the ASDC SED tool, including only sources with Synchrotron peak $\nu_{\rm peak} > 10^{15}$ Hz (equivalent to $\nu_{\rm peak} > 4$ eV). Results. We have assembled the largest and most complete catalog of HSP blazars to date, which includes 1691 sources. A number of population properties, such as infrared colours, synchrotron peak, redshift distributions, and -ray spectral properties, have been used to characterise the sample and maximize completeness. We also derived the radio logN-logS distribution. This catalog has already been used to provide seeds to discover new very high energy objects within Fermi-LAT data and to look for the counterparts of neutrino and ultra high energy cosmic ray sources, showing its potential for the identification of promising high-energy gamma-ray sources and multi-messenger targets.

We explore the correlation of $\gamma$-ray emitting blazars with IceCube neutrinos by using three very recently completed, and independently built, catalogues and the latest neutrino lists. We introduce a new observable, namely the number of neutrino events with at least one $\gamma$-ray counterpart, $N_{\nu}$. In all three catalogues we consistently observe a positive fluctuation of $N_{\nu}$ with respect to the mean random expectation at a significance level of $0.4 - 1.3$ per cent. This applies only to extreme blazars, namely strong, very high energy $\gamma$-ray sources of the high energy peaked type, and implies a model-independent fraction of the current IceCube signal $\sim 10 - 20$ per cent. An investigation of the hybrid photon -- neutrino spectral energy distributions of the most likely candidates reveals a set of $\approx 5$ such sources, which could be linked to the corresponding IceCube neutrinos. Other types of blazars, when testable, give null correlation results. Although we could not perform a similar correlation study for Galactic sources, we have also identified two (further) strong Galactic $\gamma$-ray sources as most probable counterparts of IceCube neutrinos through their hybrid spectral energy distributions. We have reasons to believe that our blazar results are not constrained by the $\gamma$-ray samples but by the neutrino statistics, which means that the detection of more astrophysical neutrinos could turn this first hint into a discovery.

Blazars are the dominant type of extragalactic sources at microwave and at $\gamma$-ray energies. In the most energetic part of the electromagnetic spectrum (E>100GeV) a large fraction of high Galactic latitude sources are blazars of the High Synchrotron Peaked (HSP) type, that is BL Lac objects with synchrotron power peaking in the UV or in the X-ray band. HSP blazars are remarkably rare, with only a few hundreds of them expected to be above the sensitivity limits of currently available surveys. To find these very uncommon objects, we have devised a method that combines ALLWISE survey data with multi-frequency selection criteria. The sample was defined starting from a primary list of infrared colour-colour selected sources from the ALLWISE all sky survey database, and applying further restrictions on IR-radio and IR-X-ray flux ratios. Using a polynomial fit to the multi-frequency data (radio to X-ray) we estimated synchrotron peak frequencies and fluxes of each object. We assembled a sample including 992 sources, which is currently the largest existing list of confirmed and candidates HSP blazars. All objects are expected to radiate up to the highest $\gamma$-ray photon energies. In fact, 299 of these are confirmed emitters of GeV $\gamma$-ray photons (based on Fermi-LAT catalogues), and 36 have already been detected in the TeV band. The majority of sources in the sample are within reach of the upcoming Cherenkov Telescope Array (CTA), and many may be detectable even by the current generation of Cherenkov telescopes during flaring episodes. The sample includes 425 previously known blazars, 151 new identifications, and 416 HSP candidates (mostly faint sources) for which no optical spectra is available yet. The full 1WHSP catalogue is on-line at http://www.asdc.asi.it/1whsp/ providing a direct link to the SED building tool where multifrequency data can be easily visualised.