{"id":"c7e9231d-3da5-4a2c-b5ae-a8f60acb5039","arxiv_id":"2506.23588","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"WR16's circumstellar shell is detected in ASKAP radio data, giving a shell age near 9500 years and an expansion speed near 280 km/s.","lead":"Astronomers report the first radio detection at 943.5 MHz of the shell of gas around the Wolf-Rayet star WR16, measuring its flux, size, age, and expansion speed. The results support the idea that WR16 went through a violent eruptive phase before becoming a Wolf-Rayet star.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Shell age and expansion velocity rest on an unquantified center fit and a constant-velocity assumption; a biased fit of the bright NW rim could shift these beyond the quoted uncertainties.","rationale":"The paper's central claim is the first radio-continuum detection of WR16's shell, with the derived age, expansion velocity, and LBV-phase interpretation following from a geometric model. The detection itself is credible: the EMU image shows a ring matching Hα and WISE, and the paper reports a 10% flux uncertainty, though the extraction method is underdescribed. However, the quantitative conclusions demanded by the abstract rest entirely on the §3.3 kinematic construction. The age is the ratio of the star's projected offset from the shell center to its tangential velocity; the expansion velocity is the shell radius divided by that age. If either the center or the velocity is wrong, both age and expansion velocity change. The paper does not quantify the center-fitting uncertainty, and the bright NW rim—attributed to the star's motion—is precisely the region most likely to bias a circular fit toward the star, shrinking the offset and age. A shift of 10 arcsec in the center changes the age by 22%, beyond the quoted 14% error. The constant-velocity assumption is also unverified for a runaway WN8 star. The reader identified the same weak point, and I agree. The detection and thermal-origin inference survive this concern, so the verdict remains conditional, pending a robustness check of the fitted center.","tokens_in":10566,"tokens_out":14572,"duration_ms":153011,"concrete_test":"Re-fit a circle to the 943.5-MHz EMU image with the bright NW quadrant (the ~90° sector centred on the star's position) masked, and repeat on the WISE 22-µm image as an independent tracer. If the best-fit center shifts by more than ~10 arcsec (or if the WISE center differs from the EMU center by that amount), the derived age and expansion velocity are not robust and must be recomputed with a proper center uncertainty; this would undermine the LBV-phase conclusion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In §3.3, the shell's age (~9500 yr) and average expansion velocity (~280 km/s) are derived from the 44.81-arcsec separation between WR16 and the 'geometric center' of the 8.42-arcmin fitted circle, divided by the current tangential peculiar velocity. This derivation requires (i) that the fitted circle's center is the true expansion origin and (ii) that WR16's transverse velocity has been constant over the inferred ~10^4-yr lifetime. Neither condition is demonstrated. The paper gives no uncertainty for the fitted center, and the circular fit is not described. The bright NW rim, which the paper attributes to the star's motion (§3.5), could bias a naive fit toward the star, reducing the measured offset and thus the age. A center displacement of ~10 arcsec (~22% of the offset) changes the age by ~22%, exceeding the quoted 14% (1300-yr) uncertainty. The constant-velocity assumption is also untested: WR16 is a WN8 runaway, and a binary kick or cluster ejection could have altered its PM direction over 9500 yr. Since the LBV-versus-RSG conclusion depends directly on the expansion velocity, this kinematic model is the weakest load-bearing pillar of the derived claims.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a 943.5 MHz ASKAP EMU detection of an extended shell around the Wolf-Rayet star WR16, with a measured shell flux density of 72.2±7.2 mJy and a stellar flux density of 0.35±0.04 mJy. Combining the EMU stellar flux with archival ATCA measurements yields a spectral index α = +0.74±0.02, interpreted as thermal free-free emission. The authors adopt a Gaia-based distance of 2.28±0.09 kpc, derive a shell diameter of 5.57±0.22 pc, use Gaia proper motions to estimate the shell age (~9500±1300 yr) and average expansion velocity (280±40 km/s), and conclude that the shell likely originated in an LBV phase. They also derive a mass-loss rate of 1.753×10^-5 M_sun/yr and a lower limit on ionizing photons N_UV > 1.406×10^47 s^-1, and report a candidate large-scale CO wind-blown bubble around the system.","tokens_in":10829,"tokens_out":4687,"duration_ms":51474,"significance":"If the detection and the kinematic interpretation hold, this is the first radio-continuum detection of the shell around WR16, a useful addition to the small sample of resolved circumstellar shells around Wolf-Rayet stars. The paper combines multi-wavelength archival data (ATCA, WISE, SuperCOSMOS, NANTEN CO, Gaia) with new ASKAP EMU imaging, and the astrometric use of Gaia proper motions to estimate a shell age is a valuable approach. The derived LBV-phase interpretation is physically interesting. However, the central claims rest on two pillars that need strengthening: the shell flux measurement is not described and its uncertainty is adopted rather than measured, and the age/expansion velocity rest on an unquantified geometric-center assumption and a constant-velocity assumption. The thermal nature of the shell is explicitly admitted to be inferred, not spectroscopically confirmed. These limitations do not invalidate the detection, but they require revision before the derived quantities can be considered robust.","major_comments":[{"comment":"The shell flux density measurement is not described. The paper does not state the aperture or integration region used to measure the 72.2±7.2 mJy shell flux, how background emission was subtracted, or how the bright north-west rim and the central point source were treated. The 10% uncertainty is adopted from Smeaton et al. (2024b) and Filipović et al. (2024) rather than derived from the present image noise and systematic effects. Since the shell detection is the central claim, please specify the measurement method and provide a detailed error budget, including the effect of the ~37 µJy/beam rms and possible contamination from the brighter NW rim.","section":"§3.1, Table 1"},{"comment":"The age (~9500±1300 yr) and average expansion velocity (280±40 km/s) are derived assuming that the shell originated at the geometric center of the fitted 8.42-arcmin circle and that WR16's current transverse peculiar velocity (50.7±6.9 km/s) has been constant over ~10^4 yr. The circular fit is not described and no uncertainty is given for the fitted center. A center displacement of ~10 arcsec, which is plausible given the bright NW rim that the paper itself attributes to the star's motion, changes the age by ~22%, exceeding the quoted 14% uncertainty. The constant-velocity assumption is also untested; a binary kick or cluster ejection could have changed the PM direction over 9500 yr. Please quantify the fit uncertainty and assess the sensitivity of the age and expansion velocity to these assumptions, since the LBV-versus-RSG conclusion depends directly on them.","section":"§3.3"},{"comment":"The thermal nature of the shell is inferred only from its morphology in Hα and WISE 22 µm images; the text admits that the shell has not been observed at other radio frequencies and that thermal origin cannot be confirmed from a spectral index. This inference is then used to support the physical interpretation. Please state explicitly that the shell's thermal nature is a working hypothesis, and discuss whether any other emission mechanism (e.g., non-thermal or bremsstrahlung from a shock) could produce the observed morphology. In addition, the mass-loss rate and N_UV lower limit assume that the stellar flux at 943.5 MHz is entirely free-free wind emission; please discuss variability and possible contamination from the shell or from unrelated background sources.","section":"§3.1, §3.4"},{"comment":"The statement that recalculating Leitherer et al. (1995) mass-loss rate with a distance of 2.28 kpc gives 2.188×10^-5 M_sun/yr is not consistent with the Wright-Barlow scaling used in the same section: if the same observed flux is used, the mass-loss rate scales as d^{3/2}, which would give ≈3.8×10^-5 M_sun/yr for the original value, not 2.19×10^-5. If instead the calculation uses the new 943.5 MHz flux of 0.35 mJy, that should be stated explicitly. Please provide the full calculation with the numerical values substituted so that the comparison is reproducible.","section":"§3.4"}],"minor_comments":[{"comment":"There are several typos: 'supported by the its appearance' in the abstract, 'W olf-Rayet' and 'A TCA' spacing in the text, and 'WR16 are released' in Section 3.4 should be 'WR16 releases' or 'WR16 emits'.","section":"Abstract, §1"},{"comment":"The reference to 'Teleios (Filipović et al. submitted)' is incomplete; either provide a full citation or remove the name.","section":"§3.5"},{"comment":"The paper states that N_UV > 1.406×10^47 s^-1 is 'in concordance with the canonical value for WN8 stars of 1.259×10^49 s^-1', but the lower limit is two orders of magnitude below the canonical value and therefore does not actually constrain it. Please rephrase to avoid overinterpretation.","section":"§3.4"},{"comment":"The column headers 'S944', 'S2.40', 'S4.80', and 'S8.64' are not defined in the table caption; please add explicit frequency labels or a note explaining the notation.","section":"Table 1"},{"comment":"The CO bubble association is presented as plausible but the text says only that the objects 'appear to be physically associated'; please mark this clearly as a tentative suggestion and consider adding a quantitative morphological or kinematical test.","section":"§3.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of PASA and the detection is potentially interesting. The main concern is that the derived age, expansion velocity, and the LBV conclusion rely on untested kinematic assumptions and an unquantified geometric fit. These are fixable with additional analysis and should be addressed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know: this is a legit first radio detection of the WR16 shell at 943.5 MHz, with a clean stellar spectral index and a new CO bubble. The derived age (~9500 yr), expansion velocity (~280 km/s), and the LBV-phase conclusion are real numbers, but they rest on a kinematic assumption that is only partly flagged.\n\nThe 943.5 MHz shell flux (72.2±7.2 mJy) is new, and the paper's star measurement combines sensibly with the archival ATCA points to give α=+0.74±0.02, consistent with a thermal wind. The CO bubble in Section 3.2 is a nice independent find. The mass-loss rate and ionizing photon limit follow from standard Wright-Barlow formulas, and the paper honestly notes the shell's thermal nature is inferred from IR/Hα rather than from a spectral index. The 2D projection caveat on age/expansion velocity is also stated.\n\nSoft spots, in proportion. The shell flux extraction is not described: no aperture, no method, just a 10% error adopted from other papers. For the headline measurement, that's thin. The bigger issue is the shell's age and expansion velocity. Section 3.3 assumes the shell originated at the geometric center of the 8.42-arcmin circle and that WR16's peculiar velocity has been constant for ~10^4 yr. The circular fit itself is not described and no uncertainty is given for the center. A bright NW rim could bias the fit, and a center offset of ~10 arcsec (~22% of the 44.8 arcsec star-center separation) would move the age by ~20%, more than the quoted 14% uncertainty. The constant-velocity assumption is plausible for a runaway but untested; a kick or cluster ejection could have changed the PM direction. These caveats don't kill the detection, but they mean the LBV-vs-RSG conclusion is model-dependent, not measured.\n\nMinor: the N_UV lower limit is two orders of magnitude below the canonical WN8 value; as a lower bound that's fine, but \"in concordance\" overstates the agreement.\n\nWho it's for: radio observers working on WR nebulae, and anyone compiling EMU detections of evolved massive stars. It deserves a serious referee; a good referee will ask for the flux extraction and circular fit details, and for a sensitivity test on the center offset. With those added, the kinematic numbers can stand with proper caveats.","headline":"First radio detection of WR16's shell is plausible and worth publishing, but the kinematic age, expansion velocity, and LBV conclusion rest on an unquantified center fit that a referee should probe.","tokens_in":11480,"tokens_out":2863,"would_cite":true,"duration_ms":27670,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"ASKAP EMU data give the first radio-continuum detection of the shell surrounding Wolf-Rayet star WR16 at 943.5 MHz, with flux density 72.2±7.2 mJy, and place its outburst at roughly 9,500 years ago.","keywords":["Wolf-Rayet star","WR16","circumstellar shell","radio continuum","thermal emission","proper motion","LBV phase","mass-loss rate"],"falsifier":"Measure the shell at a second radio frequency near 2 GHz with comparable angular resolution and fit its own spectral index: a value clearly below the flat free-free range would rule out the thermal-shell interpretation. Independently, re-fit the shell's centre with the bright north-west rim masked and compare with the Gaia proper-motion vector; a significantly shifted centre would alter the ~9,500-year age and 280 km s$^{-1}$ expansion velocity.","tokens_in":10385,"feed_emoji":"📡","tokens_out":14235,"duration_ms":129214,"temperature":0.7,"pith_summary":"The paper reports the first radio-continuum detection of the ring-like circumstellar shell around the Wolf-Rayet star WR16 (a massive, late-stage star), made at 943.5 MHz in the ASKAP EMU survey. The shell has a flux density of 72.2±7.2 mJy, and the star's spectral index of $\\alpha = +0.74\\pm 0.02$ between 943.5 MHz and 8.64 GHz marks its emission as thermal; the authors argue the shell is also thermal because it appears in near-infrared and H$\\alpha$ images. Using the Gaia parallax, they place WR16 at 2.28±0.09 kpc, making the 8.42-arcmin shell 5.57±0.22 pc across, and the Gaia proper motions give a north-west peculiar velocity of 50.7±6.9 km s$^{-1}$. Because that velocity vector passes through the shell's geometric centre, the star must have crossed that centre about 9500±1300 years ago, and the shell's radius then implies an average expansion velocity of 280±40 km s$^{-1}$. This expansion speed, between those of typical LBV shells and of Eta Carinae-like ejecta, is used to conclude that the star's preceding transitional phase was a Luminous Blue Variable epoch rather than a red supergiant one.","feed_headline":"ASKAP EMU detects a 9,500-year-old shell around WR16","feed_subtitle":"The 5.6-pc shell dates a past LBV outburst, not a red supergiant phase.","key_machinery":"The load-bearing mechanism is the measured geometry: an extended, nearly circular radio shell fitted with an angular diameter of 8.42 arcmin at 943.5 MHz, overlaid with the Gaia proper-motion vector of the star. The star's peculiar velocity (50.7±6.9 km s$^{-1}$) points north-west and passes through the geometric centre of the circle; the offset between star and centre (44.81 arcsec, or 0.49 pc) divided by that velocity gives the crossing time of about 9500 years, and the shell radius (2.79 pc) divided by that time gives the average expansion speed of 280±40 km s$^{-1}$. The thermal interpretation rests on the spectral index $\\alpha = +0.74 \\pm 0.02$ computed from the new EMU point plus earlier ATCA points, with the shell's thermal nature inferred from its near-infrared and H$\\alpha$ morphology; the Wright-Barlow free-free relations then convert the observed flux into the mass-loss rate and ionising-photon lower limit.","core_discovery":"We present the first radio-continuum detection of the circumstellar shell around the WN8-type Wolf-Rayet star WR16 at 943.5 MHz using ASKAP EMU data. The shell flux is 72.2±7.2 mJy; combined with archival ATCA measurements at 2.4, 4.8, and 8.64 GHz, the star's spectral index is $\\alpha = +0.74 \\pm 0.02$, indicating thermal free-free emission, and the shell's appearance in near-infrared and H$\\alpha$ suggests it shares that thermal origin. From the Gaia DR3 parallax and proper motions we derive a distance of 2.28±0.09 kpc, a shell diameter of 5.57±0.22 pc, and a peculiar tangential velocity of 50.7±6.9 km s$^{-1}$ directed north-west, toward the brighter part of the shell. Because the proper-motion vector passes through the geometric centre of the fitted 8.42-arcmin circle, the shell's origin is placed at that centre, giving an age of roughly 9500±1300 yr and an average expansion velocity of 280±40 km s$^{-1}$, which supports an LBV rather than RSG transitional phase. The 943.5 MHz flux also yields a mass-loss rate of $1.753\\times 10^{-5}$ $M_\\odot$ yr$^{-1}$ and a lower limit on ionising photons of $N_{\\rm UV} > 1.406\\times 10^{47}$ s$^{-1}$.","pith_inferences":["An optically thin thermal shell should have a spectral index near zero, so observing the shell alone at a second radio frequency would directly test the thermal-shell claim, which currently rests on the star's index and the shell's IR/H$\\alpha$ appearance.","The geometric-centre-plus-proper-motion method could be applied to a sample of runaway Wolf-Rayet stars with circular nebulae; a clustering of derived expansion velocities near LBV values would give independent support for the LBV-before-WR interpretation.","If WR16 is indeed approaching its bright north-west rim, spatially resolved spectroscopy of that rim should show shocked or broadened H$\\alpha$ emission; detecting such emission would confirm the interaction and give a direct, geometry-independent expansion measure."],"forward_implications":["If the age and expansion velocity hold, WR16's last major eruption occurred about 9,500 years ago and has been expanding at roughly 280 km s$^{-1}$, placing it between canonical LBV shell speeds near 50 km s$^{-1}$ and Eta Carinae-type ejecta near 600 km s$^{-1}$.","WR16 is moving toward the bright north-west rim of its shell, so the authors expect that rim to continue brightening as the star compresses it.","The derived mass-loss rate of $1.753\\times10^{-5}$ $M_\\odot$ yr$^{-1}$ is close to the canonical WN8 value, and the implied ionising-photon rate exceeds $1.4\\times10^{47}$ s$^{-1}$, consistent with a thermal radio source.","The ~30 pc CO bubble expanding at ~7 km s$^{-1}$ with a ~4 Myr dynamical timescale is identified as a wind-blown bubble from WR16's main-sequence progenitor, linking the small inner shell to the star's earlier evolutionary stage.","EMU's sensitivity to diffuse, low-surface-brightness radio emission should reveal similar extended thermal shells around other Wolf-Rayet stars that have not been seen at centimetre wavelengths."],"supporting_citations":[{"why":"Supplies the EMU survey whose 943.5 MHz image is the new detection of the shell.","marker":"Norris et al. (2021)"},{"why":"Provides earlier ATCA flux densities at 4.8 and 8.64 GHz used in the spectral index fit and in the mass-loss comparison.","marker":"Leitherer et al. (1995)"},{"why":"Provides the 2.4 GHz ATCA flux density used in the spectral index fit.","marker":"Chapman et al. (1999)"},{"why":"Supplies the parallax and proper motions that set the distance, peculiar velocity, and shell age.","marker":"Gaia Collaboration (2020)"},{"why":"Gives the free-free emission equations and canonical spectral behaviour used to compute the mass-loss rate and ionising-photon limit.","marker":"Wright & Barlow (1975)"},{"why":"Previous study of WR16's nebulosity; supplies the angular-size comparison and the peculiar-velocity method used here.","marker":"Cichowolski et al. (2020)"},{"why":"Provides the WN8 classification and canonical mass-loss and ionising-photon values used as benchmarks.","marker":"Crowther (2007)"},{"why":"Supplies WR16's coordinates and wind terminal velocity of 630 km s$^{-1}$, used to argue the shell started fast and slowed.","marker":"Toalá et al. (2015)"}],"fun_headline_variants":["First radio detection of WR16's shell reveals its 9500-year age","WR16's shell: a 9500-year-old LBV outburst traced by radio","ASKAP spots WR16's shell, hinting at LBV not red supergiant","Radio shell around WR16 points to a past LBV eruption"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derived age and expansion velocity assume the shell formed at the geometric centre of the fitted 8.42-arcmin circle and that WR16's current proper-motion velocity has been constant for about 9,500 years; if the bright north-west rim biases the fitted centre or the star's direction or speed changed, both numbers change.","fun_headline_variants_meta":{"raw":{"variants":["First radio detection of WR16's shell reveals its 9500-year age","WR16's shell: a 9500-year-old LBV outburst traced by radio","ASKAP spots WR16's shell, hinting at LBV not red supergiant","Radio shell around WR16 points to a past LBV eruption"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000244,"raw_usage":{"total_tokens":1751,"prompt_tokens":1382,"completion_tokens":369,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":998,"completion_tokens_details":{"reasoning_tokens":284}},"tokens_in":998,"tokens_out":369,"duration_ms":4358,"temperature":1.0,"reasoning_tokens":284,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:37:40.206884+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the shell at a second radio frequency near 2 GHz with comparable angular resolution and fit its own spectral index: a value clearly below the flat free-free range would rule out the thermal-shell interpretation. Independently, re-fit the shell's centre with the bright north-west rim masked and compare with the Gaia proper-motion vector; a significantly shifted centre would alter the ~9,500-year age and 280 km s$^{-1}$ expansion velocity.","supporting_citations":[{"cited_title":"P., Marvil, J., Collier, J","cited_arxiv_id":null,"evidence_quote":"Supplies the EMU survey whose 943.5 MHz image is the new detection of the shell."},{"cited_title":"M., & Koribalski, B","cited_arxiv_id":null,"evidence_quote":"Provides earlier ATCA flux densities at 4.8 and 8.64 GHz used in the spectral index fit and in the mass-loss comparison."},{"cited_title":"U., Suad, L","cited_arxiv_id":null,"evidence_quote":"Previous study of WR16's nebulosity; supplies the angular-size comparison and the peculiar-velocity method used here."}],"review_version":1}