{"id":"70bcb0a4-545a-4eb1-a24c-dcb046576855","arxiv_id":"2608.13243","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Numerical simulations show that the habitable zones of 36 Oph and 70 Oph are dynamically stable, while gamma Leo's are not.","lead":"This paper simulates whether three nearby binary star systems could host habitable planets. It finds that 36 Oph and 70 Oph are promising targets, while gamma Leo is dynamically hostile.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"1 Myr integrations do not establish the 'permanently habitable' claim over stellar lifetimes; a convergence test with longer runs is needed.","rationale":"The reader's conditional verdict is appropriate, but the weakest assumption they identified is not the one most likely to break the central claim. The 36 Oph orbital solution is indeed unpublished and conflicts with a published low-eccentricity solution, and a sensitivity test would be valuable. However, the published alternative orbit has e ≈ 0.34 with a similar or larger period, which gives a much larger periastron distance; that configuration is dynamically more benign for the HZ, so it would probably strengthen rather than overturn the 'promising target' conclusion. The more load-bearing gap is the mismatch between 1 Myr simulations and the manuscript's 'permanently habitable' and 'to the present day' language. For 36 Oph in particular, only ~2,000 binary periods are integrated, no analytic stability check is provided for e > 0.8, and no convergence test is reported. This directly threatens the positive claim for 36 Oph, which is the system with the most uncertain and extreme input orbit. The gamma Leo negative result has an internal Kepler inconsistency (P = 554 yr, a = 67 AU, M_tot = 3.21 Msun), but recomputing a self-consistent orbit still leaves the binary periastron inside or near the red-giant HZs, so that conclusion is likely robust. The paper otherwise uses a standard N-body tool (REBOUND/IAS15) and a standard flux-based HZ prescription, and the 70 Oph result is well supported by the Holman & Wiegert-style estimate. Since the timescale issue is addressable by longer integrations, a conditional acceptance remains the right call.","tokens_in":17825,"tokens_out":23705,"duration_ms":256828,"concrete_test":"Extend the 36 Oph simulations (both 0-degree and 45-degree inclination cases, around both stars) to at least 10 Myr, and preferably 100 Myr, using the same posterior sampling. Compare the cumulative ejection and uninhabitable fractions, and the surviving permanently habitable boundaries, against the 1 Myr results. If the fractions plateau and the PHZ boundaries converge, the timescale concern is resolved; if fractions grow or the PHZ shrinks materially, the 'can support permanently habitable planets' claim must be weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that 36 Oph and 70 Oph can support permanently habitable planets rests on 1 Myr simulations (§3.1). This is too short to support the stronger language in §4.1.1 that planets 'could very likely have survived the early development of the system and would continue to be habitable to the present day.' For 36 Oph, the binary period is ~500 yr, so the integrations cover only ~2,000 periastron passages. The paper does not apply an analytic stability criterion to 36 Oph: the Holman & Wiegert-style check in §3.2.2 is used only for 70 Oph, and the authors note the formula is restricted to e < 0.8. High-e binaries can exhibit secular chaos and slow diffusion, so a test particle that survives 2,000 periastron passages could still be ejected on longer timescales. The PHZ in this paper is defined as survival for the simulation duration, so the 'permanently habitable' conclusion is only as strong as the 1 Myr runs. The reader's alternative-orbit concern is legitimate, but the published low-e alternative (e ≈ 0.34) would likely widen the stable region rather than eliminate it, making the timescale issue the more direct threat to the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents N-body simulations of massless test particles placed in the habitable zones (HZs) of the three nearby stellar binaries 36 Oph, 70 Oph, and gamma Leo. Using the REBOUND IAS15 integrator, the authors inject 100 particles per star, linearly spaced in semi-major axis across the HZ, at mutual inclinations of 0 degrees (coplanar) and 45 degrees relative to the binary plane, and integrate for 1 Myr. Binary orbital parameters are drawn from adopted solutions: Giovinazzi et al. (2026, in prep) for 36 Oph, Li et al. (2026) for 70 Oph, and Romanenko & Kiselev (2014) for gamma Leo. Particles are classified as dynamically unstable, uninhabitable (instantaneous flux outside the HZ at any time), or habitable (always within the HZ). The authors find that 36 Oph and 70 Oph can support habitable planets in their HZs, with coplanar orbits more favorable than 45-degree-inclined orbits, while gamma Leo's HZs are dynamically unstable, with all injected particles ejected. They recommend 36 Oph and 70 Oph as targets for the Habitable Worlds Observatory and SHERA, and argue that gamma Leo is unlikely to host habitable planets. The paper also addresses the stability of the known planet gamma Leo Ab and the unconfirmed ~1340-day candidate around gamma Leo A.","tokens_in":18016,"tokens_out":14663,"duration_ms":141410,"significance":"If the conclusions hold, the paper provides a practical framework for prioritizing nearby binaries for future habitable-planet searches, and it identifies two concrete systems (36 Oph and 70 Oph) with apparently stable HZs. The simulation setup is conventional and the paper's logic is largely transparent: it uses publicly available integration software (REBOUND/IAS15), a standard flux-based HZ definition, and a consistency check against the Holman & Wiegert (1999) criterion for 70 Oph. The paper does not fit any free parameters; its outcomes are genuine predictions from literature-based inputs, which is a strength. However, the central 'permanently habitable' claim is stronger than the 1-Myr integrations can support, particularly for 36 Oph, and the 36 Oph result depends on an unpublished orbital solution whose eccentricity is dramatically different from a published alternative. These issues are fixable with additional robustness tests and more qualified language, but they presently overstate the certainty of the conclusions.","major_comments":[{"comment":"The 1 Myr integration time does not support the 'permanently habitable' language in the abstract and in §4.1.1 ('would continue to be habitable to the present day'). The PHZ is explicitly defined in §3.2 as survival for the simulation duration, not for the stellar lifetime. For 36 Oph, the binary period is ~500 yr (Table 1), so the runs include only ~2000 periastron passages, and no analytic stability criterion is applied because the Holman & Wiegert (1999) formula used for 70 Oph (§3.2.2, Eq. 4) is restricted to e<0.8. High-e binaries can exhibit slow secular chaos on timescales beyond the 1 Myr window; the authors should either run longer integrations (e.g., 10-100 Myr), apply a high-e stability boundary, or explicitly qualify the conclusions as 'stable over 1 Myr' rather than 'permanently habitable.'","section":"§3.1, §3.2, §4.1.1"},{"comment":"The 36 Oph result rests on the orbit of Giovinazzi et al. (2026, in prep), which is not publicly available, and its eccentricity e=0.90 differs dramatically from the published Izmailov & Khovritchev (2025) solution (e=0.34) discussed in Appendix A.1. Section 4.1.1 asserts that future observations will 'refine, rather than lead to a dramatically different posterior,' but this is an unsupported assertion, not a robustness test. The authors should run the 36 Oph suite with the alternative published orbit (or with a sample that includes its range), or explicitly state that the adopted high-e solution is the conservative choice because lower eccentricity should broaden the stable region, and base the robustness claim on that directional argument.","section":"§2.2, §4.1.1, Appendix A.1"},{"comment":"The binary orbital initial conditions for 36 Oph and 70 Oph are not fully specified: Tables 1 and 2 list P, a, and e, but omit the inclination i, argument of periastron ω, and longitude of ascending node Ω. Since the test-particle inclinations are defined relative to the binary plane (§3.1), and since the geometry of close encounters in the 45-degree runs can depend on ω at octupole order (important for e=0.9), the simulations cannot be reproduced from the information given. The full adopted orbital elements (or explicit fixed values) should be reported for both systems.","section":"§3.1.1, §3.1.2, Tables 1-2"},{"comment":"The binary orbital parameters are sampled from independent normal distributions (§3.1.1-§3.1.3), ignoring covariances between a and e. For astrometric orbits these parameters are strongly anti-correlated through the periastron distance q=a(1-e); independent sampling can generate unphysical orbital pairs and overstate the spread in q, which directly controls the perturbation strength. The authors already sample posteriors for the orbit plots (Figures 2, 4, 6); they should use the joint posterior draws, or at least demonstrate that the results are insensitive to the ignored correlations.","section":"§3.1.1-§3.1.3"}],"minor_comments":[{"comment":"The paper states in §3.1.3 that the ~1340-day candidate is neglected, but §3.2.3 reports a suite that includes it; please clarify which statement applies to the primary analysis and describe the candidate's initial conditions in the secondary suite.","section":"§3.1.3 and §3.2.3"},{"comment":"The abstract's statistic 'become uninhabitable ... only 1.5% - 1.8% of the time' is ambiguous; the text in §3.2.1 reports a rate of planets becoming uninhabitable, so the abstract should say 'of the planets' or define the temporal-average interpretation explicitly.","section":"Abstract and §3.2"},{"comment":"The caption's phrase 'apparent overlap ... largely a consequence of visualization' is contradicted by §4.1.3, which states that the binary orbits 'physically intersect with the habitable zone'; the caption should be corrected for gamma Leo.","section":"Figure 5 caption and §4.1.3"},{"comment":"The statement that the secondary star's insolation is negligible for binaries 'separated by several hundred AU' does not literally apply to 70 Oph (a=23 AU); the justification should be quantified for each system (flux ratio at the HZ) rather than citing the several-hundred-AU distance criterion.","section":"§2.1 and §3.2"},{"comment":"The paper would benefit from a data availability statement and a repository containing the simulation scripts and the adopted posterior samples, particularly for the two in-preparation orbital solutions.","section":"General"},{"comment":"Table 3 lists the mass of gamma Leo Ab as sourced to Takeda (2023), while Section 2.4 states that mass estimates are adopted from Han et al. (2010); please reconcile the source attribution.","section":"Table 3"},{"comment":"For reproducibility, please state the IAS15 accuracy parameter (e.g., epsilon) and the output cadence used for the flux classification.","section":"§3.1"},{"comment":"There is a typographical double period after '0.300-0.308 AU for 70 Oph B'; please fix the punctuation.","section":"§3.2.2"}],"recommendation":"major_revision","confidential_remarks":"The paper's primary new result for 36 Oph relies on an orbit from a paper in preparation by co-authors (Giovinazzi et al. 2026). This is not necessarily a problem, but the editor should require that the orbit's posterior samples be made available, or that a robustness test with the published alternative orbit be added, before publication. The 'permanently habitable' overclaim based on 1 Myr integrations is the main scientific issue and is raised as a major comment; it applies to the abstract and conclusions and needs to be addressed before the paper can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one. The paper does something concrete: it takes three nearby, bright binaries and runs standard IAS15/REBOUND suites to ask which ones could still host a habitable-zone planet after a gigayear's worth of perturbations. The headline result — 36 Oph and 70 Oph are viable targets, gamma Leo is not — is exactly the kind of input HWO/SHERA target selection needs. The simulations are conventional but careful: the HZ flux includes eccentricity effects, the PHZ definition is explicit, and the known giant planet around gamma Leo A is included rather than ignored. No one is fitting to the conclusion; the orbital elements are external inputs, and the outputs are genuine predictions. The 70 Oph result is the most solid: the binary's full orbit is observed, and the Holman-Wiegert/Quarles stability check lines up with the numerics.\n\nWhere it frays: the word 'permanently'. The PHZ is defined as survival for the 1 Myr simulation, which covers many Kozai-Lidov cycles but only about two thousand periastron passages for 36 Oph's 500-year orbit. High-eccentricity secular chaos can eject test particles on much longer timescales, and the paper does not apply any analytic stability criterion to 36 Oph because Holman-Wiegert stops at e=0.8. So the claim in §4.1.1 that planets 'could very likely have survived... and would continue to be habitable to the present day' is stronger than the runs support. A 10 Myr convergence suite, or at least a stated fallback for slow diffusion, is needed before that language can stand.\n\nThe bigger soft spot is 36 Oph's input orbit. The paper adopts an in-prep solution with e=0.90 and explicitly notes the published Izmailov & Khovritchev solution has e=0.34. That is a factor-of-three difference in periastron distance, which directly controls secular perturbation strength. The authors assert future observations will refine rather than overturn their posterior, but they do not run the obvious test with the published low-e solution. It is quite possible a low-e run would widen the stable region and keep the target recommendation intact — which would make the test easy and persuasive. As written, the 36 Oph result is hostage to an unreleased orbit.\n\nMinor points: the gamma Leo coplanar case is excluded with a reasonable geometric justification, and the 1340-day candidate is actually included in one suite and ejected, so the dismissal is a simulation result rather than an arbitrary omission. No code or data are released, which is a shame but not a blocker at this stage.\n\nFor a reader: this is a subfield paper, not a methods breakthrough. Its value is that it gives mission planners a concrete, reproducible argument for where to look. I'd send it to peer review with a request for longer integrations and the 36 Oph sensitivity run, and I'd expect a solid paper after that revision.","headline":"Useful target-list study with standard methods, but the 'permanently habitable' language outruns the 1 Myr integrations and the 36 Oph conclusion rests on an unpublished high-e orbit.","tokens_in":18611,"tokens_out":3305,"would_cite":true,"duration_ms":33626,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Simulations show 36 Oph and 70 Oph can host permanently habitable planets, while gamma Leo cannot.","keywords":["binary stars","habitable zone","exoplanet dynamics","N-body simulations","permanently habitable zone","Kozai-Lidov cycles","nearby stellar multiples","target selection"],"falsifier":"Measure the 36 Oph binary orbit with the next Gaia data release; if its eccentricity turns out to be about 0.34 rather than 0.90, rerun the paper's 1 Myr test-particle simulations. If habitable-zone planets then become frequently ejected or uninhabitable, the promising-target conclusion for 36 Oph collapses, while if they remain stable the claim is robust to the orbit ambiguity. Separately, a confirmed long-lived planet inside gamma Leo's habitable zone would directly contradict the paper's exclusion of that system.","tokens_in":17616,"feed_emoji":"🪐","tokens_out":7479,"duration_ms":63658,"temperature":0.7,"pith_summary":"This paper asks whether three bright, nearby binary systems could host planets that stay habitable over long timescales. Using 1-million-year N-body simulations with test planets spread across each star's habitable zone, it finds that 36 Oph and 70 Oph can keep most such planets both dynamically stable and within liquid-water flux bounds, with low ejection rates even in the least favorable cases. The same simulations find that the high-eccentricity red-giant binary gamma Leo ejects every test planet in its habitable zone, so that system is unlikely to host any habitable world. If the conclusions hold, 36 Oph and 70 Oph become priority targets for the next generation of habitable-planet searches.","feed_headline":"Two nearby binaries can host habitable planets","feed_subtitle":"Simulations rank 36 Oph and 70 Oph as top targets while gamma Leo's eccentric orbit rules it out","key_machinery":"The load-bearing machinery is a suite of N-body integrations in which massless test particles are inserted into each star's habitable zone, computed from the single-star flux limits of 1.7 and 0.3 times Earth's insolation, and followed for 1 Myr with a high-accuracy integrator. The mechanism that decides the outcome is secular gravitational forcing by the binary companion: it drives oscillations in each planet's eccentricity, and when those oscillations push the instantaneous flux outside the habitable boundaries the planet is classified as uninhabitable, while collisions or unbound orbits are classified as dynamically unstable. The paper uses the Kozai-Lidov timescale, evaluated from the binary and planet orbital periods, to confirm that the 1 Myr integration covers many eccentricity cycles for the 45-degree misaligned cases.","core_discovery":"The paper's central claim is that the K-dwarf binaries 36 Oph and 70 Oph can support permanently habitable planets around both components, whereas the red-giant binary gamma Leo cannot. In the two favorable systems, planets placed in the habitable zone on orbits coplanar with the binary become uninhabitable only 1.5%-1.8% of the time, while planets started 45 degrees off the binary plane undergo larger eccentricity oscillations and become uninhabitable 4.8%-5.4% of the time; ejection rates remain low. For gamma Leo, the binary's eccentricity of 0.90 makes the stellar orbits cross the habitable zones, and every simulated test planet is ejected, eliminating the habitable zone around both stars under the adopted orbital solution. The paper therefore recommends 36 Oph and 70 Oph as targets for the Habitable Worlds Observatory and SHERA, and argues that the unconfirmed roughly 1340-day planet candidate around gamma Leo A is likely not a real planet.","pith_inferences":["The authors do not explore, but a lower true eccentricity for 36 Oph, as in the alternative published solution, would likely make its habitable zone even more stable because the binary's perturbing effect scales with eccentricity; the paper's conclusion would then survive a downward revision, though the quantitative uninhabitable percentages would change.","The paper's habitability criterion is conservative in one direction and optimistic in another: a planet that briefly leaves the flux limits might retain liquid water through climate inertia or oceans, while the test-particle treatment ignores planet-planet interactions and unseen companions that could destabilize real systems.","The same selection method could be extended to estimate how long a planet must remain in the habitable zone to be considered promising, using the simulation output to compute residence-time distributions rather than a binary habitable or uninhabitable flag.","If future astrometry finds a habitable-zone planet in gamma Leo, the adopted orbit would need revision; the paper's exclusion claim is tied to the current 0.90-eccentricity solution."],"forward_implications":["36 Oph and 70 Oph should be added to the shortlist of nearby systems for HWO and SHERA follow-up, with the inner parts of their habitable zones identified as the most survivable.","A planet discovered in the 70 Oph habitable zone by astrometry would most plausibly be near-coplanar with the binary, so measuring its orbit would constrain its inclination and hence its mass.","Gamma Leo can be deprioritized for habitable-planet searches, and the roughly 1340-day radial-velocity signal around gamma Leo A is more likely stellar activity than a planet.","The same simulation recipe can be applied to other nearby binaries to screen targets before committing expensive observing time."],"supporting_citations":[{"why":"Supplies the adopted 36 Oph binary orbit, including the high eccentricity of 0.90, and the stellar masses.","marker":"Giovinazzi et al. (2026, in prep)"},{"why":"Supplies the adopted 70 Oph binary orbit and precise stellar masses used as simulation input.","marker":"Li et al. (2026)"},{"why":"Provides the flux formula used to define habitable-zone boundaries and to judge habitability over time.","marker":"Bolmont et al. (2016)"},{"why":"Supplies the Kozai-Lidov timescale used to verify that the 1 Myr simulations cover many eccentricity cycles.","marker":"Antognini (2015)"},{"why":"Provides the empirical stability limit against which the 70 Oph dynamical stability results are checked.","marker":"Holman & Wiegert (1999)"},{"why":"Provides the gamma Leo binary orbital solution with eccentricity 0.90 that drives the simulation outcome.","marker":"Romanenko & Kiselev (2014)"},{"why":"Supplies the parameters for gamma Leo Ab and the unconfirmed 1340-day planet candidate.","marker":"Han et al. (2010)"},{"why":"Defines the permanently habitable zone criterion used to classify simulated planets.","marker":"Eggl et al. (2013)"},{"why":"Provides the N-body integration architecture used to run the simulations.","marker":"Rein & Liu (2012)"},{"why":"Justifies treating each star's habitable zone as single-star-like in these wide binaries.","marker":"Kaltenegger & Haghighipour (2013)"}],"fun_headline_variants":["Two nearby binaries can host stable habitable planets","36 Oph and 70 Oph could support permanently habitable worlds","gamma Leo's high-eccentricity orbit rules out habitable planets","Simulations find habitable zones in 36 Oph and 70 Oph survive","36 Oph and 70 Oph: stable habitats possible, gamma Leo not"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results for 36 Oph and 70 Oph rest on the adopted binary orbital solutions being accurate, and the 36 Oph solution in particular is a new, not-yet-published orbit with an eccentricity of about 0.90 that differs sharply from an alternative published value of about 0.34.","fun_headline_variants_meta":{"raw":{"variants":["Two nearby binaries can host stable habitable planets","36 Oph and 70 Oph could support permanently habitable worlds","gamma Leo's high-eccentricity orbit rules out habitable planets","Simulations find habitable zones in 36 Oph and 70 Oph survive","36 Oph and 70 Oph: stable habitats possible, gamma Leo not"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001039,"raw_usage":{"total_tokens":4426,"prompt_tokens":1057,"completion_tokens":3369,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":673,"completion_tokens_details":{"reasoning_tokens":3273}},"tokens_in":673,"tokens_out":3369,"duration_ms":23089,"temperature":1.0,"reasoning_tokens":3273,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:29:55.970527+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the 36 Oph binary orbit with the next Gaia data release; if its eccentricity turns out to be about 0.34 rather than 0.90, rerun the paper's 1 Myr test-particle simulations. If habitable-zone planets then become frequently ejected or uninhabitable, the promising-target conclusion for 36 Oph collapses, while if they remain stable the claim is robust to the orbit ambiguity. Separately, a confirmed long-lived planet inside gamma Leo's habitable zone would directly contradict the paper's exclusion of that system.","supporting_citations":[{"cited_title":"2026, in prep, AAS","cited_arxiv_id":null,"evidence_quote":"Supplies the adopted 36 Oph binary orbit, including the high eccentricity of 0.90, and the stellar masses."},{"cited_title":"2013, ApJ, 764, 130, doi: 10.1088/0004-637X/764/2/130","cited_arxiv_id":null,"evidence_quote":"Defines the permanently habitable zone criterion used to classify simulated planets."}],"review_version":1}