Erratum to: Is it possible to measure new general relativistic third-body effects on the orbit of Mercury with BepiColombo?

Recently, Will calculated an additional contribution to the Mercury’s precession of the longitude of perihelion ϖ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varpi $$\end{document} of the order of ϖ˙W≃0.22milliarcseconds per century\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\dot{\varpi }}_\text {W}\simeq 0.22~\text {milliarcseconds per century}$$\end{document} (mascty-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {mas cty}^{-1}$$\end{document}). It is partly a direct consequence of certain 1pN third-body accelerations entering the planetary equations of motion, and partly an indirect, mixed effect due to the simultaneous interplay of the standard 1pN pointlike acceleration of the primary with the Newtonian N-body acceleration, to the quadrupole order, in the analytical calculation of the secular perihelion precession with the Gauss equations. We critically discuss the actual measurability of the mixed effects with respect to direct ones. The current uncertainties in either the magnitude of the Sun’s angular momentum S⊙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_\odot $$\end{document} and the orientation of its spin axis S^⊙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\hat{S}}}_\odot $$\end{document} impact the precessions ϖ˙J2⊙,ϖ˙LT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\dot{\varpi }}_{J_2^\odot },~{\dot{\varpi }}_\text {LT}$$\end{document} induced by the Sun’s quadrupole mass moment and angular momentum via the Lense–Thirring effect to a level which makes almost impossible to measure ϖ˙W\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\dot{\varpi }}_\text {W}$$\end{document}, even in the hypothesis that it comes entirely from the aforementioned 1pN third-body accelerations. On the other hand, from the point of view of the Lense–Thirring effect itself, the mismodeled quadrupolar precession δϖ˙J2⊙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta {\dot{\varpi }}_{J_2^\odot }$$\end{document} due to the uncertainties in S^⊙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\hat{S}}}_\odot $$\end{document} corresponds to a bias of ≃9%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\simeq 9\%$$\end{document} of the relativistic one. The resulting simulated mismodeled range and range-rate times series of BepiColombo are at about the per cent level of the nominal gravitomagnetic ones.


Introduction
Recently, Will [18] calculated a new general relativistic contributioṅ W 0.22 mas cty −1 (1) to the secular precession of the longitude of the perihelion of Mercury arising from the other planets of our solar system a e-mail: lorenzo.iorio@libero.it up to Saturn. A similar scenario, but with the perturbing body moving in an inner orbit with respect to the test particle, was treated in Yamada and Asada [19]. The precession of Eq.
(1) is, partly, a direct consequence of some post-Newtonian accelerations of order O c −2 (1 pN) induced by a distant, pointlike body X; see 1 [a] Cross in Eq. (4) of Will [18]. On the other hand, a mixed, indirect contribution, allegedly of the same order of magnitude of the direct ones, comes also from the interplay between the standard Newtonian third-body 2 acceleration, which, to the quadrupole order, is and the usual 1pN pointlike acceleration due to only the primary's mass in the perturbative calculation by means of the Gauss equations inasmuch the same way as in the case of the Newtonian acceleration due to the quadrupole mass moment of the primary and Eq. (3) [9,17]. In particular, the largest contributioṅ W max 0.16 mas cty −1 (4) to the new precession of Eq. (1) is due to the direct and mixed effects which do not depend on the velocity v X of the distant perturber. Will [18] did not display the direct and indirect contributions to Eq. (4) separately, so that it is not possible to establish the weights of both the effects. Actually, it may have its importance in view of the fact that, as explained below, the mixed effects may be unobservable. The direct acceleration in Eq. (4) of Will [18] which contains v X gives rise to a de Sitter-like precession which is about 0.4 times smaller than Eq. (4) [18]. In Appendix B, we offer our contribution by analytically working out the direct precessions induced by all the accelerations entering [a] Cross in Eq. (4) of Will [18] without making any simplifying assumptions concerning the orbital configuration of both the perturbed test particle and the distant pointlike perturber X. For Mercury, we find a total pN third-body perihelion precession induced by the planets from Venus to Saturn which amounts to 0.15 mas cty −1 , which disagrees with Eq. (1). In particular, the total direct precession due to the first two accelerations entering Eq. (4) of Will [18] amounts to 0.087 mas cty −1 instead of 3 Eq. (4).
In view of the fact that the largest 1pN contribution to the Mercury's perihelion precession oḟ is rescaled in terms of the PPN parameters β, γ , which are equal to 1 in general relativity, Will [18] argues that, since the forthcoming BepiColombo mission is expected to improve out knowledge of β, γ to the 10 −6 level [6,14], then it would be likely possible to measure Eq. (1). Indeed, the resulting theoretical mismodeling in Eq. (5) would be as little as δ˙ GR 0.03 mas cty −1 .
More specifically, Will [18] in the Abstract writes: "At a few parts in 10 −6 of the leading general relativistic precession of 42.98 arcseconds per century, these effects are likely to be detectable by the BepiColombo mission". Furthermore, Will [18] at pag. 191101-4 writes: "If BebiColombo can reach a part per million accuracy in measuring the perihelion advance, […] it will measure, for the first time, relativistic effects on Mercury's orbit arising from the planets that surround it." Conversely, if one is interested in determining the Sun's quadrupole mass moment and angular momentum through their precessions, Eq. (1) would act as a systematic bias on them. Will [18] at pag. 191101-4 writes about his new effects: "[…] their existence and cross-correlations may play a role […] in measurements of the contributions to Mercury's perihelion advance arising from the solar quadrupole moment and frame dragging that will be carried out using data from BepiColombo". In this Communication, we will show that measuring Eq. (1), or our smaller result in Appendix B, is unlikely, mainly because of the uncertainties in the magnitude of the Sun's angular momentum entering the gravitomagnetic apsidal rate of change and in the spatial orientation of the Sun's spin axis affecting especially the precession induced by the solar quadrupole mass moment. As a byproduct, our results will be useful in assessing the impact of the latter source of systematic uncertainty on the possible measurement of the Lense-Thirring effect itself with BepiColombo. Finally, our exact calculation of the direct precessions have a general validity, and can be fruitfully applied in several astronomical and astrophysical scenarios like, e.g., exoplanets or the stellar system orbiting the supermassive black hole in the Galactic Center characterized by arbitrary eccentricities and inclinations.

Our analysis
As a general remark, we note that the indirect, mixed effects, which arise from the simultaneous interplay of at least two accelerations A, B in the calculation of the averaged precessions of the Keplerian orbital elements with the Gauss equations [9,17], are likely undetectable in practical data reductions. Indeed, as far as our case is concerned in which A is, say, Eq. (3) and B is Eq. (2), data analysts of virtually all groups scattered around the world routinely model the Newtonian N -body interactions and the 1PN pointlike acceleration due to the primary to the best of our current knowledge of the parameters entering them which, of course, is necessarily imperfect. Thus, the actual output of data reductions like residuals of, say, ranges, range-rates, etc. would not show the indirect, mixed effects in full. They could only contain negligible signatures, if any, due to the mismodeling in the planetary masses and in the PPN parameters β, γ in terms of which the 1PN point particle acceleration is expressed. Instead, at least in principle, the observables' residuals should fully display the direct effects (unless they have been somewhat removed in the estimation of, say, the initial state vectors) induced by some new accelerations, like those of [a] Cross in Eq. (4) of Will [18] which, perhaps, may still not be included in the dynamical models fit to the observations by some groups. Otherwise, one should not model both Eqs. (2) and (3) at all, and subtract their theoretically computed signals from the resulting huge residuals. It does not seem certainly viable. Even from the point of view of a covariance analysis, while it would be possible, in principle, to explicitly solve for and estimate dedicated scaling parameter(s) accounting for every single acceleration entering the equations of motion, this could not be done for the indirect, mixed effects. In the following analysis, we will treat Eq. (1) as if it were a potentially measurable effect, irrespectively of its origin.
In addition to the well known 1PN pointlike precession of Eq. (5) due to solely the primary's mass, there are other two further effects affecting the perihelion of Mercury which should be regarded as serious sources of potential systematic uncertainties: they are due to the first even zonal harmonic J 2 of the multipolar expansion of the Sun's Newtonian gravitational potential, and the general relativistic gravitomagnetic field of the Sun induced by its angular momentum S . Their precessions depend not only on the size of J 2 , S , but also on the orientation of the Sun's spin axisŜ in space which must enter the error budget as well. Their exact expressions, valid in any coordinate system and for arbitrary orbital configurations, are [7] The Sun's quadrupole mass moment and angular momentum are currently known to the level of accuracy listed in Table 1 along with the nominal values of the precessions of Eqs. (7), (8). It can be noted that, if, on the one hand, it could be hoped that the expected determinations of J 2 by BepiColombo to the 10 −10 level [1,6,14] may be accurate enough to make Eq. (1) at least larger than the mismodelled J 2 -induced precession, on the other hand, a lingering 6% uncertainty in S would imply an a priori theoretical uncertainty in the Lense-Thirring precession of Eq. (8) as large as 0.13 mas cty corresponding to 58% of Eq. (1) and 86% of our result in Table 2.
As announced before, also the current uncertainties in the Carrington elements parameterizingŜ play a crucial role in view of the resulting mismodeling in Eq. (7). Indeed, a standard Root-Sum-Square (RSS) calculation of the error in˙ J 2 due to the uncertainties in i , , treated as two independent variables, yields Furthermore, Figs. 1, and 2 straightforwardly depict Eq. (7) as function of J 2 , i , as independent variables allowed to vary within their ranges of assumed uncertainties [2,6]. Their full range of variation is about twice Eq. (9). Instead, as shown by Table 1 Relevant Sun's physical parameters along with the most recent uncertainties for some of them appeared in the literature, and nominal quadrupolar and Lense-Thirring perihelion precessions for Mercury. As far as S is concerned, the values quoted for its size and uncertainty were obtained by calculating the mean and the standard deviation of the figures quoted in Table 1 the Lense-Thirring precession is not significantly impacted by the uncertainty in the Sun's spin axis orientation. From the point of view of a possible measurement of the Lense-Thirring effect, Eq. (9) corresponds to a 9% uncertainty in the gravitomagnetic precession. Figure 3 shows the impact of the uncertainties in the Carrington elements on the direct BepiColombo observables, i.e. range and range-rate. It can be noticed that the resulting mismodeled signatures amount to 1 − 1.5% of the nominal Lense-Thirring ones. Schettino et al. [13], with dedicated covariance analyses performed with simulated data of BepiColombo, detailed the practical difficulty of satisfactorily separating J 2 from S , and the impact of S itself in estimating of J 2 in various scenarios.

Conclusions
The overall post-Newtonian third-body precession of the longitude of the perihelion of Mercury recently calculated by Will [18] amounts to˙ W 0.22 mas cty −1 ; according to Will [18], it should be measurable by the forthcoming Bepi-Colombo mission. If, on the one hand, a determination of J 2 at the 5 × 10 −10 level, expected from BepiColombo, may reduce the mismodeling in the quadrupolar perihelion pre- Table 2 Doubly averaged 1pN third-body perihelion precessions of Mercury, in mas cty −1 , induced by Venus, Earth, Mars, Jupiter, Saturn via Eqs. (B1)-(B7). The resulting total precession amounts to 0.15 mas cty −1 ; in particular, Eqs. (B1)-(B4) yield a combined overall precession of 0.087 mas cty −1 , contrary to 0.16 mas cty −1 claimed by Will [18]. The discrepancy seems to be due to the indirect, mixed effects  on the other hand, the uncertainties inŜ would yield δ˙ J 2 81%˙ W = 9%˙ LT . Furthermore, the current 6% uncertainty in S would cause a further bias as large as δ˙ LT 58%˙ W . It seems that the indirect contributions to˙ W arising from the mixing of the Newtonian Nbody term with the 1pN pointlike acceleration of the Sun in the perturbative analytical calculation, which may not be measurable, amounts to about 0.07 mas cty −1 . Indeed, our calculation returns 0.15 mas cty −1 for the total direct post-Newtonian perihelion precession of Mercury induced by the other planets from Venus to Saturn, making, thus, even more pessimistic the perspective of measuring it. The simulated Earth-Mercury range and range-rate time series due to the imperfect knowledge ofŜ are about at a per cent level of the nominal Lense-Thirring signatures. Finally, we not that our exact calculation for such kind of general relativistic precessions are valid for any orbital configuration of both the test particle and the third body. Thus, they can be applied also to other astronomical and astrophysical natural laboratories characterized by large eccentricities and inclinations like, e.g., several exoplanetary systems and the stars orbiting the supermassive black hole in Sgr A * in which the coplanarity condition is not fulfilled.  Table 1 J 2 : Dimensionless quadrupole mass moment of the primary S : Magnitude of the angular momentum of the primarŷ S : Unit vector of the spin axis of the primary : Longitude of the ascending node of the Sun's equatorial plane with respect to the Vernal Equinox along the Ecliptic. One of the Carrington elements i : Inclination of the Sun's equatorial plane to the plane of the Ecliptic. One of the Carrington elementŝ Sun's spin axis unit vector in terms of the Carrington elements r : Distance of the test particle from the primarŷ r : Unit vector of the position vector of the test particle v : Velocity vector of the test particle a : Semimajor axis By assuming e X = 0, I = I X , = X [18], the resulting full shift per revolution of the test particle turns out to be (B6) In the case of Mercury, the discrepancy between the full precession and the coplanarity-based approximated one, from which Eq. (B6) was derived, amounts to −0.03 mas cty −1 for X=Venus. The third line of Eq. (4) of Will [18] provides us with the following 1pN "gravitomagnetic" acceleration of order O (G) due to the velocity v X of the third body Its exact, doubly averaged perihelion precession turns out to be d dt = 2μ X n X b c 2 a X 1 − e 2 X cos I X + sin I X tan I 2 cos .
(B8) For = 0, I = I X , e X = 0, Eq. (B8) agrees with the precession which can be inferred from the fourth term of Eq. (1) of Will [18] by taking the ratio of it to the orbital period P b of the perturbed test particle. The numerical discrepancy between Eq. (B8) and the approximated expression by Will is negligible; indeed, in the case of Mercury perturbed by Venus, they differ by just 2 × 10 −5 mas cty −1 yielding botḣ = 0.014 mas cty −1 . The total contribution of all planets from Venus to Saturn to the Mercury's precession of Eq. (B8) amounts to˙ = 0.06 mas cty −1 .
See Table 2 for a detailed overview of the contributions of the planets from Venus to Saturn to the Mercury's direct 1pN third-body perihelion precession. It can be noted that, while the total "gravitomagnetic" effect arising from the third line of Eq. (4) of Will [18] agrees with the results by Will [18] himself, our total precession due to the first two lines of Eq. (4) of Will [18] is about half than that claimed by Will [18]. Such a discrepancy seems to be attributable to the indirect, mixed effects.