Showing 1–50 of 73 results for author: Cholis, I

Robustness of the Galactic Center Excess Morphology Against Masking

Authors: Yi-Ming Zhong, Ilias Cholis

Abstract: The Galactic Center Excess (GCE) remains an enduring mystery, with leading explanations being annihilating dark matter or an unresolved population of millisecond pulsars. Analyzing the morphology of the GCE provides critical clues to identify its exact origin. We investigate the robustness of the inferred GCE morphology against the effects of masking, an important step in the analysis where the ga… ▽ More The Galactic Center Excess (GCE) remains an enduring mystery, with leading explanations being annihilating dark matter or an unresolved population of millisecond pulsars. Analyzing the morphology of the GCE provides critical clues to identify its exact origin. We investigate the robustness of the inferred GCE morphology against the effects of masking, an important step in the analysis where the gamma-ray emission from point sources and the galactic disk are excluded. Using different masks constructed from Fermi point source catalogs and a wavelet method, we find that the GCE morphology, particularly its ellipticity and cuspiness, is relatively independent of the choice of mask for energies above 2-3 GeV. The GCE morphology systematically favors an approximately spherical shape, as expected for dark matter annihilation. Compared to various stellar bulge profiles, a spherical dark matter annihilation profile better fits the data across different masks and galactic diffuse emission backgrounds, except for the stellar bulge profile from Coleman et al. (2020), which provides a similar fit to the data. Modeling the GCE with two components, one from dark matter annihilation and one tracing the Coleman Bulge, we find this two-component model outperforms any single component or combinations of dark matter annihilation and other stellar bulge profiles. Uncertainty remains about the exact fraction contributed by each component across different background models and masks. However, when the Coleman Bulge dominates, its corresponding spectrum lacks characteristics typically associated with millisecond pulsars, suggesting that it mostly models the emission from other sources instead of the GCE that is still present and spherically symmetric. △ Less

Submitted 11 June, 2024; v1 submitted 4 January, 2024; originally announced January 2024.

Comments: v2 as accepted by PRD, masks available at https://github.com/ymzhong/gce_mask

Observing Signals of Spectral Features in the Cosmic-Ray Positrons and Electrons from Milky Way Pulsars

Authors: Ilias Cholis, Thressay Hoover

Abstract: The Alpha Magnetic Spectrometer (AMS-02) has provided unprecedented precision measurements of the electron and positron cosmic-ray fluxes and the positron fraction spectrum. At the higher energies, sources as energetic local pulsars, may contribute to both cosmic-ray species. The discreteness of the source population, can result in features both on the positron fraction measurement and in the resp… ▽ More The Alpha Magnetic Spectrometer (AMS-02) has provided unprecedented precision measurements of the electron and positron cosmic-ray fluxes and the positron fraction spectrum. At the higher energies, sources as energetic local pulsars, may contribute to both cosmic-ray species. The discreteness of the source population, can result in features both on the positron fraction measurement and in the respective electron and positron spectra. For the latter, those would coincide in energy and would contrast predictions of smooth spectra as from particle dark matter. In this work, using a library of pulsar population models for the local part of the Milky Way, we perform a power-spectrum analysis on the cosmic-ray positron fraction. We also develop a technique to cross-correlate the electron and positron fluxes. We show that both such analyses, can be used to search statistically for the presence of spectral wiggles in the cosmic-ray data. For a significant fraction of our pulsar simulations, those techniques are already sensitive enough to give a signal for the presence of those features above the regular noise, with forthcoming observations making them even more sensitive. Finally, by cross-correlating the AMS-02 electron and positron spectra, we find an intriguing first hint for a positive correlation between them, of the kind expected by a population of local pulsars. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: 16 pages, 10 figures, 2 tables

Revisiting GeV-scale annihilating dark matter with the AMS-02 positron fraction

Authors: Iason Krommydas, Ilias Cholis

Abstract: Antimatter cosmic-rays are used to probe new phenomena in physics, including dark matter annihilation. We use the cosmic-ray positron fraction spectrum by the Alpha Magnetic Spectrometer, to search for such an annihilation signal in the Galaxy. We focus on dark matter with mass between 5 and 120 GeV, producing high-energy electrons and positrons. In these cosmic-ray energies the interplay of multi… ▽ More Antimatter cosmic-rays are used to probe new phenomena in physics, including dark matter annihilation. We use the cosmic-ray positron fraction spectrum by the Alpha Magnetic Spectrometer, to search for such an annihilation signal in the Galaxy. We focus on dark matter with mass between 5 and 120 GeV, producing high-energy electrons and positrons. In these cosmic-ray energies the interplay of multiple astrophysical sources and phenomena, makes this search highly sensitive to the underlying astrophysical background assumptions. We use a vast public library of astrophysical models for the cosmic-ray positron fraction background, to derive robust upper limits on the dark matter's annihilation cross section for a number of annihilation channels. This library accounts for different types of cosmic-ray sources and uncertainties on their distribution in space and time. Also, it accounts for uncertainties on those sources' output, their injected into the interstellar medium cosmic-ray spectra and for uncertainties on cosmic-ray propagation. For any given dark matter particle mass and annihilation channel, upper limits on the annihilation cross section are given by bands that stretch a full order of magnitude in its value. Our work provides weaker limits compared to earlier results, that are however robust to all the relevant astrophysical uncertainties. Between 5 and 15 GeV, we find indications for a possible excess flux of cosmic-ray electrons and positrons. That excess is found for most, but not all of our astrophysical background parameter space, and its significance can vary appreciably. Further scrutiny is necessary to improve the understanding of these lower energy cosmic rays. Finally, we note that even if an excess signal is found in these energies, the current background uncertainties do not allow us to accurately deduce its underlying particle properties. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: 21 pages, 8 figures, 1 table, zenodo link: https://zenodo.org/record/7178634

A Phantom Menace: On the Morphology of the Galactic Center Excess

Authors: Samuel D. McDermott, Yi-Ming Zhong, Ilias Cholis

Abstract: The characteristics of the Galactic Center Excess (GCE) emission observed in gamma-ray energies -- especially the morphology of the GCE -- remain a hotly debated subject. The manner in which the dominant diffuse gamma-ray background is modeled has been claimed to have a determining effect on the preferred morphology. In this work, we compare two distinct approaches to the galactic diffuse gamma-ra… ▽ More The characteristics of the Galactic Center Excess (GCE) emission observed in gamma-ray energies -- especially the morphology of the GCE -- remain a hotly debated subject. The manner in which the dominant diffuse gamma-ray background is modeled has been claimed to have a determining effect on the preferred morphology. In this work, we compare two distinct approaches to the galactic diffuse gamma-ray emission background: the first approach models this emission through templates calculated from a sequence of well-defined astrophysical assumptions, while the second approach divides surrogates for the background gamma-ray emission into cylindrical galactocentric rings with free independent normalizations. At the latitudes that we focus on, we find that the former approach works better, and that the overall best fit is obtained for an astrophysically motivated fit when the GCE follows the morphology expected of dark matter annihilation. Quantitatively, the improvement compared to the best ring-based fits is roughly 6500 in the chi^2 and roughly 4000 in the log of the Bayesian evidence. △ Less

Submitted 6 October, 2023; v1 submitted 31 August, 2022; originally announced September 2022.

Comments: 10 pages (including the appendix), code available at https://github.com/samueldmcdermott/gcepy

Report number: FERMILAB-PUB-22-650-T

Journal ref: MNRAS, Vol 522, 1, L21-25 (2023)

Possible counterpart signal of the Fermi bubbles at the cosmic-ray positrons

Authors: Ilias Cholis, Iason Krommydas

Abstract: The inner galaxy has hosted cosmic-ray burst events including those responsible for the gamma-ray Fermi bubbles and the eROSITA bubbles in X-rays. In this work, we study the AMS-02 positron fraction and find three features around 12, 21 and 48 GeV of which the lowest energy has a 1.4 to 4.9-$σ$ significance, depending on astrophysical background assumptions. Using background simulations that expla… ▽ More The inner galaxy has hosted cosmic-ray burst events including those responsible for the gamma-ray Fermi bubbles and the eROSITA bubbles in X-rays. In this work, we study the AMS-02 positron fraction and find three features around 12, 21 and 48 GeV of which the lowest energy has a 1.4 to 4.9-$σ$ significance, depending on astrophysical background assumptions. Using background simulations that explain the cosmic-ray positron fraction, positron flux and electron plus positron flux, by primary, secondary cosmic rays and cosmic rays from local pulsars, we test these spectral features as originating from electron/positron burst events from the inner galaxy. We find the 12 GeV feature, to be explained by an event of age $τ\simeq 3 - 10$ Myr; in agreement with the proposed age of the Fermi bubbles. Furthermore, the energy in cosmic-ray electrons and positrons propagating along the galactic disk and not within the Fermi bubbles volume, is estimated to be $10^{51.5}-10^{57.5}$ ergs, or $O(10^{-4}) -O(1)$ the cosmic-ray energy causing the Fermi bubbles. We advocate that these positron fraction features, are the counterpart signals of the Fermi bubbles, or of substructures in them, or of the eROSITA bubbles. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: 10 pages, 3 figures and 1 table, Zenodo link: https://zenodo.org/record/6995268#.YvsoJS ZNp8

Constraining the Charge-, Time- and Rigidity-Dependence of Cosmic-Ray Solar Modulation with AMS-02 Observations during Solar Cycle 24

Authors: Ilias Cholis, Ian McKinnon

Abstract: Our basic theoretical understanding of the sources of cosmic rays and their propagation through the interstellar medium is hindered by the Sun, that through the solar wind affects the observed cosmic-ray spectra. This effect is known as solar modulation. Recently released cosmic-ray observations from the Alpha Magnetic Spectrometer (AMS-02) and publicly available measurements of the solar wind pro… ▽ More Our basic theoretical understanding of the sources of cosmic rays and their propagation through the interstellar medium is hindered by the Sun, that through the solar wind affects the observed cosmic-ray spectra. This effect is known as solar modulation. Recently released cosmic-ray observations from the Alpha Magnetic Spectrometer (AMS-02) and publicly available measurements of the solar wind properties from the Advanced Composition Explorer and the Wilcox observatory allow us to test the analytical modeling of the time-, charge- and rigidity-dependence of solar modulation. We rely on associating measurements on the local heliospheric magnetic field and the heliospheric current sheet's tilt angle, to model the time-dependence and amplitude of cosmic-ray solar modulation. We find evidence for the solar modulation's charge- and rigidity-dependence during the era of solar cycle 24. Our analytic prescription to model solar modulation can explain well the large-scale time-evolution of positively charged cosmic-ray fluxes in the range of rigidities from 1 to 10 GV. We also find that cosmic-ray electron fluxes measured during the first years of cycle 24 are less trivial to explain, due to the complex and rapidly evolving structure of the Heliosphere's magnetic field that they experienced as they propagated inwards. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 13 pages, 7 figures, 2 tables

Snowmass2021 Cosmic Frontier White Paper: Puzzling Excesses in Dark Matter Searches and How to Resolve Them

Authors: Rebecca K. Leane, Seodong Shin, Liang Yang, Govinda Adhikari, Haider Alhazmi, Tsuguo Aramaki, Daniel Baxter, Francesca Calore, Regina Caputo, Ilias Cholis, Tansu Daylan, Mattia Di Mauro, Philip von Doetinchem, Ke Han, Dan Hooper, Shunsaku Horiuchi, Doojin Kim, Kyoungchul Kong, Rafael F. Lang, Qing Lin, Tim Linden, Jianglai Liu, Oscar Macias, Siddharth Mishra-Sharma, Alexander Murphy , et al. (14 additional authors not shown)

Abstract: Intriguing signals with excesses over expected backgrounds have been observed in many astrophysical and terrestrial settings, which could potentially have a dark matter origin. Astrophysical excesses include the Galactic Center GeV gamma-ray excess detected by the Fermi Gamma-Ray Space Telescope, the AMS antiproton and positron excesses, and the 511 and 3.5 keV X-ray lines. Direct detection excess… ▽ More Intriguing signals with excesses over expected backgrounds have been observed in many astrophysical and terrestrial settings, which could potentially have a dark matter origin. Astrophysical excesses include the Galactic Center GeV gamma-ray excess detected by the Fermi Gamma-Ray Space Telescope, the AMS antiproton and positron excesses, and the 511 and 3.5 keV X-ray lines. Direct detection excesses include the DAMA/LIBRA annual modulation signal, the XENON1T excess, and low-threshold excesses in solid state detectors. We discuss avenues to resolve these excesses, with actions the field can take over the next several years. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 57 pages, solicited white paper submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)

Snowmass2021 Cosmic Frontier: Synergies between dark matter searches and multiwavelength/multimessenger astrophysics

Authors: Shin'ichiro Ando, Sebastian Baum, Michael Boylan-Kolchin, Esra Bulbul, Michael Burgess, Ilias Cholis, Philip von Doetinchem, JiJi Fan, Patrick J. Harding, Shunsaku Horiuchi, Rebecca K. Leane, Oscar Macias, Katie Mack, Kohta Murase, Lina Necib, Ibles Olcina, Laura Olivera-Nieto, Jong-Chul Park, Kerstin Perez, Marco Regis, Nicholas L. Rodd, Carsten Rott, Kuver Sinha, Volodymyr Takhistov, Yun-Tse Tsai , et al. (1 additional authors not shown)

Abstract: This whitepaper focuses on the astrophysical systematics which are encountered in dark matter searches. Oftentimes in indirect and also in direct dark matter searches, astrophysical systematics are a major limiting factor to sensitivity to dark matter. Just as there are many forms of dark matter searches, there are many forms of backgrounds. We attempt to cover the major systematics arising in dar… ▽ More This whitepaper focuses on the astrophysical systematics which are encountered in dark matter searches. Oftentimes in indirect and also in direct dark matter searches, astrophysical systematics are a major limiting factor to sensitivity to dark matter. Just as there are many forms of dark matter searches, there are many forms of backgrounds. We attempt to cover the major systematics arising in dark matter searches using photons -- radio and gamma rays -- to cosmic rays, neutrinos and gravitational waves. Examples include astrophysical sources of cosmic messengers and their interactions which can mimic dark matter signatures. In turn, these depend on commensurate studies in understanding the cosmic environment -- gas distributions, magnetic field configurations -- as well as relevant nuclear astrophysics. We also cover the astrophysics governing celestial bodies and galaxies used to probe dark matter, from black holes to dwarf galaxies. Finally, we cover astrophysical backgrounds related to probing the dark matter distribution and kinematics, which impact a wide range of dark matter studies. In the future, the rise of multi-messenger astronomy, and novel analysis methods to exploit it for dark matter, will offer various strategic ways to continue to enhance our understanding of astrophysical backgrounds to deliver improved sensitivity to dark matter. △ Less

Submitted 13 March, 2022; originally announced March 2022.

Comments: contribution to Snowmass 2021

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

The Return of the Templates: Revisiting the Galactic Center Excess with Multi-Messenger Observations

Authors: Ilias Cholis, Yi-Ming Zhong, Samuel D. McDermott, Joseph P. Surdutovich

Abstract: The Galactic center excess (GCE) remains one of the most intriguing discoveries from the Fermi Large Area Telescope (LAT) observations. We revisit the characteristics of the GCE by first producing a new set of high-resolution galactic diffuse gamma-ray emission templates, which are ultimately due to cosmic-ray interactions with the interstellar medium. Using multi-messenger observations we constra… ▽ More The Galactic center excess (GCE) remains one of the most intriguing discoveries from the Fermi Large Area Telescope (LAT) observations. We revisit the characteristics of the GCE by first producing a new set of high-resolution galactic diffuse gamma-ray emission templates, which are ultimately due to cosmic-ray interactions with the interstellar medium. Using multi-messenger observations we constrain the properties of the galactic diffuse emission. The broad properties of the GCE that we find in this work are qualitatively unchanged despite the introduction of this new set of templates, though its quantitative features appear mildly different than those obtained in previous analyses. In particular, we find a high-energy tail at higher significance than previously reported. This tail is very prominent in the northern hemisphere, and less so in the southern hemisphere. This strongly affects one prominent interpretation of the excess: known millisecond pulsars are incapable of producing this high-energy emission, even in the relatively softer southern hemisphere, and are therefore disfavored as the sole explanation of the GCE. The annihilation of dark matter particles of mass $40^{+10}_{-7}$ GeV (95$\%$ CL) to $b$ quarks with a cross-section of $σv = 1.4^{+0.6}_{-0.3} \times 10^{-26}$ cm$^{3}$s$^{-1}$ provides a good fit to the excess especially in the relatively cleaner southern sky. Dark matter of the same mass range annihilating to $b$ quarks or heavier dark matter particles annihilating to heavier Standard Model bosons can combine with millisecond pulsars to provide a good fit to the southern hemisphere emission. As part of this paper, we make publicly available all of our templates and the data covariance matrix we have generated to account for systematic uncertainties.[abridged] △ Less

Submitted 20 May, 2022; v1 submitted 17 December, 2021; originally announced December 2021.

Comments: v1: 41 pages, 23 figures and 9 tables. v2: 42 pages 24 figures and 9 tables. Reflects PRD final version with an additional figure. Diffuse flux maps available at https://zenodo.org/record/6423495#.YogWjS-ZPGI

Report number: FERMILAB-PUB-21-709-T

Journal ref: Phys. Rev. D 105, 103023 (2022)

Utilizing cosmic-ray positron and electron observations to probe the averaged properties of Milky Way pulsars

Authors: Ilias Cholis, Iason Krommydas

Abstract: Pulsars have long been studied in the electromagnetic spectrum. Their environments are rich in high-energy cosmic-ray electrons and positrons likely enriching the interstellar medium with such particles. In this work we use recent cosmic-ray observations from the AMS-02, CALET and DAMPE collaborations to study the averaged properties of the local Milky Way pulsar population. We perform simulations… ▽ More Pulsars have long been studied in the electromagnetic spectrum. Their environments are rich in high-energy cosmic-ray electrons and positrons likely enriching the interstellar medium with such particles. In this work we use recent cosmic-ray observations from the AMS-02, CALET and DAMPE collaborations to study the averaged properties of the local Milky Way pulsar population. We perform simulations of the local Milky Way pulsar population, for interstellar medium assumptions in agreement with a range of cosmic-ray nuclei measurements. Each such simulation contains $\sim 10^{4}$ pulsars of unique age, location, initial spin-down power and cosmic-ray electron/positron spectra. We produce more than $7\times 10^{3}$ such Milky Way pulsar simulations. We account for and study i) the pulsars' birth rates and the stochastic nature of their birth, ii) their initial spin-down power distribution, iii) their time evolution in terms of their braking index and characteristic spin-down timescale, iv) the fraction of spin-down power going to cosmic-ray electrons and positrons and v) their propagation through the interstellar medium and the Heliosphere. We find that pulsars of ages $\sim 10^{5}-10^{7}$ yr, have a braking index that on average has to be 3 or larger. Given that electromagnetic spectrum observations of young pulsars find braking indices lower than 3, our work provides strong hints that pulsars' braking index increases on average as they age, allowing them to retain some of their rotational energy. Moreover, we find that pulsars have relatively uniform properties as sources of cosmic-ray electrons and positrons in terms of the spectra they produce and likely release O($10\%$) of their rotational energy to cosmic-rays in the ISM. Finally, we find at $\simeq$12 GeV positrons a spectral feature that suggests a new subpopulation of positron sources contributing at these energies. △ Less

Submitted 10 November, 2021; originally announced November 2021.

Comments: 20 pages, 13 figures and 6 tables. Simulations are available at https://zenodo.org/record/5659004#.YYqnbi-ZN0s

Can Thorne-$\dot{\textrm{Z}}$ytkow Objects source GW190814-type events?

Authors: Ilias Cholis, Konstantinos Kritos, David Garfinkle

Abstract: The LIGO-Virgo collaboration reported in their third run the coalescence event GW190814 involving a 2.6 $M_{\odot}$ object with a 23 $M_{\odot}$ black hole. In this letter we study the conditions under which Thorne-$\dot{\textrm{Z}}$ytkow objects (T$\dot{\textrm{Z}}$Os) can be connected to that type of events. We evaluate first the rate of appearance of T$\dot{\textrm{Z}}$Os in the local Universe.… ▽ More The LIGO-Virgo collaboration reported in their third run the coalescence event GW190814 involving a 2.6 $M_{\odot}$ object with a 23 $M_{\odot}$ black hole. In this letter we study the conditions under which Thorne-$\dot{\textrm{Z}}$ytkow objects (T$\dot{\textrm{Z}}$Os) can be connected to that type of events. We evaluate first the rate of appearance of T$\dot{\textrm{Z}}$Os in the local Universe. Under the assumption that T$\dot{\textrm{Z}}$Os eventually become low mass gap black holes we evaluate how those black holes end up in binaries with other stellar mass black holes and compare to the reported rate for GW190814-type of events (1-23 $\textrm{Gpc}^{-3} \textrm{yr}^{-1}$). We find that T$\dot{\textrm{Z}}$Os in dense stellar clusters can not explain the LIGO-Virgo rate without a T$\dot{\textrm{Z}}$Os population in the field providing a dominant contribution. We also find that T$\dot{\textrm{Z}}$Os formed within hierarchical triple systems in the field with the third more distant star being the progenitor of a stellar mass black hole, may be able to give a rate comparable to that of GW190814-type events. In that case, future observations should discover mergers between stellar mass and low mass gap black holes, with the lower mass spanning the entire low mass gap range. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Comments: 5 pages and 2 figures

Black Holes merging with Low Mass Gap Objects inside Globular Clusters

Authors: Konstantinos Kritos, Ilias Cholis

Abstract: Recently, the LIGO-Virgo collaborations have reported the coalescence of a binary involving a black hole and a low-mass gap object (LMGO) with mass in the range $\sim2.5-5M_\odot$. Such detections, challenge our understanding of the black hole and neutron star mass spectrum, as well as how such binaries evolve especially if isolated. In this work we study the dynamical formation of compact object… ▽ More Recently, the LIGO-Virgo collaborations have reported the coalescence of a binary involving a black hole and a low-mass gap object (LMGO) with mass in the range $\sim2.5-5M_\odot$. Such detections, challenge our understanding of the black hole and neutron star mass spectrum, as well as how such binaries evolve especially if isolated. In this work we study the dynamical formation of compact object pairs, via multiple binary-single exchanges that occur at the cores of globular clusters. We start with a population of binary star systems, which interact with single compact objects as first generation black holes and LMGOs. We evaluate the rate of exchange interactions leading to the formation of compact object binaries. Our calculations include all possible types of binary-single exchange interactions and also the interactions of individual stars with compact object binaries that can evolve their orbital properties, leading to their eventual merger. We perform our calculations for the full range of the observed Milky Way globular cluster environments. We find that the exchanges are efficient in forming hard compact object binaries at the cores of dense astrophysical stellar environments. Furthermore, if the population size of the LMGOs is related to that of neutron stars, the inferred merger rate density of black hole-LMGO binaries inside globular clusters in the local Universe is estimated to be about $0.1 \, \text{Gpc}^{-3}\text{yr}^{-1}$. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 18 pages, 10 figures and 3 appendices

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

Evaluating the Merger Rate of Binary Black Holes from Direct Captures and Third-Body Soft Interactions Using the Milky Way Globular Clusters

Authors: Konstantinos Kritos, Ilias Cholis

Abstract: The multitude of binary black hole coalescence detections in gravitational waves has renewed our interest on environments that can be the cradle of these mergers. In this work we study merger rates of binary black holes in globular clusters that are among the most dense stellar environments and a natural place for the creation of black hole binaries. To model these systems with all their variation… ▽ More The multitude of binary black hole coalescence detections in gravitational waves has renewed our interest on environments that can be the cradle of these mergers. In this work we study merger rates of binary black holes in globular clusters that are among the most dense stellar environments and a natural place for the creation of black hole binaries. To model these systems with all their variations we rely on the observational properties of the known Milky Way globular clusters. We consider direct capture events between black holes, as well as soft interactions of black hole binaries with stars as third bodies that accelerate the evolution of these binaries. We find that binary black holes from direct captures merge at an averaged rate of $0.3-5 \times 10^{-11}$ yr$^{-1}$ per cluster. Third body soft interactions are a much more prominent channel giving an averaged rate of $2-4 \times 10^{-10}$ yr$^{-1}$ per cluster. Those rates in globular clusters can lead to a cumulative merger rate of about 100 mergers per year up to redshift of 1, i.e. a significant fraction of the detectable in the near future binary black hole coalescence events. Further observations of cluster properties both in terms of their masses, profile properties, velocity dispersion of stars and their cosmological distribution, will allow us to better constrain the contribution of these environments to the detectable coalescence events rate. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: 16 pages, 6 figures and 2 appendices

Journal ref: Phys. Rev. D 102, 083016 (2020)

Constraining the Charge-Sign and Rigidity-Dependence of Solar Modulation

Authors: Ilias Cholis, Dan Hooper, Tim Linden

Abstract: Our ability to identify the sources of cosmic rays and understand how these particles propagate through the interstellar medium is hindered by the combined effects of the solar wind and its embedded magnetic field, collectively known as solar modulation. In this paper, we build upon our previous work to model and constrain the effects of solar modulation on the cosmic-ray spectrum, using data from… ▽ More Our ability to identify the sources of cosmic rays and understand how these particles propagate through the interstellar medium is hindered by the combined effects of the solar wind and its embedded magnetic field, collectively known as solar modulation. In this paper, we build upon our previous work to model and constrain the effects of solar modulation on the cosmic-ray spectrum, using data from AMS-02 and BESS Polar II collected between 2007 and 2012, during which the heliospheric magnetic field was in a state of negative polarity. Our model uses measurements of the heliospheric magnetic field and the tilt angle of the heliospheric current sheet to accurately predict the effects of solar modulation as a function of time, charge, and rigidity. By incorporating data from a period of negative polarity, we have been able to robustly observe and constrain the charge-dependent effects of solar modulation. △ Less

Submitted 1 July, 2020; originally announced July 2020.

Comments: 9 pages, 8 figures, 2 appendices

Report number: FERMILAB-PUB-20-271-AE-T

Anti-Deuterons and Anti-Helium Nuclei from Annihilating Dark Matter

Authors: Ilias Cholis, Tim Linden, Dan Hooper

Abstract: Recent studies of the cosmic-ray antiproton-to-proton ratio have identified an excess of $\sim$10-20 GeV antiprotons relative to the predictions of standard astrophysical models. Intriguingly, the properties of this excess are consistent with the same range of dark matter models that can account for the long-standing excess of $γ$-rays observed from the Galactic Center. Such dark matter candidates… ▽ More Recent studies of the cosmic-ray antiproton-to-proton ratio have identified an excess of $\sim$10-20 GeV antiprotons relative to the predictions of standard astrophysical models. Intriguingly, the properties of this excess are consistent with the same range of dark matter models that can account for the long-standing excess of $γ$-rays observed from the Galactic Center. Such dark matter candidates can also produce significant fluxes of anti-deuterium and anti-helium nuclei. Here we study the production and transport of such particles, both from astrophysical processes as well as from dark matter annihilation. Importantly, in the case of AMS-02, we find that Alfvénic reacceleration (i.e., diffusion in momentum space) can boost the expected number of $\bar{\rm d}$ and ${}^{3}\overline{\textrm{He}}$ events from annihilating dark matter by an order of magnitude or more. For relatively large values of the Alfvén speed, and for dark matter candidates that are capable of producing the antiproton and $γ$-ray excesses, we expect annihilations to produce a few anti-deuteron events and about one anti-helium event in six years of AMS-02 data. This is particularly interesting in light of recent reports from the AMS-02 Collaboration describing the detection of a number of anti-helium candidate events. △ Less

Submitted 26 August, 2020; v1 submitted 23 January, 2020; originally announced January 2020.

Comments: 9 pages, 7 figures, 3 appendices, in v2 minor changes to reflect journal version

Report number: FERMILAB-PUB-20-021-A

Journal ref: Phys. Rev. D 102, 103019 (2020)

A New Mask for An Old Suspect: Testing the Sensitivity of the Galactic Center Excess to the Point Source Mask

Authors: Yi-Ming Zhong, Samuel D. McDermott, Ilias Cholis, Patrick J. Fox

Abstract: The Fermi-LAT collaboration has recently released a new point source catalog, referred to as 4FGL. For the first time, we perform a template fit using information from this new catalog and find that the Galactic center excess is still present. On the other hand, we find that a wavelet-based search for point sources is highly sensitive to the use of the 4FGL catalog: no excess of bright regions on… ▽ More The Fermi-LAT collaboration has recently released a new point source catalog, referred to as 4FGL. For the first time, we perform a template fit using information from this new catalog and find that the Galactic center excess is still present. On the other hand, we find that a wavelet-based search for point sources is highly sensitive to the use of the 4FGL catalog: no excess of bright regions on small angular scales is apparent when we mask out 4FGL point sources. We postulate that the 4FGL catalog contains the large majority of bright point sources that have previously been suggested to account for the excess in gamma rays detected at the Galactic center in Fermi-LAT data. Furthermore, after identifying which bright sources have no known counterpart, we place constraints on the luminosity function necessary for point sources to explain the smooth emission seen in the template fit. △ Less

Submitted 27 November, 2019; originally announced November 2019.

Comments: 9 pages, 8 figures

Report number: Fermilab-Pub-19-575-T

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

Electromagnetic probes of primordial black holes as dark matter

Authors: Y. Ali-Haimoud, S. Clesse, J. Garcia-Bellido, A. Kashlinsky, L. Wyrzykowski, A. Achucarro, L. Amendola, J. Annis, A. Arbey, R. G. Arendt, F. Atrio-Barandela, N. Bellomo, K. Belotsky, J-L. Bernal, S. Bird, V. Bozza, C. Byrnes, S. Calchi Novati, F. Calore, B. J. Carr, J. Chluba, I. Cholis, A. Cieplak, P. Cole, I. Dalianis , et al. (69 additional authors not shown)

Abstract: The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence wit… ▽ More The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing. △ Less

Submitted 12 March, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: Science whitepaper submitted to the Astro2020 Decadal Survey

Tests of General Relativity and Fundamental Physics with Space-based Gravitational Wave Detectors

Authors: Emanuele Berti, Enrico Barausse, Ilias Cholis, Juan Garcia-Bellido, Kelly Holley-Bockelmann, Scott A. Hughes, Bernard Kelly, Ely D. Kovetz, Tyson B. Littenberg, Jeffrey Livas, Guido Mueller, Priya Natarajan, David H. Shoemaker, Deirdre Shoemaker, Jeremy D. Schnittman, Michele Vallisneri, Nicolas Yunes

Abstract: Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as… ▽ More Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe. △ Less

Submitted 7 March, 2019; originally announced March 2019.

Comments: White Paper submitted on March 7, 2019 to Astro2020 (2020 Decadal Survey on Astronomy and Astrophysics)

A Robust Excess in the Cosmic-Ray Antiproton Spectrum: Implications for Annihilating Dark Matter

Authors: Ilias Cholis, Tim Linden, Dan Hooper

Abstract: An excess of $\sim$10-20 GeV cosmic-ray antiprotons has been identified in the spectrum reported by the AMS-02 Collaboration. The systematic uncertainties associated with this signal, however, have made it difficult to interpret these results. In this paper, we revisit the uncertainties associated with the time, charge and energy-dependent effects of solar modulation, the antiproton production cro… ▽ More An excess of $\sim$10-20 GeV cosmic-ray antiprotons has been identified in the spectrum reported by the AMS-02 Collaboration. The systematic uncertainties associated with this signal, however, have made it difficult to interpret these results. In this paper, we revisit the uncertainties associated with the time, charge and energy-dependent effects of solar modulation, the antiproton production cross section, and interstellar cosmic-ray propagation. After accounting for these uncertainties, we confirm the presence of a 4.7$σ$ antiproton excess, consistent with that arising from a $m_χ \approx 64-88$ GeV dark matter particle annihilating to $b\bar{b}$ with a cross section of $σv \simeq (0.8-5.2) \times 10^{-26}$ cm$^{3}$/s. If we allow for the stochastic acceleration of secondary antiprotons in supernova remnants, the data continues to favor a similar range of dark matter models ($m_χ\approx 46-94$ GeV, $σv \approx (0.7-3.8)\times 10^{-26}$ cm$^3/$s) with a significance of 3.3$σ$. The same range of dark matter models that are favored to explain the antiproton excess can also accommodate the excess of GeV-scale gamma rays observed from the Galactic Center. △ Less

Submitted 16 August, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: 13 pages, 9 figures and one appendix, changes in v2 in agreement with published version, small revisions in text, results and conclusions unchanged

Report number: FERMILAB-PUB-19-091-A

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

Bounds on Ultra-Light Hidden-Photon Dark Matter from 21cm at Cosmic Dawn

Authors: Ely D. Kovetz, Ilias Cholis, David E. Kaplan

Abstract: Ultra-light hidden-photon dark matter produces an oscillating electric field in the early Universe plasma, which in turn induces an electric current in its ionized component whose dissipation results in heat transfer from the dark matter to the plasma. This will affect the global 21cm signal from the Dark Ages and Cosmic Dawn. In this work we focus on the latter, in light of the reported detection… ▽ More Ultra-light hidden-photon dark matter produces an oscillating electric field in the early Universe plasma, which in turn induces an electric current in its ionized component whose dissipation results in heat transfer from the dark matter to the plasma. This will affect the global 21cm signal from the Dark Ages and Cosmic Dawn. In this work we focus on the latter, in light of the reported detection by the EDGES collaboration of an absorption signal at frequencies corresponding to redshift z~17. By measuring the 21cm global signal, a limit can be placed on the amount of gas heating, and thus the kinetic mixing strength $\varepsilon$ between the hidden and ordinary photons can be constrained. Our inferred 21cm bounds on $\varepsilon$ in the mass range $10^{-23}\,{\rm eV}\lesssim m_χ\lesssim10^{-13}\,{\rm eV}$ are the strongest to date. △ Less

Submitted 4 September, 2018; originally announced September 2018.

Comments: 5 pages, 3 figures

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

Where do the AMS-02 anti-helium events come from?

Authors: Vivian Poulin, Pierre Salati, Ilias Cholis, Marc Kamionkowski, Joseph Silk

Abstract: We discuss the origin of the anti-helium-3 and -4 events possibly detected by AMS-02. Using up-to-date semi-analytical tools, we show that spallation from primary hydrogen and helium nuclei onto the ISM predicts a $\overline{{}^3{\rm He}}$ flux typically one to two orders of magnitude below the sensitivity of AMS-02 after 5 years, and a $\overline{{}^4{\rm He}}$ flux roughly 5 orders of magnitude… ▽ More We discuss the origin of the anti-helium-3 and -4 events possibly detected by AMS-02. Using up-to-date semi-analytical tools, we show that spallation from primary hydrogen and helium nuclei onto the ISM predicts a $\overline{{}^3{\rm He}}$ flux typically one to two orders of magnitude below the sensitivity of AMS-02 after 5 years, and a $\overline{{}^4{\rm He}}$ flux roughly 5 orders of magnitude below the AMS-02 sensitivity. We argue that dark matter annihilations face similar difficulties in explaining this event. We then entertain the possibility that these events originate from anti-matter-dominated regions in the form of anti-clouds or anti-stars. In the case of anti-clouds, we show how the isotopic ratio of anti-helium nuclei might suggest that BBN has happened in an inhomogeneous manner, resulting in anti-regions with a anti-baryon-to-photon ratio $\barη\simeq10^{-3}η$. We discuss properties of these regions, as well as relevant constraints on the presence of anti-clouds in our Galaxy. We present constraints from the survival of anti-clouds in the Milky-Way and in the early Universe, as well as from CMB, gamma-ray and cosmic-ray observations. In particular, these require the anti-clouds to be almost free of normal matter. We also discuss an alternative where anti-domains are dominated by surviving anti-stars. We suggest that part of the unindentified sources in the 3FGL catalog can originate from anti-clouds or anti-stars. AMS-02 and GAPS data could further probe this scenario. △ Less

Submitted 23 March, 2019; v1 submitted 27 August, 2018; originally announced August 2018.

Comments: 16 pages, 4 figures. Comments welcome! v2: Discussion about constraints on anti-matter in our Galaxy and in the early Universe greatly improved. A number of typos corrected. Conclusions unaffected. Some Refs added. v3: some minor text corrections/improvements. updated to match published version

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

Studying the Milky Way Pulsar Population with Cosmic-Ray Leptons

Authors: Ilias Cholis, Tanvi Karwal, Marc Kamionkowski

Abstract: Recent measurements of cosmic-ray electron and positron spectra at energies from a GeV to 5 TeV, as well as radio, X-ray and a wide range of gamma-ray observations of pulsar-wind nebulae, indicate that pulsars are significant sources of high-energy cosmic-ray electrons and positrons. Here we calculate the local cosmic-ray $e^\pm$ energy spectra from pulsars taking into account models for (a) the d… ▽ More Recent measurements of cosmic-ray electron and positron spectra at energies from a GeV to 5 TeV, as well as radio, X-ray and a wide range of gamma-ray observations of pulsar-wind nebulae, indicate that pulsars are significant sources of high-energy cosmic-ray electrons and positrons. Here we calculate the local cosmic-ray $e^\pm$ energy spectra from pulsars taking into account models for (a) the distribution of the pulsars spin-down properties; (b) the cosmic-ray source spectra; and (c) the physics of cosmic-ray propagation. We then use the measured cosmic-ray fluxes from AMS-02, CALET and DAMPE to constrain the space of pulsar and cosmic-ray-propagation models and in particular, local cosmic-ray diffusion and energy losses, the pulsars' energy-loss time-dependence, and the injected $e^{\pm}$ spectra. We find that the lower estimates for the local $e^{\pm}$ energy losses are inconsistent with the data. We also find that pulsar braking indexes of 2.5 or less for sources with ages more than 10 kyr are strongly disfavored. Moreover, the cosmic-ray data are consistent with a wide range of assumptions on the $e^{\pm}$ injection spectral properties and on the distribution of initial spin-down powers. Above a TeV in energy, we find that pulsars can easily explain the observed change in the $e^{+} + e^{-}$ spectral slope. These conclusions are valid as long as pulsars contribute $\gtrsim10\%$ of the observed cosmic-ray $e^\pm$ at energies $\gtrsim100$ GeV. △ Less

Submitted 13 July, 2018; originally announced July 2018.

Comments: 14 pages, 7 figures and 4 tables

Journal ref: Phys. Rev. D 98, 063008 (2018)

Analyzing the Gamma-ray Sky with Wavelets

Authors: Bhaskaran Balaji, Ilias Cholis, Patrick J. Fox, Samuel D. McDermott

Abstract: We analyze the gamma-ray sky at energies of 0.5 to 50 GeV using the undecimated wavelet transform on the sphere. Focusing on the inner $60^{\circ} \times 60^{\circ}$ of the sky, we identify and characterize four separate residuals beyond the expected Milky Way diffuse emission. We detect the \textit{Fermi} Bubbles, finding compelling evidence that they are diffuse in nature and contain very little… ▽ More We analyze the gamma-ray sky at energies of 0.5 to 50 GeV using the undecimated wavelet transform on the sphere. Focusing on the inner $60^{\circ} \times 60^{\circ}$ of the sky, we identify and characterize four separate residuals beyond the expected Milky Way diffuse emission. We detect the \textit{Fermi} Bubbles, finding compelling evidence that they are diffuse in nature and contain very little small-scale structure. We detect the "cocoon" inside the Southern Bubble, and we also identify its northern counterpart above 2 GeV. The Northern Cocoon lies along the same axis but is $\sim 30 \%$ dimmer than the southern one. We characterize the Galactic center excess, which we find extends up to $20^{\circ}$ in $|b|$. At latitudes $|b| \leq 5^{\circ}$ we find evidence for power in small angular scales that could be the result of point-source contributions, but for $|b| \geq 5^{\circ}$ the Galactic center excess is dominantly diffuse in its nature. Our findings show that either the Galactic center excess and {\it Fermi} Bubbles connect smoothly or that the Bubbles brighten significantly below $15^\circ$ in latitude. We find that the Galactic center excess appears off-center by a few degrees towards negative $\ell$. Additionally, we find and characterize two emissions along the Galactic disk centered at $\ell \simeq +25^{\circ}$ and $-20^{\circ}$. These emissions are significantly more elongated along the Galactic disk than the Galactic center excess. △ Less

Submitted 5 March, 2018; originally announced March 2018.

Comments: 23 pages, 25 figures, 4 tables, 3 appendices

Report number: FERMILAB-PUB-18-057-A-T

Journal ref: Phys. Rev. D 98, 043009 (2018)

Limits on Runaway Growth of Intermediate Mass Black Holes from Advanced LIGO

Authors: Ely D. Kovetz, Ilias Cholis, Marc Kamionkowski, Joseph Silk

Abstract: There is growing evidence that intermediate-mass black holes (IMBHs), defined here as having a mass in the range M=500-10^5 Msun, are present in the dense centers of certain globular clusters (GCs). Gravitational waves (GWs) from their mergers with other IMBHs or with stellar BHs in the cluster are mostly emitted in frequencies <10 Hz, which unfortunately is out of reach for current ground-based o… ▽ More There is growing evidence that intermediate-mass black holes (IMBHs), defined here as having a mass in the range M=500-10^5 Msun, are present in the dense centers of certain globular clusters (GCs). Gravitational waves (GWs) from their mergers with other IMBHs or with stellar BHs in the cluster are mostly emitted in frequencies <10 Hz, which unfortunately is out of reach for current ground-based observatories such as advanced LIGO (aLIGO). Nevertheless, we show that aLIGO measurements can be used to efficiently probe one of the possible formation mechanisms of IMBHs in GCs, namely a runaway merger process of stellar seed BHs. In this case, aLIGO will be sensitive to the lower-mass rungs of the merger ladder, ranging from the seed BH mass to masses >~50-300 Msun, where the background from standard mergers is expected to be very low. Assuming this generic IMBH formation scenario, we calculate the mass functions that correspond to the limiting cases of possible merger trees. Based on estimates for the number density of GCs and taking into account the instrumental sensitivity, we show that current observations do not effectively limit the occupation fraction f_occ of IMBHs formed by runaway mergers of stellar BHs in GCs. However, we find that a six-year run of aLIGO at design sensitivity will be able to probe down to f_occ<3% at a 99.9% confidence level, either finding evidence for this formation mechanism, or necessitating others if the fraction of GCs that harbor IMBHs is higher. △ Less

Submitted 1 March, 2018; originally announced March 2018.

Comments: 13 pages, 5 figures, 2 tables

Journal ref: Phys. Rev. D 97, 123003 (2018)

Features in the Spectrum of Cosmic-Ray Positrons from Pulsars

Authors: Ilias Cholis, Tanvi Karwal, Marc Kamionkowski

Abstract: Pulsars have been invoked to explain the origin of recently observed high-energy Galactic cosmic-ray positrons. Since the positron propagation distance decreases with energy, the number of pulsars that can contribute to the observed positrons decreases from $O(10^3)$ for positron energies $E\gtrsim10$ GeV to only a few for $E \gtrsim 500$ GeV. Thus, if pulsars explain these positrons, the positron… ▽ More Pulsars have been invoked to explain the origin of recently observed high-energy Galactic cosmic-ray positrons. Since the positron propagation distance decreases with energy, the number of pulsars that can contribute to the observed positrons decreases from $O(10^3)$ for positron energies $E\gtrsim10$ GeV to only a few for $E \gtrsim 500$ GeV. Thus, if pulsars explain these positrons, the positron energy spectrum should become increasingly bumpy at higher energies. Here we present a power-spectrum analysis that can be applied to seek such spectral features in the energy spectrum for cosmic-ray positrons and for the energy spectrum of the combined electron/positron flux. We account for uncertainties in the pulsar distribution by generating hundreds of simulated spectra from pulsar distributions consistent with current observational constraints. Although the current AMS-02 data do not exhibit evidence for spectral features, we find that such features would be detectable in $\simeq 10\%$ of our simulations, with 20 years of AMS-02 data or three years of DAMPE measurements on the electron-plus-positron flux. △ Less

Submitted 17 July, 2018; v1 submitted 30 November, 2017; originally announced December 2017.

Comments: 9 pages, 6 Figures, 3 Tables, 5 Appendices; version 2: citations added, small chances to reflect published version including added appendix E

Journal ref: Phys. Rev. D 97, 123011 (2018)

TeV Gamma Rays From Galactic Center Pulsars

Authors: Dan Hooper, Ilias Cholis, Tim Linden

Abstract: Measurements of the nearby pulsars Geminga and B0656+14 by the HAWC and Milagro telescopes have revealed the presence of bright TeV-emitting halos surrounding these objects. If young and middle-aged pulsars near the Galactic Center transfer a similar fraction of their energy into TeV photons, then these sources could plausibly dominate the emission that is observed by HESS and other ground-based t… ▽ More Measurements of the nearby pulsars Geminga and B0656+14 by the HAWC and Milagro telescopes have revealed the presence of bright TeV-emitting halos surrounding these objects. If young and middle-aged pulsars near the Galactic Center transfer a similar fraction of their energy into TeV photons, then these sources could plausibly dominate the emission that is observed by HESS and other ground-based telescopes from the innermost $\sim$$10^2$ parsecs of the Milky Way. In particular, both the spectral shape and the angular extent of this emission is consistent with TeV halos produced by a population of pulsars, although the reported correlation of this emission with the distribution of molecular gas suggests that diffuse hadronic processes also must contribute. The overall flux of this emission requires a birth rate of ~100-1000 neutron stars per Myr near the Galactic Center, in good agreement with recent estimates. △ Less

Submitted 11 December, 2017; v1 submitted 25 May, 2017; originally announced May 2017.

Comments: 13 pages, 2 figure

Report number: FERMILAB-PUB-17-173-A

Using HAWC to Discover Invisible Pulsars

Authors: Tim Linden, Katie Auchettl, Joseph Bramante, Ilias Cholis, Ke Fang, Dan Hooper, Tanvi Karwal, Shirley Weishi Li

Abstract: Observations by HAWC and Milagro have detected bright and spatially extended TeV gamma-ray sources surrounding the Geminga and Monogem pulsars. We argue that these observations, along with a substantial population of other extended TeV sources coincident with pulsar wind nebulae, constitute a new morphological class of spatially extended TeV halos. We show that HAWCs wide field-of-view unlocks an… ▽ More Observations by HAWC and Milagro have detected bright and spatially extended TeV gamma-ray sources surrounding the Geminga and Monogem pulsars. We argue that these observations, along with a substantial population of other extended TeV sources coincident with pulsar wind nebulae, constitute a new morphological class of spatially extended TeV halos. We show that HAWCs wide field-of-view unlocks an expansive parameter space of TeV halos not observable by atmospheric Cherenkov telescopes. Under the assumption that Geminga and Monogem are typical middle-aged pulsars, we show that ten-year HAWC observations should eventually observe 37$^{+17}_{-13}$ middle-aged TeV halos that correspond to pulsars whose radio emission is not beamed towards Earth. Depending on the extrapolation of the TeV halo efficiency to young pulsars, HAWC could detect more than 100 TeV halos from mis-aligned pulsars. These pulsars have historically been difficult to detect with existing multiwavelength observations. TeV halos will constitute a significant fraction of all HAWC sources, allowing follow-up observations to efficiently find pulsar wind nebulae and thermal pulsar emission. The observation and subsequent multi-wavelength follow-up of TeV halos will have significant implications for our understanding of pulsar beam geometries, the evolution of PWN, the diffusion of cosmic-rays near energetic pulsars, and the contribution of pulsars to the cosmic-ray positron excess. △ Less

Submitted 28 March, 2017; originally announced March 2017.

Comments: 14 pages, 1 figure, To Be Submitted to PRD

Journal ref: Phys. Rev. D 96, 103016 (2017)

HAWC Observations Strongly Favor Pulsar Interpretations of the Cosmic-Ray Positron Excess

Authors: Dan Hooper, Ilias Cholis, Tim Linden, Ke Fang

Abstract: Recent measurements of the Geminga and B0656+14 pulsars by the gamma-ray telescope HAWC (along with earlier measurements by Milagro) indicate that these objects generate significant fluxes of very high-energy electrons. In this paper, we use the very high-energy gamma-ray intensity and spectrum of these pulsars to calculate and constrain their expected contributions to the local cosmic-ray positro… ▽ More Recent measurements of the Geminga and B0656+14 pulsars by the gamma-ray telescope HAWC (along with earlier measurements by Milagro) indicate that these objects generate significant fluxes of very high-energy electrons. In this paper, we use the very high-energy gamma-ray intensity and spectrum of these pulsars to calculate and constrain their expected contributions to the local cosmic-ray positron spectrum. Among models that are capable of reproducing the observed characteristics of the gamma-ray emission, we find that pulsars invariably produce a flux of high-energy positrons that is similar in spectrum and magnitude to the positron fraction measured by PAMELA and AMS-02. In light of this result, we conclude that it is very likely that pulsars provide the dominant contribution to the long perplexing cosmic-ray positron excess. △ Less

Submitted 27 February, 2017; originally announced February 2017.

Comments: 13 pages, 6 figures, To Be Submitted to JCAP

Journal ref: Phys. Rev. D 96, 103013 (2017)

Possible Evidence for the Stochastic Acceleration of Secondary Antiprotons by Supernova Remnants

Authors: Ilias Cholis, Dan Hooper, Tim Linden

Abstract: The antiproton-to-proton ratio in the cosmic-ray spectrum is a sensitive probe of new physics. Using recent measurements of the cosmic-ray antiproton and proton fluxes in the energy range of 1-1000 GeV, we study the contribution to the $\bar{p}/p$ ratio from secondary antiprotons that are produced and subsequently accelerated within individual supernova remnants. We consider several well-motivated… ▽ More The antiproton-to-proton ratio in the cosmic-ray spectrum is a sensitive probe of new physics. Using recent measurements of the cosmic-ray antiproton and proton fluxes in the energy range of 1-1000 GeV, we study the contribution to the $\bar{p}/p$ ratio from secondary antiprotons that are produced and subsequently accelerated within individual supernova remnants. We consider several well-motivated models for cosmic-ray propagation in the interstellar medium and marginalize our results over the uncertainties related to the antiproton production cross section and the time-, charge-, and energy-dependent effects of solar modulation. We find that the increase in the $\bar{p}/p$ ratio observed at rigidities above $\sim$ 100 GV cannot be accounted for within the context of conventional cosmic-ray propagation models, but is consistent with scenarios in which cosmic-ray antiprotons are produced and subsequently accelerated by shocks within a given supernova remnant. In light of this, the acceleration of secondary cosmic rays in supernova remnants is predicted to substantially contribute to the cosmic-ray positron spectrum, accounting for a significant fraction of the observed positron excess. △ Less

Submitted 15 June, 2017; v1 submitted 16 January, 2017; originally announced January 2017.

Comments: 5 pages, 3 figures and 4 tables. Changes in v2: extended discussion on the method, title and some text revisions, all in agreement with final published PRD version. Results unchanged

Report number: FERMILAB-PUB-17-010-A

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

The Black Hole Mass Function from Gravitational Wave Measurements

Authors: Ely D. Kovetz, Ilias Cholis, Patrick C. Breysse, Marc Kamionkowski

Abstract: We examine how future gravitational-wave measurements from merging black holes (BHs) can be used to infer the shape of the black-hole mass function, with important implications for the study of star formation and evolution and the properties of binary BHs. We model the mass function as a power law, inherited from the stellar initial mass function, and introduce lower and upper mass cutoff paramete… ▽ More We examine how future gravitational-wave measurements from merging black holes (BHs) can be used to infer the shape of the black-hole mass function, with important implications for the study of star formation and evolution and the properties of binary BHs. We model the mass function as a power law, inherited from the stellar initial mass function, and introduce lower and upper mass cutoff parameterizations in order to probe the minimum and maximum BH masses allowed by stellar evolution, respectively. We initially focus on the heavier BH in each binary, to minimize model dependence. Taking into account the experimental noise, the mass measurement errors and the uncertainty in the redshift-dependence of the merger rate, we show that the mass function parameters, as well as the total rate of merger events, can be measured to <10% accuracy within a few years of advanced LIGO observations at its design sensitivity. This can be used to address important open questions such as the upper limit on the stellar mass which allows for BH formation and to confirm or refute the currently observed mass gap between neutron stars and BHs. In order to glean information on the progenitors of the merging BH binaries, we then advocate the study of the two-dimensional mass distribution to constrain parameters that describe the two-body system, such as the mass ratio between the two BHs, in addition to the merger rate and mass function parameters. We argue that several years of data collection can efficiently probe models of binary formation, and show, as an example, that the hypothesis that some gravitational-wave events may involve primordial black holes can be tested. Finally, we point out that in order to maximize the constraining power of the data, it may be worthwhile to lower the signal-to-noise threshold imposed on each candidate event and amass a larger statistical ensemble of BH mergers. △ Less

Submitted 3 November, 2016; originally announced November 2016.

Comments: 12 pages, 7 figures; To be submitted to PRD; Comments welcome

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

On the Gravitational Wave Background from Black Hole Binaries after the First LIGO Detections

Authors: Ilias Cholis

Abstract: The detection of gravitational waves from the merger of binary black holes by the LIGO Collaboration has opened a new window to astrophysics. With the sensitivities of ground based detectors in the coming years we can only detect the local black hole binary mergers. The integrated merger rate can instead be probed by the gravitational-wave background, the incoherent superposition of the released e… ▽ More The detection of gravitational waves from the merger of binary black holes by the LIGO Collaboration has opened a new window to astrophysics. With the sensitivities of ground based detectors in the coming years we can only detect the local black hole binary mergers. The integrated merger rate can instead be probed by the gravitational-wave background, the incoherent superposition of the released energy in gravitational waves during binary-black-hole coalescence. Through that, the properties of the binary black holes can be studied. In this work we show that by measuring the energy density $Ω_{GW}$ (in units of the cosmic critical density) of the gravitational-wave background, we can search for the rare $\sim 100 M_{\odot}$ massive black holes formed in the Universe. In addition, we can answer how often the least massive BHs of mass $> 3 M_{\odot}$ form. Finally, if there are multiple channels for the formation of binary black holes and if any of them predicts a narrow mass range for the black holes, then the gravitational-wave-background spectrum may have features that with the future Einstein Telescope can be detected. △ Less

Submitted 28 April, 2017; v1 submitted 12 September, 2016; originally announced September 2016.

Comments: 7 pages, 5 figures; version 2: citations added, minor chances

Stochastic Gravitational-Wave Background due to Primordial Binary Black Hole Mergers

Authors: Vuk Mandic, Simeon Bird, Ilias Cholis

Abstract: Recent Advanced LIGO detections of binary black hole mergers have prompted multiple studies investigating the possibility that the heavy GW150914 binary system was of primordial origin, and hence could be evidence for dark matter in the form of black holes. We compute the stochastic background arising from the incoherent superposition of such primordial binary black hole systems in the universe an… ▽ More Recent Advanced LIGO detections of binary black hole mergers have prompted multiple studies investigating the possibility that the heavy GW150914 binary system was of primordial origin, and hence could be evidence for dark matter in the form of black holes. We compute the stochastic background arising from the incoherent superposition of such primordial binary black hole systems in the universe and compare it to the similar background spectrum due to binary black hole systems of stellar origin. We investigate the possibility of detecting this background with future gravitational wave detectors, and discuss the possibility of using the stochastic gravitational-wave background measurement to constrain the dark matter component in the form of black holes. △ Less

Submitted 23 August, 2016; originally announced August 2016.

Comments: 6 pages, 5 figures

Report number: LIGO-P1600241

Journal ref: Phys. Rev. Lett. 117, 201102 (2016)

Orbital eccentricities in primordial black holes binaries

Authors: Ilias Cholis, Ely D. Kovetz, Yacine Ali-Haïmoud, Simeon Bird, Marc Kamionkowski, Julian B. Muñoz, Alvise Raccanelli

Abstract: It was recently suggested that the merger of $\sim30\,M_\odot$ primordial black holes (PBHs) may provide a significant number of events in gravitational-wave observatories over the next decade, if they make up an appreciable fraction of the dark matter. Here we show that measurement of the eccentricities of the inspiralling binary black holes can be used to distinguish these binaries from those pr… ▽ More It was recently suggested that the merger of $\sim30\,M_\odot$ primordial black holes (PBHs) may provide a significant number of events in gravitational-wave observatories over the next decade, if they make up an appreciable fraction of the dark matter. Here we show that measurement of the eccentricities of the inspiralling binary black holes can be used to distinguish these binaries from those produced by more traditional astrophysical mechanisms. These PBH binaries are formed on highly eccentric orbits and can then merge on timescales that in some cases are years or less, retaining some eccentricity in the last seconds before the merger. This is to be contrasted with massive-stellar-binary, globular-cluster, or other astrophysical origins for binary black holes (BBHs) in which the orbits have very effectively circularized by the time the BBH enters the observable LIGO window. Here we discuss the features of the gravitational-wave signals that indicate this eccentricity and forecast the sensitivity of LIGO and the Einstein Telescope to such effects. We show that if PBHs make up the dark matter, then roughly one event should have a detectable eccentricity given LIGO's expected sensitivity and observing time of six years. The Einstein Telescope should see $O(10)$ such events after ten years. △ Less

Submitted 23 June, 2016; originally announced June 2016.

Comments: 13 pages, 8 figures, 3 tables and 1 appendix

Journal ref: Phys. Rev. D 94, 084013 (2016)

Determining the progenitors of merging black-hole binaries

Authors: Alvise Raccanelli, Ely D. Kovetz, Simeon Bird, Ilias Cholis, Julian B. Munoz

Abstract: We investigate a possible method for determining the progenitors of black hole (BH) mergers observed via their gravitational wave (GW) signal. We argue that measurements of the cross-correlation of the GW events with overlapping galaxy catalogs may provide an additional tool in determining if BH mergers trace the stellar mass of the Universe, as would be expected from mergers of the endpoints of s… ▽ More We investigate a possible method for determining the progenitors of black hole (BH) mergers observed via their gravitational wave (GW) signal. We argue that measurements of the cross-correlation of the GW events with overlapping galaxy catalogs may provide an additional tool in determining if BH mergers trace the stellar mass of the Universe, as would be expected from mergers of the endpoints of stellar evolution. If on the other hand the BHs are of primordial origin, as has been recently suggested, their merging would be preferentially hosted by lower biased objects, and thus have a lower cross-correlation with luminous galaxies. Here we forecast the expected precision of the cross-correlation measurement for current and future GW detectors such as LIGO and the Einstein Telescope. We then predict how well these instruments can distinguish the model that identifies high-mass BH-BH mergers as the merger of primordial black holes that constitute the dark matter in the Universe from more traditional astrophysical sources. △ Less

Submitted 4 May, 2016; originally announced May 2016.

Comments: 9 pages, 6 figures

Journal ref: Phys. Rev. D 94, 023516 (2016)

Did LIGO detect dark matter?

Authors: Simeon Bird, Ilias Cholis, Julian B. Muñoz, Yacine Ali-Haïmoud, Marc Kamionkowski, Ely D. Kovetz, Alvise Raccanelli, Adam G. Riess

Abstract: We consider the possibility that the black-hole (BH) binary detected by LIGO may be a signature of dark matter. Interestingly enough, there remains a window for masses $20\,M_\odot \lesssim M_{\rm bh} \lesssim 100\, M_\odot$ where primordial black holes (PBHs) may constitute the dark matter. If two BHs in a galactic halo pass sufficiently close, they radiate enough energy in gravitational waves to… ▽ More We consider the possibility that the black-hole (BH) binary detected by LIGO may be a signature of dark matter. Interestingly enough, there remains a window for masses $20\,M_\odot \lesssim M_{\rm bh} \lesssim 100\, M_\odot$ where primordial black holes (PBHs) may constitute the dark matter. If two BHs in a galactic halo pass sufficiently close, they radiate enough energy in gravitational waves to become gravitationally bound. The bound BHs will rapidly spiral inward due to emission of gravitational radiation and ultimately merge. Uncertainties in the rate for such events arise from our imprecise knowledge of the phase-space structure of galactic halos on the smallest scales. Still, reasonable estimates span a range that overlaps the $2-53$ Gpc$^{-3}$ yr$^{-1}$ rate estimated from GW150914, thus raising the possibility that LIGO has detected PBH dark matter. PBH mergers are likely to be distributed spatially more like dark matter than luminous matter and have no optical nor neutrino counterparts. They may be distinguished from mergers of BHs from more traditional astrophysical sources through the observed mass spectrum, their high ellipticities, or their stochastic gravitational wave background. Next generation experiments will be invaluable in performing these tests. △ Less

Submitted 30 May, 2016; v1 submitted 1 March, 2016; originally announced March 2016.

Comments: 5 pages, 2 figures, updated to match version published in PRL

Journal ref: Phys. Rev. Lett. 116, 201301 (2016)

Wavelet-Based Techniques for the Gamma-Ray Sky

Authors: Samuel D. McDermott, Patrick J. Fox, Ilias Cholis, Samuel K. Lee

Abstract: We demonstrate how the image analysis technique of wavelet decomposition can be applied to the gamma-ray sky to separate emission on different angular scales. New structures on scales that differ from the scales of the conventional astrophysical foreground and background uncertainties can be robustly extracted, allowing a model-independent characterization with no presumption of exact signal morph… ▽ More We demonstrate how the image analysis technique of wavelet decomposition can be applied to the gamma-ray sky to separate emission on different angular scales. New structures on scales that differ from the scales of the conventional astrophysical foreground and background uncertainties can be robustly extracted, allowing a model-independent characterization with no presumption of exact signal morphology. As a test case, we generate mock gamma-ray data to demonstrate our ability to extract extended signals without assuming a fixed spatial template. For some point source luminosity functions, our technique also allows us to differentiate a diffuse signal in gamma-rays from dark matter annihilation and extended gamma-ray point source populations in a data-driven way. △ Less

Submitted 30 November, 2015; originally announced December 2015.

Comments: 18+8 pages, 7 figures

Report number: FERMILAB-PUB-15-504-T, YITP-SB-15-43

A Predictive Analytic Model for the Solar Modulation of Cosmic Rays

Authors: Ilias Cholis, Dan Hooper, Tim Linden

Abstract: An important factor limiting our ability to understand the production and propagation of cosmic rays pertains to the effects of heliospheric forces, commonly known as solar modulation. The solar wind is capable of generating time and charge-dependent effects on the spectrum and intensity of low energy ($\lsim$ 10 GeV) cosmic rays reaching Earth. Previous analytic treatments of solar modulation hav… ▽ More An important factor limiting our ability to understand the production and propagation of cosmic rays pertains to the effects of heliospheric forces, commonly known as solar modulation. The solar wind is capable of generating time and charge-dependent effects on the spectrum and intensity of low energy ($\lsim$ 10 GeV) cosmic rays reaching Earth. Previous analytic treatments of solar modulation have utilized the force-field approximation, in which a simple potential is adopted whose amplitude is selected to best fit the cosmic-ray data taken over a given period of time. Making use of recently available cosmic-ray data from the Voyager 1 spacecraft, along with measurements of the heliospheric magnetic field and solar wind, we construct a time, charge and rigidity-dependent model of solar modulation that can be directly compared to data from a variety of cosmic-ray experiments. We provide a simple analytic formula that can be easily utilized in a variety of applications, allowing us to better predict the effects of solar modulation and reduce the number of free parameters involved in cosmic ray propagation models. △ Less

Submitted 25 November, 2016; v1 submitted 4 November, 2015; originally announced November 2015.

Comments: 15 pages, 7 figures and 2 tables, changes in v2: in agreement with PRD published version, small revisions in text, results and conclusions unchanged, changes in v3: fixed typo with regards to Eqs. 2 & 12 and extended Table II to include data up to 2016

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

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

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

The GeV Excess Shining Through: Background Systematics for the Inner Galaxy Analysis

Authors: Francesca Calore, Ilias Cholis, Christoph Weniger

Abstract: Recently, a spatially extended excess of gamma rays collected by the Fermi-LAT from the inner region of the Milky Way has been detected by different groups and with increasingly sophisticated techniques. Yet, any final conclusion about the morphology and spectral properties of such an extended diffuse emission are subject to a number of potentially critical uncertainties, related to the high densi… ▽ More Recently, a spatially extended excess of gamma rays collected by the Fermi-LAT from the inner region of the Milky Way has been detected by different groups and with increasingly sophisticated techniques. Yet, any final conclusion about the morphology and spectral properties of such an extended diffuse emission are subject to a number of potentially critical uncertainties, related to the high density of cosmic rays, gas, magnetic fields and abundance of point sources. We will present a thorough study of the systematic uncertainties related to the modelling of diffuse background and to the propagation of cosmic rays in the inner part of our Galaxy. We will test a large set of models for the Galactic diffuse emission, generated by varying the propagation parameters within extreme conditions. By using those models in the analysis of Fermi-LAT data as Galactic foreground, we will show that the gamma-ray excess survives and we will quantify the uncertainties affecting the excess morphology and energy spectrum. △ Less

Submitted 10 February, 2015; originally announced February 2015.

Comments: 2014 Fermi Symposium proceedings - eConf C14102.1 7 pages, 4 figures

Report number: FERMILAB-CONF-15-037-A

A Tale of Tails: Dark Matter Interpretations of the Fermi GeV Excess in Light of Background Model Systematics

Authors: Francesca Calore, Ilias Cholis, Christopher McCabe, Christoph Weniger

Abstract: Several groups have identified an extended excess of gamma rays over the modeled foreground and background emissions towards the Galactic center (GC) based on observations with the Fermi Large Area Telescope. This excess emission is compatible in morphology and spectrum with a telltale sign from dark matter (DM) annihilation. Here, we present a critical reassessment of DM interpretations of the GC… ▽ More Several groups have identified an extended excess of gamma rays over the modeled foreground and background emissions towards the Galactic center (GC) based on observations with the Fermi Large Area Telescope. This excess emission is compatible in morphology and spectrum with a telltale sign from dark matter (DM) annihilation. Here, we present a critical reassessment of DM interpretations of the GC signal in light of the foreground and background uncertainties that some of us recently outlaid in Calore et al. 2014. We find that a much larger number of DM models fits the gamma-ray data than previously noted. In particular: (1) In the case of DM annihilation into $\bar{b}b$, we find that even large DM masses up to $m_χ\simeq$ 74 GeV are allowed with a $p$-value $> 0.05$. (2) Surprisingly, annihilation into non-relativistic hh gives a good fit to the data. (3) The inverse Compton emission from $μ^+μ^-$ with $m_χ\sim$ 60-70 GeV can also account for the excess at higher latitudes, $|b|>2^\circ$, both in its spectrum and morphology. We also present novel constraints on a large number of mixed annihilation channels, including cascade annihilation involving hidden sector mediators. Finally, we show that the current limits from dwarf spheroidal observations are not in tension with a DM interpretation when uncertainties on the DM halo profile are accounted for. △ Less

Submitted 26 April, 2015; v1 submitted 17 November, 2014; originally announced November 2014.

Comments: 18 pages, 9 figures, 4 tables. It matches the published version. Figure 1 corrected

Report number: FERMILAB-PUB-14-477-A

Journal ref: Phys. Rev. D 91, 063003 (2015)

Background model systematics for the Fermi GeV excess

Authors: Francesca Calore, Ilias Cholis, Christoph Weniger

Abstract: The possible gamma-ray excess in the inner Galaxy and the Galactic center (GC) suggested by Fermi-LAT observations has triggered a large number of studies. It has been interpreted as a variety of different phenomena such as a signal from WIMP dark matter annihilation, gamma-ray emission from a population of millisecond pulsars, or emission from cosmic rays injected in a sequence of burst-like even… ▽ More The possible gamma-ray excess in the inner Galaxy and the Galactic center (GC) suggested by Fermi-LAT observations has triggered a large number of studies. It has been interpreted as a variety of different phenomena such as a signal from WIMP dark matter annihilation, gamma-ray emission from a population of millisecond pulsars, or emission from cosmic rays injected in a sequence of burst-like events or continuously at the GC. We present the first comprehensive study of model systematics coming from the Galactic diffuse emission in the inner part of our Galaxy and their impact on the inferred properties of the excess emission at Galactic latitudes $2^\circ<|b|<20^\circ$ and 300 MeV to 500 GeV. We study both theoretical and empirical model systematics, which we deduce from a large range of Galactic diffuse emission models and a principal component analysis of residuals in numerous test regions along the Galactic plane. We show that the hypothesis of an extended spherical excess emission with a uniform energy spectrum is compatible with the Fermi-LAT data in our region of interest at $95\%$ CL. Assuming that this excess is the extended counterpart of the one seen in the inner few degrees of the Galaxy, we derive a lower limit of $10.0^\circ$ ($95\%$ CL) on its extension away from the GC. We show that, in light of the large correlated uncertainties that affect the subtraction of the Galactic diffuse emission in the relevant regions, the energy spectrum of the excess is equally compatible with both a simple broken power-law of break energy $2.1\pm0.2$ GeV, and with spectra predicted by the self-annihilation of dark matter, implying in the case of $\bar{b}b$ final states a dark matter mass of $49^{+6.4}_{-5.4}$ GeV. △ Less

Submitted 29 August, 2014; originally announced September 2014.

Comments: 65 pages, 28 figures, 7 tables

Report number: FERMILAB-PUB-14-289-A

A Critical Reevaluation of Radio Constraints on Annihilating Dark Matter

Authors: Ilias Cholis, Dan Hooper, Tim Linden

Abstract: A number of groups have employed radio observations of the Galactic center to derive stringent constraints on the annihilation cross section of weakly interacting dark matter. In this paper, we show that electron energy losses in this region are likely to be dominated by inverse Compton scattering on the interstellar radiation field, rather than by synchrotron, considerably relaxing the constraint… ▽ More A number of groups have employed radio observations of the Galactic center to derive stringent constraints on the annihilation cross section of weakly interacting dark matter. In this paper, we show that electron energy losses in this region are likely to be dominated by inverse Compton scattering on the interstellar radiation field, rather than by synchrotron, considerably relaxing the constraints on the dark matter annihilation cross section compared to previous works. Strong convective winds, which are well motivated by recent observations, may also significantly weaken synchrotron constraints. After taking these factors into account, we find that radio constraints on annihilating dark matter are orders of magnitude less stringent than previously reported, and are generally weaker than those derived from current gamma-ray observations. △ Less

Submitted 5 February, 2016; v1 submitted 26 August, 2014; originally announced August 2014.

Comments: 6 pages, 3 Figures, arXiv version updated to reflect published version

Journal ref: Phys. Rev. D 91, 083507 (2015)

Challenges in Explaining the Galactic Center Gamma-Ray Excess with Millisecond Pulsars

Authors: Ilias Cholis, Dan Hooper, Tim Linden

Abstract: Millisecond pulsars have been discussed as a possible source of the gamma-ray excess observed from the region surrounding the Galactic Center. With this in mind, we use the observed population of bright low-mass X-ray binaries to estimate the number of millisecond pulsars in the Inner Galaxy. This calculation suggests that only ~1-5% of the excess is produced by millisecond pulsars. We also use th… ▽ More Millisecond pulsars have been discussed as a possible source of the gamma-ray excess observed from the region surrounding the Galactic Center. With this in mind, we use the observed population of bright low-mass X-ray binaries to estimate the number of millisecond pulsars in the Inner Galaxy. This calculation suggests that only ~1-5% of the excess is produced by millisecond pulsars. We also use the luminosity function derived from local measurements of millisecond pulsars, along with the number of point sources resolved by Fermi, to calculate an upper limit for the diffuse emission from such a population. While this limit is compatible with the millisecond pulsar population implied by the number of low-mass X-ray binaries, it strongly excludes the possibility that most of the excess originates from such objects. △ Less

Submitted 26 February, 2016; v1 submitted 21 July, 2014; originally announced July 2014.

Comments: 5pages, 2 figures, changes in v2: in agreement with JCAP published version, Fig. 1 updated, small revisions in text, results and conclusions unchanged

Report number: FERMILAB-PUB-14-240-A

Journal ref: JCAP 1506 (2015) 06, 043

A New Determination of the Spectra and Luminosity Function of Gamma-Ray Millisecond Pulsars

Authors: Ilias Cholis, Dan Hooper, Tim Linden

Abstract: In this article, we revisit the gamma-ray emission observed from millisecond pulsars and globular clusters. Based on 5.6 years of data from the Fermi Gamma-Ray Space Telescope, we report gamma-ray spectra for 61 millisecond pulsars, finding most to be well fit by a power-law with an exponential cutoff, producing to a spectral peak near ~1-2 GeV (in $E^2 dN/dE$ units). Additionally, while most glob… ▽ More In this article, we revisit the gamma-ray emission observed from millisecond pulsars and globular clusters. Based on 5.6 years of data from the Fermi Gamma-Ray Space Telescope, we report gamma-ray spectra for 61 millisecond pulsars, finding most to be well fit by a power-law with an exponential cutoff, producing to a spectral peak near ~1-2 GeV (in $E^2 dN/dE$ units). Additionally, while most globular clusters exhibit a similar spectral shape, we identify a few with significantly softer spectra. We also determine the gamma-ray luminosity function of millisecond pulsars using the population found in the nearby field of the Milky Way, and within the globular cluster 47 Tucanae. We find that the gamma-ray emission observed from globular clusters is dominated by a relatively small number of bright millisecond pulsars, and that low-luminosity pulsars account for only a small fraction of the total flux. Our results also suggest that the gamma-ray emission from millisecond pulsars is more isotropic and less strongly beamed than the emission at X-ray wavelengths. Furthermore, the observed distribution of apparent gamma-ray efficiencies provides support for the slot gap or the outer gap models over those in which the gamma-ray emission originates from regions close to the neutron star's magnetic poles (polar cap models). △ Less

Submitted 21 July, 2014; originally announced July 2014.

Comments: 28 pages, 21 figures, 4 tables

Report number: FERMILAB-PUB-14-236-A

Indirect Detection Analysis: Wino Dark Matter Case Study

Authors: Andrzej Hryczuk, Ilias Cholis, Roberto Iengo, Maryam Tavakoli, Piero Ullio

Abstract: We perform a multichannel analysis of the indirect signals for the Wino Dark Matter, including one-loop electroweak and Sommerfeld enhancement corrections. We derive limits from cosmic ray antiprotons and positrons, from continuum galactic and extragalactic diffuse gamma-ray spectra, from the absence of gamma-ray line features at the galactic center above 500 GeV in energy, from gamma-rays toward… ▽ More We perform a multichannel analysis of the indirect signals for the Wino Dark Matter, including one-loop electroweak and Sommerfeld enhancement corrections. We derive limits from cosmic ray antiprotons and positrons, from continuum galactic and extragalactic diffuse gamma-ray spectra, from the absence of gamma-ray line features at the galactic center above 500 GeV in energy, from gamma-rays toward nearby dwarf spheroidal galaxies and galaxy clusters, and from CMB power-spectra. Additionally, we show the future prospects for neutrino observations toward the inner Galaxy and from antideuteron searches. For each of these indirect detection probes we include and discuss the relevance of the most important astrophysical uncertainties that can impact the strength of the derived limits. We find that the Wino as a dark matter candidate is excluded in the mass range bellow ~800 GeV from antiprotons and between 1.8 and 3.5 TeV from the absence of a gamma-ray line feature toward the galactic center. Limits from other indirect detection probes confirm the main bulk of the excluded mass ranges. △ Less

Submitted 28 August, 2014; v1 submitted 23 January, 2014; originally announced January 2014.

Comments: 38 pages, 17 figures, 4 tables; minor corrections, fig. 17 added discussing the dependence of the limits on the Wino contribution to the whole DM, matches the published version

Report number: FERMILAB-PUB-14-010-A, TUM-HEP 927/14

Journal ref: JCAP 07(2014)031

Cosmic Neutrino Pevatrons: A Brand New Pathway to Astronomy, Astrophysics, and Particle Physics

Authors: Luis A. Anchordoqui, Vernon Barger, Ilias Cholis, Haim Goldberg, Dan Hooper, Alexander Kusenko, John G. Learned, Danny Marfatia, Sandip Pakvasa, Thomas C. Paul, Thomas J. Weiler

Abstract: The announcement by the IceCube Collaboration of the observation of 28 cosmic neutrino candidates has been greeted with a great deal of justified excitement. The data reported so far depart by 4.3σfrom the expected atmospheric neutrino background, which raises the obvious question: "Where in the Cosmos are these neutrinos coming from?" We review the many possibilities which have been explored in t… ▽ More The announcement by the IceCube Collaboration of the observation of 28 cosmic neutrino candidates has been greeted with a great deal of justified excitement. The data reported so far depart by 4.3σfrom the expected atmospheric neutrino background, which raises the obvious question: "Where in the Cosmos are these neutrinos coming from?" We review the many possibilities which have been explored in the literature to address this question, including origins at either Galactic or extragalactic celestial objects. For completeness, we also briefly discuss new physics processes which may either explain or be constrained by IceCube data. △ Less

Submitted 31 January, 2014; v1 submitted 23 December, 2013; originally announced December 2013.

Comments: This is a review article solicited for the inaugural edition of Journal of High Energy Astrophysics (JHEAp). Matching version accepted for publication

Report number: FERMILAB-PUB-13-541-A

Journal ref: JHEAp1:1,2014

Constraining the origin of the rising cosmic ray positron fraction with the boron-to-carbon ratio

Authors: Ilias Cholis, Dan Hooper

Abstract: The rapid rise in the cosmic ray positron fraction above 10 GeV, as measured by PAMELA and AMS, suggests the existence of nearby primary sources of high energy positrons, such as pulsars or annihilating/decaying dark matter. In contrast, the spectrum of secondary positrons produced through the collisions of cosmic rays in the interstellar medium is predicted to fall rapidly with energy, and thus i… ▽ More The rapid rise in the cosmic ray positron fraction above 10 GeV, as measured by PAMELA and AMS, suggests the existence of nearby primary sources of high energy positrons, such as pulsars or annihilating/decaying dark matter. In contrast, the spectrum of secondary positrons produced through the collisions of cosmic rays in the interstellar medium is predicted to fall rapidly with energy, and thus is unable to account for the observed rise. It has been proposed, however, that secondary positrons could be produced and then accelerated in nearby supernova remnants, potentially explaining the observed rise, without the need of primary positron sources. Yet, if secondary positrons are accelerated in such shocks, other secondary cosmic ray species (such as boron nuclei, and antiprotons) will also be accelerated, leading to rises in the boron-to-carbon and antiproton-to-proton ratios. The measurements of the boron-to-carbon ratio by the PAMELA and AMS collaborations, however, show no sign of such a rise. With this new data in hand, we revisit the secondary acceleration scenario for the rising positron fraction. Assuming that the same supernova remnants accelerate both light nuclei (protons, helium) and heavier cosmic ray species, we find that no more than ~25% of the observed rise in the positron fraction can result from this mechanism (at the 95% confidence level) △ Less

Submitted 4 March, 2014; v1 submitted 10 December, 2013; originally announced December 2013.

Comments: 9 pages, 5 figures; changes in v2: in agreement with PRD published version, small revisions in text, results and conclusions unchanged

Report number: FERMILAB-PUB-13-547-A

Journal ref: Phys. Rev. D 89, 043013 (2014)

Dissecting the Gamma-Ray Background in Search of Dark Matter

Authors: Ilias Cholis, Dan Hooper, Samuel D. McDermott

Abstract: Several classes of astrophysical sources contribute to the approximately isotropic gamma-ray background measured by the Fermi Gamma-Ray Space Telescope. In this paper, we use Fermi's catalog of gamma-ray sources (along with corresponding source catalogs at infrared and radio wavelengths) to build and constrain a model for the contributions to the extragalactic gamma-ray background from astrophysic… ▽ More Several classes of astrophysical sources contribute to the approximately isotropic gamma-ray background measured by the Fermi Gamma-Ray Space Telescope. In this paper, we use Fermi's catalog of gamma-ray sources (along with corresponding source catalogs at infrared and radio wavelengths) to build and constrain a model for the contributions to the extragalactic gamma-ray background from astrophysical sources, including radio galaxies, star-forming galaxies, and blazars. We then combine our model with Fermi's measurement of the gamma-ray background to derive constraints on the dark matter annihilation cross section, including contributions from both extragalactic and galactic halos and subhalos. The resulting constraints are competitive with the strongest current constraints from the Galactic Center and dwarf spheroidal galaxies. As Fermi continues to measure the gamma-ray emission from a greater number of astrophysical sources, it will become possible to more tightly constrain the astrophysical contributions to the extragalactic gamma-ray background. We project that with 10 years of data, Fermi's measurement of this background combined with the improved constraints on the astrophysical source contributions will yield a sensitivity to dark matter annihilations that exceeds the strongest current constraints by a factor of ~ 5 - 10. △ Less

Submitted 2 December, 2013; originally announced December 2013.

Comments: 15 pages, 17 figures

Report number: FERMILAB-PUB-13-546-A, MCTP-13-40

Constraints on dark matter annihilations from diffuse gamma-ray emission in the Galaxy

Authors: Maryam Tavakoli, Ilias Cholis, Carmelo Evoli, Piero Ullio

Abstract: Recent advances in gamma-ray cosmic ray, infrared and radio astronomy have allowed us to develop a significantly better understanding of the galactic medium properties in the last few years. In this work using the DRAGON code, that numerically solves the CR propagation equation and calculating gamma-ray emissivities in a 2-dimensional grid enclosing the Galaxy, we study in a self consistent manner… ▽ More Recent advances in gamma-ray cosmic ray, infrared and radio astronomy have allowed us to develop a significantly better understanding of the galactic medium properties in the last few years. In this work using the DRAGON code, that numerically solves the CR propagation equation and calculating gamma-ray emissivities in a 2-dimensional grid enclosing the Galaxy, we study in a self consistent manner models for the galactic diffuse gamma-ray emission. Our models are cross-checked to both the available CR and gamma-ray data. We address the extend to which dark matter annihilations in the Galaxy can contribute to the diffuse gamma-ray flux towards different directions on the sky. Moreover we discuss the impact that astrophysical uncertainties of non DM nature, have on the derived gamma-ray limits. Such uncertainties are related to the diffusion properties on the Galaxy, the interstellar gas and the interstellar radiation field energy densities. Light ~10 GeV dark matter annihilating dominantly to hadrons is more strongly constrained by gamma-ray observations towards the inner parts of the Galaxy and influenced the most by assumptions of the gas distribution; while TeV scale DM annihilating dominantly to leptons has its tightest constraints from observations towards the galactic center avoiding the galactic disk plane, with the main astrophysical uncertainty being the radiation field energy density. In addition, we present a method of deriving constraints on the dark matter distribution profile from the diffuse gamma-ray spectra. These results critically depend on the assumed mass of the dark matter particles and the type of its end annihilation products. △ Less

Submitted 4 March, 2014; v1 submitted 19 August, 2013; originally announced August 2013.

Comments: 34 pages, 11 figures, 3 appendices; v2 in agreement with JCAP published version, small revisions in text, results and conclusions unchanged

Report number: FERMILAB-PUB-13-320-A

Journal ref: JCAP 1401 (2014) 017