Analysis of supernova neutrino fluxes
and neutron star properties
Andrea Gallo Rosso Gran Sasso Science Institute Astroparticule et Cosmologie (APC) Advisors: F. Vissani and C. Volpe 5th April 2019
Analysis of supernova neutrino fluxes and neutron star properties - - PowerPoint PPT Presentation
Analysis of supernova neutrino fluxes and neutron star properties Andrea Gallo Rosso Gran Sasso Science Institute Astroparticule et Cosmologie (APC) Advisors: F. Vissani and C. Volpe 5 th April 2019 List of publications A. Gallo Rosso et al .
Andrea Gallo Rosso Gran Sasso Science Institute Astroparticule et Cosmologie (APC) Advisors: F. Vissani and C. Volpe 5th April 2019
…
Theoretical Papers Theoretical Reviews XENON collaboration
1
…
Theoretical Papers Theoretical Reviews XENON collaboration
1
SN 1987A
SN H SN I SN II Si SN Ia He SN IIb SN IIL SN IIF SN IIpec SN IIP SN IIn SN Ib SN Ic
(linear) (faint) (peculiar) (plateau) (narrow) no yes yes no rich poor He dominant H dominant
CORE COLLAPSE THERMONUCLEAR
3
CORE-COLLAPSE SUPERNOVA EXPLOSION
4
DELAYED EXPLOSION
5
DELAYED EXPLOSION
5
DELAYED EXPLOSION
5
DELAYED EXPLOSION
5
NEUTRINO MESSENGERS
S.P. Mikheyev and A.Y. Smirnov, Sov. J. Nucl. Phys. 42 (1985).
6
WHAT CAN WE LEARN FROM SUPERNOVA NEUTRINOS?
7
SN 1987A: THE ONLY NEUTRINO SIGNAL (SO FAR)
8
Kamiokande-II
∼ 10×
Super-Kamiokande
∼ 10×
Hyper-Kamiokande
MANY DETECTION CHANNELS — ENERGY, TIME, FLAVOR
9
νe νe νx=νμ,νμ,ντ,ντ
10 20 30 40 0.0 0.5 1.0 1.5 2.0 2.5 Neutrino energy [MeV] Time integrated flux [1010 MeV-1 cm-2]
TIME INTEGRATED FLUX (FLUENCE) Total energy E ⇔ normalization Mean energy ⟨E⟩ ⇔ 1st moment Pinching α ⇔ width 3 PARAMETERS × 3 SPECIES = 9 D.O.F.
10
TIME INTEGRATED FLUX (FLUENCE) Total energy E ⇔ normalization Mean energy ⟨E⟩ ⇔ 1st moment Pinching α ⇔ width 3 PARAMETERS × 3 SPECIES = 9 D.O.F.
10
NUMBER OF PARAMETERS ARBITRARILY REDUCED
with MSW transformation w/o equipartition (Etot ̸= Ei/6) Etot known up to 13% but for spectral shape (i.e. pinching) fully known
2Lu et al. Phys. Rev. D 94, 023006 (2016).
11
DIFFICULTY IN RECONSTRUCTING THE BINDING ENERGY
12
WHAT CAN WE LEARN FROM SUPERNOVA NEUTRINOS?
without any usual assumptions?
the determination of key properties?
13
REFERENCE PAPERS
14
REFERENCE PAPER
no.12, 006
15
SUPERNOVA PARAMETERS
DETECTORS
16
TIME-INTEGRATED FLUXES (FLUENCES)
d F0
i
d Eν = Ei 4πD2 (αi + 1)(αi+1) Γ(αi + 1) Eαi ⟨Ei⟩αi+2 exp [ −(αi + 1) E ⟨Ei⟩ ]
νe νe νx=νμ,νμ,ντ,ντ 10 20 30 40 0.0 0.5 1.0 1.5 2.0 2.5 Neutrino energy [MeV] Time integrated flux [1010 MeV-1 cm-2]
4M.T. Keil et al., Astrophys. J. 590 (2003).
17
NEUTRINO FLAVOR TRANSFORMATIONS IN SUPERNOVAE
{ Fνe = F0
x
Fνe = |Ue1|2 · F0
νe + (1 − |Ue1|2) · F0 x
18
TOTAL ENERGIES
i = 0.5 × 1053 erg
MEAN ENERGIES
PINCHING PARAMETERS
i = 2.5
Teun Hocks, Measuring
19
— HYPER-KAMIOKANDE EXTRACTED EVENTS —
Expected distribution Extracted events 10 20 30 40 50 60 2 4 6 8 10
×103
Ee [MeV]
νe + p → e+ + n (76 × 103 expected events)
True distribution Extracted events 5 10 15 20 25 30 35 100 200 300 400 500 600 Ke [MeV]
ν + e− → ν + e− (4 × 103 expected events)
20
ν + 16O → ν + X + γ
֒ → NCR = IBD + ES + OS
1 2 5 10 15 20 25 3
Energy [MeV] Signal [a.u.]
21
NEUTRINO-OXYGEN CROSS SECTION σOS(Eν) ≈ κ · σ0 · (Eν/MeV − 15)4 [5]
10th parameter κ
5J.F. Beacom and P. Vogel, PRD 58 (1998) 053010.
22
LIKELIHOODS Lj (param.) ∝
Nbin
∏
i=1
νni
i
ni e−νi with j = IBD, ES LNCR (param.) ∝ exp [ −(nNCR − NNCR)2 2NNCR − (κ − 1)2 2σ2
κ
] 3 ANALYSES IBD → L = LIBD IBD + ES → L = LIBD × LES IBD + ES + NCR → L = LIBD × LES × LNCR
23
PRIOR 0.2 × 1053 erg ≤ Ei ≤ 1.0 × 1053 erg 5.0 MeV ≤ ⟨Ei⟩ ≤ 30 MeV 1.5 ≤ αi ≤ 3.5 0.8 ≤ κ ≤ 1.2 CONDITION log L ≥ log Lmax − 1 2Adof,CL with ∫ A χ2
dof(z)dz = C.L. 24
Comparison with Minakata et al. (2008): Good agreement
14.6 14.8 15.0 15.2 15.4 15.6 0.3 0.4 0.5 0.6 0.7
〈E(νe)〉 [MeV] ℰ(νe) [1053 erg]
2σ
log L(Pi) ≥ log Lmax − 1 2Adof,CL with ∫ A χ2
dof(z) dz = CL 25
— THE IMPORTANCE OF MANY DETECTION CHANNELS — Degeneracy broken for νe and νx total energies
26
— THE IMPORTANCE OF MANY DETECTION CHANNELS — Degeneracy broken for νe and νx mean energies
27
IBD
True spectrum Set of values P1 Set of values P2 10 20 30 40 50 50 100 150 200 e+ energy [MeV] Events spectrum [MeV-1] ES (νe+νx) True spectrum Set of values P1 Set of values P2 10 20 30 40 50 10 20 30 40 e– energy [MeV] Events spectrum [MeV-1]
P1 P2 E(νe) [1052erg] 6.65 2.94 E(νx) [1052erg] 2 10 ⟨E(νe)⟩ [MeV] 12.8 13.5 ⟨E(νx)⟩ [MeV] 9.3 11.9 α(νe) 2.08 2.08 α(νx) 2.16 2.16
28
— STILL SOME RESIDUAL UNCERTAINTIES — νe species undetermined and almost all pinching parameters α
29
— TOTAL NEUTRINO EMITTED ENERGY —
Prior IBD IBD+ES IBD+ES+NCR 2 3 4 5 6 0.00 0.02 0.04 0.06 0.08 0.10
ℰB reconstructed [1053 erg]
[1053 erg]
IBD
3.40 ± 0.86 25.1
IBD+ES
3.27 ± 0.37 11.2
IBD+ES+NCR
3.18 ± 0.35 11.0
Super-Kamiokande
IBD IBD+ES IBD+ES+NCR 2.0 2.5 3.0 3.5 4.0 4.5 5.0 0.00 0.02 0.04 0.06 0.08 0.10
ℰB reconstructed [1053erg] Probability
[1053 erg]
IBD
3.64 ± 0.79 21.7
IBD+ES
3.10 ± 0.16 5.3
IBD+ES+NCR
3.07 ± 0.18 5.8
Hyper-Kamiokande
30
— TOTAL NEUTRINO EMITTED ENERGY IN EQUIPARTITION (Etot = Ei/6) —
IBD IBD+ES IBD+ES+NCR 2.6 2.8 3.0 3.2 3.4 3.6 3.8 0.00 0.01 0.02 0.03 0.04
ℰB reconstructed [1053erg] Probability
EB
[1053 erg]
IBD
3.15 ± 0.25 7.9
IBD+ES
3.06 ± 0.10 3.4
IBD+ES+NCR
3.023 ± 0.095 3.1
Super-Kamiokande
IBD IBD+ES IBD+ES+NCR 2.90 2.95 3.00 3.05 3.10 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07
ℰB reconstructed [1053erg] Probability
EB
[1053 erg]
IBD
3.13 ± 0.23 7.4
IBD+ES
3.130 ± 3.035 0.89
IBD+ES+NCR
3.015 ± 0.021 0.68
Hyper-Kamiokande
31
EQUATION OF STATE [7] (0.60 ± 0.05)β 1 − β/2 = EB Mc2
Fit T VII Inc Buch T IV WFF1 WFF2 WFF3 N4 AP4 AP3 MS0 MS1 GM3 ENG PAL6 GS1 GS2 PCL2 PS
0.10 0.15 0.20 0.25 0.30 0.35 0.05 0.10 0.15 0.20
Compactness β = GM/Rc2 ℰB/M
32
M = √ EBR 0.6 G [√ 1 + ϵ2 − ϵ ] with ϵ = 1 4 √ EBG 0.6 R c4 Total energy Etot + Equation of state
ℰB ± 10% (SK) ℰB ± 1% (HK)
8 10 12 14 16 0.0 0.5 1.0 1.5 2.0 2.5 3.0 Radius [km] Gravitational mass [M⊙]
8An estimation of the baryonic mass may also be needed.
33
Solving Tolman-Oppenheimer-Volkoff (TOV) equations Standard and f(R) gravity: Ricci scalar R → f(R) = R + αGR2
10 12 14 0.0 0.1 0.2 0.3 0.4 B/M c2 GR = 1 = 2 = 20
12 14 16
16 APR SLy
— SENSITIVE TO EoS AND POTENTIALLY TO EXTENDED GR —
34
— QUASI-STATIC COOLING (6 ÷ 10 s window) —
10−4 10−3 10−2 10−1 100 101 102 Lν (1052 erg s−1)
νe ¯ νe νx
4 6 8 10 12 14 16 e
ν (MeV)
10−1 100 101 102 tpost-bounce (s) 101 102
R (km)
Rpns Rν
10−1 100 101 102 tpost-bounce (s) 2 4 6 8 10 12 Tν−sphere (MeV)
Almost constant behavior
L.F. Roberts and S. Reddy (2017) arXiv:1612.03860.
35
FERMI-DIRAC BLACK BODY Pinching parameter η(α) L = −24π2c (hc)3 Li4(−eη) R2 [⟨E ⟩ F2(η) F3(η) ]4 RADIUS RECONSTRUCTION
36
— CORRELATION BETWEEN R AND α — νe species νx species L ∝ Li4 ( −eη(α)) F2[η(α)] F3[η(α)]
37
— RADII RECONSTRUCTED @ 10 kpc —
νe species νe species νx species 5 10 15 20 25 30 35 0.00 0.05 0.10 0.15
R reconstructed [km] Probability [km-1]
R R∗ [km]
% νe 11.9 19 ± 19 100 νe 11.5 7 ± 4 56 νx 11.4 11 ± 6 55 Default α ∈ [2.1, 3.5]
νe species νe species νx species 5 10 15 20 25 30 35 0.00 0.05 0.10 0.15
R reconstructed [km] Probability [km-1]
38
— RADII RECONSTRUCTED @ 2 kpc —
νe species νe species νx species 5 10 15 20 25 30 35 0.00 0.02 0.04 0.06 0.08 0.10 0.12
R reconstructed [km] Probability [km-1]
νe species νe species νx species 5 10 15 20 25 30 35 40 0.00 0.05 0.10 0.15 0.20 0.25 0.30
R reconstructed [km] Probability [km-1]
39
USING R TO CONSTRAIN α
40
41
42
DEVELOPMENTS AND IMPROVEMENTS
43