skip to main content

Title: Approximate Bayesian Computation applied to the Diffuse Gamma-Ray Sky

Many sources contribute to the diffuse gamma-ray background (DGRB), including star forming galaxies, active galactic nuclei, and cosmic ray interactions in the Milky Way. Exotic sources, such as dark matter annihilation, may also make some contribution. The photon counts-in-pixels distribution is a powerful tool for analysing the DGRB and determining the relative contributions of different sources. However, including photon energy information in a likelihood analysis of the counts-in-pixels distribution quickly becomes computationally intractable as the number of source types and energy bins increase. Here, we apply the likelihood-free method of approximate Bayesian computation (ABC) to the problem. We consider a mock analysis that includes contributions from dark matter annihilation in Galactic subhaloes as well as astrophysical backgrounds. We show that our results using ABC are consistent with the exact likelihood when energy information is discarded, and that significantly tighter parameter constraints can be obtained with ABC when energy information is included. ABC presents a powerful tool for analysing the DGRB and understanding its varied origins.

; ;
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 2326-2336
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. Ultra-high-energy (UHE) photons are an important tool for studying the high-energy Universe. A plausible source of photons with exa-eV (EeV) energy is provided by UHE cosmic rays (UHECRs) undergoing the Greisen–Zatsepin–Kuzmin process (Greisen 1966; Zatsepin & Kuzmin 1966) or pair production process (Blumenthal 1970) on a cosmic background radiation. In this context, the EeV photons can be a probe of both UHECR mass composition and the distribution of their sources (Gelmini, Kalashev & Semikoz 2008; Hooper, Taylor & Sarkar 2011). At the same time, the possible flux of photons produced by UHE protons in the vicinity of their sources by pion photoproduction or inelastic nuclear collisions would be noticeable only for relatively near sources, as the attenuation length of UHE photons is smaller than that of UHE protons; see, for example, Bhattacharjee & Sigl (2000) for a review. There also exists a class of so-called top-down models of UHECR generation that efficiently produce the UHE photons, for instance by the decay of heavy dark-matter particles (Berezinsky, Kachelriess & Vilenkin 1997; Kuzmin & Rubakov 1998) or by the radiation from cosmic strings (Berezinsky, Blasi & Vilenkin 1998). The search for the UHE photons was shown to be the most sensitive methodmore »of indirect detection of heavy dark matter (Kalashev & Kuznetsov 2016, 2017; Kuznetsov 2017; Kachelriess, Kalashev & Kuznetsov 2018; Alcantara, Anchordoqui & Soriano 2019). Another fundamental physics scenario that could be tested with UHE photons (Fairbairn, Rashba & Troitsky 2011) is the photon mixing with axion-like particles (Raffelt & Stodolsky 1988), which could be responsible for the correlation of UHECR events with BL Lac type objects observed by the High Resolution Fly’s Eye (HiRes) experiment (Gorbunov et al. 2004; Abbasi et al. 2006). In most of these scenarios, a clustering of photon arrival directions, rather than diffuse distribution, is expected, so point-source searches can be a suitable test for photon - axion-like particle mixing models. Finally, UHE photons could also be used as a probe for the models of Lorentz-invariance violation (Coleman & Glashow 1999; Galaverni & Sigl 2008; Maccione, Liberati & Sigl 2010; Rubtsov, Satunin & Sibiryakov 2012, 2014). The Telescope Array (TA; Tokuno et al. 2012; Abu-Zayyad et al. 2013c) is the largest cosmic ray experiment in the Northern Hemisphere. It is located at 39.3° N, 112.9° W in Utah, USA. The observatory includes a surface detector array (SD) and 38 fluorescence telescopes grouped into three stations. The SD consists of 507 stations that contain plastic scintillators, each with an area of 3 m2 (SD stations). The stations are placed in the square grid with 1.2 km spacing and cover an area of ∼700 km2. The TA SD is capable of detecting extensive air showers (EASs) in the atmosphere caused by cosmic particles of EeV and higher energies. The TA SD has been operating since 2008 May. A hadron-induced EAS significantly differs from an EAS induced by a photon because the depth of the shower maximum Xmax for a photon shower is larger, and a photon shower contains fewer muons and has a more curved front (see Risse & Homola 2007 for a review). The TA SD stations are sensitive to both muon and electromagnetic components of the shower and therefore can be triggered by both hadron-induced and photon-induced EAS events. In the present study, we use 9 yr of TA SD data for a blind search for point sources of UHE photons. We utilize the statistics of the SD data, which benefit from a high duty cycle. The full Monte Carlo (MC) simulation of proton-induced and photon-induced EAS events allows us to perform the photon search up to the highest accessible energies, E ≳ 1020 eV. As the main tool for the present photon search, we use a multivariate analysis based on a number of SD parameters that make it possible to distinguish between photon and hadron primaries. While searches for diffuse UHE photons were performed by several EAS experiments, including Haverah Park (Ave et al. 2000), AGASA (Shinozaki et al. 2002; Risse et al. 2005), Yakutsk (Rubtsov et al. 2006; Glushkov et al. 2007, 2010), Pierre Auger (Abraham et al. 2007, 2008a; Bleve 2016; Aab et al. 2017c) and TA (Abu-Zayyad et al. 2013b; Abbasi et al. 2019a), the search for point sources of UHE photons has been done only by the Pierre Auger Observatory (Aab et al. 2014, 2017a). The latter searches were based on hybrid data and were limited to the 1017.3 < E < 1018.5 eV energy range. In the present paper, we use the TA SD data alone. We perform the searches in five energy ranges: E > 1018, E > 1018.5, E > 1019, E > 1019.5 and E > 1020 eV. We find no significant evidence of photon point sources in all energy ranges and we set the point-source flux upper limits from each direction in the TA field of view (FOV). The search for unspecified neutral particles was also previously performed by the TA (Abbasi et al. 2015). The limit on the point-source flux of neutral particles obtained in that work is close to the present photon point-source flux limits.« less
  2. Abstract

    The thermal Sunyaev–Zel’dovich (tSZ) effect is a powerful tool with the potential for constraining directly the properties of the hot gas that dominates dark matter halos because it measures pressure and thus thermal energy density. Studying this hot component of the circumgalactic medium (CGM) is important because it is strongly impacted by star formation and active galactic nucleus (AGN) activity in galaxies, participating in the feedback loop that regulates star and black hole mass growth in galaxies. We study the tSZ effect across a wide halo-mass range using three cosmological hydrodynamical simulations: Illustris-TNG, EAGLE, and FIRE-2. Specifically, we present the scaling relation between the tSZ signal and halo mass and the (mass-weighted) radial profiles of gas density, temperature, and pressure for all three simulations. The analysis includes comparisons to Planck tSZ observations and to the thermal pressure profile inferred from the Atacama Cosmology Telescope (ACT) measurements. We compare these tSZ data to simulations to interpret the measurements in terms of feedback and accretion processes in the CGM. We also identify as-yet unobserved potential signatures of these processes that may be visible in future measurements, which will have the capability of measuring tSZ signals to even lower masses. We alsomore »perform internal comparisons between runs with different physical assumptions. We conclude (1) there is strong evidence for the impact of feedback atR500, but that this impact decreases by 5R500, and (2) the thermodynamic profiles of the CGM are highly dependent on the implemented model, such as cosmic-ray or AGN feedback prescriptions.

    « less

    Since the discovery of an excess in gamma rays in the direction of M31, its cause has been unclear. Published interpretations focus on dark matter or stellar related origins. Studies of a similar excess in the Milky Way centre motivate a correlation of the spatial morphology of the signal with the distribution of stellar mass in M31. However, a robust determination of the best theory for the observed excess emission is challenging due to uncertainties in the astrophysical gamma-ray foreground model. We perform a spectro-morphological analysis of the M31 gamma-ray excess using state-of-the-art templates for the distribution of stellar mass in M31 and novel astrophysical foreground models for its sky region. We construct maps for the old stellar populations of M31 based on data from the PAndAS survey and carefully remove the foreground stars. We also produce improved astrophysical foreground models via novel image inpainting techniques based on machine learning methods. Our stellar maps, mimicking the location of a population of millisecond pulsars in the bulge of M31, reach a 5.4σ significance, making them as strongly favoured as the simple phenomenological models usually considered in the literature, e.g. disc-like templates. This detection is robust to generous variations of themore »astrophysical foreground model. Once the stellar templates are included in the astrophysical model, we show that the dark matter annihilation interpretation of the signal is unwarranted. We demonstrate that about one million unresolved millisecond pulsars naturally explain the observed gamma-ray luminosity per stellar mass, energy spectrum, and stellar bulge-to-disc flux ratio.

    « less

    Ongoing and planned weak lensing (WL) surveys are becoming deep enough to contain information on angular scales down to a few arcmin. To fully extract information from these small scales, we must capture non-Gaussian features in the cosmological WL signal while accurately accounting for baryonic effects. In this work, we account for baryonic physics via a baryonic correction model that modifies the matter distribution in dark matter-only N-body simulations, mimicking the effects of galaxy formation and feedback. We implement this model in a large suite of ray-tracing simulations, spanning a grid of cosmological models in Ωm−σ8 space. We then develop a convolutional neural network (CNN) architecture to learn and constrain cosmological and baryonic parameters simultaneously from the simulated WL convergence maps. We find that in a Hyper-Suprime Cam-like survey, our CNN achieves a 1.7× tighter constraint in Ωm−σ8 space (1σ area) than the power spectrum and 2.1× tighter than the peak counts, showing that the CNN can efficiently extract non-Gaussian cosmological information even while marginalizing over baryonic effects. When we combine our CNN with the power spectrum, the baryonic effects degrade the constraint in Ωm−σ8 space by a factor of 2.4, compared to the much worse degradation by amore »factor of 4.7 or 3.7 from either method alone.

    « less
  5. Context. Inferences about dark matter, dark energy, and the missing baryons all depend on the accuracy of our model of large-scale structure evolution. In particular, with cosmological simulations in our model of the Universe, we trace the growth of structure, and visualize the build-up of bigger structures from smaller ones and of gaseous filaments connecting galaxy clusters. Aims. Here we aim to reveal the complexity of the large-scale structure assembly process in great detail and on scales from tens of kiloparsecs up to more than 10 Mpc with new sensitive large-scale observations from the latest generation of instruments. We also aim to compare our findings with expectations from our cosmological model. Methods. We used dedicated SRG/eROSITA performance verification (PV) X-ray, ASKAP/EMU Early Science radio, and DECam optical observations of a ~15 deg 2 region around the nearby interacting galaxy cluster system A3391/95 to study the warm-hot gas in cluster outskirts and filaments, the surrounding large-scale structure and its formation process, the morphological complexity in the inner parts of the clusters, and the (re-)acceleration of plasma. We also used complementary Sunyaev-Zeldovich (SZ) effect data from the Planck survey and custom-made Galactic total (neutral plus molecular) hydrogen column density maps based onmore »the HI4PI and IRAS surveys. We relate the observations to expectations from cosmological hydrodynamic simulations from the Magneticum suite. Results. We trace the irregular morphology of warm and hot gas of the main clusters from their centers out to well beyond their characteristic radii, r 200 . Between the two main cluster systems, we observe an emission bridge on large scale and with good spatial resolution. This bridge includes a known galaxy group but this can only partially explain the emission. Most gas in the bridge appears hot, but thanks to eROSITA’s unique soft response and large field of view, we discover some tantalizing hints for warm, truly primordial filamentary gas connecting the clusters. Several matter clumps physically surrounding the system are detected. For the “Northern Clump,” we provide evidence that it is falling towards A3391 from the X-ray hot gas morphology and radio lobe structure of its central AGN. Moreover, the shapes of these X-ray and radio structures appear to be formed by gas well beyond the virial radius, r 100 , of A3391, thereby providing an indirect way of probing the gas in this elusive environment. Many of the extended sources in the field detected by eROSITA are also known clusters or new clusters in the background, including a known SZ cluster at redshift z = 1. We find roughly an order of magnitude more cluster candidates than the SPT and ACT surveys together in the same area. We discover an emission filament north of the virial radius of A3391 connecting to the Northern Clump. Furthermore, the absorption-corrected eROSITA surface brightness map shows that this emission filament extends south of A3395 and beyond an extended X-ray-emitting object (the “Little Southern Clump”) towards another galaxy cluster, all at the same redshift. The total projected length of this continuous warm-hot emission filament is 15 Mpc, running almost 4 degrees across the entire eROSITA PV observation field. The Northern and Southern Filament are each detected at >4 σ . The Planck SZ map additionally appears to support the presence of both new filaments. Furthermore, the DECam galaxy density map shows galaxy overdensities in the same regions. Overall, the new datasets provide impressive confirmation of the theoretically expected structure formation processes on the individual system level, including the surrounding warm-hot intergalactic medium distribution; the similarities of features found in a similar system in the Magneticum simulation are striking. Our spatially resolved findings show that baryons indeed reside in large-scale warm-hot gas filaments with a clumpy structure.« less