Dark Matter Density Profile Research Articles

Closing the window on high-mass strongly interacting dark matter

New and improved constraints are placed on the spin-independent interaction cross section of dark matter with regular matter by refining two methods. First, dark matter–cosmic ray interactions are considered, wherein cosmic ray protons collide with dark matter to contribute to the gamma ray sky. This constraint is developed using the NFW, Moore, and Einasto dark matter density profiles and new data from the Fermi gamma ray space telescope. Second, the Earth capture scenario is considered, wherein particles that are captured self-annihilate at Earth's center, thus adding to its internal heat flow. The constraint presented here is based on analysis of the drift time of dark matter particles through Earth, modeled as a core composed of iron and a mantle composed of oxygen with linear density gradients between layers. An analysis of the cosmic ray constraint (which rules out dark matter–regular matter interaction cross sections greater than its value) shows that it overlaps considerably with the Earth drift time constraint (which rules out cross sections smaller than its value). Most significantly, the new cosmic ray constraint excludes the last remaining open region for strongly interacting dark matter particles up to a mass of about 1017 GeV when combined with other exclusions, closing the window for strongly interacting dark matter.

The impact of the phase-space density on the indirect detection of dark matter

We study the indirect detection of dark matter when the local dark matter velocity distribution depends upon position, as expected for the Milky Way and its dwarf spheroidal satellites, and the annihilation cross-section is not purely s-wave. Using a phase-space distribution consistent with the dark matter density profile, we present estimates of cosmic and gamma-ray fluxes from dark matter annihilations. The expectations for the indirect detection of dark matter can differ significantly from the usual calculation that assumes that the velocity of the dark matter particles follows a Maxwell-Boltzmann distribution.

Clustering analysis of the morphology of the 130 GeV gamma-ray feature

Recent observations indicating the existence of a monochromatic gamma-ray line with energy ~130 GeV in the Fermi-LAT data have attracted great interest due to the possibility that the line feature stems from the annihilation of dark matter particles. Many studies examining the robustness of the putative line-signal have concentrated on its spectral attributes. Here, we study the morphological features of the gamma-ray line photons, which can be used to differentiate a putative dark matter signal from astrophysical backgrounds or instrumental artifacts. Photons stemming from dark matter annihilation will produce events tracing a specific morphology, with a statistical clustering that can be calculated based on models of the dark matter density profile in the inner Galaxy. We apply the DBSCAN clustering algorithm to Fermi gamma-ray data, and show that we can rule out the possibility that 1 (2, 4) or fewer point-like sources produce the observed morphology for the line photons at a 99% (95%, 90%) confidence level. Our study strongly disfavors the main astrophysical background envisioned to produce a line feature at energies above 100 GeV: cold pulsar winds. It is highly unlikely that 4 or more such objects have exactly the same monochromatic cosmic-ray energy needed to produce a gamma-ray line, to within instrumental energy resolution. Furthermore, we show that the larger photon statistics expected with Air Cherenkov Telescopes such as H.E.S.S.-II will allow for extraordinarily stringent morphological tests of the origin of the "line photons".

Orbital anisotropy in cosmological haloes revisited

The velocity anisotropy of particles inside dark matter (DM) haloes is an important physical quantity, which is required for the accurate modelling of mass profiles of galaxies and clusters of galaxies. It is typically measured using the ratio of the radial-to-tangential velocity dispersions at a given distance from the halo centre. However, this measure is insufficient to describe the dynamics of realistic haloes, which are typically quite elongated. Studying the velocity distribution in massive DM haloes in cosmological simulations, we find that in the inner parts of the haloes the local velocity ellipsoids are strongly aligned with the major axis of the halo, the alignment being stronger for more relaxed haloes. In the outer regions of the haloes, the alignment becomes gradually weaker and the orientation is more random. These two distinct regions of different degree of the alignment coincide with two characteristic regimes of the DM density profile: shallower and steeper than \rho r^{-2}. This alignment of the local velocity ellipsoids requires reinterpretation of features found in measurements based on the spherically averaged ratio of the radial-to-tangential velocity dispersions. In particular, we show that the velocity distribution in the central halo regions is highly anisotropic. For cluster-size haloes with mass 10^{14}-10^{15} h^-1 Msun, the velocity anisotropy along the major axis is nearly independent of radius and is equal to \beta=1-\sigma^2_{perp}/\sigma^2_{radial}=0.4, which is significantly larger than the previously estimated spherically averaged velocity anisotropy. The alignment of density and velocity anisotropies, and the radial trends may also have some implications for the mass modelling based on kinematical data of such objects as galaxy clusters or dwarf spheroidals, where the orbital anisotropy is a key element in an unbiased mass inference.

Constraints on dark matter annihilation by radio observations of M31

We used radio observations of the neighbour galaxy M31 in order to put constraints on dark matter particle mass and annihilation cross section. Dark matter annihilation in M31 halo produces highly energetic leptons, which emit synchrotron radiation on radio frequencies in the galactic magnetic field. We predicted expected radio fluxes for the two annihilation channels: \chi\chi -> bb* and \chi\chi -> \tau^+\tau^-. We then compared them with available data on the central radio emission of M31 as observed by four radio surveys: VLSS (74 MHz), WENSS (325 MHz), NVSS (1400 MHz) and GB6 (4850 MHz). Assuming a standard NFW dark matter density profile and a conservative magnetic field distribution inside the Andromeda galaxy, we find that the thermal relic annihilation cross section <\sigma v> = 3*10^{-26} cm^3/s or higher are only allowed for WIMP masses greater than 100 GeV and 55 GeV for annihilation into bb* and \tau^+\tau^- respectively. Taking into account potential uncertainties in the distributions of DM density and magnetic field, the mentioned WIMP limiting masses can be as low as 23 GeV for both channels, and as high as 280 and 130 GeV for annihilation into bb* and \tau^+\tau^- respectively. These mass values exceed the best up-to-day known constraints from Fermi gamma observations: 40 GeV and 19 GeV respectively [A.Geringer-Sameth and S.M.Koushiappas, Phys. Rev. Lett. 107, 241303 (2011)]. Precise measurements of the magnetic field in the relevant region and better reconstruction of the DM density profile of M31 will be able to reduce the uncertainties of our exclusion limits.

Scaling relations of the slightly self-interacting cold dark matter in galaxies and clusters

Recent observations in galaxies and clusters indicate dark matter density profiles exhibit core-like structures which contradict to the numerical simulation results of collisionless cold dark matter. The idea of self-interacting cold dark matter (SICDM) has been invoked to solve the discrepancies between the observations and numerical simulations. In this article, I derive some important scaling relations in galaxies and clusters by using the long-range SICDM model. These scaling relations give good agreements with the empirical fittings from observational data in galaxies and clusters if the dark matter particles are only slightly self-interacting. Also, there may exist a universal critical optical depth $\tau_c$ that characterizes the core-like structures. These results generally support the idea of SICDM to tackle the long-lasting dark matter problem.

OBSERVATIONAL EVIDENCE FOR DARK MATTER INTERACTING THROUGH A YUKAWA POTENTIAL

Recent observations in galaxies and clusters indicate dark matter density profiles exhibit core-like structures which contradict to the numerical simulation results of collisionless cold dark matter. On the other hand, it has been shown that cold dark matter particles interacting through a Yukawa potential could naturally explain the cores in dwarf galaxies. In this article, I use the Yukawa Potential interacting dark matter model to derive two simple scaling relations on the galactic and cluster scales respectively, which give excellent agreements with observations. Also, in our model, the masses of the force carrier and dark matter particle can be constrained by the observational data.

Gamma rays from cosmic-ray proton scattering in AGN jets: the intra-cluster gas vastly outshines dark matter

Active Galactic Nuclei (AGN) host powerful jets containing high-energy electrons and protons. The astrophysical environment where AGNs and their jets are found is characterized by large concentrations of both dark matter (DM) and intra-cluster medium (ICM) gas. As the high-energy jet particles transverse the DM and the ICM, elastic and inelastic scattering processes generically lead to the production of final-state photons. As first envisioned by Bloom and Wells (1998), and as more recently pointed out by us and others, the scattering of electrons off of DM could lead to a potentially detectable gamma-ray signal, with the parton-level contribution from protons offering dimmer perspectives. Recently, Chang et al. argued that taking into account photons from hadronization and showering, the actual photon flux is substantially increased. Here, we point out that the proton-jets have to be highly collimated, contrary to what predicted by simple blob-geometry jet-models sometimes employed in these studies, otherwise they would produce a very large flux of photons from inelaystic collisions with ICM nucleons, which would outshine by many orders of magnitude the signal from DM, for almost any reasonable ICM and DM density profiles.

THE DARK HALO—SPHEROID CONSPIRACY AND THE ORIGIN OF ELLIPTICAL GALAXIES

Dynamical modeling and strong lensing data indicate that the total density profiles of early-type galaxies are close to isothermal, i.e., rho_tot ~ r^gamma with gamma approx -2. To understand the origin of this universal slope we study a set of simulated spheroids formed in isolated binary mergers as well as the formation within the cosmological framework. The total stellar plus dark matter density profiles can always be described by a power law with an index of gamma approx -2.1 with a tendency toward steeper slopes for more compact, lower-mass ellipticals. In the binary mergers the amount of gas involved in the merger determines the precise steepness of the slope. This agrees with results from the cosmological simulations where ellipticals with steeper slopes have a higher fraction of stars formed in situ. Each gas-poor merger event evolves the slope toward gamma ~ -2, once this slope is reached further merger events do not change it anymore. All our ellipticals have flat intrinsic combined stellar and dark matter velocity dispersion profiles. We conclude that flat velocity dispersion profiles and total density distributions with a slope of gamma ~ -2 for the combined system of stars and dark matter act as a natural attractor. The variety of complex formation histories as present in cosmological simulations, including major as well as minor merger events, is essential to generate the full range of observed density slopes seen for present-day elliptical galaxies.

THE DENSITY PROFILES OF MASSIVE, RELAXED GALAXY CLUSTERS. II. SEPARATING LUMINOUS AND DARK MATTER IN CLUSTER CORES

We present stellar and dark matter (DM) density profiles for a sample of seven massive, relaxed galaxy clusters derived from strong and weak gravitational lensing and resolved stellar kinematic observations within the centrally-located brightest cluster galaxies (BCGs). In Paper I of the series, we demonstrated that the total density profile derived from these data, which span 3 decades in radius, is consistent with numerical DM-only simulations at radii >~ 5-10 kpc, despite the significant contribution of stellar material in the core. Here we decompose the inner mass profiles of these clusters into stellar and dark components. Parametrizing the DM density profile as a power law rho_DM ~ r^{-\beta} on small scales, we find a mean slope <\beta> = 0.50 +- 0.10 (random) +0.14-0.13 (systematic). Alternatively, cored Navarro-Frenk-White (NFW) profiles with <log r_core/kpc> = 1.14 +- 0.13 (random) +0.14-0.22 (systematic) provide an equally good description. These density profiles are significantly shallower than canonical NFW models at radii <~ 30 kpc, comparable to the effective radii of the BCGs. The inner DM profile is correlated with the distribution of stars in the BCG, suggesting a connection between the inner halo and the assembly of stars in the central galaxy. The stellar mass-to-light ratio inferred from lensing and stellar dynamics is consistent with that inferred using stellar population synthesis models if a Salpeter initial mass function is adopted. We compare these results to theories describing the interaction between baryons and DM in cluster cores, including adiabatic contraction models and the possible effects of galaxy mergers and active galactic nucleus feedback, and evaluate possible signatures of alternative DM candidates.

TIDAL STIRRING OF DISKY DWARFS WITH SHALLOW DARK MATTER DENSITY PROFILES: ENHANCED TRANSFORMATION INTO DWARF SPHEROIDALS

(Abridged) The origin of dSphs in the Local Group (LG) remains an enigma. The tidal stirring model posits that late-type, rotationally-supported dwarfs resembling present-day dwarf irregular (dIrr) galaxies can transform into dSphs via interactions with Milky Way-sized hosts. Using collisionless N-body simulations, we investigate for the first time how tidal stirring depends on the dark matter (DM) density distribution in the central stellar region of the progenitor disky dwarf. Specifically, we explore various asymptotic inner slopes gamma of the dwarf DM density profiles (rho \propto r^{-gamma} as r -> 0). For a given orbit inside the primary, rotationally-supported dwarfs embedded in DM halos with core-like density distributions (gamma = 0.2) and mild density cusps (gamma = 0.6) demonstrate a substantially enhanced likelihood and efficiency of transformation into dSphs compared to their counterparts with steeper DM density profiles (gamma = 1). Such shallow DM distributions are akin to those of observed dIrrs, highlighting tidal stirring as a plausible model for the LG morphology-density relation. When gamma <1, a single pericentric passage can induce dSph formation and disky dwarfs on low-eccentricity or large-pericenter orbits are able to transform into dSphs; these new results allow the tidal stirring model to explain the existence of virtually all known dSphs across a wide range of distances from their hosts. A subset of rotationally-supported dwarfs with gamma <1 are eventually disrupted by the primary; those that survive as dSphs are generally on orbits that are biased towards lower eccentricities and/or larger pericenters relative to those of typical CDM satellites. The latter could explain the rather peculiar orbits of several classic LG dSphs such as Fornax, Leo I, Tucana, and Cetus.

MEASURING DARK MATTER PROFILES NON-PARAMETRICALLY IN DWARF SPHEROIDALS: AN APPLICATION TO DRACO

We introduce a novel implementation of orbit-based (or Schwarzschild) modeling that allows dark matter density profiles to be calculated non-parametrically in nearby galaxies. Our models require no assumptions to be made about velocity anisotropy or the dark matter profile. The technique can be applied to any dispersion-supported stellar system, and we demonstrate its use by studying the Local Group dwarf spheroidal (dSph) galaxy Draco. We use existing kinematic data at larger radii and also present 12 new radial velocities within the central 13 pc obtained with the VIRUS-W integral field spectrograph on the 2.7m telescope at McDonald Observatory. Our non-parametric Schwarzschild models find strong evidence that the dark matter profile in Draco is cuspy for 20 < r < 700 pc. The profile for r > 20 pc is well-fit by a power law with slope \alpha=-1.0 +/- 0.2, consistent with predictions from Cold Dark Matter (CDM) simulations. Our models confirm that, despite its low baryon content relative to other dSphs, Draco lives in a massive halo.

Dark matter chaos in the Solar system

Abstract We study the capture of Galactic dark matter particles in the Solar system produced by rotation of Jupiter. It is shown that the capture cross-section is much larger than the area of the Jupiter orbit being inversely diverging at small particle energy. We show that the dynamics of captured particles is chaotic and is well described by a simple symplectic dark map. This dark map description allows us to simulate the scattering and dynamics of 1014 dark matter particles during the lifetime of the Solar system and to determine the dark matter density profile as a function of distance from the Sun. The mass of captured dark matter in the radius of the Neptune orbit is estimated to be 2 × 1015g. The radial density of captured dark matter is found to be approximately constant behind the Jupiter orbit being similar to the density profile found in galaxies.

Characterization of dark-matter-induced anisotropies in the diffuse gamma-ray background

The Fermi-LAT collaboration has recently reported the detection of angular power above the photon noise level in the diffuse gamma-ray background between 1 and 50 GeV. Such signal can be used to constrain a possible contribution from Dark-Matter-induced photons. We estimate the intensity and features of the angular power spectrum (APS) of this potential Dark Matter (DM) signal, for both decaying and annihilating DM candidates, by constructing template all-sky gamma-ray maps for the emission produced in the galactic halo and its substructures, as well as in extragalactic (sub)halos. The DM distribution is given by state-of-the-art N-body simulations of cosmic structure formation, namely Millennium-II for extragalactic (sub)halos, and Aquarius for the galactic halo and its subhalos. We use a hybrid method of extrapolation to account for (sub)structures that are below the resolution limit of the simulations, allowing us to estimate the total emission all the way down to the minimal self-bound halo mass. We describe in detail the features appearing in the APS of our template maps and we estimate the effect of various uncertainties such as the value of the minimal halo mass, the fraction of substructures hosted in a halo and the shape of the DM density profile. Our results indicate that the fluctuation APS of the DM-induced emission is of the same order as the Fermi-LAT APS, suggesting that one can constrain this hypothetical emission from the comparison with the measured anisotropy. We also quantify the uncertainties affecting our results, finding "theoretical error bands" spanning more than two orders of magnitude and dominated (for a given particle physics model) by our lack of knowledge of the abundance of low-mass (sub)halos.

Fitting functions for dark matter density profiles

We present a unified parameterization of the fitting functions suitable for density profiles of dark matter haloes or elliptical galaxies. A notable feature is that the classical Einasto profile appears naturally as the continuous limiting case of the cored subfamily amongst the double power-law profiles of Zhao (1996). Based on this, we also argue that there is basically no qualitative difference between halo models well-fitted by the Einasto profile and the standard NFW model. This may even be the case quantitatively unless the resolutions of simulations and the precisions of fittings are sufficiently high to make meaningful distinction possible.

Evolution of the first supernovae in protogalaxies: Dynamics of mixing of heavy elements

The paper considers the evolution of the supernova envelopes produced by Population III stars with masses of $M_*\sim 25-200 M_\odot$ located in non-rotating protogalaxies with masses of $M\sim 10^7 M_\odot$ at redshifts $z=12$, with dark-matter density profiles in the form of modified isothermal spheres. The supernova explosion occurs in the ionization zone formed by a single parent star. The properties of the distribution of heavy elements (metals) produced by the parent star are investigated, as well as the efficiency with which they are mixed with the primordial gas in the supernova envelope. In supernovae with high energies ($E\simgt 5\times 10^{52}$ erg), an appreciable fraction of the gas can be ejected from the protogalaxy, but nearly all the heavy elements remain in the protogalaxy. In explosions with lower energies ($E\simlt 3\times 10^{52}$ erg), essentially no gas and heavy elements are lost from the protogalaxy: during the first one to threemillion years, the gas and heavy elements are actively carried from the central region of the protogalaxy ($r\sim 0.1 r_{vir}$, where $r_{vir}$ is the virial radius of the protogalaxy), but an appreciable fraction of the mass of metals subsequently returns when the hot cavity cools and the envelope collapses. Supernovae with high energies ($E\simgt 5\times 10^{52}$ erg) are characterized by a very low efficiency of mixing of metals; their heavy elements are located in the small volume occupied by the disrupted envelope (in a volume comparable with that of the entire envelope), with most of the metals remaining inside the hot, rarified cavity of the envelope. (abridged)

A MODEL-INDEPENDENT METHOD TO STUDY DARK MATTER INDUCED LEPTONS AND GAMMA RAYS

By using recent data, we directly determine the dark matter (DM)-induced e± spectrum at the source from experimental measurements at the earth, without reference to specific particle physics models. The DM-induced gamma rays emitted via inverse Compton scattering are then obtained in a model-independent way. However, the results depend on the choice of the astrophysical e± background, which is not reliably known. Nevertheless, we calculate, as an illustration, the fluxes of gamma rays from the Fornax cluster in the decaying DM scenario with various astrophysical e± backgrounds. Without any assumptions on details of the DM model, the predictions turn out to be either in disagreement with or only marginally below the upper limits measured recently by the Fermi-LAT Collaboration. In addition, these DM-induced ICS gamma rays in the GeV range are shown to be almost independent of choices of cosmic ray propagation model and of DM density profile, when a given astrophysical e± background is assumed. This provides a strong constraint on decaying DM scenario as the gamma rays may be produced in other processes besides inverse Compton scattering, such as the bremsstrahlung and neutral pion decays.

Semi-empirical catalog of early-type galaxy-halo systems: dark matter density profiles, halo contraction and dark matter annihilation strength

With Sloan Digital Sky Survey galaxy data and halo data from up-to-date N-body simulations within the ΛCDM framework we construct a semi-empirical catalog (SEC) of early-type galaxy-halo systems by making a self-consistent bivariate statistical match of stellar mass (M⋆) and velocity dispersion (σ) with halo virial mass (Mvir) as demonstrated here for the first time. We then assign stellar mass profile and velocity dispersion profile parameters to each system in the SEC using their observed correlations with M⋆ and σ. Simultaneously, we solve for dark matter density profile of each halo using the spherical Jeans equation. The resulting dark matter density profiles deviate in general from the dissipationless profile of Navarro-Frenk-White or Einasto and their mean inner density slope and concentration vary systematically with Mvir. Statistical tests of the distribution of profiles at fixed Mvir rule out the null hypothesis that it follows the distribution predicted by dissipationless N-body simulations for Mvir≲1013.5 – 14.5 M⊙. These dark matter profiles imply that dark matter density is, on average, enhanced significantly in the inner region of halos with Mvir≲1013.5 – 14.5 M⊙ supporting halo contraction. The main characteristics of halo contraction are: (1) the mean dark matter density within the effective radius has increased by a factor varying systematically up to ≈ 3–4 at Mvir = 1012 M⊙, and (2) the inner density slope has a mean of ⟨α⟩ ≈ 1.3 with ρdm(r)∝r−α and a halo-to-halo rms scatter of rms(α) ∼ 0.4–0.5 for 1012 M⊙≲Mvir≲1013 – 14 M⊙ steeper than the NFW profile (α = 1). Based on our results we predict that halos of nearby elliptical and lenticular galaxies can, in principle, be promising targets for γ-ray emission from dark matter annihilation.

Illuminating the 130 GeV gamma line with continuum photons

Abstract There is evidence for a 130 GeV γ-ray line at the Galactic Center in the Fermi Large Area Telescope data. Dark matter candidates that explain this feature should also annihilate to Standard Model particles, resulting in a continuous spectrum of photons. To study this continuum, we analyze the Fermi data down to 5 GeV, restricted to the inner 3◦ of the Galaxy. We place a strong bound on the ratio of continuum photons to monochromatic line photons that is independent of uncertainties in the dark matter density profile. The derived constraints exclude neutralino dark matter as an explanation for the line.

Resolution of dark matter problem in f(T) gravity

In this paper, we attempt to resolve the dark matter problem in f(T) gravity. Specifically, from our model we successfully obtain the flat rotation curves of galaxies containing dark matter. Further, we obtain the density profile of dark matter in galaxies. Comparison of our analytical results shows that our torsion-based toy model for dark matter is in good agreement with empirical data-based models. It shows that we can address the dark matter as an effect of torsion of the space.