{"id":"5008d3ce-5e41-4b0c-b2d6-b6845e59e8d1","arxiv_id":"1908.09670","paper_version":3,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The 1pN N-body perihelion precession of Mercury is 30.4 microarcseconds per century, five times smaller than an earlier published value.","lead":"A new calculation of tiny general-relativistic effects on Mercury's orbit shows they are even smaller than previously thought. The result makes it harder to measure these effects with the BepiColombo mission.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central value rests on Eqs. (1)–(3), imported from Will (2018) without verifying they are the complete 1PN N-body acceleration; if a same-order term is missing, all quoted rates shift.","rationale":"The reader's weakest assumption—validity and completeness of Eqs. (1)–(3)—is exactly the load-bearing point. The paper provides strong internal evidence: the numerical and analytical routes agree, and the correction to Iorio (2018) is traced to a specific equation. However, the physical input is an imported set of accelerations whose derivation is not reproduced. The paper explicitly identifies them as a 'particular case' of the full EIH equations, but it does not demonstrate that the neglected EIH terms cancel in the Sun-Mercury-other-planets system or are numerically small at the claimed precision. Since both computational methods share the same input, agreement between them cannot establish completeness. The central quantitative claim therefore stands only conditionally on the completeness of the Will equations. The proposed check—an independent reduction of the full EIH equations to the same hierarchy—would settle the issue directly. This is not an accusation of error; it is a standard verification step for a result that corrects a previously published value by a factor of five.","tokens_in":9099,"tokens_out":22416,"duration_ms":224307,"concrete_test":"Re-derive the 1PN acceleration of a test particle near a static primary due to a distant mass X from Poisson & Will (2014), Eq. (9.127), in the limit m -> 0, subtracting the acceleration of the primary due to X and keeping all terms linear in M_X and of order c^{-2} at leading order in r/r_X. Compare term-by-term with Eqs. (1)–(3). If extra terms survive, add them to the Mercury integrations of Sect. 2.1, re-fit the 1-century slopes, and check whether \\dot\\varpi shifts by more than the EPM2017 formal error (8 μas/cty).","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's headline perihelion rate (30.4 μas/cty) and the companion rates for I, Ω, and ε all follow from the three accelerations in Eqs. (1)–(3), which are quoted from Will (2018) as 'a particular case' of the full Einstein-Infeld-Hoffmann equations (Poisson & Will 2014, Eq. 9.127). Both the numerical integration (Sect. 2.1) and the doubly averaged analytical calculation (Sect. 2.2) use exactly these three terms, so their mutual agreement tests only internal consistency, not completeness. The full EIH equations contain additional same-order terms in the test-particle limit, including pieces proportional to v_X^2, (n_X · v_X)^2, and the acceleration of X; after subtracting the primary's acceleration these may not cancel. The paper does not show that such terms vanish or are numerically negligible for the actual configuration, where Venus and Earth have r/r_X ~ 0.5, so the 'distant body' hierarchy is only marginally satisfied. If any omitted term contributes at the level of a few μas/cty or more, the corrected perihelion value and the comparison with EPM2017 errors would change, although the factor-of-five reduction relative to Iorio (2018) would likely survive.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes the first-order post-Newtonian (1pN) orbital effects on Mercury's Keplerian elements that arise from the accelerations (1)–(3), adopted from Will (2018) as a 'particular case' of the full Einstein-Infeld-Hoffmann equations, with the perturbing planets taken to be Venus through Saturn. The author finds secular rates for the inclination, node, perihelion, and mean longitude at epoch (−4.3, 18.2, 30.4, and 271.4 μas cty−1, respectively), no secular rates for the semimajor axis and eccentricity, and a perihelion precession about five times smaller than the value reported in Iorio (2018), which is traced to an error in Eq. (B5) of that paper. The results are obtained both by numerical integration of the barycentric equations and by doubly averaged Gauss equations, and are compared with formal uncertainties derived from EPM2017 ephemerides, leading to a pessimistic conclusion about measurability with BepiColombo.","tokens_in":9338,"tokens_out":28886,"duration_ms":282746,"significance":"If the adopted acceleration model is accepted as complete, the paper is a valuable self-correction of a published value and provides the first quantitative estimates for all of Mercury's orbital elements under the Will (2018) accelerations. The numerical and analytical methods cross-validate each other, there are no fitted parameters, and the comparison with published ephemeris uncertainties is concrete. The paper also explicitly identifies the origin of the previous error, which is a strength for the literature. The main risk is that the three accelerations may not represent the complete 1pn N-body acceleration in the adopted test-particle limit, and the numerical slopes lack formal uncertainties; both issues affect the central quantitative claims.","major_comments":[{"comment":"The paper adopts the three accelerations as 'a particular case' of the full 1PN Einstein-Infeld-Hoffmann equations, but it does not demonstrate that these three terms are the complete 1PN N-body acceleration in the Sun-centered test-particle limit. The full EIH equations contain additional contributions proportional to v_X^2, (r_X_hat · v_X)^2, and the acceleration of X; after subtracting the primary's acceleration, these yield terms of the same order as Eq. (1) for the actual configuration, e.g., r/r_X ≈ 0.5 for Venus and v_X^2/c^2 ≈ (GM_S/r_X)/c^2. The numerical integration (Section 2.1) and the analytical calculation (Section 2.2) both use exactly Eqs. (1)–(3), so their agreement is a check of internal consistency, not of completeness. Since the abstract and Section 3 present the results as the 1pN N-body precessions and compare them with EPM2017 errors, the quoted rates (30.4 μas cty−1 for ϖ, and the companion rates) are not yet established as the full 1pN N-body effect. Please either include the omitted terms or state explicitly that the results refer only to the subset in Eqs. (1)–(3) and justify numerically that the omitted terms are negligible at the μas cty−1 level.","section":"§1, Eqs. (1)–(3)"},{"comment":"The secular rates are obtained by fitting a linear trend to difference time series, but the fitted slopes are reported without formal uncertainties or residual statistics. This matters because the rates are at the same level as the EPM2017 formal errors used in Section 3, and because the analytical formulas are only leading order in e (Mercury's e ≈ 0.21). Please report, for each element, the 1σ slope uncertainty from the fit, the rms residual of the difference time series, and the numerical-versus-analytical difference; this would make the 'excellent agreement' claim in Section 2.2 quantitative and support the three-significant-figure rates in the abstract.","section":"§2.1 and Fig. 1"}],"minor_comments":[{"comment":"The description of the numerical integration is not sufficient for reproducibility; please specify the integrator, time step, tolerance, and the procedure used to compute the difference time series, and state whether the same HORIZONS initial conditions are used in both runs.","section":"§2.1"},{"comment":"The manuscript contains typographical artifacts from the LaTeX source ('orbit al' in the running title, 'pontlike' for 'pointlike', and several 'doteq' symbols); please clean these up in the published version.","section":"Title and text"},{"comment":"The phrase 'at the same level of, or larger by one order of magnitude than' is too vague; for each element, the predicted rate should be compared numerically with the quoted EPM2017-derived uncertainty, e.g., 30.4 vs 8 μas cty−1 for ϖ and 18.2 vs 24 μas cty−1 for Ω.","section":"Abstract and §3"},{"comment":"The estimate δn_b^obs = 20 mas cty−1 uses δμ_obs = 1e10 m^3 s−2 and δa_obs; please make the propagation formula explicit so that the reader can reproduce this factor.","section":"§3"},{"comment":"The proposed linear combinations of supplementary advances for several planets are only sketched; a concrete example, even for a simplified two-planet scenario, would help establish that the 1pN N-body effects can be separated from the J2 and Lense-Thirring precessions.","section":"§4"}],"recommendation":"major_revision","confidential_remarks":"The completeness issue in Eqs. (1)–(3) is likely to be the main point of contention. If the authors can show that the omitted EIH terms are negligible at the required accuracy, or that the paper's scope is explicitly restricted to Will's three terms, the paper may become acceptable. The self-correction of Iorio (2018) is transparent and a strength of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper is a straightforward, honest self-correction, not a big discovery. Iorio finds that his own earlier value for the 1pN N-body perihelion precession of Mercury (150 μas/cty) was five times too large; the corrected figure is 30 μas/cty. He also provides the first secular rates for the inclination, node, and mean longitude at epoch due to these accelerations, and compares them with EPM2017 formal errors.\n\nWhat is actually new: the correction itself, traced to a specific error in Eq. (B5) of Iorio (2018); new analytical formulas for a, e, I, Ω, and ε; and numerical integrations that match the analytic results. That agreement is real evidence the formulas are right, though both methods use the same starting accelerations. The comparison with ephemeris uncertainties is careful and ends with a measured negative conclusion: the effects are at or below the level of the formal errors, so viewing them as detectable with BepiColombo is not justified.\n\nSoft spots: the calculation imports Eqs. (1)–(3) from Will (2018) without an independent check that they are the complete 1PN three-body acceleration. The stress-test note about possible missing v_X^2 terms did not land as a concrete flaw on my reading—Will’s equations are standard in this context—but the paper would have been stronger with a sentence confirming completeness. More minor: the fitted numerical slopes have no error bars; the analytic formulas are leading order in e; and the stated “excellent” agreement between numerics and analytics is not quantified. None of these undercut the main result.\n\nThis paper is for people working on solar-system tests of GR and on BepiColombo data analysis. It deserves a serious referee, mainly because it corrects a value already in the literature and prevents an overclaim from propagating further. I’d send it out.","headline":"A small but solid self-correction: the 1pN N-body Mercury perihelion precession is 30 μas/cty, not 150, with new rates for other elements, all likely too small to measure soon.","tokens_in":9865,"tokens_out":3731,"would_cite":true,"duration_ms":39688,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.10.Ce","04.80.Cc"],"model":"deepseek-v4-flash","headline":"General relativity's N-body contribution to Mercury's perihelion precession is 30.4 microarcseconds per century, five times smaller than the previously reported value.","keywords":["general relativity","post-Newtonian approximation","Mercury","perihelion precession","N-body dynamics","Keplerian orbital elements","BepiColombo","planetary ephemerides"],"falsifier":"Numerically integrate Mercury's orbit with the full 1pN N-body equations of motion of the Solar System, rather than the restricted three-body accelerations of Eqs. (1)-(3), with and without the direct 1pN third-body terms, and fit the difference in the perihelion longitude over one century; a slope significantly different from $30.4\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ would invalidate the paper's correction.","tokens_in":8887,"feed_emoji":"🪐","tokens_out":15138,"duration_ms":122964,"temperature":0.7,"pith_summary":"This paper computes, analytically and by numerical integration, the first-order post-Newtonian (1pN) effects that the other planets exert on Mercury's orbit through the restricted N-body accelerations of general relativity. Its central result is that Mercury's 1pN N-body perihelion precession is $30.4\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$, about five times smaller than the value reported in the author's own 2018 paper and much smaller than the $220\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ estimate proposed by Will (2018). The same calculation yields secular rates of $-4.3\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ for the inclination, $18.2\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ for the node, and $271.4\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ for the mean longitude at epoch. Compared with the formal uncertainties from the EPM2017 ephemerides, these signals sit at or below the currently achievable errors, so a clean detection by BepiColombo would require both improved orbit knowledge and a way to handle the competing Sun quadrupole and Lense-Thirring precessions. The paper matters because it corrects a published number that had been proposed as a target for a forthcoming space mission.","feed_headline":"Mercury's relativity shift is just 30 microarcseconds","feed_subtitle":"Five times smaller than a 2018 estimate, it sits at the edge of current measurement errors.","key_machinery":"The central machinery is a set of three 1pN $N$-body accelerations, $A_{G2}$, $A_G$, and $A_{v_X}$ (Eqs. 1-3, from Will 2018), describing how a distant third body $X$ perturbs a test particle at first post-Newtonian order, applied through the Gauss perturbative equations that turn a perturbing acceleration into time derivatives of the Keplerian orbital elements $a,e,I,\\Omega,\\varpi,\\epsilon$. The paper doubly averages these equations over the orbital periods of Mercury and of each perturbing planet to produce leading-order analytic rates, and verifies them with a one-century numerical integration. The load-bearing step is the corrected analytic formula for the perihelion rate produced by $A_G$: the correction to Eq. (B5) of Iorio (2018) removes the factor that had produced the earlier $150\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ value.","core_discovery":"The core discovery is that the 1pN $N$-body perihelion precession of Mercury is $\\dot\\varpi_{1\\mathrm{pN}}^\\mathrm{X}=30.4\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$, not the $150\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ reported in Iorio (2018), because a specific analytic expression for the velocity-dependent acceleration $A_G$ in that paper was wrong. For the same accelerations, the inclination, node, and mean-longitude-at-epoch rates are $\\dot I=-4.3\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$, $\\dot\\Omega=18.2\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$, and $\\dot\\epsilon=271.4\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$. The corrected perihelion value is even smaller than the approximate $220\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ of Will (2018), which strengthens the earlier conclusion that the effect is very hard to measure; the other rates are at or below the formal EPM2017 errors, and the largest signal, in the mean longitude at epoch, is masked by the much larger uncertainty in Mercury's mean motion arising from the solar gravitational parameter.","pith_inferences":["The same doubly averaged machinery could be carried to the Earth-Moon system or to artificial satellites around the Earth, where range measurements are far more accurate than Mercury's astrometry, potentially turning this 1pN third-body acceleration into a measurable signal.","Because the analytic formulas express each rate through the same orbital elements of the perturber, ratios such as $\\dot\\Omega/\\dot\\varpi$ or $\\dot I/\\dot\\varpi$ for a given planet are fixed predictions; a future measurement of one rate could be checked against the others without an absolute calibration.","A direct numerical test of the correction would integrate the full 1pN N-body equations of motion, rather than the restricted hierarchical approximation, and compare the extracted perihelion slope with $30.4\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$; disagreement would reopen the question.","If BepiColombo does eventually deliver a perihelion residual after subtracting the standard Schwarzschild, $J_2$, and Lense-Thirring effects, that residual should be compared with $30\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ rather than the earlier $150$ or $220\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ figures."],"forward_implications":["Mercury's 1pN N-body perihelion precession is $30.4\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$, roughly four times the formal EPM2017 error of $8\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ but likely below realistic systematic errors, so direct detection with BepiColombo is doubtful.","The inclination and node rates, $-4.3$ and $18.2\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$, are comparable to their formal uncertainties ($3$ and $24\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$), which makes them difficult to isolate as clean signals.","The largest computed effect, $\\dot\\epsilon=271.4\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$, cannot currently be measured because the uncertainty in Mercury's mean motion from the solar gravitational parameter is about $20\\,\\mathrm{mas}\\,\\mathrm{cty}^{-1}$, roughly seventy times larger.","If Mercury's orbit is later determined more accurately, linear combinations of the supplementary advances of several planets' elements, following Shapiro (1990), could cancel the competing Sun $J_2$ and Lense-Thirring precessions and extract the 1pN N-body rates.","The analytic formulas apply to any test particle around a static, spherically symmetric primary with a distant massive third body, so the same rates can be obtained for other planets or artificial satellites without new integrations."],"supporting_citations":[{"why":"Supplies the three 1pN N-body accelerations (Eqs. 1-3) and the earlier 220 μas cty^-1 perihelion estimate that this paper corrects and extends.","marker":"Will (2018)"},{"why":"The author's earlier analytical calculation whose perihelion rate of 150 μas cty^-1 is shown here to be wrong, with the error located in Eq. (B5).","marker":"Iorio (2018)"},{"why":"Provides the formal uncertainties in Mercury's secular orbital-element rates used for the measurability comparison.","marker":"Iorio (2019)"},{"why":"Releases the EPM2017 ephemerides formal errors in Mercury's orbital elements and mean longitude.","marker":"Pitjeva & Pitjev (2018)"},{"why":"Gives the full 1pN N-body equations of motion of which Eqs. (1)-(3) are a particular hierarchical case.","marker":"Poisson & Will (2014)"},{"why":"Introduces the linear-combination approach that the paper proposes to use for disentangling the 1pN N-body rates from competing precessions.","marker":"Shapiro (1990)"}],"fun_headline_variants":["Mercury's perihelion shift corrected: 30 μas/century, not 150","Einstein's Mercury shift is 30 μas/century, not 150","Mercury's perihelion shift is 30 μas/cty, 5× smaller than 2018","Mercury's perihelion drift: 30, not 150 μas/cty"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the three post-Newtonian force formulas taken from Will (2018) correctly and completely describe, to first order, how a distant planet's gravity acts on a test particle; if those formulas are wrong or incomplete, every computed rate, including the perihelion correction, would change.","fun_headline_variants_meta":{"raw":{"variants":["Mercury's perihelion shift corrected: 30 μas/century, not 150","Einstein's Mercury shift is 30 μas/century, not 150","Mercury's perihelion shift is 30 μas/cty, 5× smaller than 2018","Mercury's perihelion drift: 30, not 150 μas/cty"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001561,"raw_usage":{"total_tokens":6370,"prompt_tokens":1211,"completion_tokens":5159,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":827,"completion_tokens_details":{"reasoning_tokens":5056}},"tokens_in":827,"tokens_out":5159,"duration_ms":33262,"temperature":1.0,"reasoning_tokens":5056,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:04:03.004694+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate Mercury's orbit with the full 1pN N-body equations of motion of the Solar System, rather than the restricted three-body accelerations of Eqs. (1)-(3), with and without the direct 1pN third-body terms, and fit the difference in the perihelion longitude over one century; a slope significantly different from $30.4\\,\\mu\\mathrm{as}\\,\\mathrm{cty}^{-1}$ would invalidate the paper's correction.","supporting_citations":[{"cited_title":"M., 2018, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the three 1pN N-body accelerations (Eqs. 1-3) and the earlier 220 μas cty^-1 perihelion estimate that this paper corrects and extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The author's earlier analytical calculation whose perihelion rate of 150 μas cty^-1 is shown here to be wrong, with the error located in Eq. (B5)."},{"cited_title":"J., 157, 220","cited_arxiv_id":null,"evidence_quote":"Provides the formal uncertainties in Mercury's secular orbital-element rates used for the measurability comparison."},{"cited_title":"V., Pitjev N","cited_arxiv_id":null,"evidence_quote":"Releases the EPM2017 ephemerides formal errors in Mercury's orbital elements and mean longitude."},{"cited_title":"M., 2014, Gravity","cited_arxiv_id":null,"evidence_quote":"Gives the full 1pN N-body equations of motion of which Eqs. (1)-(3) are a particular hierarchical case."},{"cited_title":"I., 1990, in General Relativity and Gravitation, 1989, Ashby N., Bartlett D","cited_arxiv_id":null,"evidence_quote":"Introduces the linear-combination approach that the paper proposes to use for disentangling the 1pN N-body rates from competing precessions."}],"review_version":1}