Impact of substructures on predictions of dark matter annihilation signals Julien Lavalle Institute & Dept. of Theoretical Physics, Madrid Aut. Univ. & CSIC
DESYTheoryAstroparticle, Hamburg 16V2011
Outline Why should we take substructures into account?
●
A quick look at what “boost factor” means
●
Going into more technical details
●
Recent results for gammarays and cosmic rays
●
Conclusions & perspectives
●
Why should we take subhalos into account? They are predicted by the theory of structure formation, while their features (minimal scale, number density, inner profile) depend on dark matter properties: eg CDM versus WDM. ●
They are “observed” in Nbody simulations of structures on different scales (from the galaxy cluster scale down to the dwarf galaxy scale); but current resolution limit => Mres ~ 104 Msun. ●
They might be observed in Nature above scales able to accrete the baryons efficiently enough (dwarf galaxies). ●
The smallscale clustering of dark matter can increase the annihilation rate quite significantly. ●
=> Selfconsistency +++ increase discovery / exclusion potential.
Credit to the original idea
WIMPs freeze out (chemically and kinematically) in a radiationdominated universe before BBN. They feed gravitational perturbations down to their freestreaming scale. ●
Perturbations start to grow efficiently after radiationmatter equivalence (z ~ 3500): DM seeds have very large densities at that time. 10 6 Msun objects enter the nonlinear regime around z ~ 100. ●
Galaxies form at z~6 in a much less dense universe, and should contain many of these very dense cores of dark matter, called clumps. ●
●
Some of these clumps may survive tidal disruption and increase the annihilation rate.
Boost factor: a simplistic view Smooth galaxy
Clumpy galaxy
Usual assumption: simulations provide us with the net density profile function (r) => keep in mind that it comes from a fit of <>
If clumps are considered:
This allows to define the socalled boost factor:
Cosmic messengers: boost related to their spatial origin
Boost factor ? … well, in fact, boost factors Smooth galaxy
Clumpy galaxy
The volume over which the average is performed depends on the cosmic messenger!
Boost factor ? … well, in fact, boost factors
Observer
The volume over which the average is performed depends on the cosmic messenger! 1) Prompt gammarays: point a telescope to a certain direction, and average over a volume set by the angular resolution
Boost factor ? … well, in fact, boost factors
Observer
The volume over which the average is performed depends on the cosmic messenger! 1) Prompt gammarays: point a telescope to a certain direction, and average over a volume set by the angular resolution a) To the Galactic center: the smooth halo is singular, clumps have no effect, B ~ 1
Boost factor ? … well, in fact, boost factors
Observer
The volume over which the average is performed depends on the cosmic messenger! 1) Prompt gammarays: point a telescope to a certain direction, and average over a volume set by the angular resolution a) To the Galactic center: the smooth halo is singular, clumps have no effect, B ~ 1 b) To high latitudes/longitudes: the smooth halo contributes much less, B>>1
Boost factor ? … well, in fact, boost factors
Observer
The volume over which the average is performed depends on the cosmic messenger! 1) Prompt gammarays: point a telescope to a certain direction, and average over a volume set by the angular resolution a) To the Galactic center: the smooth halo is singular, clumps have no effect, B ~ 1 b) To high latitudes/longitudes: the smooth halo contributes much less, B>>1 2) Cosmic rays: stochastic motion, define energydependent propagation scale. a) Large propagation scale: if enough to feel regions close to GC, then B ~ 1 b) Small propagation scale: if we are sitting on a clump, then B>>1, otherwise B moderate
Summary pictures Gamma rays: Bergström et al (1999), assuming that clumps spatially track the smooth halo
Caution: ● Statistical meaning only ● Energy dependence of propagation depends on the species (nuclei/electrons)
Antimatter cosmic rays: Lavalle et al (2007,2008)
Going into more details: a statistical approach (1) General expression for the flux measured on Earth
The Green function encodes the propagation properties => trivial for gammarays
Subhalos: pointlike sources provided G does not vary too much over the object
Going into more details: a statistical approach (2) Define the flux pdf for subhalos
The flux pdf is completely set by:
The average subhalo flux is entirely defined (same way for variance)
Remind: the subhalo properties are fully set by its mass and its concentration (and position in the Galaxy)
Use Nbody info: Via Lactea II versus Aquarius Via Lactea II: Diemand et al (2008) Aquarius: Springel et al (2008)
http://www.mpagarching.mpg.de/aquarius/
MWlike halos with ~ 1 billion particles of ~103 M⊙ > 50,000300,000 subhalos with masses > 10 6 104.5 M⊙ Slightly different cosmologies: WMAP3 vs WMAP5 (8 = 0.74 vs 0.9) http://www.ucolick.org/~diemand/vl/index.html
Gammaray studies in: Kuhlen et al (2008) – VL2 Springel et al (2008) – AQ Overall DM
Subhalos
Mass function Subhalo mass function from Aquarius (Springel et al, 2008)
Power law mass function
= 1.9 f(>3 104 Msun) = 13.2 % Press & Schechter (1974): = 2 (theory) Aquarius (Springel et al): ~ 1.9 Via Lactea (Diemand et al): ~ 2.0 NB: resolution limit => assume scale invariance Calibrate subhalo mass content from simulations:
=> Get total number of subhalos: The subhalo mass content is determined by the minimal mass Mmin and the slope , and is calibrated from the mass fraction resolved in Nbody simulations => Extrapolation down over > 10 OM !!!
The small scale issue (see review by T. Bringmann (2009)) The free streaming scale depends on the time of kinetic (≠ chemical) decoupling of WIMPs from the primordial soup. The weaker the collision rate, the smaller the free streaming scale and the cutoff mass. Subhalo mass down to 1011103 M⊙ (SUSY). The lighter the denser. Tidal effects ? Large survival fraction (Berezinsky et al, 2008)
T. Bringmann arXiv:0903.0189
Extrapolation down to 106 Msun
Luminosity function Analytical luminosity function (NFW, concentration from Bullock et al 2001) (Lavalle et al, 2008)
Luminosity per mass decade
<= Impact of Mmin depends on : ~ 1.9: constant luminosity per mass decade ~ 2.0: luminosity dominated by small objects Similar dependence
Concentration function Concentration vs mass and location in the MW
Average subhalo luminosity vs distance to GC
Concentrations: 1) Large theoretical uncertainties (impact of cosmological inputs). 2) Tidal effects: concentrations get larger when closer to the GC (demonstrated in VL2 and Aquarius).
Average luminosity strongly affected!
Spatial distribution: a selfconsistent method Trivial cases: 1) given from Nbody analysis (still to check consistency) 2) subhalos track the host halo: dP/dV = rho(r)/MMW (i) Global fit to the Nbody simulation (eg NFW)
(ii) Adding subhalos means splitting the global fit into a smooth + clumpy components Warning !!!
often assumed = in the past
(iii) Use Nbody prescriptions: subhalo distribution cored in the center. in Via Lactea, antibiased relation: subhalo distrib r global smooth distrib
Spatial distribution: a selfconsistent method Trivial cases: 1) given from Nbody analysis (still to check consistency) 2) subhalos track the host halo: dP/dV = rho(r)/MMW (i) Global fit to the Nbody simulation (eg NFW)
(ii) Adding subhalos means splitting the global fit into a smooth + clumpy components Warning !!!
often assumed = in the past
(iii) Use Nbody prescriptions: subhalo distribution cored in the center. in Via Lactea, antibiased relation: subhalo distrib r global smooth distrib
Gammaray skymaps (see Steve Blanchet's seminar) Empirical diffuse emission model: template maps from EGRET (Cillis & Hartman 05) But EGRET is no longer a reference: our Bg = EGRET – 50%
Johannesson (Moriond 2009) – Abdo et al (2009)
Sensitivity to individual subhalos Pieri, JL, Bertone & Branchini (2009)
Galactic center: astrophysical contributions not under control, notably cosmic ray electrons. Subhalos: clean signal if located at high latitude, no counterpart at lower energies ... but have to be very massive and nearby to be observable.
N <10 objects detectable with Fermi in 5 years. Model A: 40 GeV WIMP going to bbbar Model B: 100 GeV WIMP going to WW
Boost factors for positrons and antiprotons Pieri, JL, Bertone & Branchini (2009)
See also Lavalle et al (2007,2008)
Dark Matter subhalos: energydependent boost factor < 5 (modulo variance) Positron flux
Positron fraction
Antiproton flux
Pieri, JL, Bertone & Branchini (2009) using results from Via Lactea II (Diemand et al) and Aquarius (Springel et al) see early calculations in Lavalle et al (20072008) Important features: ● 40 GeV WIMP (bbbar) excluded by antiproton constraints ● 100 GeV WIMP (WW) at the edge of tension with the antiproton data ● 100 GeV WIMP going t o e+e can fit the PAMELA data; but pulsars not included => background must be known before any claim.
Highresolution is not the end of the story: what about baryons? VL2/Aquarius + baryons from Sofue et al 09
Subhalos: Small sclae issue + more efficient tidal stripping in the disk and the bulge leading to a dark disk (cf Read et al) Galactic center: Adiabatic compression might increase the DM density, but competition with dynamical friction from SF feedback reheating the gas. => Still large uncertainties Governato et al 10: CDM + highresolution baryon physics can lead to cores
Kinematics data are available for the MW: → try to use them to improve predictions
Conclusions & perspectives Highresolution Nbody simulations, e.g. Via Lactea II (Diemand et al) and Aquarius (Springel et al), provide new insights on the subhalo phasespace. ●
Small scale issue + still large scatter in predictions when using different prescriptions
●
The prospect to observe subhalos with Fermi is weak: only a few objects are detectable in 5 years [astrophysical diffuse emission to be deeply revisited – connection with cosmic rays]. ●
The antiproton signal provides interesting constraints for < 100 GeV candidates. The local positron background is not under control. ●
Caveats: still large theoretical uncertainties due (i) to understanding of small scales (ii) to the impact of baryons and (iii) to our current (mis)understanding of the Galactic diffuse emission. Relevant to subhalos and the Galactic center. Many ongoing studies on that => discovery/exclusion potential is likely to increase in the future. ●
Complementary methods mandatory. [LHC, direct detection, multimessenger wavelengthscale astrophysical signals]. Difficult, but maybe soon ... ●
Backup
Complementarity with antimatter signals
Main arguments: ● DM annihilation provides as many particles as antiparticles ● Antimatter cosmic rays are rare because secondary products ● DMinduced antimatter CRs may have specific spectral properties
But: ● We must control the backgrounds ● Antiprotons are secondaries, what about positrons ? ● Do the natural DM particle models provide clean signatures?
Single object wandering around The game one can play: ● Assume a single DM source at any distance d to the Earth. ● Assume a WIMP mass and its annihilation final states. ● Search for the brightness necessary to fit PAMELA. ● Check against other data (gamma, antiprotons, etc.) Bringmann, Lavalle & Salati (2009)