Showing 1–50 of 91 results for author: Weniger, C

Numerical Analysis of Resonant Axion-Photon Mixing: Part I

Authors: Estanis Utrilla Ginés, Dion Noordhuis, Christoph Weniger, Samuel J. Witte

Abstract: Many present-day axion searches attempt to probe the mixing of axions and photons, which occurs in the presence of an external magnetic field. While this process is well-understood in a number of simple and idealized contexts, a strongly varying or highly inhomogeneous background can impact the efficiency and evolution of the mixing in a non-trivial manner. In an effort to develop a generalized fr… ▽ More Many present-day axion searches attempt to probe the mixing of axions and photons, which occurs in the presence of an external magnetic field. While this process is well-understood in a number of simple and idealized contexts, a strongly varying or highly inhomogeneous background can impact the efficiency and evolution of the mixing in a non-trivial manner. In an effort to develop a generalized framework for analyzing axion-photon mixing in arbitrary systems, we focus in this work on directly solving the axion-modified form of Maxwell's equations across a simulation domain with a spatially varying background. We concentrate specifically on understanding resonantly enhanced axion-photon mixing in a highly magnetized plasma, which is a key ingredient for developing precision predictions of radio signals emanating from the magnetospheres of neutron stars. After illustrating the success and accuracy of our approach for simplified limiting cases, we compare our results with a number of analytic solutions recently derived to describe mixing in these systems. We find that our numerical method demonstrates a high level of agreement with one, but only one, of the published results. Interestingly, our method also recovers the mixing between the axion and magnetosonic-t and Alfvén modes; these modes cannot escape from the regions of dense plasma, but could non-trivially alter the dynamics in certain environments. Future work will focus on extending our calculations to study resonant mixing in strongly variable backgrounds, mixing in generalized media (beyond the strong magnetic field limit), and the mixing of photons with other light bosonic fields, such as dark photons. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 15 pages, 7 figures + appendix

Simulation-based inference for stochastic gravitational wave background data analysis

Authors: James Alvey, Uddipta Bhardwaj, Valerie Domcke, Mauro Pieroni, Christoph Weniger

Abstract: The next generation of space- and ground-based facilities promise to reveal an entirely new picture of the gravitational wave sky: thousands of galactic and extragalactic binary signals, as well as stochastic gravitational wave backgrounds (SGWBs) of unresolved astrophysical and possibly cosmological signals. These will need to be disentangled to achieve the scientific goals of experiments such as… ▽ More The next generation of space- and ground-based facilities promise to reveal an entirely new picture of the gravitational wave sky: thousands of galactic and extragalactic binary signals, as well as stochastic gravitational wave backgrounds (SGWBs) of unresolved astrophysical and possibly cosmological signals. These will need to be disentangled to achieve the scientific goals of experiments such as LISA, Einstein Telescope, or Cosmic Explorer. We focus on one particular aspect of this challenge: reconstructing an SGWB from (mock) LISA data. We demonstrate that simulation-based inference (SBI) - specifically truncated marginal neural ratio estimation (TMNRE) - is a promising avenue to overcome some of the technical difficulties and compromises necessary when applying more traditional methods such as Monte Carlo Markov Chains (MCMC). To highlight this, we show that we can reproduce results from traditional methods both for a template-based and agnostic search for an SGWB. Moreover, as a demonstration of the rich potential of SBI, we consider the injection of a population of low signal-to-noise ratio supermassive black hole transient signals into the data. TMNRE can implicitly marginalize over this complicated parameter space, enabling us to directly and accurately reconstruct the stochastic (and instrumental noise) contributions. We publicly release our TMNRE implementation in the form of the code saqqara. △ Less

Submitted 2 July, 2024; v1 submitted 14 September, 2023; originally announced September 2023.

Comments: 6 pages, 2 figures, 1 table + supplemental material. Code: saqqara available at https://github.com/PEREGRINE-GW/saqqara. Matches published version

Report number: CERN-TH-2023-167

Axion Clouds around Neutron Stars

Authors: Dion Noordhuis, Anirudh Prabhu, Christoph Weniger, Samuel J. Witte

Abstract: Recent work has shown that axions can be efficiently produced via non-stationary pair plasma discharges in the polar cap region of pulsars. Here, we point out that for axion masses $10^{-9} \, {\rm eV} \lesssim m_a \lesssim 10^{-4} \, \rm eV$, a sizable fraction of the sourced axion population will be gravitationally confined to the neutron star. These axions accumulate over astrophysical timescal… ▽ More Recent work has shown that axions can be efficiently produced via non-stationary pair plasma discharges in the polar cap region of pulsars. Here, we point out that for axion masses $10^{-9} \, {\rm eV} \lesssim m_a \lesssim 10^{-4} \, \rm eV$, a sizable fraction of the sourced axion population will be gravitationally confined to the neutron star. These axions accumulate over astrophysical timescales, thereby forming a dense `axion cloud' around the star. We argue that the existence of such a cloud, with densities reaching and potentially exceeding $\mathcal{O}(10^{22}) \, {\rm GeV \, cm^{-3}}$, is a generic expectation across a wide range of parameter space. For axion masses $m_a \gtrsim 10^{-7} \, \rm eV$, energy is primarily radiated from the axion cloud via resonant axion-photon mixing, generating a number of distinctive signatures that include: a sharp line in the radio spectrum of each pulsar (located at the axion mass, and with a percent-level width), and transient events arising from the reconfiguration of charge densities in the magnetosphere. While a deeper understanding of the systematic uncertainties in these systems is required, our current estimates suggest that existing radio telescopes could improve sensitivity to the axion-photon coupling by more than an order of magnitude. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 17 pages, 14 figures

Contrastive Neural Ratio Estimation for Simulation-based Inference

Authors: Benjamin Kurt Miller, Christoph Weniger, Patrick Forré

Abstract: Likelihood-to-evidence ratio estimation is usually cast as either a binary (NRE-A) or a multiclass (NRE-B) classification task. In contrast to the binary classification framework, the current formulation of the multiclass version has an intrinsic and unknown bias term, making otherwise informative diagnostics unreliable. We propose a multiclass framework free from the bias inherent to NRE-B at opt… ▽ More Likelihood-to-evidence ratio estimation is usually cast as either a binary (NRE-A) or a multiclass (NRE-B) classification task. In contrast to the binary classification framework, the current formulation of the multiclass version has an intrinsic and unknown bias term, making otherwise informative diagnostics unreliable. We propose a multiclass framework free from the bias inherent to NRE-B at optimum, leaving us in the position to run diagnostics that practitioners depend on. It also recovers NRE-A in one corner case and NRE-B in the limiting case. For fair comparison, we benchmark the behavior of all algorithms in both familiar and novel training regimes: when jointly drawn data is unlimited, when data is fixed but prior draws are unlimited, and in the commonplace fixed data and parameters setting. Our investigations reveal that the highest performing models are distant from the competitors (NRE-A, NRE-B) in hyperparameter space. We make a recommendation for hyperparameters distinct from the previous models. We suggest two bounds on the mutual information as performance metrics for simulation-based inference methods, without the need for posterior samples, and provide experimental results. This version corrects a minor implementation error in $γ$, improving results. △ Less

Submitted 4 July, 2024; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 11 pages. 34 pages with references and supplemental material. Accepted at NeurIPS 2022. Updated version corrects code implementation error and all experiments. Code at https://github.com/bkmi/cnre

One never walks alone: the effect of the perturber population on subhalo measurements in strong gravitational lenses

Authors: Adam Coogan, Noemi Anau Montel, Konstantin Karchev, Meiert W. Grootes, Francesco Nattino, Christoph Weniger

Abstract: Analyses of extended arcs in strong gravitational lensing images to date have constrained the properties of dark matter by measuring the parameters of one or two individual subhalos. However, since such analyses are reliant on likelihood-based methods like Markov-chain Monte Carlo or nested sampling, they require various compromises to the realism of lensing models for the sake of computational tr… ▽ More Analyses of extended arcs in strong gravitational lensing images to date have constrained the properties of dark matter by measuring the parameters of one or two individual subhalos. However, since such analyses are reliant on likelihood-based methods like Markov-chain Monte Carlo or nested sampling, they require various compromises to the realism of lensing models for the sake of computational tractability, such as ignoring the numerous other subhalos and line-of-sight halos in the system, assuming a particular form for the source model and requiring the noise to have a known likelihood function. Here we show that a simulation-based inference method called truncated marginal neural ratio estimation (TMNRE) makes it possible to relax these requirements by training neural networks to directly compute marginal posteriors for subhalo parameters from lensing images. By performing a set of inference tasks on mock data, we verify the accuracy of TMNRE and show it can compute posteriors for subhalo parameters marginalized over populations of hundreds of subhalos and line-of-sight halos, as well as lens and source uncertainties. We also find the MLP Mixer network works far better for such tasks than the convolutional architectures explored in other lensing analyses. Furthermore, we show that since TMNRE learns a posterior function it enables direct statistical checks that would be extremely expensive with likelihood-based methods. Our results show that TMNRE is well-suited for analyzing complex lensing data, and that the full subhalo and line-of-sight halo population must be included when measuring the properties of individual dark matter substructures. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 13 pages, 8 figures

Novel Constraints on Axions Produced in Pulsar Polar-Cap Cascades

Authors: Dion Noordhuis, Anirudh Prabhu, Samuel J. Witte, Alexander Y. Chen, Fábio Cruz, Christoph Weniger

Abstract: Axions can be copiously produced in localized regions of neutron star magnetospheres where the ambient plasma is unable to efficiently screen the induced electric field. As these axions stream away from the neutron star they can resonantly transition into photons, generating a large broadband contribution to the neutron star's intrinsic radio flux. In this work, we develop a comprehensive end-to-e… ▽ More Axions can be copiously produced in localized regions of neutron star magnetospheres where the ambient plasma is unable to efficiently screen the induced electric field. As these axions stream away from the neutron star they can resonantly transition into photons, generating a large broadband contribution to the neutron star's intrinsic radio flux. In this work, we develop a comprehensive end-to-end framework to model this process from the initial production of axions to the final detection of radio photons, and derive constraints on the axion-photon coupling, $g_{aγγ}$, using observations of 27 nearby pulsars. We study the modeling uncertainty in the sourced axion spectrum by comparing predictions from 2.5 dimensional particle-in-cell simulations with those derived using a semi-analytic model; these results show remarkable agreement, leading to constraints on the axion-photon coupling that typically differ by a factor of no more than $\sim 2$. The limits presented here are the strongest to date for axion masses $10^{-8} \, {\rm eV} \lesssim m_a \lesssim 10^{-5} \, {\rm eV}$, and crucially do not rely on the assumption that axions are dark matter. △ Less

Submitted 25 July, 2023; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: v2: Updated to match published version. Added new SM sections on analysis and uncertainties, updated plots, and corrected minor bugs and typos. v1: 5 pages, 2 figures + Supplementary Material

Snowmass Theory Frontier: Astrophysics and Cosmology

Authors: Daniel Green, Joshua T. Ruderman, Benjamin R. Safdi, Jessie Shelton, Ana Achúcarro, Peter Adshead, Yashar Akrami, Masha Baryakhtar, Daniel Baumann, Asher Berlin, Nikita Blinov, Kimberly K. Boddy, Malte Buschmann, Giovanni Cabass, Robert Caldwell, Emanuele Castorina, Thomas Y. Chen, Xingang Chen, William Coulton, Djuna Croon, Yanou Cui, David Curtin, Francis-Yan Cyr-Racine, Christopher Dessert, Keith R. Dienes , et al. (62 additional authors not shown)

Abstract: We summarize progress made in theoretical astrophysics and cosmology over the past decade and areas of interest for the coming decade. This Report is prepared as the TF09 "Astrophysics and Cosmology" topical group summary for the Theory Frontier as part of the Snowmass 2021 process. We summarize progress made in theoretical astrophysics and cosmology over the past decade and areas of interest for the coming decade. This Report is prepared as the TF09 "Astrophysics and Cosmology" topical group summary for the Theory Frontier as part of the Snowmass 2021 process. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 57 pages

Estimating the warm dark matter mass from strong lensing images with truncated marginal neural ratio estimation

Authors: Noemi Anau Montel, Adam Coogan, Camila Correa, Konstantin Karchev, Christoph Weniger

Abstract: Precision analysis of galaxy-galaxy strong gravitational lensing images provides a unique way of characterizing small-scale dark matter halos, and could allow us to uncover the fundamental properties of dark matter's constituents. Recently, gravitational imaging techniques made it possible to detect a few heavy subhalos. However, gravitational lenses contain numerous subhalos and line-of-sight hal… ▽ More Precision analysis of galaxy-galaxy strong gravitational lensing images provides a unique way of characterizing small-scale dark matter halos, and could allow us to uncover the fundamental properties of dark matter's constituents. Recently, gravitational imaging techniques made it possible to detect a few heavy subhalos. However, gravitational lenses contain numerous subhalos and line-of-sight halos, whose subtle imprint is extremely difficult to detect individually. Existing methods for marginalizing over this large population of sub-threshold perturbers to infer population-level parameters are typically computationally expensive, or require compressing observations into hand-crafted summary statistics, such as a power spectrum of residuals. Here, we present the first analysis pipeline to combine parametric lensing models and a recently-developed neural simulation-based inference technique called truncated marginal neural ratio estimation (TMNRE) to constrain the warm dark matter halo mass function cutoff scale directly from multiple lensing images. Through a proof-of-concept application to simulated data, we show that our approach enables empirically testable inference of the dark matter cutoff mass through marginalization over a large population of realistic perturbers that would be undetectable on their own, and over lens and source parameters uncertainties. To obtain our results, we combine the signal contained in a set of images with Hubble Space Telescope resolution. Our results suggest that TMNRE can be a powerful approach to put tight constraints on the mass of warm dark matter in the multi-keV regime, which will be relevant both for existing lensing data and in the large sample of lenses that will be delivered by near-future telescopes. △ Less

Submitted 14 November, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: 15 pages, 11 figures; Updated to match version published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2022; stac3215

Astrophysical and Cosmological Probes of Dark Matter

Authors: Kimberly K. Boddy, Mariangela Lisanti, Samuel D. McDermott, Nicholas L. Rodd, Christoph Weniger, Yacine Ali-Haïmoud, Malte Buschmann, Ilias Cholis, Djuna Croon, Adrienne L. Erickcek, Vera Gluscevic, Rebecca K. Leane, Siddharth Mishra-Sharma, Julian B. Muñoz, Ethan O. Nadler, Priyamvada Natarajan, Adrian Price-Whelan, Simona Vegetti, Samuel J. Witte

Abstract: While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for t… ▽ More While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade. △ Less

Submitted 12 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021, 24 pages, 4 figures, comments welcome

Extraterrestrial Axion Search with the Breakthrough Listen Galactic Center Survey

Authors: Joshua W. Foster, Samuel J. Witte, Matthew Lawson, Tim Linden, Vishal Gajjar, Christoph Weniger, Benjamin R. Safdi

Abstract: Axion dark matter (DM) may efficiently convert to photons in the magnetospheres of neutron stars (NSs), producing nearly monochromatic radio emission. This process is resonantly triggered when the plasma frequency induced by the underlying charge distribution approximately matches the axion mass. We search for evidence of this process using archival Green Bank Telescope data collected in a survey… ▽ More Axion dark matter (DM) may efficiently convert to photons in the magnetospheres of neutron stars (NSs), producing nearly monochromatic radio emission. This process is resonantly triggered when the plasma frequency induced by the underlying charge distribution approximately matches the axion mass. We search for evidence of this process using archival Green Bank Telescope data collected in a survey of the Galactic Center in the C-Band by the Breakthrough Listen project. While Breakthrough Listen aims to find signatures of extraterrestrial life in the radio band, we show that their high-frequency resolution spectral data of the Galactic Center region is ideal for searching for axion-photon transitions generated by the population of NSs in the inner pc of the Galaxy. We use data-driven models to capture the distributions and properties of NSs in the inner Galaxy and compute the expected radio flux from each NS using state-of-the-art ray tracing simulations. We find no evidence for axion DM and set leading constraints on the axion-photon coupling, excluding values down to the level $g_{a γγ} \sim 10^{-11}$ GeV$^{-1}$ for DM axions for masses between 15 and 35 $μ$eV. △ Less

Submitted 16 February, 2022; originally announced February 2022.

Comments: 9+23 pages, 3+17 figures

Report number: MIT-CTP/5398

Truncated Marginal Neural Ratio Estimation

Authors: Benjamin Kurt Miller, Alex Cole, Patrick Forré, Gilles Louppe, Christoph Weniger

Abstract: Parametric stochastic simulators are ubiquitous in science, often featuring high-dimensional input parameters and/or an intractable likelihood. Performing Bayesian parameter inference in this context can be challenging. We present a neural simulation-based inference algorithm which simultaneously offers simulation efficiency and fast empirical posterior testability, which is unique among modern al… ▽ More Parametric stochastic simulators are ubiquitous in science, often featuring high-dimensional input parameters and/or an intractable likelihood. Performing Bayesian parameter inference in this context can be challenging. We present a neural simulation-based inference algorithm which simultaneously offers simulation efficiency and fast empirical posterior testability, which is unique among modern algorithms. Our approach is simulation efficient by simultaneously estimating low-dimensional marginal posteriors instead of the joint posterior and by proposing simulations targeted to an observation of interest via a prior suitably truncated by an indicator function. Furthermore, by estimating a locally amortized posterior our algorithm enables efficient empirical tests of the robustness of the inference results. Since scientists cannot access the ground truth, these tests are necessary for trusting inference in real-world applications. We perform experiments on a marginalized version of the simulation-based inference benchmark and two complex and narrow posteriors, highlighting the simulator efficiency of our algorithm as well as the quality of the estimated marginal posteriors. △ Less

Submitted 26 October, 2021; v1 submitted 2 July, 2021; originally announced July 2021.

Comments: 10 pages. 27 pages with references and supplemental material. Implementation of experiments at https://github.com/bkmi/tmnre/. Ready-to-use implementation of underlying algorithm at https://github.com/undark-lab/swyft/. Accepted at NeurIPS 2021

Axion-Photon Conversion in Neutron Star Magnetospheres: The Role of the Plasma in the Goldreich-Julian Model

Authors: Samuel J. Witte, Dion Noordhuis, Thomas D. P. Edwards, Christoph Weniger

Abstract: The most promising indirect search for the existence of axion dark matter uses radio telescopes to look for narrow spectral lines generated from the resonant conversion of axions in the magnetospheres of neutron stars. Unfortunately, a large list of theoretical uncertainties has prevented this search strategy from being fully accepted as robust. In this work we attempt to address major outstanding… ▽ More The most promising indirect search for the existence of axion dark matter uses radio telescopes to look for narrow spectral lines generated from the resonant conversion of axions in the magnetospheres of neutron stars. Unfortunately, a large list of theoretical uncertainties has prevented this search strategy from being fully accepted as robust. In this work we attempt to address major outstanding questions related to the role and impact of the plasma, including: $(i)$ does refraction and reflection of radio photons in the magnetosphere induce strong inhomogeneities in the flux, $(ii)$ can refraction induce premature axion-photon de-phasing, $(iii)$ to what extent do photon-plasma interactions induce a broadening of the spectral line, $(iv)$ does the flux have a strong time dependence, and $(v)$ can radio photons sourced by axions be absorbed by the plasma. We present an end-to-end analysis pipeline based on ray-tracing that exploits a state-of-the-art auto-differentiation algorithm to propagate photons from the conversion surface to asymptotically large distances. Adopting a charge symmetric Goldreich-Julian model for the magnetosphere, we show that for reasonable parameters one should expect a strong anisotropy of the signal, refraction induced axion-photon de-phasing, significant line-broadening, a variable time-dependence of the flux, and, for large enough magnetic fields, anisotropic absorption. Our simulation code is flexible enough to serve as the basis for follow-up studies with a large range of magnetosphere models. △ Less

Submitted 4 July, 2022; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: v3: Fix "?" citation. v2: Accepted for publication PRD. Convergence improved, plots updated, minor bugs and typos corrected. No major changes to results/conclusions. v1: 9 pages + Appendix

Simulation-efficient marginal posterior estimation with swyft: stop wasting your precious time

Authors: Benjamin Kurt Miller, Alex Cole, Gilles Louppe, Christoph Weniger

Abstract: We present algorithms (a) for nested neural likelihood-to-evidence ratio estimation, and (b) for simulation reuse via an inhomogeneous Poisson point process cache of parameters and corresponding simulations. Together, these algorithms enable automatic and extremely simulator efficient estimation of marginal and joint posteriors. The algorithms are applicable to a wide range of physics and astronom… ▽ More We present algorithms (a) for nested neural likelihood-to-evidence ratio estimation, and (b) for simulation reuse via an inhomogeneous Poisson point process cache of parameters and corresponding simulations. Together, these algorithms enable automatic and extremely simulator efficient estimation of marginal and joint posteriors. The algorithms are applicable to a wide range of physics and astronomy problems and typically offer an order of magnitude better simulator efficiency than traditional likelihood-based sampling methods. Our approach is an example of likelihood-free inference, thus it is also applicable to simulators which do not offer a tractable likelihood function. Simulator runs are never rejected and can be automatically reused in future analysis. As functional prototype implementation we provide the open-source software package swyft. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: Accepted at Machine Learning and the Physical Sciences at NeurIPS 2020. Package: https://github.com/undark-lab/swyft/

Transient Radio Signatures from Neutron Star Encounters with QCD Axion Miniclusters

Authors: Thomas D. P. Edwards, Bradley J. Kavanagh, Luca Visinelli, Christoph Weniger

Abstract: The QCD axion is expected to form dense structures known as axion miniclusters if the Peccei-Quinn symmetry is broken after inflation. Miniclusters that have survived until today will interact with neutron stars (NSs) in the Milky Way to produce transient radio signals from axion-photon conversion in the NS magnetosphere. We quantify the properties of these encounters and find that they occur freq… ▽ More The QCD axion is expected to form dense structures known as axion miniclusters if the Peccei-Quinn symmetry is broken after inflation. Miniclusters that have survived until today will interact with neutron stars (NSs) in the Milky Way to produce transient radio signals from axion-photon conversion in the NS magnetosphere. We quantify the properties of these encounters and find that they occur frequently ($\mathcal{O}(1-100)\,\mathrm{day}^{-1}$); last between a day and a few months; are spatially clustered towards the Galactic center; and can reach observable fluxes. These radio transients are within reach of current generation telescopes and therefore offer a promising pathway to discovering QCD axion dark matter. △ Less

Submitted 22 October, 2021; v1 submitted 10 November, 2020; originally announced November 2020.

Comments: 11 pages, 5 figures. Code is available online at https://github.com/bradkav/axion-miniclusters/ (archived at https://doi.org/10.5281/zenodo.4006128 ). Joint submission with "Stellar Disruption of Axion Miniclusters in the Milky Way" in astro-ph (arXiv:2011.05377). Updated to match published version

Journal ref: Phys. Rev. Lett. 127, 131103 (2021)

Stellar Disruption of Axion Miniclusters in the Milky Way

Authors: Bradley J. Kavanagh, Thomas D. P. Edwards, Luca Visinelli, Christoph Weniger

Abstract: Axion miniclusters are dense bound structures of dark matter axions that are predicted to form in the post-inflationary Peccei-Quinn symmetry breaking scenario. Although dense, miniclusters can easily be perturbed or even become unbound by interactions with baryonic objects such as stars. Here, we characterize the spatial distribution and properties of miniclusters in the Milky Way (MW) today afte… ▽ More Axion miniclusters are dense bound structures of dark matter axions that are predicted to form in the post-inflationary Peccei-Quinn symmetry breaking scenario. Although dense, miniclusters can easily be perturbed or even become unbound by interactions with baryonic objects such as stars. Here, we characterize the spatial distribution and properties of miniclusters in the Milky Way (MW) today after undergoing these stellar interactions throughout their lifetime. We do this by performing a suite of Monte Carlo simulations which track the miniclusters' structure and, in particular, accounts for partial disruption and mass loss through successive interactions. We consider two density profiles - Navarro-Frenk-White (NFW) and Power-law (PL) - for the individual miniclusters in order to bracket the uncertainties on the minicluster population today due to their uncertain formation history. For our fiducial analysis at the Solar position, we find a survival probability of 99% for miniclusters with PL profiles and 46% for those with NFW profiles. Our work extends previous estimates of this local survival probability to the entire MW. We find that towards the Galactic center, the survival probabilities drop drastically. Although we present results for a particular initial halo mass function, our simulations can be easily recast to different models using the provided data and code. Finally, we comment on the impact of our results on lensing, direct, and indirect detection. △ Less

Submitted 22 December, 2021; v1 submitted 10 November, 2020; originally announced November 2020.

Comments: 31 pages, 18 figures. Code available online at https://github.com/bradkav/axion-miniclusters/ (archived at https://doi.org/10.5281/zenodo.4006128 ). Joint submission w/ "Transient Radio Signatures from Neutron Star Encounters with QCD Axion Miniclusters", arXiv:2011.05378 . v2: Minor changes and references added. v3: Updated to match version published in PRD

Journal ref: Phys. Rev. D 104, 063038 (2021)

Targeted Likelihood-Free Inference of Dark Matter Substructure in Strongly-Lensed Galaxies

Authors: Adam Coogan, Konstantin Karchev, Christoph Weniger

Abstract: The analysis of optical images of galaxy-galaxy strong gravitational lensing systems can provide important information about the distribution of dark matter at small scales. However, the modeling and statistical analysis of these images is extraordinarily complex, bringing together source image and main lens reconstruction, hyper-parameter optimization, and the marginalization over small-scale str… ▽ More The analysis of optical images of galaxy-galaxy strong gravitational lensing systems can provide important information about the distribution of dark matter at small scales. However, the modeling and statistical analysis of these images is extraordinarily complex, bringing together source image and main lens reconstruction, hyper-parameter optimization, and the marginalization over small-scale structure realizations. We present here a new analysis pipeline that tackles these diverse challenges by bringing together many recent machine learning developments in one coherent approach, including variational inference, Gaussian processes, differentiable probabilistic programming, and neural likelihood-to-evidence ratio estimation. Our pipeline enables: (a) fast reconstruction of the source image and lens mass distribution, (b) variational estimation of uncertainties, (c) efficient optimization of source regularization and other hyperparameters, and (d) marginalization over stochastic model components like the distribution of substructure. We present here preliminary results that demonstrate the validity of our approach. △ Less

Submitted 27 November, 2020; v1 submitted 14 October, 2020; originally announced October 2020.

Comments: Accepted for the NeurIPS 2020 workshop Machine Learning and the Physical Sciences; 9 pages, 4 figures; minor edits for clarity since v1

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

Authors: The Cherenkov Telescope Array Consortium, :, A. Acharyya, R. Adam, C. Adams, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, L. Amati, G. Ambrosi, E. O. Angüner, L. A. Antonelli, C. Aramo, A. Araudo, T. Armstrong, F. Arqueros, K. Asano, Y. Ascasíbar, M. Ashley, C. Balazs, O. Ballester , et al. (427 additional authors not shown)

Abstract: We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models giv… ▽ More We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]" △ Less

Submitted 30 January, 2021; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 68 pages (including references) and 26 figures; text identical to the version published in JCAP

Journal ref: JCAP01(2021)057

Green Bank and Effelsberg Radio Telescope Searches for Axion Dark Matter Conversion in Neutron Star Magnetospheres

Authors: Joshua W. Foster, Yonatan Kahn, Oscar Macias, Zhiquan Sun, Ralph P. Eatough, Vladislav I. Kondratiev, Wendy M. Peters, Christoph Weniger, Benjamin R. Safdi

Abstract: Axion dark matter (DM) may convert to radio-frequency electromagnetic radiation in the strong magnetic fields around neutron stars. The radio signature of such a process would be an ultra-narrow spectral peak at a frequency determined by the mass of the axion particle. We analyze data we collected from the Robert C. Byrd Green Bank Telescope in the L-band and the Effelsberg 100-m Telescope in the… ▽ More Axion dark matter (DM) may convert to radio-frequency electromagnetic radiation in the strong magnetic fields around neutron stars. The radio signature of such a process would be an ultra-narrow spectral peak at a frequency determined by the mass of the axion particle. We analyze data we collected from the Robert C. Byrd Green Bank Telescope in the L-band and the Effelsberg 100-m Telescope in the L-Band and S-band from a number of sources expected to produce bright signals of axion-photon conversion, including the Galactic Center of the Milky Way and the nearby isolated neutron stars RX J0720.4-3125 and RX J0806.4-4123. We find no evidence for axion DM and are able to set some of the strongest constraints to-date on the existence of axion DM in the highly-motivated mass range between ~5-11 $μ$eV. △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: 7+20 pages, 2+17 figures, Supplementary Data at http://github.com/joshwfoster/RadioAxionSearch

Report number: LCTP-20-04

Journal ref: Phys. Rev. Lett. 125, 171301 (2020)

Radio Signal of Axion-Photon Conversion in Neutron Stars: A Ray Tracing Analysis

Authors: Mikaël Leroy, Marco Chianese, Thomas D. P. Edwards, Christoph Weniger

Abstract: Axion dark matter can resonantly convert into photons in the magnetospheres of neutron stars (NSs). It has recently been shown that radio observations of nearby NSs can therefore provide a highly sensitive probe of the axion parameter space. Here we extend existing calculations by performing the first three-dimensional computation of the photon flux, taking into account the isotropic phase-space d… ▽ More Axion dark matter can resonantly convert into photons in the magnetospheres of neutron stars (NSs). It has recently been shown that radio observations of nearby NSs can therefore provide a highly sensitive probe of the axion parameter space. Here we extend existing calculations by performing the first three-dimensional computation of the photon flux, taking into account the isotropic phase-space distribution of axions and the structure of the NS magnetosphere. In particular, we study the overall magnitude of the flux and its possible time variation. We find that overall signal strength is robust to our more realistic analysis. In addition, we find that the variance of the signal with respect to the NS rotation is washed out by the additional trajectories in our treatment. Nevertheless, we show that SKA observations towards J0806.4-4123 are sensitive to $g_{aγγ}\sim 3\times10^{-13}\mathrm{\,GeV^{-1}}$ at $m_a\sim 7\times10^{-6}\mathrm{\,eV}$, even when accounting for Doppler broadening. Finally, we provide the necessary code to calculate the photon flux for any given NS system. △ Less

Submitted 2 June, 2020; v1 submitted 18 December, 2019; originally announced December 2019.

Comments: 12 pages, 6 figures

Journal ref: Phys. Rev. D 101, 123003 (2020)

Differentiable Strong Lensing: Uniting Gravity and Neural Nets through Differentiable Probabilistic Programming

Authors: Marco Chianese, Adam Coogan, Paul Hofma, Sydney Otten, Christoph Weniger

Abstract: Since upcoming telescopes will observe thousands of strong lensing systems, creating fully-automated analysis pipelines for these images becomes increasingly important. In this work, we make a step towards that direction by developing the first end-to-end differentiable strong lensing pipeline. Our approach leverages and combines three important computer science developments: (a) convolutional neu… ▽ More Since upcoming telescopes will observe thousands of strong lensing systems, creating fully-automated analysis pipelines for these images becomes increasingly important. In this work, we make a step towards that direction by developing the first end-to-end differentiable strong lensing pipeline. Our approach leverages and combines three important computer science developments: (a) convolutional neural networks, (b) efficient gradient-based sampling techniques, and (c) deep probabilistic programming languages. The latter automatize parameter inference and enable the combination of generative deep neural networks and physics components in a single model. In the current work, we demonstrate that it is possible to combine a convolutional neural network trained on galaxy images as a source model with a fully-differentiable and exact implementation of gravitational lensing physics in a single probabilistic model. This does away with hyperparameter tuning for the source model, enables the simultaneous optimization of nearly one hundred source and lens parameters with gradient-based methods, and allows the use of efficient gradient-based posterior sampling techniques. These features make this automated inference pipeline potentially suitable for processing a large amount of data. By analyzing mock lensing systems with different signal-to-noise ratios, we show that lensing parameters are reconstructed with percent-level accuracy. More generally, we consider this work as one of the first steps in establishing differentiable probabilistic programming techniques in the particle astrophysics community, which have the potential to significantly accelerate and improve many complex data analysis tasks. △ Less

Submitted 2 June, 2020; v1 submitted 14 October, 2019; originally announced October 2019.

Comments: 13 pages, 9 figures. Updated to match version published in MNRAS

A frequentist analysis of three right-handed neutrinos with GAMBIT

Authors: Marcin Chrzaszcz, Marco Drewes, Tomas E. Gonzalo, Julia Harz, Suraj Krishnamurthy, Christoph Weniger

Abstract: The extension of the Standard Model by right-handed neutrinos can not only explain the active neutrino masses via the seesaw mechanism, it is also able solve a number of long standing problems in cosmology. Especially, masses below the TeV scale are of particular interest as they can lead to a plethora of signatures in experimental searches. We present the first full frequentist analysis of the ex… ▽ More The extension of the Standard Model by right-handed neutrinos can not only explain the active neutrino masses via the seesaw mechanism, it is also able solve a number of long standing problems in cosmology. Especially, masses below the TeV scale are of particular interest as they can lead to a plethora of signatures in experimental searches. We present the first full frequentist analysis of the extension of the Standard Model by three right-handed neutrinos, with masses between 60 MeV and 500 GeV, using the Global and Modular BSM (beyond the Standard Model) Inference Tool GAMBIT. Our analysis is based on the Casas-Ibarra parametrisation and includes a large range of experimental constraints: active neutrino mixing, indirect constraints from, e.g., electroweak precision observables and lepton universality, and numerous direct searches for right-handed neutrinos. To study their overall effect, we derive combined profile likelihood results for the phenomenologically most relevant parameter projections. Furthermore, we discuss the role of (marginally) statistically preferred regions in the parameter space. Finally, we explore the flavour mixing pattern of the three right-handed neutrinos for different values of the lightest neutrino mass. Our results comprise the most comprehensive assessment of the model with three right-handed neutrinos model below the TeV scale so far, and provide a robust ground for exploring the impact of future constraints or detections. △ Less

Submitted 3 July, 2020; v1 submitted 6 August, 2019; originally announced August 2019.

Comments: 56 pages, 40 figures, 22 tables

Report number: gambit-physics-2019

Journal ref: Eur.Phys.J.C 80 (2020) 6, 569

Gravitational wave probes of dark matter: challenges and opportunities

Authors: Gianfranco Bertone, Djuna Croon, Mustafa A. Amin, Kimberly K. Boddy, Bradley J. Kavanagh, Katherine J. Mack, Priyamvada Natarajan, Toby Opferkuch, Katelin Schutz, Volodymyr Takhistov, Christoph Weniger, Tien-Tien Yu

Abstract: In this white paper, we discuss the prospects for characterizing and identifying dark matter using gravitational waves, covering a wide range of dark matter candidate types and signals. We argue that present and upcoming gravitational wave probes offer unprecedented opportunities for unraveling the nature of dark matter and we identify the most urgent challenges and open problems with the aim of e… ▽ More In this white paper, we discuss the prospects for characterizing and identifying dark matter using gravitational waves, covering a wide range of dark matter candidate types and signals. We argue that present and upcoming gravitational wave probes offer unprecedented opportunities for unraveling the nature of dark matter and we identify the most urgent challenges and open problems with the aim of encouraging a strong community effort at the interface between these two exciting fields of research. △ Less

Submitted 23 September, 2019; v1 submitted 24 July, 2019; originally announced July 2019.

Comments: 17 pages, 2 figures. Comments welcome. v2: Added references and minor corrections

Journal ref: SciPost Phys. Core 3, 007 (2020)

Paleo-Detectors for Galactic Supernova Neutrinos

Authors: Sebastian Baum, Thomas D. P. Edwards, Bradley J. Kavanagh, Patrick Stengel, Andrzej K. Drukier, Katherine Freese, Maciej Górski, Christoph Weniger

Abstract: Paleo-detectors are a proposed experimental technique in which one would search for traces of recoiling nuclei in ancient minerals. Natural minerals on Earth are as old as $\mathcal{O}(1)\,$Gyr and, in many minerals, the damage tracks left by recoiling nuclei are also preserved for timescales long compared to 1 Gyr once created. Thus, even reading out relatively small target samples of order 100 g… ▽ More Paleo-detectors are a proposed experimental technique in which one would search for traces of recoiling nuclei in ancient minerals. Natural minerals on Earth are as old as $\mathcal{O}(1)\,$Gyr and, in many minerals, the damage tracks left by recoiling nuclei are also preserved for timescales long compared to 1 Gyr once created. Thus, even reading out relatively small target samples of order 100 g, paleo-detectors would allow one to search for very rare events thanks to the large exposure, $\varepsilon \sim 100\,{\rm g}\,{\rm Gyr} = 10^5\,{\rm t}\,{\rm yr}$. Here, we explore the potential of paleo-detectors to measure nuclear recoils induced by neutrinos from galactic core collapse supernovae. We find that they would not only allow for a direct measurement of the average core collapse supernova rate in the Milky Way, but would also contain information about the time-dependence of the local supernova rate over the past $\sim$1 Gyr. Since the supernova rate is thought to be directly proportional to the star formation rate, such a measurement would provide a determination of the local star formation history. We investigate the sensitivity of paleo-detectors to both a smooth time evolution and an enhancement of the core collapse supernova rate on relatively short timescales, as would be expected for a starburst period in the local group. △ Less

Submitted 13 May, 2020; v1 submitted 13 June, 2019; originally announced June 2019.

Comments: 15 pages, 5 figures. code available at https://github.com/tedwards2412/SN-paleology (archived at https://doi.org/10.5281/zenodo.3066206 ) v2: matches the published version

Report number: NORDITA-2019-060; LCTP-19-12

Journal ref: Phys. Rev. D 101, 103017 (2020)

A Unique Multi-Messenger Signal of QCD Axion Dark Matter

Authors: Thomas D. P. Edwards, Marco Chianese, Bradley J. Kavanagh, Samaya M. Nissanke, Christoph Weniger

Abstract: We propose a multi-messenger probe of QCD axion Dark Matter based on observations of black hole-neutron star binary inspirals. It is suggested that a dense Dark Matter spike may grow around intermediate mass black holes ($10^{3}-10^{5} \mathrm{\,M_{\odot}}$). The presence of such a spike produces two unique effects: a distinct phase shift in the gravitational wave strain during the inspiral and an… ▽ More We propose a multi-messenger probe of QCD axion Dark Matter based on observations of black hole-neutron star binary inspirals. It is suggested that a dense Dark Matter spike may grow around intermediate mass black holes ($10^{3}-10^{5} \mathrm{\,M_{\odot}}$). The presence of such a spike produces two unique effects: a distinct phase shift in the gravitational wave strain during the inspiral and an enhancement of the radio emission due to the resonant axion-photon conversion occurring in the neutron star magnetosphere throughout the inspiral and merger. Remarkably, the observation of the gravitational wave signal can be used to infer the Dark Matter density and, consequently, to predict the radio emission. We study the projected reach of the LISA interferometer and next-generation radio telescopes such as the Square Kilometre Array. Given a sufficiently nearby system, such observations will potentially allow for the detection of QCD axion Dark Matter in the mass range $10^{-7}\,\mathrm{eV}$ to $10^{-5}\,\mathrm{eV}$. △ Less

Submitted 22 April, 2020; v1 submitted 12 May, 2019; originally announced May 2019.

Comments: 5 pages, 3 figures. Appendix added with additional figures. Updated to published version

Journal ref: Phys. Rev. Lett. 124, 161101 (2020)

Primordial Black Holes as Silver Bullets for New Physics at the Weak Scale

Authors: Gianfranco Bertone, Adam Coogan, Daniele Gaggero, Bradley J. Kavanagh, Christoph Weniger

Abstract: Observational constraints on gamma rays produced by the annihilation of weakly interacting massive particles around primordial black holes (PBHs) imply that these two classes of Dark Matter candidates cannot coexist. We show here that the successful detection of one or more PBHs by radio searches (with the Square Kilometer Array) and gravitational waves searches (with LIGO/Virgo and the upcoming E… ▽ More Observational constraints on gamma rays produced by the annihilation of weakly interacting massive particles around primordial black holes (PBHs) imply that these two classes of Dark Matter candidates cannot coexist. We show here that the successful detection of one or more PBHs by radio searches (with the Square Kilometer Array) and gravitational waves searches (with LIGO/Virgo and the upcoming Einstein Telescope) would set extraordinarily stringent constraints on virtually all weak-scale extensions of the Standard Model with stable relics, including those predicting a WIMP abundance much smaller than that of Dark Matter. Upcoming PBHs searches have in particular the potential to rule out almost the entire parameter space of popular theories such as the minimal supersymmetric standard model and scalar singlet Dark Matter. △ Less

Submitted 7 January, 2020; v1 submitted 3 May, 2019; originally announced May 2019.

Comments: 10 pages, 3 figures. Code available at https://github.com/adam-coogan/pbhs_vs_wimps , archived at https://zenodo.org/badge/latestdoi/169754838 . v2: Matches version published in PRD

Journal ref: Phys. Rev. D 100, 123013 (2019)

Digging for Dark Matter: Spectral Analysis and Discovery Potential of Paleo-Detectors

Authors: Thomas D. P. Edwards, Bradley J. Kavanagh, Christoph Weniger, Sebastian Baum, Andrzej K. Drukier, Katherine Freese, Maciej Górski, Patrick Stengel

Abstract: Paleo-detectors are a recently proposed method for the direct detection of Dark Matter (DM). In such detectors, one would search for the persistent damage features left by DM--nucleus interactions in ancient minerals. Initial sensitivity projections have shown that paleo-detectors could probe much of the remaining Weakly Interacting Massive Particle (WIMP) parameter space. In this paper, we improv… ▽ More Paleo-detectors are a recently proposed method for the direct detection of Dark Matter (DM). In such detectors, one would search for the persistent damage features left by DM--nucleus interactions in ancient minerals. Initial sensitivity projections have shown that paleo-detectors could probe much of the remaining Weakly Interacting Massive Particle (WIMP) parameter space. In this paper, we improve upon the cut-and-count approach previously used to estimate the sensitivity by performing a full spectral analysis of the background- and DM-induced signal spectra. We consider two scenarios for the systematic errors on the background spectra: i) systematic errors on the normalization only, and ii) systematic errors on the shape of the backgrounds. We find that the projected sensitivity is rather robust to imperfect knowledge of the backgrounds. Finally, we study how well the parameters of the true WIMP model could be reconstructed in the hypothetical case of a WIMP discovery. △ Less

Submitted 26 March, 2019; v1 submitted 26 November, 2018; originally announced November 2018.

Comments: 14 pages, 5 figures, code available at https://github.com/tedwards2412/paleo_detectors/ . v2: Added additional analysis theory details, matches version published in PRD

Report number: NORDITA-2018-119; LCTP-18-26

Journal ref: Phys. Rev. D 99, 043541 (2019)

Axion global fits with Peccei-Quinn symmetry breaking before inflation using GAMBIT

Authors: Sebastian Hoof, Felix Kahlhoefer, Pat Scott, Christoph Weniger, Martin White

Abstract: We present global fits of cosmologically stable axion-like particle and QCD axion models in the mass range 0.1 neV to 10 eV. We focus on the case where the Peccei-Quinn symmetry is broken before the end of inflation, such that the initial value of the axion field can be considered to be homogeneous throughout the visible Universe. We include detailed likelihood functions from light-shining-through… ▽ More We present global fits of cosmologically stable axion-like particle and QCD axion models in the mass range 0.1 neV to 10 eV. We focus on the case where the Peccei-Quinn symmetry is broken before the end of inflation, such that the initial value of the axion field can be considered to be homogeneous throughout the visible Universe. We include detailed likelihood functions from light-shining-through-wall experiments, haloscopes, helioscopes, the axion relic density, horizontal branch stars, supernova 1987A, white dwarf cooling, and gamma-ray observations. We carry out both frequentist and Bayesian analyses, with and without the inclusion of white dwarf cooling. We explore the degree of fine-tuning present in different models and identify parameter regions where it is possible for QCD axion models to account for both the dark matter in the Universe and the cooling hints, comparing them to specific DFSZ- and KSVZ-type models. We find the most credible parameter regions, allowing us to set (prior-dependent) upper and lower bounds on the axion mass. Our analysis also suggests that QCD axions in this scenario most probably make up a non-negligible but sub-dominant component of the dark matter in the Universe. △ Less

Submitted 17 October, 2019; v1 submitted 16 October, 2018; originally announced October 2018.

Comments: 77 pages, 29 figures. Version v3: revised version of the paper due to a bug in the code (see Footnote 2); please use GAMBIT 1.3.1 or newer. Updated supplementary material available on Zenodo (version v3) at http://dx.doi.org/10.5281/zenodo.2615235

Report number: gambit-physics-2018

Journal ref: J. High Energ. Phys. 03 (2019) 191

Fundamental Physics with the Square Kilometre Array

Authors: A. Weltman, P. Bull, S. Camera, K. Kelley, H. Padmanabhan, J. Pritchard, A. Raccanelli, S. Riemer-Sørensen, L. Shao, S. Andrianomena, E. Athanassoula, D. Bacon, R. Barkana, G. Bertone, C. Bonvin, A. Bosma, M. Brüggen, C. Burigana, C. Bœhm, F. Calore, J. A. R. Cembranos, C. Clarkson, R. M. T. Connors, Á. de la Cruz-Dombriz, P. K. S. Dunsby , et al. (28 additional authors not shown)

Abstract: The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA a… ▽ More The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime. △ Less

Submitted 22 October, 2019; v1 submitted 5 October, 2018; originally announced October 2018.

Comments: 69 pages, PASA, in press

Journal ref: Publ. Astron. Soc. Aust. 37 (2020) e002

The Lyman-$α$ forest as a diagnostic of the nature of the dark matter

Authors: Antonella Garzilli, Andrii Magalich, Tom Theuns, Carlos S. Frenk, Christoph Weniger, Oleg Ruchayskiy, Alexey Boyarsky

Abstract: The observed Lyman-$α$ flux power spectrum (FPS) is suppressed on scales below $\sim~ 30~{\rm km~s}^{-1}$. This cutoff could be due to the high temperature, $T_0$, and pressure, $p_0$, of the absorbing gas or, alternatively, it could reflect the free streaming of dark matter particles in the early universe. We perform a set of very high resolution cosmological hydrodynamic simulations in which we… ▽ More The observed Lyman-$α$ flux power spectrum (FPS) is suppressed on scales below $\sim~ 30~{\rm km~s}^{-1}$. This cutoff could be due to the high temperature, $T_0$, and pressure, $p_0$, of the absorbing gas or, alternatively, it could reflect the free streaming of dark matter particles in the early universe. We perform a set of very high resolution cosmological hydrodynamic simulations in which we vary $T_0$, $p_0$ and the amplitude of the dark matter free streaming, and compare the FPS of mock spectra to the data. We show that the location of the dark matter free-streaming cutoff scales differently with redshift than the cutoff produced by thermal effects and is more pronounced at higher redshift. We, therefore, focus on a comparison to the observed FPS at $z>5$. We demonstrate that the FPS cutoff can be fit assuming cold dark matter, but it can be equally well fit assuming that the dark matter consists of $\sim 7$ keV sterile neutrinos in which case the cutoff is due primarily to the dark matter free streaming. △ Less

Submitted 7 August, 2019; v1 submitted 18 September, 2018; originally announced September 2018.

Comments: 16 pages, 11 figures

Statistical challenges in the search for dark matter

Authors: Sara Algeri, Melissa van Beekveld, Nassim Bozorgnia, Alyson Brooks, J. Alberto Casas, Jessi Cisewski-Kehe, Francis-Yan Cyr-Racine, Thomas D. P. Edwards, Fabio Iocco, Bradley J. Kavanagh, Judita Mamužić, Siddharth Mishra-Sharma, Wolfgang Rau, Roberto Ruiz de Austri, Benjamin R. Safdi, Pat Scott, Tracy R. Slatyer, Yue-Lin Sming Tsai, Aaron C. Vincent, Christoph Weniger, Jennifer Rittenhouse West, Robert L. Wolpert

Abstract: The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Albert… ▽ More The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Alberta. △ Less

Submitted 24 July, 2018; originally announced July 2018.

Comments: 32 pages, 8 figures, 331 references

Report number: IPPP/18/60

Dark Matter Model or Mass, but Not Both: Assessing Near-Future Direct Searches with Benchmark-free Forecasting

Authors: Thomas D. P. Edwards, Bradley J. Kavanagh, Christoph Weniger

Abstract: Forecasting the signal discrimination power of dark matter (DM) searches is commonly limited to a set of arbitrary benchmark points. We introduce new methods for benchmark-free forecasting that instead allow an exhaustive exploration and visualization of the phenomenological distinctiveness of DM models, based on standard hypothesis testing. Using this method, we reassess the signal discrimination… ▽ More Forecasting the signal discrimination power of dark matter (DM) searches is commonly limited to a set of arbitrary benchmark points. We introduce new methods for benchmark-free forecasting that instead allow an exhaustive exploration and visualization of the phenomenological distinctiveness of DM models, based on standard hypothesis testing. Using this method, we reassess the signal discrimination power of future liquid Xenon and Argon direct DM searches. We quantify the parameter regions where various non-relativistic effective operators, millicharged DM, and magnetic dipole DM can be discriminated, and where upper limits on the DM mass can be found. We find that including an Argon target substantially improves the prospects for reconstructing the DM properties. We also show that only in a small region with DM masses in the range 20-100 GeV and DM-nucleon cross sections a factor of a few below current bounds can near-future Xenon and Argon detectors discriminate both the DM-nucleon interaction and the DM mass simultaneously. In all other regions only one or the other can be obtained. △ Less

Submitted 13 November, 2018; v1 submitted 10 May, 2018; originally announced May 2018.

Comments: 5 pages, 4 figures + appendices. Code for the calculations throughout the paper can be found at https://github.com/tedwards2412/benchmark_free_forecasting/ . v2: Minor corrections, matches version published in PRL as "Assessing Near-Future Direct Dark Matter Searches with Benchmark-Free Forecasting"

Journal ref: Phys. Rev. Lett. 121, 181101 (2018)

swordfish: Efficient Forecasting of New Physics Searches without Monte Carlo

Authors: Thomas D. P. Edwards, Christoph Weniger

Abstract: We introduce swordfish, a Monte-Carlo-free Python package to predict expected exclusion limits, the discovery reach and expected confidence contours for a large class of experiments relevant for particle- and astrophysics. The tool is applicable to any counting experiment, supports general correlated background uncertainties, and gives exact results in both the signal- and systematics-limited regi… ▽ More We introduce swordfish, a Monte-Carlo-free Python package to predict expected exclusion limits, the discovery reach and expected confidence contours for a large class of experiments relevant for particle- and astrophysics. The tool is applicable to any counting experiment, supports general correlated background uncertainties, and gives exact results in both the signal- and systematics-limited regimes. Instead of time-intensive Monte Carlo simulations and likelihood maximization, it internally utilizes new approximation methods that are built on information geometry. Out of the box, swordfish provides straightforward methods for accurately deriving many of the common sensitivity measures. In addition, it allows one to examine experimental abilities in great detail by employing the notion of information flux. This new concept generalizes signal-to-noise ratios to situations where background uncertainties and component mixing cannot be neglected. The user interface of swordfish is designed with ease-of-use in mind, which we demonstrate by providing typical examples from indirect and direct dark matter searches as jupyter notebooks. △ Less

Submitted 14 December, 2017; originally announced December 2017.

Comments: 17 pages including appendix, 13 figures

SpecBit, DecayBit and PrecisionBit: GAMBIT modules for computing mass spectra, particle decay rates and precision observables

Authors: The GAMBIT Models Workgroup, :, Peter Athron, Csaba Balázs, Lars A. Dal, Joakim Edsjö, Ben Farmer, Tomás E. Gonzalo, Anders Kvellestad, James McKay, Antje Putze, Chris Rogan, Pat Scott, Christoph Weniger, Martin White

Abstract: We present the GAMBIT modules SpecBit, DecayBit and PrecisionBit. Together they provide a new framework for linking publicly available spectrum generators, decay codes and other precision observable calculations in a physically and statistically consistent manner. This allows users to automatically run various combinations of existing codes as if they are a single package. The modular design allow… ▽ More We present the GAMBIT modules SpecBit, DecayBit and PrecisionBit. Together they provide a new framework for linking publicly available spectrum generators, decay codes and other precision observable calculations in a physically and statistically consistent manner. This allows users to automatically run various combinations of existing codes as if they are a single package. The modular design allows software packages fulfilling the same role to be exchanged freely at runtime, with the results presented in a common format that can be easily passed to downstream dark matter, collider and flavour codes. These modules constitute an essential part of the broader GAMBIT framework, a major new software package for performing global fits. In this paper we present the observable calculations, data, and likelihood functions implemented in the three modules, as well as the conventions and assumptions used in interfacing them with external codes. We also present 3-BIT-HIT, a command-line utility for computing mass spectra, couplings, decays and precision observables in the MSSM, which shows how the three modules can be easily used independently of GAMBIT. △ Less

Submitted 17 January, 2019; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 67 pages, 9 figures, updated to match version published in EPJC

Report number: CoEPP-MN-17-7, gambit-code-2017

Journal ref: Eur. Phys. J. C (2018) 78: 22

Global fits of GUT-scale SUSY models with GAMBIT

Authors: The GAMBIT Collaboration, Peter Athron, Csaba Balázs, Torsten Bringmann, Andy Buckley, Marcin Chrząszcz, Jan Conrad, Jonathan M. Cornell, Lars A. Dal, Joakim Edsjö, Ben Farmer, Paul Jackson, Abram Krislock, Anders Kvellestad, Farvah Mahmoudi, Gregory D. Martinez, Antje Putze, Are Raklev, Christopher Rogan, Roberto Ruiz de Austri, Aldo Saavedra, Christopher Savage, Pat Scott, Nicola Serra, Christoph Weniger , et al. (1 additional authors not shown)

Abstract: We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, dir… ▽ More We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos. △ Less

Submitted 23 October, 2018; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 50 pages, 21 figures, 7 tables, v2 accepted for publication in EPJC, v3 update Zenodo link

Report number: CoEPP-MN-17-9, gambit-physics-2017

FlavBit: A GAMBIT module for computing flavour observables and likelihoods

Authors: Florian U. Bernlochner, Marcin Chrzaszcz, Lars A. Dal, Ben Farmer, Paul Jackson, Anders Kvellestad, Farvah Mahmoudi, Antje Putze, Christopher Rogan, Pat Scott, Nicola Serra, Christoph Weniger, Martin White

Abstract: Flavour physics observables are excellent probes of new physics up to very high energy scales. Here we present FlavBit, the dedicated flavour physics module of the global-fitting package GAMBIT. FlavBit includes custom implementations of various likelihood routines for a wide range of flavour observables, including detailed uncertainties and correlations associated with LHCb measurements of rare,… ▽ More Flavour physics observables are excellent probes of new physics up to very high energy scales. Here we present FlavBit, the dedicated flavour physics module of the global-fitting package GAMBIT. FlavBit includes custom implementations of various likelihood routines for a wide range of flavour observables, including detailed uncertainties and correlations associated with LHCb measurements of rare, leptonic and semileptonic decays of B and D mesons, kaons and pions. It provides a generalised interface to external theory codes such as SuperIso, allowing users to calculate flavour observables in and beyond the Standard Model, and then test them in detail against all relevant experimental data. We describe FlavBit and its constituent physics in some detail, then give examples from supersymmetry and effective field theory illustrating how it can be used both as a standalone library for flavour physics, and within GAMBIT. △ Less

Submitted 28 November, 2017; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 25 pages, 2 figures, 9 tables. v2: References added

Report number: CERN-TH-2017-167, NORDITA 2017-077, gambit-code-2017

Journal ref: Eur.Phys.J. C77 (2017) no.11, 786

Status of the scalar singlet dark matter model

Authors: The GAMBIT Collaboration, Peter Athron, Csaba Balázs, Torsten Bringmann, Andy Buckley, Marcin Chrząszcz, Jan Conrad, Jonathan M. Cornell, Lars A. Dal, Joakim Edsjö, Ben Farmer, Paul Jackson, Felix Kahlhoefer, Abram Krislock, Anders Kvellestad, James McKay, Farvah Mahmoudi, Gregory D. Martinez, Antje Putze, Are Raklev, Christopher Rogan, Aldo Saavedra, Christopher Savage, Pat Scott, Nicola Serra , et al. (2 additional authors not shown)

Abstract: One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection… ▽ More One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above $\sim$1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned. △ Less

Submitted 19 October, 2018; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 19 pages, 5 figures, matches published version. v3 updates the Zenodo URL for supplementary data (available at https://doi.org/10.5281/zenodo.801510)

Report number: CoEPP-MN-17-10, gambit-physics-2017

Journal ref: Eur. Phys. J. C 77, 568 (2017)

DarkBit: A GAMBIT module for computing dark matter observables and likelihoods

Authors: Torsten Bringmann, Jan Conrad, Jonathan M. Cornell, Lars A. Dal, Joakim Edsjö, Ben Farmer, Felix Kahlhoefer, Anders Kvellestad, Antje Putze, Christopher Savage, Pat Scott, Christoph Weniger, Martin White, Sebastian Wild

Abstract: We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages. This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well… ▽ More We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages. This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well as a dedicated likelihood calculator for current and future experiments (gamlike). This provides a general solution for studying complex particle physics models that predict dark matter annihilation to a multitude of final states. We also supply a direct detection package that models a large range of direct detection experiments (DDcalc), and provides the corresponding likelihoods for arbitrary combinations of spin-independent and spin-dependent scattering processes. Finally, we provide custom relic density routines along with interfaces to DarkSUSY, micrOMEGAs, and the neutrino telescope likelihood package nuLike. DarkBit is written in the framework of the Global And Modular Beyond the Standard Model Inference Tool (GAMBIT), providing seamless integration into a comprehensive statistical fitting framework that allows users to explore new models with both particle and astrophysics constraints, and a consistent treatment of systematic uncertainties. In this paper we describe its main functionality, provide a guide to getting started quickly, and show illustrative examples for results obtained with DarkBit (both as a standalone tool and as a GAMBIT module). This includes a quantitative comparison between two of the main dark matter codes (DarkSUSY and micrOMEGAs), and application of DarkBit's advanced direct and indirect detection routines to a simple effective dark matter model. △ Less

Submitted 8 January, 2018; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 53 pages, 9 figures, 25 tables. v2: Updated to reflect new functionality in version 1.1.0, including the addition of Xenon1T and PICO-60 likelihoods. Minor corrections and clarifications throughout. Matches published version

Report number: DESY-17-235, NORDITA-2017-076, gambit-code-2017

Journal ref: Eur.Phys.J. C77 (2017) no.12, 831

ColliderBit: a GAMBIT module for the calculation of high-energy collider observables and likelihoods

Authors: Csaba Balázs, Andy Buckley, Lars A. Dal, Ben Farmer, Paul Jackson, Abram Krislock, Anders Kvellestad, Daniel Murnane, Antje Putze, Are Raklev, Christopher Rogan, Aldo Saavedra, Pat Scott, Christoph Weniger, Martin White

Abstract: We describe ColliderBit, a new code for the calculation of high energy collider observables in theories of physics beyond the Standard Model (BSM). ColliderBit features a generic interface to BSM models, a unique parallelised Monte Carlo event generation scheme suitable for large-scale supercomputer applications, and a number of LHC analyses, covering a reasonable range of the BSM signatures curre… ▽ More We describe ColliderBit, a new code for the calculation of high energy collider observables in theories of physics beyond the Standard Model (BSM). ColliderBit features a generic interface to BSM models, a unique parallelised Monte Carlo event generation scheme suitable for large-scale supercomputer applications, and a number of LHC analyses, covering a reasonable range of the BSM signatures currently sought by ATLAS and CMS. ColliderBit also calculates likelihoods for Higgs sector observables, and LEP searches for BSM particles. These features are provided by a combination of new code unique to ColliderBit, and interfaces to existing state-of-the-art public codes. ColliderBit is both an important part of the GAMBIT framework for BSM inference, and a standalone tool for efficiently applying collider constraints to theories of new physics. △ Less

Submitted 2 April, 2019; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 31 pages, 7 figures, 17 tables. Replacement version accepted by EPJC

Report number: gambit-code-2017

A global fit of the MSSM with GAMBIT

Authors: The GAMBIT Collaboration, Peter Athron, Csaba Balázs, Torsten Bringmann, Andy Buckley, Marcin Chrząszcz, Jan Conrad, Jonathan M. Cornell, Lars A. Dal, Joakim Edsjö, Ben Farmer, Paul Jackson, Abram Krislock, Anders Kvellestad, Farvah Mahmoudi, Gregory D. Martinez, Antje Putze, Are Raklev, Christopher Rogan, Aldo Saavedra, Christopher Savage, Pat Scott, Nicola Serra, Christoph Weniger, Martin White

Abstract: We study the seven-dimensional Minimal Supersymmetric Standard Model (MSSM7) with the new GAMBIT software framework, with all parameters defined at the weak scale. Our analysis significantly extends previous weak-scale, phenomenological MSSM fits, by adding more and newer experimental analyses, improving the accuracy and detail of theoretical predictions, including dominant uncertainties from the… ▽ More We study the seven-dimensional Minimal Supersymmetric Standard Model (MSSM7) with the new GAMBIT software framework, with all parameters defined at the weak scale. Our analysis significantly extends previous weak-scale, phenomenological MSSM fits, by adding more and newer experimental analyses, improving the accuracy and detail of theoretical predictions, including dominant uncertainties from the Standard Model, the Galactic dark matter halo and the quark content of the nucleon, and employing novel and highly-efficient statistical sampling methods to scan the parameter space. We find regions of the MSSM7 that exhibit co-annihilation of neutralinos with charginos, stops and sbottoms, as well as models that undergo resonant annihilation via both light and heavy Higgs funnels. We find high-likelihood models with light charginos, stops and sbottoms that have the potential to be within the future reach of the LHC. Large parts of our preferred parameter regions will also be accessible to the next generation of direct and indirect dark matter searches, making prospects for discovery in the near future rather good. △ Less

Submitted 22 October, 2018; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 31 pages, 14 figures; v2 matches version accepted by EPJC; v3 has corrected link to supplementary data on Zenodo

Report number: CERN-TH-2017-169, CoEPP-MN-17-11, NORDITA 2017-081, gambit-physics-2017

Journal ref: Eur.Phys.J. C77 (2017) no.12, 879

GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool

Authors: The GAMBIT Collaboration, Peter Athron, Csaba Balazs, Torsten Bringmann, Andy Buckley, Marcin Chrząszcz, Jan Conrad, Jonathan M. Cornell, Lars A. Dal, Hugh Dickinson, Joakim Edsjö, Ben Farmer, Tomás E. Gonzalo, Paul Jackson, Abram Krislock, Anders Kvellestad, Johan Lundberg, James McKay, Farvah Mahmoudi, Gregory D. Martinez, Antje Putze, Are Raklev, Joachim Ripken, Christopher Rogan, Aldo Saavedra , et al. (7 additional authors not shown)

Abstract: We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-Standard-Model Inference Tool. GAMBIT combines extensive calculations of observables and likelihoods in particle and astroparticle physics with a hierarchical model database, advanced tools for automatically building analyses of essentially any model, a flexible and powerful system for interfacing to exter… ▽ More We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-Standard-Model Inference Tool. GAMBIT combines extensive calculations of observables and likelihoods in particle and astroparticle physics with a hierarchical model database, advanced tools for automatically building analyses of essentially any model, a flexible and powerful system for interfacing to external codes, a suite of different statistical methods and parameter scanning algorithms, and a host of other utilities designed to make scans faster, safer and more easily-extendible than in the past. Here we give a detailed description of the framework, its design and motivation, and the current models and other specific components presently implemented in GAMBIT. Accompanying papers deal with individual modules and present first GAMBIT results. GAMBIT can be downloaded from gambit.hepforge.org. △ Less

Submitted 3 November, 2017; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 84 pages, 8 Figures. Code available at http://gambit.hepforge.org. v2: Main paper now matches version accepted by EPJC. This version also adds an addendum submitted to EPJC on Mathematica support in GAMBIT v1.1

Report number: CERN-TH-2017-166, CoEPP-MN-17-6, NORDITA 2017-074, DESY 17-236, gambit-code-2017

A Fresh Approach to Forecasting in Astroparticle Physics and Dark Matter Searches

Authors: Thomas D. P. Edwards, Christoph Weniger

Abstract: We present a toolbox of new techniques and concepts for the efficient forecasting of experimental sensitivities. These are applicable to a large range of scenarios in (astro-)particle physics, and based on the Fisher information formalism. Fisher information provides an answer to the question what is the maximum extractable information from a given observation?. It is a common tool for the forecas… ▽ More We present a toolbox of new techniques and concepts for the efficient forecasting of experimental sensitivities. These are applicable to a large range of scenarios in (astro-)particle physics, and based on the Fisher information formalism. Fisher information provides an answer to the question what is the maximum extractable information from a given observation?. It is a common tool for the forecasting of experimental sensitivities in many branches of science, but rarely used in astroparticle physics or searches for particle dark matter. After briefly reviewing the Fisher information matrix of general Poisson likelihoods, we propose very compact expressions for estimating expected exclusion and discovery limits (equivalent counts method). We demonstrate by comparison with Monte Carlo results that they remain surprisingly accurate even deep in the Poisson regime. We show how correlated background systematics can be efficiently accounted for by a treatment based on Gaussian random fields. Finally, we introduce the novel concept of Fisher information flux. It can be thought of as a generalization of the commonly used signal-to-noise ratio, while accounting for the non-local properties and saturation effects of background and instrumental uncertainties. It is a powerful and flexible tool ready to be used as core concept for informed strategy development in astroparticle physics and searches for particle dark matter. △ Less

Submitted 6 February, 2018; v1 submitted 18 April, 2017; originally announced April 2017.

Comments: 33 pages, 12 figures

Search for right-handed neutrinos from dark matter annihilation with gamma-rays

Authors: Miguel D. Campos, Farinaldo S. Queiroz, Carlos E. Yaguna, Christoph Weniger

Abstract: Several extensions of the Standard Model contain right-handed (sterile) neutrinos in the GeV-TeV mass range. Due to their mixing with the active neutrinos, they may give rise to novel effects in cosmology, neutrino physics, and collider searches. In addition, right-handed neutrinos can also appear as final states from dark matter annihilations, with important implications for dark matter indirect… ▽ More Several extensions of the Standard Model contain right-handed (sterile) neutrinos in the GeV-TeV mass range. Due to their mixing with the active neutrinos, they may give rise to novel effects in cosmology, neutrino physics, and collider searches. In addition, right-handed neutrinos can also appear as final states from dark matter annihilations, with important implications for dark matter indirect detection searches. In this paper, we use current data from the Fermi Large Area Telescope (6-year observation of dwarf spheroidal galaxies) and H.E.S.S. (10-year observation of the Galactic center) to constrain the annihilation of dark matter into right-handed neutrinos. We consider right-handed neutrino with masses between 10 GeV and 1 TeV, including both two-body and three-body decays, to derive bounds on the dark matter annihilation rate, $\langle σv\rangle$, as a function of the dark matter mass. Our results show, in particular, that the thermal dark matter annihilation cross section, $3\times 10^{-26}\, {\rm cm^3 s^{-1}} $, into right-handed neutrinos is excluded for dark matter masses smaller than $200$ GeV. △ Less

Submitted 13 October, 2017; v1 submitted 20 February, 2017; originally announced February 2017.

Comments: 19 pages, 7 figures. It matches the version published in JCAP

Cold dark matter plus not-so-clumpy dark relics

Authors: Roberta Diamanti, Shin'ichiro Ando, Stefano Gariazzo, Olga Mena, Christoph Weniger

Abstract: Various particle physics models suggest that, besides the (nearly) cold dark matter that accounts for current observations, additional but sub-dominant dark relics might exist. These could be warm, hot, or even contribute as dark radiation. We present here a comprehensive study of two-component dark matter scenarios, where the first component is assumed to be cold, and the second is a non-cold the… ▽ More Various particle physics models suggest that, besides the (nearly) cold dark matter that accounts for current observations, additional but sub-dominant dark relics might exist. These could be warm, hot, or even contribute as dark radiation. We present here a comprehensive study of two-component dark matter scenarios, where the first component is assumed to be cold, and the second is a non-cold thermal relic. Considering the cases where the non-cold dark matter species could be either a fermion or a boson, we derive consistent upper limits on the non-cold dark relic energy density for a very large range of velocity dispersions, covering the entire range from dark radiation to cold dark matter. To this end, we employ the latest Planck Cosmic Microwave Background data, the recent BOSS DR11 and other Baryon Acoustic Oscillation measurements, and also constraints on the number of Milky Way satellites, the latter of which provides a measure of the suppression of the matter power spectrum at the smallest scales due to the free-streaming of the non-cold dark matter component. We present the results on the fraction $f_{\rm ncdm}$ of non-cold dark matter with respect to the total dark matter for different ranges of the non-cold dark matter masses. We find that the 2$σ$ limits for non-cold dark matter particles with masses in the range 1--10 keV are $f_{\rm ncdm}\leq0.29$ (0.23) for fermions (bosons), and for masses in the 10--100 keV range they are $f_{\rm ncdm}\leq0.43$ (0.45), respectively. △ Less

Submitted 6 June, 2017; v1 submitted 11 January, 2017; originally announced January 2017.

Comments: 16 pages, 6 figures; matches the version accepted for publication

Novel Spectral Features in MeV Gamma Rays from Dark Matter

Authors: Torsten Bringmann, Ahmad Galea, Andrzej Hryczuk, Christoph Weniger

Abstract: Astrophysical searches for gamma rays are one of the main strategies to probe the annihilation or decay of dark matter particles. We present a new class of distinct sub-GeV spectral features that generically appear in kinematical situations where the available center-of-mass energy in such processes is just above threshold to produce excited meson states. Using a Fisher forecast with realistic ast… ▽ More Astrophysical searches for gamma rays are one of the main strategies to probe the annihilation or decay of dark matter particles. We present a new class of distinct sub-GeV spectral features that generically appear in kinematical situations where the available center-of-mass energy in such processes is just above threshold to produce excited meson states. Using a Fisher forecast with realistic astrophysical backgrounds, we demonstrate that for upcoming experiments like e-ASTROGAM and ComPair these signals can turn out to be the smoking gun in the search for particle dark matter. △ Less

Submitted 16 February, 2017; v1 submitted 14 October, 2016; originally announced October 2016.

Comments: 7 pages revtex4-1, 5 figures; minor updates in sensitivity figures and discussion. Matches published version

Journal ref: Phys. Rev. D 95, 043002 (2017)

Effective Field Theory of Dark Matter: a Global Analysis

Authors: Sebastian Liem, Gianfranco Bertone, Francesca Calore, Roberto Ruiz de Austri, Tim M. P. Tait, Roberto Trotta, Christoph Weniger

Abstract: We present global fits of an effective field theory description of real, and complex scalar dark matter candidates. We simultaneously take into account all possible dimension 6 operators consisting of dark matter bilinears and gauge invariant combinations of quark and gluon fields. We derive constraints on the free model parameters for both the real (five parameters) and complex (seven) scalar dar… ▽ More We present global fits of an effective field theory description of real, and complex scalar dark matter candidates. We simultaneously take into account all possible dimension 6 operators consisting of dark matter bilinears and gauge invariant combinations of quark and gluon fields. We derive constraints on the free model parameters for both the real (five parameters) and complex (seven) scalar dark matter models obtained by combining Planck data on the cosmic microwave background, direct detection limits from LUX, and indirect detection limits from the Fermi Large Area Telescope. We find that for real scalars indirect dark matter searches disfavour a dark matter particle mass below 100 GeV. For the complex scalar dark matter particle current data have a limited impact due to the presence of operators that lead to p-wave annihilation, and also do not contribute to the spin-independent scattering cross- section. Although current data are not informative enough to strongly constrain the theory parameter space, we demonstrate the power of our formalism to reconstruct the theoretical parameters compatible with an actual dark matter detection, by assuming that the excess of gamma rays observed by the Fermi Large Area Telescope towards the Galactic centre is entirely due to dark matter annihilations. Please note that the excess can very well be due to astrophysical sources such as millisecond pulsars. We find that scalar dark matter interacting via effective field theory operators can in principle explain the Galactic centre excess, but that such interpretation is in strong tension with the non-detection of gamma rays from dwarf galaxies in the real scalar case. In the complex scalar case there is enough freedom to relieve the tension. △ Less

Submitted 18 March, 2016; originally announced March 2016.

Comments: 20 pages, 5 figures

Gamma-ray Limits on Neutrino Lines

Authors: Farinaldo S. Queiroz, Carlos E. Yaguna, Christoph Weniger

Abstract: MMonochromatic neutrinos from dark matter annihilations ($χχ\to ν\barν$) are always produced in association with a gamma-ray spectrum generated by electroweak bremsstrahlung. Consequently, these neutrino lines can be searched for not only with neutrino detectors but also indirectly with gamma-ray telescopes. Here, we derive limits on the dark matter annihilation cross section into neutrinos based… ▽ More MMonochromatic neutrinos from dark matter annihilations ($χχ\to ν\barν$) are always produced in association with a gamma-ray spectrum generated by electroweak bremsstrahlung. Consequently, these neutrino lines can be searched for not only with neutrino detectors but also indirectly with gamma-ray telescopes. Here, we derive limits on the dark matter annihilation cross section into neutrinos based on recent Fermi-LAT and HESS data. We find that, for dark matter masses above 200 GeV, gamma-ray data actually set the most stringent constraints on neutrino lines from dark matter annihilation and, therefore, an upper bound on the dark matter total annihilation cross section. In addition, we point out that gamma-ray telescopes, unlike neutrino detectors, have the potential to distinguish the flavor of the final state neutrino. Our results indicate that we have already entered into a new era where gamma-ray telescopes are more sensitive than neutrino detectors to neutrino lines from dark matter annihilation. △ Less

Submitted 5 July, 2016; v1 submitted 18 February, 2016; originally announced February 2016.

Comments: 5 pages, 3 figures. It matches published version

Journal ref: JCAP 1605 (2016) no.05, 050

Radio detection prospects for a bulge population of millisecond pulsars as suggested by Fermi LAT observations of the inner Galaxy

Authors: Francesca Calore, Mattia Di Mauro, Fiorenza Donato, Jason W. T. Hessels, Christoph Weniger

Abstract: Analogously to globular clusters, the dense stellar environment of the Galactic center has been proposed to host a large population of as-yet undetected millisecond pulsars (MSPs). Recently, this hypothesis found support in the analysis of gamma rays from the inner Galaxy seen by the Large Area Telescope (LAT) aboard the Fermi satellite, which revealed a possible excess of diffuse GeV photons in t… ▽ More Analogously to globular clusters, the dense stellar environment of the Galactic center has been proposed to host a large population of as-yet undetected millisecond pulsars (MSPs). Recently, this hypothesis found support in the analysis of gamma rays from the inner Galaxy seen by the Large Area Telescope (LAT) aboard the Fermi satellite, which revealed a possible excess of diffuse GeV photons in the inner 15 deg about the Galactic center (Fermi GeV excess). The excess can be interpreted as the collective emission of thousands of MSPs in the Galactic bulge, with a spherical distribution that strongly peaks towards the Galactic center. In order to fully establish the MSP interpretation, it is essential to find corroborating evidence in multi-wavelength searches, most notably through the detection of radio pulsation from individual bulge MSPs. Based on globular cluster observations and the gamma-ray emission from the inner Galaxy, we investigate the prospects for detecting MSPs in the Galactic bulge. While previous pulsar surveys failed to identify this population, we demonstrate that, in the upcoming years, new large-area surveys with focus on regions a few degrees north or south of the Galactic center should lead to the detection of dozens of bulge MSPs. Additionally, we show that, in the near future, deep targeted searches of unassociated Fermi sources should be able to detect the first few MSPs in the bulge. The prospects for these deep searches are enhanced by a tentative gamma-ray/radio correlation that we infer from high-latitude gamma-ray MSPs. Such detections would constitute the first clear discoveries of field MSPs in the Galactic bulge, with far-reaching implications for gamma-ray observations, the formation history of the central Milky Way and strategy optimization for future radio observations. △ Less

Submitted 12 September, 2016; v1 submitted 21 December, 2015; originally announced December 2015.

Comments: 24 pages, 17 figures, 5 tables. Minor clarifications. Matches version published in ApJ

Journal ref: Astrophys.J. 827 (2016) no.2, 143

Combined analysis of effective Higgs portal dark matter models

Authors: Ankit Beniwal, Filip Rajec, Christopher Savage, Pat Scott, Christoph Weniger, Martin White, Anthony G. Williams

Abstract: We combine and extend the analyses of effective scalar, vector, Majorana and Dirac fermion Higgs portal models of dark matter (DM), in which DM couples to the Standard Model (SM) Higgs boson via an operator of the form $\mathcal{O}_{\textrm{DM}}\, H^\dagger H$. For the fermion models, we take an admixture of scalar $\overlineψ ψ$ and pseudoscalar $\overlineψ iγ_5 ψ$ interaction terms. For each mod… ▽ More We combine and extend the analyses of effective scalar, vector, Majorana and Dirac fermion Higgs portal models of dark matter (DM), in which DM couples to the Standard Model (SM) Higgs boson via an operator of the form $\mathcal{O}_{\textrm{DM}}\, H^\dagger H$. For the fermion models, we take an admixture of scalar $\overlineψ ψ$ and pseudoscalar $\overlineψ iγ_5 ψ$ interaction terms. For each model, we apply constraints on the parameter space based on the Planck measured DM relic density and the LHC limits on the Higgs invisible branching ratio. For the first time, we perform a consistent study of the indirect detection prospects for these models based on the WMAP7/Planck observations of the cosmic microwave background, a combined analysis of 15 dwarf spheroidal galaxies by Fermi-LAT and the upcoming Cherenkov Telescope Array (CTA). We also perform a correct treatment of the momentum-dependent direct search cross section that arises from the pseudoscalar interaction term in the fermionic DM theories. We find, in line with previous studies, that current and future direct search experiments such as LUX and XENON1T can exclude much of the parameter space, and we demonstrate that a joint observation in both indirect and direct searches is possible for high mass weakly interacting massive particles. In the case of a pure pseudoscalar interaction of a fermionic DM candidate, future gamma-ray searches are the only class of experiment capable of probing the high mass range of the theory. △ Less

Submitted 14 June, 2016; v1 submitted 20 December, 2015; originally announced December 2015.

Comments: 53 pages, 16 figures, REVTeX-4.1 processed with BibTeX and pdfLaTeX. v2: added references; submitted to Phys. Rev. D. v3: added references and corrected few typos; matches the version accepted for publication in Phys. Rev. D. v4: matches the version published in Phys. Rev. D

Report number: ADP-15-52/T954

Journal ref: Phys. Rev. D 93, 115016 (2016)

Global analysis of the pMSSM in light of the Fermi GeV excess: prospects for the LHC Run-II and astroparticle experiments

Authors: Gianfranco Bertone, Francesca Calore, Sascha Caron, Roberto Ruiz de Austri, Jong Soo Kim, Roberto Trotta, Christoph Weniger

Abstract: We present a new global fit of the 19-dimensional phenomenological Minimal Supersymmetric Standard Model (pMSSM-19) that comply with all the latest experimental results from dark matter indirect, direct and accelerator dark matter searches. We show that the model provides a satisfactory explanation of the excess of gamma-rays from the Galactic centre observed by the Fermi~Large Area Telescope, ass… ▽ More We present a new global fit of the 19-dimensional phenomenological Minimal Supersymmetric Standard Model (pMSSM-19) that comply with all the latest experimental results from dark matter indirect, direct and accelerator dark matter searches. We show that the model provides a satisfactory explanation of the excess of gamma-rays from the Galactic centre observed by the Fermi~Large Area Telescope, assuming that it is produced by the annihilation of neutralinos in the Milky Way halo. We identify two regions that pass all the constraints: the first corresponds to neutralinos with a mass ~80-100 GeV annihilating into WW with a branching ratio of 95% ; the second to heavier neutralinos, with mass ~180-200 GeV annihilating into t tbar with a branching ratio of 87%. We show that neutralinos compatible with the Galactic centre GeV excess will soon be within the reach of LHC run-II -- notably through searches for charginos and neutralinos, squarks and light smuons -- and of Xenon1T, thanks to its unprecedented sensitivity to spin-dependent cross-section off neutrons. △ Less

Submitted 18 April, 2016; v1 submitted 24 July, 2015; originally announced July 2015.

Comments: Minor changes following referee reports. Main conclusions unchanged. Matches version published in JCAP

Journal ref: JCAP04(2016)037

The Galactic Center GeV Excess from a Series of Leptonic Cosmic-Ray Outbursts

Authors: Ilias Cholis, Carmelo Evoli, Francesca Calore, Tim Linden, Christoph Weniger, Dan Hooper

Abstract: It has been proposed that a recent outburst of cosmic-ray electrons could account for the excess of GeV-scale gamma rays observed from the region surrounding the Galactic Center. After studying this possibility in some detail, we identify scenarios in which a series of leptonic cosmic-ray outbursts could plausibly generate the observed excess. The morphology of the emission observed outside of… ▽ More It has been proposed that a recent outburst of cosmic-ray electrons could account for the excess of GeV-scale gamma rays observed from the region surrounding the Galactic Center. After studying this possibility in some detail, we identify scenarios in which a series of leptonic cosmic-ray outbursts could plausibly generate the observed excess. The morphology of the emission observed outside of $\sim1^{\circ}-2^{\circ}$ from the Galactic Center can be accommodated with two outbursts, one which took place approximately $\sim10^6$ years ago, and another (injecting only about 10\% as much energy as the first) about $\sim10^5$ years ago. The emission observed from the innermost $\sim1^{\circ}-2^{\circ}$ requires one or more additional recent outbursts and/or a contribution from a centrally concentrated population of unresolved millisecond pulsars. In order to produce a spectrum that is compatible with the measured excess (whose shape is approximately uniform over the region of the excess), the electrons from the older outburst must be injected with significantly greater average energy than those injected more recently, enabling their spectra to be similar after $\sim10^6$ years of energy losses. △ Less

Submitted 16 June, 2015; originally announced June 2015.

Comments: 28 pages, 7 figures, 3 tables, 1 appendix

Report number: FERMILAB-PUB-15-255-A