arXiv:1508.03133 Jonas Lippuner Luke Roberts MICRA 2015 Stockholm, - - PowerPoint PPT Presentation

Parameter study of r-process lanthanide production and heating rates in kilonovae with SkyNet

arXiv:1508.03133

Jonas Lippuner Luke Roberts

MICRA 2015 Stockholm, August 17 − 21, 2105

Solar System Abundances

1st peak 2nd peak 3rd peak log10 of relative abundance (Si = 106) Mass number A even A

• dd A

Big Bang Nucleosynthesis Stars (α-elements) iron-peak −2 2 4 6 8 10 25 50 75 100 125 150 175 200 225

Data credit: Katharina Lodders, ApJ 591, 1220 (2003) 2 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

The r-process

τn ≪ τβ− ∼ 10 ms – 10 s Astrophysical sites: supernovae (A ≲ 120) NS–NS/BH merger (up to A ∼ 250)

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

Closed neutron shell

3 Jonas Lippuner

Kilonova

▸ Radioactively powered transient after r-process nucleosynthesis ▸ Triple coincidence: Gravitational wave signal, short GRB, kilonova ▸ Possibly observed after GRB060614 and GRB130603B ▸ Observational signature:

▸ Heating rate → luminosity ▸ Amount of lanthanides and actinides → opacity

4 Jonas Lippuner

Impact of lanthanides and actinides

Lanthanides (between 2nd and 3rd peak) and actinides (beyond 3rd peak) have

• pen f-shells → very high line opacities

lanthanide-free lanthanide-rich

• pacity

∼ 0.1cm2 g−1 ∼ 10cm2 g−1 luminosity bright dim timescale ∼ day ∼ week band blue / optical red / infrared

5 Jonas Lippuner

Parameter study

Parameters 0.01 ≤ Ye ≤ 0.50 initial electron fraction 1kB baryon−1 ≤ s ≤ 100kB baryon−1 initial specific entropy 0.1ms ≤ τ ≤ 500ms expansion time scale Density profile ρ(t,τ) = ⎧ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎩ ρ0e−t/τ t ≤ 3τ ρ0 (3τ te )

3

t ≥ 3τ

t = 3τ Density ρ/ρ0 Time t/τ 10−9 10−6 10−3 1 10−3 10−2 10−1 1 101 102 103

Initial conditions

▸ Choose initial temperature T0 = 6GK ▸ Find ρ0 by solving for NSE at T0 and Ye that produces specified s

6 Jonas Lippuner

SkyNet

▸ General-purpose nuclear reaction network ▸ ∼8000 isotopes, ∼110,000 nuclear reactions ▸ Evolves temperature and entropy based on nuclear reactions ▸ Input: ρ(t), initial composition, initial entropy or temperature ▸ Open source (soon)

7 Jonas Lippuner

SkyNet: Additional Features

Science

▸ Helmholtz equation of state (EOS) ▸ Calculate nuclear statistical equilibrium (NSE) ▸ Calculate inverse rates from detailed balance to be consistent with NSE ▸ NSE evolution mode

Code

▸ Object-oriented C++11 ▸ Python bindings ▸ Uses REACLIB rates, but can easily be extended ▸ Convenient HDF5 output ▸ Make movie with chart of nuclides

8 Jonas Lippuner

Final abundances vs. electron fraction

s = 10kB baryon−1 τ = 7.1ms Relative final abundance Mass number A Ye = 0.01 Ye = 0.13 Ye = 0.25 Lanthanides Solar r-process 10−10 10−9 10−8 10−7 10−6 10−5 10−4 10−3 10−2 50 100 150 200 250

9 Jonas Lippuner

Final abundances vs. entropy

Ye = 0.13 τ = 7.1ms Relative final abundance Mass number A skB = 1 skB = 10 skB = 100 Lanthanides Solar r-process 10−10 10−9 10−8 10−7 10−6 10−5 10−4 10−3 10−2 50 100 150 200 250

10 Jonas Lippuner

Example light curves

s = 10kB baryon−1 τ = 7.1ms M = 0.01M⊙ Luminosity, heating rate [erg s−1] Time [day] Ye = 0.01 Ye = 0.13 Ye = 0.25 Luminosity Heating rate 1039 1040 1041 1042 5 10 15

11 Jonas Lippuner

Light curves vs. electron fraction

s = 10kB baryon−1 τ = 1ms 6days 1day 1600K 6000K log XLa+Ac, tp/3 − 5, Teff/3000 − 4.5 log Peak luminosity Electron fraction Ye Lanthanide and actinide mass fraction XLa+Ac Peak time tp [day] Peak effective temperature Teff [K] Peak Luminosity [erg s−1] −5 −4 −3 −2 −1 0.0 0.1 0.2 0.3 0.4 0.5 1038 1039 1040 1041 1042

12 Jonas Lippuner

Y Y Ye

e e slices

0.1 1 10 100 500 Ye = 0.01 τ [ms] Ye = 0.25 Ye = 0.50 1 10 100 0.1 1 10 100 500 s [kB baryon−1] τ [ms] 1 10 100 s [kB baryon−1] 1 10 100 s [kB baryon−1] 10−5 10−4 10−3 10−2 10−1 0.3 final XLa+Ac 1037 1038 1039 1040 1041 1042 Mǫ at 1 day [erg s−1]

13 Jonas Lippuner

τ τ τ slices

1 10 100 τ = 0.10 ms s [kB baryon−1] τ = 7.1 ms τ = 500 ms 0.01 0.1 0.2 0.3 0.4 0.5 1 10 100 Ye s [kB baryon−1] 0.01 0.1 0.2 0.3 0.4 0.5 Ye 0.01 0.1 0.2 0.3 0.4 0.5 Ye 10−5 10−4 10−3 10−2 10−1 0.3 final XLa+Ac 1037 1038 1039 1040 1041 1042 Mǫ at 1 day [erg s−1]

14 Jonas Lippuner

Summary and conclusions

▸ Lanthanides and actinides have big impact on kilonova light curve ▸ Lanthanide-rich for Ye ≲ 0.22 − 0.30, except:

▸ High s and small τ: neutron-rich freeze-out ▸ Low s and large τ: restarted r-process at high Ye

▸ Uniform heating rate for Ye ≲ 0.4

15 Jonas Lippuner

Additional slides

Additional slides

16 Jonas Lippuner

Lanthanides and actinides

Credit: http://www.sciencegeek.net

17 Jonas Lippuner

Lanthanides and actinides

Credit: http://www.sciencegeek.net

18 Jonas Lippuner

Lanthanides and actinides

Credit: User:DePiep (Wikipedia)

19 Jonas Lippuner

Nucleosynthesis: Overview

▸ Creating nuclides out of nucleons (protons and neutrons) ▸ Cost: overcoming Coulomb barrier (but there is a loophole) ▸ Reward: ∼ 8 MeV per nucleon binding energy ▸ Sources

▸ Big Bang ▸ Stellar fusion ▸ Supernovae ▸ Neutron star mergers (NS–NS and NS–BH)

20 Jonas Lippuner

Big Bang Nucleosynthesis (BBN)

▸ ∼ 100 s after Big Bang, T ∼ 1 GK

H

4He

D/H

3He/H 7Li/H

75.2% 24.8% 3 × 10−5 1 × 10−5 1.5 × 10−10

by mass by mass predicted 5 × 10−10, 4σ discrepancy Credit: NASA / CXC / M. Weiss, note that such big nuclides where not produced in BBN 21 Jonas Lippuner

Stellar Fusion

▸ p–p Coulomb barrier: ∼ 1 MeV ∼ 1010 K ▸ Sun’s core temperature: ∼ 1.6 × 107 K ▸ Produces mostly α-elements: 4He, 12C, 16O, 20Ne, 24Mg, 28Si, 32S, 36Ar, 40Ca, 44Ti∗, 48Cr∗, 52Fe∗, 56Ni∗ (∗ = unstable)

A = 56 ⟨Binding energy⟩/A [MeV] Mass number A average binding energy per nucleon of stable nuclides α-element 1 2 3 4 5 6 7 8 9 10 25 50 75 100 125 150 175 200 225

22 Jonas Lippuner

Stellar Fusion

▸ p–p Coulomb barrier: ∼ 1 MeV ∼ 1010 K ▸ Sun’s core temperature: ∼ 1.6 × 107 K Quantum tunneling! ▸ Produces mostly α-elements: 4He, 12C, 16O, 20Ne, 24Mg, 28Si, 32S, 36Ar, 40Ca, 44Ti∗, 48Cr∗, 52Fe∗, 56Ni∗ (∗ = unstable)

A = 56 ⟨Binding energy⟩/A [MeV] Mass number A average binding energy per nucleon of stable nuclides α-element 1 2 3 4 5 6 7 8 9 10 25 50 75 100 125 150 175 200 225

22 Jonas Lippuner

Supernovae

▸ Thermonuclear (Type Ia)

▸ Thermonuclear explosion of one or two

(?) white dwarf(ves)

▸ Produces mostly iron-peak elements

through nuclear statistical equilibrium (NSE): Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn

▸ Core-collapse (Type II, Ib/c)

▸ Fluorine (ν knocks out a p from 20Ne) ▸ Produces some iron-peak elements ▸ Ejects s-process elements ▸ Maybe r-process

Credit: S. Richers, P. M¨

• sta

23 Jonas Lippuner

Neutron Star Mergers

▸ r-process in neutron-rich ejecta ▸ Ejecta: ∼ 10−3 − 10−2 M⊙

Credit: R. Haas 24 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

s-process Path

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

25 Jonas Lippuner

Double Peaks due to Closed Neutron Shells

neutron drip line

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

26 Jonas Lippuner

Double Peaks due to Closed Neutron Shells

neutron drip line

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

26 Jonas Lippuner

Double Peaks due to Closed Neutron Shells

neutron drip line

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

26 Jonas Lippuner

Double Peaks due to Closed Neutron Shells

neutron drip line

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

26 Jonas Lippuner

Solar System Abundances

r s N = 50 r s N = 82 r s N = 126 log10 of relative abundance (Si = 106) Mass number A even A

• dd A

−2 2 4 6 8 10 25 50 75 100 125 150 175 200 225

Data credit: Katharina Lodders, ApJ 591, 1220 (2003) 27 Jonas Lippuner

s-process vs. r-process

s-process r-process mechanism neutron capture, β− decay No Coulomb barrier! τn 102 − 105 yr ≪ τβ− τβ− ≪ τn 0.01 − 10 s site inside massive stars supernovae? NS-NS/BH mergers? neutron source

13C + 4He → 16O + n 22Ne + 4He → 25Mg + n

neutrino driven wind tidal ejecta of NS material path valley of stability neutron drip line peaks∗ A = 88, 138, 208 strontium, barium, lead A = 80, 130, 194 selenium, xenon, platinum

∗ due to closed neutron shells at N = 50, 82, 126

28 Jonas Lippuner

Shielding

neutron drip line

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

29 Jonas Lippuner

Shielding

s-process neutron drip line

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

29 Jonas Lippuner

Shielding

s-process neutron drip line

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

29 Jonas Lippuner

Shielding

s-process r-process neutron drip line

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

29 Jonas Lippuner

Shielding

s-process only r-process only both s and r neutron drip line

65Cu 66Zn 67Zn 68Zn 70Zn 69Ga 71Ga 70Ge 72Ge 73Ge 74Ge 76Ge 75As 74Se 76Se 77Se 78Se 80Se 82Se 79Br 81Br 78Kr 80Kr 82Kr 83Kr 84Kr 86Kr 85Rb 87Rb 84Sr 86Sr 87Sr 88Sr 89Y 90Zr 91Zr 92Zr

N = 50

29 Jonas Lippuner

SkyNet

▸ General purpose nuclear reaction network ▸ Inputs:

▸ List of nuclides (∼ 8000) ▸ List of nuclear reactions and rates (∼ 100,000) ▸ Initial composition, initial entropy / temperature ▸ Density vs. time

▸ Outputs:

▸ Composition vs. time ▸ Temperature, entropy, heating rate vs. time

30 Jonas Lippuner

SkyNet: Additional Features

Science

▸ Helmholtz equation of state (EOS) ▸ Calculate nuclear statistical equilibrium (NSE) ▸ Calculate inverse rates from detailed balance to be consistent with NSE ▸ NSE evolution mode

Code

▸ Object-oriented C++11 (with CMake build system) ▸ Python bindings (with SWIG) ▸ Support for different matrix solver packages

▸ Sparse: Intel MKL, Trilinos (KLU, UMFPACK, SuperLU), Pardiso ▸ Dense: LAPACK, Trilinos (LAPACK), Armadillo

▸ Convenient HDF5 output ▸ Make movie with chart of nuclides ▸ Open source (soon)

31 Jonas Lippuner

Nuclear Reaction Network

Consider reaction [j] + [k] → [m] (1) cross section = σ = # of reactions per target [j] per second flux of projectiles [k] = R/(Vnj) nkv = r njnkv , (2) and so r = R V = σvnjnk = # of reactions per second per volume, (3) where R = # of reactions per second, V = volume, nj,k = number density of species [j], [k], v = relative speed between [j] and [k].

32 Jonas Lippuner

Nuclear Reaction Network

In general rj,k = ∫ σ(∥v j − v k∥)∥v j − v k∥d3njd3nk, (4) using Boltzmann distribution rj,k = njnk⟨σv⟩j,k = njnk ( 8 µπ )

1/2

(kBT)−3/2 ∫

∞

Eσ(E)e−E/(kB T)dE, (5) where µ = reduced mass = mjmk mj + mk , T = temperature, kB = Boltzmann constant. Note that ⟨σv⟩j,k = ⟨σv⟩j,k(T).

33 Jonas Lippuner

Nuclear Reaction Network

Define abundance Yi = ni nB = # of species [i] # of baryons , (6) where nB is baryon number density, then for [j] + [k] → [m] ˙ Ym = rj,kV # of baryons = rj,k nB = YjnBYknB⟨σv⟩j,k nB = YjYkλj,k, (7) where λj,k = nB⟨σv⟩j,k = NAρ⟨σv⟩j,k(T) = λj,k(T,ρ), (8) where NA is Avogadro’s number, and ρ is the mass density. And, of course ˙ Yj = ˙ Yk = − ˙ Ym. (9)

34 Jonas Lippuner

Nuclear Reaction Network

In general ˙ Yi = ∑

α

Nα

i λα(T,ρ) ∏ m∈Rα

Y

∣Nα

m ∣

m

, (10) where α = index running over all reactions, Nα

i = # of species [i] destroyed/created in α,

Rα = set of reactants of α. Example: ˙ Y4He = ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ decay

4He → 2d

ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ producing reaction p + 7Li → 2 4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ destroying reaction n + p + 2 4He → 7Li + 3He + ⋯ (11)

35 Jonas Lippuner

Nuclear Reaction Network

In general ˙ Yi = ∑

α

Nα

i λα(T,ρ) ∏ m∈Rα

Y

∣Nα

m ∣

m

, (10) where α = index running over all reactions, Nα

i = # of species [i] destroyed/created in α,

Rα = set of reactants of α. Example: ˙ Y4He = − λ4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ decay

4He → 2d

ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ producing reaction p + 7Li → 2 4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ destroying reaction n + p + 2 4He → 7Li + 3He + ⋯ (11)

35 Jonas Lippuner

Nuclear Reaction Network

In general ˙ Yi = ∑

α

Nα

i λα(T,ρ) ∏ m∈Rα

Y

∣Nα

m ∣

m

, (10) where α = index running over all reactions, Nα

i = # of species [i] destroyed/created in α,

Rα = set of reactants of α. Example: ˙ Y4He = − λ4HeY4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ decay

4He → 2d

ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ producing reaction p + 7Li → 2 4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ destroying reaction n + p + 2 4He → 7Li + 3He + ⋯ (11)

35 Jonas Lippuner

Nuclear Reaction Network

In general ˙ Yi = ∑

α

Nα

i λα(T,ρ) ∏ m∈Rα

Y

∣Nα

m ∣

m

, (10) where α = index running over all reactions, Nα

i = # of species [i] destroyed/created in α,

Rα = set of reactants of α. Example: ˙ Y4He = − λ4HeY4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ decay

4He → 2d

+ 2λp,7Li ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ producing reaction p + 7Li → 2 4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ destroying reaction n + p + 2 4He → 7Li + 3He + ⋯ (11)

35 Jonas Lippuner

Nuclear Reaction Network

In general ˙ Yi = ∑

α

Nα

i λα(T,ρ) ∏ m∈Rα

Y

∣Nα

m ∣

m

, (10) where α = index running over all reactions, Nα

i = # of species [i] destroyed/created in α,

Rα = set of reactants of α. Example: ˙ Y4He = − λ4HeY4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ decay

4He → 2d

+ 2λp,7LiYpY7Li ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ producing reaction p + 7Li → 2 4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ destroying reaction n + p + 2 4He → 7Li + 3He + ⋯ (11)

35 Jonas Lippuner

Nuclear Reaction Network

In general ˙ Yi = ∑

α

Nα

i λα(T,ρ) ∏ m∈Rα

Y

∣Nα

m ∣

m

, (10) where α = index running over all reactions, Nα

i = # of species [i] destroyed/created in α,

Rα = set of reactants of α. Example: ˙ Y4He = − λ4HeY4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ decay

4He → 2d

+ 2λp,7LiYpY7Li ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ producing reaction p + 7Li → 2 4He − 2λn,p,2 4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ destroying reaction n + p + 2 4He → 7Li + 3He + ⋯ (11)

35 Jonas Lippuner

Nuclear Reaction Network

In general ˙ Yi = ∑

α

Nα

i λα(T,ρ) ∏ m∈Rα

Y

∣Nα

m ∣

m

, (10) where α = index running over all reactions, Nα

i = # of species [i] destroyed/created in α,

Rα = set of reactants of α. Example: ˙ Y4He = − λ4HeY4He ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ decay

4He → 2d

+ 2λp,7LiYpY7Li ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ producing reaction p + 7Li → 2 4He − 2λn,p,2 4HeYnYpY 2

4He

ÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜÜ destroying reaction n + p + 2 4He → 7Li + 3He + ⋯ (11)

35 Jonas Lippuner

Time Stepping Method

Given Y , T, and ρ, we have a (big) system of coupled ODEs: ˙ Y = ˙ Y (Y ,T,ρ). (12) Use implicit method due to stiffness ˙ Y (t + ∆t) = Y (t + ∆t) − Y (t) ∆t . (13) Use Newton–Raphson to find unknown x = Y (t + ∆t) by finding root of F(x) = x − c ∆t − ˙ Y (x,T,ρ) = 0, (14) where c = Y (t) is a known constant. Update guess with xn+1 = xn − [JF(xn)]−1 F(xn), (15) where the (sparse) Jacobian is (JF)ij = ∂Fi ∂Yj = δij ∆t − ∂ ˙ Yi ∂Yj . (16)

36 Jonas Lippuner

Time Stepping Method

inputs: ρ(t), initial T or s, initial Y calculate ∆t from previous ∆Y calculate reaction rates from T, ρ(t + ∆t) initial guess x0 → Y calculate ˙ Y (xn), J(xn) for guess xn, xn+1 → xn − J−1 ( xn−Y

∆t

− ˙ Y ) converged? n > 10? Y → xn+1, calculate ∆s from ∆Y , s → s + ∆s, T → EOS(s, ρ(t + ∆t), Y ) t → t + ∆t ∆t → ∆t/2 yes no no yes

37 Jonas Lippuner

Nuclear Statistical Equilibrium (NSE)

Let species [i] have Ni neutrons and Zi protons. NSE means the reactions [i] ↔ Ni n + Zi p (17) are in equilibrium so µi = Niµn + Ziµp, (18) where µx are chemical potentials. Abundance is Yi = eµi /(kB T) Gi(T) nB (mikBT 2π̵ h2 )

3/2

, (19) where Gi(T) is internal partition function. Now find µn and µp such that 1 = ∑

i

AiYi and Ye = ∑

i

ZiYi, (20) where Ai = Ni + Zi.

38 Jonas Lippuner

Nuclear Statistical Equilibrium (NSE)

T = 2.5 GK ρ = 1.0 × 107 g cm−3 Ye = 0.50 T = 7.0 GK ρ = 2.2 × 108 g cm−3 Ye = 0.051 T = 6.9 GK ρ = 7.8 × 106 g cm−3 Ye = 0.22

1 10−2 10−4 10−6 10−8 10−10 10−12 Abundance N Z N Z N Z Abundance Mass number A 10−12 10−10 10−8 10−6 10−4 10−2 1 20 40 60 80 100 Mass number A 20 40 60 80 100 Mass number A 20 40 60 80 100

39 Jonas Lippuner

Neutron Star–Black Hole Merger

• 1. Full GR neutron star–black hole

merger simulation (right), up to 10’s of ms

Francois Foucart (LBL), Foucart et al.,

• Phys. Rev. D 90, 024026 (2014)
• 2. Follow ejecta (few × 10−2M⊙) in

SPH simulation, up to 10 s, get many ρ(t) histories,

Matt Duez (WSU)

• 3. Nucleosynthesis with SkyNet,

extend ρ ∝ t−3, 7841 isotopes, 95,467 reactions

JL with Luke Roberts (Caltech) Figure credit: F. Foucart 40 Jonas Lippuner

Neutron Star–Black Hole Merger: Results

Relative final abundance Mass number A SkyNet Ye = 0.051 SkyNet Ye = 0.110 SkyNet Ye = 0.219 Observed solar r-process (scaled) 10−10 10−9 10−8 10−7 10−6 10−5 10−4 10−3 50 100 150 200 250

41 Jonas Lippuner

Parameter Study: Running SkyNet with Python

1

#!/usr/bin/env python

2

from SkyNet import *

3

import numpy as np

4 5

init_T9 = 6.0; Ye = 0.01; s = 10.0; tau_ms = 7.1

6

• pts = NetworkOptions()

7

• pts.ConvergenceCriterion = NetworkConvergenceCriterion.Mass

8

• pts.MassDeviationThreshold = 1.0E-10

9

• pts.IsSelfHeating = True

10 11

nuclib = NuclideLibrary.CreateFromWebnucleoXML(SkyNetRoot + "/data/webnucleo_nuc_v2.0.xml")

12

helm = HelmholtzEOS(SkyNetRoot + "/data/helm_table.dat")

13

weakReacs = REACLIBReactionLibrary(SkyNetRoot + "/data/reaclib",

14

ReactionType.Weak, False, "Weak reactions", nuclib, opts)

15

strongReacs = REACLIBReactionLibrary(SkyNetRoot + "/data/reaclib",

16

ReactionType.Strong, True, "Strong reactions", nuclib, opts)

17

neutFissReacs = REACLIBReactionLibrary(SkyNetRoot + "/data/netsu_panov_symmetric_0neut",

18

ReactionType.Strong, False, "Symmetric neutron induced fission reactis", nuclib, opts)

19

spontFissReacs = REACLIBReactionLibrary(SkyNetRoot + "/data/netsu_sfis_Roberts2010rates",

20

ReactionType.Strong, False, "Spontaneous fission reactions", nuclib, opts)

21 22

net = ReactionNetwork(nuclib, [weakReacs, strongReacs, neutFissReacs, spontFissReacs], helm, opts)

23

nse = NSE.CalcFromTemperatureAndEntropy(init_T9, s, Ye, net.GetNuclideLibrary(), helm)

24

density_profile = ExpTMinus3(nse.Rho(), tau_ms / 1000.0)

25

• utput = net.EvolveSelfHeatingWithInitialTemperature(nse.Y(), 0.0, 1.0E9, init_T9, density_profile,

26

"Ye_%4.2f_s_%07.3f_tau_%07.3f" % (Ye, s, tau_ms)); final_y = output.FinalYVsA()

27

np.savetxt("Y_vs_A", np.column_stack((np.arange(len(final_y)), final_y)), "%6i %30.20E")

42 Jonas Lippuner