With ten years of operation and an exceptional dataset, the Fermi-Large Area Telescope allows us to unveil the detailed composition of the extragalactic $\gamma$-ray sky above $\rm 100\,MeV$. In this paper, we derive the intrinsic source-count distribution (logN-logS) of extragalactic sources (i.e., blazars) at $|b|>20^{\circ}$ via the efficiency correction method. With this approach, we are able to measure the distribution down to a photon flux of $\sim10^{-10}$${\rm \, ph \,\, cm^{-2} \, s^{-1}}$ and to an energy flux of $\sim10^{-12}$${\rm \, erg \,\, cm^{-2} \, s^{-1}}$. In both cases the logN-logS becomes flatter at low fluxes. Moreover, we show that this logN-logS is representative of the blazar population (assuming the majority of unassociated sources are blazars) and allows us to constrain its evolution quite effectively. Among recently proposed evolutionary models, we find that the Pure Density Evolution (PDE) model best describes the evolutionary properties of the blazar population and that their integrated emission accounts for $\sim 50^{+10}_{-5}\%$ of the total extragalactic \gm-ray background.

High synchrotron peak (HSP; $\nu_{sy}^{pk} > 10^{15}$ Hz) BL Lacs are some of the most extreme accelerators in the Universe. Those found at high redshifts ($z>1$) challenge our understanding of blazar evolution models and are crucial for cosmological measurements of the Extragalactic Background Light. In this paper, we study a high-$z$ BL Lac, 4FGL J2146.5-1344, detected to be at $z$=1.34 using the photometric dropout technique. We collected multi-wavelength data for this source from optical up to $\gamma$-rays, in order to study its spectral energy distribution (SED). In particular, this source was observed for the first time with NuSTAR, which accurately measures the synchrotron emission of this blazar up to 50 keV. Despite being classified as an HSP BL Lac object, the modeling of the SED reveals that this source likely belongs to the "masquerading BL Lac" class, which comprises of FSRQs appearing as disguised BL Lac objects.

With bolometric luminosities exceeding $10^{48}$ erg s$^{-1}$, powerful jets and supermassive black holes at their center, MeV blazars are some of the most extreme sources in the Universe. Recently, the Fermi-Large Area Telescope detected five new $\gamma$-ray emitting MeV blazars beyond redshift $z=3.1$. With the goal of precisely characterizing the jet properties of these extreme sources, we started a multiwavelength campaign to follow them up with joint NuSTAR, Swift and SARA observations. We observe six high-redshift quasars, four of them belonging to the new $\gamma$-ray emitting MeV blazars. Thorough X-ray analysis reveals spectral flattening at soft X-ray for three of these objects. The source NVSS J151002$+$570243 also shows a peculiar re-hardening of the X-ray spectrum at energies $E>6\,\rm keV$. Adopting a one-zone leptonic emission model, this combination of hard X-rays and $\gamma$-rays enables us to determine the location of the Inverse Compton peak and to accurately constrain the jet characteristics. In the context of the jet-accretion disk connection, we find that all six sources have jet powers exceeding accretion disk luminosity, seemingly validating this positive correlation even beyond $z>3$. Our six sources are found to have $10^9 \rm M_{\odot}$ black holes, further raising the space density of supermassive black holes in the redshift bin $z=[3,4]$.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

The Hubble constant $H_{0}$ and matter density $\Omega_{m}$ of the Universe are measured using the latest $\gamma$-ray attenuation results from Fermi-LAT and Cherenkov telescopes. This methodology is based upon the fact that the extragalactic background light supplies opacity for very high energy photons via photon-photon interaction. The amount of $\gamma$-ray attenuation along the line of sight depends on the expansion rate and matter content of the Universe. This novel strategy results in a value of $H_{0}=67.4_{-6.2}^{+6.0}$~km~s$^{-1}$~Mpc$^{-1}$ and $\Omega_{m}=0.14_{-0.07}^{+0.06}$. These estimates are independent and complementary to those based on the distance ladder, cosmic microwave background (CMB), clustering with weak lensing, and strong lensing data. We also produce a joint likelihood analysis of our results from $\gamma$ rays and these from more mature methodologies, excluding the CMB, yielding a combined value of $H_{0}=66.6\pm 1.6$~km~s$^{-1}$~Mpc$^{-1}$ and $\Omega_{m}=0.29\pm 0.02$.

The recently released 105-month {\it Swift}-Burst Alert Telescope (BAT) all-sky hard X-ray survey catalog presents an opportunity to study astrophysical objects detected in the deepest look at the entire hard X-ray (14$-$195 keV) sky. Here we report the results of a multifrequency study of 146 blazars from this catalog, quadrupling the number compared to past studies, by utilizing recent data from the {\it Fermi}-Large Area Telescope (LAT), Swift-BAT, and archival measurements. In our $\gamma$-ray analysis of $\sim$10 years of the LAT data, 101 are found as $\gamma$-ray emitters, whereas, 45 remains LAT undetected. We model the broadband spectral energy distributions with a synchrotron-inverse Compton radiative model. On average, BAT detected sources host massive black holes ($M_{\rm bh}\sim10^9$ M$_{\odot}$) and luminous accretion disks ($L_{\rm d}\sim10^{46}$ erg s$^{-1}$). At high-redshifts ($z>2$), BAT blazars host more powerful jets with luminous accretion disks compared to those detected only with the {\it Fermi}-LAT. We find good agreement in the black hole masses derived from the single-epoch optical spectroscopic measurements and standard accretion disk modeling approaches. Other physical properties of BAT blazars are similar to those known for {\it Fermi}-LAT detected objects.

MeV blazars are the most luminous persistent sources in the Universe and emit most of their energy in the MeV band. These objects display very large jet powers and accretion luminosities and are known to host black holes with a mass often exceeding $10^9 M_{\odot}$. An MeV survey, performed by a new generation MeV telescope which will bridge the entire energy and sensitivity gap between the current generation of hard X-ray and gamma-ray instruments, will detect $>$1000 MeV blazars up to a redshift of $z=5-6$. Here we show that this would allow us: 1) to probe the formation and growth mechanisms of supermassive black holes at high redshifts, 2) to pinpoint the location of the emission region in powerful blazars, 3) to determine how accretion and black hole spin interplay to power the jet.

We present the analysis of a sample of 35 candidate Compton thick (CT-) active galactic nuclei (AGNs) selected in the nearby Universe (average redshift <z>~0.03) with the Swift-BAT 100-month survey. All sources have available NuSTAR data, thus allowing us to constrain with unprecedented quality important spectral parameters such as the obscuring torus line-of-sight column density (N_{H, z}), the average torus column density (N_{H, tor}) and the torus covering factor (f_c). We compare the best-fit results obtained with the widely used MyTorus (Murphy et al. 2009) model with those of the recently published borus02 model (Balokovic et al. 2018) used in the same geometrical configuration of MyTorus (i.e., with f_c=0.5). We find a remarkable agreement between the two, although with increasing dispersion in N_{H, z} moving towards higher column densities. We then use borus02 to measure f_c. High-f_c sources have, on average, smaller offset between N_{H, z} and N_{H, tor} than low-f_c ones. Therefore, low f_c values can be linked to a "patchy torus" scenario, where the AGN is seen through an over-dense region in the torus, while high-f_c objects are more likely to be obscured by a more uniform gas distribution. Finally, we find potential evidence of an inverse trend between f_c and the AGN 2-10 keV luminosity, i.e., sources with higher f_c values have on average lower luminosities.

High-redshift ($z>2$) blazars are the most powerful members of the blazar family. Yet, only a handful of them have both X-ray and $\gamma$-ray detection, thereby making it difficult to characterize the energetics of the most luminous jets. Here, we report, for the first time, the Fermi-Large Area Telescope detection of the significant $\gamma$-ray emission from the high-redshift blazar DA 193 ($z=2.363$). Its time-averaged $\gamma$-ray spectrum is soft ($\gamma$-ray photon index = $2.9\pm0.1$) and together with a relatively flat hard X-ray spectrum (14$-$195 keV photon index = $1.5\pm0.4$), DA 193 presents a case to study a typical high-redshift blazar with inverse Compton peak being located at MeV energies. An intense GeV flare was observed from this object in the first week of 2018 January, a phenomenon rarely observed from high-redshift sources. What makes this event a rare one is the observation of an extremely hard $\gamma$-ray spectrum (photon index = $1.7\pm0.2$), which is somewhat unexpected since high-redshift blazars typically exhibit a steep falling spectrum at GeV energies. The results of our multi-frequency campaign, including both space- (Fermi, NuSTAR, and Swift) and ground-based (Steward and Nordic Optical Telescope) observatories, are presented and this peculiar $\gamma$-ray flare is studied within the framework of a single-zone leptonic emission scenario.

We present the combined NuSTATR and XMM-Newton 0.6-79 keV spectral analysis of a Seyfert 2 galaxy, NGC 1358, which we selected as a candidate Compton thick (CT-) active galactic nucleus (AGN) on the basis of previous Swift/BAT and Chandra studies. According to our analysis, NGC 1358 is confirmed to be a CT-AGN using physical motivated models, at >3 $\sigma$ confidence level. Our best-fit shows that the column density along the 'line-of-sight' of the obscuring material surrounding the accreting super-massive black hole is N$\rm _H$ = [1.96--2.80] $\times$ 10$^{24}$ cm$^{-2}$. The high-quality data from NuSTAR gives the best constraints on the spectral shape above $\sim$10 keV to date on NGC 1358. Moreover, by combining NuSTAR and XMM-Newton data, we find that the obscuring torus has a low covering factor ($f_c$ <0.17), and the obscuring material is distributed in clumps, rather than uniformly. We also derive an estimate of NGC 1358's Eddington ratio, finding it to be $\lambda_{\rm Edd}$ $\sim$$4.7_{-0.3}^{+0.3}$ $\times$ 10$^{-2}$, which is in acceptable agreement with previous measurements. Finally, we find no evidence of short-term variability, over a $\sim$100 ks time-span, in terms of both 'line-of-sight' column density and flux.

We present the 2-100 keV spectral analysis of 30 candidate Compton thick (CT-) active galactic nuclei (AGN) selected in the Swift-BAT 100-month survey. The average redshift of these objects is $\langle z\rangle\sim$0.03 and they all lie within $\sim$500 Mpc. We used the MyTorus (Murphy et al. 2009) model to perform X-ray spectral fitting both without and with the contribution of the NuSTAR data in the 3-50 keV energy range. When the NuSTAR data are added to the fit, 14 out of 30 of these objects (47% of the whole sample) have intrinsic absorption N$_{\rm H}<10^{24}$ cm$^{-2}$ at the $>$3$\sigma$ confidence level, i.e., they are re-classified from Compton thick to Compton thin. Consequently, we infer an overall observed fraction of CT-AGN with respect to the whole AGN population lower than the one reported in previous works, and as low as $\sim$4%. We find evidence that this over-estimation of N$_{\rm H}$ is likely due to the low quality of a subsample of spectra, either in the 2-10 keV band or in the Swift-BAT one.

We present the results of a multi-frequency, time-averaged analysis of blazars included in the Candidate Gamma-ray Blazar Survey catalog. Our sample consists of 324 $\gamma$-ray detected ($\gamma$-ray loud) and 191 non $\gamma$-ray detected ($\gamma$-ray quiet) blazars, and we consider all the data up to 2016 April 1. We find that both the $\gamma$-ray loud and the $\gamma$-ray quiet blazar populations occupy similar regions in the WISE color-color diagram, and in the radio and X-ray bands $\gamma$-ray loud sources are brighter. A simple one-zone synchrotron inverse-Compton emission model is applied to derive the physical properties of both populations. We find that the central black hole mass and the accretion disk luminosity ($L_{\rm disk}$) computed from the modeling of the optical-UV emission with a Shakura-Sunyaev disk reasonably matches with that estimated from the optical spectroscopic emission-line information. A significantly larger Doppler boosting in the $\gamma$-ray loud blazars is noted, and their jets are more radiatively efficient. On the other hand, the $\gamma$-ray quiet objects are more MeV-peaked, thus could be potential targets for next-generation MeV missions. Our results confirm the earlier findings about the accretion-jet connection in blazars; however, many of the $\gamma$-ray quiet blazars tend to deviate from the recent claim that the jet power exceeds $L_{\rm disk}$ in blazars. A broadband study, considering a larger set of $\gamma$-ray quiet objects and also including BL Lacs, will be needed to confirm/reject this hypothesis and also to verify the evolution of the powerful high-redshift blazars into their low-power nearby counterparts.

Based on current models of the Cosmic X-ray Background (CXB), heavily obscured Active Galactic Nuclei (AGN) are expected to make up ~10% of the peak emission of the CXB and ~20% of the total population of AGN, yet few of these sources have been recorded and characterized in current surveys. Here we present the Chandra follow-up observation of 14 AGN detected by Swift-BAT. For five sources in the sample, NuSTAR observations in the 3-80 keV band are also available. The X-ray spectral fitting over the 0.3-150 keV energy range allows us to determine the main X-ray spectral parameters, such as the photon index and the intrinsic absorption, of these objects, and to make hypotheses on the physical structures responsible for the observed spectra. We find that 13 of the 14 objects are absorbed AGN, and one is a candidate Compton thick AGN, having intrinsic absorption NH>1E24 cm^{-2}. Finally, we verified that the use of NuSTAR observations is strategic to strongly constrain the properties of obscured AGN, since the best-fit values we obtained for parameters such as the power-law photon index and the intrinsic absorption NH changed sometimes significantly fitting the spectra with and without the use of NuSTAR data.

The most powerful sources among the blazar family are MeV blazars. Often detected at $z>2$, they usually display high X- and \gm-ray luminosities, larger-than-average jet powers and black hole masses $\gtrsim 10^9 M_{\odot}$. In the present work we perform a multiwavelength study of three high redshift blazars: 3FGL J0325.5+2223 ($z=2.06$), 3FGL J0449.0+1121 ($z= 2.15$), and 3FGL J0453.2$-$2808 ($z=2.56$), analysing quasi simultaneous data from GROND, \swift-UVOT and XRT, \nustar, and \fermi-LAT. Our main focus is on the hard X-ray band recently unveiled by \nustar~(3$-$79 keV) where these objects show a hard spectrum which enables us to constrain the inverse Compton peak and the jet power. We found that all three targets resemble the most powerful blazars, with the synchrotron peak located in the sub-millimeter range and the inverse Compton peak in the MeV range, and therefore belong to the MeV blazar class. Using a simple one zone leptonic emission model to reproduce the spectral energy distributions, we conclude that a simple combination of synchrotron and accretion disk emission reproduces the infrared-optical spectra while the X-ray to \gm-ray part is well reproduced by the inverse Compton scattering of low energy photons supplied by the broad line region. The black hole masses for each of the three sources are calculated to be $\gtrsim 4 \times 10^{8} M_{\odot}$. The three studied sources have jet power at the level of, or beyond, the accretion luminosity.

MeV blazars are a sub--population of the blazar family, exhibiting larger--than--average jet powers, accretion luminosities and black hole masses. Because of their extremely hard X--ray continua, these objects are best studied in the X-ray domain. Here, we report on the discovery by the $Fermi$ Large Area Telescope and subsequent follow-up observations with $NuSTAR$, $Swift$ and GROND of a new member of the MeV blazar family: PMN J0641$-$0320. Our optical spectroscopy provides confirmation that this is a flat--spectrum radio quasar located at a redshift of $z=1.196$. Its very hard $NuSTAR$ spectrum (power--law photon index of $\sim$1 up to $\sim$80 keV) indicates that the emission is produced via inverse Compton scattering off photons coming from outside the jet.The overall spectral energy distribution of PMN J0641$-$0320 is typical of powerful blazars and by reproducing it with a simple one-zone leptonic emission model we find the emission region to be located either inside the broad line region or within the dusty torus.