1. Introduction
Classical novae are eruption events on the surfaces of accreting white dwarf (WD) binaries (Gallagher & Starrfield Reference Gallagher and Starrfield1978; Bode & Evans Reference Bode and Evans2012; Warner Reference Warner1995). In a semi-detached binary system, Roche lobe overflow from a late type main sequence star accumulates on the WD surface, igniting a thermonuclear fusion reaction (see reviews by Starrfield, Truran, & Sparks Reference Starrfield, Truran and Sparks1978; Wolf et al. Reference Wolf, Bildsten, Brooks and Paxton2013). As a result, the temperature and density at the base of the accreted envelope rise, leading to a thermonuclear runaway that causes
$10^{-6}$
–
$10^{-4}\rm M_\odot$
of the accreted mass to be ejected along with heavier materials dredged up in the deeper layers of the WD (Gehrz et al. Reference Gehrz, Truran, Williams and Starrfield1998; José & Hernanz Reference José and Hernanz2007). The ejecta flows out at velocities ranging from 500 to 5000
$\rm km\:s^{-1}$
from the WD surface (Gallagher & Starrfield Reference Gallagher and Starrfield1978; Aydi et al. Reference Aydi2020a) in a roughly spherical shell (Seaquist & Palimaka Reference Seaquist and Palimaka1977) with accompanied asymmetries such as jets, rings, clumping, bipolarities (e.g., Sokoloski, Rupen, & Mioduszewski Reference Sokoloski, Rupen and Mioduszewski2008; O’Brien et al. Reference O’Brien, Davis, Bode, Eyres, Porter and Schilizzi2001; Weston et al. Reference Weston2016). This eruption event enriches the interstellar medium (Romano et al. Reference Romano, Matteucci, Zhang, Papadopoulos and Ivison2017) and has a typical energy of
$10^{38}$
–
$10^{43}$
$\rm ergs$
.
From weeks to years post initial optical brightening at mass ejection (Knigge, Baraffe, & Patterson Reference Knigge, Baraffe and Patterson2011; Shafter Reference Shafter2017), the remnants of the accreted envelope on the WD surface undergo nuclear burning (Shara, Prialnik, & Shaviv Reference Shara, Prialnik and Shaviv1977). As the resulting luminous ionising radiation from the WD diffuses through the expanding ejecta, the peak of the spectral energy distribution of this photospheric emission from within the envelope shifts to shorter wavelengths (Gallagher & Code Reference Gallagher1974), upto the supersoft X-ray regime through UV (see reviews by Chomiuk et al. Reference Chomiuk2021b; Della Valle & Izzo Reference Della Valle and Izzo2020; Poggiani Reference Poggiani2018). Non-thermal GeV gamma ray emission is produced due to relativistic acceleration of charged particles in hadronic or leptonic scenarios that play out between the primary and hypothesised delayed mass ejection in novae (see discussion by Chomiuk et al. Reference Chomiuk2014; Metzger et al. Reference Metzger2015; Vurm & Metzger Reference Vurm and Metzger2018; Acciari et al. Reference Acciari2022). Gamma-ray emission begins hours to days post-eruption and declines on a time-scale of days to weeks.
Radio emission from novae is produced by either free-free thermal radiation powered by the central WD (Seaquist & Palimaka Reference Seaquist and Palimaka1977; Hjellming et al. Reference Hjellming, Wade, Vandenberg and Newell1979; Taylor & Seaquist Reference Taylor and Seaquist1984; Bode & Evans Reference Bode and Evans2012), or non-thermal radiation from additional emission sources such as shocks (Kantharia Reference Kantharia2012). At first, the dense ejecta stay optically thick at all radio frequencies, and we detect optically thick thermal emission from the shell’s surface, resulting in a radio emission peak. As the ejecta density falls, the radio photosphere recedes through the nova shell, and we detect a mix of optically thick and optically thin thermal emission from deeper layers of the shell. When the photosphere within the envelope begins to shrink due to decreasing ejecta density, the lightcurve turns over and drops until the ejecta become optically thin throughout, resulting in relatively brighter emission at lower frequencies (Bode & Evans Reference Bode and Evans2012). Thermal emission traces the warm, ionised, expanding ejecta as it remains optically thick to thermal emission at radio frequencies over days to years. Non thermal emission is produced by electron acceleration at shock fronts that emit synchrotron radiation (Vlasov, Vurm, & Metzger Reference Vlasov, Vurm and Metzger2016; Weston et al. Reference Weston2016). Shock heated thermal gas also contributes to thermal emission, although likely in a negligible way (Metzger et al. Reference Metzger2014; Sokolovsky Reference Sokolovsky2023).
Radio lightcurves of novae typically brighten after the optical peak on timescales ranging from a few weeks to a few years, with single or double peaks following from high to low frequencies (Chomiuk et al. Reference Chomiuk2021b). Following the spectral development of novae at high (Bode & Evans Reference Bode and Evans2012; Ribeiro et al. Reference Ribeiro2014) and low frequencies (Kantharia Reference Kantharia2012) shows that the radio emission from novae is consistent with multi-component emission mechanisms: with a free-free thermal component dominant at higher frequencies and a non-thermal synchrotron component, subject to variable absorption dominant at lower frequencies (Steinberg & Metzger Reference Steinberg and Metzger2020; de Ruiter et al. Reference de Ruiter2023). This multi-component emission mechanism is also supported in the population study of brightness temperatures in radio detected novae by Chomiuk et al. (Reference Chomiuk2021b).
Radio studies of novae have helped to model the expanding ejecta and its morphology (e.g., Hjellming et al. Reference Hjellming, Wade, Vandenberg and Newell1979; Weston Reference Weston2016; Taylor et al. Reference Taylor, Hjellming, Seaquist and Gehrz1988), yielding measurements for distances (e.g., V959 Mon; Linford et al. Reference Linford2015), shell masses (e.g., T Pyx; Nelson et al. Reference Nelson2014) and kinetic energies (e.g., Nova Aquilae; Snijders et al. Reference Snijders1987). Following the discovery of GeV gamma-ray emission from novae (e.g., Abdo et al. Reference Abdo2010; Ackermann et al. Reference Ackermann2014; Gordon et al. Reference Gordon2021a; Franckowiak et al. Reference Franckowiak, Jean, Wood, Cheung and Buson2018), an increasing number are being targeted at radio wavelengths in large-scale (e.g., Chomiuk et al. Reference Chomiuk2021b; Bode, Seaquist, & Evans Reference Bode, Seaquist and Evans1987) and individual follow-up campaigns focused on novae with high optical brightness (e.g., V1500 Cyg; Hjellming et al. Reference Hjellming, Wade, Vandenberg and Newell1979), gamma-ray observations (e.g., V959 Mon; Chomiuk et al. Reference Chomiuk2014, V1324 Sco; Finzell et al. Reference Finzell2018), and recurrent nova status (e.g., U Sco; Chomiuk et al. Reference Chomiuk2021b).
There is still a gap between nova theory (Yaron et al. Reference Yaron, Prialnik, Shara and Kovetz2005) and observation (see review by Roy et al. Reference Roy2012), particularly in terms of the observed and theoretical ejecta masses in different novae sub-types. Similarly, the mechanism of mass ejection is not completely understood, and often models assuming a single impulsive burst (Prialnik Reference Prialnik1986) in a spherical outflow is too simple a geometry to explain observed asymmetries (Weston et al. Reference Weston2016). Radio observations are ideal tracers of the ejecta mass because the receding radio photospere through the expanding ejecta samples the entire ejecta mass profile by interacting with it. The main reason for the uncertainty in predicting accurate ejecta masses is that the observed population of novae is biased towards the optically brightest, and, recently, gamma-ray loud novae, hence, not representative of the whole population. Furthermore, only a small subset of novae have been observed in the radio bands, mostly as multi-wavelength follow-up. To increase the accuracy of ejecta mass prediction and comparison with theoretical estimates, we need unbiased radio observations.
The Australian Square Kilometre Array Pathfinder (ASKAP) is a widefield radio telescope that images large areas of the radio sky at frequencies between 0.8 and 1.7 GHz (Hotan et al. Reference Hotan2021a). The surveying capability of ASKAP allows for unbiased high cadence radio observations of novae at currently undersampled, lower radio frequencies. In this paper, we use the datasets from multi-epoch widefield surveys, such as the Variable and Slow Transients Pilot Survey (VAST; Murphy et al. Reference Murphy2021) and Rapid ASKAP Continuum survey (RACS; McConnell et al. Reference McConnell2020), to investigate the utility of ASKAP for detecting novae. In Section 2, we describe the nova sample, data acquisition, relevant methodology, and search methods. In Section 3, we report the four novae detected in our search. In Section 4, we discuss the novae detectability in current and future ASKAP surveys; and conclude in Section 5.
Surveys conducted using ASKAP, included in the search for radio emission from novae in this work and the corresponding estimated integration times of the observations,
$t_{int}$
. Figure 1 shows the survey sky coverages.
$^{\rm a}$
Murphy et al. (Reference Murphy2021)
$^{\rm b}$
McConnell et al. (Reference McConnell2020)
$^{\rm c}$
Meyer (Reference Meyer2009)
$^{\rm d}$
Norris et al. (Reference Norris2021)
$^{\rm e}$
Allison et al. (Reference Allison2022)
$^{\rm f}$
Dickey et al. (Reference Dickey2013)
$^{\rm g}$
Gaensler et al. (Reference Gaensler, Landecker, Taylor and POSSUM2010)
$^{\rm h}$
Koribalski et al. (Reference Koribalski2020)
$^{\rm i}$
15 observation epochs from the VAST survey have been included in this search.
2. Data and methodology
2.1. ASKAP data
ASKAP is a widefield radio survey instrument made up of an array of 36 prime focus antennas, located at the CSIRO Murchison Radio-astronomy Observatory in Western Australia (Johnston et al. Reference Johnston2007; Hotan et al. Reference Hotan2021b). Over 2020–2022 a number of pilot surveys were carried out with ASKAP, and we have used these as the basis of the work presented here. From 2023, these surveys will continue as ASKAP commences full operations.
Positions of novae in both samples in equatorial coordinates in relation to the ASKAP survey footprints. Circle markers show the distribution of Sample 1. Star markers show the distribution of Sample 2. Diamond markers show the positions of novae having radio emission detected at their positions, labeled by source name. Filled-in background regions depict the ASKAP pilot surveys searched for radio emission from novae, as labeled. The gray-shaded background shows the Milky Way dust distribution map from the Planck Collaboration et al. (2014), where the light to dark shading corresponds to extinction magnitudes,
$A_V$
, ranging from
$3.61\times10^{-3}$
to
$1.36\times10^{3}$
mag respectively.
We searched for novae radio counterparts in ASKAP data from RACS, VAST, and other surveys listed in Table 1. Except for RACS, all surveys listed in Table 1 have been in the pilot phase so far, and will be continued and expanded in ASKAP’s full surveying mode starting in 2023. RACS covers the entire sky south of declination
$\sim$
$40^{\circ}$
(McConnell et al. Reference McConnell2020). In this search, we used the first epoch of RACS at 0.888 GHz, RACS-low, conducted in 2019 (declination of
$-80^{\circ}$
to
$41^{\circ}$
; McConnell et al. Reference McConnell2020). We also used the second epoch of RACS at 1.367 GHz, RACS-mid, conducted in 2021 (declination of
$-90^{\circ}$
to
$49^{\circ}$
; Duchesne et al. submitted) to get an additional measurement for two of our candidates—V1369 Cen and YZ Ret. The VAST project conducted
$\sim$
162 h of observations in its pilot survey phase. These were spread over multiple epochs, with the goal of detecting transient astronomical phenomena on timescales ranging from 10 s to several years (Murphy et al. Reference Murphy2013). We used data from 13 epochs of the VAST pilot survey Phase 1 (VAST-P1) and two epochs of the Phase 2 (VAST-P2 (low)) observations at 0.888 GHz from 2019 August to 2021 August (Murphy et al. Reference Murphy2021). The VAST and RACS surveys have integration times of 12–15 min over a total field of view of
$\sim$
5000
$\rm deg^2$
and
$\sim$
34000
$\rm deg^2$
respectively. All the other pilot surveys listed in Table 1 have much longer integration times,
$t_{int}$
, for each field over smaller regions.
ATCA and ASKAP observations for ASKAP-detected novae in this targeted search. Column 1 gives the GCVS ID of the nova detected; Column 2 gives the GCVS class (Duerbeck Reference Duerbeck1981) or variability type of the novae; Columns 3 and 4 give the optical coordinates of the nova; Column 5 gives the nova discovery date,
$t_0$
, in the optical band; Column 6 gives the radio observation date; Column 7 gives the time of the observations in days post-discovery in optical,
$\Delta T$
; Column 8 gives the radio telescope or survey used for the observations; Column 9 gives the central frequency of the observations,
$\nu$
; and Column 10 gives the catalogued integrated flux,
$S_{int}$
, for detections, and
$3\sigma$
limits for non-detections.
2.2. Sample selection
We constructed two samples of novae to search for in the ASKAP surveys. For our first sample, we selected novae from the list maintained by Bill Gray under Project Pluto.Footnote a We included classical novae; objects that were classified as fast novae (NA), slow novae (NB), extremely slow novae (NC), recurrent novae (NR), and those that have not yet been assigned a confirmed classification (for example, NC/ZAND, UG/N?). We excluded objects classified as symbiotic Z Andromedae systems (Z And class), U Gem variables or dwarf novae, X-ray nova systems and variable star types—Mira, UV Ceti, FU Ori and S Dor class. As of 2021 July, there are 470 identified classical novae, recurrent novae, and previously unclassified candidate novae systems discovered at optical wavelengths in the Milky Way. Four hundred and forty of those 470 sources occured within the region covered by one or more ASKAP surveys. These 440 objects constitute Sample 1 of our targeted search.
For our second sample we used the dataset of radio observations compiled by Chomiuk et al. (Reference Chomiuk2021b, Reference Chomiuk2022). This consists of 36 classical or not yet confirmed recurrent novae discovered after 1980, and well-sampled at radio wavelengths. Most of the data is from published and archival unpublished very large array (VLA) datasets, as well as new observations made with the Australian Telescope Compact Array (ATCA). Figure 1 shows the positions of novae from Samples 1 and 2 in relation to the ASKAP pilot survey footprints.Footnote b For each nova in the samples above, we recorded the discovery date from literature as
$t_0$
(refer to Table 2), and plot all times relative to that
$t_0$
as
$\Delta T$
in Figure 2. We additionally recorded the distances, optical variability type and previous multi-wavelength observations from relevant literature for each nova in Section 3.2.
2.3. Search method
We searched for radio counterparts by crossmatching (search radius 5′′) the proper motion corrected optical positions from SIMBAD (corrected to date: 2020 July 1) for novae in the two samples with sources catalogued in the surveys in Table 1. We extracted radio data, images, and flux density estimations at the optical nova positions for VAST using VAST Tools.Footnote c We obtained the source finder output catalogues for the other surveys from the CSIRO ASKAP Science Data Archive (CASDA).Footnote d The signal to noise ratio for detections in the main catalogues is
$5\sigma$
. We also searched for radio counterparts at a lower signal to noise ratio of
$3\sigma$
by measuring the flux density and background noise at each nova position in the extracted image cutouts.
We examined the positional offsets between our candidate radio counterparts and their matched proper-motion-corrected optical progenitor positions to determine whether our potential candidates were genuine counterparts. We also looked for previous radio emission in surveys such as the Very Large Array Sky Survey (Gordon et al. Reference Gordon2021b; Lacy et al. Reference Lacy2020), the National Radio Astronomy Observatory VLA Sky Survey (Condon et al. Reference Condon1998), the Sydney University Molonglo Sky Survey (Mauch et al. Reference Mauch2003), the second epoch Molonglo Galactic Plane Survey (Murphy et al. Reference Murphy2007), The GMRT 150 MHz All-Sky Radio Survey alternative data release 1 (Intema et al. Reference Intema, Jagannathan, Mooley and Frail2017) and The Galactic and Extra-Galactic All-Sky MWA (Hurley-Walker et al. Reference Hurley-Walker2017). Where radio data were available, we checked to see if the potential candidate’s radio flux density in the ASKAP observations was consistent with expected nova evolution (i) in time, i.e., radio lightcurves from earlier radio observations, and (ii) in frequency, i.e., spectral energy distributions from concurrent multi-frequency observations per the evolution of the radio spectral index,
$\alpha$
(defined as
$S_{\nu} \propto \nu^{\alpha}$
). Bode & Evans (Reference Bode and Evans2012) suggests a spectral index,
$\alpha$
, approaching 2 for rising optically thick thermal bremsstrahlung and
$\alpha\sim -0.1$
for declining optically thin thermal emission. A steep spectrum of
$\alpha\sim-0.75$
is typical for optically thin synchrotron emission (observed in for example, V3885 Sgr; Körding et al. Reference Körding, Knigge, Tzioumis and Fender2011, GK Persei; Reynolds & Chevalier Reference Reynolds and Chevalier1984, RS Oph; Kantharia et al. Reference Kantharia2007 and V1723 Aql; Weston et al. Reference Weston2016). The late-time transition from rising to a falling spectrum is characterised by one or two spectral breaks (Weston et al. Reference Weston2016).
We quantified the false-positive rate for chance alignment between optical nova positions and unrelated sources using a Monte Carlo analysis. We took the sample of 470 novae described in Section 2.2 and offset their positions by a random angle and a random distance between 1′-10′ and calculated the number of detections within a crossmatch radius of
$5^{\prime\prime}$
for each random offset of the original nova distribution in the RACS survey. Unlike the original distribution, each of the 100 offset distributions with 470 sources produced
$\leq1$
match (22 distributions had 0 matches, 78 distributions had 1 match) to the 2.1 million radio sources, and no distribution had five positive matches in RACS. This implies a false positive rate of
$4.8 \times 10^{-4}$
, indicating that the majority of the novae we detected in RACS are real.
Radio lightcurves for V5668 Sgr, V1369 Cen, YZ Ret, and RR Tel. For each source, the top panels show the optical light curve created using V-band or visual AAVSO measurements (Kafka Reference Kafka2021). The bottom panel for: V5668 Sgr shows the radio lightcurve created using ASKAP data observed with the RACS survey (
$\alpha_{1.26-1.74\rm\,GHz} = 1.65\pm0.55$
;
$\alpha_{1.74-7.0\rm\,GHz} = 0.07\pm0.02$
for the grey stripe) and observations from Chomiuk et al. (Reference Chomiuk2022); V1369 Cen shows the radio lightcurve created using ASKAP and ATCA observations from this work highlighted in the blue striped region, and observations from Chomiuk et al. (Reference Chomiuk2022) (
$\alpha_{2.9-9.0\rm\,GHz}= -1.00\pm0.05$
for the grey stripe); YZ Ret shows the radio lightcurve (time in linear space) created using ASKAP and ATCA observations from this work –
$\alpha_{1.7-5.5\rm\,GHz} = 0.85\pm0.32$
;
$\alpha_{5.5-9.0\rm\,GHz} = -0.15$
(2022 Feb 6);
$\alpha_{5.5-21.2\rm\,GHz}=-0.94\pm0.10$
(2022 Mar 10) and
$\alpha_{5.5-16.7\rm\,GHz}= -1.00\pm0.15$
(2022 Mar 18); and RR Tel shows the radio lightcurve created using ASKAP data observed with VAST and RACS surveys and observations from Seaquist (Reference Seaquist1977), Wright & Allen (Reference Wright and Allen1978), Purton et al. (Reference Purton, Feldman, Marsh, Allen and Wright1982), Ivison et al. (Reference Ivison, Seaquist, Schwarz, Hughes and Bode1995), Dickey et al. (Reference Dickey, Weston, Sokoloski, Vrtilek and McCollough2021), Mauch et al. (Reference Mauch2003). Downward triangle markers are 3
$\sigma$
detection limits.
3. Results
We found candidate radio counterparts for six novae from Sample 1. All of these were detected in RACS, and two (RR Tel and V5590 Sgr) were also detected in VAST. Two of these six novae (V1369 Cen and V5668 Sgr) are also in Sample 2. We found no detections in the WALLABY, GASKAP, FLASH, DINGO, EMU, and POSSUM surveys.
After further investigation we confirmed radio counterparts for four objects: the extremely slow symbiotic nova RR Tel and the classical novae V1369 Cen, V5668 Sgr, and YZ Ret. All confirmed radio counterparts were spatially consistent with their proper motion-corrected optical progenitor positions within astrometric errors. The conservative astrometric error for a weak source (
$6\sigma$
) in RACS is of the order
$\sim 2^{\prime\prime}.5$
(see McConnell et al. Reference McConnell2020). Two of the counterparts—those for V1280 Sco and V5590 Sgr—are uncertain, pending further multifrequency radio observations. These are discussed in Appendix A.
Our ASKAP and ATCA observations are summarised in Table 2, and discussed below.
3.1. Follow-up ATCA observations
We observed the radio counterpart of ASKAP-detected nova V1369 Cen on 2022 February 6, with the ATCA (Project code C3363; PI T. Murphy), measuring the flux density at a wide range of frequencies to determine its evolution phase. In addition, based on a tentative
$\geq3\sigma$
detection (Lenc, private communication) in preliminary RACS-mid data, we observed YZ Ret on 2022 February 6, March 10, and March 18 to confirm the ASKAP detection. All our observations used 2048 MHz-wide bands centered on 2.1, 5.5, 9.0, 16.7, and 21.2 GHz, except the last epoch that had a zoom mode observation bandwidth of 58 MHz in the L-band.
We reduced the visibility data using standard routines in MIRIAD (Sault, Teuben, & Wright Reference Sault, Teuben, Wright, Shaw, Payne and Hayes1995). We used a combination of manual and automatic RFI flagging before calibration, conducted with MIRIAD tasks uvflag and pgflag, respectively. We used PKS 1934
$-$
63 to determine the bandpass response for all frequency bands in all epochs, except the 16.7/21.2 GHz bands. For these two exceptions, we used PKS 0727
$-$
115 as the bandpass calibrator for the first epoch and PKS 1921
$-$
293 for the second epoch. We used PKS 1934
$-$
63 to calibrate the flux density scale for all epochs and frequency bands for both sources. For YZ Ret, we used PKS 0334
$-$
546 to calibrate the time-variable complex gains for all epochs and frequency bands. For V1369 Cen, we used IVS B1338
$-$
582 to calibrate the gains for the 2.9/5.5/9.0-GHz bands and IVS B1325
$-$
558 for the 16.7/21.2-GHz bands. After calibration, we inverted the visibilities using a robust weighting of 0.5 and then used the CLEAN algorithm (Clark Reference Clark1980) with 250 iterations to the target source field using standard MIRIAD tasks INVERT, CLEAN and RESTOR to obtain the final images.
3.2. Notes on individual objects
We discuss each source briefly below. Figure 2, shows the lightcurve for each object, with all available radio data.
3.2.1. V5668 Sgr
V5668 Sgr is a slow-evolving dust-forming classical nova (Harvey et al. Reference Harvey2018; Chomiuk et al. Reference Chomiuk2021b) discovered by J. Seach on 2015 March 15.634 UT, at a distance of
$2.8\pm0.5$
kpc (Gordon et al. Reference Gordon2021a). It is well-sampled at centimetre radio wavelengths with the VLA and ATCA (Chomiuk et al. Reference Chomiuk2021b), and has near-simultaneous detections in the RACS survey with ASKAP (see Figure 2 top panel left image). Ejecta morphology and dust formation have been studied for V5668 Sgr at mid- (Gehrz et al. Reference Gehrz2015) and near-infrared (Banerjee et al. Reference Banerjee, Srivastava, Ashok and Venkataraman2016), X-ray (Gehrz et al. Reference Gehrz2018) and millimetre wavelengths (Diaz et al. Reference Diaz, Abraham, Ribeiro, Beaklini and Takeda2018). The radio counterpart of V5668 Sgr in RACS is offset by 2′′.8 from its optical position.
We measured a 0.888 GHz flux density of
$1.56\pm0.06$
mJy for V5668 Sgr (1505 d post-discovery in optical). This value agrees with the unbroken power-law extrapolation of flux densities measured in near-simultaneous, higher frequency radio observations with the ATCA at 1500.5 d by Chomiuk et al. (Reference Chomiuk2021b). For the near-simultaneous multifrequency ASKAP and ATCA detections at this epoch, the peak radio flux density is somewhere between 1.26 and 5 GHz (
$\alpha_{1.26-1.74\rm\,GHz} = 1.65\pm0.55$
;
$\alpha_{1.74-7.0\rm\,GHz} = 0.07\pm0.02$
), which indicates that the emission is a mix of optically thick and optically thin thermal emission and at the point of spectral turnover optical depth,
$\tau=1$
. This characteristic indicates that the classical nova V5668 Sgr was in its declining radio evolution phase at the time of our observations.
3.2.2. V1369 Cen
V1369 Cen was discovered in an eruption on 2013 December 2.69 UT by J. Seach at a distance of
$1.0\pm0.1$
kpc (Gordon et al. Reference Gordon2021a). It was detected with the ATCA between 2013 December 5 and 2014 May 27 in multi-frequency observations (5.5–34 GHz) by Bannister et al. (Reference Bannister, Chomiuk, Rupen and Sokoloski2014), Chomiuk et al. (Reference Chomiuk2021b). It is characterised as a fast radio classical nova by Chomiuk et al. (Reference Chomiuk2021b). Despite multiwavelength studies of the object at near-IR (Mason et al. Reference Mason2021), UV (Shore et al. Reference Shore2014), gamma ray (Cheung et al. Reference Cheung2016), and X-ray wavelengths (Gordon et al. Reference Gordon2021a), radio observations of V1369 Cen around its thermal maximum are limited. The radio counterpart of V1369 Cen in RACS is offset by 0′′.9 from its optical position.
Hubble Flow model fit to V1369 Cen radio detections. Minimization of
$\chi^2$
between 0.9, 5, and 9 GHz radio detections, given by filled-in circles, and models with ejecta masses ranging between
$5\times10^{-6}\: \rm M_\odot$
and
$1\times10^{-3}\: \rm M_\odot$
, results in the fit in the teal and yellow solid line fits respectively.
We measured a 0.888 GHz flux density of
$2.04\pm0.04$
mJy for V1369 Cen. We also measured a
$5\sigma$
detection for the source in the RACS-mid survey (1.36 GHz). For our 2022 February 6 ATCA follow-up observation, we measured a
$\geq10\sigma$
detection at 2.1/5.5/9 GHz and non-detections at 16.7/21.2 GHz (Gulati et al. Reference Gulati, Murphy, Wang, Leung and Lenc2022a). In this epoch, the steep spectral index,
$\alpha_{2.9-9.0\rm\,GHz}= -1.00\pm0.05$
, is indicative of synchrotron emission. Peaking of synchrotron radiation, dominant at lower frequencies due to decreasing free–free emission, suggests that the ejecta are optically thin across the observed bands. These characteristics indicate that classical nova V1369 Cen was in its declining radio evolution phase at the time of our observations (see Figure 2 top panel right image).
We fit a Hubble Flow modelFootnote e (described in Hjellming et al. Reference Hjellming, Wade, Vandenberg and Newell1979; Hjellming Reference Hjellming, Taylor and Paredes1996; Seaquist & Palimaka Reference Seaquist and Palimaka1977; Seaquist et al. Reference Seaquist1980) to multi-frequency radio detections (0.9, 5 and 9 GHz) for V1369 Cen obtained in this work as well as previous works via
$\chi^2$
Footnote f minimisation between the data and models with ejecta masses ranging between
$5\times10^{-6}\: \rm M_\odot$
and
$1\times10^{-3}\: \rm M_\odot$
(see Figure 3). We assumed that the mass was ejected at the time of the optical discovery. We used model parameters previously derived in literature for V1369 Cen, such as velocities of the inner edge,
$v_1=1150\, \rm km\:s^{-1}$
, and outer edge,
$v_2=1700\,\rm km\:s^{-1}$
, of the expanding shell (Chomiuk et al. Reference Chomiuk2021b), a constant electron temperature
$T_e=10^4$
K (Osterbrock Reference Osterbrock1974), and assumed pure hydrogen composition (
$\rho$
/
$N_e$
=
$1.67 \times 10^{-24}$
g). The best fit ejecta mass was derived from this fit (
$\chi^2=1234$
, DOF=8) is
$1.65\pm0.17\times 10^{-4}\: \rm M_\odot$
. When compared to theoretical predictions by Yaron et al. (Reference Yaron, Prialnik, Shara and Kovetz2005), our derived ejecta mass places the mass of the WD that hosts V1369 Cen’s explosion in the range 1–1.25
$\rm M_\odot$
.
3.2.3. YZ Ret
YZ Ret is one of the brightest classical novae (
$V_{peak}=3.7$
mag), identified by R. H. McNaught (McNaught Reference McNaught2020) on 2020 July 15.590 UT at a distance of
$2.7\substack{+0.4 \\ -0.3}$
kpc (Bailer-Jones et al. Reference Bailer-Jones, Rybizki, Fouesneau, Mantelet and Andrae2018). It is associated with known cataclysmic variable MGAB-V207 and had not been reported as radio-loud until this search.
Following previous multi-wavelength observations (Aydi et al. Reference Aydi2020b; Li et al. Reference Li2020; Drake et al. Reference Drake2020; Sokolovsky et al. Reference Sokolovsky2020; Izzo et al. Reference Izzo2020), we measured a
$\geq14\sigma$
detection for YZ Ret on 2022 February 6 at frequencies 2.9/5.5/9 GHz with the ATCA (Gulati et al. Reference Gulati2022b). In this epoch, the peak radio flux density is somewhere between 2.9 and 9 GHz (
$\alpha_{1.7-5.5\rm\,GHz} = 0.85\pm0.32$
;
$\alpha_{5.5-9.0\rm\,GHz} = -0.15$
), which indicates that the emission is a mix of optically thick and optically thin thermal emission. Using relations derived by Chomiuk et al. (Reference Chomiuk2021b), we calculated the size of the radio photosphere to be
$0^{\prime\prime}.65\pm 0^{\prime\prime}.1$
and the corresponding surface brightness temperature to be
$250\pm80$
K. The surface brightness temperature was calculated at the peak of the 5 GHz observations and an outer edge expanding ejecta velocity,
$v_{2}$
, of 2700
$\,{\rm km\,s}^{-1}$
(Aydi et al. Reference Aydi2020b). YZ Ret’s surface brightness temperature is within the conservative temperature range for thermal radio emission (
$\lesssim$
$10^4$
K), supporting the presence of optically thick thermal emission suggested in the spectral analysis of multi-frequency observations in this epoch. However, due to a lack of radio observations prior to our observations, this surface brightness temperature is not calculated at a known thermal maximum, and is insufficient to rule out the presence of non-thermal radio emission from YZ Ret. Thermal emission mechanisms are also suggested for YZ Ret in its study at X-ray wavelengths (Sokolovsky et al. Reference Sokolovsky2022).
For our 2022 March 10 follow-up observation, we measured a
$\geq$
$40\sigma$
detection at 5.5/9 GHz and a
$\geq7\sigma$
detection at 16.7/21.2 GHz; for our 2022 March 18 observation, we measured a
$\geq$
$45\sigma$
detection at 5.5/9 GHz, a
$11\sigma$
detection at 16.7 GHz, and no detection at 21.2 GHz (Gulati et al. Reference Gulati2022c). We measured steep spectral indices of
$\alpha_{5.5-21.2\rm\,GHz}=-0.94\pm0.10$
and
$\alpha_{5.5-16.7\rm\,GHz}= -1.00\pm0.15$
for the second and third epochs respectively. The fading of emission at higher frequencies, and the steep spectral index observed at lower frequencies indicates that free-free thermal emission component dominant at higher frequencies has decayed and only the synchrotron emission dominant at lower frequencies is being detected. This suggests that the ejecta are optically thin across the observed bands. Furthermore, the steepening of the spectral index from the second to third epoch suggests energy losses by the synchrotron (Pacholczyk Reference Pacholczyk1973). These characteristics indicate that the classical nova YZ Ret was in its declining radio evolution phase at the time of our observations (see Figure 2 bottom panel left image).
3.2.4. RR Tel
RR Tel is a symbiotic nova that evolves at an extremely slow rate (Nussbaumer & Dumm Reference Nussbaumer and Dumm1997) over a period of years to centuries. The analysis of X-ray observations for the nova by González-Riestra, Selvelli, & Cassatella (Reference González-Riestra, Selvelli and Cassatella2013) also suggest slow evolution, with the outflow from the WD being a direct result of the eruption even 65 yr after the burst in 1944. The radio counterparts of RR Tel in VAST and RACS are offset by 0′′.6
$-$
1′′.1 from its optical position.
Over a year of ASKAP radio observations, RR Tel does not show any significant increase or decrease, as shown by the measured 0.888 GHz flux densities ranging from 9 to 10 mJy (see Figure 2 bottom panel right image). Nonetheless, when compared to RR Tel detections with the ATCA in 2007 at a 8.6 GHz flux density of
$16.2\pm2.3 $
mJy and a 4.8 GHz flux density of
$17.6\pm0.8$
mJy (Dickey et al. Reference Dickey, Weston, Sokoloski, Vrtilek and McCollough2021), and previous radio observations at frequencies ranging from 5.0 to 8.9 GHz (Seaquist Reference Seaquist1977; Wright & Allen Reference Wright and Allen1978; Purton et al. Reference Purton, Feldman, Marsh, Allen and Wright1982; Ivison et al. Reference Ivison, Seaquist, Schwarz, Hughes and Bode1995; Mocanu et al. Reference Mocanu2013), there is a visible trend of prolonged fading in the radio band.
Summary of the properties of six slowly evolving radio classical novae from Sample 2 observed near-simultaneously with ASKAP and ATCA/VLA
$^{\rm a}$
. Column 1 gives the GCVS ID of the nova that has simultaneous radio observations with ASKAP; Column 2 gives the number of days between discovery in optical and the ASKAP observation,
$\Delta T$
; Column 3 gives the number of days between ATCA or VLA observations at
$t_{V}$
days post-discovery, and ASKAP observations at
$t_{A}$
days post-discovery; Column 4 gives the observed ASKAP flux density,
$S_{A}$
; Column 5 gives the estimated flux density at 0.888 GHz from the spectral fitting,
$S_{est}$
; Column 6 gives the spectral index,
$\alpha$
, for
$\sim$
1–9 GHz; and Column 7 gives the nova evolution phase at the time of observation.
aATCA/VLA data available in Chomiuk et al. (Reference Chomiuk2022).
bThe spectral indices for V5668 Sgr correspond to 1.26–1.74 and 1.74–7 GHz respectively.
Simulated novae light curves at 1 GHz using the Hubble Flow model with assumed pure hydrogen composition. The model used average parameters derived in literature: maximum ejecta velocity,
$v_2$
, from 500 to 3000
$\rm km\: s^{-1}$
, ejected mass,
$M_ej$
, from
$5\times10^{-6}$
to
$5\times10^{-4}\: \rm M_\odot$
, the ratio of minimum to maximum ejecta velocity,
$v_1/v_2$
, as 0.5, temperature of ejected material,
$T_e$
, as
$10^4$
K. The legend describes the combination of
$v_2$
and
$M_{ej}$
used for each curve. The nova distance is given in each plot. The VAST/RACS and the GASKAP 5
$\sigma$
sensitivity limits are given by the black and red dotted-dashed lines, respectively.
4. Detectability of classical novae with ASKAP
Our results from a targeted search demonstrated we can detect novae in a single epoch of the RACS survey, and in the VAST pilot surveys. In this section, we use previous novae radio studies in conjunction with ASKAP data to investigate the detectability of classical novae in ASKAP’s full surveying mode (across a frequency range of 0.888–1.65 GHz).
We selected six novae from Sample 2 that had near-simultaneous ASKAP and ATCA/VLA multi-frequency observations. Five of the six novae detected in radio at higher frequencies—V357 Mus, V906 Car, V5855 Sgr, V5856 Sgr, and V5667 Sgr—were undetected with ASKAP at 0.888 GHz (see Table 3). The spectral index scaling of near-simultaneous multi-frequency observations yields flux density estimates,
$S_{est}$
, in the range
$\sim$
0.2–0.3 mJy at 0.888 GHz, which are below ASKAP’s VAST/RACS
$5\sigma$
sensitivity threshold of 1.5 mJy (median RMS 0.326 mJy
$\rm{beam^{-1}}$
). Due to lack of flux-density sensitivity for the limited radio observations at ASKAP frequencies, we simulated radio lightcurves for novae at 1 GHz using average parameters for ejected mass and velocity from Chomiuk, Metzger, & Shen (Reference Chomiuk, Metzger and Shen2021a) in the Hubble Flow model. These simulated curves are presented in Figure 4 in relation with RACS/VAST and a deeper Galactic survey, GASKAP’s 5
$\sigma$
sensitivity limits.
As shown in the plots, the greater the mass ejected, the longer the nova eruption would be visible above the VAST or GASKAP 5
$\sigma$
sensitivity limits. At 0.8 kpc, novae with
$M_{ej}$
in the range
$5\times 10^{-6}\rm \:M_\odot$
are not visible above the VAST detection limit. As the distance of the explosion increases to 3 kpc, only explosions with
$M_{ej}\gtrsim 5\times 10^{-4}\rm \:M_\odot$
are detectable above the VAST sensitivity limit. With a 5
$\sigma$
sensitivity limit of 1.5 mJy in 15 min survey observations (RACS and VAST), we can detect radio emission up to a distance of 4 kpc if the ejecta mass is in the range
$10^{-3}\rm \:M_\odot$
, and upto 1 kpc if ejecta mass is in the range
$10^{-5}$
–
$10^{-3}\rm \:M_\odot$
. Since, the mass ejected is inversely proportional to the mass of the WD hosting the explosion, the detectability limit on the ejecta mass puts an upper limit of
$1.25\:\rm M_\odot$
and
$1\:\rm M_\odot$
on the masses of the white dwarfs that host ASKAP detectable novae explosions at 1 and 4 kpc respectively (Yaron et al. Reference Yaron, Prialnik, Shara and Kovetz2005). YZ Ret’s inferred WD mass of
$0.98\pm0.23\:\rm M_\odot$
in König et al. (Reference König2022) at 2.7 kpc is consistent with this ASKAP detectability limit. A higher 5
$\sigma$
sensitivity limit of 0.7 mJy in the GASKAP survey, focused on the Galactic plane, will be useful for detecting radio emission from nova explosions hosted on higher mass host white dwarfs (
$\gtrsim$
$1\,\rm M_\odot$
) that eject material
$\lesssim5\times10^{-6}\:\rm M_\odot$
.
Figure 4 also shows that as ejecta velocity,
$v_2$
, increases, the radio thermal maximum peaks earlier and spends less time above the 5
$\sigma$
sensitivity limits (see Table B.1). At 0.8 kpc, novae with
$M_{ej}$
in the range
$5\times 10^{-4}\rm \:M_\odot$
spend
$\sim$
3000–4000 d above the RACS/VAST 5
$\sigma$
limit where the greater number of days correspond to lower velocities (
$v_{2}\approx 500\: \rm km\:s^{-1}$
). As the distance of explosion increases to 3 kpc, for the same
$M_{ej}$
, the time spent above the RACS/VAST 5
$\sigma$
limit is
$\sim$
700–2000 d. For faster novae (
$v_{2}\approx 3000\,\rm km\:s^{-1}$
), we would need one epoch of radio sky surveying per 2 months since the minimum time detectable above the higher sensitivity GASKAP 5
$\sigma$
limit in the curves simulated for the lowest
$M_{ej}\approx 5\times 10^{-6}\rm \:M_\odot$
is 72 d (see Table B.1). VAST surveys, with their multi-epoch surveying strategy covering 265
$\rm deg^2$
on the Galactic centre, would be the best ASKAP survey to capture fast novae. The Galactic centre has high dust concentration along its line of sight (see Figure 1), implying pronounced reddening and higher optical extinction (Planck Collaboration et al. 2014). This region is, therefore, particularly suited for radio searches of novae with VAST.
Along with nova ejecta characteristics, and survey sensitivities, novae will of course only be detectable if they are within the sky coverage of a given survey. From the current pilot surveys listed in Table 1, the three surveys that cover the highest percentage of all previous nova eruptions are RACS (93.6%), VAST (31.9%) and GASKAP (14.1%). Single epoch ASKAP detections, from all-sky surveys such as RACS, will be useful for multi-frequency radio follow-up observations with other telescopes such as the ATCA, as demonstrated for the case of YZ Ret in Section 3.
Because of the irregular cadence for pilot survey observations used in this work, multiwavelength data was required to determine radio counterparts. Future ASKAP surveys will increase the number of novae with radio lightcurves, allowing spectral and temporal modeling. Radio detections could also overcome extinction effects and gaps in optical novae monitoring, which account for a large discrepancy of
$\sim$
30–50 novae between observed and theoretical discovery rates per year (Shafter Reference Shafter2017; De et al. Reference De2021). For example, VVV-WIT-01 was classified as a classical nova using only radio and IR multiwavelength analysis, with no detection in optical (Lucas et al. Reference Lucas2020). The availability of a statistically robust sample of novae observed in radio would also aid in more accurate ejecta mass predictions. While our targeted search did not uncover radio rebrightening from recurrent novae, data from routine ASKAP surveys will allow these experiments to enable insights about their ejecta masses and spectral signatures (for e.g., RS Oph; Cheung et al. Reference Cheung2022).
We have primarily considered the discovery rate of novae due to thermal emission, which can be faint at low radio frequencies (
$\lesssim$
1 GHz) even when not absorbed, as evidenced by the ASKAP non-detections for all novae in their optically thick thermal emission phase in Table 3. Furthermore, spectral index analysis for our detected novae in Section 3 reveals that when the ejecta is optically thin, we identify non-thermal synchrotron emission for V1369 Cen and YZ Ret. This is consistent with previous observations of novae systems at low frequencies (Kantharia Reference Kantharia2012; Weston et al. Reference Weston2016; de Ruiter et al. Reference de Ruiter2023). As long as the thermally emitting material is not optically thick, synchrotron emission can be brighter at lower frequencies, and is present in 25% of the radio novae examined in Chomiuk et al. (Reference Chomiuk2021b), with a hypothesised upper limit of 100%. ASKAP observations of radio emission from novae at lower frequencies could help place strict limitations on the fraction of novae that generate synchrotron radiation. This will advance our understanding of the luminosity and timing of non-thermal components in novae explosions and aid in understanding how frequently these contribute to the observed diversity in radio lightcurves and, thereby, reflect the shock energetics of novae. Low-frequency VAST lightcurves, combined with high-frequency observations, could prove useful for studying different regions and physics of the nova system (see review by Kantharia Reference Kantharia2012).
5. Conclusions
We have completed the first systematic search for radio counterparts of 440 classical novae previously identified in optical at centimetre wavelengths using RACS, VAST and other ASKAP pilot surveys. We identified three known radio classical novae, and one with no previous radio detection. The detection of V1369 Cen, V5668 Sgr, and YZ Ret, all in their declining optically thin phase, and RR Tel in its ongoing slow variability phase, demonstrates ASKAP’s sensitivity to classical nova detection. We conducted ATCA observations of V1369 Cen and YZ Ret which indicate the presence of optically thin synchrotron component dominant at lower frequencies from the computed spectral indices. The multi-frequency observations also helped fit an ejecta mass of
$1.65 \pm0.17\times 10^{-4}\: \rm M_\odot$
for V1369 Cen; and strengthen the hypothesis for optically thick thermal radio emission from YZ Ret in the first observed epoch.
The surveying nature of ASKAP observations eliminate the need for large-scale radio campaigns across decades as follow-up observations can be conducted for those detectable with ASKAP, for example, monitoring detections for novae with host WD masses in the range 0.4–1.25
$\rm M_\odot$
, that erupt within a distance of 1 kpc in 1.5 mJy (5
$\sigma$
) flux density sensitivity surveys. Single-epoch detections of previously identified optical novae with observed false-positive rates of
$4.8 \times 10^{-4}$
, combined with multi-frequency follow-up strategy using other radio telescopes, could enable observations of novae near thermal radio maximum, as well as potential easier discovery of recurrent novae. Multi epoch VAST surveys with a region covering the optically extinct Galactic centre line of sight is well suited for radio searches of highly reddened novae with ASKAP. Detection and classification of the radio emission mechanism at ASKAP frequencies could aid the study of relative timing and luminosity of the components in the multi-component emission mechanism in novae and help understand underlying shock energetics.
The extension and application of this search method to surveys conducted at higher ASKAP frequencies of 1.36 and 1.65 GHz with a more complete novae list that includes progenitors discovered in UV/IR searches (De et al. Reference De2021) and from space based optical telescopes (Mandel et al. Reference Mandel, Shara, Zurek, Conroy and van Dokkum2023) will help better understand the scope of novae detectability with ASKAP. The position of the survey region on the sky, survey cadence, and survey flux density sensitivity are essential factors to consider when designing campaigns for targeting novae in the future. These observations will provide an opportunity to improve nova population statistics and study associated nova physics. To better understand the implications of these results, an unbiased search for a year’s worth of ASKAP data can be conducted when the full surveys begin.
Acknowledgement
We thank the referee for their helpful comments. AG, JL and JP are supported by Australian Government Research Training Program Scholarships. DK is supported by NSF grant AST-1816492. RS acknowledges grant number 12073029 from the National Natural Science Foundation of China (NSFC). Parts of this research were conducted by the Australian Research Council Centre of Excellence for Gravitational Wave Discovery (OzGrav), project number CE170100004.
This scientific work uses data obtained from Inyarrimanha Ilgari Bundara/the Murchison Radio-astronomy Observatory. We acknowledge the Wajarri Yamaji People as the Traditional Owners and native title holders of the Observatory site. CSIRO’s Australian Square Kilometre Array Pathfinder radio telescope is part of the Australia Telescope National Facility. Operation of ASKAP is funded by the Australian Government with support from the National Collaborative Research Infrastructure Strategy. ASKAP uses the resources of the Pawsey Supercomputing Research Centre. Establishment of ASKAP, Inyarrimanha Ilgari Bundara, the CSIRO Murchison Radio-astronomy Observatory and the Pawsey Supercomputing Research Centre are initiatives of the Australian Government, with support from the Government of Western Australia and the Science and Industry Endowment Fund. The Australia Telescope Compact Array is part of the Australia Telescope National Facility which is funded by the Australian Government for operation as a National Facility managed by CSIRO. We acknowledge the Gomeroi people as the Traditional Owners of the Observatory site.
This research has made use of the SIMBAD database (Wenger et al. Reference Wenger2000), operated at CDS, Strasbourg, France, and Software: Astropy (Astropy Collaboration et al. Reference Astropy2013), MATPLOTLIB (Hunter Reference Hunter2007), NumPy (Harris et al. Reference Harris2020), VAST Tools (Murphy et al. Reference Murphy2021), MIRIAD (Sault et al. Reference Sault, Teuben, Wright, Shaw, Payne and Hayes1995).
Data Availability
The ASKAP data used in this paper (RACS and VAST) can be accessed through the CSIRO ASKAP Science Data Archive (CASDA) under project codes AS110 and AS107. The ATCA data used in this paper can be accessed through the Australia Telescope Online Archive (ATOA) under project code C3363.
Appendix A. Uncertain Nova Counterparts
In Section 3, we ruled two out of the six nova radio counterpart candidates as uncertain. For each candidate we provide an explanation; further observations would be required to confirm these as counterparts.
V1280 Sco is an extremely slow evolving classical nova that was discovered in an outburst on 2007 February 4.8 by Yamaoka et al. (Reference Yamaoka, Nakamura, Nakano, Sakurai and Kadota2007) at a distance of
$1.1\pm0.5$
kpc (Ness et al. Reference Ness2009). Its matched radio source in RACS detected at a
$6.5\sigma$
significance flux density of
$1.62\pm0.04$
mJy is 3′′.3 offset from its optical position at 4469 d post-discovery (
$t_\textrm{RACS}$
). An ejecta expansion velocity of
$\sim$
$350\pm160$
$\rm km\:s^{-1}$
(Chesneau et al. Reference Chesneau2012) gives an ejecta radius of
$0^{\prime\prime}.8\pm 0^{\prime\prime}.5$
. The lack of previously confirmed radio observations, a single epoch detection in RACS, inconsistency of ejecta radius with optical radio offset (2′′.5 differential) make it difficult to assert whether V1280 Sco’s radio counterpart is related to it, hence, we ruled its radio counterpart as uncertain.
V5590 Sgr is a candidate Type D symbiotic nova candidate (Mróz et al. Reference Mróz2014) announced in 2012 (Ayani, Fujii, & Maeno Reference Ayani, Fujii and Maeno2012; Nakano et al. Reference Nakano2012). Its matched radio source in VAST and RACS was not detected in all epochs of our search (see Figure A.1). Despite its spatial consistency with the nova optical position (1′′.5 offset), the near-detection limit 3
$-$
5
$\sigma$
detections, lack of previous multiwavelength observations and confirmed radio detections make it difficult to assert whether its radio counterpart is related to the nova or is some unrelated background source, hence, we ruled its radio counterpart as uncertain.
The ASKAP radio lightcurve at V5590 Sgr’s position. Filled-in circles indicate catalogued flux density within 5′′ of the nova optical position, downward triangles are
$3\sigma$
limits.
Appendix B. Simulated Radio Lightcurve Charcteristics
In Table B.1, we give the nova ejecta parameters from Figure 4 along with corresponding lightcurve characteristics.
Simulated nova radio lightcurve characteristics at 1 GHz. Column 1 gives the ejected mass,
$M_{ej}$
; Column 2 gives the maximum ejecta velocity,
$v_2$
; Column 3 gives the maximum radio flux density,
$S_{peak}$
; Column 4 gives the time of
$S_{peak}$
,
$t_{peak}$
; Columns 5 and 6 give the time spent above VAST and GASKAP flux density sensitivity limits,
$t_{VAST}$
and
$t_{GASKAP}$
.
Open access
