FleCSPH is a smoothed particle hydrodynamics simulation tool, based on the compile-time configurable framework FleCSI. The asynchronous distributed tree topology combined with a fast multipole method allows FleCSPH to efficiently compute hydrodynamics and long range particle-particle interactions. FleCSPH provides initial data generators, particle relaxation techniques, and standard evolution drivers, which can be easily modified and extended to user-specific setups. Data input/output uses the H5part format, compatible with modern visualization software.

The mergers of compact binaries with at least one neutron star component are the potential leading sites of the production and ejection of $r$-process elements. Discoveries of galactic binary pulsars, short gamma-ray bursts, and gravitational-wave detections have all been constraining the rate of these events, while the gravitational wave plus broadband electromagnetic coverage of binary neutron star merger (GW170817) has also placed constraints on the properties (mass and composition) of the merger ejecta. But uncertainties and ambiguities in modeling the optical and infrared emission make it difficult to definitively measure the distribution of heavy isotopes in these mergers. In contrast, gamma rays emitted in the decay of these neutron-rich ejecta may provide a more direct measurement of the yields. We calculate the gamma production in remnants of neutron star mergers, considering two epochs: a kilonova epoch, lasting about two weeks, and a much later epoch of tens and hundreds of thousands of years after the merger. For the kilonova epoch, when the expanding ejecta is still only partially transparent to gamma radiation, we use 3D radiative transport simulations to produce the spectra. We show that the gamma-ray spectra associated with beta- and alpha-decay provide a fingerprint of the ejecta properties and, for a sufficiently nearby remnant, may be detectable, even for old remnants. We compare our gamma spectra with the potential detection limits of next generation detectors, including the Lunar Occultation Explorer (LOX), the All-sky Medium Energy Gamma-ray Observatory (AMEGO), and the Compton Spectrometer and Imager (COSI). We show that fission models can be discriminated via the presence of short-lived fission fragments in the remnant spectra.

The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g. in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide "curtaining", which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent P Cygni features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances.

Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters -- such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key parameters on explosive nucleosynthesis and final explosion yields. Post-bounce temperature and density evolution of ZAMS 15, 20, and 25 solar mass progenitors are post-processed through the Nucleosynthesis Grid (NuGrid) nuclear network to obtain detailed explosive yields. In particular, this study focuses on radio-isotopes that are of particular interest to the next generation of gamma-ray astronomical observations; 43K, 47Ca, 44Sc, 47Sc, 48V, 48Cr, 51Cr, 52Mn, 59Fe, 56Co, 57Co, 57Ni. These nuclides may be key in advancing our understanding of the inner workings of core-collapse supernovae by probing the parameters of the explosion engine. We find that the isotopes that are strong indicators of explosion energy are 43K, 47Ca, 44Sc, 47Sc, and 59Fe, those that are dependent on the progenitor structure are 48V, 51Cr, and 57Co, and those that probe neither are 48Cr, 52Mn, 57Ni, and 56Co. We discuss prospects of observing these radionuclides in supernova remnants.

We explore the impact of pulsar electromagnetic dipole and fallback accretion emission on the luminosity of a suite of kilonova models. The pulsar models are varied over pulsar magnetic field strength, pulsar lifetime, ejecta mass, and elemental abundances; the fallback models are varied over fallback accretion rate and ejecta mass. For the abundances, we use Fe and Nd as representatives of the wind and dynamical ejecta, respectively. We simulate radiative transfer in the ejecta in either 1D spherical or 2D cylindrical spatial geometry. For the grid of 1D simulations, the mass fraction of Nd is 0, $10^{-4}$, or $10^{-3}$ and the rest is Fe. Our models that fit the bolometric luminosity of AT 2017gfo (the kilonova associated with the first neutron star merger discovered in gravitational waves, GW170817) do not simultaneously fit the B, V, and I time evolution. However, we find that the trends of the evolution in B and V magnitudes are better matched by the fallback model relative to the pulsar model, implying the time dependence of the remnant source influences the color evolution. Further exploration of the parameter space and model deficiencies is needed before we can describe AT 2017gfo with a remnant source.

The merger of neutron star binaries is believed to eject a wide range of heavy elements into the universe. By observing the emission from this ejecta, scientists can probe the ejecta properties (mass, velocity and composition distributions). The emission (a.k.a. kilonova) is powered by the radioactive decay of the heavy isotopes produced in the merger and this emission is reprocessed by atomic opacities to optical and infra-red wavelengths. Understanding the ejecta properties requires calculating the dependence of this emission on these opacities. The strong lines in the optical and infra-red in lanthanide opacities have been shown to significantly alter the light-curves and spectra in these wavelength bands, arguing that the emission in these wavelengths can probe the composition of this ejecta. Here we study variations in the kilonova emission by varying individual lanthanide (and the actinide uranium) concentrations in the ejecta. The broad forest of lanthanide lines makes it difficult to determine the exact fraction of individual lanthanides. Nd is an exception. Its opacities above 1 micron are higher than other lanthanides and observations of kilonovae can potentially probe increased abundances of Nd. Similarly, at early times when the ejecta is still hot (first day), the U opacity is strong in the 0.2-1 micron wavelength range and kilonova observations may also be able to constrain these abundances.

Pair-Instability Supernovae (PISNe) may signal the deaths of extremely massive stars in the local Universe or massive primordial stars after the end of the Cosmic Dark Ages. Hydrodynamic simulations of these explosions, performed in 1D, 2D, and 3D geometry, have revealed the strong dependence of mixing in the PISN ejecta on dimensionality. This chemical rearrangement is mainly driven by Rayleigh-Taylor instabilities that start to grow shortly after the collapse of the carbon-oxygen core. We investigate the effects of such mixing on the spectroscopic evolution of PISNe by post-processing explosion profiles with the radiation diffusion-equilibrium code SNEC and the implicit Monte Carlo-discrete diffusion Monte Carlo (IMC-DDMC) radiation transport code SuperNu. The first 3D radiation transport calculation of a PISN explosion is presented yielding viewing angle-dependent synthetic spectra and lightcurves. We find that while 2D and 3D mixing does not significantly affect the lightcurves of PISNe, their spectroscopic and color evolution is impacted. Strong features of intermediate mass elements dominated by silicon, magnesium and oxygen appear at different phases and reach different intensities depending on the extent of mixing in the silicon/oxygen interface of the PISN ejecta. On the other hand, we do not find a significant dependence of PISN lightcurves and spectra on viewing angle. Our results showcase the capabilities of SuperNu to handle 3D radiation transport and highlight the importance of modeling time-series of spectra in identifying PISNe with future missions.

We investigate $^{60}$Fe in massive stars and core-collapse supernovae focussing on uncertainties that influence its production in 15, 20 and 25 $M_\odot$ stars at solar metallicity. We find that the $^{60}$Fe yield is a monotonic increasing function of the uncertain $^{59}$Fe$(n,\gamma)^{60}$Fe cross section and that a factor of 10 reduction in the reaction rate results in a factor 8-10 reduction in the $^{60}$Fe yield; while a factor of 10 increase in the rate increases the yield by a factor 4-7. We find that none of the 189 simulations we have performed are consistent with a core-collapse supernova triggering the formation of the Solar System, and that only models using $^{59}$Fe$(n,\gamma)^{60}$Fe cross section that is less than or equal to that from NON-SMOKER can reproduce the observed $^{60}$Fe/$^{26}$Al line flux ratio in the diffuse ISM. We examine the prospects of detecting old core-collapse supernova remnants (SNRs) in the Milky Way from their $\gamma$-ray emission from the decay of $^{60}$Fe, finding that the next generation of gamma-ray missions could be able to discover up to $\sim100$ such old SNRs as well as measure the $^{60}$Fe yields of a handful of known Galactic SNRs. We also predict the X-ray spectrum that is produced by atomic transitions in $^{60}$Co following its ionization by internal conversion and give theoretical X-ray line fluxes as a function of remnant age as well as the Doppler and fine-structure line broadening effects. The X-ray emission presents an interesting prospect for addressing the missing SNR problem with future X-ray missions.

Neutron star mergers offer unique conditions for the creation of the heavy elements and additionally provide a testbed for our understanding of this synthesis known as the $r$-process. We have performed dynamical nucleosynthesis calculations and identified a single isotope, $^{254}$Cf, which has a particularly high impact on the brightness of electromagnetic transients associated with mergers on the order of 15 to 250 days. This is due to the anomalously long half-life of this isotope and the efficiency of fission thermalization compared to other nuclear channels. We estimate the fission fragment yield of this nucleus and outline the astrophysical conditions under which $^{254}$Cf has the greatest impact to the light curve. Future observations in the middle-IR which are bright during this regime could indicate the production of actinide nucleosynthesis.

We present the results of an investigation of the dredge-up and mixing during the merger of two white dwarfs with different chemical compositions by conducting hydrodynamic simulations of binary mergers for three representative mass ratios. In all the simulations, the total mass of the two white dwarfs is $\lesssim1.0~{\rm M_\odot}$. Mergers involving a CO and a He white dwarf have been suggested as a possible formation channel for R Coronae Borealis type stars, and we are interested in testing if such mergers lead to conditions and outcomes in agreement with observations. Even if the conditions during the merger and subsequent nucleosynthesis favor the production of $^{18}{\mathrm O}$, the merger must avoid dredging up large amounts of $^{16}{\mathrm O}$, or else it will be difficult to produce sufficient $^{18}{\mathrm O}$ to explain the oxygen ratio observed to be of order unity. We performed a total of 9 simulations using two different grid-based hydrodynamics codes using fixed and adaptive meshes, and one smooth particle hydrodynamics (SPH) code. We find that in most of the simulations, $>10^{-2}~{\rm M_\odot}$ of $^{16}{\mathrm O}$ is indeed dredged up during the merger. However, in SPH simulations where the accretor is a hybrid He/CO white dwarf with a $\sim 0.1~{\rm M_\odot}$ layer of helium on top, we find that no $^{16}{\mathrm O}$ is being dredged up, while in the $q=0.8$ simulation $<10^{-4}~{\rm M_\odot}$ of $^{16}{\mathrm O}$ has been brought up, making a WD binary consisting of a hybrid CO/He WD and a companion He WD an excellent candidate for the progenitor of RCB stars.

Core-collapse supernova science is now entering an era where engine models are beginning to make both qualitative and quantitative predictions. Although the evidence in support of the convective engine for core-collapse supernova continues to grow, it is difficult to place quantitative constraints on this engine. Some studies have made specific predictions for the remnant distribution from the convective engine, but the results differ between different groups. Here we use a broad parameterization for the supernova engine to understand the differences between distinct studies. With this broader set of models, we place error bars on the remnant mass and basic yields from the uncertainties in the explosive engine. We find that, even with only 3 progenitors and a narrow range of explosion energies, we can produce a wide range of remnant masses and nucleosynthetic yields.

We report the discovery and monitoring of the near-infrared counterpart (AT2017gfo) of a binary neutron-star merger event detected as a gravitational wave source by Advanced LIGO/Virgo (GW170817) and as a short gamma-ray burst by Fermi/GBM and Integral/SPI-ACS (GRB170817A). The evolution of the transient light is consistent with predictions for the behaviour of a "kilonova/macronova", powered by the radioactive decay of massive neutron-rich nuclides created via r-process nucleosynthesis in the neutron-star ejecta. In particular, evidence for this scenario is found from broad features seen in Hubble Space Telescope infrared spectroscopy, similar to those predicted for lanthanide dominated ejecta, and the much slower evolution in the near-infrared Ks-band compared to the optical. This indicates that the late-time light is dominated by high-opacity lanthanide-rich ejecta, suggesting nucleosynthesis to the 3rd r-process peak (atomic masses A~195). This discovery confirms that neutron-star mergers produce kilo-/macronovae and that they are at least a major - if not the dominant - site of rapid neutron capture nucleosynthesis in the universe.

With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the binary neutron star merger GW170817. The bright, rapidly fading ultraviolet emission indicates a high mass ($\approx0.03$ solar masses) wind-driven outflow with moderate electron fraction ($Y_{e}\approx0.27$). Combined with the X-ray limits, we favor an observer viewing angle of $\approx 30^{\circ}$ away from the orbital rotation axis, which avoids both obscuration from the heaviest elements in the orbital plane and a direct view of any ultra-relativistic, highly collimated ejecta (a gamma-ray burst afterglow).

To assess the effectiveness of optical emission as a probe of spatial asymmetry in core-collapse supernovae (CCSNe), we apply the radiative transfer software, \supernu, to a unimodal CCSN model. The \snsph\ radiation-hydrodynamics software was used to simulate an asymmetric explosion of a 16 M$_{\odot}$ ZAMS binary star. The ejecta has 3.36 M$_{\odot}$ with 0.024 M$_{\odot}$ of radioactive $^{56}$Ni, with unipolar asymmetry along the z-axis. For 96 discrete angular views, we find the ratio between maximum and minimum peak total luminosities is $\sim$1.36. The brightest light curves emerge from views orthogonal to the z-axis. Multigroup spectra from UV to IR are obtained. We find a shift in wavelength with viewing angle in a near-IR Ca II emission feature, consistent with Ca being mostly in the unimode. We compare emission from the grey gamma-ray transfer in \supernu\ and from the detailed gamma-ray transfer code \maverick. Relative to the optical light curves, the brightness of the gamma-ray emission is more monotonic with respect to viewing angle. UBVRI broad-band light curves are also calculated. Parallel with the unimode, the U and B bands have excess luminosity at $\gtrsim 10$ days post-explosion, due to $^{56}$Ni on the unimode. We compare our CCSN model with SN 2002ap, which is thought to have a similar ejecta morphology.

The electromagnetic transients accompanying compact binary mergers ($\gamma$-ray bursts, afterglows and 'macronovae') are crucial to pinpoint the sky location of gravitational wave sources. Macronovae are caused by the radioactivity from freshly synthesised heavy elements, e.g. from dynamic ejecta and various types of winds. We study macronova signatures by using multi-dimensional radiative transfer calculations. We employ the radiative transfer code SuperNu and state-of-the-art LTE opacities for a few representative elements from the wind and dynamical ejecta (Cr, Pd, Se, Te, Br, Zr, Sm, Ce, Nd, U) to calculate synthetic light curves and spectra for a range of ejecta morphologies. The radioactive power of the resulting macronova is calculated with the detailed input of decay products. We assess the detection prospects for our most complex models, based on the portion of viewing angles that are sufficiently bright, at different cosmological redshifts ($z$). The brighter emission from the wind is unobscured by the lanthanides (or actinides) in some of the models, permitting non-zero detection probabilities for redshifts up to $z=0.07$. We also find the nuclear mass model and the resulting radioactive heating rate are crucial for the detectability. While for the most pessimistic heating rate (from the FRDM model) no reasonable increase in the ejecta mass or velocity, or wind mass or velocity, can possibly make the light curves agree with the observed nIR excess after GRB130603B, a more optimistic heating rate (from the Duflo-Zuker model) leads to good agreement. We conclude that future reliable macronova observations would constrain nuclear heating rates, and consequently help constrain nuclear mass models.

Currently there is great interest in the outcomes and astrophysical implications of mergers of double degenerate binaries. In a commonly adopted approximation, the components of such binaries are represented by polytropes with an index n=3/2. We present detailed comparisons of stellar mass-transfer and merger simulations of polytropic binaries that have been carried out using two very different numerical algorithms --- a finite-volume "grid" code and a smoothed-particle hydrodynamics (SPH) code. We find that there is agreement in both the ultimate outcomes of the evolutions and the intermediate stages if the initial conditions for each code are chosen to match as closely as possible. We find that even with closely matching initial setups, the time it takes to reach a concordant evolution differs between the two codes because the initial depth of contact cannot be matched exactly. There is a general tendency for SPH to yield higher mass transfer rates and faster evolution to the final outcome. We also present comparisons of simulations calculated from two different energy equations: in one series we assume a polytropic equation of state and in the other series an ideal gas equation of state. In the latter series of simulations an atmosphere forms around the accretor which can exchange angular momentum and cause a more rapid loss of orbital angular momentum. In the simulations presented here, the effect of the ideal equation of state is to de-stabilize the binary in both SPH and grid simulations, but the effect is more pronounced in the grid code.

The SuperNovae Analysis aPplication (SNAP) is a new tool for the analysis of SN observations and validation of SN models. SNAP consists of an open source relational database with (a) observational light curve, (b) theoretical light curve, and (c) correlation table sets, statistical comparison software, and a web interface available to the community. The theoretical models are intended to span a gridded range of parameter space. The goal is to have users to upload new SN models or new SN observations and run the comparison software to determine correlations via the web site. There are looming problems on the horizon that SNAP begins to solve. Namely, large surveys will discover thousands of SNe annually. Frequently, the parameter space of a new SN event is unbounded. SNAP will be a resource to constrain parameters and determine if an event needs follow-up without spending resources to create new light curve models from scratch. Secondly, there is not a rapidly available, systematic way to determine degeneracies between parameters or even what physics is needed to model a realistic SNe. The correlations made within the SNAP system begin to solve these problems.

With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before the explosion. Based on current observations of SLSNe, the PISN origin has been disfavoured for a number of reasons. Many PISN models provide overly broad light curves and too reddened spectra, because of massive ejecta and a high amount of nickel. In the current study we re-examine PISN properties using progenitor models computed with the GENEC code. We calculate supernova explosions with FLASH and light curve evolution with the radiation hydrodynamics code STELLA. We find that high-mass models (200 and 250 solar masses) at relatively high metallicity (Z=0.001) do not retain hydrogen in the outer layers and produce relatively fast evolving PISNe Type I and might be suitable to explain some SLSNe. We also investigate uncertainties in light curve modelling due to codes, opacities, the nickel-bubble effect and progenitor structure and composition.

The recent discovery of the unprecedentedly superluminous transient ASASSN-15lh (or SN 2015L) with its UV-bright secondary peak challenges all the power-input models that have been proposed for superluminous supernovae. Here we examine some of the few viable interpretations of ASASSN-15lh in the context of a stellar explosion, involving combinations of one or more power inputs. We model the lightcurve of ASASSN-15lh with a hybrid model that includes contributions from magnetar spin-down energy and hydrogen-poor circumstellar interaction. We also investigate models of pure circumstellar interaction with a massive hydrogen-deficient shell and discuss the lack of interaction features in the observed spectra. We find that, as a supernova ASASSN-15lh can be best modeled by the energetic core-collapse of a ~40 Msun star interacting with a hydrogen-poor shell of ~20 Msun. The circumstellar shell and progenitor mass are consistent with a rapidly rotating pulsational pair-instability supernova progenitor as required for strong interaction following the final supernova explosion. Additional energy injection by a magnetar with initial period of 1-2 ms and magnetic field of 0.1-1 x 10^14 G may supply the excess luminosity required to overcome the deficit in single-component models, but this requires more fine-tuning and extreme parameters for the magnetar, as well as the assumption of efficient conversion of magnetar energy into radiation. We thus favor a single-input model where the reverse shock formed in a strong SN ejecta-CSM interaction following a very powerful core-collapse SN explosion can supply the luminosity needed to reproduce the late-time UV-bright plateau.

Population III supernovae have been of growing interest of late for their potential to directly probe the properties of the first stars, particularly the most energetic events that are visible near the edge of the observable universe. But until now, hypernovae, the unusually energetic Type Ib/c supernovae that are sometimes associated with gamma-ray bursts, have been overlooked as cosmic beacons at the highest redshifts. In this, the latest of a series of studies on Population III supernovae, we present numerical simulations of 25 - 50 M$_{\odot}$ hypernovae and their light curves done with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible at z = 10 - 15 to the James Webb Space Telescope (JWST) and z = 4 - 5 to the Wide-Field Infrared Survey Telescope (WFIRST), tracing star formation rates in the first galaxies and at the end of cosmological reionization. If, however, the hypernova crashes into a dense shell ejected by its progenitor, it is expected that a superluminous event will occur that may be seen at z ~ 20, in the first generation of stars.

The era of the universe's first (Population III) stars is essentially unconstrained by observation. Ultra-luminous and massive stars from this time altered the chemistry of the cosmos, provided the radiative scaffolding to support the formation of the first protogalaxies, and facilitated the creation and growth of now-supermassive black holes. Unfortunately, because these stars lie literally at the edge of the observable universe, they will remain beyond the reach of even the next generation of telescopes such as the James Webb Space Telescope and the Thirty-Meter Telescope. In this paper, we provide a primer to supernovae modeling and the first stars to make our discussion accessible to those new to or outside our field. We review recent work of the Los Alamos Supernova Light Curve Project and Brigham Young University to explore the possibility of probing this era through observations of the spectacular deaths of the first stars. We find that many such brilliant supernova explosions will be observable as far back as $\sim 99$% of the universe's current age, tracing primordial star formation rates and the locations of their protogalaxies on the sky. The observation of Population III supernovae will be among the most spectacular discoveries in observational astronomy in the coming decade.

Mass loss remains one of the primary uncertainties in stellar evolution. In the most massive stars, mass loss dictates the circumstellar medium and can significantly alter the fate of the star. Mass loss is caused by a variety of wind mechanisms and also through binary interactions. Supernovae are excellent probes of this mass loss, both the circumstellar material and the reduced mass of the hydrogen-rich envelope. In this paper, we focus on the effects of reducing the hydrogen-envelope mass on the supernova light curve, studying both the shock breakout and peak light curve emission for a wide variety of mass loss scenarios. Even though the trends of this mass loss will be masked somewhat by variations caused by different progenitors, explosion energies, and circumstellar media, these trends have significant effects on the supernova light-curves that should be seen in supernova surveys. We conclude with a comparison of our results to a few key observations.

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}$.

Population III stars that die as pair-instability supernovae are usually thought to fall in the mass range of 140 - 260 M$_{\odot}$. But several lines of work have now shown that rotation can build up the He cores needed to encounter the pair instability at stellar masses as low as 90 $_{\odot}$. Depending on the slope of the initial mass function of Population III stars, there could be 4 - 5 times as many stars from 90 - 140 $_{\odot}$ in the primordial universe than in the usually accepted range. We present numerical simulations of the pair-instability explosions of such stars performed with the MESA, FLASH and RAGE codes. We find that they will be visible to supernova factories such as Pan-STARRS and LSST in the optical out to z $\sim$ 1 - 2 and to JWST and the 30 m-class telescopes in the NIR out to $z \sim$ 7 - 10. Such explosions will thus probe the stellar populations of the first galaxies and cosmic star formation rates in the era of cosmological reionization. These supernovae are also easily distinguished from more massive pair-instability explosions, underscoring the fact that there is far greater variety to the light curves of these events than previously understood.

Supernovae are Nature's high-energy, high density laboratory experiments, reaching densities in excess of nuclear densities and temperatures above 10MeV. Astronomers have built up a suite of diagnostics to study these supernovae. If we can utilize these diagnostics, and tie them together with a theoretical understanding of supernova physics, we can use these cosmic explosions to study the nature of matter at these extreme densities and temperatures. Capitalizing on these diagnostics will require understanding a wide range of additional physics. Here we review the diagnostics and the physics needed to use them to learn about the supernova engine, and ultimate nuclear physics.

Population III supernovae have been the focus of growing attention because of their potential to directly probe the properties of the first stars, particularly the most energetic events that can be seen at the edge of the observable universe. But until now pair-pulsation supernovae, in which explosive thermonuclear burning in massive stars fails to unbind them but can eject their outer layers into space, have been overlooked as cosmic beacons at the earliest redshifts. These shells can later collide and, like Type IIn supernovae, produce superluminous events in the UV at high redshifts that could be detected in the near infrared today. We present numerical simulations of a 110 M$_{\odot}$ pair-pulsation explosion done with the Los Alamos radiation hydrodynamics code RAGE. We find that collisions between consecutive pair pulsations are visible in the near infrared out to z $\sim$ 15 - 20 and can probe the earliest stellar populations at cosmic dawn.

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.

One of the leading contenders for the origin of supermassive black holes at $z \gtrsim$ 7 is catastrophic baryon collapse in atomically-cooled halos at $z \sim$ 15. In this scenario, a few protogalaxies form in the presence of strong Lyman-Werner UV backgrounds that quench H$_2$ formation in their constituent halos, preventing them from forming stars or blowing heavy elements into the intergalactic medium prior to formation. At masses of 10$^ 8$ \Ms\ and virial temperatures of 10$^4$ K, gas in these halos rapidly cools by H lines, in some cases forming 10$^4$ - 10$^6$ \Ms\ Pop III stars and, a short time later, the seeds of supermassive black holes. Instead of collapsing directly to black holes some of these stars died in the most energetic thermonuclear explosions in the universe. We have modeled the explosions of such stars in the dense cores of line-cooled protogalaxies in the presence of cosmological flows. In stark contrast to the explosions in diffuse regions in previous simulations, these SNe briefly engulf the protogalaxy but then collapse back into its dark matter potential. Fallback drives turbulence that efficiently distributes metals throughout the interior of the halo and fuels the rapid growth of nascent black holes at its center. The accompanying starburst and x-ray emission from these line-cooled galaxies easily distinguish them from more slowly evolving neighbors and might reveal the birthplaces of supermassive black holes on the sky.

The first stars are the key to the formation of primitive galaxies, early cosmological reionization and chemical enrichment, and the origin of supermassive black holes. Unfortunately, in spite of their extreme luminosities, individual Population III stars will likely remain beyond the reach of direct observation for decades to come. However, their properties could be revealed by their supernova explosions, which may soon be detected by a new generation of NIR observatories such as JWST and WFIRST. We present light curves and spectra for Pop III pair-instability supernovae calculated with the Los Alamos radiation hydrodynamics code RAGE. Our numerical simulations account for the interaction of the blast with realistic circumstellar envelopes, the opacity of the envelope, and Lyman absorption by the neutral IGM at high redshift, all of which are crucial to computing the NIR signatures of the first cosmic explosions. We find that JWST will detect pair-instability supernovae out to z > 30, WFIRST will detect them in all-sky surveys out to z ~ 15 - 20 and LSST and Pan-STARRS will find them at z ~ 7 - 8. The discovery of these ancient explosions will probe the first stellar populations and reveal the existence of primitive galaxies that might not otherwise have been detected.

Supermassive primordial stars are expected to form in a small fraction of massive protogalaxies in the early universe, and are generally conceived of as the progenitors of the seeds of supermassive black holes (BHs). Supermassive stars with masses of ~55,000 M_Sun, however, have been found to explode and completely disrupt in a supernova (SN) with an energy of up to ~10^55 erg instead of collapsing to a BH. Such events, ~10,000 times more energetic than typical SNe today, would be among the biggest explosions in the history of the universe. Here we present a simulation of such a SN in two stages. Using the RAGE radiation hydrodynamics code we first evolve the explosion from an early stage through the breakout of the shock from the surface of the star until the blast wave has propagated out to several parsecs from the explosion site, which lies deep within an atomic cooling dark matter (DM) halo at z ~ 15. Then, using the GADGET cosmological hydrodynamics code we evolve the explosion out to several kiloparsecs from the explosion site, far into the low-density intergalactic medium. The host DM halo, with a total mass of 4 x 10^7 M_Sun, much more massive than typical primordial star-forming halos, is completely evacuated of high density gas after < 10 Myr, although dense metal-enriched gas recollapses into the halo, where it will likely form second-generation stars with metallicities of ~ 0.05 Z_Sun after > 70 Myr. The chemical signature of supermassive star explosions may be found in such long-lived second-generation stars today.

The first primitive galaxies formed from accretion and mergers by z ~ 15, and were primarily responsible for cosmological reionization and the chemical enrichment of the early cosmos. But a few of these galaxies may have formed in the presence of strong Lyman-Werner UV fluxes that sterilized them of H_2, preventing them from forming stars or expelling heavy elements into the IGM prior to assembly. At masses of 10^8 Ms and virial temperatures of 10^4 K, these halos began to rapidly cool by atomic lines, perhaps forming 10^4 - 10^6 Ms Pop III stars and, later, the seeds of supermassive black holes. We have modeled the explosion of a supermassive Pop III star in the dense core of a line-cooled protogalaxy with the ZEUS-MP code. We find that the supernova (SN) expands to a radius of ~ 1 kpc, briefly engulfing the entire galaxy, but then collapses back into the potential well of the dark matter. Fallback fully mixes the interior of the protogalaxy with metals, igniting a violent starburst and fueling the rapid growth of a massive black hole at its center. The starburst would populate the protogalaxy with stars in greater numbers and at higher metallicities than in more slowly-evolving, nearby halos. The SN remnant becomes a strong synchrotron source that can be observed with eVLA and eMERLIN and has a unique signature that easily distinguishes it from less energetic SN remnants. Such explosions, and their attendant starbursts, may well have marked the birthplaces of supermassive black holes on the sky.

(Abridged) Some supernovae (SNe) may be powered by the interaction of the SN ejecta with a large amount of circumstellar matter (CSM). Here we outline a method to measure the mass of the optically thick CSM around such SNe. We present observations of SN2010jl, including the first detection of a SN using NuSTAR. The total radiated luminosity of SN2010jl is extreme, at least 9e50 erg. By modeling the visible-light data, we robustly show that the mass of the circumstellar material within ~1e16 cm of the progenitor was in excess of 10 solar masses, likely ejected tens of years prior to the SN explosion. Our modeling suggests that the shock velocity during shock breakout was ~6000 km/s, decelerating to ~2600 km/s about two years after maximum light. Our late-time NuSTAR+XMM spectra of the SN presumably provide the first direct measurement of SN shock velocity two years after the SN maximum light -- measured to be in the range of 2000 to 4500 km/s if the ions and electrons are in equilibrium, and >~2000 km/s if they are not in equilibrium. This measurement is in agreement with the shock velocity predicted by our modeling of the optical data. We also show that the mean radial density distribution of the CSM roughly follows an r^-2 law. A possible explanation for the massive CSM with a wind-like profile is that they are the result of multiple pulsational pair instability events prior to the SN explosion, separated from each other by years.

The detection of Pop III supernovae could directly probe the primordial IMF for the first time, unveiling the properties of the first galaxies, early chemical enrichment and reionization, and the seeds of supermassive black holes. Growing evidence that some Pop III stars were less massive than 100 solar masses may complicate prospects for their detection, because even though they would have been more plentiful they would have died as core-collapse supernovae, with far less luminosity than pair-instability explosions. This picture greatly improves if the SN shock collides with a dense circumstellar shell ejected during a prior violent LBV type eruption. Such collisions can turn even dim SNe into extremely bright ones whose luminosities can rival those of pair-instability SNe. We present simulations of Pop III Type IIn SN light curves and spectra performed with the Los Alamos RAGE and SPECTRUM codes. Taking into account Lyman-alpha absorption in the early universe and cosmological redshifting, we find that 40 solar mass Pop III Type IIn SNe will be visible out to z ~ 20 with JWST and out to z ~ 7 with WFIRST. Thus, even low mass Pop III SNe can be used to probe the primeval universe.

We have entered the era of explosive transient astronomy, in which upcoming real-time surveys like the Large Synoptic Survey Telescope (LSST), the Palomar Transient Factory (PTF) and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) will detect supernovae in unprecedented numbers. Future telescopes such as the James Webb Space Telescope may discover supernovae from the earliest stars in the universe and reveal their masses. The observational signatures of these astrophysical transients are the key to unveiling their central engines, the environments in which they occur, and to what precision they will pinpoint cosmic acceleration and the nature of dark energy. We present a new method for modeling supernova light curves and spectra with the radiation hydrodynamics code RAGE coupled with detailed monochromatic opacities in the SPECTRUM code. We include a suite of tests that demonstrate how the improved physics is indispensable to modeling shock breakout and light curves.

Observations with the Swift UVOT have unambiguously uncovered for the first time a long-lived, UV "plateau" in a Type II-P supernova (SN). Although this flattening in slope is hinted at in a few other SNe, due to its proximity and minimal line-of-sight extinction, SN 2012aw has afforded the first opportunity to clearly observe this UV plateau. The observations of SN 2012aw revealed all Swift UV and u-band lightcurves initially declined rapidly, but 27 days after explosion the light curves flattened. Some possible sources of the UV plateau are: the same thermal process that cause the optical plateau, heating from radioactive decay, or a combination of both processes.

The first stars ended the cosmic Dark Ages and created the first heavy elements necessary for the formation of planets and life. The properties of these stars remain uncertain, and it may be decades before individual Pop III stars are directly observed. Their masses, however, can be inferred from their supernova explosions, which may soon be found in both deep-field surveys by JWST and in all-sky surveys by WFIRST. We have performed radiation hydrodynamical simulations of the near infrared signals of Pop III pair-instability supernovae in realistic circumstellar environments with Lyman absorption by the neutral intergalactic medium. We find that JWST and WFIRST will detect these explosions out to z ~ 30 and 20, respectively, unveiling the first generation of stars in the universe.

A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the results of five 3-dimensional hydrodynamic simulations of the merger of a double white dwarf system where the total mass is 0.9 Mdot and the initial mass ratio (q) varies between 0.5 and 0.99. We identify in simulations with $q\lesssim0.7$ a feature around the merged stars where the temperatures and densities are suitable for forming 18O. However, more 16O is being dredged-up from the C- and O-rich accretor during the merger than the amount of 18O that is produced. Therefore, on a dynamical time scale over which our hydrodynamics simulation runs, a 16O/18O ratio of ~2000 in the "best" case is found. If the conditions found in the hydrodynamic simulations persist for 10^6 seconds the oxygen ratio drops to 16 in one case studied, while in a hundred years it drops to ~4 in another case studied, consistent with the observed values in RCB stars. Therefore, the merger of two white dwarfs remains a strong candidate for the formation of these enigmatic stars.

We have developed a self-consistent-field technique similar to the one described by Hachisu, Eriguchi, & Nomoto (1986b) that can be used to construct detailed force-balanced models of synchronously rotating, double white dwarf (DWD) binaries that have a wide range of total masses, mass ratios, and separations. In addition to providing a computational tool that can be used to provide quiet initial starts for dynamical studies of the onset of mass transfer in DWD systems, we show that this SCF technique can be used to construct model sequences that mimic the last portion of the detached inspiral phase of DWD binary evolutions, and semi-detached model sequences that mimic a phase of conservative mass transfer.