Radio emission of SN1993J: the complete picture. I. Re-analysis of all the available VLBI data

We have performed a complete re-calibration and re-analysis of all the available VLBI observations of supernova SN1993J, following an homogeneous and well-defined methodology. Observations of SN1993J at 69 epochs, spanning 13 years, were performed by two teams, which used different strategies and analysis tools. The results obtained by each group are similar, but their conclusions on the supernova expansion and the shape and evolution of the emitting region differ significantly. From our analysis of the combined set of observations, we have obtained an expansion curve with unprecedented time resolution and coverage. We find that the data from both teams are compatible when analyzed with the same methodology. One expansion index ($m_3 = 0.87 \pm 0.02$) is enough to model the expansion observed at 1.7\,GHz, while two expansion indices ($m_1 = 0.933\pm0.010$ and $m_2 = 0.796\pm0.005$), separated by a break time, $t_{br} = 390\pm30$ days, are needed to model the data, at frequencies higher than 1.7\,GHz, up to day 4000 after explosion. We thus confirm the wavelength dependence of the size of the emitting region reported by one of the groups. We also find that all sizes measured at epochs later than day 4000 after explosion are systematically smaller than our model predictions. We estimate the fractional shell width ($0.31 \pm 0.02$, average of all epochs and frequencies) and the level of opacity to the radio emission by the ejecta. We find evidence of a spectral-index radial gradient in the supernova shell, which is indicative of a frequency-dependent ejecta opacity. Finally, we study the distribution and evolution of the azimuthal anisotropies (hot spots) found around the radio shell during the expansion. These anisotropies have intensities of $\sim 20$% of the mean flux density of the shell, and appear to systematically evolve during the expansion.Comment: 13 pages, 9 figures, accepted for publication in A&

Limit to the radio emission from a putative central compact source in SN1993J

SN1993J in M81 is the best studied young radio-luminous supernova in the Northern Hemisphere. We recently reported results from the analysis of a complete set of VLBI observations of this supernova at 1.7, 2.3, 5.0, and 8.4 GHz, covering a time baseline of more than one decade. Those reported results were focused on the kinematics of the expanding shock, the particulars of its evolving non-thermal emission, the density profile of the circumstellar medium, and the evolving free-free opacity by the supernova ejecta. In the present paper, we complete our analysis by performing a search for any possible signal from a compact source (i.e., a stellar-mass black hole or a young pulsar nebula) at the center of the expanding shell. We have performed a stacking of all our VLBI images at each frequency, after subtraction of our best-fit shell model at each epoch, and measured the peak intensity in the stacked residual image. Given the large amount of available global VLBI observations, the stacking of all the residual images allows us to put upper limits to the eventual emission of a putative compact central source at the level of $\sim102$ $\mu$Jy at 5 GHz (or, more conservatively, $\sim192$ $\mu$Jy, if we make a further correction for the ejecta opacity) and somewhat larger at other wavelengths.Comment: 4 pages, 3 figures. Accepted for publication in A&

Jet precession in the active nucleus of M81. Ongoing VLBI monitoring

In a recent publication, we reported results of a multi-frequency VLBI campaign of observations of the Active Galactic Nucleus (AGN) in galaxy M\,81, phase-referenced to the supernova SN\,1993J. We were able to extract precise information on the relative astrometry of the AGN radio emission at different epochs and frequencies. We found strong evidence of precession in the AGN jet (i.e., a systematic evolution in the jet inclination at each frequency) coupled to changes in the overall flux density at the different frequencies. In these proceedings, we summarise the main contents of our previous publication and we report on (preliminary) new results from our follow-up VLBI observations, now phase-referenced to the young supernova SN2008iz. We also briefly discuss how these results match the picture of our previously-reported precession model.Comment: 6 pages, 3 figures. Proceedings of EVN meeting 201

Differential Astrometry over 15 degrees

We observed the pair of radio sources 1150+812 and 1803+784 in November 1993 with a VLBI array, simultaneously recording at 8.4 and 2.3 GHz. We determined the angular separation between the two sources with submilliarcsecond accuracy by using differential techniques. This result demonstrates the feasibility of high precision differential astrometry for radio sources separated in the sky by almost 15 degrees, and opens the avenue to its application to larger samples of radio sources.Comment: 6 pages, latex2e, 2 figures, To appear in the proceedings of the EVN/JIVE Symposium No. 4, New Astronomy Reviews (eds. Garret, M. Campbell, R., and Gurvits, L.

Highlights of the expansion of SN1993J

SN 1993J in M 81, a powerful radio supernova high in the northern sky, has been carefully monitored with VLBI throughout its lifetime. Its extremely circular shell-like radio structure has expanded over 15 years in a rather self-similar way in reasonable agreement with Chevalier’s model. An extension of Chevalier’s model simultaneously accounts for all the light-curve and VLBI results. Beyond the first year a single deceleration parameter, mtrue = 0.87± 0.02, characterizes the shock expansion. However, at short wavelengths and beyond 5 years, an enhanced deceleration, mshortwavelengths = 0.79±0.01, is measured. This enhancement is interpreted as due to a combination of effects ranging from varying free-free opacity in the supernova ejecta at the short wavelengths to the radial shape of the intensity of the magnetic fields in the emitting region. The final radio-loud stage is characterized by an abrupt decrease of radio emission. This large flux-density decay rate can be explained as due to the supernova shock surpassing the outer boundary of the circumstellar medium generated by the supernova progenitor. Presently, the supernova expansion beyond the circumstellar/interstellar boundary is rather radio silent, although there are indications that the supernova remnant could be studied by LOFAR and eMERLIN

1.6 GHz VLBI Observations of SN 1979C: almost-free expansion

We report on 1.6 GHz Very-Long-Baseline-Interferometry (VLBI) observations of supernova SN 1979C made on 18 November 2002. We derive a model-dependent supernova size. We also present a reanalysis of VLBI observations made by us on June 1999 and by other authors on February 2005. We conclude that, contrary to our earlier claim of strong deceleration in the expansion, SN 1979C has been undergoing almost-free expansion ($m = 0.91\pm0.09$; $R \propto t^m$) for over 25 years.Comment: 4 pages, 4 figures; submitted to A&A on 14 May 2009. Accepted on 7 Jul 200

Radio astrometry with chromatic AGN core positions

Aims: The effect of frequency-dependent AGN core positions (``core-shifts'') on radio Very Long Baseline Interferometry (VLBI) global astrometry measurements is investigated. Methods: The basic equations relating to VLBI astrometry are reviewed, including the effects of source structure. A power-law representation of core-shifts, based on both observations and theoretical considerations of jet conditions, is incorporated. Results: It is shown that, in the presence of core-shifts, phase and group-delay astrometry measurements yield different positions. For a core displacement from the jet base parametrized by Delta x (lambda) = k lambda^beta group delays measure a ``reduced'' core-shift of (1-beta) Delta x (lambda). For the astrophysically-significant case of beta = 1, group delays measure no shift at all, giving the position of the jet base. At 8.4 GHz an estimated typical offset between phase and group-delay positions of ~170 uas is smaller than the current ~250 uas precision of group-delay positions of the sources used to define the ICRF; however, this effect must be taken into account for future measurements planned with improved accuracy when comparing with optical positions of AGN to be obtained with the GAIA mission.Comment: Accepted for publication in Astronomy & Astrophysics, 4 page