Dark Matter Halo Profile Research Articles

The SLUGGS Survey: The Inner Dark Matter Density Slope of the Massive Elliptical Galaxy NGC 1407

Abstract We investigate the dark matter density profile of the massive elliptical galaxy, NGC 1407, by constructing spherically symmetric Jeans models of its field star and globular cluster systems. Two major challenges in such models are the degeneracy between the stellar mass and the dark matter halo profiles, and the degeneracy between the orbital anisotropy of the tracer population and the total mass causing the observed motions. We address the first issue by using new measurements of the mass-to-light ratio profile from stellar population constraints that include a radially varying initial mass function. To mitigate the mass–anisotropy degeneracy, we make use of multiple kinematic tracers, including two subpopulations of globular clusters in addition to the galaxy’s field stars. We create a hierarchical Bayesian model that addresses several often-neglected systematic uncertainties, such as the statistical weight given to various data sets and the adopted distance. After sampling the posterior probability distribution with a Markov chain Monte Carlo method, we find evidence for a central cusp with a log slope of , with the quantified systematic uncertainty dominated by choice of anisotropy profile. This is lower than expected for dark matter halos that have undergone adiabatic contraction, supporting inferences from gravitational lensing that some process has suppressed the steepening of halos in massive galaxies. We also confirm radially biased orbits for the metal-rich globular clusters and tangentially biased orbits for the metal-poor globular clusters, which remains a puzzling finding for an accretion-dominated halo.

The impact of baryonic discs on the shapes and profiles of self-interacting dark matter haloes

We employ isolated N-body simulations to study the response of self-interacting dark matter (SIDM) halos in the presence of the baryonic potentials. Dark matter self-interactions lead to kinematic thermalization in the inner halo, resulting in a tight correlation between the dark matter and baryon distributions. A deep baryonic potential shortens the phase of SIDM core expansion and triggers core contraction. This effect can be further enhanced by a large self-scattering cross section. We find the final SIDM density profile is sensitive to the baryonic concentration and the strength of dark matter self-interactions. Assuming a spherical initial halo, we also study evolution of the SIDM halo shape together with the density profile. The halo shape at later epochs deviates from spherical symmetry due to the influence of the non-spherical disc potential, and its significance depends on the baryonic contribution to the total gravitational potential, relative to the dark matter one. In addition, we construct a multi-component model for the Milky Way, including an SIDM halo, a stellar disc and a bulge, and show it is consistent with observations from stellar kinematics and streams.

Halo mass and weak galaxy-galaxy lensing profiles in rescaled cosmological N-body simulations

We investigate 3D density and weak lensing profiles of dark matter haloes predicted by a cosmology-rescaling algorithm for $N$-body simulations. We extend the rescaling method of Angulo & White (2010) and Angulo & Hilbert (2015) to improve its performance on intra-halo scales by using models for the concentration-mass-redshift relation based on excursion set theory. The accuracy of the method is tested with numerical simulations carried out with different cosmological parameters. We find that predictions for median density profiles are more accurate than $\sim 5\,\%$ for haloes with masses of $10^{12.0} - 10^{14.5} h^{-1}\,M_{\odot}$ for radii $0.05 < r/r_{200\text{m}} < 0.5$, and for cosmologies with $\Omega_\text{m} \in [0.15,\,0.40]$ and $\sigma_8 \in [0.6,\,1.0]$. For larger radii, $0.5 < r/r_{200\text{m}} < 5$, the accuracy degrades to $\sim20\,\%$, due to inaccurate modelling of the cosmological and redshift dependence of the splashback radius. For changes in cosmology allowed by current data, the residuals decrease to $\lesssim2\,\%$ up to scales twice the virial radius. We illustrate the usefulness of the method by estimating the mean halo mass of a mock galaxy group sample. We find that the algorithm's accuracy is sufficient for current data. Improvements in the algorithm, particularly in the modelling of baryons, are likely required for interpreting future (dark energy task force stage IV) experiments.

What the Milky Way’s dwarfs tell us about the Galactic Center extended gamma-ray excess

The Milky Way's Galactic Center harbors a gamma-ray excess that is a candidate signal of annihilating dark matter. Dwarf galaxies remain predominantly dark in their expected commensurate emission. In this work we quantify the degree of consistency between these two observations through a joint likelihood analysis. In doing so we incorporate Milky Way dark matter halo profile uncertainties, as well as an accounting of diffuse gamma-ray emission uncertainties in dark matter annihilation models for the Galactic Center extended gamma-ray excess (GCE) detected by the Fermi Gamma-Ray Space Telescope. The preferred range of annihilation rates and masses expands when including these unknowns. Even so, using two recent determinations of the Milky Way halo's local density leaves the GCE preferred region of single-channel dark matter annihilation models to be in strong tension with annihilation searches in combined dwarf galaxy analyses. A third, higher Milky Way density determination, alleviates this tension. Our joint likelihood analysis allows us to quantify this inconsistency. We provide a set of tools for testing dark matter annihilation models' consistency within this combined data set. As an example, we test a representative inverse Compton sourced self-interacting dark matter model, which is consistent with both the GCE and dwarfs.

Robustness of dark matter constraints and interplay with collider searches for New Physics

We study the implications of dark matter searches, together with collider constraints, on the phenomenological MSSM with neutralino dark matter and focus on the consequences of the related uncertainties in some detail. We consider, inter alia, the latest results from AMS-02, Fermi-LAT and XENON1T. In particular, we examine the impact of the choice of the dark matter halo profile, as well as the propagation model for cosmic rays, for dark matter indirect detection and show that the constraints on the MSSM differ by one to two orders of magnitude depending on the astrophysical hypotheses. On the other hand, our limited knowledge of the local relic density in the vicinity of the Earth and the velocity of Earth in the dark matter halo leads to a factor 3 in the exclusion limits obtained by direct detection experiments. We identified the astrophysical models leading to the most conservative and the most stringent constraints and for each case studied the complementarities with the latest LHC measurements and limits from Higgs, SUSY and monojet searches. We show that combining all data from dark matter searches and colliders, a large fraction of our supersymmetric sample could be probed. Whereas the direct detection constraints are rather robust under the astrophysical assumptions, the uncertainties related to indirect detection can have an important impact on the number of the excluded points.

Search for neutrinos from dark matter self-annihilations in the center of the Milky Way with 3 years of IceCube/DeepCore

We present a search for a neutrino signal from dark matter self-annihilations in the Milky Way using the IceCube Neutrino Observatory (IceCube). In 1005 days of data we found no significant excess of neutrinos over the background of neutrinos produced in atmospheric air showers from cosmic ray interactions. We derive upper limits on the velocity averaged product of the dark matter self-annihilation cross section and the relative velocity of the dark matter particles langle sigma _{text {A}}vrangle . Upper limits are set for dark matter particle candidate masses ranging from 10 GeV up to 1 TeV while considering annihilation through multiple channels. This work sets the most stringent limit on a neutrino signal from dark matter with mass between 10 and 100 GeV, with a limit of 1.18cdot 10^{-23}text { cm}^3text {s}^{-1} for 100 GeV dark matter particles self-annihilating via tau ^+tau ^- to neutrinos (assuming the Navarro–Frenk–White dark matter halo profile).

Density profile of dark matter haloes and galaxies in the horizon–agn simulation: the impact of AGN feedback

Using a suite of three large cosmological hydrodynamical simulations, Horizon-AGN, Horizon-noAGN (no AGN feedback) and Horizon-DM (no baryons), we investigate how a typical sub-grid model for AGN feedback affects the evolution of the inner density profiles of massive dark matter haloes and galaxies. Based on direct object-to-object comparisons, we find that the integrated inner mass and density slope differences between objects formed in these three simulations (hereafter, H_AGN, H_noAGN and H_DM) significantly evolve with time. More specifically, at high redshift (z~5), the mean central density profiles of H_AGN and H_noAGN dark matter haloes tend to be much steeper than their H_DM counterparts owing to the rapidly growing baryonic component and ensuing adiabatic contraction. By z~1.5, these mean halo density profiles in H_AGN have flattened, pummelled by powerful AGN activity ("quasar mode"): the integrated inner mass difference gaps with H_noAGN haloes have widened, and those with H_DM haloes have narrowed. Fast forward 9.5 billion years, down to z=0, and the trend reverses: H_AGN halo mean density profiles drift back to a more cusped shape as AGN feedback efficiency dwindles ("radio mode"), and the gaps in integrated central mass difference with H_noAGN and H_DM close and broaden respectively. On the galaxy side, the story differs noticeably. Averaged stellar profile central densities and inner slopes are monotonically reduced by AGN activity as a function of cosmic time, resulting in better agreement with local observations. As both dark matter and stellar inner density profiles respond quite sensitively to the presence of a central AGN, there is hope that future observational determinations of these quantities can be used constrain AGN feedback models.

Unbiased constraints on ultralight axion mass from dwarf spheroidal galaxies

It has been suggested that the internal dynamics of dwarf spheroidal galaxies (dSphs) can be used to test whether or not ultralight axions with $m_a\sim 10^{-22}\text{eV}$ are a preferred dark matter candidate. However, comparisons to theoretical predictions tend to be inconclusive for the simple reason that while most cosmological models consider only dark matter, one observes only baryons. Here we use realistic kinematic mock data catalogs of Milky Way dSph's to show that the "mass-anisotropy degeneracy" in the Jeans equations leads to biased bounds on the axion mass in galaxies with unknown dark matter halo profiles. In galaxies with multiple chemodynamical components this bias can be partly removed by modelling the mass enclosed within each subpopulation. However, analysis of the mock data reveals that the least-biased constraints on the axion mass result from fitting the luminosity-averaged velocity dispersion of the individual chemodynamical components directly. Applying our analysis to two dSph's with reported stellar subcomponents, Fornax and Sculptor, and assuming that the halo profile has not been acted on by baryons, yields core radii $r_{c}>1.5$ kpc and $r_c> 1.2$ kpc respectively, and $m_a<0.4\times 10^{-22}\text{eV}$ at 97.5\% confidence. These bounds are in tension with the number of observed satellites derived from simple (but conservative) estimates of the subhalo mass function in Milky Way-like galaxies. We discuss how baryonic feedback might affect our results, and the impact of such a small axion mass on the growth of structures in the Universe.

Group galaxy number density profiles far out: Is the ‘one-halo’ term NFW out to &amp;gt;10 virial radii?

Abstract While the density profiles (DPs) of Lambda cold dark matter haloes obey the Navarro, Frenk &amp; White (NFW) law out to roughly one virial radius, rvir, the structure of their outer parts is still poorly understood, because the one-halo term describing the halo itself is dominated by the two-halo term representing the other haloes picked up. Using a semi-analytical model, we measure the real-space one-halo number DP of groups out to 20rvir by assigning each galaxy to its nearest group above mass Ma, in units of the group rvir. If Ma is small (large), the outer DP of groups falls rapidly (slowly). We find that there is an optimal Ma for which the stacked DP resembles the NFW model to 0.1 dex accuracy out to 13 virial radii. We find similar long-range NFW surface DPs (out to 10rvir) in the Sloan Digital Sky Survey observations using a galaxy assignment scheme that combines the non-linear virialized regions of groups with their linear outer parts. The optimal Ma scales as the minimum mass of the groups that are stacked to a power 0.25–0.3. The NFW model thus does not solely originate from violent relaxation. Moreover, populating haloes with galaxies using halo occupation distribution models must proceed out to much larger radii than usually done.

Discovery of a new extragalactic population of energetic particles

We report the discovery of a statistically significant hardening in the Fermi-LAT $\gamma$-ray spectrum of Centaurus A's core, with the spectral index hardening from $\Gamma_{1}=2.73 \pm 0.02$ to $\Gamma_{1}=2.29 \pm 0.07$ at a break energy of ($2.6 \pm 0.3$) GeV. Using a likelihood analysis, we find no evidence for flux variability in Cen A's core lightcurve above or below the spectral break when considering the entire 8 year period. Interestingly, however, the first $\sim3.5$ years of the low energy lightcurve shows evidence of flux variability at the $\sim3.5 \sigma$ confidence level. To understand the origin of this spectral break, we assume that the low energy component below the break feature originates from leptons in Centaurus A's radio jet and we investigate the possibility that the high energy component above the spectral break is due to an additional source of very high energy particles near the core of Cen A. We show for the first time that the observed $\gamma$-ray spectrum of an Active Galactic Nucleus is compatible with either a very large localized enhancement (referred to as a spike) in the dark matter halo profile or a population of millisecond pulsars. Our work constitutes the first robust indication that new $\gamma$-ray production mechanisms can explain the emission from active galaxies and could provide tantalizing first evidence for the clustering of heavy dark matter particles around black holes.

H.E.S.S. Limits on Linelike Dark Matter Signatures in the 100GeV to 2TeV Energy Range Close to the Galactic Center.

No significant excess associated with dark matter annihilations was found in the energy range of 100GeV to 2TeV and upper limits on the gamma-ray flux and the velocity weighted annihilation cross section are derived adopting an Einasto dark matter halo profile. Expected limits for present and future large statistics H.E.S.S. observations are also given.

All-flavour search for neutrinos from dark matter annihilations in the Milky Way with IceCube/DeepCore

We present the first IceCube search for a signal of dark matter annihilations in the Milky Way using all-flavour neutrino-induced particle cascades. The analysis focuses on the DeepCore sub-detector of IceCube, and uses the surrounding IceCube strings as a veto region in order to select starting events in the DeepCore volume. We use 329 live-days of data from IceCube operating in its 86-string configuration during 2011-2012. No neutrino excess is found, the final result being compatible with the background-only hypothesis. From this null result, we derive upper limits on the velocity-averaged self-annihilation cross-section, < \sigma_A v >, for dark matter candidate masses ranging from 30 GeV up to 10 TeV, assuming both a cuspy and a flat-cored dark matter halo profile. For dark matter masses between 200 GeV and 10 TeV, the results improve on all previous IceCube results on < \sigma_A v >, reaching a level of 10^{-23} cm^3 s^-1, depending on the annihilation channel assumed, for a cusped NFW profile. The analysis demonstrates that all-flavour searches are competitive with muon channel searches despite the intrinsically worse angular resolution of cascades compared to muon tracks in IceCube.

The Energetic Particle Population in Centaurus A

AbstractWe report a significant hardening of theFermi-LAT gamma-ray spectrum from the core of Cen A at E &gt; 2.4 GeV, suggesting there is a source of high energy particles in the core of Cen A which is in addition to the jet component. We show that the observed gamma-ray spectrum is compatible with either a spike in the dark matter halo profile or a population of millisecond pulsars. This work gives a strong indication of new gamma-ray production mechanisms in active galactic nuclei and could even provide evidence for the clustering of heavy dark matter particles around black holes.

Dissipative dark matter and the rotation curves of dwarf galaxies

There is ample evidence from rotation curves that dark matter halos around disk galaxies have nontrivial dynamics. Of particular significance are: a) the cored dark matter profile of disk galaxies, b) correlations of the shape of rotation curves with baryonic properties, and c) Tully-Fisher relations. Dark matter halos around disk galaxies may have nontrivial dynamics if dark matter is strongly self interacting and dissipative. Multicomponent hidden sector dark matter featuring a massless `dark photon' (from an unbroken dark U(1) gauge interaction) which kinetically mixes with the ordinary photon provides a concrete example of such dark matter. The kinetic mixing interaction facilitates halo heating by enabling ordinary supernovae to be a source of these `dark photons'. Dark matter halos can expand and contract in response to the heating and cooling processes, but for a sufficiently isolated halo could have evolved to a steady state or `equilibrium' configuration where heating and cooling rates locally balance. This dynamics allows the dark matter density profile to be related to the distribution of ordinary supernovae in the disk of a given galaxy. In a previous paper a simple and predictive formula was derived encoding this relation. Here we improve on previous work by modelling the supernovae distribution via the measured UV and Hα fluxes, and compare the resulting dark matter halo profiles with the rotation curve data for each dwarf galaxy in the LITTLE THINGS sample. The dissipative dark matter concept is further developed and some conclusions drawn.

NIHAO IX: the role of gas inflows and outflows in driving the contraction and expansion of cold dark matter haloes

We use ∼100 cosmological galaxy formation ‘zoom-in’ simulations using the smoothed particle hydrodynamics code gasoline to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M200 ∼ 1010 M⊙) to Milky Way (M200 ∼ 1012 M⊙) masses. Our simulations exhibit a wide range of halo responses, primarily varying with mass, from expansion to contraction, with up to factor ∼10 changes in the enclosed dark matter mass at 1 per cent of the virial radius. Confirming previous studies, the halo response is correlated with the integrated efficiency of star formation: ϵSF ≡ (Mstar/M200)/(Ωb/Ωm). In addition, we report a new correlation with the compactness of the stellar system: ϵR ≡ r1/2/R200. We provide an analytic formula depending on ϵSF and ϵR for the response of cold dark matter haloes to baryonic processes. An observationally testable prediction is that, at fixed mass, larger galaxies experience more halo expansion, while the smaller galaxies more halo contraction. This diversity of dark halo response is captured by a toy model consisting of cycles of adiabatic inflow (causing contraction) and impulsive gas outflow (causing expansion). For net outflow, or equal inflow and outflow fractions, f, the overall effect is expansion, with more expansion with larger f. For net inflow, contraction occurs for small f (large radii), while expansion occurs for large f (small radii), recovering the phenomenology seen in our simulations. These regularities in the galaxy formation process provide a step towards a fully predictive model for the structure of cold dark matter haloes.

Strong lensing in the inner halo of galaxy clusters

We present an axially symmetric formula to calculate the probability of finding gravitational arcs in galaxy clusters, being induced by their massive dark matter haloes, as a function of clusters redshifts and virial masses. The formula includes the ellipticity of the clusters dark matter potential by using a pseudo-elliptical approximation. The probabilities are calculated and compared for two dark-matter halo profiles, the Navarro, Frenk and White (NFW) and the Non-Singular-Isothermal-Sphere (NSIS). We demonstrate the power of our formulation through a Kolmogorov-Smirnov (KS) test on the strong lensing statistics of an X-ray bright sample of low redshift Abell clusters. This KS test allows to establish limits on the values of the concentration parameter for the NFW profile ($c_\Delta$) and the core radius for the NSIS profile (\rc), which are related to the lowest cluster redshift ($z_{\rm cut}$) where strong arcs can be observed. For NFW dark matter profiles, we infer cluster haloes with concentrations that are consistent to those predicted by $\Lambda$CDM simulations. As for NSIS dark matter profiles, we find only upper limits for the clusters core radii and thus do not rule out a purely SIS model. For alternative mass profiles, our formulation provides constraints through $z_{\rm cut}$ on the parameters that control the concentration of mass in the inner region of the clusters haloes. We find that $z_{\rm cut}$ is expected to lie in the 0.0--0.2 redshift, highlighting the need to include very low-$z$ clusters in samples to study the clusters mass profiles.

Wide binaries in ultrafaint galaxies: a window on to dark matter on the smallest scales

Abstract We carry out controlled N-body simulations that follow the dynamical evolution of binary stars in the dark matter (DM) haloes of ultrafaint dwarf spheroidals (dSphs). We find that wide binaries with semimajor axes a ≳ at tend to be quickly disrupted by the tidal field of the halo. In smooth potentials the truncation scale, at, is mainly governed by (i) the mass enclosed within the dwarf half-light radius (Rh) and (ii) the slope of the DM halo profile at R ≈ Rh, and is largely independent of the initial eccentricity distribution of the binary systems and the anisotropy of the stellar orbits about the galactic potential. For the reported velocity dispersion and half-light radius of Segue I, the closest ultrafaint, our models predict at values that are a factor of 2–3 smaller in cuspy haloes than in any of the cored models considered here. Using mock observations of Segue I we show that measuring the projected two-point correlation function of stellar pairs with sub-arcsecond resolution may provide a useful tool to constrain the amount and distribution of DM in the smallest and most DM-dominated galaxies.

The outer profile of dark matter haloes: an analytical approach

A steepening feature in the outer density profiles of dark matter halos indicating the splashback radius has drawn much attention recently. Possible observational detections have even been made for galaxy clusters. Theoretically, Adhikari et al. have estimated the location of the splashback radius by computing the secondary infall trajectory of a dark matter shell through a growing dark matter halo with an NFW profile. However, since they imposed a shape of the halo profile rather than computing it consistently from the trajectories of the dark matter shells, they could not provide the full shape of the dark matter profile around the splashback radius. We improve on this by extending the self-similar spherical collapse model of Fillmore \& Goldreich to a $\Lambda$CDM universe. This allows us to compute the dark matter halo profile and the trajectories simultaneously from the mass accretion history. Our results on the splashback location agree qualitatively with Adhikari et al. but with small quantitative differences at large mass accretion rates. We present new fitting formulae for the splashback radius $R_{\rm sp}$ in various forms, including the ratios of $R_{\rm sp} / R_{\rm 200c}$ and $R_{\rm sp} / R_{\rm 200m}$. Numerical simulations have made the puzzling discovery that the splashback radius scales well with $R_{\rm 200m}$ but not with $R_{\rm 200c}$. We trace the origin of this to be the correlated increase of $\Omega_{\rm m}$ and the average halo mass accretion rate with an increasing redshift.

Bright gamma-ray Galactic Center excess and dark dwarfs: Strong tension for dark matter annihilation despite Milky Way halo profile and diffuse emission uncertainties

We incorporate Milky Way dark matter halo profile uncertainties, as well as an accounting of diffuse gamma-ray emission uncertainties in dark matter annihilation models for the Galactic Center Extended gamma-ray excess (GCE) detected by the Fermi Gamma Ray Space Telescope. The range of particle annihilation rate and masses expand when including these unknowns. However, two of the most precise empirical determinations of the Milky Way halo's local density and density profile leave the signal region to be in considerable tension with dark matter annihilation searches from combined dwarf galaxy analyses for single-channel dark matter annihilation models. The GCE and dwarf tension can be alleviated if: one, the halo is very highly concentrated or strongly contracted; two, the dark matter annihilation signal differentiates between dwarfs and the GC; or, three, local stellar density measures are found to be significantly lower, like that from recent stellar counts, increasing the local dark matter density.

Studying dark matter haloes with Schwarzschild modeling

AbstractWe present a study exploring the reliability of recovering the mass and velocity anisotropy profiles of dark matter haloes using the Schwarzschild method of modeling galaxies.