We report on the substructure of visibility functions in the delay domain of PSRs B0329+54, B0823+26, B0834+06, B1933+16 and B0833-45 (Vela) observed with earth-earth and RadioAstron space-earth two-element interferometers at frequencies of 324 MHz and 1668 MHz. All visibility functions display unresolved spikes distributed over a range of delays. They are due to band-limited scintillation noise and related to the scattering time. The envelopes for each but the Vela pulsar are well fit by a single Lorentzian which we interpret as being indicative of isotropic scattering on the plane of the sky due to a thin scattering screen between the pulsar and us. In contrast, the envelope for the Vela pulsar needs to be mostly fit by at least two Lorentzians, a narrow and a broad one at the same zero delay. We interpret this characteristic as indicative of anisotropic scattering due to more complex structure of scattering screens in the supernova remnant. The possibility of describing the delay visibility functions by Lorentzians is likely a general property of pulsars and offers a new way of describing scattering parameters of the intervening interstellar medium. Furthermore, for all our pulsars, the unresolved spikes in visibility functions of similar projected baselines were well correlated indicating that the telescopes are located in the same diffraction spot. The correlation vanished for visibilities from largely different baselines, when some radio telescopes are not in the same spot.

We report on VLBI observations of the fast and blue optical transient (FBOT), AT 2018cow. At ~62 Mpc, AT 2018cow is the first relatively nearby FBOT. The nature of AT 2018cow is not clear, although various hypotheses from a tidal disruption event to different kinds of supernovae have been suggested. It had a very fast rise time (3.5 d) and an almost featureless blue spectrum although high photospheric velocities (40,000 km s$^{-1}$) were suggested early on. The X-ray luminosity was very high, ~$1.4 \times 10^{43}$ erg s$^{-1}$, larger than those of ordinary SNe, and more consistent with those of SNe associated with gamma-ray bursts. Variable hard X-ray emission hints at a long-lived "central engine." It was also fairly radio luminous, with a peak 8.4-GHz spectral luminosity of ~$4 \times 10^{28}$ erg s$^{-1}$ Hz$^{-1}$, allowing us to make VLBI observations at ages between 22 and 287 d. We do not resolve AT 2018cow. Assuming a circularly symmetric source, our observations constrain the average apparent expansion velocity to be <0.49$c$ by t = 98 d (3$\sigma$ limit). We also constrain the proper motion of AT 2018cow to be <0.51$c$. Since the radio emission generally traces the fastest ejecta, our observations make the presence of a long-lived relativistic jet with a lifetime of more than one month very unlikely.

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

The largest Key Science Program of the RadioAstron space VLBI mission is a survey of active galactic nuclei (AGN). The main goal of the survey is to measure and study the brightness of AGN cores in order to better understand the physics of their emission while taking interstellar scattering into consideration. In this paper we present detection statistics for observations on ground-space baselines of a complete sample of radio-strong AGN at the wavelengths of 18, 6, and 1.3 cm. Two-thirds of them are indeed detected by RadioAstron and are found to contain extremely compact, tens to hundreds of $\mu$as structures within their cores.

We discuss our VLA and VLBI observations of supernova 1986J, which is characterized by a compact radio-bright component within the expanding shell of ejecta. No other supernova (SN) has such a central component at cm wavelengths. The central component therefore provides a unique probe of the propagation of radio signals at cm wavelengths through the ejecta of a young SN. Such a probe is important in the context of Fast Radio Bursts (FRB), which, in many models, are thought to be produced by young magnetars or neutron stars. The FRB signals would therefore have to propagate through the expanding SN ejecta. Our observations of the Type II SN 1986J show that the ejecta would remain opaque to cm-wave emission like FRBs for at least several decades after the explosion, and by the time the ejecta have become transparent, the contribution of the ejecta to the dispersion measure is likely small.

We report on our efforts to test the Einstein Equivalence Principle by measuring the gravitational redshift with the VLBI spacecraft RadioAstron, in an eccentric orbit around Earth with geocentric distances as small as $\sim$ 7,000 km and up to 350,000 km. The spacecraft and its ground stations are each equipped with stable hydrogen maser frequency standards, and measurements of the redshifted downlink carrier frequencies were obtained at both 8.4 and 15 GHz between 2012 and 2017. Over the course of the $\sim$ 9 d orbit, the gravitational redshift between the spacecraft and the ground stations varies between $6.8 \times 10^{-10}$ and $0.6 \times 10^{-10}$. Since the clock offset between the masers is difficult to estimate independently of the gravitational redshift, only the variation of the gravitational redshift is considered for this analysis. We obtain a preliminary estimate of the fractional deviation of the gravitational redshift from prediction of $\epsilon = -0.016 \pm 0.003_{\rm stat} \pm 0.030_{\rm syst}$ with the systematic uncertainty likely being dominated by unmodelled effects including the error in accounting for the non-relativistic Doppler shift. This result is consistent with zero within the uncertainties. For the first time, the gravitational redshift has been probed over such large distances in the vicinity of Earth. About three orders of magnitude more accurate estimates may be possible with RadioAstron using existing data from dedicated interleaved observations combining uplink and downlink modes of operation.

We present the first extensive radio to gamma-ray observations of a fast-rising blue optical transient (FBOT), AT2018cow, over its first ~100 days. AT2018cow rose over a few days to a peak luminosity $L_{pk}\sim4\times 10^{44}$ erg/s exceeding those of superluminous supernovae (SNe), before declining as $\propto t^{-2}$. Initial spectra at $\lesssim 15$ days were mostly featureless and indicated large expansion velocities v~0.1c and temperatures reaching 30000 K. Later spectra revealed a persistent optically-thick photosphere and the emergence of H and He emission features with v~sim 4000 km/s with no evidence for ejecta cooling. Our broad-band monitoring revealed a hard X-ray spectral component at $E\ge 10$ keV, in addition to luminous and highly variable soft X-rays, with properties unprecedented among astronomical transients. An abrupt change in the X-ray decay rate and variability appears to accompany the change in optical spectral properties. AT2018cow showed bright radio emission consistent with the interaction of a blastwave with $v_{sh}$~0.1c with a dense environment ($\dot M\sim10^{-3}-10^{-4}\,M_{\odot}yr^{-1}$ for $v_w=1000$ km\s). While these properties exclude Ni-powered transients, our multi-wavelength analysis instead indicates that AT2018cow harbored a "central engine", either a compact object (magnetar or black hole) or an embedded internal shock produced by interaction with a compact, dense circumstellar medium. The engine released $\sim10^{50}-10^{51.5}$ erg over $\sim10^3-10^5$ s and resides within low-mass fast-moving material with equatorial-polar density asymmetry ($M_{ej,fast}\lesssim0.3\,\rm{M_{\odot}}$). Successful SNe from low-mass H-rich stars (like electron-capture SNe) or failed explosions from blue supergiants satisfy these constraints. Intermediate-mass black-holes are disfavored by the large environmental density probed by the radio observations.

We discuss the possibility of obtaining Fast Radio Bursts (FRBs) from the interior of supernovae, in particular SN 1986J. Young neutron stars are involved in many of the possible scenarios for the origin of FRBs, and it has been suggested that the high dispersion measures observed in FRBs might be produced by the ionized material in the ejecta of associated supernovae. Using VLA and VLBI measurements of the Type IIn SN 1986J, which has a central compact component not so far seen in other supernovae, we can directly observe for the first time radio signals which originate in the interior of a young (~30 yr old) supernova. We show that at age 30 yr, any FRB signal at ~1 GHz would still be largely absorbed by the ejecta. By the time the ejecta have expanded so that a 1-GHz signal would be visible, the internal dispersion measure due to the SN ejecta would be below the values typically seen for FRBs. The high dispersion measures seen for the FRBs detected so far could of course be due to propagation through the intergalactic medium provided that the FRBs are at distances much larger than that of SN 1986J, which is 10 Mpc. We conclude that if FRBs originate in Type II SNe/SNRs, they would likely not become visible till 60 ~ 200 yr after the SN explosion.

We present an approach to testing the gravitational redshift effect using the RadioAstron satellite. The experiment is based on a modification of the Gravity Probe A scheme of nonrelativistic Doppler compensation and benefits from the highly eccentric orbit and ultra-stable atomic hydrogen maser frequency standard of the RadioAstron satellite. Using the presented techniques we expect to reach an accuracy of the gravitational redshift test of order $10^{-5}$, a magnitude better than that of Gravity Probe A. Data processing is ongoing, our preliminary results agree with the validity of the Einstein Equivalence Principle.

We report on VLA measurements between 1 and 45 GHz of the evolving radio spectral energy distribution (SED) of SN 1986J, made in conjunction with VLBI imaging. The SED of SN 1986J is unique among supernovae, and shows an inversion point and a high-frequency turnover. Both are due to the central component seen in the VLBI images, and both are progressing downward in frequency with time. The optically-thin spectral index of the central component is almost the same as that of the shell. We fit a simple model to the evolving SED consisting of an optically-thin shell and a partly-absorbed central component. The evolution of the SED is consistent with that of a homologously expanding system. Both components are fading, but the shell more rapidly. We conclude that the central component is physically inside the expanding shell, and not a surface hot-spot central only in projection. Our observations are consistent with the central component being due to interaction of the shock with the dense and highly-structured circumstellar medium that resulted from a period of common-envelope evolution of the progenitor. However a young pulsar-wind nebula or emission from an accreting black hole can also not be ruled out at this point.

We present a new 5-GHz global-VLBI image of supernova 1986J, observed in 2014 at $t=31.6$ yr after the explosion, and compare it to previous images to show the evolution of the supernova. Our new image has a dynamic range of ~100 and a background rms noise level of 5.9 $\mu$Jy beam$^{-1}$. There is no significant linear polarization, with the image peak being $<$3% polarized. The latest image is dominated by the compact central component, whose flux density is now comparable to that of the extended supernova shell. This central component is marginally resolved with a FWHM width of $900_{-500}^{+100} \; \mu$as, corresponding to a radius of $r_{\rm comp}=6.7 _{-3.7}^{+0.7} \times 10^{16}$ cm for a distance of 10 Mpc. Using VLBI observations between 2002 and 2014, we measured the proper motions of both the central component and a hot-spot to the NE in the shell relative to the quasar 3C66A. The central component is stationary to within the uncertainty of 12 $\mu$as yr$^{-1}$, corresponding to 570 km s$^{-1}$. Our observations argue in favor of the central component being located near the physical center of SN 1986J. The shell hot-spot had a mean velocity of 2810+-750 km s$^{-1}$ to the NE, which is consistent with it taking part in the homologous expansion of the shell seen earlier. The shell emission is evolving in a non-selfsimilar fashion, with the brightest emission shifting inwards within the structure, and with only relatively faint emission being seen near the outer edge and presumed forward shock. An animation is available.

Observational constraints of relativistic jets from black holes has largely come from the most powerful and extended jets\cite{Jorstad05,Asada14}, leaving the nature of the low luminosity jets a mystery\cite{Falcke04}. M81* is one of the nearest low-luminosity jets, which underwent an extremely large radio flare in 2011, allowing us to study compact core emission with unprecedented sensitivity and linear resolution. Utilizing a multi-wavelength campaign, we were able to track the flare as it re-brightened and became optically thick. Simultaneous X-ray observations indicated the radio re-brightening was preceded by a low energy X-ray flare at least $t_{\rm delay}>12\ {\rm days}$ prior. Associating the time delay between the two bands as the cooling time in a synchrotron flare\cite{Urry97,Bai03}, we find the magnetic field strength was $1.9<B<9.2\ {\rm G}$, which is consistent with magnetic field estimate from spectral-energy distribution modeling\cite{Kellerman81}, $B<10.2\ {\rm G}$. In addition, VLBA observations at 23 GHz clearly illustrate a discrete knot moving mildly relativistically at $v_{\rm app}/c=0.51\pm0.17$ associated with the initial radio flare. The observations indicate radial jet motions for the first time in M81*. This has profound implications for jet production, as it means radial motion can be observed in even the lowest-luminosity AGN, but at slower velocities and smaller radial extents ($\approx10^4\ R_{\rm G}$).

We have resolved the scatter-broadened image of PSR B0329+54 and detected substructure within it. These results are not influenced by any extended structure of a source but instead are directly attributed to the interstellar medium. We obtained these results at 324 MHz with the ground-space interferometer RadioAstron which included the space radio telescope (SRT), ground-based Westerbork Synthesis Radio Telescope and 64-m Kalyazin Radio Telescope on baseline projections up to 330,000 km in 2013 November 22 and 2014 January 1 to 2. At short 15,000 to 35,000 km ground-space baseline projections the visibility amplitude decreases with baseline length providing a direct measurement of the size of the scattering disk of 4.8$\pm$0.8 mas. At longer baselines no visibility detections from the scattering disk would be expected. However, significant detections were obtained with visibility amplitudes of 3 to 5% of the maximum scattered around a mean and approximately constant up to 330,000 km. These visibilities reflect substructure from scattering in the interstellar medium and offer a new probe of ionized interstellar material. The size of the diffraction spot near Earth is 17,000$\pm$3,000 km. With the assumption of turbulent irregularities in the plasma of the interstellar medium, we estimate that the effective scattering screen is located 0.6$\pm$0.1 of the distance from Earth toward the pulsar.

A test of a cornerstone of general relativity, the gravitational redshift effect, is currently being conducted with the RadioAstron spacecraft, which is on a highly eccentric orbit around Earth. Using ground radio telescopes to record the spacecraft signal, synchronized to its ultra-stable on-board H-maser, we can probe the varying flow of time on board with unprecedented accuracy. The observations performed so far, currently being analyzed, have already allowed us to measure the effect with a relative accuracy of $4\times10^{-4}$. We expect to reach $2.5\times10^{-5}$ with additional observations in 2016, an improvement of almost a magnitude over the 40-year old result of the GP-A mission.

We discovered fine-scale structure within the scattering disk of PSR B0329+54 in observations with the RadioAstron ground-space radio interferometer. Here, we describe this phenomenon, characterize it with averages and correlation functions, and interpret it as the result of decorrelation of the impulse-response function of interstellar scattering between the widely-separated antennas. This instrument included the 10-m Space Radio Telescope, the 110-m Green Bank Telescope, the 14x25-m Westerbork Synthesis Radio Telescope, and the 64-m Kalyazin Radio Telescope. The observations were performed at 324 MHz, on baselines of up to 235,000 km in November 2012 and January 2014. In the delay domain, on long baselines the interferometric visibility consists of many discrete spikes within a limited range of delays. On short baselines it consists of a sharp spike surrounded by lower spikes. The average envelope of correlations of the visibility function show two exponential scales, with characteristic delays of $\tau_1=4.1\pm 0.3\ \mu{\rm s}$ and $\tau_2=23\pm 3\ \mu{\rm s}$, indicating the presence of two scales of scattering in the interstellar medium. These two scales are present in the pulse-broadening function. The longer scale contains 0.38 times the scattered power of the shorter one. We suggest that the longer tail arises from highly-scattered paths, possibly from anisotropic scattering or from substructure at large angles.

We report on third epoch VLBI observations of the radio-bright supernova SN 2011dh located in the nearby (7.8 Mpc) galaxy M51. The observations took place at $t=453$ d after the explosion and at a frequency of 8.4 GHz. We obtained a fairly well resolved image of the shell of SN 2011dh, making it one of only six recent supernovae for which resolved images of the ejecta are available. SN 2011dh has a relatively clear shell morphology, being almost circular in outline, although there may be some asymmetry in brightness around the ridge. By fitting a spherical shell model directly to the visibility measurements we determine the angular radius of SN 2011dh's radio emission to be $636 \pm 29$ $\mu$as. At a distance of 7.8 Mpc, this angular radius corresponds to a linear radius of $(7.4 \pm 0.3) \times 10^{16}$ cm and an average expansion velocity since the explosion of $19000^{+2800}_{-2400}$ kms$^{-1}$. We combine our VLBI measurements of SN 2011dh's radius with values determined from the radio spectral energy distribution under the assumption of a synchrotron-self-absorbed spectrum, and find all the radii are consistent with a power-law evolution, with $R \sim t^{0.97\pm0.01}$, implying almost free expansion over the period $t=4$ d to 453 d.

We review the radio very long baseline interferometry (VLBI) observations of the guide star, IM Peg, and three compact extragalactic reference sources, 3C 454.3, B2250+194, and B2252+172, made in support of the NASA/Stanford gyroscope relativity mission, GP-B. The main goal of the observations was the determination of the proper motion of IM Peg relative to the distant universe. VLBI observations made between 1997 and 2005 yield a proper motion of IM Peg of -20.83 $\pm$ 0.09 mas yr$^{-1}$ in RA and -27.27 $\pm$ 0.09 mas yr$^{-1}$ in dec, in a celestial reference frame of extragalactic radio galaxies and quasars virtually identical to the International Celestial Reference Frame 2 (ICRF2). They also yield a parallax for IM Peg of 10.37 $\pm$ 0.07 mas, corresponding to a distance of 96.4 $\pm$ 0.7 pc. The uncertainties are standard errors with statistical and estimated systematic contributions added in quadrature. These results met the pre-launch requirements of the GP-B mission to not discernibly degrade the estimates of the geodetic and frame-dragging effects. The paper also reports on a 1$\sigma$ upper limit on the magnitude of the components of the proper motion of the 'core' of 3C 454.3 relative to the ICRF2 of 46 and 56 $\mu$as yr$^{-1}$ in RA and dec, respectively, and presents densely sampled ellipses of the parallax and the orbit of the giant of the binary system. It further gives a sequence of images of the flickering radio emission relative to the disk of the giant. For a 'movie of a star,' see http://www.yorku.ca/bartel/impeg.mpg .

Gamma Ray Bursts (GRBs) are characterized by ultra-relativistic outflows, while supernovae are generally characterized by non-relativistic ejecta. GRB afterglows decelerate rapidly usually within days, because their low-mass ejecta rapidly sweep up a comparatively larger mass of circumstellar material. However supernovae, with heavy ejecta, can be in nearly free expansion for centuries. Supernovae were thought to have non-relativistic outflows except for few relativistic ones accompanied by GRBs. This clear division was blurred by SN 2009bb, the first supernova with a relativistic outflow without an observed GRB. Yet the ejecta from SN 2009bb was baryon loaded, and in nearly-free expansion for a year, unlike GRBs. We report the first supernova discovered without a GRB, but with rapidly decelerating mildly relativistic ejecta, SN 2012ap. We discovered a bright and rapidly evolving radio counterpart driven by the circumstellar interaction of the relativistic ejecta. However, we did not find any coincident GRB with an isotropic fluence of more than a sixth of the fluence from GRB 980425. This shows for the first time that central engines in type Ic supernovae, even without an observed GRB, can produce both relativistic and rapidly decelerating outflows like GRBs.

A unique test of general relativity is possible with the space radio telescope RadioAstron. The ultra-stable on-board hydrogen maser frequency standard and the highly eccentric orbit make RadioAstron an ideal instrument for probing the gravitational redshift effect. Large gravitational potential variation, occurring on the time scale of $\sim$24 hr, causes large variation of the on-board H-maser clock rate, which can be detected via comparison with frequency standards installed at various ground radio astronomical observatories. The experiment requires specific on-board hardware operating modes and support from ground radio telescopes capable of tracking the spacecraft continuously and equipped with 8.4 or 15 GHz receivers. Our preliminary estimates show that $\sim$30 hr of the space radio telescope's observational time are required to reach $\sim 2\times10^{-5}$ accuracy in the test, which would constitute a factor of 10 improvement over the currently achieved best result.

Accreting black holes are observed to launch relativistic, collimated jets of matter and radiation. In some sources, discrete ejections have been detected with highly relativistic velocities. These particular sources typically have very high mass accretion rates, while sources lower knot velocities are predominantly associated with black holes with relatively low mass accretion rates. We quantify this behavior by examining knot velocity with respect to X-ray luminosity, a proxy for mass accretion rate onto the black hole. We find a positive correlation between the mass-scaled X-ray luminosity and jet knot velocity. In addition, we find evidence that the jet velocity is also a function of polar angle, supporting the "spine-sheath" model of jet production. Our results reveal a fundamental aspect of how accretion shapes mechanical feedback from black holes into their host environments.

RadioAstron space-ground VLBI observations of the pulsar B0950+08, conducted with the 10-m space radio telescope in conjunction with the Arecibo 300-m telescope and Westerbork Synthesis Radio Telescope at a frequency of 324 MHz, were analyzed in order to investigate plasma inhomogeneities in the direction of this nearby pulsar. The observations were conducted at a spacecraft distance of 330,000 km, resulting in a projected baseline of 220,000 km, providing the greatest angular resolution ever achieved at meter wavelengths. Our analysis is based on fundamental behavior of structure and coherence functions. We find that the pulsar shows scintillation on two frequency scales, both much less than the observing frequency; but modulation is less than 100%. We infer that the scattering is weak, but a refracting wedge disperses the scintillation pattern. The refraction angle of this "cosmic prism" is measured as theta_0=1.1 - 4.4 mas, with the refraction direction being approximately perpendicular to the observer velocity. We show that the observed parameters of scintillation effects indicate that two plasma layers lie along the line of sight to the pulsar, at distances of 4.4 - 16.4 pc and 26 - 170 pc, and traveling in different directions relative to the line of sight. Spectra of turbulence for the two layers are found to follow a power law with the indices gamma_1 = gamma_2 = 3.00 +/- 0.08, significantly different from the index expected for a Kolmogorov spectrum of turbulence, gamma=11/3.

We report on a VLA survey for late-time radio emission from 59 supernovae (SNe) of Type I b/c, which have been associated with long-duration gamma-ray bursts (GRBs). An "off-axis" GRB burst (i.e. whose relativistic jet points away from us) is expected to have late-time radio emission even in the absence of significant prompt gamma-ray emission. From our sample, we detected only SN 2003gk with an 8.4-GHz flux density of $2260 \pm 130 \,\mu$Jy. Our subsequent VLBI observations of SN 2003gk, at an age of $\sim$8 yr, allowed us to determine its radius to be $(2.4 \pm 0.4) \times 10^{17}$ cm, or $94 \pm 15$ light days. This radius rules out relativistic expansion as expected for an off-axis GRB jet, and instead suggests an expansion speed of $\sim 10\:000$ km s$^{-1}$ typical for non-relativistic core-collapse supernovae. We attribute the late-onset radio emission to interaction of the ejecta with a dense shell caused by episodic mass-loss from the progenitor. In addition, we present new calculations for the expected radio lightcurves from GRB jets at various angles to the line of sight, and compare these to our observed limits on the flux densities of the remainder of our SN sample. From this comparison we can say that only a fraction of broadlined Type I b/c SNe have a radio-bright jet similar to those seen for GRB afterglows at cosmological distances. However, we also find that for a reasonable range of parameters, as might be representative of the actual population of GRB events rather than the detected bright ones, the radio emission from the GRB jets can be quite faint, and that at present, radio observations do not place strong constraints on off-axis GRB jets.

We present the various science cases for building Band 1 receivers as part of ALMA's ongoing Development Program. We describe the new frequency range for Band 1 of 35-52 GHz, a range chosen to maximize the receiver suite's scientific impact. We first describe two key science drivers: 1) the evolution of grains in protoplanetary disks and debris disks, and 2) molecular gas in galaxies during the era of re-ionization. Studies of these topics with Band 1 receivers will significantly expand ALMA's Level 1 Science Goals. In addition, we describe a host of other exciting continuum and line science cases that require ALMA's high sensitivity and angular resolution. For example, ALMA Band 1 continuum data will probe the Sunyaev-Zel'dovich Effect in galaxy clusters, Very Small Grains and spinning dust, ionized jets from young stars, spatial and flaring studies of Sgr A*, the acceleration sites of solar flares, pulsar wind nebulae, radio supernovae, and X-ray binaries. Furthermore, ALMA Band 1 line data will probe chemical differentiation in cloud cores, complex carbon chain molecules, extragalactic radio recombination lines, masers, magnetic fields through Zeeman effect measurements, molecular outflows from young stars, the co-evolution of star formation and active galactic nuclei, and the molecular content of galaxies at z ~ 3. ALMA provides similar to better sensitivities than the JVLA over 35-50 GHz, with differences increasing with frequency. ALMA's smaller antennas and shorter baselines, greater number of baselines, and single-dish capabilities, however, give it a significant edge for observing extended emission, making wide-field maps (mosaics), or attaining high image fidelity, as required by the described science cases.

We report on sensitive phase-referenced and gated 1.4-GHz VLBI radio observations of the pulsar PSR J0205+6449 in the young pulsar-wind nebula 3C 58, made in 2007 and 2010. We employed a novel technique where the ~105-m Green Bank telescope is used simultaneously to obtain single-dish data used to determine the pulsar's period as well as to obtain the VLBI data, allowing the VLBI correlation to be gated synchronously with the pulse to increase the signal-to-noise. The high timing noise of this young pulsar precludes the determination of the proper motion from the pulsar timing. We derive the position of the pulsar accurate at the milliarcsecond level, which is consistent with a re-determined position from the Chandra X-ray observations. We reject the original tentative optical identification of the pulsar by Shearer and Neustroev (2008), but rather identify a different optical counterpart on their images, with R-band magnitude ~24. We also determine an accurate proper motion for PSR J0205+6449 of (2.3 +- 0.3) mas/yr, corresponding to a projected velocity of only (35 +- 6) km/s for a distance of 3.2 kpc, at p.a. -38 deg. This projected velocity is quite low compared to the velocity dispersion of known pulsars of ~200 km/s. Our measured proper motion does not suggest any particular kinematic age for the pulsar.

We present X-ray, UV/optical, and radio observations of the stripped-envelope, core-collapse supernova (SN) 2011ei, one of the least luminous SNe IIb or Ib observed to date. Our observations begin with a discovery within 1 day of explosion and span several months afterward. Early optical spectra exhibit broad, Type II-like hydrogen Balmer profiles that subside rapidly and are replaced by Type Ib-like He-rich features on the timescale of one week. High-cadence monitoring of this transition suggests that absorption attributable to a high velocity (> 12,000 km/s) H-rich shell is not rare in Type Ib events. Radio observations imply a shock velocity of v = 0.13c and a progenitor star mass-loss rate of 1.4 x 10^{-5} Msun yr^{-1} (assuming wind velocity v_w=10^3 km/s). This is consistent with independent constraints from deep X-ray observations with Swift-XRT and Chandra. Overall, the multi-wavelength properties of SN 2011ei are consistent with the explosion of a lower-mass (3-4 Msun), compact (R* <= 1x10^{11} cm), He core star. The star retained a thin hydrogen envelope at the time of explosion, and was embedded in an inhomogeneous circumstellar wind suggestive of modest episodic mass-loss. We conclude that SN 2011ei's rapid spectral metamorphosis is indicative of time-dependent classifications that bias estimates of explosion rates for Type IIb and Ib objects, and that important information about a progenitor star's evolutionary state and mass-loss immediately prior to SN explosion can be inferred from timely multi-wavelength observations.

We have studied the statistical and polarization properties of giant pulses (GPs) emitted by the millisecond pulsar B1937+21, with high sensitivity and time resolution. The observations were made in June 2005 with the 100-m Robert C. Byrd Green Bank Telescope at S-band (2052-2116 MHz) using the Mk5A VLBI recording system, with formal time resolution of 16 ns. The total observing time was about 4.5 hours; the rate of detection of GPs was about 130 per hour at the average longitudes of the main pulse (MPGPs) and 60 per hour at the interpulse (IPGPs). While the average profile shows well-defined polarization behavior, with regular evolution of the linear polarization position angle (PA), GPs exhibit random properties, occasionally having high linear or circular polarization. Neither MPGPs nor IPGPs show a preferred PA. The cumulative probability distribution (CPD) of GP pulse energy was constructed down to the level where GPs merge with regular pulses and noise. For both MPGPs and IPGPs, the CPD follows a power law with a break, the power index changing from -2.4 at high energy to -1.6 for low energy. Pulse smearing due to scattering masks the intrinsic shape and duration of the detected GPs. The smearing time varied during the observing session within a range of a few hundred nanoseconds. The measured polarization and statistical properties of GPs impose strong constraints on physical models of GPs. Some of these properties support a model in which GPs are generated by the electric discharge caused by magnetic reconnection of field lines connecting the opposite magnetic poles of a neutron star.

We present a set of compelling science cases for the ALMA Band 1 receiver suite. For these cases, we assume in tandem the updated nominal Band 1 frequency range of 35-50 GHz with a likely extension up to 52 GHz; together these frequencies optimize the Band 1 science return. The scope of the science cases ranges from nearby stars to the re-ionization edge of the Universe. Two cases provide additional leverage on the present ALMA Level One Science Goals and are seen as particularly powerful motivations for building the Band 1 Receiver suite: (1) detailing the evolution of grains in protoplanetary disks, as a complement to the gas kinematics, requires continuum observations out to ~35 GHz (~9mm); and (2) detecting CO 3-2 line emission from galaxies like the Milky Way during the epoch of re-ionization, i.e., 6 < z < 10, also requires Band 1 receiver coverage. The range of Band 1 science is wide, however, and includes studies of very small dust grains in the ISM, pulsar wind nebulae, radio supernovae, X-ray binaries, the Galactic Center (i.e., Sgr A*), dense cloud cores, complex carbon-chain molecules, masers, magnetic fields in the dense ISM, jets and outflows from young stars, distant galaxies, and galaxy clusters (i.e., the Sunyaev-Zel'dovich Effect). A comparison of ALMA and the Jansky VLA (JVLA) at the same frequencies of Band 1 finds similar sensitivity performance at 40-50 GHz, with a slight edge for ALMA at higher frequencies (e.g., within a factor of 2 for continuum observations). With its larger number of instantaneous baselines, however, ALMA Band 1data will have greater fidelity than those from the JVLA at similar frequencies.

We used 8.4 GHz VLBI images obtained at up to 35 epochs between 1997 and 2005 to examine the radio structures of the main reference source, 3C 454.3, and two secondary reference sources, B2250+194 and B2252+172, for the guide star for the NASA/Stanford relativity mission Gravity Probe B (GP-B). For one epoch in 2004 May, we also obtained images at 5.0 and 15.4 GHz. The 35 8.4 GHz images for quasar 3C 454.3 confirm a complex, evolving, core-jet structure. We identified at each epoch a component, C1, near the easternmost edge of the core region. Simulations of the core region showed that C1 is located, on average, 0.18 +- 0.06 mas west of the unresolved "core" identified in 43 GHz images. We also identified in 3C 454.3 at 8.4 GHz several additional components which moved away from C1 with proper motions ranging in magnitude between 0.9c and 5c. The detailed motions of the components exhibit two distinct bends in the jet axis located ~3 and ~5.5 mas west of C1. The spectra between 5.0 and 15.4 GHz for the "moving" components are steeper than that for C1. The 8.4 GHz images of B2250+194 and B2252+172, in contrast to those of 3C 454.3, reveal compact structures. The spectrum between 5.0 and 15.4 GHz for B2250+194 is inverted while that for B2252+172 is flat. Based on its position near the easternmost edge of the 8.4 GHz radio structure, close spatial association with the 43 GHz core, and relatively flat spectrum, we believe 3C 454.3 component C1 to be the best choice for the ultimate reference point for the GP-B guide star. The compact structures and inverted to flat spectra of B2250+194 and B2252+172 make these objects valuable secondary reference sources

We describe the NASA/Stanford gyroscope relativity mission, Gravity Probe B (GP-B), and provide an overview of the following series of six astrometric and astrophysical papers that report on our radio observations and analyses made in support of this mission. The main goal of this 8.5 year program of differential VLBI astrometry was to determine the proper motion of the guide star of the GP-B mission, the RS CVn binary IM Pegasi (IM Peg; HR 8703). This proper motion is determined with respect to compact, extragalactic reference sources. The results are: -20.833 +- 0.090 mas/yr and -27.267 +- 0.095 mas/yr for, respectively, the right ascension and declination, in local Cartesian coordinates, of IM Peg's proper motion, and 10.370 +- 0.074 mas (i.e., 96.43 +- 0.69 pc) for its parallax (and distance). Each quoted uncertainty is meant to represent an ~70% confidence interval that includes the estimated contribution from systematic error. These results are accurate enough not to discernibly degrade the GP-B estimates of its gyroscopes' relativistic precessions: the frame-dragging and geodetic effects.

We made VLBI observations at 8.4 GHz between 1997 and 2005 to estimate the coordinates of the "core" component of the superluminal quasar, 3C 454.3, the ultimate reference point in the distant universe for the NASA/Stanford Gyroscope Relativity Mission, Gravity Probe B. These coordinates are determined relative to those of the brightness peaks of two other compact extragalactic sources, B2250+194 and B2252+172, nearby on the sky, and within a celestial reference frame (CRF), defined by a large suite of compact extragalactic radio sources, and nearly identical to the International Celestial Reference Frame 2 (ICRF2). We find that B2250+194 and B2252+172 are stationary relative to each other, and also in the CRF, to within 1-sigma upper limits of 15 and 30 micro-arcsec/yr in RA and decl., respectively. The core of 3C 454.3 appears to jitter in its position along the jet direction over ~0.2 mas, likely due to activity close to the putative supermassive black hole nearby, but on average is stationary in the CRF within 1-sigma upper limits on its proper motion of 39 micro-arcsec/yr (1.0c) and 30 micro-arcsec/yr (0.8c) in RA and decl., respectively, for the period 2002 - 2005. Our corresponding limit over the longer interval, 1998 - 2005, of more importance to GP-B, is 46 and 56 micro-arcsec/yr in RA and decl., respectively. Some of 3C 454.3's jet components show significantly superluminal motion with speeds of up to ~200 micro-arcsec/yr or 5c in the CRF. The core of 3C 454.3 thus provides for Gravity Probe B a sufficiently stable reference in the distant universe.

We present the principal astrometric results of the very-long-baseline interferometry (VLBI) program undertaken in support of the Gravity Probe B (GP-B) relativity mission. VLBI observations of the GP-B guide star, the RS CVn binary IM Pegasi (HR 8703), yielded positions at 35 epochs between 1997 and 2005. We discuss the statistical assumptions behind these results and our methods for estimating the systematic errors. We find the proper motion of IM Peg in an extragalactic reference frame closely related to the International Celestial Reference Frame 2 (ICRF2) to be -20.83 +- 0.03 +- 0.09 mas/yr in right ascension and -27.27 +- 0.03 +- 0.09 mas/yr in declination. For each component the first uncertainty is the statistical standard error and the second is the total standard error (SE) including plausible systematic errors. We also obtain a parallax of 10.37 +- 0.07 mas (distance: 96.4 +- 0.7 pc), for which there is no evidence of any significant contribution of systematic error. Our parameter estimates for the ~25-day-period orbital motion of the stellar radio emission have SEs corresponding to ~0.10 mas on the sky in each coordinate. The total SE of our estimate of IM Peg's proper motion is ~30% smaller than the accuracy goal set by the GP-B project before launch: 0.14 mas/yr for each coordinate of IM Peg's proper motion. Our results ensure that the uncertainty in IM Peg's proper motion makes only a very small contribution to the uncertainty of the GP-B relativity tests.

We present a physical interpretation for the locations of the sources of radio emission in IM Pegasi (IM Peg, HR 8703), the guide star for the NASA/Stanford relativity mission Gravity Probe B. This emission is seen in each of our 35 epochs of 8.4-GHz VLBI observations taken from 1997 to 2005. We found that the mean position of the radio emission is at or near the projected center of the primary to within about 27% of its radius, identifying this active star as the radio emitter. The positions of the radio brightness peaks are scattered across the disk of the primary and slightly beyond, preferentially along an axis with position angle, p.a. = (-38 +- 8) deg, which is closely aligned with the sky projections of the orbit normal (p.a. = -49.5 +- 8.6 deg) and the expected spin axis of the primary. Comparison with simulations suggests that brightness peaks are 3.6 (+0.4,-0.7) times more likely to occur (per unit surface area) near the pole regions of the primary (|latitude| >= 70 deg) than near the equator (|latitude| <= 20 deg), and to also occur close to the surface with ~2/3 of them at altitudes not higher than 25% of the radius of the primary.

We present measurements of the total radio flux density as well as very-long-baseline interferometry (VLBI) images of the star, IM Pegasi, which was used as the guide star for the NASA/Stanford relativity mission Gravity Probe B. We obtained flux densities and images from 35 sessions of observations at 8.4 GHz (wavelength = 3.6 cm) between 1997 January and 2005 July. The observations were accurately phase-referenced to several extragalactic reference sources, and we present the images in a star-centered frame, aligned by the position of the star as derived from our fits to its orbital motion, parallax, and proper motion. Both the flux density and the morphology of IM Peg are variable. For most sessions, the emission region has a single-peaked structure, but 25% of the time, we observed a two-peaked (and on one occasion perhaps a three-peaked) structure. On average, the emission region is elongated by 1.4 +- 0.4 mas (FWHM), with the average direction of elongation being close to that of the sky projection of the orbit normal. The average length of the emission region is approximately equal to the diameter of the primary star. No significant correlation with the orbital phase is found for either the flux density or the direction of elongation, and no preference for any particular longitude on the star is shown by the emission region.

When VLBI observations are used to determine the position or motion of a radio source relative to reference sources nearby on the sky, the astrometric information is usually obtained via: (i) phase-referenced maps; or (ii) parametric model fits to measured fringe phases or multiband delays. In this paper we describe a "merged" analysis technique which combines some of the most important advantages of these other two approaches. In particular, our merged technique combines the superior model-correction capabilities of parametric model fits with the ability of phase-referenced maps to yield astrometric measurements of sources that are too weak to be used in parametric model fits. We compare the results from this merged technique with the results from phase-referenced maps and from parametric model fits in the analysis of astrometric VLBI observations of the radio-bright star IM Pegasi (HR 8703) and the radio source B2252+172 nearby on the sky. In these studies we use central-core components of radio sources 3C 454.3 and B2250+194 as our positional references. We obtain astrometric results for IM Peg with our merged technique even when the source is too weak to be used in parametric model fits, and we find that our merged technique yields superior astrometric results to the phase-referenced mapping technique. We used our merged technique to estimate the proper motion and other astrometric parameters of IM Peg in support of the NASA/Stanford Gravity Probe B mission.

We report on phase-referenced VLBI radio observations of the Type IIb supernova 2011dh, at times t = 83 days and 179 days after the explosion and at frequencies, respectively, of 22.2 and 8.4 GHz. We detected SN 2011dh at both epochs. At the first epoch only an upper limit on SN 2011dh's angular size was obtained, but at the second epoch, we determine the angular radius SN 2011dh's radio emission to be 0.25 +- 0.08 mas by fitting a spherical shell model directly to the visibility measurements. At a distance of 8.4 Mpc this angular radius corresponds to a time-averaged (since t=0) expansion velocity of the forward shock of 21000 +- 7000 km/s. Our measured values of the radius of the emission region are in excellent agreement with those derived from fitting synchrotron self-absorbed models to the radio spectral energy distribution, providing strong confirmation for the latter method of estimating the radius. We find that SN 2011dh's radius evolves in a power-law fashion, with R proportional to t^(0.92 +- 0.10).

We present the most recent VLBI images of SN 1993J, taken at 1.7 GHz on 2010 March 5-6, along with a discussion of its evolution with time. The new image is the latest in a sequence covering almost the entire lifetime of the supernova. For these latest observations we used an "in beam calibrator" technique, and obtained a background rms brightness of 3.7 micro-Jy/beam. The supernova shell remains quite circular in outline. Modulations in brightness are seen around the rim which evolve relatively slowly, having remained generally similar over the last several years of observation. We determine the outer radius of the supernova using visibility-plane model-fitting. The supernova has slowed down to around 30% of its original expansion velocity, and continues to expand with radius approximately proportional to t^0.8, however, deviations from a strict power-law evolution are seen. We do not find any clear-cut evidence for systematically frequency-dependent evolution, suggesting that the radii as determined from visibility-plane model-fitting continue to provide reasonable estimates of the physical outer shock-front radius.

We report on VLBI, as well as VLA radio observations of the Type Ib/c supernova 2009bb. The high radio luminosity of this supernova seems to require relativistic outflow, implying that the early radio emission was "engine-driven", that is driven by collimated outflow from a compact object, even though no gamma-ray emission was seen. The radio light curve shows a general decline, with a "bump" near t = 52 d, seen most prominently at 5 GHz. The lightcurve bump could be either engine-driven, or it might represent the turn-on of the normal radio emission from a supernova, driven by interaction with the CSM rather than by the engine. We undertook VLBI observations to resolve SN 2009bb's relativistic outflow. Our observations constrain the angular outer radius at an age of 85 d to be <0.64 mas, corresponding to <4 x 10^17 cm and an average apparent expansion speed of <1.74c. This result is consistent with the moderately relativistic ejecta speeds implied by the radio luminosity and spectrum.

We report on sensitive new 1.4-GHz VLA radio observations of the pulsar wind nebula G21.5-0.9, powered by PSR J1833-1034, and its environs. Our observations were targeted at searching for the radio counterpart of the shell-like structure seen surrounding the pulsar wind nebula in X-rays. Some such radio emission might be expected as the ejecta from the <~ 1000 yr old supernova expand and interact with the surrounding medium. We find, however, no radio emission from the shell, and can place a conservative 3-sigma upper limit on its 1-GHz surface brightness of 7 x 10^-22 W/m^2/Hz/sr, comparable to the lowest limits obtained for radio emission from shells around other pulsar-wind nebulae. Our widefield radio image also shows the presence of two extended objects of low-surface brightness. We re-examine previous 327-MHz images, on which both the new objects are visible. We identify the first, G21.64-0.84, as a new shell-type supernova remnant, with a diameter of ~13' and an unusual double-shell structure. The second, G21.45-0.59, ~1' in diameter, is likely an HII region.

We report on sensitive 1.4-GHz VLBI radio observations of the unusually luminous supernova remnant SNR 4449-1 in the galaxy NGC 4449, which gave us the first well-resolved image of this object. The remnant's radio morphology consists of two approximately parallel bright ridges, suggesting similarities to the barrel shape seen for many older Galactic supernova remnants or possibly to SN 1987A. The angular extent of the remnant is 65 x 40 mas, corresponding to (3.7 x 2.3) x 10^{18} (D/3.8 Mpc) cm. We also present a new, high signal-to-noise optical spectrum. By comparing the remnant's linear size to the maximum velocities measured from optical lines, as well as using constraints from historical images, we conclude that the supernova explosion occurred between ~1905 and 1961, likely around 1940. The age of the remnant is therefore likely ~70 yr. We find that SNR 4449-1's shock wave is likely still interacting with the circumstellar rather than interstellar medium.

We present new VLBI images of supernova 1986J, taken at 5, 8.4 and 22 GHz between t=22 to 25 yr after the explosion. The shell expands as t^(0.69+-0.03). We estimate the progenitor's mass-loss rate at (4 ~ 10) * 10^-5 Msol/yr (for v_w = 10 km/s). Two bright spots are seen in the images. The first, in the northeast, is now fading. The second, very near the center of the projected shell and unique to SN1986J, is still brightening relative to the shell, and now dominates the VLBI images. It is marginally resolved at 22 GHz (diameter ~0.3 mas; ~5 * 10^16 cm at 10 Mpc). The integrated VLA spectrum of SN1986J shows an inversion point and a high-frequency turnover, both progressing downward in frequency and due to the central bright spot. The optically-thin spectral index of the central bright spot is indistinguishable from that of the shell. The small proper motion of 1500+-1500 km/s of the central bright spot is consistent with our previous interpretation of it as being associated with the expected black-hole or neutron-star remnant. Now, an alternate scenario seems also plausible, where the central bright spot, like the northeast one, results when the shock front impacts on a condensation within the circumstellar medium (CSM). The condensation would have to be so dense as to be opaque at cm wavelengths (~1000x denser than the average corresponding CSM) and fortuitously close to the center of the projected shell. We include a movie of the evolution of SN1986J at 5 GHz from t=0 to 25 yr.

We report on two epochs of very-long-baseline interferometry (VLBI) observations of the Type Ib/c supernova SN 2008D, which was associated with the X-ray outburst XRF 080109. At our first epoch, at t = 30 days after the explosion, we observed at 22 and 8.4 GHz, and at our second, at t = 133 days, at 8.4 and 5.0 GHz. The VLBI observations allow us to accurately measure the source's size and position at each epoch, and thus constrain its expansion velocity and proper motion. We find the source at best marginally resolved at both epochs, allowing us to place a 3sigma upper limit of ~0.75c on the expansion velocity of a circular source. For an elongated source, our measurements are compatible with mildly relativistic expansion. However, our 3sigma upper limit on the proper motion is 4 micro-arcsec/day, corresponding to an apparent velocity of <0.6c, and is consistent with a stationary flux centroid. This limit rules out a relativistic jet such as an gamma-ray burst jet away from the line of sight, which would be expected to show apparent proper motion of >c. Taken together, our measurements argue against the presence of any long-lived relativistic outflow in SN 2008D. On the other hand, our measurements are consistent with the nonrelativistic expansion velocities of <30,000 km/s and small proper motions (<500 km/s) seen in typical supernovae.

[ABRIDGED VERSION] The future of cm and m-wave astronomy lies with the Square Kilometre Array (SKA), a telescope under development by a consortium of 17 countries. The SKA will be 50 times more sensitive than any existing radio facility. A majority of the key science for the SKA will be addressed through large-area imaging of the Universe at frequencies from 300 MHz to a few GHz. The Australian SKA Pathfinder (ASKAP) is aimed squarely in this frequency range, and achieves instantaneous wide-area imaging through the development and deployment of phase-array feed systems on parabolic reflectors. This large field-of-view makes ASKAP an unprecedented synoptic telescope poised to achieve substantial advances in SKA key science. The central core of ASKAP will be located at the Murchison Radio Observatory in inland Western Australia, one of the most radio-quiet locations on the Earth and one of the sites selected by the international community as a potential location for the SKA. Following an introductory description of ASKAP, this document contains 7 chapters describing specific science programmes for ASKAP. The combination of location, technological innovation and scientific program will ensure that ASKAP will be a world-leading radio astronomy facility, closely aligned with the scientific and technical direction of the SKA. A brief summary chapter emphasizes the point, and considers discovery space.

We review VLBI observations of supernovae over the last quarter century and discuss the prospect of imaging future supernovae with space VLBI in the context of VSOP-2. From thousands of discovered supernovae, most of them at cosmological distances, ~50 have been detected at radio wavelengths, most of them in relatively nearby galaxies. All of the radio supernovae are Type II or Ib/c, which originate from the explosion of massive progenitor stars. Of these, 12 were observed with VLBI and four of them, SN 1979C, SN 1986J, SN 1993J, and SN 1987A, could be imaged in detail, the former three with VLBI. In addition, supernovae or young supernova remnants were discovered at radio wavelengths in highly dust-obscured galaxies, such as M82, Arp 299, and Arp 220, and some of them could also be imaged in detail. Four of the supernovae so far observed were sufficiently bright to be detectable with VSOP-2. With VSOP-2 the expansion of supernovae can be monitored and investiated with unsurpassed angular resolution, starting as early as the time of the supernova's transition from its opaque to transparent stage. Such studies can reveal, in a movie, the aftermath of a supernova explosion shortly after shock breakout.

VLBI observations at 5 GHz have revealed that supernova 1979C, in the galaxy M100 in the Virgo cluster, has shell structure. The shell is approximately circular with the 50% contour deviating from a circle by an rms of no more than 8% of the radius. The brightness distribution along the ridge may vary. The position of the center of the shell is at R.A.=12h 22m 58.66758s and decl.=15^o 47' 51.9695" (J2000), with a standard error of 0.8 mas in each of the coordinates. No isolated central component is visible above a flux density limit of ~0.15 mJy which corresponds to an upper limit of ~15 times the corresponding spectral luminosity of the Crab Nebula. The radio lightcurve is clearly falling again and the radio spectrum is now flatter than at earlier times. SN 1979C is only the fourth, and oldest (optically identified) supernova of which a detailed image could be obtained.

We report on new VLT optical spectroscopic and multi-wavelength archival observations of SN1996cr, a previously identified ULX known as Circinus Galaxy X-2. Our optical spectrum confirms SN1996cr as a bona fide type IIn SN, while archival imaging isolates its explosion date to between 1995-02-28 and 1996-03-16. SN1996cr is one of the closest SNe (~3.8 Mpc) in the last several decades and in terms of flux ranks among the brightest radio and X-ray SNe ever detected. The wealth of optical, X-ray, and radio observations that exist for this source provide relatively detailed constraints on its post-explosion expansion and progenitor history, including an preliminary angular size constaint from VLBI. The archival X-ray and radio data imply that the progenitor of SN1996cr evacuated a large cavity just prior to exploding: the blast wave likely expanded for ~1-2 yrs before eventually striking the dense circumstellar material which surrounds SN1996cr. The X-ray and radio emission, which trace the progenitor mass-loss rate, have respectively risen by a factor of ~2 and remained roughly constant over the past 7 yr. This behavior is reminiscent of the late rise of SN1987A, but 1000 times more luminous and much more rapid to onset. Complex Oxygen line emission in the optical spectrum further hints at a possible concentric shell or ring-like structure. The discovery of SN1996cr suggests that a substantial fraction of the closest SNe observed in the last several decades have occurred in wind-blown bubbles. An Interplanetary Network position allows us to reject a tentative GRB association with BATSE 4B960202. [Abridged]

We report on new 5-GHz VLA radio observations of the pulsar-powered supernova remnant G21.5-0.9. These observations have allowed us to make a high-quality radio image of this remnant with a resolution of ~0.7". It has a filamentary structure similar to that seen in the Crab Nebula. Radio structure suggestive of the torus seen around the Crab pulsar is tentatively identified. We also compared the new image with one taken ~15 yr earlier at 1.5 GHz, both to find the expansion speed of the remnant and to make a spectral index image. Between 1991 and 2006, we find that the average expansion rate of the remnant is 0.11 +/- 0.02 %/year, corresponding, for a distance of 5 kpc, to a speed of 910 +/- 160 km/s wrt. the centre of the nebula. Assuming undecelerated expansion, this expansion speed implies that the age of G21.5-0.9 is 870 (+200,-150) yr, which makes PSR J1833-1034 one of the youngest, if not the youngest, known pulsars in the Galaxy.

We compare the angular expansion velocities, determined with VLBI, with the linear expansion velocities measured from optical spectra for supernova 1993J in the galaxy M81, over the period from 7 d to ~9 yr after shock breakout. We estimate the distance to SN 1993J using the Expanding Shock Front Method (ESM). We find the best distance estimate is obtained by fitting the angular velocity of a point halfway between the contact surface and outer shock front to the maximum observed hydrogen gas velocity. We obtain a direct, geometric, distance estimate for M81 of D=3.96+-0.05+-0.29 Mpc with statistical and systematic error contributions, respectively, corresponding to a total standard error of $+-0.29 Mpc. The upper limit of 4.25 Mpc corresponds to the hydrogen gas with the highest observed velocity reaching no farther out than the contact surface a few days after shock breakout. The lower limit of 3.67 Mpc corresponds to this hydrogen gas reaching as far out as the forward shock for the whole period, which would mean that Rayleigh-Taylor fingers have grown to the forward shock already a few days after shock breakout. Our distance estimate is 9+-13 % larger than that of 3.63+-0.34 Mpc from the HST Key Project, which is near our lower limit but within the errors.

We report on our second-epoch VLBI and VLA observations of the Type Ib/c supernova 2001em, five years after the explosion. It was suggested that SN 2001em might be a jet-driven gamma ray burst (GRB), with the jet oriented near the plane of the sky, which would entail relativistic expansion or motion. Our VLBI image shows that SN 2001em is still unresolved five years after the explosion. For a distance of 83 Mpc (H_0 = 70 km/s/Mpc), the nominal expansion velocity is 5800 +/- 10,000 km/s, and the proper motion is 33,000 +/- 34,000 km/s. Our values are inconsistent with either relativistic expansion or motion, but are consistent with the non-relativistic expansion speeds and small proper motions seen in other supernovae. In particular these values are consistent with radio emission from SN 2001em being due to normal, non-relativistic supernova ejecta interacting with the circumstellar medium. Our VLA observations show a power-law decay in flux density since the time of the peak in the 8.4 GHz radio lightcurve in ~2003.

VLBI observations of the Crab pulsar with the 64-m radio telescope at Kalyazin (Russia) and the 46-m radio telescope of the Algonquin Radio Observatory (Canada) at 2.2 GHz and single-dish observations of the millisecond pulsar B1937+21 with the GBT (USA) at 2.1 GHz were conducted to probe the interstellar medium and study the properties of giant pulses. The VLBI data were processed with a dedicated software correlator, which allowed us to obtain the visibility of single giant pulses. Two frequency scales of 50 and 450 kHz were found in the diffraction spectra of giant pulses from the Crab pulsar. The location of the scattering region was estimated to be close to the outer edge of the nebula. No correlation was found between the power spectra of giant pulses at left- and right-hand circular polarization. We explain this lack of correlation through the influence of the strong magnetic field on circularly polarized emission in the region close to the Crab pulsar. Combining the measurement of the decorrelation bandwidth with that of the scattering time of giant pulses for B1937+21, we found three frequency scales of 1.7, 3.8, and 16.5 MHz. The scattering time of giant pulses of B1937+21 at 2.1 GHz was found to be 40+-4 ns. We obtained an upper limit of the intrinsic width of giant pulses from B1937+21 of less than 8 ns. The frequency dependences of the scattering times for the Crab pulsar and PSR B1937+21 were found to be different. They are characterized by exponents of -3.5 and -4.2, respectively. We attribute the difference to the large influence of scattering in the Crab nebula.

The second fastest millisecond pulsar, B1937+21, is one of several pulsars known to emit giant pulses (GPs). GPs are characterized by their huge energy, power-law cumulative energy distribution, and particular longitudes of occurrence. All these characteristics are different from those of regular pulses. Here, we present a study of GPs from our observations of the pulsar B1937+21 with the GBT at 2.1 GHz in both left and right circular polarization with a time resolution of 8 ns. The Mark5 data acquisition system was used for the first time in single-dish observations with the GBT. This allowed us to obtain continuous and uniform recording for 7.5 hours with a data rate of 512 Mbps. As a result, more than 6000 GPs were found above a detection threshold of 200 Jy. We report on instantaneous spectra of GPs, as well as on a comparison with scintillation spectra of regular emission, on the distribution of GP energies, and on polarization properties of GPs.

We report on VLBI observations of the Type I b/c supernova 2001em, three years after the explosion. It had been suggested that SN 2001em might be a jet-driven gamma ray burst (GRB), with the jet oriented far from the line of sight so that the GRB would not be visible from earth. To test this conjecture, we determined the size of SN 2001em. It is only marginally resolved at our resolution of ~0.9 mas. The 3-sigma upper limit on the major axis angular size of the radio source was 0.59 mas (FWHM of an elliptical Gaussian), corresponding to a one-sided apparent expansion velocity of 70,000 km/s for a distance of 80 Mpc. No low-brightness jet was seen in our image to a level of 4% of the peak brightness. If we assume instead a spherical shell geometry typical of a supernova, we find the angular radius of SN 2001em was 0.17 (-0.10,+0.06) mas, implying an isotropic expansion velocity of 20,000 (-12,000, +7,000) km/s, which is comparable to the expansion velocities of supernova shells. Our observations, therefore, are not consistent with a relativistically expanding radio source in SN 2001em, but are instead consistent with a supernova shell origin of SN 2001's radio emission.

We present four new, high-resolution VLA radio images of the Crab nebula, taken between 2001 Feb. and Apr. The radio images show systematic variability in the Crab's radio emission throughout the region near the pulsar. The principal geometry of the variable features is that of elliptical ripples very similar to the optical wisps. The radio wisps are seen to move systematically outward with projected speeds of up to 0.3c. Comparing the new radio images to our earlier ones from 1998 and 2000, we show there are also more slowly moving features somewhat farther away from the pulsar. In particular, there is a prominent moving feature to the northwest of the pulsar which has a projected speed of order 10,000 km/s. Striation is seen throughout the nebula, suggesting the presence of wave-like disturbances propagating through the synchrotron bubble. The radio images were taken simultaneously with HST optical observations. Comparing the radio to the optical images, we find that the radio wisps are sometimes displaced from the optical ones or have no optical counterparts. We also find that some optical wisps in particular, the brightest optical wisps near the pulsar, do not seem to have radio counterparts. In the exterior of the nebula, by contrast, there is generally a good correspondence between the radio and optical features.

We report on VLBI observations of M81*, the northwest-southeast oriented nuclear core-jet source of the spiral galaxy M81, at five different frequencies between 1.7 and 14.8 GHz. By phase referencing to supernova 1993J we can accurately locate the emission region of M81* in the galaxy's reference frame. Although the emission region's size decreases with increasing frequency while the brightness peak moves to the southwest, the emission region seems sharply bounded to the southwest at all frequencies. We argue that the core must be located between the brightness peak at our highest frequency (14.8 GHz) and the sharp bound to the southwest. This narrowly constrains the location of the core, or the purported black hole in the center of the galaxy, to be within a region of +/-0.2 mas or +/-800 AU (at a distance of ~4 Mpc). This range includes the core position that we determined earlier by finding the most stationary point in the brightness distribution of M81* at only a single frequency. This independent constraint therefore strongly confirms our earlier core position. Our observations also confirm that M81* is a core-jet source, with a one-sided jet that extends to the northeast from the core, on average curved somewhat to the east, with a radio spectrum that is flat or inverted near the core and steep at the distant end. The brightness peak is unambiguously identified with the variable jet rather than the core, which indicates limitations in determining the proper motion of nearby galaxies and in refining the extragalactic reference frame.

Giant radio pulses of the millisecond pulsar B1937+21 were recorded with the S2 VLBI system at 1.65 GHz with NASA/JPL's 70-m radio telescope at Tidbinbilla, Australia. These pulses have been observed as strong as 65000 Jy with widths <= 15 ns, corresponding to a brightness temperature T_b >= 5 x 10^39 K, the highest observed in the universe. The vast majority of these pulses occur in a 5.8 mcs and 8.2 mcs window at the very trailing edges of the regular main pulse and interpulse profiles, respectively. Giant pulses occur in general with a single spike. Only in one case out of 309 was the structure clearly more complex. The cumulative distribution is fit by a power law with index -1.40 +/- 0.01 with a low-energy but no high-energy cutoff. We estimate that giant pulses occur frequently but are only rarely detected. When corrected for the directivity factor, 25 giant pulses are estimated to be generated in one neutron star revolution alone. The intensities of the giant pulses of the main pulses and interpulses are not correlated with each other nor with the intensities or energies of the main pulses and interpulses themselves. Their radiation energy density can exceed 300 times the plasma energy density at the surface of the neutron star and can even exceed the magnetic field energy density at that surface. We therefore do not think that the generation of giant pulses is linked to the plasma mechanisms in the magnetosphere. Instead we suggest that it is directly related to discharges in the polar cap region of the pulsar.

We present two sequences of VLBI images of SN1993J in M81, with 24 images at 8.4 GHz and 19 images at 5.0 GHz, These sequences, from 50 d to ~9 yr after shock breakout, show the evolution of the expanding radio shell of an exploded star in detail. The images are all phase-referenced to the stable reference point of the core of M81, allowing us to display them relative to the supernova explosion center. At 50 d, SN1993J is almost unresolved with a radius of 520 AU. The shell structure becomes discernible at 175 d. The brightness of the ridge of the projected shell is not uniform, but rather varies by a factor of two, having a distinct maximum to the south-east and a minimum to the west. Over the next ~350 d, this pattern appears to rotate counter-clockwise. After two years, the structure becomes more complex with hot spots developing in the east, south, and west. The pattern of modulation continues to change, and after five years the three hot spots have shifted somewhat. After nine years, the radio shell has expanded to a radius of 19,000 AU. The brightness in the center of the images is lower than expected for an optically thin, spherical shell. Absorption in the center is favored over a thinner shell in the back and/or front. Allowing for absorption, we find that the thickness of the shell is (25+/-3)% of its outer radius. We find no compact source in the central region and conclude that any pulsar nebula in the center of SN 1993J is either much fainter than the Crab or affected by remaining significant internal radio absorption.

VLBI astrometry allows coordinates of celestial radio sources to be determined with submilliarcsecond accuracy. In particular, with differential VLBI astrometry, the standard errors of relative positions and proper motions can be as small as 10 microarcseconds and 10 microarcseconds per year, respectively. I will review astrometric observations of AGN, M81*, a supernova, pulsars, and radio stars, including the guide star, HR 8703, for the NASA/Stanford gyroscope relativity mission, Gravity Probe B.

We present VLBI images of the RS CVn binary star HR 5110 (=BH CVn; HD 118216), obtained from observations made at 8.4 GHz on 1994 May 29/30 in support of the NASA/Stanford Gravity Probe B project. Our images show an emission region with a core-halo morphology. The core was 0.39 +/- 0.09 mas (FWHM) in size, or 66% +/- 20% of the 0.6 +/- 0.1 mas diameter of the chromospherically active K subgiant star in the binary system. The halo was 1.95 +/- 0.22 mas (FWHM) in size, or 1.8 +/- 0.2 times the 1.1 +/- 0.1 mas separation of the centers of the K and F stars. The core increased significantly in brightness over the course of the observations and seems to have been the site of flare activity that generated an increase in the total flux density of ~200% in 12 hours. The fractional circular polarization simultaneously decreased from ~10% to 2.5%.

We report on a study of microstructure and its quasi-periodicities of three pulsars at 1.65 GHz with the S2 VLBI system at a resolution of 62.5 ns, by far the highest for any such statistical study yet. For PSR B1929+10 we found in the average cross-correlation function (CCF) broad microstructure with a characteristic timescale of 95+-10 mcs and confirmed microstructure with characteristic timescales between 100 and 450 mcs for PSRs B0950+08 and B1133+16. On a finer scale PSRs B0950+08, B1133+16 (component II) and B1929+10 show narrow microstructure with a characteristic timescale in the CCFs of ~10 mcs, the shortest found in the average CCF or autocorrelation function (ACF) for any pulsar, apart perhaps for the Crab pulsar. Histograms of microstructure widths are skewed heavily toward shorter timescales but display a sharp cutoff. The shortest micropulses have widths between 2 and 7 mcs. No nanopulses or unresolved pulse spikes were detected. Cross-power spectra of single pulses show a large range of complexity with single spectral features representing classic quasi-periodicities and broad and overlapping features with essentially no periodicities at all. Significant differences were found for the two components of PSR B1133+16 in every aspect of our statistical analysis of micropulses and their quasi-periodicities. Asymmetries in the magnetosphere and the hollow cone of emission above the polar cap of the neutron star may be responsible for these differences.

We report on VLBI observations of supernova 1986J in the spiral galaxy NGC 891 at two new epochs, 1990 July and 1999 February, t=7.4 and 15.9 yr after the explosion, and on a comprehensive analysis of these and earlier observations from t~4 yr after the explosion date, which we estimate to be 1983.2 +/- 1.1. The source is a shell or composite, and continues to show a complex morphology with large brightness modulations along the ridge and with protrusions. The supernova is moderately to strongly decelerated. The average outer radius expands as t^(0.71 +/- 0.11), and the expansion velocity has slowed to 6000 km/s at t=15.9 yr from an extrapolated 20,000 km/s at t=0.25 yr. The structure changes significantly with time, showing that the evolution is not self-similar. The shell structure is best visible at the latest epoch, when the protrusions have diminished somewhat in prominence and a new, compact component has appeared. The radio spectrum shows a clear inversion above 10 GHz. This might be related to a pulsar nebula becoming visible through the debris of the explosion. The radio flux density between 1.5 and 23 GHz decreases strongly with time, with the flux density proportional to t^(-2.94 +/- 0.24) between t~15 to 19 yr. This decrease is much more rapid than that found in earlier measurements up to t~6 yr.

Phase-referenced VLBI observations of supernova 1993J at 24 epochs, from 50 days after shock breakout to the present, allowed us to determine the coordinates of the explosion center relative to the quasi-stationary core of the host galaxy M81 with an accuracy of 45 micro-arcsec, and to determine the nominal proper motion of the geometric center of the radio shell with an accuracy of 9micro-arcsec/yr. The uncertainties correspond to 160 AU for the position and 160 km/s for the proper motion at the distance of the source of 3.63 Mpc. After correcting for the expected galactic proper motion of the supernova around the core of M81 using HI rotation curves, we obtain a peculiar proper motion of the radio shell center of only 320 +/- 160 km/s to the south, which limits any possible one-sided expansion of the shell. We also find that the shell is highly circular, the outer contours in fact being circular to within 3%. Combining our proper motion values with the degree of circular symmetry, we find that the expansion of the shockfront from the explosion center is isotropic to within 5.5% in the plane of the sky. This is a more fundamental result on isotropic expansion than can be derived from the circularity of the images alone. The brightness of the radio shell, however, varies along the ridge and systematically changes with time. The degree of isotropy in the expansion of the shockfront contrasts with the asymmetries and polarization found in optical spectral lines. Asymmetric density distributions in the ejecta or more likely in the circumstellar medium, are favored to reconcile the radio and optical results. We see no sign of any disk-like density distribution of the circumstellar material, with the average axis ratio of the radio shell of SN1993J being less than 1.04.

The nucleus of M81 was observed at 8.3 GHz with VLBI at 20 epochs over 4.5yrs, with a linear resolution at the source of about 2000 AU or 0.01 pc. Phase-referenced mapping with respect to the geometric center of supernova 1993J enabled us to find, with a standard error of about 600 AU, a stationary point in the source. We identify this point as the location of the core and the putative black hole at the gravitational center of the galaxy. The upper bound on the core's average velocity on the sky is <730 km/s. A short, one-sided jet extends from the core. The orientation of the jet varies smoothly, with timescales of about 1yr and an rms of 6degrees about the mean of 50degrees. Occasionally the jet appears to bend to the east. The length of the jet is only about 1 mas (3,600 AU), and varies with an rms of about 20%. The inferred speeds are below 0.08c. The total flux density of the core-jet varies erratically, changing on occasion by a factor of two over a few weeks, without any significant changes in the source size and orientation. The inferred velocity of the plasma flow is >0.25c. The results are consistent with a model in which plasma condensations with short lifetimes are ejected relativistically from the core on a timescale of less than a few weeks, and then travel along a tube whose pattern and geometry are also variable but only on a timescale of about one year. The central engine of M81 has qualitative similarities to those of powerful the AGN of radio galaxies and quasars, and may also represent in power and size a scaled-up version of the largely hidden nucleus in our own Galaxy.

M81 is the nearest galaxy with an active galactic nucleus, with the possibly exception of Cen A. The nucleus is exceptionally compact. Here we present results from new VLA and VLBI observations of the nucleus of M81. The VLA shows that its flux density and spectral index are variable by 50% on a timescale of weeks. The VLBI imaging results show that there is structure on scales <0.5 mas, and that this structure is variable, also on a timescale of weeks.