{"id":"68387aa8-6448-4591-bc50-aa59063d0f38","arxiv_id":"1908.06331","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Massive circumbinary planet simulations match the secular analytic predictions for stationary and critical tilts, including new crescent retrograde librating orbits.","lead":"Three-body simulations of a nonzero-mass planet orbiting an eccentric binary confirm that the planet's stationary tilt and the threshold tilt for libration depend only on the binary eccentricity and the planet-to-binary angular momentum ratio. The results validate and extend the secular analytic model, and they reveal new crescent-shaped librating orbits for retrograde planets at high angular momentum.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative validation may depend on using initial eb and j even though both evolve strongly in the massive-planet runs; a cycle-averaged comparison would settle it.","rationale":"The paper is a solid numerical study: it uses an independent N-body integrator, systematically maps a wide parameter space, and confirms the qualitative predictions of Martin & Lubow (2019), including the new crescent retrograde librating orbits. The reader's identified weakest assumption, the quadrupole truncation, is explicitly tested in Fig. 7 and is shown to fail only very close to the binary, so it does not constitute a decisive objection. My concern is adjacent but more specific: the analytic formulas are parametrized by fixed eb and j, while the simulations of massive planets show large amplitude evolution in eb (and therefore in j). The paper compares to formulas evaluated at the initial values without demonstrating that this is the correct comparison. This does not invalidate the paper, because the observed agreement may still be genuine, but it is the least secure link in the central quantitative claim. A cycle-averaged re-analysis would settle whether the agreement is robust or parameterization-dependent. Since the qualitative conclusions and the main parameter map are well supported, I do not recommend changing the accept verdict, but the reproducibility concern should be addressed in a revision.","tokens_in":18519,"tokens_out":17981,"duration_ms":180552,"concrete_test":"For Model C2 or F3, compute the cycle-averaged binary eccentricity <eb> and the cycle-averaged angular-momentum ratio <j> over one full libration period from the simulation, then re-evaluate Eqs. (5), (10), and (11) with these averaged values and compare residuals against the simulation-derived is and imin. If the residuals grow substantially relative to the initial-parameter comparison, the paper's use of initial eb and j is load-bearing; if the residuals remain small, the initial-parameterization is validated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is that the three-body simulations agree with Eqs. (5), (10), and (11), with stationary and minimum libration tilts depending only on eb and j. Those analytic formulas are derived for a secular quadrupole Hamiltonian in which eb and j are fixed parameters. In the massive-planet simulations used to test this claim, however, eb changes substantially over a libration cycle: Fig. 2 shows eb approaching 1 for Models D2, F2, and related high-j cases, and j = L_p/L_b is not constant because L_b is proportional to sqrt(1 - e_b^2). The paper compares the simulation-derived angles with the analytic formulas evaluated at the initial eb and the initial j, but it does not justify why the initial values, rather than the instantaneous or cycle-averaged values, are the correct parameters for a system whose binary eccentricity is evolving. If the matching formulas should use the time-averaged eb or j, the claimed agreement could be partly an artifact of the chosen parameterization. The quadrupole truncation itself is acknowledged and tested in Fig. 7, but this separate parameter-evolution issue is not addressed. The paper also does not specify how the numerical stationary and critical angles are extracted from the phase-plane plots, which makes the comparison difficult to reproduce independently.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the three-body dynamics of a nonzero-mass, initially circular circumbinary planet around an eccentric binary using REBOUND N-body simulations. It maps prograde and retrograde circulating and librating orbits, determines the stationary inclination and the critical minimum (and maximum) libration angles, and compares the results with the analytic secular quadrupole formulas of Martin & Lubow (2019), namely Eq. (5) for the stationary tilt and Eqs. (10) and (11) for the minimum libration tilt. A broad parameter grid is covered, including binary eccentricities eb = 0.2, 0.5, 0.8, mass fractions fb = 0.5, 0.1, planet masses up to ~0.1 mb, and orbital radii 5 ab and 20 ab. The paper reports very good quantitative agreement with the analytic predictions and highlights new crescent-shaped retrograde librating orbits for high-angular-momentum planets. It closes with an application to the HD 106906 system.","tokens_in":18653,"tokens_out":7107,"duration_ms":83426,"significance":"If the claimed agreement holds, the paper provides a valuable numerical validation of the secular quadrupole framework for massive circumbinary planets and shows that the key angles depend only on eb and the angular momentum ratio j, rather than on the binary mass fraction or other parameters. The comparison is a genuine prediction test: the analytic formulas are evaluated without any fitted parameters, and the N-body integrations are independent of the analytic model. The paper also documents a new morphology of retrograde librating orbits and connects the results to circumbinary disc evolution and to a concrete observed system. The strengths are the breadth of the parameter survey, the parameter-free nature of the comparison, and the identification of qualitative features (crescent orbits, the jcr threshold, and the crossing of retrograde stationary curves) that are directly testable.","major_comments":[{"comment":"The analytic formulas (5), (10), and (11) treat eb and j as fixed parameters of a secular Hamiltonian. In the massive-planet simulations used for the comparison, however, eb evolves strongly: Fig. 2 shows eb approaching 1 for Models D2, F2, and related high-j cases, and j = Lp/Lb then changes because Lb is proportional to sqrt(1 - eb^2). The paper does not state whether the dots in Fig. 6 and Fig. 10 were placed using the initial, instantaneous, or cycle-averaged values of eb and j, nor does it justify why the chosen convention is the correct one. Since the central claim of \"very good quantitative agreement\" is a test of Eqs. (5), (10), and (11), this is load-bearing. Please specify the definition of j used for each dot and either provide a comparison using cycle-averaged values of eb and j or give an explicit argument for why the initial values are the exact parameters of the comparison.","section":"§3.1, Fig. 6 and Fig. 2"},{"comment":"The numerical extraction of is, imin, and imax is not reproducible from the text. The paper states that these angles were determined from the simulations, but it does not describe the algorithm: how the libration center is identified in the i-phi plane, how the separatrix boundary is interpolated from the finite grid of initial inclinations, or how the green/blue/magenta boundaries are assigned. Without this information, the reader cannot judge whether the quoted agreement is within the numerical uncertainty of the phase-plane classification. Please describe the extraction procedure and provide at least a grid-resolution uncertainty estimate for the plotted numerical points.","section":"§3.1 and §3.2, Figs. 6-10"}],"minor_comments":[{"comment":"The sentence describing the green and red curves in Figs. 8 and 9 is inconsistent with the definitions of Eqs. (10) and (11): the text labels Eq. (10) as green for chi < 0 and Eq. (11) as red for chi > 0, whereas the preceding paragraph defines Eq. (10) for chi > 0 and Eq. (11) for chi < 0. The figure captions use the correct assignment; please fix the text.","section":"§3.2, text after Eq. (11)"},{"comment":"The ratio j is central to the comparison, but the paper never explicitly defines it in terms of the model parameters. Please give the formula for j, for example j = (mp/mb) sqrt(a/ab) with the appropriate reduced-mass factors, so that the reader can reproduce the values used in Figs. 6 and 10.","section":"§2.1 and §3.1"},{"comment":"All simulations start with the same initial phase angle phi = 90 deg. Because the libration/circulation separatrix is classified along this single line in phase space, a short convergence test with a different initial phase angle would help confirm that the reported critical angles are independent of the initial phase, especially for the non-nested crescent orbits.","section":"§2.1, initial conditions"},{"comment":"The simulation-derived points are shown without error bars or point-size information. Given the finite inclination grid used to classify orbits, an uncertainty estimate on each numerical angle would make the claimed \"very good quantitative agreement\" more quantitative.","section":"Figs. 6, 7, and 10"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the central comparison is a genuine, parameter-free test of the analytic model. The main concern is that the quantitative comparison is not fully specified: the convention for eb and j and the numerical extraction of the angles need to be documented and, if necessary, corrected to cycle-averaged values. These are local but consequential issues, and I expect they can be fixed with a revised analysis. The self-citation to Martin & Lubow (2019) is appropriate and is not a circularity problem, since the simulations are independent of that model."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a useful, clearly written numerical study that does what it says—maps the orbital dynamics of a massive, circular circumbinary planet around an eccentric binary. The genuinely new pieces are the suite of three-body simulations with planet masses up to about 0.1 times the binary mass, the systematic extraction of stationary and critical inclinations, and the discovery of crescent-shaped retrograde librating orbits that appear when the planet-to-binary angular momentum ratio j is below the analytic threshold jcr.\n\nWhat the paper does well: the simulations are straightforward and documented (REBOUND, WHfast), with a sensible parameter grid over binary eccentricity, mass fraction, planet mass, and semi-major axis. They recover the test-particle limit of Doolin & Blundell, so the code is benchmarked. The comparison with the analytic formulas from Martin & Lubow is extensive—Figures 6, 8, 9, and 10 cover many models—and the agreement for low-j cases is good. The crescent orbits are a real qualitative finding, and the explanation via the absence of retrograde stationary states below jcr is coherent and consistent. The paper also correctly connects these results to circumbinary disc evolution and misaligned planets. The self-citation is not a problem here because the N-body integrations are an independent test of the analytic model.\n\nSoft spots: the stress-test note has teeth. The analytic formulas treat eb and j as fixed parameters, but in the massive-planet runs both evolve significantly over a libration cycle (see Figure 2, especially Models D2, E2, F2). The paper compares simulation-derived angles to the formulas evaluated at the initial eb and j without justifying why the initial values are the relevant ones. This matters most for the high-j models, where the claimed \"very good quantitative agreement\" is less clean. It is not a fatal flaw, but it weakens the headline validation claim. Also, the paper never states how the numerical stationary and critical angles are extracted from the phase-plane plots, which makes independent replication harder than it should be. Minor: some orbits are removed as \"unstable\" but the stability criterion is not specified.\n\nBottom line: for someone working on circumbinary planets or misaligned discs, this is a solid and citable contribution. It deserves peer review and likely publication after a minor revision that addresses the parameter-evolution issue and details the measurement procedure. I'd send it to a good referee.","headline":"Solid numerical extension of circumbinary planet dynamics to nonzero planet masses, broadly confirming the self-authored analytic model, with a real caveat about using initial eb and j in the high-mass comparisons.","tokens_in":19262,"tokens_out":2946,"would_cite":true,"duration_ms":27381,"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":"For massive circumbinary planets, stationary tilt and libration are set by binary eccentricity and angular momentum ratio alone.","keywords":["circumbinary planets","secular dynamics","orbital inclination","libration","three-body simulations","binary eccentricity","angular momentum ratio","quadrupole approximation"],"falsifier":"Run direct three-body integrations at a planet separation of $2.0\\,a_b$ (inside the unstable boundary but where quadrupole error should be larger) or at angular momentum ratios $j > 2$, measuring $i_s$ and $i_{\\rm min}$ and comparing to Eqs. (5), (10), and (11). A systematic offset larger than the scatter in Figs. 6–10, or a measurable dependence of these angles on binary mass fraction at fixed $e_b$ and $j$, would falsify the claim that $e_b$ and $j$ are the only controlling parameters.","tokens_in":18244,"feed_emoji":"🪐","tokens_out":7601,"duration_ms":69961,"temperature":0.7,"pith_summary":"This paper asks what determines the tilt of a massive planet orbiting an eccentric binary star on a circular, misaligned path. Through three-body simulations, it argues that two dimensionless numbers — the binary eccentricity $e_b$ and the ratio $j$ of planet-to-binary angular momentum — completely set the stationary tilt and the critical tilt separating libration from circulation. The simulations reproduce the analytic secular quadrupole formulas of Farago & Laskar (2010) and Martin & Lubow (2019) across a wide grid of masses, separations, and binary mass fractions. A sympathetic reader cares because observed circumbinary planets are nearly coplanar only by selection bias, and these formulas predict where tilted, librating planets should live.","feed_headline":"Two numbers control a circumbinary planet's tilt","feed_subtitle":"Simulations confirm analytic formulas: binary eccentricity and angular momentum ratio set tilt and libration.","key_machinery":"The argument is carried by the secular quadrupole approximation for the binary potential, which reduces the three-body secular problem to equations for the planet's inclination $i$ and nodal phase $\\varphi$ relative to the binary's eccentricity and angular-momentum vectors. Within this model, the load-bearing formulas are Eq. (5) for the stationary tilt $\\cos i_s$, Eq. (9) defining the branch parameter $\\chi$, and Eqs. (10)–(11) for the minimum libration tilt $\\cos i_{\\rm min}$. The comparison is made by plotting phase-space trajectories in the $i\\cos\\varphi$–$i\\sin\\varphi$ plane and reading off the centers and boundaries of libration islands.","core_discovery":"The paper's central claim is that for a nonzero-mass planet on an initially circular circumbinary orbit, the stationary planet-to-binary tilt $i_s$ and the minimum libration tilt $i_{\\rm min}$ depend only on the binary eccentricity $e_b$ and the planet-to-binary angular momentum ratio $j = |\\mathbf{l}_p|/|\\mathbf{l}_b|$; binary mass fraction and planet mass enter only through $j$. Direct N-body integrations spanning $e_b = 0.2, 0.5, 0.8$, binary mass fractions $f_b = 0.1$ and $0.5$, planet masses $m_p = 0.001\\,m_b$ to $0.116\\,m_b$, and separations $5\\,a_b$ to $20\\,a_b$ agree quantitatively with Eqs. (5) and (9)–(11). The paper also maps a previously unexplored regime: for $j$ below $j_{\\rm cr} = (1+4e_b^2)/(2+3e_b^2)$, no noncoplanar retrograde stationary state exists, and instead retrograde orbits librate in crescent-shaped paths of nonzero phase extent.","pith_inferences":["If the quadrupole map holds beyond the tested grid, the near-coplanarity of all Kepler-detected circumbinary planets is likely a transit selection effect; wide-orbit binaries should host planets with a continuum of tilts up to and including polar orbits.","The crescent retrograde orbits found for $j < j_{\\rm cr}$ suggest a dynamical boundary that could mark a zone of instability or bifurcation; a dedicated long-time stability map in that region would test this extension.","The same ratio-symmetric structure might extend to eccentric planet orbits or octupole-order terms, where the stationary tilt would acquire additional dependence on the planet's eccentricity and the binary mass ratio."],"forward_implications":["A measured circumbinary planet tilt, combined with known binary eccentricity, fixes the angular momentum ratio $j$ and thereby constrains the planet's mass and semi-major axis.","The analytic boundary between circulating and librating orbits lets observers predict which systems should show tilted, precessing planets rather than coplanar ones.","For retrograde planets with $j < j_{\\rm cr}$, the absence of a stationary state implies crescent-shaped librations with nonzero phase extent; such orbits may have distinct stability and observational signatures.","Because the same quadrupole formulas describe massive discs, the verified scalings transfer to circumbinary disc alignment, predicting the final tilt of debris discs and planets formed in them."],"supporting_citations":[{"why":"Derives the secular quadrupole approximation and the analytic precession and libration behavior for test particles, supplying the basis for the formulas.","marker":"Farago & Laskar 2010"},{"why":"Provides Eqs. (5), (9)–(11) for the stationary and minimum libration tilts with nonzero angular momentum ratio $j$.","marker":"Martin & Lubow 2019"},{"why":"Supplies the massless-particle numerical baseline that the paper verifies against and extends to massive planets.","marker":"Doolin & Blundell 2011"},{"why":"Supplies the WHfast symplectic integrator used for the three-body simulations.","marker":"Rein & Tamayo 2015"}],"fun_headline_variants":["Two numbers alone set a circumbinary planet's tilt","Binary eccentricity and momentum ratio dictate planet tilt","Retrograde circumbinary orbits lose stationary tilt at low j","Simulations confirm analytic tilt formulas for circumbinary planets","Tilt and libration of circumbinary planets follow two parameters"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the quadrupole truncation of the binary potential is accurate enough in the tested region that Eqs. (5), (10), and (11) serve as a faithful benchmark; the paper acknowledges the approximation degrades close to the binary, where it already sees deviations.","fun_headline_variants_meta":{"raw":{"variants":["Two numbers alone set a circumbinary planet's tilt","Binary eccentricity and momentum ratio dictate planet tilt","Retrograde circumbinary orbits lose stationary tilt at low j","Simulations confirm analytic tilt formulas for circumbinary planets","Tilt and libration of circumbinary planets follow two parameters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000347,"raw_usage":{"total_tokens":1937,"prompt_tokens":1017,"completion_tokens":920,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":633,"completion_tokens_details":{"reasoning_tokens":837}},"tokens_in":633,"tokens_out":920,"duration_ms":8427,"temperature":1.0,"reasoning_tokens":837,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:48:50.406388+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run direct three-body integrations at a planet separation of $2.0\\,a_b$ (inside the unstable boundary but where quadrupole error should be larger) or at angular momentum ratios $j > 2$, measuring $i_s$ and $i_{\\rm min}$ and comparing to Eqs. (5), (10), and (11). A systematic offset larger than the scatter in Figs. 6–10, or a measurable dependence of these angles on binary mass fraction at fixed $e_b$ and $j$, would falsify the claim that $e_b$ and $j$ are the only controlling parameters.","supporting_citations":[],"review_version":1}