The DJEHUTY project is an intensive effort at the Lawrence Livermore National Laboratory (LLNL) to produce a general purpose 3-D stellar structure and evolution code to study dynamic processes in whole stars.

We present 30 optical spectra and 49 photometric epochs sampling the first 517 days after discovery of supernova (SN) 1999em, and derive its distance through the expanding photosphere method (EPM). SN 1999em is shown to be a Type II-plateau (II-P) event, with a photometric plateau lasting until about 100 days after explosion. We identify the dominant ions responsible for most of the absorption features seen in the optical portion of the spectrum during the plateau phase. We find the distance to SN 1999em to be D = 8.2 +/- 0.6 Mpc, with an explosion date of 5.3 +/- 1.4 days before discovery. We examine potential sources of systematic error in EPM-derived distances, and find the most significant to result from uncertainty in the theoretical modeling of the flux distribution emitted by the SN photosphere (i.e., the ``flux dilution factor''). We compare previously derived EPM distances to 5 SNe II in galaxies for which a recently revised Cepheid distance exists from the HST Key Project and find D(Cepheids) / D(EPM) = 0.96 +/- 0.09. Finally, we investigate the possible use of SNe II-P as standard candles and find that for 8 photometrically confirmed SNe II-P with previously derived EPM distances and SN 1999em, the mean plateau absolute brightness is M_V(plateau) = -16.4 +/- 0.6 mag, implying that distances good to ~30% (1-sigma) may be possible without the need for a complete EPM analysis. At M_V(plateau) = -15.9 +/- 0.2 mag, SN 1999em is somewhat fainter than the average SN II-P. The general consistency of absolute SNe II-P brightness during the plateau suggests that the standard candle assumption may allow SNe II-P to be viable cosmological beacons at z > 2.

We present optical and IR spectroscopy of the first two months of evolution of the Type II SN 1999em. We combine these data with high-quality optical/IR photometry beginning only three days after shock breakout, in order to study the performance of the ``Expanding Photosphere Method'' (EPM) in the determination of distances. With this purpose we develop a technique to measure accurate photospheric velocities by cross-correlating observed and model spectra. The application of this technique to SN 1999em shows that we can reach an average uncertainty of 11% in velocity from an individual spectrum. Our analysis shows that EPM is quite robust to the effects of dust. In particular, the distances derived from the VI filters change by only 7% when the adopted visual extinction in the host galaxy is varied by 0.45 mag. The superb time sampling of the BVIZJHK light-curves of SN 1999em permits us to study the internal consistency of EPM and test the dilution factors computed from atmosphere models for Type II plateau supernovae. We find that, in the first week since explosion, the EPM distances are up to 50% lower than the average, possibly due the presence of circumstellar material. Over the following 65 days, on the other hand, our tests lend strong credence to the atmosphere models, and confirm previous claims that EPM can produce consistent distances without having to craft specific models to each supernova. This is particularly true for the VI filters which yield distances with an internal consistency of 4%. From the whole set of BVIZJHK photometry, we obtain an average distance of 7.5+/-0.5 Mpc, where the quoted uncertainty (7%) is a conservative estimate of the internal precision of the method obtained from the analysis of the first 70 days of the supernova evolution.

Currently popular models for Type Ia supernovae (SNe Ia) fall into two general classes. The first comprises explosions of nearly pure carbon/oxygen (C/O) white dwarfs at the Chandrasekhar limit which ignite near their centers. The second consists of lower-mass C/O cores which are ignited by the detonation of an accreted surface helium layer. Explosions of the latter type produce copious Fe, Co and Ni K-alpha emission from 56Ni and 56Co decay in the detonated surface layers, emission which is much weaker from Chandrasekhar-mass models. The presence of this emission provides a simple and unambiguous discriminant between these two models for SNe Ia. Both mechanisms may produce 0.1-0.6 solar masses of 56Ni, making them bright gamma-ray line emitters. The time to maximum brightness of 56Ni decay lines is distinctly shorter in the sub-Chandrasekhar mass class of model (approximately 15 days) than in the Chandrasekhar mass model (approximately 30 days), making gamma-ray line evolution another direct test of the explosion mechanism. It should just be possible to detect K-shell emission from a sub-Chandrasekhar explosion from SNe Ia as far away as the Virgo cluster with the XMM Observatory. A 1 to 2 square meter X-ray telescope such as the proposed Con-X Observatory could observe K-alpha emission from sub-Chandrasekhar mass SNe Ia in the Virgo cluster, providing not just a detection, but high-accuracy flux and kinematic information.

Observations of type Ia supernovae at high redshifts have become an important tool for studying the geometry of the universe. The relation between the duration of the peak phase of a type Ia supernova's lightcurve and its luminosity (broader is brighter) forms the cornerstone of this measurement, yet it is a purely empirical relation. In this paper we show that the relation is a natural consequence of the radiation transport in type Ia supernovae and suggest constraints on the nature of the explosions which arise from our interpretation of the observed relation. The principle parameter underlying the relation is the mass of radioactive 56Ni produced in the explosion. The relation is shown to be relatively insensitive to most other parameters in Chandrasekhar-mass explosions. Cosmological results are thus unlikely to suffer from systematic effects stemming from evolution in the explosions' progenitors.

We have developed an implicit, multi-group, time-dependent, spherical neutrino transport code based on the Feautrier variables, the tangent-ray method, and accelerated ${\bf \Lambda}$ iteration. The code achieves high angular resolution, is good to O($v/c$), is equivalent to a Boltzmann solver (without gravitational redshifts), and solves the transport equation at all optical depths with precision. In this paper, we present our formulation of the relevant numerics and microphysics and explore protoneutron star atmospheres for snapshot post-bounce models. Our major focus is on spectra, neutrino-matter heating rates, Eddington factors, angular distributions, and phase-space occupancies. In addition, we investigate the influence on neutrino spectra and heating of final-state electron blocking, stimulated absorption, velocity terms in the transport equation, neutrino-nucleon scattering asymmetry, and weak magnetism and recoil effects. Furthermore, we compare the emergent spectra and heating rates obtained using full transport with those obtained using representative flux-limited transport formulations to gauge their accuracy and viability. Finally, we derive useful formulae for the neutrino source strength due to nucleon-nucleon bremsstrahlung and determine bremsstrahlung's influence on the emergent $\nu_{\mu}$ and $\nu_{\tau}$ neutrino spectra.

Recently a Type Ic supernova, SN 1998bw, was discovered coincident with a gamma-ray burst, GRB 980425. The supernova had unusual radio, optical, and spectroscopic properties. We explore here models based upon helium stars in the range 9 - 14 Msun and carbon-oxygen stars 6 - 11 Msun which experience unusually energetic explosions (kinetic energy 0.5 - 2.8 x 10**52 erg). Bolometric light curves and multi-band photometry are calculated and compared favorably with observations. Both LTE and non-LTE spectra are also calculated for the model that agrees best with the light curve, a carbon-oxygen core of 6 Msun exploded with a kinetic energy of 2.2 x 10**52 erg. We also examine potential mechanisms for producing the observed gamma-ray burst (GRB) - shock break-out and relativistic shock deceleration in circumstellar material. For spherically symmetric models, both fail to produce a GRB of even the low luminosity inferred for GRB 980425. The high explosion energies required to understand the supernova are in contrast to what is expected for such massive stars and may indicate that a new sort of explosion has been identified, possibly the consequence of a collapsar (Woosley 1993, 1996). Indeed a more likely explanation for what was seen is a highly asymmetric explosion in which the GRB was produced by a relativistic jet, perhaps viewed obliquely, and only a fraction of the total stellar mass was ejected, the remainder accreting into a black hole. The ejected mass (but not the 56Ni mass), explosion energy, and velocities may then be smaller. Other Type Ic supernova - GRB coincidences are suggested.

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.

We examine the basic physics of type Ia supernova (SNe Ia) light curves with a view toward interpreting the relations between peak luminosity, peak width, and late-time slope in terms of the properties of the underlying explosion models. We develop an analytic solution of the radiation transport problem in the dynamic diffusion regime and show that it reproduces bolometric light curves produced by more detailed calculations under the assumption of a constant extinction coefficient. This model is used to derive the thermal conditions in the interior of SNe Ia, to address the issue of time dependence in the modeling of SNe Ia light curves and spectra and to show that these are intrinsically time-dependent phenomena. The analytic model is then used to study the sensitivity of light curves to various properties of supernova explosions. We show that the dominant opacity arises from line transitions and discuss the nature of line opacities in expanding media and the appropriateness of various mean opacities used in light curve calculations.

We have used observations gathered at CTIO to measure distances by the Expanding Photosphere Method (EPM) to 5 Type II supernovae. These supernovae lie at redshifts from cz = 1100 km/s to z = 5500 km/s, and increase to 18 the number of distances measured using EPM. We compare distances derived to 11 Type II supernovae with distances to their host galaxies measured using the Tully-Fisher method. We find that the Tully-Fisher distances average 11%+-7% smaller. The comparison shows no significant evidence of any large distance-dependent bias in the Tully-Fisher distances. We employ the sample of EPM distances from 4.5~Mpc to 180~Mpc to derive a value for the Hubble Constant. We find H_0 = 73 +- 6 (statistical) +- 7 (systematic}) km/s/Mpc.

We use the Expanding Photosphere Method to determine distances to 10 type II supernovae. The effects of asymmetries, extinction, and flux dilution are explored. Using empirical evidence and time-independent, spherical models which treat H and He in non-LTE, we show that blackbody corrections caused by flux dilution are small for type II supernovae in the infrared, and in the optical when their color temperatures are less than 6000~K. The extinction to a type II-P supernova can be estimated from its light curve: the uncertainty introduced into a distance measurement due to extinction is usually less than 10\%. Correcting for extinction and flux dilution we derive distances to 10 supernovae: SN 1968L, SN 1969L, SN 1970G, SN 1973R, SN 1979C, SN 1980K, SN 1987A, SN 1988A, SN 1990E, and SN 1990ae. The distance measurements span a wide range, 50 kpc to 120 Mpc, which is unique among the methods for establishing the extragalactic distance scale. The distances measured to SN 1970G in M101 and SN 1987A in the LMC are in good agreement with distances determined from Cepheid variable stars. Our distance to the Virgo Cluster, 22 +- 3 Mpc, is larger than recent distances estimates made using surface brightness fluctuations, planetary nebula luminosity functions, and the Tully-Fisher method. Using the distances determined from these type II supernovae we derive a value of $H_0 = 60 \pm 10$ km sec$^{-1}$Mpc$^{-1}$. This value is subject to errors caused by local deviations in the Hubble flow, but will soon be improved by applying the Expanding Photosphere Method to several distant type II supernovae.