{"id":"4c822ff3-75d5-42ba-99df-1ecb8464f1eb","arxiv_id":"1908.03468","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Symplectic correctors improve energy conservation but not the accuracy of secular precession frequencies in 20 Myr Solar System integrations, because the {B,{A,B}} shadow-Hamiltonian term causes a non-oscillatory periastron drift.","lead":"This numerical-methods paper shows that symplectic correctors, which dramatically reduce energy error in planetary N-body simulations, do not improve the accuracy of the secular precession frequencies that drive long-term Solar System evolution. It traces the failure to a specific error term in the integrator's shadow Hamiltonian and shows that alternative integrators which remove this term perform much better.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Reference-solution accuracy is the key assumption; indirect evidence supports it, but a direct high-precision check would settle it.","rationale":"I agree with the reader that the IAS15 reference is the weakest assumption. The paper's Appendix A provides an indirect but sensible argument; the distinct, non-collapsing round-off noise at small timesteps indicates the tested integrators' errors dominate over any IAS15 error. The reproduction of Laskar et al. (2011) frequencies to within 0.5% for g2 and better for other modes shows the model captures the relevant physics, though it is not a direct accuracy bound. The central mechanism—{B,{A,B}} causing secular periastron drift—has independent support from the composition-operator tests, which match the observed one-step and cumulative errors in simpler systems. Therefore, even if the IAS15 reference had a small bias, the qualitative conclusion that WHC does not improve secular frequencies over long timescales is backed by these simpler tests and by the presence of the {B,{A,B}} term in WHC's shadow Hamiltonian. A direct high-precision reference check would remove the residual uncertainty and is cheap enough to run. I do not see grounds to lower the confidence or change the ACCEPT verdict; the concern is noted but does not move the verdict.","tokens_in":18513,"tokens_out":23311,"duration_ms":232468,"concrete_test":"Rerun the 20 Myr Solar System integration with IAS15 at a stricter tolerance (e.g., 1e-13 instead of the default) and recompute the g-mode frequencies; if they shift by more than 1e-12 arcsec/yr, the reference is not converged and the relative WH/WHC errors in Fig. 2 could be distorted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline result—that WHC does not improve secular frequency accuracy over WH for the inner Solar System—is established by comparing each integrator's measured g-mode frequencies to an IAS15 reference. If IAS15 itself had a secular-frequency bias at the level of the differences being measured (e.g., ~1e-7 arcsec/yr in g1 at dt = 8 d), the relative ranking of WH and WHC could be an artifact of the reference rather than an intrinsic property of the integrators. Appendix A argues indirectly: because the error curves for different integrators remain distinct and non-collapsing at small timesteps, IAS15 must be more accurate than the differences between the tested integrators. This is plausible, and the agreement with Laskar et al. (2011) adds confidence, but it is not a direct error bound on IAS15's secular frequencies. The explanatory mechanism—the {B,{A,B}} shadow-Hamiltonian term producing secular periastron drift—is independently verified in one-planet GR and two-planet tests using the composition-operator predictions, so the mechanistic claim is less exposed. The quantitative frequency-error plots in Fig. 2, and hence the precise claim of 'no improvement,' rest on the reference assumption, making it load-bearing though adequately defended.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper performs 20 Myr integrations of the Solar System with REBOUND using the WH, WHC, WHCKL, and SABACL4 symplectic integrators and an IAS15 reference solution. It compares energy error, semi-major axis error, periastron error, and the g-mode secular frequencies. The central finding is that adding high-order symplectic correctors to WH reduces energy error by orders of magnitude but does not improve the accuracy of the inner-planet secular frequencies, and can even slightly degrade g1. Using the BCH expansion of the shadow Hamiltonian, the authors identify the {B,{A,B}} term as responsible: it contributes negligibly to the energy but produces a non-oscillatory, secularly growing periastron error. They develop a composition-operator method to isolate individual shadow-Hamiltonian terms without computing Poisson brackets, verify the mechanism on one-planet GR and two-planet Mercury-Saturn problems, and show that WHCKL and SABACL4, which remove this term, do improve secular-frequency accuracy. The paper concludes that energy conservation alone is an insufficient convergence metric for secularly evolving systems.","tokens_in":18633,"tokens_out":38787,"duration_ms":339849,"significance":"This is a substantive and counterintuitive result for long-term N-body simulation: a three-order-of-magnitude improvement in energy conservation does not translate into improved secular-frequency accuracy for WHC. The paper's main strength is that the empirical convergence study is paired with a parameter-free, falsifiable mechanistic explanation. The composition-operator construction is simple, uses only already implemented evolution operators, and makes quantitative predictions that are independently confirmed on reduced test problems before being used to interpret the full Solar System runs. The manuscript is also reproducible in practice, as it relies on the open-source REBOUND and REBOUNDx packages. If the conclusions hold, the paper establishes an important methodological lesson: the convergence metric must match the physical quantities of interest, and integral-of-motion diagnostics can be misleading. The authors are appropriately cautious about the physical meaning of ultra-precise secular frequencies and about the role of chaos.","major_comments":[],"minor_comments":[{"comment":"Equation (22) is missing the identity subtraction: the correct relation is phi[{A,{A,B}}]_t(y0) = Id(y0) + t^{-2}(C^AAB_t(y0) - Id(y0)) + O(t^2). As printed, the formula contains a spurious O(t^{-2}) term and is inconsistent with Eq. (19). Although the subsequent applications use the composition directly rather than this formula, the typo should be corrected.","section":"Eq. 22"},{"comment":"The error terms in Eqs. (23) and (24) appear to be misstated: since H(y0)-H(C^AAB_{alpha t}(y0)) is of order t^3, the remainder should be O(t^4), not O(t^2). As written, the claimed error term dominates the leading contribution, which would make the approximation useless. Please verify and correct the order notation.","section":"Eqs. 23 and 24"},{"comment":"In Eq. (17), the expansion of exp(t(LA+LB)) contains the term (t^3/6)(LA+LB)^3, not (LA+LB)^6. The exponent 6 appears to be a typographical error; as written it spoils the Taylor expansion.","section":"Eq. 17"},{"comment":"The abstract states that symplectic correctors do not improve the accuracy of secular frequencies compared with WH, but Fig. 2 and Sec. 4.6 show that for the outer Solar System modes g5-g8 the correctors do improve accuracy. The abstract should explicitly say 'for the inner Solar System', as Sec. 5 does, to avoid overgeneralization.","section":"Abstract and Sec. 5"},{"comment":"The argument in Appendix A that IAS15 can be taken as the true solution rules out a dominant random round-off error in the reference solution, but it does not directly exclude a systematic bias in IAS15's secular frequencies at the level of the differences being measured in Fig. 2. The manuscript would be strengthened by stating the IAS15 tolerance settings and by reporting a direct check, e.g., comparing g-mode frequencies obtained with IAS15 at two different error tolerances.","section":"Appendix A"},{"comment":"There are two minor presentation issues: the figure legend uses 'SABA4CL' while the text uses 'SABACL4' for the same method, and the statement in Sec. 2.6 that a standard Fourier transform determines the frequency 'to within one Nyquist frequency' should read 'to within one frequency bin (1/T)'.","section":"Fig. 2 and Sec. 2.6"}],"recommendation":"minor_revision","confidential_remarks":"The central claim is sound and the paper is well within the journal's scope. The only substantive caveat is the IAS15 reference-solution assumption, but the independent mechanism tests and the convergence of many integrators to the same frequencies make this adequately supported for publication. The requested changes are local typographical and clarity fixes. I do not see any novelty-disclosure or citation-pattern concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should read this one. Rein, Brown, and Tamayo show something counterintuitive that holds up: adding symplectic correctors to Wisdom-Holman makes the energy error three orders of magnitude smaller but does not improve the secular frequencies in the inner Solar System over 20 Myr. They also explain why, which is what makes it a real paper rather than a curiosity.\n\nThe new piece is the empirical demonstration plus a composition-operator trick that lets you compute the leading shadow-Hamiltonian error term on any observable without deriving Poisson brackets by hand. The BCH/shadow-Hamiltonian framework is established, but the application to secular periapsis drift is new. The simple tests — one planet with GR, two planets — confirm the predicted secular growth from the {B,{A,B}} term, and WHCKL/SABACL4, which remove that term, do much better. No fitting is involved; the error terms come from the BCH expansion and the simulations test the prediction. That is solid.\n\nThe main soft spot is the assumption that IAS15 is the true solution. It is load-bearing for the quantitative frequency-error plots. The defense in Appendix A is indirect: at small timesteps the error curves for different integrators do not collapse onto each other, which would happen if IAS15 were the noisy one. That is reasonable but not a direct error bound. I do not think it sinks anything; the mechanism is verified independently in the simple tests, and the agreement with Laskar et al. (2011) adds confidence.\n\nMinor things: Eq. 22 appears to drop the identity and the -Id from Eq. 21 (a typo, not a substantive error), and no code or scripts are shipped, though REBOUND is public and the runs are reproducible. The citation pattern is clean; self-citations to REBOUND and prior integrator papers are appropriate.\n\nBottom line: the central claim holds, the diagnostic tool is useful, and the recommendation about choosing validation metrics for long-term integrations is one people should actually follow. Send it out for review; it deserves a serious referee and should get published with minor revisions.","headline":"A solid methods paper: it convincingly shows energy error alone misjudges symplectic integrators for secular dynamics, and the WHC result should change how the field validates long integrations.","tokens_in":19210,"tokens_out":2352,"would_cite":true,"duration_ms":25312,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65P10","70F15"],"pacs":[],"model":"deepseek-v4-flash","headline":"Symplectic correctors cut energy error 1000x but not secular frequencies","keywords":["symplectic integrators","Wisdom-Holman method","secular frequencies","shadow Hamiltonian","periastron precession","energy conservation","convergence study","Solar System N-body simulations"],"falsifier":"In the same 20-Myr inner Solar System setup, replace IAS15 with an independent high-accuracy reference, such as an even tighter-tolerance IAS15 or a different high-order integrator, and compare WH versus WHC at timesteps around 8 to 20 days: if WHC's g1 or g2 error becomes clearly smaller than WH's once the reference changes, the central paradox disappears.","tokens_in":18242,"feed_emoji":"🪐","tokens_out":6370,"duration_ms":57924,"temperature":0.7,"pith_summary":"This paper asks which accuracy metric actually matters for long integrations of planetary systems, and answers that energy conservation can be deeply misleading. In 20-million-year Solar System integrations, the Wisdom-Holman integrator with symplectic correctors reduces the energy error by three orders of magnitude compared with standard Wisdom-Holman, yet reproduces the inner planets' secular precession frequencies no better, and sometimes slightly worse. The cause is a single term in the integrator's shadow Hamiltonian, {B,{A,B}}, whose errors do not oscillate away but accumulate secularly as an artificial periastron precession. Methods that remove that term, the modified-kernel Wisdom-Holman method and the SABACL4 family, do improve secular frequencies substantially. A reader should care because stability and long-term evolution studies that validate simulations by energy conservation alone may be overestimating their accuracy.","feed_headline":"Correctors cut energy error 1000x but not secular frequencies","feed_subtitle":"A 20-Myr Solar System test traces the paradox to one shadow-Hamiltonian term driving artificial periastron precession.","key_machinery":"The machinery is the shadow Hamiltonian of the split integrator plus a composition-operator method that evaluates the effect of an individual Poisson-bracket term on any observable without computing derivatives. For WH, the shadow Hamiltonian contains $t^{2}$ {A,{A,B}} and $t^{2}$ {B,{A,B}} at order $t^{3}$. The paper constructs operators C^AAB_t and C^BAB_t from the already-available Kepler and kick evolutions, which reproduce to leading order the evolution generated by {A,{A,B}} and {B,{A,B}}. These operators let the authors compute periastron and phase errors due to each term, show that the {B,{A,B}} errors accumulate linearly, and predict the crossover timescale without running long integrations.","core_discovery":"The central discovery is that the error of a symplectic integrator in a secularly evolving planetary system is controlled not by the size of the shadow-Hamiltonian terms but by whether a given term produces oscillatory or secular errors. For the Wisdom-Holman split H=A+B, the leading correction {A,{A,B}} (order ϵ $dt^{2}$) oscillates and averages out in the periastron, so removing it with correctors does not help long-term secular accuracy. The smaller term {B,{A,B}} (order $ϵ^{2}$ $dt^{2}$) generates a non-oscillatory, linearly growing periastron error that dominates after roughly $10^{5}$ to $10^{6}$ orbits; once it dominates, standard WH and corrected WHC have the same secular-frequency error despite a $10^{3}$ difference in energy error. Integrators that explicitly remove {B,{A,B}}, namely WHCKL and SABACL4, restore the expected gain in secular accuracy.","pith_inferences":["The same composition-operator technique could be applied to other observables, such as resonant angles, tidal or general-relativistic precession rates, to identify non-oscillatory error terms that energy checks miss.","For exoplanet stability surveys that use WH-family integrators, this suggests a validation protocol: check convergence of the system's secular frequencies, not just energy, because artificial periastron precession could bias resonance boundaries.","A testable extension would vary the planet mass ratio epsilon and orbital period in a two-planet system to map where the {B,{A,B}} crossover time falls, giving a practical rule for when correctors are worthwhile.","The framework implies that labels like 'higher-order' or 'corrector-enhanced' do not guarantee long-term accuracy unless the specific non-oscillatory term in the shadow Hamiltonian is removed."],"forward_implications":["Energy conservation alone is not a sufficient convergence test for secularly evolving planetary systems; round-off-level energy error does not mean the trajectory has converged.","For inner Solar System secular frequencies over about 20 Myr, WHC buys nothing over WH near their recommended timestep, despite a 10^3 smaller energy error.","WHCKL and SABACL4, which remove the {B,{A,B}} term, deliver substantially better periastron and secular-frequency accuracy.","With such methods at a 10-day timestep, about eight steps per Mercury period, Mercury's secular frequency g1 can be determined to 10^-9 arc-seconds per year in this model.","Over very long timescales the error budget is a sum of oscillatory discretization error, secularly growing discretization error, round-off error, and chaotic divergence; the composition operators identify which term dominates."],"supporting_citations":[{"why":"Introduces the WH operator splitting and the Kepler-plus-perturbation Hamiltonian split on which the whole analysis is built.","marker":"Wisdom & Holman (1991)"},{"why":"Introduces symplectic correctors (WHC) and the modified-kernel method (WHCKL); removing {A,{A,B}} versus {B,{A,B}} is the central comparison.","marker":"Wisdom et al. (1996)"},{"why":"Introduces the SABA integrator family whose correctors remove the {B,{A,B}} term, one of the two better-behaved methods tested.","marker":"Laskar & Robutel (2001)"},{"why":"Provides the implementations of WHC, WHCKL, and SABACL4 used in the convergence study.","marker":"Rein et al. (2019b)"},{"why":"Supplies the IAS15 integrator used as the reference true solution of the model.","marker":"Rein & Spiegel (2015)"},{"why":"Provides the published Solar System secular frequencies used as a consistency check for the 20-Myr integrations.","marker":"Laskar et al. (2011)"},{"why":"Provides the proper-mode transformation and modified Fourier algorithm used to extract the secular frequencies.","marker":"Laskar (1990)"},{"why":"Supplies the 1/r^3 general-relativity potential included in the model.","marker":"Nobili & Roxburgh (1986)"}],"fun_headline_variants":["Energy error 1000x smaller, but secular frequencies unchanged: why?","Shadow term causes artificial periastron precession in symplectic sims","The real error culprit in long-term N-body: not energy but a small term","Correctors don't fix secular drift: it's a specific shadow term","One subtle error term dominates secular inaccuracy in symplectic methods"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The reference 'true solution' is the IAS15 integration; if IAS15 itself had appreciable secular-frequency error, the measured ordering of the methods could be distorted.","fun_headline_variants_meta":{"raw":{"variants":["Energy error 1000x smaller, but secular frequencies unchanged: why?","Shadow term causes artificial periastron precession in symplectic sims","The real error culprit in long-term N-body: not energy but a small term","Correctors don't fix secular drift: it's a specific shadow term","One subtle error term dominates secular inaccuracy in symplectic methods"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00095,"raw_usage":{"total_tokens":4074,"prompt_tokens":988,"completion_tokens":3086,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":604,"completion_tokens_details":{"reasoning_tokens":2987}},"tokens_in":604,"tokens_out":3086,"duration_ms":22922,"temperature":1.0,"reasoning_tokens":2987,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:12:48.653192+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In the same 20-Myr inner Solar System setup, replace IAS15 with an independent high-accuracy reference, such as an even tighter-tolerance IAS15 or a different high-order integrator, and compare WH versus WHC at timesteps around 8 to 20 days: if WHC's g1 or g2 error becomes clearly smaller than WH's once the reference changes, the central paradox disappears.","supporting_citations":[{"cited_title":"1996, Fields Insti- tute Communications, Vol","cited_arxiv_id":null,"evidence_quote":"Introduces symplectic correctors (WHC) and the modified-kernel method (WHCKL); removing {A,{A,B}} versus {B,{A,B}} is the central comparison."},{"cited_title":"& Robutel, P","cited_arxiv_id":null,"evidence_quote":"Introduces the SABA integrator family whose correctors remove the {B,{A,B}} term, one of the two better-behaved methods tested."},{"cited_title":"2011, A&A, 532, A89","cited_arxiv_id":null,"evidence_quote":"Provides the published Solar System secular frequencies used as a consistency check for the 20-Myr integrations."},{"cited_title":"& Roxburgh, I","cited_arxiv_id":null,"evidence_quote":"Supplies the 1/r^3 general-relativity potential included in the model."}],"review_version":1}