Null Gravitational Redshift by a Reissner-Nordstr\"{o}m Black Hole in the Strong Field Limit

The gravitational shift of electromagnetic frequency in the strong field limit is usually investigated under the common scenario, where the light receiver is far away from the central body while the emitter is in the strong-field region of the lens. In this paper, the gravitational frequency shift of light caused by a Reissner-Nordstr\"{o}m (RN) black hole is studied numerically in the traditional strong-field scenario, as well as in the scenario where both the light emission and reception events happen in the strong-field region of the black hole. In order to obtain the numerical results of the gravitational redshift, we first derive the exact null equations of motion in the RN geometry in harmonic coordinates. For a given light observer, a new numerical technique is proposed in the integration of the geodesic equations to determine the spatial position of the emitter, considering the fact that their spatial positions are not always known simultaneously. Our work might be helpful to the related observations for probing strong gravity.


Introduction
The study on gravitational frequency shift of light signals is a important topic in modern astronomy since its original prediction by general relativity (GR) [1]. Acting as a crucial experimental verification of the local position invariance [2] of Einstein equivalence principle, null gravitational shift (especially redshift) has been applied extensively from groundbased [3][4][5][6][7][8] to space-borne-based [9][10][11] measurements and has received its high-accuracy consistency (at the 7 × 10 −5 level) with GR's prediction. In contrast to these efforts devoted to the weak-field counterparts in the solar system, investigations of gravitational shift of spectrum lines produced a e-mail: hgs@usc.edu.cn (corresponding author) b e-mail: lilinmath@usc.edu.cn in the strong-field region of a compact object (e.g., a neutron star or black hole) beyond the solar system were also performed with a fascinating prospect (see [12][13][14][15][16][17][18][19][20], and references therein). In particular, Müller and Wold [17] adopted the ray-tracing technique to numerically investigate the gravitational redshift effect of emission lines from both the weakand strong-field regions of a Kerr black hole, which was compared with the observations of galaxy Mrk 110 [21]. Actually, probing strong gravitational fields using the gravitational shift tests along with other tests (e.g., Refs. [22][23][24][25][26][27][28][29][30][31]) plays a significant role in the examinations of general relativity as well as alternative theories of gravity [2,32,33]. Two reasons are responsible for this. The first one is that strong-field regions in our Universe are common, however, most of actual tests have been restricted in our solar system [33]. It is necessary and naturally desired to grasp the degree to which the consistency of the predictions of gravitational theories including GR with the astronomical observations can achieve. A second one is that gravitational shift tests are classical and relative feasible, and they have become one of the most convenient chooses [16]. Recently, the most direct evidence of a binary black hole system in the first direct detections of gravitational waves was reported [34]. Very recently, the images of a supermassive black hole (i.e., the one at the center of M87 Galaxy) observed directly by the Event Horizon Telescope have also been presented for the first time [35,36]. There is no doubt that these two events as well as other observations (e.g, Ref. [37]) may promote further the theoretical studies of gravitational shifts of frequency of light in the strong field of a black hole.
Historically, considerations of gravitational shift in the strong field limit were usually done under the common scenario, where the light receiver was far away from the central body while the light emitter was in the strong-field region of the lens, corresponding to the famous gravitational redshift phenomenon. Under this scenario, the gravitational shift approximately depends only on the emitter's position. However, further consideration is necessary if both the light emission and reception events happen in the strong-field region for probing strong gravity in future. This is because the gravitational shift for this case depends not only on the position of the receiver but also on that of the emitter, which, however, may not be always known simultaneously. This paper gives a new numerical approach to tackle the issue how to obtain the emitter's position (and thus the gravitational shift) for a given observer and a given received direction of the light curve.
In this work, based on the exact equations of motion of light in harmonic coordinates, we study the gravitational shift of light induced by a Reissner-Nordström black hole numerically in the traditional strong-field scenario, as well as in the scenario where both the light emitter and the observer are located in the strong field of the lens. Since the emitter is connected to the receiver by a light-like world line, our focus is the determination of the unknown spatial position of one of the emitter and the receiver by integrating the null geodesics in the strong field limit, when the other's position is given. For this purpose, a new numerical technique is proposed. Two things should be pointed out. First, both the issue and the method of this work are different from that in the numerical investigation [17] and theoretical considerations [19,20]. Second, we note that rapid progress in the precision to measure the gravitational shift of frequency has been made during the last two decades. For example, the gravitational frequency shift of radio photons was measured with a high precision (∼ 10 −14 ) in the Cassini experiment in 2002 [38][39][40]. There might be a possibility to detect or constrain the intrinsic electric charge of a RN black hole by measuring its component contribution to the gravitational redshift, although it may be very weak up to now. This paper is organized as follows. In Section 2, the analytical form of the null gravitational redshift in Reissner-Nordström spacetime in the strong field limit is presented in harmonic coordinates. In Section 3, we first derive the exact null equations of motion in this spacetime, and then integrate them via a new numerical technique to calculate the spacial position of the light emitter for a given observer. The numerical results of the gravitational redshift are thus obtained for both the traditional and non-traditional strong-field scenarios. A summary is given in Section 4. In what follows, we use natural units in which G = c = 1.

Null Gravitational shift caused by a RN black hole in harmonic coordinates
Consider the gravitational frequency shift of light signals in the Reissner-Nordström geometry. The Reissner-Nordström metric in harmonic coordinates (X 0 , X 1 , X 2 , X 3 ) (defined as (T, X, Y, Z)) can be written as [41] where δ i j denotes Kronecker delta, i and j run over the values 1, 2, 3, X X X ·dX X X ≡ X 1 dX 1 +X 2 dX 2 +X 3 dX 3 , and R ≡ |X X X| = X 2 1 + X 2 2 + X 2 3 . M and Q are the rest mass and electrical charge of the gravitational source, respectively, with the relation M 2 ≥ Q 2 to avoid the naked singularity of the black hole. Figure 1 shows the schematic diagram for the propagation of light in the gravitational field of a RN black hole. The spacial coordinates of the light emitter (marked by A) and the receiver (marked by B) are denoted as A is connected to B by a null world line (red line). We denote the coordinate times of the emission and reception events as T A and T B respectively. For simplicity, we discuss the special case in which both the light emitter and the receiver are assumed to be static in the background's rest frame (T, X, Y, Z).
The gravitational redshift of frequency of light is conventionally defined by [17,[42][43][44] where τ A and τ B denote the proper times of the light emitter and the receiver, respectively. A negative redshift means a blueshift. For the case of the emitter and the receiver being static, dτ A dT A , dT B dτ B , and dT A dT B in the Reissner-Nordström spacetime read as follows: where Thus, the exact form of the gravitational redshift of frequency induced by a RN source in harmonic coordinates is easily rewritten as follows: In the conventional strong-field case, the light observer is far away from the central body (R B → +∞), and the red- However, this is not the case in our non-traditional scenario where both the light emitter and the observer are located in the strong-field region. Notice that the positions of the receiver Fig. 1 Schematic diagram for the light propagation in the gravitational field of a Reissner-Nordström black hole. The central body is assumed to be located at the origin O(0, 0, 0) of a three-dimensional Cartesian coordinate system (X, Y, Z). and the emitter may not be always known at the same time in some cases, which is different from the assumption that they are always known simultaneously given in previous work (see [45][46][47], and references therein). Since they are connected by a light-like geodesic, our focus is the determination of the position of one of the light emitter and the observer via integrating the exact null geodesics, once the spacial position of the other is given. Without loss of generality, the observer's position is assumed to be known while the emitter's position unknown in the following calculations.

Numerical calculations in the strong field limit
In this section, we first derive the exact null equations of motion in the RN spacetime. A new numerical technique, which is different from the Kerr Black Hole Ray Tracer technique [17], is then adopted to compute the gravitational redshift z RN in the strong field of the RN black hole for both the traditional and non-traditional strong-field scenarios.

Exact null geodesics in the RN geometry
Based on the harmonic metric given in Equation (1), we obtain the exact form of the equations of motion of light in the RN geometry via calculating the nonvanishing Christoffel symbols tediously but straightforwardly as follows: where a dot denotes the derivative with respect to the affine parameter ξ along the geodesic [48,49]. It can be seen that, dropping the electrical charge (Q = 0) in Equations (7) -(10), we can get the exact null geodesics of the Schwarzschild spacetime in harmonic coordinates which can be reduced to the result for the equatorial propagation [50] in the second post-Minkowskian approximation.

Basics of numerical simulations
For the convenience of discussion, three assumptions are made. First, we follow Wucknitz and Sperhake's idea [49] and assume the trajectory parameter ξ to take the dimension of length. Second, since the receiver B (X B , Y B , Z B ) for the traditional scenario is far enough away from the gravitational source, we assume |X B | to be large enough, with X B Y B < 0 and X B Z B < 0. We finally assume the light propagating vector k k k at the position of the receiver B to be approximately parallel to the X-axis of the coordinate system, which is not difficult to be realized by adjusting the coordinate frame. Notice that the light emitter A (X A , Y A , Z A ) is close to the black hole but beyond its gravitational radius, and that the receiver for the non-traditional scenario is denoted byB (XB, YB, ZB) and is located at some place (between A and B) of the same null world line.
With the consideration of the reversibility of light paths, there are three steps to determine the spatial coordinates of the light emitter A for calculating the redshift z RN : First, for convenience, we make our numerical integrating direction be reverse to the direction of the light propagation shown in Figure 1, namely, computing from the ob- Second, for a given observer B and a given X A , we need to integrate Equations (7) -(10) numerically to obtain Y A and Z A . Notice that for a given X B in the starting conditions, the adjustment of the coordinates Y B and Z B will make us get any desired position for the emitter A. Since there is only one specific value of ξ corresponding to one point of the photon's world line, we can obtain the numerical value of the affine parameter ξ A for the emitter via performing a definite integral forẊ over ξ as follows: with the affine parameter ξ B of the observer being given to be X B in the asymptotically flat spacetime [49,51]. It is worth to point out that there might be several numerical values of ξ A to correspond to a given X A in the strong field limit, and we take the one whose value is closet to X A in the following calculations for physical reasons. After the determination of ξ A , we can immediately get the numerical values of Y A and Z A from the following definite integrals: Finally, knowing the desired spatial position of the light emitter A(X A , Y A , Z A ), the gravitational redshift of light propagating from the emitter A to the receiver B in the field of a given RN black hole can be calculated according to Equation (6). We emphasize that the position of the observer B is the place where the starting conditions are given for both the traditional and non-traditional strong-field cases, and that the numerical process above is also applied to calculating the gravitational redshift of light propagating from the emitter A to the receiverB caused by the same central body.
In the numerical calculations, the computation domain is , and we use X max to replace the infinity +∞ for a large enough X max . The initial and boundary conditions for both the traditional and non-traditional strong-field cases are given as follows:

The traditional case
We first consider the redshift z RN as the function of the positions of the emitter and the receiver. Since the event horizon of a RN black hole is located at R H = M 2 − Q 2 in harmonic coordinates, the numerical values of Y B and Z B can not be too small for avoiding the swallow of light by the black hole. Table 1 gives the numerical values of the redshift z RN for variable Y B , Z B , and X A in units of M. We can see that the magnitues of z RN are much larger than nowadays precesion (∼ 10 −14 ) of the astronomical measurements of the shift. For the convenience of display, Figure 2 shows the propagating trajectories of light in the strong field of a RN black hole, plotted for different Y B and Z B . We next discuss the contribution of the electrical charge Q to the frequency shift. Table 2 gives the numerical values of the redshift z RN for variable Q and X A in units of M. It indicates that it is possible to detect the charge-produced contribution to the gravitational redshift in the high-accuracy observations, although the electrical charge of the black hole may be weak (such as Q = 0.01M).
With respect to the gravitational redshift as a function of the rest mass M of the gravitational lens for a given observer, Table 3 presents the numerical values of z RN for different rest mass M and the emitter's position (X A ) in units of M .

The non-traditional case
The numerical results above are based on the scenario where the light receiver B is located far enough from the black hole. We now consider the gravitational redshift of light coming from the emitter A (X A , Y A , Z A ) to the receiverB (XB, YB, ZB) (X B XB < X A ). With the starting conditions (i.e., Equations (18) -(25)), the spacial positions of both the light emitter A and the receiverB can be determined via the process in Subsection 3.3.1, if X A and XB are given. Note that in the scenario where both the light emitter and the receiver are located in the strong-field region, the spacial point (X B , Y B , Z B ) is just regarded as the position of the starting conditions rather than a reception point. In units of M, Table 4 shows the numerical values of the gravitational redshift z RN of light propagating from A toB for different values of Y B , Z B , and X A , with X B = −1.0 × 10 10 and XB = 20 as an example.

Summary
In this work, the gravitational redshift effect of light in the Reissner-Nordström geometry has been studied numerically, in both the traditional strong-field scenario and the non-traditional scenario where both the light emission and reception events happen in the strong-field region. In order to obtain the numerical results of the gravitational redshift, the exact null geodesics in the Reissner-Nordström spacetime have been derived and solved via a new numerical technique to determine the spatial position of the emitter for a given light observer. Since the positions of the emission and reception events are not always known simultaneously, this work might be helpful to the related astronomical measurements for probing strong gravity in future. Table 1 The numerical values of the redshift z RN for different Y B , Z B , and X A (in units of M). As an example, we set X B = −1.0 × 10 10 and Q = 0.05.   Table 3 The  Table 4 The numerical values of z RN for light propagating from the emitter A to the receiverB for different Y B , Z B , and X A (in units of M).