{"id":"0ab1fd80-15a6-4447-98a7-5425e90b75d5","arxiv_id":"2506.19779","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Quiescent millimeter emission of the young M dwarf AD Leo is 2-7 times brighter than a 1D thermal chromospheric model predicts, indicating quasi-steady non-thermal emission, and an 18-second double-hump frequency-rising flare was detected.","lead":"Astronomers measured the millimeter light from the young star AD Leo and found it up to seven times brighter than models of its hot atmosphere predict. The extra brightness points to steady non-thermal processes, and the team also caught a rare fast double-peaked flare.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Quiescent non-thermal claim hinges on an unconstrained 1D thermal baseline at mm-forming heights.","rationale":"The paper makes a valuable observational contribution: a clean NOEMA detection of AD Leo's quiescent mm spectrum at three frequencies and a well-characterized double-hump flare, with data products released on Zenodo. The central interpretive claim, however, is the quiescent non-thermal excess. That claim is only as strong as the thermal baseline. The 1D PHOENIX model is constrained by optical/near-UV lines that form lower than the mm continuum, so the mm prediction is not directly tested by the same data that validate the high-energy part of the SED. The reader's CONDITIONAL verdict correctly captures this: the non-thermal interpretation is plausible but not uniquely required. I agree with the reader's weakest-assumption identification. The concrete test proposed here—a sensitivity grid of allowed chromospheric temperature structures—would directly determine whether any thermal model consistent with the UVES data can account for the observed mm flux. If it cannot, the non-thermal claim would be substantially strengthened; if it can, the abstract's certainty should be reduced. Since the reader already conditioned the verdict on this model dependence, no further adjustment is needed.","tokens_in":16154,"tokens_out":12612,"duration_ms":134255,"concrete_test":"Generate a grid of 1D PHOENIX chromospheric models with modified temperature structures at column masses above the Balmer/Fe I formation region (e.g., varying T by ±0.2-0.5 dex for log m < -4 g cm^-2), each refit to the UVES line fluxes within their measurement errors, and compute the 84-102 GHz continuum for each. If any member of this grid reproduces the observed quiescent S(ν) within 1σ of Table 2, the non-thermal excess is not required by the current data. If instead all grid models that fit the line data underpredict the observed flux by >5σ, the non-thermal interpretation is substantially supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the quiescent mm-S(ν) of AD Leo contains a quasi-steady non-thermal component rests on the accuracy of the 1D PHOENIX chromospheric model (Sec. 3.2.2) as the thermal baseline. The temperature structure of this model was tuned to reproduce Balmer and Fe I lines in the UVES spectrum, which form in the middle/lower chromosphere. However, the 84-102 GHz free-free continuum is expected to become optically thick at much higher layers (small column mass, upper chromosphere/transition region) that these lines do not directly constrain. The observed excess grows toward lower frequencies (Fig. 3b inset), i.e., toward higher formation heights, exactly where the model is unconstrained. The authors' counterargument—that scaling the model by an 11% filling factor of 7 MK active regions (Sec. 3.2.2) does not remove the excess—is not a decisive test: that scaling is an ad hoc multiplication of the whole SED, not a self-consistent radiative-transfer calculation with a physically motivated inhomogeneous temperature/density structure. A modest increase in the temperature and density of the upper chromosphere in active regions could plausibly raise the thermal mm continuum by the needed factor 2-7 without violating the Balmer/Fe I constraints. The paper itself concedes this possibility in Sec. 3.2.2: 'This persistent ΔS/S_ch may imply the limitations of 1D models, that cannot factor in the surface inhomogeneities.' Thus the data do not uniquely force the non-thermal interpretation; the abstract states the conclusion more firmly than the modeling permits.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports NOEMA Band-1 (82–106 GHz) observations of the young M3.5V dwarf AD Leo, deriving a quiescent millimeter SED with detections at 84.3, 88.2, and 101.7 GHz. The authors measure a steep spectral index δ = −2.7 ± 1.3 and a brightness-temperature spectral index α_mm that is about three times larger than the extrapolated I-branch scaling law. Comparing the SED with a purely photospheric PHOENIX model and with a UVES-constrained 1D chromospheric model, they find a factor 2–7 excess over the chromospheric model, which they interpret as quasi-steady non-thermal emission powered by supra-thermal electrons. They also report an ~18 s double-hump flare with second-scale spectral-index variability and a frequency-rising behavior above 50% of peak flux, interpreted as multiple injections of accelerated electrons.","tokens_in":16381,"tokens_out":5445,"duration_ms":54816,"significance":"If the non-thermal interpretation is correct, the paper would provide the first evidence that young active M dwarfs host a quasi-steady non-thermal millimeter component, challenging the usual assumption that quiescent cool-star mm emission is purely thermal chromospheric emission. The flare observation is also a rare, well-characterized second-scale mm flare with a frequency-rising spectrum. The data analysis is careful in several respects: explicit imaging reliability criteria, multi-epoch stability checks for quiescence, a high-significance flare detection, and confirmation of variability by phase-only self-calibration. The main risk is that the non-thermal quiescent claim rests on a 1D thermal baseline that is not directly constrained at the heights where the millimeter continuum forms; this needs to be addressed before the central conclusion can be accepted.","major_comments":[{"comment":"The inference of quasi-steady non-thermal quiescent emission is not uniquely forced by the data. The thermal baseline S_ch^mod is tied to a 1D temperature structure constrained by Balmer and Fe I lines in UVES data, which form in the lower/middle chromosphere, whereas the 84–102 GHz free-free continuum becomes optically thick in the upper chromosphere/transition region, a height range not directly constrained by those lines. The residual ΔS/S_ch^mod increases toward lower frequencies (Fig. 3b inset), toward the unconstrained layers. The 11% filling-factor test with 7 MK active regions appears to be a uniform rescaling of the model SED rather than a self-consistent radiative-transfer calculation with active-region temperature and density stratification, so it cannot rule out a thermal explanation. The paper's own statement in §3.2.2, \"This persistent ΔS/S_ch may imply the limitations of 1D models, that cannot factor in the surface inhomogeneities,\" is in tension with the abstract's conclusion \"This indicates a quasi-steady non-thermal emission.\" Please either (a) compute a physically motivated active-region contribution and show that it cannot reproduce the excess, or (b) soften the conclusion to an excess over the 1D model and identify observables (polarization, spectral shape, variability) that would test the non-thermal origin.","section":"§3.2.2 and §4.1"},{"comment":"The derived spectral index δ = −2.7 ± 1.3 and the subsequent α_mm comparison are sensitive to the weakest detection: the 101.7 GHz point is only a 5σ detection. A robustness test excluding this point should be reported; without it, the reader cannot judge whether the steep spectrum and the factor-3 α_mm deviation are driven by one marginal measurement. In addition, the quoted errors appear to be thermal noise only; systematic uncertainties from amplitude calibration should be quantified and propagated into δ and α_mm.","section":"§3.1, Table 2, Fig. 2a"},{"comment":"The claim that AD Leo deviates by a factor of 3 from the α_mm–T_eff scaling law relies on extrapolating a relation calibrated on F–K dwarfs from Mohan et al. (2022) down to T_eff ≈ 3500 K, with UV Ceti as a single anchor at 2728 K. This is a long extrapolation and is not an independent test of non-thermal emission. Please quantify the extrapolation uncertainty or explicitly rephrase the result as a deviation from an extrapolated I-branch trend rather than from an established scaling law.","section":"§3.1, Fig. 2b"}],"minor_comments":[{"comment":"The sentence \"The flare resemble certain solar events\" should read \"The flare resembles certain solar events.\"","section":"Abstract"},{"comment":"The header \"T able 1\" has a stray space and should be \"Table 1.\"","section":"Table 1"},{"comment":"The range written as \"10 3 - 10 6 GHz\" should be formatted as 10^3–10^6 GHz.","section":"§2.1"},{"comment":"The downward arrow denoting the 104 GHz upper limit is not explained in the caption or the text; please clarify whether this is a non-detection from the upper sideband.","section":"Fig. 2a"},{"comment":"Please add explicit error bars or shaded uncertainty bands to the spectral-index time series in panel (c), and specify the time averaging used for the displayed points.","section":"Fig. 4"},{"comment":"The phrase \"strongly suggest the need to incorporate non-thermal emission mechanisms\" is stronger than the model-comparison evidence presented; consider aligning it with the revised, more cautious conclusion.","section":"§4.1"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the flare detection and the multi-epoch quiescent measurements are valuable and appear technically sound. My main reservation concerns the interpretation of the quiescent excess: the 1D thermal baseline is not constrained at the mm-forming heights, so the non-thermal claim needs either a stronger thermal-model test or a softened conclusion. This is fixable within the scope of the manuscript, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nRead the AD Leo NOEMA paper. The observations are the real thing: first quiescent mm SED for a young (<1 Gyr) M dwarf, three clean detections with careful imaging and reliability cuts, plus an 18 s double-hump flare with frequency-rising behavior that hasn't been seen in dMe stars. The α_mm comparison with the old-star scaling law, and the UV Ceti sanity check, are nice touches. If you work on stellar activity or mm cool-star observations, this paper is worth your time.\n\nThe soft spot is the central interpretive claim. The quiescent excess over the 1D PHOENIX chromospheric model is taken as evidence for quasi-steady non-thermal emission. But the stress-test concern is real: the UVES lines that constrained the temperature structure form in the middle/lower chromosphere, while the 84–102 GHz free-free continuum becomes optically thick in the upper chromosphere/transition region. The excess grows toward lower frequencies—i.e., higher formation heights—exactly where the model is unconstrained. The paper's counterargument, scaling the whole SED by an 11% filling factor of 7 MK active regions, is an ad hoc multiplicative factor, not a self-consistent inhomogeneous radiative-transfer calculation. A modest temperature/density bump in active-region upper chromospheres could plausibly raise the mm continuum by the needed factor 2–7 without violating the Balmer/Fe I constraints. The paper itself concedes this in Sec. 3.2.2: 'This persistent ΔS/S_ch may imply the limitations of 1D models.' So the data do not uniquely force the non-thermal interpretation; the abstract states it more firmly than the modeling allows.\n\nThe flare is a single event and lacks full Stokes calibration, but the relative polarization variability argument is reasonable and the source-born nature is supported by independent images with phase-only self-cal. The δ=-2.7±1.3 quiescent index is steep but the 101.7 GHz point is only 5σ; still, the trend is consistent.\n\nOverall: a valuable observational contribution, with an interpretation that needs to be reframed as one of two viable explanations until better thermal models (3D or empirical upper-chromosphere constraints) are available. I'd send it to a serious referee, expecting a major revision to the abstract and discussion. For my own work, I'd cite the quiescent SED and flare morphology regardless of the non-thermal interpretation.","headline":"Solid first mm SED of a young M dwarf and a novel flare, but the quiescent non-thermal claim doesn't uniquely follow from the 1D thermal baseline.","tokens_in":17042,"tokens_out":2520,"would_cite":true,"duration_ms":25266,"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":"AD Leo's quiet radio glow is up to 7 times too bright to be thermal.","keywords":["AD Leo","M dwarfs","millimeter continuum","chromospheric activity","non-thermal emission","stellar flares","supra-thermal electrons","brightness temperature spectral index"],"falsifier":"Measure AD Leo's millimeter spectrum with full Stokes polarization while simultaneously constructing a 3D magnetohydrodynamic model that reproduces its observed surface magnetic field and hot active regions; if that model reproduces the observed $2$–$7\\times$ excess and $\\alpha_{mm}$ without any supra-thermal electrons, the non-thermal quiescent claim collapses, whereas detection of strong circular polarization or a coherent spectral feature would confirm it.","tokens_in":15879,"feed_emoji":"📡","tokens_out":8726,"duration_ms":81546,"temperature":0.7,"pith_summary":"The paper tries to establish that quiescent millimeter emission from the young M dwarf AD Leo is not purely thermal. Interferometric observations near 94 GHz show a flux density that is 2–7 times higher than the spectrum predicted by a 1D chromospheric model whose temperature structure is constrained by optical and ultraviolet spectroscopy. The authors read this persistent excess, together with a steep mm brightness-temperature spectral index, as evidence of a quasi-steady non-thermal component powered by supra-thermal electrons. They also report an 18-second double-hump flare whose flux density rises with frequency above half-peak, which they attribute to multiple episodes of electron acceleration. If the interpretation holds, the standard assumption that quiescent mm emission in cool stars is chromospheric thermal emission must be relaxed for young, active M dwarfs.","feed_headline":"Quiet radio glow of a young M dwarf is 7x too bright for heat alone","feed_subtitle":"The excess points to steady supra-thermal electrons, and an 18-second flare shows the same acceleration at work.","key_machinery":"The load-bearing object is the 1D chromospheric model of AD Leo, a recalculation of an earlier model in which the temperature structure is fixed by fitting optical-UV emission lines and then used to predict the thermal millimeter continuum. Around it sit two diagnostics: the mm brightness-temperature spectral index $\\alpha_{mm}$ defined by $T_B(\\nu)\\propto\\nu^{-\\alpha_{mm}}$, which measures the chromospheric heating gradient, and the flare spectral index $\\delta$ from the ratio of lower-sideband to upper-sideband fluxes. The excess $\\Delta S/S_{\\rm ch}^{\\rm mod}$ between observed and modeled flux is the quantity whose magnitude (up to a factor 7) carries the non-thermal argument.","core_discovery":"The central discovery claim is that the quiescent 82–106 GHz spectrum of AD Leo (a ~250 Myr, M3.5V dwarf) is in excess of the thermal baseline set by a recalibrated 1D chromospheric model: the observed flux densities (155, 123 and 94 µJy at 84.3, 88.2 and 101.7 GHz) exceed the model by factors of roughly 2–7, with a spectral index $\\delta = -2.7 \\pm 1.3$. The paper argues this excess, which grows toward lower frequencies, cannot be removed by plausible active-region filling factors and implies quasi-steady non-thermal emission from supra-thermal electrons. As a corollary, the mm brightness-temperature spectral index $\\alpha_{mm}$ ($T_B(\\nu)\\propto \\nu^{-\\alpha_{mm}}$) is about three times larger than the $\\alpha_{mm}$–$T_\\mathrm{eff}$ scaling law derived for older I-branch stars, while the older M dwarf UV Ceti fits that law. The same data set contains an $\\sim18$ s double-hump flare at 86.3 and 101.7 GHz that is frequency-rising above half peak flux, with second-scale variability in Stokes I spectral index, which the authors interpret as evidence of multiple injections of accelerated electrons.","pith_inferences":["If the non-thermal quiescent interpretation is right, the same supra-thermal electron population should reveal itself in full-Stokes circular polarization or in a spectral turnover at other bands; a simultaneous 34–230 GHz campaign would test this without needing new theory.","The paper's own caveat that 1D models cannot capture surface inhomogeneities implies a direct falsification path: a 3D radiative-MHD model with AD Leo's measured surface magnetic fields could shift the inferred thermal baseline enough to absorb the excess.","The C/I branch difference suggests $\\alpha_{mm}$ may be a proxy for stellar age or rotation, not just $T_\\mathrm{eff}$; observing intermediate-age stars like $\\epsilon$ Eridani at multiple epochs could map how the scaling law re-establishes itself.","The frequency-rising flare threshold (above 50% of peak) hints at an optically thick non-thermal source whose spectrum turns over between 100 and 150 GHz; targeted high-cadence 150/230 GHz follow-up would locate the turnover."],"forward_implications":["Quiescent mm emission is not a safe thermal chromospheric thermometer for young, rapidly rotating M dwarfs; thermal model fits can understate the required heating unless a non-thermal component is included.","The $\\alpha_{mm}$–$T_\\mathrm{eff}$ scaling from old I-branch stars is not universal; young C-branch stars like AD Leo can deviate by a factor of about 3.","Frequency-rising, second-scale, double-hump mm flares are a stellar phenomenon, not just solar: AD Leo joins AU Mic and Proxima Cen with comparable luminosity and duration, but with a frequency-rising signature not seen in those stars.","Multiple humps in flare light curves map to multiple injection episodes of accelerated electrons, so mm light curves can resolve particle acceleration timing even without spatially resolved imaging.","The SPT95 minutes-long dM flares are likely stronger cousins of the weaker, shorter AU Mic/Proxima Cen/AD Leo events, implying a continuous luminosity-duration relation across dM mm flares."],"supporting_citations":[{"why":"Provides the $\\alpha_{mm}$–$T_\\mathrm{eff}$ scaling law for old I-branch stars that AD Leo's spectral index is compared against.","marker":"(Mohan et al. 2022)"},{"why":"Supplies the AD Leo chromospheric model that the paper recalculates with updated NLTE line lists to produce the thermal baseline.","marker":"(Fuhrmeister et al. 2005)"},{"why":"Supplies the photospheric model grid used for the purely photospheric comparison SED.","marker":"(Husser et al. 2013)"},{"why":"Establishes the archival SED compilation method and the earlier result that quiescent mm emission of older F-K stars matches 1D thermal models.","marker":"(Mohan et al. 2021)"},{"why":"Frames mm emission as a chromospheric activity diagnostic connecting brightness temperature to atmospheric height.","marker":"(Wedemeyer et al. 2016)"},{"why":"Provides the UV Ceti quiescent flux densities at 34 and 98 GHz used to compute the comparison $\\alpha_{mm}\\sim2.7$.","marker":"(Plant et al. 2024)"},{"why":"Supplies the SPT95 sample of dM mm flares against which the AD Leo flare's luminosity and spectral index are compared.","marker":"(Tandoi et al. 2024)"},{"why":"Provides the AU Mic flare luminosity and duration baseline, establishing the comparable-but-weaker relation for AD Leo.","marker":"(MacGregor et al. 2020)"}],"fun_headline_variants":["Young M dwarf's quiet mm glow is 7x too bright for heat alone","Quiescent mm emission of young M dwarf exceeds thermal models 7x","Supra-thermal electrons power young M dwarf's extra mm glow","Steady non-thermal mm glow from young M dwarf, not just heat"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on the 1D chromospheric model's temperature structure, fixed by optical-UV lines, being the correct predictor of the thermal millimeter continuum; if surface inhomogeneities or active-region geometry make the true thermal spectrum brighter, the observed excess could be thermal rather than non-thermal.","fun_headline_variants_meta":{"raw":{"variants":["Young M dwarf's quiet mm glow is 7x too bright for heat alone","Quiescent mm emission of young M dwarf exceeds thermal models 7x","Supra-thermal electrons power young M dwarf's extra mm glow","Steady non-thermal mm glow from young M dwarf, not just heat"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000865,"raw_usage":{"total_tokens":3899,"prompt_tokens":1241,"completion_tokens":2658,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":857,"completion_tokens_details":{"reasoning_tokens":2576}},"tokens_in":857,"tokens_out":2658,"duration_ms":21121,"temperature":1.0,"reasoning_tokens":2576,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:24:31.269740+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure AD Leo's millimeter spectrum with full Stokes polarization while simultaneously constructing a 3D magnetohydrodynamic model that reproduces its observed surface magnetic field and hot active regions; if that model reproduces the observed $2$–$7\\times$ excess and $\\alpha_{mm}$ without any supra-thermal electrons, the non-thermal quiescent claim collapses, whereas detection of strong circular polarization or a coherent spectral feature would confirm it.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the AD Leo chromospheric model that the paper recalculates with updated NLTE line lists to produce the thermal baseline."},{"cited_title":"M., Osten , R","cited_arxiv_id":null,"evidence_quote":"Provides the AU Mic flare luminosity and duration baseline, establishing the comparable-but-weaker relation for AD Leo."}],"review_version":2}