{"id":"843492c0-5273-41ce-a504-fc4f82b96e5c","arxiv_id":"2412.05002","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"The maximum number of Earth-like habitats in the Milky Way's disk is about 2.5×10^5 for 10% CO2 atmospheres (or 0.6×10^5 for 1% CO2), making such habitats rare.","lead":"This paper calculates how many planets in the Milky Way could be 'Earth-like habitats': rocky planets with Earth-like nitrogen-oxygen atmospheres in the habitable zone of complex life. It finds at most a few hundred thousand such planets, implying that Earth-like habitats are rare and that complex, animal-like life may need many candidate planets to arise.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported 'maximum' is not a true upper bound because XUV stability is evaluated at the mean HZCL distance rather than the outer edge, where incident flux is lowest and atmospheric stability is highest; this can systematically exclude numerous low-mass stars and increase the final count.","rationale":"The reader's weakest assumption identifies the atmospheric stability thresholds from Johnstone et al. (2021) as the most fragile input. That is an external validity concern: if different CO2 cooling or higher planetary masses stabilize atmospheres around late K or M dwarfs, the viable stellar sample grows. My concern is more internal and arguably more directly tied to the central claim. The paper labels its computation a 'maximum number' but evaluates the key atmospheric-stability filter at the mean HZCL distance rather than the most favorable location (the outer HZCL edge). Because XUV flux decreases with squared distance, the outer edge gives the largest possible stellar sample consistent with the authors' own thresholds. By using the mean, the calculation produces a typical or average estimate rather than a rigorous maximum. The magnitude could be large: the nominal XUV cutoff excludes all stars below about 0.44 M⊙, and the excluded low-mass stars are numerous; moving to the outer edge will allow a non-negligible fraction of them to enter. This could raise NEH by a substantial factor, potentially more than an order of magnitude, although still likely leaving EHs rare relative to the ~10¹¹ stars in the disk. The central quantitative claim (a maximum of 2.5×10⁵) would be wrong as stated. For these reasons, the paper requires a revision of the distance assumption and a recomputation before the headline result can be considered an upper bound. This is consistent with a conditional acceptance, not a flat rejection, because the overall framework and rarity conclusion may survive the correction.","tokens_in":54410,"tokens_out":4059,"duration_ms":233125,"concrete_test":"Recompute the αXUV filter (Section 5.2.1.2, Figures 12–13) using the outer HZCL boundary (Ramirez 2020 for 10% CO2, Schwieterman et al. 2019b for 1% CO2) instead of the mean distance, for each stellar mass bin and rotation percentile, keeping all other parameters unchanged. If the number of stars surviving the cutoff increases by more than 20%, or if any star with M⋆ < 0.44 M⊙ now passes the nominal FXUV = 35 erg s⁻¹ cm⁻² threshold, then the reported maximum is not an upper bound and the abstract and Table 3 require revision.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"In Section 5.2.1.1, the lower stellar limit αXUV is implemented by computing FXUV at the mean HZCL distance, d⟨HZCL⟩, defined as the average of the inner and outer HZCL boundaries (around Eq. 15). The text even notes that planets farther out receive less XUV flux, but then deliberately averages over the zone. For an upper bound, one must adopt the most favorable orbital location, i.e., the outer HZCL boundary, where the star's XUV flux is smallest. Using the mean instead of the outer edge excludes stars that could host a stable N2-O2 atmosphere on a planet near the outer boundary. This is not a minor detail: the nominal XUV cutoff already removes all stars below M⋆ ≈ 0.44 M⊙ at the mean distance; the most numerous stellar population. For a star of mass 0.35 M⊙, the XUV flux at the outer HZCL is lower by roughly (d_outer/d_mean)² compared to the mean, so a star that fails the threshold at d⟨HZCL⟩ can easily pass at d_outer. Consequently, the surviving stellar sample after αXUV is not the maximum possible set, and the headline NEH values (2.5×10⁵ and 0.6×10⁵) are not guaranteed to be maxima. This is an internal methodological inconsistency with the paper's stated goal of deriving a maximum, independent of external uncertainties in the atmospheric thresholds themselves.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies the Earth-like Habitat (EH) framework developed in Paper I to the Galactic disk. It constructs a model of present-day main-sequence stars with masses 0.1–1.25 M_sun from an IMF, a star formation history, and a Galactic disk mass model, then sequentially applies supernova-sterilization, metallicity, XUV-driven atmospheric stability, bolometric-luminosity, rocky-planet occurrence, ocean/subaerial-land, and large-moon requirements. Six model cases spanning two CO2 limits (10% and 1%) and nominal/minimum/maximum input choices yield maximum counts of 2.5^{+71.6}_{-2.4} x 10^5 and 0.6^{+27.1}_{-0.59} x 10^5 EHs, implying that roughly 10^3–10^6 rocky HZCL planets are needed per EH. The authors conclude that EHs are rare and that the Copernican Principle cannot be used to infer that complex animal-like life is common in the Galaxy.","tokens_in":54751,"tokens_out":8365,"duration_ms":87223,"significance":"If established as an upper bound, the result would be a valuable quantitative update to Eta-Earth estimates because it combines stellar evolution, Galactic environment, and atmospheric stability in a single framework. The paper has real strengths: it uses the public Mors stellar evolution code, documents literature-based parameter ranges in tables, separates quantifiable from currently unquantifiable requirements, and is unusually candid about its caveats. It also makes falsifiable statements about the expected rarity of N2-O2-dominated atmospheres in the HZCL. However, the central 'maximum' claim is not yet supported, because several load-bearing choices are not conservative in the direction of maximizing the count. The result is better described as a plausible central estimate under a specific set of assumptions rather than a demonstrated upper bound.","major_comments":[{"comment":"The XUV stability filter is applied at the mean HZCL distance d<HZCL> rather than at the outer HZCL boundary, where the incident XUV flux is smallest. The text explicitly notes that planets farther out receive less FXUV, but then averages over the zone. For a claimed maximum, a star should be counted if there exists any orbital position within the HZCL at which an N2-O2 atmosphere is stable; that position is the outer edge. This choice is not neutral: under the nominal FXUV,max = 35 erg s^-1 cm^-2 threshold, all stars below about 0.44 M_sun are removed, and a 0.35 M_sun star that fails at d<HZCL> can pass at the outer edge because the flux is lower by roughly (d_outer/d_mean)^2. The same issue applies to the upper bolometric limit in §5.2.2.1, where using the mean distance makes the Seff cutoff stricter than an outer-edge criterion. The headline NEH values are therefore not guaranteed maxima; the maximum case should evaluate atmospheric stability at the outer HZCL boundary and should apply the bolometric upper limit at the same outer location.","section":"§5.2.1.1, Eq. (15)"},{"comment":"The atmospheric stability thresholds (FXUV,max = 35 erg s^-1 cm^-2 for 10% CO2, FX,max = 5 erg s^-1 cm^-2 for the minimum case) are taken from Johnstone et al. (2021) simulations for 1 Earth-mass planets. The manuscript states this choice and notes that higher-mass planets may be more stable, but it does not provide a maximum-case analysis with the most favorable planetary mass or with CO2 cooling effects included. Since these thresholds are the main reason that M and late-K dwarfs are excluded from the sample, an upper-bound claim requires either adopting the most favorable threshold values that are still consistent with the adopted atmospheric models or demonstrating that the maximum case already brackets them. As written, the result depends on this unvalidated atmospheric assumption and is not a formal upper bound.","section":"§5.2.1.1 and Appendix A2"},{"comment":"Equations (2)–(4) multiply independently estimated fractions, but the paper itself notes that requirements may be positively correlated. If two necessary requirements are positively correlated, the fraction satisfying both is at least as large as the product of the individual fractions; treating them as independent can therefore lower the estimate and break the 'maximum' property. The paper handles the metallicity–βHZCL correlation by weighting, but it does not establish that all remaining implemented fractions are independent or that positive correlations cannot increase NEH. The text should either derive the maximum under an explicit independence assumption or replace 'maximum' with 'plausible estimate under nominal assumptions' in the abstract and conclusions.","section":"§3.2.3 and Eq. (4)"}],"minor_comments":[{"comment":"The notation M̄⊙ is used where the mean stellar mass is meant; this should be M̄⋆ to avoid confusion with the solar mass symbol.","section":"§4.2.4"},{"comment":"The piecewise main-sequence lifetime relation is taken from Westby & Conselice (2020), but the units of the numerical coefficients should be stated explicitly; as written, terms such as 7.1M^{-2.5} appear dimensionally inconsistent.","section":"Eq. (8)"},{"comment":"The asymmetric ranges quoted for NEH (e.g., 2.5^{+71.6}_{-2.4} x 10^5) are not labeled; please state explicitly whether these are 1σ errors, full min–max ranges, or something else, and make the same distinction in the tables.","section":"Abstract and §8 tables"},{"comment":"The caption and text describing the direction in which the minimum and maximum stable HZCL boundaries are shifted appear contradictory; please check whether the minimum case shifts the boundaries inward or outward and make the figure and text consistent.","section":"Figure 16 and §5.2.2.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is long but generally well organized and unusually candid about its limitations. My main concern is that the central quantitative claim is framed as a maximum, while the implementation uses several choices that are not conservative in the maximizing direction. The mean-vs-outer HZCL issue is the most concrete and can be fixed by rerunning the maximum case at the outer boundary and checking whether the headline numbers change. If the authors prefer to keep the mean-distance assumption, they should remove 'maximum' from the title and abstract and present the result as a nominal estimate under stated assumptions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe headline is that this paper gives the first concrete quantitative upper limit for Earth-like habitats in the galactic disk: about 2.5e5 and 0.6e5 planets for 10% and 1% CO2 atmospheres, implying EHs are rare relative to ~1e11 stars. It's a transparent, caveat-heavy synthesis of Paper I's formula with literature models for star formation, metallicity, XUV-driven atmospheric escape, planet occurrence, water, and moons. The step-by-step application is careful, and they run nominal/min/max cases across a wide parameter range. The qualitative claim—that not every star or rocky HZCL planet can host an Earth-like N2-O2 atmosphere—is well supported.\n\nThe main soft spot is internal to the 'maximum' claim. They evaluate XUV stability at the mean HZCL distance, not at the outer edge where incident flux is lowest. For a true upper bound, you'd take the most favorable orbit. At the nominal threshold, this choice excludes all stars below ~0.44 Msun; many of those low-mass stars could host stable atmospheres at the outer HZCL. So the headline numbers are not guaranteed maxima—they are maxima conditional on a particular orbital averaging. The paper acknowledges the simplification but doesn't address its implications for the upper-bound claim. The stress-test concern holds up on reading, and it's a genuine methodological flaw, though it doesn't necessarily overturn the rarity conclusion.\n\nSecond, the atmospheric stability thresholds (e.g., FXUV,max=35 erg/s/cm2) come from Johnstone et al. 2021 for 1-Earth-mass planets. If CO2 cooling or higher masses allow survival around late K/M dwarfs, the viable stellar sample grows substantially. That's an external uncertainty, and the paper is admirably upfront about it. Still, the numbers are model-dependent.\n\nMinor issues: the asymmetric ranges (e.g., +71.6/-2.4) are parameter ranges presented in the style of error bars, which can mislead; and no code/data are provided, which hurts reproducibility. The correlation between parameters is discussed but not resolved.\n\nWho's this for? Astrobiologists and SETI researchers who want a concrete, if conditional, estimate of Earth-like habitats. It deserves a serious referee—the topic is important, the work is extensive, and the flaws are identifiable and fixable. I'd recommend sending it out with a request to address the mean-vs-outer-edge issue and to reframe the numbers as 'plausible' rather than strict maxima.","headline":"A transparent and useful synthesis that delivers concrete numbers, but the 'maximum' label is undercut by evaluating XUV stability at the mean rather than the outer HZCL distance.","tokens_in":55295,"tokens_out":3963,"would_cite":true,"duration_ms":49179,"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":"At most 250,000 Earth-like habitats exist in the galactic disk","keywords":["Earth-like Habitats","habitable zone of complex life","N2-O2 atmosphere","atmospheric escape","galactic habitable zone","stellar XUV flux","exoplanet occurrence rate","astrobiology"],"falsifier":"A spectroscopic survey of rocky planets in the HZCL around M and late-K dwarfs that receive mean XUV surface fluxes above $35\\ \\mathrm{erg\\,s^{-1}\\,cm^{-2}}$ would test the central threshold: finding even a few long-lived N$_2$-dominated or N$_2$-O$_2$ atmospheres on such planets would contradict the assumed maximum.","tokens_in":54153,"feed_emoji":"🌍","tokens_out":12657,"duration_ms":112576,"temperature":0.7,"pith_summary":"The paper sets out to put an upper bound on how many planets in the Milky Way's disk can be Earth-like Habitats: rocky planets in the habitable zone of complex life that can hold an N$_2$-O$_2$-dominated atmosphere with only minor CO$_2$. Applying its maximum-count formula to the present-day disk-star population, with filters for galactic metallicity and supernova sterilization, stellar XUV-driven atmospheric escape, rocky-planet occurrence, water and land availability, and large moons, it finds maxima of $2.5^{+71.6}_{-2.4}\\times10^5$ planets for atmospheres with at most 10% CO$_2$ and $0.6^{+27.1}_{-0.59}\\times10^5$ for at most 1% CO$_2$. The consequence is that Earth-like habitats are rare: on average, at least roughly $10^3$ to $10^6$ rocky HZCL planets are needed for one Earth-like habitat to exist. These are maxima, the paper emphasizes, because several requirements that are not yet quantifiable, such as the origin of life, are set to unity.","feed_headline":"At most 250,000 Earth-like habitats exist in the galactic disk","feed_subtitle":"The true count is likely far lower, and any technological civilizations would be rarer still.","key_machinery":"The load-bearing object is the inequality $N_{\\mathrm{EH}} \\le N_\\star \\times \\eta_\\star \\times \\eta_{\\mathrm{EH}}$, where $N_\\star$ is the number of disk stars with masses between $0.1$ and $1.25\\,M_\\odot$, $\\eta_\\star$ is the fraction of those stars offering a habitable galactic and stellar environment, and $\\eta_{\\mathrm{EH}}$ is the fraction of their rocky planets that meet the planetary requirements. Each fraction is a product of requirement probabilities: supernova survival, metallicity threshold, XUV/X-ray atmospheric stability, bolometric-luminosity upper limit, rocky-planet occurrence, water plus subaerial land, and large-moon frequency. The physical hinge is the atmospheric stability threshold: for a one-Earth-mass planet at the middle of the HZCL, an N$_2$-dominated atmosphere with 10% CO$_2$ is assumed to survive only if the stellar XUV surface flux stays below $F_{\\mathrm{XUV,max}} = 35\\ \\mathrm{erg\\,s^{-1}\\,cm^{-2}}$ (with $F_{\\mathrm{X,max}} = 5\\ \\mathrm{erg\\,s^{-1}\\,cm^{-2}}$ in the minimum case), and lower thresholds for the 1% CO$_2$ case. Applying this threshold with stellar evolution tracks across the rotational distribution removes most M dwarfs and shifts the surviving sample toward K and G stars.","core_discovery":"The central claim is that the number of Earth-like Habitats in the galactic disk, while not known, is bounded above by a few hundred thousand, and that most rocky planets in the habitable zone of complex life do not become Earth-like habitats. The calculation reaches this bound by filtering the disk-star population stepwise: a star must lie in a galactic environment with sufficient metallicity and low enough supernova exposure; at the mean HZCL distance, the star's XUV and X-ray fluxes must be low enough for an N$_2$-O$_2$ atmosphere to survive thermal escape; and its luminosity must not yet have pushed the planet through the inner habitable-zone boundary. The paper then multiplies by the rocky-planet occurrence rate and by the requirement of surface water with subaerial land and a large moon, using literature minima and maxima for each factor. The result is a plausible maximum range of $2.5^{+71.6}_{-2.4}\\times10^5$ Earth-like Habitats for the 10% CO$_2$ case and $0.6^{+27.1}_{-0.59}\\times10^5$ for the 1% CO$_2$ case, with the statement that the actual number is likely lower, since several biological and geological requirements are not included. The authors conclude that the Copernican Principle cannot be used to infer that complex animal-like life is common in the Galaxy.","pith_inferences":["If the atmospheric-stability thresholds hold, the model predicts a sharp stellar-type signature: N$_2$-dominated atmospheres should be found preferentially around K-type and early M stars in a restricted birth-age window, while late M dwarfs should be systematically empty of them.","The same filter chain could be applied to other habitat definitions by replacing the atmospheric-stability and HZCL terms, which would shift the upper bound by orders of magnitude for, say, CO$_2$-dominated or H$_2$-dominated habitats.","The paper's spatial distribution of EHs could be integrated over any survey volume to build a target list for future atmospheric spectroscopy, so that observed atmospheric compositions would test the filter chain region by region."],"forward_implications":["At most a few hundred thousand Earth-like Habitats exist in the galactic disk, and the actual number is likely smaller.","On average, at least $\\sim 10^3$ to $10^6$ rocky planets in the habitable zone of complex life are needed for one Earth-like Habitat to exist.","Most M dwarfs, especially late M dwarfs, are strongly disfavored: for the stricter 1% CO$_2$ atmosphere case essentially no M dwarfs remain in the viable stellar sample.","Extraterrestrial intelligence, if it requires an Earth-like Habitat, would be much rarer than the population of rocky habitable-zone planets suggests, and the Copernican Principle cannot be invoked to infer that such life is common."],"supporting_citations":[{"why":"Defines the Earth-like Habitat and supplies the maximum-number formula $N_{\\mathrm{EH}} \\le N_\\star \\eta_\\star \\eta_{\\mathrm{EH}}$ that this paper applies.","marker":"Lammer et al. 2024"},{"why":"Upper-atmosphere simulations used to set the XUV and X-ray surface-flux thresholds for thermal stability of N$_2$-dominated atmospheres.","marker":"Johnstone et al. (2021a)"},{"why":"Provides the Mors stellar evolution package used to track XUV flux, bolometric luminosity, and HZCL distances over stellar mass and age.","marker":"Johnstone et al. (2021b)"},{"why":"Sets the inner habitable-zone Runaway Greenhouse boundary and the effective stellar flux values used for the upper stellar limit.","marker":"Kopparapu et al. (2014)"},{"why":"Defines the outer HZCL boundary for $p_{\\mathrm{CO_2}}=0.1$ bar and the CO$_2$ toxicity limit for the nominal 10% CO$_2$ case.","marker":"Ramirez (2020)"},{"why":"Provides HZCL boundaries for $p_{\\mathrm{CO_2}}=0.01$ bar and the CO$_2$ toxicity framework behind the 1% CO$_2$ case.","marker":"Schwieterman et al. (2019b)"},{"why":"Supplies the supernova-sterilization probability distributions in the galactic disk used for the galactic habitable zone filter.","marker":"Gowanlock et al. (2011)"},{"why":"Gives the present-day metallicity distribution functions of the disk used to apply the metallicity threshold.","marker":"Hayden et al. (2015)"},{"why":"Provides the star formation history and metallicity evolution tracks used to assign stellar birth ages and metallicities.","marker":"Snaith et al. (2015)"},{"why":"Provides the thin and thick disk mass model used to distribute stars radially and vertically in the galactic disk.","marker":"McMillan (2017)"}],"fun_headline_variants":["At most 250k Earth-like habitats in the galactic disk","Galactic disk caps Earth-like habitats at 250k","Maximum 250k Earth-like habitats in disk, likely fewer","Up to 250k Earth-like habitats in Milky Way disk"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the XUV and X-ray surface-flux thresholds that define atmospheric stability, taken from simulations for roughly one-Earth-mass N$_2$-dominated atmospheres, are representative of real planets; if additional cooling agents, higher planetary masses, or different atmospheric histories allow such atmospheres to survive around late K or M dwarfs at fluxes above $35\\ \\mathrm{erg\\,s^{-1}\\,cm^{-2}}$, the viable stellar sample grows substantially and the stated maximum is no longer an upper bound.","fun_headline_variants_meta":{"raw":{"variants":["At most 250k Earth-like habitats in the galactic disk","Galactic disk caps Earth-like habitats at 250k","Maximum 250k Earth-like habitats in disk, likely fewer","Up to 250k Earth-like habitats in Milky Way disk"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001426,"raw_usage":{"total_tokens":5899,"prompt_tokens":1236,"completion_tokens":4663,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":852,"completion_tokens_details":{"reasoning_tokens":4593}},"tokens_in":852,"tokens_out":4663,"duration_ms":35546,"temperature":1.0,"reasoning_tokens":4593,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T21:00:26.241057+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A spectroscopic survey of rocky planets in the HZCL around M and late-K dwarfs that receive mean XUV surface fluxes above $35\\ \\mathrm{erg\\,s^{-1}\\,cm^{-2}}$ would test the central threshold: finding even a few long-lived N$_2$-dominated or N$_2$-O$_2$ atmospheres on such planets would contradict the assumed maximum.","supporting_citations":[],"review_version":1}