{"id":"bddb0c8a-0533-420d-bfe2-b030773139a6","arxiv_id":"2501.12258","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"N-body simulations show proposed exomoons around Kepler-1625, Kepler-1708, and HD 23079 are mostly orbitally stable, while submoons destabilize above about 40 degrees inclination and show a secular 13:2 precession resonance feature.","lead":"This paper maps where hypothetical moons and moons-of-moons would stay in orbit around three exoplanet systems, using computer simulations with varied tilt angles. It finds that the proposed exomoon orbits are mostly stable, and that a submoon would face instability above about 40 degrees of tilt, with some hint of a slow precession resonance.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The secular-resonance identification for the submoon ridges is not demonstrated: Eq. (3) is an ad hoc 13:2 commensurability, no resonant-angle libration is shown, and Fig. 9 only matches amplitude, not frequency.","rationale":"The reader's weakest assumption (no tides, adopted exomoon parameters) is real and is acknowledged in the paper itself, but it is a scope limitation: the maps are explicitly pure Newtonian 10^5-yr integrations. The more damaging issue for the paper's novel contribution is that the secular-resonance identification is not independently established. The 13:2 commensurability in Sec. 3.2.2 is extracted from the same simulations it is used to explain, the resonant angle is never shown to librate, and the secular-theory check (Fig. 9) matches only amplitude. This is a circularity/correctness risk, not just a scope limitation. I therefore give it as the load-bearing concern. The 3-body stability results for Kepler-1625 (Fig. 3a) are much more robust and stand independently, so a full rejection is not warranted; the resonance claim should be either supported with a libration test or downgraded to 'near-commensurability'. Since the reader already issued CONDITIONAL, no verdict change is needed.","tokens_in":15915,"tokens_out":6784,"duration_ms":69971,"concrete_test":"Re-run the Kepler-1625 4-body case of Fig. 7 (a_sm = 0.00084 au, i_sm = 3.6°, MA = 0°) for at least 2000 yr and compute \\phi(t) = 13\\omega_m(t) − 2\\omega_sm(t). If \\phi circulates through 2π rather than librating about a fixed value, the secular-resonance interpretation is unsupported. As a control, repeat for an off-ridge point (e.g., a_sm = 0.00080 au, same inclination) and verify that the ridge locations in Fig. 5 align with librating \\phi, not merely with a near-13:2 frequency ratio.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's novel 4-body claim — that the elevated-eccentricity ridges in Fig. 5(a) are secular resonances where the misaligned moon and submoon precess together — rests on an unproven commensurability. Section 3.2.2 reports measured precession rates \\dot{\\omega}_m ≈ 1.22°/yr and \\dot{\\omega}_sm ≈ 8.026°/yr, whose ratio (≈6.58) is close to 13:2, and defines the resonant-angle derivative in Eq. (3) as \\dot{\\phi} = 13\\dot{\\omega}_m − 2\\dot{\\omega}_sm. With the quoted numbers, \\dot{\\phi} ≈ −0.19°/yr, so \\phi would circulate with a period of roughly 1900 yr. The paper never shows \\phi(t) librating; the 100-yr simulations cover only a small fraction of that circulation period, so a slow drift and a trapped libration are indistinguishable. The supporting secular-theory comparison (Fig. 9) is explicitly limited to matching the inclination amplitude, with the authors noting a frequency/phase difference. Thus the central resonance identification is currently a fit to the simulation data rather than an independent dynamical confirmation. This does not invalidate the 3-body stability maps, but it undercuts the strongest novel conclusion about submoon dynamics.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the rebound N-body package to map orbital stability for hierarchical 3-body (star-planet-moon) systems for Kepler-1625, Kepler-1708, and HD 23079, and for 4-body (star-planet-moon-submoon) systems for Kepler-1625 and HD 23079. The authors validate their integration setup against known von Zeipel-Lidov-Kozai behavior, identify and patch a coordinate-frame issue in rebound's hierarchical-element initialization, and produce log maximum-eccentricity maps over satellite semimajor-axis and inclination grids. They interpret elevated-eccentricity ridges in the submoon stability map as secular resonances between the host moon and submoon precession, supported by measured apsidal precession rates, a claimed 13:2 commensurability, MEGNO maps, and a secular perturbation theory comparison. The main conclusions are that the proposed Kepler-1625 exomoon parameters are mostly orbitally stable and that certain submoon configurations remain gravitationally bound over 10^5 yr.","tokens_in":16292,"tokens_out":3566,"duration_ms":39143,"significance":"If confirmed, the paper would provide a useful framework for evaluating exomoon and submoon stability with inclination, and it would identify a specific secular-resonance mechanism for submoon eccentricity excitation. The manuscript has several genuine strengths: the authors compare multiple integrators, explicitly document and patch a coordinate-frame bug in rebound, randomize mean anomalies to reduce phase bias, verify their secular-theory implementation against Murray and Dermott's Jupiter-Saturn test case, and use MEGNO as an independent chaos indicator. These elements make the 3-body stability maps and the general 4-body instability structure credible and reproducible. The main weakness is that the central novel claim, the identification of the Fig. 5(a) ridges as secular resonances, is not demonstrated with a librating resonant angle or a matched secular frequency; it rests on an approximate commensurability of fitted precession rates.","major_comments":[{"comment":"The secular-resonance identification is not demonstrated. The coefficients 13 and 2 in Eq. (3) are chosen because the measured precession-rate ratio omega_dot_sm/omega_dot_m ~ 6.58 is close to 13/2, but the paper does not show that the resonant angle phi librates or remains bounded on any timescale. With the quoted rates omega_dot_m = 1.22 deg/yr and omega_dot_sm = 8.026 deg/yr, Eq. (3) gives phi_dot ~ -0.19 deg/yr, corresponding to a circulation period of roughly 1900 years; the 100-year simulations in Fig. 7 cover only a small fraction of that cycle, so slow circulation and trapped libration are indistinguishable. The authors should show phi(t) over at least a full cycle, report the libration width or a stroboscopic map, and test whether the 13:2 choice is dynamically preferred over nearby integer combinations.","section":"Sec. 2.3, Eq. (3); Sec. 3.2.2"},{"comment":"The secular perturbation theory comparison does not confirm the resonance because it matches only the amplitude of the submoon inclination variation and not its frequency or phase, as the text and figure caption admit. Since a secular resonance is defined by a commensurability of precession frequencies, an amplitude match is insufficient evidence. The authors should compare the dominant periods or eigenfrequencies of the secular theory with the N-body result, or otherwise show that the frequency mismatch arises from a benign non-secular effect, before concluding that the ridges in Fig. 5(a) are secular resonances.","section":"Sec. 2.4; Fig. 9"},{"comment":"The paper's own cited literature limits the physical interpretation of the 4-body stability maps. Section 4 cites Kollmeier and Raymond (2018) and Rosario-Franco et al. (2020) showing that tidal migration would preferentially remove submoons, but the 10^5-year simulations in Sec. 2.1 are pure Newtonian and include no tides. The stable regions in Fig. 5 are therefore gravitational-stability regions, not necessarily regions where a submoon could survive in the real Kepler-1625 or HD 23079 systems. The conclusion that a submoon 'could be orbitally stable' should be explicitly qualified as stability in the absence of tides, with a discussion of whether the tidal migration timescale is longer than the 10^5-year simulation time or the age of the system.","section":"Sec. 4, final discussion of submoon viability"}],"minor_comments":[{"comment":"The phrase 'pursing orbital stability analyses' contains a typo; it should be 'pursuing.'","section":"Abstract"},{"comment":"In the HD 23079 4-body row, the entry '0.0298 9' is visually ambiguous; the moon semimajor axis and inclination should be separated or labeled so that the reader can identify a_m = 0.0298 au and i_m = 9 deg without guessing.","section":"Table 1"},{"comment":"The caption says 'Kepler-162' where 'Kepler-1625' is meant.","section":"Fig. 6 caption"},{"comment":"The statement that the precession-rate ratio is '~13:2' is imprecise; the quoted values give 8.026/1.22 = 6.58, which is close to but measurably different from 13/2 = 6.50, and the paper should quantify this difference and discuss whether it is within the fitting uncertainty.","section":"Sec. 3.2.2"},{"comment":"The description of the coordinate-frame bug says the second transformation was not applied, but the precise failure mode in rebound's hierarchical initialization would be clearer if the authors stated whether the bug affected the submoon's inclination relative to the moon or its longitude of ascending node, since the correction in Appendix A addresses only inclination rotations.","section":"Sec. 2.1"}],"recommendation":"major_revision","confidential_remarks":"No additional issues beyond those stated; the manuscript is within the scope of a planetary-dynamics journal, and the stability-map portion is sound enough to warrant revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe useful part of this paper is the inclined stability maps. The 3-body maps for Kepler-1625, Kepler-1708, and HD 23079 are straightforward N-body products with standard integrator choices, and the 4-body submoon maps for Kepler-1625 and HD 23079 are new. The authors also found and patched a coordinate-frame bug in rebound when extracting hierarchical elements for the submoon, which is a good sign for the numerics. The MEGNO chaos maps add a useful cross-check. I'd trust the stability boundaries at the level expected of 10^5-year integrations.\n\nThe soft spot is the resonance claim. Section 3.2.2 identifies a 13:2 secular resonance between the moon and submoon precession rates. The coefficients are chosen because the measured ratio is close to 6.5, and Eq. (3) gives phi_dot about -0.19 deg/yr, i.e., a circulation period around 1900 years. The 100-year simulations cover less than 20 degrees of phase, so a slow drift and a trapped libration are indistinguishable. The paper never shows phi(t) librating. The secular-theory comparison in Fig. 9 matches the inclination amplitude but has a frequency offset, and the authors note this. So the 'we identify these resonances as secular' conclusion is not supported. The elevated-eccentricity ridges in Fig. 5(a) are real features, but calling them a resonance is an interpretation that needs more work.\n\nAlso, the 4-body runs ignore tidal dissipation. The paper cites Kollmeier & Raymond and Rosario-Franco showing tides would remove submoons, so the stable submoon regions may be dynamically but not physically viable. The authors acknowledge this, so it's a limitation, not an oversight. The HD 23079 exomoon initial condition is arbitrary, but that's fine for a framework paper.\n\nMy overall take: this is a useful contribution to the exomoon-stability subfield, with a solid core of stability maps and a resonance interpretation that should not be repeated as established. The paper deserves a serious referee: the maps are worth publishing, and the resonance claim needs to be pushed on in review. I'd recommend the editor send it to review, with a request for the resonant-angle time series and a demonstration of libration over at least one circulation period, or a softening of the resonance language.","headline":"Useful inclined stability maps for the leading exomoon candidates; the 13:2 secular resonance claim is a fit, not a demonstrated resonance.","tokens_in":16839,"tokens_out":4195,"would_cite":true,"duration_ms":37857,"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":"This paper maps where moons and submoons around exoplanets stay bound over 100,000 years, and finds that the proposed exomoon orbits in Kepler-1625 and Kepler-1708 fall mostly in stable regions.","keywords":["exomoon stability","submoon","secular resonance","von Zeipel-Lidov-Kozai effect","N-body simulations","Kepler-1625","Kepler-1708","HD 23079"],"falsifier":"Re-run the 4-body stability maps with tidal dissipation included, using tidal quality factors comparable to those in Kollmeier & Raymond (2018), and check whether any submoon in the 20-33% moon-Hill-radius range remains bound for 100,000 years; if none do, the claimed stable regions are not physically real.","tokens_in":15703,"feed_emoji":"🌙","tokens_out":10249,"duration_ms":88029,"temperature":0.7,"pith_summary":"This paper uses 100,000-year N-body simulations to map where a moon around an exoplanet, and a moon-of-a-moon (a submoon) around that moon, can remain bound, treating orbital size and inclination as free parameters. Applying the maps to the two strongest exomoon candidates, it finds that the observationally proposed orbits of Kepler-1625b-i and Kepler-1708b-i lie mostly in stable, low-eccentricity regions. For a hypothetical submoon around Kepler-1625b-i, it identifies stable islands alongside ridges of elevated eccentricity, which it attributes to a secular resonance where the moon's and submoon's orbits precess together at a 13:2 rate ratio. The result is a general framework for screening exomoon and submoon candidates against dynamical stability before investing in follow-up observations.","feed_headline":"Proposed exomoons pass a 100,000-year stability test","feed_subtitle":"Simulations map where moons and submoons stay bound, and where a secular resonance drives eccentricity up.","key_machinery":"The central machinery is a grid of stability maps built from 100,000-year N-body integrations that vary the satellite's semimajor axis and inclination, using maximum eccentricity as a stability proxy and checking orbit crossing beyond the Hill radius. The resonance claim rests on the resonant-angle derivative $\\dot{\\phi} = 13\\dot{\\omega}_{\\rm m} - 2\\dot{\\omega}_{\\rm sm}$, computed by linear fits to the moon's and submoon's arguments of pericenter; secular perturbation theory (the classical disturbing-function expansion) reproduces the amplitude of the submoon's inclination oscillation, which anchors the identification of the ridges as secular rather than mean-motion resonances. A double coordinate rotation—tilting the submoon relative to the moon, then the moon relative to the planet—is needed to initialize the inclined hierarchical 4-body problem correctly.","core_discovery":"On its own terms, the paper establishes that the published parameter estimates for the exomoon candidate Kepler-1625b-i are dynamically viable: most of the observationally allowed semimajor-axis and inclination range falls in regions of low maximum eccentricity in 100,000-year N-body integrations, with instability only where the von Zeipel-Lidov-Kozai mechanism excites eccentricity at inclinations above roughly 40 degrees. The same holds even more strongly for Kepler-1708b-i, whose proposed orbit sits so close to its planet that maximum eccentricity stays below about 0.005 for inclinations under 40 degrees. Extending to a fourth body, the paper finds that a submoon can remain bound around Kepler-1625b-i in portions of the 20-33% moon-Hill-radius range for inclinations below about 40 degrees, and it identifies the curved ridges of elevated eccentricity in that stability map as secular resonances, confirmed by a 13:2 ratio of apsidal precession rates and by secular perturbation theory reproducing the amplitude of the submoon's inclination oscillation.","pith_inferences":["Because the integrations are purely Newtonian, the stable submoon regions are upper limits: tidal dissipation, as the paper itself notes citing Kollmeier & Raymond (2018), would likely remove submoons on shorter timescales, so a submoon detection would require either weak tidal dissipation or a young system.","The 13:2 apsidal precession resonance is a high-order secular resonance; extending the maps to longer integration times (10^6 years) would test whether the resonant ridges eventually drive the submoon to ejection, and whether other precession ratios (e.g., 11:2, 15:2) produce similar ridges.","The coordinate-rotation fix for inclined hierarchical initial conditions suggests that earlier coplanar-only stability studies may have missed inclination-driven resonances of this kind; future exomoon stability surveys should scan inclination as a matter of course."],"forward_implications":["Kepler-1625b-i remains dynamically viable: most of its observed parameter range falls in low-maximum-eccentricity, bound regions, with only high-inclination parts susceptible to von Zeipel-Lidov-Kozai excitation.","Kepler-1708b-i is even more robust: its proposed semimajor axis sits at 6-11% of the planet's Hill radius, where maximum eccentricity stays below about 0.005 for inclinations under roughly 40 degrees.","Stable submoon orbits exist around Kepler-1625b-i in the 20-33% moon-Hill-radius range for inclinations below about 40 degrees, so a moon of a moon is dynamically possible.","The elevated-eccentricity ridges in the submoon stability map are secular resonances, identified by a 13:2 ratio of apsidal precession rates and confirmed by secular perturbation theory matching the inclination amplitude.","The same framework applied to HD 23079 (which has no known exomoon) predicts a stable region similar to Kepler-1625's, which can inform future exomoon searches there."],"supporting_citations":[{"why":"Supplies the Kepler-1625 exomoon candidate parameters (semimajor axis, mass, inclination range) that define the 3-body initial conditions and the stability comparison region.","marker":"Teachey & Kipping 2018"},{"why":"Supplies the Kepler-1708 exomoon candidate parameters used for the 3-body stability maps and the observed-range overlay.","marker":"Kipping et al. 2022"},{"why":"Establishes the 20-33% Hill-radius submoon stability range and the revised prograde satellite stability limit, and provides the prior coplanar submoon stability study this work extends to inclined orbits.","marker":"Rosario-Franco et al. 2020"},{"why":"Provides the von Zeipel-Lidov-Kozai secular approximation and maximum-eccentricity formula used to design the inclination grid and to compare theoretical contours with simulation.","marker":"Naoz 2016"},{"why":"Supplies the secular perturbation theory framework and the definition of secular resonances used to interpret the 13:2 precession ratio and the ridge structures.","marker":"Murray & Dermott 2000"},{"why":"Provides tidal migration calculations for Kepler-1625 submoons, cited as the reason the gravity-only stable regions may be tidally removed.","marker":"Kollmeier & Raymond 2018"},{"why":"Supplies the N-body integrator software used for all stability maps.","marker":"Rein & Liu 2012"}],"fun_headline_variants":["Exomoon orbits hold for 100,000 years in simulations","Submoons survive around Kepler-1625's moon at low tilt","High inclination triggers exomoon instability via secular resonance","Kepler-1625b-i passes stability test, submoons possible"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that 100,000 years of gravity-only motion is a fair test of stability: the simulations omit tidal dissipation, which the paper's own cited sources say would drag submoons inward and preferentially remove them, and they assume the contested exomoon detections with their adopted masses and orbits.","fun_headline_variants_meta":{"raw":{"variants":["Exomoon orbits hold for 100,000 years in simulations","Submoons survive around Kepler-1625's moon at low tilt","High inclination triggers exomoon instability via secular resonance","Kepler-1625b-i passes stability test, submoons possible"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000289,"raw_usage":{"total_tokens":1733,"prompt_tokens":1022,"completion_tokens":711,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":636}},"tokens_in":638,"tokens_out":711,"duration_ms":7312,"temperature":1.0,"reasoning_tokens":636,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T17:20:26.078085+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the 4-body stability maps with tidal dissipation included, using tidal quality factors comparable to those in Kollmeier & Raymond (2018), and check whether any submoon in the 20-33% moon-Hill-radius range remains bound for 100,000 years; if none do, the claimed stable regions are not physically real.","supporting_citations":[{"cited_title":"D., Dermott S","cited_arxiv_id":null,"evidence_quote":"Supplies the secular perturbation theory framework and the definition of secular resonances used to interpret the 13:2 precession ratio and the ridge structures."},{"cited_title":"A., Raymond S","cited_arxiv_id":null,"evidence_quote":"Provides tidal migration calculations for Kepler-1625 submoons, cited as the reason the gravity-only stable regions may be tidally removed."}],"review_version":1}