{"id":"75946bc1-4d05-4c4d-8268-24e28b205ea1","arxiv_id":"1908.01048","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper introduces a modified Laplace-Lagrange formula, Eq. (10), that adds the Heppenheimer eccentricity factor to handle eccentric companion stars, and extends it to multiple perturbers via superposition.","lead":"This paper compares five analytic formulas for how fast a planet's orbit precesses when its host star has a companion star, then proposes a sixth formula that stitches two classic ones together. The new formula is simple enough for quick estimates of where these precession effects matter in binary-planet systems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The LLM's eccentricity factor is ad hoc and never tested in isolation: Fig. 2 omits LLM, and Fig. 3 confounds it with the multi-perturber sum.","rationale":"Agreement with the reader's weakest_assumption: the ad hoc factorization is indeed the key unproven step. My stress-test reinforces rather than changes the reader's CONDITIONAL verdict. The paper's contribution is clearly framed as an interpolation formula with known limits, and the multi-perturber superposition (Eq. 11) is a standard LL property; those parts are sound. However, the novel element—the eccentricity factor—is asserted, not derived, and the presented numerical evidence does not isolate it. In Fig. 2, the LLM is absent from the very comparison that would show its e_B dependence. In Fig. 3, the LLM is the only model that includes the giant-planet contribution via Eq. 11, so its success is confounded with model completeness. Thus the paper's central assertion of 'quantitatively correct estimates' currently rests on an untested ansatz. The condition I would attach is exactly the missing like-for-like ER3BP test. Since the reader already requested this, no verdict change is needed.","tokens_in":11036,"tokens_out":6308,"duration_ms":60485,"concrete_test":"Perform a single like-for-like numerical test in the ER3BP with one eccentric binary perturber and no giant planet: for α ∈ {0.05, 0.1} and e_B ∈ {0, 0.2, 0.4, 0.6}, extract the secular frequency from direct numerical integration of the full 3BP (as in Fig. 2) and compare LLM (Eq. 10) against LL (Eq. 6) and HEP (Eq. 1). If LLM's mean absolute relative error does not improve on LL's at large e_B, or if the error grows with α, the claim that the ad hoc factor is quantitatively correct fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the modified Laplace–Lagrange model (Eq. 10), g_LLM = 1/4 µ α^2 n b^(1)_{3/2}(α)(1-e_B^2)^(-3/2), gives 'quantitatively correct estimates' of a test particle's secular frequency. The only new ingredient relative to LL is the ad hoc Heppenheimer factor (1-e_B^2)^(-3/2); Section 4 explicitly says it is introduced 'ad hoc' and is not derived. This factor is the load-bearing assumption. If the true dependence of the secular frequency on e_B does not factor out of the Laplace coefficient's α-dependence, then Eq. 10 (and its multi-perturber sum, Eq. 11) will misestimate frequencies at large e_B and moderate α. The paper provides no clean test of this factor. Figure 2, which scans e_B in the ER3BP for α=0.05 and 0.1, does not include the LLM curve. Figure 3, which does show LLM matching a numerical R4BP reference, is not like-for-like: LLM is evaluated in the R4BP with the sum over both the giant planet and the secondary, while the competing models (HEP, AND, GIU, GEO) are evaluated in the ER3BP with only the secondary. The agreement in Fig. 3 may therefore be due to the superposition in Eq. 11 rather than to the new eccentricity factor. Higher-order secular theories (e.g., AND, GEO, Georgakarakos 2002-2009) indicate that e_B corrections do not simply multiply the circular-orbit coefficient at all orders in α, so the ad hoc ansatz has a real correctness risk. No derivation, code, or error bars are provided to support quantitative accuracy.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reviews several analytical models for the secular precession frequency of a circumstellar planet in a binary star system (Heppenheimer, Giuppone, Andrade-Ines & Eggl, Georgakarakos, and the Laplace-Lagrange theory), applies them to synthetic and observed systems, and introduces a 'modified Laplace-Lagrange model' (LLM). The LLM is defined in Eq. (10) as the standard Laplace-Lagrange frequency multiplied by the Heppenheimer eccentricity factor (1 - e_B^2)^(-3/2), and is generalized to multiple perturbers via linear superposition in Eq. (11). The central claim is that this ad hoc modification gives 'quantitatively correct estimates' of the secular frequency, supported by a comparison with numerical integrations of the restricted four-body problem for the HD 41004 system (Figure 3).","tokens_in":11293,"tokens_out":5334,"duration_ms":53042,"significance":"If the central claim holds, the paper provides a simple closed-form extension of the Laplace-Lagrange theory that accounts for eccentric perturbers and can be applied to multiple perturbers, which would be useful for parameter surveys and secular resonance studies. The paper's strengths are that no parameters are fitted to the reference data, the proposed model has the correct limiting behavior (reducing to LL for e_B=0 and to HEP for small alpha), and the numerical reference is obtained from independent direct integrations. However, the validation is incomplete: the new eccentricity factor is never tested in isolation, and the quantitative claim is not backed by a defined error metric. The contribution is modest but potentially useful if the identified gaps are addressed.","major_comments":[{"comment":"The eccentricity factor (1 - e_B^2)^(-3/2) is introduced 'ad hoc' (the paper's own wording) and is not derived from the secular perturbation equations. The two limiting cases (e_B -> 0 and alpha -> 0) do not establish that the factorized form is valid at intermediate alpha and large e_B. In fact, the higher-order model in Eq. (5) (Georgakarakos) shows e_B-dependent corrections that do not simply factor out of the alpha-dependence, so the factorization in Eq. (10) is a genuinely load-bearing assumption. Since the central claim of 'quantitatively correct estimates' depends on this factorization, the paper should either provide a derivation or a clean numerical test that isolates the factor.","section":"Section 4, Eq. (10)"},{"comment":"The comparison in Figure 3 does not isolate the new ingredient of the LLM. The LLM curve is computed in the R4BP with the sum over both perturbers (Eq. 11), while all competing models are computed in the ER3BP with only the secondary star. The unmodified Laplace-Lagrange model (Eq. 6) evaluated in the same R4BP setup (i.e., summed over both perturbers but without the eccentricity factor) is not shown. Without this benchmark, the agreement of the LLM with the reference could be attributable to the linear superposition property of the LL theory rather than to the ad hoc factor. The authors should add the conventional LL model (summed over both perturbers) to Figure 3.","section":"Section 4, Figure 3"},{"comment":"Figure 2 is the only experiment that directly varies the perturber's eccentricity, and it shows that the conventional LL model fails for e_B > 0.1. However, the LLM is not included in this figure. To support the claim that Eq. (10) improves the LL model, the authors should show the LLM curve in the ER3BP comparison (Figure 2), where the only difference from LL is the eccentricity factor. This would provide a direct, controlled test of the load-bearing assumption.","section":"Section 3, Figure 2"},{"comment":"The statement that the LLM gives 'quantitatively correct estimates' is not supported by a quantitative measure. Figure 3 shows a visible systematic overestimate for larger alpha, and no error bars, RMS differences, or relative error values are provided. The authors should present a quantitative comparison (e.g., relative error as a function of alpha) for the LLM against the numerical reference, and ideally also against the unmodified LL model in the same setup, to substantiate the central claim.","section":"Section 4, final paragraph and Figure 3"}],"minor_comments":[{"comment":"The symbol alpha is defined in Section 2.1.4 as a_P/a_B, but in Figure 3 and its description it is reused for a/a_P (the ratio of the test particle's semi-major axis to that of the giant planet). This notational conflict is confusing and should be resolved, for example by using a different symbol for the test-particle ratio.","section":"Section 2.1.4 and Figure 3"},{"comment":"The abstract claims 'better frequency estimates than the traditional model for large eccentricities of the perturber,' but the only application of the new model (Figure 3) uses e_B = 0.2, which is not 'large' in the context of the paper (Figure 2 covers up to e_B = 0.6). Either test the LLM at larger e_B or soften the claim.","section":"Abstract and Section 4"},{"comment":"The notation for the Laplace coefficient b^{(1)}_{3/2} is introduced without explicitly stating that the superscript (1) denotes the order k in Eq. (7). This should be stated for clarity, especially for readers not familiar with Laplace coefficients.","section":"Section 2.2.5, Eq. (7)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a proceedings contribution, and its ambitions are modest. However, as a journal submission, the central claim is not yet adequately supported because the new model's only new ingredient is not tested in a controlled way. The requested additions (LLM in Fig. 2, LL in Fig. 3, and quantitative errors) are straightforward and should be feasible within the manuscript's scope, so major revision rather than rejection seems appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know about this paper: it's a proceedings contribution that reviews five secular-frequency models for circumstellar planets in binaries and then proposes a modified Laplace-Lagrange model (LLM, Eq. 10) by multiplying the LL frequency by the Heppenheimer eccentricity factor (1-e_B^2)^(-3/2). A straightforward linear-superposition extension (Eq. 11) makes it applicable to multiple perturbers. The formula reduces correctly to LL at e_B=0 and to HEP at small alpha. That part is legitimate and potentially useful for parameter surveys.\n\nThe main problem is the validation. The authors never show LLM in the controlled e_B scan of Figure 2, where all other models are compared. So we never learn whether the eccentricity factor actually improves LL for a single external perturber. In Figure 3, the comparison is not like-for-like: LLM is computed in the restricted four-body problem (including the giant planet's contribution), while the competitors are computed in the ER3BP (which only includes the secondary star). The good agreement of LLM there could be entirely due to the summed giant-planet term, not the new factor. The text compounds this with a garbled sentence claiming the three-body models \"include the contribution from the giant planet\" when in fact they miss it.\n\nOther soft spots: the ad hoc nature of the factor is acknowledged, which I respect, but it means the formula is a heuristic without a derivation; the reference frequencies have no error bars; and the claim of \"quantitatively correct estimates\" goes beyond what the evidence shows. No code or data, but for a proceedings that's a minor issue.\n\nWhat's genuinely good: the review is clear, the limits are correct, and the multi-perturber sum is a nice practical tool. The numerical methods are standard and the authors are honest about the ad hoc step.\n\nWho is this for? Someone who wants a closed-form frequency estimate for a binary with a giant planet for a quick parameter survey. With a proper head-to-head test, it could be a solid minor contribution. Without it, it's an interesting suggestion.\n\nMy recommendation: send to peer review only if the authors add a like-for-like comparison of LLM vs LL over a range of e_B in the ER3BP, and fix the text inconsistency. Otherwise it's a workshop note, not a journal paper.","headline":"A short proceedings paper that honestly proposes an ad hoc eccentricity factor for Laplace-Lagrange secular frequencies, but the new formula never gets an isolated head-to-head test.","tokens_in":11921,"tokens_out":3842,"would_cite":false,"duration_ms":36391,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A modified Laplace-Lagrange model yields quantitatively correct secular precession frequencies for test particles in multi-star systems.","keywords":["binary stars","secular perturbation theory","Laplace-Lagrange model","Heppenheimer model","restricted three-body problem","restricted four-body problem","secular resonance","exoplanets in binaries"],"falsifier":"Numerically integrate the planar restricted three-body problem for an equal-mass binary with $\\alpha=0.1$ and $e_B=0.6$, extract the apsidal precession frequency by Fourier analysis of $e\\sin\\varpi$ and $e\\cos\\varpi$, and compare it with $g_{\\mathrm{LLM}}$; a deviation that grows faster than $(1-e_B^2)^{-3/2}$ with $e_B$ would falsify the factorized ansatz.","tokens_in":10714,"feed_emoji":"🪐","tokens_out":9503,"duration_ms":80684,"temperature":0.7,"pith_summary":"The paper tries to establish that one closed-form expression—the modified Laplace-Lagrange model (LLM)—can replace direct numerical integration for the secular precession rate of a test particle in a binary or multiple star system. The formula multiplies the classical Laplace-Lagrange frequency by the Heppenheimer eccentricity factor $(1-e_B^2)^{-3/2}$, so it reproduces the standard LL model for a circular companion and the Heppenheimer model for a very distant companion. In the restricted four-body system HD 41004, the summed form of the model follows the numerically measured precession frequency outside the resonance gaps, while the single-perturber three-body models drift away as the semi-major axis ratio grows. This matters because it makes fast analytical scans of secular resonances possible for observed binaries and multiple-star systems.","feed_headline":"Modified formula tracks planet precession in eccentric binaries","feed_subtitle":"Adding a Heppenheimer eccentricity factor to Laplace-Lagrange theory matches four-body numerical runs.","key_machinery":"The load-bearing object is the modified Laplace-Lagrange frequency (LLM): the classical single-perturber secular frequency $g_{\\mathrm{LL}}=\\frac{1}{4}\\mu\\alpha^2 n b_{3/2}^{(1)}(\\alpha)$ multiplied by Heppenheimer's eccentricity factor $(1-e_B^2)^{-3/2}$. The Laplace coefficient $b_{3/2}^{(1)}(\\alpha)$, defined by an integral over the perturbing body's orbital phase, carries the full dependence on the semi-major axis ratio $\\alpha$, while the added factor carries the binary's eccentricity $e_B$. Equation (11) then exploits the linear additivity of Laplace-Lagrange secular theory, summing the modified single-perturber frequencies of every massive companion to obtain a multi-perturber estimate.","core_discovery":"The central assertion is that $g_{\\mathrm{LLM}} = \\frac{1}{4}\\mu\\alpha^2 n b_{3/2}^{(1)}(\\alpha)(1-e_B^2)^{-3/2}$ provides quantitatively correct estimates of a test particle's secular frequency in a binary system, and that summing this term over perturbing companions, as in Eq. (11), extends the estimate to four-body and multi-perturber systems. The construction is an interpolation: as $e_B\\to0$ the formula reduces to the Laplace-Lagrange frequency, and as $\\alpha\\to0$ it reduces to the Heppenheimer frequency. In the HD 41004 restricted four-body test the LLM tracks the numerical reference values, with only a small systematic overestimate at larger $\\alpha$, whereas the Heppenheimer, Giuppone, Andrade-Ines, and Georgakarakos models, all built for a single external perturber, fail to reproduce the secondary star's contribution in that regime.","pith_inferences":["If the factorized eccentricity correction is more than an interpolation, a similar $(1-e_B^2)^{-3/2}$ multiplier might apply to higher-order terms in the Laplace-coefficient expansion, a possibility the paper does not test.","The reported systematic overestimate at larger $\\alpha$ suggests the true frequency is bracketed by the LLM and the small-$\\alpha$ Heppenheimer value, so the two could be used together as an error estimate.","Since the summation in Eq. (11) is linear, the same formula could in principle also treat giant planets as perturbers of a massless test particle, extending the approach from multi-star to multi-planet secular dynamics."],"forward_implications":["Because the LLM reduces to the standard Laplace-Lagrange model when $e_B\\to0$ and to the Heppenheimer model when $\\alpha\\to0$, it provides one interpolation formula that covers both classical limits.","The linear summation in Eq. (11) applies to any number of perturbing companions, so multi-star systems with additional planets can be treated without new numerical integrations.","In the HD 41004 restricted four-body system the LLM matches the numerical secular frequencies outside the resonance gaps, while the three-body models deviate increasingly for larger $\\alpha$.","The formula makes fast parameter studies of secular resonance locations, including resonances inside habitable zones, feasible for large samples of observed binary systems."],"supporting_citations":[{"why":"Supplies the Heppenheimer model and the eccentricity factor $(1-e_B^2)^{-3/2}$ that the LLM grafts onto the Laplace-Lagrange frequency.","marker":"(Heppenheimer, 1978)"},{"why":"Provides the Laplace-Lagrange secular theory, the definition of the Laplace coefficient $b_{3/2}^{(1)}$, and the additivity property used in Eq. (11).","marker":"(Murray & Dermott, 1999)"},{"why":"Defines the HD 41004 configuration and identifies the secular resonance that motivates the four-body application.","marker":"(Pilat-Lohinger et al., 2016)"},{"why":"Surveys binaries with secular resonances near habitable zones; the LLM is intended to replace this survey's numerical frequency-determination step.","marker":"(Bazsó et al., 2017)"},{"why":"Supplies the GIU second-order model used as a comparison baseline in the figures.","marker":"(Giuppone et al., 2011)"},{"why":"Supplies the AND empirical correction model used as a comparison baseline in the figures.","marker":"(Andrade-Ines & Eggl, 2017)"},{"why":"Supplies the GEO second-order model used as a comparison baseline in the figures.","marker":"(Georgakarakos, 2002)"}],"fun_headline_variants":["New analytical model improves precession in eccentric binaries","Laplace-Lagrange extension tracks planets with multiple perturbers","Refined secular frequency formula handles eccentric companions","Four-body test backs modified secular precession estimate","Sharper precession predictions for planets in binary systems"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Heppenheimer eccentricity factor $(1-e_B^2)^{-3/2}$ can be attached as a simple multiplier to the full Laplace-Lagrange frequency, even though that factor is derived only in the small-$\\alpha$ quadrupole limit; if the true eccentricity dependence does not factor this way, the LLM will misestimate at large $e_B$ and larger $\\alpha$.","fun_headline_variants_meta":{"raw":{"variants":["New analytical model improves precession in eccentric binaries","Laplace-Lagrange extension tracks planets with multiple perturbers","Refined secular frequency formula handles eccentric companions","Four-body test backs modified secular precession estimate","Sharper precession predictions for planets in binary systems"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00061,"raw_usage":{"total_tokens":2798,"prompt_tokens":860,"completion_tokens":1938,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":476,"completion_tokens_details":{"reasoning_tokens":1877}},"tokens_in":476,"tokens_out":1938,"duration_ms":15936,"temperature":1.0,"reasoning_tokens":1877,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:25:07.600274+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate the planar restricted three-body problem for an equal-mass binary with $\\alpha=0.1$ and $e_B=0.6$, extract the apsidal precession frequency by Fourier analysis of $e\\sin\\varpi$ and $e\\cos\\varpi$, and compare it with $g_{\\mathrm{LLM}}$; a deviation that grows faster than $(1-e_B^2)^{-3/2}$ with $e_B$ would falsify the factorized ansatz.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Heppenheimer model and the eccentricity factor $(1-e_B^2)^{-3/2}$ that the LLM grafts onto the Laplace-Lagrange frequency."},{"cited_title":"D., & Dermott, S","cited_arxiv_id":null,"evidence_quote":"Provides the Laplace-Lagrange secular theory, the definition of the Laplace coefficient $b_{3/2}^{(1)}$, and the additivity property used in Eq. (11)."},{"cited_title":"2016, AJ, 152, 139","cited_arxiv_id":null,"evidence_quote":"Defines the HD 41004 configuration and identifies the secular resonance that motivates the four-body application."},{"cited_title":"A., Leiva, A","cited_arxiv_id":null,"evidence_quote":"Supplies the GIU second-order model used as a comparison baseline in the figures."},{"cited_title":"Secular Orbit Evolution in Systems with a Strong External Perturber - A Simple and Accurate Model","cited_arxiv_id":"1701.03425","evidence_quote":"Supplies the AND empirical correction model used as a comparison baseline in the figures."},{"cited_title":"2002, MNRAS, 337, 559 —","cited_arxiv_id":null,"evidence_quote":"Supplies the GEO second-order model used as a comparison baseline in the figures."}],"review_version":1}