Free Access
Issue
A&A
Volume 595, November 2016
Article Number L9
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201629763
Published online 11 November 2016

© ESO, 2016

1. Introduction

Following the first detection of a supernova at radio wavelengths close to fifty years ago (SN 1970G; Gottesman et al. 1972), it was soon widely realized that multi-frequency radio light curves can yield information about the configuration, density, and prevailing conditions of the circumstellar medium (e.g. Chevalier 1982b; Weiler et al. 1986). For cases where the radio light curves are well sampled, this information in turn can be translated into mass loss rates from the progenitor as a function of time, thereby revealing the likely evolutionary stages the star must have passed through prior to explosion (e.g. Sramek et al. 1984; Williams et al. 2002; Kotak & Vink 2006). Additionally, carrying out very long baseline radio-interferometry (VLBI) on nearby (few Mpc) objects can spatially resolve the expanding ejecta, uncovering underlying asymmetries, potentially due to a binary companion (e.g. Marcaide et al. 1994; Bartel & Bietenholz 2003).

It is generally well established that the radio emission arises from a combination of synchrotron self-absorption and free-free absorption as the blast wave from the supernova encounters the surrounding medium ionized by the initial high energy flash from the explosion, but shaped by the progenitor system (Chevalier 1982a,b). At early epochs, the optical depth is large, suppressing the radio emission. As the shock moves outwards, the optical depth drops, and its frequency dependence results in the transition from optically thick to optically thin occurring first at higher frequencies. Post peak, the emission declines smoothly as a power-law, unless further dense material is encountered, perhaps from earlier mass loss episodes.

Although several hundred supernovae (SNe) have now been detected at radio wavelengths, the overwhelming majority of these are detections at early times. This can be attributed to a combination of intrinsic faintness and sensitivity of currently available instrumentation on the one hand, and the lack of circumstellar material (CSM) at large distances from the supernova that is dense enough to generate radio emission on the other.

Supernova 1978K (SN 1978K) in the late-type barred spiral galaxy NGC 1313 was only the second supernova to be detected and recognized as a supernova from its radio emission. and the first from its X-rays (Ryder et al. 1993). At peak, it would have been one of the brightest radio supernovae yet seen. The Cepheid-based distance to NGC 1313 is 4.61 ± 0.21 Mpc (Qing et al. 2015), although estimates based on the tip of the Red Giant Branch method (e.g. Jacobs et al. 2009) put it slightly closer. Regardless, SN 1978K is one of the closest supernovae undergoing ejecta-circumstellar interaction (Type IIn). It is one of a small number of SNe amenable to long term multi-wavelength monitoring, and it remains bright at all wavelengths almost four decades after explosion.

In this Letter we report new high frequency, and high spatial resolution imaging of SN 1978K with the Australia Telescope Compact Array (ATCA) and Long Baseline Array (LBA).

2. Observations and results

2.1. ATCA observations

SN 1978K was observed at 34 GHz and at 94 GHz with the ATCA in September 2014. The 34 GHz observations were performed with the array in the 6B configuration, and the 94 GHz observations with the array in the more compact H214 configuration. The observations were carried out in two bands, each of 2 GHz bandwidth, centred either on 33 and 35 GHz, or on 93 and 95 GHz, using the Compact Array Broad-band Backend (CABB: Wilson et al. 2011) with each 2 GHz band consisting of 2048 × 1 MHz channels. The parameters of each observing session are listed in Table 1.

Table 1

ATCA observations of SN 1978K.

The observations consisted of 12 min scans on SN 1978K, interleaved with 2 min scans on the phase calibrator PKS B0302-623 having flux densities of ~0.8 Jy at 34 GHz and ~0.6 Jy at 94 GHz. The sources PKS B2155-152 and PKS B1921-293 served as the bandpass calibrators at 34 GHz and at 94 GHz, respectively. Uranus was used as the primary flux density calibrator at both frequencies, with hourly pointing scans on the phase calibrator used to update the pointing model. At 94 GHz “paddle” scans preceded each phase calibrator observation to monitor the system temperatures.

The data were processed using the miriad1 package, and procedures outlined in Sect. 4.3 of the ATCA User Guide2. The 34 GHz visibility data were corrected for atmospheric opacity, while the 94 GHz visibility data had corrections applied to the antenna positions and system temperatures. After editing and calibration of the data, images at each frequency were made using a Briggs robust weighting factor of 0.5 as a best compromise between resolution and sensitivity, then cleaned down to three times the r.m.s. noise level. A clear detection of SN 1978K was obtained in each frequency, and fitting of a Gaussian point source yielded the flux densities shown in Table 1 where the uncertainty in each case is a combination of the fitting and absolute flux calibration errors. The resulting image of SN 1978K at 94 GHz is shown in Fig. 1. Core-collapse supernovae have been detected previously at these frequencies (e.g. SN 2011dh: Horesh et al. 2013), but usually only within a few days or weeks after explosion.

thumbnail Fig. 1

Continuum image of SN 1978K at 94 GHz from the ATCA on 2014 Sep. 29, showing a 10σ detection at the SN location (03h17m38.8s, − 66°3304′′ J2000). The synthesised beam is indicated in the lower-left corner and the apparent source elongation is an artifact of the underlying noise pattern.

In Fig. 2 we present the 1100 GHz spectral energy distribution (SED) of SN 1978K. The 110 GHz flux densities come from ongoing monitoring with the ATCA, the most recent epoch for which is 2013 June 7 UT. The expected flux densities at 1.4, 2.4, 4.8, and 8.6 GHz for 2014 September have been derived using the model light curve fits in Smith et al. (2007). While a power-law index of α = −0.6 (where flux density Sννα) provides a good fit to the data below 10 GHz, the power law steepens to α ~ −0.9 between the two higher frequency data points. We note that even with the spectral steepening at high frequencies, extrapolation of this power law to the regime of the Atacama Large Millimetre/submillimetre Array (ALMA) Bands 4 (144 GHz) through 7 (325 GHz) should make the continuum emission from SN 1978K easily detectable with ALMA. The Type IIn SN 1996cr (Meunier et al. 2013) in the Circinus Galaxy lies at a distance comparable to SN 1978K; it was found using the ATCA to have flux densities of 23.9 ± 6.0 mJy at 88.6 GHz at 3616 days after the explosion, and 40.3 ± 6.0 mJy at 34.5 GHz at 4842 days. The flux densities measured recently for SN 1987A at 44 GHz (Zanardo et al. 2013) and at 94 GHz (Lakićević et al. 2012) would, if it were placed at the distance of NGC 1313 instead of 51.4 kpc (Panagia 1999), correspond to ~5 μJy and ~3 μJy respectively, more than 2 orders of magnitude fainter than observed for SN 1978K.

thumbnail Fig. 2

Radio spectral energy distribution for SN 1978K in September 2014. The dashed lines mark a power law slope of − 0.6 passing through the ν< 10 GHz data, and a steeper power law with slope of − 0.9 passing through the ν > 10 GHz data. The red solid line shows an SED fit to the full dataset using the model of Weiler et al. (2002), which includes low frequency absorption due to an H ii region along the line of sight.

2.2. VLBI observations

SN 1978K was observed with the LBA on 2007 Nov. 13 and 2015 March 29 UT. These gave similar results and we discuss only the latter here because of the longer time since explosion. For this observation the array was comprised of 6 antennas: ATCA 5 × 22 m phased array; Parkes 64 m; Mopra 22 m; Ceduna 30 m; Hobart 26 m; and a single 12 m antenna of the Australian Square Kilometre Array Pathfinder (ASKAP) in Western Australia, offering baselines in excess of 3000 km. Although the observation block spanned 12 h, the source was only above the elevation limit of some antennas for ~10 h. A bandwidth of 64 MHz centred on 8.4 GHz was used in dual polarisation mode, and data recorded to disk. The sources PKS B1921-293 and PKS B0537-441 served as fringe-finders, while 3 min scans of SN 1978K were alternated with equivalent scans of the nearby source PKS B0355-669 (~0.5 Jy at 8.4 GHz) to allow phase referencing.

All the data were correlated using the DiFX software correlator (Deller et al. 2011), and the correlated data were exported to AIPS3 for editing, applying nominal sensitivities for each antenna, and iterative fringe-fitting. Imaging using standard self-calibration, robust=0 weighting, and clean techniques was carried out with the imagr task. The resultant image (Fig. 3) has a beam size of ~3 milli-arcsec and a dynamic range of ~100:1. Deconvolution of the VLBI source in the image plane using the jmfit task yields a deconvolved Gaussian size of 4.9 mas × 4.4 mas (± 0.5 mas in each axis) at a position angle of (20 ± 40)° compared to the beam position angle of 85°. The integrated flux density of 5.2 ± 0.3 mJy is just half the ATCA-measured value at 9 GHz in Fig. 2, suggesting that some of the more-extended emission has been resolved out on account of the longer baselines and sparse uv-coverage. Fitting of the visibilities with a model in the uv-plane yielded similar results.

3. Discussion

3.1. Radio SED

While the 1100 GHz radio spectrum of SN 1987A is consistent with a single power law having α = −0.74 (Zanardo et al. 2013), Fig. 2 shows this is clearly not the case for SN 1978K. A break such as this is not unusual in core collapse supernovae at early phases, where synchrotron self-absorption and/or free-free absorption results in a flattening of the spectral index at lower frequencies that are still optically-thick. However, after 36 yr the expectation would be that SN 1978K should be optically thin at all these frequencies.

Montes et al. (1997) noted a similar curvature in the 0.89 GHz spectrum of SN 1978K at ~5000 days, which they attributed to an additional time-independent free-free absorption component along the line-of-sight, namely a foreground H ii region. Using their Eq. (4) for the SED, but with model parameters from our fit to the longer time baseline data in Smith et al. (2007) we are able to reproduce the observed low frequency turnover as shown in Fig. 2 using a value for the thermal optical depth parameter K4 that is 2.5 × larger than that determined by Montes et al. (1997). This is nevertheless still consistent with the presence of an intervening H ii region, or a circumstellar ejecta nebula as argued by Chu et al. (1999) on the basis of optical spectra. Observations at even lower frequencies would help to further constrain this fit, but the GaLactic and Extragalactic All-Sky MWA Survey (GLEAM; Wayth et al. 2015) at 70230 MHz currently lacks the necessary resolution to cleanly separate SN 1978K from its host galaxy emission.

thumbnail Fig. 3

VLBI image of SN 1978K at a frequency of 8.4 GHz on 2015 March 29. The map is centered at 03:17:38.62, 66:33:03.4 (J2000). The beam shown at lower-left by the hashed ellipse is 3.9 mas × 2.8 mas. Contours are shown at − 2%, 10%, 20%, 30%,..., 90% of the peak brightness of 1.95 mJy/beam.

3.2. Expansion velocity

Were SN 1978K undergoing free expansion at a mean velocity of 10 000 km s-1, it should span nearly 40 mas and be easily resolvable by now. Our previous LBA observations at 1.4 GHz and 2.3 GHz of SN 1978K in 2003 (Smith et al. 2007) put an upper limit on the expanding remnant diameter of ~10 mas, or 0.22 pc at the distance of NGC 1313. It appears that more than a decade later SN 1978K is still barely resolved even at 8.4 GHz, with no apparent shell or ring-like structure yet visible as seen in SN 1987A (Zanardo et al. 2013) or SN 1993J (Marcaide et al. 2009a). We can, however, halve the upper limit on the diameter to <5 mas (0.11 pc), and consequently the average expansion velocity for SN 1978K since explosion is <1500 km s-1. While much lower than the 800010 000 km s-1 average expansion velocity inferred from VLBI measurements of objects like SN 1979C (Marcaide et al. 2009b) and SN 1986J (Bietenholz et al. 2010) some 25 yr after their explosions, we note the recently measured range of current expansion velocities for SN 1987A of between 1600 and 6300 km s-1 depending on azimuth (Zanardo et al. 2013). This confirms that the radio remnant of SN 1987A has also undergone significant deceleration since the initial explosion at ~35 000 km s-1 (Gaensler et al. 1997).

We can make use of this new VLBI constraint on the current diameter of the remnant of SN 1978K, together with information from the radio light curve fits, to infer how the ejecta has decelerated over time. The blast wave radius is assumed to grow with time as R ~ tm, where m = 1 for no deceleration. In the Weiler et al. (2002) parameterisation of the Chevalier (1982a) model m = −(αβ−3) / 3, where α is the radio spectral index, and β is the rate of decline in flux density during the post-peak, optically-thin phase. The SED fit shown in Fig. 2 has α = −0.85 and β = −1.51, so m = −0.78 which is at the extreme end for decelerations measured this way (Weiler et al. 2002).

thumbnail Fig. 4

Modelled evolution of the remnant diameter of SN 1978K (black line), together with the current (blue line) and past average (red line) expansion velocities, for R ~ t0.78.

Figure 4 shows a model of the evolution in diameter, current expansion velocity, and past average expansion velocity of SN 1978K for R ~ t0.78 such that the current diameter is no more than 0.1 pc. Recent optical spectroscopy of SN 1978K by Kuncarayakti et al. (2016) shows remarkable similarities with SN 1987A, and line-widths indicating a current expansion velocity ~500600 km s-1, within a factor of two of that predicted by this model. Whereas the radio emission at such late times is thought to trace the forward shock, the optical emission lines are associated with the slower-moving reverse shock. Kuncarayakti et al. (2016) inferred an extreme mass-loss rate (0.01 M yr-1) based on the Hα luminosity, leading to a massive circumstellar medium (~1 M). Although there are several caveats that need to be considered when converting Hα luminosities to mass loss rates, the values above could simultaneously account for the significant radio luminosity as well as the dramatic deceleration in the expansion. Nevertheless we note that while Type IIn supernovae are among the most luminous at radio wavelengths, SN 1978K is in fact not unusually luminous for its type (e.g. Romero-Cañizales et al. 2014).

4. Conclusions

New radio observations of the Type IIn SN 1978K at high frequencies and at high spatial resolution have led to some intriguing results:

  • 1.

    SN 1978K turns out to be >400× brighter at frequencies above 30 GHz than SN 1987A, the only other extragalactic supernova to have been detected at these frequencies more than 2 decades after explosion.

  • 2.

    A single power-law does not fit the 1100 GHz SED of SN 1978K. This may be because the intrinsic emission has a broken power law, or due to the presence of a nebular absorber along the line of sight as previously suspected from low-frequency observations. Future metre-wave and sub-millimetre observations should be able to discriminate between these two scenarios.

  • 3.

    The remnant of SN 1978K is an order of magnitude smaller than expected for the case of free expansion, and instead appears to be decelerating as R ~ t0.78. The past average expansion rate is below 1 500 km s-1, and the current expansion rate modelled to be 20% lower still.

These findings lend support to the arguments from contemporary X-ray and optical analyses that the progenitor star of SN 1978K underwent extreme mass-loss in its final centuries before explosion. This material is responsible for drastically decelerating the expansion of the remnant in much the same way as observed in SN 1987A, which shows other similarities with the older SN 1978K despite not being classed as a Type IIn event.


Acknowledgments

S.R. and R.K. acknowledge support from the Royal Society International Exchange scheme (IE140343). S.R. was a Visiting Research Fellow at Queen’s University Belfast where much of this paper was written. R.K. also acknowledges support from STFC via ST/L000709/1. We are grateful to Cormac Reynolds for invaluable assistance with reducing the LBA data. We thank F. Bauer for bringing to our attention the millimetre observations of SN 1996cr. The Australia Telescope Compact Array and Long Baseline Array are part of the Australia Telescope National Facility which is funded by the Commonwealth of Australia for operation as a National Facility managed by CSIRO. The data reported here were obtained under Programs C184 and V157 (P.I. S. Ryder).

References

  1. Bartel, N., & Bietenholz, M. F. 2003, ApJ, 591, 301 [NASA ADS] [CrossRef] [Google Scholar]
  2. Bietenholz, M. F., Bartel, N., & Rupen, M. P. 2010, ApJ, 712, 1057 [NASA ADS] [CrossRef] [Google Scholar]
  3. Chevalier, R. A. 1982a, ApJ, 258, 790 [NASA ADS] [CrossRef] [Google Scholar]
  4. Chevalier, R. A. 1982b, ApJ, 259, 302 [NASA ADS] [CrossRef] [Google Scholar]
  5. Chu, Y.-H., Caulet, A., Montes, M. J., et al. 1999, ApJ, 512, L51 [NASA ADS] [CrossRef] [Google Scholar]
  6. Deller, A. T., Brisken, W. F., Phillips, C. J., et al. 2011, PASP, 123, 275 [NASA ADS] [CrossRef] [Google Scholar]
  7. Gaensler, B. M., Manchester, R. N., Staveley-Smith, L., et al. 1997, ApJ, 479, 845 [NASA ADS] [CrossRef] [Google Scholar]
  8. Gottesman, S. T., Broderick, J. J., Brown, R. L., Balick, B., & Palmer, P. 1972, ApJ, 174, 383 [NASA ADS] [CrossRef] [Google Scholar]
  9. Horesh, A., Stockdale, C., Fox, D. B., et al. 2013, MNRAS, 436, 1258 [NASA ADS] [CrossRef] [Google Scholar]
  10. Jacobs, B. A., Rizzi, L., Tully, R. B., et al. 2009, AJ, 138, 332 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  11. Kotak, R., & Vink, J. S. 2006, A&A, 460, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Kuncarayakti, H., Maeda, K., Anderson, J. P., et al. 2016, MNRAS, 458, 2063 [NASA ADS] [CrossRef] [Google Scholar]
  13. Lakićević, M., Zanardo, G., van Loon, J. T., et al. 2012, A&A, 541, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Marcaide, J. M., Alberdi, A., Elosegui, P., et al. 1994, ApJ, 424, L25 [NASA ADS] [CrossRef] [Google Scholar]
  15. Marcaide, J. M., Martí-Vidal, I., Alberdi, A., et al. 2009a, A&A, 505, 927 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Marcaide, J. M., Martí-Vidal, I., Perez-Torres, M. A., et al. 2009b, A&A, 503, 869 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Meunier, C., Bauer, F. E., Dwarkadas, V. V., et al. 2013, MNRAS, 431, 2453 [NASA ADS] [CrossRef] [Google Scholar]
  18. Montes, M. J., Weiler, K. W., & Panagia, N. 1997, ApJ, 488, 792 [NASA ADS] [CrossRef] [Google Scholar]
  19. Panagia, N. 1999, in New Views of the Magellanic Clouds, eds. Y.-H. Chu, N. Suntzeff, J. Hesser, & D. Bohlender, IAU Symp., 190, 549 [Google Scholar]
  20. Qing, G., Wang, W., Liu, J.-F., & Yoachim, P. 2015, ApJ, 799, 19 [NASA ADS] [CrossRef] [Google Scholar]
  21. Romero-Cañizales, C., Herrero-Illana, R., Pérez-Torres, M. A., et al. 2014, MNRAS, 440, 1067 [NASA ADS] [CrossRef] [Google Scholar]
  22. Ryder, S., Staveley-Smith, L., Dopita, M., et al. 1993, ApJ, 416, 167 [NASA ADS] [CrossRef] [Google Scholar]
  23. Smith, I. A., Ryder, S. D., Böttcher, M., et al. 2007, ApJ, 669, 1130 [NASA ADS] [CrossRef] [Google Scholar]
  24. Sramek, R. A., Panagia, N., & Weiler, K. W. 1984, ApJ, 285, L59 [NASA ADS] [CrossRef] [Google Scholar]
  25. Wayth, R. B., Lenc, E., Bell, M. E., et al. 2015, PASA, 32, e025 [NASA ADS] [CrossRef] [Google Scholar]
  26. Weiler, K. W., Sramek, R. A., & Panagia, N. 1986, Science, 231, 1251 [NASA ADS] [CrossRef] [Google Scholar]
  27. Weiler, K. W., Panagia, N., Montes, M. J., & Sramek, R. A. 2002, ARA&A, 40, 387 [NASA ADS] [CrossRef] [Google Scholar]
  28. Williams, C. L., Panagia, N., Van Dyk, S. D., et al. 2002, ApJ, 581, 396 [NASA ADS] [CrossRef] [Google Scholar]
  29. Wilson, W. E., Ferris, R. H., Axtens, P., et al. 2011, MNRAS, 416, 832 [NASA ADS] [CrossRef] [Google Scholar]
  30. Zanardo, G., Staveley-Smith, L., Ng, C.-Y., et al. 2013, ApJ, 767, 98 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

ATCA observations of SN 1978K.

All Figures

thumbnail Fig. 1

Continuum image of SN 1978K at 94 GHz from the ATCA on 2014 Sep. 29, showing a 10σ detection at the SN location (03h17m38.8s, − 66°3304′′ J2000). The synthesised beam is indicated in the lower-left corner and the apparent source elongation is an artifact of the underlying noise pattern.

In the text
thumbnail Fig. 2

Radio spectral energy distribution for SN 1978K in September 2014. The dashed lines mark a power law slope of − 0.6 passing through the ν< 10 GHz data, and a steeper power law with slope of − 0.9 passing through the ν > 10 GHz data. The red solid line shows an SED fit to the full dataset using the model of Weiler et al. (2002), which includes low frequency absorption due to an H ii region along the line of sight.

In the text
thumbnail Fig. 3

VLBI image of SN 1978K at a frequency of 8.4 GHz on 2015 March 29. The map is centered at 03:17:38.62, 66:33:03.4 (J2000). The beam shown at lower-left by the hashed ellipse is 3.9 mas × 2.8 mas. Contours are shown at − 2%, 10%, 20%, 30%,..., 90% of the peak brightness of 1.95 mJy/beam.

In the text
thumbnail Fig. 4

Modelled evolution of the remnant diameter of SN 1978K (black line), together with the current (blue line) and past average (red line) expansion velocities, for R ~ t0.78.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.