The light curve of SN 1987A revisited: constraining production masses of radioactive nuclides

We revisit the evidence for the contribution of the long-lived radioactive nuclides 44Ti, 55Fe, 56Co, 57Co, and 60Co to the UVOIR light curve of SN 1987A. We show that the V-band luminosity constitutes a roughly constant fraction of the bolometric luminosity between 900 and 1900 days, and we obtain an approximate bolometric light curve out to 4334 days by scaling the late time V-band data by a constant factor where no bolometric light curve data is available. Considering the five most relevant decay chains starting at 44Ti, 55Co, 56Ni, 57Ni, and 60Co, we perform a least squares fit to the constructed composite bolometric light curve. For the nickel isotopes, we obtain best fit values of M(56Ni) = (7.1 +- 0.3) x 10^{-2} Msun and M(57Ni) = (4.1 +- 1.8) x 10^{-3} Msun. Our best fit 44Ti mass is M(44Ti) = (0.55 +- 0.17) x 10^{-4} Msun, which is in disagreement with the much higher (3.1 +- 0.8) x 10^{-4} Msun recently derived from INTEGRAL observations. The associated uncertainties far exceed the best fit values for 55Co and 60Co and, as a result, we only give upper limits on the production masses of M(55Co) < 7.2 x 10^{-3} Msun and M(60Co) < 1.7 x 10^{-4} Msun. Furthermore, we find that the leptonic channels in the decay of 57Co (internal conversion and Auger electrons) are a significant contribution and constitute up to 15.5% of the total luminosity. Consideration of the kinetic energy of these electrons is essential in lowering our best fit nickel isotope production ratio to [57Ni/56Ni]=2.5+-1.1, which is still somewhat high but is in agreement with gamma-ray observations and model predictions.Comment: 7 pages, 6 pages, 2 table

Constraints on explosive silicon burning in core-collapse supernovae from measured Ni/Fe ratios

Measurements of explosive nucleosynthesis yields in core-collapse supernovae provide tests for explosion models. We investigate constraints on explosive conditions derivable from measured amounts of nickel and iron after radioactive decays using nucleosynthesis networks with parameterized thermodynamic trajectories. The Ni/Fe ratio is for most regimes dominated by the production ratio of 58Ni/(54Fe + 56Ni), which tends to grow with higher neutron excess and with higher entropy. For SN 2012ec, a supernova that produced a Ni/Fe ratio of $3.4\pm1.2$ times solar, we find that burning of a fuel with neutron excess $\eta \approx 6\times 10^{-3}$ is required. Unless the progenitor metallicity is over 5 times solar, the only layer in the progenitor with such a neutron excess is the silicon shell. Supernovae producing large amounts of stable nickel thus suggest that this deep-lying layer can be, at least partially, ejected in the explosion. We find that common spherically symmetric models of $M_{\rm ZAMS} \lesssim 13$ Msun stars exploding with a delay time of less than one second ($M_{\rm cut} < 1.5$ Msun) are able to achieve such silicon-shell ejection. Supernovae that produce solar or sub-solar Ni/Fe ratios, such as SN 1987A, must instead have burnt and ejected only oxygen-shell material, which allows a lower limit to the mass cut to be set. Finally, we find that the extreme Ni/Fe value of 60-75 times solar derived for the Crab cannot be reproduced by any realistic-entropy burning outside the iron core, and neutrino-neutronization obtained in electron-capture models remains the only viable explanation.Comment: 13 pages, 9 figures, accepted for publication in Ap

Proton-Rich Nuclear Statistical Equilibrium

Proton-rich material in a state of nuclear statistical equilibrium (NSE) is one of the least studied regimes of nucleosynthesis. One reason for this is that after hydrogen burning, stellar evolution proceeds at conditions of equal number of neutrons and protons or at a slight degree of neutron-richness. Proton-rich nucleosynthesis in stars tends to occur only when hydrogen-rich material that accretes onto a white dwarf or neutron star explodes, or when neutrino interactions in the winds from a nascent proto-neutron star or collapsar-disk drive the matter proton-rich prior to or during the nucleosynthesis. In this paper we solve the NSE equations for a range of proton-rich thermodynamic conditions. We show that cold proton-rich NSE is qualitatively different from neutron-rich NSE. Instead of being dominated by the Fe-peak nuclei with the largest binding energy per nucleon that have a proton to nucleon ratio close to the prescribed electron fraction, NSE for proton-rich material near freeze-out temperature is mainly composed of Ni56 and free protons. Previous results of nuclear reaction network calculations rely on this non-intuitive high proton abundance, which this paper will explain. We show how the differences and especially the large fraction of free protons arises from the minimization of the free energy as a result of a delicate competition between the entropy and the nuclear binding energy.Comment: 4 pages, 7 figure

The effect of 12C + 12C rate uncertainties on s-process yields

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.Comment: Conference proceedings for the Nuclear Physics in Astrophysics IV conference, 8-12 June 2009. 4 pages, 3 figures. Accepted for publication to the Journal of Physics: Conference Serie

NuGrid stellar data set. 1. Stellar yields from H to Bi for stars with metallicities Z=0.02 and Z=0.01

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 database for applications in areas such as presolar grain studies. Our non-rotating models assume convective boundary mixing (CBM) 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 simple analytic core-collapse supernova (SN) models with two options for fallback and shock velocities. The explosions show which pre-SN yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impact 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 CBM at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the ¹³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. The entirety of our stellar nucleosynthesis profile and time evolution output are available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced

Reaction rate sensitivity of 44Ti production in massive stars and implications of a thick target yield measurement of 40Ca(alpha,gamma)44Ti

We evaluate two dominant nuclear reaction rates and their uncertainties that affect 44Ti production in explosive nucleosynthesis. Experimentally we develop thick-target yields for the 40Ca(alpha,gamma)44Ti reaction at E(alpha) = 4.13, 4.54, and 5.36 MeV using gamma-ray spectroscopy. At the highest beam energy, we also performed an activation measurement that agrees with the thick target result. From the measured yields a stellar reaction rate was developed that is smaller than current statistical-model calculations and recent experimental results, which would suggest lower 44Ti production in scenarios for the alpha-rich freeze out. Special attention has been paid to assessing realistic uncertainties of stellar rates produced from a combination of experimental and theoretical cross sections, which we use to develop a re-evaluation of the 44Ti(alpha,p)47V reaction rate. Using these we carry out a sensitivity survey of 44Ti synthesis in eight expansions representing peak temperature and density conditions drawn from a suite of recent supernova explosion models. Our results suggest that the current uncertainty in these two reaction rates could lead to as large an uncertainty in 44Ti synthesis as that produced by different treatments of stellar physics.Comment: Comments: 45 pages, 19 postscript figures Minor corrections from Referee and Proof Editors Figs 9 & 10 now in colo

Revisiting INTEGRAL/SPI observations of 44^44Ti from Cassiopeia A

T. Siegert, R. Diehl, M. G. H. Krause, and J. Greiner, “Revisiting INTEGRAL/SPI observations of 44Ti from Cassiopeia A”, Astronomy & Astrophysics, Vol. 579, July 2015. This version of record is available online at: https://doi.org/10.1051/0004-6361/201525877 Reproduced with Permission from Astronomy and Astrophysics, © ESO 2015.Context. The 340-yr old supernova remnant Cassiopeia A, located at 3.4 kpc distance, is the best-studied young core-collapse super- nova remnant. Nucleosynthesis yields in radioactive isotopes have been studied with different methods, in particular for production and ejection of 44Ti and 56Ni, which originate in the innermost regions of the supernova. 44Ti was first discovered in this remnant, but is not seen consistently in other core-collapse sources. Aims. We aim to measure radioactive 44Ti ejected in Cassiopeia A and to place constraints on velocities of these ejecta by determining X- and γ-ray line-shape parameters of the emission lines. Methods. We analyzed the observations made with the SPI spectrometer on INTEGRAL together with an improved instrumental background method, to achieve a high spectroscopic resolution that enables interpretation for a velocity constraint on 44Ti ejecta from the 1.157 MeV γ-ray line of the 44Sc decay. Results. We observe both the hard X-ray line at 78 keV and the γ-ray line at 1157 keV from the 44Ti decay chain at a combined significance of 3.8σ. Measured fluxes are (2.1 ± 0.4) × 10−5 ph cm−2 s−1 and (3.5 ± 1.2) × 10−5 ph cm−2 s−1, which corresponds to (1.5 ± 0.4) × 10−4 and (2.4 ± 0.9) × 10−4 M⊙ of 44Ti, respectively. The measured Doppler broadening of the lines implies expansion velocities of 4300 and 2200 km s−1 , respectively. By combining our results with previous studies, we determine a more precise esti- mate of ejected 44Ti of (1.37 ± 0.19) × 10−4 M⊙. Conclusions. The measurements of the two lines are consistent with previous studies. The flux in the line originating from excited 44Ca is significantly higher than the flux determined in the lines from 44Sc. Cosmic-ray acceleration within the supernova remnant may be responsible for an additional contribution to this line from nuclear de-excitation following energetic particle collisions in the remnant and swept-up material.Peer reviewe

Asymmetries in core-collapse supernovae from maps of radioactive 44Ti in CassiopeiaA

Asymmetry is required by most numerical simulations of stellar core-collapse explosions, but the form it takes differs significantly among models. The spatial distribution of radioactive 44Ti, synthesized in an exploding star near the boundary between material falling back onto the collapsing core and that ejected into the surrounding medium1, directly probes the explosion asymmetries. Cassiopeia A is a young2, nearby3, core-collapse4 remnant from which 44Ti emission has previously been detected5, 6, 7, 8 but not imaged. Asymmetries in the explosion have been indirectly inferred from a high ratio of observed 44Ti emission to estimated 56Ni emission9, from optical light echoes10, and from jet-like features seen in the X-ray11 and optical12 ejecta. Here we report spatial maps and spectral properties of the 44Ti in Cassiopeia A. This may explain the unexpected lack of correlation between the 44Ti and iron X-ray emission, the latter being visible only in shock-heated material. The observed spatial distribution rules out symmetric explosions even with a high level of convective mixing, as well as highly asymmetric bipolar explosions resulting from a fast-rotating progenitor. Instead, these observations provide strong evidence for the development of low-mode convective instabilities in core-collapse supernovae

Nucleosynthesis simulations for a wide range of nuclear production sites from NuGrid

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