{"id":"8a11b945-8a6a-496b-bbef-ed9212e8e15e","arxiv_id":"2507.11594","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"N-body simulations with tidal spin evolution show that a Moon-like satellite around K2-18b is ejected within about 10 Myr, so the 3 Gyr-old system is unlikely to host an exomoon today.","lead":"This paper simulates whether a Moon-like satellite could survive around K2-18b, a possibly habitable planet orbiting a red dwarf. Tides push any moon outward so quickly that it leaves the planet's gravitational reach within about 10 million years, far too fast to still exist in the 3 billion-year-old system.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fixed 5 h initial spin (Table 1) forces outward tidal migration; the conclusion is explicitly conditional on this assumption, and no alternative spins are tested.","rationale":"The load-bearing concern is the unvaried initial spin. The paper's mechanism is outward tidal migration driven by a 5 h spin; the Moon's initial orbital period is ~36 h, so the sign of the tidal torque is determined by this single parameter. The authors themselves limit the conclusion in Section 4 to a rapid rotator, but the abstract and title state a stronger result. Because the spin of K2-18b is not observable and formation models allow a range of spins, the central claim is conditional on an untested assumption. Other issues (single Moon mass, neglect of K2-18c, stability criterion from the same group) are secondary: they would change lifetimes by factors, not change the sign of migration. The proposed spin grid directly tests whether plausible initial conditions can yield long-lived moons; this is the minimal experiment that would settle the robustness of the headline. The reader's CONDITIONAL verdict correctly reflects that the paper supports its claim only for the assumed scenario; our stress-test does not move that verdict.","tokens_in":7797,"tokens_out":14337,"duration_ms":174093,"concrete_test":"For a representative case (e.g., M_p = 8.63 M_E, e_p = 0.12, τ_p = 10 s, Neptune-like k2 and C), rerun the same rebound/reboundx setup with initial spin periods of 2, 5, 10, 20, 30, 36, 40, 50, and 100 h, plus one retrograde case, and record the maximum Moon lifetime. If any spin other than 5 h yields a lifetime exceeding 3 Gyr (or even exceeding 100 Myr), the paper's general claim that 'any moons would be extremely unlikely' fails for a plausible initial condition. The 36 h case, where the spin is synchronous with the initial lunar orbital period, is the critical control.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The planet's initial spin period is fixed at 5 h (Table 1), far shorter than the Moon's initial orbital period (~36 h at 3R_Roche). This makes the planet's tidal bulge lead the Moon, forcing a_m to increase and the Moon to cross the stability limit (Eq. 3) within ~10 Myr. If K2-18b instead formed with a slower spin (near or below the Moon's orbital frequency) or a retrograde spin, the tidal torque would pull the Moon inward or allow a near-synchronous state, and the lifetime could be much longer. Section 4 concedes that the conclusion assumes a rapid rotator, but the abstract and title present the result without this caveat. Since the actual initial spin of K2-18b is unconstrained, this single untested parameter controls the sign of tidal migration and therefore the headline conclusion. A 'migration reversal' of the type studied by Sasaki et al. (2012) is dismissed rather than simulated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates whether K2-18b could host a detectable exomoon by integrating star–planet–moon systems with rebound and reboundx's tides_spin module under the constant time-lag tidal model. The authors vary planetary mass (271 values), initial planetary eccentricity (three values), planetary tidal time-lag (three values), and two interior parameter sets (Earth-like and Neptune-like k2 and C). Each simulation places a Luna-mass moon at 3 R_Roche and defines instability as crossing the a_crit stability boundary of Eq. (3). Across 2,439 simulations, the maximum moon lifetime is found to be about 9.1 Myr, with an extrapolated circular-orbit limit of roughly 22–24 Myr, far shorter than the ~3 Gyr system age. The paper concludes that observable moons are unlikely around K2-18b and questions moon-based habitability scenarios for short-period M-dwarf planets in general.","tokens_in":7908,"tokens_out":4845,"duration_ms":62324,"significance":"If the result holds, it provides a concrete dynamical constraint on exomoon survivability in a high-profile habitable-zone M-dwarf system and offers a template for using N-body tidal simulations to filter exomoon survey targets. The main strengths are the broad parameter sweep (2,439 integrations, three eccentricities, three time-lags, two interior models), the use of a well-documented tidal evolution module that includes spin evolution, and the fact that the lifetimes are outputs of direct integration rather than assumed inputs. The principal weakness is that the fixed 5-hour initial spin period controls the sign of tidal migration and is not varied, so the headline conclusion is conditional on an unconstrained assumption. The absence of public code or data and the neglect of the known companion K2-18c further limit the strength of the broadest generalizations.","major_comments":[{"comment":"The fixed initial planetary spin period of 5 h (Table 1) places the planet in a super-synchronous state relative to the moon's initial orbit, which forces the tidal torque to push the moon outward. This single assumption controls the sign of the tidal migration that drives the headline result. If K2-18b instead began with a slower, near-synchronous, or retrograde spin, the tidal torque could draw the moon inward, keep it near a synchronous state, or otherwise allow substantially longer survival; the migration-reversal mechanism cited from Sasaki et al. (2012) is dismissed rather than simulated. Section 4 explicitly concedes that the conclusion assumes a rapid rotator, but the abstract and title present the result without this caveat. I request additional simulations with a range of initial spin periods (e.g., near the synchronous period, longer periods, and retrograde spin), or at least an analytic sign criterion for the tidal migration, and a revised abstract and title that state the conditional nature of the claim.","section":"Section 2, Table 1; Section 4"},{"comment":"The generalization that the result \"cast[s] doubt on moon-based habitability scenarios for short-period M-dwarf planets in general\" goes beyond what the simulations establish. The study varies M_p, e_p, tau_p, k2, and C, but it holds the moon mass at a single lunar mass, fixes the initial moon semi-major axis at 3 R_Roche, fixes the initial planetary spin at 5 h, and omits the known companion K2-18c. Any of these choices could alter the tidal migration rate or the effective stability boundary. The broad statement should either be removed or supported by simulations that explore at least the moon-mass and initial-spin dimensions.","section":"Abstract; Section 4"}],"minor_comments":[{"comment":"The column header for a_m lists units of au, but the entry is \"3 R_Roche\"; please specify the corresponding au value or change the header so the units are consistent.","section":"Table 1"},{"comment":"The MCMC linear fits use only three eccentricity values per time-lag value, and the error bars are described as \"not statistically significant.\" The slopes and y-intercepts in Table 2 should be presented as descriptive linear interpolations rather than robust statistical inferences, and the circular-orbit extrapolation should be labeled as such.","section":"Section 3, Fig. 3, Table 2"},{"comment":"The simulations treat the star and moon as point masses and the planet with a constant radius of 2.61 R_E. The tidal response of the star is not modeled; a brief justification or citation for neglecting stellar tides on the 10 Myr integration timescale would improve the completeness of the methods.","section":"Section 2"},{"comment":"The statement that \"the data underlying this article will be shared on reasonable request\" is insufficient for a study with 2,439 simulations; please deposit the rebound/reboundx scripts, initial-condition files, and summary outputs in a public repository to enable reproducibility.","section":"Data Availability"},{"comment":"K2-18c is mentioned as a companion but is not included in the N-body model. A sentence quantifying the expected perturbation from K2-18c (for example, the ratio of its semi-major axis to K2-18b's Hill radius) would make the Payne et al. (2013) assumption of negligible planetary perturbations more transparent.","section":"Section 1"}],"recommendation":"major_revision","confidential_remarks":"The central N-body calculation is sound and the parameter sweep is a genuine strength, but the fixed 5-hour spin assumption is load-bearing for the headline conclusion. The authors should be encouraged to add spin-variation simulations; if that is not feasible, the abstract and title must be qualified to state that the result applies under the assumption of a rapidly rotating initial planet. The paper fits the journal's scope and is publishable after these revisions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read. The headline result is real but narrower than the title. The paper simulates, with rebound/reboundx tides_spin, a Moon-mass satellite around K2-18b across a wide grid of planet mass, eccentricity, time-lag, and two interior models, and finds maximum lifetimes of 10^4-10^7 yr, under 10 Myr even at the lowest dissipation. That is a new, concrete, K2-18b-specific result. The parameter sweep is broad, the integration scheme is standard and clearly documented, and the MCMC fit to e_p=0 is a reasonable way to extrapolate. The authors also correctly note that varying the moon mass changes lifetimes by at most a few.\n\nThe soft spot is exactly what the stress-test flags: the initial planetary spin is fixed at 5 h, which forces the tidal bulge to lead the Moon and pushes the Moon outward. That is a physically plausible choice for a forming planet, but it is not tested against slower or retrograde spins. The authors do concede in Section 4 that the conclusion assumes a rapid rotator, but the abstract and title do not. The consequence is that the paper does not actually show that K2-18b lacks moons; it shows that K2-18b lacks moons if it formed spinning fast. That makes the 'extremely unlikely' language in the abstract too strong, and the sweep to 'short-period M-dwarf planets in general' is a step beyond the evidence.\n\nAlso minor: only one moon mass (though they argue the mass dependence is weak), no perturbation from K2-18c, and no code or data release beyond 'available on request'. These are minor next to the spin issue.\n\nThe paper is worth a serious referee. The methods are solid, the result is useful for exomoon survey targeting, and the caveat is fixable by adding a spin parameter or at least an explicit discussion of why 5 h is the only relevant case. I would send to peer review, but I would require the authors to either test the spin sensitivity or soften the abstract to match the conditional nature of the claim. Cite it if you work on exomoon lifetimes; it's a clean, well-documented numerical study, just not the final word on K2-18b.","headline":"A solid, well-documented simulation that shows K2-18b likely has no moons if it formed as a rapid rotator, but the abstract oversells the result by dropping that caveat.","tokens_in":8587,"tokens_out":3098,"would_cite":true,"duration_ms":36207,"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":"Tidal forces would push any moon around K2-18b past its stability limit within ~10 million years, far under the system's 3-billion-year age, so the habitable-zone planet is moonless today.","keywords":["exomoons","K2-18b","tidal evolution","habitable zone","M-dwarf planets","N-body simulations","reboundx","astrobiology"],"falsifier":"Two concrete tests would decide the matter: a transit-timing or transit-duration search that detects a moon signal around K2-18b would directly contradict the null prediction, and a re-run of the same simulations with a slow or retrograde initial planetary spin that yields lifetimes above 3 Gyr would show the conclusion depends on the spin assumption rather than on the system itself.","tokens_in":7473,"feed_emoji":"🌙","tokens_out":12344,"duration_ms":118915,"temperature":0.7,"pith_summary":"K2-18b sits in its star's habitable zone and has shown tentative signs of biosignature gases in its atmosphere, making it a natural test case for whether a large moon could add habitable surface area. This paper asks whether a moon as massive as Earth's own Moon could actually remain in orbit there. Across 2,439 N-body simulations that include tidal spin evolution, the authors find that tides push any moon outward past the orbital stability limit within about 10 million years, roughly 300 times shorter than the system's 3-billion-year age. They conclude that K2-18b cannot host an observable moon today, and that moon-based habitability for short-period M-dwarf planets is doubtful because moons are removed before they can stabilize climate.","feed_headline":"Moons around K2-18b are doomed within 10 million years","feed_subtitle":"Tidal drag pushes any moon past the stability limit about 300 times faster than the star ages, leaving none today.","key_machinery":"The argument rides on two mechanisms working together. The first is the constant time-lag model of tides, implemented as the tides_spin module: the planet's tidal bulge lags the moon by a fixed time $\\tau_p$, transferring angular momentum and steadily evolving the moon's semi-major axis and eccentricity. The second is the critical semi-major axis $a_{\\rm crit} = 0.4031(1 - 1.123\\,e_p)\\,R_H$, with $R_H$ the planet's Hill radius, which defines the stability boundary beyond which a moon counts as lost; each simulation is scored by the moon's maximum lifetime $t_{\\max}$, the time to cross $a_{\\rm crit}$. The assumed rapid initial spin of the planet (5-hour period) sets the sign of the tidal torque so that the moon always migrates outward, and the relatively small Hill radius of this close-in planet makes the outward journey short.","core_discovery":"The paper's central claim is that exomoons are extremely unlikely to survive around K2-18b because planetary tides drive a moon's orbit outward on a timescale far shorter than the system age. In the full grid of simulations, the longest moon lifetime is about 9.1 Myr, achieved only with the most favorable parameters tested: a Neptune-like tidal response (Love number $k_2 = 0.120$), the smallest tidal time-lag $\\tau_p = 10$ s, and the lowest observed planetary eccentricity $e_p = 0.12$. Faster dissipation (larger $\\tau_p$) and higher eccentricities shorten lifetimes, and even the circular-orbit limit only extends the maximum to roughly 22–24 Myr. Since the K2-18 system is about 3 Gyr old, the authors conclude that no moon can be observed there now, and they extend the result to short-period M-dwarf planets generally: rapid tidal-driven migration removes moons before they can play any climate-stabilizing role, so moon-based habitability scenarios for such systems are doubtful.","pith_inferences":["The assumed 5-hour initial spin is load-bearing: if K2-18b instead formed with a slow or retrograde spin, the tidal torque would reverse and could pull a moon inward, letting it survive far beyond 10 Myr. The paper states this assumption but does not vary it, so the result is a statement about rapidly rotating young planets.","The ~10 Myr removal timescale is so short that even a moon formed in situ after the planet's birth would have to appear within an extremely narrow window to be observable now; this strengthens the null prediction but also means the simulations assume the moon existed at all.","The same simulation machinery could be turned into a quick screening catalog for other habitable-zone M-dwarf planets: any planet with a short period and small Hill radius is a poor exomoon target, no matter how promising its atmosphere.","If future JWST-style observations confirm biosignature gases while independent data rule out a moon, the habitability discussion for K2-18b would have to lean entirely on the planet itself rather than on moon-supported environments."],"forward_implications":["No moon is observable around K2-18b today: the ~10 Myr migration timescale is shorter than the 3 Gyr system age by a factor of about 300.","Because even the circular-planet limit gives lifetimes of only ~22–24 Myr, the null prediction is robust across the observed eccentricity range.","The result acts as a predictive filter for exomoon surveys, flagging K2-18b as a false-positive-prone target and helping prioritize other systems.","Moon-based habitability scenarios for short-period M-dwarf planets are doubtful in general, since tides remove moons before they can stabilize planetary obliquity or host subsurface oceans.","Changing the moon's mass changes the lifetime by at most a factor of a few, so the conclusion does not hinge on the assumed Luna-mass satellite."],"supporting_citations":[{"why":"Supplies the constant time-lag tidal model that the simulations integrate, defining how the lagging tidal bulge exchanges angular momentum with the moon.","marker":"Eggleton et al. (1998)"},{"why":"Provides the tides_spin module in reboundx that evolves the moon's semi-major axis and eccentricity together with the planetary spin.","marker":"Lu et al. (2023)"},{"why":"Defines the critical semi-major axis used as the stability boundary that determines each simulated moon's lifetime.","marker":"Rosario-Franco et al. (2020)"},{"why":"Supplies the observed planetary mass range and the three eccentricity values (0.12, 0.20, 0.28) that grid the simulations.","marker":"Sarkis et al. (2018)"},{"why":"Provides the planet's fixed radius of 2.61 Earth radii and the mass prior that set the physical scale tides act on.","marker":"Benneke et al. (2019)"},{"why":"Establishes the tidal-evolution equations showing that stronger dissipation speeds up semi-major axis migration, used to interpret the time-lag dependence.","marker":"Hut (1981)"},{"why":"Supports the assumed 5-hour initial spin period of the planet, which sets the outward direction of tidal migration.","marker":"Kokubo & Genda (2010)"},{"why":"Identifies inward migration and spin effects as a mechanism that could lengthen moon lifetimes, which the paper argues cannot occur under rapid-rotation assumptions.","marker":"Sasaki et al. (2012)"}],"fun_headline_variants":["Moons around K2-18b die in under 10 million years","K2-18b's moons: a 10-million-year death sentence","Tidal forces wipe out K2-18b moons in 10 million years","Exomoons can't survive K2-18b: tidal migration ends in 10 Myr"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole conclusion assumes K2-18b began as a rapid rotator with a 5-hour spin period, so tides always push a moon outward; if the planet instead started with a slow or retrograde spin, the tides could pull a moon inward and let it survive much longer.","fun_headline_variants_meta":{"raw":{"variants":["Moons around K2-18b die in under 10 million years","K2-18b's moons: a 10-million-year death sentence","Tidal forces wipe out K2-18b moons in 10 million years","Exomoons can't survive K2-18b: tidal migration ends in 10 Myr"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000644,"raw_usage":{"total_tokens":2963,"prompt_tokens":948,"completion_tokens":2015,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":1935}},"tokens_in":564,"tokens_out":2015,"duration_ms":16524,"temperature":1.0,"reasoning_tokens":1935,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:08:18.472288+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Two concrete tests would decide the matter: a transit-timing or transit-duration search that detects a moon signal around K2-18b would directly contradict the null prediction, and a re-run of the same simulations with a slow or retrograde initial planetary spin that yields lifetimes above 3 Gyr would show the conclusion depends on the spin assumption rather than on the system itself.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the tidal-evolution equations showing that stronger dissipation speeds up semi-major axis migration, used to interpret the time-lag dependence."},{"cited_title":"W., O'Brien D","cited_arxiv_id":null,"evidence_quote":"Identifies inward migration and spin effects as a mechanism that could lengthen moon lifetimes, which the paper argues cannot occur under rapid-rotation assumptions."}],"review_version":1}