Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 1
Gravitationslinsen Rotationskurven
Indirekter Nachweis der DM
( Annihilation der DM in Materie-Antimaterie)
Direkter Nachweis der DM
( Elastische Streuung an Kernen)
Nachweismethoden der DM
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 2
Gravitationslinsen
ART: Die Ausbreitung von Licht ändert sich
beim Durchgang durch
ein Gravitationsfeld
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 3
Gravitationslinsen
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 4
Colliding Clusters Shed Light on Dark Matter
Observations with bullet cluster:
•Chandra X-ray telescope shows distribution of hot gas
•Hubble Space Telescope and others show distribution of dark matter from weak gravitational lensing
•Distributions are clearly different after collision->
dark matter is weakly interacting!
Rot: sichtbares Gas
Blau: dunkle Materie aus Gravitations- potential
dunkel
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 5
Simulation der “Colliding Clusters”
http://www.sciam.com/
August 22, 2006
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 6
Center of the Coma Cluster by
Hubble space telescope ©Dubinski
Discovery of DM in 1933 Zwicky, Fritz (1898-1974
Zwicky notes in 1933 that
outlying galaxies in Coma cluster moving much faster than mass calculated for the visible
galaxies would indicate
DM attracts galaxies with more force->
higher speed.
But still bound!
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 7
Dunkle Materie im Universum
Die Rotationskurven von
Spiralgalaxien sind weitgehend flach, während die leuchtende Materie eine abfallende Kurve erwarten lässt. Erklärung: dunkle Materie.
Spiralgalaxien bestehen aus einem zentralen Klumpen und einer sehr dünnen Scheibe
leuchtender Materie, welche von einem nahezu sphährischen,
sehr ausgedehnten Halo umgeben
ist.
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 8
Messung der Masse durch Newtons Gravitationsgesetz v=ωr
v1/r
mv
2/r=GmM/r
2Milchstraße
Cygnus Perseus
Orion
Sagittarius
Scutum Crux
Norma
Sun (8 kpc from center
)
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 9
Do we have Dark Matter in our Galaxy?
Rotationcurve Solarsystem
rotation curve Milky Way
1/ r
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 10
Estimate of DM density
DM density falls off like 1/r
2for v=const.
Averaged DM density “1 WIMP/coffee cup”
(for 100 GeV WIMP)
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 11
• Für Ensemble wechselwirkender Systeme im mechanischen Gleichgewicht gilt
• Für N Galaxien also N(N-1)/2 Teilchenpaaren
Für N groß: und
0 2 E
Kin E
Pot
2 0 ) 1
(
22
r N m
G N v
m N
N 1 N m
2 m
2 N m M 2 rGv2Erwarte also für ´Gas` gravitativ wechselwirkender Teilchen M r !
Aber dann v
2M/r = konst -> flat rot. curve
Virialsatz
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 12
Kandidaten der DM
Problem: max. 4% der Gesamtenergie
des Univ. in Baryonen nach CMB und BBN.
Sichtbar nur 0.5%, d.h. 3.5% in obigen Kandidaten möglich. Rest der DM muss aus nicht-baryonischen Materie bestehen.
Probleme:
ν < 0.7% aus WMAP Daten
kombiniert mit Dichtekorrelationen der Galaxien.
•Für kosmische Strings keine Vorhersagekraft.
•Abweichungen von Newtons
Gravitationsgesetz nicht plausibel.
•WIMPS ergeben nach Virialtheorem flache Rotationskurven.
In Supersymmetrie sind die WIMPS Supersymmetrische Partner der CMB d.h. Spin ½ Photonen (Photinos genannt).
†
†
?
?
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 13
• 95% of the energy of the Universe is non-baryonic
23% in the form of Cold Dark Matter
• Dark Matter enhanced in Galaxies and Clusters of Galaxies but DM widely distributed in halo->
DM must consist of weakly interacting and massive particles -> WIMP’s
• Annihilation with <σv>=2.10
-26cm
3/s, if thermal relic
From CMB + SN1a + surveys
DM halo profile of galaxy cluster from weak lensing
If it is not dark It does not matter
What is known about Dark Matter?
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 14
How do particles annihilate?
e+ ~~
e-
q q
LEP collider:
e+e- annihilation
~~ ~~
e e e
photon annihilation
In CM: Eq=Ee
monoenergetic quarks
from monoenergetic leptons Quarks fragment into jets, mostly light mesons:π+,π-,π0 π0 decays 100% in 2 photons
So as many photons as charged particles from annihilation
On averaged: 37 photons pro annihilation into quarks at LEP
Spectral shape VERY WELL MEASURED
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 15
Indirect Dark Matter Searches
Annihilation products from dark matter annihilation:
Gamma rays
(EGRET, FERMI)
Positrons (PAMELA)
Antiprotons (PAMELA) e+ + e-
(ATIC, FERMI, HESS, PAMELA)
Neutrinos (Icecube, no results yet)
e-, p drown in cosmic rays?
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 16
Diffuse gamma rays
Great advantage of pointing to the source and propagation is
„straightforward“ without dependence on magnetic field and diffusion,
which plagues charged particles.
Astrophysical point sources can be
pinpointed and subtracted.
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 17
Propagation of charged cosmic rays (CR)
Present models use isotropic
propagation, i.e. same diffusion constant in halo and disc.
This does not allow for
significant convection, since CR‘s do not return to disc->
too little secondary production from CR hitting gas in disc
HOWEVER, significant
convection observed by ROSAT CRs propagation can be
described by diffusion and convection, very much like a drop of ink inside streaming water (with water
velocity=convection velocity)
Radiaactive clocks like
10Be
determine time from source to Sun (10
7yrs) Need slow
diffusion in disc, but particles in halo drift to outer space with convection
With convection little flux of
charged particles from DMA,
since particles drift away.
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 18
Present models: isotropic propagation
Is this right?
Isotropic propagation leads to
“propagation enhancement”:
of charged particles: trapping of charged
particles in “leaky” Galaxy for a long time->
Flux of gamma rays from DMA
Flux of antiprotons in such propagation models, Although we KNOW from LEP that fragmentation gives many more photons than antiprotons
Not nessarily!
CONVECTION = negligible with isotropic
propagation in contrast to observation
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 19
Thermische Geschichte der WIMPS
Thermal equilibrium abundance Actual abundance
T=M/22
C o m o vi n g n u m b er d en si ty
x=m/T
Jungmann,Kamionkowski, Griest, PR 1995
WMAP -> h
2=0.1130.009 ->
<v>=2.10
-26cm
3/s
DM nimmt wieder zu in Galaxien:
1 WIMP/Kaffeetasse 10
5<ρ>.
DMA (ρ
2) fängt wieder an.
T>>M: f+f->M+M; M+M->f+f T<M: M+M->f+f
T=M/22: M decoupled, stable density
(wenn Annihilationrate
Expansions- rate, i.e. =< v>n (x
fr) H(x
fr) !)
Annihilation in leichtere Teilchen, wie
Quarks und Leptonen -> 0’s -> Gammas!
Einzige Annahme: WIMP = thermisches
Relikt, d.h. im thermischen Bad des
frühen Universums erzeugt.
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 20
DM Annihilation in Supersymmetrie
Dominant
+ A b bbar quark pair B-Fragmentation bekannt!
Daher Spektren der Positronen, Gammas und Antiprotonen bekannt!
f
f
f f
f f
Z Z W
W
0f ~
A Z
Galaxie = Super B-Fabrik mit Rate 10
40x B-Fabrik
≈37
gammas
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 21
Usual astrophysicist’s search strategies
Particle physicist: get rid of model
dependence by DATA DRIVEN calibration
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 22
R
Sun disc
Basic principle for indirect dark matter searches
R
Sun
bulge
disc
From rotation curve:
Forces: mv
2/r=GmM/r
2or M/r=const.for v=cons.
and (M/r)/r
2 1/r
2for flat rotation curve
Expect highest DM density IN CENTRE OF GALAXY
IF FLUX AND SHAPE MEASURED IN ONE DIRECTION, THEN FLUX AND SHAPE FIXED IN ALL (=180) SKY DIRECTIONS!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
R
1THIS IS AN INCREDIBLE CONSTRAINT, LIKE SAYING I VERIFY
THE EXCESS AND WIMP MASS WITH 180 INDEPENDENT MEAS.
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 23
Instrumental parameters:
Energy range: 0.02-30 GeV Energy resolution: ~20%
Effective area: 1500 cm
2Angular resol.: <0.5
0Data taking: 1991-1994 Main results:
Catalogue of point sources Excess in diffuse gamma rays
EGRET (Energetic Gamma Ray Experiment Telescope.) Data publicly available from NASA archive
EGRET excess
Hunter et al. 1997
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 24
Woher erwartet man Untergrund?
Quarks fromWIMPS
Quarks in protons
Background from nuclear interactions (mainly p+p-> π0 + X -> + X inverse Compton scattering (e-+ -> e- + )
Bremsstrahlung (e- + N -> e- + + N)
Shape of background KNOWN if Cosmic Ray spectra of p and e- known
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 25
Untergrund + DM Annihilation beschreiben Daten
W. de Boer et al., 2005
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 26
Analyse der EGRET Daten in 6 Himmelsrichtungen
A: inner Galaxy (l=±300, |b|<50) B: Galactic plane avoiding A
C: Outer Galaxy
D: low latitude (10-200)
E: intermediate lat. (20-600) F: Galactic poles (60-900)
A: inner Galaxy B: outer disc C: outer Galaxy
D: low latitude E: intermediate lat. F: galactic poles
Total 2 for all regions :28/36 Prob.= 0.8 Excess above background > 10σ.
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 27
Fits für 180 statt 6 Regionen
180 regions:
8
0in longitude 45 bins
4 bins in latitude 0
0<|b|<5
05
0<|b|<10
010
0<|b|<20
020
0<|b|<90
0 4x45=180 bins
bulge
disk
sun
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 28
Halo density on scale of 300 kpc
(from normalization factors in 180 sky directions)
Sideview Topview
Cored isothermal profile with scale 4 kpc
Total mass: O(10
12) solar masses
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 29
Halo density on scale of 30 kpc
Sideview Topview
Enhancement of inner (outer) ring over 1/r
2profile 6 (8).
Mass in rings 0.3 (3)% of total DM
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 30
The Milky Way and its 13 satellite galaxies
Canis Major
Tidal force ΔF
G 1/r
3Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 31
Tidal streams of dark matter from CM and Sgt
CM
Sgt Sun
From David Law, Caltech
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 32
Canis Major Dwarf orbits from N-body simulations to fit visible ring of stars at 13 and 18 kpc
Canis Major leaves at 13 kpc tidal stream of
gas(10
6M
☉from 21 cm line), stars (10
8M
☉,visible), dark matter (10
10M
☉, EGRET)
Movie from Nicolas Martin, Rodrigo Ibata
http://astro.u-strasbg.fr/images_ri/canm-e.html
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 33
Core of Canis Major Dwarf
just below Galactic disc
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 34
Tidal disruption of Sagittarius
Movie from Kathryn Johnston (Wesleyan University )
http://astsun.astro.virginia.edu/~mfs4n/sgr/
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 35
N-body simulation from Canis-Major dwarf galaxy
prograde retrograde
O bs er ve d st ar s R=13 kpc,φ=-20
0,ε=0.8
Canis Major (b=-15
0)
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 36
EGRET Excess predicts shape of rotation curve!
Outer Ring Inner Ring
bulge
totalDM
1/r2 halo
disk
Rotation Curve
Normalize to solar velocity of 220 km/s
R0=8.3 kpc
R0=7.0
v
R/R0 Inner
rotation curve
Outer RC
Black hole at centre:
R0=8.00.4 kpc
Sofue &Honma
Note 1: Absolute value of rotation curve depends on distances.
But chance of slope can ONLY
be explained by ringlike structure.
Note 2: fact that shape of DM halo can describe shape of RC implies
that EGRET excess has exactly right
intensity to deliver grav. potential!
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 37
Gas flaring in the Milky Way
no ring
with ring
P M W Kalberla, L Dedes, J Kerp and U Haud, http://arxiv.org/abs/0704.3925
Gas flaring needs EGRET ring with mass of 2.10
10M
☉!
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 38
Enhancement of inner (outer) ring over 1/r
2profile 6 (8).
Mass in rings 0.3 (3)% of total DM
Inner Ring coincides with ring of dust and H
2->
gravitational potential well!
H
24 kpc coincides with ring of neutral hydrogen molecules!
H+H->H
2in presence of dust->
grav. potential well at 4-5 kpc.
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 39
Probleme die durch DM Annihilation gelöst werden
Big Bang
Spektren der Gamma- strahlung für Untergrund und DMA
Teilchenphysik Kosmologie
Astroteilchenphysik
23%DM, Hubble Annihilation
Strukturformation Kosmische Strahlung (Gammastrahlen)
Astronomie
Rotationskurve
Ringförmiger Struktur von Sternen bei 14 kpc Ringförmiger Struktur von Wasserstoff bei 4 kpc
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 40
What about
Supersymmetry?
Assume mSUGRA
5 parameters: m
0, m
1/2, tanb, A, sign μ
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 41
Fundamental questions of modern physics
Particle physics Cosmology
What is the origin of mass?
Why forces different strength?
Why hydrogen atom neutral?
What is Dark Matter?
What is Dark Energy?
Why no antimatter?
How did galaxies form?
Magic solution: SUPERSYMMETRY
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 42
Motivation of SUSY in Particle Physics
1.Unification with Gravity
2.Unification of gauge couplings
3.Solution of the hierarchy problem
4.Higgs mechanism by radiative corrections 5.No quadratic divergencies,
i.e. theory valid to high energies 6.Dark matter in the Universe
7.Superstrings
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 43
Expected SUSY mass spectra in mSUGRA
mSUGRA: common masses m
0and m
1/2for spin 0 and spin ½ particles
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 44
Gauge unification perfect with SUSY spectrum from EGRET
With SUSY spectrum from EGRET + WMAP data and start values of couplings from final LEP data perfect gauge coupling unification!
U pd at e fr om A m al di , dB , F ür st en au , PL B 26 0 19 91
SM SUSY
Also b->s and g-2 agree within 2σ with SUSY spectrum from EGRET
NO FREE PARAMETER
WdB, C. Sander,PLB585(2004).
e-Print: hep-ph/0307049
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 45
• Was wissen wir über Dunkle Materie?
massive Teilchen
23% der Energie des Universums
schwache Wechselwirkung mit Materie Annihilation mit <σv>=2.10
-26cm
3/s
• Annihilation in Quarkpaare ->
Überschuss in galaktischen Gammastrahlen Tatsächlich beobachtet (EGRET Satellit)
• WIMP Masse 50-100 GeV aus Spektrum
• Verteilung der Dunklen Materie
• Data konsistent mit Supersymmetrie
Dunkle Materie, was ist das?
From CMB + SN1a
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 46
New gamma ray data from FERMI Satellite
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 47
FERMI EGRET
FERMI requires factor 3 smaller boost factor then EGRET,
But fit including DM improves also for FERMI fit significantly
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 48
Clustering of DM -> boosts annihilation rate
Clumps with M
min-> dominant contribution -> MANY clumps in given direction -> same
boostfaktor in all directions Annihilation SQUARE of DM density
Clustersize: ≈ Solarsystem?
M
min 10
-8-10
-6M סּ ?
Steeply falling mass spectrum.
Boost factor <
2>/< >
2 20-2000
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 49
CR spectra
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 50
Resurs Dk1 Satellite
300 - 600 km
Bottom Scintillator Transition
Radiation Detector
(removed for tech.reasons)
Time of Flight Counters
Silicon
Tracker and Permanent Magnet
Si-W
Electromagnetic Calorimeter
Neutron Detector Anticoincidence
Shield 1.2 m
20.5 cm
2sr
~450 kg
~10 T
The PAMELA Satellite Experiment (launched July 2006)
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 51
Positron fraction
PAMELA, positron and antiproton measurements
Positrons: excess Galp rop Pam ela
Nature 458:60,2009,arXiv:0810.4995
Antiprotons: NO excess
Antiproton/proton ratio
+prelim. new data, Boezio, Pamela-WS 2009 (O. Adriani et. al., PRL (2009)[0810.4994])
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 52
Interpretations
Many possibilities:
Background from hadronic showers
with large electromagnetic component -> a
p->0 astrophysical sources
pulsars -> a
pulsar positron acceleration in SNR -> a
sec locality of sources -> a
SNR dark matter annihilation -> a
DMA leptophilic?
bound states?
Kaluza-Klein
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 53
Truth?
Depends on whom you ask!
My assumption:
|Data>= a
p->0|Background> + a
DMA|DMA>
+ a
sec|SNR> + a
local|SNR(x)> + a
pulsar|Pulsar>
Unitarity must be fulfilled. However,
each component has enough uncertainty
to saturate observations…
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 54
Cosmic ray spectra
Lipari, PAMELA Workshop, 2009
1 TeV E
-2.7E
-3.3E
-3.0e- mainly from SNR
e+ mainly p+p e p+p 3p+p+X
3 orders of magn.
<2 orders of magn.
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 55
dP/P
2~ 0.004 2.5 TV, p rejection = 10
-5(ECAL +TRD); Δx=10µm; Δt=100ps
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 56
2009 2010
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 57
AMS to be launched in 2010
AMS
Space S huttle
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 58
Simulation of AMS on ISS
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 59
The AMS superconducting Magnet at CERN (2008)
59 Coils
He Tank
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010
60
60Magnet inside vacuum tank
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 61
Hope to suppress background with AMS
because of excellent particle identification Should provide excellent data for
indirect DM searches
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 62
Neutralino Annihilation channels
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 63
Neutralino-Quark elastic scattering
Wirkungsquerschnitte sehr klein, weil Higgs nur an Masse koppelt, aber u,d Quarks praktisch keine Masse haben.
Sehr empfindlich für s-Quark Anteil im Nukleon.
Squark Austausch sehr klein, wenn Squark schwer
Z-Austausch klein, wenn Neutralino hauptsächlich Bino ist (Bino koppelt nur an elektrische Ladung)
Zusätzlich geringer Impulsübertrag bei Streuung, so weit von Masse des ausgetauschten Teilchens -> Unterdrückung
Spin independent Spin dependent
Resultat: ( N) 10 Größenordnungen kleiner als ( ) (Annihilation)
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 64
Direkter Nachweis von WIMPs
Wir gehen davon aus, dass
DM Neutralino oder WIMP ist.
Es ist kalte DM, d.h. Impuls<<Masse (oder E
2=p
2+m
2 m
2, da p=mv mit
v 10
-3c und m 100 GeV
Geschwindigkeitsverteilung der WIMPs in einem Gravitationsfeld folgt wie
bei Gas in der Atmosphäre
Maxwell-Boltzmann-Verteilung e
-Ekin/kTmit häufigster Wert v=270 km/h
χ χ
Science Voisinage réseau.ico
E
R~ E
kin(1 - cos) Neutralino kann wegen R-Paritätserhaltung
NUR elastische Streuung an Kernen durchführen
Streuung von nicht-relativ. Teilchen meist
koherent, d.h. Wellenlänge des einlaufenden
Teilchens hat de Broglie Wellenlänge =h/p
größer als Kernradius, so es kann einzelne
Kerne nicht auflösen und Rückstoß wird an
den gesamten Kern abgegeben. Wirkungs-
querschnitt A
2(A= Anzahl der Nukleonen)
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 65
0
0
WIMPs elastically scatter off nuclei => nuclear recoils Measure recoil energy spectrum in target
Direct Detection of WIMPs
Spin independent Number of nuclei
2(coherent scattering on all nuclei!)
Spin dependent
Spin dependent and indep.
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 66
Principles of WIMP detection
• Elastic scattering of a WIMP on a nucleus inside a detector
• The recoil energy of a nucleus with mass
For
• This recoil can be detected in some ways :
Electric charges released (ionization detector)
Flashes of light produced (scintillation detector)
Vibrations produced (phonon detector)
m
N10
3v
c
2 2
(max) 2
2( )
recoil x N
N
E v m m
m m
10
610
recoil N
E
m keV
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 67
Direkter Nachweis von WIMPs
Berechnung des Streuwirkungsquerschnitt an einem Kern kompliziert:
Koherente Streuung am ganzen Kern meistens dominant, aber bei Streuung kann auch Drehimpuls eine Rolle spielen Dann wird abhängig vom Spin S der Kerne im Detektormaterial. Spin S ist gegeben durch Differenz der Nukleonen mit Spin up und Spin down.
Koherenz geht verloren bei Stößen mit hohem Impuls- übertrag q, also wenn die Wellenlänge klein gegenüber Kernradius R ist oder Kohärenzbedingung q · R « 1
Impulstransfer q = A ·10
-3GeV
Kernradius R~ 1.14 fm · A ⅓
R ~ 7 GeV
-1· A ⅓
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 68
Direkter Nachweis von WIMPs
Koherenzbedingung meistens nur erfüllt für Kerne bis A=50, d.h. perfekt für Neutralinomassen von ca. 50 GeV, denn bei gleicher Kern und WIMP Masse wird q max, weil dann reduzierte Masse = M · M
N/(M
N+ M ) maximal wird.
Wenn Koherenzbedingung nicht erfüllt, dann Kernmassenverteilung wichtig, wird beschrieben durch Formfaktor (Fouriertransformierte der Massenverteilung)
Bei sehr leichten Kernen wird Verstärkung durch Koherenz
der Streuung A
2gering und spinabh. Streuung wird wichtig
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 69
Neutralino-Quark elastic scattering
scalar interaction
5 5
( ) ( ) ( ) ( ) ....
q q
L f qq d
q
q
spin-dep. interaction
• The other terms are velocity-dependent contributions and can be neglected in the non-relativistic limit for the direct detection.
• The axial vector currents are proportional to spin operators in the non-relativistic limit.
Effective
Lagrangian
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 70
Direkter Nachweis von WIMPs
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 71
Direct detection event rates
Jodi Cooley, SMU, CDMS Collaboration
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 72
Detection challenges
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 73
Background Rejection
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 74
Shielding
Underground +
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 75
Direct Dark Matter Detection
CRESST ROSEBUD CUORICINO
DAMAZEPLIN I UKDM NaI LIBRA
CRESST II ROSEBUD CDMSEDELWEISS
XENON
ZEPLIN II,III,IV HDMSGENIUS
IGEXMAJORANA DRIFT (TPC)
E
RPhonons
Ionization Scintillation
Large spread of technologies:
varies the systematic errors, important if positive signal!
All techniques have equally aggressive projections for future performance But different methods for improving sensitivity
L. Baudis
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 76
WIMP Searches Worldwide
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 77
Diskutiere nur 4 Beispiele:
Edelweiss und CDMS (Halbleiterdetektoren:
Ionisation und Wärme)
DAMA/Libra (Szintillator) XENON (Flüssigkeit:
Ionisation und Szintillation)
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 78
Der Edelweiss Detektor
Messprinzip eines Halbleiter-Bolometers. Kommt es zu einem
elastischen Stoß eines WIMP-Teilchens mit einem Atomkern des Germanium-Kristalls führt der Kern-Rückstoß zu einer
Temperaturerhöhung des Kristalls, die über ein Thermometer
registriert wird. Gleichzeitig ionisiert der Ge-Kern das Material
in seiner Umgebung, was zu einem Ladungssignal führt, das an
den Oberflächenelektroden ausgelesen wird.
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 79
Array von
Phasenübergangs- Thermometern
Schnelle (großflächige) Auslese
von Phononen DM-Suche mit Tieftemperatur-Kalorimetern / CDMS
Si oder
Ge Einkristall
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 80 Rückstoß-Energie(keV)
Elektron-Rückstöße
Kern-Rückstöße
Ionisations-Energieschwelle
0 0.5
1 1.5
0 50 100 150 200
Kalibration mit
252Cf
Kalibration eines Ge-Bolometers durch Bestrahlung mit einer
252Cf-Neutronenquelle: Deutlich erkennbar sind zwei
Ereignispopulationen, die durch das Verhältnis von Ionisations- zu Rückstoß-Energie separiert
werden können. Die auf das Ionisationssignal angelegte Energieschwelle (grüne Kurve) entspricht einer Rückstoßenergie von 3.5keV. Die Bänder
beschreiben die Bereiche, in denen 90% der Elektron- bzw.
Kern-Rückstöße liegen.
Kalibration
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 81
Edelweiss Experiment
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 82
CDM detectors
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 83
Ionization measurement in CDMS
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 84
SQUID: Superconducting Quantum Interference Device zur Messung von minimalen Änderungen der magnetischen
Feldstärke (bis 10
-14T !)
Phonon measurement in CDMS
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 85
CDMS in Soudan mine in Minnesota (USA)
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 86
Fiducial Volume removes edges
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 87
-Flüssiges Xe als Detektormaterial (LXe)
-hohe Dichte gute Selbstabschirmung kompakte
Detektoren XENON
-hohe Massenzahl
-niedrige Energieschwelle der Rückstoßenergie
-gute Ionisations- und Szintillationseigenschaften -Betriebstemperatur
„leicht“ zu halten (180 K)
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 88
Noble liquids
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 89
Ionization and Scintillation in Xe
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 90
Double Phase Detector Concept
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 91
The XENON10 Experiment (10 kg)
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 92
XENON10 Backgrounds
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 93
Cross section limits
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 94
Annual Modulation as unique signature
June June
Dec Dec
95 97 99 101 103 105
-0.5 -0.1 0.3 0.7 1.1 1.5
±2%
0 25 50 75 100 125
-0.5 -0.1 0.3 0.7 1.1 1.5
Background WIMP Signal
June June Dec
Annual modulation: v, so signal in June larger than
in December due to motion of earth around sun (5-9% effect).
June
v
0galactic center
Sun 230 km/s Dec.
L. Baudis, CAPP2003
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 95
• DAMA NaI-1 to 4: 58k kg.day
• DAMA NaI-5 to 7: 50k kg.day
• Full substitution of electronics and DAQ in 2000
The data favor the presence of a modulated signal with the
proper features at the 6.3 σ C.L.
0
0cos with t =152.5, T=1.00 y A t t
Running conditions stable at level < 1%
DAMA/NaI 1 to 7: Riv.N.Cim 26 n.1. (2003) 1-73
Schael, EPS2003
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 96
Daten bis 2008 Daten bis 2008
Modulation nur in 2-6 keV
Region -> leichte WIMPs
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 97
Comparison with indirect searches and SUSY
Note: N
90%CL=n
<
90%CLv>N
targetTo get
90%CLone has to assume v and n
:
v assume Maxwellian
and NO corotation of DM halo n
: assume DM mass from
rotation curve to be completely diffuse.
Theory: x-section can be
order of magnitude lower due to matrix element uncertainties Conclusion: can easily move up exp. limits by order of magn.
and move down theory by order of magnitude.
CDMS
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 98
Large uncertainties in direct scattering x-section
Ellis, Olive, Savage, arXiv:0801.3656
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 99
Clustering of DM
An artist picture of what we should see if our eyes were sensitive to 3 GeV gamma rays and we are flying
with 220 km/s (=speed of sun) through the DM halo
Diemand et al., Nature 2005, astro-ph/0501589:
The earth passes through a dark matter mini-halo every 10,000 years, an encounter which lasts for about 50
years, therefore most of the time the earth is within an UNDERDENSE region of dark matter.
Consequently the averaged DM density on a large scale (from the rotation curve) has very little to do with the
LOCAL DM density!
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 100
Es gibt interessante Hinweise für Teilchencharakter der DM:
a)Überschuss an Gammastrahlung von EGRET gemessen (aber FERMI misst weniger und Konsistenz mit
Antiprotonenfluss steht nach aus, abhängig vom Propagationsmodell)
b) Jährliche Modulation der Signale in Libra/DAMA (aber inkonsistent mit anderen Experimenten)
c)Überschüsse in Positronen (PAMELA Satellit)
(aber Pulsare oder andere Quellen bieten gute Lösung)
Zusammenfassung
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 101
Zukunft
Sind Beobachtungen konsistent mit SUPERSYMMETRIE?
LHC Experimente werden ab 2011 klären ob dies stimmt.
Antwort: Wenn ja,dann hat WIMP Eigenschaften vermutlich
ähnlich einem Spin ½ Photons, d.h.
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 102
1. „The dark ages“: Vom Begriff her entsteht der Eindruck, dass mit der Rekombination das Universum quasi schlagartig dunkel wurde.
Tatsächlich muss es jedoch noch für eine ganze Weile sehr hell und heiß gewesen sein. Von t = 380.000 yr (Rekombination) mit T ≈ 3.000 K (weißglühend) bis zur Rotglut (T ≈ 750 K) bei t ≈ 4 Myr war das
Universum von sichtbarer Strahlung erfüllt. Allerdings dauerte es dann
≈ 200 Myr, bis die ersten Sterne leuchteten.
2. Neutrinomasse: Die durchschnittliche Neutrinomasse beträgt (aus
WMAP-Messungen) m
ν< 0.23 eV. Da Elektronneutrinos vermutlich die geringste Neutrinomasse besitzen, müsste deren Masse deutlich unter 0.23 eV liegen. Im Großexperiment KATRIN soll die Masse der
Elektronneutrinos bzw. deren Obergrenze bestimmt werden, wobei die Nachweisgrenze von KATRIN bei 0.2 eV liegen soll. Wenn die
Auswertung der WMAP-Daten korrekt ist, wäre damit KATRIN überflüssig, oder?
A: eine unabhängige Bestätigung, dass die Neutrinomassen tatsächlich so klein sind, ist immer gut.
Fragen
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 103 3. Polarisation der CMB: Die Polarisation setzt m. E. voraus, dass Elektronen in der
LSS in lokalen Bereichen keine stochastische Bewegung ausgeführt haben bzw.
deren Spins nicht isotrop verteilt waren, da ansonsten die CMB-Polarisation
„random“ sein müsste. Was wären Ursachen für großräumige Bereiche von Anisotropien der Elektronen-Flüsse?
A: Die relative Bewegung der Photonen besitzt durch die CMB Anisotropie (vor allem Quadrupolasymmetrie) eine bevorzugte Richtung gegenüber Elektronen, wodurch eine Polarization entsteht.
4. Annihilation von Materie/Antimaterie: Protonen und Antiprotonen wurden bis auf wenige 10-10 durch Annihilation in Photonen umgewandelt. Derselbe Prozess hat für Elektronen und Positronen stattgefunden. Erstaunlich ist, dass offenbar
exakt der identische winzige Anteil η an Elektronen „übriggeblieben“ ist, wie der der Protonen, denn sonst wäre das Universum nicht elektrisch neutral. Woher kommt die identische Asymmetrie für Protonen/Elektronen und deren
Antiteilchen? (Klar: Im Urknall war das Universum auch elektrisch neutral, aber warum ist die Asymmetrie identisch?)
A: Man geht davon aus, dass es eine B-L Symmetrie gibt, d.h. B-L=konstant. Hier ist B die Baryonzahl und L die Leptonzahl. Diese Symmetrie erzeugt oder vernichtet immer gleich viele Leptonen und Baryonen. B-L ist in allen bekannten
Wechselwirkungen erhalten (und von den einfachsten GUT's vorhergesagt).
Fragen
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 104
Fragen
5. Der Urknall: Vor der inflationären Phase war auf kleinstem Raum immense Energie (= Masse) konzentriert. Der Schwarzschild-Radius dieser Masse war jedoch wesentlich größer als die Ausdehnung der Massenkonzentration. (Beispiel: Bereits für m
≈ 10 μg ist der Schwarzschildradius rc = Planklänge lP). Also hätte es eigentlich bei einer derartig hohen Massenkonzentration, wie sie bei τP vorlag, gar nicht zu einem Big Bang kommen dürfen. (Aus einem Schwarzen Loch entweicht nichts!) Oder aber, die Energie/Masse, die sich im Urknall ausgebreitet hat, ist erst während des Urknalls entstanden (Umwandlung „falsches Vakuum“ in Energie?). Dann müsste lokal die
Massenkonzentration immer kleiner als die kritische „Schwarzschild-Masse“ gewesen sein, d.h. bereits zur Zeit der Quantenfluktuationen dürften recht kleine Raumbereiche in der Regel nicht in kausalem Kontakt mit Nachbarbereichen gestanden haben. Die
daraus resultierenden Irregularitäten wurden dann „eingefroren“ und sind heute in der CMB nachweisbar.
Oder aber, ganz einfach: Die gesamte Masse des Universums war bei t = τP in einem Raumbereich lP konzentriert, der Schwarzschildradius dieser Masse entsprach aber bereits seiner heutigen Dimension, d.h., der Big Bang lief in einem Schwarzen Loch ab.
Aber dann: Woher kommt diese Masse?
A: gute Frage. Universum so groß wegen Inflation, die nach einer Symmetriebrechung entstand, z.B. die Brechung einer GUT Symmetrie in die bekannten Kraefte. Bei der Symmetriebrechung entstehen Higgsfelder, die die Austauschteilchen Masse geben und so die Kraft ausschalten, aber gleichzeitig durch die Vakuumenergie Inflation
hervorrufen und die freiwerdende Energie in Masse umwandeln.
D.h. vor der Inflation war noch keine Masse vorhanden und Gesamtenergie null.
Es ist jedoch nicht ausgeschlossen, dass Urknall in einem SL stattfand (siehe nächste Folie).
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 105
Größe und Dichte eines schwarzen Loches.
Radius eines SL:
R = 2GM/c
2, d.h.
wächst mit Masse!
Masse unseres
Universums, die kritische Dichte von 10
-29g/cm
3(10
23M
☼) entspricht, liegt auf diese Linie, d.h. es ist nicht
ausgeschlossen, dass wir in einem SL leben.
J. Luminet
Wim de Boer, Karlsruhe Kosmologie VL, 29.01.2010 106
Fragen
6. Kosmische Zeitskala: Eine Zeitskala ist abhängig von der Stärke des Gravitationsfeldes in dem Bereich, in dem die Zeit gemessen wird. Aufgrund der
extremen Massenkonzentration im frühen Universum müsste für große z eine andere (verzögerte) Zeitskala gelten als heute. Oder könnte eine solche Zeitdilatation lediglich ein „äußerer“ Beobachter feststellen? (den es natürlich grundsätzlich nicht geben
kann.)
A: die unterschiedlichen Zeitskalen können nur gemessen werden von zwei
Beobachtern, die „Frequenzen von Gammastrahlen“ miteinander vergleichen. Daher praktisch schwierig.
7). Energie der Neutrinos aus der Entkopplung: Die Energie der Photonen aus der LSS skaliert mit 1/S (λS) Da Neutrinos Ruhmasse besitzen, müsste deren gesamte Energie bei der Entkopplung nahezu vollständig Ekin sein ( 2.5 – 3.5 MeV). Die Neutrinos müssten dann kinetische Energie verlieren, also mit der Zeit langsamer werden. Wird diese Energie dem Raum
übertragen (Energieerhaltung!), also z.B. durch Zunahme der
Vakuumenergie? Oder müssen wir Neutrinos hier quantenmechanisch betrachten und wie bei Photonen der Neutrinoenergie eine „Frequenz“
zuordnen, die S abnimmt?
A: relativistische Materie geht mit 1/S4, nicht relativistische mit 1/S3. Daher werden die Neutrinos, wenn sie relativ. sind, erst mit 1/S4 skalieren /wie Photonen) und bei Temp. T<mν als 1/S3 (wie Teilchen).