{"id":"53c22f1f-5fd5-4eac-86df-20720af60bc0","arxiv_id":"2608.13330","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"ASKAP monitoring at 887.5 MHz found three classical novae whose steep early radio flares are better matched by a broken-power-law synchrotron model than by thermal free-free emission.","lead":"ASKAP's routine radio monitoring of the Galactic plane caught three classical novae with steep, early radio flares that fit shock-powered synchrotron models better than thermal emission models, and all three were also detected in gamma rays. A smart generalist might read this as evidence that low-frequency survey monitoring can identify shock-dominated novae and connect radio to gamma-ray activity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Dominant-synchrotron claim not settled by single-frequency light-curve shape: hand-fixed broken power-law indices and poor reduced chi2 (2.9–37.8) make the model-selection argument inconclusive without alternative density models.","rationale":"The paper makes a genuinely useful observational contribution: a homogeneous, well-sampled, single-frequency radio sample of three gamma-ray-detected novae, with a statistically significant association (p ≈ 0.004) between radio and Fermi-LAT detection. That part of the paper is robust. The concern is specifically the abstract's dominant-synchrotron claim. It is supported by (a) the model comparison, which is loaded by hand-fixed k, m, n; (b) poor absolute fits, especially V1723 Sco; and (c) the absence of a clumpy thermal free-free comparison. However, the independent evidence—Fermi-LAT detections for all three, VLA spectral index α = 0.5 for V1723 Sco, and luminous radio emission ≥ 10^20 erg/s/Hz—is genuinely consistent with synchrotron emission, so the paper should not be rejected. The correct verdict is CONDITIONAL: the central claim should be softened or the alternative models tested. The reader's weakest_assumption exactly identified the hand-fixed indices and the loaded model-selection, and my analysis agrees. The concrete test I propose would settle whether the decline slope really requires the broken power-law model or whether a clumpy thermal model fits equally well; that is the decisive check.","tokens_in":18434,"tokens_out":7716,"duration_ms":59151,"concrete_test":"Perform a model-selection test treating k, m, n as free parameters (recompute AIC/BIC with the full 8-parameter broken power-law model) and add a clumpy/aspherical thermal free-free model with comparable parameter count. If the broken power-law synchrotron model no longer yields a decisive AIC/BIC improvement (e.g., Δ < 10 against a clumpy thermal model, or after parameter-count correction), then the 'dominant synchrotron' claim should be softened.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim—that all three novae are dominated by synchrotron emission at 887.5 MHz—rests on the modified synchrotron model beating the thermal and single power-law models in AIC/BIC. But the decisive feature of that model is a broken power-law density profile (Eq. 5) whose indices k, m, n are not fit; they are fixed by hand from the observed light-curve slopes (Section 3.2, Table 6). Thus the comparison is partly circular: the model is tuned to the very shape it is supposed to explain. The absolute fit quality is also poor: reduced chi2 are 2.9, 3.7, and 37.8 for V6598 Sgr, V1716 Sco, and V1723 Sco, respectively, and the paper states 'none of the models has a formally acceptable fit' and that V1723 Sco 'has no acceptable fit with any of the models.' The steep decline (S_nu ~ t^-3.3 to t^-3.6) is the key discriminator, but the paper never fits a clumpy or aspherical thermal free-free model—which can also produce steep, early radio light curves in novae (V1723 Aql, V959 Mon)—so the dominant-mechanism conclusion is not established by the light-curve comparison alone. Independent Fermi-LAT and X-ray evidence supports shock activity and makes synchrotron plausible, but does not quantify the radio emission fraction at 887.5 MHz.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript searches the ASKAP VAST 887.5 MHz survey for radio emission from classical novae, cross-matching 43 optically discovered novae erupting between 2021 September and 2025 November within the Galactic footprint and identifying three with significant radio emission: V6598 Sgr, V1716 Sco, and V1723 Sco. The authors fit the radio light curves with three models: a thermal free-free model, a standard synchrotron model with a wind-like density profile, and a modified synchrotron model incorporating a broken power-law density profile (Eq. 5). They report that the modified synchrotron model is preferred by AIC/BIC for all three novae, that the broken power-law density profile fits better than a standard wind profile, and that all three radio-detected novae are among the six Fermi-LAT-detected novae in the sample, with a Fisher-type p≈0.004 for the association. The paper concludes that single-frequency light-curve shape can identify shock-driven synchrotron-dominated novae.","tokens_in":18906,"tokens_out":5220,"duration_ms":47537,"significance":"If the central claim holds, the paper provides a method to identify non-thermal radio emission from classical novae using single-frequency light-curve morphology rather than multi-frequency spectra or brightness-temperature measurements, which would be valuable for interpretation of current and future synoptic radio surveys such as SKA. The measured 887.5 MHz light curves of three gamma-ray-detected novae are a useful data product, and the statistical association between ASKAP detections and Fermi-LAT detections, though based on small numbers, is suggestive. The paper honestly reports the poor reduced-chi2 values and parameter degeneracies. However, the central claim that all three novae are dominated by synchrotron emission is not established by the analysis as presented; the broken power-law indices are fixed from the data and no clumpy or aspherical thermal alternative is tested.","major_comments":[{"comment":"The broken power-law indices k, m, n are not fitted; they are chosen by hand from the observed rise and decay slopes, as stated in the text ('we fixed them at certain values depending on the steepness of temporal evolution'). Because the modified synchrotron model's superior AIC/BIC is largely attributable to these hand-tuned slopes, the model-selection comparison is partially circular: the model is engineered to reproduce the very light-curve shape it is then preferred for. The analysis would need to treat k, m, n as free parameters, or marginalize over a physically motivated prior, to support the claim that the broken power-law density profile is genuinely required.","section":"Section 3.2, Eq. (5), Table 6"},{"comment":"The reduced chi2 values are 2.9 (V6598 Sgr), 3.7 (V1716 Sco), and 37.8 (V1723 Sco), and the text states that 'none of the models has a formally acceptable fit' and that V1723 Sco 'has no acceptable fit with any of the models.' Despite this, the abstract asserts that all three novae 'show evidence of non-thermal synchrotron emission as the dominant emission mechanism at this frequency.' This is an overstatement; for V1723 Sco the paper explicitly does not identify a preferred model, so the abstract should be revised to reflect the uncertainty.","section":"Table 7 and Section 4.3"},{"comment":"The thermal comparison model is a single, spherically symmetric, uniformly filled shell. Novae with aspherical or clumpy thermal ejecta can produce steep, early single-frequency radio light curves (e.g., V1723 Aql, V959 Mon), and the paper does not fit any such thermal model. Without comparing against a clumpy or aspherical thermal free-free model, the steep observed decline S_nu ~ t^{-3.3} to t^{-3.6} does not uniquely require synchrotron emission, so the dominant-mechanism conclusion is not established by light-curve shape alone.","section":"Section 3.2 and Section 5"}],"minor_comments":[{"comment":"The text states that V6598 Sgr radio observations ran from 2022 November 14, but Table 1 lists observations beginning 2023 August 4; similarly V1716 Sco is said to be observed from 2022 November 19, while Table 1 starts 2023 May 21. Please reconcile the dates, perhaps by clarifying whether the earlier dates refer to the start of the VAST survey monitoring campaign rather than the first detection.","section":"Section 2"},{"comment":"The p-value is computed from a 2x2 contingency table, but the comparison is post-hoc and the 95% confidence interval for the detection fraction in gamma-ray novae (12% to 88%) is very broad; the text should state more explicitly that the association is tentative rather than definitive.","section":"Section 5, Figure 10"},{"comment":"The fixed values of epsilon_e = epsilon_B = 0.01 and p = 2.5 are standard but arbitrary; since the fitted physical parameters scale strongly with these choices, please briefly justify or discuss their influence on the derived energetics.","section":"Section 3.2 and Table 6"},{"comment":"The abstract and conclusion should be harmonized: the conclusion says V1723 Sco is 'not well described by either model,' while the abstract claims all three novae are dominated by synchrotron emission; please adjust both to match the stated uncertainty.","section":"Abstract and Section 6"}],"recommendation":"major_revision","confidential_remarks":"The manuscript presents a potentially useful dataset and a reasonable exploratory analysis, but the central claim is not supported as stated. The model-selection step is partly circular because the broken power-law indices are fixed from the observed light-curve slopes, and the absolute fit quality is formally unacceptable for all sources, especially V1723 Sco. The abstract overstates the conclusion. I recommend major revision with the following specific requests: (1) treat k, m, n as free parameters or explicitly demonstrate robustness over a plausible range; (2) fit at least one aspherical or clumpy thermal free-free model as a control; (3) revise the abstract so that it does not claim a dominant mechanism for V1723 Sco; (4) discuss the small-number gamma-ray association more cautiously. The paper is not beyond repair, but the load-bearing evidence needs strengthening."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper gives you three well-sampled single-frequency nova light curves from the ASKAP VAST survey and a clean statistical statement that radio detections cluster with Fermi-LAT detections. The modeling is where I'd pump the brakes. The abstract says all three novae show synchrotron as the dominant mechanism, but the reduced chi2 values are bad (2.9, 3.7, 37.8), and V1723 Sco explicitly has no acceptable fit. For the two better sources, the 'modified synchrotron' model wins AIC/BIC partly because the broken-power-law indices k, m, n are fixed by hand from the observed rise and decay slopes before fitting. That doesn't make the conclusion wrong, but it is not an independent test.\n\nWhat is genuinely new: this is the first systematic single-frequency modeling of a homogeneous survey sample of novae, and it demonstrates that steep, fast radio flares are common enough to be caught at 887.5 MHz with 2-week cadence. The gamma-ray association, 3/6 vs 0/26, with p~0.004, is the strongest quantitative result in the paper, though the binomial uncertainty is large (95% CI 12-88% for the detection fraction among gamma-ray novae). I also credit the authors for including distance scaling tables and for being upfront that none of the models has a formally acceptable fit. The corner plots show they know the parameters are degenerate.\n\nWeak spots: the dominant-mechanism claim is load-bearing and not supported by absolute fit quality. The paper never fits a clumpy or aspherical thermal free-free model, which can produce steep early light curves in novae (V1723 Aql, V959 Mon). So the light-curve shape alone does not rule out thermal emission. And the gamma-ray/radio association, while suggestive, is a 2x2 table with three radio detections; the p-value is fragile to one event. Neither of these should sink the paper, but the abstract should be softened to 'consistent with' or 'suggestive of' rather than 'show evidence of dominant.' The current phrasing overstates what single-frequency light-curve fitting can establish.\n\nWho this is for: observers working on nova radio surveys, and theorists who care about shock physics in novae. It deserves serious peer review—the data release and the statistical association are useful, and the model-selection issue is a good forcing function for the field to think about what single-frequency curves can actually deliver. I would recommend acceptance after major revision: soften the claim, add or at least discuss alternative thermal morphologies, and make the choice of k, m, n more transparent (e.g., test a grid or treat them parametrically). The observational core is solid.","headline":"A useful survey sample and a suggestive radio/gamma-ray link, but the 'dominant synchrotron' claim outruns the fits because the key model features are hand-set to match the light curves.","tokens_in":19405,"tokens_out":2254,"would_cite":true,"duration_ms":22320,"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":"Three classical novae detected at 887.5 MHz show non-thermal synchrotron emission as their dominant radio mechanism, argued from single-frequency light-curve fits.","keywords":["classical novae","non-thermal synchrotron emission","free-free thermal emission","radio light curves","ASKAP VAST survey","shock-driven particle acceleration","gamma-ray novae","broken power-law density profile"],"falsifier":"Observe any of the three novae at two or more radio frequencies through the rise, peak, and decline; if the spectral index is consistent with free-free emission (roughly $-0.1$ to $+2$) rather than synchrotron with self-absorption or optically thin synchrotron (roughly $-0.5$ or steeper), the dominant-mechanism claim fails. A clumpy thermal model that reproduces the steep light curves without any synchrotron component would also falsify the conclusion.","tokens_in":18261,"feed_emoji":"📡","tokens_out":9305,"duration_ms":81795,"temperature":0.7,"pith_summary":"The paper seeks to establish that radio light curves at a single frequency, 887.5 MHz, can identify the dominant emission mechanism in classical novae without multi-frequency spectra. Cross-matching 43 optically discovered novae with ASKAP VAST survey data, the authors find three significant radio sources: V6598 Sgr, V1716 Sco, and V1723 Sco. Fitting thermal free-free, standard synchrotron, and a synchrotron model with a broken power-law density profile, they conclude that non-thermal synchrotron emission dominates all three and that a broken power-law density profile fits better than a standard wind profile. The three radio-detected novae are exactly the gamma-ray-detected ones in the sample, so the paper links radio synchrotron emission to shock-driven particle acceleration and GeV gamma-ray production. If correct, this makes wide-field single-frequency surveys a viable route to discovering and physically classifying shock-driven novae.","feed_headline":"Three novae show synchrotron shocks at one radio frequency","feed_subtitle":"One 887.5 MHz survey frequency separates shock-driven novae from thermal ones and links radio to gamma rays.","key_machinery":"The load-bearing object is a broken power-law radial density profile for the circumbinary material that the nova shock encounters, written as $\\rho \\propto (R/R_0)^{-k} \\left(\\frac{1}{2}\\left[1 + (R/R_0)^{1/n}\\right]\\right)^{(k-m)n}$, where $k$ sets the inner slope, $m$ the outer slope, and $n$ the smoothness of the transition at the break radius $R_0$. The paper fixes $k$, $m$, and $n$ by hand from the observed rise and decay slopes of each light curve rather than fitting them, then runs MCMC over five physical parameters: filling factor $f$, wind velocity $V_{\\rm wind}$, mass-loss rate $\\dot{M}$, explosion energy $E$, and ejecta mass $M_{\\rm ej}$. This density profile is what allows the model to reproduce the steep rise and rapid post-peak fall of these flares, which a single power-law profile cannot reproduce.","core_discovery":"The central claim is that non-thermal synchrotron radiation from shock-accelerated electrons dominates the 887.5 MHz radio emission of V6598 Sgr, V1716 Sco, and V1723 Sco, and that this can be inferred from the shape of a single-frequency light curve when interpreted with a physically motivated model. The thermal free-free model under-predicts the observed flux and rises too slowly; the standard synchrotron model over-predicts and cannot reproduce the sharp, short-lived flares; the modified synchrotron model with a broken power-law density profile yields the lowest reduced chi-square, AIC, and BIC for V6598 Sgr and V1716 Sco and better captures the steep evolution of V1723 Sco, although none of the fits is formally acceptable for that source. The paper also reports that all three radio-detected novae lie within the six Fermi-LAT gamma-ray-detected novae in the survey footprint, a coincidence with $p \\approx 0.004$, supporting the same shock population producing both radio and gamma-ray emission.","pith_inferences":["If single-frequency light-curve shape reliably identifies synchrotron-dominated novae, then archival survey data could be re-mined to build a larger radio-selected nova population, independent of gamma-ray telescope exposure.","The paper does not fit the density-profile indices $k$, $m$, and $n$; they are set by hand from the same light curves they are used to explain, so a clumpy or aspherical thermal ejecta model with a variable filling factor might also reproduce the steep light curves and would be a direct challenge to the synchrotron conclusion.","With sharper distance measurements for these systems, the degenerate fitted parameters (mass-loss rate, explosion energy, ejecta mass) would become physically informative, potentially separating the three novae into distinct shock regimes.","A joint radio-gamma analysis of a larger nova sample could test whether the $p \\approx 0.004$ association persists; if it does, low-frequency radio surveys could serve as an unbiased finder for shock-powered novae."],"forward_implications":["A single 887.5 MHz light curve, modeled with a broken power-law density profile, can separate shock-dominated novae from thermal ones, so future wide-field surveys can classify novae without multi-frequency follow-up.","Because all three ASKAP-detected novae are among the six Fermi-LAT gamma-ray novae in the footprint, radio detection at this frequency selects the same shock-powered population that produces GeV gamma rays, and the radio detection fraction among gamma-ray novae is 50 percent versus zero otherwise.","The fitted radio luminosities, around or above $10^{20}\\,\\mathrm{erg\\,s^{-1}\\,Hz^{-1}}$, place the three novae among the most luminous synchrotron-emitting novae known, implying efficient particle acceleration.","The steep post-peak decays, with flux falling as roughly $t^{-3.3}$ to $t^{-3.6}$, are much faster than the $t^{-2}$ expected for freely expanding thermal ejecta, indicating that the shock encounters a rapidly thinning circumbinary medium.","The success of the model with regular two-week monitoring provides a template for future radio surveys, such as the Square Kilometre Array, where targeted follow-up of large source samples will not be feasible."],"supporting_citations":[{"why":"provides the non-thermal synchrotron shock emission model at the base of the fits.","marker":"Chevalier (1998)"},{"why":"supplies the standard synchrotron model and parameter framework that the paper modifies with a broken power-law density.","marker":"Nyamai et al. (2023)"},{"why":"supplies the thermal free-free emission model used as the competing baseline.","marker":"Hjellming et al. (1979)"},{"why":"describes the internal and external shock geometry and absorption that motivate the synchrotron interpretation.","marker":"Vlasov et al. (2016)"},{"why":"is the earlier VAST-based nova study whose methods and comparison sample ground the survey analysis.","marker":"Gulati et al. (2023)"},{"why":"is the multi-frequency radio census of classical novae that defines the luminous synchrotron comparison class.","marker":"Chomiuk et al. (2021b)"},{"why":"establishes Fermi-LAT gamma-ray detections of classical novae, the baseline for the radio-gamma association.","marker":"Ackermann et al. (2014)"},{"why":"supports shock activity and gamma-ray behavior for V1723 Sco.","marker":"Fauverge et al. (2026)"},{"why":"provides the adopted distances for the three novae used in flux scaling and parameter estimation.","marker":"Schaefer (2025)"}],"fun_headline_variants":["Single-frequency radio reveals shock-driven novae","Novae show synchrotron shocks at 887.5 MHz","Three novae's radio flares point to shocks","Radio light curves distinguish nova emission mechanisms"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the light-curve shape alone, fitted with a broken power-law density model whose slope indices are chosen by hand, can attribute the dominant emission mechanism to non-thermal synchrotron rather than thermal free-free or a different density structure.","fun_headline_variants_meta":{"raw":{"variants":["Single-frequency radio reveals shock-driven novae","Novae show synchrotron shocks at 887.5 MHz","Three novae's radio flares point to shocks","Radio light curves distinguish nova emission mechanisms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000198,"raw_usage":{"total_tokens":1423,"prompt_tokens":1052,"completion_tokens":371,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":668,"completion_tokens_details":{"reasoning_tokens":310}},"tokens_in":668,"tokens_out":371,"duration_ms":3985,"temperature":1.0,"reasoning_tokens":310,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:30:57.186045+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe any of the three novae at two or more radio frequencies through the rise, peak, and decline; if the spectral index is consistent with free-free emission (roughly $-0.1$ to $+2$) rather than synchrotron with self-absorption or optically thin synchrotron (roughly $-0.5$ or steeper), the dominant-mechanism claim fails. A clumpy thermal model that reproduces the steep light curves without any synchrotron component would also falsify the conclusion.","supporting_citations":[{"cited_title":"M., Linford, J","cited_arxiv_id":null,"evidence_quote":"supplies the standard synchrotron model and parameter framework that the paper modifies with a broken power-law density."},{"cited_title":"M., Wade, C","cited_arxiv_id":null,"evidence_quote":"supplies the thermal free-free emission model used as the competing baseline."},{"cited_title":"L., et al","cited_arxiv_id":null,"evidence_quote":"is the earlier VAST-based nova study whose methods and comparison sample ground the survey analysis."},{"cited_title":"2014, Science, 345, 554","cited_arxiv_id":null,"evidence_quote":"establishes Fermi-LAT gamma-ray detections of classical novae, the baseline for the radio-gamma association."},{"cited_title":"2026, A&A, 705, A19","cited_arxiv_id":null,"evidence_quote":"supports shock activity and gamma-ray behavior for V1723 Sco."}],"review_version":1}