Rapid neutron capture process (r-process) elements have been detected in a large fraction of metal-poor halo stars, with abundances relative to iron (Fe) that vary by over two orders of magnitude. This scatter is reduced to less than a factor of 3 in younger Galactic disc stars. The large scatter of r-process elements in the early Galaxy suggests that the r-process is made by rare events, like compact binary mergers and rare sub-classes of supernovae. Although being rare, neutron star mergers alone have difficulties to explain the observed enhancement of r-process elements in the lowest metallicity stars compared to Fe. The supernovae producing the two neutron stars already provide a substantial Fe abundance where the r-process ejecta from the merger would be injected. In this work we investigate another complementary scenario, where the r-process occurs in neutron star-black hole mergers in addition to neutron star mergers. Neutron star-black hole mergers would eject similar amounts of r-process matter as neutron star mergers, but only the neutron star progenitor would have produced Fe. Furthermore, a reduced efficiency of Fe production from single stars significantly alters the age-metallicity relation, which shifts the onset of r-process production to lower metallicities. We use the high-resolution [(20 pc)3/cell] inhomogeneous chemical evolution tool `ICE' to study the outcomes of these effects. In our simulations, an adequate combination of neutron star mergers and neutron star-black hole mergers qualitatively reproduces the observed r-process abundances in the Galaxy.

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.

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

Among extinct radioactivities present in meteorites, 60Fe (t1/2 = 1.49 Myr) plays a key role as a high-resolution chronometer, a heat source in planetesimals, and a fingerprint of the astrophysical setting of solar system formation. A critical issue with 60Fe is that it could have been heterogeneously distributed in the protoplanetary disk, calling into question the efficiency of mixing in the solar nebula or the timing of 60Fe injection relative to planetesimal formation. If this were the case, one would expect meteorites that did not incorporate 60Fe (either because of late injection or incomplete mixing) to show 60Ni deficits (from lack of 60Fe decay) and collateral effects on other neutron-rich isotopes of Fe and Ni (coproduced with 60Fe in core-collapse supernovae and AGB-stars). Here, we show that measured iron meteorites and chondrites have Fe and Ni isotopic compositions identical to Earth. This demonstrates that 60Fe must have been injected into the protosolar nebula and mixed to less than 10 % heterogeneity before formation of planetary bodies.

With presently known input physics and computer simulations in 1D, a self-consistent treatment of core collapse supernovae does not yet lead to successful explosions, while 2D models show some promise. Thus, there are strong indications that the delayed neutrino mechanism works combined with a multi-D convection treatment for unstable layers. On the other hand there is a need to provide correct nucleosynthesis abundances for the progressing field of galactic evolution and observations of low metallicity stars. The innermost ejecta is directly affected by the explosion mechanism, i.e. most strongly the yields of Fe-group nuclei for which an induced piston or thermal bomb treatment will not provide the correct yields because the effect of neutrino interactions is not included. We apply parameterized variations to the neutrino scattering cross sections and alternatively, parameterized variations are applied to the neutrino absorption cross sections on nucleons in the ``gain region''. We find that both measures lead to similar results, causing explosions and a Ye>0.5 in the innermost ejected layers, due to the combined effect of a short weak interaction time scale and a negligible electron degeneracy, unveiling the proton-neutron mass difference. We include all weak interactions (electron and positron capture, beta-decay, neutrino and antineutrino capture on nuclei, and neutrino and antineutrino capture on nucleons) and present first nucleosynthesis results for these innermost ejected layers to discuss how they improve predictions for Fe-group nuclei. The proton-rich environment results in enhanced abundances of 45Sc, 49Ti, and 64Zn as requested by chemical evolution studies and observations of low metallicity stars as well as appreciable production of nuclei in the mass range up to A=80.