The origin of high-energy galactic cosmic rays is yet to be understood, but some galactic cosmic-ray accelerators can accelerate cosmic rays up to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding material or radiation, resulting in the production of gamma-rays and neutrinos. To optimize for the detection of such associated production of gamma-rays and neutrinos for a given source morphology and spectrum, a multimessenger analysis that combines gamma-rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework with IceCube Maximum Likelihood Analysis software and HAWC Accelerated Likelihood to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 yr of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that, for five of the sources, the gamma-ray emission observed by HAWC cannot be produced purely from hadronic interactions. We report the limit for the fraction of gamma-rays produced by hadronic interactions for these five sources.
A Compound Poisson Generator Approach to Point-source Inference in Astrophysics
Abstract The identification and description of point sources is one of the oldest problems in astronomy, yet even today the correct statistical treatment for point sources remains one of the field’s hardest problems. For dim or crowded sources, likelihood-based inference methods are required to estimate the uncertainty on the characteristics of the source population. In this work, a new parametric likelihood is constructed for this problem using compound Poisson generator (CPG) functionals that incorporate instrumental effects from first principles. We demonstrate that the CPG approach exhibits a number of advantages over non-Poissonian template fitting (NPTF)—an existing method—in a series of test scenarios in the context of X-ray astronomy. These demonstrations show that the effect of the point-spread function, effective area, and choice of point-source spatial distribution cannot, generally, be factorized as they are in NPTF, while the new CPG construction is validated in these scenarios. Separately, an examination of the diffuse-flux emission limit is used to show that most simple choices of priors on the standard parameterization of the population model can result in unexpected biases: when a model comprising both a point-source population and diffuse component is applied to this limit, nearly all observed flux will be assigned to either the population or to the diffuse component. A new parameterization is presented for these priors that properly estimates the uncertainties in this limit. In this choice of priors, CPG correctly identifies that the fraction of flux assigned to the population model cannot be constrained by the data.
more »
« less
- Award ID(s):
- 2019786
- PAR ID:
- 10343316
- Date Published:
- Journal Name:
- The Astrophysical Journal Supplement Series
- Volume:
- 260
- Issue:
- 2
- ISSN:
- 0067-0049
- Page Range / eLocation ID:
- 29
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
Abstract γ -ray observations of the Cygnus Cocoon, an extended source surrounding the Cygnus X star-forming region, suggest the presence of a cosmic-ray accelerator reaching energies up to a few PeV. The very-high-energy (VHE; 0.1–100 TeV) γ -ray emission may be explained by the interaction of cosmic-ray hadrons with matter inside the Cocoon, but an origin of inverse Compton radiation by relativistic electrons cannot be ruled out. Inverse Compton γ -rays at VHE are accompanied by synchrotron radiation peaked in X-rays. Hence, X-ray observations may probe the electron population and magnetic field of the source. We observed 11 fields in or near the Cygnus Cocoon with the Neil Gehrels Swift Observatory’s X-Ray Telescope (Swift-XRT) totaling 110 ks. We fit the fields to a Galactic and extragalactic background model and performed a log-likelihood ratio test for an additional diffuse component. We found no significant additional emission and established upper limits in each field. By assuming that the X-ray intensity traces the TeV intensity and follows a dN / dE ∝ E − 2.5 spectrum, we obtained a 90% upper limit of F X < 8.7 × 10 −11 erg cm −2 s −1 or <5.2 × 10 −11 erg cm −2 s −1 on the X-ray flux of the entire Cygnus Cocoon between 2 and 10 keV depending on the choice of hydrogen column density model for the absorption. The obtained upper limits suggest that no more than one-quarter of the γ -ray flux at 1 TeV is produced by inverse Compton scattering, when assuming an equipartition magnetic field of ∼20 μ G.more » « less
-
Abstract The majority of astrophysical neutrinos have undetermined origins. The IceCube Neutrino Observatory has observed astrophysical neutrinos but has not yet identified their sources. Blazars are promising source candidates, but previous searches for neutrino emission from populations of blazars detected in ≳GeV gamma rays have not observed any significant neutrino excess. Recent findings in multimessenger astronomy indicate that high-energy photons, coproduced with high-energy neutrinos, are likely to be absorbed and reemitted at lower energies. Thus, lower-energy photons may be better indicators of TeV–PeV neutrino production. This paper presents the first time-integrated stacking search for astrophysical neutrino emission from MeV-detected blazars in the first Fermi Large Area Telescope low energy (1FLE) catalog using ten years of IceCube muon–neutrino data. The results of this analysis are found to be consistent with a background-only hypothesis. Assuming an E −2 neutrino spectrum and proportionality between the blazars MeV gamma-ray fluxes and TeV–PeV neutrino flux, the upper limit on the 1FLE blazar energy-scaled neutrino flux is determined to be 1.64 × 10 −12 TeV cm −2 s −1 at 90% confidence level. This upper limit is approximately 1% of IceCube’s diffuse muon–neutrino flux measurement.more » « less
-
null (Ed.)The Bayesian formulation of inverse problems is attractive for three primary reasons: it provides a clear modelling framework; it allows for principled learning of hyperparameters; and it can provide uncertainty quantification. The posterior distribution may in principle be sampled by means of MCMC or SMC methods, but for many problems it is computationally infeasible to do so. In this situation maximum a posteriori (MAP) estimators are often sought. Whilst these are relatively cheap to compute, and have an attractive variational formulation, a key drawback is their lack of invariance under change of parameterization; it is important to study MAP estimators, however, because they provide a link with classical optimization approaches to inverse problems and the Bayesian link may be used to improve upon classical optimization approaches. The lack of invariance of MAP estimators under change of parameterization is a particularly significant issue when hierarchical priors are employed to learn hyperparameters. In this paper we study the effect of the choice of parameterization on MAP estimators when a conditionally Gaussian hierarchical prior distribution is employed. Specifically we consider the centred parameterization, the natural parameterization in which the unknown state is solved for directly, and the noncentred parameterization, which works with a whitened Gaussian as the unknown state variable, and arises naturally when considering dimension-robust MCMC algorithms; MAP estimation is well-defined in the nonparametric setting only for the noncentred parameterization. However, we show that MAP estimates based on the noncentred parameterization are not consistent as estimators of hyperparameters; conversely, we show that limits of finite-dimensional centred MAP estimators are consistent as the dimension tends to infinity. We also consider empirical Bayesian hyperparameter estimation, show consistency of these estimates, and demonstrate that they are more robust with respect to noise than centred MAP estimates. An underpinning concept throughout is that hyperparameters may only be recovered up to measure equivalence, a well-known phenomenon in the context of the Ornstein–Uhlenbeck process. The applicability of the results is demonstrated concretely with the study of hierarchical Whittle–Matérn and ARD priors.more » « less
-
more » « less
ABSTRACT We present the largest low frequency (120 MHz) arcminute resolution image of the radio synchrotron background (RSB) to date, and its corresponding angular power spectrum of anisotropies (APS) with angular scales ranging from 3° to 0.3 arcmin. We show that the RSB around the north celestial pole has a significant excess anisotropy power at all scales over a model of unclustered point sources based on source counts of known source classes. This anisotropy excess, which does not seem attributable to the diffuse Galactic emission, could be linked to the surface brightness excess of the RSB. To better understand the information contained within the measured APS, we model the RSB varying the brightness distribution, size, and angular clustering of potential sources. We show that the observed APS could be produced by a population of faint clustered point sources only if the clustering is extreme and the size of the Gaussian clusters is ≲1 arcmin. We also show that the observed APS could be produced by a population of faint diffuse sources with sizes ≲1 arcmin, and this is supported by features present in our image. Both of these cases would also cause an associated surface brightness excess. These classes of sources are in a parameter space not well probed by even the deepest radio surveys to date.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3847/1538-4365/ac5cb7
