We report on our observing campaign of the compact binary merger GW190814, detected by the Advanced LIGO and Advanced Virgo detectors on August 14th, 2019. This signal has the best localisation of any observed gravitational wave (GW) source, with a 90% probability area of 18.5 deg$^2$, and an estimated distance of ~ 240 Mpc. We obtained wide-field observations with the Deca-Degree Optical Transient Imager (DDOTI) covering 88% of the probability area down to a limiting magnitude of $w$ = 19.9 AB. Nearby galaxies within the high probability region were targeted with the Lowell Discovery Telescope (LDT), whereas promising candidate counterparts were characterized through multi-colour photometry with the Reionization and Transients InfraRed (RATIR) and spectroscopy with the Gran Telescopio de Canarias (GTC). We use our optical and near-infrared limits in conjunction with the upper limits obtained by the community to constrain the possible electromagnetic counterparts associated with the merger. A gamma-ray burst seen along its jet's axis is disfavoured by the multi-wavelength dataset, whereas the presence of a burst seen at larger viewing angles is not well constrained. Although our observations are not sensitive to a kilonova similar to AT2017gfo, we can rule out high-mass (> 0.1 M$_{\odot}$) fast-moving (mean velocity >= 0.3c) wind ejecta for a possible kilonova associated with this merger.

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.

LOX is a lunar-orbiting astrophysics mission that will probe the cosmos at MeV energies. It is guided by open questions regarding thermonuclear, or Type-Ia, supernovae (SNeIa) and will characterize these inherently radioactive objects by enabling a systematic survey of SNeIa at gamma-ray energies for the first time. Astronomical investigations from lunar orbit afford new opportunities to advance our understanding of the cosmos. The foundation of LOX is an observational approach well suited to the all-sky monitoring demands of supernova investigations and time-domain astronomy. Its inherently wide field-of-view and continuous all-sky monitoring provides an innovative way of addressing decadal survey questions at MeV energies (0.1-10 MeV). The LOX approach achieves high sensitivity with a simple, high-heritage instrument design that eliminates the need for complex, position-sensitive detectors, kinematic event reconstruction, masks, or other insensitive detector mass, while also mitigating technology development, implementation complexity, and their associated costs. LOX can be realized within existing programs, like Explorer.

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.

The electromagnetic observations of GW170817 were able to dramatically increase our understanding of neutron star mergers beyond what we learned from gravitational waves alone. These observations provided insight on all aspects of the merger from the nature of the gamma-ray burst to the characteristics of the ejected material. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called kilonovae or macronovae. Characteristics of the ejecta include large velocity gradients, relative to supernovae, and the presence of heavy $r$-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound-bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-binned opacities that preserve the integral of the opacity over frequency. Advantages of this area-preserving approach over the traditional expansion-opacity formalism include the ability to pre-calculate opacity tables that are independent of the type of hydrodynamic expansion and that eliminate the computational expense of calculating opacities within radiation-transport simulations. Tabular opacities are generated for all 14 lanthanides as well as a representative actinide element, uranium. We demonstrate that spectral simulations produced with the line-binned opacities agree well with results produced with the more accurate continuous Monte Carlo Sobolev approach, as well as with the commonly used expansion-opacity formalism. Additional investigations illustrate the convergence of opacity with respect to the number of included lines, and elucidate sensitivities to different atomic physics approximations, such as fully and semi-relativistic approaches.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

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.

Gravitational wave observations need accompanying electromagnetic signals to accurately determine the sky positions of the sources. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called macronovae. Characteristics of the ejecta include large velocity gradients and the presence of heavy $r$-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound-bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-smeared opacities that preserve the integral of the opacity over frequency, which is an alternative approach to the expansion-opacity formalism. The use of area-preserving line profiles produces frequency-dependent opacities that are one to two orders of magnitude greater than those obtained with the use of expansion opacities in the Sobolev approximation. Opacities are calculated for four $r$-process elements (cerium, neodymium, samarium and uranium) using fully and semi-relativistic methods, as well as different amounts of configuration interaction in the atomic structure, in order to test the sensitivity of the emission to the underlying atomic physics. We determine the effect of these new opacities on simulated spectra and broad-band light curves by applying a multi-dimensional ray-trace method to the ejecta predicted from 3-dimensional merger calculations. With our substantially larger opacities the simulations yield slightly lower luminosities that peak in the mid-IR. These results suggest that those radioactively powered transients that are related to the very heaviest $r$-process material are more difficult to observe than previously believed.

The hydrogen and helium accreted by X-ray bursting neutron stars is periodically consumed in runaway thermonuclear reactions that cause the entire surface to glow brightly in X-rays for a few seconds. With models of the emission, the mass and radius of the neutron star can be inferred from the observations. By simultaneously probing neutron star masses and radii, X-ray bursts are one of the strongest diagnostics of the nature of matter at extremely high densities. Accurate determinations of these parameters are difficult, however, due to the highly non-ideal nature of the atmospheres where X-ray bursts occur. Observations from X-ray telescopes such as RXTE and NuStar can potentially place strong constraints on nuclear matter once uncertainties in atmosphere models have been reduced. Here we discuss current progress on modeling atmospheres of X-ray bursting neutron stars and some of the challenges still to be overcome.

We provide a set of stellar evolution and nucleosynthesis calculations that applies established physics assumptions simultaneously to low- and intermediate-mass and massive star models. Our goal is to provide an internally consistent and comprehensive nuclear production and yield data base for applications in areas such as pre-solar grain studies. Our non-rotating models assume convective boundary mixing where it has been adopted before. We include 8 (12) initial masses for $Z = 0.01$ ($0.02$). Models are followed either until the end of the asymptotic giant branch phase or the end of Si burning, complemented by a simple analytic core-collapse supernova models with two options for fallback and shock velocities. The explosions show which pre-supernova yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impacts the light elements and the $s$ and $p$ process. For low- and intermediate-mass models our stellar yields from H to Bi include the effect of convective boundary mixing at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the $^{13}$C pocket. All post-processing nucleosynthesis calculations use the same nuclear reaction rate network and nuclear physics input. We provide a discussion of the nuclear production across the entire mass range organized by element group. All our stellar nucleosynthesis profile and time evolution output is available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced.

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.

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.

[Shortened] The 12C + 12C fusion reaction has been the subject of considerable experimental efforts to constrain uncertainties at temperatures relevant for stellar nucleosynthesis. In order to investigate the effect of an enhanced carbon burning rate on massive star structure and nucleosynthesis, new stellar evolution models and their yields are presented exploring the impact of three different 12C + 12C reaction rates. Non-rotating stellar models were generated using the Geneva Stellar Evolution Code and were later post-processed with the NuGrid Multi-zone Post-Processing Network tool. The enhanced rate causes core carbon burning to be ignited more promptly and at lower temperature. This reduces the neutrino losses, which increases the core carbon burning lifetime. An increased carbon burning rate also increases the upper initial mass limit for which a star exhibits a convective carbon core. Carbon shell burning is also affected, with fewer convective-shell episodes and convection zones that tend to be larger in mass. Consequently, the chance of an overlap between the ashes of carbon core burning and the following carbon shell convection zones is increased, which can cause a portion of the ashes of carbon core burning to be included in the carbon shell. Therefore, during the supernova explosion, the ejecta will be enriched by s-process nuclides synthesized from the carbon core s process. The yields were used to estimate the weak s-process component in order to compare with the solar system abundance distribution. The enhanced rate models were found to produce a significant proportion of Kr, Sr, Y, Zr, Mo, Ru, Pd and Cd in the weak component, which is primarily the signature of the carbon-core s process. Consequently, it is shown that the production of isotopes in the Kr-Sr region can be used to constrain the 12C + 12C rate using the current branching ratio for a- and p-exit channels.

The contribution by massive stars (M > 9 solar masses) to the weak s-process component of the solar system abundances is primarily due to the 22Ne neutron source, which is activated near the end of helium-core burning. The residual 22Ne left over from helium-core burning is then reignited during carbon burning, initiating further s-processing that modifies the isotopic distribution. This modification is sensitive to the stellar structure and the carbon burning reaction rate. Recent work on the 12C + 12C reaction suggests that resonances located within the Gamow peak may exist, causing a strong increase in the astrophysical S-factor and consequently the reaction rate. To investigate the effect of an increased rate, 25 solar mass stellar models with three different carbon burning rates, at solar metallicity, were generated using the Geneva Stellar Evolution Code (GENEC) with nucleosynthesis post-processing calculated using the NuGrid Multi-zone Post-Processing Network code (MPPNP). The strongest rate caused carbon burning to occur in a large convective core rather than a radiative one. The presence of this large convective core leads to an overlap with the subsequent convective carbon-shell, significantly altering the initial composition of the carbon-shell. In addition, an enhanced rate causes carbon-shell burning episodes to ignite earlier in the evolution of the star, igniting the 22Ne source at lower temperatures and reducing the neutron density.

We compare the yields of Ti44 and Ni56 produced from post-processing the thermodynamic trajectories from three different core-collapse models -- a Cassiopeia A progenitor, a double shock hypernova progenitor, and a rotating 2D explosion -- with the yields from exponential and power-law trajectories. The peak temperatures and densities achieved in these core-collapse models span several of the distinct nucleosynthesis regions we identify, resulting in different trends in the Ti44 and Ni56 yields for different mass elements. The Ti44 and Ni56 mass fraction profiles from the exponential and power-law profiles generally explain the tendencies of the post-processed yields, depending on which regions are traversed by the model. We find integrated yields of Ti44 and Ni56 from the exponential and power-law trajectories are generally within a factor 2 or less of the post-process yields. We also analyze the influence of specific nuclear reactions on the Ti44 and Ni56 abundance evolution. Reactions that affect all yields globally are the 3a, p(e-,nu)n and n(e+,nubar)p. The rest of the reactions are ranked according to their degree of impact on the synthesis of Ti44. The primary ones include Ti44(a,p)V47, Ca40(a,g)Ti44, V45(p,g)Cr46, Ca40(a,p)Sc43, F17(a,p)Ne20, Na21(a,p)Mg24, Sc41(p,g)Ti42, Sc43(p,g)Ti44, Ti44(p,g)V45, and Ni57(p,g)Cu58, along with numerous weak reactions. Our analysis suggests that not all Ti44 need be produced in an a-rich freeze-out in core-collapse events, and that reaction rate equilibria in combination with timescale effects for the expansion profile may account for the paucity of Ti44 observed in supernovae remnants.

The merger of two white dwarfs (a.k.a. double degenerate merger) has often been cited as a potential progenitor of type Ia supernovae. Here we combine population synthesis, merger and explosion models with radiation-hydrodynamics light-curve models to study the implications of such a progenitor scenario on the observed type Ia supernova population. Our standard model, assuming double degenerate mergers do produce thermonuclear explosions, produces supernova light-curves that are broader than the observed type Ia sample. In addition, we discuss how the shock breakout and spectral features of these double degenerate progenitors will differ from the canonical bare Chandrasekhar-massed explosion models. We conclude with a discussion of how one might reconcile these differences with current observations.

The slow neutron capture process in massive stars (the weak s-process) produces most of the s-only isotopes in the mass region 60 < A < 90. The nuclear reaction rates used in simulations of this process have a profound effect on the final s-process yields. We generated 1D stellar models of a 25 solar mass star varying the 12C + 12C rate by a factor of 10 and calculated full nucleosynthesis using the post-processing code PPN. Increasing or decreasing the rate by a factor of 10 affects the convective history and nucleosynthesis, and consequently the final yields.

Astronomers have proposed a number of mechanisms to produce supernova explosions. Although many of these mechanisms are now not considered primary engines behind supernovae, they do produce transients that will be observed by upcoming ground-based surveys and NASA satellites. Here we present the first radiation-hydrodynamics calculations of the spectra and light curves from three of these "failed" supernovae: supernovae with considerable fallback, accretion induced collapse of white dwarfs, and energetic helium flashes (also known as type .Ia supernovae).

We investigate the physical conditions where 44Ti and 56Ni are created in core-collapse supernovae. In this preliminary work we use a series of post-processing network calculations with parametrized expansion profiles that are representative of the wide range of temperatures, densities and electron-to-baryon ratios found in 3D supernova simulations. Critical flows that affect the final yields of 44Ti and 56Ni are assessed.

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation. We also stress the dynamical importance of including shocks in the simulations. These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae. In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase.

We review some of the uncertainties in calculating nucleosynthetic yields, focusing on the explosion mechanism. Current yield calculations tend to either use a piston, energy injection, or enhancement of neutrino opacities to drive an explosion. We show that the energy injection, or more accurately, an entropy injection mechanism is best-suited to mimic our current understanding of the convection-enhanced supernova engine. The enhanced neutrino-opacity technique is in qualitative disagreement with simulations of core-collapse supernovae and will likely produce errors in the yields. But piston-driven explosions are the most discrepant. Piston-driven explosion severely underestimate the amount of fallback, leading to order-of-magnitude errors in the yields of heavy elements. To obtain yields accurate to the factor of a few level, we must use entropy or energy injection and this has become the NuGrid collaboration approach.

The "collapsar" engine for gamma-ray bursts invokes as its energy source the failure of a normal supernova and the formation of a black hole. Here we present the results of the first three-dimensional simulation of the collapse of a massive star down to a black hole, including the subsequent accretion and explosion. The explosion differs significantly from the axisymmetric scenario obtained in two-dimensional simulations; this has important consequences for the nucleosynthetic yields. We compare the nucleosynthetic yields to those of hypernovae. Calculating yields from three-dimensional explosions requires new strategies in post-process nucleosynthesis; we discuss NuGrid's plan for three-dimensional yields.

Simulations of nucleosynthesis in astrophysical environments are at the intersection of nuclear physics reaction rate research and astrophysical applications, for example in the area of galactic chemical evolution or near-field cosmology. Unfortunately, at present the available yields for such applications are based on heterogeneous assumptions between the various contributing nuclear production sites, both in terms of modeling the thermodynamic environment itself as well as the choice of specifc nuclear reaction rates and compilations. On the other side, new nuclear reaction rate determinations are often taking a long time to be included in astrophysical applications. The NuGrid project addresses these issues by providing a set of codes and a framework in which these codes interact. In this contribution we describe the motivation, goals and first results of the NuGrid project. At the core is a new and evolving post-processing nuclesoynthesis code (PPN) that can follow quiescent and explosive nucleosynthesis following multi-zone 1D-stellar evolution as well as multi-zone hydrodynamic input, including explosions. First results are available in the areas of AGB and massive stars.

The s-process production in massive stars at very low metallicities is expected to be negligible due to the low abundance of the neutron source 22Ne, to primary neutron poisons and decreasing iron seed abundances. However, recent models of massive stars including the effects of rotation show that a strong production of 22Ne is possible in the helium core, as a consequence of the primary nitrogen production (observed in halo metal poor stars). Using the PPN post-processing code, we studied the impact of this primary 22Ne on the s process. We find a large production of s elements between strontium and barium, starting with the amount of primary 22Ne predicted by stellar models. There are several key reaction rate uncertainties influencing the s-process efficiency. Among them, 17O(alpha,gamma) may play a crucial role strongly influencing the s process efficiency, or it may play a negligible role, according to the rate used in the calculations. We also report on the development of a new parallel (MPI) post-processing code (MPPNP) designed to follow the complete nucleosynthesis in stars on highly resolved grids. We present here the first post-processing run from the ZAMS up to the end of helium burning for a 15 solar mass model.

We examine observed heavy element abundances in the Cassiopeia A supernova remnant as a constraint on the nature of the Cas A supernova. We compare bulk abundances from 1D and 3D explosion models and spatial distribution of elements in 3D models with those derived from X-ray observations. We also examine the cospatial production of 26Al with other species. We find that the most reliable indicator of the presence of 26Al in unmixed ejecta is a very low S/Si ratio (~0.05). Production of N in O/S/Si-rich regions is also indicative. The biologically important element P is produced at its highest abundance in the same regions. Proxies should be detectable in supernova ejecta with high spatial resolution multiwavelength observations.

Many nucleosynthesis and mixing processes of low-mass stars as they evolve from the Main Sequence to the thermal-pulse Asymptotic Giant Branch phase (TP-AGB) are well understood (although of course important physics components, e.g. rotation, magnetic fields, gravity wave mixing, remain poorly known). Nevertheless, in the last years presolar grain measurements with high resolution have presented new puzzling problems and strong constraints on nucleosynthesis processes in stars. The goal of the NuGrid collaboration is to present uniform yields for a large range of masses and metallicities, including low$-$mass stars and massive stars and their explosions. Here we present the first calculations of stellar evolution and high-resolution, post-processing simulations of an AGB star with an initial mass of 2 M_sun and solar-like metallicity (Z=0.01), based on the post-processing code PPN. In particular, we analyze the formation and evolution of the radiative 13C-pocket between the 17th TP and the 18th TP. The s-process nucleosynthesis profile of a sample of heavy isotopes is also discussed, before the next convective TP occurrence.

The nucleosynthetic yield from a supernova explosion depends upon a variety of effects: progenitor evolution, explosion process, details of the nuclear network, and nuclear rates. Especially in studies of integrated stellar yields, simplifications reduce these uncertainties. But nature is much more complex, and to actually study nuclear rates, we will have to understand the full, complex set of processes involved in nucleosynthesis. Here we discuss a few of these complexities and detail how the NuGrid collaboration will address them.

The currently-favored model for long-duration gamma-ray bursts (GRBs) invokes explosions from the collapse of a massive star down to a black hole: either directly or through fallback. Those GRBs forming via fallback will produce much less radioactive nickel, and hence it has been argued (without any real calculation) that these systems produce dim supernovae. These fallback black-hole GRBs have been recently been argued as possible progenitors of a newly discovered set of GRBs lacking any associated supernovae. Here we present the first ever radiation-hydrodynamics calculations of the light-curves produced in the hypernova explosion by a delayed-fallback gamma-ray burst. We find that the bolometric light-curve is dominated by shock-deposited energy, not the decay of radioactive elements. As such, observations of such bursts actually probe the density in the progenitor wind more than it does the production of radioactive nickel.

Although there is strong support for the collapsar engine as the power source of long-duration gamma-ray bursts (GRBs), we still do not definitively know the progenitor of these explosions. Here we review the current set of progenitor scenarios for long-duration GRBs and the observational constraints on these scenarios. Examining these, we find that single-star models cannot be the only progenitor for long-duration GRBs. Several binary progenitors can match the solid observational constraints and also have the potential to match the trends we are currently seeing in the observations. Type Ib/c supernovae are also likely to be produced primarily in binaries; we discuss the relationship between the progenitors of these explosions and those of the long-duration GRBs.

Recent observations of the Galactic center in high-energy gamma-rays (above 0.1TeV) have opened up new ways to study this region, from understanding the emission source of these high-energy photons to constraining the environment in which they are formed. We present a revised theoretical density model of the inner 5pc surrounding Sgr A* based on the fact that the underlying structure of this region is dominated by the winds from the Wolf-Rayet stars orbiting Sgr A*. An ideal probe and application of this density structure is this high energy gamma-ray emission. We assume a proton-scattering model for the production of these gamma-rays and then determine first whether such a model is consistent with the observations and second whether we can use these observations to further constrain the density distribution in the Galactic center.

The conditions for the leading r-process site candidate, neutrino-driven winds, can not be reproduced self-consistently in current supernova models. For that reason, we investigate an alternate model involving the mass ejected by fallback in a supernova explosion, through hydrodynamic and nucleosynthesis calculations. The nucleosynthetic products of this ejected material produces r-process elements, including those in the vicinity of the elusive 3rd peak at mass number 195. Trans-iron element production beyond the second peak is made possible by a rapid (<1ms) freezeout of alpha particles which leaves behind a large nucleon (including protons!) to r-process seed ratio. This rapid phase is followed by a relatively long (>15 ms) simmering phase at ~2e9 K, which is the thermodynamic consequence of the hydrodynamic trajectory of the turbulent flows in the fallback outburst. During the slow phase high mass elements beyond the second peak are first made through rapid capture of both protons and neutrons. The flow stays close to valley of stability during this phase. After freeze-out of protons the remaining neutrons cause a shift out to short-lived isotopes as is typical for the r-process. A low electron fraction isn't required in this model, however, the detailed final distribution is sensitive to the electron fraction. Our simulations suggest that supernova fallback is a viable alternative scenario for the r-process.

The collapsar engine behind long-duration gamma-ray bursts extracts the energy released from the rapid accretion of a collapsing star onto a stellar-massed black hole. In a collapsing star, this black hole can form in two ways: the direct collapse of the stellar core into a black hole and the delayed collapse of a black hole caused by fallback in a weak supernova explosion. In the case of a delayed-collapse black hole, the strong collapsar-driven explosion overtakes the weak supernova explosion before shock breakout, and it is very difficult to distinguish this black hole formation scenario from the direct collapse scenario. However, the delayed-collapse mechanism, with its double explosion, produces explosive nucleosynthetic yields that are very different from the direct collapse scenario. We present 1-dimensional studies of the nucleosynthetic yields from both black hole formation scenarios, deriving differences and trends in their nucleosynthetic yields.

We present 39 nights of optical photometry, 34 nights of infrared photometry, and 4 nights of optical spectroscopy of the Type Ia SN 1999ac. This supernova was discovered two weeks before maximum light, and observations were begun shortly thereafter. At early times its spectra resembled the unusual SN 1999aa and were characterized by very high velocities in the Ca II H and K lines, but very low velocities in the Si II 6355 A line. The optical photometry showed a slow rise to peak brightness but, quite peculiarly, was followed by a more rapid decline from maximum. Thus, the B- and V-band light curves cannot be characterized by a single stretch factor. We argue that the best measure of the nature of this object is not the decline rate parameter Delta m_15 (B). The B-V colors were unusual from 30 to 90 days after maximum light in that they evolved to bluer values at a much slower rate than normal Type Ia supernovae. The spectra and bolometric light curve indicate that this event was similar to the spectroscopically peculiar slow decliner SN 1999aa.

We compare a suite of 3D explosion calculations and stellar models incorporating advanced physics with observational constraints on the progenitor of Cassiopeia A. We consider binary and single stars from 16 to 40 solar masses with a range of explosion energies and geometries. The parameter space allowed by observations of nitrogen rich high velocity ejecta, ejecta mass, compact remnant mass, and 44Ti and 56Ni abundances individually and as an ensemble is considered. A progenitor of 15-25 solar masses which loses its hydrogen envelope to a binary interaction and undergoes an energetic explosion can match all the observational constraints.

Multi-dimensional simulations of the neutrino-driven mechanism behind core-collapse supernovae have long shown that the explosions from this mechanism would be asymmetric. Recently, detailed core-collapse simulations have shown that the explosion may be strongest in a single direction. We present a suite of simulations modeling these ``single-lobe'' supernova explosions of a 15 solar mass red supergiant star, focusing on the effect these asymmetries have on the gamma-ray emission and the mixing in the explosion. We discuss how asymmetries in the explosion mechanism might explain many of the observed ``asymmetries'' of supernovae, focusing on features of both supernova 1987A and the Cas A supernova remnant. In particular, we show that single-lobe explosions provide a promising solution to the redshifted iron lines of supernova 1987A. We also show that the extent of mixing for explosive burning products depends sensitively on the angular profile of the velocity asymmetry and can be much more extensive than previously assumed.

The possible impact Sgr A East is having on the Galactic center has fueled speculation concerning its age and the energetics of the supernova explosion that produced it. We have carried out the first in-depth analysis of the remnant's evolution and its various interactions: with the stellar winds flowing out from the inner ~2 pc, with the supermassive black hole, Sgr A*, and with the 50 km/s molecular cloud behind and to the East of the nucleus. We have found that a rather "standard" supernova explosion with energy ~1.5e51 ergs is sufficient to create the remnant we see today, and that the latter is probably only ~1,700 years old. The X-ray Ridge between ~9" and 15" to the NE of Sgr A* appears to be the product of the current interaction between the remaining supernova ejecta and the outflowing winds. Perhaps surprisingly, we have also found that the passage of the remnant across the black hole would have enhanced the accretion rate onto the central object by less than a factor 2. Such a small increase cannot explain the current Fe fluorescence observed from the molecular cloud Sgr B2; this fluorescence would have required an increase in Sgr A*'s luminosity by 6 orders of magnitude several hundred years ago. Instead, we have uncovered what appears to be a more plausible scenario for this transient irradiation--the interaction between the expanding remnant and the 50 km/s molecular cloud. The first impact would have occurred about 1,200 years after the explosion, producing a 2-200 keV luminosity of ~1e39 ergs/s. During the intervening 300-400 years, the dissipation of kinetic energy subsided considerably, leading to the much lower luminosity (~1e36 ergs/s at 2-10 keV) we see today.

An explanation for the large random velocities of pulsars is presented. Like many other models, we propose that the momentum imparted to the star is given at birth. The ultimate source of energy is provided by the intense optically thick neutrino flux that is responsible for radiating the proto-neutron star's gravitational binding energy during the Kelvin-Helmholtz phase. The central feature of the kick mechanism is a radiative-driven magnetoacoustic instability, which we refer to as ``neutrino bubbles.'' Identical in nature to the photon bubble instability, the neutrino bubble instability requires the presence of an equilibrium radiative flux as well as a coherent steady background magnetic field. Over regions of large magnetic flux densities, the neutrino bubble instability is allowed to grow on dynamical timescales ~ 1ms, potentially leading to large luminosity enhancements and density fluctuations. Local luminosity enhancements, which preferentially occur over regions of strong magnetic field, lead to a net global asymmetry in the neutrino emission and the young neutron star is propelled in the direction opposite to these regions. For favorable values of magnetic field structure, size, and strength as well as neutrino bubble saturation amplitude, momentum kicks in excess of 1000 km/s can be achieved. Since the neutrino-powered kick is delivered over the duration of the Kelvin-Helmholtz time ~ a few seconds, one expects spin-kick alignment from this neutrino bubble powered model.

We present a comprehensive study of the rho Ophiuchi molecular cloud that addresses aspects of the physical structure and condition of the molecular cloud and its photodissociation region (PDR) by combining far-infrared and submillimeter-wave observations with a wide range of angular scale and resolution. We present 40'x40' maps in CO(4-3) and [CI](3P1-3P0) line emission from the Antarctic Submillimeter Telescope and Remote Observatory (AST/RO), and pointed observations in the CO(7-6) and [CI](3P2-3P1) lines. Within, smaller spectral line maps of the cloud in CO, CS, HCO+ and their rare isotopomers are made at the Heinrich Hertz Submillimeter Telescope Observatory (HHT) in Arizona. Comparison with far-infrared and submillimeter continuum emission, and near-infrared H2 emission allows clearer determination of the physical and chemical structure of the rho Oph photodissociation region (PDR). The excitation conditions needed to produce the observed HCO+ and [OI] emission directly imply inhomogeneous structure. Strong chemical gradients are observed in HCO+ and CS; the former is ascribed to a local enhancement in the H2 ionization rate, the latter is principally due to shocks. The distribution of [CI] is very similar to C18O, and generally consistent with illumination from the 'far' side of the cloud. A notable exception is found at the the western edge of the cloud, where UV photons create a PDR viewed `edge-on'. The abundance of atomic carbon is accurately modeled using a radiation field that decreases with increasing projected distance from the exciting star HD147889. In contrast to conclusions of other studies, we find that no non-equilibrium chemistry is needed to enhance the atomic carbon abundance.

The light curves of Type Ia Supernovae (SN Ia) are powered by gamma-rays emitted by the decay of radioactive elements such as $^{56}$Ni and its decay products. These gamma-rays are downscattered,absorbed, and eventually reprocessed into the optical emission which makes up the bulk of all supernova observations. Detection of the gamma-rays that escape the expanding star provide the only direct means to study this power source for SN Ia light curves. Unfortunately, disagreements between calculations for the gamma-ray lines have made it difficult to interpret any gamma-ray observations. Here we present a detailed comparison of the major gamma-ray line transport codes for a series of 1-dimensional Ia models. Discrepancies in past results were due to errors in the codes, and the corrected versions of the seven different codes yield very similar results. This convergence of the simulation results allows us to infer more reliable information from the current set of gamma-ray observations of SNe Ia. The observations of SNe 1986G, 1991T and 1998bu are consistent with explosion models based on their classification: sub-luminous, super-luminous and normally-luminous respectively.

We present 3-dimensional SPH simulations of supernova explosions from 100 seconds to 1 year after core-bounce. By extending our modelling efforts to a 3-dimensional hydrodynamics treatment, we are able to investigate the effects of explosion asymmetries on mixing and gamma-ray line emergence in supernovae. A series of initial explosion conditions are implemented, including jet-like and equatorial asymmetries of varying degree. For comparison, symmetric explosion models are also calculated. A series of time slices from the explosion evolution are further analyzed using a 3-dimensional Monte Carlo gamma-ray transport code. The emergent hard X- and gamma-ray spectra are calculated as a function of both viewing angle and time, including trends in the gamma-ray line profiles. We find significant differences in the velocity distribution of radioactive nickel between the symmetric and asymmetric explosion models. The effects of this spatial distribution change are reflected in the overall high energy spectrum, as well as in the individual gamma-ray line profiles.