-
We study the motion of a charged test particle in the most general static and spherically symmetric spacetime, whose line element is described by
$ {\rm{d}}s^2=-A(r){\rm{d}}t^2+B(r){\rm{d}}r^2+C(r)({\rm{d}}\theta^2+\sin^2\theta{\rm{d}}\phi^2), $
(1) where
$ (t,r,\theta,\phi) $ are the coordinates, and the metric functions depend on r only. For the electromagnetic field in this spacetime, we assume that only a pure electric field described by the potential$ \left[{\cal{A}}_0\left(r\right),0,0,0\right] $ exists. However, in some spacetime solutions with electric fields, no such potential but the electric field itself$ E(r) $ is given. Fortunately, for the two methods presented in this paper, we only require the series expansion of$ {\cal{A}}_0(r) $ near the orbit radius; therefore, the electric field$ E(r) $ is also sufficient. Owing to the symmetry of the spacetime and field, without losing any generality, we can always investigate the equatorial motion of charged test particles in such spacetimes.The dynamics of a point particle possessing non-zero mass m and specific charge
$ \hat{q}\equiv q/m $ can be described by the Lorentz equation [29]:$ \frac{{\rm{d}}^2x^\rho}{{\rm{d}}\tau^2}+\Gamma^{\rho}{}_{\mu\nu}\frac{{\rm{d}}x^\mu}{{\rm{d}}\tau}\frac{{\rm{d}}x^\nu}{{\rm{d}}\tau}=\hat{q}{\cal{F}} ^{\rho}{}_{\mu}\frac{{\rm{d}}x^\mu}{{\rm{d}}\tau}, $
(2) where
$ F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu $ , and τ is the proper time of the test particle. After integrating Eq. (2) on the equatorial plane$ (\theta= {\pi}/{2}) $ , we obtain three first-order differential equations [13]:$ \dot{t}=\frac{\hat{E}+\hat{q}{\cal{A}}_0}{A}, $
(3a) $ \dot{\phi}=\dfrac{\hat{L}}{C}, $
(3b) $ \dot{r}^2=\left[\frac{(\hat{E}+\hat{q}{\cal{A}}_0)^2}{AB}-\frac{1}{B}\right]-\frac{\hat{L}^2}{BC}, $
(3c) where integration constants
$ \hat{E} $ and$ \hat{L} $ correspond to the energy and angular momentum per unit mass of the particle respectively, and$ \dot{{}} $ denotes the derivative with respect to proper time τ.To analyze the radial motion of the particle in more detail, we first rewrite Eq. (3c) in an alternative form:
$ B(r)C(r)\dot{r}^2=\left[\frac{(\hat{E}+\hat{q}{\cal{A}}_0(r))^2}{A(r)}-1\right]C(r)-\hat{L}^2. $
(4) The zeros and singularities of
$ B(r) $ and$ C(r) $ generally correspond to event horizons or singularities of the spacetime, and in this paper, we focus on region$ B(r)C(r)>0 $ . Thus, we can consider of the first term on the right-hand side of Eq. (4) as a type of effective potential, against which$ \hat{L}^2 $ can be compared with$ V(r)=\left[\frac{(\hat{E}+\hat{q}{\cal{A}}_0(r))^2}{A(r)}-1\right]C(r). $
(5) Subsequently, the motion of the particle is only allowed when
$ V(r)\geq \hat{L}^2 $ . We assume that this inequality has nonempty interval solution$ [r_1,\; r_2] $ where$ V(r_1)=V(r_2)=\hat{L}^2. $
(6) Clearly,
$ r_1 $ and$ r_2 $ are the radii of the periapsis and apoapsis. For simplicity, we also assume that this interval contains only one local maximum of$ V(r) $ at$ r=r_c $ . For later usage, we denote value$ V(r_c) $ as$ \hat{L}_c^2 $ , i.e.,$ \frac{{\rm{d}}V(r)}{{\rm{d}}r}\bigg|_{r=r_c}=0,\; V(r_c)\equiv \hat{L}_c^2. $
(7) The existence of such
$ r_c $ is always possible as long as$ \hat{L}^2 $ is sufficiently close to$ \hat{L}_c^2 $ . The physical interpretation of$ r_c $ is straightforward: it represents the radius of the circular orbit of the particle with energy$ \hat{E} $ and angular momentum$ \hat{L}_c $ . Note that, except for very simple spacetime and electric potential (e.g., the RN spacetime), solving the algebraic equation (7) to obtain the closed form solution of$ r_c $ is often difficult. Therefore, in the following sections where$ r_c $ is used, we assume that$ r_c $ is already solved, numerically if necessary.With the above consideration, the PS can be defined as
$ \alpha= 2\int_{r_1}^{r_2}\left|\frac{{\rm{d}}\phi}{{\rm{d}}r}\right|{\rm{d}}r -2\pi, $
(8) which, after using Eq. (3), becomes
$ \alpha=2\left(\int_{r_1}^{r_c}+\int_{r_c}^{r_2}\right)\sqrt{\frac{B(r)}{C(r)}} \frac{\hat{L} \; {\rm{d}} r}{\sqrt{V(r)-\hat{L}^2}}-2\pi. $
(9) To reveal the effect of electric interaction on the PS, our task becomes to systematically solve the PS in Eq. (9) using appropriate approximations. We present a perturbative approach based on NC approximation, and then in the next section, a PN method is used.
For the NC approximation, the first step for the computation of the integrals in Eq. (9) is to utilize the deviation of the particle's orbit from a perfect circle as a measure of perturbation for series expansion. This deviation can be quantified by parameter a defined as
$ a\equiv\frac{\hat{L}_c}{\hat{L}}-1. $
(10) When
$ r_1 $ and$ r_2 $ are close, i.e., the orbit is close to a circle, the above defined a will be a small quantity.To use a, we extend the method developed in Refs. [10, 13, 30] to the case with electric interaction. We first propose a change of variables in Eq. (9) from r to ξ, which are linked by
$ \xi\equiv \hat{L}_c P\left(\frac{1}{r}\right), $
(11) where function
$ P(1/r) $ is defined as$ P\left(\frac{1}{r}\right)\equiv \frac{1}{\sqrt{V(r)}}-\frac{1}{\hat{L}_c}. $
(12) From Eq. (7), we observe that function
$ P(x) $ exhibits opposite monotonic behavior on the two sides of$ x=1/r_c $ . This implies that its inverse function has two different branches when$ x\geqslant 0 $ . We use$ \omega_+ $ and$ \omega_- $ to denote the branch of the inverse function of$ P(x) $ for$ r>r_c $ and$ r< r_c $ respectively. In other words, using Eq. (11) we obtain$ r= \left\{\begin{array}{l}1/\omega_-\left( \xi/\hat{L}_c\right)\ \text{for}\; r<r_c,\\ 1/\omega_+\left( \xi/\hat{L}_c\right)\ \text{for}\; r\geqslant r_c. \end{array}\right. $
(13) Substituting this change of variables, the terms in the integral (9) become
$ \frac{V(r)}{\hat{L}^2}-1 \to \frac{(2+a+\xi)(a-\xi)}{(\xi+1)^2}, $
(14a) $ r_{1,2} \to \hat{L}_cP\left(\frac{1}{r_{1,2}}\right)=a, $
(14b) $ r_c \to \hat{L}_cP\left(\frac{1}{r_{c}}\right)=0, $
(14c) and consequently the PS in Eq. (9) can be rewritten as
$ \alpha=\int_0^a \frac{y(a,\,\xi)}{\sqrt{a-\xi}} \; {\rm{d}} \xi -2\pi, $
(15) where factor
$ y(a,\,\xi) $ of the integrand is$ y(a,\,\xi)= 2\sum\limits_{s=+,-}s\sqrt{\frac{B(1 / \omega_s)}{C(1 / \omega_s)}} \frac{\xi+1}{\hat{L}_c \sqrt{2+a+\xi}} \frac{-\omega_s^{\prime}}{\omega_s^2}. $
(16) An evident advantage of this change of variables is that it transforms two different integral limits
$ r_1 $ and$ r_2 $ into the same value a and therefore simplifies the computation.Because we are focusing on the
$ a\ll 1 $ case, the next step is naturally expanding the function$ y(\xi) $ for small ξ. Carrying out this expansion, we can show that it generally assumes the form$ y(a,\,\xi)=\sum\limits_{n=0}^{\infty}y_n(a) \xi^{n-\frac12}. $
(17) Here,
$ y_n(a) $ can be determined once the metric functions and potential$ {\cal{A}}_0 $ are known around$ r_c $ . Substituting this back into Eq. (15), and using the integral formula$ \int_0^a \frac{\xi^{n-\frac12}}{\sqrt{a-\xi}} \; {\rm{d}}\xi=\frac{(2n-1)!!}{(2n)!!}\pi a^n, $
(18) the PS is computed as
$ \alpha= \sum\limits_{n=0}^{\infty}\frac{(2n-1)!!}{(2n)!!}\pi y_n(a) a^n-2\pi. $
(19) Because we have assumed that deviation parameter a is small, in principle, we can further expand
$ y_n(a) $ for small a and reorganize the PS into a true power series:$ \alpha_{{\rm{NC}}}= \sum\limits_{n=0}^{\infty}t_n a^n =\sum\limits_{n=0}^{\infty}t_n \left(\frac{\hat{L}_c}{\hat{L}}-1\right)^n . $
(20) These coefficients
$ t_n $ are now independent of a and, similar to$ y_n $ , can be completely determined when the metric functions and electric potential are fixed. In the following, we show this determination process.We assume that the metric functions and the electric potential
$ {\cal{A}}_0 $ have the following expansions around$ r_c $ $ A(r)=\sum\limits_{n=0}^{\infty}a_{cn}(r-r_c)^n, $
(21a) $ B(r)=\sum\limits_{n=0}^{\infty}b_{cn}(r-r_c)^n, $
(21b) $ C(r)=\sum\limits_{n=0}^{\infty}c_{cn}(r-r_c)^n, $
(21c) $ {\cal{A}}_0(r)=\sum\limits_{n=0}^{\infty}h_{cn}(r-r_c)^n. $
(21d) The definition of
$ r_c $ in Eq. (7) immediately implies a constraint between these coefficients:$ \begin{aligned}[b] \hat{E}+\hat{q}h_{c0}=\;&\frac{a_{c0}c_{c0}}{a_{c0}c_{c1}-a_{c1}c_{c0}} \\ &\times \left(\sqrt{\frac{c_{c1}}{c_{c0}^2}\left(a_{c0}c_{c1}-a_{c1}c_{c0}\right)+\hat{q}^2h_{c1}^2}-\hat{q}h_{c1}\right). \end{aligned} $
(22) We will use this relation occasionally to simplify
$ t_n $ and$ T_n $ below. Critical$ \hat{L}_c $ in Eq. (20) can also be expressed in terms of these coefficients using again Eq. (7):$ \hat{L}_c=\sqrt{c_{c0} \left[\frac{\left(\hat{E}+\hat{q} h_{c0}\right){}^2}{a_{c0}}-1\right]}. $
(23) Substituting the expansions in Eq. (21) into Eq. (12), taking its inverse function
$ \omega_\pm $ and then substituting it into (16), expanding for small ξ, and eventually expanding$ y_n(a) $ again for small a, we an determine coefficients$ t_n $ in Eq. (20). Here, for simplicity, we list only the first two:$\begin{aligned}\\[-5pt] t_0= \frac{ 2\pi b_{c0}^{1/2} \hat{L}_c }{c_{c0}^{1/2}\left(-T_2\right)^{1/2}},\end{aligned} $ (24a) $ \begin{aligned}[b] t_1=\;& -\frac{\pi }{4 b_{c0}^{3/2} c_{c0}^{5/2} (-T_2)^{7/2}}\left\{\left[\left(-3 b_{c0}^2 c_{c1}^2+4 b_{c0}^2 c_{c0} c_{c2}-4 b_{c2} b_{c0} c_{c0}^2+2 b_{c1} b_{c0} c_{c0} c_{c1}+b_1^2 c_{c0}^2\right) T_2^2\right.\right. \\ &\left.\left.+\left(6 b_{c0} b_{c1} c_{c0}^2-6 b_{c0}^2 c_{c0} c_{c1}\right) T_2 T_3 -15 b_{c0}^2 c_{c0}^2 T_3^2+12 b_{c0}^2 c_{c0}^2 T_2 T_4\right]\hat{L}_c^3 -8 b_{c0}^2 c_{c0}^2 T_2^3 \hat{L}_c \right\}, \end{aligned} $
(24b) where
$ T_n $ represents the n-th derivative of the potential at$ r_c $ , i.e.,$ T_n\equiv V^{(n)}(r_c)/n! $ , and the first two of them are$ \begin{aligned}[b] T_2=\;&\frac{\left(\hat{E}+\hat{q}h_{c0}\right)^2}{a_{c0}^3}\left(a_{c0}^2 c_{c2}-a_{c2} a_{c0} c_{c0}-a_{c1} a_{c0} c_{c1}+a_{c1}^2 c_{c0}\right) \\ &+\frac{2 \hat{q} \left(\hat{E}+\hat{q}h_{c0}\right)}{a_{c0}^2}\left(-a_{c1} c_{c0} h_{c1}+a_{c0} c_{c1} h_{c1}+a_{c0} c_{c0} h_{c2}\right)+\frac{ \hat{q}^2 c_{c0} h_{c1}^2}{a_{c0}}-c_{c2} , \end{aligned}$
(25a) $\begin{aligned}[b] T_3=\frac{\left(\hat{E}+\hat{q}h_{c0}\right)^2}{a_{c0}^4}\left[a_{c0} a_{c1}^2 c_{c1}+a_{c0} a_{c1} \left(2 a_{c2} c_{c0}-a_{c0} c_{c2}\right)-a_{c1}^3 c_{c0}+a_{c0}^2 \left(-a_{c3} c_{c0}-a_{c2} c_{c1}+a_{c0} c_3\right)\right] \end{aligned} $
$\begin{aligned}[b] &+\frac{2 \hat{q} \left(\hat{E}+\hat{q}h_{c0}\right)}{a_{c0}^3}\left[a_{c0} a_{c0} \left(c_{c2} h_{c1}+c_{c1} h_{c2}+c_{c0} h_{c3}\right)-a_{c0} a_{c2} c_{c0} h_{c1}-a_{c0} a_{c1} \left(c_{c1} h_{c1}+c_{c0} h_{c2}\right)\right.\left.+a_{c1}^2 c_{c0} h_{c1}\right]\\ & +\frac{\hat{q}^2 h_{c1} }{a_{c0}^2}\left(-a_{c1} c_{c0} h_{c1}+a_{c0} c_{c1} h_{c1}+2 a_{c0} c_{c0} h_{c2}\right)-c_{c3} .\end{aligned} $
(25b) The higher order
$ t_n\; (n=2,\,3,\cdots) $ and$ T_n\; (n=4,\,5,\cdots) $ can also be obtained easily using an algebraic system.In obtaining the PS in Eq. (20), although we used small a or equivalently small eccentricity e approximation, we need not assume that the orbit size itself is large compared with the characteristic length scale (e.g. mass M) of the spacetime. The latter is the key assumption for the PN method, which can also be used to compute the corresponding PS. In next section, we show how the PN method can be used to consider the electric interaction and to determine the PS when the orbit size is large.
-
Many previous studies computed the PS using the PN method. However, most of them, if not all, focused only on the lowest order(s) result, and some used this method very loosely and the results are large. Moreover, only a few of them used the PN method to solve the gravitational and electric interactions simultaneously. In this section, we systematically develop a PN method that not only yields the PS to a high order of the semilatus rectum p but also considers the electric interaction. Some of the techniques we use here are similar to those in Refs. [10, 15, 31]. In the last section, we further show that the PN method result agrees perfectly with the result obtained using the NC method when the orbit size is large and the eccentricity is small.
For the PN method, we first assume that the orbit can be described by the relation
$ \dfrac{1}{r}=\dfrac{1}{p}\left[1+e\cos \psi\left(\phi\right)\right] , $
(26) where function
$ \psi(\phi) $ describes the deviation of the orbit from ellipse. Note that the periapsis and apoapsis of the orbit correspond to$ \psi=0 $ and$ \psi=\pi $ respectively, and radii of the periapsis$ r_1 $ and apoapsis$ r_2 $ satisfy$ \frac{1}{r_1}=\frac{1+e}{p},\quad \frac{1}{r_2}=\frac{1-e}{p}. $
(27) The radius of the orbit evolves for a full period when ψ changes from 0 to
$ 2\pi $ ; therefore, the PS using this description becomes$ \alpha_{{\rm{PN}}}= \int_0^{2\pi}\frac{{\rm{d}}\phi}{{\rm{d}}\psi}{\rm{d}}\psi-2\pi . $
(28) $ {\rm{d}}\phi/{\rm{d}}\psi $ can be expressed from Eq. (26) as$ \frac{{\rm{d}}\phi}{{\rm{d}}\psi}=-\frac{e}{p} \frac{\sin \psi}{{\rm{d}}u / {\rm{d}}\phi} , $
(29) whence we have defined
$ u=1/r $ . Term$ {\rm{d}}u/{\rm{d}}\phi $ in Eq. (29) can be transformed using Eq. (3) to$ \begin{aligned}[b] \left(\frac{{\rm{d}}u}{{\rm{d}}\phi}\right)^2 =\;&u^4\frac{\dot{r}^2}{\dot{\phi}^2} =u^4\left[\frac{(\hat{E}+\hat{q}{\cal{A}}_0(\frac{1}{u}))^2C(\frac{1}{u})^2}{A(\frac{1}{u})B(\frac{1}{u})\hat{L}^2}-\frac{C(\frac{1}{u})^2}{B(\frac{1}{u})\hat{L}^2}\right]\\&-\frac{u^4C(\frac{1}{u})}{B(\frac{1}{u})} \equiv F(u) \end{aligned} $
(30) where we have defined the right-hand side as
$ F(u) $ . Next, we show that$ F(u) $ can be expressed as a serial function of ψ, and after substituting it back into (28), the integral over ψ can be performed to determine the PS.The key to accomplishing this is the PN assumption, i.e., we assume that p is large and therefore u is always a small quantity. This enables us to expand
$ F(u) $ as a power series of u:$ F(u)=\sum\limits_{n=0}^\infty f_nu^n, $
(31) where coefficients
$ f_n $ can be determined from Eq. (30) using the asymptotic expansion of the metric functions and$ {\cal{A}}_0 $ . Assuming that they are of the form$ \begin{aligned}[b] A(r) =1+\sum\limits_{n=1} \frac{a_n}{r^n},&\ B(r) =1+\sum\limits_{n=1} \frac{b_n}{r^n}, \\ \frac{C(r)}{r^2}=1+\sum\limits_{n=1} \frac{c_n}{r^n},&\ {\cal{A}}_0 =\sum\limits_{n=1} \frac{\mathfrak{q}_n}{r^n} , \end{aligned} $
(32) we can then determine
$ f_n $ order by order as functions of coefficients$ a_n,\; b_n,\; c_n $ , and$ \mathfrak{q}_n $ . Here, we only illustrate the first few of them:$ f_0= \frac{ \hat{E}^2-1}{ \hat{L}^2}, $
(33a) $ f_1= \frac{ \left(b_1 -2 c_1\right) \left( 1-\hat{E}^2\right)- a_1 \hat{E}^2-2 \hat{q} \mathfrak{q}_1 \hat{E}}{\hat{L}^2}, $
(33b) $ \begin{aligned}[b] f_2=\;& -1+\frac{1}{\hat{L}^2}\left\{\left(a_1 b_1-2 a_1 c_1+a_1^2-a_2\right)\hat{E}^2\right. \\ &\left.+\left(-2 b_1 c_1+b_1^2-b_2+c_1^2+2 c_2\right) \left(\hat{E}^2-1\right)\right. \\ &\left.+ 2 \hat{q} \hat{E}\left[ \left(-a_1-b_1+2 c_1\right)\mathfrak{q}_1 +\mathfrak{q}_2\right] + \hat{q}^2 \mathfrak{q}_1^2\right\}. \end{aligned}$
(33c) Repetitively substituting Eq. (26) as u, all u instances on the right hand side of Eq. (31) can be completely replaced by variables
$ p,\, e $ , and ψ. After collecting terms proportional to$ \cos\psi $ and$ \sin^{2n}\psi $ , we show that$ F(u) $ always assumes the form [10]$ \begin{aligned}[b] F(u)=\;&G_0(\hat{E},\hat{L},p,e)+G_0^\prime(\hat{E},\hat{L},p,e)\cos \psi \\ &+\sum\limits_{n=1}^\infty G_n(\hat{E},\hat{L},p,e,\cos \psi)\sin^{2n} \psi ,\end{aligned} $
(34) where
$ G_0,\; G_0^\prime,\; G_n $ are linear combinations of$ f_n $ with coefficients being power series functions of$ \hat{E},\; 1/\hat{L},\; e $ , and$ 1/p $ . Their exact forms can be dertemined without difficulty although the algebra is too tedious to show here.Because at the apoapsis and periapsis, by definition (30) we have
$ F(u)=0 $ , using the condition (27), we immediately obtain$ F(\psi=0)=F(\psi=\pi)=0 $ , i.e.,$ G_0(\hat{E},\hat{L},p,e)=0,~\; G_0^\prime(\hat{E},\hat{L},p,e)=0. $
(35) These two power series conditions effectively establish two relations between the kinetic variables
$ (\hat{E},\, \hat{L}) $ and$ (p,\, e) $ . These conditions can be solved perturbatively when p is large, enabling us to express$ (\hat{E},\, \hat{L}) $ in terms of$ (p,\, e) $ $ \begin{aligned}[b] \hat{E}=\;&1-\frac{\left(1-e ^2\right) \left(2 \hat{q} \mathfrak{q}_1 -a_1 \right)}{4 }\frac{1}{p} \\ &+\frac{\left(1-e ^2\right)^2 \left(2 \hat{q} \mathfrak{q}_1 -a_1 \right) \left(2 \hat{q} \mathfrak{q}_1 -3 a_1 +4 c_1 \right)}{32 }\frac{1}{p^2}+{\cal{O}}\left(p\right)^3, \end{aligned} $
(36a) $ \begin{aligned}[b] \hat{L}^2=\;&\frac{ m \left(2 \hat{q} \mathfrak{q}_1 -a_1 \right)}{2 }p \\ &+\frac{1}{4 }\left\{\left[-3 a_1 c_1+3 a_1^2-4 a_2+e^2 \left(a_1^2-a_1 c_1\right)\right]\right. \\ &\left.+\hat{q} \left[-5 a_1 \mathfrak{q}_1+6 c_1 \mathfrak{q}_1+8 \mathfrak{q}_2+e^2 \left(2 c_1 \mathfrak{q}_1-3 a_1 \mathfrak{q}_1\right)\right]\right. \\ &\left.+2\hat{q}^2 \mathfrak{q}_1^2\left( e^2+1\right)\right\}+{\cal{O}}\left(p\right)^{-1}. \end{aligned}$
(36b) In the final expression for the PS, these expressions can aid us in eliminating the dependence on unnecessary kinetic variables. Note that other methods can be used to obtain the relation between
$ (\hat{E},\, \hat{L}) $ and$ (p,\, e) $ . Substituting Eq. (27) into Eq. (6), we obtain$ V\left(\frac{p}{1+e}\right)=\hat{L}^2=V\left(\frac{p}{1-e}\right). $
(37) Clearly, this is more concise than Eq. (35), and we can show that it is essentially equivalent to Eq. (35).
Now substituting Eq. (35) into Eq. (34) and further into Eq. (29), we find
$ \begin{aligned}[b]\\[-8pt] \frac{{\rm{d}}\phi}{{\rm{d}}\psi}=\;&\frac{e}{p}\left[\sum\limits_{n=1}^\infty G_n(\hat{E},\hat{L},p,e,\cos \psi)\sin^{2n-2} \psi\right]^{-\frac{1}{2}} =\frac{1}{\sqrt{-f_2}}-\frac{f_3 (e \cos (\psi)+3)}{2 \sqrt{-f_2} }\frac{1}{p}+\frac{1}{16 \left(-f_2\right){}^{3/2} }\left[3f_3^2 \left(e ^2+18\right)-12 f_2 f_4 \left(e ^2+4\right)\right.\\ &\left.+ \left(36 f_3^2-32 f_2 f_4\right) e \cos (\psi)+ \left(3 f_3^2-4 f_2 f_4\right) e ^2 \cos (2 \psi)\right]\frac{1}{p^2}+{\cal{O}}[p^{-3}] ,\end{aligned} $ (38) where, in the second step, we have substituted the first few
$ G_n\; (n=1,2,\cdots) $ and performed the large p expansion again. Substituting into Eq. (28) and noting that terms proportional to$ \cos(k \psi)\; (k=1,2,\cdots) $ will not survive the integration, the result for the PN PS is simply determined as$ \begin{aligned}[b] \alpha_{{\rm{PN}}}=\;&\frac{2 \pi }{\sqrt{-f_2}}+\frac{3 \pi f_3}{\left(-f_2\right){}^{3/2} p}\\ &+\frac{3 \pi \left[f_3^2 \left(e ^2+18\right)-4 f_2 f_4 \left(e ^2+4\right)\right]}{8 \left(-f_2\right){}^{5/2} p^2}+{\cal{O}}\left(p\right)^{-3}. \end{aligned} $
(39) However, this result is still not a true series of p because of the dependence of
$ f_n $ on$ (\hat{E},\, \hat{L}) $ and equivalently on$ (p,\, e) $ . Substituting Eq. (36) as$ f_n $ in Eqs. (33) and re-expanding for large p, we can obtain a true power series form of p for$ \alpha_{{\rm{PN}}} $ $ \alpha_{{\rm{PN}}}=\sum\limits_{n=1}^{\infty} \left(\sum\limits_{j=0}^{n-1}d_{n,j}e^{2j}\right)\frac{1}{p^n}, $
(40) with the coefficient for order
$ p^{-n} $ , a polynomial of order$ 2(n-1) $ of eccentricity e. The first few coefficients for$ d_{n,j} $ are$ d_{1,0}= \frac{\pi}{ a_1 -2 \hat{q} \mathfrak{q}_1 }\left\{a_1 b_1+a_1 c_1-2 a_1^2+2 a_2+2\hat{q} \left[\left(2a_1-b_1-c_1\right)\mathfrak{q}_1-2\mathfrak{q}_2\right] -2\hat{q}^2\mathfrak{q}_1^2 \right\}, $
(41a) $ \begin{aligned}[b] d_{2,0}=\;&\frac{\pi}{\left( a_1 -2 \hat{q} \mathfrak{q}_1\right)^2} \left\{ \frac{1}{2} a_1^2 b_1 c_1- a_1^3 b_1-\frac{1}{4} a_1^2 b_1^2+ a_1^2 b_2+ a_2 a_1 b_1-4 a_1^3 c_1-2 a_1^2 c_1^2+4 a_1^2 c_2\right. +4 a_2 a_1 c_1+5 a_1^4-8 a_2 a_1^2+6 a_3 a_1\\ &- a_2^2+\hat{q} [ (-2 a_1 b_1 c_1+4 a_1^2 b_1+ a_1 b_1^2-4 a_1 b_2 -2 a_2 b_1+16 a_1^2 c_1+ a_1 c_1^2-16 a_1 c_2-8 a_2 c_1-20 a_1^3+31 a_2 a_1-12 a_3)\mathfrak{q}_1\\ &+(-2 a_1 b_1 -8 a_1 c_1+9 a_1^2+4 a_2)\mathfrak{q}_2-12 a_1\mathfrak{q}_3]+\hat{q}^2[(-5 a_1 b_1-20 a_1 c_1+25 a_1^2-20 a_2 +2 b_1 c_1 \\ &- b_1^2+4 b_2- c_1^2+16 c_2)\mathfrak{q}_1^2+(-30 a_1 +4 b_1 +16 c_1 )\mathfrak{q}_1\mathfrak{q}_2-4 \mathfrak{q}_2^2 ]+\hat{q}^3[(-11a_1 +2b_1+8c_1 )\mathfrak{q}_1^3+20\mathfrak{q}_1^2\mathfrak{q}_2] +\hat{q}^4\mathfrak{q}_1^4 \},\\[-7pt]\end{aligned} \tag{41b} $
(41b) $ \begin{aligned}[b] d_{2,1}=\;&\frac{\pi}{\left( a_1 -2 \hat{q} \mathfrak{q}_1\right)^2} \left\{ -\frac{1}{4} a_1^2 b_1 c_1+\frac{1}{2} a_1^2 b_2-\frac{1}{8}a_1^2 b_1^2+ a_1^3 c_1-\frac{5}{8} a_1^2 c_1^2+\frac{1}{2} a_1^2 c_2- a_2 a_1 c_1 \right. \\ &\left.+\hat{q} \left[\left( a_1 b_1 c_1+\frac{1}{2} a_1 b_1^2-2 a_1 b_2-4 a_1^2 c_1+\frac{5}{2} a_1 c_1^2-2 a_1 c_2+2 a_2 c_1+ a_2 a_1\right)\mathfrak{q}_1\right.\right. \\&\left.\left.+\left(-a_1^2+2a_1c_1\right)\mathfrak{q}_2\right]+ \hat{q}^2\left[\left(5 a_1 c_1- a_1^2-2 a_2- b_1 c_1-\frac{1}{2} b_1^2+2 b_2-20 c_1^2+2 c_2\right)\mathfrak{q}_1^2\right.\right. \\ &\left.\left.+\left(2a_1-4c_1\right)\mathfrak{q}_1\mathfrak{q}_2\right]+\hat{q}^3\left[\left(3a_1-2c_1\right)\mathfrak{q}_1^3\right]-2\hat{q}^4\mathfrak{q}_1^4\right\}.\end{aligned} $
(41c) Before we analyze the implication of result (40) in the next section, we compare the two PSs,
$ \alpha_{{\rm{NC}}} $ in Eq. (20) obtained using NC approximation, and$ \alpha_{{\rm{PN}}} $ . Clearly, these two PSs are only comparable when the orbit is both large and near-circular.Under this assumption, we now rewrite the PS
$ \alpha_{{\rm{NC}}} $ using PN kinetic variables$ (p,\, e) $ . Hence, we must express$ (\hat{E},\; \hat{L},\; r_c,\; a) $ in terms of$ (p,\, e) $ and other spacetime parameters. First, note that the first two of them,$ (\hat{E},\, \hat{L}) $ , have already been linked to$ (p,\, e) $ using Eq. (36). For$ r_c $ , although the defining Eq. (7) is not always solvable in closed form, it can always be solved perturbatively for large p. Substituting this into Eqs. (36) and asymptotic (32), the solution to$ r_c $ as functions of$ (p,\, e) $ can be obtained using the method of undetermined coefficients:$ \begin{aligned}[b] r_c=\;&\frac{p}{1-e^2}+ \left\{ -a_1^2 c_1+a_1 c_2+a_2 c_1+a_1^3-2 a_2 a_1+a_3\right.\\ &\left.+\hat{q}\left[\left(2 a_1 c_1 -2 a_1^2 +2 a_2 -2 c_2\right) \mathfrak{q}_1+\left(2 a_1-2 c_1\right) \mathfrak{q}_2-2 \mathfrak{q}_3\right]\right.\\ &\left.+\hat{q}^2\left[\left(a_1 -c_1\right) \mathfrak{q}_1^2-2 \mathfrak{q}_1 \mathfrak{q}_2\right]\right\} \frac{e^2}{(a_1-2 \mathfrak{q}_1)p}+{\cal{O}}(p)^{-2}. \end{aligned} $
(42) Similarly, by substituting Eqs. (36), (32), and (42) into
$ \hat{L}_c $ in Eq. (7) and further into Eq. (10), small parameter a can also be expressed using$ (p,\, e) $ as$\begin{aligned}[b] a=\;&\frac{1}{\sqrt{1-e^2}}-1+\frac{e^2}{\sqrt{1-e^2}\left(a_1-2 \hat{q} \mathfrak{q}_1\right)}\left\{-a_1 c_1+a_1^2-a_2\right.\\ &\left.+\hat{q}\left[\left(-2 a_1 +2 c_1\right) \mathfrak{q}_1+2 \mathfrak{q}_2\right]+\hat{q}^2\mathfrak{q}_1^2\right\}\frac{1}{p}+{\cal{O}}\left(\frac{1}{p}\right)^2. \end{aligned} $
(43) Eqs. (42) and (43) enable us to express all terms that depend on
$ r_c $ and a explicitly in$ \alpha_{{\rm{NC}}} $ in Eq. (20) in terms of the asymptotic expansion coefficients$ a_n,\,b_n,\,c_n,\,\mathfrak{q}_n $ and kinetic variables$ (p,\,e) $ . The expansion coefficients of the metric and electric potential functions around$ r_c $ , i.e.,$ a_{cn},\,b_{cn},\,c_{cn},\,h_{cn} $ , can also be expressed using$ a_n,\,b_n,\,c_n,\,\mathfrak{q}_n $ and variables$ (p,\,e) $ . Hence, we need only to substitute the asymptotic forms (32) and Eq. (42) into the expansion (21) and then further expand these functions for large p. Although$ a_{cn},\,b_{cn},\,c_{cn},\,h_{cn} $ can be solved easily, their expressions are too lengthy to present here. Finally, substituting all these coefficients, as well as Eqs. (42) and (43), into Eq. (20), we can express$ \alpha_{{\rm{NC}}} $ as power series of$ 1/p $ . Not surprisingly, we find that the result agrees perfectly with the PN formula (40). This agreement demonstrates the correctness of results from both methods. Moreover, for the coefficient of each fixed order of$ p^{-n} $ , its dependence on the eccentricity e is automatically in a polynomial form, indicating that if the small a expansion in Eq. (20) is performed to high enough order, the rewritten PN result to low orders will be valid even for large e. -
In this section, we implement the two aforementioned methods to determine PSs within several specific spacetimes and investigate the effect of the electric interaction and other spacetime parameters on these shifts.
-
The RN spacetime is the simplest and cleanest for the analysis of the PSs of charged particles. However, as we show next, it still captures the main characteristic of the electric interaction effect on the PS. The line element and electric potential of the RN spacetime are given by
$ A(r)=\frac{1}{B(r)}=1-\frac{2M}{r}+\frac{Q^2}{r^2},\; C(r)=r^2, $
(44a) $ {\cal{A}}_0(r)=-\frac{Q}{r}. $
(44b) Next, we analyze
$ \alpha_{{\rm{NC}}} $ in Eq. (20) and$ \alpha_{{\rm{PN}}} $ in Eq. (40) in this spacetime. First, substituting Eq. (44) into Eq. (7), we see that$ r_c $ for the RN spacetime is determined by a quartic polynomial:$\begin{aligned}[b] r_c^4 \left(\hat{E}^2-1\right)+r_c^3 \left(-3 \hat{E}^2 M-\hat{q}\hat{E} Q+4 M\right) \end{aligned}$
$\begin{aligned}[b] &+2 r_c^2 \left( \hat{E}^2 Q^2+2\hat{q} \hat{E} M Q-2 M^2- Q^2\right) \\ &+r_c\left(-3 \hat{q}\hat{E} Q^3+4 M Q^2-\hat{q}^2 M Q^2\right)+Q^4 \left(\hat{q}^2-1\right)=0 . \end{aligned}$
(45) Although its solution can still be expressed in a closed form, it is too lengthy to show here [13]. When
$ r_c $ is determined, using$ \hat{L}_c $ in Eqs. (7) and (10), we immediately determine small parameter a:$ a=\frac{r_c}{\hat{L}}\left[\frac{( r_c \hat{E} -\hat{q} Q)^2}{r_c (r_c-2 M)+Q^2}-1\right]^{{1}/{2}}-1. $
(46) Using Eq. (24), coefficients
$ t_n $ for RN spacetime can be determined without much difficulty. After combining with a, we finally obtain$ \alpha_{{\rm{NC}}} $ in RN spacetime:$ \begin{aligned}[b] \alpha_{{\rm{RN,NC}}}=\;&2\pi r_c\left[\frac{( r_c \hat{E} -\hat{q} Q)^2}{r_c (r_c-2 M)+Q^2}-1\right]\left[r_c (r_c-2 M)+Q^2\right] / \left\{ \left(1-\hat{E}^2\right)r_c^6 -6 M \left(1- \hat{E}^2\right)r_c^5-3 \left(\hat{E}^2-1\right) \left(4 M^2+Q^2\right)r_c^4\right.\\ &\left.+ \left[-2 Q \left(-8 \hat{E}^2 M Q+ \hat{q}\hat{E} Q^2-4 \hat{q} \hat{E} M^2+\hat{q}^2 M Q\right)-8 M^3-6 M Q^2\right]r_c^3 +3 Q^2 \left[Q \left(-2 \hat{E}^2 Q-4 \hat{q}\hat{E} M +\hat{q}^2 Q\right)\right.\right.\\ &\left.\left.+ 4 M^2+Q^2\right] r_c^2 +6 Q^4 \left(\hat{q}\hat{E} Q- M\right)r_c+Q^6 \left(1-\hat{q}^2\right)\right\}^{\frac{1}{2}}-2\pi +{\cal{O}}\left(a\right)^1, \end{aligned} $ (47) Higher orders can also be determined but are too lengthy to show here. Using a computer algebraic system, we computed this result to the tenth order of a.
To determine the PN PS in RN spacetime, we need only expand the metric functions and electric potential (44) asymptotically and then substitute the coefficients into Eq. (40). These steps are very simple, and the result to order
$ p^{-2} $ is$ \begin{aligned}[b] \alpha_{{\rm{RN,PN}}}=\;&\pi \left( 6 -\hat{Q}^2-6 \hat{q} \hat{Q} +\hat{q}^2 \hat{Q}^2 \right)\frac{M}{(1-\hat{q}\hat{Q} )p}\\&+\frac{\pi }{4 }\left\{\left[6 \left(18+e ^2\right)-2 \hat{Q}^2 \left(24+e ^2\right)-\hat{Q}^4\right]\right. \\ &-2 \hat{q} \hat{Q} \left[6 \left(18+e ^2\right)- \hat{Q}^2 \left(37+3 e ^2\right)\right]\\&+2 \hat{q}^2 \hat{Q}^2 \left[ \left(66+e ^2\right)-2 \hat{Q}^2 \left(6+e ^2\right)\right] \\ &\left.-2 \hat{q}^3 \hat{Q}^3 \left(13-3 e ^2\right)+\hat{q}^4 \hat{Q}^4 \left(1-2 e ^2\right)\right\}\\&\times\left[\frac{M}{ (1-\hat{q}\hat{Q} )p}\right]^2+{\cal{O}}\left(\frac{M}{p}\right)^3, \end{aligned} $
(48) where
$ \hat{Q}\equiv Q/M $ denotes the charge-to-mass ratio of the central object. As indnicated in the last section, when p is large and e is small, this PS should coincide with the NC result$ \alpha_{\rm{RN,NC}} $ in Eq. (47). We have converted all$ (\hat{E},\, \hat{L}) $ instances in$ \alpha_{\rm{RN,NC}} $ , including those in$ r_c $ and a to$ (p,\, e) $ , and the result agrees perfectly with (48), which is naturally small e expanded. When truncated to first order and set to neutral test particle, Eq. (48) agrees with results in Refs. [14–17, 19–23]. For nonzero$ \hat q $ , Eq. (48) to the first order agrees with Eq. (28) of Ref. [24] and Eq. (42) of Ref. [28]. However, the results in Refs. [25, 26] cannot be determined from our work.To check the correctness of the PSs in Eqs. (47) and (48), Fig. 1 plots the difference between the series result
$ \alpha_{\rm{RN,NC}} $ (Fig. 1(a) and 1(b)) and$ \alpha_{\rm{RN,PN}} $ (Fig. 1(c) and 1(d)) truncated to order$ \bar{n} $ and PS$ \alpha_{{\rm{num}}} $ computed using the numerical integration of the definition (8). The latter is performed to very high accuracy and can therefore be considered the true value of the PS. In both plots, we fix$ |\hat{Q}\hat{q}|\ll 1 $ to ensure that the orbit is a bound one and the PS is well defined. For the NC PS plotted in Fig. 1(a) and 1(b), although$ \alpha_{\rm{RN,NC}} $ only explicitly depends on$ (\hat{E},\, \hat{L}) $ , we still convert them to$ (p,\, e) $ using Eq. (36) to make a comparison with Fig. 1(c) and 1(d) for PN PS$ \alpha_{\rm{RN,PN}} $ . Note that in Fig. 1(a) and 1(c), p is varied, whereas in Fig. 1(b) and 1(d), e is changed.Figure 1. (color online) Differences between the numerical and analytic results in Eqs. (47) and (48) truncated to order
$ \bar{n} $ , with parameters fixed at$ \hat{Q}=1/2,\hat{q}=-1/5,p=500M,e=1/10 $ , except for the running one.All the plots in Fig. 1 show that as the truncation order increases, both the NC and PN results approach the true value of the PS exponentially. This shows that both methods are effective as the series order increases. Moreover, for
$ \alpha_{\rm{RN,NC}} $ Fig.1(a) shows that when e is small, the NC approximation almost functions equally well for large p and small p, and Fig. 1(b) shows that this method is very sensitive to e. The smaller the e, i.e, the more circular the orbit is, the more accurate the$ \alpha_{\rm{RN,NC}} $ . In contrast, we observe that the characteristics of plots Fig.1(c) and 1(d) are the opposite of those of Fig.1(a) and 1(b): the accuracy of$ \alpha_{\rm{RN,PN}} $ is not sensitive to e but increases rapidly as p increases. Again, this corresponds with the expectation of a PN approximation of the orbit. The abovementioned characteristics were also observed in Ref. [10], where the effect of the spacetime spin was studied.Effect of
$ \hat{q} $ and$ \hat{Q} $ on PSWith the correctness of the PS, particularly the PN PS (48), verified, we can now use it to analyze the physics contained in this formula. A few points must be noted. First, this PS depends on parameters
$ q,\; m,\; Q,\; M $ only through the two charge-to-mass ratios$ \hat q $ and$ \hat Q $ , indicating a certain redundancy in these parameters. Therefore, in the following, when analyzing the effects of these parameters, we need only focus on these two ratios. Second, the charge q appears only in the form of$ (qQ)^n $ , showing that q influences the PS only through the Coulomb interaction. Q also manifests in terms that are not a product with q; therefore, it can affect the PS through the gravitational channel. Third, because this result is post-Newtonian, it is more accurate in the asymptotic regions, where gravity follows the universal gravitational law. In such regions, the gravitational attraction and the Coulomb force, which depends on the sign of$ qQ $ , will strengthen or weaken each other, resulting in a net interaction proportional to$ (mM-qQ) $ under the natural unit system we are using. Therefore, for the orbit to be bounded such that a study of the PS is meaningful, this total force must be attractive, which means that$ mM-qQ>0 $ , i.e,$ 1-\hat{q}\hat{Q}>0 $ . This point is reflected in the denominator of each order in Eq. (48).To fully study the effect of parameters
$ (m,\; q,\; Q,\; M,\; p) $ on the PS using Eq. (48), we must first determine the boundary of the parameter space in which this series result is convergent. Inspecting Eq. (48) more closely, we observe that this formula contains only the power series of the three quantities in the left-hand sides of the inequalities in Eq. (49), as well as their product series. For the total series to converge, the necessary and sufficient condition is that the sizes of all these three quantities are less than one, i.e.,$ \begin{aligned}[b] \left[\frac{1}{(1-\hat q\hat Q)}\frac{M}{p}\right]^{1/2}&\leq 1, \\ \left[\frac{\hat Q^2}{(1-\hat q\hat Q)}\frac{M}{p}\right]^{1/2}&\leq 1, \\ \left[\frac{\hat q^2\hat Q^2}{(1-\hat q\hat Q)}\frac{M}{p}\right]^{1/2}&\leq 1. \end{aligned} $
(49) These conditions can be further simplified to
$ 1-\hat{q}\hat{Q}\geq \frac{M}{p},\; 1-\hat{q}\hat{Q}\geq\hat{Q}^2\frac{M}{p},\; 1-\hat{q}\hat{Q}\geq\hat{q}^2\hat{Q}^2\frac{M}{p}. $
(50) We plot a parameter space spanned by
$ (\hat{q}, \hat{Q}) $ in Fig. 2 to show the allowed regions bounded by condition (50) and requirement$ 1-\hat{q}\hat{Q}>0 $ . This region is primarily concentrated around the$ \hat{q} $ and$ \hat{Q} $ axes. In the$ \hat{Q} $ direction, it is completely bounded for given$ M/p $ , whereas in the$ \hat{q} $ direction, the region is not limited. All three conditions in Eq. (50) are effective in some part of this parameter space. In our following analyses of the effects of various parameters, we limit the parameter ranges according to this figure. Note that because the PS depends on$ \hat q $ and$ \hat Q $ only through$ \hat q\hat Q $ and$ \hat Q^2 $ , in principle, we can further limit our analysis to the case of$ \hat Q\geq0 $ for general$ \hat q $ . The case with$ \hat Q<0 $ is deducible by switching the sign of$ \hat q $ but maintaining that of$ \hat Q $ . However, to be as straightforward as possible, we make the statements in the following applicable to any sign choices of$ \hat q $ and$ \hat Q $ .Figure 2. (color online) Allowed parameter space for PN PS
$ \alpha_{\rm{RN,PN}} $ to be valid (gray area). The red, blue, and brown boundaries are due to the first, second, and third conditions in Eq. (50), respectively. We fix$ p=10M $ in this plot to separate the curves visually. A more relativistic value of p will only change this figure quantitatively.In the following analysis, we examine the influence of variables
$ \hat q $ and$ \hat Q $ on the PS. We focus on the leading order(s) of$ M/p $ . When$ \hat q\hat Q\ll 1 $ , we can expand the denominator of Eq. (48), and the PS to order$ (\hat{q}\hat{Q})^2 $ becomes$\begin{aligned}[b]\alpha_{\rm{RN,PN}}=\;& \frac{\pi M}{p}\left[6-\hat{Q}^2-\hat{q}\hat{Q}^3 +\left(1-\hat{Q}^2\right) \hat{q}^2\hat{Q}^2 \right]\\ &+{\cal{O}}\left[\frac{M}{p},(\hat{q}\hat{Q})^3\right].\end{aligned}$
(51) When
$ \hat q\hat Q=0 $ , the electrostatic effect no longer exists, and the Schwarzschild result to the leading order is recovered. For fixed nonzero$ \hat Q $ and$ 0<|\hat q\hat Q|\ll 1 $ , the first order term$ (\hat q\hat Q)^1\hat Q^2 $ dominates the second order$ (\hat q\hat Q)^2 $ term. Comparing to neutral particles, we observe that a small electrostatic attractive (or repulsive) force due to small$ |\hat q| $ with$ {\rm{sign}}(\hat q\hat Q)=-1 $ (or$ {\rm{sign}}(\hat q\hat Q)= +1 $ ) increases (or decreases) the PS, similar to the effect of larger (or smaller) M to the PS of neutral particles in Schwarzschild spacetime. However, when$ |\hat q| $ exceeds$ |\hat Q|/2 $ , Eq. (51) indicates that the$ \hat q^2\hat Q^2 $ term becomes larger than$ \hat q\hat Q^3 $ and, therefore, the PS will increase as$ |\hat q| $ increases. In Fig. 3, we plot the dependence of the PS on$ \hat q $ for three typical$ \hat Q $ values ($ \hat Q= $ 1/10, 1/5, 2/5) using the red curves. For very small$ \hat Q $ , the dependence is visually suppressed by the largeness of the PS at large$ \hat q $ and$ \hat Q $ . If we fix$ \hat q $ as nonzero and vary$ \hat Q $ to but still maintain$ 0<|\hat q\hat Q|\ll 1 $ , as$ \hat Q $ deviates from zero, the PS (51) to the leading orders of$ \hat Q $ will behave asFigure 3. (color online) Dependence of
$ \alpha_{\rm{RN,PN}} $ for small$ \hat q\hat Q $ . The red curves correspond to$ \hat Q=1/10,\,1/5,\,2/5 $ and the solid and dashed blue curves correspond to$ \hat{q}=-7/10,2/5 $ and$ \hat{q}=3/2 $ , respectively. We select$ p=20M,\; e=4/5,\; |\hat{Q}|<2/5,\; $ $ |\hat{q}|<5/3 $ for the qualitative features to be visible.$ \alpha_{\rm{RN,PN}}= \frac{\pi M}{p}\left[6-(1-\hat q^2)\hat{Q}^2+{\cal{O}}(\hat Q)^3\right]+{\cal{O}}\left[\frac{M}{p},(\hat{q}\hat{Q})^3\right]. $
(52) This means that, in this region, the increase in
$ \hat Q $ decreases (or increases) the PS if$ |\hat q| $ is smaller (or larger) than 1. This dependence on$ \hat Q $ is shown by the blue curves in Fig. 3.In Fig. 4, we plot the orbits for several choices of
$ \hat q $ and$ \hat Q $ to illustrate the dependence of the PS on them. The starting points of these orbits are all set at the positive x-axis, which are also their apoapses. We observe that the PS decreases first and then increases as$ \hat q $ increases for fixed$ \hat Q $ . For fixed$ |\hat q|<1 $ (or$ |\hat q|>1 $ ), it decreases (or increases) as$ |\hat Q| $ increases. These agree perfectly with the analysis above. -
The EMD theory is a captivating theoretical framework that combines general relativity, electromagnetism, and the dilaton field minimally, providing a useful starting point for unifying all the fundamental interactions. For the static and spherical solution in this theory, its line element and electric field are described by [32, 33]
$ {\rm{d}}s^{2}=-f_+f_-^{\gamma}{\rm{d}}t^{2}+f_+^{-1}f_-^{-\gamma}{\rm{d}}r^{2}+r^{2}f_-^{1-\gamma}{\rm{d}}\Omega^{2}, $
(53a) $ f_{\pm}=1-\frac{r_{\pm}}{r},\quad\gamma=\frac{1-\eta_{1}\lambda^{2}}{1+\eta_{1}\lambda^{2}}, $
(53b) $ {\cal{A}}_0(r)=-\frac{Q}{r}, $
(53c) where
$ r_\pm $ represents the locations of the outer and inner event horizons, respectively,$ \eta_1=\pm1 $ depending on the dilaton/antidilaton nature of the scalar field, and λ is the real (anti)dilaton-Maxwell coupling constant. These parameters can be linked to the mass and charge of the spacetime by$ 2M=r_+ +\gamma r_-,\quad2Q^2=\eta_2(1+\gamma)r_+r_- $
(54) or equivalently
$ \begin{aligned}[b] r_+=\; & M+ \sqrt{ M^2-\frac{2\eta_2 \gamma Q^2}{1+\gamma}}, \\ r_-=\;&\frac{1}{\gamma}\left( M- \sqrt{M^2-\frac{2\eta_2 \gamma Q^2}{1+\gamma}}\right).\end{aligned} $
(55) Here,
$ \eta_2=\pm1 $ for Maxwell and anti-Maxwell fields, respectively. When$ \gamma=\eta_2=1 $ , metric (53a) reduces to the normal RN spacetime.Substituting the asymptotic expansion coefficients of the metric and the electric potential functions in Eq. (53a) into Eq. (40), we determine the PS of charged test objects in the PN approximation in the EMD gravity to the leading order of
$ 1/p $ as$ \begin{aligned}[b]\alpha_{{\rm{EMD,PN}}}=& \pi \Bigg\{ 6 -\eta_2 \hat{Q}^2- \hat{q} \hat{Q}\\&\times\left[6+\frac{\gamma-1}{\gamma}\left(1-\sqrt{1-\frac{2\eta_2\gamma \hat{Q}^2}{ \left(\gamma+1\right)}}\right)\right] \\&+\hat{q}^2 \hat{Q}^2 \Bigg\}\frac{M}{(1-\hat{q}\hat{Q} )p}+{\cal{O}}\left(\frac{M}{p}\right)^2. \end{aligned} $
(56) The Maxwell-gravity coupling constant
$ \eta_2 $ also determines the sign of its contribution to the PS through gravity (the$ \eta_2\hat{Q}^2 $ term). This sign choice also affects the PS through the electric interaction term proportional to$ \hat{q}\hat{Q} $ , but only weakly when$ \hat{Q} $ is small. Parameter γ, which is related to the strength λ of the Maxwell-(anti-)dilaton coupling, determines the amount of deviation of the electric interaction from the standard RN case. In limits$ \eta_2=\gamma=1 $ , this PS agrees with that of the RN case in Eq. (48).We observe that, by inspecting the higher order terms of the PS (56), the coefficient of order
$ 1/p^{-n} $ contains a factor of$ \eta_2\left(1+\gamma\right) $ up to order$- {n}/{2}$ . Therefore, for the total PS to converge, condition$ p>1/\sqrt{|1+\gamma|} $ should be fulfilled. In Fig. 5, we show the effect of Maxwell-(anti)dilaton parameter γ on the PS for different signs of$ \hat{q} $ and$ \eta_2 $ , with all other continuous parameters fixed. For the region of γ, we observe that the PSs are well defined for both$ \eta_2 $ values, and the effects of γ are always opposite for different$ \eta_2 $ and signs of$ \hat q $ . Near the location of$ \gamma=1 $ , the PS decreases (or increases) as γ increases for normal Maxwell field, i.e.,$ \eta_2=+1 $ (or anti-Maxwell field i.e.$ \eta_2=-1 $ ), for positive$ \hat q $ . This means that stronger Maxwell-dilaton coupling will result in a smaller (or larger) PS if the electromagnetic field is Maxwell (or anti-Maxwell) for positive particle charges. In comparison, the effect of γ when$ \gamma\lesssim 2 $ is the opposite to its effect when$ \gamma\gtrsim 1 $ . When$ |\gamma|\gg 1 $ , then its effect on the PS diminishes as$ |\gamma| $ increases, as shown by Eq. (56).Figure 5. (color online) Dependence of
$ \alpha_{\rm{EMD,PN}} $ on γ with$ p=500M,\; e=1/10,\; \hat{Q}=1/2,\; \hat{q}=\pm 1/5 $ . For$ r_+,\; r_- $ in Eq. (55) to exist, we require$ \gamma<-2 $ for$ \eta_2=+1 $ and$ \gamma>-2/3 $ for$ \eta_2=-1 $ . The blue and red curves correspond to$ \eta_2=+1 $ and$ \eta_2=-1 $ , respectively. -
A Lorentzian traversable wormhole with a charge is described by the line element (1) with metric functions and field strength [34]:
$ \begin{aligned}[b] A(r)=\;&\left(1+\frac{Q^2}{r^2}\right), \\ B(r)=\;&\left(1-\frac{s(r)}{r}+\frac{Q^2}{r^2}\right)^{-1}, \end{aligned}$
(57a) $ \begin{aligned}[b] C(r)=\;&r^2, \\ {\cal{F}}_{tr}=\;&-\frac{Q}{r^{2}}\sqrt{A(r)B(r)},\end{aligned} $
(57b) where
$ s(r) $ is the shape function,$ s(r)=s_{0}^{\frac{2\beta}{2\beta+1}}r^{\frac{1}{2\beta+1}}, $
(58) and
$ \beta<-1/2 $ . For the PN method to be effective, exponent$ \dfrac{2\beta}{2\beta+1} $ must be an integer. Therefore, we fix$ \beta=-1 $ and, consequently,$ s(r)=s_0^2/r $ . This implies that$ s_0 $ plays the role of an additional charge, although only partially, through the metric function$ B(r) $ . Parameter$ s_0 $ in this spacetime must satisfy$ s_0^2>Q^2 $ to maintain the wormhole.The asymptotic form of the metric functions (57b) can be obtained easily. For potential
$ {\cal{A}}_0 $ , we can obtain its asymptotic expansion by expanding and then integrating (57d):$ {\cal{A}}_0(r)=-\frac{Q}{r}-\frac{Qs_0^2}{6 r^3}+\frac{Qs_0^2\left(4Q^2-3s_0^2\right)}{40r^5}+{\cal{O}}\left(\frac{Q}{r}\right)^7. $
(59) Substituting these expansions into Eq. (40), the PS in this charged wormhole spacetime becomes
$\begin{aligned}[b] \alpha_{{\rm{CW,PN}}}=\;&\pi\left(\frac{1}{\hat{q}}-\hat{q}\right)\frac{Q}{p}-\frac{\pi}{4}\left[24+4e^2+\frac{1}{\hat{q}^2}-\hat{q}^2\left(1-2e^2\right)\right.\\ &\left.-2\frac{s_{0}^{2}}{Q^2}\left(4+e^2\right)\right]\frac{Q^2}{p^2}+{\cal{O}}\left(\frac{Q}{p}\right)^3. \end{aligned} $
(60) When
$ s_0 $ is set to zero, this agrees with the PS (48) in the RN case with$ M=0 $ . Therefore, the PS in this spacetime can be considered as that of a massless RN spacetime with some modifications from parameter$ s_0 $ . Moreover, the fundamental scale in this spacetime can be selected as the charge Q.Among the two terms in the leading order, the term proportional to
$ Q/\hat{q} $ is due to the gravitational effect of the spacetime charge, whereas the term proportional to$ -\hat{q}Q $ is still caused by the electric interaction. The extra parameter$ s_0 $ only appears starting from the$ (Q/p)^2 $ order and affects the PS only through gravitation (no multiplication with$ \hat{q} $ ). This also agrees with the fact that$ s_0 $ does not appear in the leading order of the electric potential in Eq. (59). Similar to the previous analysis after Eq. (48) for the convergence requirement in the RN spacetime, the$ \alpha_{{\rm{CW,PN}}} $ here also consists of four distinct series and the convergence of the PS requires$ |Q|<p,\; \frac{|Q|}{p}\leq |\hat{q}| \leq \frac{p}{|Q|},\; |s_0|<p. $
(61) Additionally, differentiating Eq. (3c) with respect to τ again, we find that, in this spacetime,
$ \begin{aligned}[b] \ddot{r}=\;&\frac{\hat{q}Q}{r^4\left(Q^2+r^2\right)}\sqrt{\frac{Q^2+r^2}{Q^2+r^2-s_0^2}}\left[\left(Q^2+r^2\right)^2-Q^2 s_0^2\right]\dot{t} \\ &+\frac{Q^2 \left(Q^2+r^2-s_0^2\right)\dot{t}^2}{r^5}+\frac{\left(s_0^2-Q^2\right)\dot{r}^2}{r\left(Q^2+r^2-s_0^2\right)} \\ &+\frac{\left(Q^2+r^2-s_0^2\right)\dot{\phi}^2}{r}. \end{aligned} $
(62) Because
$ Q^2+r^2-s_0^2>0 $ due to$ B(r)>0 $ in the observable region, the above equation shows that if$ qQ>0 $ ,$\ddot{r}$ will always be positive, and no bounded orbits or well-defined PS exists. Therefore, we also require that$ qQ<0 $ in this spacetime to study its PS. This is also consistent with the instinct that when no gravitational attraction exists from the mass ($ M=0 $ ),$ qQ $ must be negative for the electrostatic interaction to be attractive and form a bound orbit.In Fig. 6, we plot the dependence of the PS on new parameter
$ s_0 $ (in scale of Q) according to Eq. (60) for a few values of$ \hat q $ . For all$ \hat q $ , the PS increases as$ s_0/Q $ increases. This is generally consistent with the observation in the RN spacetime in Fig. 3 that as charge$ Q^2 $ increases, the PS also increases. For each$ -1\lesssim\hat q<0 $ , a value of$ s_0/Q $ below which the PS becomes negative exists. However, this characteristic is not present in the regular RN case because the PS is mostly positive for reasonable parameter choices. This critical value of$ s_0 $ can also be solved from Eq. (60) to the leading two orders asFigure 6. (color online) Dependence of
$ \alpha_{\rm{CW,PN}} $ on$ s_0 $ with fixed$ Q=1,\; p=20Q,\; e=1/2 $ .$ s_{0}^2 = \frac{2Q^2}{4+e^2}\left(\hat{q}-\frac{1}{\hat{q}}\right)\frac{p}{Q}+{\cal{O}}\left(Q\right)^2. $
(63) From this, we observe that for positive Q, because
$ qQ $ must be negative, the existence of such critical$ s_0 $ requires that$ 1\lesssim\hat q<0 $ , which is consistent with the observations in Fig. 6.In Fig. 7, we use
$ \hat q=-1/3 $ , which corresponds to the red dotted curve in Fig. 6, and three typical values of$ s_0/Q\; (2, 7, 8) $ to plot the corresponding orbits. The first of these according to Fig. 6 corresponds to the orbit with a negative PS, whereas the latter two have positive PSs. Fig. 7 verifies this, and we have checked that the PS in these orbits agrees quantitatively with the value in Fig. 6 and Eq. (60). -
In this section, we apply the PS in the RN spacetime to the observed PS of Mercury around the Sun and S2 around Sgr A* to constrain the charges of these objects.
-
The MESSENGER spacecraft measured an uncertainty of
$ 9\times 10^{-4}\; ''/ $ cty for the Schwarzschild-like precession of Mercury around the Sun [35]. If we associate this uncertainty to the electric effect of solar charge$ \hat Q_\odot $ and Mercury charge$\hat q$ ☿, then this uncertainty can be used to constrain the charges in the parameter space spanned by them. The result is shown in Fig. 8, where the allowed region of$ \hat q$ ☿ and$ \hat Q_\odot $ are enclosed by four segments of blue boundaries and two vertical red boundaries. The blue boundaries are due to the leading order of Eq. (48) because, in this case, p☿/$ M_{\odot}=3.8\times 10^7 $ , which makes the first order sufficient for the estimation of the PS. In other words, we treat this uncertainty as$ \begin{aligned}[b] \Delta\alpha =\;&\pi \left( 6 -\hat{Q}^2-6 \hat{q} \hat{Q} +\hat{q}^2 \hat{Q}^2 \right)\frac{M}{(1-\hat{q}\hat{Q} )p}-\frac{6\pi M}{p} \\ =\;&\left(\hat{q}^2 -1\right)\hat{Q}^2 \frac{\pi M}{(1-\hat{q}\hat{Q} )p},\end{aligned} $
(64) where, in the last term of the first line, we deduce the Schwarzschild contribution. Therefore, we observe that these blue boundaries are oddly symmetric, meaning that the change of signs
$ (\hat q\to-\hat q,\,\hat Q\to -\hat Q) $ will not change the extra PS. Moreover, we observe that these boundaries are generally not evenly symmetric, i.e., when only one charge changes sign, the value of$ \Delta\alpha $ changes. However, in the current case, since the extra PS is very small, the total value of the allowed$ \hat q\hat Q $ is already very small, and this nonsymmetry cannot be recognized in the plot using bare eyes. The red boundaries on the solar charge at two ends of$ \hat Q_\odot $ are obtained from the tight constraint in Ref. [36].The constrained value of the solar charge is much tighter than the previously reported value using the PS data, if we assume the same charge for Mercury (see Eqs. (42) and (43) of Ref. [24]).
-
The precession of S2 around Srg A
$ ^{*} $ was measured recently by the Gravity group [9] to yield ratio f of the measured value to the standard Schwarschild-precession. A value of$ f=1.10\pm 0.19 $ was obtained. If we associate the deviation of the central value from 1 to the electric effect of the Sgr A* and S2 charges, then similar to the Mercury precession, we can use this deviation to determine the allowed region in the parameter space of$ \hat q_{\rm{S2}} $ and$ \hat Q_{\rm{Sgr\ A}^*} $ .The result is shown in Fig. 9, where the blue boundaries indicate the leading order of Eq. (48). In this case,
$ p/M_{\rm{Sgr\ A}^*}=5\times 10^3 $ [40] is sufficiently large such that the first order of the PS is sufficient for estimating$ \hat q_{\rm{S2}} $ and$ \hat Q_{\rm{Sgr\ A}^*} $ . In this allowed parameter space, because the maximal$ \hat q\hat Q $ value can be as large as$ 0.8 $ , from Eq. (64) and comparing to the plot in Fig. 8, its nonsymmetry under sign change$ (\hat q\to-\hat q,\,\hat Q\to\hat Q) $ or$ (\hat q\to\hat q,\,\hat Q\to-\hat Q) $ is very clear. The red boundary on the Sgr A* charge is due to the constraint from the shadow size of the Sgr A* [41] (see [42–44] for other bounds on this charge). To our best knowledge, the allowed region in Fig. 9 is the first combined constraint on charge$ \hat q_{\rm{S2}} $ and$ \hat Q_{\rm{Sgr\ A}^*} $ using the S2 PS data. -
For a minimally coupled EMS gravity, the spacetime metric and electric field are described by [47]
$ \begin{aligned}[b] A(r)=\;&\frac{1}{B(r)}=\left[\frac{r_+\left(\dfrac{r-r_-}{r-r_+}\right)^{n/2}-r_-\left(\dfrac{r-r_+}{r-r_-}\right)^{n/2}}{r_+-r_-}\right]^{-2},\\C(r)=\;&A^{-1}(r)\left(r-r_+\right)\left(r-r_-\right),\end{aligned} $
$ \begin{aligned}[b] {\cal{F}}_{tr}=\;&-\frac{Q}{C(r)}, \end{aligned} $
(A1) where
$ n\in(1/2,1] $ , and$ r_{\pm}=\frac{M\pm\sqrt{M^{2}-Q^{2}}}{n} $
(A2) are the location of the outer and inner event horizons, and n is related to the scalar charge.
Substituting the asymptotics of these functions into Eq. (40), the PS in the EMS field in the PN approximation is found as
$ \begin{aligned}[b] \\[-8pt]\alpha_{{\rm{EMs,PN}}}=\;&\pi \left( 6 -\hat{Q}^2-6 \hat{q} \hat{Q} +\hat{q}^2 \hat{Q}^2 \right)\frac{M}{(1-\hat{q}\hat{Q} )p}+\frac{\pi }{4 n^2}\left[\left(92n^2+24n-8+e^2 \left\{-16n^2+24n-2\right\}\right.\right. \\ &\left.\left.- \hat{Q}^2 \left\{4\left[13n^2+n-2\right]+2e ^2\left[2n-1\right]\right\}-n^2\hat{Q}^4\right)-2 \hat{q} \hat{Q} \left(92n^2+24n-8 \right.\right. \\ &\left.\left.-2e^2\left\{8n^2-12n+1\right\}- \hat{Q}^2 \left\{43n^2+2n-8-e^2\left[3n^2+2n-2\right] \right\}\right)\right. \\ &\left.+2 \hat{q}^2 \hat{Q}^2 \left(56n^2+14n-4-e^2\left\{12n^2-14n+1-\hat{Q}^2\left[16n^2-4+e^2\left(3n^2-1\right)\right]\right\} \right)\right. \\ &\left.-2 \hat{q}^3 \hat{Q}^3 \left(11n^2+2n-e^2\left\{5n^2-2n\right\}\right)+\hat{q}^4 \hat{Q}^4 n^2\left(1-2 e ^2\right)\right]\left[\frac{M}{ (1-\hat{q}\hat{Q} )p}\right]^2+{\cal{O}}\left(\frac{M}{p}\right)^3. \end{aligned} $ (A3) To the leading order, this agrees with the RN result (48). Therefore, distinguishing this spacetime from the RN one using the current observation data would be difficult.
-
The charged Horndeski spacetime is described by the metric and potential functions [48, 49]:
$ A(r)=1-\frac{2M}r+\frac{Q^2}{4r^2}-\frac{Q^4}{192r^4}, $
(A4a) $ B(r)=\frac{1}{A(r)}\left(1-\frac{Q^2}{8r^2}\right)^2, C(r)=r^2, $
(A4b) $ {\cal{A}}_0(r)=-\frac {Q}{r}+\frac{Q^3}{24r^3}. $
(A4c) The expansion of these functions are easy. Subsequently, substituting this into Eq. (40), the PN PS in this spacetime is determined as
$ \begin{aligned}[b] \alpha_{{\rm{CH,PN}}} =& \frac{\pi}{4} \left( 24 -\hat{Q}^2-24 \hat{q} \hat{Q} +4\hat{q}^2 \hat{Q}^2 \right)\frac{M}{\left(1-\hat{q}\hat{Q} \right)p} \\ &+\frac{\pi }{64 }\left\{\left[96 \left(18+e ^2\right)-16 \hat{Q}^2 \left(13+e ^2\right)-\hat{Q}^4\right]\right. \\ &\left.-8 \hat{q} \hat{Q} \left[24 \left(18+e ^2\right)- \hat{Q}^2 \left(43+5 e ^2\right)\right] \right.\end{aligned}$
$ \begin{aligned}[b] &\left.+8 \hat{q}^2 \hat{Q}^2 \left[4\left(66+e ^2\right)- \hat{Q}^2 \left(16+3e ^2\right)\right]\right.\\ &\left. -32 \hat{q}^3 \hat{Q}^3 \left(13-3 e ^2\right) +16\hat{q}^4 \hat{Q}^4 \left(1-2 e ^2\right) \right\} \\ & \times \left[\frac{M}{ (1-\hat{q}\hat{Q} )p}\right]^2+{\cal{O}}\left(\frac{M}{p}\right)^3 . \end{aligned} $
(A5) To the leading order, this is different from the RN result in the gravitational contribution of
$ \hat{Q} $ by a factor of a quarter. -
The charge black-bounce RN spacetime is described by the RN metric functions with r replaced by
$ \sqrt{r^2+l^2} $ , where l is a length scale (usually associated with the Planck length), i.e., the new metric functions. The electric potential becomes [50]$ A(r)= \frac{1}{B(r)}=1-\frac{2M}{\sqrt{r^2+l^2}}+\frac{Q^2}{r^2+l^2}, $
(A6a) $ C(r)= r^2+l^2, $
(A6b) $ {\cal{A}}_0(r)= -\frac{Q}{\sqrt{r^2+l^2}}. $
(A6c) The asymptotic expansions of these functions are also simple, and after substituting the coefficients into Eq. (40), we can directly obtain the PS in this spacetime. To the second order, this is
$ \alpha_{{\rm{BB,PN}}}= \alpha_{{\rm{RN,PN}}}+\frac{\pi \left(2+e^2\right)l^2}{2p^2}+{\cal{O}}\left(\frac{M}{p}\right)^3, $
(A7) where the
$ \alpha_{{\rm{RN,PN}}} $ is the PS of the RN spacetime in Eq. (48). The extra scale parameter l only appears from the second order of$ 1/p $ ; therefore, the PS in this spacetime is also difficult to distinguish from that in the regular RN spacetime. Moreover, to the$ {\cal{O}}(M/p)^2 $ order, the scale parameter l does not participate in any electric interaction effect to the PS. We also checked that this will change beginning from the$ {\cal{O}}(M/p)^3 $ order, i.e., terms proportional to$ l^2\hat q\hat Q $ and$ l^2\hat q^2\hat Q^2 $ exist in the PS at this order. The result (A7), after setting$ \hat q=0 $ , reduces to the PS of neutral particles in black-bounce RN spacetime [23].
Periapsis shift in spherically symmetric spacetimes and effects of electric interactions
- Received Date: 2024-03-21
- Available Online: 2024-08-15
Abstract: The periapsis shift of charged test particles in arbitrary static and spherically symmetric charged spacetimes are studied. Two perturbative methods, the near-circular approximation and post-Newtonian methods, are developed and shown to be very accurate when the results are determined to high orders. The near-circular approximation method is more precise when eccentricity e of the orbit is small, whereas the post-Newtonian method is more effective when orbit semilatus rectum p is large. Results from these two methods are shown to agree when both e is small and p is large. These results are then applied to the Reissner-Nordström spacetime, the Einstein-Maxwell-dilation gravity, and a charged wormhole spacetime. The effects of various parameters on the periapsis shift, particularly that of the electrostatic interaction, are carefully studied. The periapsis shift data of the solar-Mercury are then used to constrain the charges of the Sun and Mercury, and the data of the Sgr A*-S2 periapsis shift are used to determine, for the first time using this method, the constraints of the charges of Sgr A* and S2.