The gamma -ray sky >100 MeV is dominated by the diffuse emissions from interactions of cosmic rays with the interstellar gas and radiation fields of the Milky Way. Observations of these diffuse ...emissions provide a tool to study cosmic-ray origin and propagation, and the interstellar medium. We present measurements from the first 21 months of the Fermi Large Area Telescope (Fermi-LAT) mission and compare with models of the diffuse gamma -ray emission generated using the GALPROP code. The models are fitted to cosmic-ray data and incorporate astrophysical input for the distribution of cosmic-ray sources, interstellar gas, and radiation fields. To assess uncertainties associated with the astrophysical input, a grid of models is created by varying within observational limits the distribution of cosmic-ray sources, the size of the cosmic-ray confinement volume (halo), and the distribution of interstellar gas. An all-sky maximum-likelihood fit is used to determine the X sub(CO) factor, the ratio between integrated CO-line intensity and H sub(2) column density, the fluxes and spectra of the gamma -ray point sources from the first Fermi-LAT catalog, and the intensity and spectrum of the isotropic background including residual cosmic rays that were misclassified as gamma -rays, all of which have some dependency on the assumed diffuse emission model. The models are compared on the basis of their maximum-likelihood ratios as well as spectra, longitude, and latitude profiles. We also provide residual maps for the data following subtraction of the diffuse emission models. The models are consistent with the data at high and intermediate latitudes but underpredict the data in the inner Galaxy for energies above a few GeV. Possible explanations for this discrepancy are discussed, including the contribution by undetected point-source populations and spectral variations of cosmic rays throughout the Galaxy. In the outer Galaxy, we find that the data prefer models with a flatter distribution of cosmic-ray sources, a larger cosmic-ray halo, or greater gas density than is usually assumed. Our results in the outer Galaxy are consistent with other Fermi-LAT studies of this region that used different analysis methods than employed in this paper.
The second catalog of active galactic nuclei (AGNs) detected by the Fermi Large Area Telescope (LAT) in two years of scientific operation is presented. The second LAT AGN catalog (2LAC) includes 1017 ...Delta *g-ray sources located at high Galactic latitudes (|b| > 10?) that are detected with a test statistic (TS) greater than 25 and associated statistically with AGNs. However, some of these are affected by analysis issues and some are associated with multiple AGNs. Consequently, we define a Clean Sample which includes 886 AGNs, comprising 395 BL Lacertae objects (BL Lac objects), 310 flat-spectrum radio quasars (FSRQs), 157 candidate blazars of unknown type (i.e., with broadband blazar characteristics but with no optical spectral measurement yet), 8 misaligned AGNs, 4 narrow-line Seyfert 1 (NLS1s), 10 AGNs of other types, and 2 starburst galaxies. Where possible, the blazars have been further classified based on their spectral energy distributions (SEDs) as archival radio, optical, and X-ray data permit. While almost all FSRQs have a synchrotron-peak frequency <1014 Hz, about half of the BL Lac objects have a synchrotron-peak frequency >1015 Hz. The 2LAC represents a significant improvement relative to the first LAT AGN catalog (1LAC), with 52% more associated sources. The full characterization of the newly detected sources will require more broadband data. Various properties, such as Delta *g-ray fluxes and photon power-law spectral indices, redshifts, Delta *g-ray luminosities, variability, and archival radio luminosities and their correlations are presented and discussed for the different blazar classes. The general trends observed in 1LAC are confirmed.
We have performed an analysis of the diffuse gamma-ray emission with the Fermi Large Area Telescope (LAT) in the Milky Way halo region, searching for a signal from dark matter annihilation or decay. ...In the absence of a robust dark matter signal, constraints are presented. We consider both gamma rays produced directly in the dark matter annihilation/decay and produced by inverse Compton scattering of the e super(+)/e super(-) produced in the annihilation/decay. Conservative limits are derived requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission using the GALPROP code. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the index of the injection cosmic-ray electron spectrum, and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as the annihilation of dark matter.
Context.
Many questions concerning the nature of astrophysical counterparts of high-energy neutrinos remain unanswered. There is increasing evidence of a connection between blazar jets and neutrino ...events, with the flare of the
γ
-ray blazar TXS 0506+056 in spatial and temporal proximity of IC 170922A representing one of the most outstanding associations of high-energy neutrinos with astrophysical sources reported so far.
Aims.
With the purpose of characterising potential blazar counterparts to high-energy neutrinos, we analysed the parsec-scale regions of
γ
-ray blazars in spatial coincidence with high-energy neutrinos, detected by the IceCube Observatory. Specifically, we intended to investigate peculiar radio properties of the candidate counterparts related to the neutrino production, such as radio flares coincident with the neutrino detection or features in jet morphology (limb brightening, transverse structures).
Methods.
We collected multi-frequency, very-long-baseline interferometry (VLBI) follow-up observations of candidate counterparts of four high-energy neutrino events detected by IceCube between January 2019 and November 2020, with a focus on
γ
-ray-associated objects. We analysed their radio characteristics soon after the neutrino arrival in comparison with archival VLBI observations and low-frequency radio observations. We discussed our results with respect to previous statistical works and studies on the case of TXS 0506+056.
Results.
We identified and analysed five potential neutrino-emitting blazars in detail. Our results suggest an enhanced state of activity for one source, PKS 1725+123. However, the lack of adequate monitoring prior to the neutrino events was a limitation in tracing radio activity and morphological changes in all the sources.
Conclusions.
We suggest that PKS 1725+123 is a promising neutrino source candidate. For the other sources, our results alone do not reveal a strong connection between the radio activity state at the neutrino arrival. A larger number of VLBI and multi-wavelength follow-up observations of neutrino events are now essential to our understanding of the neutrino production mechanisms in astrophysical sources.
Abstract
Flat-spectrum radio quasars (FSRQs) are the most luminous blazars at GeV energies but only rarely emit detectable fluxes of TeV gamma rays, typically during bright GeV flares. We explore the ...gamma-ray variability and spectral characteristics of three FSRQs that have been observed at GeV and TeV energies by Fermi-LAT and VERITAS, making use of almost 100 hr of VERITAS observations spread over 10 yr: 3C 279, PKS 1222+216, and Ton 599. We explain the GeV flux distributions of the sources in terms of a model derived from a stochastic differential equation describing fluctuations in the magnetic field in the accretion disk and estimate the timescales of magnetic flux accumulation and stochastic instabilities in their accretion disks. We identify distinct flares using a procedure based on Bayesian blocks and analyze their daily and subdaily variability and gamma-ray energy spectra. Using observations from VERITAS, as well as Fermi, Swift, and the Steward Observatory, we model the broadband spectral energy distributions of PKS 1222+216 and Ton 599 during very high energy (VHE)–detected flares in 2014 and 2017, respectively, strongly constraining the jet Doppler factors and gamma-ray emission region locations during these events. Finally, we place theoretical constraints on the potential production of PeV-scale neutrinos during these VHE flares.
Classical novae are shock-powered multi-wavelength transients triggered by a thermonuclear runaway on an accreting white dwarf. V1674 Her is the fastest nova ever recorded (time to declined by two ...magnitudes is t_2=1.1 d) that challenges our understanding of shock formation in novae. We investigate the physical mechanisms behind nova emission from GeV gamma-rays to cm-band radio using coordinated Fermi-LAT, NuSTAR, Swift and VLA observations supported by optical photometry. Fermi-LAT detected short-lived (18 h) 0.1-100 GeV emission from V1674 Her that appeared 6 h after the eruption began; this was at a level of (1.6 +/- 0.4)x10^-6 photons cm^-2 s^-1. Eleven days later, simultaneous NuSTAR and Swift X-ray observations revealed optically thin thermal plasma shock-heated to kT_shock = 4 keV. The lack of a detectable 6.7 keV Fe K_alpha emission suggests super-solar CNO abundances. The radio emission from V1674 Her was consistent with thermal emission at early times and synchrotron at late times. The radio spectrum steeply rising with frequency may be a result of either free-free absorption of synchrotron and thermal emission by unshocked outer regions of the nova shell or the Razin-Tsytovich effect attenuating synchrotron emission in dense plasma. The development of the shock inside the ejecta is unaffected by the extraordinarily rapid evolution and the intermediate polar host of this nova.
The Fermi Large Area Telescope (LAT) First Source Catalog (1FGL) provided spatial, spectral, and temporal properties for a large number of gamma -ray sources using a uniform analysis method. After ...correlating with the most-complete catalogs of source types known to emit gamma rays, 630 of these sources are "unassociated" (i.e., have no obvious counterparts at other wavelengths). Here, we employ two statistical analyses of the primary gamma -ray characteristics for these unassociated sources in an effort to correlate their y-ray properties with the active galactic nucleus (AGN) and pulsar populations in 1FGL. Based on the correlation results, we classify 221 AGN-like and 134 pulsar-like sources in the 1FGL unassociated sources. The results of these source "classifications" appear to match the expected source distributions, especially at high Galactic latitudes. While useful for planning future multiwavelength follow-up observations, these analyses use limited inputs, and their predictions should not be considered equivalent to "probable source classes" for these sources. We discuss multiwavelength results and catalog cross-correlations to date, and provide new source associations for 229 Fermi-LAT sources that had no association listed in the 1FGL catalog. By validating the source classifications against these new associations, we find that the new association matches the predicted source class in ~80% of the sources.
Context. Globular clusters with their large populations of millisecond pulsars (MSPs) are believed to be potential emitters of high-energy gamma-ray emission. The observation of this emission ...provides a powerful tool to assess the millisecond pulsar population of a cluster, is essential for understanding the importance of binary systems for the evolution of globular clusters, and provides complementary insights into magnetospheric emission processes. Aims. Our goal is to constrain the millisecond pulsar populations in globular clusters from analysis of gamma-ray observations. Methods. We use 546 days of continuous sky-survey observations obtained with the Large Area Telescope aboard the Fermi Gamma-ray Space Telescope to study the gamma-ray emission towards 13 globular clusters. Results. Steady point-like high-energy gamma-ray emission has been significantly detected towards 8 globular clusters. Five of them (47 Tucanae, Omega Cen, NGC 6388, Terzan 5, and M 28) show hard spectral power indices (0.7 < Γ < 1.4) and clear evidence for an exponential cut-off in the range 1.0 − 2.6 GeV, which is the characteristic signature of magnetospheric emission from MSPs. Three of them (M 62, NGC 6440 and NGC 6652) also show hard spectral indices (1.0 < Γ < 1.7), however the presence of an exponential cut-off can not be unambiguously established. Three of them (Omega Cen, NGC 6388, NGC 6652) have no known radio or X-ray MSPs yet still exhibit MSP spectral properties. From the observed gamma-ray luminosities, we estimate the total number of MSPs that is expected to be present in these globular clusters. We show that our estimates of the MSP population correlate with the stellar encounter rate and we estimate 2600 − 4700 MSPs in Galactic globular clusters, commensurate with previous estimates. Conclusions. The observation of high-energy gamma-ray emission from globular clusters thus provides a reliable independent method to assess their millisecond pulsar populations.
Context. Cosmic rays (CRs) can be studied through the galaxy-wide gamma-ray emission that they generate when propagating in the interstellar medium. The comparison of the diffuse signals from ...different systems may inform us about the key parameters in CR acceleration and transport. Aims. We aim to determine and compare the properties of the cosmic-ray-induced gamma-ray emission of several Local Group galaxies. Methods. We use 2 years of nearly continuous sky-survey observations obtained with the Large Area Telescope aboard the Fermi Gamma-ray Space Telescope to search for gamma-ray emission from M 31 and M 33. We compare the results with those for the Large Magellanic Cloud, the Small Magellanic Cloud, the Milky Way, and the starburst galaxies M 82 and NGC 253. Results. We detect a gamma-ray signal at 5σ significance in the energy range 200 MeV–20 GeV that is consistent with originating from M 31. The integral photon flux above 100 MeV amounts to (9.1 ± 1.9stat ± 1.0sys) × 10-9 ph cm-2 s-1. We find no evidence for emission from M 33 and derive an upper limit on the photon flux >100 MeV of 5.1 × 10-9 ph cm-2 s-1 (2σ). Comparing these results to the properties of other Local Group galaxies, we find indications of a correlation between star formation rate and gamma-ray luminosity that also holds for the starburst galaxies. Conclusions. The gamma-ray luminosity of M 31 is about half that of the Milky Way, which implies that the ratio between the average CR densities in M 31 and the Milky Way amounts to ξ = 0.35 ± 0.25. The observed correlation between gamma-ray luminosity and star formation rate suggests that the flux of M 33 is not far below the current upper limit from the LAT observations.