{"id":"97b7f07b-b23c-4e89-a2e7-504f9364e4c3","arxiv_id":"2501.14511","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Including a constant turbulent pressure in 1D O-star atmosphere models reproduces 2D-simulated density stratifications and lowers the spectroscopic mass discrepancy.","lead":"Astronomers run 1D atmosphere models for massive O stars and compare them with 2D radiation-hydrodynamic simulations. They find that adding a turbulent pressure term in the hydrostatic equation lets the 1D models match the 2D densities and changes inferred stellar masses, potentially resolving a long-standing mass discrepancy.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fitted v_turb values inherit a known 2D temperature/opacity bias: Section 2.2 sets energy and Planck means equal to the flux mean, and Section 4.3 concedes this biases the 2D gas temperature; since the hydrostatic density fit (Eqs.","rationale":"The reader's weakest assumption is that the unweighted 2D averages from the LTE/opacity-approximated RHD simulations are a reliable target; my stress test converges on that point and sharpens it. The paper itself supplies the key evidence: Section 2.2 states that energy and Planck means are set equal to the flux mean, and Section 4.3 explicitly says this very likely biases the 2D gas temperature. Since the 1D hydrostatic solution depends on temperature through Eq. 4, a biased T(r) is degenerate with the fitted v_turb in setting the pressure scale height. This is not an external disagreement with the 2D methodology; it is an internal admission that the target of the fit is suspect in exactly the quantity most coupled to the fitted parameter. The mass-discrepancy conclusion is downstream of v_turb and inherits the same systematic. I do not see a fatal internal inconsistency: the paper is honest about the approximations and frames the mass-discrepancy resolution as potential. The requested check—refitting after correcting the 2D temperature bias—would settle whether v_turb absorbs a known systematic. If the check shows v_turb stable to within the grid step, the paper's quantitative conclusions survive; otherwise they are conditional on the 2D opacity treatment. Because the reader already assigned CONDITIONAL with high confidence and the identified concern reinforces rather than overturns that verdict, no verdict change is needed.","tokens_in":26507,"tokens_out":15167,"duration_ms":148848,"concrete_test":"Quantitative sensitivity test with corrected 2D temperature: re-extract the 2D-averaged density profile from Debnath et al. (2024) after recomputing the gas temperature with separate Planck and energy means (or, if a full rerun is unavailable, with a post-hoc correction that removes the forced T_gas ≈ T_rad equality in the wind), then repeat the Section 3.1 PoWR fitting using a chi-squared metric in log rho over x in [-0.1, 0.3]. If the best-fit v_turb shifts by more than the 25/50 km/s grid step, or if a v_turb = 0 model with re-optimized Mdot and beta reaches residuals within the 2D temporal scatter, the central claim is not robust to the known opacity-mean approximation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central inference—that 1D models need a constant turbulent pressure to reproduce the 2D-averaged densities and that this can resolve the mass discrepancy—is only as strong as the fidelity of the 2D-averaged target profiles. In Section 2.2 the authors state that the RHD models assume LTE level populations and, 'for practical reasons, both the energy- and Planck-mean opacities are assumed to be equal to this flux mean, which very likely significantly overestimates heating and cooling effects in the stellar wind.' Section 4.3 reiterates that this 'very likely results in an over-efficient heating and cooling, which in practice forces the gas and radiation temperatures to the same value, the net effect being a higher gas temperature.' The density stratification that the 1D fits target is dynamically coupled to this biased T(r) through the equation of state: in the hydrostatic equation (Eq. 3) with P = ρ(a_s^2(T)+v_turb^2) (Eq. 4), a systematically high 2D temperature produces a larger pressure scale height that the 1D model can mimic by increasing v_turb. Hence the fitted values v_turb = 35/88/106 km/s (Table 1) and the +0.1..+0.3 dex Mdot offsets are not cleanly attributable to physical turbulence. The mass-discrepancy argument (Section 4.5, Fig. 9) rests on the same v_turb values: the mechanism raises the inferred surface gravity at fixed line wings, but if v_turb partly compensates for a 2D temperature bias, the inferred mass shift is not established. The by-eye fitting and the unweighted-vs-mass-weighted averaging choice (Section 4.2, Appendix B) add further uncertainty, but the load-bearing premise is the reliability of the 2D average density/temperature profiles.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper compares 1D PoWR atmosphere models with laterally and temporally averaged 2D RHD simulation profiles from Debnath et al. (2024) for three O-type supergiant models (O8, O4, O2). The authors add a turbulent pressure term Pturb = rho v_turb^2 to the hydrostatic equation (Eq. 4), and then vary v_turb, the mass-loss rate, the wind beta exponent, and the hydrostatic/wind connection point to reproduce the 2D-averaged density, velocity, and temperature stratifications. Their main results are that a constant turbulent velocity (35, 88, and 106 km/s for O8, O4, O2) is sufficient to reproduce the 2D-averaged density profiles, that the 1D mass-loss rates must be increased by roughly 0.1 to 0.3 dex relative to the 2D values, and that the spectral line wings of the v_turb>0 models can be reproduced by lower-gravity v_turb=0 models. The latter effect is used to argue that turbulent pressure may resolve the long-standing mass discrepancy between spectroscopic and evolutionary masses. The paper also reports that turbulent pressure smooths the radiative acceleration 'dip' in the quasi-hydrostatic domain, with potential consequences for theoretical mass-loss rates.","tokens_in":26879,"tokens_out":6269,"duration_ms":59500,"significance":"If the central result holds, this paper is a useful step toward embedding insights from multi-dimensional radiation-hydrodynamic simulations into practical 1D spectral modeling: it proposes a concrete parameterization (constant v_turb in the hydrostatic equation), documents its spectral consequences, and connects the turbulent pressure to the mass-discrepancy problem. The authors are transparent about the approximations in both frameworks, and the analytic estimates in Eqs. (7) and (8) and the several appendices are valuable. However, the quantitative support for the central claim is currently incomplete. The fits to the 2D averages are performed by visual inspection with coarse parameter steps, and the 2D target profiles themselves carry a known temperature/opacity bias (Sections 2.2 and 4.3) that is degenerate with the inferred turbulent pressure. The mass-discrepancy conclusion in Section 4.5 therefore needs a sensitivity analysis before it can be considered robust.","major_comments":[{"comment":"The 'best-fit' values v_turb = 35/88/106 km/s and the mass-loss adjustments are obtained by visual inspection, with v_turb varied in steps of 25 or 50 km/s and log(Mdot) in steps of 0.25 dex, and no figure of merit is reported. Since the derived surface gravities, masses, and the mass-discrepancy ratios in Table 2 and Fig. 9 depend directly on these fitted values, please provide a quantitative residual statistic (e.g., rms or chi-square in log rho over the fitted radial range), report confidence intervals, and discuss the degeneracy among v_turb, Mdot, and the connection point.","section":"Section 3.1, Table 1"},{"comment":"The 2D benchmark models assume the energy- and Planck-mean opacities are equal to the flux mean, which, as the authors themselves state, likely overestimates heating and cooling and forces gas and radiation temperatures to the same value, yielding a higher gas temperature. Because v_turb enters the hydrostatic equation through Eq. (4), a systematically high 2D temperature produces a larger pressure scale height that a 1D model can mimic by increasing v_turb. The fitted values v_turb = 35/88/106 km/s and the mass-loss offsets are therefore not unambiguously attributable to physical turbulent pressure, and the mass-discrepancy argument in Section 4.5 inherits this ambiguity. Please quantify this sensitivity, for example by repeating the fits to 2D averages with a corrected or reduced gas temperature, and report how v_turb and the inferred masses change.","section":"Sections 2.2, 4.3, Eq. (4)"},{"comment":"For the O2 model, the lower-gravity comparison model did not converge, and Table 2 explicitly lists log g0 < 3.45, Teff0 < 40.6 kK, and M0 < 33.60 Msun as upper limits. Nevertheless, Fig. 9 presents the O2 no-turbulence case as a point in the Mevol/Mspec comparison, and the text in Section 4.5 discusses the trend without distinguishing limits from measurements. Please plot upper limits with appropriate symbols and separate the measured ratios from the bound in the mass-discrepancy claim, since the largest model in the sample is currently the least constrained.","section":"Sections 4.4, 4.5, Table 2, Fig. 9"}],"minor_comments":[{"comment":"The velocity symbol is typeset as '3' in many places (e.g., Eq. 1 and Figure captions), which appears to be a font/encoding problem; please ensure the variable v is consistently and correctly rendered.","section":"Throughout"},{"comment":"The text contains typos: 'O4 and 02 models' should be 'O4 and O2 models', and the phrase 'one might now in turn expert offsets' in Section 4.2 should read 'one might in turn expect offsets'.","section":"Section 3.1"},{"comment":"The phrase 'assuming mu ~ 0.6 as of a fully ionized plasma' should be 'assuming mu ~ 0.6 for a fully ionized plasma'.","section":"Section 4.4"},{"comment":"The evolutionary masses are selected by eye to the nearest 0.5 Msun in the HRD. For the O8 model, where Mevol/Mspec is close to 1.0 even with turbulence, please state the interpolation uncertainty explicitly or provide a more objective estimate of Mevol.","section":"Appendix G"}],"recommendation":"major_revision","confidential_remarks":"The paper's benchmark 2D models come from the same collaboration (Debnath et al. 2024), which is disclosed and is not by itself a reason for concern. However, the revised version should explicitly justify why these particular simulations are the appropriate target and ideally include a comparison with an independent multi-D calculation or with observations, so that the central claim does not rest solely on same-group models."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, know this: the paper is a prototype study, not a calibrated correction to 1D atmosphere modeling. The authors say so themselves. What's new is the first systematic 1D-versus-2D comparison for O stars, and the demonstration that a constant turbulent pressure in the hydrostatic equation reproduces the 2D-averaged density stratifications for O8, O4, and O2 supergiants. That result is real and clearly shown. The spectral consequence — adding v_turb mimics a lower surface gravity — follows directly from their Eq. 4 and is convincingly demonstrated in the model-model comparisons. The paper is also honest about its own limitations: by-eye fitting in 25–50 km/s and 0.25 dex steps, the non-converged O2 lower-gravity model, by-eye interpolation of evolutionary masses, and no observed spectra.\n\nThe main soft spot is the load-bearing premise. The target profiles come from the Debnath et al. (2024) 2D RHD simulations, which share authors with this paper and carry a known temperature bias: setting energy and Planck mean opacities equal to the flux mean very likely overestimates heating and cooling and forces gas and radiation temperatures together. Because the hydrostatic density gradient depends on temperature through the sound speed, a systematically high 2D temperature can be partly mimicked by a larger v_turb. The authors flag this in Section 4.3 but do not quantify how much of the fitted 35/88/106 km/s is physical turbulence versus compensation for the 2D temperature bias. The mass-discrepancy conclusion, which inherits those v_turb values, is therefore a motivated suggestion, not a demonstrated solution. The by-eye fitting adds further uncertainty, but it is not fatal; these are addressable issues.\n\nFor whom: anyone building or using 1D unified models for hot massive stars, and anyone who interprets spectroscopic masses of O stars. I would send it to a referee — it deserves serious review — but the referee should ask for a quantitative fitting metric, error bars on the derived masses, and ideally at least one observed O-star spectrum. A sensitivity test of the 2D temperature bias on v_turb would also be needed before I'd take the mass-discrepancy claim at face value.","headline":"A useful, honestly-flagged prototype: the constant-turbulence term in the hydrostatic equation reproduces the 2D-averaged density profiles, but the fitted v_turb values and the mass-discrepancy conclusion are not yet cleanly separated from the known temperature bias in the 2D target models.","tokens_in":27573,"tokens_out":3465,"would_cite":true,"duration_ms":30094,"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":"Standard 1D O-star atmosphere models need a turbulent pressure term to match multi-dimensional simulations, and doing so can reconcile spectroscopic and evolutionary masses.","keywords":["massive stars","O stars","stellar atmospheres","turbulent pressure","radiation-hydrodynamics","mass discrepancy","stellar winds","hydrostatic equilibrium"],"falsifier":"Compute the same 2D models with non-LTE level populations and a proper frequency-dependent treatment of the three opacity means instead of flux-limited diffusion with Planck- and energy-means set equal to the flux mean; if the averaged density stratification then no longer requires a turbulent pressure term for a 1D fit, the paper's central claim collapses. A cheaper observational check is to measure surface gravities of O supergiants with known dynamical masses: if including turbulent pressure pushes the inferred log g above the dynamical value, the correction over-shoots.","tokens_in":26241,"feed_emoji":"🌟","tokens_out":8661,"duration_ms":74086,"temperature":0.7,"pith_summary":"Standard one-dimensional models of O-star atmospheres leave out a physically real ingredient: the pressure of the turbulent, radiation-driven motions that multi-dimensional simulations find just below the photosphere. This paper asks whether that omission matters for how we read O-star spectra, and shows that it does. By fitting 1D models to the averaged density, velocity, and temperature profiles of 2D radiation-hydrodynamic simulations of O8, O4, and O2 supergiants, it finds that a single constant turbulent velocity term in the hydrostatic equation is enough to reproduce the 2D density structure, once the wind connection point is adjusted and the mass-loss rate is raised by roughly 0.2 dex for the earlier types. The same term changes the predicted line wings in the same way a lower surface gravity would, so spectroscopic masses derived without it come out too low. If the paper is right, the standard 1D analysis pipeline has been underestimating O-star masses, and the long-standing mismatch between spectroscopic and evolutionary masses has a straightforward mechanical explanation.","feed_headline":"Turbulent pressure reproduces 2D O-star atmospheres in 1D","feed_subtitle":"Adding a constant turbulence term to the hydrostatic equation raises inferred spectroscopic masses toward evolutionary values.","key_machinery":"The mechanism is a modified hydrostatic equilibrium. In the subsonic layers the pressure is $P = \\rho a_s^2$ with an effective sound speed $a_s^2 = k_B T/(\\mu m_H) + v_{\\rm turb}^2$, so the turbulent velocity $v_{\\rm turb}$ adds a pressure $P_{\\rm turb} = \\rho v_{\\rm turb}^2$; the density follows from $dP/dr = -\\rho (g - a_{\\rm rad})$, where $a_{\\rm rad}$ is the radiative acceleration from the co-moving-frame radiative transfer. Inserting $v_{\\rm turb}$ flattens the density gradient exactly as a lower gravity would, and it changes where the hydrostatic solution joins the prescribed $\\beta$-law wind velocity field, the standard analytic form $v(r) \\propto (1 - R_*/r)^\\beta$. The paper shows that a single constant $v_{\\rm turb}$ per model is sufficient, even though the 2D simulations exhibit a depth-dependent turbulence, because the deep layers where the turbulence varies are too optically thick to affect the emergent spectrum.","core_discovery":"The central claim is that a constant turbulent pressure in the hydrostatic equation of a 1D expanding atmosphere model reproduces the density stratification of the averaged 2D simulations well enough to change spectral diagnostics. The paper derives best-fit turbulent velocities of 35, 88, and 106 km/s for its O8, O4, and O2 models, respectively, close to the density-weighted values extracted from the 2D runs. With turbulence included, the 1D model's effective gravity drops, its density scale height grows, and the connection point between the quasi-hydrostatic photosphere and the beta-law wind must move outward; for the earlier O stars a roughly 0.2 dex higher mass-loss rate is also needed. In the synthetic spectra, the turbulent model produces narrower H-zeta and He ii line wings than the zero-turbulence model, mimicking the spectral signature of a lower surface gravity. The paper's conclusion is that spectroscopic masses inferred without turbulent pressure are systematically low, which offers an explanation for the mass discrepancy between spectroscopic and evolutionary masses.","pith_inferences":["Beyond the paper: the same correction should be applied to B supergiants and other stars near the Eddington limit, where sub-photospheric turbulence is expected to be strongest, since the fitted turbulent velocities increase steadily from the O8 to the O2 model.","Beyond the paper: replacing the constant $v_{\\rm turb}$ with a depth-dependent profile would likely improve the sub-photospheric density fit and may connect the hydrostatic pressure term to the photospheric macroturbulence already seen in observations.","Beyond the paper: an observational test can be built from existing high-resolution spectra of O supergiants with independently known masses; if turbulent-pressure models require surface gravities above the independently measured values, the proposed correction is too large.","Beyond the paper: the mapping from 2D to 1D depends on how the average is taken, and using mass-weighted averages for velocity while keeping unweighted averages for density could remove part of the wind-onset mismatch the paper reports."],"forward_implications":["O-star atmosphere codes should include turbulent pressure in the hydrostatic solution, because without it the modelled density stratification deviates from the structure that multi-dimensional simulations predict.","Spectroscopic masses derived from 1D models that ignore turbulent pressure are underestimated, and including the term shifts the evolutionary-to-spectroscopic mass ratio toward one.","Mass-loss rates inferred from 1D fits to early O stars are likely about 0.2 dex too low, with correspondingly larger photospheric radii and lower effective temperatures at fixed luminosity.","Turbulent pressure smooths the radiative acceleration profile in the quasi-hydrostatic region and removes the pre-wind dip, so future hydrodynamically consistent mass-loss predictions from 1D models will change once this term is included.","A single constant turbulent velocity plus an adjusted hydrostatic-to-wind connection point is enough for the density profile, but reproducing the full 2D wind velocity field still requires improved velocity-law prescriptions."],"supporting_citations":[{"why":"It supplies the three 2D radiation-hydrodynamic O-star models and their averaged profiles that the 1D fits must reproduce.","marker":"Debnath et al. (2024)"},{"why":"It establishes the multi-dimensional radiation-hydrodynamic module and box-in-a-star setup used in those simulations.","marker":"Moens et al. (2022b)"},{"why":"It provides the hybrid opacity method that accelerates the 2D wind and shapes the turbulent structure being fitted.","marker":"Poniatowski et al. (2022)"},{"why":"It underpins the 1D model atmosphere code used here, including the non-LTE treatment and the iron super-level approach.","marker":"Gräfener et al. (2002)"},{"why":"It supplies the radiative acceleration calculation inside the 1D code that enters the modified hydrostatic equation.","marker":"Sander et al. (2015)"},{"why":"It defines the mass discrepancy problem that the paper argues turbulent pressure can partially resolve.","marker":"Herrero et al. (1992)"},{"why":"It provides the evolutionary tracks used to assign evolutionary masses for the mass-ratio comparison.","marker":"Ekström et al. (2012)"},{"why":"It offers the observed O-star sample, including HD 46223, where evolutionary masses exceed spectroscopic masses.","marker":"Markova et al. (2018)"},{"why":"It shows a precedent for including turbulent pressure in hydrostatic 1D atmosphere fits, for B supergiants.","marker":"Weßmayer et al. (2022)"},{"why":"It motivates the beta-law velocity prescription and the standard beta equals 0.8 choice that the paper varies in its fits.","marker":"Pauldrach et al. (1986)"}],"fun_headline_variants":["Turbulent pressure reproduces 2D O-star density in 1D","Adding constant turbulence to 1D matches 2D O-star stratifications","1D turbulence reproduces 2D O-star density and addresses mass gap","Turbulent pressure in 1D mimics 2D O-star structure","Turbulent pressure aligns 1D and 2D O-star atmospheres"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the unweighted averages of the 2D radiation-hydrodynamic simulations, which use LTE populations, flux-limited diffusion, and flux-mean opacities, give a physically reliable target for fitting 1D models.","fun_headline_variants_meta":{"raw":{"variants":["Turbulent pressure reproduces 2D O-star density in 1D","Adding constant turbulence to 1D matches 2D O-star stratifications","1D turbulence reproduces 2D O-star density and addresses mass gap","Turbulent pressure in 1D mimics 2D O-star structure","Turbulent pressure aligns 1D and 2D O-star atmospheres"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001148,"raw_usage":{"total_tokens":4844,"prompt_tokens":1113,"completion_tokens":3731,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":729,"completion_tokens_details":{"reasoning_tokens":3628}},"tokens_in":729,"tokens_out":3731,"duration_ms":23240,"temperature":1.0,"reasoning_tokens":3628,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:04:53.068060+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same 2D models with non-LTE level populations and a proper frequency-dependent treatment of the three opacity means instead of flux-limited diffusion with Planck- and energy-means set equal to the flux mean; if the averaged density stratification then no longer requires a turbulent pressure term for a 1D fit, the paper's central claim collapses. A cheaper observational check is to measure surface gravities of O supergiants with known dynamical masses: if including turbulent pressure pushes the inferred log g above the dynamical value, the correction over-shoots.","supporting_citations":[{"cited_title":"O., Moens , N., et al","cited_arxiv_id":null,"evidence_quote":"It supplies the three 2D radiation-hydrodynamic O-star models and their averaged profiles that the 1D fits must reproduce."},{"cited_title":"G., Kee , N","cited_arxiv_id":null,"evidence_quote":"It provides the hybrid opacity method that accelerates the 2D wind and shapes the turbulent structure being fitted."},{"cited_title":"2018, , 613, A12","cited_arxiv_id":null,"evidence_quote":"It offers the observed O-star sample, including HD 46223, where evolutionary masses exceed spectroscopic masses."}],"review_version":1}