{"id":"431b0c14-3811-4144-b589-6ff0cc64d0d7","arxiv_id":"2607.07571","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":5,"one_line_summary":"Synthetic SKA-Mid observations of simulated MHD and photoevaporative disk winds show that free-free emission is detectable in hours and stacked hydrogen recombination lines are spectrally resolvable in ~10 hours.","lead":"This paper uses simulations of photoevaporative and MHD disk winds to predict that the SKA telescope can detect and spatially resolve free-free emission and hydrogen recombination lines from protoplanetary disks. A smart generalist might read it to understand what the upcoming SKA observatory will be able to measure regarding disk dispersal and planet formation.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The MHD model's ~100x higher density (attributed by the authors to modeling approach, not physics) drives the most striking results — high SNR and clear line-profile distinctions — yet this density contrast may be an artifact of the isothermal assumption.","rationale":"The reader correctly identified model inconsistency as a concern, particularly the isothermal MHD assumption, but did not trace it to the specific mechanism by which it undermines the central claim. The ~100x density contrast in the MHD model is not a minor caveat — it is the proximate cause of the MHD_B4 model's high SNR (31 vs 6–12) and its distinct line wings. The authors themselves flag this as likely a modeling artifact. Without this density contrast, the MHD models would look much more like the PE/MT models in both flux and profile, collapsing the distinguishing power that the paper highlights. Additionally, the inclination dependence is underemphasized in the summary: the paper's own Figure 4 shows that at i=40°, the line profiles converge across models, with PE sometimes broader than magnetized models. The face-on case, where distinctions are clearest, is the least observationally common. The paper is a competent feasibility study and the methodology (mocassin, ProDiMo, SKA sensitivity calculator) is sound. The issue is that the comparative conclusions — which are the most scientifically interesting part — are built on a model comparison that the authors themselves acknowledge is not apples-to-apples, and the specific parameter driving the distinction (wind density) is the one most affected by the modeling differences. The verdict should remain CONDITIONAL; the paper is suitable as an SKA science case contribution, but the claim about distinguishing wind mechanisms should be more heavily qualified, particularly regarding the inclination dependence and the density contrast's origin.","tokens_in":25565,"tokens_out":2568,"duration_ms":165507,"concrete_test":"Re-run the MHD wind model with a non-isothermal thermal structure (or, as a minimal check, re-scale the MHD wind density to match the n_H ≈ 5×10^6 cm^-3 range of the PE/MT models) and recompute the stacked Hα SNR and line profiles. If the MHD SNR drops below ~10 and the FWHM difference with PE models falls below ~10 km/s at i=40°, the claim that SKA can distinguish wind mechanisms via line profiles is not supported by these models.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim that SKA can distinguish PE from MHD winds rests on two pillars: (1) stacked Hα lines detectable for all models (SNR ≥ 6), and (2) line profiles distinguishable between mechanisms. Both pillars depend heavily on the MHD_B4 model's properties. The MHD models show n_H ≈ 4–8×10^8 cm^-3, roughly two orders of magnitude higher than all other models (n_H ≈ 5×10^6 cm^-3; Section 4.1.2). The authors explicitly state this 'huge difference is likely a consequence of the different modeling approaches.' This density contrast is what gives MHD_B4 its SNR of 31 (vs 6–12 for others) and its broad, distinct line wings at face-on inclination. If the isothermal MHD model overestimates the wind density — plausible, since the thermal structure determines the ionization fraction and thus the emitting gas mass — then the MHD line fluxes would drop substantially, potentially below detection threshold, and the line-profile distinctions would diminish. Separately, the paper's own results at i=40° (Fig. 4, right panel) show that PE and MHD_B6 profiles become nearly indistinguishable, with PE actually showing broader FWHM than MT_B4. The text acknowledges 'identifying the dominant launching mechanism will be challenging' at moderate inclinations. Yet the summary (Section 7) presents line-width discrimination as a key finding. The face-on case is the most favorable and least realistic; at typical inclinations the distinguishing power weakens considerably, and this is not reflected in the headline claims.","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","summary":"This paper presents synthetic SKA-Mid (AA4) observations of free-free continuum and hydrogen radio recombination lines (RRLs) arising from photoevaporative (PE), magneto-thermal (MT), and MHD disk wind models. The authors post-process existing simulations with mocassin (for free-free) and ProDiMo (for Hα) and use the official SKA sensitivity calculator to predict detectability at 140 pc. The main claims are: (1) SKA-Mid can spatially resolve free-free emission with ~10 h integration (peak SNR 4–13 depending on band/inclination); (2) stacking 12 Hα lines in Band 5b yields detections for all models with ~10 h integration (peak SNR ≥ 6); and (3) line profiles (FWHM, wing structure) can distinguish PE from magnetically driven winds. The paper is written as a chapter for the 'Advancing Astrophysics with the SKA – II' proceedings.","tokens_in":25825,"tokens_out":1575,"duration_ms":218902,"significance":"The paper provides a timely and useful forecast for SKA capabilities in a scientifically important area—disk wind diagnostics—that is currently poorly constrained observationally. The use of established radiative transfer codes (mocassin, ProDiMo) and the official SKA sensitivity calculator lends credibility to the quantitative predictions. The extension of the ProDiMo hydrogen excitation model to n=200 (19,900 transitions) is a concrete technical contribution. The stacking strategy for Band 5b RRLs is a practical and falsifiable observational proposal. The comparison of multiple wind models (PE, MT, MHD) within a single synthetic-observational framework is valuable for the community, even given the model inhomogeneities discussed below.","major_comments":[{"comment":"Section 4.1.2 and Figure 4: The central claim that SKA can distinguish PE from MHD winds via line-profile FWHM is heavily dependent on the MHD_B4 model, which has n_H ≈ 4–8×10^8 cm^-3, roughly 100× higher than all other models (n_H ≈ 5×10^6 cm^-3). The authors state this 'huge difference is likely a consequence of the different modeling approaches.' The MHD model (Lesur 2021) is isothermal and independent of stellar properties, meaning the thermal structure that determines ionization (and thus free-free and Hα emissivity) is not self-consistently computed. If the isothermal assumption overestimates the wind density, the MHD_B4 fluxes (SNR=31) and its broad, distinct line wings would be substantially reduced. The paper should more explicitly quantify or bound this sensitivity: e.g., by how much would the MHD_B4 density need to decrease before its stacked Hα SNR falls below the detection阈值","section":null},{"comment":"Section 4.1.2, Figure 4 (right panel), and Section 7: At i=40°, the PE and MHD_B6 line profiles become nearly indistinguishable, and the PE model actually shows a broader FWHM than MT_B4. The text acknowledges that 'identifying the dominant launching mechanism will be challenging' at moderate inclinations. However, the Summary (Section 7) presents line-width discrimination as a key finding without this caveat. The summary should be revised to reflect that the distinguishing power is strongly inclination-dependent and is most robust only for the face-on case (which is the least typical geometry).","section":null},{"comment":"Section 4.1.1: The models use different computational domains (MT: 0.5–15 au; PE/MHD: 0.3–60 au) and different surface density profiles (MHD has a shallower profile). The authors note this 'somewhat limits a direct comparison.' However, the inner radius difference (0.5 au vs. 0.3 au) is particularly relevant for the MHD models, where the inner disk edge may be the origin of high-density, high-velocity flow. The text itself notes (Section 4.1.2) that the MT model's larger inner radius 'might have an impact on the line profiles and fluxes' and that line wings 'could be more pronounced in reality.' This is a load-bearing caveat for the line-profile comparison claims and should be more prominently discussed, ideally with a brief estimate of how the 0.5 au inner boundary affects the MT model's predicted FWHM relative to the MHD models.","section":null}],"minor_comments":[{"comment":"Section 4.2: The text states 'in one hour we can detect free-free with a SNR of 5' but the preceding discussion focuses on 10 h and 100 h integrations. It would help to clarify which band, model, and inclination this 1-hour SNR=5 refers to.","section":null},{"comment":"Figure 5: The y-axis label 'peak line flux [Jy]' spans 10^2 to 10^6, which seems unusually large for a protoplanetary disk at 140 pc. Please verify the units (should these be μJy or mJy?).","section":null},{"comment":"Section 4.2: The RFI exclusion range (10.7–12.7 GHz) is mentioned but the number of surviving lines used for stacking (stated as 12) should be cross-checked against the total number of Hα lines in Band 5b shown in Figure 5.","section":null},{"comment":"Figure 7: The y-axis label 'Flux (μJy/beam)' appears to have a formatting issue (the μ symbol renders as a box in some readers). Please verify the encoding.","section":null},{"comment":"Section 4.1.2: The accretion luminosity is quoted as L_accr ≈ 0.3 L_☉ in the figure caption for Figure 2, but Section 4.1.1 states L_accr = 2.6×10^-2 L_☉ for the free-free models. Please clarify which value applies to which set of models.","section":null},{"comment":"Section 6.2: The GRAVITY+ synergy section focuses almost entirely on massive YSOs (MYSOs), while the rest of the paper concerns classical T Tauri stars. A brief statement on whether GRAVITY+ can also access T Tauri systems would improve coherence.","section":null},{"comment":"The paper would benefit from a concise table summarizing all model parameters (domain, surface density profile, L_accr, L_X, β, code used) to help the reader track the differences across models.","section":null}],"recommendation":"minor_revision","confidential_remarks":"This is a proceedings chapter (AASKAII), so the bar for novelty and completeness should be calibrated accordingly. The core methodology is sound and the predictions are useful for the community. The main concern—the MHD model's isothermal assumption and its impact on the density contrast that drives the strongest detections—is real but can be addressed by more careful framing and caveating rather than new calculations. The authors should be asked to temper the summary claims to match what the models actually demonstrate, given the acknowledged inhomogeneities. No concerns about citation patterns or scope fit."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and constructive report. The referee raises three major comments, all of which concern the robustness of our line-profile discrimination between PE and MHD winds. We agree that these caveats deserve more prominent treatment and will revise the manuscript accordingly. Two of the three comments can be fully addressed through revised text and a quantitative estimate; for the third, we can provide a bounding argument but cannot fully resolve the underlying model inhomogeneity without new simulations that are beyond the scope of this proceedings contribution.","responses":[{"response":"The referee is correct that the MHD_B4 model's high density (n_H ≈ 4–8×10^8 cm^-3) is the primary driver of its high stacked Hα SNR (31), and that this density is a consequence of the isothermal assumption in the Lesur (2021) model, which does not self-consistently compute the thermal structure. We agree that this sensitivity should be quantified. We can provide a bounding estimate as follows. The stacked Hα peak flux scales approximately linearly with the emission measure, i.e., as n_e^2 × V (for optically thin emission) or roughly as n_e × V (for optically thick emission). The Hα emitting region in MHD_B4 is compact and partially optically thick, so the scaling is intermediate. Taking the conservative (optically thin) case, the SNR scales as n_H^2. The detection threshold is SNR ≈ 6, so the MHD_B4 density would need to decrease by a factor of sqrt(31/6) ≈ 2.3 for the stacked SNR to fall to the detection threshold. In the optically thick limit, the scaling is linear in n_H, giving a required decrease factor of 31/6 ≈ 5.2. Thus, the MHD_B4 density would need to decrease by a factor of roughly 2–5 (depending on optical depth regime) before its stacked Hα detection becomes marginal. We note that even a factor of 5 reduction would still leave n_H ≈ 10^8 cm^-3, well above the other models, and the broad line wings (driven by the wind velocity structure, not density) would be unaffected. We will add this quantitative bounding argument to Section 4.1.2 and explicitly state that the MHD_B4 flux predictions carry a systematic uncertainty tied to the isothermal assumption. We will also add a sentence noting that the line-profile FWHM distinction between PE and magnetically driven winds is driven by the wind velocity field, which is a more robust prediction of the MHD launching机制","revision_made":"partial","referee_comment":"Section 4.1.2 and Figure 4: The central claim that SKA can distinguish PE from MHD winds via line-profile FWHM is heavily dependent on the MHD_B4 model, which has n_H roughly 100x higher than all other models. The isothermal assumption may overestimate the wind density. The paper should quantify or bound this sensitivity: e.g., by how much would the MHD_B4 density need to decrease before its stacked Hα SNR falls below detection threshold?"},{"response":"We fully agree. The body text (Section 4.1.2) already acknowledges that 'identifying the dominant launching mechanism will be challenging' at moderate inclinations and that Keplerian broadening washes out the wind-velocity differences. However, the Summary (Section 7) does not carry this caveat and presents line-width discrimination as a general finding. This is an oversight. We will revise the Summary to explicitly state that the line-profile discrimination between PE and magnetically driven winds is most robust for face-on or near-face-on inclinations, and that at moderate inclinations (i ≳ 40°) Keplerian broadening reduces the distinguishing power, though spectral features such as the shoulder at ±20 km/s in the MHD_B6 profile may still provide diagnostic information if the SNR is sufficient.","revision_made":"yes","referee_comment":"Section 4.1.2, Figure 4 (right panel), and Section 7: At i=40°, the PE and MHD_B6 line profiles become nearly indistinguishable, and the PE model actually shows a broader FWHM than MT_B4. The Summary presents line-width discrimination as a key finding without this caveat. The summary should be revised to reflect that the distinguishing power is strongly inclination-dependent and is most robust only for the face-on case."},{"response":"The referee correctly identifies the inner radius difference as a load-bearing caveat for the line-profile comparison. The text already notes that the MT model's larger inner radius 'might have an impact on the line profiles and fluxes' and that line wings 'could be more pronounced in reality,' but we agree this deserves more prominent treatment and at least a rough quantitative estimate. We can provide the following: the MHD_B4 model's FWHM of ≈110 km/s is dominated by the high-velocity flow launched from the innermost disk region (r < 1 au), where the Keplerian velocity is v_K ≈ 30 km/s at 1 au and increases as r^-1/2. At 0.3 au, v_K ≈ 55 km/s; at 0.5 au, v_K ≈ 42 km/s. The wind velocity in MHD models is typically a fraction of the local Keplerian speed, so the inner boundary at 0.5 au (MT) vs. 0.3 au (MHD) implies a difference in maximum launch velocity of roughly 30%. If the MT model's inner radius were extended from 0.5 au to 0.3 au, we would expect the high-velocity wings to be more pronounced, potentially increasing the FWHM by on the order of 10–20 km/s. This would not close the gap with MHD_B4 (FWHM ≈ 110 km/s vs. MT_B4 ≈ 40 km/s face-on), but it could make the MT_B4 profile more similar to the MHD_B6 profile. We will add this estimate to Section 4.1.1 and move the caveat to a more prominent position, making clear that the MT model's FWHM should be considered a lower bound. We cannot, however, fully resolve this issue without rerunning the MT simulations with a smaller inner radius, which is beyond the scope of this proceedings paper.","revision_made":"partial","referee_comment":"Section 4.1.1: The models use different computational domains (MT: 0.5–15 au; PE/MHD: 0.3–60 au) and different surface density profiles. The inner radius difference (0.5 au vs. 0.3 au) is particularly relevant for the MHD models, where the inner disk edge may be the origin of high-density, high-velocity flow. This caveat should be more prominently discussed, ideally with a brief estimate of how the 0.5 au inner boundary affects the MT model's predicted FWHM relative to the MHD models."}],"tokens_in":25611,"tokens_out":1605,"duration_ms":156110,"standing_objections":["The MHD_B4 model's density structure is a direct output of the isothermal Lesur (2021) model and cannot be self-consistently improved without new MHD simulations that include thermal physics. We can bound the sensitivity (as described above) but cannot eliminate the systematic uncertainty within this paper."]},"desk_editor":{"model":"glm-5.2","letter":"This is a competent feasibility study for SKA-Mid observations of ionized gas in protoplanetary disks. The two things to know: (1) the paper shows SKA-Mid AA4 can spatially resolve free-free emission at 140 pc with ~10 hours and detect stacked Hα recombination lines in Band 5b for all tested wind models, and (2) the claim that line profiles can distinguish photoevaporative from MHD winds is real but fragile — it works best face-on and depends heavily on one model's density structure.","headline":"Solid SKA feasibility study; the MHD density contrast does real work but the paper is honest about it","tokens_in":26558,"tokens_out":175,"would_cite":false,"duration_ms":63616,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"SKA can resolve disk winds and distinguish wind mechanisms","keywords":[],"falsifier":"If real protoplanetary disk winds at 140 pc produce free-free emission below the predicted flux levels, or if recombination line profiles from photoevaporative and MHD winds are more similar than the models suggest (e.g., due to Keplerian broadening washing out wind-velocity differences at moderate inclinations), SKA observations would not cleanly distinguish the mechanisms.","tokens_in":25638,"feed_emoji":"📡","tokens_out":1154,"duration_ms":231542,"temperature":0.7,"pith_summary":"This paper asks whether the upcoming SKA-Mid telescope, in its full AA4 configuration, can detect and spatially resolve the ionized gas component of protoplanetary disk winds at nearby star-forming distances (140 pc). The authors simulate free-free continuum emission and hydrogen radio recombination lines from three classes of wind models—pure photoevaporative, magneto-thermal, and magnetohydrodynamic—and generate synthetic SKA observations to test detectability. The central claim is that SKA-Mid can spatially resolve free-free emission from disk winds with roughly 10 hours of integration (peak signal-to-noise of 4–13 depending on band and inclination), and that stacking the twelve hydrogen recombination lines within Band 5b yields spectrally resolved detections for all tested wind models at similar integration times (peak SNR of 6–31). The paper further argues that the line profiles and continuum spectral slopes differ enough between thermally driven and magnetically driven winds that SKA observations could distinguish the dominant wind-launching mechanism, with magnetically driven winds producing faster, denser, and more compact emission.","feed_headline":"SKA can resolve ionized disk winds in 10 hours","feed_subtitle":"Synthetic observations show free-free emission and stacked hydrogen lines distinguish photoevaporative from magnetic winds at 140 pc","key_machinery":"The central objects are synthetic SKA-Mid observations of free-free continuum emission and hydrogen radio recombination lines (Hα at cm wavelengths), generated by post-processing three classes of disk wind simulations: photoevaporative (PE), magneto-thermal (MT), and magnetohydrodynamic (MHD). The discriminant between wind mechanisms is the combination of continuum spectral slope and recombination line profile shape.","core_discovery":"The paper's central result is a feasibility demonstration: synthetic observations show that SKA-Mid AA4 can both spatially resolve free-free continuum emission from protoplanetary disk winds at 140 pc and spectrally resolve stacked hydrogen recombination lines in Band 5b with approximately 10 hours of integration. The discriminant power between wind mechanisms comes from two observables: the free-free spectral index, which is steeper (approximately 1.97, fully optically thick) for magneto-thermal winds versus shallower (approximately 1.13, partially optically thick) for photoevaporative winds, and the recombination line widths, which are broader for magnetically driven winds (FWHM up to 110+","pith_inferences":["The claim that line profiles distinguish wind mechanisms depends on the wind models accurately representing real disk physics. The models use inconsistent computational domains and disk surface density profiles, and the MHD model is isothermal, so the discriminant power may be over- or under-estimated relative to real disks with self-consistent thermal structure.","If the MHD model's high density (two orders of magnitude above other models) is an artifact of the semi-analytic approach rather than a physical prediction, the recombination line fluxes for magnetically driven winds could be significantly lower in reality, making detection harder than simulated.","The stacking strategy assumes all twelve Band 5b recombination lines originate from the same spatial region; if real disks have more complex ionization structure than the models, stacking may not improve SNR as predicted.","Radio frequency interference from satellite mega-constellations in the 10.7–12.7 GHz range removes several lines from the stacking analysis, and growing RFI could further degrade the effective sensitivity by the time SKA reaches full operation."],"forward_implications":["If SKA achieves the predicted sensitivity, it would provide the first systematic survey of ionized gas in protoplanetary disks at cm wavelengths, accessing a wind component that optical/IR observations probe only indirectly.","The ability to distinguish photoevaporative from MHD winds via line profiles would constrain the relative contributions of thermal versus magnetic disk dispersal mechanisms, which is currently a major open question in disk evolution theory.","A Band 5b survey of a star-forming region like Ophiuchus could be completed in approximately 60 pointings, making population-level studies of disk wind incidence feasible.","Multi-epoch observations could separate variable non-thermal (gyrosynchrotron) emission from steady free-free emission, isolating the wind signal.","Combined with JWST infrared spectroscopy and GRAVITY+ near-IR interferometry, SKA would bridge spatial scales from sub-au wind-launching regions to tens-of-au extended outflows."],"fun_headline_variants":["SKA-Mid can spatially resolve ionized disk winds at 140 pc","Free-free spectral index distinguishes photoevaporative from magnetic disk winds","SKA synthetic observations resolve stacked hydrogen lines in 10 hours","SKAO can detect magneto-thermal winds via steep free-free spectral indices","Synthetic SKA data shows 10-hour feasibility for ionized disk wind detection"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The synthetic observations rely on wind models that are not directly comparable to each other: they cover different radial ranges, use different disk surface density profiles, and the MHD model assumes an isothermal disk independent of stellar properties. The claim that SKA can distinguish wind mechanisms depends on these models faithfully representing the density and velocity structure of real winds.","fun_headline_variants_meta":{"raw":{"variants":["SKA-Mid can spatially resolve ionized disk winds at 140 pc","Free-free spectral index distinguishes photoevaporative from magnetic disk winds","SKA synthetic observations resolve stacked hydrogen lines in 10 hours","SKAO can detect magneto-thermal winds via steep free-free spectral indices","Synthetic SKA data shows 10-hour feasibility for ionized disk wind detection"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":657,"prompt_tokens":559,"completion_tokens":98,"prompt_tokens_details":null},"tokens_in":559,"tokens_out":98,"duration_ms":49569,"temperature":1.0,"reasoning_tokens":null,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-09T06:34:11.154517+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If real protoplanetary disk winds at 140 pc produce free-free emission below the predicted flux levels, or if recombination line profiles from photoevaporative and MHD winds are more similar than the models suggest (e.g., due to Keplerian broadening washing out wind-velocity differences at moderate inclinations), SKA observations would not cleanly distinguish the mechanisms.","supporting_citations":[],"review_version":1}