skip to main content

Title: Evidence for a high-energy tail in the gamma-ray spectra of globular clusters
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 versus 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-yr 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 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 ± 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.  more » « less
Award ID(s):
1914409 1908960
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
5161 to 5176
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    We use Fermi-LAT data to analyse the faint gamma-ray source located at the centre of the Sagittarius (Sgr) dwarf spheroidal galaxy. In the 4FGL-DR3 catalogue, this source is associated with the globular cluster, M54. We investigate the spectral energy distribution and spatial extension of this source, with the goal of testing two hypotheses: (1) the emission is due to millisecond pulsars within M54, or (2) the emission is due to annihilating dark matter from the Sgr halo. For the pulsar interpretation, we consider a two-component model which describes both the lower-energy magnetospheric emission and possible high-energy emission arising from inverse Compton scattering. We find that this source has a point-like morphology at low energies, consistent with magnetospheric emission, and find no evidence for a higher-energy component. For the dark matter interpretation, we find the signal favours a dark matter mass of mχ = 29.6 ± 5.8 GeV and an annihilation cross section of $\sigma v = (2.1 \pm 0.59) \times 10^{-26} \, \text{cm}^3$ s−1 for the $b \bar{b}$ channel (or mχ = 8.3 ± 3.8 GeV and $\sigma v = (0.90 \pm 0.25) \times 10^{-26} \, \text{cm}^3$ s−1 for the τ+τ− channel), when adopting a J-factor of $J=10^{19.6} \, \text{GeV}^2 \, \text{cm}^{-5}$. This parameter space is consistent with gamma-ray constraints from other dwarf galaxies and with dark matter interpretations of the Galactic Centre Gamma-Ray Excess.

    more » « less
  2. Abstract

    Terzan 5 is a rich globular cluster within the galactic bulge containing 39 known millisecond pulsars, the largest known population of any globular cluster. These faint pulsars do not have sufficient signal-to-noise ratio (S/N) to measure reliable flux density or polarization information from individual observations in general. We combined over 5.2 days of archival data, at 1500 and 2000 MHz, taken with the Green Bank Telescope over the past 12 years. We created high-S/N profiles for 32 of the pulsars and determined precise rotation measures (RMs) for 28. We used the RMs, pulsar positions, and dispersion measures to map the projected parallel component of the Galactic magnetic field toward the cluster. The 〈B∣∣〉 shows a rough gradient of ∼6 nG arcsec−1(∼160 nG pc−1) or, fractionally, a change of ∼20% in the R.A. direction across the cluster, implying Galactic magnetic field variability at sub-parsec scales. We also measured average flux densitiesSνfor the pulsars, ranging from ∼10μJy to ∼2 mJy, and an average spectral indexα= −1.35, whereSννα. This spectral index is flatter than most known pulsars, likely a selection effect due to the high frequencies used in pulsar searches to mitigate dispersion and scattering. We used flux densities from each observation to constrain the scintillation properties toward the cluster, finding strong refractive modulation on timescales of months. The inferred pulsar luminosity function is roughly power law, with slope(dlogN)/(dlogL)=1at the high-luminosity end. At the low-luminosity end, there are incompleteness effects, implying that Terzan 5 contains many more pulsars.

    more » « less
  3. Abstract

    Some of the most energetic pulsars exhibit rotation-modulatedγ-ray emission in the 0.1–100 GeV band. The luminosity of this emission is typically 0.1%–10% of the pulsar spin-down power (γ-ray efficiency), implying that a significant fraction of the available electromagnetic energy is dissipated in the magnetosphere and reradiated as high-energy photons. To investigate this phenomenon we model a pulsar magnetosphere using 3D particle-in-cell simulations with strong synchrotron cooling. We particularly focus on the dynamics of the equatorial current sheet where magnetic reconnection and energy dissipation take place. Our simulations demonstrate that a fraction of the spin-down power dissipated in the magnetospheric current sheet is controlled by the rate of magnetic reconnection at microphysical plasma scales and only depends on the pulsar inclination angle. We demonstrate that the maximum energy and the distribution function of accelerated pairs is controlled by the available magnetic energy per particle near the current sheet, the magnetization parameter. The shape and the extent of the plasma distribution is imprinted in the observed synchrotron emission, in particular, in the peak and the cutoff of the observed spectrum. We study how the strength of synchrotron cooling affects the observed variety of spectral shapes. Our conclusions naturally explain why pulsars with higher spin-down power have wider spectral shapes and, as a result, lowerγ-ray efficiency.

    more » « less
  4. ABSTRACT The afterglow emission from gamma-ray bursts (GRBs) is believed to originate from a relativistic blast wave driven into the circumburst medium. Although the afterglow emission from radio up to X-ray frequencies is thought to originate from synchrotron radiation emitted by relativistic, non-thermal electrons accelerated by the blast wave, the origin of the emission at high energies (HE; ≳GeV) remains uncertain. The recent detection of sub-TeV emission from GRB 190114C by the Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) raises further debate on what powers the very high energy (VHE; ≳300 GeV) emission. Here, we explore the inverse Compton scenario as a candidate for the HE and VHE emissions, considering two sources of seed photons for scattering: synchrotron photons from the blast wave (synchrotron self-Compton or SSC) and isotropic photon fields external to the blast wave (external Compton). For each case, we compute the multiwavelength afterglow spectra and light curves. We find that SSC will dominate particle cooling and the GeV emission, unless a dense ambient infrared photon field, typical of star-forming regions, is present. Additionally, considering the extragalactic background light attenuation, we discuss the detectability of VHE afterglows by existing and future gamma-ray instruments for a wide range of model parameters. Studying GRB 190114C, we find that its afterglow emission in the Fermi-Large Area Telescope (LAT) band is synchrotron dominated. The late-time Fermi-LAT measurement (i.e. t ∼ 104 s), and the MAGIC observation also set an upper limit on the energy density of a putative external infrared photon field (i.e. ${\lesssim} 3\times 10^{-9}\, {\rm erg\, cm^{-3}}$), making the inverse Compton dominant in the sub-TeV energies. 
    more » « less
  5. Abstract

    The number of millisecond pulsars (MSPs) observed in Milky Way globular clusters has increased explosively in recent years, but the underlying population is still uncertain due to observational biases. We use state-of-the-artN-body simulations to study the evolution of MSP populations in dense star clusters. These cluster models span a wide range in initial conditions, including different initial masses, metallicities, and virial radii, which nearly cover the full range of properties exhibited by the population of globular clusters in the Milky Way. We demonstrate how different initial cluster properties affect the number of MSPs, for which we provide scaling relations as a function of cluster age and mass. As an application, we use our formulae to estimate the number of MSPs delivered to the Galactic center from inspiralling globular clusters to probe the origin of the Galactic-center gamma-ray excess detected by Fermi. We predict about 400 MSPs in the Galactic center from disrupted globular clusters, which can potentially explain most of the observed gamma-ray excess.

    more » « less