MNRAS 000,1–19 (2021) Preprint 28 September 2021 Compiled using MNRAS LATEX style file v3.0 Evidence for a high-energy tail in the gamma-ray spectra of globular clusters Deheng Song,1¢ Oscar Macias,1,2,3y Shunsaku Horiuchi,1z Roland M. Crocker4 and David M. Nataf5 1Center for Neutrino Physics, Department of Physics, Virginia Tech, Blacksburg, VA 24061, USA 2Kavli Institute for the Physics and Mathematics of the Universe (WPI),University of Tokyo, Kashiwa, Chiba 277-8583, Japan 3GRAPPA Institute, Institute of Physics, University of Amsterdam, 1098 XH Amsterdam, The Netherlands 4Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia 5Center for Astrophysical Sciences and Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD 21218, USA Accepted XXX. Received YYY; in original form ZZZ 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 rela- tivistic 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 Fermi-LAT data from the directions of 157 Milky Way globular clusters. We find a mildly statistically significant (3.8f) 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 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.2f–has a power index of 2.79 ± 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. Key words: gamma-rays:general – globular clusters: general – pulsars:general 1 INTRODUCTION two GCs (Freire et al. 2011; Johnson et al. 2013), providing further support for this scenario. Over two dozen globular clusters (GCs) have been detected in W rays The high-energy emission from MSPs emerges from the pri- arXiv:2102.00061v3 [astro-ph.HE] 26 Sep 2021 in Fermi Large Area Telescope (LAT) data (Abdo et al. 2009a, 2010; mary electrons accelerated by them and subsequent radiation by Kong et al. 2010; Tam et al. 2011; Zhou et al. 2015; Zhang et al. secondary, relativistic electrons and positrons (4±) pair created in 2016). The millisecond pulsar (MSP) populations of those GCs are their magnetospheres. In particular, Harding et al.(2005) studied the believed to be the main source of this W-ray emission. In particular, curvature radiation (CR) of primary electrons within MSP magne- MSPs have been firmly established as W-ray sources (Verbunt et al. tospheres with a focus on GeV-scale emission. Bednarek & Sitarek 1996; Abdo et al. 2009b; Espinoza et al. 2013; Abdo et al. 2013) (2007) then considered a scenario in which 4±, injected by MSPs, and a large fraction of them have been discovered in GCs (Camilo gradually diffuse through a GC, up-scattering ambient photons, & Rasio 2005). Recently, Fermi-LAT detected W-ray pulsations in and thus producing inverse Compton (IC) W-ray emission in the GeV−TeV energy range. Venter et al.(2009) calculated the CR and IC spectra for an ensemble of MSPs in the GCs 47 Tucanae and ¢ [email protected] Terzan 5. Cheng et al.(2010) found that the spectra of 47 Tucanae y [email protected] and seven other GCs can be explained by IC alone, invoking back- z [email protected] ground photons from the cosmic microwave background (CMB) © 2021 The Authors 2 Song et al. or Galactic infrared/stellar radiation. For a review of the observa- we consider the photon field energy densities, the stellar encounter tions and models about the W-ray emission from globular clusters, rate, and the metallicity of the GCs. Unlike previous studies of cor- see Tam et al.(2016). In general, the GeV emission mechanism of relations of the GC W-ray emissions, we consider also the upper MSPs remains in contention with CR, IC, and synchrotron radiation limits placed by null detections, which we implement via a Kendall all proposed (Harding 2021). g coefficient test statistic. Furthermore, we also look for evidence of Here, motivated by the increasing number of GCs detected IC from the spectra of GCs. For the first time, we implement a uni- in W rays, we perform a collective statistical study of their prop- versal two-component model to study the spectra of W-ray-detected erties in order to gain insight into the nature of their high-energy GCs. The two-component model comprises a CR component, which emission. Our particular aim is to investigate the importance of the is spectrally curved, plus an IC component modeled as a power law contribution of IC emission to the overall W-ray emission of GCs. in the energy range of interest. Relations between the detected GC gamma-ray luminosities The remainder of the paper is as follows: In Section2 we dis- and properties of GC can be used to probe the origins of W-ray cuss the choice of GC samples and the data analysis procedure. emission and their underlying sources. For example, correlations Section3 presents the methodology and results of our correlation with the photon field energy density of GCs could unveil the poten- analysis. Section4 describes the spectral analysis method and re- tial contribution from IC, and correlations with the stellar encounter ports the 4± injection efficiency in the GCs. We discuss the impli- rate and metallicity could provide insight into the dynamical for- cations of our results in Section5 and conclude in Section6. mation of MPSs in GCs. Previous work here includes a study by Abdo et al.(2010) that reported a correlation between the W-ray luminosity !W and the stellar encounter rate of eight GCs. Hui et al. 2 W-RAY EMISSION FROM GLOBULAR CLUSTERS (2011) studied a group of 15 W-ray emitting GCs with 2 years of In this section, we describe our choice of GC sample and the GCs’ W- Fermi data and found a positive correlation between !W and, re- spectively, encounter rate, metallicity, and Galactic photon energy ray-related parameters. The Fermi data analysis process is reported density. Hooper & Linden(2016) studied 25 GCs using 85 months as well. For GCs with a W-ray counterpart in the 4FGL, we update of Fermi data, and found that the W-ray luminosity function of MSPs their spectral parameters through a maximum likelihood procedure. in GCs is consistent with that applying to MSPs detected in the field. For those not detected in the 4FGL, we estimate their 95% C.L. Lloyd et al.(2018) studied W-ray emission from high-latitude GCs W-ray upper limits. and its connection to their X-ray emission. de Menezes et al.(2019) reanalysed 9 years of Fermi data and found 23 W-ray emitting GCs; 2.1 Globular cluster sample they found that the metallicity only mildly contributes to !W while a very high encounter rate seemed to reduce the !W from GCs. We consider the Harris(1996) catalog (2010 edition), which con- In parallel, modeling of GCs’ observed broadband spectral tains identifications and basic parameters for 157 GCs in the Milky energy distributions provides a handle on their CR and IC emissions. Way. Here, we reanalyse publicly available Fermi-LAT data from Recently, Kopp et al.(2013) and Ndiyavala et al.(2019) modelled the direction of all GCs in the Harris(1996) catalog. Figure1 shows the multiwavelength emission from MSPs considering a potential the spatial distribution of the GCs. The top panel shows the pro- CR origin for GeV and IC emissions for TeV W rays, as well as jected direction of the GCs on the celestial plane while the bottom synchrotron radiation for the radio and X-ray wavebands. These two panels display their 3D coordinates. The GCs which are de- models are successful in explaining the multiwavelength spectra tected in the 4FGL are marked by red stars while null detections are of Terzan 5. However, Terzan 5 is the only GC (perhaps) detected indicated by green circles. Most W-ray GCs are near the Sun (yellow above TeV energies (H. E. S. S. Collaboration et al. 2011). Detailed circle) or located in the Galactic bulge (assumed to be sphere of 3 spectral modelling similar to that presented by Kopp et al.(2013) kpc radius, grey circular area). and Ndiyavala et al.(2019) is difficult for other GCs at present due The origin of the W-ray emission from GCs can be studied by to a lack of TeV W-ray data. Although Fermi-LAT is sensitive to comparing its dependency on GC properties. IC emission is sensi- ± W-rays of up to ≈ 1 TeV, the photon count statistics at the highest tive to the ambient photon field on which the 4 scatter. Hui et al. energies are very low. (2011) reported a positive correlation between the W-ray luminosity In the recently published Fermi-LAT fourth source cata- !W and the photon field energy density at the cluster location, indi- log (Abdollahi et al.