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.

The optical/UV light curves of SN 1987A are analyzed with the multi-energy group radiation hydrodynamics code STELLA. The calculated monochromatic and bolometric light curves are compared with observations shortly after shock breakout, during the early plateau, through the broad second maximum, and during the earliest phase of the radioactive tail. We have concentrated on a progenitor model calculated by Nomoto & Hashimoto and Saio, Nomoto, & Kato, which assumes that 14 solar masses of the stellar mass is ejected. Using this model, we have updated constraints on the explosion energy and the extent of mixing in the ejecta. In particular, we determine the most likely range of E/M (explosion energy over ejecta mass) and R_0 (radius of the progenitor). In general, our best models have energies in the range E = (1.1 +/- 0.3) x 10^{51} ergs, and the agreement is better than in earlier, flux-limited diffusion calculations for the same explosion energy. Our modeled B and V fluxes compare well with observations, while the flux in U undershoots after about 10 days by a factor of a few, presumably due to NLTE and line transfer effects. We also compare our results with IUE observations, and a very good quantitative agreement is found for the first days, and for one IUE band (2500-3000 A) as long as for 3 months. We point out that the V flux estimated by McNaught & Zoltowski should probably be revised to a lower value.Comment: 27 pages AASTeX v.4.0 + 35 postscript figures. ApJ, accepte

Aims.We present synthetic bolometric and broad-band UBVRI light curves of SNe Ia for four selected 3D deflagration models of thermonuclear supernovae. Methods.The light curves are computed with the 1D hydro code stella, which models (multi-group time-dependent) non-equilibrium radiative transfer inside SN ejecta. Angle-averaged results from 3D hydrodynamical explosion simulations with the composition determined in a nucleosynthetic postprocessing step served as the input to the radiative transfer model.Results.The predicted model light curves do agree reasonably well with the observed ones for SNe Ia in the range of low to normal luminosities, although the underlying hydrodynamical explosion models produced only a modest amount of radioactive Ni(i.e. ~0.24–0.42 ) and relatively low kinetic energy in the explosion (less than erg). The evolution of predicted B and V fluxes in the model with a Nimass of 0.42 follows the observed decline rate after the maximum very well, although the behavior of fluxes in other filters deviates somewhat from observations, and the bolometric decline rate is a bit slow. The material velocity at the photospheric level is on the order of 104 km s-1 for all models. Using our models, we check the validity of Arnett's rule, relating the peak luminosity to the power of the deposited radioactive heating, and we also check the accuracy of the procedure for extracting the Nimass from the observed light curves.Conclusions.We find that the comparison between theoretical light curves and observations provides a useful tool to validate SN Ia models. The steps necessary for improving the agreement between theory and observations are set out.

Various types of expansions in series of Chebyshev-Hermite polynomials currently used in astrophysics for weakly non-normal distributions are compared, namely the Gram-Charlier, Gauss-Hermite and Edgeworth expansions. It is shown that the Gram-Charlier series is most suspect because of its poor convergence properties. The Gauss-Hermite expansion is better but it has no intrinsic measure of accuracy. The best results are achieved with the asymptotic Edgeworth expansion. We draw attention to the form of this expansion found by Petrov for arbitrary order of the asymptotic parameter and present a simple algorithm realizing Petrov's prescription for the Edgeworth expansion. The results are illustrated by examples similar to the problems arising when fitting spectral line profiles of galaxies, supernovae, or other stars, and for the case of approximating the probability distribution of peculiar velocities in the cosmic string model of structure formation.Comment: 13 pages with 11 eps figures, aa.cls + graphics packages, Submitted to Astronomy & Astrophysics Supplement on July 16, 1997, accepted November 10, 199