SN1993J VLBI (I): The Center of the Explosion and a Limit on Anisotropic Expansion

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.Comment: 21 pages, LaTex + 5 Figures (encapulsated PostScript), Accepted for Publication in the Astrophysical Journa

The Distribution of High and Low Redshift Type Ia Supernovae

The distribution of high redshift Type Ia supernovae (SNe Ia) with respect to projected distance from the center of the host galaxy is studied and compared to the distribution of local SNe. The distribution of high-z SNe Ia is found to be similar to the local sample of SNe Ia discovered with CCDs, but different than the sample discovered photographically. This is shown to be due to the Shaw effect. These results have implications for the use of SNe Ia to determine cosmological parameters if the local sample of supernovae used to calibrate the light curve decline relationships is drawn from a sample discovered photographically. A K-S test shows that the probability that the high redshift SNe of the Supernova Cosmology Project are drawn from the same distribution as the low redshift calibrators of Riess et al. is 0.1%. This is a potential problem because photographically discovered SNe are preferentially discovered farther away from the galaxy nucleus, where SNe show a lower scatter in absolute magnitude, and are on average 0.3 magnitudes fainter than SNe located closer to the center of their host galaxy. This raises questions about whether or not the calibration SNe sample the full range of parameters potentially present in high redshift SNe Ia. The limited data available suggest that the calibration process is adequate; however, it would be preferable if high redshift SNe and the low redshift SNe used to calibrate them were drawn from the same sample, as subtle differences may be important. Data are also presented which suggest that the seeming anti-Malmquist trend noticed by Tammann et al.(1996, 1998) for SNe Ia in galaxies with Cepheid distances may be due to the location of the SNe in their host galaxies.Comment: 16 pages, 2 figures Accepted for publication in the Astrophysical Journa

Thermonuclear Burning Regimes and the Use of SNe Ia in Cosmology

The calculations of the light curves of thermonuclear supernovae are carried out by a method of multi-group radiation hydrodynamics. The effects of spectral lines and expansion opacity are taken into account. The predictions for UBVI fluxes are given. The values of rise time for B and V bands found in our calculations are in good agreement with the observed values. We explain why our results for the rise time have more solid physical justification than those obtained by other authors. It is shown that small variations in the chemical composition of the ejecta, produced in the explosions with different regimes of nuclear burning, can influence drastically the light curve decline in the B band and, to a lesser extent, in the V band. We argue that recent results on positive cosmological constant Lambda, found from the high redshift supernova observations, could be wrong in the case of possible variations of the preferred mode of nuclear burning in the earlier Universe.Comment: 20 pages, 5 figures, presented at the conference "Astronomy at the Eve of the New Century", Puschino, May 17-22, 1999. A few references and a table added, typos correcte

Extinction and the Radial Distribution of Supernova Properties in Their Parent Galaxies

We use a Monte Carlo technique and assumed spatial distributions of dust and supernova (SN) progenitors in a simple model of a characteristic SN--producing disk galaxy to explore the effects of extinction on the radial distributions of SN properties in their parent galaxies. The model extinction distributions and projected radial number distributions are presented for various SN types. Even though the model has no core-collapse SNe within three kpc of the center, a considerable fraction of the core-collapse SNe are projected into the inner regions of inclined parent galaxies owing to their small vertical scale height. The model predicts that because of extinction, SNe projected into the central regions should on average appear dimmer and have a much larger magnitude scatter than those in the outer regions. In particular, the model predicts a strong deficit of bright core-collapse events inside a projected radius of a few kpc. Such a deficit is found to be present in the observations. It is a natural consequence of the characteristic spatial distributions of dust and core-collapse SNe in galaxies, and it leads us to offer an alternative to the conventional interpretation of the Shaw effect.Comment: 14 pages, 4 figure

A comparative modeling of supernova 1993J

The light curve of Supernova (SN) 1993J is calculated using two approaches to radiation transport as exemplified by the two computer codes, STELLA and EDDINGTON. Particular attention is paid to shock breakout and the photometry in the U, B, and V bands during the first 120 days. The hydrodynamical model, the explosion of a 13 Msun star which had lost most of its hydrogenic envelope to a companion, is the same in each calculation. The comparison elucidates differences between the approaches and also serves to validate the results of both. STELLA includes implicit hydrodynamics and is able to model supernova evolution at early times, before the expansion is homologous. STELLA also employs multi-group photonics and is able to follow the radiation as it decouples from the matter. EDDINGTON uses a different algorithm for integrating the transport equation, assumes homologous expansion, and uses a finer frequency resolution. Good agreement is achieved between the two codes only when compatible physical assumptions are made about the opacity. A new result for SN 1993J is a prediction of the continuum spectrum near the shock breakout (calculated by STELLA) which is superior to the results of other standard single energy group hydrocodes such as VISPHOT or TITAN. Based on the results of our independent codes, we discuss the uncertainties involved in the current time dependent models of supernova light curves.Comment: 43 pages with 22 eps figures, aaspp4.sty + epsf.sty, Accepted by ApJ, to appear in March 20, 1998 issue, Vol. 49

SN 1993J VLBI (III): The Evolution of the Radio Shell

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.Comment: Accepted for publication in the Astrophysical Journal. 34 Pages, 14 Figures (Figure 3 in color

A Parameter Study of Type II Supernova Light Curves Using 6 M_odot He Cores

Results of numerical calculations of Type II supernova light curves are presented. The model progenitor stars have 6 $M{_\odot}$ cores and various envelopes, originating from a numerically evolved 20 $M{_\odot}$ star. Five parameters that affect the light curves are examined: the ejected mass, the progenitor radius, the explosion energy, the $^{56}$Ni mass, and the extent of $^{56}$Ni mixing. The following affects have been found: 1) the larger the progenitor radius the brighter the early--time light curve, with little affect on the late--time light curve, 2) the larger the envelope mass the fainter the early light curve and the flatter the slope of the late light curve, 3) the larger the explosion energy the brighter the early light curve and the steeper the slope of the late light curve, 4) the larger the $^{56}$Ni mass the brighter the overall light curve after 20 to 50 days, with no affect on the early light curve, 5) the more extensive the $^{56}$Ni mixing the brighter the early light curve and the steeper the late light curve. The primary parameters affecting the light curve shape are the progenitor radius and the ejected mass. The secondary parameters are the explosion energy, $^{56}$Ni mass and $^{56}$Ni mixing. I find that while in principle the general shape and absolute magnitude of a light curve indicate a unique set of parameters, in practice it is difficult to avoid some ambiguity in the parameters. I find that the nickel--powered diffusion wave and the recombination of helium produce a prominent secondary peak in all our calculations. The feature is less prominent when compositional mixing, both $^{56}$Ni mixing and mixing between the hydrogen and helium layers, occurs. The model photospheric temperatures and velocities are presented, for comparison to observation.Comment: 39 pages, 15 figures. Astrophysical Journal (Accepted, Dec. 20, 2004