The discovery of 150 - 300 M$_{\odot}$ stars in the Local Group and pair-instability supernova candidates at low redshifts has excited interest in this exotic explosion mechanism. Realistic light curves for pair-instability supernovae at near-solar metallicities are key to identifying and properly interpreting these events as more are found. We have modeled pair-instability supernovae of 150 - 500 M$_{\odot}$ Z $\sim$ 0.1 - 0.4 Z$_{\odot}$ stars. These stars lose up to 80% of their mass to strong line-driven winds and explode as bare He cores. We find that their light curves and spectra are quite different from those of Population III pair-instability explosions, which therefore cannot be used as templates for low-redshift events. Although non-zero metallicity pair-instability supernovae are generally dimmer than their Population III counterparts, in some cases they will be bright enough to be detected at the earliest epochs at which they can occur, the formation of the first galaxies at $z \sim$ 10 - 15. Others can masquerade as dim, short duration supernovae that are only visible in the local universe and that under the right conditions could be hidden in a wide variety of supernova classes. We also report for the first time that some pair-instability explosions can create black holes with masses of $\sim$ 100 M$_{\odot}$.

The existence of supermassive black holes as early as z ~ 7 is one of the great unsolved problems in cosmological structure formation. One leading theory argues that they are born during catastrophic baryon collapse in z ~ 15 protogalaxies in strong Lyman-Werner UV backgrounds. Atomic line cooling in such galaxies fragments baryons into massive clumps that are thought to directly collapse to 10^4 - 10^5 solar-mass black holes. We have now discovered that some of these fragments can instead become supermassive stars that eventually explode as pair-instability supernovae with energies of ~ 10^55 erg, the most energetic explosions in the universe. We have calculated light curves and spectra for supermassive Pop III PI SNe with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible in NIR all-sky surveys by Euclid out to z ~ 10 - 15 and by WFIRST and WISH out to z ~ 15 - 20, perhaps revealing the birthplaces of the first quasars.

Understanding the properties of Pop III stars is prerequisite to elucidating the nature of primeval galaxies, the chemical enrichment and reionization of the early IGM, and the origin of supermassive black holes. While the primordial IMF remains unknown, recent evidence from numerical simulations and stellar archaeology suggests that some Pop III stars may have had lower masses than previously thought, 15 - 50 \Ms in addition to 50 - 500 \Ms. The detection of Pop III supernovae by JWST, WFIRST or the TMT could directly probe the primordial IMF for the first time. We present numerical simulations of 15 - 40 \Ms Pop III core-collapse SNe done with the Los Alamos radiation hydrodynamics code RAGE. We find that they will be visible in the earliest galaxies out to z ~ 10 - 15, tracing their star formation rates and in some cases revealing their positions on the sky. Since the central engines of Pop III and solar-metallicity core-collapse SNe are quite similar, future detection of any Type II supernovae by next-generation NIR instruments will in general be limited to this epoch.

The observational signatures of the first cosmic explosions and their chemical imprint on second-generation stars both crucially depend on how heavy elements mix within the star at the earliest stages of the blast. We present numerical simulations of the early evolution of Population III pair-instability supernovae with the new adaptive mesh refinement code CASTRO. In stark contrast to 15 - 40 Msun core-collapse primordial supernovae, we find no mixing in most 150 - 250 Msun pair-instability supernovae out to times well after breakout from the surface of the star. This may be the key to determining the mass of the progenitor of a primeval supernova, because vigorous mixing will cause emission lines from heavy metals such as Fe and Ni to appear much sooner in the light curves of core-collapse supernovae than in those of pair-instability explosions. Our results also imply that unlike low-mass Pop III supernovae, whose collective metal yields can be directly compared to the chemical abundances of extremely metal-poor stars, further detailed numerical simulations will be required to determine the nucleosynthetic imprint of very massive Pop III stars on their direct descendants.

We present multidimensional simulations of the post-explosion hydrodynamics in three different 15 solar mass supernova models with zero, 10^{-4} solar metallicity, and solar metallicities. We follow the growth of the Rayleigh-Taylor instability that mixes together the stellar layers in the wake of the explosion. Models are initialized with spherically symmetric explosions and perturbations are seeded by the grid. Calculations are performed in two-dimensional axisymmetric and three-dimensional Cartesian coordinates using the new Eulerian hydrodynamics code, CASTRO. We find as in previous work, that Rayleigh-Taylor perturbations initially grow faster in 3D than in 2D. As the Rayleigh-Taylor fingers interact with one another, mixing proceeds to a greater degree in 3D than in 2D, reducing the local Atwood number and slowing the growth rate of the instability in 3D relative to 2D. By the time mixing has stopped, the width of the mixed region is similar in 2D and 3D simulations provided the Rayleigh-Taylor fingers show significant interaction. Our results imply that 2D simulations of light curves and nucleosynthesis in supernovae (SNe) that die as red giants may capture the features of an initially spherically symmetric explosion in far less computational time than required by a full 3D simulation. However, capturing large departures from spherical symmetry requires a significantly perturbed explosion. Large scale asymmetries cannot develop through an inverse cascade of merging Rayleigh-Taylor structures; they must arise from asymmetries in the initial explosion.

We present a new code, CASTRO, that solves the multicomponent compressible hydrodynamic equations for astrophysical flows including self-gravity, nuclear reactions and radiation. CASTRO uses an Eulerian grid and incorporates adaptive mesh refinement (AMR). Our approach to AMR uses a nested hierarchy of logically-rectangular grids with simultaneous refinement in both space and time. The radiation component of CASTRO will be described in detail in the next paper, Part II, of this series.

The inclusion of rotationally-induced mixing in stellar evolution can alter the structure and composition of presupernova stars. We survey the effects of progenitor rotation on nucleosynthetic yields in Population III and II supernovae using the new adaptive mesh refinement (AMR) code CASTRO. We examine spherical explosions in 15, 25 and 40 solar mass stars at Z = 0 and 10^-4 solar metallicity with three explosion energies and two rotation rates. Rotation in the Z = 0 models resulted in primary nitrogen production and a stronger hydrogen burning shell which led all models to die as red supergiants. On the other hand, the Z=10^-4 solar metallicity models that included rotation ended their lives as compact blue stars. Because of their extended structure, the hydrodynamics favors more mixing and less fallback in the metal free stars than the Z = 10^-4 models. As expected, higher energy explosions produce more enrichment and less fallback than do lower energy explosions, and less massive stars produce more enrichment and leave behind smaller remnants than do more massive stars. We compare our nucleosynthetic yields to the chemical abundances in the three most iron-poor stars yet found and reproduce the abundance pattern of one, HE 0557-4840, with a zero metallicity 15 solar mass, 2.4 x 10^51 erg supernova. A Salpeter IMF averaged integration of our yields for Z=0 models with explosion energies of 2.4x10^51 ergs or less is in good agreement with the abundances observed in larger samples of extremely metal-poor stars, provided 15 solar mass stars are included. Since the abundance patterns of extremely metal-poor stars likely arise from a representative sample of progenitors, our yields suggest that low-mass supernovae contributed the bulk of the metals to the early universe.

Two-dimensional simulations of mixing and fall back in non-rotating massive stars have been carried out using realistic initial models for the presupernova star and assuming standard spherically symmetric explosions of 1.2*10^51ergs. Stars of 15 and 25 solar masses with both primordial and solar composition were modeled. The zero metallicity supernova progenitors were compact blue stars and the amount of Rayleigh-Taylor induced mixing in them was greatly reduced compared with what was seen in the red supergiants with solar metallicity. The compact zero-metal stars also experienced more fallback than their solar metallicity counterparts. As a result, the ejected nucleosynthesis from the two populations was very different. For the simple explosion model assumed, low metallicity stars ejected ejected too little iron and intermediate mass elements even to explain the abundance patterns in the most iron-poor stars found to date, suggesting that some important ingredient is missing. Rotation is likely to alter these conclusions by producing a greater fraction of red supergiants among Population III stars. The velocities of the heavy elements in all models considered - both red and blue supergiants - were less than observed in SN 1987A, suggesting that at least occasionally, asymmetric aspects of the explosion mechanism and fallback play a major role in mixing.