{"id":"ec126e11-fe9a-489f-b68d-fc0c95cc630f","arxiv_id":"2505.09148","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"ASKAP detections of 48 magnetic massive stars, including 14 new, confirm the CBO radio-luminosity scaling relation with best-fit Lrad = 10^-8.0 L_CBO^0.87 and yield 9 MRP candidates.","lead":"This paper reports ASKAP radio detections of 48 magnetic massive stars, 14 of which were previously undetected at radio wavelengths, and identifies 9 new candidate 'Main-sequence Radio Pulse Emitter' stars. The expanded sample supports the centrifugal-breakout scaling law for radio luminosity, with a slightly lower efficiency than earlier estimates.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed slope (0.87) and lower efficiency may be artifacts of estimating Lrad for new stars from single-band ~1 GHz flux; §4.1.2 and Appendix 2 show 1 GHz/trapezoid estimates can be systematically low by ~0.47 dex, comparable to the fit scatter.","rationale":"The reader's weakest assumption identifies exactly the load-bearing issue: the trapezoidal spectral shape used to convert single-band ASKAP flux densities into integrated incoherent radio luminosities. My stress-test sharpens this into a concrete threat to the headline numbers: the fit mixes two estimators (wideband-based for old stars, 1 GHz-based for new stars), Figure 4 demonstrates a luminosity-dependent offset between them, and Appendix 2 shows the assumed shape can be wrong by 0.47 dex, equal to the fitted scatter. The qualitative conclusion that a CBO-related correlation exists is well supported by the expanded sample and by the independent support of previous studies; the paper is transparent about its assumptions and makes the data public. But the specific slope 0.87 and the claimed lower efficiency should remain provisional. Because the reader already assigned CONDITIONAL and this is the same concern, no verdict change is needed.","tokens_in":39255,"tokens_out":9249,"duration_ms":95778,"concrete_test":"Using the 23 stars common to §4.1.2, re-fit the Lrad–LCBO relation twice: once with all Lrad set to the Shultz et al. (2022) values and once with all Lrad set to the ASKAP-only values, then compare with the mixed-estimator fit of Fig. 6. If the slope moves by more than its 1σ interval, or if the ASKAP-only fit already gives 0.87, the claimed lower efficiency is an artifact of estimator inhomogeneity. As a complement, recompute Lrad for the newly detected stars with archival 5 GHz measurements to test the trapezoid assumption directly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is the best-fit Lrad = 10^-8.0 L_CBO^0.87 (§4.1.3), read as a slightly lower radio-production efficiency than Owocki et al. (2022). This claim is not secure because the Lrad values entering the fit are estimated with two different methods whose difference is luminosity-dependent. For stars already in Shultz et al. (2022), the paper mostly adopts their wideband/multi-frequency luminosities; for newly added stars, §4.1.1 computes Lrad from the minimum ASKAP flux density near 0.9 GHz assuming the trapezoidal spectrum (flat 0.9–30 GHz, zero outside 0.6–100 GHz). Figure 4 then shows that 1 GHz-based ASKAP luminosities fall below the Shultz et al. values by a factor that grows with Lrad, i.e. the ratio Lrad,ASKAP/Lrad,Shultz+2022 decreases as Lrad increases. Since Lrad and LCBO are positively correlated, this introduces a systematic, luminosity-dependent offset between the two subsets rather than random scatter. Appendix 2 independently quantifies the failure of the trapezoid: for HD 142184 the assumed-shape integration underestimates the directly measured wideband luminosity by 0.47 dex, essentially the same as the 'true uncertainty' of 0.47 dex returned by the MCMC fit in §4.1.3. Therefore the fitted slope 0.87 (1σ [0.80,0.94]) and the claimed lower efficiency could be produced, or at least significantly shifted, by the estimator mismatch; the data do not yet demonstrate a deviation from the original slope-unity relation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the VAST-MeMeS project, an ASKAP survey-based search for non-thermal radio emission from magnetic massive stars, using cross-matching of Selavy source catalogues and forced Stokes I/V photometry. It reports radio detections of 48 magnetic hot stars, 14 of which have no prior radio detection, and identifies nine Main-sequence Radio Pulse emitter candidates. Combining these with previously known radio-bright stars yields an expanded sample of 70 objects, from which the authors re-derive the correlation between incoherent radio luminosity Lrad and centrifugal-breakout luminosity LCBO, obtaining a best fit of log Lrad = -8.0 + 0.87 log LCBO. This is interpreted as a slightly lower radio-production efficiency than the slope-unity relation reported by Owocki et al. (2022). The paper also argues that incoherent radio spectra can extend to lower frequencies than previously assumed, and uses partial correlation coefficients to argue against a direct role of effective temperature in driving the Lrad-LCBO relation.","tokens_in":39544,"tokens_out":5707,"duration_ms":56316,"significance":"If the central quantitative result were secure, this would be an important contribution: it approximately doubles the sample of radio-bright magnetic hot stars, extends the tested range of LCBO, and provides new southern-hemisphere detections and MRP candidates that will be valuable for follow-up. The paper is careful in its catalogue construction, uses publicly available ASKAP data, reports a Monte-Carlo-based cross-match reliability, and explicitly quantifies several systematic limitations in the appendices. However, the headline claim of a slope below unity and a slightly lower efficiency is not yet robust, because the luminosities entering the fit are estimated with two different methods whose difference is luminosity-dependent. The paper's own Appendix 2 shows that the assumed spectral shape can underestimate Lrad by 0.47 dex for a well-observed star, an amount equal to the inferred intrinsic scatter; this makes the specific slope and intercept values in Section 4.1.3 vulnerable to estimator bias.","major_comments":[{"comment":"The central quantitative claim, the fitted Lrad = 10^-8.0 L_CBO^0.87 relation, is not secure because of a luminosity-dependent mismatch between the luminosity estimators for the two subsamples. For stars already in Shultz et al. (2022), Section 4.1.2 adopts their wideband/multi-frequency luminosities, while for newly added stars Section 4.1.1 computes Lrad from a single ASKAP flux density near 0.9 GHz under the assumed trapezoidal spectrum. Figure 4 shows that the ratio Lrad,ASKAP/Lrad,Shultz+2022 decreases with increasing Lrad, so the new low-luminosity and previously known high-luminosity stars enter the fit with systematically different biases. Appendix 2 independently shows that for HD 142184 the assumed trapezoid underestimates the directly measured wideband luminosity by 0.47 dex, which is the same size as the 'true uncertainty' returned by the MCMC in Section 4.1.3. The authors should demonstrate that the fitted slope and intercept are robust to a homogeneous luminosity estimation, for example by fitting only stars with multi-frequency spectral coverage, or by applying the single-band/trapezoid estimator to all stars, or by explicitly modeling the spectral-shape bias as a function of Lrad.","section":"4.1.2, 4.1.3, Appendix 2"},{"comment":"The error bars on Lrad in Table 5 and Figure 6 include only flux-density measurement uncertainties, as the table footnote states, and the fit is performed with is_weighted=False, meaning the quoted 1σ intervals on the slope [0.80, 0.94] and intercept [-8.1, -7.8] do not propagate the dominant systematic uncertainties: spectral shape, rotational-phase variability, and the estimator mismatch described above. The derived 'true uncertainty' of 0.47 dex should therefore be interpreted as absorbing these systematics rather than as a measurement of astrophysical scatter. A sensitivity analysis, such as repeating the fit with the Shultz et al. (2022) luminosities for all stars or with a broader set of spectral-shape priors, would be needed before the claimed deviation from the slope-unity relation can be taken as evidence for a lower radio-production efficiency.","section":"4.1.3, Table 5"}],"minor_comments":[{"comment":"The third summary item states the best-fit relation as Lrad = 10^-8.52 L_CBO^0.88, which is the Keszthelyi et al. (2024) relation before the conversion described in Section 4.1; the paper's own best fit is Lrad = 10^-8.0 L_CBO^0.87. This should be corrected to avoid confusing the paper's result with the earlier untranslated relation.","section":"Section 6 (Summary)"},{"comment":"The caption begins 'he stellar magnetospheric parameters', which is missing the leading 'T', and the use of '-1' for missing values could be misread as a physical placeholder; a clearer notation such as '...' or 'N/A' would be preferable.","section":"Table 5 caption"},{"comment":"The sentence noting that Lrad,ASKAP/Lrad,Shultz+2022 decreases with Lrad is an important observation, but it is not quantified; reporting a rank correlation coefficient and its uncertainty would strengthen the point and make the luminosity dependence more transparent.","section":"4.1.2"},{"comment":"Appendix 2 discusses the two stars with turn-over observed at both ends of the spectrum, but Figure A1 is not cross-referenced anywhere in the main text; adding an explicit reference near Section 4.1.1 or 5.1 would help readers connect the spectral-shape caveat to the luminosity estimates.","section":"Appendix 2 and Figure A1"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"One line: this is a genuinely useful observational sample paper, but the revised slope 0.87 is not secure because the two luminosity estimators used are not consistent, and the discrepancy is luminosity-dependent.\n\nWhat's actually new: fourteen first-time radio detections, nine MRP candidates (three already confirmed by follow-up), and the first sub-GHz flux density measurements for many stars. The cross-matching is careful, with Monte Carlo reliability radii, and the data are in CASDA; a competent group could reproduce the catalogue. The paper also does a sensible partial-correlation analysis showing LCBO, not Teff, drives the correlation.\n\nThe soft spot is the one the stress test flags. For stars already in Shultz et al. (2022) they adopt the wideband-based luminosities; for new stars they estimate Lrad from a single ~1 GHz flux density and a trapezoidal spectrum. Figure 4 shows the ASKAP-based luminosities fall below the Shultz values by a factor that grows with Lrad. That means the two subsets are not offset by a constant; the offset is correlated with the independent variable. Appendix 2 shows for HD 142184 the trapezoid integration underestimates the directly measured wideband luminosity by 0.47 dex—the same as the fit's intrinsic scatter. So the fitted slope 0.87, and the claim of slightly lower efficiency, could be produced by the estimator mismatch. The authors acknowledge that Lrad uncertainties exclude spectral shape and rotational phase, and they explicitly note the spectral index depends on luminosity, but they don't propagate that into the fit. The post-hoc exclusion of HD 101412 is also a bit quick; unlike HD 148937 (no CM), there is no physical justification, and including it would change the fit.\n\nNone of this kills the main result. The correlation between Lrad and LCBO is visually robust across 70 stars, and the sample expansion is a real contribution. But the precise slope and intercept should be treated as provisional. If I were the editor I would send this to a competent referee, and if I were the referee I would ask for a fit that includes a systematic error term or a reanalysis using only stars with wideband spectra. The paper is for exactly the community that needs more radio-bright magnetic massive stars: observers and magnetospheric theorists. It deserves publication after those revisions.","headline":"Genueninely useful sample expansion for radio-bright magnetic massive stars, but the revised CBO slope is not secure because the two L_rad estimators differ in a luminosity-dependent way.","tokens_in":40207,"tokens_out":3158,"would_cite":true,"duration_ms":31967,"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":"ASKAP data on 70 radio-bright magnetic hot stars confirm that radio luminosity scales with centrifugal-breakout luminosity at slightly lower efficiency than previously reported.","keywords":["magnetic massive stars","non-thermal radio emission","centrifugal breakout","ASKAP","radio luminosity","main-sequence radio pulse emitters","stellar magnetospheres","incoherent gyrosynchrotron emission"],"falsifier":"Observe a sample of the newly detected stars simultaneously from roughly 0.3 to 30 GHz; if their flux densities decline steeply below 1 GHz instead of staying flat, the trapezoid integration overestimates $L_{\\mathrm{rad}}$ and the reported slope of 0.87 and intercept of $-8.0$ would need revision.","tokens_in":38987,"feed_emoji":"📡","tokens_out":8654,"duration_ms":73395,"temperature":0.7,"pith_summary":"This paper uses survey data from the Australian SKA Pathfinder to grow the sample of radio-bright magnetic massive stars from 47 to 70, adding 14 stars never before detected in radio. With the expanded sample, the authors test the claim that the incoherent radio luminosity of a magnetic hot star is set by the power released in centrifugal breakouts, the episodic escapes of wind plasma trapped in the star's magnetosphere. They find the scaling still holds, with a best-fit relation $L_{\\mathrm{rad}} = 10^{-8.0} L_{\\mathrm{CBO}}^{0.87}$, i.e. a slightly lower radio-production efficiency than the earlier $L_{\\mathrm{rad}} = 10^{-8} L_{\\mathrm{CBO}}$. They also show that ASKAP detections at about 1 GHz imply incoherent spectra can stay flat to lower frequencies than previously assumed, and they identify nine candidate main-sequence radio pulse emitters. Because the same scaling links radio output to fundamental stellar parameters, a robust relation gives a calibrated way to infer magnetospheric physics from a single radio measurement.","feed_headline":"Radio power of 70 magnetic hot stars tracks centrifugal breakouts","feed_subtitle":"An expanded ASKAP sample upholds the centrifugal-breakout scaling law at slightly lower radio efficiency.","key_machinery":"The load-bearing quantity is the centrifugal-breakout luminosity $L_{\\mathrm{CBO}} = B_{\\mathrm{eq}}^2 R_*^4 \\Omega^2 / v_{\\mathrm{orb}}$, the power released when plasma trapped in the co-rotating magnetosphere periodically breaks open the field lines; the theory assumes the field behaves like a monopole at the reconnection site. The paper combines this with a prescription for turning sparse radio measurements into luminosities: a trapezoidal spectrum that is flat between about 0.9 GHz and 30 GHz and vanishes at 0.6 and 100 GHz, calibrated on the ten stars with wideband spectra. Cross-matching 761 known magnetic hot stars against ASKAP point-source catalogues at 98% reliability radii, and forced fitting in Stokes I and V images, supplies the flux densities; the fit itself is a Markov Chain Monte Carlo line fit in log-log space.","core_discovery":"The central claim is that the empirical correlation between incoherent radio luminosity $L_{\\mathrm{rad}}$ and centrifugal-breakout luminosity $L_{\\mathrm{CBO}}$ survives a roughly 50% expansion of the known radio-bright magnetic hot star population. Fitting the 68 usable stars (after excluding the non-CM star HD 148937 and the overluminous outlier HD 101412) gives $L_{\\mathrm{rad}} = 10^{-8.0} L_{\\mathrm{CBO}}^{0.87}$ with 1$\\sigma$ ranges [0.80, 0.94] on the slope and [$-8.1$, $-7.8$] on the intercept; the intercept agrees with the $L_{\\mathrm{rad}} = 10^{-8} L_{\\mathrm{CBO}}$ reported earlier, while the sub-unity slope indicates a slightly lower efficiency at the high-luminosity end. The paper further reports that for several stars the ~1 GHz ASKAP flux density exceeds what the previously assumed 1.5-30 GHz flat spectrum would predict, and that a full-spectrum integration for HD 142184 changes its luminosity by 0.47 dex relative to the earlier trapezoidal estimate. Partial-correlation analysis finds that $L_{\\mathrm{rad}}$ remains strongly correlated with $L_{\\mathrm{CBO}}$ after removing effective temperature, with no significant residual dependence on temperature.","pith_inferences":["If flat sub-GHz spectra are the norm rather than the exception, then single-epoch low-frequency surveys could eventually provide nearly unbiased radio luminosities for hundreds of magnetic hot stars, turning the current 70-star sample into a population-scale test of the CBO relation.","The sub-unity slope (0.87) implies the radio-production efficiency drops as $L_{\\mathrm{CBO}}$ grows; one plausible but unmodelled cause is stronger free-free absorption in the denser winds of the more luminous systems, which would flatten the observed relation artificially.","The HD 142184 case suggests that published $L_{\\mathrm{rad}}$ values built on the old trapezoid assumption may be systematically biased; if similar biases affect many stars, re-deriving the fit from true broadband spectra could change both the slope and the scatter of the scaling relation.","HD 101412, with its unusually steep inferred spectrum and overluminosity, is the clearest test case: simultaneous multi-frequency observations would show whether a cool star can genuinely violate the CBO scaling or whether rotational sampling created an artefact."],"forward_implications":["The CBO mechanism remains viable across OBA spectral types: the relation built on 47 stars still describes a sample of 70, extending the range of $L_{\\mathrm{CBO}}$ by an order of magnitude.","Incoherent radio spectra can stay flat below 1 GHz, so sub-GHz survey data can contribute directly to measuring radio luminosities rather than being dismissed as dominated by coherent emission.","Producing radio luminosity does not appear to depend on stellar temperature once $L_{\\mathrm{CBO}}$ is accounted for; the correlation between $L_{\\mathrm{rad}}$ and $T_{\\mathrm{eff}}$ seen in the plot is an indirect effect of the mass-radius-temperature relation.","Nine stars meet the circular-polarisation or variability criteria for main-sequence radio pulse emitters, widening the hunting ground for electron-cyclotron-maser emission.","Single-frequency luminosity estimates can be off by factors of two or more: five common stars show ASKAP luminosities at least double the earlier published values, arguing for wideband spectral campaigns."],"supporting_citations":[{"why":"Provides the centrifugal-breakout model and the original $L_{\\mathrm{rad}} = 10^{-8} L_{\\mathrm{CBO}}$ scaling that the expanded sample is tested against.","marker":"Owocki et al. (2022)"},{"why":"Supplies the prior 47-star sample, the trapezoidal spectral integration strategy, and the radio luminosities the new ASKAP estimates are compared to.","marker":"Shultz et al. (2022)"},{"why":"Established the empirical radio luminosity relation with stellar parameters and the ~1 GHz spectral turnover that the new data extend.","marker":"Leto et al. (2021)"},{"why":"Provides the most recent independent fit ($L_{\\mathrm{rad}} = 10^{-8.52} L_{\\mathrm{CBO}}^{0.88}$) and the 650 MHz detection of HD 55522, whose $L_{\\mathrm{CBO}}$ definition differs by a factor of four.","marker":"Keszthelyi et al. (2025)"},{"why":"Established centrifugal breakout as the plasma transport mechanism in centrifugal magnetospheres from H$\\alpha$ observations.","marker":"Shultz et al. (2020)"},{"why":"Supplies the Monte Carlo cross-matching method and ASKAP catalogue filtering used to match stars to radio sources.","marker":"Driessen et al. (2024)"},{"why":"Defines the minimum flux-density-gradient condition used as one of the MRP candidate criteria.","marker":"Das et al. (2022c)"},{"why":"Provides earlier ~1 GHz detections and candidate identifications whose high circular polarisation indicates coherent emission for HD 61556 and HD 105382.","marker":"Pritchard et al. (2021)"},{"why":"Reports follow-up ATCA observations confirming three of the MRP candidates identified here.","marker":"Das et al. (2025)"}],"fun_headline_variants":["Magnetic massive stars' radio output follows centrifugal breakouts","ASKAP finds 48 magnetic stars with radio glow, 14 new","Centrifugal-breakout radio scaling holds for 68 magnetic stars","Magnetic star radio emission less efficient but still matches breakouts"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The weakest link is the assumption that the newly detected stars have the same flat, trapezoidal incoherent radio spectrum as the ten well-observed stars, so that flux densities measured only near 1 GHz can stand in for the full 0.6-100 GHz luminosity; if the true spectra turn over below 1 GHz or have different shapes, the derived $L_{\\mathrm{rad}}$ values and the fitted $L_{\\mathrm{rad}}$-$L_{\\mathrm{CBO}}$ relation would be biased.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic massive stars' radio output follows centrifugal breakouts","ASKAP finds 48 magnetic stars with radio glow, 14 new","Centrifugal-breakout radio scaling holds for 68 magnetic stars","Magnetic star radio emission less efficient but still matches breakouts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000816,"raw_usage":{"total_tokens":3681,"prompt_tokens":1154,"completion_tokens":2527,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":770,"completion_tokens_details":{"reasoning_tokens":2453}},"tokens_in":770,"tokens_out":2527,"duration_ms":23379,"temperature":1.0,"reasoning_tokens":2453,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:39:06.479462+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe a sample of the newly detected stars simultaneously from roughly 0.3 to 30 GHz; if their flux densities decline steeply below 1 GHz instead of staying flat, the trapezoid integration overestimates $L_{\\mathrm{rad}}$ and the reported slope of 0.87 and intercept of $-8.0$ would need revision.","supporting_citations":[{"cited_title":"E., Owocki , S., Rivinius , T., et al","cited_arxiv_id":null,"evidence_quote":"Established centrifugal breakout as the plasma transport mechanism in centrifugal magnetospheres from H$\\alpha$ observations."},{"cited_title":"2021, , 502, 5438","cited_arxiv_id":null,"evidence_quote":"Provides earlier ~1 GHz detections and candidate identifications whose high circular polarisation indicates coherent emission for HD 61556 and HD 105382."}],"review_version":1}