ABSTRACT We present the third Fermi Large Area Telescope (LAT) source catalog (3FGL) of sources in the 100 MeV-300 GeV range. Based on the first 4 yr of science data from the Fermi Gamma-ray Space ...Telescope mission, it is the deepest yet in this energy range. Relative to the Second Fermi LAT catalog, the 3FGL catalog incorporates twice as much data, as well as a number of analysis improvements, including improved calibrations at the event reconstruction level, an updated model for Galactic diffuse γ-ray emission, a refined procedure for source detection, and improved methods for associating LAT sources with potential counterparts at other wavelengths. The 3FGL catalog includes 3033 sources above significance, with source location regions, spectral properties, and monthly light curves for each. Of these, 78 are flagged as potentially being due to imperfections in the model for Galactic diffuse emission. Twenty-five sources are modeled explicitly as spatially extended, and overall 238 sources are considered as identified based on angular extent or correlated variability (periodic or otherwise) observed at other wavelengths. For 1010 sources we have not found plausible counterparts at other wavelengths. More than 1100 of the identified or associated sources are active galaxies of the blazar class; several other classes of non-blazar active galaxies are also represented in the 3FGL. Pulsars represent the largest Galactic source class. From source counts of Galactic sources we estimate that the contribution of unresolved sources to the Galactic diffuse emission is ∼3% at 1 GeV.
We analyze the prompt emission of GRB 100724B and GRB 160509A, two of the brightest gamma-ray bursts (GRBs) observed by Fermi at MeV energies but surprisingly faint at 100 MeV energies. Time-resolved ...spectroscopy reveals a sharp high-energy cutoff at energies Ec ∼ 20-60 MeV for GRB 100724B and Ec ∼ 80-150 MeV for GRB 160509A. We first characterize phenomenologically the cutoff and its time evolution. We then fit the data to two models where the high-energy cutoff arises from intrinsic opacity to pair production within the source (τγγ): (i) a Band spectrum with τγγ from the internal-shocks-motivated model of Granot et al. (2008) and (ii) the photospheric model of Gill & Thompson (2014). Alternative explanations for the cutoff, such as an intrinsic cutoff in the emitting electron energy distribution, appear to be less natural. Both models provide a good fit to the data with very reasonable physical parameters, providing an estimate of bulk Lorentz factors in the range Γ ∼ 100-400, on the lower end of what is generally observed in Fermi GRBs. Surprisingly, their lower cutoff energies Ec compared to other Fermi/LAT GRBs arise not only predominantly from the lower Lorentz factors, but also at a comparable level from differences in variability time, luminosity, and high-energy photon index. Finally, particularly low Ec values may prevent detection by Fermi/LAT, thus introducing a bias in the Fermi/LAT GRB sample against GRBs with low Lorentz factors or variability times.
ABSTRACT
Indirect detection strategies of particle dark matter (DM) in Dwarf spheroidal satellite galaxies (dSphs) typically entail searching for annihilation signals above the astrophysical ...background. To robustly compare model predictions with the observed fluxes of product particles, most analyses of astrophysical data – which are generally frequentist – rely on estimating the abundance of DM by calculating the so-called J factor. This quantity is usually inferred from the kinematic properties of the stellar population of a dSph using the Jeans equation, commonly by means of Bayesian techniques that entail the presence (and additional systematic uncertainty) of prior choice. Here, extending earlier work, we develop a scheme to derive the profile likelihood for J factors of dwarf spheroidals for models with five or more free parameters. We validate our method on a publicly available simulation suite, released by the Gaia Challenge, finding satisfactory statistical properties for bias and probability coverage. We present the profile likelihood function and maximum likelihood estimates for the J-factor of 10 dSphs. As an illustration, we apply these profile likelihoods to recently published analyses of γ-ray data with the Fermi Large Area Telescope to derive new, consistent upper limits on the DM annihilation cross-section. We do this for a subset of systems, generally referred to as classical dwarfs. The implications of these findings for DM searches are discussed, together with future improvements and extensions of this technique.
ABSTRACT
Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior ...probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing 12 photo-z algorithms applied to mock data produced for The Rubin Observatory Legacy Survey of Space and Time Dark Energy Science Collaboration. By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. We report systematic biases and overall over/underbreadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performance metrics.
Abstract
TeV photons from extragalactic sources are absorbed in the intergalactic medium and initiate electromagnetic cascades. These cascades offer a unique tool to probe the properties of the ...universe at cosmological scales. We present a new Monte Carlo code dedicated to the physics of such cascades. This code has been tested against both published results and analytical approximations, and is made publicly available. Using this numerical tool, we investigate the main cascade properties (spectrum, halo extension and time delays), and study in detail their dependence on the physical parameters (extragalactic magnetic field, extragalactic background light, source redshift, source spectrum and beaming emission). The limitations of analytical solutions are emphasized. In particular, analytical approximations account only for the first generation of photons and higher branches of the cascade tree are neglected.
Abstract
Line-of-sight integrals of the squared density, commonly called the J-factor, are essential for inferring dark matter (DM) annihilation signals. The J-factors of DM-dominated dwarf ...spheroidal satellite galaxies (dSphs) have typically been derived using Bayesian techniques, which for small data samples implies that a choice of priors constitutes a non-negligible systematic uncertainty. Here we report the development of a new fully frequentist approach to construct the profile likelihood of the J-factor. Using stellar kinematic data from several classical and ultra-faint dSphs, we derive the maximum likelihood value for the J-factor and its confidence intervals. We validate this method, in particular its bias and coverage, using simulated data from the Gaia Challenge. We find that the method possesses good statistical properties. The J-factors and their uncertainties are generally in good agreement with the Bayesian-derived values, with the largest deviations restricted to the systems with the smallest kinematic data sets. We discuss improvements, extensions, and future applications of this technique.
ABSTRACT The origin of the extragalactic γ-ray background (EGB) has been debated for some time. The EGB comprises the γ-ray emission from resolved and unresolved extragalactic sources, such as ...blazars, star-forming galaxies, and radio galaxies, as well as radiation from truly diffuse processes. This Letter focuses on the blazar source class, the most numerous detected population, and presents an updated luminosity function and spectral energy distribution model consistent with the blazar observations performed by the Fermi-Large Area Telescope (LAT). We show that blazars account for 50 % of the EGB photons (>0.1 GeV), and that Fermi-LAT has already resolved ∼70% of this contribution. Blazars, and in particular hard-spectrum sources such as BL Lacs, are responsible for most of the EGB emission above 100 GeV. We find that the extragalactic background light, which attenuates blazars' high-energy emission, is responsible for the high-energy cutoff observed in the EGB spectrum. Finally, we show that blazars, star-forming galaxies, and radio galaxies can naturally account for the amplitude and spectral shape of the background in the 0.1-820 GeV range, leaving only modest room for other contributions. This allows us to set competitive constraints on the dark matter annihilation cross section.
The launch of the gamma-ray telescope Fermi Large Area Telescope (Fermi-LAT) started a pivotal period in indirect detection of dark matter. By outperforming expectations, for the first time a robust ...and stringent test of the paradigm of weakly interacting massive particles (WIMPs) is within reach. In this paper, we discuss astrophysical targets for WIMP detection and the challenges they present, review the analysis tools which have been employed to tackle these challenges, and summarize the status of constraints on and the claimed detections in the WIMP parameter space. Methods and results will be discussed in comparison to Imaging Air Cherenkov Telescopes. We also provide an outlook on short term and longer term developments.
Spatial extension is an important characteristic for correctly associating gamma-ray-emitting sources with their counterparts at other wavelengths and for obtaining an unbiased model of their ...spectra. We present a new method for quantifying the spatial extension of sources detected by the Large Area Telescope (LAT), the primary science instrument on the Fermi Gamma-ray Space Telescope (Fermi). We perform a series of Monte Carlo simulations to validate this tool and calculate the LAT threshold for detecting the spatial extension of sources. We then test all sources in the second Fermi-LAT catalog for extension. We report the detection of seven new spatially extended sources.
Abstract
The Vera Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to process ∼10
6
transient detections per night. For precision measurements of cosmological parameters and ...rates, it is critical to understand the detection efficiency, magnitude limits, artifact contamination levels, and biases in the selection and photometry. Here we rigorously test the LSST Difference Image Analysis (DIA) pipeline using simulated images from the Rubin Observatory LSST Dark Energy Science Collaboration Data Challenge (DC2) simulation for the Wide-Fast-Deep survey area. DC2 is the first large-scale (300 deg
2
) image simulation of a transient survey that includes realistic cadence, variable observing conditions, and CCD image artifacts. We analyze ∼15 deg
2
of DC2 over a 5 yr time span in which artificial point sources from Type Ia supernova (SNIa) light curves have been overlaid onto the images. The magnitude limits per filter are
u
= 23.66 mag,
g
= 24.69 mag,
r
= 24.06 mag,
i
= 23.45 mag,
z
= 22.54 mag, and
y
= 21.62 mag. The artifact contamination levels are ∼90% of all detections, corresponding to ∼1000 artifacts deg
–2
in
g
band, and falling to 300 deg
–2
in
y
band. The photometry has biases <1% for magnitudes 19.5 <
m
< 23. Our DIA performance on simulated images is similar to that of the Dark Energy Survey difference-imaging pipeline on real images. We also characterize DC2 image properties to produce catalog-level simulations needed for distance bias corrections. We find good agreement between DC2 data and simulations for distributions of signal-to-noise ratio, redshift, and fitted light-curve properties. Applying a realistic SNIa cosmology analysis for redshifts
z
< 1, we recover the input cosmology parameters to within statistical uncertainties.