Evolutionary Map of the Universe: Detection and Analysis of the Shell Surrounding the Runaway Wolf-Rayet Star WR16

TL;DR

This study reports the first radio detection of the circumstellar shell around Wolf-Rayet star WR16, analyzing its properties, origin, and implications for stellar evolution using multi-frequency observations and Gaia data.

Contribution

It provides the first radio continuum detection of WR16's shell, estimates its physical parameters, and discusses its evolutionary stage based on multi-wavelength data and Gaia measurements.

Findings

Detected the shell at 943.5 MHz with 72.2 mJy flux.

Estimated shell size to be 5.57 pc at 2.28 kpc distance.

Determined the shell's age to be approximately 9500 years.

Abstract

We present the first radio--continuum detection of the circumstellar shell around the well-known WN8 type Wolf-Rayet star WR16 at 943.5\,MHz using the \ac{ASKAP} \ac{EMU} survey. At this frequency, the shell has a measured flux density of 72.2$\pm$7.2\,mJy. Using previous \ac{ATCA} measurements at 2.4, 4.8, and 8.64~GHz, as well as the \ac{EMU} observations of the star itself, we determine a spectral index of $\alpha\,=\,+0.74\pm0.02$, indicating thermal emission. We propose that the shell and star both exhibit thermal emission, supported by the its appearance in near-infrared and H$\alpha$ observations. The latest \textit{Gaia} parallax is used to determine a distance of 2.28$\pm$0.09\,kpc. This star is well-known for its surrounding circular nebulosity, and using the distance and an angular diameter of 8\farcm42, we determine the shell size to be 5.57$\pm$0.22~pc. We use the \textit{Gaia} \acp{PM} of WR16 to determine peculiar velocities of the star as $V_{\alpha}(pec) =$ --45.3$\pm$5.4\,\kms\ and $V_{\delta}(pec) =$ 22.8$\pm$4.7\,\kms, which indicates that the star is moving in a north-west direction, and translates to a peculiar tangential velocity to be 50.7$\pm$6.9\,\kms. We also use these \acp{PM} to determine the shell's origin, estimate an age of $\sim 9500\pm 1300$\,yr, and determine its average expansion velocity to be $280\pm40$\,\kms. This average expansion velocity suggests that the previous transitional phase is a \ac{LBV} phase, rather than a \ac{RSG} phase. We also use the measured flux at 943.5~MHz to determine a mass-loss rate of $1.753\times 10^{-5}~M_\odot~$yr$^{-1}$, and use this to determine a lower-limit on ionising photons of $N_{UV} > 1.406\times 10^{47}~$s$^{-1}$.