Search for a command to run...
Context . The Galactic-Centre γ -ray excess (GCE) is a spatially extended surplus of γ -rays around Sgr A* observed by Fermi -LAT that exceeds predictions from standard cosmic-ray interactions in the inner Galaxy. The GCE has been mainly attributed either to dark-matter annihilation or to an unresolved population of millisecond pulsars (MSPs) dynamically delivered to the inner Galaxy by globular clusters (GCs). Aims . We reassess the MSP interpretation by following a fully dynamical framework to explore how neutron stars (NSs) produced in GCs are stripped and deposited into the central kiloparsec. We also quantify the resulting γ -ray emission to compare it with Fermi -LAT measurements. Methods . We ran high-resolution direct N -body simulations of GCs in a time-varying Milky Way (MW) potential, capturing the internal cluster dynamics and external tides. We modelled two channels: (i) present-day observed MW GCs on their orbits and (ii) a population of early, now-destroyed clusters whose debris was accreted during their evolution into the inner Galaxy. The simulations provide the full phase-space distribution of deposited NSs to central 1 kpc from these both sources. We converted the NS count to an MSP count using an empirically calibrated efficiency, adopting a sample-averaged ratio between observed and predicted MSP counts and NSs that are bound to clusters. As a result, we were able to generate mock sky maps and cumulative flux profiles assuming representative per-pulsar luminosities. Results . We find that MSPs inferred from the observed GCs already supply a substantial γ -ray signal, including disrupted clusters increases both amplitude and central concentration. By taking the observational ratio of MSPs to the total number of NSs in observed Galactic stellar systems as a basis and adopting an average luminosity of ⟨ L ⟩∼8 × 10 33 erg s −1 , we were able to reproduce the GCE quite well. As a free parameter to achieve better agreement with the observed flux, we had to increase the number of NSs originating from previously disrupted stellar systems by approximately a factor of 2. The deposited NSs from destroyed clusters exhibit an axisymmetric morphology with pronounced over-densities in the Galactic plane and perpendicular to it. Conclusions . The results of our modelling favour an MSP origin of the GalC γ -ray excess over dark matter annihilation, primarily because the combined contribution of MSPs delivered by surviving and disrupted GCs naturally reproduces both the amplitude and concentration of the observed signal under reasonable assumptions.