Arxiv:0707.3260V3 [Astro-Ph] 18 Sep 2007 O Akmte.Silms Ftewr Sn Hs The- These Using Work the of Need Most the Still Without Observations Matter
Total Page:16
File Type:pdf, Size:1020Kb
Cosmology of fermionic dark matter Tillmann Boeckel and J¨urgen Schaffner-Bielich Institut f¨ur Theoretische Physik / Astrophysik, J.W. Goethe Universit¨at Max-von-Laue-Straße 1, D-60438 Frankfurt am Main, Germany (Dated: October 24, 2018) We explore a model for a fermionic dark matter particle family which decouples from the rest of the particles when at least all standard model particles are in equilibrium. We calculate the allowed ranges for mass and chemical potential to be compatible with big bang nucleosynthesis 0 0 (BBN) calculations and WMAP-data for a flat universe with dark energy (ΩΛ = 0.72, ΩM = 0.27, h = 0.7). Futhermore we estimate the free streaming length for fermions and antifermions to allow comparison to large scale structure data (LSS). We find that for dark matter decoupling when all standard model particles are present even the least restrictive combined BBN calculation and WMAP results allow us to constrain the initial dark matter chemical potential to a highest value of 6.3 times the dark matter temperature. In this case the resulting mass range is at most 1.8 eV ≤ s m ≤ 53 eV, where the upper bound scales linearly with geff (TDec). From LSS we find that similar to ordinary warm dark matter models the particle mass has to be larger than ∼ 500 eV (meaning s 3 geff (TDec) > 10 ) to be compatible with observations of the Ly α forest at high redshift, but still the dark matter chemical potential over temperature ratio can exceed unity. I. INTRODUCTION ories has been done on galactic rotation curves, stellar velocities in dwarf galaxies and the Pioneer anomaly. It In recent years observations and the high precision is unclear whether they even succeed in explaining the analysis of the cosmic microwave background radiation observations at these scales so they still have to prove [1], of large scale structure [2], of the abundancies of being an viable cosmological alternative to dark matter. light elements in the universe [3] and of the evolution of For example Ref. [15] find that a simple MONDian model the Hubble parameter by observation of the luminosity- with a low matter density gives a rather bad fit to the redshift relation of supernovae of type Ia [4], [5] have been power spectrum of the CMBR anisotropies when includ- strongly indicating that the universe is flat and mainly ing the third acoustic peak extracted from BOOMERanG consists of two unknown components, dark matter and and WMAP first year data. A quite interesting recent ap- dark energy. Dark matter appears to be consisting of proach to explain at least galactic rotation curves without non-relativistic particles that only interact gravitation- the need for exotic dark matter or a new theory of grav- ally and perhaps by weak interactions. In addition there itation has been done by Cooperstock et al. [16] using a has been a recent independent support for the dark mat- axially symmetric metric and linearized general relativity ter paradigm by combined x-ray and weak lensing ob- for disk galaxies. Although the mentioned different ap- servation of the aftermath of a collision of two galaxies proaches may be able to remove the need for dark matter clusters [6]. on and below galactic scales it is highly unclear whether they can do so for the CMB, large scale structure or weak The nature of the second unknown component, dark lensing observations. All in all there is at present no con- energy, is even less clear, it is mostly characterized by vincing alternative to the dark matter paradigm that can having a negative pressure that allows to explain the ob- explain all the mentioned independent observations. served late time exponential expansion of the universe. In the simplest models dark energy is just a constant The ΛCDM model is successful in explaining most of arXiv:0707.3260v3 [astro-ph] 18 Sep 2007 vacuum energy. Combined with a sizeble amount of cold the observations although it currently predicts an exces- dark matter this constitues the so called ΛCDM or con- sive amount of low mass dark matter halos orbiting larger cordance model. galaxies in both linear structure formation approaches For an alternative to dark energy using a large scale as well as in numerical n-body simulations, as found by local inhomogeneity see e.g. [7, 8]. Some authors have Ref. [17, 18]. Another problem of the CDM model is the even tried to unify dark matter and dark energy within prediction of cuspy inner density profiles of dark matter so called Chaplygin gas models [9, 10, 11, 12]. This ap- halos in dwarf and spiral galaxies, that was not confirmed proach faces the difficulty of having a speed of sound by observations (see e.g. Ref. [19, 20, 21]). Instead it was increasing at late times which makes it probably unable found that the core density of dark matter must to be to reproduce the matter power spectrum of density fluc- rather constant. It is unclear whether these problems can tuations on smaller scales. be resolved in more realistic simulations/models that in- Also there has been the phenomenological MOND ap- clude for example more effects from the baryonic matter, proach by Milgrom [13] using a modified Newtonian force or if these actually show that the dark matter properties law and the more recent TeVeS theory by Bekenstein [14] are different from the standard CDM case. For example that tried to explain the observations without the need Ref. [22] proposed self-interacting dark matter as a solu- for dark matter. Still most of the work using these the- tion to the cuspy-core problem, but Ref. [23] found that 2 the interesting parameter range in which self-interacting at the GUT scale and at the time of BBN. The tempera- dark matter could flatten halo cores can probably be ex- ture of equilibrium particles therefore is shifted upwards cluded by other observations. Another solution could be compared to the temperature of decoupled particles (like warm dark matter (WDM) as studied by Ref. [24], where our dark matter particles) leading to the temperature the core densities are reduced by the higher velocity dis- relation: persion of WDM. 1/3 It is usually assumed that the dark matter constituent, gs (T ) whether it is bosonic (a candidate would be the axion) or eff Dec T = TDM s (3) fermionic (for example the neutralino, the gravitino or a geff (T ) ! massive sterile neutrino), initially has a negligible chem- ical potential (or to be more exact that the condition The equilibrium relations µ = −µ¯ and T = T¯ will still T ≫ µ − m is fullfilled). The reason for this assumption be fullfilled as long as both fermions and antifermions is the expectation that the matter-antimatter asymmetry are ultra-relativistic (where barred quantities refer to an- in the universe is of the same magnitude for all particle tiparticles), although they are decoupled and no inter- species. So spin statistics are mostly ignored for dark actions take place anymore. The transition from ultra- matter and the usual approach for cold dark matter is relativistic to non-relativistic behaviour of both fluids then a free Boltzmann gas. In contrast to that we will needs to be treated numerically, meaning that we can use a free Fermi gas in our approach to find constraints on start our numerical calculation at relatively late times the allowed matter-antimatter-asymmetry through the (i.e. after BBN) and simply add the energy density, related chemical potentials and the resulting mass win- pressure etc. of the underlying photon, neutrino and dow of a thermal fermionic dark matter relic from cosmo- baryon fluids where their behaviour and magnitude is logical observations. We assume an additional conserved known from standard cosmology, WMAP-data and BBN- quantum number beyond those of the standard model to calculations. Neutrinos are treated as massless for sim- keep the particle stable, but this is common in many dark plicity. matter models. We will also explore possible implications for structure formation for such a particle species. III. ANALYTIC LIMITS FOR ENERGY II. INITIAL ASSUMPTIONS AND DENSITY AND EQUATION OF STATE MOTIVATION The energy density and equation of state (EoS) of a We assume a thermal distribution of a fermionic parti- 1 non-interacting Fermi gas with non-zero chemical poten- cle family with 4 degrees of freedom (Spin 2 ) was created tial can be calculated analytically in the ultra-relativistic asymmetrically (i.e. a surplus of particles over antiparti- as well as in the non-relativistic limit, an outline of the cles or the other way around) at an early stage somewhere 15 derivations will be given below, a more detailed descrip- around the GUT-scale (approx. 10 GeV). The assym- tion can be found in [25]. metric process then became ineffective and the surplus The energy density ǫ+ = ǫ+¯ǫ and pressure p+ = p+¯p and of fermions was conserved leading to a non-zero chem- net number density n− = n − n¯ of a free non-interacting ical potential of these particles. We then assume ther- Fermi gas in chemical equilibrium are given in natural mal decoupling at a later stage while it was still ultra- units (c = ~ = kB = 1) by the following integrals over relativistic. Thus the particles will be unaffected by the the momentum ~k of the particle: following numerous phase transitions and the subsequent reduction of the number of equilibrium particle species ∞ and their degrees of freedom gi. Assuming overall en- + g 3~ ~ 1 tropy conservation from the dark matter decoupling scale ǫ = d k E(k) − (2π)3 exp( E µ )+1 up to today demands: Z0 T ∞ gs a3T 3 = gs a3T 3 (1) 1 eff Dec eff later + d3~k E(~k) (4) exp( E−µ¯ )+1 where Z0 T T 3 7 T 3 gs = g i + g i (2) eff i T 8 i T i=bosons i=fermions ∞ X X g k2 1 p+ = d3~k which allows to write the entropy using just one overall (2π)3 ~ E−µ 3E(k) exp( T )+1 temperature (usually the photon temperature).