Search

We present a comprehensive multi-messenger study of NGC 1068, the prototype Seyfert II galaxy recently associated with high-energy IceCube neutrinos. Various aspects of the source, including its nuclear activity, jet, outflow, and starburst region, are analyzed in detail using a multi-wavelength approach and relevant luminosities are derived. We then explore its gamma-ray and neutrino emissions and investigate potential mechanisms underlying these phenomena and their relations with the different astrophysical components to try to understand which one is responsible for the IceCube neutrinos. By first using simple order-of-magnitude arguments and then applying specific theoretical models, we infer that only the region close to the accretion disc around the supermassive black hole has both the right density of X-ray photons needed to provide the targets for protons to sustain neutrino production and of optical/infrared photons required to absorb the associated but unobserved gamma rays. We conclude by highlighting ongoing efforts to constrain a possible broad connection between neutrinos and active galactic nuclei, as well as future synergies between astronomical and neutrino facilities., Comment: 18 pages, 4 figures, review paper accepted for publication in Nature Astronomy; authors' version before editorial review. Revised version matches the title of the published paper; no other changes were made to the text

Recent results from numerical simulations and models of galaxy formation suggest that recently discovered ultra-faint compact stellar systems (UFCSs) in the halo of the Milky Way (MW) may be some of the smallest and faintest galaxies. If this is the case, these systems would be attractive targets for indirect searches of weakly interacting massive particle (WIMP) dark matter (DM) annihilation due to their relative proximity and high expected DM content. In this study, we analyze 14.3 years of gamma-ray data collected by the Fermi-LAT coincident with 26 UFCSs. No significant excess gamma-ray emission is detected, and we present gamma-ray flux upper limits for these systems. Assuming that the UFCSs are dark-matter-dominated galaxies consistent with being among the faintest and least massive MW dwarf spheroidal (dSphs) satellite galaxies, we derive the projected sensitivity for a dark matter annihilation signal. We find that observations of UFCSs have the potential to yield some of the most powerful constraints on DM annihilation, with sensitivity comparable to observations of known dSphs and the Galactic center. This result emphasizes the importance of precise kinematic studies of UFCSs to empirically determine their DM content., Comment: 11 pages, 4 figures, 1 table, submitted to ApJL

Active galactic nuclei (AGNs) can power relativistic jets, which are called blazars when pointed close to our line of sight. Depending on the presence or absence of emission lines in their optical spectra, blazars are categorized into flat spectrum radio quasars (FSRQs) or BL Lacertae (BL Lacs) objects. According to the 'blazar sequence', as synchrotron peak frequency ($\nu^{sy}_{pk}$) shifts to higher energies, the synchrotron peak luminosity decreases. This means that BL Lacs as luminous as FSRQs, and with synchrotron peak frequencies $\nu^{sy}_{pk}>10^{15}$ Hz, should not exist. Detected as a high-synchrotron peak (HSP; $\nu^{sy}_{pk}>10^{15}$ Hz) BL Lac, 4FGL J1520.8-0348 shows high gamma-ray luminosity ($L_{\gamma}>10^{46}\,\rm erg~s^{-1}$), being at a high redshift of $z=$1.46. Since it is an outlier in the 'blazar sequence', the process of its jet acceleration and power may be different from bona fide BL Lacs. In this work, we constrain its spectral energy distribution (SED) by modeling the multi-wavelength data from infrared to $\gamma$-ray regime. Simultaneous X-ray data was obtained from X-ray Multi-Mirror Mission and Nuclear Spectroscopic Telescope Array to constrain the synchrotron emission and underlying electron distribution. On undertaking the SED modeling of the source, including the effect of extragalactic background light, we conclude that the source is more likely to be a 'blue FSRQ' or 'masquerading BL Lac' where the BL Lac is actually a FSRQ in disguise.

We present the multi-epoch analysis of 13 variable, nearby (z<0.1), Compton-thin (22

Blazars present highly variable $\gamma$-ray emission. This variability, which can range from a few minutes to several years, is also observed at other wavelengths across the entire electromagnetic spectrum. We make use of the first 12 years of data from the Fermi Large Area Telescope (LAT), complemented with multiwavelength (MWL) archival data from different observatories and facilities in radio, infrared and optical bands, to study the possible periodic emission from 19 blazars previously claimed as periodic candidates. A periodicity analysis is performed with a pipeline for periodicity searches. Moreover, we study the cross-correlations between the $\gamma$-ray and MWL light curves. Additionally, we use the fractional variability and the structure function to evaluate the variability timescales. We find five blazars showing hints of periodic modulation with $\geq$3.0$\sigma$ ($\approx$0$\sigma$ post-trials), with periods ranging from 1.2 to 4 years, both in their $\gamma$-ray and MWL emission. The results provide clues for understanding the physical mechanisms generating the observed periodicity., Comment: 21 pages, 5 figures

Blazars display variable emission across the entire electromagnetic spectrum, with timescales that can range from a few minutes to several years. Our recent work has shown that a sample of five blazars exhibit hints of periodicity with a global significance $\gtrsim2\,\sigma$ at $\gamma$-ray energies, in the range of 0.1~GeV$<$E$<$800~GeV. In this work, we study their multiwavelength (MWL) emission, covering the X-ray, ultraviolet, optical, and radio bands. We show that three of these blazars present similar periodic patterns in the optical and radio bands. Additionally, fluxes in the different bands of the five blazars are correlated, suggesting a co-spatial origin. Moreover, we detect a long-term ($\approx$10 year) rising trend in the light curves of PG~1553+113, and we use it to infer possible constraints on the binary black hole hypothesis., Comment: 20 pages, 10 figures, 7 tables

The extragalactic high-energy $\gamma$-ray sky is dominated by blazars, which are active galactic nuclei with their jets pointing towards us. Distance measurements are of fundamental importance yet for some of these sources are challenging because any spectral signature from the host galaxy may be outshone by the non-thermal emission from the jet. In this paper, we present a method to constrain redshifts for these sources that relies only on data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. This method takes advantage of the signatures that the pair-production interaction between photons with energies larger than approximately 10 GeV and the extragalactic background light leaves on $\gamma$-ray spectra. We find upper limits for the distances of 303 $\gamma$-ray blazars, classified as 157 BL Lacertae objects, 145 of uncertain class, and 1 flat-spectrum-radio quasar, whose redshifts are otherwise unknown. These derivations can be useful for planning observations with imaging atmospheric Cherenkov telescopes and also for testing theories of supermassive black hole evolution. Our results are applied to estimate the detectability of these blazars with the future Cherenkov Telescope Array, finding that at least 21 of them could be studied in a reasonable exposure of 20 h., Comment: 13 pages, 7 figures, 1 table; Accepted by MNRAS

We present $\gamma$-ray optical-depth calculations from a recently published extragalactic background light (EBL) model built from multiwavelength galaxy data from the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (HST/CANDELS). CANDELS gathers one of the deepest and most complete observations of stellar and dust emissions in galaxies. This model resulted in a robust derivation of the evolving EBL spectral energy distribution up to $z\sim 6$, including the far-infrared peak. Therefore, the optical depths derived from this model will be useful for determining the attenuation of $\gamma$-ray photons coming from high-redshift sources, such as those detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope, and for multi-TeV photons that will be detected from nearby sources by the future Cherenkov Telescope Array. From these newly calculated optical depths, we derive the cosmic $\gamma$-ray horizon and also measure the expansion rate and matter content of the Universe including an assessment of the impact of the EBL uncertainties. We find $H_{0}=61.9$ $^{+2.9}_{-2.4}$ km s$^{-1}$ Mpc$^{-1}$ when fixing $\Omega_{m}=0.32$, and $H_{0}=65.6$ $^{+5.6}_{-5.0}$ km s$^{-1}$ Mpc$^{-1}$ and $\Omega_{m}=0.19\pm 0.07$, when exploring these two parameters simultaneously., Comment: 11 pages, 8 figures, 1 tables; Accepted by MNRAS

We report the discovery of GRB 221009A, the highest flux gamma-ray burst ever observed by the Fermi Gamma-ray Burst Monitor (GBM). This GRB has continuous prompt emission lasting more than 600 seconds which smoothly transitions to afterglow visible in the GBM energy range (8 keV--40 MeV), and total energetics higher than any other burst in the GBM sample. By using a variety of new and existing analysis techniques we probe the spectral and temporal evolution of GRB 221009A. We find no emission prior to the GBM trigger time (t0; 2022 October 9 at 13:16:59.99 UTC), indicating that this is the time of prompt emission onset. The triggering pulse exhibits distinct spectral and temporal properties suggestive of the thermal, photospheric emission of shock-breakout, with significant emission up to $\sim$15 MeV. We characterize the onset of external shock at t0+600 s and find evidence of a plateau region in the early-afterglow phase which transitions to a slope consistent with Swift-XRT afterglow measurements. We place the total energetics of GRB 221009A in context with the rest of the GBM sample and find that this GRB has the highest total isotropic-equivalent energy ($\textrm{E}_{\gamma,\textrm{iso}}=1.0\times10^{55}$ erg) and second highest isotropic-equivalent luminosity ($\textrm{L}_{\gamma,\textrm{iso}}=9.9\times10^{53}$ erg/s) based on redshift of z = 0.151. These extreme energetics are what allowed us to observe the continuously emitting central engine of GBM from the beginning of the prompt emission phase through the onset of early afterglow., Comment: 26 pages 7 figures - accepted for publication in ApJL

Including millimeter-wave (mm-wave) data in multi-wavelength studies of the variability of active galactic nuclei (AGN) can provide insights into AGN physics that are not easily accessible at other wavelengths. We demonstrate in this work the potential of cosmic microwave background (CMB) telescopes to provide long-term, high-cadence mm-wave AGN monitoring over large fractions of sky. We report on a pilot study using data from the SPTpol instrument on the South Pole Telescope (SPT), which was designed to observe the CMB at arcminute and larger angular scales. Between 2013 and 2016, SPTpol was used primarily to observe a single 500 deg^2 field, covering the entire field several times per day with detectors sensitive to radiation in bands centered at 95 and 150 GHz. We use SPT 150 GHz observations to create AGN light curves, and we compare these mm-wave light curves to those at other wavelengths, in particular gamma-ray and optical. In this Letter, we focus on a single source, PKS 2326-502, which has extensive, day-timescale monitoring data in gamma-ray, optical, and now mm-wave between 2013 and 2016. We find PKS 2326-502 to be in a flaring state in the first two years of this monitoring, and we present a search for evidence of correlated variability between mm-wave, optical R band, and gamma-ray observations. This pilot study is paving the way for AGN monitoring with current and upcoming CMB experiments such as SPT-3G, Simons Observatory, and CMB-S4, including multi-wavelength studies with facilities such as VRO-LSST., Comment: 9 pages, 3 figures, accepted to Astrophysical Journal Letters

We present the analysis of multiepoch observations of a set of 12 variable, Compton-thin, local (z<0.1) active galactic nuclei (AGN) selected from the 100-month BAT catalog. We analyze all available X-ray data from \chandra, \xmm, and \nustar, adding up to a total of 53 individual observations. This corresponds to between 3 and 7 observations per source, probing variability timescales between a few days and $\sim 20$~yr. All sources have at least one \nustar observation, ensuring high-energy coverage, which allows us to disentangle the line-of-sight and reflection components in the X-ray spectra. For each source, we model all available spectra simultaneously, using the physical torus models \myt, \bor, and \uxc. The simultaneous fitting, along with the high-energy coverage, allows us to place tight constraints on torus parameters such as the torus covering factor, inclination angle, and average column density. We also estimate the line-of-sight column density ($N_{\rm H}$) for each individual observation. Within the 12 sources, we detect clear line-of-sight $N_{\rm H}$ variability in 5, non-variability in 5, and for 2 of them it is not possible to fully disentangle intrinsic-luminosity and $N_{\rm H}$ variability. We observe large differences between the average values of line-of-sight $N_{\rm H}$ (or $N_{\rm H}$ of the obscurer) and the average $N_{\rm H}$ of the torus (or $N_{\rm H}$ of the reflector), for each source, by a factor between $\sim2$ to $>100$. This behavior, which suggests a physical disconnect between the absorber and the reflector, is more extreme in sources that present $N_{\rm H}$ variability. $N_{\rm H}$-variable AGN also tend to present larger obscuration and broader cloud distributions than their non-variable counterparts. We observe that large changes in obscuration only occur at long timescales, and use this to place tentative lower limits on torus cloud sizes., Comment: 16 pages, 6 figures, 3 tables, 3 Appendix (19 pages, 2 figures, 11 tables). Submitted to A&A

The Fermi Large Area Telescope (LAT) light curve repository (LCR) is a publicly available, continually updated library of gamma-ray light curves of variable Fermi-LAT sources generated over multiple timescales. The Fermi-LAT LCR aims to provide publication-quality light curves binned on timescales of 3 days, 7 days, and 30 days for 1525 sources deemed variable in the source catalog of the first 10 years of Fermi-LAT observations. The repository consists of light curves generated through full likelihood analyses that model the sources and the surrounding region, providing fluxes and photon indices for each time bin. The LCR is intended as a resource for the time-domain and multi-messenger communities by allowing users to quickly search LAT data to identify correlated variability and flaring emission episodes from gamma-ray sources. We describe the sample selection and analysis employed by the LCR and provide an overview of the associated data access portal., Comment: Accepted for publication in ApJ Supplement Series

Active galactic Nuclei (AGNs) with their relativistic jets pointed toward the observer, are a class of luminous gamma-ray sources commonly known as blazars. The study of this source class is essential to unveil the physical processes powering these extreme jets, to understand their cosmic evolution, as well as to indirectly probe the extragalactic background light. To do so, however, one needs to correctly classify and measure a redshift for a large sample of these sources. The Third Fermi-LAT Catalog of High-Energy Sources (3FHL) contains $1212$ blazars detected at energies greater than 10 GeV. However, $\sim$25% of these sources are unclassified and $\sim$56% lack redshift information. To increase the optical completeness of blazars in the 3FHL catalog, we devised an optical spectroscopic follow-up campaign using 4 $m$ and 8 $m$ telescopes. In this paper, we present the results of the last part of this campaign, where we observed 5 blazars using the 8 $m$ Gemini-S telescope in Chile. We report all 5 sources to be classified as BL Lacs, a redshift lower limit for 2 sources, and featureless spectra for the remaining 3 sources. We also performed a one-zone leptonic fit to the two sources with the redshift lower limits., Comment: 10 pages, 3 Figures, 5 Tables, accepted for publication in The Astronomical Journal

The Third Catalog of Hard Fermi Large Area Telescope Sources (3FHL) reports the detection of 1556 objects at E > 10 GeV. However, 177 sources remain unassociated and 23 are associated with a ROSAT X-ray detection of unknown origin. Pointed X-ray observations were conducted on 30 of these unassociated and unknown sources with Swift-XRT. A bright X-ray source counterpart was detected in 21 out of 30 fields. In five of these 21 fields, we detected more than one X-ray counterpart, totaling 26 X-ray sources analyzed. Multiwavelength data was compiled for each X-ray source detected. We find that 21 out of the 26 X-ray sources detected display the multiwavelength properties of blazars, while one X-ray source displays the characteristics of a Galactic source. Using trained decision tree, random forest, and support vector machine models, we predict all 21 blazar counterpart candidates to be BL Lacertae objects (BL Lacs). This is in agreement with BL Lacs being the most populous source class in the 3FHL., Comment: 15 pages, 11 figures, interactive plot link see, https://authortools.aas.org/AAS40568/figNint.html

It is well known that blazars can show variability on a wide range of time scales. This behavior can include periodicity in their $\gamma$-ray emission, whose clear detection remains an ongoing challenge, partly due to the inherent stochasticity of the processes involved and also the lack of adequately-well sampled light curves. We report on a systematic search for periodicity in a sample of 24 blazars, using twelve well-established methods applied to Fermi-LAT data. The sample comprises the most promising candidates selected from a previous study, extending the light curves from nine to twelve years and broadening the energy range analyzed from $>$1 GeV to $>$0.1 GeV. These improvements allow us to build a sample of blazars that display a period detected with global significance $\gtrsim3\,\sigma$. Specifically, these sources are PKS 0454$-$234, S5 0716+714, OJ 014, PG 1553+113, and PKS 2155$-$304. Periodic $\gamma$-ray emission may be an indication of modulation of the jet power, particle composition, or geometry but most likely originates in the accretion disk, possibly indicating the presence of a supermassive black hole binary system., Comment: 26 pages, 8 figures, 7 tables

We derive the most up-to-date Swift-Burst Alert Telescope (BAT) blazar luminosity function in the 14-195 keV range, making use of a clean sample of 118 blazars detected in the BAT 105-month survey catalog, with newly obtained redshifts from the BAT AGN Spectroscopic Survey (BASS). We determine the best-fit X-ray luminosity function for the whole blazar population, as well as for Flat Spectrum Radio Quasars (FSRQs) alone. The main results are: (1) at any redshift, BAT detects the most luminous blazars, above any possible break in their luminosity distribution, which means we cannot differentiate between density and luminosity evolution; (2) the whole blazar population, dominated by FSRQs, evolves positively up to redshift z~4.3, confirming earlier results and implying lower number densities of blazars at higher redshifts than previously estimated. The contribution of this source class to the Cosmic X-ray Background at 14-195 keV can range from 5-18%, while possibly accounting for 100% of the MeV background. We also derived the average 14 keV-10 GeV SED for BAT blazars, which allows us to predict the number counts of sources in the MeV range, as well as the expected number of high-energy (>100 TeV) neutrinos. A mission like COSI, will detect 40 MeV blazars and 2 coincident neutrinos. Finally, taking into account beaming selection effects, the distribution and properties of the parent population of these extragalactic jets are derived. We find that the distribution of viewing angles is quite narrow, with most sources aligned within < 5{\deg} of the line of sight. Moreover, the average Lorentz factor, <$\Gamma$>= 8-12, is lower than previously suggested for these powerful sources., Comment: Accepted for publication in the Astrophysical Journal; 33 pages; 8 Tables; 16 Figures

BL Lacertae (BL Lac) objects are a subclass of blazar, distinguished by their featureless optical spectrum. The featureless spectrum presents a challenge in measuring the redshift of the BL Lacs. In this paper, we measure the redshift of BL Lacs using the photometric dropout technique. The space-based telescope \emph{Swift} and the ground-based SARA telescopes are employed to provide magnitudes in the $uvw2,\ uvm2,\ uvw1,\ u,\ b,\ v,\ g',\ r',\ i',\ z'$ filters. We observe 60 BL Lacs and report reliable redshift upper limits for 41 of them. We discover three new high-$z$ BL Lacs ($z>1.3$) at $1.74_{-0.08}^{+0.05}$, $1.88_{-0.03}^{+0.07}$, and $2.10_{-0.04}^{+0.03}$, bringing the number of high-$z$ BL Lacs found by this method up to 19. Discussions are made on the implications for the blazar sequence, the \emph{Fermi} blazar divide, and the gamma-ray horizon based on an analysis of the 4LAC catalog and all high-$z$ BL Lacs found with the photo-$z$ technique., Comment: 16pages, 11 figures

After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95\% credible limit on the GWB characteristic strain of $1.0\times10^{-14}$ at 1 yr$^{-1}$, which scales as the observing time span $t_{\mathrm{obs}}^{-13/6}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models., Comment: 3 figures in the main text. 3 figures and 8 tables are in the supplementary material

We present the analysis of simultaneous NuSTAR and XMM-Newton data of 8 Compton-thick (CT-) active galactic nuclei (AGN) candidates selected in the Swift-Burst Alert Telescope (BAT) 100 month survey. This work is part of an ongoing effort to find and characterize all CT-AGN in the local ($z\leq$0.05) Universe. We used two physically motivated models, MYTorus and borus02, to characterize the sources in the sample, finding 5 of them to be confirmed CT-AGN. These results represent an increase of $\sim19$% over the previous NuSTAR-confirmed, BAT-selected CT-AGN at $z\leq0.05$, bringing the total number to 32. This corresponds to an observed fraction of $\sim 8$\% of all AGN within this volume-limited sample, although it increases to $20\pm5$% when limiting the sample to $z\leq0.01$. Out of a sample of 48 CT-AGN candidates, selected using BAT and soft (0.3$-$10 keV) X-ray data, only 24 are confirmed as CT-AGN with the addition of the NuSTAR data. This highlights the importance of NuSTAR when classifying local obscured AGN. We also note that most of the sources in our full sample of 48 Seyfert 2 galaxies with NuSTAR data have significantly different line-of-sight and average torus column densities, favouring a patchy torus scenario., Comment: 31 pages, 13 figures, accepted for publication in ApJ

The Large Area Telescope (LAT), the primary instrument for the Fermi Gamma-ray Space Telescope (Fermi) mission, is an imaging, wide field-of-view, high-energy gamma-ray telescope, covering the energy range from 30 MeV to more than 300 GeV. We describe the performance of the instrument at the 10-year milestone. LAT performance remains well within the specifications defined during the planning phase, validating the design choices and supporting the compelling case to extend the duration of the Fermi mission. The details provided here will be useful when designing the next generation of high-energy gamma-ray observatories., Comment: 60 pages, 28 figures. Published in ApJS

As a follow-up of the optical spectroscopic campaign aimed at achieving completeness in the Third Catalog of Hard Fermi-LAT Sources (3FHL), we present here the results of a sample of 28 blazars of uncertain type observed using the 4m telescope at Cerro Tololo Inter-American Observatory (CTIO) in Chile. Out of these 28 sources, we find that 25 are BL Lacertae objects (BL Lacs) and 3 are Flat Spectrum Radio Quasars (FSRQs). We measure redshifts or lower limits for 16 of these blazars, whereas it is observed that the 12 remaining blazars have featureless optical spectra. These results are part of a more extended campaign of optical spectroscopy follow-up of 3FHL blazars, where until now 51 blazars of uncertain type have been classified into BL Lac and FSRQ categories. Further, this campaign has resulted in redshift measurements and lower limits for 15 of these sources. Our results contribute towards attaining a complete sample of blazars above 10 GeV, which then will be crucial in extending our knowledge on blazar emission mechanisms and the extragalactic background light.

We present a catalog of central engine properties, i.e., black hole mass (mbh) and accretion luminosity (ld), for a sample of 1077 blazars detected with the Fermi~Large Area Telescope. This includes broad emission line systems and blazars whose optical spectra lack emission lines but dominated by the absorption features arising from the host galaxy. The average mbh for the sample is $\langle \log~M_{{\rm BH,all}~msun} \rangle=8.60$ and there are evidences suggesting the association of more massive black holes with absorption line systems. Our results indicate a bi-modality of ld~in Eddington units (ld/ledd) with broad line objects tend to have a higher accretion rate (ld/ledd$>$0.01). We have found that ld/ledd and Compton dominance (CD, the ratio of the inverse Compton to synchrotron peak luminosities) are positively correlated at $>$5$\sigma$ confidence level, suggesting that the latter can be used to reveal the state of accretion in blazars. Based on this result, we propose a CD based classification scheme. Sources with CD$>$1 can be classified as High-Compton Dominated or HCD blazars, whereas, that with CD$\lesssim$1 are Low-Compton Dominated (LCD) objects. This scheme is analogous to that based on the mass accretion rate proposed in the literature, however, it overcomes the limitation imposed by the difficulty in measuring ld and mbh for objects with quasi-featureless spectra. We conclude that the overall physical properties of Fermi blazars are likely to be controlled by the accretion rate in Eddington units. The catalog is made public at http://www.ucm.es/blazars/engines and Zenodo., Comment: to appear in the ApJ Supplement series

We present the first Fermi - Large Area Telescope (LAT) solar flare catalog covering the 24 th solar cycle. This catalog contains 45 Fermi -LAT solar flares (FLSFs) with emission in the gamma-ray energy band (30 MeV - 10 GeV) detected with a significance greater than 5 sigma over the years 2010-2018. A subsample containing 37 of these flares exhibit delayed emission beyond the prompt-impulsive hard X-ray phase with 21 flares showing delayed emission lasting more than two hours. No prompt-impulsive emission is detected in four of these flares. We also present in this catalog the observations of GeV emission from 3 flares originating from Active Regions located behind the limb (BTL) of the visible solar disk. We report the light curves, spectra, best proton index and localization (when possible) for all the FLSFs. The gamma-ray spectra is consistent with the decay of pions produced by >300 MeV protons. This work contains the largest sample of high-energy gamma-ray flares ever reported and provides the unique opportunity to perform population studies on the different phases of the flare and thus allowing to open a new window in solar physics., Comment: Published in ApJS

The obscuration observed in active galactic nuclei (AGN) is mainly caused by dust and gas distributed in a torus-like structure surrounding the supermassive black hole (SMBH). However, properties of the obscuring torus of the AGN in X-ray have not been fully investigated yet due to the lack of high-quality data and proper models. In this work, we perform a broadband X-ray spectral analysis of a large, unbiased sample of obscured AGN (with line-of-sight column density 23$\le$log(NH)$\le$24) in the nearby universe which has high-quality archival NuSTAR data. The source spectra are analyzed using the recently developed borus02 model, which enables us to accurately characterize the physical and geometrical properties of AGN obscuring tori. We also compare our results obtained from the unbiased Compton thin AGN with those of Compton-thick AGN. We find that Compton thin and Compton-thick AGN may possess similar tori, whose average column density is Compton thick (N$\rm _{H,tor,ave}$ $\sim$1.4$\times$10$^{24}$ cm$^{-2}$), but they are observed through different (under-dense or over-dense) regions of the tori. We also find that the obscuring torus medium is significantly inhomogeneous, with the torus average column densities significantly different from their line-of-sight column densities (for most of the sources in the sample). The average torus covering factor of sources in our unbiased sample is c$_f$=0.67, suggesting that the fraction of unobscured AGN is $\sim$33%. We develop a new method to measure the intrinsic line-of-sight column density distribution of AGN in the nearby universe, which we find the result is in good agreement with the constraints from recent population synthesis models., Comment: 16 pages, 14 figures, 7 tables; accepted by A&A

One of the major challenges in studying the cosmic evolution of relativistic jets is the identification of the high-redshift ($z>3$) BL Lacertae objects, a class of jetted active galactic nuclei characterized by their quasi-featureless optical spectra. Here we report the identification of the first $\gamma$-ray emitting BL Lac object, 4FGL~J1219.0+3653 (J1219), beyond $z=3$, i.e., within the first two billion years of the age of the Universe. The optical and near-infrared spectra of J1219 taken from 10.4 m Gran Telescopio Canarias exhibit no emission lines down to an equivalent width of $\sim$3.5 A supporting its BL Lac nature. The detection of a strong Lyman-$\alpha$ break at $\sim$5570 A, on the other hand, confirms that J2119 is indeed a high-redshift ($z\sim3.59$) quasar. Based on the prediction of a recent BL Lac evolution model, J1219 is one of the only two such objects expected to be present within the comoving volume at $z=3.5$. Future identifications of more $z>3$ $\gamma$-ray emitting BL Lac sources, therefore, will be crucial to verify the theories of their cosmic evolution., Comment: Accepted for publication in The Astrophysical Journal Letters

We report on the investigation of a very high energy (VHE), Galactic gamma-ray source recently discovered at >50GeV using the Large Area Telescope (LAT) on board the Fermi Gamma-Ray Space Telescope. This object, 2FHL J1703.4-4145, displays a very hard >50GeV spectrum with a photon index ~1.2 in the 2FHL catalog and, as such, is one of the most extreme sources in the 2FHL sub-sample of Galactic objects. A detailed analysis of the available multi-wavelength data shows that this source is located on the western edge of the supernova remnant (SNR) G344.7--0.1, along with extended TeV source, HESS J1702-420. The observations and the spectral energy distribution modeling support a scenario where this gamma-ray source is the byproduct of the interaction between the SNR shock and the dense surrounding medium, with escaping cosmic rays (CRs) diffusing into the dense environment and interacting with a large local cloud, generating the observed TeV emission. If confirmed, an interaction between the SNR CRs and a nearby cloud would make 2FHL J1703.4-4145 another promising candidate for efficient particle acceleration of the 2FHL Galactic sample, following the first candidate from our previous investigation of a likely shock-cloud interaction occurring on the West edge of the Vela SNR., Comment: 9 pages, 7 figures, Submitted to ApJ June 15, 2020. Accepted for publication Oct 2, 2020

The uncharted territory of the high-redshift ($z\gtrsim3$) Universe holds the key to understand the evolution of quasars. In an attempt to identify the most extreme members of the quasar population, i.e., blazars, we have carried out a multi-wavelength study of a large sample of radio-loud quasars beyond $z=3$. Our sample consists of 9 $\gamma$-ray detected blazars and 133 candidate blazars selected based on the flatness of their soft X-ray spectra (0.3$-$10 keV photon index $\leq1.75$), including 15 with NuSTAR observations. The application of the likelihood profile stacking technique reveals that the high-redshift blazars are faint $\gamma$-ray emitters with steep spectra. The high-redshift blazars host massive black holes ($\langle \log~M_{\rm BH,~M{\odot}} \rangle>9$) and luminous accretion disks ($\langle L_{\rm disk} \rangle>10^{46}$ erg s$^{-1}$). Their broadband spectral energy distributions are found to be dominated by high-energy radiation indicating their jets to be among the most luminous ones. Focusing on the sources exhibiting resolved X-ray jets (as observed with the Chandra satellite), we find the bulk Lorentz factor to be larger with respect to other $z>3$ blazars, indicating faster moving jets. We conclude that the presented list of the high-redshift blazars may act as a reservoir for follow-up observations, e.g., with NuSTAR, to understand the evolution of relativistic jets at the dawn of the Universe., Comment: To appear in the Astrophysical journal. Full tables and plots are provided

We report a detailed study of an unidentified gamma-ray source located in the region of the compact stellar cluster NGC 3603. This is a star-forming region (SFR) powered by a massive cluster of OB stars. A dedicated analysis of about 10 years of data from 10 GeV - 1 TeV, provided by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope, yields the detection of a pointlike source at a significance of 9$\sigma$. The source photon spectrum can be described by a power-law model with best fit spectral index of $2.35 \pm 0.03$. In addition, the analysis of a deep Chandra image in the 0.5 - 7 keV band reliably rules out an extragalactic origin for the gamma rays. We also conclude that the broadband spectral energy distribution of the point source can be explained well with both leptonic and hadronic models. No firm evidence of association with any other classes of known gamma-ray emitters is found, therefore we speculate that 4FGL J1115.1-6118 is a gamma-ray emitting SFR., Comment: 16 pages, 9 figures, 7 tables. Submitted to ApJ

In modeling the X-ray spectra of active galactic nuclei (AGNs), the inclination angle is a parameter that can play an important role in analyzing the X-ray spectra of AGN, but it has never been studied in detail. We present a broadband X-ray spectral analysis of the joint NuSTAR-XMM-Newton observations of 13 sources with [OIII] measured inclinations determined by Fischer et al. 2013. By freezing the inclination angles at the [OIII] measured values when modeling the observations, the spectra are well fitted and the geometrical properties of the obscuring structure of the AGNs are slightly better constrained than those fitted when the inclination angles are left free to vary. We also test if one could freeze the inclinations at other specific angles in fitting the AGN X-ray spectra as commonly did in the literatures. We find that one should always let the inclination angle free to vary in modeling the X-ray spectra of AGNs, while fixing the inclination angle at [OIII] measured values and fixing the inclination angle at 60$^\circ$ also present correct fits of the sources in our sample.Correlations between the covering factor and the average column density of the obscuring torus with respect to the Eddington ratio are also measured, suggesting that the distribution of the material in the obscuring torus is regulated by the Eddington ratio, which is in agreement with previous studies. In addition, no geometrical correlation is found between the narrow line region of the AGN and the obscuring torus, suggesting that the geometry might be more complex than what is assumed in the simplistic unified model., Comment: 22 pages, 10 figures, accepted to ApJ

A majority of the $\gamma$-ray emission from star-forming galaxies is generated by the interaction of high-energy cosmic rays with the interstellar gas and radiation fields. Star-forming galaxies are expected to contribute to both the extragalactic $\gamma$-ray background and the IceCube astrophysical neutrino flux. Using roughly 10\,years of $\gamma$-ray data taken by the {\it Fermi} Large Area Telescope, in this study we constrain the $\gamma$-ray properties of star-forming galaxies. We report the detection of 11 bona-fide $\gamma$-ray emitting galaxies and 2 candidates. Moreover, we show that the cumulative $\gamma$-ray emission of below-threshold galaxies is also significantly detected at $\sim$5\,$\sigma$ confidence. The $\gamma$-ray luminosity of resolved and unresolved galaxies is found to correlate with the total (8-1000\,$\mu$m) infrared luminosity as previously determined. Above 1\,GeV, the spectral energy distribution of resolved and unresolved galaxies is found to be compatible with a power law with a photon index of $\approx2.2-2.3$. Finally, we find that star-forming galaxies account for roughly 5\,\% and 3\,\% of the extragalactic $\gamma$-ray background and the IceCube neutrino flux, respectively., Comment: Accepted for publication in The Astrophysical Journal

We use nine years of gamma-ray data provided by the Fermi Large Area Telescope (LAT) to systematically study the light curves of more than two thousand active galactic nuclei (AGN) included in recent Fermi-LAT catalogs. Ten different techniques are used, which are organized in an automatic periodicity-search pipeline, in order to search for evidence of periodic emission in gamma rays. Understanding the processes behind this puzzling phenomenon will provide a better view about the astrophysical nature of these extragalactic sources. However, the observation of temporal patterns in gamma-ray light curves of AGN is still challenging. Despite the fact that there have been efforts on characterizing the temporal emission of some individual sources, a systematic search for periodicities by means of a full likelihood analysis applied to large samples of sources was missing. Our analysis finds 11 AGN, of which 9 are identified for the first time, showing periodicity at more than 4sigma in at least four algorithms. These findings will help in solving questions related to the astrophysical origin of this periodic behavior., Comment: 16 pages, 5 figures, 4 tables. Accepted by ApJ

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]$., Comment: 17 pages, 7 figures, 8 tables, 1 appendix, accepted for publication in ApJ

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process., Comment: 24 pages, 9 figures, Accepted to ApJ

LOX is a lunar-orbiting astrophysics mission that will probe the cosmos at MeV energies. It is guided by open questions regarding thermonuclear, or Type-Ia, supernovae (SNeIa) and will characterize these inherently radioactive objects by enabling a systematic survey of SNeIa at gamma-ray energies for the first time. Astronomical investigations from lunar orbit afford new opportunities to advance our understanding of the cosmos. The foundation of LOX is an observational approach well suited to the all-sky monitoring demands of supernova investigations and time-domain astronomy. Its inherently wide field-of-view and continuous all-sky monitoring provides an innovative way of addressing decadal survey questions at MeV energies (0.1-10 MeV). The LOX approach achieves high sensitivity with a simple, high-heritage instrument design that eliminates the need for complex, position-sensitive detectors, kinematic event reconstruction, masks, or other insensitive detector mass, while also mitigating technology development, implementation complexity, and their associated costs. LOX can be realized within existing programs, like Explorer., Comment: Astro2020 APC White Paper, Submitted to the National Academy of Sciences Decadal Survey on Astronomy and Astrophysics (Astro2020); 15 pages

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., Comment: 16 pages, 7 figures, 6 tables, to appear in the Astrophysical Journal. Full data tables and SED plots are included in the arXiv uploaded files

We present the results of a detailed multi-wavelength study of a sample of sixteen narrow-line Seyfert 1 (NLSy1) galaxies known in $\gamma$-rays so far. This includes a focused X-ray spectral analysis, to compare them with a more general NLSy1 population, and a broadband spectral energy distribution (SED) modeling considering the leptonic radiative processes to understand their similarity/dissimilarity with blazars. The X-ray spectra of $\gamma$-NLSy1 galaxies exhibit similarities both with typical NLSy1 sources below $\lesssim$2 keV and blazars at higher energies. We also find weak evidences of excess absorption in the soft X-ray spectra of 3 sources and a tentative detection of Fe line in 1H~0323+342. A comparison of the broadband parameters derived from the SED modeling reveals the similarity of $\gamma$-NLSy1 objects with blazars, in particular more with flat spectrum radio quasars. The $\gamma$-NLS1 galaxies host relatively low power jets with small bulk Lorentz factors with respect to blazars and this explains the paucity of such sources in the $\gamma$-ray catalogs and identification of an even fewer parent population members. Many of the observed features of these intriguing objects suggest they host low-mass black holes. The accretion rate in Eddington units can be considered as the driving factor to explain the physical properties of $\gamma$-NLSy1 sources, similar to that which is known for general blazar population., Comment: 34 pages, 11 figures, 8 tables, to appear in the Astrophysical journal

The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail., Comment: 12 pages, 3 figures. Accepted for publication in MNRAS

We report the first results obtained from our campaign to characterize the intranight-optical variability (INOV) properties of {\it Fermi}~detected blazars, using the observations from the recently commissioned 1.3 m J C Bhattacharya telescope (JCBT). During the first run, we were able to observe 17 blazars in the Bessel $R$ filter for $\sim$137 hrs. Using $C$ and scaled $F$-statistics, we quantify the extent of INOV and derive the duty cycle (DC) which is the fraction of time during which a source exhibits a substantial flux variability. We find a high DC of 40\% for BL Lac objects and the flat spectrum radio quasars are relatively less variable (DC $\sim15$\%). However, when estimated for blazars sub-classes, a high DC of $\sim$59\% is found in low synchrotron peaked (LSP) blazars, whereas, intermediate and high synchrotron peaked objects have a low DC of $\sim$11\% and 13\%, respectively. We find evidences about the association of the high amplitude INOV with the $\gamma$-ray flaring state. We also notice a high polarization during the elevated INOV states (for the sources that have polarimetric data available), thus supporting the jet based origin of the observed variability. We plan to enlarge the sample and utilize the time availability from the small telescopes, such as 1.3 m JCBT, to strengthen/verify the results obtained in this work and those existed in the literature., Comment: 9 figures, 5 tables, published in the Astrophysical Journal

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., Comment: 14 pages, 9 figures, 2 tables, to appear in the Astrophysical Journal

We present the 2--78 keV spectral analysis of the deep NuSTAR and XMM-Newton observation of a nearby Seyfert 2 galaxy, ESO 116-G018, which is selected as a candidate Compton-thick (CT-) active galactic nucleus (AGN) based on a previous Chandra-Swift/BAT study. Through our analysis, the source is for the first time confirmed to be a CT-AGN at a $>$3 $\sigma$ confidence level, with the 'line-of-sight' column density N$_{\rm H,Z}$ = [2.46-2.76] $\times$ $10^{24}$ cm$^{-2}$. The 'global average' column density of the obscuring torus is N$_{\rm H,S}$ = [0.46-0.62] $\times$ $10^{24}$ cm$^{-2}$, which suggests a clumpy, rather than uniform, distribution of the obscuring material surrounding the accreting supermassive black hole. The excellent-quality data given by the combined NuSTAR and XMM-Newton observations enable us to produce a strong constraint on the covering factor of the torus of ESO 116-G018, which is found to be $f_c$ = [0.13-0.15]. We also estimate the bolometric luminosity from the broad-band X-ray spectrum, being L$\rm _{bol}$ = [2.57-3.41] $\times$ 10$^{44}$ erg s$^{-1}$., Comment: 11 pages, 5 figures

The gamma-ray sky has been observed with unprecedented accuracy in the last decade by the Fermi large area telescope (LAT), allowing us to resolve and understand the high-energy Universe. The nature of the remaining unresolved emission (unresolved gamma-ray background, UGRB) below the LAT source detection threshold can be uncovered by characterizing the amplitude and angular scale of the UGRB fluctuation field. This work presents a measurement of the UGRB autocorrelation angular power spectrum based on eight years of Fermi LAT Pass 8 data products. The analysis is designed to be robust against contamination from resolved sources and noise systematics. The sensitivity to subthreshold sources is greatly enhanced with respect to previous measurements. We find evidence (with $\sim$3.7$\sigma$ significance) that the scenario in which two classes of sources contribute to the UGRB signal is favored over a single class. A double power law with exponential cutoff can explain the anisotropy energy spectrum well, with photon indices of the two populations being 2.55 $\pm$ 0.23 and 1.86 $\pm$ 0.15., Comment: 8 + 12 (SOM) pages, Contact author: Michela Negro (michela.negro@to.infn.it). (Typos in Tab. VII in previous versions)

The light emitted by all galaxies over the history of the Universe produces the extragalactic background light (EBL) at ultraviolet, optical, and infrared wavelengths. The EBL is a source of opacity for $\gamma$ rays via photon-photon interactions, leaving an imprint in the spectra of distant $\gamma$-ray sources. We measure this attenuation using {739} active galaxies and one gamma-ray burst detected by the {\it Fermi} Large Area Telescope. This allows us to reconstruct the evolution of the EBL and determine the star-formation history of the Universe over 90\% of cosmic time. Our star-formation history is consistent with independent measurements from galaxy surveys, peaking at redshift $z\sim2$. Upper limits of the EBL at the epoch of re-ionization suggest a turnover in the abundance of faint galaxies at $z\sim 6$., Comment: Published on Science. This is the authors' version of the manuscript

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., Comment: 12 pages, 6 figures

The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region., Comment: Accepted for publication in the ApJ, Corresponding author: Nahee Park (VERITAS Collaboration), John W. Hewitt (Fermi-LAT Collaboration), Ignacio Taboada (HAWC Collaboration), 30 pages, 12 figures