{"id":"45fc9ec4-3a85-416e-9a4c-b4459043408b","arxiv_id":"2505.20558","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A new CHZ2 metric, integrating the Bayesian likelihood of a 2 Gyr continuous habitable zone outside a coronagraph inner working angle, ranks the 164 HWO target stars.","lead":"This paper computes a Bayesian metric for how much of a star's habitable zone has stayed continuously habitable for 2 billion years, and applies it to the 164 stars on NASA's target list for the future Habitable Worlds Observatory. It ranks which stars offer the largest accessible zones where life could have had Earth-like time to alter an atmosphere.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"CHZ2 metric ranking is dominated by the assumed 2 Gyr CHZ definition and the IWA cutoff; the 3-4 Gyr peak for late-F/early-G stars is not demonstrated to be robust to plausible changes in these assumptions.","rationale":"The reader identified the 2 Gyr CHZ definition as the weakest assumption, and I agree that the biological timescale is a fundamental interpretational assumption. However, I see an even more immediate technical load-bearing concern: the paper's quantitative claims (3-4 Gyr peak, F3/6600 K cutoff) are not accompanied by any robustness analysis or uncertainty propagation for the CHZ2 metric. The paper itself notes the IWA is assumed and that the metric assumes all planets outside IWA are observable, and the age uncertainties are large. The comparison with SPORES validates masses and ages, but not the CHZ2 metric, which is the central output. Thus, the central claim holds only if the metric is insensitive to plausible variations in the CHZ duration, IWA, and input stellar properties. The paper is a useful demonstration of a method, and the code is open-source, but the headline scientific claims are not yet shown to be robust. A sensitivity analysis would settle this. I partially agree with the reader because they focused on the biological timescale (which is a conceptual assumption), while I focus on the metric's robustness to technical assumptions that are within the paper's control. Both are load-bearing, but the technical one is more directly testable and more central to the paper's quantitative claims.","tokens_in":17875,"tokens_out":1677,"duration_ms":13953,"concrete_test":"Recompute the CHZ2 metric for the 164 EMSL stars under three perturbations: (1) CHZ duration of 1 Gyr and 4 Gyr instead of 2 Gyr; (2) IWA of 60 mas and 100 mas instead of 83 mas; (3) resampling Teff, R*, [Fe/H] from the SPORES uncertainties to propagate errors into the metric. If the 3-4 Gyr peak and the F3/6600 K cutoff persist across all perturbations (or if the peak shifts but the qualitative ranking is preserved), the central claim is robust. If the peak moves or the cutoff disappears, the headline claim should be weakened to a model-dependent suggestion.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that the CHZ2 metric peaks at 3-4 Gyr for late-F/early-G dwarfs and drops for stars hotter than ~6600 K. This claim depends on the fixed-duration CHZ2 definition (Section 1) and on the metric's integration only over radii outside the IWA (Section 2.5). The paper does not propagate uncertainties from stellar properties (Teff, R*, [Fe/H]) into the CHZ2 metric values, and it does not test how the metric changes with the CHZ duration or with the IWA choice. The paper itself notes the CHZ2 definition is tied to the GOE timescale (2 Gyr) and the IWA is assumed (83 mas). If the CHZ duration is, e.g., 1 Gyr or 4 Gyr (relevant since the paper's own age uncertainties are >1 Gyr), or if the IWA is 60 mas or 100 mas, the ranking and the claimed peak could shift. The paper also does not report the metric uncertainties, so it is unclear whether the 3-4 Gyr peak and the F3/6600 K cutoff are significant or just noise. The comparison with SPORES is only for masses and ages, not for the CHZ2 metric itself, so the metric robustness is untested. This is a load-bearing concern because the headline quantitative claims (3-4 Gyr peak, F3/6600 K cutoff) are the main novel results and their scientific value depends on their stability.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a Bayesian method for computing the likelihood that a given orbital radius around a star lies in the 2 Gyr continuous habitable zone (CHZ2), then integrates this likelihood outside a fixed inner working angle to define a single 'CHZ2 metric' for HWO target prioritization. The method is applied to the 164 stars in the EMSL using Tycho stellar evolution tracks interpolated with kiauhoku, SPORES input properties, and Kopparapu et al. (2014) plus Turbet et al. (2023) habitable zone prescriptions. The headline result is that the CHZ2 metric peaks at 3–4 Gyr for late-F and early-G dwarfs, and that stars earlier than roughly F3 or hotter than roughly 6600 K are unlikely to have a CHZ2 at the time of observation. The paper also reports MCMC-derived masses and ages for the sample, compares them with the SPORES catalog, and releases an open-source module for computing CHZ2 posteriors on arbitrary stellar grids.","tokens_in":18101,"tokens_out":4081,"duration_ms":42759,"significance":"If the central claims are robust, this paper provides a practical, physically motivated ranking tool for the most observationally expensive decision facing HWO: which ~25 systems to observe for biosignatures. The application to the actual EMSL sample, the detailed MCMC setup with convergence diagnostics, the comparison of derived masses and ages with SPORES, and the release of reproducible code within kiauhoku are genuine strengths. The paper is also honest about several limitations, notably modeling difficulties for low-mass K/M dwarfs and the need for future comparisons across stellar grids and HZ prescriptions. However, the headline quantitative claims rest on untested assumptions (fixed 2 Gyr CHZ duration, fixed 83 mas IWA, OWA set to infinity) and on metric values reported without uncertainties. Because the 3–4 Gyr peak and the F3/6600 K cutoff are the paper's main new results, the lack of robustness and uncertainty quantification is load-bearing rather than cosmetic. The stress-test concern lands: the manuscript does not demonstrate that the ranking and the claimed peak are stable to plausible changes in the adopted CHZ duration or IWA.","major_comments":[{"comment":"The CHZ2 metric is reported in Table 2 as a bare point estimate with no uncertainty, even though it is an integral of a posterior distribution. The headline claims about a 3–4 Gyr peak and a drop near 6600 K are made by visual inspection of Figure 5, and without per-star metric uncertainties or a statistical comparison of temperature/age bins, the reader cannot tell whether these features are significant or are scatter. Please provide metric uncertainties (e.g., by propagating the MCMC posterior) and a quantitative test of the peak and cutoff.","section":"§2.5, Eq. (3); Table 2; Figure 5"},{"comment":"The metric's lower integration limit is fixed by assuming IWA = 83 mas for all stars, converted to a physical radius using distance, but no sensitivity test is presented. The IWA directly controls which part of the CHZ2 is counted: for a nearby star the physical IWA radius may lie well inside the CHZ2, while for a more distant star it may exclude a substantial fraction of the CHZ2. A sensitivity test over a plausible IWA range (e.g., 60–100 mas) is needed to show that the target ranking and the claimed age/Teff trends are stable.","section":"§2.5, Eq. (3)"},{"comment":"The outer working angle is set to infinity, with the text stating the OWA is 'likely to be well beyond the OHZ in all cases,' but this assumption is never tested. For early-F dwarfs, which have wide HZs extending to larger orbital radii, an unlimited outer integration limit may preferentially inflate their CHZ2 metrics relative to later-type stars, potentially biasing the very peak the paper claims. The authors should either adopt a plausible finite OWA or quantify how much of the CHZ2 metric lies beyond a realistic outer working angle.","section":"§2.5"},{"comment":"The choice of a 2 Gyr continuous-habitable-zone duration, taken from Truitt et al. (2020) and tied to the Great Oxidation Event timescale, is the biological pivot of the entire prioritization metric, yet no alternative CHZ durations are explored. Given that the paper's own age uncertainties frequently exceed 1 Gyr (e.g., HD 72905: 2.25+0.79−2.21 Gyr in Table 2), it is plausible that a 1 Gyr or 4 Gyr CHZ duration would shift the age at which the metric peaks and possibly the spectral-type cutoff. A sensitivity analysis varying the CHZ duration is essential to support the central claim.","section":"§1 and §2.3"},{"comment":"The abstract states without qualification that stars earlier than ~F3 or hotter than ~6600 K are unlikely to have a CHZ2, while Section 4.1 itself labels this result 'tentative' and notes it is based on only ~10 stars with a median metric of 0.053. The abstract should carry the same qualification, and the paper should quantify how many of these stars have metric values consistent with zero once uncertainties are included, which again requires the uncertainty estimates missing from Table 2.","section":"Abstract and §4.1"}],"minor_comments":[{"comment":"In the sentence describing the Turbet et al. (2023) prescription, 'compliment' should be 'complement'.","section":"§2.3"},{"comment":"The caption and table would benefit from stating the units of the CHZ2 metric; the text describes it as a proxy for the width of the CHZ2, so specifying that it is an effective AU width would aid the reader.","section":"Table 2"},{"comment":"The color scale for Teff is not described in the caption; please state the range of the color bar and how the colors map to temperature.","section":"Figure 5"},{"comment":"The age prior in Equation (1) is written with a generic Age variable, but Table 1 reports ages in Log yr; please clarify whether the prior is applied in linear or logarithmic units.","section":"§2.4, Eq. (1)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a useful application of an existing Bayesian CHZ framework to the HWO target list, and the release of code is commendable. The main concern is that the two headline claims (the 3–4 Gyr peak and the F3/6600 K cutoff) are presented in the abstract as results while the body of the paper repeatedly qualifies them as tentative and dependent on assumptions that are not tested. I believe the authors can address this with sensitivity analyses and uncertainty propagation, which is why I recommend major revision rather than rejection. The fit to the journal is appropriate for astro-ph.EP."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a solid, readable application paper, not a fundamental breakthrough. The new thing is an integrated CHZ2 ranking metric for all 164 EMSL stars, built from the authors' earlier Bayesian method plus a water-condensation inner-HZ limit, with stellar ages, masses, and CHZ2 posteriors sampled simultaneously inside the open-source kiauhoku code. Credit is due: the method is described in enough detail to reproduce, the code and model grid are released, the mass and age comparisons to SPORES are honest and mostly consistent, and the authors flag several weaknesses themselves. The F3/6600 K cutoff is explicitly labeled tentative, the poor low-mass-star fits get their own section, and the assumption that everything outside the IWA is observable is discussed rather than hidden. I also do not see circularity: the metric is not fitted to the peak or cutoff; those emerge from stellar evolution and the adopted HZ prescriptions.\n\nThe soft spots are real but proportionate. Most important: the CHZ2 metrics are reported as point values with no uncertainties, so we cannot tell whether the differences between ranked stars matter. The paper also never tests sensitivity to the two inputs that most directly shape the result: the fixed 2 Gyr CHZ duration and the 83 mas IWA. If the duration were 1 or 4 Gyr, or if the IWA were 60 or 100 mas, the ranking and even the claimed 3–4 Gyr peak and early-F cutoff could shift. The authors openly tie the 2 Gyr to Earth's Great Oxidation Event, but since their own age uncertainties are often larger than 1 Gyr, a robustness check is not optional. The OWA=infinity approximation and the neglect of exoplanet brightness/contrast are also assumptions the paper acknowledges but does not quantify.\n\nNone of this sinks the central qualitative result: late-F and early-G stars around 3–4 Gyr being strong candidates is a natural consequence of stellar evolution and HZ physics, and the paper says it cautiously. But the quantitative ranking and the sharp cutoff need more support than this draft provides. The stress-test note overstates the fragility a bit, because the authors do flag many of the limitations; it does not overstate the missing uncertainty propagation and sensitivity checks.\n\nWho gets value: mission planners and target-selection studies, plus anyone building habitability-prioritization metrics. The paper deserves a serious referee. I would send it to review and ask for metric uncertainties and a sensitivity analysis over CHZ duration and IWA; both are doable and would materially harden the claims.","headline":"A useful, clearly written application paper that packages a Bayesian continuous-habitable-zone calculation into a target-ranking metric for all 164 EMSL stars; the headline peak and cutoff are plausible but not stress-tested, and the missing uncertainties and sensitivity analysis are the real soft spots.","tokens_in":18714,"tokens_out":1903,"would_cite":true,"duration_ms":24784,"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":"For a future direct-imaging telescope hunting Earth-like worlds, the best targets are late-F and early-G dwarfs around 3-4 Gyr old, while hotter stars rarely keep any orbit continuously habitable for 2 Gyr.","keywords":["exoplanets","habitable zone","continuous habitable zone","direct imaging","target prioritization","stellar evolution","Bayesian inference","habitable worlds observatory"],"falsifier":"Find a star hotter than about 6600 K with an age near 3-4 Gyr whose Bayesian CHZ2 posterior, computed with the same pipeline, gives a high probability of a wide continuous habitable zone at its current age; that single counterexample would directly contradict the claimed hot-star cut-off. Alternatively, rerunning the pipeline with a 3D GCM-based habitable-zone prescription and seeing whether the 3-4 Gyr peak and the F3/6600 K boundary persist would test whether the result is an artifact of the adopted 1D habitable-zone prescriptions.","tokens_in":17589,"feed_emoji":"🔭","tokens_out":7369,"duration_ms":68025,"temperature":0.7,"pith_summary":"The paper proposes a way to rank stars for a future direct-imaging mission (HWO) by how much orbital area has remained continuously habitable for 2 billion years, the time it took Earth's life to oxygenate the atmosphere. Applying a Bayesian posterior over stellar ages, masses, and the evolving habitable zone to the 164-star HWO target list, it computes a single \"CHZ2 metric\" per star: the integrated likelihood that an orbit outside the coronagraph's inner working angle is in the 2 Gyr continuous habitable zone. The central result is that this metric peaks at 3-4 Gyr for late-F and early-G dwarfs, and drops steeply for stars earlier than roughly F3 or hotter than roughly 6600 K, suggesting those hotter stars are unlikely to host a continuously habitable planet at the time of observation. If correct, this gives mission planners a quantitative way to spend scarce observing time on systems most likely to show a detectable biosphere.","feed_headline":"3-4 Gyr-old F and G dwarfs top habitability rankings","feed_subtitle":"Hotter stars can't hold a 2-Gyr continuous habitable zone, so a new metric says telescope time should skip them.","key_machinery":"The load-bearing object is the CHZ2 metric, defined as $\\int_{\\mathrm{IWA}}^{\\infty} P(\\mathrm{CHZ}_2)\\,dr$, where $P(\\mathrm{CHZ}_2)$ is the Bayesian posterior probability that an orbit at radius $r$ lies within the 2 Gyr continuous habitable zone. The CHZ2 itself comes from the fixed-duration definition: an orbit is counted only if it stays inside the evolving habitable zone for 2 Gyr, using the runaway-greenhouse inner edge, the maximum-greenhouse outer edge, and a water-condensation inner limit. Posterior probabilities are generated by Markov-chain Monte Carlo sampling over a grid of stellar evolution tracks, with effective temperature, radius, and metallicity as observed inputs, so stellar age and the CHZ2 boundaries are inferred together rather than separately.","core_discovery":"The paper's central claim is that continuous habitability, measured by the fixed-duration 2 Gyr continuous habitable zone (CHZ2), is not uniformly distributed across the target list: late-F and early-G dwarfs with ages near 3-4 Gyr carry the widest and most likely CHZ2 zones, with metric values around 0.7-0.9, while stars earlier than about F3, hotter than about 6600 K, or more massive than about 1.3 solar masses appear unlikely to have any CHZ2 at the present epoch. This is established by simultaneously sampling each star's age, mass, and CHZ2 boundaries from stellar evolution tracks, converting the sampled boundaries into a likelihood curve for each orbital radius, and integrating that curve outside a fixed 83 mas inner working angle to define the CHZ2 metric. The paper also reports that the metric is lower for subgiants and the oldest stars, and that K dwarfs, though long-lived, score lower because their habitable zones are physically narrow.","pith_inferences":["A natural extension would replace the fixed 83 mas inner working angle with per-star brightness and contrast limits, which would likely shift rankings toward cooler stars because F-star habitable zones are dimmer and farther out.","If the 3-4 Gyr peak holds for a larger sample, it implies that the expected yield of biosignature-detection surveys is maximized by weighting targets by age, not just by spectral type or proximity.","The same Bayesian machinery could incorporate a planet-occurrence prior, turning the CHZ2 metric into an expected-yield estimate rather than a geometric prioritization.","Because the cut-off near 6600 K is based on about ten stars, the sharpest test of this paper's thesis is simply to apply the pipeline to the larger HWO Tier 2 and Tier 3 catalogs."],"forward_implications":["Target selection for HWO can be ranked by CHZ2 metric, with late-F and early-G dwarfs near 3-4 Gyr receiving the highest priority.","Stars earlier than roughly F3 or hotter than roughly 6600 K can be deprioritized, since the analysis suggests they will not host a 2 Gyr continuous habitable zone today.","The metric can be folded into an existing target list as a complement to exposure-time or brightness criteria, distinguishing between systems with similar observational cost.","The open-source implementation allows the same calculation to be rerun with other stellar model grids or other habitable-zone prescriptions, so the ranking can be tested against model choices.","For stars younger than 2 Gyr, the metric is zero by construction, and for very old stars it contracts, so the peak near solar age is a direct prediction of stellar evolution plus the 2 Gyr requirement."],"supporting_citations":[{"why":"Defines the fixed-duration 2 Gyr continuous habitable zone (CHZ2) that the metric is built on.","marker":"Truitt et al. 2020"},{"why":"Supplies the runaway-greenhouse inner habitable zone and maximum-greenhouse outer habitable zone prescriptions used to evolve the HZ over time.","marker":"Kopparapu et al. 2014"},{"why":"Adds the water-condensation limit that sets the initial inner edge of the habitable zone at the zero-age main sequence.","marker":"Turbet et al. 2023"},{"why":"Establishes the Bayesian CHZ posterior likelihood method that the paper extends to whole stars rather than individual planets.","marker":"Ware et al. 2022"},{"why":"Provides the 164-star HWO target list and the 83 mas inner working angle assumption used in the metric.","marker":"Mamajek & Stapelfeldt 2024"},{"why":"Supplies the uniformly derived effective temperatures, radii, and metallicities used as observed inputs and comparison ages and masses.","marker":"Harada et al. 2024"},{"why":"Provides the kiauhoku grid-interpolation tool that allows simultaneous sampling of stellar age, mass, and CHZ2 boundaries.","marker":"Claytor et al. 2020"}],"fun_headline_variants":["3-4 Gyr F and G dwarfs get best CHZ2 scores","New metric for HWO targets: continuous habitable zone","Hotter than 6600 K stars lack continuous habitable zones","Prioritize late-F and early-G dwarfs for life search"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole ranking rests on the assumption that staying in the habitable zone for 2 billion years, the time Earth took to oxygenate its atmosphere, is a good proxy for a planet being able to produce detectable life.","fun_headline_variants_meta":{"raw":{"variants":["3-4 Gyr F and G dwarfs get best CHZ2 scores","New metric for HWO targets: continuous habitable zone","Hotter than 6600 K stars lack continuous habitable zones","Prioritize late-F and early-G dwarfs for life search"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000213,"raw_usage":{"total_tokens":1482,"prompt_tokens":1067,"completion_tokens":415,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":683,"completion_tokens_details":{"reasoning_tokens":339}},"tokens_in":683,"tokens_out":415,"duration_ms":4601,"temperature":1.0,"reasoning_tokens":339,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:51:31.936057+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Find a star hotter than about 6600 K with an age near 3-4 Gyr whose Bayesian CHZ2 posterior, computed with the same pipeline, gives a high probability of a wide continuous habitable zone at its current age; that single counterexample would directly contradict the claimed hot-star cut-off. Alternatively, rerunning the pipeline with a 3D GCM-based habitable-zone prescription and seeing whether the 3-4 Gyr peak and the F3/6600 K boundary persist would test whether the result is an artifact of the adopted 1D habitable-zone prescriptions.","supporting_citations":[],"review_version":1}