{"id":"b76ebe61-d13e-40ee-8046-98cb0b67de69","arxiv_id":"2505.16458","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"1D hydrodynamically consistent atmosphere models reproduce the average density, velocity, and mass-loss behavior of 3D Wolf-Rayet wind simulations within the simulations' own time variability.","lead":"This paper compares 1D stellar atmosphere models that compute the wind self-consistently with 3D radiation-hydrodynamic simulations of Wolf-Rayet stars. It finds that the 1D models reproduce the average density structure and can match mass-loss rates when small adjustments are made, supporting 1D models for spectral analysis.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 3D LTE/flux-mean opacity treatment could bias the benchmark mass-loss rates, not just the temperatures, so the claimed 0.2 dex 1D–3D agreement may be partly coincidental.","rationale":"The paper's core result is the 1D–3D structure comparison, and the density match is genuinely impressive: with identical L, M, R, and tuned τ_Ross boundaries, the 1D PoWR_HD models reproduce the 3D average density over several orders of magnitude. The mass-loss offset is small. However, the statement that differences are reconciled within the 3D dispersion depends on the 3D models being physically reliable. The LTE-like closure in the 3D benchmark (Sect. 2.2, Sect. 3.1.1) is the strongest unquantified systematic: the authors invoke it to explain the temperature mismatch, but they do not estimate its effect on the mass-loss rate. Since Ṁ is set near the sonic point, where the gas temperature enters the critical-point condition, a biased 3D temperature can bias the target Ṁ. The proposed fixed-structure post-processing test would isolate this effect. The reader's concern about spurious Teff offsets is part of the same issue; I agree partially, but I place the weight on the mass-loss comparison because it is the headline quantitative result. Other potential concerns (non-monotonic velocity handling for Γ2, unquantified error bars on the 3D dispersion, sensitivity to vDop) are acknowledged in the paper and are addressable, but they are less direct threats to the central claim. I therefore do not change the reader's CONDITIONAL verdict; I reinforce the condition and specify a test that would either resolve it or expose a larger problem.","tokens_in":22225,"tokens_out":16513,"duration_ms":145323,"concrete_test":"For the Γ3 model, take the time-averaged 3D density and velocity profiles from Moens et al. (2022a) and feed them into the PoWR radiative transfer in a fixed-structure mode to compute the non-LTE gas temperature and the corresponding radiative acceleration. Then compare the sonic-point condition (Eq. 4) obtained with the non-LTE temperature against that obtained with the 3D LTE (T_gas = T_rad) temperature. If the implied mass-loss rate shifts by more than ~0.1 dex, the LTE closure materially contributes to the 1D–3D Ṁ offset; if the shift is negligible, the approximation can be safely confined to the optically thin temperature diagnostics.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that 1D PoWR_HD models reproduce the averaged 3D density structure with mass-loss rates only about 0.2 dex higher, and that the offset is reconciled by 3D dispersion. This rests on the validity of the 3D benchmark. The paper is transparent that the 3D models approximate the energy and Planck mean opacities by the flux mean, forcing T_gas = T_rad (Sect. 2.2 and Sect. 3.1.1). The authors invoke this to explain the temperature mismatch, but they do not bound its effect on the mass-loss rate. The wind is launched around the sonic point, where the gas temperature enters the sound speed and the critical-point condition in Eq. (4); a systematically high 3D gas temperature can shift the sonic point and hence the predicted 3D mass-loss rate. Table 1 shows 1D–3D offsets of 0.18, 0.13, and 0.03 dex in log Ṁ. If a significant part of the 3D Ṁ is distorted by the LTE-like closure, the claimed '0.2 dex overprediction' and its reconciliation by 3D dispersion would be an artifact of the benchmark rather than a genuine property of the 1D models. The paper explicitly leaves the optically thin temperature mismatch open, so a quantitative test of the LTE effect on mass loss is needed before the headline agreement can be considered robust.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper compares 1D, stationary, non-LTE PoWR_HD models with the time- and density-averaged stratifications of three 3D radiation-hydrodynamic simulations of classical WR stars, the Γ2, Γ3, and Γ4 models of Moens et al. (2022a). In the PoWR_HD branch, the mass-loss rate and velocity field are iterated to satisfy the momentum equation (Eq. 4) using the CMF radiative acceleration, so both are predictions rather than inputs. Keeping the 3D stellar parameters (M = 10 M_sun, log(L/L_sun) = 5.47–5.74, Rc = 1 R_sun) and tuning the inner Rosseland optical depth to approximately match the 3D models, the authors find that the 1D density stratification closely tracks the averaged 3D density; that the 1D mass-loss rates are higher by 0.03–0.18 dex; that the 1D terminal velocities, and velocities at 6 R*, are higher; and that the 1D effective temperatures are lower by 0.14–0.22 dex. They then test modified inputs: fixing the 1D mass-loss rate to the 3D value requires either a mass increase (to 11.3 M_sun for Γ2 and Γ3, to 10.2 M_sun for Γ4) or a luminosity decrease (0.08 dex for Γ3, 0.02 dex for Γ4; no such solution was found for Γ2). Additional variations of the Doppler velocity, turbulent velocity, and clumping contrast are explored, and synthetic UV/optical spectra are computed for all variants, with the C IV 1550 Å P-Cygni profile identified as the main wind diagnostic.","tokens_in":22552,"tokens_out":22419,"duration_ms":158924,"significance":"If the central comparison is robust, this paper provides a valuable validation of 1D hydrodynamically consistent modeling for classical WR stars: a 1D stationary non-LTE code reproduces the averaged density and, to a lesser degree, the velocity structure of 3D time-dependent simulations at roughly the level of the 3D run-to-run dispersion (~0.2 dex in mass-loss rate). The comparison is non-trivial because the 1D mass-loss rates and terminal velocities are genuine outputs of the momentum equation (Eq. 4) and are not fitted to the 3D results; the luminosity and mass adjustments of Sect. 3.1.2 are therefore sensitivity tests rather than circular refits. The paper is transparent about the known weaknesses of both frameworks, including the flux-mean closure of the 3D models, the boundary spike in the 3D average velocity, the non-monotonic 1D velocity solution for Γ2, and the lack of a converged Γ1 model. The spectral synthesis section makes the model differences falsifiable through specific line diagnostics (C IV 1550 Å, N V 4604 Å, N III 4640 Å).","major_comments":[{"comment":"The headline quantitative claim — that 1D PoWR_HD models overproduce the 3D mass-loss rates by at most 0.2 dex and that this is reconciled by physically small adjustments — rests on the 3D reference values, yet the paper does not assess whether the LTE-like closure of the 3D benchmark (energy and Planck mean opacities replaced by the flux mean, forcing T_gas = T_rad; Sect. 2.2 and 3.1.1) biases those reference mass-loss rates. The temperature bias is invoked to explain the 1D–3D temperature differences, but the gas temperature also enters the sound speed and the critical-point condition of Eq. (4), so a systematically distorted T_gas in the wind-launching region could shift the emergent 3D mass flux. Since Fig. 1 shows that the 1D and 3D temperature profiles agree reasonably well inside r ≈ 1.6 Rc, the authors may be able to argue that the launch region is largely insensitive to the closure, but that argument is not made and no quantitative bound is given. Please either locate the sonic point and estimate the implied uncertainty in the 3D mass-loss rate, quantify the sensitivity of the 3D mass flux to the closure with an explicit test, or rephrase the agreement as being relative to the published 3D benchmark with its stated approximation.","section":"Sect. 2.2 and 3.1.1; Table 1; Eq. (4)"},{"comment":"The concluding statement that 'discrepancies can largely be mitigated by small adjustments of the mass (< 2%)' is inconsistent with Table 2, which lists required masses of 11.3 M_sun for both the Γ2 and Γ3 models — a 13% increase over the nominal 10 M_sun — while only Γ4 requires the roughly 2% change. The quantitative claim about the size of the mass adjustment should be corrected or made model-specific, since the magnitude of the adjustment is part of the paper's reconciliation argument.","section":"Conclusions vs. Table 2"}],"minor_comments":[{"comment":"The conclusion that v(6 R*) is overestimated 'by up to ~400 km/s' is not supported by Table 1, which shows differences of 210, 600, and 500 km/s for Γ2, Γ3, and Γ4, respectively; the maximum is about 600 km/s.","section":"Conclusions vs. Table 1"},{"comment":"The claim that effective temperatures are lower 'by 0.2 dex (~30 kK)' holds only for Γ4; Table 1 shows differences of 0.14, 0.17, and 0.22 dex (about 19, 25, and 30 kK), so the range should be stated.","section":"Conclusions vs. Table 1"},{"comment":"The scaling in Eq. (5) is derived from a small set of Γ3 test models (three Doppler-velocity values, Fig. B.1); the text should state this limited basis explicitly when presenting the scaling, since the next sentence contrasts it with the OB-star scaling of Björklund et al. (2021).","section":"Sect. 3.1.3; Eq. (5); Fig. B.1"},{"comment":"Please fix the typos: 'the the hydrodynamically-consistent' (Sect. 2), 'Noteably' (Sect. 3.1.1), 'preform' (Sect. 2.2), and 'disussed' (Sect. 2.1).","section":"Various sections"},{"comment":"The Table 2 note 'All fundamental parameters are scaled at τ = 2/3' should read 'defined at' or 'evaluated at'; the Fig. 3 caption 'including a turbulence v_turb' would be clearer as 'including different turbulent velocities v_turb'.","section":"Table 2 note; Fig. 3 caption"},{"comment":"For the Γ2 case, Sect. 3.1.2 states that no solution with decreased luminosity could be found; a sentence explaining whether this reflects a physical limit (e.g., the reduced radiative acceleration no longer sustaining the assumed mass-loss rate) or a numerical convergence issue would help the reader interpret the difference from the Γ3 and Γ4 cases.","section":"Sect. 3.1.2"}],"recommendation":"major_revision","confidential_remarks":"The principal risk to the headline result remains the unquantified sensitivity of the 3D benchmark mass-loss rates to the flux-mean closure; a targeted test or an explicit sensitivity argument in the launch region would resolve it. The numerical inconsistencies in the conclusions (mass adjustment '<2%' versus the 13% required for Γ2 and Γ3, and '~400 km/s' versus the 600 km/s maximum velocity difference) suggest that the summary statements need a careful audit against Tables 1 and 2. The paper fits A&A's scope well and the citation practice to the prior 3D and 1D work is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hey,\n\nTwo things to know: this is the first direct comparison of 1D hydrodynamically consistent PoWR_HD models with 3D radiation-hydrodynamic simulations for classical WR winds driven by the hot iron bump, and the headline agreement (density matched, mass-loss within ~0.2 dex) is credible but rests on a 3D benchmark with an unquantified LTE closure. The paper is worth reading, just not without that caveat.\n\nThe new content is solid. The authors take the three cWR models from Moens et al. (2022a) and run 1D PoWR_HD models with identical input. The density match is genuinely good, and the mass-loss offsets are small: 0.18, 0.13, and 0.03 dex in log Mdot. The parameter exploration is useful: lowering the Doppler velocity from 100 to 50 km/s helps for the more Eddington-limited cases, and they extract a scaling d log Mdot / d log v_Dop ~ 0.26, though that slope comes from only a couple of test points. The synthetic spectra comparison adds practical value for quantitative spectroscopy.\n\nThe main soft spot is the benchmark. The 3D models approximate the energy and Planck mean opacities by the flux mean, forcing T_gas = T_rad. The authors invoke this to explain the ~30 kK temperature mismatch, but they never bound its effect on the 3D mass-loss rate. The wind is launched near the sonic point, where gas temperature enters the critical-point condition; a systematically high gas temperature could shift the sonic point and change Mdot. If that happens, the claimed \"0.2 dex overprediction\" becomes partly a comparison artifact rather than a property of the 1D models. They note the 3D time-dependent dispersion is ~0.2 dex, which covers the offsets, but a dispersion around a biased mean does not remove the bias. This is a caveat on the reference data, not a flaw in the 1D code, and it doesn't sink the paper, but it should be quantified in a revision.\n\nMinor issues: the Eq. (5) scaling is based on few models, PoWR_HD is closed so reproducibility is limited, and the optically thin temperature mismatch is left unresolved. The mass/luminosity adjustments in Sect. 3.1.2 are honest sensitivity tests, but they're tuned to the 3D values, so they don't strengthen the predictive claim.\n\nOverall, a careful, transparent paper. It deserves a serious referee and will be a reference for WR atmosphere modeling. I'd encourage the authors to estimate the LTE effect on Mdot before publication.\n\nSend it to peer review.","headline":"The paper is a useful, honest 1D-3D WR comparison, but the headline 0.2 dex agreement depends on an unquantified LTE closure in the 3D benchmark.","tokens_in":23171,"tokens_out":6350,"would_cite":true,"duration_ms":45446,"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":"One-dimensional hydrodynamically consistent atmosphere models reproduce the average 3D density structure of classical Wolf-Rayet star simulations, with mass-loss rates only about 0.2 dex higher; small mass or luminosity adjustments within…","keywords":["Wolf-Rayet stars","stellar atmospheres","radiation hydrodynamics","mass loss","stellar winds","non-LTE line transfer","iron opacity bump","1D-3D model comparison"],"falsifier":"Recompute the Γ3 simulation with the same 3D set-up but with a proper non-LTE treatment of the energy- and Planck-mean opacities (relaxing the flux-mean approximation), then compare the averaged density and temperature stratifications with the PoWR_HD 1D solution; if the density agreement degrades or the temperature discrepancy reverses, the claimed 1D-3D consistency is an artifact of the flux-mean assumption.","tokens_in":22033,"feed_emoji":"🌟","tokens_out":8191,"duration_ms":57518,"temperature":0.7,"pith_summary":"The paper asks whether a one-dimensional, hydrodynamically consistent atmosphere model can capture what three-dimensional, time-dependent radiation-hydrodynamic simulations reveal about Wolf-Rayet winds launched by the hot iron opacity bump. Using the PoWR_HD branch, the authors show that the 1D models reproduce the time-averaged density stratification of the 3D simulations rather well, with mass-loss rates typically only about 0.2 dex higher. They further show that small adjustments in stellar mass (under 2%) or luminosity (under 0.1 dex), all within the scatter of the time-dependent 3D simulations, bring the mass-loss rates into agreement. The remaining systematic differences - 1D models more radially extended, higher terminal velocities, lower effective temperatures - are traced to methodological choices such as the Doppler velocity in the co-moving-frame transfer and the flux-mean opacity approximation in the 3D code. The upshot is that 1D hydrodynamically consistent models are a viable tool for quantitative spectral analysis of classical WR stars, with quantified uncertainty from the 3D dispersion.","feed_headline":"1D models match 3D Wolf-Rayet wind structure","feed_subtitle":"Mass-loss rates differ by at most 0.2 dex; adjusting mass or luminosity within the 3D scatter brings 1D and 3D models into agreement.","key_machinery":"The central object is the hydrodynamically-consistent branch of the PoWR model atmosphere code (PoWR_HD), which integrates the stationary 1D equation of motion from the sonic (critical) point - using the radiative acceleration obtained from detailed non-LTE co-moving-frame radiative transfer - inwards and outwards, iterating the mass-loss rate and velocity field until flux consistency and conserved Rosseland optical depth are achieved. The comparison benchmark consists of density-weighted, time-averaged radial stratifications from the 3D box-in-a-star radiation-hydrodynamic simulations of Moens et al. (2022a), which use hybrid opacities (OPAL Rosseland mean near the core and CAK-like line forces with a Gayley line-strength cutoff in the supersonic regime). The paper also varies, within the 1D framework, the Doppler velocity used for the CMF opacity profiles, the turbulent pressure in the hydrodynamic solution, and the clumping density contrast, and computes synthetic UV and optical spectra to gauge the observational impact of the stratification differences.","core_discovery":"The central discovery is that the averaged density profile of a 3D radiation-hydrodynamic simulation of a classical Wolf-Rayet star with a wind launched at the hot iron bump can be reproduced by a 1D stationary, spherically symmetric, non-LTE PoWR_HD model that solves the hydrodynamic equation of motion with the radiative acceleration computed from co-moving-frame radiative transfer. For the three models Γ2, Γ3, and Γ4 from Moens et al. (2022a), the 1D models match the 3D density structure at the same inner boundary, while predicting mass-loss rates about 0.2 dex higher; adjusting the stellar mass to about 10.2-11.3 M\\odot or lowering the luminosity by 0.02-0.08 dex yields the same mass-loss rate as the 3D models. The authors interpret these adjustments as consistent with the mass-loss and luminosity dispersion of the time-dependent simulations. The 1D models launch the wind slightly further out and reach higher velocities through the hot iron bump, resulting in higher terminal velocities and lower effective temperatures; these differences propagate into synthetic spectra, most visibly in UV P-Cygni lines such as C IV 1550 Å.","pith_inferences":["If the flux-mean approximation in the 3D models is responsible for the systematically higher optically-thin gas temperatures, then the 1D non-LTE temperature stratification may better represent the real outer wind, and the 3D temperature profiles could be recalibrated against the 1D solution.","The empirically found scaling ∂log Ṁ/∂log v_Dop ≈ 0.26 for WR models suggests that the Doppler velocity acts as a tunable broadening parameter that mimics the effect of unresolved microturbulence on the radiative force; this could be tested against 3D simulations with different velocity dispersion.","The failure to converge the lower-luminosity Γ1 model indicates that the current 1D turbulent-pressure treatment cannot capture the inflated, agitated atmosphere of hot stripped stars; extending PoWR_HD with a depth-dependent turbulent pressure may bridge the gap between classical WR stars and hot (sub)dwarfs."],"forward_implications":["Hydrodynamically consistent 1D PoWR_HD models can be used for quantitative spectral analysis of classical WR stars, providing mass-loss rates, terminal velocities, and stellar-mass estimates with an uncertainty set by the 3D dispersion.","A mass adjustment of under 2% or a luminosity adjustment of under 0.1 dex reconciles the 1D and 3D mass-loss rates, meaning the 3D averages are compatible with 1D solutions within their inherent time variability.","For stars close to the Eddington limit, reducing the Doppler velocity from 100 km s^{-1} to 50 km s^{-1} improves agreement in mass-loss rate, effective temperature, and outer-wind velocity, although the optimal value is regime-dependent.","The failure of β-law models to reproduce the outer-wind deceleration seen in the 3D simulations confirms that hydrodynamically consistent 1D models are preferable for WR wind structure.","Models sharing the same transformed mass-loss rate produce similar spectra, extending the transformed-rate scaling to hydrodynamically consistent 1D models."],"supporting_citations":[{"why":"Supplies the three 3D classical WR reference models (Γ2, Γ3, Γ4) whose averaged stratifications are the comparison target.","marker":"Moens et al. (2022a)"},{"why":"Introduces and details the PoWR_HD hydrodynamically-consistent branch used to compute the 1D models.","marker":"Sander et al. (2017, 2023)"},{"why":"Establishes the 1D versus multi-dimensional comparison methodology and the turbulent-pressure treatment for O stars, extended here to WR stars.","marker":"González-Torà et al. (2025)"},{"why":"Provides the hybrid opacity and CAK-like force-multiplier approach used in the 3D simulations.","marker":"Poniatowski et al. (2022)"}],"fun_headline_variants":["1D Wolf-Rayet models reproduce 3D wind density","Wolf-Rayet 1D models align with 3D simulations, mass-loss close","1D vs 3D Wolf-Rayet: density matches, mass-loss tweakable","1D Wolf-Rayet winds match 3D density, mass-loss within 0.2 dex"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole comparison rests on the 3D simulations being trustworthy enough that their averaged density and temperature profiles are meaningful, despite their LTE-like flux-mean opacity approximation that forces gas and radiation temperatures to be equal.","fun_headline_variants_meta":{"raw":{"variants":["1D Wolf-Rayet models reproduce 3D wind density","Wolf-Rayet 1D models align with 3D simulations, mass-loss close","1D vs 3D Wolf-Rayet: density matches, mass-loss tweakable","1D Wolf-Rayet winds match 3D density, mass-loss within 0.2 dex"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001303,"raw_usage":{"total_tokens":5418,"prompt_tokens":1150,"completion_tokens":4268,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":766,"completion_tokens_details":{"reasoning_tokens":4178}},"tokens_in":766,"tokens_out":4268,"duration_ms":23813,"temperature":1.0,"reasoning_tokens":4178,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T15:00:19.367984+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the Γ3 simulation with the same 3D set-up but with a proper non-LTE treatment of the energy- and Planck-mean opacities (relaxing the flux-mean approximation), then compare the averaged density and temperature stratifications with the PoWR_HD 1D solution; if the density agreement degrades or the temperature discrepancy reverses, the claimed 1D-3D consistency is an artifact of the flux-mean assumption.","supporting_citations":[{"cited_title":"G., Kee, N","cited_arxiv_id":null,"evidence_quote":"Provides the hybrid opacity and CAK-like force-multiplier approach used in the 3D simulations."}],"review_version":1}