Dark Matters

51 downloads 241 Views 6MB Size Report
17 Apr 2009 ... Nucleosynethesis. Rotation Curves. Weak Lensing. CMB. BAO. Structure Formation. Evidence for Dark Matter: Overwhelming ...
WestCoast UCR 17 April 2009

Dark Matters

Graham Kribs University of Oregon

Evidence for Dark Matter: Overwhelming Nucleosynethesis Rotation Curves Structure Formation

Weak Lensing CMB BAO

Dark Matter: New Physics Not about: • • • • •

Tuning Fine-fining Beauty Aesthetics ...

It’s real new physics, waiting to be discovered!

We don’t know:

Mass of Dark Matter Composition of Dark Matter Interactions of Dark Matter

We do know: 1) Dark matter is dark 2) ρDM ≈ 5 ρmatter (averaging over Universe) 3) Dark matter is cold (well, not hot)

What to make of ρDM ≈ 5ρmatter? 1) Coincidence (axions, etc.) 2) Dark matter “number” related to baryon number, asymmetry! 3) Weak scale inspired coincidence (WIMPs)

1 pb 2 Ωh ≈ 0.1 ———

andidate electron events (162,000) with energy over 50 GeV are plotted as a stogram with the horizontal axis showing the sum of the ‘weighted energy action’ (F values as defined below) in the last two BGO layers and the ower width (root mean squared, r.m.s.) in the first two layers. The shower dth is calculated as

Figure 2 | ATIC-1 and ATIC-2 spectra at balloon altitude, showing good agreement with each other. The measured primary electron flux (scaled by E3) at flight altitude is shown for ATIC-1 (open squares) and ATIC-2 (filled circles). The errors are one standard deviation. Both balloon flights were from McMurdo, Antarctica, and circumnavigated that continent. ATIC-1 was a test flight in 2000–01 and the usable data correspond to an exposure of 2-4 keV 0.61 m2 sr days. ATIC-2 was a science flight in 2002–03 with an exposure of 2.47 m2 sr days. To eliminate edge effects, we restrict the incident zenith DAMA/NaI (0.29 ton!yr) DAMA/LIBRA (0.53 ton!yr) angle to be less than ,37u (cos h $ 0.8), use only the central 80% of the SiM (targetouter mass = 87.3crystals kg) (targetlayer. mass = 232.8 kg) these and eliminate events in the in each BGO Within limits, the electron detection efficiency above 60 GeV is 84% essentially independent of energy. The effective acceptance was determined as a function of particle energy considering the trigger efficiency, trajectory reconstruction efficiency and the geometrical restrictions. The effective acceptance of the instrument increases from 0.075 m2 sr at 20 GeV to 0.15 m2 sr for E . 60 GeV. Above 100 GeV, a total of 1,724 electron events were observed, with the highest energy event at 2.3 TeV. The total background is also shown in the figure as the open triangles and is a combination of unresolved protons, unidentified c-rays and atmospheric secondary electrons produced in the material (,4.5 g cm22) above the instrument. ATIC becomes background limited for electrons only above several teraelectronvolts.

i~1

2

Ei ðXi { Xc Þ =

n X

Ei

i~1

+

here Xc is the coordinate of the energy centre, Xi is the coordinate of the 0.2 ntre of the ith crystal and Ei is the energy deposited in the ith crystal. The F lue is calculated as Fn ~ðEn =SumÞ hr:m:s:i2 where En is the energy deposit BGO layer n,0.1Sum isMuller the& total energy deposit in all BGO layers and Ær.m.s.æ Tang 1987 fers to layer n (ref. 12). Each event is also fitted to an electromagnetic MASS 1989 scade profile to estimate the starting point and the depth of the cascade TS93 HEAT94+95 aximum. An event is accepted if the cascade starts above the first BGO CAPRICE94 yer, which eliminates many protons (,75%) but passes most electrons AMS98 ,90%). Next a diagonal cut in r.m.s. and F is determined for each energy n and used 0.02 to isolateHEAT00 the electrons. This removes most of the protons (2 in Clem & Evenson 2007 4 remain) and retains 84% of the electrons12. The selected electrons are PAMELA own as the dotted histogram. 0.01 0.1

1

10

100

Energy (GeV)

part from the calculated curve. They show an excess electron flux p to about 650 GeV,positron abovefraction whichwith theother spectrum drops rapidly, with a FIG. 3: PAMELA experimental data. The positron fraction the PAMELA experimentline compared with other recentIn experimental data[24, 29, 30, turnmeasured to the by ‘general’ spectrum at ,800 GeV. particular, over 31, 32, 33, 34, 35]. One to standard error bars are shown. not visible, they lie inside the e energy range 300 800 deviation GeV we observe 210 Ifelectrons, whereas data points. ALPROP predicts only 140 events, an excess of about six standard viations. a source-on/source-off method determining a showerUsing tail catcher scintillator (S4) and a neutron detector. for The ToF system provides 15 a fast signal, for acquisition and measures the standard time-of-flightdeviaand ionizagnificance’ wetriggering obtaintheandata excess of roughly four tion energy losses (dE/dx) of traversing section particles. 4). It also allows down-going particles to ons (Supplementary Information reliably identified. Multiple tracks, produced in interactions above the spectrometer, Databe recently became available from the Polar Patrol Balloon were rejected by requiring that only one strip of the top ToF scintillator (S1 and S2) layers Antarctic) flight of the BETS detector. Although of lower statistical registered an energy deposition (’hit’). Similarly no hits were permitted in either top scintilecision, results from the PPB-BETS calorimeter16 also indicate a lators of the AC system (CARD and CAT). The central part of the PAMELA apparatus is ossible structure and agree with the ATIC results (see Fig. 3), giving 12 dded confidence to the conclusion that this feature is real. We varied the source injection parameters in the GALPROP code try to reproduce the data points at 500 to 700 GeV. This required a ard injection spectrum which could not reproduce the drop in flux

Residuals (cpd/kg/keV)

hr:m:s:i ~

n X

Time (day)

2-5 keV

1,000 Residuals (cpd/kg/keV)

0.3

2

Ee3.0dN/dEe (m−2 s−1 sr−1 GeV2)

0.4

+

-

Positron fraction !(e ) / (!(e )+ !(e ))

Some Intriguing Hints...

DAMA/NaI (0.29 ton!yr) (target mass = 87.3 kg)

DAMA/LIBRA (0.53 ton!yr) (target mass = 232.8 kg)

100

10

10

100 Energy (GeV)

1,000

Time (day)

Main Focus: PAMELA antiproton fraction -3

0.4

!10

0.4 0.3

0.35

0.2

0.3

Donato 2001 (D, " =500MV) Simon 1998 (LBM, " =500MV) Ptuskin 2006 (PD, " =550MV) PAMELA

+

+

-

Positron fraction !(e ) / (!(e )+ !(e ))

positron fraction

0.02

0.25 Muller & Tang 1987

p/p

0.1

MASS 1989 TS93 HEAT94+95

0.15

CAPRICE94 AMS98

0.1

HEAT00 Clem & Evenson 2007

0.05

PAMELA

0.01 0.1

0.2

1

10

100

Energy (GeV)

0 1

10 kinetic energy (GeV)

102

G. 3: PAMELA positron fraction with other experimental data. The positron fraction

FIG. 3: The antiproton-to-proton flux ratio obtained in this work com PAMELA

asured by the PAMELA experiment compared with other recent experimental data[24, 29, 30,

32, 33, 34, 35]. One standard deviation error bars are shown. If not visible, they lie inside the

a points.

These hints (signals & no signals) strongly disfavor the MSSM neutralino

These hints (positive and negative) strongly disfavor the MSSM neutralino Motivating various theorists to theorize...

Results

Post-PAMELA Which DM can fitpapers... the data? M.Pospelov and A.Ritz, 0810.1502: Secluded DM - A.Nelson and C.Spitzer, 0810.5167: Slightly Non-Minimal DM - Y.Nomura and J.Thaler, 0810.5397: DM through the Axion Portal - R.Harnik and G.Kribs, 0810.5557: Dirac DM - D.Feldman, Z.Liu, P.Nath, 0810.5762: Hidden Sector - T.Hambye, 0811.0172: Hidden Vector - Yin, Yuan, Liu, Zhang, Bi, Zhu, 0811.0176: Leptonically decaying DM - K.Ishiwata, S.Matsumoto, T.Moroi, 0811.0250: Superparticle DM - Y.Bai and Z.Han, 0811.0387: sUED DM - P.Fox, E.Poppitz, 0811.0399: Leptophilic DM - C.Chen, F.Takahashi, T.T.Yanagida, 0811.0477: Hidden-Gauge-Boson DM - K.Hamaguchi, E.Nakamura, S.Shirai, T.T.Yanagida, 0811.0737: Decaying DM in Composite Messenger - E.Ponton, L.Randall, 0811.1029: Singlet DM - A.Ibarra, D.Tran, 0811.1555: Decaying DM - S.Baek, P.Ko, 0811.1646: U(1) Lmu-Ltau DM - C.Chen, F.Takahashi, T.T.Yanagida, 0811.3357: Decaying Hidden-Gauge-Boson DM I.Cholis, G.Dobler, D.Finkbeiner, L.Goodenough, N.Weiner, 0811.3641: 700+ GeV WIMP - E.Nardi, F.Sannino, A.Strumia, 0811.4153: Decaying DM in TechniColor - K.Zurek, 0811.4429: Multicomponent DM - M.Ibe, H.Murayama, T.T.Yanagida, 0812.0072: Breit-Wigner enhancement of DM annihilation - E.Chun, J.-C.Park, 0812.0308: sub-GeV hidden U(1) in GMSB - M.Lattanzi, J.Silk, 0812.0360: Sommerfeld enhancement in cold substructures - M.Pospelov, M.Trott, 0812.0432: super-WIMPs decays DM - Zhang, Bi, Liu, Liu, Yin, Yuan, Zhu, 0812.0522: Discrimination with SR and IC - Liu, Yin, Zhu, 0812.0964: DMnu from GC - M.Pohl, 0812.1174: electrons from DM - J.Hisano, M.Kawasaki, K.Kohri, K.Nakayama, 0812.0219: DMnu from GC A.Arvanitaki, S.Dimopoulos, S.Dubovsky, P.Graham, R.Harnik, S.Rajendran, 0812.2075: Decaying DM in GUTs - R.Allahverdi, B.Dutta, K.Richardson-McDaniel, Y.Santoso, 0812.2196: SuSy B-L DM- S.Hamaguchi, K.Shirai, T.T.Yanagida, 0812.2374: Hidden-Fermion DM decays - D.Hooper, A.Stebbins, K.Zurek, 0812.3202: Nearby DM clump - C.Delaunay, P.Fox, G.Perez, 0812.3331: DMnu from Earth - Park, Shu, 0901.0720: SplitUED DM - .Gogoladze, R.Khalid, Q.Shafi, H.Yuksel, 0901.0923: cMSSM DM with additions - Q.H.Cao, E.Ma, G.Shaughnessy, 0901.1334: Dark Matter: the leptonic connection - E.Nezri, M.Tytgat, G.Vertongen, 0901.2556: Inert Doublet DM - C.-H.Chen, C.-Q.Geng, D.Zhuridov, 0901.2681: Fermionic decaying DM J.Mardon, Y.Nomura, D.Stolarski, J.Thaler, 0901.2926: Cascade annihilations (light non-abelian new bosons) - P.Meade, M.Papucci, T.Volansky, 0901.2925: DM sees the light - D.Phalen, A.Pierce, N.Weiner, 0901.3165: New Heavy Lepton - T.Banks, J.-F.Fortin, 0901.3578: Pyrma baryons - Goh, Hall, Kumar, 0902.0814: Leptonic Higgs - K.Bae, J.-H. Huh, J.Kim, B.Kyae, R.Viollier, 0812.3511: electrophilic axion from flipped-SU(5) with extra spontaneously broken symmetries and a two component DM with Z2 parity - ...

(from Cirelli; Moriond 09)

Four basic explanations: 1) PAMELA screwed up, e.g., misidentified positrons 4 (proton flux 10 times positron flux) 2) Point sources (pulsars, supernovae?) dumping more energy than previously thought into positrons 3) Dark matter annihilation 4) Dark matter decay

Four basic explanations: 1) PAMELA screwed up, e.g., misidentified positrons 4 (proton flux 10 times positron flux) 2) Point sources (pulsars, supernovae?) dumping more energy than previously thought into positrons 3) Dark matter annihilation 4) Dark matter decay

Interlude: Direct Detection of WIMPs

One of the most striking constraints on WIMPs is direct detection:

WIMP

Nucleus

WIMP

Nucleus

If the WIMP-nucleon coupling is coherent w.r.t. mass Generally expressed in terms of an effective nucleon cross section: σn

2 µ (n,D)

I(n) σ(N) ≈ ————————— ≈ ———— µ2(N,D) I(N) A4 [for M ≥ m(Nuc)]

That allows experiments with different materials to be compared (for c.w.r.t.m.) against each other

10-40

Direct Detection Bounds

10-41

σn

10-42 -43 10

Xenon 0706 CDMS 0802

10-44 -45 10

1995

2000

2005

2010

(adapted from Gaitskell 2004)

Original (1980s) hope for WIMPs

Acquire • mass from EWSB • coupling to SM through EW interactions e.g. fourth geneneration neutrino that acquires Dirac mass with ν4R

Such WIMPs have true Weak Interactions: Vector interactions to SM with GF strength: — — µ ν4 γ ν4 q γµ q ———————————— v2 Leads to WIMP-nucleon scattering cross section σn ≈ 0.1 pb!

(for ≈ 100 GeV WIMP)

s g . e e

-37 10 4 !41

Baltz Gondolo 2004 Roszkowski et al. 2007 95% CL Roszkowski et al. 2007 68% CL CDMS II 1T+2T Ge Reanalysis XENON10 2007 CDMS II 2008 Ge CDMS II Ge combined

2

Spin!independent cross section [cm ]

10

10

10

10

!42

Ruled Out by Direct Detection Bounds

!43

!44

10

1

2

10 WIMP mass [GeV/c2]

3

10

FIG. 4: Spin-independent WIMP-nucleon cross-section upper limits (90% C.L.) versus WIMP mass. The upper curve (dash-dot) is the result of a re-analysis [17] of our previously published data. The upper solid line is the limit from this

Direct Detection Suggests: • WIMPs coupling to quarks with sub-weak interaction strength (vector “g” ≤ 0.01; or Higgs exchange; or ...) or • WIMPs couple to leptons, not quarks (or gluons) Evades direct detection constraints. or • Not thermal freezeout (forget about 1 pb scale) or • ...

WIMPs Coupling to Leptons

New Dirac Fermion: D (Harnik, GK) • New Dirac fermion D neutral under SM gauge group • Global U(1)D conserved • Interactions with SM through higher dimensional operators -- effective theory! • Choose operators we want; discard those we don’t.

Higher Dimensional Operators _ DDff ————— 2 Λ We want: f = lepton (e) vector interaction no quarks

D _ D

f f

2 1/Λ

Thermally averaged cross section

= σ0 + σ2v2 + ... Dirac fermion:

2/Λ4 M ~

(not velocity suppressed)

D _ D

l l

2 1/Λ

Matching with thermal relic abundance: Ωh2 ≈ 0.1 (1 pb/) IV.

($.# 1

/0 /0 1 ! +0 +0

($."

!!!"#$%&

!

D _ D

(!.#

(!."

leptons only!

(".#

("." !""

!#"

$""

$#" %"" '(()*+,-

%#"

&""

A.l

POSITRON DIRA

Backgrou 2 1/Λ

Determining the b l is of utmost impor ratio excess does, in culation of the bac from the Galprop c ated as secondary p other protons and l the inputs to the c and parameter dep the uncertainties in Harnik,both GK sign propagate consistent way to t

Now onto indirect detection implications...

To make indirect DM annihilation predictions... Astrophysics Propagation: - diffusion - energy loss Backgrounds: - secondary production - pulsars (neglected) Abundance - average density - local clumpiness “BOOST factor”

Particle Physics - annihilation rate - annihilation channels Effective theory allows: + - annihilation into e e , - no other collider constraints (M > 100 GeV)!

Propagation Diffusion: charged particles move in galactic magnetic field modeled as a random walk. Models (convection; reacceleration); and fit to observations (B/C ratio...) Energy loss: light charged particles lose energy by synchrotron radiation & inverse Compton off starlight and CMB

We used “Galprop” [Moskolenko & Strong]. It requires as input, among other things: - proton (nuclei) spectrum - electron spectrum

For dark matter... Source

Propagated

while diffusion randomizes the direction. Typical scale for energy loss is ~ kpc, i.e., our galactic neighborhood.

Secondary production p

p, nuclei

+ π ->

+ e

... K+->μ+-> ... ->e+

Position and antiproton spectrum determined essentially by: input proton spectrum + propagation

PAMELA proton data

! ! y 2.75 r Proton flux * E a

in

e r P

m i l

(statistical errors only)

Power-law fit: ~ E-# # ~ 2.76

• Proton of primary origin • Diffusive shock-wave acceleration in SNRs • Local spectrum: injection spectrum ! galactic propagation

LBM

Local primary spectral shape: " study of particle acceleration mechanism

! Elena Vannuccini ! SSI 2008 !

Vannuccini SSI 08

na

i

Pr

im l e

!

! ry

Proton flux * E2.75

(statistical errors only)

Power-law fit: ~ E-# # ~ 2.76

• Proton of primary origin • Diffusive shock-wave acceleration in SNRs • Local spectrum: injection spectrum ! galactic propagation

LBM

Local primary spectral shape: " study of particle acceleration mechanism

! Elena Vannuccini ! SSI 2008 !

Galprop

Absolute secondary positron spectrum

Propagation: “best” of DC, DR, DRB; astro-ph/0502406 (Lionetto, Morselli, Zdravkovic)

-BETS data with an index of −2.77±0.21.

Observed

+

+ e)

flux

Tang 1984 Golden et al. 1984 Boezio et al. 2000 DuVernois et al. 2001 Torii et al. 2001 Aguilar et al. 2002 Kobayashi et al. 2002 Chang et al. 2005 This work

104 2

E3J (electrons m!2 s!1 sr!1 GeV )

(e

103

1.82"10!4(E/100GeV)!3.05

102

101 100

101

102 Electron Energy (GeV)

103

104

(from BETS paper)

Tang 1984 Golden et al. 1984 Boezio et al. 2000 DuVernois et al. 2001 Torii et al. 2001 Aguilar et al. 2002 Kobayashi et al. 2002 Chang et al. 2005 This work

4

2

E3J (electrons m!2 s!1 sr!1 GeV )

10

103

Our approach (Oct 08)...

1.82"10!4(E/100GeV)!3.05

102

101 100

101

102 Electron Energy (GeV)

103

104

Fig. 12. Electron energy spectrum observed with PPB-BETS (solid circles) in comparison with the energy spectra of BETS (solid squares) ant the other observations. The dash line shows the best fit power-law function of the combined spectrum of PPB-BETS and BETS with an index of −3.05±0.05.

21

Absolute electron spectrum

(and use PAMELA flux ratio at 5 GeV to normalize e spectrum)

Tang 1984 Golden et al. 1984 Boezio et al. 2000 DuVernois et al. 2001 Torii et al. 2001 Aguilar et al. 2002 Kobayashi et al. 2002 Chang et al. 2005 This work

4

2

E3J (electrons m!2 s!1 sr!1 GeV )

10

103

1.82"10!4(E/100GeV)!3.05

102

101 100

101

102 Electron Energy (GeV)

103

104

Fig. 12. Electron energy spectrum observed with PPB-BETS (solid circles) in comparison with the energy spectra of BETS (solid squares) ant the other observations. The dash line shows the best fit power-law function of the combined spectrum of PPB-BETS and BETS with an index of −3.05±0.05.

21

Absolute electron spectrum

Which data to use?

PB-BETS data with an index of −2.77±0.21.

Electron data 5-100 GeV Tang 1984 Golden et al. 1984 Boezio et al. 2000 DuVernois et al. 2001 Torii et al. 2001 Aguilar et al. 2002 Kobayashi et al. 2002 Chang et al. 2005 This work

2

E3J (electrons m!2 s!1 sr!1 GeV )

104

103

1.82"10!4(E/100GeV)!3.05

102

101 100

101

102 Electron Energy (GeV)

103

104

12. Electron energy spectrum observed with PPB-BETS (solid circles) in com-

PB-BETS data with an index of −2.77±0.21.

Electron data 5-100 GeV Tang 1984 Golden et al. 1984 Boezio et al. 2000 DuVernois et al. 2001 Torii et al. 2001 Aguilar et al. 2002 Kobayashi Experiment et al. 2002 Chang et AMS-01 al. 2005 [29] This work

2

E3J (electrons m!2 s!1 sr!1 GeV )

104

power law index α 3.15 ± 0.04 ATIC [30]!3.05 3.14 ± 0.08 1.82"10!4(E/100GeV) BETS [31, 32] 3.05 ± 0.05 CAPRICE [33] 3.47 ± 0.34 HEAT [34] 2.82 ± 0.16 MASS [35] 2.89 ± 0.10

103

102

101 100

101

12. Electron energy spectrum

TABLE I: Our weighted least-squares best fit to the flux, Φe− (E) ∝ E −α , measured by the various exp The BETS best fit was taken from [32]; their best 2 lower energy 3data between 10 the 10 10 104 to 100 GeV is 3. [31]; the error is assumed to be 1σ. The MASS best Electron Energy (GeV) from [35]; the error is assumed to be 1σ. We empha our reported errors for the other experiments are p (95% CL) with(solid regard circles) to fitting in data (wit observedtistical with PPB-BETS comto a power law, and do not necessarily reflect the i

PB-BETS data with an index of −2.77±0.21.

Electron data 5-100 GeV Tang 1984 Golden et al. 1984 Boezio et al. 2000 DuVernois et al. 2001 Torii et al. 2001 Aguilar et al. 2002 Kobayashi Experiment et al. 2002 Chang et AMS-01 al. 2005 [29] This work

2

E3J (electrons m!2 s!1 sr!1 GeV )

104

power law index α 3.15 ± 0.04 ATIC [30]!3.05 3.14 ± 0.08 1.82"10!4(E/100GeV) BETS [31, 32] 3.05 ± 0.05 -2.8 CAPRICE [33] 3.47 ± 0.34 HEAT [34] 2.82 ± 0.16 -3.15 MASS [35] 2.89 ± 0.10

103

E E

102

-3.5 E

101 100

101

TABLE I: Our weighted least-squares best fit to the flux, Φe− (E) ∝ E −α , measured by the various exp The BETS best fit was taken from [32]; their best 2 lower energy 3data between 10 the 10 10 104 to 100 GeV is 3. [31]; the error is assumed to be 1σ. The MASS best Electron Energy (GeV) from [35]; the error is assumed to be 1σ. We empha our reported errors for the other experiments are p -2.8 -3.5 tistical (95% CL) with regard to fitting in data (wit observed with PPB-BETS (solid circles) comto a power law, and do not necessarily reflect the i

Interpretation:

12. Electron energy spectrum

E

> Φ(e) > E

Preliminary PAMELA Electron data

~ E- 3.24±0.04

Adriani; talk given in Waseda 14 Jan 09

PB-BETS data with an index of −2.77±0.21.

Electron data 5-100 GeV Tang 1984 Golden et al. 1984 Boezio et al. 2000 DuVernois et al. 2001 Torii et al. 2001 Aguilar et al. 2002 Kobayashi Experiment et al. 2002 Chang et AMS-01 al. 2005 [29] This work

2

E3J (electrons m!2 s!1 sr!1 GeV )

104

power law index α 3.15 ± 0.04 ATIC [30]!3.05 3.14 ± 0.08 1.82"10!4(E/100GeV) BETS [31, 32] 3.05 ± 0.05 -2.8 CAPRICE [33] 3.47 ± 0.34 HEAT [34] 2.82 ± 0.16 -3.15 MASS [35] 2.89 ± 0.10

103

E E

102

-3.5 E

101 100

101

TABLE I: Our weighted least-squares best fit to the flux, Φe− (E) ∝ E −α , measured by the various exp The BETS best fit was taken from [32]; their best 2 lower energy 3data between 10 the 10 10 104 to 100 GeV is 3. [31]; the error is assumed to be 1σ. The MASS best Electron Energy (GeV) from [35]; the error is assumed to be 1σ. We empha our reported errors for the other experiments are p -2.8 -3.5 tistical (95% CL) with regard to fitting in data (wit observed with PPB-BETS (solid circles) comto a power law, and do not necessarily reflect the i

Interpretation:

12. Electron energy spectrum

E

> Φ(e) > E

Secondary production ONLY (no dark matter): Propagation models: convection reaccelaration reaccelaration+break

e /(e + e )

M PAMELA0.15 B ratio data

0.1 0.07

+

-3.5 E

0.1 0.07

Ratio derived 0.05 from: • e- data

0.05

+

+

+



e /(e + e )

0.15

0.2



0.2

E-3.15

• e+ predicted 0.03 secondary

0.03 -2.8 E

10

15 20

30

50

Energy (GeV)

70

100

Harnik, GK

a in

Tang 1984 Golden et al. 1984 Boezio et al. 2000 DuVernois et al. 2001 Torii et al. 2001 Aguilar et al. 2002 Kobayashi et al. 2002 Chang et al. 2005 This work

4

10 2

(statistical errors only)

E3J (electrons m!2 s!1 sr!1 GeV )

lim

!! y r Proton flux * E2.75 Power-law fit: ~ E-# # ~ 2.76

• Proton of primary origin • Diffusive shock-wave acceleration in SNRs • Local spectrum: injection spectrum ! galactic propagation

Local primary spectral shape: " study of particle acceleration mechanism

LBM

103

DM Annihilation

1.82"10!4(E/100GeV)!3.05

102

101 100

101

102 Electron Energy (GeV)

103

104

! Elena Vannuccini ! SSI 2008 !

Galprop

Fig. 12. Electron energy spectrum observed with PPB-BETS (solid circles) in comparison with the energy spectra of BETS (solid squares) ant the other observations. The dash line shows the best fit power-law function of the combined spectrum of PPB-BETS and BETS with an index of −3.05±0.05.

Galprop

21

Absolute secondary positron spectrum

Absolute electron spectrum

Spectrum from DM

M=100 GeV Dirac Dark Matter 0.2

+

0.1 0.07

+



e /(e + e )

M = 100 GeV 0.15 B = 1 (no boost)

0.05

0.03

100 Annihilation

to e+e- only

10

20

30

50

Energy (GeV)

100

200

Harnik, GK

Dirac Dark Matter: BOOST

BOOST = 1

BOOST = 5

for

for

α

2 M

M = 100 GeV -3.5 Φe- ≈ E ; ρlocal = 0.3 GeV/cm3 M = 100 GeV Φe- ≈ E-3.15; ρlocal = 0.3 GeV/cm3

be automatic. The relic abundance of an exa calculated before in Ref. [46] u channel) scalar exchange. Left scalars give rise to a four-fermio be Fierz transformed into our effe (5) with cL = (YL g ! )2 /2 and cR is the hypercharge of the Standa g ! is the hypercharge coupling. mass of the exchanged scalar Λ to immediately re-evaluate Fig. 1 of the physical scalar states that operators. This is shown in Fig. -3.15 of the leptonic The dominance efor two reasons. First, the four-fe 3pr Standard Model fermion are local much like KK dark matter [2, 3 the right-handed leptons. Secon scale as 1/m4f , even a modest hier are lighter than squarks will ove dominant annihilation channel to handed leptons. Hence, a Dirac annihilation to charged leptons. is clear: relatively light sleptons, about 200-400 GeV, are an ines obtain a thermal relic abundance ogy and a positron signal consist A pure Dirac bino-eigenstate Z. This eliminates one source o

M=100,200,400 GeV Dirac Dark Matter

0.2

M (GeV) 100 200 400

B 5 20 80

Φ ρ

0.1

+

+



e /(e + e )

0.5

0.05

10

20 30

50

100

Energy (GeV)

200

400

≈E ; = 0.3 GeV/cm

Lesson: Need enhancement FIG. 7: Same Fig. 6, for M = 100, 200, 400 GeV. The DCfor “heavy” DM model was used for propagation, and annihilation was assumed only (M into e » e .100 GeV) not “light” DM (M ≈ 100 GeV) + −

+ e

Feature in

e?

+

500

s ! e"#e$ ATIC 08 HESS 08

! !

200

!

E3 dN#dE !GeV2 m$2 s$1 sr$1 "

! !

! ! !

! !

!

!

!

! !

! !

" !

" "

100

! !

" "

50

" "

!

"

20

10

ADDGHR (adapted) 50

100

500 Energy !GeV"

1000

5000

Fermi/GLAST feature: FSR radiation D

f

_ D

- γ f

Size of signal *dependent* on DM density profile of galactic center

Mathev, Perelstein, Spray, umber ofBirkedal, signal (blue) and background (red) events at hep-ph/0507194 the GLAST space ion area A = 10000 cm2 , exposure time T = 2 years, field of view ∆Ω =

ependent of energy. The effective acceptance was determined as a ction of particle energy considering the trigger efficiency, trajectory onstruction efficiency and the geometrical restrictions. The effective eptance of the instrument increases from 0.075 m2 sr at 20 GeV to 5 m2 sr for E . 60 GeV. Above 100 GeV, a total of 1,724 electron events e observed, with the highest energy event at 2.3 TeV. The total kground is also shown in the figure as the open triangles and is a mbination of unresolved protons, unidentified c-rays and atmospheric ondary electrons produced in the material (,4.5 g cm22) above the trument. ATIC becomes background limited for electrons only above eral teraelectronvolts.

ATIC/PPB-BETS bump @ 500 GeV?

Ee3.0dN/dEe (m−2 s−1 sr−1 GeV2)

1,000

(Cannot be explained without BOOST factor or Sommerfeld or ...)

100

10

10

100 Energy (GeV)

1,000

ure 3 | ATIC results showing agreement with previous data at lower rgy and with the imaging calorimeter PPB-BETS at higher energy. The ctron differential energy spectrum measured by ATIC (scaled by E3) at the of the atmosphere (red filled circles) is compared with previous ervations from the Alpha Magnetic Spectrometer AMS (green stars)31, AT (open black triangles)30, BETS (open blue circles)32, PPB-BETS (blue sses)16 and emulsion chambers (black open diamonds)4,8,9, with certainties of one standard deviation. The GALPROP code calculates a wer-law spectral index of 23.2 in the low-energy region (solid curve)14. he dashed curve is the solar modulated electron spectrum and shows that

Fermi/GLAST will be definitive...

Dirac Bino as Dirac Dark Matter

Dirac Bino? Yes, the gauginos in low energy supersymmetry could acquire Dirac masses, not Majorana masses. This occurs when the gauginos are paired with new fermions (beyond MSSM) in the adjoint representation. Dirac masses are automatic in supersymmetric models that preserve an R-symmetry.

Interpretating D as a Dirac Bino Resolve the 4-fermion vertex as D _ D

f f

D _ D

f ~

f

f

The dominance of leptonic annihilation results automatically given YeR=1 and some mild hierarchy, m~l < m~q (and, dim-5 Higgs operator is absent)

Matching thermal relic abundance, = 1 pb 4"" 9

58 58 9 ! 18 18

#""

.5167*)+,-..++/0123

!

&""

%""

$""

!"" !""

!#"

$"" $#" %"" '()*+,-..++/0123

%#"

&""

Harnik, GK

Absence of Direct Detection (and no Antiproton Annihilation signal) Implies

D q

~ q

q D

cannot be large. Rough estimates of direct detection bounds suggest m~q > 1.5 TeV for first generation, right-handed squarks.

Also constraint on Higgsino content of LSP from direct detection D _ D

Z

f f

Our estimate, from (g’v)2/μ2 < 0.01, µ > 600 GeV

Sketch of Spectrum

~ dR

1.5 TeV

~ H

600 GeV

~ e

200 GeV 100 GeV

B

~~ µ,τ

~L ~ uL,d

~ uR

What’s unusual? R-symmetric models have very interesting flavor properties; much of the supersymmetric flavor problem is absent. With a Bino lighter than selectron, different from smuon/stau, suggests observable LFV (work in progress) In any case, sleptons become the prime targets for LHC...

Summary • Strong constraints and hints towards the particle nature of DM; MSSM in trouble. _ _• One pair D,D with DDll operator can: - thermally produce Ωh2 ≈ 0.1 relic abundance - avoid direct detection - explain PAMELA ratio with minimal boost factor - explain WMAP haze, right scale for DAMA, ... - can arise in R-symmetric (extended) SUSY • Collider implications of “unusual” dark matter candidates is only just beginning... stay tuned!