Arxiv:2106.15290V1 [Gr-Qc] 29 Jun 2021 the Aim of This Work Is to Account Fully for the Contri- Different Types of Object Following the Same Geodesic
Total Page:16
File Type:pdf, Size:1020Kb
Constraints on Equivalence Principle Violation from Gamma Ray Bursts D. J. Bartlett,1, ∗ D. Bergsdal,2, 3 H. Desmond,1 P. G. Ferreira,1 and J. Jasche3 1Astrophysics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford, OX1 3RH, UK 2Particle and Astroparticle Physics, Department of Physics, KTH Royal Institute of Technology, SE 106 91, Stockholm, Sweden 3The Oskar Klein Centre, Department of Physics, Stockholm University, AlbaNova University Centre, SE 106 91 Stockholm, Sweden Theories of gravity that obey the Weak Equivalence Principle have the same Parametrised Post- Newtonian parameter γ for all particles at all energies. The large Shapiro time delays of extragalactic sources allow us to put tight constraints on differences in γ between photons of different frequencies from spectral lag data, since a non-zero ∆γ would result in a frequency-dependent arrival time. The majority of previous constraints have assumed that the Shapiro time delay is dominated by a few local massive objects, although this is a poor approximation for distant sources. In this work we consider the cosmological context of these sources by developing a source-by-source, Monte Carlo-based forward model for the Shapiro time delays by combining constrained realisations of the local density field using the BORG algorithm with unconstrained large-scale modes. Propagating uncertainties in the density field reconstruction and marginalising over an empirical model describing other contributions to the time delay, we use spectral lag data of Gamma Ray Bursts from the BATSE satellite to constrain ∆γ < 3:4 × 10−15 at 1σ confidence between photon energies of 25 keV and 325 keV. I. INTRODUCTION they experience time delays identically. This can be achieved either when the source emits light at varying A basic feature of metric theories of gravity is that frequencies which can be independently measured [6], or objects on either timelike or null geodesics experience when it emits other types of energy in addition to light, apparent time delays due to motion through regions of such as gravitational waves [7]. Within the PPN frame- varying spacetime curvature, caused by the difference work, this constrains the difference between their γ fac- between proper and observer time induced by a grav- tors. itational field. This effect, first derived by Shapiro in In general metric theories the time delay is propor- 1964 [1], has now been accurately measured in the Solar tional to the integral of the fluctuation in Newtonian po- System: the best constraint on the fractional deviation tential φ along the line of sight to the source; within the of the time delay from the prediction of General Rel- PPN framework it has a prefactor 1 + γ where γ = 1 ativity (GR) is currently at the 10−5 level, using data in GR. Predicting time delays is therefore equivalent to from the radio link with the Cassini spacecraft [2]. In determining δφ. In typical analyses, δφ φ is modelled the Parametrised Post-Newtonian (PPN) framework the as arising from one or a few isolated sources' near the line time delay is governed by the light-bending parameter γ of sight that are believed to be predominantly respon- [3]. For sources at greater distance, however, the lack of sible for sourcing the potential. However, as φ 1=r knowledge of the time of emission of a signal precludes one should expect it to be sensitive to the distribution∼ of direct measurement of the delay. distant sources, and indeed the contribution to φ from a Nevertheless, the time delay effect has found use in spherical shell at distance r from a test point diverges as testing various aspects of the theory of gravity and the r . This casts doubt on the multiple point masses standard model of cosmology. This is done by comparing approximation! 1 [8, 9]. time delays, either between nearby geodesics or between arXiv:2106.15290v1 [gr-qc] 29 Jun 2021 The aim of this work is to account fully for the contri- different types of object following the same geodesic. The butions to the time delay from all mass in the non-linear former may be achieved by comparing the time of recep- cosmological density field. We derive the contribution tion of photons originating from a common source but from local structures using constrained density fields gen- traversing different paths due to gravitational lensing. erated by the BORG algorithm [10{14] and combine this When the source is a time-varying quasar and the lens with an unconstrained contribution from distant sources a single massive elliptical galaxy or cluster, this method to produce a Monte-Carlo based source-by-source for- has been used to constrain the Hubble parameter, H0, ward model for the expected Shapiro time delay. on which time delay distances depend [4, 5]. Conversely, We apply our method to high-energy astrophysical comparing time delays between different objects travel- sources that have previously been used to test the ling along the same geodesic allows the Weak Equivalence Equivalence Principle for photons of varying wavelength. Principle (WEP) to be tested by investigating whether Marginalising over uncertainties in the density field re- construction and parameters describing non-Equivalence- Principle-violating contributions, we compare our pre- ∗ [email protected] dictions to the observations via a Markov Chain Monte 2 Carlo (MCMC) algorithm. We find our constraints are which can be parametrised as a difference in the PPN pa- 40 times tighter than literature results, illustrating the rameter γ between two wavelengths. Thus by comparing benefit∼ of using complete mass distributions when study- the measured time delays with the line-of-sight integral ing non-local relativistic effects such as time delays in a in Equation 2, one can put constraints on how γ changes cosmological setting. with photon frequency, and thus on any violation of the In Section II we discuss time delays in the context of WEP. metric theories of gravity, including the importance of long wavelength modes in calculating the Shapiro time delay. We describe the Gamma Ray Burst (GRB) obser- B. Time delay angular power spectrum vations in Section III and detail our inference methods in Section IV. The results are presented in Section V. The angular power spectrum of the gravitational time A discussion of systematic uncertainties, future applica- delay is [8] tions of our methods and a comparison to the literature Z Z rs 2 are presented in Section VI. We conclude in Section VII. 2 2 C` = dk k Pφ (k) dr D (r) j` (kr) ; (4) π 0 II. TIME DELAY DIFFERENCES IN METRIC such that the mean-squared time delay is THEORIES OF GRAVITY X 2` + 1 t2 = C ; (5) grav 4π ` A. Equivalence principle violation ` where Pφ is the power spectrum of the potential. For Consider a perturbed FRW-like metric in the Newto- k 1=rs, the spherical Bessel function can be expanded nian gauge, as 2 2 2 2 pπ ds = (1 + 2γδφ) dt + a (t) (1 2δφ) dr ; (1) j (x) = x` + x2+` ; (6) − − ` `+1 3 2 Γ ` + 2 O where we have introduced the PPN parameter γ to allow for deviations from General Relativity (where γ = 1). so that For a massless particle, the gravitational time delay to 4 pπ Z dk C = [T (k)]2 A kns−1k2` a source at distance rs is given by (1 + γ) tgrav, where ` 3 `+1 3 φ s π 2 Γ ` + 2 k (7) Z rs Z 2 ! ` 2+2` tgrav = dr δφ0(r)D(r); (2) dr D (r) r + k ; − 0 O where D(r) is the linear growth factor, δφ is the poten- 0 where Tφ is the potential transfer function and we assume tial fluctuation evaluated using the present-day matter a nearly scale-invariant primordial power spectrum field and we have used the weak field limit δφ 1. For P ns−1 a set of point masses φ = i GMi=ri, but in a general R (k) = Ask : (8) density field δφ must be found by solving the Poisson P equation On super-Hubble scales, Tφ (k) constant, so the small- est power of k in the integral has∼ an exponent n 2+2`. s − 2δφ(r) = 4πGδρ(r); (3) For ` 1, this is ns, which, given that ns 0:97 r [15], means≥ that, for≥ ` 1, the integral does not∼ di- ≥ where δρ is the total matter density fluctuation. This is verge at small k. However for ` = 0, the exponent is most readily solved in Fourier space. Thus determining ns 2 < 1, and so the integral diverges. We therefore − − time delays amounts to mapping out the three dimen- see that the monopole diverges due to the contribution sional density field to at least the redshift of the source. from large scales. This is not to say that there is an We describe how we do this in Section IV A. infinite time delay for a given universe, but that its vari- There exist several formulations of the Equivalence ance across all possible universes is infinite. Thus when Principle, of varying strength. The WEP, our focus here, we forward model the time delays in Section IV, we will states that all freely-falling test objects follow the same only be able to predict angular fluctuations about the trajectories given the same initial conditions, irrespec- mean Shapiro time delay at a given redshift, and not its tive of their composition or structure. This requires that absolute value. all objects from a given source experience the same time delay, regardless of their composition or energy. Multiple astrophysical sources have been detected with C. Other contributions to the time delay precise timings of photons of different wavelengths.