Galactic Shapiro delay to the Crab pulsar and limit on weak equivalence principle violation

We calculate the total galactic Shapiro delay to the Crab pulsar by including the contributions from the dark matter as well as baryonic matter along the line of sight. The total delay due to dark matter potential is about 3.4 days. For baryonic matter, we included the contributions from both the bulge and the disk, which are approximately 0.12 and 0.32 days respectively. The total delay from all the matter distribution is therefore 3.84 days. We also calculate the limit on violations of Weak equivalence principle by using observations of “nano-shot” giant pulses from the Crab pulsar with time-delay <0.4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$<0.4$$\end{document} ns, as well as using time differences between radio and optical photons observed from this pulsar. Using the former, we obtain a limit on violation of Weak equivalence principle in terms of the PPN parameter Δγ<2.41×10-15\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varDelta \gamma < 2.41\times 10^{-15}$$\end{document}. From the time-difference between simultaneous optical and radio observations, we get Δγ<1.54×10-9\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varDelta \gamma < 1.54\times 10^{-9}$$\end{document}. We also point out differences in our calculation of Shapiro delay and that from two recent papers (Yang and Zhang, Phys Rev D 94(10):101501, 2016; Zhang and Gong, Astrophys J 837:134, 2017), which used the same observations to obtain a corresponding limit on Δγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varDelta \gamma $$\end{document}.


Introduction
In 1964, Shapiro [3] pointed out that the round-trip time of an electromagnetic pulse to the inner planets of our solar system experiences a delay due to the non-zero gravitational potential of the Sun. This delay is referred to in the literature as "Shapiro delay" and has been measured precisely in the solar system for more than five decades, enabling very stringent tests of general relativity (GR) [4] and also an astrophysical probe to measure neutron star masses in binary systems [5]. Following the detection of neutrinos from SN 1987A [6,7], it was pointed out that the neutrinos also encountered a Shapiro a e-mail: shntn05@gmail.com delay of about 1-6 months due to the gravitational potential of the intervening matter along the line of sight [8,9]. From the first measured GW signal GW150914 [10], one can deduce that the Shapiro delay for GWs is frequency independent [11,12]. The most recent electromagnetic observations seen in association with GW170817, show that gravitational waves also experience the same Shapiro delay as photons to about O(10 −8 ) [13][14][15][16][17].
Wei et al. [18] (following earlier suggestions in [19]), have pointed out that one can use the relative time difference between different astrophysics messengers seen across a broad swath of frequencies, to constrain energy-dependent or inter-messenger Shapiro delay violations, enabling us to set a stringent limit on any violations of weak equivalence principle. 1 WEP violation is predicted in many quantum gravity models [20] and the equality of Shapiro delays across the electromagnetic spectrum can be used to test for such violations. The violation of WEP is usually parameterized in terms of the post-Newtonian parameter Δγ [4]. This technique has been applied to EM observations from a wide variety of extra-galactic astrophysical objects such as FRBs [18,21], blazars [22], GRBs [23,24], etc. A complete summary of limits on WEP principle can be found in Wu et al. [25]. Most recently, two independent groups have using the Crab pulsar to obtain limit on violations of WEP [1,2]. Yang and Zhang [1] (hereafter Y16) have used "nano-shot" giant pulses from the Crab pulsar with time-delay between different energies of about 0.4 ns [26] to set the most stringent limit on the violation of equivalence principle of Δγ < (0.6 − 1.8) × 10 −15 . The Crab pulsar has been simultaneously timed at radio, optical, X-Ray, and γ -ray wavelengths. Zhang and Gong [2] (hereafter Z16) have used the arrival time difference between various combinations of the above observations to obtain a limit on Δγ < (2.63 − 4.01) × 10 −9 . In this work, we calculate the total Shapiro delay to the Crab pulsar (by incorporating both the dark matter and baryonic contribution) using the same method used to calculate Shapiro delay to a variety of galactic sources in our previous works [27][28][29]. From the calculated Shapiro delay, we then obtain constraints on WEP using the same observations as in Y16 and Z16.

Estimated Shapiro delay to Crab pulsar
The Crab pulsar (PSR B0531+21) is located at RA = 05 h 34 m 32 s and Dec = 22 • 52.1 at a distance of 2.2 Kpc [30]. To calculate the Shapiro delay to this pulsar, we follow the same procedure as used for the delay calculation to PSR 1937B+21 [29]. We provide a brief synopsis of the calculation. More details can be found in Refs. [27][28][29]. We assume static symmetric geometry and posit the Schwarzschild metric to model the gravitational potential of the dark matter distribution. The coefficients of the metric are obtained in terms of the density profile and mass distribution by solving Einstein's equations. For the dark matter distribution we use the NFW profile [31] and mass-halo concentration relation from Klypin et al. [32]. With these assumptions, the total Shapiro delay due to the dark matter potential turns out to be 3.4 days for a distance of 2.2 kpc and its variation with distance in the vicinity of the Crab pulsar is shown in Fig. 1. We note that the Shapiro delay for alternate dark matter density profiles has been calculated in our previous works [27,28]. Sensitivty to alternate baryonic mass profiles can be found in Ref. [29].
To calculate the Shapiro delay from the baryonic matter, we sum the contributions from both the bulge and the disk and posit spherical symmetry for the mass distributions of both of them. We assume a Hernquist profile [33] Fig. 1 The Shapiro delay due to the dark matter potential along the line of sight as a function of distance in the vicinity of the Crab pulsar and Nagai [35] profile for the mass of the disk, with the total mass equal to M disk = 5 × 10 10 M [36]. Therefore, the total mass due to the baryonic component is equal to 6.5 × 10 10 M . From the total mass, we can calculate the coefficients of the Schwarzschild metric and obtain the total Shapiro delay for both the bulge and the disk. This delay turns out to be 0.12 days for the bulge and 0.32 days for the disk. The total Shapiro delay to the Crab pulsar after summing all these contributions turns out to be 3.84 days. Y16 have used the total gravitational potential of the Milky way and with their assumptions, one obtains a Shapiro delay of (5.14-15.42) days. Z16 have assumed an NFW profile for the dark matter halo and a Miyamoto-Nagai disk for the baryonic component and have used the Milky way parameters for these from Gomez et al. [37]. From their value of the gravitational potential, the inferred Shapiro delay is about 1.98 days. Therefore, the delay which we have calculated by integrating the geodesics for the Schwarzschild metric is inbetween the values obtained by Y16 and Z16, but is of the same order of magnitude.
We now emphasize some of the key differences between our calculations and those done by Y16 and Z16. We have calculated the Shapiro delay by solving the geodesic deviation equation perturbatively by positing pressureless, static, and spherically symmetric system. To calculate the coefficients of the Schwarzschild metric, we sum the contributions from both the dark matter and the baryonic matter. Both Y16 and Z16 have assumed the point mass approximation to calculate the Shapiro delay and used the original formula derived by Shapiro in 1964 [3], where one gets the logarithm. This is valid only for the case of a first order perturbation theory in the PPN formalism. For the case of dark matter in addition to the baryonic matter, there are two small parameters: usual 2GM/Rc 2 and 2v 2 /c 2 , which is due to the constancy of the asymptotic rotation speed of spiral galaxies and these are of the same order. Therefore, the Shapiro delay calculation should be done more carefully as in this work, and the difference between the two cases can be found in Kahya and Woodard [38] and more details in Soussa and Woodard [39].
Therefore, since our calculation is a fully general relativistic one and does not assume a point source for the gravitating mass and assume a varying density, our result should be more accurate than what was done in Y16 and Z16. Furthermore, Y16 have used the galactic rotation curve from Irrgang et al. [40] extending upto 10 kpc. to calculate the contribution of the Milky way galaxy. Since the Crab is located at a distance of about 3 Kpc, incorporating the total potential upto 10 Kpc would result in an overestimate of the Shapiro delay. Moreover, Z16 have not included the contribution due to the bulge. Therefore their calculation would be a slight overestimate of the true value.

Constraints on WEP
Once the Shapiro delay for a given mass distribution is calculated along a line of sight to the Crab pulsar, if photons of different frequencies/energies arrive from the same source within a time interval (Δt), after traversing the Cosmos, one can constrain the violations of WEP in terms of the PPN parameter Δγ and the calculated Shapiro delay T shapiro [18]: We note that Δt includes contributions from the intrinsic time delay, violation of Lorentz invariance induced delays, delay due to possible non-zero photon mass, and an additional delay due to cold plasma dispersion along the line of sight (valid only for radio waves) [18,23]. If we assume that all the other time delays are zero or negligible, we can estimate a limit on Δγ [18] If we use the nano-second pulse from Crab with flux exceeding 2 MJy [26], which had Δt < 0.4 ns, we get a corresponding limit on Δγ < 2.41 × 10 −15 . We note that the radio pulse observations include a correction for the dispersion induced delay (Δt DM ). However, an additional assumption made in obtaining a limit on Δγ from the Crab nano-shot observations is that the PPN γ does not have the same 1 ν 2 dependence on frequency as Δt DM [41]. We do however know from extragalactic observations of GRBs that the PPN γ is independent of energy over a range of eV to GeV to within O10 −8 [12]. The most precise time difference from multi-frequency monitoring of the Crab pulsar is between the radio and optical wavelengths equal to 255 ± 21µs [42], where the radio data is obtained from the Nançay radio telescope and optical data from the S-CAM3 imager on the OGS telescope in Tenerife. Using this observed value of Δt, we get Δγ < 1.54 × 10 −9 . Z16 have also used other combinations of multi-wavelength observations to set limits on WEP, but these are less stringent than those obtained using radio and optical observations. Therefore, we only report results on the violation of WEP using radio and optical observations.

Conclusions
We have calculated the total Shapiro delay to the Crab pulsar due to the gravitational potential along the line of sight. We included the contributions from both the dark and baryonic matter. The total delay from the gravitational potential of the dark matter distribution to the Crab pulsar is 3.4 days. The delay from the baryonic matter is equal to 0.12 days and 0.32 days for the bulge and disk respectively. Therefore, the total Shapiro delay is equal to 3.84 days. We also reviewed the differences in our calculations of Shapiro delay and those from other groups [1,2], which also estimated this delay. Using this value for Shapiro delay, we then used the same multiwavelength observations of the Crab pulsar as in Y16 and Z16 [1,2] to obtain limits on violation of weak equivalence principle in terms of the PPN parameter γ . Similar to Y16, we use observations of "nano-shot" giant pulses with time-delay < 0.4 ns resulting in Δγ < 2.41 × 10 −15 . We then follow Z16, and use the time-differences between radio and optical photons from a multi-wavelength observing campaign of the Crab pulsar and obtain Δγ < 1.54 × 10 −9 .