Point Source Population Research Articles

Are there correlations in the HAWC and IceCube high energy skymaps outside the Galactic plane?

We use publicly available data to perform a search for correlations of high energy neutrino candidate events detected by IceCube and high-energy photons seen by the HAWC Collaboration. Our search is focused on unveiling such correlations outside of the Galactic plane. This search is sensitive to correlations in the neutrino candidate and photon skymaps which would arise from a population of unidentified point sources. We find no evidence for such a correlation, but suggest strategies for improvements with new datasets. Published by the American Physical Society 2024

Gamma-rays and neutrinos from giant molecular cloud populations in the galactic plane

The recent IceCube detection of significant neutrino flux from the inner Galactic plane has provided us valuable insights on the spectrum of cosmic rays in our Galaxy. This flux can be produced either by a population of Galactic point sources or by diffused emission from cosmic ray interactions with the interstellar medium or by a mixture of both. In this work, we compute diffused gamma-ray and neutrino fluxes produced by a population of giant molecular clouds (GMCs) in our Galaxy, assuming different parametrizations of the Galactic diffused cosmic ray distribution. In particular, we take into account two main cases: (I) constant cosmic ray luminosity in our Galaxy, and (II) space-dependent cosmic ray luminosity, based on the supernovae distribution in our Galaxy. For Case-I, we found that the neutrino flux from GMCs is a factor of ∼ 10 below compared to π 0 and KRA γ best-fitted models of IceCube observations at 105 GeV. Instead, for Case-II the model can explain up to ∼ 90 % of the neutrino flux at that energy. Moreover, for this last scenario IceCube detector could be able to detect neutrino events from the Galactic centre regions. We then calculated gamma-ray and neutrino fluxes from individual GMCs and noticed that several current and future Cherenkov telescopes and neutrino observatories have the right sensitivities to study these objects. In particular, very neutrino-bright region such as Aquila Rift is favourable for detection by the IceCube-Gen2 observatory.

The high energy X-ray probe: resolved X-ray populations in extragalactic environments

We construct simulated galaxy data sets based on the High Energy X-ray Probe (HEX-P) mission concept to demonstrate the significant advances in galaxy science that will be yielded by the HEX-P observatory. The combination of high spatial resolution imaging (<20 arcsec FWHM), broad spectral coverage (0.2–80 keV), and sensitivity superior to current facilities (e.g., XMM-Newton and NuSTAR) will enable HEX-P to detect hard (4–25 keV) X-ray emission from resolved point-source populations within ∼800 galaxies and integrated emission from ∼6,000 galaxies out to 100 Mpc. These galaxies cover wide ranges of galaxy types (e.g., normal, starburst, and passive galaxies) and properties (e.g., metallicities and star-formation histories). In such galaxies, HEX-P will: 1) provide unique information about X-ray binary populations, including accretor demographics (black hole and neutron stars), distributions of accretion states and state transition cadences; 2) place order-of-magnitude more stringent constraints on inverse Compton emission associated with particle acceleration in starburst environments; and 3) put into clear context the contributions from X-ray emitting populations to both ionizing the surrounding interstellar medium in low-metallicity galaxies and heating the intergalactic medium in the z > 8 Universe.

Testing the surface brightness fluctuation method on dwarf galaxies in the COSMOS field

ABSTRACT Dwarf galaxies are important tracers of small-scale cosmological structure, yet much of our knowledge about these systems comes from the limited sample of dwarf galaxies within the Local Group. To make a comprehensive inventory of dwarf populations in the local Universe, we require effective methods for deriving distance estimates for large numbers of faint, low surface brightness objects. Here we test the surface brightness fluctuation (SBF) method, traditionally applied to brighter early-type galaxies, on a sample of 20 nearby dwarf galaxies detected in the Cosmological Evolution Survey (COSMOS) field. These objects are partially resolved in the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) images and have confirmed redshift distances in the range 17–130 Mpc. We discuss the many model choices required in applying the SBF method and explore how these affect the final distance estimates. Amongst other variations on the method, when applying the SBF method, we alter the standard equation to include a term accounting for the power spectrum of the background, greatly improving our results. For the most robust modelling choices, we find a roughly Gaussian SBF signal that correlates linearly with distance out to distances of 50–100 Mpc, but with only a fraction of the power expected. At larger distances, there is excess power relative to that predicted, probably from undetected point sources. Overall, obtaining accurate SBF distances to faint, irregular galaxies remains challenging, but may yet prove possible with the inclusion of more information about galaxy properties and point source populations, and the use of more advanced techniques.

The X-ray point source population hosted by globular clusters in the elliptical galaxy NGC 4261

ABSTRACT Utilizing archival Chandra X-ray Observatory data and Hubble Space Telescope globular cluster catalogues, we probe the time-domain properties of the low mass X-ray binary population in the elliptical galaxy NGC 4261. Of the 98 unique X-ray sources identified in this study, 62 sources are within the optical field of view and, of those, 33 per cent are aligned with an optical cluster counterpart. We find twenty X-ray sources coincident with globular clusters; two are previously discovered ultra-luminous X-ray sources (ULXs) and eighteen are low mass X-ray binaries (GCLMXBs) with LX < 1039 erg s−1. ULXs are a heterogeneous class of extremely bright X-ray binaries (LX > 1039 erg s−1) and ULXs located in globular clusters (GCULXs) and may be indicators of black holes. Identifying these unusually X-ray bright sources and measuring their optical properties can provide valuable constraints on the progenitors of gravitational wave sources. We compare observations of these sources to the twenty previously studied GCULXs from five other early-type galaxies, and find that GCULXs in NGC 4261 are of similar colour and luminosity and do not significantly deviate from the rest of the sample in terms of distance from the galaxy centre or X-ray luminosity. Both the GCULX and GCLMXB populations of NGC 4261 show long-term variability; the former may have implications for fast radio bursts originating in globular clusters and the latter will likely introduce additional scatter into the low mass end of GCLMXB X-ray luminosity functions.

Observation of high-energy neutrinos from the Galactic plane.

The origin of high-energy cosmic rays, atomic nuclei that continuously impact Earth's atmosphere, is unknown. Because of deflection by interstellar magnetic fields, cosmic rays produced within the Milky Way arrive at Earth from random directions. However, cosmic rays interact with matter near their sources and during propagation, which produces high-energy neutrinos. We searched for neutrino emission using machine learning techniques applied to 10 years of data from the IceCube Neutrino Observatory. By comparing diffuse emission models to a background-only hypothesis, we identified neutrino emission from the Galactic plane at the 4.5σ level of significance. The signal is consistent with diffuse emission of neutrinos from the Milky Way but could also arise from a population of unresolved point sources.

Diffuse sources, clustering, and the excess anisotropy of the radio synchrotron background

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.

A Compound Poisson Generator Approach to Point-source Inference in Astrophysics

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.

Luminosity functions consistent with a pulsar-dominated Galactic Center excess

A new population of millisecond pulsars is a long-standing proposed explanation for the excess of GeV-scale gamma rays emanating from the region surrounding the center of the Milky Way (the “Galactic Center excess”). We examine several simple parameterizations of possible luminosity functions for this population, as well as several benchmark luminosity functions proposed in the literature, and compare the predicted populations of resolved point sources to the Fermi 4FGL-DR2 point source catalog and a sub-population recently identified using wavelet-based methods. We provide general results that can be used to translate upper limits on the number of resolved point sources associated with the excess, and the fraction of the flux in the excess that can be attributed to resolved sources, into limits on the luminosity function parameter space. We discuss a number of important systematic uncertainties, including in the detection threshold model and the total flux attributed to the excess. We delineate regions of parameter space (containing existing benchmark models) where there is no apparent tension with current data, and the number of total pulsars needed to explain the excess is in the range of \U0001d4aa(104-5). In the future, lowered point source detection thresholds could be achieved either by new analysis methods or new data. An order-of-magnitude reduction in the sensitivity threshold (which may already be achieved by novel analyses probing sub-threshold source populations) could hope to resolve more than 30% of the flux of the excess even in pessimistic scenarios.

Neural simulation-based inference approach for characterizing the Galactic Center γ -ray excess

The nature of the Fermi gamma-ray Galactic Center Excess (GCE) has remained a persistent mystery for over a decade. Although the excess is broadly compatible with emission expected due to dark matter annihilation, an explanation in terms of a population of unresolved astrophysical point sources e.g., millisecond pulsars, remains viable. The effort to uncover the origin of the GCE is hampered in particular by an incomplete understanding of diffuse emission of Galactic origin. This can lead to spurious features that make it difficult to robustly differentiate smooth emission, as expected for a dark matter origin, from more "clumpy" emission expected for a population of relatively bright, unresolved point sources. We use recent advancements in the field of simulation-based inference, in particular density estimation techniques using normalizing flows, in order to characterize the contribution of modeled components, including unresolved point source populations, to the GCE. Compared to traditional techniques based on the statistical distribution of photon counts, our machine learning-based method is able to utilize more of the information contained in a given model of the Galactic Center emission, and in particular can perform posterior parameter estimation while accounting for pixel-to-pixel spatial correlations in the gamma-ray map. This makes the method demonstrably more resilient to certain forms of model misspecification. On application to Fermi data, the method generically attributes a smaller fraction of the GCE flux to unresolved point sources when compared to traditional approaches. We nevertheless infer such a contribution to make up a non-negligible fraction of the GCE across all analysis variations considered, with at least $38^{+9}_{-19}\%$ of the excess attributed to unresolved point sources in our baseline analysis.

Statistical strong lensing

Context.Existing samples of strong lenses have been assembled by giving priority to sample size, but this is often at the cost of a complex selection function. However, with the advent of the next generation of wide-field photometric surveys, it might become possible to identify subsets of the lens population with well-defined selection criteria, trading sample size for completeness.Aims.There are two main advantages of working with a complete sample of lenses. First, such completeness makes possible to recover the properties of the general population of galaxies, of which strong lenses are a biased subset. Second, the relative number of lenses and non-detections can be used to further constrain models of galaxy structure. The present work illustrates how to carry out a statistical strong lensing analysis that takes advantage of these features.Methods.I introduce a general formalism for the statistical analysis of a sample of strong lenses with known selection function, and then test it on simulated data. The simulation consists of a population of 105galaxies with an axisymmetric power-law density profile, a population of background point sources, and a subset of ∼103strong lenses, which form a complete sample above an observational cut.Results.The method allows the user to recover the distribution of the galaxy population in Einstein radius and mass density slope in an unbiased way. The number of non-lenses helps to constrain the model when magnification data are not available.Conclusions.Complete samples of lenses are a powerful asset with which to turn precise strong lensing measurements into accurate statements on the properties of the general galaxy population.

A Chandra Virgo cluster survey of spiral galaxies – I. Introduction to the survey and a new ULX sample

ABSTRACT We present an analysis of the ultraluminous X-ray source (ULX) population in 75 Virgo cluster late-type galaxies, including all those with a star formation rate ≳1 M⊙ yr−1 and a representative sample of the less star forming ones. This study is based on 110 observations obtained over 20 yr with the Chandra X-ray Observatory Advanced Camera for Imaging Spectroscopy. As part of a Large Chandra Program, new observations were obtained for 52 of these 75 galaxies. The data are complete to a sensitivity of ≈1039 erg s−1, with a typical detection limit of ≈3 × 1038 erg s−1 for the majority of the sources. The catalogue contains about 80 ULXs (0.3–10 keV luminosity >1039 erg s−1), and provides their location, observed flux, de-absorbed luminosity, and (for the 25 most luminous ones) simple X-ray spectral properties. We discuss the ULX luminosity function in relation to the mass and star formation rate of the sample galaxies. We show that recent models of low-mass plus high-mass X-ray binary populations (scaling with stellar mass and star formation rate, respectively) are mostly consistent with our observational results. We tentatively identify the most luminous X-ray source in the sample (a source in IC 3322A with LX ≈ 6 × 1040 erg s−1) as a recent supernova or its young remnant. The properties of the sample galaxies (morphologies, stellar masses, star formation rates, total X-ray luminosities from their point-source population) are also summarized.

The Origin of Intergalactic Light in Compact Groups of Galaxies

We investigate the origin of intergalactic light (IGL) in close groups of galaxies. IGL is hypothesized to be the byproduct of interaction and merger within compact groups. Comparing the X-ray point source population in our sample of compact groups that have intergalactic light with compact groups without IGL, we find marginal evidence for a small increase in ultra-luminous X-ray sources (ULXs). There is also a significant bias towards lower luminosity high mass X-ray binaries (HMXRBs). We interpret this as an indication that groups with visible IGL represent a later evolutionary phase than other compact groups. They have galaxies characterized by quenching of star formation (lower star formation rate (SFR) inferred from lower HMXRB luminosity) after stellar material has been removed from the galaxies into the intergalactic medium, which is the source of the IGL. We conclude that the presence of an increased fraction of ULXs is due to past interaction and mergers within groups that have IGL.

The Fermi–LAT Galactic Center Excess: Evidence of Annihilating Dark Matter?

The center of the Galaxy is one of the prime targets in the search for a signal of annihilating (or decaying) dark matter. If such a signal were to be detected, it would shed light on one of the biggest mysteries in physics today: What is dark matter? Fundamental properties of the particle nature of dark matter, such as its mass, annihilation cross section, and annihilation final states, could be measured for the first time. Several experiments have searched for such a signal, and some have measured excesses that are compatible with it. A long-standing and compelling excess is observed in γ-rays by the Fermi Large Area Telescope ( Fermi–LAT). This excess is consistent with a dark matter particle with a mass of approximately 50 (up to ∼200) GeV annihilating with a velocity-averaged cross section of ∼10−26 cm3 s−1. Although a dark matter origin of the excess remains viable, other interpretations are possible. In particular, there is some evidence that the excess is produced by a population of unresolved point sources of γ-rays—for example, millisecond pulsars. In this article, I review the current status of the observation of the Fermi–LAT Galactic center excess, the possible interpretations of the excess, the evidence and counterevidence for each, and the prospects for resolving its origin with future measurements.

Foreground mismodeling and the point source explanation of the Fermi Galactic Center excess

The Fermi Large Area Telescope has observed an excess of ~GeV energy gamma rays from the center of the Milky Way, which may arise from near-thermal dark matter annihilation. Firmly establishing the dark matter origin for this excess is however complicated by challenges in modeling diffuse cosmic-ray foregrounds as well as unresolved astrophysical sources, such as millisecond pulsars. Non-Poissonian Template Fitting (NPTF) is one statistical technique that has previously been used to show that at least some fraction of the GeV excess is likely due to a population of dim point sources. These results were recently called into question by Leane and Slatyer (2019), who showed that a synthetic dark matter annihilation signal injected on top of the real Fermi data is not recovered by the NPTF procedure. In this work, we perform a dedicated study of the Fermi data and explicitly show that the central result of Leane and Slatyer (2019) is likely driven by the fact that their choice of model for the Galactic foreground emission does not provide a sufficiently good description of the data. We repeat the NPTF analyses using a state-of-the-art model for diffuse gamma-ray emission in the Milky Way and introduce a novel statistical procedure, based on spherical-harmonic marginalization, to provide an improved description of the Galactic diffuse emission in a data-driven fashion. With these improvements, we find that the NPTF results continue to robustly favor the interpretation that the Galactic Center excess is due, in part, to unresolved astrophysical point sources across the analysis variations that we have explored.

The C-Band All-Sky Survey: total intensity point-source detection over the northern sky

ABSTRACT We present a point-source detection algorithm that employs the second-order Spherical Mexican Hat wavelet filter (SMHW2), and use it on C-Band All-Sky Survey (C-BASS) northern intensity data to produce a catalogue of point sources. This catalogue allows us to cross-check the C-BASS flux-density scale against existing source surveys, and provides the basis for a source mask that will be used in subsequent C-BASS and cosmic microwave background (CMB) analyses. The SMHW2 allows us to filter the entire sky at once, avoiding complications from edge effects arising when filtering small sky patches. The algorithm is validated against a set of Monte Carlo simulations, consisting of diffuse emission, instrumental noise, and various point-source populations. The simulated source populations are successfully recovered. The SMHW2 detection algorithm is used to produce a $4.76\, \mathrm{GHz}$ northern sky source catalogue in total intensity, containing 1784 sources and covering declinations δ ≥ −10°. The C-BASS catalogue is matched with the Green Bank 6 cm (GB6) and Parkes-MIT-NRAO (PMN) catalogues over their areas of common sky coverage. From this we estimate the 90 per cent completeness level to be approximately $610\, \mathrm{mJy}$, with a corresponding reliability of 98 per cent, when masking the brightest 30 per cent of the diffuse emission in the C-BASS northern sky map. We find the C-BASS and GB6 flux-density scales to be consistent with one another to within approximately 4 per cent.

A Search for Neutrino Point-source Populations in 7 yr of IceCube Data with Neutrino-count Statistics

Abstract The presence of a population of point sources in a data set modifies the underlying neutrino-count statistics from the Poisson distribution. This deviation can be exactly quantified using the non-Poissonian template fitting technique, and in this work we present the first application of this approach to the IceCube high-energy neutrino data set. Using this method, we search in 7 yr of IceCube data for point-source populations correlated with the disk of the Milky Way, the Fermi bubbles, the Schlegel, Finkbeiner, and Davis dust map, or with the isotropic extragalactic sky. No evidence for such a population is found in the data using this technique, and in the absence of a signal, we establish constraints on population models with source-count distribution functions that can be described by a power law with a single break. The derived limits can be interpreted in the context of many possible source classes. In order to enhance the flexibility of the results, we publish the full posterior from our analysis, which can be used to establish limits on specific population models that would contribute to the observed IceCube neutrino flux.

The ASKAP EMU Early Science Project: radio continuum survey of the Small Magellanic Cloud

ABSTRACT We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observed at 960 MHz (4489 sources) and 1320 MHz (5954 sources) with a bandwidth of 192 MHz and beam sizes of 30.0 × 30.0 arcsec2 and 16.3 × 15.1 arcsec2, respectively. The median root mean square (RMS) noise values are 186 $\mu$Jy beam−1 (960 MHz) and 165 $\mu$Jy beam−1 (1320 MHz). To create point source catalogues, we use these two source lists, together with the previously published Molonglo Observatory Synthesis Telescope (MOST) and the Australia Telescope Compact Array (ATCA) point source catalogues to estimate spectral indices for the whole population of radio point sources found in the survey region. Combining our ASKAP catalogues with these radio continuum surveys, we found 7736 point-like sources in common over an area of 30 deg2. In addition, we report the detection of two new, low surface brightness supernova remnant candidates in the SMC. The high sensitivity of the new ASKAP ESP survey also enabled us to detect the bright end of the SMC planetary nebula sample, with 22 out of 102 optically known planetary nebulae showing point-like radio continuum emission. Lastly, we present several morphologically interesting background radio galaxies.

Constraining Galactic dark matter with gamma-ray pixel counts statistics

Gamma-ray searches for new physics such as dark matter are often driven by investigating the composition of the extragalactic gamma-ray background (EGB). Classic approaches to EGB decomposition manifest in resolving individual point sources and dissecting the intensity spectrum of the remaining unresolved component. Furthermore, statistical methods have recently been proven to outperform the sensitivity of classic source detection algorithms in finding point-source populations in the unresolved flux regime. In this article, we employ the 1-point photon count statistics of eight years of Fermi-LAT data to resolve the population of extragalactic point sources and to decompose the diffuse isotropic background contribution for Galactic latitudes |b|>30 deg. We use three adjacent energy bins between 1 and 10 GeV. For the first time, we extend the analysis to incorporate a potential contribution from annihilating dark matter smoothly distributed in the Galaxy. We investigate the sensitivity reach of 1-point statistics for constraining the thermally-averaged self-annihilation cross section <\sigma v> of dark matter, using different template models for the Galactic foreground emission. Given the official Fermi-LAT interstellar emission model, we set upper bounds on the DM self-annihilation cross section <\sigma v> that are comparable with the constraints obtained by other indirect detection methods, in particular by the stacking analysis of several dwarf spheroidal galaxies.

Interplanetary scintillation studies with the Murchison Widefield Array III: comparison of source counts and densities for radio sources and their sub-arcsecond components at 162 MHz

We use Murchison Widefield Array observations of interplanetary scintillation (IPS) to determine the source counts of point ($<$0.3 arcsecond extent) sources and of all sources with some subarcsecond structure, at 162 MHz. We have developed the methodology to derive these counts directly from the IPS observables, while taking into account changes in sensitivity across the survey area. The counts of sources with compact structure follow the behaviour of the dominant source population above $\sim$3 Jy but below this they show Euclidean behaviour. We compare our counts to those predicted by simulations and find a good agreement for our counts of sources with compact structure, but significant disagreement for point source counts. Using low radio frequency SEDs from the GLEAM survey, we classify point sources as Compact Steep-Spectrum (CSS), flat spectrum, or peaked. If we consider the CSS sources to be the more evolved counterparts of the peaked sources, the two categories combined comprise approximately 80% of the point source population. We calculate densities of potential calibrators brighter than 0.4 Jy at low frequencies and find 0.2 sources per square degrees for point sources, rising to 0.7 sources per square degree if sources with more complex arcsecond structure are included. We extrapolate to estimate 4.6 sources per square degrees at 0.04 Jy. We find that a peaked spectrum is an excellent predictor for compactness at low frequencies, increasing the number of good calibrators by a factor of three compared to the usual flat spectrum criterion.