Effect of electric interaction on the deflection and gravitational lensing in the strong field limit

The deflection angle $\Delta\phi$ of charged signals in general charged spacetime in the strong field limit is analyzed in this work using a perturbative method generalized from the neutral signal case. The solved $\Delta\phi$ naturally contains the finite distance effect and takes a quasi-power series form with a logarithmic divergence at the leading order. The coefficients of the series contain both the gravitational and electric contributions. Using the Reissner-Nordstr\"{o}m spacetime as an example, we found that an electric repulsion (or attraction) tends to decrease (or increase) the critical impact parameter $b_c$. If the repulsion is strong enough, then $b_c$ can shrink to zero and the critical particle sphere $r_{0c}$ will disappear. These results are applied to the gravitational lensing of charge signal, from which we solved the image positions, their magnifications and time delays. It is found that in general, the electric repulsion (or attraction) will decrease (or increase) the image apparent angles, the black hole shadow size as well as their magnifications but increase (or decrease) the time delay.


I. INTRODUCTION
Deflection of light in gravitational field was one of the most important evidence for the establishment of General Relativity (GR) [1].Based on this, the gravitational lensing (GL) phenomenon has been developed into a powerful tool in astronomy.GL can not only link the properties of the source, the lens and the messengers to the observables [2][3][4][5][6][7][8][9], but also be used to test gravitational theories beyond GR [10,11].
With the discovery of supernova neutrinos [12,13] and blazer neutrinos [14,15], the gravitational wave [16][17][18][19][20], and the more historical cosmic rays [21], people become more and more interested in the deflection and GL of massive signals.Among these, cosmic rays with ultrahigh energies are known to be composed of charged massive particles of protons and heavier nuclei [22].Therefore, not only are they massive but carry charge and experience electric interaction if the spacetime is also charged.
Previously, there have been a large amount of studies on the trajectory of charged signals in electrically or magnetically charged spacetimes, including the Reissner-Nordström (RN) spacetime in the former case [23][24][25][26], and some particular spacetime background supplemented by a weak magnetic field in the latter [27][28][29].However, most of these are either concerned with the circular motion of the signal or the general features of the motion.Concentrating on the deflection and GL of charged signals, then they have been studied in the weak field limit in arbitrary static and spherically symmetric (SSS) spacetimes using perturbative method [30] and in particular spherically or axially symmetric spacetimes using Gauss-Bonnet theorem method [31][32][33][34].It was shown that in general, the repulsion (or attraction) between small signal and lens charges will increase (or decrease) the deflection angle at the leading order, i.e., at the order O(M/b) where M and b are the mass of the lens and impact parameter of the signal respectively.For a gravitational-electric dual lensing in the weak field limit, this implies that in the repulsion (or attraction) case the impact parameters, the apparent angles, the magnifications and total travel times of the images will be smaller (or larger) [30].
However, these results are more or less intuitively expected because in the weak field limit, gravity resembles its Newtonian limit and therefore both the gravitational and electric potentials take the inverse square form.Consequently the physics due to both interactions can be expected in the same way, except a possible choice of signs of the electric interaction.In contrast, in the strong field limit (SFL), GR deviates most from its Newtonian limit but the Coulomb potential keeps its form.Then one would naturally anticipate more interesting interplay or competition between these two kinds of interactions and more nontrivial results are expected.
The signal deflection and GL in the SFL have been studied for a relatively long time.Virbhadra et al. first pointed out the existence of the relativistic image sequences [35,36].Bozza et al. systematically studied the GL in the SFL for null signals in SSS spacetimes and in equatorial plane of stationary and axially symmetric spacetimes [37][38][39], and then many authors followed in particular spacetimes [40][41][42][43][44][45][46][47][48][49][50].In Ref. [51,52] we were able to develop a perturbative method applicable to arbitrary SSS spacetime to study the deflection, GL and time delay for neutral timelike signals.It was found that in the SFL, i.e. as the impact parameter b approaches its critical value b c , the deflection takes a simple quasi-power series form and the time delay has a simple physical interpretation.
In this work, we would like to examine how the electric interaction will affect these results, both qualitatively and quantitatively.We will show that our method not only works for charged signals in arbitrary SSS charged spacetime, but it can automatically take into account the finite distance effect of the source and observer.The paper is organized as follows.In Sec.II we lay out the general equations involved in the problem.In Sec.III, the perturbative method to compute the deflection angle and total travel time with electric interaction is presented.In Sec.IV, we take the RN spacetime as an example to show how the formalism works in a particular spacetime.The corresponding GL and black hole (BH) shadow are analyzed in Sec.V, where the images' apparent angles, their magnifications, time delays as well as the BH shadow sizes are solved.Sec.VI closed the work by a discussion.Throughout the paper we work in the units where G = c = 1/(4πε 0 ) = 1.

II. METRIC, MOTION EQUATION AND PARTICLE SPHERE WITH ELECTRIC INTERACTION
The most general SSS metric can be described by the line element ds 2 = −A(r)dt 2 + B(r)dr 2 + C(r)(dθ 2 + sin 2 θdφ 2 ), (2) where t, r, θ, φ are coordinates and A(r), B(r) and C(r) are functions of the radial coordinate only.
When the signal is charged and the electromagnetic force is taken into account, the signal will not follow the geodesic but the Lorentz equation [53] where q and m are the charge and mass of the signal and τ is the proper time.The electromagnetic field strength F µν is given by where A µ is the four potential of the electromagnetic field.For simplicity, in this work we will only consider the effect of an electric field but not the magnetic field.
That is, we assume A 0 (r) = 0 and A i = 0 (i = 1, 2, 3).Without losing any generality, for such a SSS spacetime and electric potential we can always assume that the trajectory lies in the equatorial plane so that θ(τ ) = π/2, and then Eq. ( 3) after the first integrals reduces to three separate equations FIG. 1: The trajectory deflection and GL in the SFL.
The source and detector are located at coordinates (r s , φ s ) and (r d , φ d ) respectively.The apparent angle in this example is where E and L are integration constants that can be interpreted as the energy and the angular momentum of the particle respectively.Comparing to the case without the electric interaction (setting q to zero), we observe that the changes in these equations happen in and only in terms containing E: all E is replaced by the E + qA 0 (r), i.e., the electric potential energy manifests itself in this way.
In asymptotically flat spacetimes, on which this work will concentrate, L and E can be related to the impact parameter b and asymptotic velocity v of the signal Moreover, L can also be related to the shortest approach r 0 of the trajectory, which is defined through ṙ| r=r0 = 0. Using Eq. (5c), the equation determining r 0 becomes Eqs. (6a) and ( 7) establish the following correspondence between b and r 0 , and this allows us in principle to express the deflection angle and time delay in either of them.For a signal from a source located at radius r s to an observer at r d (see Fig. 1), from Eqs. (5) For general SSS spacetime metric and potential A 0 (r), these integrals usually can not be worked out analytically even for the simplest RN case.Therefore, some special techniques are required to approximate them.The technique we will use in this work is actually a generalization of the perturbative method developed in Ref. [51,52] for neutral signals to the current case, in which the electric interaction is now taken into account.
In the SFL, the spacetime allows the existence of some critical value r 0c of the closest approach r 0 .When r 0 → r + 0c , the trajectory will circulate many times around the center before propagating to the observer; while below it, the signal will enter the event horizon and therefore not reach the observer anymore.Previously, the existence of such r 0c has been known for many BH spacetimes for photons (the photon sphere) [54], and for massive signals too (the particle sphere (PS)) [51,55].In this work, because of the presence of the electric interaction, the existence and size of r 0c might both be affected.Therefore, next, we derive the equation determining r 0c in this new case first.This process will also help us to better understand its dependence on various parameters of the spacetime and signal.
We first re-write Eq. (5c) into the following form Then we can define the second term on the left hand side as the effective potential When A = 1/B and A 0 = 0, as in many well-known SSS chargeless spacetimes, this reduces to the effective potential in the corresponding spacetimes [55].The terms proportional to q 2 A 2 0 and qEA 0 in Eq. ( 11) are due to the electric interaction.From Eq. (10) we see that the allowed range of r for a signal is when V eff (r) is smaller than the right hand side of Eq. (10), i.e., The equal sign of this equation is actually equivalent to Eq. ( 7) and therefore also determines the closest radius Now for many interesting spacetimes, their Eq.( 10) also possesses a critical behavior.As certain parameter(s) in V eff (r) and E 2 /(2AB) varies (e.g., L or equivalently b, q and Q, etc), if these two quantities become barely intersecting and eventually detach from each other (see Fig. 2), then at this point, the r 0 , which is defined by the equal sign of Eq. ( 12), ceases to exist and the signal coming from larger r will be able to continue to smaller r.This critical behavior can be thought to happen at the critical value of the above mentioned parameters.But on the other hand, since this detaching point in the radial direction is also the location of r 0 , it can also be thought as the critical value of r 0 , which we will denote as r 0c .Mathematically, r 0c can be determined from the conditions And we can also supplement the condition for r > r 0c (14) to ensure the signal can come from large r.Once the spacetime metric is specified, Eq. ( 13) will fix the value of r 0c completely.Using Eq. ( 8), the corresponding critical impact parameter b c becomes

III. PERTURBATIVE EXPANSION IN THE SFL
To compute the deflection angle and total travel time defined in Eq. ( 9) in the SFL, in the following we will extend a perturbative method developed previously in Ref. [51,52] to the case with electric interaction.It turns out that the methodology and general form of ∆φ and ∆t will also work in this case after some modification of the critical r 0c and the expansion coefficients in Eq. (40).
In this method, we first define a function p(x) inspired by Eq. ( 15) (16) Using Eq. ( 8), it is clear then Denoting the inverse function of p(x) as w(x), the above means Using the function w(x), we can define a change of variables in the integrals (9) from r to ξ, which are connected by With this, the integral limits and integrands of Eq. ( 9) change in the following way where and in Eq. (20b) we have defined η s,d , and in Eq. (20e) and henceforth we set Collecting terms in Eq. ( 20) together and using Eq. ( 6) and ( 22), Eq. ( 9) becomes In the SFL, i.e., r 0 → r + 0c , b → b + c and a → 0 + , the integrands of Eq. ( 23) can be perturbatively expanded for small ξ and therefore it allows us to find an effective approximation of ∆φ and ∆t.Among the four factors of each of the integrands, the last one is the simplest and can be kept since its product with a power series of ξ can be directly integrated.The third factor has a simple Taylor expansion Denoting the product of the first two factors of Eq. (23a) as f (ξ) and those of Eq. (23b) as g(ξ), we have we then are able to show that f (ξ) and g(ξ) have the following form of expansions Note in these expansions, the series expansion of w(x) is needed, and it can always be obtained using the Lagrange Inversion Theorem from the series form of its inverse function p(x).In Eq. ( 27), the initial summation index is -1 and the powers are half integers because q (x) has a singularity at x = 0. We can indeed collectively denote the product of the first three factors in Eq. (23a) as y(ξ, a) and that of Eq. (23b) as z(ξ, a) Then their expansions are given by the product of expansions ( 24) and (27a), and ( 24) and (27b).After some re-organization, they become respectively where y n,m and z n,m denote coefficients of corresponding powers of ξ and a.Note that these y n,m and z n,m can be completely determined once the metric functions are fixed.The first few of them for general SSS spacetime with a critical r 0c are given in Eqs. ( 40) and (43).Substituting Eq. ( 29) into Eq.( 23), ∆φ and ∆t are written as series of integrals of We emphasize that the integrals in the above series can always be carried out and the results are some elementary functions presented in Eqs.(A1) of Ref. [51].Using these formulas, ∆φ and ∆t are finally found to be A few remarks are in order here.First of all, for both ∆φ and ∆t, we observe that when r s and r d are not infinite, there is only one divergence proportional to ln a in the SFL a → 0 + , which is only contained by the k = m = 0 term.Secondly, we can further expand other functions involving a, i.e. ln( √ η i + √ η i − a) and 1 − a/η i etc., in the limit a → 0 + , and the results should become quasi-power series of a with the coefficient of a n containing one ln a.
where C n and D n , and C n and D n are constants determined by the initial conditions of the trajectory and metric function parameters.Thirdly, from Eq.(31a), it is seen that to the O(a) 0 order, the m = 0, j = k terms contribute dominantly to ∆φ so that it becomes where we have identified The travel time (31b) can also be similarly expanded again for small a so that to the leading orders it becomes where we see that At this point, it is appropriate to point out the difference between the perturbative method and result in this work and the case without electric interaction [51,52].Superficially, the steps from Eq. ( 17) to the result (38) are very similar to the case without electric interaction.However, indeed there exist a fundamental difference in the definition of p(x) in Eq. ( 15), where the electric potential explicitly contributes in the current case but not so in Ref. [51,52].Consequently, its inverse function w(x) is also changed by A 0 (r) and this further changed all formulas involving the expansions of terms containing w(ξ/b c ).More specifically, the coefficients f n and g n in Eq. ( 27) and consequently y n,m and z n,m in Eqs. ( 29)- (31), as well as the C n , D n , C n , D n in Eqs. ( 32)-( 38), all contain contribution from the electric potential.We now show their dependence on A 0 (r), as well as other metric functions, more explicitly.
Assuming the metric functions and the electric poten-tial have the following series expansions near r 0c where a n , b n , c n and a 0n (n = 0, 1, • • • ) are the coefficients, and substituting them into Eqs.( 16), we can solve the series p(x).Inverting it to find w(x) and then further substituting into Eq.( 28), we can compute y n,m and z n,m to any desired order.The first two of y n,m , i.e. y −1,0 in Eq. ( 34) which fixes the divergence of ∆φ and y 0,0 in Eq. ( 35) which dominates the constant term of ∆φ, are where b c is given in Eq. ( 15) and T 2 and T 3 are related to the first few coefficients in Eq. ( 39) by Similarly, the first few z n,m are found to be The higher order terms of y n,m and z n,m can also be computed but are too long to present here.It is seen from these equations that a 0n explicitly appear in all of the y n,m and z n,m and therefore will affect how the deflection angle ∆φ and total travel time ∆t diverge.In Sec.
IV, we will study the dependence of ∆φ on the electric interaction using RN spacetime as an example.It is also interesting to note from Eqs. ( 41) and ( 42) by factoring out an E 2 from the curly brackets that both T 2 , T 3 , and indeed all orders of y n,m and z n,m and consequently the entire ∆φ and ∆t, depend on (q, m, E) only through the ratios q/E and v = √ E 2 − m 2 /E.That is, there is one degree of degeneracy among these three parameters.This fact was also observed in the weak deflection case in Ref. [30,34].

IV. THE RN SPACETIME CASE
In this section, We apply our result to the RN spacetime, which is the simplest SSS spacetime that allows both electric and gravitational interactions.The metric functions and the electric potential in the RN spacetime are given by We will examine the validity of the results in Secs.II and III, particularly concentrating on how the electric interaction will affect the PS radius r 0c defined in Eq. ( 13), b c defined in (15) and the deflection ∆φ in Eq. (33).
A. The critical r0c and bc Substituting Eq. ( 44) into Eq.( 13) and simplifying the result, it is not difficult to find that r 0c satisfies the following quartic equation Denoting the coefficients of r n 0c by g n for n = 0, • • • , 4, this equation can be solved to yield an explicit formula of the only physical root of r 0c (see Eq. (A3) of Ref. [56] and [57]) where , Although Eq. ( 46) is not very transparent to see the effects of various parameters such as M, Q, E, m and q on the above r 0c , it is however not difficult to demonstrate by dividing Eq. ( 45) by m 2 M 4 that r 0c /M depends on them only through the ratios q/m ≡ q, Q/M ≡ Q and In other words, Moreover, there are a few limits of Eq. ( 46) that one can check.The first is its neutral particle limit, which can be obtained by letting q → 0 and then the result agrees with Eq. ( 18) of Ref. [56].The second is its Schwarzschild limit reached by letting Q → 0 This agrees with Eq. ( 18) of Ref. [58].The third limit is the ultra-relativistic limit of the charged signal, which is simply the limit E → ∞ while holding m finite or equivalently v → 1.The result is Finally, one can also attempt the v → 0 limit of r 0c , in which case the Eq. ( 45) becomes cubic and a root simpler than Eq. ( 46) can be obtained (see Eq. ( 26) of Ref. [56] for the neutral case).
In both the limits ( 48) and ( 49), the dependence of r 0c on the charge-mass ratio q of the charged signal disappears.The disappearance in the former is quite easy to understand, since in a neutral spacetime like the Schwarzschild one, the signal although charged will not experience any electric interaction.While for the second disappearance, one can show by dividing Eq. ( 45) by E 2 M 4 that r 0c will depend on q only through q/(E/m) = q√ 1 − v 2 .That is, Eq. ( 47) can be further transformed into Therefore as v → 1, q effectively drops out from r 0c .The r 0c solved from Eq. ( 45) is plotted in Fig. 3 (a) as a function of q/m and v for a typical Q/M .It is seen that for any fixed v ≤ 1 and 0 < Q/M ≤ 1, as q increases to some positive maximal value q max , r 0c decreases monotonically if v is large or otherwise it could increase again when q increases towards q max .Only for q ≤ q max , a critical r 0c is allowed.This implies that for any given Q (|Q| ≤ M ), as the repulsion between q and Q increases, the critical r 0c may cease to exist.For cases with q > q max , signals incoming along any direction will not be captured.This q max is actually determined by the requirement that even when L = 0 = b, there still exists a closest approach of radius and its value should be equal to r 0c .Solving r 0 from L = 0 using Eq. ( 7), the positiveness of r 0 yields the q max for the RN spacetime to be and at this boundary, the value or r 0c is simply Note that since v < 1 this value is always larger than the RN BH exterior horizon radius r The boundary (51) and its corresponding value of r 0c is also plotted in Fig. 3 (a) using the red curve.We also note that in the limit v → 0, q max approaches mM/Q and r 0c (q, v → 0) approaches an infinite value, except in the extremal RN case r 0c (q = q max , Q = M, v) = M .For any r 0c , there exists a corresponding critical value of the impact parameter, b c .For RN spacetime, using Eq. ( 15), b c is found to be This b c is also important because it can be directly linked to the angular size θ sh of the BH shadow observed by the remote observer at radius r d , through relation θ sh ≈ b c /r d (see Eq. ( 63)) when E and r d are large.Using Eq. ( 53), the b c corresponding to Fig. 3 (a) is plotted in Fig. 3 (b).It is seen that unlike r 0c , b c for any fixed Q and v decreases monotonically as q increases, and reaches 0 at q max .While for its dependence on v, then as v decreases to zero, b c increases monotonically to infinity.
It is more interesting to study the effect of Q on b c than that of q/m, because q/m only affects the electric deflection while Q affects both the electric and the gravitational deflections.Previous works have implied that in the weak field limit increasing |Q| would decrease b c gravitationally for neutral signals (see Eq. (122) of Ref. [56]), and the desired b to reach the same observer will be increased if the signal is charged and qQ < 0 (see by Eq. (3.3) of Ref. [30]).These mean that when q is fixed and qQ < 0, a nonzero Q will cause a competition between its gravitational and electric effects on b c , which is worthy to investigate.FIG.4: b c as a function of Q/M and q/m for v = 0.9c.The red curve is the maximal Q max given in Eq. (54) for a fixed q.
In Fig. 4 we plotted the b c as a function of Q/M and q/m.It is seen that for fixed q/m, b c monotonically decreases as Q increases, as long as it is positive so that qQ > 0. For q/m 2.3, when Q increases to a certain point Q max , the b c decreases to zero.The critical points where b c approaches zero is also where r 0c ceases to exist and therefore we can work out this Q max directly from Eq. ( 51).Replacing q max by q in this equation, we find which is shown as the red curve in Fig. 4. Using relation (65), b c approaching zero corresponds to the point that the BH shadow angular size shrinks to zero.This will be verified in Fig. 6.While for small q/m (0 < q/m < 2.3), b c only decreases to a finite value when Q reaches its extreme value of Q ext = M .On the other hand, for Q < 0 so that qQ < 0, we also see from Fig. 4 that b c also increases as Q decreases by a small amount from zero (more apparent for large q).This suggests that for the chosen ranges of parameters q and v, the electric effect on b c is stronger than that of the gravitational one for small negative Q, because we know that Q's gravitational effect alone would decrease b c [30].With Q further decreases to more negative value, we see that for small q/m, b c decreases again, implying the gravitational effect exceeds the electric effect again.While for larger q/m, b c keeps increasing for the entire range of Q from 0 to −Q ext .This will also be confirmed in Fig. 6 by ultra-high energy protons.The above features for small and large negative Q are intuitively understandable after a quantitative comparison of the two effects on ∆φ in the weak field limit.The gravitational and electric contributions of Q to ∆φ in this limit are respectively about the size where c and c are some positive constants (see Eq. (4.4) of Ref. [30]).Consequently for a small but fixed q, when Q is very small (or large) the gravitational effect on the deflection is always weaker (or stronger) than the electric one.Also from these contributions in Eq. ( 55), we see that if q is very large, such as in the case of q/m 2.3 in Fig. 4, the charge Q might not be able to reach the desired value for the gravitational effect to overcome the electric one, before it reaches its extreme value −M .Finally, it is also instructive to comment on the case when q exceeds the q max in Eq. ( 51) for a fixed Q or Q exceeds the Q max in Eq. ( 54) for a fixed q.The vanishing of r 0c or b c reaching zero do not mean that the signal will be able to travel to the very central region of the BH.Indeed in these cases the closet approach r 0 of the trajectory solvable from Eq. ( 7) still exists and can be shown to be well beyond the RN BH outer horizon.However, the critical behavior of the effective potential (11) detaching the right hand side of Eq. ( 12) will not happen.A key noticeable feature following this is that the charged signal will always only experience a finite amount of angular deflection regardless how small is b.

B. Deflection angle ∆φ
After obtaining r 0c , we can continue to expand the RN metric and electric potential functions in Eq. ( 44) and obtain Reading off the coefficients a n , b n , c n and a 0n in Eq. ( 39) from the above, and substituting them into Eqs.( 40)-( 42), the y −1,0 , y 0,0 and higher order y n,m 's can be obtained.Further substituting into Eq.( 33), the perturbative ∆φ in the RN spacetime in the SFL for charged signal is found immediately where the coefficients are Higher order y n,0 were also computed but they are numerically less important and too long to be shown here.Other parameters are the same as in (a).
To check the validity of this ∆φ, we can define a truncated ∆φ m (m = 0, 2, • • • ) by including only terms to order m in the summation of n in Eq. ( 57) and compare its value with a ∆φ num obtained by numerically integrating Eq. (9a).In Fig. 5 (a) we chose some exemplary parameters M = 1, Q = M/2, q/m = 1, v = 9/10 and plot the ∆φ m and ∆φ num as functions of a = 1 − b c /b.We see that as the truncation order increases, the perturbative result converges to the numerical value in general, and more rapidly so for smaller a.This is in accord with the fact that as a → 0, ∆φ diverges as ∼ ln a and the a-independent terms composing D 0 in Eq. ( 35) become less important to ∆φ.
To study the effect of the charges Q and q and the electric interaction on the deflection in the SFL, in Fig. 5 (b) we plot ∆φ using Eq.(57) for a fixed impact parameter.We chose a small parameter range of Q ∈ (0.5±0.0005)M and q ∈ (1 ± 0.002)m in order to make sure that for the fixed b, the parameter a = 1 − b c /b is still smaller than 1.4×10 −4 so that the approximation (57) is always valid.It is seen that as Q or q increases, ∆φ decreases monotonically in this range, suggesting that a stronger repulsion results in a smaller deflection.This actually is understandable from Fig. 4 that larger Q and q when qQ > 0 result in a smaller b c and therefore the trajectory with fixed b experiences weaker field and smaller deflection.

V. GL IN THE SFL AND BH SHADOW BY CHARGED SIGNALS
To reveal the effect of the electric interaction on the GL in the SFL, we will have to solve the GL equation for this case.Formally, this equation and its solution process are the same as in the case of neutral particle which was studied in Ref. [51], except now all parameters or coefficients have to be updated to include the electric interaction.Therefore in this work, we will do a parallel analysis of the GL but concentrate on the electric effect this time.
When one has a ∆φ that takes the finite distance effect of the source and detector into account, establishing the GL is particularly simple.For a source located at (r s , φ s ) and detector at (r d , φ d ) (0 ≤ φ d < φ s < 2π) (see Fig. 1), the change of the angular coordinate ∆φ satisfies where integer |n| = 1, 2, • • • is the looping number of the trajectory around the center and n > 0 and n < 0 corresponds to the anti-clockwise and the clockwise looping directions respectively.Substituting Eq. ( 33) for ∆φ into Eq.( 60), it becomes where we have denoted ∆φ sd ≡ φ s − φ d .From this, we can easily solve b that allows the signal to reach the detector in terms of other quantities Formally this is the same as Eq. ( 38) of Ref. [51], however the parameters b c , C 0 and D 0 are now given by the updated values in Eqs. ( 15), ( 34) and (35).
A. The relativistic image and BH shadow locations Corresponding to this series of b n are two series of images in the SFL, one series on each side of the lens, and their apparent angles are given by the following formula [30] Substituting Eqs. ( 62) and (15), this becomes Clearly if all other parameters are fixed but |n| increases, θ n will decrease monotonically.At |n| → ∞, this yields the angular size θ sh of the BH shadow formed by charged signals.That is, Clearly, the effect of first term inside the arcsin function in Eq. ( 64), including parameters sign(n) and ∆ sd , becomes irrelevant to the shadow size.
To see the effect of the electric interaction on θ n and θ sh , in Fig. 6 we plot them for the SgrA* SMBH by assuming it is a RN BH and the shadow is formed by cosmic protons.We take the proton energy to be 10 19  [eV] and the SgrA* SMBH mass M = 4.30 × 10 6 M and source/detector distance r s = r d = 8.28 [kpc] to be the distance to us [59].For the BH charge, we use the magnetically induced charge as the rough unit [30].Here we take the typical order of B mag ≈10 [G].We plot the first relativistic image θ 1 and then the BH shadow θ sh for several Q.We chose ∆φ sd = π so that the source, lens and detector are aligned.From Fig. 6 (a) it is clear that for our fixed q, as Q increases positively from zero (the red solid curve), the shadow size decreases monotonically, until Q reaches Q max at which point b c and consequently θ sh becomes zero.This corresponds to the large q case in Fig. 4, for which a maximal Q exists.For the chosen parameter values (v, q/m, M etc.) in this figure, we can actually work out the critical Q max = 5862Q mag from Eq. ( 54).Also, in this case, one can see from the green dash-dotted curve that for negative Q, the shadow size will increase, which again is in accord with Fig. 4. Because of the same electrical repulsion effect, we see from Fig. 6 (b) that if charge Q of the BH is fixed to a positive value, decreasing q/m to its value of typical heavier nuclei, i.e., half its value for proton, will increase the shadow size.From both zoom-in's in Fig. 6 we see that all the relativistic images are very packed, and close to the BH shadow location.This is a general feature of all such lensed images in the SFL.

B. Magnification
The magnification of the images are defined as where β is the angle of source against the lens-detector axis.To carry out the differentiation in above, we have to link β to quantities in Eq. ( 64).This can be achieved by considering the following geometrical relation Formally, this agrees with Eq. ( 44) of Ref. [51] (after a typo corrected there), but again the coefficients C 0 , D 0 and b c are the ones containing the electric interaction.
To see the electric effect on µ n more clearly, in Fig. 7 we have plotted the magnification as a function of Q and q.It is seen that comparing to neutral particles, the magnification receives a correction that is also determined by sign(qQ).When sign(qQ) is positive (or negative), the magnification is decreased (or increased).Therefore qualitatively the effect of the electric interaction to the magnification here in the SFL is the same as that in the weak field limit [30].

C. Time delay
With the impact parameter b n of the n-th images known in Eq. ( 62), we can find the total travel time FIG.6: The Einstein ring θ 1 (∆φ sd = π) in the SFL using Eq. ( 64) and BH shadow θ sh using Eq.(65) for RN spacetime: (a) Formed by proton with energy 10 19  [eV] for several Q.From outer to inner: Formed by protons with energy 10 19 [eV] and heavier nuclei with the same velocity and Q = 10 3 Q mag .From inner to outer: proton and heavier nuclei with half q/m of proton.Each line actually contains two separate lines: the outer θ 1 and inner θ sh respectively (see the zoom-in).
for the n-th trajectory using Eq. ( 36) FIG. 7: The magnification µ n using Eq. ( 69) for n = 1, ∆φ sd = 5π/6.q/e is the ratio of particle charge to elementary charge from −1 to 1.Other parameters are the same as in Fig. 6.
Using this, we can obtain the more observationally interesting quantity, the time delay ∆ 2 t n,m between the n-th and m-th images where Eq. ( 62) for b n was substituted.
We now show that similar to the case of neutral particles [52], Eq. ( 71) also has a simple interpretation as the amount of local time to move around one full circle of the PS times the redshift factor from the PS to the observer, and then times the difference of the looping numbers of the two trajectories.For the first factor of Eq. (71), using Eqs.( 35), (37), (40a) and (43a), it becomes Now for a SSS spacetime, it is always possible to scale the metric function C(r) = r 2 and therefore its expansion at r 0c yields c 0 = r 2 0c .Consequently the numerator 2π √ c 0 of Eq. ( 72) is actually the circumference of the PS.For the first denominator of this equation, we can show that it is nothing but the local speed v l of the particle's motion around the PS.Since in the SFL the signal circulates the PS, v l is given by where γ sr = 1/ 1 − v 2 l is the gamma factor due to special relativity.Substituting Eqs.(5b), (6a), (15), (39) and c 0 = r 2 0c into Eq.( 73) and after some simple algebra, one can solve v l as which is exactly the first denominator of Eq. ( 72).Finally, the factor 1/ √ a 0 = 1/ A(r 0c ) in Eq. ( 72) is the redshift factor from the PS to an asymptotic observer in an SSS spacetime described by the metric (2).Combining these three parts, therefore the first factor of Eq. ( 71) is the time observed by the asymptotic observer for the signal to loop one circle of the PS.The second factor of Eq. ( 71) is the difference between looping number of the n-th and m-th trajectories.Putting these factors together, we obtain the claimed interpretation of the formula (71).That is, In the above time delay, usually E, m and q of the signal are measurable.And if we assume that the spacetime is RN type, then it is seen from Eqs. ( 45), (56a) and (56d) that r 0c , a 0 and a 00 will be fixed once M and Q are known.Since mass M of a lens can often be known using other astronomical means, from Eq. ( 75) we immediately see that the charge Q can be constrained by the measurement of ∆ 2 t n,m .In Fig. 8, we plot the PS radius r 0c , the redshift factor 1/ √ a 0 and the time delay ∆ 2 t m+1,m as functions of Q by again assuming the Sgr A* SMBH is a RN one.Note that ∆ 2 t m+1,m actually does not depend on the value of ∆φ sd or m as long as m ≥ 1 or m ≤ −2.
It is seen from the lower subplot of Fig. 8 (a) that for charged signal with q > 0, the time delay increases when Q deviates from zero.This is actually a combine effect of a decreasing PS radius and increasing redshift factor (see the upper subplot of Fig. 8 (a)).While for neutral signals, although all quantities seem flat in Fig. 8 (a), they are not truly invariant.They appear quite flat because in this case there is no electric interaction and therefore Q only affect them gravitationally and then it takes a much larger Q to make a comparable difference.
As one can see from Fig. 8 (b), only when Q reaches about 0.25M ≈ 7.6 × 10 10 Q mag , the time delay variation can reach the same order as the charged case in Fig. 8 (a).This comparison suggests that although the time delay of the neutral signal can also be used to constrain the lens charge as suggested in Ref. [52], that of the charged signal is more sensitive when Q is small, which is generally expected by many astronomers.

VI. CONCLUSIONS
In this work we studied the deflection, GL and time delay of charged signal in a charged spacetime in the SFL using a perturbative method.The effect of the extra electric interaction was carefully analyzed.It is found that as the trajectory approaches the PS.Perturbatively, a small electric repulsion (or attraction) will decrease (or increase) the critical impact parameter and therefore decreases (or increases) the deflection.For any fixed Q, in general there is a maximal q max given in Eq. ( 51) (qQ > 0) beyond which the PS ceases to exist and b c shrinks to zero.While for fixed q, the electric effect of Q on b c will be weaker (or stronger) than its gravitational one when Q is small (or large).
For the GL in the SFL, similar to the neutral signal case, there exists one series of very packed and weakly magnified images on each side of the lens.Their angular positions and magnifications are given by Eqs. ( 64) and (69).The BH shadow with electric interaction is still given by a simple formula Eq. (65).The electric repulsion (or attraction) tends to decrease (or increase) the angular sizes of these images and that of the BH shadow.For a fixed q/m of the charged signal, the BH shadow size shrinks to zero when Q is beyond Q max (qQ > 0).
The time delay between two images given by Eq. ( 75) has an intuitive and yet quantitatively precise explanation, as the PS circumference dividing the local velocity and then multiplying the redshift factor from the PS to the observer and the difference of the loop numbers of the two trajectories.It is shown that comparing to neutral signals, the time delay of charged signal is much more sensitive to the spacetime charge Q when it is small.
We emphasize that most of the results in this work, except Sec.IV, are completely general to all SSS charged spacetimes.Therefore besides the RN one, it is equally simple to apply them to other kinds of charged space-times, such as the Gibbons-Maeda-Garfinkle-Horowitz-Strominger (GMGHS) [61][62][63] and charged Horndeski [64] spacetimes.We also expect that it is straightforward to extend the perturbative method in this work to the deflection and GL in the equatorial plane of stationary and axially symmetric spacetimes, with or without the electric interactions.
we see then the change of the angular coordinate ∆φ and the total travel time ∆t become respectively ∆φ =

FIG. 2 :
FIG. 2: The effective potential V eff (r) and E 2 /(2AB) as functions of r for different values of parameters in them.(a) The case with a regular trajectory bouncing back at r 0 .(b) The case that the trajectory experiencing the critical r 0c .(c) The supercritical trajectory case, for which the r 0 does not exist anymore.

FIG. 3 :
FIG. 3: (a) and (b) are respectively r 0c and b c as functions of q/m and v for Q/M = 1/2.

FIG. 5 :
FIG. 5: (a) Comparison between the truncated ∆φ m for m = 0, 2, 4, 6 and the numerical integration result for a = 1 − b c /b.(b) ∆φ as a function of Q/M and q/m.Other parameters are the same as in (a).

FIG. 8 :
FIG. 8: (a): Upper: the PS radius (solid curves) and the redshift factor (dashed curves) as functions of Q for protons (red curves) with energy 10 19 [eV] and neutral signal (blue curves) with Sgr A* as the lens.Lower: the corresponding time delay |∆ 2 t m,m+1 |. (b): the case for neutral signal for a much larger range of Q.