{"id":"1a0a68e0-1d6c-4b44-8afd-eb891d51b490","arxiv_id":"2505.12124","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"For a Sun-like star with Jupiter-like and Neptune-like planets, dynamics changes sharply at a mutual inclination near 30-40 degrees, where the planets' precession frequencies become equal.","lead":"This paper maps how two giant planets orbiting a Sun-like star behave when their orbits are tilted and stretched. It finds that once the mutual tilt exceeds roughly 30 to 40 degrees, the motion switches from a calm, predictable regime to a wilder one driven by inclination changes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The regime boundary rests on SCC>0.05 spectral detections over 10 Myr, and the paper's own Sec. 3.3 admits that in high-(e,i) cases identifying the fundamental frequencies is more difficult; a threshold or integration-time artifact could masquerade as the g1=g2 transition.","rationale":"The paper's central claim requires that the frequencies detected with SCC>0.05 over 10 Myr are the true fundamental frequencies, and that the g1=g2 crossing is the controlling bifurcation. The reader identified this as the weakest assumption, and I agree: Sec. 3.3 explicitly admits that in high-(e,i) regimes spectral identification is difficult and every peak above a fixed threshold is recorded. Because the transition is inferred from the appearance of combination frequencies and the crossing of g1 and g2, a threshold or integration-window artifact could produce the same qualitative signature without a true change of regime. The paper also does not quantify uncertainties or test convergence, so the sharpness of the 30-40 degree boundary is not established. This is a genuine but addressable concern: the qualitative picture is consistent with prior work, and the paper's numerical maps are rich and carefully produced, but the quantitative regime boundary needs robustness checks before the claim is accepted as definitive. The reader's CONDITIONAL verdict is therefore appropriate, and my stress test does not move it.","tokens_in":16090,"tokens_out":4843,"duration_ms":49217,"concrete_test":"Choose a representative configuration from Fig. 10 (e.g., a_N=4 au, e_i=0.05, initial omega_i=90 deg) and repeat the frequency analysis for initial mutual inclinations i=20, 28, 32, 36, and 44 deg using integration times T=1, 10, and 100 Myr and SCC thresholds 0.01, 0.05, and 0.1, recovering frequencies from (k,h) and (q,p) with the same pipeline. Also analyze a synthetic quasi-periodic signal with known frequencies and amplitudes comparable to the system's modes. If the inferred ic (the g1=g2 crossing) shifts by more than a few degrees with T or SCC, or if combination frequencies with SCC>0.05 appear below ic for longer T, then the regime boundary is a detection artifact rather than a robust dynamical transition. In addition, verify for the same runs that mutual inclination remains constant for i0.05.' A fixed SCC cutoff and finite integration window can systematically drop weak fundamental peaks and retain combination peaks as the dynamics become more complex, producing an apparent emergence of combination frequencies above ic even if no qualitative change occurs at exactly g1=g2. The paper gives no error bars on frequencies, no convergence test in integration time, and no validation of the spectral pipeline on a synthetic signal with known frequencies. It also asserts that mutual inclination is nearly constant for i i_c combinations of frequencies appear, the mutual inclination changes, and the dynamics is dominated by the secular resonance or by the vZLK mechanism. The paper also documents two families of pseudo-equilibrium configurations: low-inclination equilibria at Δϖ = 0° or 180°, and high-inclination equilibria at ω1 = ω2 = K·90°.","tokens_in":16367,"tokens_out":2202,"duration_ms":24417,"significance":"If correct, the proposed 30–40 degree boundary would provide a simple, physically interpretable organizing principle for two-planet secular dynamics, connecting the g1 = g2 resonance with the onset of vZLK-like evolution, and it would strengthen earlier stability studies that place the unstable/inclined boundary near 30–40 degrees. The paper is valuable for its numerical atlases (Figs. 1–2, 4), its explicit use of direct integrations, and its careful labeling of pseudo-equilibria and finite-time limitations. The semi-analytical Rs topology calculation (Figs. 5–7) is independent of planetary masses and provides a useful geometric characterization. However, the central regime boundary rests on spectral detections with a fixed SCC > 0.05 threshold over 10 Myr, and the paper itself notes (Sec. 3.3) that in high-(e,i) cases identifying fundamental frequencies is difficult; without convergence tests or uncertainty quantification, the sharpness and exact location of the boundary remain not fully established.","major_comments":[{"comment":"The central claim that a sharp two-regime boundary occurs at the g1 = g2 crossing rests entirely on frequencies extracted from 10 Myr integrations with a fixed spectral correlation coefficient threshold (SCC > 0.05). The paper states in Sec. 3.3 that in high-(e,i) cases 'several frequencies appear in the spectra... so we simply record all those with SCC>0.05.' There is no convergence test in integration time, no test of the spectral pipeline on a synthetic signal with known frequencies, and no uncertainty estimate for the detected frequencies. A fixed SCC cutoff can systematically drop weak fundamental peaks and retain combination peaks as the dynamics complexify, producing an apparent emergence of combination frequencies above i_c even if no qualitative change occurs exactly at g1 = g2. Since the g1 = g2 crossing is the definition of the boundary, this is load-bearing and needs to be addressed.","section":"Sec. 3.3, Figs. 8–10"},{"comment":"The claimed critical inclination range 30 < i_c < 40 is inferred from a finite set of initial conditions and displayed as frequency-versus-inclination curves without uncertainty quantification. The paper states that 'in all experiments we performed for the range of aN we studied is located between 30 and 40 degrees' (Sec. 3.3), but no statistical or systematic error bars are given, and the number of experiments is not specified. Moreover, Fig. 11 shows that i_c varies with a_N and eccentricity, so the statement that 'in general' i_c is between 30 and 40 needs a quantitative characterization of the range of variation, especially for a_N/a_J near the chaotic boundary, before it can support a universal regime split.","section":"Sec. 3.3, Figs. 8 and 11"},{"comment":"The equilibrium-topology switch in Figs. 5–6 occurs at mutual inclinations below 30 degrees (e.g., 15 < i < 30 in the text describing Fig. 5), which is distinct from the g1 = g2 boundary at 30–40 degrees used for the regime split. The paper says the g1 = g2 crossing is 'not equal but linked to' the change in the distribution of equilibrium points (Sec. 3.3), but no dynamical mechanism or quantitative relation is provided. This leaves the central claim with two separate, potentially inconsistent critical inclinations; a demonstration of how the equilibrium switch and the frequency crossing are related is needed.","section":"Figs. 5–6 and Sec. 3.3"},{"comment":"The identification of the vZLK onset with the condition g2 = f is based on the same finite-time spectral analysis, and the paper reports that after g2 = f the frequency g2 'disappears or is merged with f' and f shows a discontinuity. These observations are qualitative and depend on the SCC threshold. Since a secular resonance and a vZLK regime have distinct phase-space signatures, the paper should verify at least one representative high-inclination case with longer integrations and/or with a direct phase-space diagnostic (e.g., libration of ωN in the appropriate variables) to confirm that the spectral transition corresponds to a true dynamical transition rather than to a loss of spectral resolution.","section":"Sec. 3.3, Fig. 8"}],"minor_comments":[{"comment":"There are several typographical and grammatical issues, e.g., 'an a Neptune-like planet', 'with1< aN <20 au' (missing space), 'the dynamics is analogue a the classic secular model', and 'intial' in Fig. 9. These should be corrected.","section":"Abstract and throughout"},{"comment":"The notation i for mutual inclination is introduced informally in Sec. 2 and used extensively, but a formal definition relating i to i1 and i2 in the invariable Laplace plane appears only in Sec. 3.2. It would help to define the mutual inclination explicitly when the variables are first introduced.","section":"Sec. 3.1"},{"comment":"Fig. 6 is described as showing the 'relative minimal mutual inclination' for the onset of the new equilibria, but the caption and text do not explain how this curve is computed from the Rs topography (e.g., what threshold in Rs is used to identify the onset). A brief description in the caption would improve reproducibility.","section":"Fig. 6"},{"comment":"The MEGNO maps use a time span of 5000 orbital revolutions of the test Neptune; the paper notes this may be insufficient to detect chaos in some regions. It would be useful to state the MEGNO threshold used to classify regular versus chaotic, since the color scale is not quantified in the caption.","section":"Sec. 2.1, Fig. 4"},{"comment":"The term 'SCC' is defined in the text, but the numerical threshold 'SCC>0.05' is used without discussion of its sensitivity. A brief statement of how the threshold was chosen and how results change for, say, SCC>0.1 or SCC>0.01 would be valuable.","section":"Sec. 3.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid numerical exploration, but the central claim about the sharp 30–40 degree boundary needs stronger support than a fixed SCC threshold over a fixed 10 Myr window. The authors should be encouraged to add convergence tests and uncertainty estimates; without those, the boundary could be an artifact of the spectral pipeline. I would not reject the paper, because the underlying integrations and atlases are useful and the regime split is plausible, but the central claim currently rests on a fragile inference."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, readable numerical survey of two-planet dynamics in a parameter regime that has been less charted than the hierarchical case. The headline result—two secular regimes split near 30–40 degrees mutual inclination—is not new; Michtchenko et al. 2006, Funk et al. 2011, Volpi et al. 2019, and Mastroianni & Efthymiopoulos 2023 all point the same way. What is new is the systematic non-hierarchical, comparable-mass mapping and the explicit identification of g1=g2 as the dividing line, which is a useful concrete criterion.\n\nThe paper does real work well. The MMR atlas of ~1300 resonances is a handy reference. The pseudo-equilibrium analysis is clearly labeled and honestly restricted to the omega subspace. The frequency analysis is careful to flag the high-(e,i) cases where spectral identification gets messy. The authors also take the time to warn against conflating apsidal libration in quasi-coplanar systems with a genuine secular resonance—good pedagogy.\n\nThe soft spots are real but not disabling. No code or data are released, and the frequency extraction has no error bars, no convergence test in integration time, and no validation against synthetic signals with known frequencies. The SCC>0.05 cutoff is arbitrary enough that the appearance of combination frequencies just above i_c could be partly a numerical effect; the paper acknowledges the difficulty in Sec. 3.3 but does not address it. More substantively, the pseudo-equilibrium topology switch in Figs. 5–6 happens at i < 30°, while the claimed g1=g2 boundary is at 30–40°; the paper calls the two \"linked\" but never shows the connection. That gap weakens the claim that g1=g2 is the actual regime boundary rather than a correlate.\n\nI don't think any of this sinks the central picture. For a first-pass map of a messy regime, the conclusions are plausible and fit the prior literature. But the paper would be stronger with a couple of longer integrations and a synthetic-signal test of the spectral pipeline, plus a more nuanced statement of what is new relative to the cited work. I'd send it to a good referee—the numeric survey is a service to the community—but I'd push for the UQ and a sharper novelty statement before accepting.","headline":"A workmanlike non-hierarchical survey of two-planet secular dynamics; the two-regime picture is not new, but the explicit g1=g2 criterion and the broad parameter map make it a useful reference—worth peer review, with requests for uncertainty quantification and a sharper novelty statement.","tokens_in":16953,"tokens_out":3234,"would_cite":true,"duration_ms":32502,"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 two comparable-mass giant planets, the secular dynamics has two regimes separated by a critical mutual inclination between about 30° and 40°, marked by the secular resonance g1 = g2.","keywords":["planetary dynamics","secular dynamics","resonances","vZLK mechanism","fundamental frequencies","mutual inclination","mean-motion resonances","two-planet systems"],"falsifier":"Integrate the working system (Sun-like star, $m_J=0.001$, $m_N=m_J/10$, $a_J=8$ au, $a_N$ from 4 to 16 au, $e\\simeq0.05$–$0.2$) for 10 Myr at mutual inclinations of 25°, 30°, 35°, 40°, and 45° and measure the fundamental frequencies. If the $g_1=g_2$ crossing is not found between 30° and 40° for a configuration with no MMR overlap, or if the mutual inclination varies strongly for some stable case with $i<30°$, the claimed regime boundary fails. A direct look at the frequency map would settle it.","tokens_in":2502,"feed_emoji":"🪐","tokens_out":2575,"duration_ms":96107,"temperature":0.7,"pith_summary":"The paper aims to establish a clean organizing principle for the long-term evolution of a star with two giant planets of comparable mass: the dynamics splits into two regimes according to the initial mutual inclination. Below a critical inclination $i_c$, located between about 30° and 40° and identified with the secular resonance $g_1=g_2$, the system behaves like the classic low-eccentricity, low-inclination secular model: three fundamental frequencies describe the motion and the mutual inclination stays nearly constant. Above $i_c$, combinations of those frequencies appear, the mutual inclination can change substantially, and either the secular resonance or the von Zeipel–Lidov–Kozai mechanism drives the evolution. The same dividing line appears in the structure of secular equilibria, which switch from $\\Delta\\varpi=0^\\circ/180^\\circ$ at low inclination to $\\omega_1=\\omega_2$ equal to integer multiples of 90° at high inclination. A reader should care because this offers a rule of thumb for which of two very different dynamical descriptions applies to a given inclined giant-planet pair.","feed_headline":"At 30–40° tilt, two-planet dynamics flips regime","feed_subtitle":"Same boundary separates calm secular motion from Kozai-driven orbital change.","key_machinery":"The load-bearing tool is frequency analysis of long numerical integrations: the authors integrate the exact equations for 10 Myr, extract periodicities with a spectral correlation coefficient (SCC) cutoff, and track the fundamental frequencies $g_1$, $g_2$, $f$ as functions of initial mutual inclination, eccentricity, semi-major axis, and mass ratio. The crossing $g_1=g_2$ is the event that defines the regime boundary, and the crossing $g_2=f$ marks the onset of the vZLK regime. Supporting this is a semi-analytical secular Hamiltonian $\\mathcal{R}_s(\\omega_1,\\omega_2)$ obtained by numerical double averaging over the mean anomalies; its extrema locate the pseudo-equilibrium configurations and show the geometric switch from $\\Delta\\varpi$-type to $\\omega_i=K\\,90^\\circ$-type equilibria. The assumption of an invariable Laplace reference plane (nodes opposed, $\\Delta\\Omega=180^\\circ$) reduces the system to two degrees of freedom while still allowing three fundamental frequencies in an inertial frame.","core_discovery":"The central claim is that the regime boundary is the condition $g_1=g_2$, i.e., the equality of the two secular precession rates of the planetary pericenters, which the authors observe to occur in their explored configurations only when the initial mutual inclination lies between roughly $i_c=30^\\circ$ and $40^\\circ$. Below $i_c$, three well-defined fundamental frequencies exist — $g_1$, $g_2$ associated with the pericenters and eccentricities, and $f$ associated with the common line of nodes and inclinations — and the mutual inclination is essentially conserved. Above $i_c$, spectral analysis shows combinations such as $f+(f-g_2)$, the identities of the $g_i$ become entangled, and the mutual inclination changes; the vZLK mechanism begins when $g_2=f$, with $\\omega_N$ librating. The authors further find that the equilibrium configurations of the secular Hamiltonian change at the same transition: at low inclination, minima occur for $\\Delta\\varpi=0^\\circ$ or $180^\\circ$ regardless of the individual $\\omega_i$, while at high inclination they occur only for $\\omega_1=\\omega_2=K\\,90^\\circ$. This is established through numerical integrations of the full equations of motion over 10 Myr, spectral decomposition with a spectral correlation coefficient threshold, and a semi-analytical secular Hamiltonian computed by double averaging.","pith_inferences":["A testable extension of the paper's picture: for a given pair of giant planets, the critical inclination should be measurable from observed orbital precession rates as the mutual inclination at which the apsidal precession rates cross; if real systems near $i\\simeq 35^\\circ$ show no such frequency crossing, the boundary would need revision.","The two regimes suggest a practical classification for exoplanet stability: systems with mutual inclination below $i_c$ can be modeled with quasi-constant-inclination secular theory, while those above it require vZLK-capable integration; this could guide choices of which systems to flag for strong mutual perturbations.","The equilibrium switch to $\\omega_i=K\\,90^\\circ$ implies that high-inclination giant-planet pairs may preferentially settle with apses aligned or perpendicular to the mutual line of nodes, a statistical preference that could be searched for in observed multi-planet systems with measured mutual inclinations.","The paper excludes secular evolution inside mean-motion resonances; extending the frequency analysis to resonant populations might reveal that the $g_1=g_2$ boundary shifts inside resonance, connecting the MMR atlas to the secular regime split."],"forward_implications":["Below the critical mutual inclination $i_c$, a two-giant-planet system keeps three well-separated fundamental frequencies, the mutual inclination stays nearly constant, and the eccentricity dynamics reduces to forced plus free modes; above $i_c$, frequency combinations appear and the mutual inclination can change substantially over time.","The locus of secular equilibria shifts from $\\Delta\\varpi=0^\\circ$ or $180^\\circ$ at low mutual inclination to $\\omega_1=\\omega_2=K\\,90^\\circ$ (apses aligned or perpendicular to the line of nodes) at high mutual inclination, so the orientation of the apsidal lines relative to the nodes becomes dynamically meaningful only in the high-inclination regime.","The von Zeipel–Lidov–Kozai mechanism sets in only above the critical inclination, and the onset inclination is higher for an exterior Neptune-like planet and higher still when the mass ratio $m_N/m_J$ is smaller.","Systems with mutual inclination above roughly 30° are predisposed to large changes in eccentricity and inclination, and for an exterior planet this can lead to instability, consistent with earlier stability studies that place the onset of inclination-driven instability around 30°–40°.","The atlas of roughly 1300 mean-motion resonances shows that for a Neptune-like planet between about 4 and 16 au with mutual inclination below about 30°, resonance overlap dominates, so the secular regime split applies primarily in the white, MMR-free regions of the atlas."],"supporting_citations":[{"why":"Supplies the 3D secular framework, the reduction to two degrees of freedom in the Laplace plane, and the domain structure that the paper's regimes are compared with.","marker":"Michtchenko et al. (2006)"},{"why":"Establishes the numerical secular disturbing function method and the apsidal Mode I/II equilibria that the low-inclination results extend.","marker":"Michtchenko and Malhotra (2004)"},{"why":"Provides the hierarchical quadrupole/octupole results with $\\Delta\\varpi=0^\\circ/180^\\circ$ equilibria that serve as the low-inclination baseline.","marker":"Lee and Peale (2003)"},{"why":"Provides the semianalytical resonant Hamiltonian model used to compute the ~1300 mean-motion resonance domains and overlap maps.","marker":"Gallardo et al. (2021)"},{"why":"Gives analytical expressions for the fundamental secular frequencies, motivating the identification of $g_1$, $g_2$, and $f$.","marker":"Libert and Henrard (2008)"},{"why":"Supplies the frequency-analysis method with spectral correlation coefficients used to extract the fundamental frequencies from the integrations.","marker":"Ferraz-Mello (1981)"},{"why":"Applies time-frequency analysis to libration, providing the basis for interpreting the detected periodicities as dynamical modes.","marker":"Gallardo and Ferraz-Mello (1997)"},{"why":"Independent stability study locating inclination-driven instabilities around 30°–40°, supporting the paper's critical-inclination finding.","marker":"Funk et al. (2011)"},{"why":"Identifies planar-like and Lidov-Kozai regimes in hierarchical three-body systems with high mutual inclination, whose regime terminology is adopted.","marker":"Mastroianni and Efthymiopoulos (2023)"},{"why":"Provides the MEGNO chaos indicator used to build the regularity maps that connect resonance overlap, chaotic regions, and the secular regimes.","marker":"Cincotta and Simó (2000)"}],"fun_headline_variants":["Two-planet systems flip dynamics at 30–40° mutual tilt","Critical tilt flips two-planet dynamics between 30–40°","Secular resonance sets boundary at 30–40° tilt","Two-planet orbits shift regime above ~35° mutual tilt","Inclined two-planet systems show two regimes, split at 30–40°"],"cache_read_input_tokens":18944,"weakest_assumption_plain":"The classification into two regimes rests on the assumption that the frequencies measured in 10-Myr integrations with spectral correlation coefficient above 0.05 are the system's true fundamental frequencies, and that below the critical inclination the mutual inclination is constant enough for the regime split to be sharp.","fun_headline_variants_meta":{"raw":{"variants":["Two-planet systems flip dynamics at 30–40° mutual tilt","Critical tilt flips two-planet dynamics between 30–40°","Secular resonance sets boundary at 30–40° tilt","Two-planet orbits shift regime above ~35° mutual tilt","Inclined two-planet systems show two regimes, split at 30–40°"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000634,"raw_usage":{"total_tokens":3007,"prompt_tokens":1111,"completion_tokens":1896,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":727,"completion_tokens_details":{"reasoning_tokens":1800}},"tokens_in":727,"tokens_out":1896,"duration_ms":12992,"temperature":1.0,"reasoning_tokens":1800,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:40:27.916277+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Integrate the working system (Sun-like star, $m_J=0.001$, $m_N=m_J/10$, $a_J=8$ au, $a_N$ from 4 to 16 au, $e\\simeq0.05$–$0.2$) for 10 Myr at mutual inclinations of 25°, 30°, 35°, 40°, and 45° and measure the fundamental frequencies. If the $g_1=g_2$ crossing is not found between 30° and 40° for a configuration with no MMR overlap, or if the mutual inclination varies strongly for some stable case with $i<30°$, the claimed regime boundary fails. A direct look at the frequency map would settle it.","supporting_citations":[],"review_version":1}