{"id":"ee99d60b-7b47-4230-8d27-e932dff97cd2","arxiv_id":"2607.28012","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A hydrodynamically consistent PoWR-HD grid yields a continuous mass-loss formula with a Γe kink and iron bistability jumps that matches the Arches O-to-WNh transition rate.","lead":"Researchers built 178 self-consistent wind models and turned them into a continuous mass-loss recipe for stars from about 20 to 500 solar masses at solar metallicity. The recipe tracks how winds strengthen near the Eddington limit and is checked against a model-independent rate in the Arches Cluster, so evolution codes can use it with clearer physical grounding.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified beyond the already-flagged clumping/turbulence priors; external anchors keep the central claim intact.","rationale":"The paper’s central claim is a usable, continuous solar-Z Mdot recipe (Eq. 5 plus auxiliary T*–Teff and Γe-ionisation corrections) that captures the Mdot–Γe kink and Fe bistability and is suitable for evolution in roughly the 20–500 M⊙ range. That claim is supported by (i) a large PoWR-HD grid with self-consistent hydro, (ii) a tight match to the model-independent Arches transition rate at τF,sonic~1, and (iii) ZAMS agreement with the Pauli et al. (2025) empirical Mdot–Γe trend. The reader correctly isolates prescribed clumping and vturb as the weakest structural assumption; App. A already shows ~0.2 dex typical uncertainty from vturb corrections, worst at cool T*. I do not elevate this to a verdict change: the Arches and Pauli anchors fix the absolute scale at the transition and along the hot ZAMS, so the priors set a precision floor rather than an unconstrained offset. Under-predicted v∞ (η~0.4 vs ~0.6) is disclosed and largely orthogonal to the Mdot recipe. Implementation caveats (iterative T*/Γe corrections; blurred kink at fixed Teff) are also disclosed and do not break internal consistency. Verdict remains ACCEPT; no stronger load-bearing failure mode identified.","tokens_in":39132,"tokens_out":774,"duration_ms":51189,"concrete_test":"Recompute the Arches-window subset of Fig. 8 (5.7<log L/L⊙<6.3, 25<Teff<35 kK) at fixed stellar parameters but with Dcl,max=4 and Dcl,max=25 (and, separately, vturb=30 and 100 km/s). If log Mdot at τF,sonic=1 moves by more than ~0.15 dex from −5.15, or if the kink in Γe-space shifts enough that Γe,ref in Eq. 5 is no longer consistent with τF,sonic~1, the absolute calibration of the recipe is prior-dominated rather than hydro-anchored.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing soft spot remains the one the reader named: every hydro solution (hence every coefficient in Eq. 5) is computed with prescribed micro-clumping (Dcl=1→10, onset τcl=0.1) and a fixed radially constant vturb=70.71 km/s (§2; App. A). Absolute Mdot and the detailed shape of the high-Γe branch are therefore prior-dependent at the ~0.1–0.2 dex level that App. A already quantifies for vturb, especially at cool T*. That said, the concern does not overturn the central claim. The recipe is anchored at the optically thin-to-thick transition by a model-independent observable (Arches log Mdot_trans≈−5.2), and Fig. 8 shows the models yield log Mdot(τF,sonic=1)≈−5.15±0.1 across varied X, M, and vturb in the relevant L–Teff window. The ZAMS comparison to Pauli et al. (2025) further constrains the absolute scale over a range of Γe. Prescribed clumping/turbulence is therefore a real systematic floor on precision and on cool-supergiant applicability, not an internal contradiction or a free absolute scale.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The authors present a continuous theoretical mass-loss recipe for massive and very massive stars at Z=0.02, derived from a grid of 178 hydrodynamically consistent PoWR-HD non-LTE wind models spanning M⋆≈16–500 M⊙, log(L⋆/L⊙)=5.5–6.8, T⋆=12–50 kK, and X=0.01–0.9. They confirm a kink in the Ṁ–Γe relation (shallow slope ∼2.8 below the kink; steep ∼10 above), linked to τF,sonic∼1 and η∼0.4, and two Fe-driven bistability features near T⋆≈25 and 17 kK. These behaviours, plus explicit L⋆ and X scalings, are captured in a single fitting formula (Eq. 5) with auxiliary T⋆–Teff(τR=2/3) and incomplete-ionisation Γe corrections for evolutionary implementation. Absolute scale is tested against the model-independent Arches O-to-WNh transition mass-loss rate (predicted log Ṁ≈−5.15 vs ≈−5.2) and against the Pauli et al. (2025) empirical Ṁ–Γe trend on the ZAMS.","tokens_in":39428,"tokens_out":1618,"duration_ms":49532,"significance":"Mass loss is the dominant uncertainty for VMS evolution, pair-instability boundaries, and heavy black-hole formation. A hydrodynamically self-consistent, continuous recipe that joins optically thin O-star and optically thick WNh regimes—and is anchored to the Arches transition point rather than only to internal model normalisations—is of clear practical value for stellar evolution. Strengths include the large documented grid, explicit coefficients with uncertainties and quoted RMSE (0.12 dex on Ṁ), auxiliary implementation relations, and two external anchors (Arches transition; Pauli et al. 2025) not used to set the hydro solutions. Terminal-velocity and ionising-flux predictions add further utility. The work is a natural and useful extension of Paper I and of the η-framework of Sabhahit et al. (2022).","major_comments":[{"comment":"§2 and Appendix A: Every hydrodynamic solution (hence every coefficient in Eq. 5) is computed with prescribed micro-clumping (Dcl=1→10, onset τcl=0.1) and a fixed radially constant vturb=70.71 km s−1. Appendix A already shows that vturb variations can shift Ṁ at the ∼0.2 dex level (especially at cool T⋆). The Arches τF,sonic=1 anchor (Fig. 8) remains persuasive and is partly robust to vturb, but the manuscript does not quantify how changing the clumping law (factor or onset) would move the absolute Ṁ(Γe) curve or the high-Γe slope. For a recipe intended for evolution codes, please add a short, explicit systematic-uncertainty floor (or a limited clumping sensitivity test) so users know the precision claimed beyond the 0.12 dex fit RMSE.","section":"§2, Appendix A, Eq. (5)"},{"comment":"§3.3–3.4 and Fig. 6: The models systematically under-predict terminal velocities relative to Arches transition objects (η≈0.4 at the kink vs the η≈0.6 used in the Vink & Gräfener 2012 argument). The paper correctly notes that the recipe is parametrized in Γe rather than η, so Ṁ at the kink is still well matched. However, evolutionary applications often need both Ṁ and v∞ (wind momentum, mechanical feedback). Please state more clearly whether v∞ from the grid (or a recommended scaling) should be used with Eq. 5, and flag the v∞ under-prediction as a separate limitation rather than only as an explanation for low η.","section":"§3.1, §3.4, Fig. 6"},{"comment":"§3.3 and §4.3: Implementation requires iterating Eqs. (5), (6), and (9) from structure-code Teff(τR=2/3) and Γe,high-T. The text describes the procedure but does not demonstrate numerical stability or provide a minimal worked example (e.g., one ZAMS row from Table 2 with intermediate T⋆, Γe, and Ṁ iterates). Given the steep high-Γe branch, a 10% Γe error is stated to cost ∼0.5 dex in Ṁ; a short convergence note or pseudocode would make the “suitable for stellar evolution” claim load-bearing rather than aspirational.","section":"§3.3, Table 2, §4.3"}],"minor_comments":[{"comment":"Abstract and Conclusions quote suitability for 20–500 M⊙ (abstract) and 30–500 M⊙ (conclusions). Align the stated mass range.","section":"Abstract, §5"},{"comment":"Fig. 1 caption refers to “best-fit relations described in Sect. 3.3” before that section; a forward reference is fine but ensure all solid-line sequences are uniquely identifiable against Table 1 sequence indices (as in Fig. 2).","section":"Fig. 1"},{"comment":"Eq. (5): T2 is written as depending on Γe/Γe,ref while T1 is fixed; a one-sentence physical motivation for allowing only the cooler dip centre to shift would help readers.","section":"§3.3, Eq. (5)"},{"comment":"Fig. 9 vs Fig. 1: the apparent high-Γe slope difference at fixed Teff versus fixed T⋆ is important for users; consider a brief callout in the Fig. 9 caption pointing to the T⋆–Teff correction in §3.3.","section":"Fig. 9, §4.2"},{"comment":"Table C.1 is very long; if journal length is an issue, consider moving the bulk to CDS/Zenodo and retaining a representative subset plus a machine-readable full table statement in the text.","section":"Appendix C"},{"comment":"Minor typography: “3 turb”, “3∞”, “PoWR hd/PoWRhd” appear in several spellings; unify notation for vturb, v∞, and PoWR-HD across text and tables.","section":"Throughout"},{"comment":"References to Vink et al. (2026, submitted) and related in-press items should be updated or marked consistently at proof stage.","section":"§3.2, §4.2"}],"recommendation":"minor_revision","confidential_remarks":"Solid, publishable recipe paper; the Arches and Pauli anchors keep the absolute scale from being purely internal. I do not see a reject-level or major-revision-level flaw. The clumping/turbulence prior is the main systematic floor and is already partly quantified for vturb; asking for a clearer user-facing uncertainty statement and a short implementation/v∞ clarification is enough. Scope fits A&A well as Paper II of an ongoing series."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The thing worth knowing is that this is a practical recipe paper, not a discovery of the kink. They ran 178 hydrodynamically consistent PoWR-HD models spanning L, M, T*, and X at Z=0.02, confirmed the shallow-to-steep Mdot–Γe break (~2.8 to ~10) and the two Fe bistability dips across that space, and packaged it as one continuous log-sum-exp fit (Eq. 5) plus the auxiliary T*–Teff and Γe-ionisation corrections that evolution codes actually need.\n\nWhat is new is the breadth and the implementable form. Paper I and Vink et al. (2011) already had the kink at fixed L/X; here they show it survives when you vary luminosity and hydrogen fraction, add explicit L and X scalings, and give coefficients with uncertainties and 0.12 dex RMSE. The external checks are real strengths: at τF,sonic≈1 they get log Mdot≈−5.15, matching the model-independent Arches transition rate ≈−5.2, and the ZAMS track sits close to the Pauli et al. (2025) empirical Mdot–Γe relation. That keeps circularity modest for a fitting paper.\n\nThe soft spot is the one the reader flagged and the stress-test correctly left standing: clumping (Dcl=1→10, onset τcl=0.1) and a fixed vturb=70.71 km/s are prescribed, not predicted. Every hydro solution and therefore every coefficient in Eq. 5 carries that prior. Appendix A quantifies ~0.2 dex shifts from vturb, worst at cool T*. That is a genuine systematic floor on absolute precision and on cool-supergiant use, not an internal contradiction and not enough to free the absolute scale, because the Arches and Pauli anchors still hold. Terminal velocities look systematically low relative to Arches transition objects; they note it and it does not break the Mdot recipe.\n\nMath and tables are clean, citations are appropriate (they credit the MC kink and earlier bistability work), and the auxiliary relations for structure codes are thoughtful. This is for people who put mass-loss into MESA/GENEC-type calculations of 20–500 Msun stars, black-hole mass spectra, and pair-instability boundaries. It deserves a serious referee. I would cite it when I need a continuous solar-Z VMS wind recipe, and I would send it to peer review.","headline":"Solid, usable solar-Z mass-loss recipe from a large PoWR-HD grid; kink and bistability confirmed across L/X, anchored at Arches and Pauli, with the main systematic being prescribed clumping/turbulence.","tokens_in":40149,"tokens_out":619,"would_cite":true,"duration_ms":12562,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"A single continuous mass-loss recipe captures the kink from optically thin O-star winds to optically thick WNh winds at solar metallicity.","keywords":["very massive stars","mass loss","stellar winds","Eddington parameter","Wolf-Rayet stars","hydrodynamical atmosphere models","bistability","O stars"],"falsifier":"Obtain independent mass-loss rates and Eddington parameters for Of/WNh transition stars in other young massive clusters; if the rates at the spectral thin-to-thick transition systematically miss log(Mdot) ≈ −5.15 at the matching luminosity, the absolute calibration fails.","tokens_in":39965,"feed_emoji":"💨","tokens_out":956,"duration_ms":32429,"temperature":0.7,"pith_summary":"The fates of the most massive stars are set largely by how much mass their radiatively driven winds remove. This paper builds a grid of 178 hydrodynamically consistent non-LTE wind models and turns the results into one continuous fitting formula for mass-loss rate that depends on the Eddington parameter, luminosity, temperature, and surface hydrogen abundance. The formula encodes a shallow-to-steep kink when winds become optically thick, plus two iron bistability features, and it matches the model-independent transition mass-loss rate in the Arches Cluster. Applied on the zero-age main sequence it also tracks recent empirical mass-loss versus Eddington trends. The result is a drop-in, empirically anchored recipe for stellar evolution calculations from roughly 20 to 500 solar masses at solar metallicity.","feed_headline":"One formula spans O-star to WNh mass loss","feed_subtitle":"178 hydro models yield a continuous recipe that matches the Arches transition rate","key_machinery":"A log-sum-exponential bridge between low- and high-Eddington power laws (the paper’s Eq. 5), fed by hydrodynamically consistent PoWR-HD solutions that predict mass-loss rate and velocity structure instead of assuming a prescribed beta-law.","core_discovery":"Across 178 PoWR-HD models at Z=0.02, mass loss versus classical Eddington parameter shows a kink: a shallow scaling near 2.8 for optically thin O-star winds becomes a steep scaling near 10 once winds are optically thick. The kink occurs where the flux-weighted optical depth at the sonic point is order unity and the wind efficiency is about 0.4. One continuous fit captures that kink, two iron ionisation bistability dips, and explicit luminosity and hydrogen-abundance scalings, and it recovers the Arches transition rate of about log(Mdot) = −5.2.","pith_inferences":["Because the kink tracks sonic-point optical depth of order unity, lower-metallicity hydro grids should shift the kink in Eddington space and change the initial-mass threshold for WNh-type winds.","Once depth-dependent turbulence from multi-D envelope simulations is solved inside the same hydro code, absolute rates may move by tenths of a dex, especially at cool temperatures.","Agreement on the ZAMS with empirical Mdot–Γe slopes implies that, at solar metallicity, the dominant remaining uncertainty for very massive star endpoints is less the main-sequence wind law than binary mass transfer and post-main-sequence mass loss."],"forward_implications":["Stellar evolution codes can replace separate O-star and VMS wind recipes with one continuous function over ~20–500 solar masses at Z=0.02.","Above the kink, steeper high-Eddington mass loss will change retained mass, black-hole masses, and pair-instability outcomes relative to older shallow recipes.","Below the transition the new rates lie below classic Vink et al. (2001) rates, reducing main-sequence mass loss for ordinary O stars.","Auxiliary relations convert surface Teff and high-T Eddington estimates into the inner-boundary quantities the fit needs, enabling direct use in structure codes.","The same grid supplies terminal velocities and ionising photon rates as secondary outputs across the parameter space."],"fun_headline_variants":["Mass-loss kink links thin O-star winds to thick WNh outflows","178 PoWR-HD models yield continuous O-to-WNh mass-loss recipe","Eddington scaling steepens from 2.8 to 10 at optical-depth unity","New fit captures mass-loss kink plus two iron bistability jumps","Solar-Z recipe recovers Arches transition rate near log Mdot -5.2"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"Clumping and turbulent velocity are fixed by hand in every model rather than predicted, so the entire fitted recipe inherits those prescribed choices.","fun_headline_variants_meta":{"raw":{"variants":["Mass-loss kink links thin O-star winds to thick WNh outflows","178 PoWR-HD models yield continuous O-to-WNh mass-loss recipe","Eddington scaling steepens from 2.8 to 10 at optical-depth unity","New fit captures mass-loss kink plus two iron bistability jumps","Solar-Z recipe recovers Arches transition rate near log Mdot -5.2"]},"model":"grok-4.5","effort":"low","cost_usd":0.006988,"raw_usage":{"total_tokens":1869,"prompt_tokens":1002,"num_sources_used":0,"completion_tokens":110,"cost_in_usd_ticks":69884000,"prompt_tokens_details":{"text_tokens":1002,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":757,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":1002,"tokens_out":110,"duration_ms":13720,"temperature":1.0,"reasoning_tokens":757,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T20:26:34.050822+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Obtain independent mass-loss rates and Eddington parameters for Of/WNh transition stars in other young massive clusters; if the rates at the spectral thin-to-thick transition systematically miss log(Mdot) ≈ −5.15 at the matching luminosity, the absolute calibration fails.","supporting_citations":[],"review_version":1}