Electromagnetic counterparts and r-process Tsvi Piran The Hebrew - - PowerPoint PPT Presentation

electromagnetic counterparts and r process
SMART_READER_LITE
LIVE PREVIEW

Electromagnetic counterparts and r-process Tsvi Piran The Hebrew - - PowerPoint PPT Presentation

Electromagnetic counterparts and r-process Tsvi Piran The Hebrew University Kenta Hotokezaka, Ehud Nakar Kyoto - Nov 2016 Outline 1. The Li-Paczynski Macronova (kilonova) 2. GRBs 060614/050709 and their Macronove 3. Plutonium 4. Dwarf


slide-1
SLIDE 1

Electromagnetic counterparts and r-process

Tsvi Piran

The Hebrew University Kenta Hotokezaka, Ehud Nakar

Kyoto - Nov 2016

slide-2
SLIDE 2

Outline

1. The Li-Paczynski Macronova (kilonova) 2. GRBs 060614/050709 and their Macronove 3. Plutonium 4. Dwarf Galaxies 5. The cocoon’ s macronova - the strongest EM counterpart? 6. Limits on magnetars from radio flares 7 . * The energy deposition rate 8. Conclusions

slide-3
SLIDE 3
  • Radioactive decay of the

neutron rich matter.

  • Eradioactive ≈ 0.001 Mc2 ≈

1050 erg

  • A weak short Supernova

like event.

  • 1. Macronova*(Li & Paczynski 1997)

Bohdan Paczynski

*Also called Kilonova

slide-4
SLIDE 4
  • Radioactive decay of the

neutron rich matter.

  • Eradioactive ≈ 0.001 Mc2 ≈

1050 erg

  • A weak short Supernova

like event.

  • 1. Macronova*(Li & Paczynski 1997)

Bohdan Paczynski

*Also called Kilonova Hektanova

slide-5
SLIDE 5
  • Radioactive decay of the

neutron rich matter.

  • Eradioactive ≈ 0.001 Mc2 ≈

1050 erg

  • A weak short Supernova

like event.

  • 1. Macronova*(Li & Paczynski 1997)

Bohdan Paczynski

*Also called Kilonova Hektanova Decanova

slide-6
SLIDE 6
  • After a second dE/dt∝t-1.3 (Freiburghaus+

1999; Korobkin + 2013)

Radioactive Decay*

Korobkin + 13; Rosswog, Korobkin + 13

slide-7
SLIDE 7

Photons escape from this region luminosity Decrease due to radioactive decay Increase as we see a large fraction of the matter. time The light curve depends on 1. mass 2. velocity 3.

  • pacity
slide-8
SLIDE 8

Photons escape from this region luminosity Decrease due to radioactive decay Increase as we see a large fraction of the matter. time The light curve depends on 1. mass 2. velocity 3.

  • pacity
slide-9
SLIDE 9
  • S. Rosswog, … Following Davies + 1994
slide-10
SLIDE 10

Lanthanides dominate the opacity

(Kassen & Barnes 13, Tanaka & Hotokezaka 13) )

κ= 10cm2/gm tmax ∝κ1/2 => l o n g e r Lmax ∝κ-0.65 => weaker

T ∝ κ-0.4 => redder

slide-11
SLIDE 11

Lanthanides dominate the opacity

(Kassen & Barnes 13, Tanaka & Hotokezaka 13) )

κ= 10cm2/gm tmax ∝κ1/2 => l o n g e r Lmax ∝κ-0.65 => weaker

T ∝ κ-0.4 => redder

slide-12
SLIDE 12

Lanthanides dominate the opacity

(Kassen & Barnes 13, Tanaka & Hotokezaka 13) )

κ= 10cm2/gm tmax ∝κ1/2 => l o n g e r Lmax ∝κ-0.65 => weaker

T ∝ κ-0.4 => redder

1 days 10

slide-13
SLIDE 13

Lanthanides dominate the opacity

(Kassen & Barnes 13, Tanaka & Hotokezaka 13) )

κ= 10cm2/gm tmax ∝κ1/2 => l o n g e r Lmax ∝κ-0.65 => weaker

T ∝ κ-0.4 => redder

1 days 10 1040 1041

slide-14
SLIDE 14

Lanthanides dominate the opacity

(Kassen & Barnes 13, Tanaka & Hotokezaka 13) )

κ= 10cm2/gm tmax ∝κ1/2 => l o n g e r Lmax ∝κ-0.65 => weaker

T ∝ κ-0.4 => redder

1 days 10 1040 1041

uv or optical -> IR

slide-15
SLIDE 15

GRB130603B @ 9 days AB

(6.6 days at the source frame)

nIR HST image (Tanvir + 13) V

slide-16
SLIDE 16

Swift Tanvir + 13, Berger + 13

Macronova?

slide-17
SLIDE 17

If correct

Confirmaiton of the GRB neutron star merger model (Eichler, Livio, TP & Schramm 1989). Confirmation of the Li-Paczynski Macronova (Li-Paczynski 1997). Confirmation that compact binary mergers are the source of heavy (A>130) r-process material: Gold, Silver, Platinum, Plotonium, Uranium etc…(Lattimer & Schramm, 75).

slide-18
SLIDE 18

The rate of Short GRBs Macronova and r- process

About 1/3 of Swift short (<2sec) GRBs are Collapsars The rate of non-Collapsar short GRBs (sGRbs) is 4.1+2.3-1.9 Gpc-3 yr-1 (depending on the assumed minimal luminosity). A LIGO detection rate of 3-100 per year (0.1-3 coinciding with a sGRB)* A typical time delay of ~3 Gyr after SFR=> an initial separation of ~2 x 1011 cm But selection effects? Maybe consistent with p(τ)~1/τ With beaming of ~30 and mass ejection of 0.02 Msun - compatible with R-process nucleosynthesis for A>110 elements.

slide-19
SLIDE 19

GRB 060614

Yang et al., 2015 Need M≃0.1M⨀ => BH-NS ?

slide-20
SLIDE 20

GRB 050709

Jin et al., 2016 Need M≃0.05M⨀ => BH-NS ?

!! !" !# !$ !% !& !' ! " $ !( "( #( ) )( !"# !"#$ %$#&'()*" *+,- .+/0- 123.4 5678.9 :')#; < = > !! !" !# !$ !% !& !' ! " $ !( "( #( ) )(

!"! !"# !"$ !"% !"& '"! '"# # ( ) '! !"#$ %&'()*+ ,-.+/ *+,-.,/01 2'!!"345 !"! !"# !"$ !"% !"& '"! '"# # ( ) '!

slide-21
SLIDE 21

Are Macronova Frequent?

There are 3 (6) possible (nearby) historical candidates with a good enough data In 3/3 (3/ 6) there are possible Macronovae

slide-22
SLIDE 22

R-Process

0.1 1 10 100 1000 10000 0.0001 0.001 0.01 0.1 1 R(z=0) [Myr-1] Mej [Msun]

Short GRBs Galactic NS2 LIGO/Virgo limit Advanced (5yr) R-element mass (A>90) Macronova candidate

slide-23
SLIDE 23

R-Process

0.1 1 10 100 1000 10000 0.0001 0.001 0.01 0.1 1 R(z=0) [Myr-1] Mej [Msun]

Short GRBs Galactic NS2 LIGO/Virgo limit Advanced (5yr) R-element mass (A>90) Macronova candidate

*

slide-24
SLIDE 24

0.1 1 10 100 1000 10000 0.0001 0.001 0.01 0.1 1 R(z=0) [Myr-1] Mej [Msun]

Short GRBs Galactic NS2 LIGO/Virgo limit Advanced (5yr) R-element mass (A>90) Macronova candidate

Wanderman Piran 15

*

slide-25
SLIDE 25

0.1 1 10 100 1000 10000 0.0001 0.001 0.01 0.1 1 R(z=0) [Myr-1] Mej [Msun]

Short GRBs Galactic NS2 LIGO/Virgo limit Advanced (5yr) R-element mass (A>90) Macronova candidate

Wanderman Piran 15

*

slide-26
SLIDE 26

0.1 1 10 100 1000 10000 0.0001 0.001 0.01 0.1 1 R(z=0) [Myr-1] Mej [Msun]

Short GRBs Galactic NS2 LIGO/Virgo limit Advanced (5yr) R-element mass (A>90) Macronova candidate

Can we break the yield - rate degeneracy?

slide-27
SLIDE 27

0.1 1 10 100 1000 10000 0.0001 0.001 0.01 0.1 1 R(z=0) [Myr-1] Mej [Msun]

Short GRBs Galactic NS2 LIGO/Virgo limit Advanced (5yr) R-element mass (A>90) Macronova candidate

*

Can we break the yield - rate degeneracy?

slide-28
SLIDE 28

t t Rare Events Frequent events

Radioactive Elements

slide-29
SLIDE 29

High 244Pu at the early solar system =>

244Pu Radioactive decay time ~ 100 Myear

A nearby event near solar system Mixing time < 150 Myr Large fluctuations possible => Event rate is low Lack of Cu => 10 Myr < Mixing length

slide-30
SLIDE 30

Tissot + 16t

slide-31
SLIDE 31

The early solar system

244Pu (half life 81Myr)

Wallner + 14

slide-32
SLIDE 32

Rare and “massive” events

slide-33
SLIDE 33

r-process material in Dwarf Galaxies

(Beniamini+ 16a,b)

slide-34
SLIDE 34
slide-35
SLIDE 35
slide-36
SLIDE 36

The Secret Signatures of GRB cocoons

Nakar & TP ApJ 16 in press

From Mizuta

slide-37
SLIDE 37

The idea in a single picture

slide-38
SLIDE 38

The Jet drills a hole in the star

Zhang, Woosley & MacFadyen 2004

slide-39
SLIDE 39

Jet breakout

(Bromberg Nakar, TP, Sari 11 ApJ 2011)

The engine must be active until the jet’ s head breaks out!*

slide-40
SLIDE 40

A prediction of the Collapsar model

T90 = Te-TB

Observed duration Engine time Break out time

slide-41
SLIDE 41

A prediction of the Collapsar model

T90 = Te-TB

T90

Observed duration Engine time Break out time

dN(T90)/dt

slide-42
SLIDE 42

A prediction of the Collapsar model

TB

T90 = Te-TB

T90

Observed duration Engine time Break out time

dN(T90)/dt

slide-43
SLIDE 43

Short Long

?

T90

slide-44
SLIDE 44

Short Long

?

dlog(N)/dT90 T90

slide-45
SLIDE 45

A second look

(Bromberg Nakar, TP & Sari, 2011) T90 d N / d T

9

slide-46
SLIDE 46

A second look

(Bromberg Nakar, TP & Sari, 2011)

A direct observational proof of the Collapsar model.

T90 d N / d T

9

slide-47
SLIDE 47

Short (Non-Collapsars)

Collapsars

slide-48
SLIDE 48

Short (Non-Collapsars)

Collapsars

slide-49
SLIDE 49
slide-50
SLIDE 50
slide-51
SLIDE 51

Swift Short (Non- Collapsars) GRBs

Collapsars

slide-52
SLIDE 52

Swift Short (Non- Collapsars) GRBs

Collapsars

slide-53
SLIDE 53

Swift Short (Non- Collapsars) GRBs

Collapsars Short Swift GRBs with T90>0.7sec are not “short”!

slide-54
SLIDE 54

EGRB≈Eejecta≈Ec

Macronova + Radio flare

slide-55
SLIDE 55

Cocoon’ s structure

slide-56
SLIDE 56

3D simulation

4Msun, R*=4x1010cm. Lj =1051erg/s, θ=8ο Using Pluto with high resolution ΔR=107cm. Credit: Ore Gottlieb

slide-57
SLIDE 57

3D simulation

4Msun, R*=4x1010cm. Lj =1051erg/s, θ=8ο Using Pluto with high resolution ΔR=107cm. Credit: Ore Gottlieb

slide-58
SLIDE 58

2D simulation 110sec after breakout

4Msun, R*=4x1010cm. Lj =1051erg/s, θ=8ο Using Pluto with high resolution ΔR=107cm. Credit: Ore Gottlieb

Jet Star Wide angle Γ≈10 material

slide-59
SLIDE 59

The cocoons

Harrison, Goetlieb and Nakar in prep, 2016

slide-60
SLIDE 60

Emission component

Newtonian Cocoon - cooling (photospheric) emission Newtonian cocoon - macronova Relativistic Jet cocoon - cooling (photospheric) emission Relativistic Jet cocoon - afterglow

slide-61
SLIDE 61

The cocoons

Harrison, Goetlieb & Nakar in prep, 2016

Light “relativistic” Jet cocoon Heavy “Newtonian” stellar cocoon Jet

slide-62
SLIDE 62

Cocoon Dynamics

R Rθ

Stellar Envelope

L=Ecc/R

slide-63
SLIDE 63

η3 η2 η1

Full mixing

η=1

Partial mixing No mixing

slide-64
SLIDE 64

Partial Mixing

Harrison, Goetlieb and Nakar in prep, 2016

Light “relativistic” Jet cocoon Heavy “Newtonian” Stellar cocoon

slide-65
SLIDE 65

2D simulation 110sec after breakout

Jet Star Wide angle Γ≈10 material

4MO, R*=4x1010cm. Lj =1051erg/s, θ=8ο Using Pluto with high resolution ΔR=107cm. Credit: Ore Gottlieb

.

slide-66
SLIDE 66

Γβ Energy per Log interval Newtonian Relativistic

slide-67
SLIDE 67

Short GRBs

Nagakura et al. 2014; Murguia-Berthier et al. 2014, 2016

From Hotokezaka & TP 2015

slide-68
SLIDE 68
slide-69
SLIDE 69

SGRB cocoon signatures

+ =>

  • Rel. Cocoon cooling
  • Rel. Cocoon Afterglow,

scaling from the regular SGRB afterglow This is a wide angle signal 0.5 rad is stronger than typical SGRB orphan afterglow

slide-70
SLIDE 70

Macronova cocoon signature

Heating due to radioactive decay Blue signal at around 0.5-1 day! Brighter or comparable to the classical Macronova

slide-71
SLIDE 71

Summary

Cocoons are the forgotten cousins in the GRB story. They carry a comparable amount of energy to the GRB and are wider than the GRBs. Short GRBs have their own cocoons whose signatures might be the best EM counterpart to

slide-72
SLIDE 72

The radio - flare (Nakar & Piran 2011) Testing the Macronova interpretation

A long lasting radio flare due to the interaction of the ejecta with surrounding matter may follow the macronova.

slide-73
SLIDE 73

The radio - flare (Nakar & Piran 2011) Testing the Macronova interpretation

A long lasting radio flare due to the interaction of the ejecta with surrounding matter may follow the macronova.

slide-74
SLIDE 74

The radio - flare (Nakar & Piran 2011) Testing the Macronova interpretation

A long lasting radio flare due to the interaction of the ejecta with surrounding matter may follow the macronova.

slide-75
SLIDE 75

The radio - flare (Nakar & Piran 2011) Testing the Macronova interpretation

A long lasting radio flare due to the interaction of the ejecta with surrounding matter may follow the macronova.

Supernova -> Supernova remnant GRB -> Afterglow Macronova -> Radio Flare

slide-76
SLIDE 76

Search for the flare from GRB 130603B by the EVLA

slide-77
SLIDE 77

Search for the flare from GRB 130603B by the EVLA

slide-78
SLIDE 78

Search for the flare from GRB 130603B by the EVLA

slide-79
SLIDE 79

Radio limits on Magnetars

Horesh + 16 060614 130603B

slide-80
SLIDE 80

Do GRBs need magnetars?

Quasars eject magnetic jets. => GRBs also have magnetic jets => Mangetars But quasars produce magnetic jets without magnetars

slide-81
SLIDE 81

Where?

Prompt? Afterglow? Is impossible to have both from the same magnetar?

slide-82
SLIDE 82

If a magnetar did this What did that?

slide-83
SLIDE 83

If a magnetar did this What did that?

slide-84
SLIDE 84

Energy Generation

Hotokezaka, Sari & TP + 16

N+n N+p e νe γ

tf = 2⇡3 G2

F

~7 m5

ec4 ≈ 104sec

˙ E = ✏e mec2 tf ✓ t tF ◆−α 1 ⌧ ∝ 1 E Z d3pe Z d3pν Relativistic 1 ⌧ ∝ E5 → ↵ = 6/5 Newtonian 1 ⌧ ∝ E7/2 → ↵ = 9/7

GF

tf = 2⇡3 G2

F

~7 m5

ec4 ≈ 104sec

˙ E = ✏e mec2 tf ✓ t tF ◆−α 1 ⌧ ∝ d dE Z d3pe Z d3pν

E3 or E3/2

E3

Z Z Relativistic 1 ⌧ ∝ E5 → ↵ = 6/5 Newtonian 1 ⌧ ∝ E7/2 → ↵ = 9/7

slide-85
SLIDE 85

Efficiency

Hotokezaka, Wajano +…TP 16; Barnes +

Photon losses: The ejecta becomes

  • ptically thin to gamma-rays long

before it becomes optically thin to

  • ptical/IR photons => photon

leakage during the macronova peak (Hotokezaka + 16) Electron losses: Unlike previous believes not all the electrons energy is deposited (Barnes + 16)

slide-86
SLIDE 86

Summary

The nIR flare that followed the short GRB 130603B could have been a

  • Macronova. If so than:

✓Short GRBs arise from mergers. ✓Gold and other A>130 elemets are

produced in mergers. (But large mej). A radio flare may confirm this! A second & third Macronovae suggest a BH-NS merger

244Pu suggests that R-process production

is in rare events. Cocoon produces a short bright macronova We wait for the sGRB-GW coincidence

slide-87
SLIDE 87

Summary

The nIR flare that followed the short GRB 130603B could have been a

  • Macronova. If so than:

✓Short GRBs arise from mergers. ✓Gold and other A>130 elemets are

produced in mergers. (But large mej). A radio flare may confirm this! A second & third Macronovae suggest a BH-NS merger

244Pu suggests that R-process production

is in rare events. Cocoon produces a short bright macronova We wait for the sGRB-GW coincidence

slide-88
SLIDE 88

Summary

The nIR flare that followed the short GRB 130603B could have been a

  • Macronova. If so than:

✓Short GRBs arise from mergers. ✓Gold and other A>130 elemets are

produced in mergers. (But large mej). A radio flare may confirm this! A second & third Macronovae suggest a BH-NS merger

244Pu suggests that R-process production

is in rare events. Cocoon produces a short bright macronova We wait for the sGRB-GW coincidence