Showing 1–50 of 55 results for author: Macias, O

Multi-messenger modeling of the Monogem pulsar halo

Authors: Youyou Li, Oscar Macias, Shinichiro Ando, Jacco Vink

Abstract: The High-Altitude Water Cherenkov Telescope (HAWC) has detected TeV halos associated with two nearby pulsars/pulsar wind nebulae (PWN) -- Geminga and B0656+14. These TeV halos extend up to tens of pc from the central accelerators, indicating that the diffusion of ultrarelativistic electrons and positrons in the interstellar medium has been suppressed by two orders of magnitude. Although Geminga an… ▽ More The High-Altitude Water Cherenkov Telescope (HAWC) has detected TeV halos associated with two nearby pulsars/pulsar wind nebulae (PWN) -- Geminga and B0656+14. These TeV halos extend up to tens of pc from the central accelerators, indicating that the diffusion of ultrarelativistic electrons and positrons in the interstellar medium has been suppressed by two orders of magnitude. Although Geminga and B0656+14 are at similar distances and in the same field of view, they have distinct histories. Notably, B0656+14 probably still resides within its parent supernova remnant, the Monogem Ring, which can be observed in X-rays. In this work, we perform high-resolution simulations of the propagation and emission of relativistic lepton pairs around B0656+14 using a two-zone diffusion model using the GALPROP numerical code. We compared the predicted inverse-Compton spectrum to the observations made by HAWC and Fermi-LAT and found physically plausible model parameters that resulted in a good fit to the data. Additionally, we estimated the contribution of this TeV-halo to the positron flux observed on Earth and found it to be smaller than 10\% of the measured flux. We conclude that future observations of the TeV halo and its synchrotron emission counterpart in radio and X-ray frequencies will be crucial to distinguish between various possible models. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Dark Matter Line Searches with the Cherenkov Telescope Array

Authors: S. Abe, J. Abhir, A. Abhishek, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, L. Angel, C. Aramo, C. Arcaro, T. T. H. Arnesen, L. Arrabito, K. Asano, Y. Ascasibar, J. Aschersleben, H. Ashkar , et al. (540 additional authors not shown)

Abstract: Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of sele… ▽ More Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 43 pages JCAP style (excluding author list and references), 19 figures

Robust inference of the Galactic centre gamma-ray excess spatial properties

Authors: Deheng Song, Christopher Eckner, Chris Gordon, Francesca Calore, Oscar Macias, Kevork N. Abazajian, Shunsaku Horiuchi, Manoj Kaplinghat, Martin Pohl

Abstract: The gamma-ray Fermi-LAT Galactic centre excess (GCE) has puzzled scientists for over 15 years. Despite ongoing debates about its properties, and especially its spatial distribution, its nature remains elusive. We scrutinize how the estimated spatial morphology of this excess depends on models for the Galactic diffuse emission, focusing particularly on the extent to which the Galactic plane and poi… ▽ More The gamma-ray Fermi-LAT Galactic centre excess (GCE) has puzzled scientists for over 15 years. Despite ongoing debates about its properties, and especially its spatial distribution, its nature remains elusive. We scrutinize how the estimated spatial morphology of this excess depends on models for the Galactic diffuse emission, focusing particularly on the extent to which the Galactic plane and point sources are masked. Our main aim is to compare a spherically symmetric morphology - potentially arising from the annihilation of dark matter (DM) particles - with a boxy morphology - expected if faint unresolved sources in the Galactic bulge dominate the excess emission. Recent claims favouring a DM-motivated template for the GCE are shown to rely on a specific Galactic bulge template, which performs worse than other templates for the Galactic bulge. We find that a non-parametric model of the Galactic bulge derived from the VVV survey results in a significantly better fit for the GCE than DM-motivated templates. This result is independent of whether a GALPROP-based model or a more non-parametric ring-based model is used to describe the diffuse Galactic emission. This conclusion remains true even when additional freedom is added in the background models, allowing for non-parametric modulation of the model components and substantially improving the fit quality. When adopted, optimized background models provide robust results in terms of preference for a boxy bulge morphology of the GCE, regardless of the mask applied to the Galactic plane. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 21 pages, 11 figures

Report number: LAPTH-009/24

Journal ref: 2024MNRAS.530.4395S

Chasing Gravitational Waves with the Cherenkov Telescope Array

Authors: Jarred Gershon Green, Alessandro Carosi, Lara Nava, Barbara Patricelli, Fabian Schüssler, Monica Seglar-Arroyo, Cta Consortium, :, Kazuki Abe, Shotaro Abe, Atreya Acharyya, Remi Adam, Arnau Aguasca-Cabot, Ivan Agudo, Jorge Alfaro, Nuria Alvarez-Crespo, Rafael Alves Batista, Jean-Philippe Amans, Elena Amato, Filippo Ambrosino, Ekrem Oguzhan Angüner, Lucio Angelo Antonelli, Carla Aramo, Cornelia Arcaro, Luisa Arrabito , et al. (545 additional authors not shown)

Abstract: The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very… ▽ More The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA. △ Less

Submitted 5 February, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/30

TeV gamma-ray sensitivity to velocity-dependent dark matter models in the Galactic Center

Authors: Alessandro Montanari, Oscar Macias, Emmanuel Moulin

Abstract: The center of the Milky Way is a prime site to search for signals of dark matter (DM) annihilation due to its proximity and expected high concentration of DM. The amplification of the dispersion velocity of DM particles in the Galactic center (GC), caused by baryonic contraction and feedback, makes this particular region of the sky an even more promising target for exploring velocity-dependent DM… ▽ More The center of the Milky Way is a prime site to search for signals of dark matter (DM) annihilation due to its proximity and expected high concentration of DM. The amplification of the dispersion velocity of DM particles in the Galactic center (GC), caused by baryonic contraction and feedback, makes this particular region of the sky an even more promising target for exploring velocity-dependent DM models. Here we demonstrate that current GC observations with the H.E.S.S. telescope, presently the most sensitive TeV-scale gamma-ray telescope in operation in this region of the sky, set the strongest constraints on velocity-dependent annihilating DM particles with masses above 200 GeV. For p-wave annihilations, they improve the current constraints by a factor of $\sim$4 for a DM mass of 1 TeV. For the spatial distribution of DM, we use the results of the latest FIRE-2 zoom cosmological simulation of Milky Way-size halos. In addition, we utilize the newest version of the GALPROP cosmic-ray propagation framework to simulate the Galactic diffuse gamma-ray emission in the GC. We have found that p-wave (d-wave) DM particles with a mass of approximately 1.7 TeV and annihilating into the $W^+$$W^-$ channel exhibit a velocity-weighted annihilation cross-section upper limit of 4.6$\times$ 10$^{-22}$ cm$^3$s$^{-1}$ (9.2$\times$10$^{-17}$ cm$^3$s$^{- 1}$) at a 95\% confidence level. This is about 460 (2$\times$ 10$^{6}$) times greater than the thermal relic cross-section for p-wave (d-wave) DM models. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 16 pages, 9 figures, 1 table, 3 Appendices. Accepted for publication in Phys. Rev. D

Report number: Vol. 108, Iss. 8-15 October 2023

Journal ref: Phys. Rev. D 129, 108, 083027, 2023

Prospects for $γ$-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, K. Abe, S. Abe, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, E. O. Angüner, L. A. Antonelli, C. Aramo, M. Araya, C. Arcaro, L. Arrabito, K. Asano, Y. Ascasíbar, J. Aschersleben , et al. (542 additional authors not shown)

Abstract: Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster med… ▽ More Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $α_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $α_{\rm CRp}$ down to about $Δα_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $τ_χ>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 93 pages (including author list, appendix and references), 143 figures. Submitted to JCAP

The Giant Radio Array for Neutrino Detection (GRAND) Collaboration -- Contributions to the 38th International Cosmic Ray Conference (ICRC 2023)

Authors: GRAND Collaboration, Rafael Alves Batista, Aurélien Benoit-Lévy, Teresa Bister, Mauricio Bustamante, Yiren Chen, LingMei Cheng, Simon Chiche, Jean-Marc Colley, Pablo Correa, Nicoleta Cucu Laurenciu, Zigao Dai, Beatriz de Errico, Sijbrand de Jong, João R. T. de Mello Neto, Krijn D. de Vries, Peter B. Denton, Valentin Deocoene, Kaikai Duan, Bohao Duan, Ralph Engel, Yizhong Fan, Arsène Ferrière, QuanBu Gou, Junhua Gu , et al. (74 additional authors not shown)

Abstract: The Giant Radio Array for Neutrino Detection (GRAND) is an envisioned observatory of ultra-high-energy particles of cosmic origin, with energies in excess of 100 PeV. GRAND uses large surface arrays of autonomous radio-detection units to look for the radio emission from extensive air showers that are triggered by the interaction of ultra-high-energy cosmic rays, gamma rays, and neutrinos in the at… ▽ More The Giant Radio Array for Neutrino Detection (GRAND) is an envisioned observatory of ultra-high-energy particles of cosmic origin, with energies in excess of 100 PeV. GRAND uses large surface arrays of autonomous radio-detection units to look for the radio emission from extensive air showers that are triggered by the interaction of ultra-high-energy cosmic rays, gamma rays, and neutrinos in the atmosphere or underground. In particular, for ultra-high-energy neutrinos, the future final phase of GRAND aims to be sensitive enough to discover them in spite of their plausibly tiny flux. Presently, three prototype GRAND radio arrays are in operation: GRANDProto300, in China, GRAND@Auger, in Argentina, and GRAND@Nancay, in France. Their goals are to field-test the design of the radio-detection units, understand the radio background to which they are exposed, and develop tools for diagnostic, data gathering, and data analysis. This list of contributions to the 38th International Cosmic Ray Conference (ICRC 2023) presents an overview of GRAND, in its present and future incarnations, and a look at the first data collected by GRANDProto13, the first phase of GRANDProto300. △ Less

Submitted 27 July, 2023; originally announced August 2023.

Comments: Note: To access the list of contributions, please follow the "HTML" link that can be found on the arXiv page

Sensitivity of the Cherenkov Telescope Array to spectral signatures of hadronic PeVatrons with application to Galactic Supernova Remnants

Authors: The Cherenkov Telescope Array Consortium, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Aloisio, N. Álvarez Crespo, R. Alves Batista, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, C. Aramo, C. Arcaro, T. Armstrong, K. Asano, Y. Ascasibar, J. Aschersleben, M. Backes, A. Baktash, C. Balazs, M. Balbo , et al. (334 additional authors not shown)

Abstract: The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The pote… ▽ More The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The potential to search for hadronic PeVatrons with the Cherenkov Telescope Array (CTA) is assessed. The focus is on the usage of very high energy $γ$-ray spectral signatures for the identification of PeVatrons. Assuming that SNRs can accelerate CRs up to knee energies, the number of Galactic SNRs which can be identified as PeVatrons with CTA is estimated within a model for the evolution of SNRs. Additionally, the potential of a follow-up observation strategy under moonlight conditions for PeVatron searches is investigated. Statistical methods for the identification of PeVatrons are introduced, and realistic Monte--Carlo simulations of the response of the CTA observatory to the emission spectra from hadronic PeVatrons are performed. Based on simulations of a simplified model for the evolution for SNRs, the detection of a $γ$-ray signal from in average 9 Galactic PeVatron SNRs is expected to result from the scan of the Galactic plane with CTA after 10 hours of exposure. CTA is also shown to have excellent potential to confirm these sources as PeVatrons in deep observations with $\mathcal{O}(100)$ hours of exposure per source. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 34 pages, 16 figures, Accepted for publication in Astroparticle Physics

Prospective Dark Matter Annihilation Signals From the Sagittarius Dwarf Spheroidal

Authors: Thomas A. A. Venville, Alan R. Duffy, Roland M. Crocker, Oscar Macias, Thor Tepper-García

Abstract: The Sagittarius Dwarf Spheroidal galaxy (Sgr) is investigated as a target for DM annihilation searches utilising J-factor distributions calculated directly from a high-resolution hydrodynamic simulation of the infall and tidal disruption of Sgr around the Milky Way. In contrast to past studies, the simulation incorporates DM, stellar and gaseous components for both the Milky Way and the Sgr progen… ▽ More The Sagittarius Dwarf Spheroidal galaxy (Sgr) is investigated as a target for DM annihilation searches utilising J-factor distributions calculated directly from a high-resolution hydrodynamic simulation of the infall and tidal disruption of Sgr around the Milky Way. In contrast to past studies, the simulation incorporates DM, stellar and gaseous components for both the Milky Way and the Sgr progenitor galaxy. The simulated distributions account for significant tidal disruption affecting the DM density profile. Our estimate of the J-factor value for Sgr, $J_{\text{Sgr}}=1.48\times 10^{10}$ M$_\odot^2$ kpc$^{-5}$ ($6.46\times10^{16}\ \text{GeV}\ \text{cm}^{-5}$), is significantly lower than found in prior studies. This value, while formally a lower limit, is likely close to the true J-factor value for Sgr. It implies a DM cross-section incompatibly large in comparison with existing constraints would be required to attribute recently observed $γ$-ray emission from Sgr to DM annihilation. We also calculate a J-factor value using a NFW profile fitted to the simulated DM density distribution to facilitate comparison with past studies. This NFW J-factor value supports the conclusion that most past studies have overestimated the dark matter density of Sgr on small scales. This, together with the fact that the Sgr has recently been shown to emit $γ$-rays of astrophysical origin, complicate the use of Sgr in indirect DM detection searches. △ Less

Submitted 13 November, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: 15 pages, 15 figures. Accepted for publication in MNRAS. Accepted manuscript uploaded

Gamma-ray emission from the Sagittarius Dwarf Spheroidal galaxy due to millisecond pulsars

Authors: Roland M. Crocker, Oscar Macias, Dougal Mackey, Mark R. Krumholz, Shin'ichiro Ando, Shunsaku Horiuchi, Matthew G. Baring, Chris Gordon, Thomas Venville, Alan R. Duffy, Rui-Zhi Yang, Felix Aharonian, J. A. Hinton, Deheng Song, Ashley J. Ruiter, Miroslav D. Filipović

Abstract: The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show vi… ▽ More The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show via a spatial template analysis that much of the gamma-ray emission associated to the brightest region of substructure -- the so-called cocoon -- is likely due to the Sagittarius dwarf spheroidal (Sgr dSph) galaxy. This large Milky Way satellite is viewed through the Fermi Bubbles from the position of the Solar System. As a tidally and ram-pressure stripped remnant, the Sgr dSph has no on-going star formation, but we nevertheless demonstrate that the dwarf's millisecond pulsar (MSP) population can plausibly supply the gamma-ray signal that our analysis associates to its stellar template. The measured spectrum is naturally explained by inverse Compton scattering of cosmic microwave background photons by high-energy electron-positron pairs injected by MSPs belonging to the Sgr dSph, combined with these objects' magnetospheric emission. This finding plausibly suggests that MSPs produce significant gamma-ray emission amongst old stellar populations, potentially confounding indirect dark matter searches in regions such as the Galactic Centre, the Andromeda galaxy, and other massive Milky Way dwarf spheroidals. △ Less

Submitted 8 September, 2022; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: Updated to match version accepted for publication in Nature Astronomy (2022). 14 pages main text, 3 main figures, 7 extended data figures. For published version of the paper, see https://rdcu.be/cUZ6X

Journal ref: Nature Astronomy (2022)

The Andromeda Gamma-Ray Excess: Background Systematics of the Millisecond Pulsars and Dark Matter Interpretations

Authors: Fabian Zimmer, Oscar Macias, Shin'ichiro Ando, Roland M. Crocker, Shunsaku Horiuchi

Abstract: Since the discovery of an excess in gamma rays in the direction of M31, its cause has been unclear. Published interpretations focus on a dark matter or stellar related origin. Studies of a similar excess in the Milky Way center motivate a correlation of the spatial morphology of the signal with the distribution of stellar mass in M31. However, a robust determination of the best theory for the obse… ▽ More Since the discovery of an excess in gamma rays in the direction of M31, its cause has been unclear. Published interpretations focus on a dark matter or stellar related origin. Studies of a similar excess in the Milky Way center motivate a correlation of the spatial morphology of the signal with the distribution of stellar mass in M31. However, a robust determination of the best theory for the observed excess emission is very challenging due to large uncertainties in the astrophysical gamma-ray foreground model. Here we perform a spectro-morphological analysis of the M31 gamma-ray excess using state-of-the-art templates for the distribution of stellar mass in M31 and novel astrophysical foreground models for its sky region. We construct maps for the old stellar populations of M31 based on observational data from the PAndAS survey and carefully remove the foreground stars. We also produce improved astrophysical foreground models by using novel image inpainting techniques based on machine learning methods. We find that our stellar maps, taken as a proxy for the location of a putative population of millisecond pulsars in the bulge of M31, reach a statistical significance of $5.4σ$, making them as strongly favoured as the simple phenomenological models usually considered in the literature, e.g., a disk-like template with uniform brightness. Our detection of the stellar templates is robust to generous variations of the astrophysical foreground model. Once the stellar templates are included in the astrophysical model, we show that the dark matter annihilation interpretation of the signal is unwarranted. Using the results of a binary population synthesis model we demonstrate that a population of about one million unresolved MSPs could naturally explain the observed gamma-ray luminosity per stellar mass, energy spectrum, and stellar bulge-to-disk flux ratio. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: 15 pages, 11 figures, comments are welcome

Assessing the Impact of Hydrogen Absorption on the Characteristics of the Galactic Center Excess

Authors: Martin Pohl, Oscar Macias, Phaedra Coleman, Chris Gordon

Abstract: We present a new reconstruction of the distribution of atomic hydrogen in the inner Galaxy that is based on explicit radiation-transport modelling of line and continuum emission and a gas-flow model in the barred Galaxy that provides distance resolution for lines of sight toward the Galactic Center. The main benefits of the new gas model are, a), the ability to reproduce the negative line signals… ▽ More We present a new reconstruction of the distribution of atomic hydrogen in the inner Galaxy that is based on explicit radiation-transport modelling of line and continuum emission and a gas-flow model in the barred Galaxy that provides distance resolution for lines of sight toward the Galactic Center. The main benefits of the new gas model are, a), the ability to reproduce the negative line signals seen with the H$I$4PI survey and, b), the accounting for gas that primarily manifests itself through absorption. We apply the new model of Galactic atomic hydrogen to an analysis of the diffuse gamma-ray emission from the inner Galaxy, for which an excess at a few GeV was reported that may be related to dark matter. We find with high significance an improved fit to the diffuse gamma-ray emission observed with the Fermi-LAT, if our new H$I$ model is used to estimate the cosmic-ray induced diffuse gamma-ray emission. The fit still requires a nuclear bulge at high significance. Once this is included there is no evidence for a dark-matter signal, be it cuspy or cored. But an additional so-called boxy bulge is still favoured by the data. This finding is robust under the variation of various parameters, for example the excitation temperature of atomic hydrogen, and a number of tests for systematic issues. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: ApJ in press, 19 pages

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

Dark matter constraints with stacked gamma rays scales with the number of galaxies

Authors: Daiki Hashimoto, Atsushi J. Nishizawa, Masahiro Takada, Oscar Macias

Abstract: Low-surface-brightness galaxies (LSBGs) are interesting targets for searches of dark matter emission due to their low baryonic content. However, predicting their expected dark matter emissivities is difficult because of observational challenges in their distance measurements. Here we present a stacking method that makes use of catalogs of LSBGs and maps of unresolved gamma-ray emission measured by… ▽ More Low-surface-brightness galaxies (LSBGs) are interesting targets for searches of dark matter emission due to their low baryonic content. However, predicting their expected dark matter emissivities is difficult because of observational challenges in their distance measurements. Here we present a stacking method that makes use of catalogs of LSBGs and maps of unresolved gamma-ray emission measured by the Fermi Gamma-Ray Space Telescope. We show that, for relatively large number of LSBGs, individual distance measurements to the LSBGs are not necessary, instead the overall distance distribution of the population is sufficient in order to impose dark matter constraints. Further, we demonstrate that the effect of the covariance between two galaxies located closely -- at an angular distance comparable to the size of the Fermi point spread function -- is negligibly small. As a case in point, we apply our pipeline to a sample of 800 faint LSBGs discovered by Hyper Suprime-Cam and find that, the 95 per cent confidence level upper limits on the dark matter annihilation cross-section scales with inverse of the number LSBGs. In light of this linear dependence with the number of objects, we argue this methodology could provide extremely powerful limits if it is applied to the more than 105 LSBGs readily available with the Legacy Survey of Space and Time. △ Less

Submitted 29 April, 2022; v1 submitted 18 September, 2021; originally announced September 2021.

Comments: 12 pages, 7 figures, submitted to PRD

Measuring the smearing of the Galactic 511 keV signal: positron propagation or supernova kicks?

Authors: Thomas Siegert, Roland M. Crocker, Oscar Macias, Fiona H. Panther, Francesca Calore, Deheng Song, Shunsaku Horiuchi

Abstract: We use 15 years of $γ$-ray data from INTEGRAL/SPI in a refined investigation of the morphology of the Galactic bulge positron annihilation signal. Our spatial analysis confirms that the signal traces the old stellar population in the bulge and reveals for the first time that it traces the boxy bulge and nuclear stellar bulge. Using a 3D smoothing kernel, we find that the signal is smeared out over… ▽ More We use 15 years of $γ$-ray data from INTEGRAL/SPI in a refined investigation of the morphology of the Galactic bulge positron annihilation signal. Our spatial analysis confirms that the signal traces the old stellar population in the bulge and reveals for the first time that it traces the boxy bulge and nuclear stellar bulge. Using a 3D smoothing kernel, we find that the signal is smeared out over a characteristic length scale of $150 \pm 50\,$pc, suggesting either annihilation in situ at astrophysical sources kicked at formation or positron propagation away from sources. The former is disfavoured by its requiring kick velocities different between the Galactic nucleus ($\gtrsim 50\,\mathrm{km\,s^{-1}}$) and wider bulge ($\lesssim 15\,\mathrm{km\,s^{-1}}$) source. Positron propagation prior to annihilation can explain the overall phenomenology of the 511 keV signal for positrons injection energies $\lesssim 1.4\,$MeV, suggesting a nucleosynthesis origin. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 5 pages, 4 figures, 2 tables; submitted with MNRAS Letters

Millisecond Pulsars from Accretion Induced Collapse as the Origin of the Galactic Centre Gamma-ray Excess Signal

Authors: Anuj Gautam, Roland M. Crocker, Lilia Ferrario, Ashley J. Ruiter, Harrison Ploeg, Chris Gordon, Oscar Macias

Abstract: Gamma-ray data from the Fermi-Large Area Telescope reveal an unexplained, apparently diffuse, signal from the Galactic bulge that peaks near 2 GeV with an approximately spherical intensity profile $\propto r^{-2.4}$ that extends to angular radial scales of at least 10 degrees, possibly to 20 degrees. The origin of this "Galactic Centre Excess" (GCE) has been debated with proposed sources prominent… ▽ More Gamma-ray data from the Fermi-Large Area Telescope reveal an unexplained, apparently diffuse, signal from the Galactic bulge that peaks near 2 GeV with an approximately spherical intensity profile $\propto r^{-2.4}$ that extends to angular radial scales of at least 10 degrees, possibly to 20 degrees. The origin of this "Galactic Centre Excess" (GCE) has been debated with proposed sources prominently including self-annihilating dark matter and a hitherto undetected population of millisecond pulsars (MSPs). However, the conventional channel for the generation of MSPs has been found to predict too many low mass X-ray binary (LMXB) systems and, because of the expected large natal kicks, may not accommodate the close spatial correspondence between the GCE signal and stars in the bulge. Here we report a binary population synthesis forward model that demonstrates that an MSP population arising from the accretion induced collapse (AIC) of O-Ne white dwarfs in Galactic bulge binaries can naturally reproduce the morphology, spectral shape, and intensity of the GCE signal while also obeying LMXB constraints. Synchrotron emission from MSP-launched cosmic ray electrons and positrons may simultaneously explain the mysterious, microwave "haze" from the inner Galaxy △ Less

Submitted 2 May, 2022; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: 4 pages main text, 3 main figures, 6 pages Methods, 33 pages supplementary material. Updated to match version accepted for publication in Nature Astronomy (2022). For published version of the paper, see https://rdcu.be/cMjKY

Journal ref: Nature Astronomy (2022)

Cherenkov Telescope Array sensitivity to the putative millisecond pulsar population responsible for the Galactic center excess

Authors: Oscar Macias, Harm van Leijen, Deheng Song, Shin'ichiro Ando, Shunsaku Horiuchi, Roland M. Crocker

Abstract: The leading explanation of the $\textit{Fermi}$ Galactic center $γ$-ray excess is the extended emission from a unresolved population of millisecond pulsars (MSPs) in the Galactic bulge. Such a population would, along with the prompt $γ$ rays, also inject large quantities of electrons/positrons ($e^\pm$) into the interstellar medium. These $e^\pm$ could potentially inverse-Compton (IC) scatter ambi… ▽ More The leading explanation of the $\textit{Fermi}$ Galactic center $γ$-ray excess is the extended emission from a unresolved population of millisecond pulsars (MSPs) in the Galactic bulge. Such a population would, along with the prompt $γ$ rays, also inject large quantities of electrons/positrons ($e^\pm$) into the interstellar medium. These $e^\pm$ could potentially inverse-Compton (IC) scatter ambient photons into $γ$ rays that fall within the sensitivity range of the upcoming Cherenkov Telescope Array (CTA). In this article, we examine the detection potential of CTA to this signature by making a realistic estimation of the systematic uncertainties on the Galactic diffuse emission model at TeV-scale $γ$-ray energies. We forecast that, in the event that $e^\pm$ injection spectra are harder than $E^{-2}$, CTA has the potential to robustly discover the IC signature of a putative Galactic bulge MSP population sufficient to explain the GCE for $e^\pm$ injection efficiencies in the range $\approx 2.9-74.1\%$, or higher, depending on the level of mismodeling of the Galactic diffuse emission components. On the other hand, for spectra softer than $E^{-2.5}$, a reliable CTA detection would require an unphysically large $e^\pm$ injection efficiency of $\gtrsim 158\%$. However, even this pessimistic conclusion may be avoided in the plausible event that MSP observational and/or modeling uncertainties can be reduced. We further find that, in the event that an IC signal were detected, CTA can successfully discriminate between an MSP and a dark matter origin for the radiating $e^\pm$. △ Less

Submitted 26 June, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 21+6 pages, 12+5 figures. Updated to match version soon to be published by MNRAS

Journal ref: MNRAS 506, 1741-1760 (2021)

Evidence for a high-energy tail in the gamma-ray spectra of globular clusters

Authors: Deheng Song, Oscar Macias, Shunsaku Horiuchi, Roland M. Crocker, David M. Nataf

Abstract: Millisecond pulsars are very likely the main source of gamma-ray emission from globular clusters. However, the relative contributions of two separate emission processes--curvature radiation from millisecond pulsar magnetospheres vs. inverse Compton emission from relativistic pairs launched into the globular cluster environment by millisecond pulsars--have long been unclear. To address this, we sea… ▽ More Millisecond pulsars are very likely the main source of gamma-ray emission from globular clusters. However, the relative contributions of two separate emission processes--curvature radiation from millisecond pulsar magnetospheres vs. inverse Compton emission from relativistic pairs launched into the globular cluster environment by millisecond pulsars--have long been unclear. To address this, we search for evidence of inverse Compton emission in 8-year $\textit{Fermi}$-LAT data from the directions of 157 Milky Way globular clusters. We find a mildly statistically significant (3.8$σ$) correlation between the measured globular cluster gamma-ray luminosities and their photon field energy densities. However, this may also be explained by a hidden correlation between the photon field densities and the stellar encounter rates of globular clusters. Analysed $\textit{in toto}$, we demonstrate that the gamma-ray emission of globular clusters can be resolved spectrally into two components: i) an exponentially cut-off power law and ii) a pure power law. The latter component--which we uncover at a significance of 8.2$σ$--has a power index of 2.79 $\pm$ 0.25. It is most naturally interpreted as inverse Compton emission by cosmic-ray electrons and positrons injected by millisecond pulsars. We find the luminosity of this power-law component is comparable to, or slightly smaller than, the luminosity of the curved component, suggesting the fraction of millisecond pulsar spin-down luminosity into relativistic leptons is similar to the fraction of the spin-down luminosity into prompt magnetospheric radiation. △ Less

Submitted 26 September, 2021; v1 submitted 29 January, 2021; originally announced February 2021.

Comments: 16+9 pages, 13+4 figures, 4+1 tables. match the published version in MNRAS

Journal ref: MNRAS 507, 5161-5176 (2021)

Comparing the Galactic Bulge and Galactic Disk Millisecond Pulsars

Authors: Harrison Ploeg, Chris Gordon, Roland Crocker, Oscar Macias

Abstract: The Galactic Center Excess (GCE) is an extended gamma-ray source in the central region of the Galaxy found in Fermi Large Area Telescope (Fermi-LAT) data. One of the leading explanations for the GCE is an unresolved population of millisecond pulsars (MSPs) in the Galactic bulge. Due to differing star formation histories it is expected that the MSPs in the Galactic bulge are older and therefore dim… ▽ More The Galactic Center Excess (GCE) is an extended gamma-ray source in the central region of the Galaxy found in Fermi Large Area Telescope (Fermi-LAT) data. One of the leading explanations for the GCE is an unresolved population of millisecond pulsars (MSPs) in the Galactic bulge. Due to differing star formation histories it is expected that the MSPs in the Galactic bulge are older and therefore dimmer than those in the Galactic disk. Additionally, correlations between the spectral parameters of the MSPs and the spin-down rate of the corresponding neutron stars have been observed. This implies that the bulge MSPs may be spectrally different from the disk MSPs. We perform detailed modelling of the MSPs from formation until observation. Although we confirm the correlations, we do not find they are sufficiently large to significantly differentiate the spectra of the bulge MSPs and disk MSPs when the uncertainties are accounted for. Our results demonstrate that the population of MSPs that can explain the gamma-ray signal from the resolved MSPs in the Galactic disk and the unresolved MSPs in the boxy bulge and nuclear bulge can consistently be described as arising from a common evolutionary trajectory for some subset of astrophysical sources common to all these different environments. We do not require that there is anything unusual about inner Galaxy MSPs to explain the GCE. Additionally, we use a more accurate geometry for the distribution of bulge MSPs and incorporate dispersion measure estimates of the MSPs' distances. We find that the elongated boxy bulge morphology means that some the bulge MSPs are closer to us and so easier to resolve. We identify three resolved MSPs that have significant probabilities of belonging to the bulge population. △ Less

Submitted 14 June, 2021; v1 submitted 25 August, 2020; originally announced August 2020.

Comments: 37 pages, 20 figures, V3: Minor grammatical corrections/clarifications. Reflects version to appear in JCAP. V4: Minor typos corrected

Journal ref: JCAP12(2020)035

Prospects for Heavy WIMP Dark Matter with CTA: the Wino and Higgsino

Authors: Lucia Rinchiuso, Oscar Macias, Emmanuel Moulin, Nicholas L. Rodd, Tracy R. Slatyer

Abstract: TeV-scale particles that couple to the standard model through the weak force represent a compelling class of dark matter candidates. The search for such Weakly Interacting Massive Particles (WIMPs) has already spanned multiple decades, and whilst it has yet to provide any definitive evidence for their existence, viable parameter space remains. In this paper, we show that the upcoming Cherenkov Tel… ▽ More TeV-scale particles that couple to the standard model through the weak force represent a compelling class of dark matter candidates. The search for such Weakly Interacting Massive Particles (WIMPs) has already spanned multiple decades, and whilst it has yet to provide any definitive evidence for their existence, viable parameter space remains. In this paper, we show that the upcoming Cherenkov Telescope Array (CTA) has significant sensitivity to uncharted parameter space at the TeV mass scale. To do so, we focus on two prototypical dark matter candidates, the Wino and Higgsino. Sensitivity forecasts for both models are performed including the irreducible background from misidentified cosmic rays, as well as a range of estimates for the Galactic emissions at TeV energies. For each candidate, we find substantial expected improvements over existing bounds from current imaging atmospheric Cherenkov telescopes. In detail, for the Wino we find a sensitivity improvement of roughly an order of magnitude in $\langle σv \rangle$, whereas for the Higgsino we demonstrate that CTA has the potential to become the first experiment that has sensitivity to the thermal candidate. Taken together, these enhanced sensitivities demonstrates the discovery potential for dark matter at CTA in the 1-100 TeV mass range. △ Less

Submitted 15 December, 2020; v1 submitted 3 August, 2020; originally announced August 2020.

Comments: 20 pages, 8 figures, Accepted for publication in Phys. Rev. D

Report number: MIT-CTP 5120, IRFU-20-13

Journal ref: Phys. Rev. D 103, 023011 (2021)

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

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

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

Submitted 31 March, 2020; originally announced April 2020.

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

Report number: LCTP-20-04

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

Strong constraints on thermal relic dark matter from Fermi-LAT observations of the Galactic Center

Authors: Kevork N. Abazajian, Shunsaku Horiuchi, Manoj Kaplinghat, Ryan E. Keeley, Oscar Macias

Abstract: The extended excess toward the Galactic Center (GC) in gamma rays inferred from Fermi-LAT observations has been interpreted as being due to dark matter (DM) annihilation. Here, we perform new likelihood analyses of the GC and show that, when including templates for the stellar galactic and nuclear bulges, the GC shows no significant detection of a DM annihilation template, even after generous vari… ▽ More The extended excess toward the Galactic Center (GC) in gamma rays inferred from Fermi-LAT observations has been interpreted as being due to dark matter (DM) annihilation. Here, we perform new likelihood analyses of the GC and show that, when including templates for the stellar galactic and nuclear bulges, the GC shows no significant detection of a DM annihilation template, even after generous variations in the Galactic diffuse emission models and a wide range of DM halo profiles. We include Galactic diffuse emission models with combinations of three-dimensional inverse Compton maps, variations of interstellar gas maps, and a central source of electrons. For the DM profile, we include both spherical and ellipsoidal DM morphologies and a range of radial profiles from steep cusps to kiloparsec-sized cores, motivated in part by hydrodynamical simulations. Our derived upper limits on the dark matter annihilation flux place strong constraints on DM properties. In the case of the pure $b$-quark annihilation channel, our limits on the annihilation cross section are more stringent than those from the Milky Way dwarfs up to DM masses of approximately TeV and rule out the thermal relic cross section up to approximately 300 GeV. Better understanding of the DM profile, as well as the Fermi-LAT data at its highest energies, would further improve the sensitivity to DM properties. △ Less

Submitted 25 January, 2021; v1 submitted 23 March, 2020; originally announced March 2020.

Comments: 23 pages, 14 figures. V2: more IC rings included, profile likelihood corrected, more tests added; main results and conclusions unchanged. V3: corrected typo found in Tab. II. Matches version published in PRD

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

Foreground Mismodeling and the Point Source Explanation of the Fermi Galactic Center Excess

Authors: Malte Buschmann, Nicholas L. Rodd, Benjamin R. Safdi, Laura J. Chang, Siddharth Mishra-Sharma, Mariangela Lisanti, Oscar Macias

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

Submitted 27 February, 2020; originally announced February 2020.

Comments: 27 pages, 18 figures, comments welcome

Report number: LCTP-20-02

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

Cross Correlation of the Extragalactic Gamma-ray Background with Thermal Sunyaev-Zel'dovich Effect in the Cosmic Microwave Background

Authors: Masato Shirasaki, Oscar Macias, Shin'ichiro Ando, Shunsaku Horiuchi, Naoki Yoshida

Abstract: Cosmic rays in galaxy clusters are unique probes of energetic processes operating with large-scale structures in the Universe. Precise measurements of cosmic rays in galaxy clusters are essential for improving our understanding of non-thermal components in the intracluster-medium (ICM) as well as the accuracy of cluster mass estimates in cosmological analyses. In this paper, we perform a cross-cor… ▽ More Cosmic rays in galaxy clusters are unique probes of energetic processes operating with large-scale structures in the Universe. Precise measurements of cosmic rays in galaxy clusters are essential for improving our understanding of non-thermal components in the intracluster-medium (ICM) as well as the accuracy of cluster mass estimates in cosmological analyses. In this paper, we perform a cross-correlation analysis with the extragalactic gamma-ray background and the thermal Sunyaev-Zel'dovich (tSZ) effect in the cosmic microwave background. The expected cross-correlation signal would contain rich information about the cosmic-ray-induced gamma-ray emission in the most massive galaxy clusters at $z\sim0.1-0.2$. We analyze the gamma-ray background map with 8 years of data taken by the Large Area Telescope onboard Fermi satellite and the publicly available tSZ map by Planck. We confirm that the measured cross-correlation is consistent with a null detection, and thus it enables us to put the tightest constraint on the acceleration efficiency of cosmic ray protons at shocks in and around galaxy clusters. We find the acceleration efficiency must be below 5\% with a $2σ$ confidence level when the hydrostatic mass bias of clusters is assumed to be 30\%, and our result is not significantly affected by the assumed value of the hydrostatic mass bias. Our constraint implies that the non-thermal cosmic-ray pressure in the ICM can introduce only a $\le 3\%$ level of the hydrostatic mass bias, highlighting that cosmic rays alone do not account for the mass bias inferred by the Planck analyses. Finally, we discuss future detectability prospects of cosmic-ray-induced gamma rays from the Perseus cluster for the Cherenkov Telescope Array. △ Less

Submitted 5 May, 2020; v1 submitted 26 November, 2019; originally announced November 2019.

Comments: 19 pages, 15 figures, 1 table. Accepted for publication in Phys. Rev. D

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

Maximum Entropy Estimation of the Galactic Bulge Morphology via the VVV Red Clump

Authors: Brenna Coleman, Dylan Paterson, Chris Gordon, Oscar Macias, Harrison Ploeg

Abstract: The abundance and narrow magnitude dispersion of Red Clump (RC) stars make them a popular candidate for mapping the morphology of the bulge region of the Milky Way. Using an estimate of the RC's intrinsic luminosity function, we extracted the three-dimensional density distribution of the RC from deep photometric catalogues of the VISTA Variables in the Via Lactea (VVV) survey. We used maximum en… ▽ More The abundance and narrow magnitude dispersion of Red Clump (RC) stars make them a popular candidate for mapping the morphology of the bulge region of the Milky Way. Using an estimate of the RC's intrinsic luminosity function, we extracted the three-dimensional density distribution of the RC from deep photometric catalogues of the VISTA Variables in the Via Lactea (VVV) survey. We used maximum entropy based deconvolution to extract the spatial distribution of the bulge from Ks-band star counts. We obtained our extrapolated non-parametric model of the bulge over the inner 40 by 40 degrees squared region of the Galactic centre. Our reconstruction also naturally matches onto a parametric fit to the bulge outside the VVV region and inpaints overcrowded and high extinction regions. We found a range of bulge properties consistent with other recent investigations based on the VVV data. In particular, we estimated the bulge mass to be in the range 13 to 17 billion solar masses, the X-component to be between 18% and 25% of the bulge mass, and the bulge angle with respect to the Sun-Galactic centre line to be between 18 and 32 degrees. Studies of the Fermi Large Area Telescope (LAT) gamma-ray Galactic centre excess suggests that the excess may be traced by Galactic bulge distributed sources. We applied our deconvolved density in a template fitting analysis of this Fermi-LAT GeV excess and found an improvement in the fit compared to previous parametric based templates. △ Less

Submitted 10 December, 2020; v1 submitted 12 November, 2019; originally announced November 2019.

Comments: 25 pages, 27 figures, minor typo corrected

Journal ref: 2020MNRAS.495.3350C

Probing heavy dark matter decays with multi-messenger astrophysical data

Authors: Koji Ishiwata, Oscar Macias, Shin'ichiro Ando, Makoto Arimoto

Abstract: We set conservative constraints on decaying dark matter particles with masses spanning a very wide range ($10^4-10^{16}$ GeV). For this we use multimessenger observations of cosmic-ray (CR) protons/antiprotons, electrons/positrons, neutrinos/antineutrinos and gamma rays. Focusing on decays into the $\bar{b}b$ channel, we simulate the spectra of dark matter yields by using the Dokshitzer-Gribov-Lip… ▽ More We set conservative constraints on decaying dark matter particles with masses spanning a very wide range ($10^4-10^{16}$ GeV). For this we use multimessenger observations of cosmic-ray (CR) protons/antiprotons, electrons/positrons, neutrinos/antineutrinos and gamma rays. Focusing on decays into the $\bar{b}b$ channel, we simulate the spectra of dark matter yields by using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equations and the $\texttt{Pythia}$ package. We then propagate the CRs of dark matter origin till Earth by using the state-of-the-art numerical frameworks $\texttt{CRPropa}$, $\texttt{GALPROP}$ and $\texttt{HelMod}$ for the solution of the CR transport equation in the extragalactic, Galactic region and the heliosphere, respectively. Conservative limits are obtained by requiring that the predicted dark matter spectra at Earth be less than the observed CR spectra. Overall, we exclude dark matter lifetimes of $10^{28}$ s or shorter for all the masses investigated in this work. The most stringent constraints reach $10^{30}$ s for very heavy dark matter particles with masses in the range $10^{11}-10^{14}$ GeV. △ Less

Submitted 3 January, 2020; v1 submitted 26 July, 2019; originally announced July 2019.

Comments: 22 pages, 10 figures, V2: Minor changes to match JCAP published version

Report number: KANAZAWA-19-06

Journal ref: JCAP01(2020)003

Constraining dark matter annihilation with HSC Low Surface Brightness Galaxies

Authors: Daiki Hashimoto, Oscar Macias, Atsushi J. Nishizawa, Kohei Hayashi, Masahiro Takada, Masato Shirasaki, Shin'ichiro Ando

Abstract: Searches for dark matter annihilation signals have been carried out in a number of target regions such as the Galactic Center and Milky Way dwarf spheroidal galaxies (dSphs), among a few others. Here we propose low surface brightness galaxies (LSBGs) asnovel targets for the indirect detection of dark matter emission. In particular, LSBGs are known to have very large dark matter contents and be les… ▽ More Searches for dark matter annihilation signals have been carried out in a number of target regions such as the Galactic Center and Milky Way dwarf spheroidal galaxies (dSphs), among a few others. Here we propose low surface brightness galaxies (LSBGs) asnovel targets for the indirect detection of dark matter emission. In particular, LSBGs are known to have very large dark matter contents and be less contaminated by extragalactic gamma-ray sources (e.g., blazars) compared to star forming galaxies. We report on an analysis that uses eight LSBGs (detected by Subaru Hyper Suprime-Cam survey data) with known redshifts to conduct a search for gamma-ray emission at the positions of these new objects in Fermi Large Area Telescope data. We found no excesses of gamma-ray emission and set constraints on the dark matter annihilation cross-section. We exclude (at the 95% C.L.) dark matter scenarios predicting a cross-section higher than 10^-23[cm^3/s] for dark matter particles of mass 10 GeV self-annihilating in the b_b channel. Although this constraint is weaker than the ones reported in recent studies using other targets, we note that in the near future, the number of detections of new LSBGs will increase by a few orders of magnitude. We forecast that with the use of the full catalog of soon-to-be-detected LSBGs the constraint will reach cross-section sensitivities of ~ 3*10^-25 [cm^3/s] for dark matter particles with masses less than 10 GeV. △ Less

Submitted 8 January, 2020; v1 submitted 16 June, 2019; originally announced June 2019.

Comments: 19 pages, 6 figures, accepted to JCAP

Discovery prospects of dwarf spheroidal galaxies for indirect dark matter searches

Authors: Shin'ichiro Ando, Bradley J. Kavanagh, Oscar Macias, Tiago Alves, Siebren Broersen, Stijn Delnoij, Thomas Goldman, Jim Groefsema, Jorinde Kleverlaan, Jaïr Lenssen, Toon Muskens, Liam X. Palma Visser, Ebo Peerbooms, Bram van der Linden, Sill Verberne

Abstract: We study the prospects for the Large Synoptic Survey Telescope (LSST) to find new dwarf spheroidal galaxies in the Milky Way. Adopting models of Milky-Way halo substructure and phenomenological prescriptions connecting subhalos and satellite galaxies, we obtain surface brightness distributions of $V$-band magnitude that lead us to predict that LSST will discover tens to hundreds of dwarf spheroida… ▽ More We study the prospects for the Large Synoptic Survey Telescope (LSST) to find new dwarf spheroidal galaxies in the Milky Way. Adopting models of Milky-Way halo substructure and phenomenological prescriptions connecting subhalos and satellite galaxies, we obtain surface brightness distributions of $V$-band magnitude that lead us to predict that LSST will discover tens to hundreds of dwarf spheroidal galaxies above its sensitivity. The soon-to-be-discovered dwarfs will be interesting targets for indirect searches of dark matter annihilation yields. We forecast the distribution function of gamma-ray emission from dark matter annihilation in these objects, and discuss the detectability of these signals at both Fermi Large Area Telescope (LAT) and Cherenkov Telescope Array (CTA). By combining information from the predicted dwarf galaxies, we obtain an expected sensitivity to the annihilation cross section $\langle σv \rangle$ of $10^{-26}$ cm$^3$ s$^{-1}$ (for dark matter particles of mass 10 GeV with Fermi-LAT) and $5\times 10^{-24}$ cm$^3$ s$^{-1}$ (for dark matter particles of mass 500 GeV with CTA). We find that the current uncertainties in the mass measurement of the Milky-Way halo are relatively minor compared with the Poisson errors associated to drawing the most promising dwarfs from the underlying flux distribution. △ Less

Submitted 14 October, 2019; v1 submitted 17 May, 2019; originally announced May 2019.

Comments: 13 pages (+ references), 6 figures. v2: minor corrections, matches version published in JCAP

Journal ref: JCAP 10 (2019) 040

Inverse Compton emission from millisecond pulsars in the Galactic bulge

Authors: Deheng Song, Oscar Macias, Shunsaku Horiuchi

Abstract: Analyses of Fermi Gamma-Ray Space Telescope data have revealed a source of excess diffuse gamma rays towards the Galactic center that extends up to roughly $\pm$20 degrees in latitude. The leading theory postulates that this GeV excess is the aggregate emission from a large number of faint millisecond pulsars (MSPs). The electrons and positrons ($e^\pm$) injected by this population could produce d… ▽ More Analyses of Fermi Gamma-Ray Space Telescope data have revealed a source of excess diffuse gamma rays towards the Galactic center that extends up to roughly $\pm$20 degrees in latitude. The leading theory postulates that this GeV excess is the aggregate emission from a large number of faint millisecond pulsars (MSPs). The electrons and positrons ($e^\pm$) injected by this population could produce detectable inverse-Compton (IC) emissions by up-scattering ambient photons to gamma-ray energies. In this work, we calculate such IC emissions using GALPROP. A triaxial three-dimensional model of the bulge stars obtained from a fit to infrared data is used as a tracer of the putative MSP population. This model is compared against one in which the MSPs are spatially distributed as a Navarro-Frenk-White squared profile. We show that the resulting spectra for both models are indistinguishable, but that their spatial morphologies have salient recognizable features. The IC component above $\sim$TeV energies carries information on the spatial morphology of the injected $e^\pm$. Such differences could potentially be used by future high-energy gamma-ray detectors such as the Cherenkov Telescope Array to provide a viable multiwavelength handle for the MSP origin of the GeV excess. △ Less

Submitted 18 July, 2019; v1 submitted 21 January, 2019; originally announced January 2019.

Comments: 13 pages, 8 figures, 3 tables. Match the version published in PRD

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

Strong Evidence that the Galactic Bulge is Shining in Gamma Rays

Authors: Oscar Macias, Shunsaku Horiuchi, Manoj Kaplinghat, Chris Gordon, Roland M. Crocker, David M. Nataf

Abstract: There is growing evidence that the Galactic Center Excess identified in the $\textit{Fermi}$-LAT gamma-ray data arises from a population of faint astrophysical sources. We provide compelling supporting evidence by showing that the morphology of the excess traces the stellar over-density of the Galactic bulge. By adopting a template of the bulge stars obtained from a triaxial 3D fit to the diffuse… ▽ More There is growing evidence that the Galactic Center Excess identified in the $\textit{Fermi}$-LAT gamma-ray data arises from a population of faint astrophysical sources. We provide compelling supporting evidence by showing that the morphology of the excess traces the stellar over-density of the Galactic bulge. By adopting a template of the bulge stars obtained from a triaxial 3D fit to the diffuse near-infrared emission, we show that it is detected at high significance. The significance deteriorates when either the position or the orientation of the template is artificially shifted, supporting the correlation of the gamma-ray data with the Galactic bulge. In deriving these results, we have used more sophisticated templates at low-latitudes for the $\textit{Fermi}$ bubbles compared to previous work and the three-dimensional Inverse Compton (IC) maps recently released by the ${\tt GALPROP}$ team. Our results provide strong constraints on Millisecond Pulsar (MSP) formation scenarios proposed to explain the excess. We find that an $\textit{admixture formation}$ scenario, in which some of the relevant binaries are $\textit{primordial}$ and the rest are formed $\textit{dynamically}$, is preferred over a primordial-only formation scenario at $7.6σ$ confidence level. Our detailed morphological analysis also disfavors models of the disrupted globular clusters scenario that predict a spherically symmetric distribution of MSPs in the Galactic bulge. For the first time, we report evidence of a high energy tail in the nuclear bulge spectrum that could be the result of IC emission from electrons and positrons injected by a population of MSPs and star formation activity from the same site. △ Less

Submitted 20 September, 2019; v1 submitted 12 January, 2019; originally announced January 2019.

Comments: 21 pages, 13 figures, V2: Minor changes to match submitted version, V3: matches JCAP published version

Journal ref: JCAP09(2019)042

Measurement of redshift dependent cross correlation of HSC clusters and Fermi $γ$ rays

Authors: Daiki Hashimoto, Atsushi J. Nishizawa, Masato Shirasaki, Oscar Macias, Shunsaku Horiuchi, Hiroyuki Tashiro, Masamune Oguri

Abstract: The cross-correlation study of the unresolved $γ$-ray background (UGRB) with galaxy clusters has a potential to reveal the nature of the UGRB. In this paper, we perform a cross-correlation analysis between $γ$-ray data by the Fermi Large Area Telescope (Fermi-LAT) and a galaxy cluster catalogue from the Subaru Hyper Suprime-Cam (HSC) survey. The Subaru HSC cluster catalogue provides a wide and hom… ▽ More The cross-correlation study of the unresolved $γ$-ray background (UGRB) with galaxy clusters has a potential to reveal the nature of the UGRB. In this paper, we perform a cross-correlation analysis between $γ$-ray data by the Fermi Large Area Telescope (Fermi-LAT) and a galaxy cluster catalogue from the Subaru Hyper Suprime-Cam (HSC) survey. The Subaru HSC cluster catalogue provides a wide and homogeneous large-scale structure distribution out to the high redshift at $z=1.1$, which has not been accessible in previous cross-correlation studies. We conduct the cross-correlation analysis not only for clusters in the all redshift range ($0.1 < z < 1.1$) of the survey, but also for subsamples of clusters divided into redshift bins, the low redshift bin ($0.1 < z < 0.6$) and the high redshift bin ($0.6 < z < 1.1$), to utilize the wide redshift coverage of the cluster catalogue. We find the evidence of the cross-correlation signals with the significance of 2.0-2.3$σ$ for all redshift and low-redshift cluster samples. On the other hand, for high-redshift clusters, we find the signal with weaker significance level (1.6-1.9$σ$). We also compare the observed cross-correlation functions with predictions of a theoretical model in which the UGRB originates from $γ$-ray emitters such as blazars, star-forming galaxies and radio galaxies. We find that the detected signal is consistent with the model prediction. △ Less

Submitted 30 January, 2019; v1 submitted 21 May, 2018; originally announced May 2018.

Comments: 11 pages, 24 figures, accepted by MNRAS

The correlation of extragalactic $γ$-rays with cosmic matter density distributions from weak-gravitational lensing

Authors: Masato Shirasaki, Oscar Macias, Shunsaku Horiuchi, Naoki Yoshida, Chien-Hsiu Lee, Atsushi J. Nishizawa

Abstract: The extragalactic $γ$-ray background (EGB) arises from the accumulation of $γ$-ray emissions from resolved and unresolved extragalactic sources as well as diffuse processes. It is important to study the statistical properties of the EGB in the context of cosmological structure formation. Known astrophysical $γ$-ray sources such as blazars, star-forming galaxies, and radio galaxies are expected to… ▽ More The extragalactic $γ$-ray background (EGB) arises from the accumulation of $γ$-ray emissions from resolved and unresolved extragalactic sources as well as diffuse processes. It is important to study the statistical properties of the EGB in the context of cosmological structure formation. Known astrophysical $γ$-ray sources such as blazars, star-forming galaxies, and radio galaxies are expected to trace the underlying cosmic matter density distribution. We explore the correlation of the EGB from Fermi-LAT data with the large-scale matter density distribution from the Subaru Hyper Suprime-Cam (HSC) SSP survey. We reconstruct an unbiased surface matter density distribution $κ$ at $z<1$ by applying weak-gravitational lensing analysis to the first-year HSC data. We then calculate the $γ- κ$ cross-correlation. Our measurements are consistent with a null detection, but a weak correlation is found at angular scales of 30-60 arcmin, especially when distant source galaxies at $z > 1$ are used for the lensing $κ$ reconstruction. The large-scale correlation suggests strong clustering of high-redshift $γ$-ray sources such as blazars. However, the inferred bias factor of $4-5$ is larger by about a factor of two than results from other clustering analyses. The final HSC data covering 1,400 squared degrees will play an essential role to determine accurately the blazar bias at $z > 0.5$. △ Less

Submitted 4 June, 2018; v1 submitted 27 February, 2018; originally announced February 2018.

Comments: 17 pages, 8 figures, 4 tables, accepted for publication in Physical Review D

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

Consistency Between the Luminosity Function of Resolved Millisecond Pulsars and the Galactic Center Excess

Authors: Harrison Ploeg, Chris Gordon, Roland Crocker, Oscar Macias

Abstract: Fermi Large Area Telescope data reveal an excess of GeV gamma rays from the direction of the Galactic Center and bulge. Several explanations have been proposed for this excess including an unresolved population of millisecond pulsars (MSPs) and self-annihilating dark matter. It has been claimed that a key discriminant for or against the MSP explanation can be extracted from the properties of the l… ▽ More Fermi Large Area Telescope data reveal an excess of GeV gamma rays from the direction of the Galactic Center and bulge. Several explanations have been proposed for this excess including an unresolved population of millisecond pulsars (MSPs) and self-annihilating dark matter. It has been claimed that a key discriminant for or against the MSP explanation can be extracted from the properties of the luminosity function describing this source population. Specifically, is the luminosity function of the putative MSPs in the Galactic Center consistent with that characterizing the resolved MSPs in the Galactic disk? To investigate this we have used a Bayesian Markov Chain Monte Carlo to evaluate the posterior distribution of the parameters of the MSP luminosity function describing both resolved MSPs and the Galactic Center excess. At variance with some other claims, our analysis reveals that, within current uncertainties, both data sets can be well fit with the same luminosity function. △ Less

Submitted 11 August, 2017; v1 submitted 2 May, 2017; originally announced May 2017.

Comments: 37 pages, 18 figures; V2: updated references and extended discussion slightly. V3: Additional material added, conclusions unchanged, matches version to appear in JCAP

Journal ref: JCAP08(2017)015

Galactic Bulge Preferred Over Dark Matter for the Galactic Center Gamma-Ray Excess

Authors: Oscar Macias, Chris Gordon, Roland M. Crocker, Brenna Coleman, Dylan Paterson, Shunsaku Horiuchi, Martin Pohl

Abstract: An anomalous gamma-ray excess emission has been found in Fermi Large Area Telescope data covering the centre of the Galaxy. Several theories have been proposed for this `Galactic Centre Excess'. They include self-annihilation of dark matter particles, an unresolved population of millisecond pulsars, an unresolved population of young pulsars, or a series of burst events. Here we report on a new a… ▽ More An anomalous gamma-ray excess emission has been found in Fermi Large Area Telescope data covering the centre of the Galaxy. Several theories have been proposed for this `Galactic Centre Excess'. They include self-annihilation of dark matter particles, an unresolved population of millisecond pulsars, an unresolved population of young pulsars, or a series of burst events. Here we report on a new analysis that exploits hydrodynamical modelling to register the position of interstellar gas associated with diffuse Galactic gamma-ray emission. We find evidence that the Galactic Centre Excess gamma rays are statistically better described by the stellar over-density in the Galactic bulge and the nuclear stellar bulge, rather than a spherical excess. Given its non-spherical nature, we argue that the Galactic Centre Excess is not a dark matter phenomenon but rather associated with the stellar population of the Galactic bulge and nuclear bulge. △ Less

Submitted 25 September, 2018; v1 submitted 20 November, 2016; originally announced November 2016.

Comments: 21 pages, 7 figures, V2: Minor modifications, main conclusions unchanged. V3: version submitted to Nature Astronomy. V4: matches Nature Astronomy published version. The hydrodynamical interstellar maps used in this work are available at https://github.com/chrisgordon1/galactic_bulge

Journal ref: Nature Astronomy (2018)

Supernova Physics at DUNE

Authors: Artur Ankowski, John Beacom, Omar Benhar, Sun Chen, John Cherry, Yanou Cui, Alexander Friedland, Ines Gil-Botella, Alireza Haghighat, Shunsaku Horiuchi, Patrick Huber, James Kneller, Ranjan Laha, Shirley Li, Jonathan Link, Alessandro Lovato, Oscar Macias, Camillo Mariani, Anthony Mezzacappa, Evan O'Connor, Erin O'Sullivan, Andre Rubbia, Kate Scholberg, Tatsu Takeuchi

Abstract: The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's lo… ▽ More The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors. △ Less

Submitted 28 August, 2016; originally announced August 2016.

Comments: Summary of workshop "Supernova Physics at DUNE" held at Virginia Tech; for more details, see, http://cnp.phys.vt.edu/SNatDUNE/

Cosmological constraints on dark matter annihilation and decay: Cross-correlation analysis of the extragalactic $γ$-ray background and cosmic shear

Authors: Masato Shirasaki, Oscar Macias, Shunsaku Horiuchi, Satoshi Shirai, Naoki Yoshida

Abstract: We derive constraints on dark matter (DM) annihilation cross section and decay lifetime from cross-correlation analyses of the data from Fermi-LAT and weak lensing surveys that cover a wide area of $\sim660$ squared degrees in total. We improve upon our previous analyses by using an updated extragalactic $γ$-ray background data reprocessed with the Fermi Pass 8 pipeline, and by using well-calibrat… ▽ More We derive constraints on dark matter (DM) annihilation cross section and decay lifetime from cross-correlation analyses of the data from Fermi-LAT and weak lensing surveys that cover a wide area of $\sim660$ squared degrees in total. We improve upon our previous analyses by using an updated extragalactic $γ$-ray background data reprocessed with the Fermi Pass 8 pipeline, and by using well-calibrated shape measurements of about twelve million galaxies in the Canada-France-Hawaii Lensing Survey (CFHTLenS) and Red-Cluster-Sequence Lensing Survey (RCSLenS). We generate a large set of full-sky mock catalogs from cosmological $N$-body simulations and use them to estimate statistical errors accurately. The measured cross correlation is consistent with null detection, which is then used to place strong cosmological constraints on annihilating and decaying DM. For leptophilic DM, the constraints are improved by a factor of $\sim100$ in the mass range of O(1) TeV when including contributions from secondary $γ$ rays due to the inverse-Compton upscattering of background photons. Annihilation cross-sections of $\langle σv \rangle \sim 10^{-23}\, {\rm cm}^3/{\rm s}$ are excluded for TeV-scale DM depending on channel. Lifetimes of $\sim 10^{25}$ sec are also excluded for the decaying TeV-scale DM. Finally, we apply this analysis to wino DM and exclude the wino mass around 200 GeV. These constraints will be further tightened, and all the interesting wino DM parameter region can be tested, by using data from future wide-field cosmology surveys. △ Less

Submitted 27 September, 2016; v1 submitted 7 July, 2016; originally announced July 2016.

Comments: 19 pages, 10 figures, 2 tables, accepted for publication in Phys. Rev. D

Report number: DESY 16-112

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

The $\textit{Fermi}$-LAT gamma-ray excess at the Galactic Center in the singlet-doublet fermion dark matter model

Authors: Shunsaku Horiuchi, Oscar Macias, Diego Restrepo, Andres Rivera, Hamish Silverwood, Oscar Zapata

Abstract: The singlet-doublet fermion dark matter model (SDFDM) provides a good DM candidate as well as the possibility of generating neutrino masses radiatively. The search and identification of DM requires the combined effort of both indirect and direct DM detection experiments in addition to the LHC. Remarkably, an excess of GeV gamma rays from the Galactic Center (GCE) has been measured with the \textit… ▽ More The singlet-doublet fermion dark matter model (SDFDM) provides a good DM candidate as well as the possibility of generating neutrino masses radiatively. The search and identification of DM requires the combined effort of both indirect and direct DM detection experiments in addition to the LHC. Remarkably, an excess of GeV gamma rays from the Galactic Center (GCE) has been measured with the \textit{Fermi} Large Area Telescope (LAT) which appears to be robust with respect to changes in the diffuse galactic background modeling. Although several astrophysical explanations have been proposed, DM remains a simple and well motivated alternative. In this work, we examine the sensitivities of dark matter searches in the SDFDM scenario using $\textit{Fermi}$-LAT, CTA, IceCube/DeepCore, LUX, PICO and LHC with an emphasis on exploring the regions of the parameter space that can account for the GCE. We find that DM particles present in this model with masses close to $\sim 99$ GeV and $\sim (173-190)$ GeV annihilating predominantly into the $W^+W^-$ channel and $t\bar{t}$ channel respectively, provide an acceptable fit to the GCE while being consistent with different current experimental bounds. We also find that much of the obtained parameter space can be ruled out by future direct search experiments like LZ and XENON-1T, in case of null results by these detectors. Interestingly, we show that the most recent data by LUX is starting to probe the best fit region in the SDFDM model. △ Less

Submitted 15 March, 2016; v1 submitted 15 February, 2016; originally announced February 2016.

Comments: 21 pages, 7 figures. v2: references added. Version accepted by JCAP

Journal ref: JCAP03 (2016) 048

The Spatial Morphology of the Secondary Emission in the Galactic Center Gamma-Ray Excess

Authors: Thomas Lacroix, Oscar Macias, Chris Gordon, Paolo Panci, Celine Boehm, Joseph Silk

Abstract: Excess GeV gamma rays from the Galactic Center (GC) have been measured with the Fermi Large Area Telescope (LAT). The presence of the GC excess (GCE) appears to be robust with respect to changes in the diffuse galactic background modeling. The three main proposals for the GCE are an unresolved population of millisecond pulsars (MSPs), outbursts of cosmic rays from the GC region, and self-annihilat… ▽ More Excess GeV gamma rays from the Galactic Center (GC) have been measured with the Fermi Large Area Telescope (LAT). The presence of the GC excess (GCE) appears to be robust with respect to changes in the diffuse galactic background modeling. The three main proposals for the GCE are an unresolved population of millisecond pulsars (MSPs), outbursts of cosmic rays from the GC region, and self-annihilating dark matter (DM). The injection of secondary electrons and positrons into the interstellar medium (ISM) by an unresolved population of MSPs or DM annihilations can lead to observable gamma-ray emission via inverse Compton scattering or bremsstrahlung. Here we investigate how to determine whether secondaries are important in a model for the GCE. We develop a method of testing model fit which accounts for the different spatial morphologies of the secondary emission. We examine several models which give secondary emission and illustrate a case where a broadband analysis is not sufficient to determine the need for secondary emission. △ Less

Submitted 20 April, 2016; v1 submitted 6 December, 2015; originally announced December 2015.

Comments: 12 pages, 3 figures. Fixed convergence problem in MSP case. Matches version accepted by PRD

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

Can Cosmic Rays Interacting With Molecular Clouds Explain The Galactic Center Gamma-Ray Excess?

Authors: Chris Gordon, Oscar Macias

Abstract: The Fermi Large Area Telescope data appear to have an excess of gamma rays from the inner 150 pc of the Galactic Center (GC). The main explanations proposed for this are: an unresolved population of millisecond pulsars (MSPs), dark matter (DM) annihilation, and cosmic rays interacting with molecular clouds at the GC. In this conference proceeding article we highlight some of the cosmic-ray results… ▽ More The Fermi Large Area Telescope data appear to have an excess of gamma rays from the inner 150 pc of the Galactic Center (GC). The main explanations proposed for this are: an unresolved population of millisecond pulsars (MSPs), dark matter (DM) annihilation, and cosmic rays interacting with molecular clouds at the GC. In this conference proceeding article we highlight some of the cosmic-ray results of [1]. The MSPs and DM explanations were modeled with spatial templates well fitted by the square of a generalized Navarro-Frenk-White (NFW) profile with inner slope $γ=1.2$. The cosmic-ray option was modeled with a 20-cm continuum emission Galactic Ridge template. A template based on the HESS residuals were shown to give similar results. The gamma-ray excess was found to be best fit by a combination of the generalized NFW squared template and a Galactic Ridge template. We also found the spectra of each template was not significantly affected in the combined fit and is consistent with previous single template fits. Due to the ridge like morphology of the molecular gas at the GC, the cosmic-ray option was not able to explain all the excess gamma rays, indicating that DM annihilation and/or an unresolved population of MSPs is also needed. △ Less

Submitted 28 October, 2014; originally announced October 2014.

Comments: 9 pages, 3 figures, conference proceedings for Cosmic Rays and the InterStellar Medium - CRISM 2014

Cosmic-Ray Models of the Ridge-Like Excess of Gamma Rays in the Galactic Center

Authors: Oscar Macias, Chris Gordon, Roland Crocker, Stefano Profumo

Abstract: The High-Energy Stereoscopic System (HESS) has detected diffuse TeV emission correlated with the distribution of molecular gas along the Ridge at the Galactic Center. Diffuse, non-thermal emission is also seen by the Fermi large area telescope (Fermi-LAT) in the GeV range and by radio telescopes in the GHz range. Additionally, there is a distinct, spherically symmetric excess of gamma rays seen by… ▽ More The High-Energy Stereoscopic System (HESS) has detected diffuse TeV emission correlated with the distribution of molecular gas along the Ridge at the Galactic Center. Diffuse, non-thermal emission is also seen by the Fermi large area telescope (Fermi-LAT) in the GeV range and by radio telescopes in the GHz range. Additionally, there is a distinct, spherically symmetric excess of gamma rays seen by Fermi-LAT in the GeV range. A cosmic ray flare, occurring in the Galactic Center, $10^4$ years ago has been proposed to explain the TeV Ridge. An alternative, steady-state model explaining all three data sets (TeV, GeV, and radio) invokes purely leptonic processes. We show that the flare model from the Galactic Center also provides an acceptable fit to the GeV and radio data, provided the diffusion coefficient is energy independent. However, if Kolmogorov-type turbulence is assumed for the diffusion coefficient, we find that two flares are needed, one for the TeV data (occurring approximately $10^4 $ years ago) and an older one for the GeV data (approximately $10^5$ years old). We find that the flare models we investigate do not fit the spherically symmetric GeV excess as well as the usual generalized Navarro-Frenk-White spatial profile, but are better suited to explaining the Ridge. We also show that a range of single-zone, steady-state models are able to explain all three spectral data sets. Large gas densities equal to the volumetric average in the region can be accommodated by an energy independent diffusion or streaming based steady-state model. Additionally, we investigate how the flare and steady-state models may be distinguished with future gamma-ray data looking for a spatial dependence of the gamma-ray spectral index. △ Less

Submitted 10 June, 2015; v1 submitted 7 October, 2014; originally announced October 2014.

Comments: V3: Minor changes made. Matches version published by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 2015 451 (2): 1833-1847

Multimessenger Search for Sources of Gravitational Waves and High-Energy Neutrinos: Results for Initial LIGO-Virgo and IceCube

Authors: The IceCube Collaboration, The LIGO Scientific Collaboration, The Virgo Collaboration, M. G. Aartsen, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, T. Anderson, C. Arguelles, T. C. Arlen, J. Auffenberg, X. Bai, S. W. Barwick, V. Baum, J. J. Beatty, J. Becker Tjus, K. -H. Becker, S. BenZvi, P. Berghaus, D. Berley, E. Bernardini, A. Bernhard , et al. (1166 additional authors not shown)

Abstract: We report the results of a multimessenger search for coincident signals from the LIGO and Virgo gravitational-wave observatories and the partially completed IceCube high-energy neutrino detector, including periods of joint operation between 2007-2010. These include parts of the 2005-2007 run and the 2009-2010 run for LIGO-Virgo, and IceCube's observation periods with 22, 59 and 79 strings. We find… ▽ More We report the results of a multimessenger search for coincident signals from the LIGO and Virgo gravitational-wave observatories and the partially completed IceCube high-energy neutrino detector, including periods of joint operation between 2007-2010. These include parts of the 2005-2007 run and the 2009-2010 run for LIGO-Virgo, and IceCube's observation periods with 22, 59 and 79 strings. We find no significant coincident events, and use the search results to derive upper limits on the rate of joint sources for a range of source emission parameters. For the optimistic assumption of gravitational-wave emission energy of $10^{-2}$\,M$_\odot$c$^2$ at $\sim 150$\,Hz with $\sim 60$\,ms duration, and high-energy neutrino emission of $10^{51}$\,erg comparable to the isotropic gamma-ray energy of gamma-ray bursts, we limit the source rate below $1.6 \times 10^{-2}$\,Mpc$^{-3}$yr$^{-1}$. We also examine how combining information from gravitational waves and neutrinos will aid discovery in the advanced gravitational-wave detector era. △ Less

Submitted 9 October, 2014; v1 submitted 3 July, 2014; originally announced July 2014.

Journal ref: Phys. Rev. D 90, 102002 (2014)

Multipole analysis of IceCube data to search for dark matter accumulated in the Galactic halo

Authors: IceCube Collaboration, M. G. Aartsen, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, T. Anderson, C. Arguelles, T. C. Arlen, J. Auffenberg, X. Bai, S. W. Barwick, V. Baum, J. J. Beatty, J. Becker Tjus, K. -H. Becker, S. BenZvi, P. Berghaus, D. Berley, E. Bernardini, A. Bernhard, D. Z. Besson, G. Binder , et al. (277 additional authors not shown)

Abstract: Dark matter which is bound in the Galactic halo might self-annihilate and produce a flux of stable final state particles, e.g. high energy neutrinos. These neutrinos can be detected with IceCube, a cubic-kilometer sized Cherenkov detector. Given IceCube's large field of view, a characteristic anisotropy of the additional neutrino flux is expected. In this paper we describe a multipole method to se… ▽ More Dark matter which is bound in the Galactic halo might self-annihilate and produce a flux of stable final state particles, e.g. high energy neutrinos. These neutrinos can be detected with IceCube, a cubic-kilometer sized Cherenkov detector. Given IceCube's large field of view, a characteristic anisotropy of the additional neutrino flux is expected. In this paper we describe a multipole method to search for such a large-scale anisotropy in IceCube data. This method uses the expansion coefficients of a multipole expansion of neutrino arrival directions and incorporates signal-specific weights for each expansion coefficient. We apply the technique to a high-purity muon neutrino sample from the Northern Hemisphere. The final result is compatible with the null-hypothesis. As no signal was observed, we present limits on the self-annihilation cross-section averaged over the relative velocity distribution $<σv>$ down to $1.9\cdot 10^{-23}\,\mathrm{cm}^3\mathrm{s}^{-1}$ for a dark matter particle mass of $700\,\mathrm{GeV}$ to $1000\,\mathrm{GeV}$ and direct annihilation into $ν\barν$. The resulting exclusion limits come close to exclusion limits from $γ$-ray experiments, that focus on the outer Galactic halo, for high dark matter masses of a few TeV and hard annihilation channels. △ Less

Submitted 22 October, 2014; v1 submitted 26 June, 2014; originally announced June 2014.

Comments: 17 pages, 18 figures, 5 tables

Journal ref: European Physical Journal C 01/2015; 75(20)

Searches for Extended and Point-like Neutrino Sources with Four Years of IceCube Data

Authors: IceCube Collaboration, M. G. Aartsen, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, T. Anderson, C. Arguelles, T. C. Arlen, J. Auffenberg, X. Bai, S. W. Barwick, V. Baum, J. J. Beatty, J. Becker Tjus, K. -H. Becker, S. BenZvi, P. Berghaus, D. Berley, E. Bernardini, A. Bernhard, D. Z. Besson, G. Binder , et al. (279 additional authors not shown)

Abstract: We present results on searches for point-like sources of neutrinos using four years of IceCube data, including the first year of data from the completed 86-string detector. The total livetime of the combined dataset is 1,373 days. For an E$^{-2}$ spectrum the median sensitivity at 90\% C.L. is $\sim 10^{-12}$ TeV$^{-1}$cm$^{-2}$s$^{-1}$ for energies between 1 TeV$-$1 PeV in the northern sky and… ▽ More We present results on searches for point-like sources of neutrinos using four years of IceCube data, including the first year of data from the completed 86-string detector. The total livetime of the combined dataset is 1,373 days. For an E$^{-2}$ spectrum the median sensitivity at 90\% C.L. is $\sim 10^{-12}$ TeV$^{-1}$cm$^{-2}$s$^{-1}$ for energies between 1 TeV$-$1 PeV in the northern sky and $\sim 10^{-11}$ TeV$^{-1}$cm$^{-2}$s$^{-1}$ for energies between 100 TeV $-$ 100 PeV in the southern sky. The sensitivity has improved from both the additional year of data and the introduction of improved reconstructions compared to previous publications. In addition, we present the first results from an all-sky search for extended sources of neutrinos. We update results of searches for neutrino emission from stacked catalogs of sources, and test five new catalogs; two of Galactic supernova remnants and three of active galactic nuclei. In all cases, the data are compatible with the background-only hypothesis, and upper limits on the flux of muon neutrinos are reported for the sources considered. △ Less

Submitted 25 June, 2014; originally announced June 2014.

Comments: 36 pages, 15 figures. Submitted to the Astrophysical Journal

Journal ref: Astrophys.J. 796:109,2014

Observation of High-Energy Astrophysical Neutrinos in Three Years of IceCube Data

Authors: M. G. Aartsen, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, D. Altmann, T. Anderson, C. Arguelles, T. C. Arlen, J. Auffenberg, X. Bai, S. W. Barwick, V. Baum, J. J. Beatty, J. Becker Tjus, K. -H. Becker, S. BenZvi, P. Berghaus, D. Berley, E. Bernardini, A. Bernhard, D. Z. Besson, G. Binder, D. Bindig , et al. (272 additional authors not shown)

Abstract: A search for high-energy neutrinos interacting within the IceCube detector between 2010 and 2012 provided the first evidence for a high-energy neutrino flux of extraterrestrial origin. Results from an analysis using the same methods with a third year (2012-2013) of data from the complete IceCube detector are consistent with the previously reported astrophysical flux in the 100 TeV - PeV range at t… ▽ More A search for high-energy neutrinos interacting within the IceCube detector between 2010 and 2012 provided the first evidence for a high-energy neutrino flux of extraterrestrial origin. Results from an analysis using the same methods with a third year (2012-2013) of data from the complete IceCube detector are consistent with the previously reported astrophysical flux in the 100 TeV - PeV range at the level of $10^{-8}\, \mathrm{GeV}\, \mathrm{cm}^{-2}\, \mathrm{s}^{-1}\, \mathrm{sr}^{-1}$ per flavor and reject a purely atmospheric explanation for the combined 3-year data at $5.7 σ$. The data are consistent with expectations for equal fluxes of all three neutrino flavors and with isotropic arrival directions, suggesting either numerous or spatially extended sources. The three-year dataset, with a livetime of 988 days, contains a total of 37 neutrino candidate events with deposited energies ranging from 30 to 2000 TeV. The 2000 TeV event is the highest-energy neutrino interaction ever observed. △ Less

Submitted 2 July, 2014; v1 submitted 21 May, 2014; originally announced May 2014.

Comments: 8 pages, 5 figures. Accepted by PRL. The event catalog, event displays, and other data tables are included after the final page of the article. Changed from the initial submission to reflect referee comments, expanding the section on atmospheric backgrounds, and fixes offsets of up to 0.9 seconds in reported event times. Address correspondence to: J. Feintzeig, C. Kopper, N. Whitehorn

Journal ref: Phys. Rev. Lett. 113, 101101 (2014)

Search for non-relativistic Magnetic Monopoles with IceCube

Authors: IceCube Collaboration, M. G. Aartsen, R. Abbasi, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, D. Altmann, C. Arguelles, T. C. Arlen, J. Auffenberg, X. Bai, M. Baker, S. W. Barwick, V. Baum, R. Bay, J. J. Beatty, J. Becker Tjus, K. -H. Becker, M. L. Benabderrahmane, S. BenZvi, P. Berghaus, D. Berley, E. Bernardini, A. Bernhard , et al. (266 additional authors not shown)

Abstract: The IceCube Neutrino Observatory is a large Cherenkov detector instrumenting $1\,\mathrm{km}^3$ of Antarctic ice. The detector can be used to search for signatures of particle physics beyond the Standard Model. Here, we describe the search for non-relativistic, magnetic monopoles as remnants of the GUT (Grand Unified Theory) era shortly after the Big Bang. These monopoles may catalyze the decay of… ▽ More The IceCube Neutrino Observatory is a large Cherenkov detector instrumenting $1\,\mathrm{km}^3$ of Antarctic ice. The detector can be used to search for signatures of particle physics beyond the Standard Model. Here, we describe the search for non-relativistic, magnetic monopoles as remnants of the GUT (Grand Unified Theory) era shortly after the Big Bang. These monopoles may catalyze the decay of nucleons via the Rubakov-Callan effect with a cross section suggested to be in the range of $10^{-27}\,\mathrm{cm^2}$ to $10^{-21}\,\mathrm{cm^2}$. In IceCube, the Cherenkov light from nucleon decays along the monopole trajectory would produce a characteristic hit pattern. This paper presents the results of an analysis of first data taken from May 2011 until May 2012 with a dedicated slow-particle trigger for DeepCore, a subdetector of IceCube. A second analysis provides better sensitivity for the brightest non-relativistic monopoles using data taken from May 2009 until May 2010. In both analyses no monopole signal was observed. For catalysis cross sections of $10^{-22}\,(10^{-24})\,\mathrm{cm^2}$ the flux of non-relativistic GUT monopoles is constrained up to a level of $Φ_{90} \le 10^{-18}\,(10^{-17})\,\mathrm{cm^{-2}s^{-1}sr^{-1}}$ at a 90% confidence level, which is three orders of magnitude below the Parker bound. The limits assume a dominant decay of the proton into a positron and a neutral pion. These results improve the current best experimental limits by one to two orders of magnitude, for a wide range of assumed speeds and catalysis cross sections. △ Less

Submitted 11 June, 2014; v1 submitted 14 February, 2014; originally announced February 2014.

Comments: 20 pages, 20 figures

Journal ref: Eur.Phys.J. C74:2938,2014

The Contribution of Cosmic Rays Interacting With Molecular Clouds to the Galactic Center Gamma-Ray Excess

Authors: Oscar Macias, Chris Gordon

Abstract: The Fermi-LAT data appear to have an excess of gamma rays from the inner 150 pc of the Galactic Center. The main explanations proposed for this are: an unresolved population of millisecond pulsars (MSPs), dark matter (DM) annihilation, and nonthermal bremsstrahlung produced by a population of electrons interacting with neutral gas in molecular clouds. The first two options have spatial templates w… ▽ More The Fermi-LAT data appear to have an excess of gamma rays from the inner 150 pc of the Galactic Center. The main explanations proposed for this are: an unresolved population of millisecond pulsars (MSPs), dark matter (DM) annihilation, and nonthermal bremsstrahlung produced by a population of electrons interacting with neutral gas in molecular clouds. The first two options have spatial templates well fitted by the square of a generalized Navarro-Frenk-White (NFW) profile with inner slope $γ=1.2$. We model the third option with a 20-cm continuum emission Galactic Ridge template. A template based on the HESS residuals is shown to give similar results. The gamma-ray excess is found to be best fit by a combination of the generalized NFW squared template and a Galactic Ridge template. We also find the spectra of each template is not significantly affected in the combined fit and is consistent with previous single template fits. That is the generalized NFW squared spectrum can be fit by either of order 1000 unresolved MSPs or DM with mass around 30 GeV, a thermal cross section, and mainly annihilating to $b \bar{b}$ quarks. While the Galactic Ridge continues to have a spectrum consistent with a population of nonthermal electrons whose spectrum also provides a good fit to synchrotron emission measurements. We also show that the current DM fit may be hard to test, even with 10 years of Fermi-LAT data, especially if there is a mixture of DM and MSPs contributing to the signal, in which case the implied DM cross section will be suppressed. △ Less

Submitted 11 March, 2014; v1 submitted 23 December, 2013; originally announced December 2013.

Comments: 14 pages, 6 figures, minor changes, matches PRD published version

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

Search for neutrino-induced particle showers with IceCube-40

Authors: IceCube Collaboration, M. G. Aartsen, R. Abbasi, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, D. Altmann, C. Arguelles, T. C. Arlen, J. Auffenberg, X. Bai, M. Baker, S. W. Barwick, V. Baum, R. Bay, J. J. Beatty, J. Becker Tjus, K. -H. Becker, S. BenZvi, P. Berghaus, D. Berley, E. Bernardini, A. Bernhard, D. Z. Besson , et al. (265 additional authors not shown)

Abstract: We report on the search for neutrino-induced particle-showers, so-called cascades, in the IceCube-40 detector. The data for this search was collected between April 2008 and May 2009 when the first 40 IceCube strings were deployed and operational. Three complementary searches were performed, each optimized for different energy regimes. The analysis with the lowest energy threshold (2 TeV) targeted… ▽ More We report on the search for neutrino-induced particle-showers, so-called cascades, in the IceCube-40 detector. The data for this search was collected between April 2008 and May 2009 when the first 40 IceCube strings were deployed and operational. Three complementary searches were performed, each optimized for different energy regimes. The analysis with the lowest energy threshold (2 TeV) targeted atmospheric neutrinos. A total of 67 events were found, consistent with the expectation of 41 atmospheric muons and 30 atmospheric neutrino events. The two other analyses targeted a harder, astrophysical neutrino flux. The analysis with an intermediate threshold of 25 TeV lead to the observation of 14 cascade-like events, again consistent with the prediction of 3.0 atmospheric neutrino and 7.7 atmospheric muon events. We hence set an upper limit of $E^2 Φ_{lim} \leq 7.46\times10^{-8}\,\mathrm{GeV sr^{-1} s^{-1} cm^{-2}}$ (90% C.L.) on the diffuse flux from astrophysical neutrinos of all neutrino flavors, applicable to the energy range 25 TeV to 5 PeV, assuming an $E_ν^{-2}$ spectrum and a neutrino flavor ratio of 1:1:1 at the Earth. The third analysis utilized a larger and optimized sample of atmospheric muon background simulation, leading to a higher energy threshold of 100 TeV. Three events were found over a background prediction of 0.04 atmospheric muon events and 0.21 events from the flux of conventional and prompt atmospheric neutrinos. Including systematic errors this corresponds to a $2.7σ$ excess with respect to the background-only hypothesis. Our observation of neutrino event candidates above 100 TeV complements IceCube's recently observed evidence for high-energy astrophysical neutrinos. △ Less

Submitted 30 November, 2013; originally announced December 2013.

Comments: 19 pages, 13 figures, submitted to Phys. Rev. D

Journal ref: Phys.Rev. D89:102001,2014

Search for a diffuse flux of astrophysical muon neutrinos with the IceCube 59-string configuration

Authors: The IceCube Collaboration, M. G. Aartsen, R. Abbasi, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, D. Altmann, C. Arguelles, T. C. Arlen, J. Auffenberg, X. Bai, M. Baker, S. W. Barwick, V. Baum, R. Bay, J. J. Beatty, J. Becker Tjus, K. -H. Becker, S. BenZvi, P. Berghaus, D. Berley, E. Bernardini, A. Bernhard, D. Z. Besson , et al. (265 additional authors not shown)

Abstract: A search for high-energy neutrinos was performed using data collected by the IceCube Neutrino Observatory from May 2009 to May 2010, when the array was running in its 59-string configuration. The data sample was optimized to contain muon neutrino induced events with a background contamination of atmospheric muons of less than 1%. These data, which are dominated by atmospheric neutrinos, are analyz… ▽ More A search for high-energy neutrinos was performed using data collected by the IceCube Neutrino Observatory from May 2009 to May 2010, when the array was running in its 59-string configuration. The data sample was optimized to contain muon neutrino induced events with a background contamination of atmospheric muons of less than 1%. These data, which are dominated by atmospheric neutrinos, are analyzed with a global likelihood fit to search for possible contributions of prompt atmospheric and astrophysical neutrinos, neither of which have yet been identified. Such signals are expected to follow a harder energy spectrum than conventional atmospheric neutrinos. In addition, the zenith angle distribution differs for astrophysical and atmospheric signals. A global fit of the reconstructed energies and directions of observed events is performed, including possible neutrino flux contributions for an astrophysical signal and atmospheric backgrounds as well as systematic uncertainties of the experiment and theoretical predictions. The best fit yields an astrophysical signal flux for $ν_μ+ \barν_μ$ of $E^2 \cdot Φ(E) = 0.25 \cdot 10^{-8} \textrm{GeV} \textrm{cm}^{-2} \textrm{s}^{-1} \textrm{sr}^{-1}$, and a zero prompt component. Although the sensitivity of this analysis for astrophysical neutrinos surpasses the Waxman and Bahcall upper bound, the experimental limit at 90% confidence level is a factor of 1.5 above at a flux of $E^2 \cdot Φ(E) = 1.44 \cdot 10^{-8} \textrm{GeV} \textrm{cm}^{-2} \textrm{s}^{-1} \textrm{sr}^{-1}$. △ Less

Submitted 27 November, 2013; originally announced November 2013.

Comments: 18 pages, 13 figures. Submitted to Phys. Rev. D

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

Energy Reconstruction Methods in the IceCube Neutrino Telescope

Authors: IceCube Collaboration, M. G. Aartsen, R. Abbasi, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, D. Altmann, C. Arguelles, J. Auffenberg, X. Bai, M. Baker, S. W. Barwick, V. Baum, R. Bay, J. J. Beatty, J. Becker Tjus, K. -H. Becker, S. BenZvi, P. Berghaus, D. Berley, E. Bernardini, A. Bernhard, D. Z. Besson, G. Binder , et al. (263 additional authors not shown)

Abstract: Accurate measurement of neutrino energies is essential to many of the scientific goals of large-volume neutrino telescopes. The fundamental observable in such detectors is the Cherenkov light produced by the transit through a medium of charged particles created in neutrino interactions. The amount of light emitted is proportional to the deposited energy, which is approximately equal to the neutrin… ▽ More Accurate measurement of neutrino energies is essential to many of the scientific goals of large-volume neutrino telescopes. The fundamental observable in such detectors is the Cherenkov light produced by the transit through a medium of charged particles created in neutrino interactions. The amount of light emitted is proportional to the deposited energy, which is approximately equal to the neutrino energy for $ν_e$ and $ν_μ$ charged-current interactions and can be used to set a lower bound on neutrino energies and to measure neutrino spectra statistically in other channels. Here we describe methods and performance of reconstructing charged-particle energies and topologies from the observed Cherenkov light yield, including techniques to measure the energies of uncontained muon tracks, achieving average uncertainties in electromagnetic-equivalent deposited energy of $\sim 15\%$ above 10 TeV. △ Less

Submitted 10 February, 2014; v1 submitted 19 November, 2013; originally announced November 2013.

Comments: 21 pages, 23 figures. New version reflects referee comments. Accepted by J. Instrumentation

Journal ref: JINST 9 (2014), P03009