{"id":"2c2e78c5-7756-4d15-bd17-ef769dde2354","arxiv_id":"2506.15067","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The LMC supernova remnant Veliki is measured to be one of the largest known, with a flat radio spectral index alpha=-0.26±0.02, interpreted as an old, fully radiative remnant with high shock compression and thermal bremsstrahlung.","lead":"New radio images of the Large Magellanic Cloud supernova remnant J0450.4-7050, nicknamed Veliki, show it is even larger than previously measured and has an unusually flat radio spectrum. The authors argue it is an old, radiative remnant whose shock compression and thermal emission explain its brightness.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The flat integrated spectral index is not confirmed by the paper's own high-frequency data, and the radiative/compression interpretation rests on this unsecured value.","rationale":"The reader identified the same load-bearing premise: the integrated spectral index. My reading sharpens it: the issue is not only that the heterogeneous points have conventional uncertainties, but that the paper's own high-frequency observations and the proposed two-component model point in opposite directions. The observational results (morphology, polarisation, environment) remain useful, and the size claim is separately caveated by the authors in Sec. 4.1. A conditional acceptance requiring the spectral-index analysis to be redone—or the interpretation to be downgraded—is the right verdict, so I do not move the reader's decision. Table 1 actually lists 19 flux measurements (17 dof), not 17 points; this correction does not change the assessment.","tokens_in":24286,"tokens_out":14044,"duration_ms":150304,"concrete_test":"Refit the Table 1 flux densities as three independent power laws in the sub-bands 88–200 MHz, 408–1295 MHz, and 2300–8850 MHz, propagating both the tabulated statistical errors and the stated 10–20% systematic errors. If the sub-band slopes are mutually inconsistent at >2σ, the single value α=-0.26±0.02 is not supported by the data, and the derived compression ratio and thermal fraction must be revised or removed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims—r≈6.8 and a 58.6% thermal bremsstrahlung contribution—are derived directly from α=-0.26±0.02 (Sec. 4.4.3, 4.4.4), but that slope is not robustly established. Sec. 3.3.1 fits 19 heterogeneous flux points (17 dof) with assigned 10–20% errors and finds reduced χ²=0.53, meaning the errors are generous and the formal ±0.02 underestimates systematics. Sec. 3.3.2's same-beam spectral index map over 888–1300 MHz gives α=-0.46±0.36; the dismissal of this as a 'smaller frequency range' artifact is not valid, since a single power law has the same slope in any sub-range—only the uncertainty changes. Table 1 itself is not a single power law: 944→1295 MHz gives α≈-0.75 and 888→1295 MHz α≈-0.33, while 2300→8850 MHz gives α≈-0.3 to -0.5. Moreover, a two-component spectrum with 58.6% thermal emission at 1 GHz would flatten toward α≈-0.1 to -0.2 at higher frequencies, opposite to the observed 888–1300 MHz and 4.75–8.85 GHz behavior. The stated formula α_total=α_th+α_nt (Sec. 4.4.3) is also incorrect for composite spectra, which combine by flux weighting. The flat slope therefore looks like an average of incompatible sub-band slopes rather than a secure physical spectral index.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents new ASKAP and MeerKAT radio-continuum observations of the LMC supernova remnant J0450.4-7050 (nicknamed Veliki), together with archival MWA, MOST, optical, infrared, X-ray, H I, and CO data. The authors report a revised physical size of about 150 by 81 pc, an integrated radio spectral index of alpha = -0.26 +/- 0.02, low fractional polarization, and a roughly uniform spectral index map. Based on these results, they argue that Veliki is an old, predominantly radiative SNR whose flat spectrum arises from a high shock compression ratio (r about 6.8) combined with a thermal bremsstrahlung contribution of about 58.6% at 1 GHz. The interpretation is placed in the context of the LMC's low-metallicity, density-structured environment, and alternative scenarios such as a pulsar wind nebula, second-order Fermi acceleration, and molecular cloud interaction are considered and rejected.","tokens_in":24651,"tokens_out":7842,"duration_ms":73533,"significance":"The new high-resolution imaging and multi-wavelength comparison are a useful observational contribution to LMC SNR studies, and the measured size, morphology, polarimetry, and environmental analysis are of interest. If the flat integrated spectral index were robust, the paper would present a significant challenge to simple DSA expectations for evolved remnants. The authors are candid about several caveats, including the uncertain nature of the southern filaments and the heterogeneous flux-density measurements. However, the central quantitative conclusions rest on a spectral index whose robustness is not established and on a two-component spectral model that is mis-specified as written. As it stands, the radiative-shock/thermal-bremsstrahlung interpretation is not supported by the evidence presented.","major_comments":[{"comment":"The quoted integrated spectral index alpha = -0.26 +/- 0.02 is not robustly established. The fit combines 17 flux-density points with adopted 10-20% uncertainties and yields reduced chi^2 = 0.53, which means the errors are sufficiently generous that the formal +/- 0.02 underestimates systematic uncertainties. More directly, the data in Table 1 are not described by a single power law: the 944 to 1295 MHz pair gives alpha about -0.75, the 888 to 1295 MHz pair gives alpha about -0.33, and the 2300 to 8850 MHz range gives alpha about -0.3 to -0.5. The paper needs to quantify this scatter, for example with per-sub-band fits or an explicit curvature test, before using alpha = -0.26 as the foundation of the physical model.","section":"Sec. 3.3.1 and Table 1"},{"comment":"The dismissal of the spectral index map's average value alpha = -0.46 +/- 0.36 as an artifact of the 'smaller frequency range' is invalid: a single power law has the same slope in any frequency sub-range, with only the uncertainty changing. The discrepancy between the map and the integrated value therefore needs a quantitative explanation, such as calibration offsets between the ASKAP and MeerKAT images, missing extended flux at high resolution, or genuine spectral curvature.","section":"Sec. 3.3.2"},{"comment":"The two-component model is mis-specified as written. The statement 'alpha_total = alpha_thermal + alpha_non-thermal' is not the correct combination law for a composite spectrum; the effective spectral index is a flux-weighted average, alpha_eff = (S_th alpha_th + S_nt alpha_nt) / (S_th + S_nt). Because the 58.6% thermal fraction is derived from this model, it must be recomputed with the proper formula, or the intended formula must be stated unambiguously, and the uncertainty in the derived thermal fraction should be reported.","section":"Sec. 4.4.3"},{"comment":"The thermal-bremsstrahlung interpretation is internally inconsistent with the frequency dependence of the data. With a 58.6% thermal contribution at 1 GHz and alpha_th = -0.1, the spectrum should flatten at higher frequencies, approaching alpha about -0.1 as the thermal component dominates; instead, the 888-1300 MHz and 4.75-8.85 GHz sub-band slopes are steeper, approximately -0.3 to -0.75. This is not merely an absence of observable curvature; it is the opposite trend from the model's prediction.","section":"Sec. 4.4.3 and Sec. 4.4.4"},{"comment":"The revised size of 150 by 81 pc is presented as a headline result even though the authors state that the southern filaments may be leaked ionising radiation rather than physical shell material. Since the 'one of the largest SNRs' claim depends on including these filaments, the paper should either adopt a conservative shell size excluding ambiguous structures or provide a quantitative criterion, such as radio-optical morphological correspondence, for including them.","section":"Sec. 4.1 and Sec. 3.1"}],"minor_comments":[{"comment":"The frequency range '88-8850 GHz' should be '88-8850 MHz'.","section":"Sec. 4.4.3"},{"comment":"The molecular cloud is called PGCC G272.62-35.35 in Sec. 3.1 and PGCC G282.62-35.35 in Sec. 4.2; the designation should be made consistent.","section":"Sec. 3.1 and Sec. 4.2"},{"comment":"The surface brightness unit should be W m^-2 Hz^-1 sr^-1, not 'W m^-1 Hz^-2 sr^-1' as currently printed.","section":"Fig. 9 caption"},{"comment":"The phrase 'one of the lowest average radio spectral indices' is ambiguous; since alpha is negative, 'flattest' or 'least negative' would be clearer.","section":"Abstract and Sec. 5"},{"comment":"The sign-convention footnotes are helpful, but the use of alpha = 0.26 in the compression-ratio equation may confuse readers; adopting one sign convention throughout would improve readability.","section":"Secs. 3.5 and 4.4.4"}],"recommendation":"major_revision","confidential_remarks":"The observational dataset is valuable and the size/morphology/polarimetry portions are likely publishable, but the spectral-index analysis and the physical interpretation built on it need substantial reworking. I recommend major revision rather than rejection because the data products are useful and the main technical issues appear fixable within the scope of a revised manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a useful new radio imaging study of a known LMC SNR, with real new information in the images, but the paper's central interpretation—the unusually flat integrated spectral index and everything built on it—does not survive close reading of their own data.\n\nWhat is actually new: ASKAP and MeerKAT resolve faint north-south filaments that extend the remnant to 150×81 pc, making it one of the largest known LMC SNRs. The multi-frequency flux table (17 measurements from 88 to 8850 MHz) is a useful community resource, and the optical/Hα/[O iii] context is nicely integrated with the radio. Credit where due: the authors state in Sec. 4.1 that the southern filaments may be leaked ionizing radiation rather than physical shell. That is the right kind of honesty.\n\nWhere it gets soft: the α = -0.26 ± 0.02 is built from heterogeneous flux points with assumed 10–20% errors. The reduced χ² of 0.53 tells you the error bars are generous, so the ±0.02 is a formal fit error, not a real uncertainty. More telling, their own same-beam spectral index map over 888–1300 MHz gives α = -0.46 ± 0.36, and the dismissal of this as a 'smaller frequency range artifact' does not hold: a single power law has the same slope in any sub-range, only the uncertainty changes. And the Table 1 points are not a single power law: 944→1295 MHz gives roughly -0.75, while 2300→8850 MHz gives about -0.3 to -0.5. The flat slope looks like an average of incompatible sub-band slopes.\n\nWorse, the thermal-decomposition argument in Sec. 4.4.3 writes α_total = α_th + α_nth. That is not how composite spectra work; flux densities add, so the effective index is flux-weighted. They then assume α_th = -0.1 and α_nth = -0.5, which forces a 58.6% thermal fraction at 1 GHz. This is circular and uses an incorrect formula. The compression ratio r ≈ 6.8 derived from α = -0.26 inherits the same problem.\n\nThis is not a desk reject: the imaging is solid, the object is genuinely interesting, and the size measurement (with the caveat noted) is a real update. But the flat-spectrum claim and the radiative-SNR interpretation need to be reworked before this can be trusted. Either quantify the systematics honestly, fit a two-component model by adding flux densities, and treat the sub-band slopes as possible evidence of curvature, or drop the strong physical claims.\n\nRecommendation: send it to peer review, but with explicit instruction that the spectral index analysis requires major revision.","headline":"New MeerKAT/ASKAP imaging reveals a larger, still interesting LMC SNR, but the paper's flat-spectrum and radiative-phase interpretation rests on a spectral index that their own data contradict.","tokens_in":25296,"tokens_out":3469,"would_cite":true,"duration_ms":31944,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.85.Bh","98.38.Mz"],"model":"deepseek-v4-flash","headline":"New radio maps enlarge supernova remnant Veliki to 150×81 pc and reveal an unusually flat radio spectrum pointing to a fully radiative shock.","keywords":["ISM: supernova remnants","supernovae: general","supernovae: individual: J0450-7050","radio continuum: general","Large Magellanic Cloud","spectral index","diffuse shock acceleration","thermal bremsstrahlung"],"falsifier":"Re-fit the integrated spectrum after treating each survey's absolute flux calibration as a correlated systematic error, and also produce a matched-resolution spectral-index map that includes the 88–200 MHz low-frequency points; if the slope moves from $\\alpha = -0.26$ to $\\alpha \\lesssim -0.35$ under either test, or the low-frequency map is uniformly steeper than the integrated value, the flat-spectrum case and the derived $r \\approx 6.8$ and 58.6% thermal fraction would be refuted.","tokens_in":24130,"feed_emoji":"📡","tokens_out":12794,"duration_ms":112085,"temperature":0.7,"pith_summary":"This paper reports new high-resolution radio-continuum observations of the Large Magellanic Cloud supernova remnant J0450.4–7050, nicknamed Veliki, and re-examines its multi-frequency properties. It claims that Veliki is substantially larger than previously measured — 150×81 pc — making it one of the largest known supernova remnants, and that its integrated radio spectrum is unusually flat, $\\alpha = -0.26 \\pm 0.02$, with little spatial variation. The paper argues that the most likely explanation is an old, predominantly radiative remnant whose shocks have a high compression ratio $r \\approx 6.8$, flattening the non-thermal spectrum, combined with a thermal bremsstrahlung contribution of about 58.6% at 1 GHz. If this interpretation is right, Veliki becomes a rare example of a giant, radiatively cooled supernova remnant and a direct test of diffuse shock acceleration outside the standard adiabatic shock limit.","feed_headline":"New radio maps reveal Veliki as one of the largest supernova remnants","feed_subtitle":"Its flat radio spectrum hints at a fully radiative shock with compression far above the standard limit.","key_machinery":"The load-bearing machinery is the radio spectral index itself plus two theoretical links. First, the integrated value $\\alpha = -0.26 \\pm 0.02$, derived from a 17-point power-law fit to flux densities spanning 88–8850 MHz, fixes the global emission law $S \\propto \\nu^{\\alpha}$. Second, the diffuse-shock-acceleration (DSA) compression-ratio formula $\\alpha = 3/(2(r-1))$ converts that index into $r \\approx 6.8$, and a two-component spectral model (non-thermal synchrotron with $\\alpha = -0.5$ plus optically thin thermal bremsstrahlung with $\\alpha = -0.1$) assigns a 58.6% thermal fraction at 1 GHz. Radiative-shock theory, where the compression ratio can approach the square of the isothermal Mach number, supplies the physical justification for $r > 4$. These pieces, not any single image, carry the interpretive claim that Veliki is a fully radiative supernova remnant.","core_discovery":"On the paper's own terms, the discovery is that Veliki is one of the largest known supernova remnants and its radio spectrum is flatter than standard shock acceleration predicts. Using 17 flux-density measurements spanning 88–8850 MHz, the authors obtain an integrated spectral index $\\alpha = -0.26 \\pm 0.02$, which corresponds, through the diffuse-shock-acceleration relation $\\alpha = 3/(2(r-1))$, to a shock compression ratio $r \\approx 6.8$ — well above the strong-shock limit $r = 4$ for an ideal adiabatic gas. They further model the spectrum as non-thermal synchrotron plus optically thin thermal bremsstrahlung and find that about 58.6% of the 1 GHz flux would need to be thermal to reproduce the flat index. Combined with a bright [SII]/H$\\alpha$ shell and a soft X-ray interior, the paper concludes that Veliki is most likely a fully radiative supernova remnant, with the flat spectrum produced by the high compression ratio and thermal contamination, and its large size and high surface brightness explained by delayed cooling in the low-metallicity LMC environment and possibly a higher-than-normal explosion energy (about $8.6\\times10^{51}$ erg).","pith_inferences":["Beyond the paper, the flat-index claim could be checked independently by re-fitting the same 17 flux points with per-survey flux-scale nuisance parameters; if a correlated-error fit yields a steeper index, the radiative-shock and thermal-fraction story would lose its foundation.","Beyond the paper, the paper's dismissal of the steeper spectral-index map ($\\alpha = -0.46 \\pm 0.36$) as a frequency-range artifact could be tested by producing a matched-resolution map that includes the 88–200 MHz low-frequency points; a genuinely uniform flat map would strengthen the thermal-contamination story, while a steep map would point to calibration or resolution issues.","Beyond the paper, the southern radio filaments without optical counterparts may be ionised leakage rather than true shell material; targeted optical and infrared follow-up could decide whether Veliki's 150 pc extent is real or an upper limit, which would also change its position as a surface-brightness outlier.","Beyond the paper, if the molecular cloud near the north-western rim is truly interacting, future very-high-energy gamma-ray observations should detect hadronic emission; Veliki is a natural target for such observations once the LMC is covered at sufficient sensitivity."],"forward_implications":["Veliki's size and surface brightness place it outside the normal $\\Sigma$–D evolutionary tracks, so if the interpretation is right it becomes a rare giant remnant formed by a delayed radiative transition in the low-metallicity LMC environment.","A compression ratio $r \\approx 6.8$ means the shock has crossed the adiabatic $r = 4$ limit, making Veliki one of the clearest cases of diffuse shock acceleration operating in a fully radiative shock.","The 58.6% thermal fraction at 1 GHz should become visible as spectral flattening at higher frequencies, with any low-frequency turnover lying below 88 MHz; microwave and infrared measurements could test this directly.","If the fitted age of about 43,000 years and explosion energy of about $8.6\\times10^{51}$ erg are correct, Veliki is a high-energy, evolved remnant rather than a typical middle-aged supernova remnant."],"supporting_citations":[{"why":"Supplies the previous radio spectral index ($\\alpha = -0.43 \\pm 0.06$), the 102×75 pc size estimate, and several catalogued flux densities against which the new results are compared.","marker":"Čajko et al. 2009"},{"why":"Provides the original 4800 MHz Parkes detection and flux density that first identified the source as MC 11.","marker":"McGee et al. 1972"},{"why":"Supplies the 408 MHz flux density, the largest outlier in the spectral-index fit that is used to test for spectral curvature.","marker":"Clarke et al. 1976"},{"why":"Contributes the 2300 MHz catalogued flux density that enters the integrated spectral-index fit.","marker":"Filipović et al. 1996"},{"why":"Contributes the 2700 MHz catalogued flux density that enters the integrated spectral-index fit.","marker":"Filipović et al. 1998"},{"why":"Provides the standard diffuse-shock-acceleration prediction $\\alpha = -0.5$, the baseline from which Veliki's flat index deviates.","marker":"Bell 1978"},{"why":"Supplies the compression-ratio relation $\\alpha = 3/(2(r-1))$ and the radiative-shock theory used to derive $r \\approx 6.8$.","marker":"Urošević 2014"},{"why":"Supplies the two-component thermal/non-thermal spectral model used to estimate a 58.6% thermal contribution at 1 GHz.","marker":"Onić et al. 2012"},{"why":"Supplies the X-ray/optical age, expansion velocities, and ambient density that anchor the radiative-phase interpretation.","marker":"Williams et al. 2004"},{"why":"Supplies the forward-shock evolution model that yields the ~43,000 yr age and ~$8.6\\times10^{51}$ erg explosion energy.","marker":"Leahy and Williams 2017"}],"fun_headline_variants":["Veliki: one of the largest SNRs, with a shock compression of 6.8","Radio view reveals Veliki as a giant SNR with an unexpectedly flat spectrum","Veliki's flat radio spectrum hints at shock compression beyond the standard limit","Meet Veliki: a supernova remnant that exceeds the shock compression limit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the single number $\\alpha = -0.26 \\pm 0.02$ accurately describes the whole remnant's radio emission, even though it is fitted to heterogeneous 88–8850 MHz flux densities with largely conventional 10–20% uncertainties and one strong outlier; if that number is off, the compression ratio and thermal fraction built on it collapse.","fun_headline_variants_meta":{"raw":{"variants":["Veliki: one of the largest SNRs, with a shock compression of 6.8","Radio view reveals Veliki as a giant SNR with an unexpectedly flat spectrum","Veliki's flat radio spectrum hints at shock compression beyond the standard limit","Meet Veliki: a supernova remnant that exceeds the shock compression limit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000834,"raw_usage":{"total_tokens":3676,"prompt_tokens":1021,"completion_tokens":2655,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":637,"completion_tokens_details":{"reasoning_tokens":2571}},"tokens_in":637,"tokens_out":2655,"duration_ms":23544,"temperature":1.0,"reasoning_tokens":2571,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:45:50.195565+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit the integrated spectrum after treating each survey's absolute flux calibration as a correlated systematic error, and also produce a matched-resolution spectral-index map that includes the 88–200 MHz low-frequency points; if the slope moves from $\\alpha = -0.26$ to $\\alpha \\lesssim -0.35$ under either test, or the low-frequency map is uniformly steeper than the integrated value, the flat-spectrum case and the derived $r \\approx 6.8$ and 58.6% thermal fraction would be refuted.","supporting_citations":[{"cited_title":"X., Brooks, J","cited_arxiv_id":null,"evidence_quote":"Provides the original 4800 MHz Parkes detection and flux density that first identified the source as MC 11."},{"cited_title":"M., Chu, Y","cited_arxiv_id":null,"evidence_quote":"Supplies the X-ray/optical age, expansion velocities, and ambient density that anchor the radiative-phase interpretation."}],"review_version":2}