We explore the sensitivity of the nucleosynthesis of intermediate mass elements (28 < A < 80) in supernovae derived from massive stars to the nuclear reaction rates employed in the model. Two standard sources of reaction rate data (Woosley et al. 1978; and Thielemann et al. 1987) are employed in pairs of calculations that are otherwise identical. Both include as a common backbone the experimental reactions rates of Caughlan & Fowler (1988). Two stellar models are calculated for each of two main sequence masses: 15 and 25 solar masses. Each star is evolved from core hydrogen burning to a presupernova state carrying an appropriately large reaction network, then exploded using a piston near the edge of the iron core as described by Woosley & Weaver (1995). The final stellar yields from the models calculated with the two rate sets are compared and found to differ in most cases by less than a factor of two over the entire range of nuclei studied. Reasons for the major discrepancies are discussed in detail along with the physics underlying the two reaction rate sets employed. The nucleosynthesis results are relatively robust and less sensitive than might be expected to uncertainties in nuclear reaction rates, though they are sensitive to the stellar model employed.

Ten years later, astronomers are still puzzled by the stellar evolution that produced SN 1987A --- a blue supergiant. In single star models, the new OPAL opacities make blue solutions more difficult to achieve, though still possible for certain choices of convection physics. We also consider rotation, which has the desirable effect of producing large surface enhancements of nitrogen and helium, but the undesirable effect of increasing the helium-core mass at the expense of the envelope. The latter makes blue solutions more difficult. Still, we seek a model that occurs with high probability in the LMC and for which the time-scale for making the last transition from red to blue, $\sim$ 20,000 years, has a physical interpretation --- the Kelvin-Helmholtz time of the helium core. Single star models satisfy both criteria and might yet prove to be the correct explanation for Sk -69 202, provided new rotational or convection physics can simultaneously give a blue star and explain the ring structure. Some speculations on how this might be achieved are presented and some aspects of binary models briefly discussed.

Using recently calculated models for massive stellar evolution and supernovae coupled to a model for Galactic chemical evolution, neutron star and black hole birth functions (number of neutron stars and black holes as a function of their mass) are determined for the Milky Way Galaxy. For those stars that explode as Type II supernovae, the models give birth functions that are bimodal with peaks at 1.27 and 1.76 M\sun \ and average masses within those peaks of 1.28 and 1.73 M\sun. For those stars that explode as Type Ib there is a narrower spread of remnant masses, the average being 1.32 M\sun, and less evidence for bimodality. These values will be increased, especially in the more massive Type II supernovae, if significant accretion continues during the initial launching of the shock, and the number of heavier neutron stars could be depleted by black hole formation. The principal reason for the dichotomy in remnant masses for Type II is the difference in the presupernova structure of stars above and below 19 M\sun, the mass separating stars that burn carbon convectively from those that produce less carbon and burn radiatively. The Type Ib's and the lower mass group of the Type II's compare favorably with measured neutron star masses, and in particular to the Thorsett et al. (1993) determination of the average neutron star mass in 17 systems; 1.35 $\pm$ 0.27 M\sun. Variations in the exponent of a Salpeter initial mass function are shown not to affect the locations of the two peaks in the distribution function, but do affect their relative amplitudes. Sources of uncertainty, in particular placement of the mass cut and sensitivity to the explosion energy, are discussed, and estimates of the total number of neutron stars and black holes in the Galaxy are given. Accretion induced collapse should give a unique gravitational mass of 1.27

Using recently calculated yields for Type II supernovae, along with models for chemical evolution and the distribution of mass in the interstellar medium, the current abundances and spatial distributions of two key gamma-ray radioactivities, $^{26}$Al and $^{60}$Fe, are determined. The estimated steady state production rates are 2.0 $\pm$ 1.0 M\sun \ Myr$^{-1}$ for $^{26}$Al and 0.75 $\pm$ 0.4 M\sun \ Myr$^{-1}$ for $^{60}$Fe. This corresponds to 2.2 $\pm$ 1.1 M\sun \ of $^{26}$Al and 1.7 $\pm$ 0.9 M\sun \ of $^{60}$Fe in the present interstellar medium. Sources of uncertainty are discussed, one of the more important being the current rate of core collapse supernovae in the Galaxy. Our simple model gives three per century, but reasonable changes in the star formation rate could easily accommodate a core collapse rate one-half as large, and thus one-half the yields. When these stellar and chemical evolution results are mapped into a three dimensional model of the Galaxy, the calculated 1809 keV gamma-ray flux map is consistent with the {\it Compton Gamma Ray Observatory} observations of a steep decline in the flux outside a longitude of $\pm$ 50$^\circ$ from the Galactic center, and the slight flux enhancements observed in the vicinity of spiral arms. Other potential stellar sources of $^{26}$Al and $^{60}$Fe are mentioned, especially the possibility of $^{60}$Fe synthesis in Type Ia supernovae. Predictions for the $^{60}$Fe mass distribution, total mass, and flux map are given.

Using the output from a grid of 60 Type II supernova models (Woosley \& Weaver 1994) of varying mass (11 $\ltaprx$ M/M\sun $\ltaprx$ 40) and metallicity (0, 10$^{-4}$, 0.01, 0.1, and 1 Z\sol), the chemical evolution of 76 stable isotopes, from hydrogen to zinc, is calculated. The chemical evolution calculation employs a simple dynamical model for the Galaxy (infall with a 4 billion year $e$-folding time scale onto a exponential disk and 1/r$^2$ bulge), and standard evolution parameters, such as a Salpeter initial mass function and a quadratic Schmidt star formation rate. The theoretical results are compared in detail with observed stellar abundances in stars with metallicities in the range -3.0 $\ltaprx$ [Fe/H] $\ltaprx$ 0.0 dex. While our discussion focuses on the solar neighborhood where there are the most observations, the supernovae rates, an intrinsically Galactic quantity, are also discussed.