Galactic Cosmic rays as a dark matter probe Carmelo Evoli - - PowerPoint PPT Presentation

galactic cosmic rays as a dark matter probe
SMART_READER_LITE
LIVE PREVIEW

Galactic Cosmic rays as a dark matter probe Carmelo Evoli - - PowerPoint PPT Presentation

Galactic Cosmic rays as a dark matter probe Carmelo Evoli (Universitt Hamburg) Oslo University - 22th of April 2015 (Hess) Radiazione altamente penetrante A 100 years old discovery! Primi studi sistematici Victor Hess


slide-1
SLIDE 1

Galactic Cosmic rays as a dark matter probe

Carmelo Evoli (Universität Hamburg)

Oslo University - 22th of April 2015

slide-2
SLIDE 2

A 100 years old discovery!

(Hess)

  • Radiazione altamente penetrante
  • Primi studi sistematici

Victor Hess

  • W. Kolhorster
  • W. Bothe
slide-3
SLIDE 3

High-energy photons or charged particles?

slide-4
SLIDE 4

High-energy photons or charged particles?

Bruno Rossi in his laboratory in Florence Uber die Eigenschaften der durchdringenden Korpuskularstrahlung usw. 16"i Mit zunehmender Dicke der Streuschicht mmmt zun~chst die H~ufigkeit der Koinzidenzen zu, wegen der Zunahme der erzeugten Sekund~irteilchen. Bald aber macht sich die Wiederabsorption der Sekund~rtefichen in der Streuschicht bemerkbar; die Intensit~t der austretenden Sekund~rstrahlung (und damit die H~ufigkeit der dreifachen Koinzidenzen) s~eigt daher immer langsamer an, geht dutch ein Maximum und setzt sich endlich (bei geniigend dicken Bleischichten) ins Gleichgewicht mit der Prim~rs~rahlung. Die Stelle des Maximums wird im wesentlichen dutch die Hs der Sekund~r- strahlung bestimmt. Nun erreichen unsere Kurven das Maximum zwischen 10 und 20 g/era2; wir schlie~en daraus, dal~ das mittlere Durchdringungs- vermSgen der in Blei erzeugten Sekund~rstrahlen "con der GrSl3enordnung 10 g/era 2 ist.

Jo

g5 go

l] , ,f,,I l i I I 1 i t

70 gO 30 ~-0 ~0 60 70 80 80 700 g/cm g

  • Fig. 8.

Was die Absolutzahl der drei~achen Koinzidenzen betrifft, so betr~gt ihr Maximalwert bei den Messungen mit der Bleischicht in 14,6 cm Ent- fernung etwa 4%, bei den Messungen mit der Bleischicht in 1,2 cm Ent- fernung etwa 8% der Anzahl der zweifachen Koinzidenzen zwischen C 1 und C 2 (oder C 1 und C3). Sieht man yon der Absolutzahl der Koinzidenzen ab (deren Abh~ngigkeit yon der Lage der Streuschicbt A l~l~t sich leich~ durch geometrische Uberlegungen deuten), so zeigen die Kurven I und II keinen wesentlich verschiedenen Verlauf, obwohl nach der Auffassung von Heisenberg die scheinbare H~rte der Sekund~rstrahlung vom Winkel zwischen den beiden Bahnzweigen (und somit vom Abstand des Verzweigungs-

151

O~ber die Eigenschaften der durchdringenden Korpuskularstrahlung im Meeresniveau.

Von Bruno Rossi in Florenz, .~reetri. Mit 16 Abbildungen. (Eingegangen am 24. Februar 1933.) Die Absorbierbarkeit der durchdringenden Korpuskularstrahlung wurde bis zu einer Absorberdicke von 101 em Blei untersueht ; die Versuehsresultate werden im Zusammenhang mit der Frage nach dem Wesen der Ult.rastrahlung und naeh ihrer Energie diskutiert. -- Das Entstehen einer Sekund~rstrahlung in der von den Ultrakorpuskularstrahlen durchsetzten Materie wurde naehgewiesen. Die wiehtigsten Eigent0anliehkeiten dieser Erseheinung und der Einflul3 der Se- kund~rstrahlung auf die Ultrastrahlungsph~nomene wurden untersueht.

  • I. Zweek und allgemeine Methode der Untersuchung.
  • 1. I~ach den Versuchen yon Bothe und KolhSrster 1) und vom Ver-

fasser 2) well3 man, dal~ im Meeresniveau die Ultrastrahlungserscheinungen yon einer durchdringenden Korpuskularstrahlung hervorgerufen werden. Um die Eigenschaften dieser Korpuskularstrahlung eingehend zu unter- suchen, babe ich Ende vorigen Jahres eine Versuchsreihe unternommen, die insbesondere eine genauere Messung ihres mittleren Durchdringungs- vermSgens 'und Aufkl~rung tlber dio in der Materie hervorgerufenen Se- kund~rerscheinungen bezweckte. Einige der gewonnenen Resultate sind schon kurz in deutschen Zeit- schriften verSffentlicht worden~). Vollst~ndigere Berichte tiber die Einzel- versuehe sind in itMienischen Zeitschriften erschienen~). Die vorliegende Arbeit ist einer zusammenfassenden Darstellung und Diskussion aller Er- gebnisse gewidmet.

  • 2. S~mtliche Versuehe wurden mit der Koinzidenzmethode durch-

gefi~hrt, die sieh bereits bei den friiheren Versuchen yon Bothe und Kol- h5rster und vom Yerfasser f~r die Untersuehung der durchdringenden Korpusknlarstr~hlung bew~hrt hatte. Wie bekannt werden bei dieser Methode die gleichzeitigen Aussehl~ge yon zwei oder mehreren Geiger-Mi]llerschen Z~hlrohren beobachtet, die entweder den Durchgang einer und derselben Korpuskel dutch alle be- nutzten Z~hlrohre, oder auch den gleichzeitigen Durchgang einer Korpuskel 1) W. Bothe u. W. Kolh5rster, ZS. f. Phys. 56, 751, ]929. 3) B. Rossi. ebenda 68, 64, 1931. 3) ]3. Rossi, Naturwissenseh. 20, 65, 1932; Phys. ZS. 33, 304, 1932. 4) B. Rossi, Rend. Lineei 15, 734, 1932; B. Rossi u. B, CrinS, ebenda 15, 741, 1932; B. Rossi, Rieerca Seientif. 3, I, Nr. 7--8, 1932; 3, II, Nr. 7--8, 1932.

The CR telescope used by Bruno Rossi during the expedition in Eritrea

slide-5
SLIDE 5

A unique particle physics laboratory

Carl D. Anderson 6mm absorption layer

@ London’s Westminster Abbey, adjacent to Newton’s grave.

The first anti-matter evidence was found in the cosmic radiation in 1933.

63 MeV 23 MeV

B=15 kG

slide-6
SLIDE 6

Today

slide-7
SLIDE 7

Cosmic-ray flux

  • Almost a perfect power-law
  • ver 12 energy decades.
slide-8
SLIDE 8

Cosmic-ray flux

  • Almost a perfect power-law
  • ver 12 energy decades.
  • Observed at energy higher than

terrestrial laboratories!

LHC

Ivo Karlović

slide-9
SLIDE 9

Cosmic-ray flux

  • Almost a perfect power-law
  • ver 12 energy decades.
  • Observed at energy higher than

terrestrial laboratories!

  • Direct measurements versus

air-cascade reconstructions.

1/cm2/s 1/km2/century

slide-10
SLIDE 10

Cosmic-ray flux

  • Almost a perfect power-law
  • ver 12 energy decades.
  • Observed at energy higher than

terrestrial laboratories!

  • Direct measurements versus

air-cascade reconstructions.

  • Anti-matter component.
slide-11
SLIDE 11

Cosmic-ray flux

  • Almost a perfect power-law
  • ver 12 energy decades.
  • Observed at energy higher than

terrestrial laboratories!

  • Direct measurements versus

air-cascade reconstructions.

  • Anti-matter component.
  • Transition from galactic to

extra-galactic?

slide-12
SLIDE 12

Cosmic-ray flux

  • Almost a perfect power-law
  • ver 12 energy decades.
  • Observed at energy higher than

terrestrial laboratories!

  • Direct measurements versus

air-cascade reconstructions.

  • Anti-matter component.
  • Transition from galactic to

extra-galactic?

  • Energy density in equipartition

with starlight, turbulent gas motions and magnetic fields.

slide-13
SLIDE 13

The SN paradigm

LSN ∼ RSNEkin ∼ 3 × 1041 erg/s

hadronic:

  • r leptonic:

Fritz Zwicky

TeV emission X-ray contours

Aharonian et al., Nature, 2007

slide-14
SLIDE 14

The pion-bump as hadronic signature

Energy (eV)

8

10

9

10

10

10

11

10

12

10

)

  • 1

s

  • 2

dN/dE (erg cm

2

E

  • 12

10

  • 11

10

  • 10

10

)

  • 1

s

  • 2

dN/dE (erg cm

2

E

  • 12

10

  • 11

10

  • 10

10

Best-fit broken power law Fermi-LAT VERITAS (30) MAGIC (29) AGILE (31)

  • decay

π Bremsstrahlung Bremsstrahlung with Break

IC 443

Energy (eV)

8

10

9

10

10

10

11

10

12

10

1

Best-fit broken power law Fermi-LAT AGILE (19)

  • decay

π Bremsstrahlung Bremsstrahlung with Break

W44

A B

FERMI collaboration, Science, 2013

log scale!

mπ0 2

≈ E−δ

pion bump symmetric symmetric

mπ0 2

E2

γ nγ

≈ E2−δ

very steep rise!

100 MeV

slide-15
SLIDE 15

The pion-bump as hadronic signature

Energy (eV)

8

10

9

10

10

10

11

10

12

10

)

  • 1

s

  • 2

dN/dE (erg cm

2

E

  • 12

10

  • 11

10

  • 10

10

)

  • 1

s

  • 2

dN/dE (erg cm

2

E

  • 12

10

  • 11

10

  • 10

10

Best-fit broken power law Fermi-LAT VERITAS (30) MAGIC (29) AGILE (31)

  • decay

π Bremsstrahlung Bremsstrahlung with Break

IC 443

Energy (eV)

8

10

9

10

10

10

11

10

12

10

1

Best-fit broken power law Fermi-LAT AGILE (19)

  • decay

π Bremsstrahlung Bremsstrahlung with Break

W44

A B

FERMI collaboration, Science, 2013

log scale!

mπ0 2

≈ E−δ

pion bump symmetric symmetric

mπ0 2

E2

γ nγ

≈ E2−δ

very steep rise!

100 MeV

Do SNRs accelerate ENOUGH protons? Do they accelerate protons up to the knee?

slide-16
SLIDE 16

Cosmic-ray composition

Solar System Cosmic Rays

>> Galaxy size!

Primary Secondary

cτesc = X(E) ¯ nISMµ ∼ 103 kpc

slide-17
SLIDE 17

Cosmic-ray clocks

0.1 0.2 0.3 0.4 0.5 0.6 ISOMAX TOF ISOMAX CK ACE Ulysses Voyager 1, 2 IMP-7/8 ISEE-3 100 10–1 10–2 101

Ekin (GeV/n)

10Be/ 9Be

Diffusive halo model Leaky box model Ulysses

  • 129. Hams T, et al. Astrophys. J. 611:892 (2004)
  • 130. Streitmatter RE, Stephens SA. Adv. Space Res.
slide-18
SLIDE 18

you are here

L = 1-10 kpc

Galactic Propagation

slide-19
SLIDE 19

you are here

L = 1-10 kpc

Galactic Propagation

high isotropy expected!

slide-20
SLIDE 20

The diffusion equation:

Ginzburg & Syrovatsky, 1964

∂Ni ∂t −∇·(D∇−vc)Ni + ∂ ∂p ⇣ ˙ p− p 3∇·vc ⌘ Ni − ∂ ∂pp2Dpp ∂ ∂p Ni p2 = Qi(p,r,z)+∑

j>i

cβngas(r,z)σi jN j −cβngasσin(Ek)Ni

CR Diffusion in the MW

Source term:!

  • assumed to trace the SNR in the Galaxy!
  • assumed the same power-law everywhere
slide-21
SLIDE 21

The diffusion equation:

Ginzburg & Syrovatsky, 1964

∂Ni ∂t −∇·(D∇−vc)Ni + ∂ ∂p ⇣ ˙ p− p 3∇·vc ⌘ Ni − ∂ ∂pp2Dpp ∂ ∂p Ni p2 = Qi(p,r,z)+∑

j>i

cβngas(r,z)σi jN j −cβngasσin(Ek)Ni

Spallation cross-section:

  • appearance of nucleus i due to spallation of nucleus j

CR Diffusion in the MW

slide-22
SLIDE 22

The diffusion equation:

Ginzburg & Syrovatsky, 1964

∂Ni ∂t −∇·(D∇−vc)Ni + ∂ ∂p ⇣ ˙ p− p 3∇·vc ⌘ Ni − ∂ ∂pp2Dpp ∂ ∂p Ni p2 = Qi(p,r,z)+∑

j>i

cβngas(r,z)σi jN j −cβngasσin(Ek)Ni

Spallation cross-section:!

  • appearance of nucleus i due to spallation of nucleus j !
  • total inelastic cross-section: disappearance of nucleus i

CR Diffusion in the MW

slide-23
SLIDE 23

The diffusion equation:

Ginzburg & Syrovatsky, 1964

∂Ni ∂t −∇·(D∇−vc)Ni + ∂ ∂p ⇣ ˙ p− p 3∇·vc ⌘ Ni − ∂ ∂pp2Dpp ∂ ∂p Ni p2 = Qi(p,r,z)+∑

j>i

cβngas(r,z)σi jN j −cβngasσin(Ek)Ni

CR Diffusion in the MW

Diffusion tensor:


  • D(E) = D0(ρ/ρ0)δexp(z/zt)
slide-24
SLIDE 24

The diffusion equation:

Ginzburg & Syrovatsky, 1964

∂Ni ∂t −∇·(D∇−vc)Ni + ∂ ∂p ⇣ ˙ p− p 3∇·vc ⌘ Ni − ∂ ∂pp2Dpp ∂ ∂p Ni p2 = Qi(p,r,z)+∑

j>i

cβngas(r,z)σi jN j −cβngasσin(Ek)Ni

CR Diffusion in the MW

Energy losses:!

  • ionization, Coulomb, synchrotron!
  • adiabatic convection
slide-25
SLIDE 25

The diffusion equation:

Ginzburg & Syrovatsky, 1964

∂Ni ∂t −∇·(D∇−vc)Ni + ∂ ∂p ⇣ ˙ p− p 3∇·vc ⌘ Ni − ∂ ∂pp2Dpp ∂ ∂p Ni p2 = Qi(p,r,z)+∑

j>i

cβngas(r,z)σi jN j −cβngasσin(Ek)Ni

CR Diffusion in the MW

Dpp ∝ p2v2

A

D

Reacceleration:!

slide-26
SLIDE 26

Fitting local observables

vA

D0

The best constraints on the halo scale height (L > 2 kpc) are obtained from the galactic diffuse synchrotron emission
 (G.Di Bernardo, CE, et al., JCAP, 2013)

slide-27
SLIDE 27

Fitting local observables

vA

D0

The best constraints on the halo scale height (L > 2 kpc) are obtained from the galactic diffuse synchrotron emission
 (G.Di Bernardo, CE, et al., JCAP, 2013)

anti-protons here!

slide-28
SLIDE 28

Fitting local observables

vA

D0

The best constraints on the halo scale height (L > 2 kpc) are obtained from the galactic diffuse synchrotron emission
 (G.Di Bernardo, CE, et al., JCAP, 2013)

anti-protons here!

A powerful probe of the tiny Dark Matter signal!

slide-29
SLIDE 29

AMS-02 (2011-2021)

slide-30
SLIDE 30
slide-31
SLIDE 31

Is it possible being not-local?

  • we can measure the anisotropy:

  • we can observe diffuse emissions:

δ / rncr φγ ∝ Z ncr · ngas dr

argets in the Galaxy

as

Gas ring s

H I from LAB survey CO from CFA

los

slide-32
SLIDE 32

Since CR sources are more abundant in the inner Galaxy, a dipole anisotropy is expected towards the Galactic center:

You are here!

δ

x = 3D(E)

c ∇

xnCR(E,

r, t) nCR ,

The anisotropy problem

Macro Collaboration, PRD, 2003; Super-Kamiokande Collaboration, PRD, 2007

δ = 0.5, vA = 15 δ = 0.33, vA = 30

muon detectors (e.g. SuperKamiokande)

slide-33
SLIDE 33

Lorimer (2004), based on pulsar cat. fitted to EGRET data

Strong & Mattox, A&A, 1996; Strong et al., ApJ, 2000

  • CR distribution inferred from gamma-ray data (method goes back to SAS-2/COS-B era) 


is flatter than that computed assuming the observed SNR (source) profile.

The gradient problem

Case & Bhattacharya (1996) fitted to EGRET data

slide-34
SLIDE 34

FERMI (2008-2018)

slide-35
SLIDE 35

The gradient problem in the FERMI era

  • The extremely accurate gamma ray maps that Fermi is providing are useful to

trace the CR distribution throughout all the Galaxy!

slide-36
SLIDE 36

The gradient problem in the FERMI era

!"#"$%&$'(%)*$+,"-*./+012$3 4 5 6 7 8 94 95 96 97 98 +3

:9

+/)

:9

+2;&%&(/+/

:57

<=*//*>*%?+094 4 4@A 9 9@A 5 5@A B B@A 6 6@A A

CDE+-"%" +944+F'G ! +

"

<

+H+54+12$

;

I +H+94+12$

;

I +H+6+12$

;

I +H+5+12$

;

I +H+9+12$

;

I

!"#"$%&$'(%)*$+,"-*./+012$3 4 5 6 7 8 94 95 96 97 98 +3

:9

+/)

:9

+2;&%&(/+/

:57

<=*//*>*%?+094 4 4@A 9 9@A 5 5@A B B@A 6 6@A A

CDE+-"%" +944+F'G ! +

"

<

:9

+/

5

+$=

58

+94 # +H+A@8+

4

+H+6+12$J+K

;

I

+H+94+12$

L1

, +H+99+12$

L1

, +H+95+12$

L1

, +H+9B+12$

L1

, +H+96+12$

L1

,

Fermi Collaboration, ApJ, 2011

FERMI detected more γ’s than a prediction based on SNR distribution and standard CR halo: 
 more CR sources, more “dark gas” or larger CR halo?

/proton

slide-37
SLIDE 37

Diffusion is more complicate than that!

  • How do the diffusion coefficient depends on turbulence?

If the turbulent field is very low:

  • In the inner galaxy, where turbulence is high, the parallel and perp. diffusion

are similar values and the perpendicular escape is the dominant one:

If the turbulent component is comparable to the regular field:

T T⊥

  • Rarm

H 2 D⊥ D 4 × 102

  • H

4 kpc −2 D⊥ D

slide-38
SLIDE 38

How to solve the gradient problem

  • In the regions where CR sources are more abundant turbulence is higher then

perpendicular escape is faster, more CR are removed.

dependence D?ðRÞ / QðRÞ, be fixed against the

slide-39
SLIDE 39

Results

Ackermann et al, ApJ 710 (2010), II quadrant analysis Ackermann et al, ApJ 726 (2011), III quadrant analysis

τ = 0.0 (no radial dependence) τ = 0.7 τ = 1.0

dependence D?ðRÞ / QðRÞ, be fixed against the

CE, D. Gaggero, D. Grasso & L. Maccione, PRL, 2012

slide-40
SLIDE 40

Anisotropy prediction

δ = 0.6 τ = 0 δ = 0.6 τ = 0.85

CE, D. Gaggero, D. Grasso & L. Maccione, PRL, 2012

phase change!

slide-41
SLIDE 41

The dark side

slide-42
SLIDE 42

The challenge!

Andromeda : a MW-like galaxy

5 1

slide-43
SLIDE 43

The importance of being WIMP

problem!
 astrophysical emissivity is much larger than Dark one

slide-44
SLIDE 44

Why antiprotons?

  • we know the background with

good accuracy

  • in a democratic WIMP model the

ratio between DM signal and background from standard astrophysical sources is usually much larger in the antiproton channel with respect to all other indirect detection methods.

slide-45
SLIDE 45

Playing with anti-protons from DM

CE, I.Cholis, D.Grasso, L.Maccione & P .Ullio, PRD, 2012, 1108.0664

Model zt (kpc)

  • D0ð1028 cm2=sÞ
  • vA ðkm=sÞ
  • dvc=dz

ðkm=s=kpcÞ 2

B=C

2

p

(GV) 2

  • p

Color in Figs. KRA 4 0.50 2.64 0:39 14.2 2.35 0.6 0.47 0.67 0.59 Red KOL 4 0.33 4.46 1. 36. 1:78=2:45 0.4 0.3 0.36 1.84 Blue THN 0.5 0.50 0.31 0:27 11.6 2.35 0.7 0.46 0.70 0.73 Green THK 10 0.50 4.75 0:15 14.1 2.35 0.7 0.55 0.69 0.62 Orange CON 4 0.6 0.97 1. 38.1 1:62=2:35 50 0.4 0.53 0.21 1.32 Gray

30% ~1 oom

THN

THK

slide-46
SLIDE 46

Varying the halo size in the range 2 - 10 kpc

vA

D0

The best constraints on the halo scale height (L > 2 kpc) are obtained from the galactic diffuse synchrotron emission
 (G.Di Bernardo, CE, D.Gaggero, D.Grasso and L.Maccione, JCAP, 2013)

slide-47
SLIDE 47

Varying the halo size in the range 2 - 10 kpc

less than 10%

4x

~ 4 kpc ~ 8 kpc

THIN halo

The ratio of the local flux obtained considering sources with distance smaller than RS to that

  • btained with RS = ∞

(see also R. Taillet & D. Maurin, A&A, 2003)

CE, I.Cholis, D.Grasso, L.Maccione & P .Ullio, PRD, 2012, 1108.0664

Much larger uncertainty in the DM fluxes!

slide-48
SLIDE 48

Unavoidable uncertainties?

wind

CE, I.Cholis, D.Grasso, L.Maccione & P .Ullio, PRD, 2012, 1108.0664

Changing diffusion conditions in the inner Galaxy gives significant effect on the DM contribution without affecting the local

  • bservables

Only a comprehensive study including local and non-local observables may succeed in reducing safely the propagation uncertainties.

anisotropic diffusion

slide-49
SLIDE 49

How to bracket the propagation uncertainties?

protons protons+B/C

slide-50
SLIDE 50

How to bracket the propagation uncertainties?

see also Bringmann&Salati 2007

slide-51
SLIDE 51

DM bounds after PAMELA data

CE, D.Gaggero & D.Grasso, 1504.05175

slide-52
SLIDE 52

AMS-02 anomaly?

CE, D.Gaggero & D.Grasso, 1504.05175