Astronomy. Astrophysics. Radio evolution of supernova SN 2008iz in M 82

Transcription

1 DOI: 1051/ / c ESO 2016 Astronomy & Astrophysics Radio evolution of supernova SN 2008iz in M 82 N. Kimani 1, K. Sendlinger 1, 2, A. Brunthaler 1, K. M. Menten 1, I. Martí-Vidal 3, C. Henkel 1, 4, H. Falcke 5, 6, T. W. B. Muxlow 7, R. J. Beswick 7, and G. C. Bower 8 1 Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, Bonn, Germany nkimani@mpifr-bonn.mpg.de 2 Argelander-Institut für Astronomie, Auf dem Hügel 71, Bonn, Germany 3 Onsala Space Observatory, Chalmers Univ. of Technology, Onsala, Sweden 4 Astron. Dept., King Abdulaziz University, PO Box 80203, Jeddah Saudi Arabia 5 Department of Astrophysics, Radboud University Nijmegen, Postbus 9010, 6500 GL Nijmegen, The Netherlands 6 ASTRON, Postbus 2, 7990 AA Dwingeloo, The Netherlands 7 Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, The university of Manchester, Oxford Road, Manchester M13 9PL, UK 8 Academia Sinica Institute of Astronomy and Astrophysics, 645 N. Aóhoku Place, Hilo, HI 96720, USA Received 26 April 2016 / Accepted 20 June 2016 ABSTRACT We report on multi-frequency Very Large Array (VLA) and Very Long Baseline Interferometry (VLBI) radio observations for a monitoring campaign of supernova SN 2008iz in the nearby irregular galaxy M 82. We fit two models to the data, a simple time power-law, S t β, and a simplified Weiler model, yielding decline indices of β = 2 ± 7 (days ) and 1.41 ± 2 (days ), respectively. The late-time radio light-curve evolution shows flux-density flares at 970 and 1400 days that are a factor of 2 and 4 higher than the expected flux, respectively. The later flare, except for being brighter, does not show signs of decline at least from results examined so far (2014 January 23; day 2167). We derive the spectral index, α, S ν α for frequencies 1.4 to 43 GHz for SN 2008iz during the period from 430 to 2167 days after the supernova explosion. The value of α shows no signs of evolution and remains steep 1 throughout the period, unlike that of SN 1993J, which started flattening at day 970. From the 4.8 and 8.4 GHz VLBI images, the supernova expansion is seen to start with a shell-like structure that becomes increasingly more asymmetric, then breaks up in the later epochs, with bright structures dominating the southern part of the ring. This structural evolution differs significantly from SN 1993J, which remains circularly symmetric over 4000 days after the explosion. The VLBI 4.8 and 8.4 GHz images are used to derive a deceleration index, m, for SN 2008iz, of 6 ± 2, and the average expansion velocity between days 73 and 1400 as (12.1 ± ) 10 3 km s 1. From the energy equipartition between magnetic field and particles, we estimate the minimum total energy in relativistic particles and the magnetic fields during the supernova expansion and also find the magnetic field amplification factor for SN 2008iz to be in the range of Key words. radio continuum: galaxies galaxies: individual: M 82 supernovae: individual: SN 2008iz 1. Introduction Radio-loud supernovae are rare events of which just a few dozen have been detected so far (Weiler et al. 2002). Most of them are relatively distant or fairly weak, making them difficult to study in great detail. To date, the best-known example is SN 1993J in M 81 (Martí-Vidal et al. 2011a,b and references therein), which has been studied extensively because of its proximity (3.63 Mpc, Freedman et al. 1994), environment (which allows for multiwavelength studies), and galaxy orientation (M 81 is seen almost face-on). The discovery of SN 2008iz (Brunthaler et al. 2009a) offers the possibility of studying another supernova at a very similar distance in great detail and to make a comparison to SN 1993J. For instance, considering the peak 5 GHz radio luminosities, SN 2008iz at L 5 GHz erg 1 Hz 1 is comparable to SN 1993J at L 5 GHz erg 1 Hz 1 (van Dyk et al. 1994). The two supernovae are also comparable in the time they took to reach their peak after the explosion, The VLBI images (FITS files) are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr ( ) or via with SN 2008iz taking 125 days, while SN 1993J took 133 days (Weiler et al. 2007). Supernova 2008iz was discovered in M 82, which is a nearby irregular (IO) galaxy forming part of the M 81 group. This galaxy harbours numerous bright radio supernova remnants with the current estimate at 50 supernovae (Muxlow et al. 1994; Beswick et al. 2006; Fenech et al. 2008; Kronberg et al. 2000; Weiler et al. 2002). Supernova 2008iz was discovered in 2009 April as a bright radio transient (Brunthaler et al. 2009a,b). The explosion date of the supernova is estimated to be 2008 February 18 ± 6 days (Brunthaler et al. 2009b; Marchili et al. 2010). The discovery was made at radio wavelengths with the Very Large Array (VLA) at 22 GHz (Brunthaler et al. 2009a) and confirmed with MERLIN (Muxlow et al. 2009; Beswick et al. 2009) and the Urumqi 25-m telescope at 5 GHz (Marchili et al. 2010). Endeavours to make detections in other astronomical windows have never been successful. For instance, there are no detections in visible light and the X-ray regimes (Brunthaler et al. 2009b). Varenius et al. (2015) also reported the Low-Frequency Array for Radio Astronomy (LOFAR) non-detection at 154 MHz. Article published by EDP Sciences A18, page 1 of 10

2 However, a detection of the near-ir counterpart for SN 2008iz at 480 days was reported by Mattila et al. (2013) after careful analyses using image-subtraction techniques and considering possible uncertainties caused by the differing filter response functions between the Gemini-North telescope and the Hubble Space Telescope. Overall, this makes it challenging to classify this supernova. However, since no type Ia supernova has yet been detected in the radio regime in M 82 (Pérez-Torres et al. 2014), we can firmly say that SN 2008iz is a core-collapse supernova, in agreement with both Brunthaler et al. (2009b) and Mattila et al. (2013). The position of SN 2008iz is estimated to be α J2000 = 09 h 55 m 51 ṣ 551 ± 0 ṣ 008, δ J2000 = ± (Brunthaler et al. 2009b), which is 2.5 (43 pc) southwest of the photometric centre of M 82 based on the 2.2 µm peak (Weiß et al. 2001). The location of the supernova is heavily obscured by dust and gas, hiding the region from direct observations at optical wavelengths (Weiß et al. 2001). Brunthaler et al. (2010) estimated the extinction towards SN 2008iz to be A υ 24.4 mag, concluding that the reason for the supernova non-detections was that it exploded behind a large dusty interstellar cloud. However, using the same data and relations, Mattila et al. (2013) obtained a total extinction of 48.9 mag, suggesting that the authors may have omitted multiplying the value of N(H 2 ) by 2 to convert to hydrogen nuclei before applying the extinction relation of Güver & Özel (2009). The near-ir detection could indicate that the supernova was located in the foreground with most of the H 2 column density behind the site of SN 2008iz (Mattila et al. 2013). Very Long Baseline Interferometry (VLBI) observations reveal the expansion of SN 2008iz to be self-similar with an expansion velocity of km s 1 over the first 430 days after the explosion, making it one of the fastest expanding radio supernovae, with an expansion index, m, of 9 ± 3 (Brunthaler et al. 2010). The modelled light curve by Marchili et al. (2010) showed that synchrotron self-absorption (SSA) is negligible, making free-free absorption (FFA) the most significant process during the supernova expansion. Radio monitoring of SN 2008iz with VLA and VLBI has been ongoing on a regular basis since its discovery. We describe the monitoring and data reduction process in Sect. 2. In Sect. 3 we present multi-frequency radio observations from 36 to 2167 days after the explosion. In Sect. 4 we discuss the variations in the spectral index α, equipartition minimum energy and magnetic fields, and the evolution of the magnetic field evolution B eq and its amplification. We present our summary in Sect Observations and data calibration The radio monitoring campaign conducted with the VLA at 1.4, 4.8, 8.4, 22, and 43 GHz to trace the evolution of the radio emission from SN 2008iz has been ongoing since the discovery of the emission in The first series of observations is archival and corresponds to the data published by (Brunthaler et al. 2009b, 2010). It consists of six observations (epochs 1 6 in Table A.1) obtained between 2008 March 24 (day 36) and 2009 September 19 (day 580). For details on the observations of these epochs, we refer to the respective publications. For the later observations, the data were taken in the standard continuum observation mode with a total bandwidth of 128 MHz, each in dual circular polarisation. During these epochs, observations were made in VLA configurations A, AB, B, and C. Flux density measurements were derived using calibrator 3C 48. The flux calibrator was observed for a total time of two minutes in each observation. The observation used a switching cycle of six minutes, spending on average one minute on the phase calibrator J and five minutes on M 82. The cycles were repeated five times over the observations, yielding an integration time of 25 min on M 82 at each frequency. The last epoch of VLA observations on 2014 January 23 was obtained during a confirmation search of a newly discovered Type Ia supernova, SN 2014j at C and K bands under the observation code TOBS0008 at configuration AB. Data were taken in the standard continuum observation mode with a total bandwidth of 128 MHz, each in dual circular polarisation. The flux calibrator 3C 48 was observed for a total time of three minutes. The observation used a switching cycle of nine minutes, spending an average time of 30 s on the phase calibrator J and eight minutes on M 82. The cycle was repeated four times over the observation, yielding an integration time of 32 min on M 82 at both frequencies. The VLA observations in this work are matched by highresolution VLBI observations at 1.6, 4.8, and 8.4 GHz. These were taken with a network consisting of the Very Long Baseline Array (VLBA), the Green Bank 100-m telescope, and the Effelsberg 100-meter telescope, leading to a maximum resolution of 4.3 milliarcsec (mas) at 1.6 GHz, mas at 4.8 GHz, and mas at 8.4 GHz. SN 2008iz was observed 15 times between 2009 October and 2013 January in dual polarisation. The total bandwidth of most datasets is 64 MHz, although two datasets have a bandwidth of 32 MHz and three have a bandwidth 260 MHz wide. The fluxes were calibrated using system temperature measurements and standard gain curves, and the AGN in M 81 (M 81*) was used as a phase reference calibrator. The calibrator J was observed several times during each observation run, each time for 90 to 120 s. SN 2008iz and M 81* were observed in turns for 45 to 55 s at 1.6 GHz and 80 s at 4.8, and 8.4 GHz during each observation run. The VLA and VLBI data were reduced using standard packages within the AIPS (Astronomical Image Processing Software) of the National Radio Astronomy Observatory (NRAO). To ascertain the late-time enhancement of flux values, which are much higher than expected from the supernova radio light curve, the VLA data was run using CASA-EVLA_pipeline1.1.3 of NRAO for all epochs after 2009 September 19. The results from both data reduction procedures agree. 3. Results Tables A.1 and A.2 contain the complete VLA and VLBI logs of radio flux measurements of SN 2008iz. The first column in each table lists the observation dates. The second column lists the number of days elapsed since explosion, which is estimated to be 2008 February 18 (±6) (Marchili et al. 2010), while the remaining columns list the integrated or peak flux density or the source size at different frequencies. To minimise the contribution of the diffuse emission in M 82 to the supernova measurements, we used only baselines longer than 30 kλ in the imaging. The AIPS task JMFIT was used to determine the flux by fitting a 2D Gaussian to the identified compact supernova source. The flux errors of the VLA observations were derived by adding in quadrature the formal errors from the 2D Gaussian, a 5% systematic error, and a small random error derived as the difference between peak and integrated flux densities. The errors of the VLBI 1.6 GHz fluxes were derived in a similar way; adding in quadrature the formal errors from 2D Gaussian fit to the sources, and a 9% systematic error. However, the supernova is resolved at 4.8 and 8.4 GHZ in the VLBI data, thus their integrated fluxes A18, page 2 of 10

3 N. Kimani et al.: Radio evolution of supernova SN 2008iz in M 82 cannot be derived in the same way as for the other observations. In Table A.2 we list the peak intensities at 4.8 and 8.4 GHz Light curve The light curve of SN 2008iz is well fit by a power law of the form ( ) β t t0 S = K 0, (1) 1 day especially in the well-sampled optically thin regime, that is, 5 and 22 GHz (see Fig. 1). K 0 is a flux scaling factor that does not change the overall behaviour of the fit, (t t 0 ) is the time elapsed after the explosion date, while β is the flux decay index. The obtained fit results (see Fig. 1) with β and K 0 as free parameters at 5 GHz ( days) are K 0 5 GHz = (6.29 ± 2.22) 10 4 mjy and β = 2 ± 7. The time range days was excluded because the radio emission from the supernova shock wave was still within the optically thick regime, while the last epochs of the data sets have their peak intensity higher than expected from the light curve. The obtained fit value at 22 GHz ( days) with both β and K 0 as free parameters yields K 0 22 GHz = (1.46 ± 3) 10 4 mjy and β = 1.18 ± 5. Considering the 5 GHz data set, which is better sampled over the time range than the one at 22 GHz, we adopted a β value of 2 ± 7 to be more representative. Fitting the 22 GHz data with our adopted β value, we obtain the scaling parameter K 0 22 GHz to be (1.73 ± 0) 10 4 mjy. A fit to all frequencies (Fig. 3) using one single β value of 2 ± 7 yields a different K 0 value for each frequency. The values of K 0 are K GHz = (1.19 ± 0) 10 5 mjy, K 0 5 GHz = (6.19 ± 7) 10 4 mjy, K GHz = (3.57 ± 0) 10 4 mjy, K 0 22 GHz = (1.73 ± ) 10 4 mjy, and K 0 43 GHz = (3.77 ± 2) 10 3 mjy. The derived scaling factor K 0 yields a normalised average spectrum that depends on frequency, as shown in Fig. 2. Although the simple power-law (see Eq. (1)) describes the optically thin regime of the light curve very well, a complete fit to the supernova light-curve model is not possible because we have data for only very few epochs in the optically thick regime. For the complete radio light curve, we fit the simplified Weiler model (Weiler et al. 2007) described as ( ν ) α ( ) β t t0 S (ν) = K 1 e τ, (2) 5 GHz 1 day where K 1 is a scaling factor in Jy and α is the spectral index of the emission in the optically thin regime. The exponent β is related to the flux density drop at late epochs, normally affected by both the spectral index, α, and the opacity, τ. The date of explosion is denoted as t 0, and t is the epoch of observation. The opacity of thermal electrons of the circumstellar medium (CSM) in Eq. (2) is modelled as ( ν ) 2.1 ( ) δ t t0 τ = K 2, (3) 5 GHz 1 day where K 2 is the scaling factor and δ is the absorption decline index related to the CSM radial density profile derived as δ = α β 3. The exponent 2.1 corresponds to the spectral dependence of FFA by thermal ionised gas in the radio regime. To constrain the simplified Weiler model, we applied the best-fit parameters of models 2 and 3 obtained by (Marchili et al. 2010) at 5 GHz (i.e. K 1 = (2.14 ± 4) 10 5 mjy, β = 1.43 ± 5, t 0 = 18 Feb. 2008, K 2 = (11.0 ± 0.7) 10 4 mjy, Flux density (mjy) MERLIN 5GHz Urumqi 5GHz 4.8GHz 22.3GHz Fig. 1. Fits to the light curve of the optically thin regime of SN 2008iz at 5 and 22 GHz. β and K 0 (see Eq. (1)) are derived from the fit to the data over the range of days for 5 GHz and days for 22 GHz. The fitted data are shown as filled symbols, while the data not used to fit are shown as open symbols. The MERLIN 5 GHz data were taken from Beswick et al. (2009) and the Urumqi 5 GHz from Marchili et al. (2010). The remaining data are tabulated in Table A.1. K 0 [mjy] Frequency [GHz] Fig. 2. Normalised average spectrum of the variation of K O with frequency for the optically thin power-law model. δ = 2.65±0), and from (Brunthaler et al. 2010) we obtained the spectral index, α = 1.08 ± 8. We fit all frequencies at the same time, as shown in Fig. 4, using β as a free parameter, for which a value of β = 1.41 ± 2 was determined. The model satisfactorily fits the flux density decline from day 500 except for the 1.4 GHz observations, which fall well below the model fit in Fig. 4. The lower flux densities at 1.4 GHz are also confirmed by a 3σ LOFAR non-detection limit of SN 2008iz at a level of 1mJy/beam at the even lower frequency of 154 MHz (Varenius et al. 2015). Plausible explanations might be FFA from a dense foreground screen along the line of sight. Other effects, such as a low-frequency cut-off caused by a Razin- Tsytovich effect, might also help explain the lower 1.4 GHz flux and LOFAR non-detection. From the same figure, the observations at 43.2 GHz drop quite fast past 1000 days. This drop is most likely related to synchrotron ageing of the emitting electrons in the shocked CSM. A18, page 3 of 10

4 GHz 8 GHz 22GHz Flux density (mjy) GHz VLBI 1.6GHz MERLIN 5GHz Urumqi 5GHz 4.8GHz 8.4GHz 22.3GHz 35GHz 43.2GHz Fig. 3. Combined multi-frequency light curve for the optically thin regime of SN 2008iz. The fitted data are shown as filled symbols, while the data not used to fit are shown as open symbols. The data sources are indicated in Fig. 1, while the 1.6 GHz VLBI data are tabulated in Table A.2.The lines represent simple power-law fits to the data at the different frequencies as colour coded in the legend. Flux density [mjy] GHz VLBI 1.6GHz MERLIN 5GHz Urumqi 5GHz 4.8GHz 8.4GHz 22.3GHz 35GHz 43.2GHz Fig. 4. Multi-frequency light curve of SN 2008iz. The MERLIN 5 GHz data were taken from Beswick et al. (2009) and the Urumqi 5 GHz from Marchili et al. (2010). The remaining data are tabulated in Tables A.1 and A.2. The black line represents a Weiler model fit to 5 GHz data obtained by Marchili et al. (2010) with the Urumqi telescope. The other lines represent light curves calculated using the equations in the text for the other frequency data sets as colour coded in the legend. Ringsize R 50% [mas] Fig. 5. Expansion curve of SN 2008iz. The two data points at 22 GHz were taken from Brunthaler et al. (2010), while the last two data points (overlaid with a black ) were excluded from the fit. telescopes, especially the Effelsberg 100 m telescope, did not take part in the observations. The images show that the supernova expansion started with a shell-like structure that became increasingly more asymmetric. Overall, we see expansion with similar velocities in all directions. During the later epochs, the emission distribution breaks up and shows substructures in the southern part of the ring. There, the ring is much brighter than in its other parts, indicating a denser surrounding medium along the southern direction. The deceleration of SN 2008iz is derived from the VLBI images convolved with a dynamic beam at 4.8 and 8.4 GHz. A dynamic beam (see Marcaide et al. 1997) scales the resolution of the later epochs to the resolution of the first epoch with respect to the size of the supernova shell size before imaging. It eliminates the systematic errors that would otherwise arise from the different resolution with which the SNR was imaged while it expanded. We obtained the size of the SNR with the AIPS task IRING, which was used to fit 2D Gaussians to the obtained light profiles. The centre of the SNR remains at the same position during the time span of our observations. We took the radius of the supernova ring to lie at 50% of the peak intensity on the outside of the bright rim since this is the real ring size times an unknown factor. This method is more reliable than taking the maximum of the ring as the radius (see a more detailed description in Brunthaler et al. 2010; Beswick et al. 2006). We fit a power law of the form R 50% = c y (t t 0 ) m (4) 3.2. Deceleration of SN 2008iz The VLBI images of SN 2008iz at 4.8 and 8.4 GHz are shown in Figs. 6 and 7, respectively. All the images are phase referenced to the M 81 core, M 81*, for consistency in their alignment. The SNR is resolved at these frequencies. We used natural weighting to reveal more information on the larger-scale structures. The size and shape of the beam used to restore the images was determined by the set of antennas (uv-coverage) that were employed for each of the considered data sets. For instance, the resolution of the images in Figs. 6e, f, and 7e and h is lower because some to the size evolution of SN 2008iz, where c y is a scaling factor, (t t 0 ) is the time elapsed after the explosion date, and m is the deceleration index. We derive values of c y = (7.5 ± 0.9) 10 6 arcsec/day and m = 6 ± 2. The errors are obtained from the post-fit covariance matrix. The size evolution and deceleration of SN 2008iz are displayed in Fig. 5. We also show two 22 GHz data points from Brunthaler et al. (2010) that were included in the analysis. The last two data points were excluded from the fit since we know from the light curve that there is a change in the surrounding medium of SN 2008iz at around 1400 days after the explosion that is not caused by internal processes. A18, page 4 of 10

5 N. Kimani et al.: Radio evolution of supernova SN 2008iz in M Peak Flux =3.73±4 mjy/b Peak Flux =0.91±3 mjy/b (h) Day Peak Flux =±3 mjy/b 0 0 (j) Day (k) Day 1314 (l) Day Peak Flux =1.65±1 mjy/b Peak Flux =2.43±5 mjy/b 0 5 (i) Day Peak Flux =1.15±3 mjy/b (g) Day Peak Flux =0.76±3 mjy/b (d) Day (f) Day 931 Peak Flux =1.73±3 mjy/b (e) Day 866 Peak Flux =0.98±3 mjy/b (c) Day 720 (b) Day (a) Day 595 Peak Flux =4.56±5 mjy/b Peak Flux =1.79±3 mjy/b Peak Flux =2.09±3 mjy/b 2.4 Peak Flux =2.61±3 mjy/b 2.8 Peak Flux =3.02±1 mjy/b (m) Day 1650 (n) Day 1806 Fig. 6. Natural weighted VLBI images of the expanding shell produced by supernova SN 2008iz at 4.8 GHz. They cover the observation periods between 2009 October and 2013 January as summarized in Table A.2. The figures show the 4.8 GHz linear grey-scale images overlaid with their contour maps that scale to the peak intensity at each epoch (listed in each panel). All image contours correspond to 6% of the peak intensity (lowest contour) and subsequent increases by 6%. The images have been convolved with the natural weighting beam, shown in the bottom left corner of each panel. 4. Discussion 4.1. Flux-density flare from day 970 A zoom into the light curve at the higher frequencies (see Fig. 8) uncovers excess emission from the supernova after t = 970 days and after t = 1400 days. We note that the extra emission after t = 970 days is detected at 8.4, 22.3, and 43.2 GHz, while the enhancement at t = 1400 days is detected at 4.8, 8.4, 22, and 35 GHz. The later flare after time 1400 days is also confirmed with our densely sampled VLBI 1.6 GHz results. The extra emission increases the total flux density of the supernova by a factor of 2 at t = 970 days, lasting for about 240 days. The second flare increases the total flux density by a factor of 4, without showing signs of flux decline so far, that is, in 2014 January. This calls for a continuation of the monitoring to determine how long the second flare lasts. A18, page 5 of 10

6 (i) Day 1125 Peak Flux =7±4 mjy/b (l) Day Peak Flux =9±1 mjy/b Peak Flux =9±3 mjy/b 8 (k) Day (h) Day Peak Flux =0.73±3 mjy/b 6 (j) Day (g) Day 1006 Peak Flux =3±2 mjy/b 4 (f) Day 931 (e) Day Peak Flux =1.02±5 mjy/b 2 6 (d) Day Peak Flux =1±3 mjy/b Peak Flux =6±3 mjy/b 6 8 (c) Day 720 (b) Day (m) Day Peak Flux =5±3 mjy/b Peak Flux =1±3 mjy/b 2 0 (a) Day 595 Peak Flux =1.34±5 mjy/b Peak Flux =8±3 mjy/b 8 Peak Flux =9±3 mjy/b Peak Flux =0.78±3 mjy/b Peak Flux =3±4 mjy/b (n) Day 1650 (o) Day 1806 Fig. 7. Natural weighted VLBI images of the expanding shell produced by supernova SN 2008iz at 8.4 GHz. They cover the observation periods between 2009 October and 2013 January as summarized in Table A.2. The figures show the 8.4 GHz linear grey-scale images overlaid with their contour maps that scale to the peak intensity at each epoch (listed in each panel). All image contours correspond to 6% of the peak intensity (lowest contour) and subsequent increases by 6%. The images have been convolved with the natural weighting beam, shown in the bottom left corner of each panel. The source is obviously optically thin at all frequencies during the flares, and its spectral index remains steep (see Sect. 4.2). Weiler et al. (2002) indicated that a supernova whose radio emission preserves its spectral index while deviating from the standard model shows evidence of a change in the conditions of the CSM that its remnant expands into. We therefore examined the changes in the density of the CSM as the explanation to the observed flares. According to Chevalier (1982), radio luminosity is A18, page 6 of 10 related to the average CSM density (ρcsm M /w) through M L w!(γ 7+12m)/4, (5) where the relativistic particle index γ = 1 2α, α = 1.08 ± 8 and expansion index m = 6 ± 2. For SN 2008iz, L ( M /w)1.62. Consequently, for the radio emission to double its

7 N. Kimani et al.: Radio evolution of supernova SN 2008iz in M Flux density [mjy] Spectral Index GHz 43GHz 35GHz Fig. 8. Zoom into our higher frequency light-curve for a period between t = 800 to t = 2400 days. The flux density enhancement at t 1400 days was detected in both VLA and VLBI results at all frequencies. flux in the first flare, the CSM density must have increased by a factor of 1.5, while for the increase by a factor of 4 in the second flare, the CSM density must have tripled. However, the VLBI images of the remnant of SN 2008iz do not show a spherically symmetric density enhancement. This suggests that the shock wave is encountering dense inhomogeneities in the CSM. The late flare (day 1400) may also be related to the transition of the front shock from the CSM bubble into the ISM. This would enhance the magnetic field and increase the emission at all frequencies. However, we obtain that the size of the shock on day 1400 is 6 pc (3.6 mas), which is very small compared to the supernova remnant sizes measured by Batejat et al. (2011) of pc in Arp 220 at the transition. If the flare were related to the transition, this would indicate a very dense ISM for the pressure balance to hold at the bubble boundary of SN 2008iz. Thus we require an ISM with a density higher than the characteristic one in Arp 220 of 10 4 cm 3 to explain the light curve of SN 2008iz at late stages, in the frame of the ISM interaction model Evolution of the spectral index Figure 9 shows the best-fit radio spectral index, α, for SN 2008iz from 430 days up to 2167 days. To obtain the radio spectral index, we fit a simple power-law spectrum of the form (S = S 0 ν α ) to our VLA data for all the epochs except for epoch 2009 April 27, which was best fit with a broken power-law of the form ( ) α ( ) δ α ν S = S 0 1 e ν ν 0 (6) ν 0 as presented by Brunthaler et al. (2010). We fit this expression to our data and left S 0 and α as free parameters. For epochs with two or three frequency data points only, that is, 2010 October 20, 2010 December 18, 2011 July 30, and 2011 August 01, a systematic error of 30% was added to the fit error to account for the low data statistics. This causes the data points for these epochs to have slightly larger error bars than the epochs with more frequency data points, whose errors are derived directly from the post-fit covariance matrix. The observed spectral index does not show signs of evolution and remains steep, that is, α 1 throughout the period. This steepness seems to persist longer than for SN 1993J, whose spectral index evolution shows flattening at all frequencies beginning at an age Fig. 9. Spectral index, α, for SN 2008iz obtained from fitting a simple power law fit to the radio continuum data at each epoch. The spectral index remains steep, near α = days (Pérez-Torres et al. 2002a; Martí-Vidal et al. 2011b; Weiler et al. 2007) 4.3. Equipartition, total minimum energy, and magnetic field For sources of radio synchrotron emission, energy equipartition between the particles and the magnetic field is usually postulated (Pacholczyk 1970). The total minimum energy content (E min ) and magnetic field (B eq ) in SN 2008iz can be expressed as described in relations 7 and 8 E min = c 13 (1 + k) 4/7 φ 3/7 R 9/7 L 4/7 (7) and B eq = 4.5 2/7 c 2/7 12 (1 + k)2/7 φ 2/7 R 6/7 L 2/7, (8) where the dimensionless parameters c 12 and c 13 are and , respectively (Pacholczyk 1970) for a spectral index α 1. The filling factor, φ, the ratio of the inner and outer radii of the emitting region, is estimated as (Marchili et al. 2010). R is the source radius in cm derived from our VLBI data R 50%, which predicts m = 6 ± 2 and c y = (7.5 ± 0.9) 10 6 arcsec/day in agreement with previous derivations from Brunthaler et al. (2010) of m = 9 ± 3. L is the integrated radio luminosity in erg 1 between the radio synchrotron cut-off frequencies ν min = 10 7 Hz and ν max = Hz. k is the ratio between the relativistic heavy particle energy density to the relativistic electron energy density. The parameter k is dependent on masses of the particles and varies from m p /m e 1 to m p /m e (Marchili et al. 2010). With a varying value of k, we cannot compute a single value of the total minimum energy and equipartition magnetic field, but can only give a range of possibilities (see Table 1). The derived E min values are in the range of erg. Our results of B eq at 2008 May 3, which range between G and 1.78 G, are comparable to the values derived by Marchili et al. (2010) at 5 GHz of G and 2.1 G and also on 2009 April 8 between 4 G and 7 G, compared to 4 G and 1 G, respectively. Figure 10 shows the magnetic field evolution as the supernova ages derived from the integrated 22.3 GHz flux densities from day 36 to 2200 after the explosion. From this figure, we note that the magnetic fields at each epoch scale with time, showing a hint of flattening at the A18, page 7 of 10

8 Table 1. SN 2008iz derived equipartition total minimum energy and magnetic fields. Date Days since Radius Luminosity B eq B eq E min E min (yy/mm/dd) (cm) (erg 1 ) (G) k = 1 (G) k = 2000 (erg) k = 1 (erg) k= /03/ /05/ /04/ /04/ /09/ /10/ /10/ /12/ /02/ /05/ /08/ /02/ /08/ /11/ /01/ Notes. The luminosity is derived from the 22.3 GHz integrated flux values, while the radius is derived from the VLBI results by assuming a constant expansion index m = 6 ± 2 and c y = (7.5 ± 0.9) 10 6 arcsec/day. time day This flattening corresponds to the flux-density flare events when the shock interacts with dense CSM. However, assuming a constant k value, the magnetic field generally decreases according to B t 1 in the full time range. Considering the equipartition between fields and particles, the B eq range between (37 266) mg on 2009 April 8 (day 416) is expected to explain the observed level of radio emission. The choice of the radius, which is cm, allows for a comparison with other supernovae at about the same point of expansion. The energy density of the magnetic field in the remnant is postulated to be lower than the kinetic energy density (i.e. ρ csm υ 2 csm/2 B 2 csm/8π), such that Magnetic field (G) Range of B eq B csm (Ṁυ w ) 1/2 r ( ) 1/2 Ṁ ( υ ) 1/2 ( w r ) 1 = 2.5 mg M yr 1 10 km s cm (9) Using a mass-loss of the progenitor of SN 2008iz of M yr 1 (Marchili et al. 2010), a standard pre-supernova wind velocity υ w of 10 km s 1, and a standard CSM density profile of ρ wind r 2, we obtain B csm = 7 mg, which is a factor of about times smaller than the equipartition field. This indicates that the magnetic field inferred for SN 2008iz cannot originate solely from compression of the existing circumstellar magnetic field, which is predicted to increase the field only by a factor of 4 (Dyson & Williams 1980). Large amplification factors of the magnetic field have also been found for other radio supernova: for SN 2001gd an amplification factor in the range of was determined (Pérez-Torres et al. 2005); the SN 1993J amplification factor of a few hundred (Fransson & Björnsson 1998; Pérez-Torres et al. 2001); for SN 1986J values in the range of (Pérez-Torres et al. 2002b); for SN 1979C in the range of (Pérez-Torres et al. 2005); and between 37 and 260 for SN 2004et (Martí-Vidal et al. 2007). Thus, if equipartition between fields and particles holds, amplification mechanisms other than compression of the circumstellar Fig GHz evolution of the magnetic field in the synchrotronemitting plasma of SN 2008iz with time. The magnetic field evolution is roughly consistent with B t 1 with signs of flattening past day 970. magnetic field need to be invoked to explain the level of radio emission from the supernova, such as turbulent amplification due to Rayleigh-Taylor instability (Chevalier 1982; Chevalier & Blondin 1995). 5. Summary We reported on multi-frequency VLA and VLBI radio data in the optically thin regime for an on-going monitoring campaign of SN 2008iz. We fitted two models to the data, a simple power-law (S t β ) and a simplified Weiler model, yielding β = 2 ± 7 (days ) and 1.41 ± 2 (days ), respectively. The light curve uncovers flux density enhancements at t = 970 days and t = 1400 days by a factor of 2 and 4, respectively. The later flare, except for being brighter, does not show signs of decline, at least from the results A18, page 8 of 10

9 N. Kimani et al.: Radio evolution of supernova SN 2008iz in M 82 examined so far (2014 January 23; day 2167). The flaring activity is attributed to increase in the number density of shocked circumstellar medium (CSM) electrons as the expanding shock wave encounters a clumpy denser medium. The late flare may also be attributed to the transition of the front shock from the CSM bubble into the interstellar medium (ISM), although a very dense ISM (>10 4 cm 3 ) would be required, which is consistent with SN 2008iz being deeply buried within the centre of M 82. The 1.4 GHz flux density values fall below the supernova model fit, a phenomenon confirmed by the 3σ LOFAR nondetection limit at 154 MHz level of 1 mjy/beam. The low flux values could be attributed to free-free absorption (FFA) from a dense foreground screen along the line of sight, or a lowfrequency cut-off caused by the Razin-Tsytovich effect. The spectrum for SN 2008iz for the period from 430 to 2167 days after the supernova explosion shows no signs of evolution and remains steep (α 1). This is different from SN 1993J, whose spectral index evolution shows α flattening of the spectrum at all frequencies, beginning at an age 970 days. From the 4.8 and 8.4 GHz VLBI images, the supernova expansion is seen to start with a shell-like structure, reflecting an expansion with similar velocities into all directions that becomes increasingly more asymmetric at later stages. Finally, in later epochs the structure breaks up, with bright structures dominating the lower part of the ring, which serves as an indication of a denser surrounding medium along the southern direction. This structural evolution differs significantly from SN 1993J, which remains circularly symmetric over 4000 days after the explosion. From the size evolution of the SN 2008iz SNR derived from the 4.8 and 8.4 GHz VLBI images, a deceleration parameter, m, of 6±2 and an expansion velocity of (12.1±) 10 3 km s 1 are derived for the time range between 73 and 1400 days. From the energy equipartition between fields and particles, we estimate the minimum total energy in relativistic particles and the magnetic fields during the supernova expansion. On 2009 April 8, E min erg, which corresponds to an average equipartition magnetic field of B min mg. We derive the average magnetic field in the circumstellar wind of B csm = 7 mg at a radius from the centre of the supernova explosion of cm. Since the supernova shock compression could only enhance the magnetic field by a factor of 4, powerful amplification mechanisms must be at play in SN 2008iz to account for the derived magnetic field amplification of that is responsible for the synchrotron radio emission. Acknowledgements. We wish to thank S. Mattila for the valuable and constructive comments. Naftali Kimani was supported for this research through a stipend from the International Max Planck Research School (IMPRS) for Astronomy and Astrophysics at the Universities of Bonn and Cologne. This work made use of the Swinburne University of Technology software correlator, developed as part of the Australian Major National Research Facilities Programme and operated under licence. More information on the software can for example be found in Deller et al. (2011). References Batejat, F., Conway, J. E., Hurley, R., et al. 2011, ApJ, 740, 95 Beswick, R. J., Riley, J. D., Marti-Vidal, I., et al. 2006, MNRAS, 369, 1221 Beswick, R. J., Muxlow, T. W. B., Pedlar, A., et al. 2009, ATel, 2060 Brunthaler, A., Menten, K. M., Reid, M. J., et al. 2009a, ATel, 2020 Brunthaler, A., Menten, K. M., Reid, M. J., et al. 2009b, A&A, 499, L17 Brunthaler, A., Martí-Vidal, I., Menten, K. M., et al. 2010, A&A, 516, A27 Chevalier, R. A. 1982, ApJ, 258, 790 Chevalier, R., & Blondin, J. M. 1995, ApJ, 444, 312 Deller, A. T., Brisken, W. F., Phillips, C. J., et al. 2011, PASP, 123, 275 Dyson, J., & Williams, D. 1980, Physics of the Interstellar Medium (Manchester: Manchester University Press) Fenech, D. M., Muxlow, T. W. B., Beswick, R. J., Pedlar, A., & Argo, M. K. 2008, MNRAS, 391, 1384 Fransson, C., & Björnsson, C.-I. 1998, ApJ, 509, 861 Freedman, W. L., Hughes, S. M., Madore, B. F., et al. 1994, ApJ, 427, 628 Güver, T., & Özel, F. 2009, MNRAS, 400, 2050 Kronberg, P. P., Sramek, R. A., Birk, G. T., et al. 2000, ApJ, 535, 706 Marcaide, J. M., Alberdi, A., Ros, E., et al. 1997, ApJ, 486, L31 Marchili, N., Martí-Vidal, I., Brunthaler, A., et al. 2010, A&A, 509, A47 Martí-Vidal, I., Marcaide, J. M., Alberdi, A., et al. 2007, A&A, 470, 1071 Martí-Vidal, I., Marcaide, J. M., Alberdi, A., et al. 2011a, A&A, 526, A142 Martí-Vidal, I., Marcaide, J. M., Alberdi, A., et al. 2011b, A&A, 526, A143 Mattila, S., Fraser, M., Smartt, S. J., et al. 2013, MNRAS, 431, 2050 Muxlow, T. W. B., Pedlar, A., Wilkinson, P. N., et al. 1994, MNRAS, 266, 455 Muxlow, T. W. B., Beswick, R. J., Pedlar, A., et al. 2009, ATel, 2073 Pacholczyk, A. G. 1970, Radio astrophysics. Nonthermal processes in galactic and extragalactic sources (San Francisco: Freeman Publishers) Pérez-Torres, M. A., Alberdi, A., & Marcaide, J. M. 2001, A&A, 374, 997 Pérez-Torres, M. A., Alberdi, A., & Marcaide, J. M. 2002a, A&A, 394, 71 Pérez-Torres, M. A., Alberdi, A., Marcaide, J. M., et al. 2002b, MNRAS, 335, L23 Pérez-Torres, M. A., Alberdi, A., Marcaide, J. M., et al. 2005, MNRAS, 360, 1055 Pérez-Torres, M. A., Lundqvist, P., Beswick, R. J., et al. 2014, ApJ, 792, 38 van Dyk, S. D., Weiler, K. W., Sramek, R. A., Rupen, M. P., & Panagia, N. 1994, ApJ, 432, L115 Varenius, E., Conway, J. E., Martí-Vidal, I., et al. 2015, A&A, 574, A114 Weiler, K. W., Panagia, N., Montes, M. J., & Sramek, R. A. 2002, ARA&A, 40, 387 Weiler, K. W., Williams, C. L., Panagia, N., et al. 2007, ApJ, 671, 1959 Weiß, A., Neininger, N., Hüttemeister, S., & Klein, U. 2001, A&A, 365, 571 A18, page 9 of 10

10 Appendix A: Additional tables Table A.1. Log of SN 2008iz integrated radio flux measurements observed with VLA. Date Days since Configuration S 1.4 GHz S 4.8 GHz S 8.4 GHz S 22.3 GHz S 35 GHz S 43.2 GHz (dd/mm/yy) (mjy) (mjy) (mjy) (mjy) (mjy) (mjy) 2008/03/24 36 C 100 ± /05/03 76 C 88.4 ± 2009/04/ B 9.2 ± 2009/04/ B 55 ± 6 30 ± 5 14 ± ± ± 2009/07/ C 21 ± 2 11 ± ± 2009/09/ B 14 ± 2 6 ± 2 5 ± ± 2009/10/ D 3 ± /10/ C 2.4 ± 1.1 ± 2010/12/ C 1 ± ± 1.7 ± 2011/02/ B 1 ± 7.7 ± ± 0.9 ± 2011/05/ AB 23 ± 5 6 ± ± ± ± 2011/07/ A 19 ± ± ± 2011/08/ A 17 ± 2 11 ± ± 2011/08/ A 2.3 ± ± 2012/02/ C 1.9 ± 0.9 ± 2012/08/ B 12 ± 2 9 ± ± 2.0 ± 2012/11/ A 26 ± ± ± 2.7 ± 2014/01/ AB 41 ± ± Notes. The data on days 24 Mar 2008 to 27 Apr 2009 are obtained from Brunthaler et al. (2009b) and Brunthaler et al. (2010). Table A.2. Log of SN 2008iz radio flux and angular size measurements obtained from the VLBI observations. Date Age S 1.6 GHz S 1.6 GHz S 4.8 GHz S 8.4 GHz Θ R50% 4.8 GHz Θ R50% 8.4 GHz (t t 0 ) (mjy) (mjy/b) (mjy/b) (mjy/b) (mas) (mas) 2009/10/ ± ± 3.02 ± 4 3 ± ± 1.7 ± 2009/12/ ± ± 2.61 ± ± ± 1.8 ± 2010/02/ ± ± 2.09 ± 3 9 ± ± 2.0 ± 2010/03/ ± ± 2010/05/ ± 9 3 ± 1.79 ± 3 8 ± ± 2.2 ± 2010/05/ ± ± 2010/06/ ± ± 2010/07/ ± ± 4.56 ± ± ± 2010/09/ ± 9 2 ± 3.73 ± 4 6 ± ± 2.5 ± 2010/11/ ± ± 1.73 ± 3 1 ± ± 2.8 ± 2011/01/ ± ± 1.15 ± ± ± 2.8 ± 2011/03/ ± ± 0.98 ± 3 1 ± ± 3.0 ± 2011/05/ ± ± 0.76 ± 3 3 ± ± 3.0 ± 2011/08/ ± ± 1.55 ± 5 7 ± ± 2011/08/ ± ± 3.78 ± ± /09/ ± ± 0.91 ± 3 9 ± ± 3.3 ± 2011/09/ ± ± 3.17 ± ± ± 2011/12/ ± ± ± 3 5 ± ± 3.7 ± 2012/08/ ± 9 3 ± 2.43 ± ± ± 4.5 ± 2013/01/ ± 9 2 ± 1.65 ± 1 9 ± ± 4.9 ± Notes. Both peak and integrated flux density values of our 1.6 GHz are displayed. The peak flux values are presented for the 4.8 GHz and 8.4 GHz measurements because the supernova SN 2008iz was resolved at these frequencies. The Θ R50% are the angular sizes derived as the radius of the supernova ring at 50% of the peak intensity on the outside of the bright rim. A18, page 10 of 10