{"id":"1f4c7031-7081-444b-b8ec-b6dfddec6b65","arxiv_id":"2505.16586","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Centimeter free-free emission in Class II disks correlates with accretion rate, and the authors infer it traces MHD winds or jets in both full and transition disks.","lead":"This paper reports a correlation between the ionized mass loss rate, inferred from centimeter free-free emission, and the accretion rate in 31 young stars with disks, and argues the emission comes from magnetohydrodynamic winds or jets. It compares full and transition disks to test whether the same accretion-outflow mechanism operates in both.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The MHD-wind/jet interpretation is not uniquely supported: unresolved cm emission could contain a roughly constant gyro-synchrotron component, which would preserve the reported Mdot_i–Mdot_acc correlation without a thermal jet.","rationale":"The reader already assigned CONDITIONAL and flagged the Reynolds-model conversion and alternative emission mechanisms as the weakest assumption. My concern overlaps with that but sharpens it: the most load-bearing issue is not merely the fixed parameter choices in Eq. 2, but the possibility that the unresolved 2 cm excess is not dominated by thermal free-free at all. Because gyro-synchrotron emission can also produce flat/positive spectral indices and may plausibly scale with accretion activity, the observed Mdot_i–Mdot_acc correlation could survive even without a jet or MHD wind. The paper explicitly acknowledges this contamination cannot be ruled out by SED fitting, so the physical interpretation is conditional on a testable assumption. This does not invalidate the empirical correlation between the 2 cm excess and accretion rate, which is well characterized, but it does mean the central physical claim should be regarded as conditional pending independent discrimination. The recommendation remains CONDITIONAL in substance; since the reader's verdict is already CONDITIONAL, the verdict should be UNCHANGED. Agreement is partial because the reader emphasized photoevaporative winds and parameter assumptions whereas I emphasize non-thermal contamination as the more likely confounder.","tokens_in":16992,"tokens_out":3642,"duration_ms":33926,"concrete_test":"Measure full-Stokes polarimetry of the 2 cm continuum for a representative subset of the full-disk sample (e.g., BP Tau, DN Tau, DO Tau, MWC 480) with VLA at 10–15 GHz. Thermal free-free is essentially unpolarized, whereas gyro-synchrotron typically shows a few to tens of percent circular polarization. If any source shows significant circular polarization (>~1–3%), the non-thermal fraction is non-negligible and the Mdot_i values from Eq. 2—and hence the MHD-wind/jet conclusion—are not secure. In that case, re-fit the Mdot_i–Mdot_acc correlation after excluding or down-weighting the polarized sources and check whether r remains above ~0.5.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central conclusion rests on interpreting the unresolved 2 cm excess as thermal free-free and converting it to an ionized mass-loss rate via the Reynolds (1986) jet model (Eq. 2). This conversion assumes a fully ionized, isothermal 10^4 K jet, fixed turnover frequency 40 GHz, opening-angle ratio 0.5, and v_jet ∝ M*^{1/2} (Eq. 3). The evidence that the emission is thermal free-free is the flat spectral index and the absence of correlations with X-ray and [O I] properties. Yet the paper itself notes in Sec. 4.1 that gyro-synchrotron emission cannot be subtracted by the SED fitting and 'may contribute to a not-determinable fraction of the flux'; it dismisses this because non-thermal contributions 'are likely to vary from target to target.' That dismissal is not secure: if the gyro-synchrotron fraction is roughly constant, or itself scales with accretion activity, the Mdot_i–Mdot_acc correlation (r=0.70, Fig. 6) would remain significant even if no collimated thermal jet is present. Moreover, since Mdot_i ∝ F_ff^{0.75}, the mass-loss correlation largely restates the already-reported F_ff–L_acc correlation (r=0.72) under a nonlinear rescaling and does not independently establish an MHD wind/jet. The transition-disk comparison is also partly built from lower limits, but the full-disk interpretation is the load-bearing piece.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes 31 Class II YSOs (18 full disks, 13 transition disks) using literature multi-wavelength photometry to isolate centimeter excess emission over the dust contribution. The authors find that the 2 cm free-free flux correlates with accretion luminosity (r=0.72), and that the ionized mass loss rate inferred from a Reynolds (1986) jet model correlates with the stellar accretion rate (r=0.70 for the combined sample). They rule out an X-ray-driven photoevaporative wind on the basis of the spectral indices and the lack of an X-ray correlation, and they argue that the emission originates from an MHD wind/jet. A comparison between full and transition disks yields consistent slopes, though most transition-disk fluxes are lower limits. The authors conclude that free-free emission in both classes of disks is likely associated with an MHD wind/jet, and that the apparent difference in outflow efficiency may be due to the biased accretion rates of the transition-disk sample.","tokens_in":17346,"tokens_out":5193,"duration_ms":44020,"significance":"If the interpretation is correct, the paper provides important evidence that MHD winds/jets are a common accretion-driving mechanism in Class II full disks, extending earlier transition-disk results to a wider sample. The robustness check that the correlation survives when the jet velocity is held fixed at 200 km/s (Sec. 4.3) is a genuine strength, as is the use of multi-wavelength SED fitting to separate free-free from dust emission. The paper is also honest about the limitations of the transition-disk comparison. However, the central conclusion depends on a model-dependent conversion from radio flux to ionized mass-loss rate and on the assumption that the non-dust emission is entirely thermal free-free; the significance for the field therefore hinges on resolving this degeneracy.","major_comments":[{"comment":"The conclusion that the centimeter emission is thermal free-free from an MHD wind/jet is not uniquely supported. The paper states that gyro-synchrotron emission 'may contribute to a not-determinable fraction of the flux' and dismisses this because such emission is 'likely to vary from target to target.' This argument is not quantitative: if the non-thermal fraction is roughly constant, or scales with accretion activity, the reported Mdot_i-Mdot_acc correlation would persist even if no collimated thermal jet is present. Please provide an explicit constraint on the non-thermal fraction (e.g., from circular polarization, spectral index at additional wavelengths, or variability) or rewrite the conclusions to present the jet interpretation as conditional on the thermal-free-free assumption.","section":"Sec. 4.1, Eq. (2)"},{"comment":"The correlation between Mdot_i and Mdot_acc is largely a nonlinear rescaling of the already-reported F_ff-L_acc correlation, because Mdot_i is proportional to F_ff^{0.75} for fixed jet parameters. The paper should explicitly acknowledge that the Mdot_i conversion mainly changes the normalization and does not independently establish a physical jet. Also, the r=0.70 value quoted in Sec. 3.1.3 is for the combined sample of full and transition disks, whereas the abstract and conclusions attribute the strong correlation specifically to full disks; the full-disk-only correlation is r=0.64±0.19 (Eq. b). Please state the subsample-specific coefficients in the abstract and avoid implying that the strongest correlation is exclusive to full disks.","section":"Sec. 3.1.3, Fig. 6; Abstract"},{"comment":"The transition-disk comparison is built on lower limits for the free-free flux in most targets, because the spectral indices are upper limits. The fitted TD-only slope (0.41±0.28) has a large uncertainty, and the statement that the slopes are 'fully consistent' is therefore weak. The conclusion that the same mechanism operates in both classes is not strongly constrained by these data. Please treat the TD correlations as tentative, or use a survival-analysis approach that accounts for censored data, and explicitly state how the lower limits affect the slope comparison.","section":"Sec. 4.2, Table 2 and Fig. 8"}],"minor_comments":[{"comment":"There is a typo: 'ionized gas closed to the star' should be 'ionized gas close to the star.'","section":"Abstract"},{"comment":"The caption uses 'Top:' twice; the second panel should be labeled 'Bottom:'.","section":"Fig. 2 caption"},{"comment":"The introduction contains the typo 'photevaporative' instead of 'photoevaporative'.","section":"Sec. 1"},{"comment":"The inclination dependence enters as (sin i)^-0.25, and for the most face-on disks (e.g., HD 135344B with i=12 deg) this factor is highly sensitive to the assumed inclination; a brief discussion of how inclination uncertainties propagate into Mdot_i would be useful.","section":"Eq. (2) and Table 1"},{"comment":"The spectral-index entries for the transition disks are presented in a nonstandard format (e.g., '0.58 -1.08 +0.98'); please clarify the notation in the table notes so that the central value and its asymmetric uncertainties are unambiguous.","section":"Table 2"}],"recommendation":"major_revision","confidential_remarks":"The reader's stress-test concern about gyro-synchrotron contamination is well founded and is the main reason for the major revision. The paper's own admission that the non-thermal fraction is 'not-determinable' makes the central claim load-bearing on an assumption that is not currently quantified. The manuscript is likely publishable after the authors either add a quantitative constraint on the non-thermal contribution or reframe the conclusions as conditional. The abstract's attribution of the r=0.70 to full disks only is also misleading and should be corrected. The sample-size and lower-limit issues are secondary but should be addressed as described in the major comments."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is an incremental but genuinely useful paper: it adds new diagnostics (no correlation with X-ray or [O I] properties) to the full-disk free-free sample from Garufi et al. (2025), and it makes the first quantitative full-disk vs transition-disk comparison. Second, the central claim—that the free-free emission traces an MHD wind/jet—is plausible but not uniquely supported; the conversion from flux to ionized mass-loss rate uses a fixed Reynolds jet model, and the resulting Mdot_i–Mdot_acc correlation largely restates the direct F_ff–L_acc correlation under a power-law transform.\n\nThe paper does several things well. The null results for X-ray and [O I] are interesting and honestly reported; if the emission were an X-ray-driven photoevaporative wind, you'd expect a correlation with L_X, and you don't see it. The check that fixing v_jet=200 km/s preserves the correlation is a good control against the stellar-mass dependence. The authors are also transparent about the transition-disk lower limits and the bias of that sample toward strong accretors; they don't overclaim the TD-FD difference.\n\nThe soft spots are real but proportionate. The stress-test note is right that gyro-synchrotron emission cannot be subtracted and could in principle keep the correlation alive if its fractional contribution is roughly constant. The paper dismisses this because non-thermal contributions 'are likely to vary from target to target,' but that is an assertion, not a test. Still, the flat spectral indices (mean 0.76) are more naturally thermal free-free, and the X-ray non-correlation is a point in the authors' favor. The bigger caveat is that Mdot_i is derived from the same flux that already correlates with L_acc, so the Mdot_i–Mdot_acc slope is not an independent measurement; it is a rescaled version of the same observable correlation. That doesn't invalidate the result, but it should be stated more explicitly. The TD comparison is explicitly provisional, with most fluxes as lower limits; the authors handle this with appropriate caution.\n\nOverall, this is a careful paper from people who know the literature. The interpretation leans on a model, but the observational correlations are solid and the controls are reasonable. The right audience is anyone working on disk winds and accretion; the paper will be a useful reference for the full-disk free-free sample and for the FD-TD comparison. It deserves a serious referee, with the main request being a clearer separation of the direct observable correlation from the model-dependent mass-loss rate.\n\nRecommendation: send it to review. The authors should be asked to quantify or better bound the non-thermal contamination, and to present the F_ff–L_acc correlation as the primary result, with Mdot_i–Mdot_acc as a derived interpretation rather than an independent confirmation.","headline":"Useful incremental study adding null diagnostics and an FD-TD comparison, but the MHD-wind interpretation rests on a model conversion that the data don't uniquely constrain.","tokens_in":17909,"tokens_out":3139,"would_cite":true,"duration_ms":26848,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["97.10.Fy","98.38.Er","97.21.+a","98.62.Mw"],"model":"deepseek-v4-flash","headline":"Radio free-free emission links accretion to outflow in 31 young stellar disks","keywords":["free-free emission","young stellar objects","protoplanetary disks","transition disks","MHD disk winds","jets","accretion rate","mass loss rate"],"falsifier":"A direct measurement that would settle the claim is a resolved image of the radio continuum around several of the full disks: if the centimeter emission is spread over a broad, roughly spherical region whose size scales with the X-ray or EUV luminosity rather than being collimated along the disk axis, the jet interpretation would fail. A second test is to observe a disk while its accretion rate drops dramatically; if the centimeter flux does not fall proportionally, the free-free emission is not tracing the accretion-driven outflow.","tokens_in":16832,"feed_emoji":"🌠","tokens_out":1368,"duration_ms":13878,"temperature":0.7,"pith_summary":"The paper argues that the centimeter-wavelength free-free emission seen toward evolved young stars with disks traces ionized gas in a jet or magnetohynamic disk wind, not an X-ray-heated photoevaporative flow. It compares 18 full disks and 13 transition disks, finds that the ionized mass loss rate inferred from the free-free emission correlates with the stellar accretion rate, and concludes that jets or MHD winds are likely a main driver of accretion in both disk classes. This matters because it pins the same accretion-outflow engine to disks at different evolutionary stages, including transition disks whose inner cavities raise questions about how material crosses the gap.","feed_headline":"Radio glow links disk accretion to outflow in 31 young stars","feed_subtitle":"Ionized gas near the star moves with the accretion rate, pointing to jets or MHD winds in full and transition disks alike.","key_machinery":"The central object is the ionized mass loss rate dot M_i inferred from the 2 cm free-free flux through the Reynolds (1986) jet model. That model converts the free-free flux, spectral index, jet velocity, opening angle, inclination, distance, and temperature into a mass loss rate. Because the jet velocity is assumed to scale as the square root of stellar mass, the paper also checks whether the accretion-outflow correlation survives fixing the jet velocity at 200 km/s, which it does, ruling out a stellar-mass-driven artifact.","core_discovery":"For a sample of 31 Class II young stellar objects with full or transition disks, the free-free emission detected at centimeter wavelengths most likely arises from ionized gas in a jet or MHD disk wind close to the star. The evidence is that the ionized mass loss rate derived from the free-free flux correlates strongly with the accretion rate onto the star (r=0.70±0.18), while the free-free emission shows no correlation with X-ray luminosity or with the [O I]6300 Å line properties. The spectral indices are flat or positive, consistent with partially optically thick free-free emission from a collimated outflow rather than the optically thin emission expected from an X-ray photoevaporative wind. The same conclusion holds for transition disks, whose free-free fluxes are mostly lower limits, and the two samples show consistent accretion-outflow correlations, with hints that the loudest accretors may have lower outflow efficiency.","pith_inferences":["If the radio free-free emission is a direct proportionality tracer of accretion-driven outflow, then centimeter monitoring of individual disks could give a time-resolved view of accretion bursts and their ejection response, which the current snapshot sample does not provide.","The hint that transition disks show shallower accretion-outflow correlations may reflect the cavity itself impeding efficient angular momentum extraction; testing this requires a transition-disk sample that is not biased toward high accretors and, ideally, resolved radio maps of the inner cavity region.","The correlation coefficients are computed on heterogeneous literature data with different measurement depths and upper limits, so the reported r-values likely change when a uniform, deeper radio survey of both disk classes is obtained."],"forward_implications":["The free-free emission in full disks, not just transition disks, is a practical radio tracer of accretion-driven outflows in evolved Class II sources.","Accretion onto the central star in both full and transition disks appears to be coupled to a jet or MHD wind, supporting the idea that MHD winds, not purely viscous processes, drive disk accretion.","The ratio of ionized mass loss to accretion is small for strong accretors and larger for weak accretors, pointing to a mass-dependent or accretion-rate-dependent efficiency of outflow launching.","The absence of correlation with [O I] line luminosity indicates free-free emission and optical forbidden lines trace different parts or phases of the outflow, so combining both tracers gives a fuller picture of the jet-wind structure."],"supporting_citations":[{"why":"Provides the multi-wavelength Taurus photometry, the free-free flux measurements, and the spectral index analysis for the full-disk sample.","marker":"Garufi et al. (2025)"},{"why":"Supplies the transition-disk free-free fluxes, spectral indices, upper limits, and the previous finding that free-free emission in transition disks is likely an MHD wind or jet.","marker":"Rota et al. (2024)"},{"why":"Defines the geometric jet model used to convert free-free flux into ionized mass loss rate.","marker":"Reynolds (1986)"},{"why":"Provides the expected free-free luminosity for X-ray and EUV driven photoevaporative winds against which the observed emission is compared.","marker":"Pascucci et al. (2012)"},{"why":"Establishes the interpretation of flat or positive free-free spectral indices as collimated ionized outflows and the jet velocity scaling with stellar mass.","marker":"Anglada et al. (2018)"}],"fun_headline_variants":["Radio glow links accretion to outflow in 31 young disks","Jets and winds revealed by radio emission in young stars","Accretion and outflow rates correlate in protoplanetary disks","Radio signals trace jet-like winds in young stellar disks","Disk outflow tied to accretion in 31 Class II young stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conversion from observed radio flux to ionized mass loss rate assumes a fully ionized Reynolds jet geometry and an ionization that is total, so if the radio emission were instead dominated by gyro-synchrotron radiation or a photoevaporative wind, the inferred mass loss rates and the correlation would not support the MHD wind or jet conclusion.","fun_headline_variants_meta":{"raw":{"variants":["Radio glow links accretion to outflow in 31 young disks","Jets and winds revealed by radio emission in young stars","Accretion and outflow rates correlate in protoplanetary disks","Radio signals trace jet-like winds in young stellar disks","Disk outflow tied to accretion in 31 Class II young stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000176,"raw_usage":{"total_tokens":1362,"prompt_tokens":1088,"completion_tokens":274,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":704,"completion_tokens_details":{"reasoning_tokens":192}},"tokens_in":704,"tokens_out":274,"duration_ms":2587,"temperature":1.0,"reasoning_tokens":192,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:57:22.967409+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct measurement that would settle the claim is a resolved image of the radio continuum around several of the full disks: if the centimeter emission is spread over a broad, roughly spherical region whose size scales with the X-ray or EUV luminosity rather than being collimated along the disk axis, the jet interpretation would fail. A second test is to observe a disk while its accretion rate drops dramatically; if the centimeter flux does not fall proportionally, the free-free emission is not tracing the accretion-driven outflow.","supporting_citations":[{"cited_title":"The centimeter emission from planet-forming disks in Taurus","cited_arxiv_id":"2501.11686","evidence_quote":"Provides the multi-wavelength Taurus photometry, the free-free flux measurements, and the spectral index analysis for the full-disk sample."},{"cited_title":"A., Meijerhof, J","cited_arxiv_id":null,"evidence_quote":"Supplies the transition-disk free-free fluxes, spectral indices, upper limits, and the previous finding that free-free emission in transition disks is likely an MHD wind or jet."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the geometric jet model used to convert free-free flux into ionized mass loss rate."},{"cited_title":"2012, ApJ, 751, L42","cited_arxiv_id":null,"evidence_quote":"Provides the expected free-free luminosity for X-ray and EUV driven photoevaporative winds against which the observed emission is compared."},{"cited_title":"F., & Carrasco-González, C","cited_arxiv_id":null,"evidence_quote":"Establishes the interpretation of flat or positive free-free spectral indices as collimated ionized outflows and the jet velocity scaling with stellar mass."}],"review_version":1}