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&

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

The Late-Time Radio Spectrum of SN1993J

We present VLA radio continuum measurements of SN1993J in M81 at the frequencies of 0.32 (P-band), 1.3 and 1.7 (L-band), 4.9 (C-band), 8.5 (X-band), and 14.9 (U-band) GHz carried out on December 17 and 21, 2000, about 2820 days after the supernova explosion. We find that a power-law spectrum, free-free absorbed by an homogeneous, or clumpy, distribution of ionized gas yields the best fit to the radio data. A combined homogeneous-clumpy model is not favored, but neither totally excluded. This result contrasts with the modeling of the early (t \lessim 230 days) radio emission from SN1993J, where a mixture of homogeneous and clumpy absorbers appeared to be necessary to adequately describe the behavior of the light curves. The radio spectrum of supernova SN1993J between 0.32 and 14.9 GHz is well characterized by $\alpha = -0.67 \pm 0.02$ ($S_\nu \propto \nu^\alpha$), typical of an optically thin radio supernova. A fit to the radio spectra of SN1993J from $\sim$ 70 up to 2820 days shows that the observed spectral index of SN1993J has been slowly evolving since $t\sim$1000 days on, with the observed spectral index changing from $\alpha \approx -1$ to $\alpha=-0.67$. This spectral evolution seems to suggest that, in addition to the radiative (synchrotron) losses, adiabatic cooling and ionization (Coulomb) losses at the lowest frequencies might be contributing significantly to the integrated electron spectrum.Comment: A&A Main Journal, in press, 6 pages. Includes minor revision

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

Radio emission of SN1993J. The complete picture: II. Simultaneous fit of expansion and radio light curves

We report on a simultaneous modelling of the expansion and radio light curves of SN1993J. We have developed a simulation code capable of generating synthetic expansion and radio light curves of supernovae by taking into consideration the evolution of the expanding shock, magnetic fields, and relativistic electrons, as well as the finite sensitivity of the interferometric arrays used in the observations. Our software successfully fits all the available radio data of SN 1993J with an standard emission model for supernovae extended with some physical considerations, as an evolution in the opacity of the ejecta material, a radial drop of the magnetic fields inside the radiating region, and a changing radial density profile of the circumstellar medium beyond day 3100 after explosion.Comment: 12 pages, 12 figures, accepted for publication in A&

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.

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