Testing astrophysical model to interpret the Auger data PhD - - PowerPoint PPT Presentation

testing astrophysical model to interpret the auger data
SMART_READER_LITE
LIVE PREVIEW

Testing astrophysical model to interpret the Auger data PhD - - PowerPoint PPT Presentation

Testing astrophysical model to interpret the Auger data PhD candidate: Antonio Condorelli Supervisors: Sergio Petrera, Denise Boncioli Admission to the third year, GSSI, 17th October 2019 Outline Restyling of the Combined fit; Application


slide-1
SLIDE 1

Admission to the third year,
 GSSI, 17th October 2019

Testing astrophysical model to interpret the Auger data

PhD candidate: Antonio Condorelli Supervisors: Sergio Petrera, Denise Boncioli

slide-2
SLIDE 2

A.Condorelli Testing astrophysical model to interpret Auger data 2

Outline

Restyling of the Combined fit; Application to a mass composition study; Interpretation of the ankle feature; Source effect mechanism; Conclusion and future perspectives;

slide-3
SLIDE 3

A.Condorelli Testing astrophysical model to interpret Auger data 3

Auger spectrum

  • What is the origin of the cosmic-ray features in the energy spectrum?
  • Energy spectrum alone remains ambiguous concerning interpretation.


 V.Verzi [PierreAugerColl.], PoS(ICRC2019)450

slide-4
SLIDE 4

A.Condorelli Testing astrophysical model to interpret Auger data 4

Combined fit of both spectrum and composition

  • A. Aab et al .,JCAP04(2017)038,arXiv:1612.07155v3
slide-5
SLIDE 5

A.Condorelli Testing astrophysical model to interpret Auger data 5

Update of the combined fit

slide-6
SLIDE 6

A.Condorelli Testing astrophysical model to interpret Auger data 6

Update of the combined fit

Working principle of the combined fit

slide-7
SLIDE 7

A.Condorelli Testing astrophysical model to interpret Auger data 7

Update of the combined fit

Working principle of the fraction fit

slide-8
SLIDE 8

A.Condorelli Testing astrophysical model to interpret Auger data 8

Working principle of the fraction fit P r e s e n t e d a t I C R C 2 1 9 !

Update of the combined fit

A.Yushkov [PierreAugerColl.], PoS(ICRC2019)482

slide-9
SLIDE 9

A.Condorelli Testing astrophysical model to interpret Auger data

17 17.5 18 18.5 19 19.5 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

H EPOS-LHC - Gumbel

17 17.5 18 18.5 19 19.5 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

He

17 17.5 18 18.5 19 19.5 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

N

17 17.5 18 18.5 19 19.5 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Fe

9

Results of fraction fit

500 600 700 800 900 1000 ]
  • 2
[g cm max X 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 D/N = 43.7/23 17.8 < lg(E/eV) < 17.9 500 600 700 800 900 1000 ]
  • 2
[g cm max X 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 D/N = 64.7/25 17.9 < lg(E/eV) < 18.0 500 600 700 800 900 1000 1100 ]
  • 2
[g cm max X 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 D/N = 31.6/26 18.0 < lg(E/eV) < 18.1 600 700 800 900 1000 1100 ]
  • 2
[g cm max X 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 D/N = 35.5/24 18.1 < lg(E/eV) < 18.2 500 600 700 800 900 1000 1100 ]
  • 2
[g cm max X 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 D/N = 32.6/24 18.2 < lg(E/eV) < 18.3 600 700 800 900 1000 1100 ]
  • 2
[g cm max X 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 D/N = 23.7/22 18.3 < lg(E/eV) < 18.4 600 700 800 900 1000 ]
  • 2
[g cm max X 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 D/N = 15.8/20 18.4 < lg(E/eV) < 18.5 600 700 800 900 1000 1100 ]
  • 2
[g cm max X 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 D/N = 32.6/22 18.5 < lg(E/eV) < 18.6 600 700 800 900 1000 1100 ]
  • 2
[g cm max X 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.22 D/N = 31.8/19 18.6 < lg(E/eV) < 18.7 600 650 700 750 800 850 900 950 1000 ]
  • 2
[g cm max X 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 D/N = 11.0/15 18.7 < lg(E/eV) < 18.8 600 650 700 750 800 850 900 950 1000 ]
  • 2
[g cm max X 0.05 0.1 0.15 0.2 0.25 D/N = 14.6/16 18.8 < lg(E/eV) < 18.9 700 800 900 1000 1100 ]
  • 2
[g cm max X 0.05 0.1 0.15 0.2 0.25 D/N = 31.4/15 18.9 < lg(E/eV) < 19.0 650 700 750 800 850 900 950 1000 ]
  • 2
[g cm max X 0.05 0.1 0.15 0.2 0.25 D/N = 18.0/13 19.0 < lg(E/eV) < 19.1 650 700 750 800 850 900 950 1000 ]
  • 2
[g cm max X 0.05 0.1 0.15 0.2 0.25 0.3 D/N = 26.3/12 19.1 < lg(E/eV) < 19.2 650 700 750 800 850 900 950 1000 ]
  • 2
[g cm max X 0.05 0.1 0.15 0.2 0.25 0.3 D/N = 11.7/12 19.2 < lg(E/eV) < 19.3 650 700 750 800 850 900 950 1000 ]
  • 2
[g cm max X 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 D/N = 24.9/9 19.3 < lg(E/eV) < 19.4 700 750 800 850 900 950 ]
  • 2
[g cm max X 0.05 0.1 0.15 0.2 0.25 0.3 D/N = 16.8/9 19.4 < lg(E/eV) < 19.5 700 750 800 850 900 ]
  • 2
[g cm max X 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 D/N = 3.3/8 19.5 < lg(E/eV) < 21.0
slide-10
SLIDE 10

Testing source effect mechanism

slide-11
SLIDE 11

Testing source effect mechanism

slide-12
SLIDE 12

A.Condorelli Testing astrophysical model to interpret Auger data 12

Auger measurements

slide-13
SLIDE 13

A.Condorelli Testing astrophysical model to interpret Auger data 13

What is the ankle?

If we assume only proton spectrum—> feature of the propagation; Transition point between galactic and extragalactic cosmic rays; Two extragalactic components; Source mechanism;

slide-14
SLIDE 14

A.Condorelli Testing astrophysical model to interpret Auger data 14

Testing a model

Consider a system in which the accelerator (also referred to as the source) is embedded in an environment in which the cosmic rays are confined by magnetic fields while interacting with the ambient radiation field. The lower the energy, the more time the nuclei have to interact before escaping, leading to a hardening of the spectrum and lightening of the composition of nuclei escaping the region surrounding the source.

Origin of the ankle in the ultrahigh energy cosmic ray spectrum and of the extragalactic protons below it, M.Unger, G. Farrar, L. Anchordoqui, arXiv:1505.02153v2

slide-15
SLIDE 15

A.Condorelli Testing astrophysical model to interpret Auger data 15

Interaction time easy to compute; Propagation and production of secondary particles; Implementation of a generic photon field; A single code for propagation inside and outside the source;

SimProp

SimProp v2r4: Monte Carlo simulation code for UHECR propagation, R.Aloisio, D.Boncioli, A.Di Matteo, A.F. Grillo, S.Petrera, F.Salamida, arXiv:1705.03729 .

What ? Why?

MonteCarlo code for propagation of particle through the Universe; Generation of a primary (proton or nucleus) and its propagation from the source to the observer; Taking into account all possible energy losses;

slide-16
SLIDE 16

A.Condorelli Testing astrophysical model to interpret Auger data 16

Interaction and escape time

The shape of the spectrum of the target photons is needed. The interaction time is given by the double integral on the cross sections and on the photon field; The escaping time is just a power law on rigidity.

Photomeson Photodisintegration Total

slide-17
SLIDE 17

A.Condorelli Testing astrophysical model to interpret Auger data 17

H He N Si Fe

slide-18
SLIDE 18

A.Condorelli Testing astrophysical model to interpret Auger data

17.5 18 18.5 19 19.5 20 20.5 log10(E)[eV]

2

10

3

10

4

10

5

10

6

10 dN/dlogE [a.u.]

Ejected Spectra

18 Origin of the ankle in the ultrahigh energy cosmic ray spectrum and of the extragalactic protons below it, M.Unger, G. Farrar, L. Anchordoqui, arXiv:1505.02153v2

slide-19
SLIDE 19

A.Condorelli Testing astrophysical model to interpret Auger data 19

Neutrino spectrum from the source

8 10 12 14 16 18 log10(E)[eV]

3

10

4

10

5

10

6

10

7

10

8

10 dN/dlogE [a.u.]

Legenda Standard configuration High Luminosity

Neutrinos

Increasing luminosity by a factor 10 Decreasing temperature by a factor 10

8 10 12 14 16 18 log10(E)[eV]

3

10

4

10

5

10

6

10

7

10

8

10 dN/dlogE [a.u.]

Legenda Standard configuration low T

Neutrinos

slide-20
SLIDE 20

A.Condorelli Testing astrophysical model to interpret Auger data 20

Application of source mechanism to SBGs

Cosmic ray transport and radiative processes in nuclei of starburst galaxies, E.Peretti, P .Blasi, F.Aharonian, G.Morlino. arXiv:1812.01996v2

slide-21
SLIDE 21

A.Condorelli Testing astrophysical model to interpret Auger data 21

Summary and future perspectives

Organizing combined fit code according to a logic structure and applications; Fit using two components. A single code for propagation inside and outside the source; How the neutrino flux (produced in the source) changes according to the source 
 parameters —> Additional observables with respect to cosmogenic neutrinos. Studying of parameter space for real sources; Studying of the diffusion process inside the source; Development of simplified analytical approach; Fit at the source. Tools Source mechanism Interpretation

slide-22
SLIDE 22

A.Condorelli Testing astrophysical model to interpret Auger data 22

Summary and future perspectives

Organizing combined fit code according to a logic structure and applications; Fit using two components. A single code for propagation inside and outside the source; How the neutrino flux (produced in the source) changes according to the source 
 parameters —> Additional observables with respect to cosmogenic neutrinos. Studying of parameter space for real sources; Studying of the diffusion process inside the source; Development of simplified analytical approach; Fit at the source. Tools Source mechanism Interpretation combined fit

Source mechanism

slide-23
SLIDE 23

A.Condorelli Testing astrophysical model to interpret Auger data 23

b a c k

  • u

p s l i d e s

slide-24
SLIDE 24

A.Condorelli Testing astrophysical model to interpret Auger data 24

Assumptions

slide-25
SLIDE 25

A.Condorelli Testing astrophysical model to interpret Auger data 25

7 −

10

6 −

10

5 −

10

4 −

10

3 −

10

2 −

10

1 −

10 Energy (eV) 2 4 6 8 10 12 14 16

9

10 × ]

  • 2

cm

  • 1

dN/dE [ev

Comparison of photon spectra

17.5 18 18.5 19 19.5 20 1 10

2

10

3

10

Time scale vs Energy

Changing the temperature

slide-26
SLIDE 26

A.Condorelli Testing astrophysical model to interpret Auger data 26

Changing the luminosity

0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 Energy (eV) 20 40 60 80 100 120 140 160

9

10 × ]

  • 2

cm

  • 1

dN/dE [ev

Comparison of photon spectra

17.5 18 18.5 19 19.5 20

1 −

10 1 10

2

10

slide-27
SLIDE 27

A.Condorelli Testing astrophysical model to interpret Auger data 27

Parameters in the fit

slide-28
SLIDE 28

A.Condorelli Testing astrophysical model to interpret Auger data 28

Fit the Auger Data

slide-29
SLIDE 29

A.Condorelli Testing astrophysical model to interpret Auger data 29

Injection

slide-30
SLIDE 30

A.Condorelli Testing astrophysical model to interpret Auger data 30

Advection and diffusion time

➤ Advection time: just the ratio between the size of the source R and the 


velocity of the wind.
 v = 500 Km/s; R = 300 pc

➤ Diffusion time is given by

Where:

➤ d = 5/3; ➤ L0 = 1 pc; ➤ R = 300 pc; ➤ (ΔB/B) = 1 ➤ if (RL > L0 ) we multiply the diffusion coefficient D0 by a factor (RL /L0)2 ➤ If (tdiff < R/c) tdiff = R/c