Post-Keplerian perturbations of the orbital time shift in binary pulsars: an analytical formulation with applications to the galactic center

We develop a general approach to analytically calculate the perturbations Δδτp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \delta \tau _\text {p}$$\end{document} of the orbital component of the change δτp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta \tau _\text {p}$$\end{document} of the times of arrival of the pulses emitted by a binary pulsar p induced by the post-Keplerian accelerations due to the mass quadrupole Q2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q_2$$\end{document}, and the post-Newtonian gravitoelectric (GE) and Lense–Thirring (LT) fields. We apply our results to the so-far still hypothetical scenario involving a pulsar orbiting the supermassive black hole in the galactic center at Sgr A∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^*$$\end{document}. We also evaluate the gravitomagnetic and quadrupolar Shapiro-like propagation delays δτprop\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta \tau _\text {prop}$$\end{document}. By assuming the orbit of the existing main sequence star S2 and a time span as long as its orbital period Pb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_\mathrm{b}$$\end{document}, we obtain ΔδτpGE≲103s,ΔδτpLT≲0.6s,ΔδτpQ2≲0.04s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| \Delta \delta \tau _\text {p}^\text {GE}\right| \lesssim 10^3~\text {s},~\left| \Delta \delta \tau _\text {p}^\text {LT}\right| \lesssim 0.6~\text {s},\left| \Delta \delta \tau _\text {p}^{Q_2}\right| \lesssim 0.04~\text {s}$$\end{document}. Faster Pb=5years\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( P_\mathrm{b}= 5~\text {years}\right) $$\end{document} and more eccentric e=0.97\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( e=0.97\right) $$\end{document} orbits would imply net shifts per revolution as large as ΔδτpGE≲10Ms,ΔδτpLT≲400s,ΔδτpQ2≲103s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| \left\langle \Delta \delta \tau _\text {p}^\text {GE}\right\rangle \right| \lesssim 10~\text {Ms},~\left| \left\langle \Delta \delta \tau _\text {p}^\text {LT}\right\rangle \right| \lesssim 400~\text {s},\left| \left\langle \Delta \delta \tau _\text {p}^{Q_2}\right\rangle \right| \lesssim 10^3~\text {s}$$\end{document}, depending on the other orbital parameters and the initial epoch. For the propagation delays, we have δτpropLT≲0.02s,δτpropQ2≲1μs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| \delta \tau _\text {prop}^\text {LT}\right| \lesssim 0.02~\text {s},~\left| \delta \tau _\text {prop}^{Q_2}\right| \lesssim 1~\upmu \text {s}$$\end{document}. The results for the mass quadrupole and the Lense–Thirring field depend, among other things, on the spatial orientation of the spin axis of the Black Hole. The expected precision in pulsar timing in Sgr A∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^*$$\end{document} is of the order of 100μs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$100~\upmu \text {s}$$\end{document}, or, perhaps, even 1–10 μs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu \text {s}$$\end{document}. Our method is, in principle, neither limited just to some particular orbital configuration nor to the dynamical effects considered in the present study.


Introduction
In a binary hosting at least one emitting pulsar 1 p, the time of arrivals τ p of the emitted radio pulses changes primarily because of the orbital motion about the common center of 1 See, e.g., [28,49] and references therein. a e-mail: lorenzo.iorio@libero.it mass caused by the gravitational tug of the unseen companion c, which can be, in principle, either a main sequence star or an astrophysical compact object like, e.g., another neutron star which does not emit or whose pulses are, for some reason, not detectable, a white dwarf or, perhaps, even a black hole [50]. Such a periodic variation δτ p ( f ) can be modeled as the ratio of the projection of the barycentric orbit r p of the pulsar p onto the line of sight to the speed of light c [12,31]. By assuming a coordinate system centered in the binary's center of mass whose reference z-axis points toward the observer along the line of sight in such a way that the reference {x, y} plane coincides with the plane of the sky, we have (1) In obtaining Eq. (1), which is somewhat the analogous of the range in Earth-Moon or Earth-planets studies [12], we used the fact that, to the Keplerian level, the barycentric semimajor axis of the pulsar A is In a purely Keplerian scenario, there is no net variation δτ p over a full orbital cycle.
In this paper, we illustrate a relatively simple and straightforward approach to analytically calculate the impact that several post-Keplerian (pK) features of motion, both Newtonian (quadrupole) and post-Newtonian (1pN static and stationary fields), have on such a key observable. As such, we will analytically calculate the corresponding net time delays per revolution δτ p ; the instantaneous shifts δτ p ( f ) will be considered as well in order to cope with systems exhibiting very long orbital periods with respect to the time spans usually adopted for data collection. Our strategy has gen-eral validity since, in principle, it can be extended to a wide range of dynamical effects, irrespectively of their physical origin, which may include, e.g., modified models of gravity as well. Furthermore, it is applicable to systems whose constituents may have arbitrary masses and orientations of their spin axes, and orbital configurations. Thus, more realistic sensitivity analyses, aimed at both reinterpreting already performed studies and designing future targeted ones, could be conducted in view of a closer correspondence with which is actually measured. We will also take into account the Shapiro-like time delays due to the propagation of the electromagnetic waves emitted by the visible pulsar(s) throughout the spacetime deformed by axisymmetric departures from spherical symmetry of the deflecting bodies [11,14,15,29,30,32,50,55].
Our results, which are not intended to replace dedicated, covariance-based real data analyses, being, instead, possible complementary companions, will be applied to the so far putative scenario involving emitting radiopulsars, not yet detected, orbiting the Supermassive Black Hole (SMBH) in the Galactic Center (GC) at Sgr A * [1, 17,20,27,33,39,54]. Moreover, we will perform also quantitative sensitivity analyses on the measurability of frame-dragging and quadrupolar-induced time delays in such a hypothesized system. In principle, our results may be applicable even to anthropogenic binaries like, e.g., those contrived in past concept studies to perform tests of fundamental physics in space [37,43], or continuously emitting transponders placed on the surface of some moons of larger astronomical bodies.
The paper is organized as follows. Section 2 details the calculational approach. The 1pN Schwarzschild-type gravitoelectric effects are calculated in Sect. 3, while Sect. 4 deals with the 1pN gravitomagnetic ones. The impact of the quadrupole mass moment of the SMBH is treated in Sect. 5. Section 6 summarizes our findings.

Outline of the proposed method
If the motion of a binary is affected by some relatively small post-Keplerian (pK) acceleration A, either Newtonian or post-Newtonian (pN) in nature, its impact on the projection of the orbit onto the line of sight can be calculated perturbatively as follows. Casotto [10] analytically worked out the instantaneous changes of the radial, transverse and out-of-plane components, r ρ , r σ , r ν of the position vector r, respectively, for the relative motion of a test particle about its primary: they are where the time derivatives dκ/dt of the Keplerian orbital elements κ are to be taken from the right-hand sides of the Gauss equations, evaluated onto the Keplerian ellipse r = p 1 + e cos f (12) and assumed as unperturbed reference trajectory; the same holds also for entering Eq. (6). The case of the mean anomaly M is more subtle and requires more care. Indeed, in the most general case encompassing the possibility that the mean motion n b is time-dependent because of some physical phenomena, it can be written as 2 [4,8,35] the Gauss equation for the variation of the mean anomaly at epoch is 3 [4,8,35] If n b is constant, as in the Keplerian case, Eq. (14) reduces to the usual form, In general, when a disturbing acceleration is present, the semimajor axis a does vary according to Eq. (7); thus, also the mean motion n b experiences a change 4 which can be calculated in terms of the true anomaly f as by means of Eqs. (7) and (13). Depending on the specific perturbation at hand, Eq. (18) does not generally vanish. Thus, the total change experienced by the mean anomaly M due to the disturbing acceleration A can be obtained: where In the literature, the contribution due to Eq. (21) has been often neglected. An alternative way to compute the perturbation of the mean anomaly with respect to Eq. (19) implies the use of the mean longitude λ and the longitude of pericenter . It turns out that 5 [46] 3 It is connected with the Gauss equation for the variation of the time of passage at pericenter t p by dη/dt = −n b dt p /dt. 4 We neglect the case μ (t). 5 The mean longitude at epoch is denoted as by Milani et al. [35], Soffel [46], Brumberg [8] and Bertotti et al. [4]. It is better suited than η at small inclinations [4].
where the Gauss equations for the variation of , are [4,8,35,46] d dt = 2 sin 2 I 2 d dt It must be remarked that, depending on the specific perturbing acceleration A under consideration, the calculation of Eq.
expressing the perturbation experienced by the binary's relative position vector r. It is It is possible to express the true anomaly as a function of time through the mean anomaly according to Brouwer and Clemence [7, p. 77] where J k (se) is the Bessel function of the first kind of order k and s max , j max are some values of the summation indices s, j adequate for the desired accuracy level. Having at disposal such analytical time series yielding the time-dependent pattern of Eq. (26) allows one to easily study some key features of it such as, e.g., its extrema along with the corresponding epochs and the values of some unknown parameters which may enter the disturbing acceleration. The net change per orbit r z can be obtained by calculating Eq. (26) with f = f 0 + 2π, and using Eq. (6) and Eqs. (19)- (21) integrated from f 0 to f 0 + 2π. In order to have the change of the times of arrival of the pulses from the binary's pulsar p, Eq. (26) and its orbit averaged expression have to be scaled by m c m −1 tot c −1 . In the following, we will look at three pK dynamical effects: the Newtonian deviation from spherical symmetry of the binary's bodies due to their quadrupole mass moments, and the velocity-dependent 1pN static (gravitoelectric) and stationary (gravitomagnetic) accelerations responsible of the time-honored anomalous Mercury perihelion precession and the Lense-Thirring frame-dragging, respectively.

The 1pN gravitoelectric effect
Let us start with the static component of the 1pN field, which, in the case of our Solar System, yields the formerly anomalous perihelion precession of Mercury of˙ =

arcsec cty
The 1pN gravitoelectric, Schwarzschild-type, acceleration of the relative motion is, in General Relativity [46], By projecting Eq. (28) onto the radial, transverse, out-ofplane unit vectorsρ,σ ,ν, its corresponding components are Here, we use the true anomaly f since it turns out computationally more convenient. The resulting net shifts per orbit of the osculating Keplerian orbital elements, obtained by integrating Eq. (6) and If, on the one hand, Eq. (33) (34) is an exact result in the sense that no a priori assumptions on e were assumed. It can be shown that, to the zeroth order in e, Eq. (34) is independent of f 0 . We will not explicitly display here the analytical expressions for the instantaneous changes κ GE ( f ) , κ = a, e, I, , ω, M GE ( f ) because of their cumbersomeness, especially as far as the mean anomaly is concerned. How- , ω can be found in Eqs. (A2.78b)-(A2.78d) of Soffel [46, p. 178]. Equations (A2.78e) to (A2.78f) of Soffel [46, p. 178] allow one to obtain the instantaneous shift of the mean anomaly, although in terms of the three anomalies f, E, M; instead, our (lengthy) expression contains only the true anomaly f . See also Eqs.
It should be noted that Eq. (35) is independent of the semimajor axis a, depending only on the shape of the orbit through e and its orientation in space through I, ω. Furthermore, Eq. (35) does depend on the initial epoch t 0 through f 0 . In the limit e → 0, Eq. (35) does not vanish, reducing to In view of its cumbersomeness, we will not display here the explicit expression of δτ GE p ( f ) whose validity was successfully checked by numerically integrating the equations of motion for a fictitious binary system, as shown by Fig. 1; see also Sect. 3.1.
We will not deal here with the Shapiro-like propagation delay since it was accurately calculated in the literature; see, e.g., [11,14] and the references therein.

The pulsar in Sgr A * and the gravitoelectric orbital time delay
An interesting, although still observationally unsupported, scenario involves the possibility that radio pulsars orbit the SMBH at the GC in Sgr A * ; in this case, the unseen companion would be the SMBH itself. Thus, in view of its huge mass, the expected time shift per orbit δτ GE p would be quite large.
By considering a hypothetical pulsar with standard mass m p = 1.4 M and, say, the same orbital parameters of the main sequence star S2 actually orbiting the galactic SMBH [19], Eq. Lower row, red curve δτ GE p (t), in s, of a hypothetical pulsar in Sgr A * obtained analytically from Eqs. (26) to (27) and the instantaneous changes of the Keplerian orbital elements, not displayed in the text, induced by Eq. (28). It turns out that the net shift after a full revolution starting at t 0 = 2003.271 amounts to δτ GE p = 1,722.6948 s. The orbital configuration of the S2 star, quoted in Table 1, was adopted for the putative pulsar in Sgr A * δτ GE p max = 1.74521212562 Such huge orbital time delays would be accurately detectable, even by assuming a pessimistic pulsar timing precision of just 100 µs [20,39]; more optimistic views point towards precisions of the order of even 1-10 µs [20,39]. Table 1 Relevant physical and orbital parameters of the S2 star and the SMBH at the GC along with their estimated uncertainties according to Table 3 of [19]; they are referred to the epoch 2000.0. D 0 is the distance to Sgr A * . The Schwarzschild radius of the SMBH is r g . = 2G M • /c 2 = 0.088 au, while the linear size of the semimajor axis of S2 is a = 1,044 au = 11,863.6 r g . We quote also the derived values of the SMBH's angular momentum and quadrupole mass moment calculated as [18,21] [22,26,40]. The dimensionless parameter χ g ≤ 1 is of the order of about 0.6 for the SMBH in Sgr A * [39]. We display also the value f 0 inferred from Eq. (27) for the true anomaly at the epoch t 0 = 2003.271 quoted in Table 5

The 1pN gravitomagnetic Lense-Thirring effect
The stationary component of the 1pN field, due to massenergy currents, is responsible of several aspects of the socalled spin-orbit coupling, or frame-dragging [16,44,45,48]. The 1pN gravitomagnetic, Lense-Thirring-type acceleration affecting the relative orbital motion of a generic binary made of two rotating bodies A, B is [2,46] In general, it iŝ i.e. the angular momenta of the two bodies are usually not aligned. Furthermore, they are neither aligned with the orbital angular momentum L, whose unit vector is given byν. Finally, also the magnitudes S A , S B are, in general, different. The radial, transverse and out-of-plane components of the gravitomagnetic acceleration, obtained by projecting Eq. (42) onto the unit vectorsρ,σ ,ν, turn out to be It is interesting to remark that, in the case of Eq. (42), both Eqs. (15) and (21) yield vanishing contributions to M LT . For previous calculations based on different approaches and formalisms, see, e.g., [2,13,25], and the references therein. Equation (26) Note that, contrary to Eq. (35), Eq. (51) does depend on the semimajor axis as a −1/2 . As Eq. (35), also Eq. (51) depends on f 0 . The instantaneous orbital time shift δτ LT p ( f ) turns out to be too unwieldy to be explicitly displayed here. Its validity was successfully checked by numerically integrating the equations of motion for a fictitious binary system, as shown by Fig. 2; see also Sect. 4.1.
The gravitomagnetic propagation time delay is treated in Sect. 4.2.  Table 1, was adopted for the putative pulsar in Sgr A * along with i • = 20.9 • , ε • = 317.9 •

The pulsar in Sgr A * and the Lense-Thirring orbital time delay
Let us, now, consider the so-far hypothetical scenario of an emitting radio pulsar orbiting the SMBH in Sgr A * [1,20]. It turns out that, in some relevant astronomical and astrophysical binary systems of interest like the one at hand, the (scaled) angular momentum S A/B of one of the bodies is usually much smaller than the other one. Let us assume that the pulsar under consideration has the same characteristics of PSR J0737-3039A. By assuming [3,36] I NS 10 38 kg m 2 (52) for the moment of inertia of a neutron star (NS), the spin of PSR J0737-3039A is S A = 2.8 × 10 40 kg m 2 s −1 .
The angular momentum of a NS of mass m NS can also be expressed in terms of the dimensionless parameter χ g > 0 as [34] Thus, Eq. (53) implies for PSR J0737-3039A. Since for the galactic SMBH it is [21] S • = χ g M 2 • G c 9.68 × 10 54 kg m 2 s −1 (56) with 6 [39] we have Thus, in this case, the dominant contribution to Eq. (42) is due to the pulsar's companion c. As far as the orientation of the SMBH's spin is concerned, we model it aŝ The angles i • , ε • are still poorly constrained [5,6,53], so that we prefer to treat them as free parameters by considering their full ranges of variation Also in this case, we assume for our putative pulsar the same orbital parameters of, say, the S2 star. By using our analytical expression for δτ LT p (t), calculated with S p → 0 in view of Eq. (58), one gets

The Lense-Thirring propagation time shift
The gravitomagnetic propagation delay for a binary with relative separation r and angular momentum S of the primary is [32,50] According to Figure 2 of [50], their unit vector K 0 agrees with ourê z , since it is directed towards the Earth. About the spin axis of the primary, identified with a BH by [50], our i • coincides with their λ • . Instead, their angle η • is reckoned from our unit vectorl, i.e. it is as if [50] set = 0. On the contrary, our angle ε • is counted from the reference x direction in the plane of the sky, whose unit vectorê x , in general, does not coincide withl. Furthermore, [50] use the symbol i for the orbital inclination angle, i.e., our I . It is important to notice that, contrary to the orbital time delay of Eq. (51), Eq. (66) is a short-term effect in the sense that there is no net shift over one orbital revolution. It is also worth As far as the putative scenario of the pulsar in the GC is concerned, the emitting neutron star is considered as the source s of the electromagnetic beam delayed by the angular momentum of the SMBH. Thus, by calculating Eq. (66) for a S2-type orbit and with S = S • , it is possible to obtain , a further non-central relative acceleration arises; it is [2] 2r 4 3μ in which the first even zonal parameter J A/B 2 is dimensionless. In the notation of [2], their J A/B 2 parameter is not dimension-less as ours, being dimensionally an area because it corresponds to our J A/B 2 R 2 A/B . Furthermore, [2] introduce an associated dimensional quadrupolar parameter I A/B , having the dimensions of a moment of inertia, which is connected to our Thus, I A/B corresponds to the dimensional quadrupolar parameter Q A/B 2 customarily adopted when astrophysical compact objects like neutron stars and black holes are considered [34,51], up to a minus sign, i.e.
Thus, Eq. (71) can be written as Projecting Eq. (71) onto the radial, transverse and out-ofplane unit vectorsρ,σ ,ν provides us with − a 4 1 − e 2 4 3μ (1 + e cos f ) 4 A straightforward consequence of Eqs. (75)-(77) is the calculation of the net quadrupole-induced shifts per revolution of the Keplerian orbital elements by means of Eq. (6) and Eqs. (19)- (21), which turn out to be a J 2 = e J 2 = 0, Also Eq. (82), as Eq. (34) for the Schwarzschild-like 1pN acceleration, is a novel result which amends the incorrect formulas widely disseminated in the literature [9,25,41,47,52]; indeed, it turns out that, in the case of Eq. (71), Eq. (21) does not vanish when integrated over a full orbital revolution. Furthermore, contrary to almost all of the other derivations existing in the literature, Eq. (82) is quite general since it holds for a two-body system with generic quadrupole mass moments arbitrarily oriented in space, and characterized by a general orbital configuration. The same remark holds also for Eqs. (78)-(81); cf. the corresponding (correct) results by [25] in the case of a test particle orbiting an oblate primary. According to Eq. (26) and Eqs. (78)-(82), the net orbitlike time change of the pulsar p after one orbital revolution is  Table 1, was adopted for the putative pulsar in Sgr A * along with i • = 146.7 • , ε • = 148.8 • (51), on the other hand, it depends on the orbital semimajor axis through a −1 . As far as δτ J 2 p ( f ) is concerned, it will not be displayed explicitly because it is far too ponderous. Also in this case, a numerical integration of the equations of motion for a fictitious binary system, displayed in Fig. 3, confirmed our analytical result for the temporal pattern of δτ J 2 p ( f ); see also Sect. 5.1. The propagation time delay is dealt with in Sect. 5.2.

The pulsar in Sgr A * and the quadrupole-induced orbital time delay
A rotating NS acquires a non-zero quadrupole moment given by [34] the absolute values of the dimensionless parameter q < 0 ranges from 0.074 to 3.507 for a variety of Equations of State (EOSs) and m NS = 1.4 M ; cf. Table 4 of [34]. It is interesting to note that [34] find the relation where the parameter α of the fit performed by [34] depends on both the mass of the neutron star and the EOS used. According to Table 7 of [34], it is for some of the EOSs adopted by [34]. In the case of PSR J0737-3039A, Eq. (84) yields According to Eq. (55) and Eq. (85), we have As a consequence of the "no-hair" or uniqueness theorems [22,26,40], the quadrupole moment of a BH is uniquely determined by its mass and spin according to [18,21] so that the quadrupole of a hypothetical emitting neutron star p orbiting the SMBH in Sgr A * can be completely neglected with respect to the quadrupole of the latter one in any practical calculation. According to our analytical expression for δτ

The quadrupole-induced propagation time shift
The propagation delay δτ J 2 prop due to the quadrupole mass moment is rather complicated to be analytically calculated; see, e.g., [29,30,32,55]. No explicit expressions analogous to the simple one of Eq. (66) for frame-dragging exist in the literature. Here, we will obtain an analytical formula for δτ J 2 prop , which will be applied to the double pulsar and the pulsar-Sgr A * systems. The approach by [55] will be adopted by adapting it to the present scenario. In the following, the subscripts d, s, o will denote the deflector, the source, and the observer, respectively. In the case of, say, the double pulsar, d is the pulsar B, while s is the currently visible pulsar A; in the pulsar-Sgr A * scenario, d is the SMBH and s is the hypothetical pulsar p orbiting it. See Fig. 4 Fig. 4 Vectors connecting the source s and the deflector d at t emi , and the observer o at t rec . The symbols adopted differ from those used in, e.g., [55]. The origin O coincides with the binary's center of mass (C.M.) is the barycentric position vector of the source s at the time of emission t emi , is the barycentric position vector of the deflector d at t emi , is the relative position vector of the source s with respect to the deflector d at t emi . Thus, to the Newtonian order, it is where m s , m d are the masses of source and deflector, respectively. Furthermore, is the barycentric position vector of the observer o at the time of reception t rec , is the position vector of the observer o at t rec with respect to the deflector d at t emi , and is the position vector of the observer o at t rec with respect to the source s at t emi . With our conventions for the coordinate axes, it is where D is the distance of the binary at t emi from us at t rec , If the pulsar-SMBH in the GC is considered, the quadrupole Shapiro-type time delay is much smaller than the orbital time shift. Indeed, by using Eqs. (124)-(125) for a S2-type orbital configuration, it turns out that  127) should be hard to be detectable in view of the expected pulsar timing precision, even in the optimistic case of 1-10 µs [20,39]. If the orbital configuration of S2 is abandoned, letting I, , ω, f, i • , ε • freely vary within their full natural ranges, we get values which can reach the 100 µs level for 0 ≤ e ≤ 0.97, 5 years ≤ P b ≤ 16 years.

Summary and conclusions
In order to perform sensitivity studies, designing suitable tests and reinterpreting existing data analyses in a way closer to the actual experimental practice in pulsar timing, we devised a method to analytically calculate the shifts δτ A p experienced by the orbital component of the time changes δτ p of a binary pulsar p due to some perturbing post-Keplerian accelerations A: Schwarzschild, Lense-Thirring and mass quadrupole. We applied it to the still hypothetical scenario encompassing an emitting neutron star which orbits the supermassive black hole in Sgr A * ; its timing precision could reach 100 µs, or, perhaps, even 1-10 µs [20,39]. The main results of the present study are resumed in Table 2. By assuming a S2-like orbital configuration and a time span as long as its orbital period, the magnitude of the post-Newtonian Schwarzschild-type gravitoelectric signature can reach δτ GE p 10 3 s. The post-Newtonian Lense-Thirring gravitomagnetic and quadrupolar effects are much smaller, amounting to at most δτ LT p 0.6 s, δτ Among other things, we also explicitly calculated an analytical formula for the Shapiro-like time delay δτ prop due to the propagation of electromagnetic waves in the field Table 2 Maximum and minimum values for the orbital and propagation time shifts δτ p , δτ prop over a full orbital revolution due to the gravitoelectric (GE), gravitomagnetic (LT) and quadrupole (Q 2 ) effects for a hypothetical pulsar-Sgr A * scenario. While the orbital time delay δτ p is cumulative over the revolutions, the propagation shift δτ prop vanishes over one orbital period. For the putative pulsar orbiting the SMBH in the GC, the orbital configuration of the main sequence S2 star was adopted. The values of the system's physical and orbital parameters corresponding to the quoted maxima and minima are not reported here: see the text for details. Only the angular momentum S • and the quadrupole Q • 2 of the SMBH, playing the role of deflector d, were taken into account. The expected timing precision for a pulsar orbiting the galactic SMBH is about 100 µs, or, perhaps, even 1-10 µs [20,39]  of a spinning oblate body, which we applied to the aforementioned binary system. As far as the Lense-Thirring and the quadrupolar effects are concerned, the Shapiro-like time shifts δτ prop are, in general, much smaller than the orbital ones δτ p which, contrary to δτ prop , are cumulative. In the case of the pulsar-Sgr A * scenario, we find, for a S2-type orbit, that the Lense-Thirring propagation delay is as little as δτ LT prop 0.02 s, while the quadrupolar one is of the order of δτ Q 2 prop 1 µs, both depending on the spin orientation of the black hole. Removing the limitation to the S2 orbital configuration yields essentially similar values for δτ LT prop , δτ Q 2 prop , even for highly eccentric and faster orbits.
Finally, we remark that our approach is general enough to be extended to arbitrary orbital geometries and symmetry axis orientations of the binary's bodies, and to whatsoever disturbing accelerations. As such, it can be applied to other binary systems as probes for, say, modified models of gravity. In principle, also man-made binaries could be considered.