The recent observation of an X-ray line at an energy of 3.5 keV mainly from galaxy clusters has initiated a discussion about whether we may have seen a possible dark matter signal. If confirmed, this ...signal could stem from a decaying sterile neutrino of a mass of 7.1 keV. Such a particle could make up all the dark matter, but it is not clear how it was produced in the early Universe. In this letter we show that it is possible to discriminate between different production mechanisms with present-day astronomical data. The most stringent constraint comes from the Lyman-α forest and seems to disfavor all but one of the main production mechanisms proposed in the literature, which is the production via decay of heavy scalar singlets. Pinning down the production mechanism will help to decide whether the X-ray signal indeed comprises an indirect detection of dark matter.
Abstract
We present a new power spectrum emulator named EuclidEmulator that estimates the nonlinear correction to the linear dark matter power spectrum depending on the six cosmological parameters ...ωb, ωm, ns, h, $w$0, and σ8. It is constructed using the uncertainty quantification software UQLab using a spectral decomposition method called polynomial chaos expansion. All steps in its construction have been tested and optimized: the large high-resolution N-body simulations carried out with PKDGRAV3 were validated using a simulation from the Euclid Flagship campaign and demonstrated to have converged up to wavenumbers $k\approx 5\, h\, {\rm Mpc}^{-1}$ for redshifts $z$ ≤ 5. The emulator is based on 100 input cosmologies simulated in boxes of (1250 Mpc/h)3 using 20483 particles. We show that by creating mock emulators it is possible to successfully predict and optimize the performance of the final emulator prior to performing any N-body simulations. The absolute accuracy of the final nonlinear power spectrum is as good as one obtained with N-body simulations, conservatively, ${\sim } 1$ per cent for $k\lesssim 1\, h\, {\rm Mpc}^{-1}$ and $z$ ≲ 1. This enables efficient forward modelling in the nonlinear regime, allowing for estimation of cosmological parameters using Markov Chain Monte Carlo methods. EuclidEmulator has been compared to HALOFIT, CosmicEmu, and NGenHalofit, and shown to be more accurate than these other approaches. This work paves a new way for optimal construction of future emulators that also consider other cosmological observables, use higher resolution input simulations, and investigate higher dimensional cosmological parameter spaces.
All commonly considered dark matter scenarios are based on hypothetical particles with small but non-zero thermal velocities and tiny interaction cross-sections. A generic consequence of these ...attributes is the suppression of small-scale matter perturbations either due to free-streaming or due to interactions with the primordial plasma. The suppression scale can vary over many orders of magnitude depending on particle candidate and production mechanism in the early Universe. While non-linear structure formation has been explored in great detail well above the suppression scale, the range around suppressed perturbations is still poorly understood. In this paper, we study structure formation in the regime of suppressed perturbations using both analytical techniques and numerical simulations. We develop simple and theoretically motivated recipes for the halo mass function, the expected number of satellites, and the halo concentrations, which are designed to work for power spectra with suppression at arbitrary scale and of arbitrary shape. As case studies, we explore warm and mixed dark matter scenarios where effects are most distinctive. Additionally, we examine the standard dark matter scenario based on weakly interacting massive particles (WIMP) and compare it to pure cold dark matter with zero primordial temperature. We find that our analytically motivated recipes are in good agreement with simulations for all investigated dark matter scenarios, and we therefore conclude that they can be used for generic cases with arbitrarily suppressed small-scale perturbations.
The dark energy dominated warm dark matter (WDM) model is a promising alternative cosmological scenario. We explore large-scale structure formation in this paradigm. We do this in two different ways: ...with the halo model approach and with the help of an ensemble of high-resolution N-body simulations. Combining these quasi-independent approaches leads to a physical understanding of the important processes which shape the formation of structures. We take a detailed look at the halo mass function, the concentrations and the linear halo bias of WDM. In all cases we find interesting deviations with respect to cold dark matter (CDM). In particular, the concentration-mass relation displays a turnover for group scale dark matter haloes, for the case of WDM particles with masses of the order of m
WDM∼ 0.25 keV. This may be interpreted as a hint for top-down structure formation on small scales. We implement our results into the halo model and find much better agreement with simulations. On small scales, the WDM halo model now performs as well as its CDM counterpart.
Abstract
The dearth of dwarf galaxies in the local Universe is hard to reconcile with the large number of low-mass haloes expected within the concordance Λ cold dark matter (ΛCDM) paradigm. In this ...paper, we perform a systematic evaluation of the uncertainties affecting the measurement of dark matter halo abundance using galaxy kinematics. Using a large sample of dwarf galaxies with spatially resolved kinematics, we derive a correction to obtain the abundance of galaxies as a function of maximum circular velocity – a direct probe of halo mass – from the line-of-sight velocity function in the Local Volume. This method provides a direct means of comparing the predictions of theoretical models and simulations (including non-standard cosmologies and novel galaxy formation physics) to the observational constraints. The new ‘galactic Vmax’ function is steeper than the line-of-sight velocity function but still shallower than the theoretical CDM expectation, implying that unaccounted baryonic physics may be necessary to reduce the predicted abundance of galaxies. Using the galactic Vmax function, we investigate the theoretical effects of feedback-powered outflows and photoevaporation of gas due to reionization. At the 3σ confidence level, we find that feedback and reionization are not effective enough to reconcile the disagreement. In the case of maximum baryonic effects, the theoretical prediction still deviates significantly from the observations for Vmax < 60 km s−1. CDM predicts at least 1.8 times more galaxies with Vmax = 50 km s−1 and 2.5 times more than observed at 30 km s−1. Recent hydrodynamic simulations seem to resolve the discrepancy but disagree with the properties of observed galaxies with spatially resolved kinematics. This abundance problem might point to the need to modify cosmological predictions at small scales.
A halo merger tree forms the essential backbone of a semi-analytic model for galaxy formation and evolution. Recent studies have pointed out that extracting merger trees from numerical simulations of ...structure formation is non-trivial; different tree building algorithms can give differing merger histories. These differences should be carefully understood before merger trees are used as input for models of galaxy formation. We investigate the impact of different halo merger trees on a semi-analytic model. We find that the z = 0 galaxy properties in our model show differences between trees when using a common parameter set. The star formation history of the universe and the properties of satellite galaxies can show marked differences between trees with different construction methods. Independently calibrating the semi-analytic model for each tree can reduce the discrepancies between the z = 0 global galaxy properties, at the cost of increasing the differences in the evolutionary histories of galaxies. Furthermore, the underlying physics implied can vary, resulting in key quantities such as the supernova feedback efficiency differing by factors of 2. Such a change alters the regimes where star formation is primarily suppressed by supernovae. Therefore, halo merger trees extracted from a common halo catalogue using different, but reliable, algorithms can result in a difference in the semi-analytic model. Given the uncertainties in galaxy formation physics, however, these differences may not necessarily be viewed as significant.
From observations at low and high redshifts, it is well known that the bulk of dark matter (DM) has to be stable or at least very long-lived. However, the possibility that a small fraction of DM is ...unstable or that all DM decays with a half-life time (
τ
) significantly longer than the age of the Universe is not ruled out. One-body decaying dark matter (DDM) consists of a minimal extension to the ΛCDM model. It causes a modification of the cosmic growth history as well as a suppression of the small-scale clustering signal, providing interesting consequences regarding the
S
8
tension, which is the observed difference in the clustering amplitude between weak-lensing (WL) and cosmic microwave background (CMB) observations. In this paper, we investigate models in which a fraction or all DM decays into radiation, focusing on the long-lived regime, that is,
τ
≳
H
0
−1
(
H
0
−1
being the Hubble time). We used WL data from the Kilo-Degree Survey (KiDS) and CMB data from
Planck
. First, we confirm that this DDM model cannot alleviate the
S
8
difference. We then show that the most constraining power for DM decay does not come from the nonlinear WL data, but from CMB via the integrated Sachs-Wolfe effect. From the CMB data alone, we obtain constraints of
τ
≥ 288 Gyr if all DM is assumed to be unstable, and we show that a maximum fraction of
f
= 0.07 is allowed to decay assuming the half-life time to be comparable to (or shorter than) one Hubble time. The constraints from the KiDS-1000 WL data are significantly weaker,
τ
≥ 60 Gyr and
f
< 0.34. Combining the CMB and WL data does not yield tighter constraints than the CMB alone, except for short half-life times, for which the maximum allowed fraction becomes
f
= 0.03. All limits are provided at the 95% confidence level.
ABSTRACT
Many non-minimal dark matter scenarios lead to oscillatory features in the matter power spectrum induced by interactions either within the dark sector or with particles from the standard ...model. Observing such dark acoustic oscillations would therefore be a major step towards understanding dark matter. We investigate what happens to oscillatory features during the process of non-linear structure formation. We show that at the level of the power spectrum, oscillations are smoothed out by non-linear mode coupling, gradually disappearing towards lower redshifts. In the halo mass function, however, the oscillatory features remain visible until the present epoch. As a consequence, dark acoustic oscillations could be detectable in observations that are either based on the halo mass function or on the high-redshift power spectrum. We investigate the effect of such features on different observables, namely the cluster mass function, the stellar-to-halo mass relation, and the Lyman α flux power spectrum. We find that oscillatory features remain visible in all of these observables, but they are very extended and of low amplitude, making it challenging to detect them as distinct features in the data.
Abstract
The number density of field galaxies per rotation velocity, referred to as the velocity function, is an intriguing statistical measure probing the smallest scales of structure formation. In ...this paper we point out that the velocity function is sensitive to small shifts in key cosmological parameters such as the amplitude of primordial perturbations (σ8) or the total matter density (Ωm). Using current data and applying conservative assumptions about baryonic effects, we show that the observed velocity function of the Local Volume favours cosmologies in tension with the measurements from Planck but in agreement with the latest findings from weak lensing surveys. While the current systematics regarding the relation between observed and true rotation velocities are potentially important, upcoming data from H i surveys as well as new insights from hydrodynamical simulations will dramatically improve the situation in the near future.