A Study Revealing Physical Attributes of Supernova Remnant in G321.3-3.9

We present a radio analysis of the recently identified supernova remnant G321.3-3.9 using archival multi-wavelength data spanning 88-2304 MHz. The source exhibits an elliptical shell-like morphology (1.3 deg x 1.7 deg) and a relatively flat non-thermal spectral index of alpha = -0.40 +/- 0.03. The distance is estimated using both the Sigma-D relation (1.6-2.9 kpc) and tentative associations with HI structures, the latter suggesting a near-side solution of 2.5-3.3 kpc, though the physical connection remains uncertain.

💡 Research Summary

This paper presents a comprehensive radio‑frequency analysis of the recently confirmed supernova remnant (SNR) G321.3‑3.9, employing archival multi‑wavelength data that span from 88 MHz to 2.3 GHz. The authors first assemble low‑frequency images from the GLEAM survey (four sub‑bands centered at 88, 118, 154, and 200 MHz) and higher‑frequency maps from CHIPASS (1.4 GHz) and the S‑Band Polarization All‑Sky Survey (S‑PASS, 2.3 GHz). After identifying and subtracting 19 unrelated point sources using two‑dimensional Gaussian fitting, all images are convolved to the lowest common resolution (the 14.4′ beam of CHIPASS) and re‑gridded onto a common spatial grid.

Using a consistent source region (an elliptical shell of 1.3° × 1.7°) and a surrounding background annulus, the authors extract integrated flux densities at each frequency. To assess systematic uncertainties arising from background selection, nine alternative background definitions are tested; the resulting covariance matrix is incorporated into a generalized least‑squares fit. The resulting radio continuum spectrum is well described by a single power law with a spectral index α = ‑0.40 ± 0.03 (statistical) ± 0.01 (systematic). This value is significantly flatter than the earlier estimate (α = ‑0.8 ± 0.2) derived from only two frequencies, indicating that the broader frequency coverage captures the true synchrotron spectrum more reliably.

Polarization data from S‑PASS are used to compute the polarized intensity (PI = √(Q² + U² − σ²)) and to map the projected magnetic‑field orientation (B‑vectors). The polarized fraction reaches ~10 % along the southern rim, and the B‑vectors display a predominantly tangential, circular pattern, consistent with magnetic‑field compression by the expanding shock front. However, the spatial mismatch between polarized and total intensity suggests that part of the polarized signal may arise from unrelated foreground/background HI structures.

Distance estimation proceeds via two independent methods. First, the empirical Σ–D relation is applied: the 1 GHz surface brightness derived from the fitted spectrum yields a physical diameter of ~56.9 pc, which, combined with the angular size, gives a nominal distance of ~2.2 kpc. Accounting for the intrinsic scatter of the Σ–D relation expands the plausible range to 1.6–2.9 kpc. Second, the authors search for associated interstellar gas using HI4PI and CO (Dame et al. 2001) data. No significant CO emission is detected, likely due to the coarse 0.5° beam and limited sensitivity. In the HI data, a partial arc‑like feature centered near l ≈ 321.5°, with a line‑of‑sight velocity of v_LSR ≈ ‑50 km s⁻¹, is spatially coincident with the SNR’s northern boundary and a central brightness depression. Using the Galactic rotation curve calibrated by VLBI parallaxes (Reid et al. 2014, 2019), this velocity translates to two kinematic distances: a near solution of 2.5–3.3 kpc and a far solution of 9.5–10.3 kpc (68 % credible intervals). Monte‑Carlo propagation of uncertainties in the velocity, solar motion, and rotation‑curve parameters is performed to obtain these intervals. The authors caution that the physical association between the HI feature and the SNR remains tentative.

In the discussion, the authors interpret the flat spectral index and moderate polarized fraction as indicative of a relatively young remnant that may still be in the transition from free‑expansion to the Sedov‑Taylor phase, or as a sign of interaction with a non‑uniform interstellar medium that modifies the electron energy distribution. The distance estimates from both methods converge on a near‑side distance of roughly 2.5–3 kpc, though the far‑side solution cannot be excluded without further evidence.

The paper concludes that G321.3‑3.9 is a large (∼100 pc) elliptical SNR with a flat radio spectrum, detectable polarized emission, and a likely distance of ~2.5–3 kpc. Future high‑resolution HI and CO observations, deeper X‑ray and optical spectroscopy, and refined magnetohydrodynamic modeling are recommended to confirm the gas association, determine the remnant’s age and explosion energy, and clarify its role in the surrounding interstellar medium.