Testing weak equivalence principle with strongly lensed cosmic transients

Current constraints on Einstein’s weak equivalence principle (WEP) utilize the observed time delay between correlated particles of astronomical sources. However, the intrinsic time delay due to particle emission time and the time delays caused by potential Lorentz-invariance violation and non-zero rest mass of photons were simply omitted in previous studies. Here we propose a robust method to test WEP using strongly lensed cosmic transients, which can naturally overcome these time delays. This can be achieved by comparing the time delays between lensed images seen in different energy bands or in gravitational waves (GWs) and their electromagnetic (EM) counterparts. The power of our method mainly depends on the timing accuracy of cosmic transient and the strong lensing time delay. If the time delay of cosmic transient can be measured with accuracy about 0.1 s (e.g. gamma-ray bursts), we show that the upper limit on the differences of the parameterized post-Newtonian parameter γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document} value is Δγ<10-7\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varDelta \gamma <10^{-7}$$\end{document} with a one-month strong lensing time delay event. This accuracy of WEP can be improved by several orders, if the lens is galaxy cluster and the strongly lensed cosmic transients have much shorter duration, such as fast radio bursts.


Introduction
The Einstein weak equivalence principle (WEP) is one of the cornerstones of general relativity. It states that the trajectories of any freely falling, uncharged test bodies are independent of their energy, internal structure, or composition [1]. In the theory of parameterized post-Newtonian (PPN), the WEP requires that the PPN parameter γ of different particles or the same type of particle with different energies (hereafter, "different particles" represents both cases) should be the same. WEP has been tested by time delay of different particles. In a e-mail: fayinwang@nju.edu.cn 1964, [2] proposed that one could use the time delays between transmission of radar pulses towards either of the inner planets and detection of the echoes. The Shapiro time delay can be formulated as t Shapiro = − 1+γ c 3 r 0 r e Ψ (r )dr , where r e and r 0 are the emitting position of the particle and the position of observer, respectively, c is the speed of light, and Ψ (r ) is the gravitational potential. If the WEP is violated, the arrival times of two different particles emitted simultaneously traveling in the same gravitational potential will be different. The relative Shapiro time delay is given by Up to presently, the observed time delays of different type particles (e.g. photons, neutrinos, or gravitational waves), or the same type of particles with different energies have been used to test WEP, such as the different arrival times of photons and neutrinos from SN1987A [3,4], the time delay of a PeVenergy neutrino event associated with a giant flare of the blazar PKS B1424-418 [5], the photons in different energy bands of gamma-ray bursts (GRBs) [6,7], polarized photons of GRBs [7], radio signals at different frequency bands of fast radio bursts (FRBs) [8][9][10] and the Crab pulsar [11], and gravitational wave (GW) sources [12][13][14][15]. However, the observed time delay should include several terms: where Δt int is the intrinsic time delay between two particles due to different emission times, Δt LIV is the time delay caused by Lorentz-invariance violation, Δt spec is the potential time delay due to photons having a rest mass, and Δt DM is the time delay contributed by the dispersion of the line-ofsight free electrons.
In the high-energy range (from keV to GeV), Δt DM can be negligible. However, most of the previous work omitted all of the other potential effects based on the following assumptions. First, the observed time delay is mainly attributed to the gravitational potential so the time delays caused by other effects are omitted. Second, the time delays caused by other effects and Δt gra have the same sign. So they can give an upper limit of the violation of WEP. These two assumptions are unreasonable. For the first assumption, in some cases, the Δt int term would dominate the observed time delay especially when the observed time delay is very small, which leads to a non-physical constraint on the WEP. In addition, Δt spec and Δt LIV are strongly correlated with Δt gra in the total observed time delay. For the second assumption, in the special case that the signs of Δt other caused by other effects and Δt gra are opposite (e.g. Δt other = − 0.99 s and Δt gra = 1 s), a much smaller time delay is obtained. This will lead to a much tighter, but incorrect constraint on WEP. Therefore, the intrinsic time delay Δt int and Δt other can severely limit the ability of this method. To avoid these effects, we propose that the strongly lensed cosmic transients can be used to test WEP. This is the first time to correct the other effects in a WEP constraint. The time delay due to strong gravitational lensing between different particles is a powerful tool, which has been used to test the Lorentz-invariance violation [16] and the speed of GW [17,18].

The method and constraints on WEP
The gravitational lens effect is a prediction of general relativity. After the first observation example of gravitational lensing, the quasar QSO 0957+561A, B [19], it has become a powerful tool in many fields of astronomy, such as probing dark matter halo, large scale structures, the Hubble constant, and parameters of universe. Generally, there will be multiple images of the source when it is strongly lensed. The differences of arrival times for images are caused by the Shapiro time delay and geometric delay due to the bending of light rays [20]. In a general model of the lens, the time delay of the images relative to the case that source, lens and image are on a straight line is where θ and β are the position vectors of the image and source, z l is the redshift of the lens, d l and d s are the angular diameter distances of the lens and source, d ls is the angular diameter distance from the lens to the source and ψ(θ ) = d ls itational potential [20]. Actually, t (θ) can be divided into 2 and t gra = 1+z l c d l d s d ls ψ(θ ), which are the geometric time delay and Shapiro time delay, respectively.
We show a strong lens of a bright cosmic transient in Fig. 1. In this figure, O and L are the observer and lens Fig. 1 The geometry of the gravitational lens considered here. Illustration of a strongly lensed cosmic transient which can be used to test WEP object. S 1 and S 2 are the two signals which are associated with each other in this transient event, for example a GW and its electromagnetic (EM) counterpart or particles at two different energy bands in a cosmic explosion. d l , d s , d ls and β are the same as those in Eq. (3). θ 1 and θ 2 are the positions of two images formed by gravitational lens. P 1 and P 2 are the trajectories of the light rays of the two images. We assume that the intrinsic time delay between S 1 and S 2 is Δt int (S 1 is earlier than S 2 ) and the first and second signals of S 1 arrive at t 11 and t 12 , respectively. Here we just assume that Δt spec and Δt LIV can be omitted, and we discuss them in the following section. If the WEP is valid, one can expect that the first and second signals of S 2 will arrive at t 21 = t 11 + (1 + z s )Δt int and t 22 = t 12 + (1 + z s )Δt int , where z s is the redshift of the source. If there is some small violation of WEP, the positions of the images have small changes δθ 11 , δθ 12 , δθ 21 , δθ 22 and the arrival times of the signals become t 11 , t 12 , t 21 , and t 22 . Similar to [21], we perform a Taylor expansion of these new arrival times and consider the first order terms Because the Fermat principle, which we assume is still valid in this case, requires that the lensed image position makes the travel time stationary, ∂t ∂θ | θ 1 and ∂t ∂θ | θ 2 are equal to 0 and all the second terms of the right parts of these equations must vanish. Then comparing t 21 − t 11 and t 22 − t 12 , the effect of the intrinsic time delay t int can be removed naturally. The difference between γ 1 and γ 2 is where α = Δt geo /Δt gra is the ratio of time delays caused by geometric effect and Shapiro delay effect. Hereafter, we define 'time delay' as t 22 − t 12 or t 21 − t 11 , which is the time difference between different particles in the same light path. The 'strong lensing time delay' is defined as t 22 − t 21 or t 12 − t 11 , which is the time delay between two different paths due to strong lensing. The value of α depends on the choice of the lens model. For the purpose of illustrating our method, we use the singular isothermal sphere (SIS) model, which has been proved to be a reliable model for lenses. In the SIS model, the distribution of stars and other mass components in galaxies are thought to be like that of particles in idea gas. The projected potential and Einstein radius of a SIS lens are where σ v is the velocity dispersion and ξ is the angular distance from the center of the lens [20]. If the lensing is strong, β < θ E , there are two images of the source at the positions θ ± = β ± θ E . Then, from Eq. (3), the time delay between the two images is which is all caused by the Shapiro time delay effect (i.e. ΔT SIS = t gra, SIS ). The difference in arrival times of two images caused by geometric time delay is Δt geo = 0, which means that the lengths of the different light trajectories of the two images P 1 and P 2 are equal. It should be pointed out that the SIS model is an idealized lens model. The real strong lensing by a galaxy is likely not to be a SIS case. For a certain strong lensing event, if there are other observations of the properties of the lens object, we can calculate the value of α in a suitable lens model. Generally, both of Δt geo and Δt gra are of the order of G M/c 3 , where M is the mass of lens [22]. Therefore, one can expect that the parameter α should be of the order of one. GRBs are promising tools to constrain the WEP for photons with different energies [6]. Due to their large luminosities, GRBs can be observed at very high redshifts [23]. Therefore, there is a much greater possibility for a GRB to be lensed by a galaxy or galaxy cluster in the universe. Due to the success of BATSE, Swift and Fermi satellites, the number of detected GRBs keeps increasing. There is a lot of work researching the potential strong lensed GRBs in several GRB catalogs. However, no such event was ever found [24,25]. Li and Li (2014) researched the potential lensing events in BATSE GRB data. They found four candidates. The second couple of GRBs, 2044 and 2368, has properties closest to a strong lensed GRB event. We use this GRB as an example, although they excluded the possibility in their work [25]. The flux ratios in four considered energy channels seem similar. In addition, the angular separation of them is Δθ = 3.88 • while the location uncertainties of them are 2.88 • and 6.06 • . The detected time delay between 2044 and 2368 is about 1.77 × 10 7 s and the time delays of the photons in energy channels of 25-60 and 60-110 keV are 0.085 ± 0.042 and 1.730 ± 0.162 s for 2044 and 2368, respectively. Assuming there is a lensed GRB event with similar time delay, we constrain the WEP with our method. From Eq. (4), the constraint on the violation of WEP between the photons in these two energy bands is Δγ ≤ 2(1 + α) 1.730−0.085 1.77×10 7 = 1.86(1 + α) × 10 −7 . If we choose the SIS lens model, α = 0, we have Δγ ≤ 1.86 × 10 −7 . For the Fermi GBM, the expected time to observe one lensed GRB is about 11 years [25]. Since the Fermi GBM has served about 9 years and will be in service for more than another 10 years, it is reasonable to expect a lensed GRB event in the operating period of the Fermi GBM.

Discussion
In this section, we discuss two potential biases. The first point is the time delay potentially caused by the non-zero rest mass of photons Δt spec and the LIV Δt LIV . Generally, these two effects are similar since both of them will cause an energydependent speed of photons. If the photon has a non-zero rest mass, then the higher-energy photons will travel faster than lower-energy photons. On the contrary, LIV will lead to the opposite effect that higher-energy photons have lower traveling speed because of the so-called vacuum dispersion effect [26][27][28][29]. Therefore, the time delay terms Δt spec and Δt LIV are both caused by the potential difference of the traveling speeds of photons with different energies.
In Fig. 1, there are two traveling paths and for each path there are two images formed by two different kinds of particles. Let us assume that the speeds of those two different particles are v 1 and v 2 and the lengths of the two traveling paths are L 1 and L 2 , respectively. Therefore, the time delays caused by this effect are L 1 /v 1 − L 1 /v 2 and L 2 /v 1 − L 2 /v 2 , respectively. Their contribution to the difference of time delays |(t 22 − t 12 ) − (t 21 − t 11 )|, which we use in Eq. (4) to constrain the WEP, is ( . Considering a strong lensing system with a strong lensing time delay, which is about 1 year, and the time delays of different particles t 22 − t 12 and t 21 − t 11 are about 1 s, the difference of the time delays |(t 22 −t 12 )−(t 21 −t 11 )| should also be of the order of 1 s. For this kind of strong lensing system, the difference between L 1 and L 2 should be of the order of 1 light year. However, the distance of a typical cosmic source should be of the order of 1 Gpc, so it has Therefore, even though the time delay caused by non-zero rest mass of photons and LIV effect is hundreds of times larger than the observed time delay between two different particles, only contributes a very small part of |(t 22 − t 12 ) − (t 21 − t 11 )|. Therefore, it is reasonable to omit the terms Δt spec and Δt LIV , which means our method can also exclude the potential effects of a non-zero rest mass of photons and LIV.
The second point is that the two traveling paths will perhaps introduce some other contributions to the observed time delay. Actually, for a typical strong lensing system, the Einstein radius is about 4 × 10 −6 ( M 10 11 M ) 0.5 ( D 1Gpc ) −0.5 [20], which is very much smaller than the relativistic beaming angle θ ∼ −1 ∼ 10 −3 where is the Lorentz factor of the GRB's jet [30]. Therefore, the difference of the specialrelativistic boost factor caused by different viewing angles can be omitted reasonably. In addition, the different paths may also lead to different effects of gravitational potential along the paths, such as the weak lensing by the large scale structure and also the effects of our Galaxy and local galaxy cluster, which may also contribute into the observed time delay. However, from Eq. (1), the difference between the Shapiro time delays for two paths relies on the total difference of the gravity potential along the whole paths which includes the gravitational potential of large scale structure, the local galaxy cluster and also our galaxy. Therefore, the effect is considered even though there is some small difference of gravitational potentials along the two different paths.

Conclusions
In this paper, we have proposed a method to constrain the violation of the WEP with strongly lensed cosmic transients. Our method does not need to make any assumption as regards the physical mechanism of the transient. Moreover, because it utilizes the difference of time delays of the two lensing images, the potential effect of intrinsic time delay Δt int , LIV time delay Δt LIV and the non-zero rest mass time delay Δt spec can be naturally removed. By analyzing the properties of time delay of gravitational lens, we find that the parameter α, which represents the ratio of time delays between two lensing images caused by geometric and Shapiro time delay effects, will be zero for the SIS lens model. Therefore, the lengths of the light ray paths for two lensing images are the same in the SIS lens model, which means that any test depending on the difference of the path length cannot work in the SIS lens model, such as testing the difference between the speeds of light and GW. Even though for a realistic lens, one should consider the parameter α when constraining the speed of the GW. Otherwise an unreasonably tighter constraint will be obtained.
Due to the number of GRBs increasing, the detection of a strong lensed GRB is promising. However, there is still no such event at present. Assuming that a strong lensed GRB with time delay between images is about 1.77 × 10 7 s and the time delays of the photons in the two energy channels are 0.085 ± 0.042 and 1.730 ± 0.162 s, one can give a constraint on WEP at about the 10 −7 level. Besides GRBs, the lensed FRBs and the GW events with their EM counterparts are also potential candidates to test the WEP using our method [31][32][33]. Interestingly, lensing of FRBs has been proposed to probe dark matter [34,35].
From Eq. (4), with the time delay measurements of a single strong gravitational lensing event one can give a tight constraint on the Δγ . It can also be found that the efficiency of the constraint on WEP is proportional to the accuracy of the time delay measurement and inversely proportional to the strong lensing time delay. A typical strong lensing by a galaxy will give multiple images with days to months strong lensing time delays [36]. Because the strong lensing time delay is proportional to the mass of lens, the time delay due to galaxy cluster lens is much longer [37], which will give much stricter constraint on Δγ . If a several-months strong lensing time delay is observed and the time delay of two images is about 0.1 s, the constraint on the violation of WEP will be up to Δγ < 10 −8 . For the measurement of the time delay between two different particles, it depends on the accuracy of the timing measurement of the transient event. The accuracies of the time measurement of FRB, GW and GRB are about 0.01 ms [38], 10 −4 ms [31] and 0.1 s [39]. Recently, the detections of GW170817 and its electromagnetic counterparts [14,40] also encouraged us to find a more reliable constraint on the WEP using strong lensing of GWs and electromagnetic counterparts. Therefore, the accuracy of the WEP can be improved by several orders of magnitude in the future, if the lens is a galaxy cluster and the strongly lensed cosmic transients have a much more precise time delay measurement, such as the GW events and FRBs.