Explosive nucleosynthesis of heavy elements: an astrophysical and nuclear physics challenge Gabriel Martínez Pinedo NUSPIN 2017 GSI, Darmstadt, June 26-29, 2017
32 30 28 34 36 38 40 42 46 44 48 26 28 60 30 32 34 36 38 40 42 44 46 48 50 52 54 56 58 50 52 54 56 58 60 62 64 66 68 70 72 70 72 74 76 74 76 78 80 88 90 82 84 98 86 88 90 110 92 112 114 94 116 118 96 98 100 138 156 158 160 162 Making Gold in Nature: r-process nucleosynthesis 10 1 10 0 10 1 10 2 Solar r abundances 10 3 100 150 200 250 62 64 66 68 78 80 82 84 86 N=82 92 94 96 100 102 104 106 108 Known mass Known half life r process waiting point (ETFSI Q) 124 126 120 122 128 130 132 134 136 N=126 r process path 152 154 140 142 144 146 148 150 164 166 168 170 172 174 176 184 180 182 178 186 188 190 N=184 The r-process requires the knowledge of the properties of extremely neutron-rich nuclei: Nuclear masses. Beta-decay half-lives. Neutron capture rates. Fission rates and yields.
r-process astrophysical sites Core-collapse supernova Explosion of massive stars (M 9 M ) Neutrino-winds from protoneutron stars. Strong sensitivity to neutrino interactions at subnuclear densities [GMP+, PRL 109, 251104 (2012)] Only intermediate mass elements are produced (A 100) [GMP+, JPG 41, 044008 (2014)] Neutron star mergers Mergers eject around 0.01 M of very neutron rich-material (Y e 0.01). Similar amount of less neutron-rich matter (Y e 0.2) ejected from accretion disk. Low frequency, high yield: consistent with astronomical observations. Observational signature: electromagnetic transient from radioactive decay of r-process nuclei
Merger channels and ejection mechanism In mergers we deal with a variety of initial configurations (netron-star neutron-star vs neutron-star black-hole) with additional variations in the mass-ratio. The evolution after the merger also allows for further variations. NSNS wind; t~100 ms dyn. ejecta; t~ 1 ms torus unbinding; t ~ 1s BH formation sgrb NSBH Rosswog, et al, Class. Quantum Gravity 34, 104001 (2017)
Evolution nucleosynthesis in mergers r-process stars once electron fermi energy drops below 10 MeV to allow for beta-decays (ρ 10 11 g cm 3 ). Important role of nuclear energy production (mainly beta decay). Energy production increases temperature to values that allow for an (n, γ) (γ, n) equilibrium for most of the trajectories. Systematic uncertainties due to variations of astrophysical conditions and nuclear input 528 trajectories Mendoza-Temis, Wu, Langanke, GMP, Bauswein, Janka, PRC 92, 055805 (2015)
Final abundances different mass models Mendoza-Temis, Wu, Langanke, GMP, Bauswein, Janka, PRC 92, 055805 (2015) abundances at 1 Gyr FRDM WS3 10-3 10-4 10-5 10-6 10-7 abundances at 1 Gyr 10-8 10-3 10-4 10-5 10-6 10-7 HFB21 DZ31 10-8 120 140 160 180 200 220 240 mass number, A 120 140 160 180 200 220 240 mass number, A Robustness astrophysical conditions, strong sensitivity to nuclear physics Second peak (A 130) sensitive to fission yields. Third peak (A 195) sensitive to masses (neutron captures) and beta-decay half-lives. Elements lighter than (A 120) are not produced. Possible contribution of the ejecta from accretion disks.
Temporal evolution (selected phases) Fission is fundamental to determine the final r-process abundances.
The role of N 130 Sn (MeV) 6 5 4 3 2 DZ31 WS3 HFB 21 FRDM 1 0 Yb (Z=70) Isotopes 120 125 130 135 140 Neutron Number Both FRDM and HFB models predict a sudden drop in neutron separation energies approaching N 130 for Z 70 (shape coexistence region).
Global beta-decay calculations Beta-decay rates determine the speed of matter flow from light100 to heavy nuclei. r-process path determined by neutron separation energies P n (%) 50 nuclei with largest impact are those with larger instantaneous half-lives. 0 Despite tremendous progress at RIB facilities Charge Number Z (RIBF at RIKEN) most of the half-lives are based on theoretical calculations. 18 15 FRDM+QRPA DF3+QRPA SM Exp. 24 26 28 30 32 Zhi et al, PRC 87, 024803 (2013) T 1/2 (ms) T1/2 (ms) 1000 100 10 1000 100 10 FRDM+QRPA DF3+QRPA HFB+QRPA SM Exp. 42 44 46 48 Charge number Z Expt. (ISOLDE data) Shell Model Expt. (RIKEN data) 12 9 1 N=82 Isotones 42 43 44 45 46 47 48 49 Charge Number Z 42 44 46 48
Global beta-decay calculations Beta-decay rates determine the speed of matter flow from light to heavy nuclei. r-process path determined by neutron separation energies nuclei with largest impact are those with larger instantaneous half-lives. Despite tremendous progress at RIB facilities (RIBF at RIKEN) most of the half-lives are based on theoretical calculations. Two microscopic calculations (GT+FF) have become available: Covariant density functional theory + QRPA (Marketin+ 2016) Skyrme finite-amplitude method (Mustonen & Engel 2016) Marketin, Huther, GMP, PRC 93, 025805 (2016) proton number 120 110 100 90 80 70 60 50 40 30 20 10-3 -2-1 0 +1 +2 +3 log 10 T 1/2 D3C /T 1/2 F RDM 0 0 20 40 60 80 100 120 140 160 180 200 220 240 neutron number Mustonen & Engel, PRC 93, 014304 (2016)
Impact on r-process abundances Shorter half-lives for Z 80 have a strong impact on the position of A 195 (Eichler+ 2015). abundances at 1 Gyr 10-2 10-3 10-4 10-5 10-6 FRDM masses solar r abundance FRDM+QRPA D3C * abundances at 1 Gyr 10-2 10-3 10-4 10-5 10-6 DZ31 masses solar r abundance FRDM+QRPA D3C * 10-7 120 140 160 180 200 220 240 mass number, A 10-7 120 140 160 180 200 220 240 mass number, A They also affect the robustness of the distribution, the shape of the 2nd peak and the amount of actinides (Wu+, in preparation)
Nucleosynthesis in black-hole accretion disk ejecta abundances at 1 Gyr Accretion disk around compact object is expected to eject material with broad Y e distribution [Fernández, Metzger, MNRAS 435, 502 (2013)] This material is expected to contribute to the production of all r-process nuclides [Wu et al, MNRAS 463, 2323 (2016)] 10-2 10-3 10-4 10-5 10-6 solar r abundance FRDM masses DZ31 masses abundances at 1 Gyr ejecta mass (10-4 M) 10-2 10-3 10-4 10-5 10-6 4 3.5 3 2.5 2 1.5 1 0.5 0 0.1 0.15 0.2 0.25 0.3 0.35 0.4 Y e,5 solar r abundance FRDM+QRPA D3C * 10-7 0 50 100 150 200 250 300 mass number, A 10-7 0 50 100 150 200 250 300 mass number, A
Introduction r process in mergers Summary Kilonova/Macronova electromagnetic transient Electromagnetic transitent from radioactive decay r-process ejecta [Li & Paczyn ski 1998] Luminosities 1000 times those of a nova [Metzger et al, 2010] Large optical opacities of Lanthanides delay the peak to timescales of a week in the red/infrared [Kasen et al, 2013] Probably observed associated to GRB 130603B Time since GRB 130603B (d) 1 0.6 m 21 10 X-ray 10 11 F606W 22 F160W 1.6 m 10 12 24 25 10 13 26 27 N 28 10 14 29 E First direct observation of an r-process event? 104 105 106 Time since GRB 130603B (s) Tanvir+, Nature 500, 457 (2013) X-ray flux (erg s 1 cm 2) 5 AB magnitude 23
Actinides affect opacities and energy production Actinides can be an important opacity source at timescales of weeks. They can substantially contribute to energy production via alpha decay. Mendoza-Temis, Wu, Langanke, GMP, Bauswein, Janka, PRC 92, 055805 (2015)
Impact on the light curve Light curve contains nuclear physics signatures that impact the luminosity up to a factor 10. 10 41 L bol (ergs s 1 ) 10 40 10 39 10 38 10 37 10 36 FRDM (Beta-decay dominates) DZ31 (Alpha-decay dominates) 0 5 10 15 20 25 30 Days Fundamental to determine the amount of material ejected in a merger. Barnes, Kasen, Wu, GMP, ApJ 829, 110 (2016); Rosswog et al, CQG 34, 104001 (2017)
Summary Neutron star mergers most likely constitute the main r process site. The combination of dynamical ejecta and disk outflow ejecta can account for the solar system r-process abundances. Dynamical ejecta of neutron star mergers produce a robust r-process abundance pattern mainly determined by the fission yields of superheavy nuclei. Role of weak interaction on Y e needs to be clarified. Ejecta from black hole accretion disks produce all r-process nuclides in all models considered. Nuclear physics is fundamental for abundance predictions and electromagnetic transient modeling. Transient detection provides a direct confirmation that r process occurs in mergers. Likely observed in GRB 130603B.