LA-UR-18-25559 Draft version June 27, 2018 Typeset using LATEX twocolumn style in AASTeX62 Californium-254 and kilonova light curves Y. Zhu,1 R. T. Wollaeger,2 N. Vassh,3 R. Surman,3, 4 T. M. Sprouse,3 M. R. Mumpower,2, 4, 5 P. Moller¨ ,4, 5 G. C. McLaughlin,1, 4 O. Korobkin,2, 4 T. Kawano,5 P. J. Jaffke,5 E. M. Holmbeck,3, 4 C. L. Fryer,2, 4 W. P. Even,2, 4, 6 A. J. Couture,2, 4 and J. Barnes7, 8 1North Carolina State University, Raleigh, NC 27695-8202, USA 2Center for Theoretical Astrophysics, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA 3University of Notre Dame, Notre Dame, Indiana 46556, USA 4Joint Institute for Nuclear Astrophysics - Center for the Evolution of the Elements, USA 5Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA 6Department of Physical Science, Southern Utah University, Cedar City, UT 84720, USA 7Department of Physics and Columbia Astrophysics Laboratory, Columbia University, New York, NY 10027, USA. 8Einstein Fellow ABSTRACT 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, 254Cf, 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 254Cf 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. Keywords: nuclear reactions, nucleosynthesis, abundances | binaries: close | stars: neutron 1. INTRODUCTION brightness of optical/near-infrared (nIR) counterparts With the first observation of two merging neutron to NSMs { kilonovae (or macronovae, Metzger 2017), stars (NSM) (Abbott et al. 2017) theoretical predictions which are powered by radioactive decays much like su- 56 are for the first time being put to the test with multi- pernovae are powered by the decay of Ni. messenger observational constraints. To understand this Before the light curves of supernovae were known to 56 event requires the combined effort and knowledge from be dominated by the Ni decay chain, it was speculated a broad range of disciplines including nuclear physics, that they might be driven by heavy, long-lived, neutron- 254 atomic physics, astrophysics, and astronomy. rich r-process nuclei, specifically Cf (e.g. Baade et al. Mergers produce extreme conditions which are 1956; Fields et al. 1956). At that time, experimen- unattainable in the laboratory. As a result, observed tal efforts established that spontaneous fission is the arXiv:1806.09724v1 [astro-ph.HE] 25 Jun 2018 signals provide a unique probe of nuclear physics and dominant decay mode of this nucleus with a half-life astrophysics. For example, we can test the idea of of 60:5 ± 0:2 days (Phillips et al. 1963) and an α-decay the production of heavy elements in mergers (Lattimer branching of 0:31 ± 0:016% (Bemis & Halperin 1968). & Schramm 1974) by comparing nucleosynthetic mod- In contrast to supernovae, kilonovae from NSMs are els to observed light curves (Kasen et al. 2017; Tanvir thought to be powered by residual radioactivity of the et al. 2017). Nuclear heating rates directly impact the r-process, and have a lanthanide-rich component. If a complete main r-process pattern is produced in mergers, actinides must be produced as well. In such case, fission Corresponding author: M. R. Mumpower may contribute substantially to nuclear heating due to [email protected] its ∼ 200 MeV energy release. Fissioning nuclei that are 2 Y. Zhu et al. likely to influence the light curve are those with half- & Takahashi 1975; Petermann et al. 2012). Neverthe- lives on the order of days, which roughly corresponds less, despite recent advances (Goriely et al. 2009; Giu- to the predicted peak timescale (Li & Paczy´nski 1998; liani et al. 2018), the treatment of spontaneous fission in Metzger 2017). the r-process is subject to large uncertainties, and it is While the potential for the late-time dominance of impossible to rule out unmeasured nuclei that are pop- 254Cfhas been noted in previous heating calculations ulated during the r-process and have decay timescales (Wanajo et al. 2014), the effect of this experimentally on the order of days. established spontaneous fission process on the late-time Turning to experimentally measured data, some of the light curve has not yet been explored. In this work we nuclei which could potentially become significantly pop- report that when experimentally known spontaneous fis- ulated in the r-process and undergo spontaneous fission sion decays are included in nucleosynthesis calculations include isotopes of Californium, Fermium, as well as for neutron-rich NSM ejecta, 254Cf and its fission daugh- 257Es and 260Md, see Fig.1. Although many isotopes of ter products are dominant contributors to the nuclear Fermium reportedly undergo sf, their experimentally es- heating at ∼ 15 − 250 days, greatly impacting late-time tablished half-lives are on the order of seconds or faster, light curves. thus they are not populated for long enough to make significant contributions to the heating. For example, 258 2. NUCLEOSYNTHESIS AND FISSION the nucleus Fm fissions very quickly; NUBASE2016 re- ports ∼ 100% sf branching with a half-life of 370 ± 14 To model r-process nucleosynthesis we use the µs. The same applies to 256Cf given its 12:3 m half-life. Portable Routines for Integrated nucleoSynthesis Mod- This leaves 257Es, 260Md, and 254Cf in this populated eling (PRISM) reaction network developed jointly at region with half-lives experimentally known to be on the the University of Notre Dame and Los Alamos National order of days. Laboratory. In this network, now at version 2:0, the We used NUBASE2016 branchings and so assume 257Es reheating of the ejecta is handled self-consistently. We to undergo only β-decay, however the possibility for this choose wind conditions that are consistent with current nucleus to decay via sf has not been experimentally ruled neutron star merger ejecta models: entropy per baryon out. Mendelevium-260 is known to dominantly undergo s=k = 40, outflow timescale τ = 20, and Y = 0:2. e sf with half-life ∼ 30 days, however the population of Our nucleosynthesis calculations contain all relevant 260Md is highly subject to the theoretical decay rates nuclear reaction channels including charged particle re- applied to its potential β-feeder 260Fm, which based on actions, neutron capture, photodissociation, β-decay, 254 FRDM2012 masses has Q − < 0. This leaves Cf as and delayed neutron emission. Fission from neutron- β a nucleus likely to influence the heating rate and light induced, β-delayed and spontaneous channels are also curve. included. The nuclear properties are based on the the- The extent of the influence of 254Cf depends on the oretical nuclear model FRDM2012 (M¨olleret al. 2015, mechanisms by which the nucleus becomes populated. 2016; Mumpower et al. 2016, 2018; M¨olleret al. 2018). Figure1 shows that we find 254Cf to be populated solely When available, evaluated data is used for masses (Wang by the β-decay of nuclei in the A = 254 isobaric chain. et al. 2017) and decay properties from NUBASE2016 As can be seen to the right of the black line in Fig.1, (Audi et al. 2017). The database contains NUBASE2016 the nuclei which β-decay to populate 254Cf are for the half-lives and branching ratios for the decay channels of most part unstudied. Alpha-feeding from 258Fm is pro- β±, electron capture, α-decay, spontaneous fission (sf), hibited by the uncertain NUBASE2016 branching data. A and the spontaneous emission of neutrons and protons. non-zero α-branching from 258Fm would seek to amplify For the inclusion of all such evaluated decay rates and the influence of 254Cf possibly by opening pathways to branching ratios, care is taken to ensure no theoretical population via α-decay chains of heavier nuclei. This decay rates interfere with experimentally established de- highlights the importance of future experimental inves- cay processes. tigations to understand the precise branchings of nuclei The choice of theoretical sf rates produce small vari- in this region. ation in the computed heating rates. In the calcula- We construct fission fragment yields of 254Cf with a tions we present we used a parameterization found to hybrid method that combines both theoretical and ex- fit within ∼ 1 − 3 orders of magnitude of the measured perimental data. For 254Cf(sf), fission fragment yields half-lives of elements with Z < 100 (Xu & Ren 2005). Y (A; Z) in both mass A and charge Z are used in order Similar heating rates are obtained using phenomenolog- to produce the most accurate estimate of the energy re- ical dependences on fission barrier height, such as (Za- lease. Our hybrid method is based on experimental data grebaev et al. 2011; Karpov et al. 2012) and (Kodama Californium-254 and kilonova light curves 3 6 4 Md 10− 254 X 254 258 Cf ] Fm Fm 10 8 6 3 − − Es 10 10− Z 2 254 Cf Cf 12 10− Abundance Bk Abundance [10 1 14 10− Cm 0 10 16 154 156 158 160 − 0.001 0.01 0.1 1 10 100 N t [days] Figure 1. Population of 254Cf. Left: abundances of nuclei at t = 0:02 days. The β-decay path into 254Cf and the potential α-decay from 258Fm are shown.