Search

Intensity scintillations of radio pulsars are known to originate from interference between waves scattered by the electron density irregularities of interstellar plasma, often leading to parabolic arcs in the two-dimensional power spectrum of the recorded dynamic spectrum. The degree of arc curvature depends on the distance to the scattering plasma and its transverse velocity with respect to the line of sight. We report the observation of annual and orbital variations in the curvature of scintillation arcs over a period of 16 yr for the bright millisecond pulsar, PSR J0437-4715. These variations are the signature of the relative transverse motions of Earth, the pulsar, and the scattering medium, which we model to obtain precise measurements of parameters of the pulsar's binary orbit and the scattering medium itself. We observe two clear scintillation arcs in most of our >5000 observations, and we show that they originate from scattering by thin screens located at distances D 1 = 89.8 ± 0.4 pc and D 2 = 124 ± 3 pc from Earth. The best-fit scattering model we derive for the brightest arc yields the pulsar's orbital inclination angle, i = 137.°1 ± 0.°3, and longitude of ascending node, Ω = 206.°3 ± 0.°4. Using scintillation arcs for precise astrometry and orbital dynamics can be superior to modeling variations in the diffractive scintillation timescale, because the arc curvature is independent of variations in the level of turbulence of interstellar plasma. This technique can be used in combination with pulsar timing to determine the full three-dimensional orbital geometries of binary pulsars and provides parameters essential for testing theories of gravity and constraining neutron star masses.

There is growing evidence that the Galactic Center Excess identified in the 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 Fermi bubbles compared to previous work and the three-dimensional Inverse Compton (IC) maps recently released by the GALPROP team. Our results provide strong constraints on Millisecond Pulsar (MSP) formation scenarios proposed to explain the excess. We find that an admixture formation scenario, in which some of the relevant binaries are primordial and the rest are formed 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.

We perform a composite likelihood analysis of subdivided regions within the central 26° × 20° of the Milky Way, with the aim of characterizing the spectrum of the gamma-ray galactic center excess in regions of varying galactocentric distance. Outside of the innermost few degrees, we find that the radial profile of the excess is background-model dependent and poorly constrained. The spectrum of the excess emission is observed to extend upwards of 10 GeV outside ∼5° in radius, but cuts off steeply between 10-20 GeV only in the innermost few degrees. If interpreted as a real feature of the excess, this radial variation in the spectrum has important implications for both astrophysical and dark matter interpretations of the galactic center excess. Single-component dark matter annihilation models face challenges in reproducing this variation; on the other hand, a population of unresolved millisecond pulsars contributing both prompt and secondary inverse Compton emission may be able to explain the spectrum as well as its spatial dependency. We show that the expected differences in the photon-count distributions of a smooth dark matter annihilation signal and an unresolved point source population are an order of magnitude smaller than the fluctuations in residuals after fitting the data, which implies that mismodeling is an important systematic effect in point source analyses aimed at resolving the gamma-ray excess.