BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Aurel Schneider

Aurel Schneider appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAhep-phastro-ph.IM

Catalog footprint

What is connected

29works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2022arXiv

A Full $w$CDM Analysis of KiDS-1000 Weak Lensing Maps using Deep Learning

We present a full forward-modeled $w$CDM analysis of the KiDS-1000 weak lensing maps using graph-convolutional neural networks (GCNN). Utilizing the $\texttt{CosmoGrid}$, a novel massive simulation suite spanning six different cosmological parameters, we generate almost one million tomographic mock surveys on the sphere. Due to the large data set size and survey area, we perform a spherical analysis while limiting our map resolution to $\texttt{HEALPix}$ $n_\mathrm{side}=512$. We marginalize over systematics such as photometric redshift errors, multiplicative calibration and additive shear bias. Furthermore, we use a map-level implementation of the non-linear intrinsic alignment model along with a novel treatment of baryonic feedback to incorporate additional astrophysical nuisance parameters. We also perform a spherical power spectrum analysis for comparison. The constraints of the cosmological parameters are generated using a likelihood free inference method called Gaussian Process Approximate Bayesian Computation (GPABC). Finally, we check that our pipeline is robust against choices of the simulation parameters. We find constraints on the degeneracy parameter of $S_8 \equiv σ_8\sqrt{Ω_M/0.3} = 0.78^{+0.06}_{-0.06}$ for our power spectrum analysis and $S_8 = 0.79^{+0.05}_{-0.05}$ for our GCNN analysis, improving the former by 16%. This is consistent with earlier analyses of the 2-point function, albeit slightly higher. Baryonic corrections generally broaden the constraints on the degeneracy parameter by about 10%. These results offer great prospects for full machine learning based analyses of on-going and future weak lensing surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Constraining baryonic feedback and cosmology with weak-lensing, X-ray, and kinematic Sunyaev-Zeldovich observations

Modern weak-lensing observations are becoming increasingly sensitive to baryonic feedback processes which are still poorly understood. So far, this challenge has been faced either by imposing scale-cuts in the data or by modelling baryonic effects with simple, one-parameter models. In this paper, we rely on a more general, seven-parameter prescription of baryonic feedback effects, which is primarily motivated by observations and has been shown to agree with a plethora of hydrodynamical simulations. By combining weak-lensing data from the Kilo-Degree Survey (KiDS-1000) with observations of gas around galaxy clusters, we are able to constrain baryonic parameters and learn more about feedback and cosmology. In particular, we use cluster gas fractions from X-ray data and gas profiles from kinematic Sunyaev-Zeldovich (kSZ) observations to provide evidence for baryonic feedback that is stronger than predicted by most hydrodynamical simulations. In terms of the matter power spectrum, we report a beyond-percent effect at wave-modes above $k\sim 0.1-0.45$ h/Mpc and a maximum suppression of $12-33$ percent at $k\sim7$ h/Mpc (68 percent confidence level). Regarding the combined parameter $Σ_8=σ_8(Ω_m/0.3)^{0.58}$, we find the known tension with the Planck satellite data to be reduced from 3.8 to 2.9 $σ$ once baryonic effects are fully included in the analysis pipeline. The tension is further decreased to 2.6 $σ$ when the weak-lensing data is combined with X-ray and kSZ observations. We conclude that, while baryonic feedback effects become more important in modern weak-lensing surveys, they are unlikely to act as the main culprit for the observed $Σ_8$-tension.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Imprints of fermionic and bosonic mixed dark matter on the 21-cm signal at cosmic dawn

The 21-cm signal from the epoch of cosmic dawn prior to reionization consists of a promising observable to gain new insights into the dark matter (DM) sector. In this paper, we investigate its potential to constrain mixed (cold + non-cold) dark matter scenarios that are characterised by the non-cold DM fraction ($f_{\rm nCDM}$) and particle mass ($m_{\rm nCDM}$). As non-cold DM species, we investigate both a fermionic (sterile neutrino) and a bosonic (ultra-light axion) particle. We show how these scenarios affect the global signal and the power spectrum using a halo-model implementation of the 21-cm signal at cosmic dawn. Next to this study, we perform an inference-based forecast study based on realistic mock power spectra from the Square Kilometre Array (SKA) telescope. Assuming inefficient, yet non-zero star-formation in minihaloes (i.e. haloes with mass below $10^8$ M$_{\odot}$), we obtain stringent constraints on both $m_{\rm nCDM}$ and $f_{\rm nCDM}$ that go well beyond current limits. Regarding the special case of $f_{\rm nCDM}\sim 1$, for example, we find a constraint of $m_{\rm nCDM}>15$ keV (thermal mass) for fermionic DM and $m_{\rm nCDM}>2\times10^{-20}$ eV for bosonic DM. For the opposite case of dominating cold DM, we find that at most one percent of the total DM abundance can be made of a hot fermionic or bosonic relic. All constraints are provided at the 95 percent confidence level.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Emulation of baryonic effects on the matter power spectrum and constraints from galaxy cluster data

Baryonic feedback effects consist of a major systematic for upcoming weak-lensing and galaxy-clustering surveys. In this paper, we present an emulator for the baryonic suppression of the matter power spectrum. The emulator is based on the baryonification model, containing seven free parameters that are connected to the gas profiles and stellar abundances in haloes. We show that with the baryonic emulator, we can not only recover the power spectra of hydro-dynamical simulations at sub-percent precision but also establish a connection between the baryonic suppression of the power spectrum and the gas and stellar fractions in haloes. This connection allows us to predict the expected deviation from a dark-matter-only power spectrum using measured X-ray gas fractions of galaxy groups and clusters. With these measurements, we constrain the suppression to exceed the percent-level at k = 0.1-0.4 h/Mpc and to reach a maximum of 20-28 percent at around k = 7 h/Mpc (68 percent confidence level). As a further step, we also perform a detailed parameter study and we present a minimum set of four baryonic parameters that are required to recover the scale and redshift dependence observed in hydro-dynamical simulations. The baryonic emulator can be found at https://github.com/sambit-giri/BCemu.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Parameter inference with non-linear galaxy clustering: accounting for theoretical uncertainties

We implement EuclidEmulator (version 1), an emulator for the non-linear correction of the matter power spectrum, into the MCMC forecasting code MontePython. We compare the performance of Halofit, HMCode, and EuclidEmulator1, both at the level of power spectrum prediction and at the level of posterior probability distributions of the cosmological parameters, for different cosmological models and different galaxy power spectrum wave number cut-offs. We confirm that the choice of the power spectrum predictor has a non-negligible effect on the computed sensitivities when doing cosmological parameter forecasting, even for a conservative wave number cut-off of $0.2\,h\,{\rm Mpc}^{-1}$. We find that EuclidEmulator1 is on average up to $17\%$ more sensitive to the cosmological parameters than the other two codes, with the most significant improvements being for the Hubble parameter of up to $42\%$ and the equation of state of dark energy of up to $26\%$, depending on the case. In addition, we point out that the choice of the power spectrum predictor contributes to the risk of computing a significantly biased mean cosmology when doing parameter estimations. For the four tested scenarios we find biases, averaged over the cosmological parameters, of between 0.5 and 2$σ$ (from below $1σ$ up to $6σ$ for individual parameters). This paper provides a proof of concept that this risk can be mitigated by taking a well-tailored theoretical uncertainty into account as this allows to reduce the bias by a factor of 2 to 5, depending on the case under consideration, while keeping posterior credibility contours small: the standard deviations are amplified by a factor of $\leq1.4$ in all cases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Baryonic effects for weak lensing. Part I. Power spectrum and covariance matrix

Baryonic feedback effects lead to a suppression of the weak lensing angular power spectrum on small scales. The poorly constrained shape and amplitude of this suppression is an important source of uncertainties for upcoming cosmological weak lensing surveys such as Euclid or LSST. In this first paper in a series of two, we use simulations to build a Euclid-like tomographic mock data-set for the cosmic shear power spectrum and the corresponding covariance matrix, which are both corrected for baryonic effects following the baryonification method of Schneider et al. (2019). In addition, we develop an emulator to obtain fast predictions of the baryonic power suppression, allowing us to perform a likelihood inference analysis for a standard $Λ$CDM cosmology with both cosmological and astrophysical parameters. Our main findings are the following: (i) ignoring baryonic effects leads to a greater than 5$σ$ bias on the cosmological parameters $Ω_m$ and $σ_8$; (ii) restricting the analysis to the largest scales, that are mostly unaffected by baryons, makes the bias disappear, but results in a blow-up of the $Ω_m$-$σ_8$ contour area by more than a factor of 10; (iii) ignoring baryonic effects on the covariance matrix does not significantly affect cosmological parameter estimates; (iv) while the baryonic suppression is mildly cosmology dependent, this effect does not noticeably modify the posterior contours. Overall, we conclude that including baryonic uncertainties in terms of nuisance parameters results in unbiased and surprisingly tight constraints on cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Baryonic effects for weak lensing. Part II. Combination with X-ray data and extended cosmologies

An accurate modelling of baryonic feedback effects is required to exploit the full potential of future weak-lensing surveys such as Euclid or LSST. In this second paper in a series of two, we combine Euclid-like mock data of the cosmic shear power spectrum with an eROSITA X-ray mock of the cluster gas fraction to run a combined likelihood analysis including both cosmological and baryonic parameters. Following the first paper of this series, the baryonic effects (based on the baryonic correction model of Schneider et al. 2019) are included in both the tomographic power spectrum and the covariance matrix. However, this time we assume the more realistic case of a $Λ$CDM cosmology with massive neutrinos, and we consider several extensions of the currently favoured cosmological model. For the standard $Λ$CDM case, we show that including X-ray data reduces the uncertainties on the sum of the neutrino mass by $\sim30$ percent, while there is only a mild improvement on other parameters such as $Ω_m$ and $σ_8$. As extensions of $Λ$CDM, we consider the cases of a dynamical dark energy model (wCDM), a $f(R)$ gravity model (fRCDM), and a mixed dark matter model ($Λ$MDM) with both a cold and a warm/hot dark matter component. We find that combining weak lensing with X-ray data only leads to a mild improvement of the constraints on the additional parameters of wCDM, while the improvement is more substantial for both fRCDM and $Λ$MDM. Ignoring baryonic effects in the analysis pipeline leads to significant false-detections of either phantom dark energy or a light subdominant dark matter component. Overall we conclude that for all cosmologies considered, a general parametrisation of baryonic effects is both necessary and sufficient to obtain tight constraints on cosmological parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A new method to quantify the effects of baryons on the matter power spectrum

Future large-scale galaxy surveys have the potential to become leading probes for cosmology provided the influence of baryons on the total mass distribution is understood well enough. As hydrodynamical simulations strongly depend on details in the feedback implementations, no unique and robust predictions for baryonic effects currently exist. In this paper we propose a baryonic correction model that modifies the density field of dark-matter-only $N$-body simulations to mimic the effects of baryons from any underlying adopted feedback recipe. The model assumes haloes to consist of 4 components: 1- hot gas in hydrostatical equilibrium, 2- ejected gas from feedback processes, 3- central galaxy stars, and 4- adiabatically relaxed dark matter, which all modify the initial dark-matter-only density profiles. These altered profiles allow to define a displacement field for particles in $N$-body simulations and to modify the total density field accordingly. The main advantage of the baryonic correction model is to connect the total matter density field to the observable distribution of gas and stars in haloes, making it possible to parametrise baryonic effects on the matter power spectrum. We show that the most crucial quantities are the mass fraction of ejected gas and its corresponding ejection radius. The former controls how strongly baryons suppress the power spectrum, while the latter provides a measure of the scale where baryonic effects become important. A comparison with X-ray and Sunyaev-Zel'dovich cluster observations suggests that baryons suppress wave modes above $k\sim0.5$ h/Mpc with a maximum suppression of 10-25 percent around $k\sim 2$ h/Mpc. More detailed observations of the gas in the outskirts of groups and clusters are required to decrease the large uncertainties of these numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Astrophysical constraints on resonantly produced sterile neutrino dark matter

Resonantly produced sterile neutrinos are considered an attractive dark matter (DM) candidate only requiring a minimal, well motivated extension to the standard model of particle physics. With a particle mass restricted to the keV range, sterile neutrinos are furthermore a prime candidate for warm DM, characterised by suppressed matter perturbations at the smallest observable scales. In this paper we take a critical look at the validity of the resonant scenario in the context of constraints from structure formation. We compare predicted and observed number of Milky-Way satellites and we introduce a new method to generalise existing Lyman-$α$ limits based on thermal relic warm DM to the case of resonant sterile neutrino DM. The tightest limits come from the Lyman-$α$ analysis, excluding the entire parameter space (at 2-$σ$ confidence level) still allowed by X-ray observations. Constraints from Milky-Way satellite counts are less stringent, leaving room for resonant sterile neutrino DM most notably around the suggested line signal at 7.1 keV.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Matter power spectrum and the challenge of percent accuracy

Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day $N$-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used $N$-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at $k\leq1$ $h\,\rm Mpc^{-1}$ and to within three percent at $k\leq10$ $h\,\rm Mpc^{-1}$. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for $k\leq 2$ $h\,\rm Mpc^{-1}$. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than $L=0.5$ $h^{-1}\rm Gpc$ to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of $M_{\rm p}=10^{9}$ $h^{-1}\rm M_{\odot}$ is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sussing Merger Trees: Stability and Convergence

Merger trees are routinely used to follow the growth and merging history of dark matter haloes and subhaloes in simulations of cosmic structure formation. Srisawat et al. (2013) compared a wide range of merger-tree-building codes. Here we test the influence of output strategies and mass resolution on tree-building. We find that, somewhat surprisingly, building the tree from more snapshots does not generally produce more complete trees; instead, it tends to short- en them. Significant improvements are seen for patching schemes which attempt to bridge over occasional dropouts in the underlying halo catalogues or schemes which combine the halo-finding and tree-building steps seamlessly. The adopted output strategy does not affec- t the average number of branches (bushiness) of the resultant merger trees. However, mass resolution has an influence on both main branch length and the bushiness. As the resolution increases, a halo with the same mass can be traced back further in time and will encounter more small progenitors during its evolutionary history. Given these results, we recommend that, for simulations intended as precursors for galaxy formation models where of order 100 or more snapshots are analysed, the tree-building routine should be integrated with the halo finder, or at the very least be able to patch over multiple adjacent snapshots.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dodelson-Widrow Production of Sterile Neutrino Dark Matter with Non-Trivial Initial Abundance

The simplest way to create sterile neutrinos in the early Universe is by their admixture to active neutrinos. However, this mechanism, connected to the Dark Matter (DM) problem by Dodelson and Widrow (DW), cannot simultaneously meet the relic abundance constraint as well as bounds from structure formation and X-rays. Nonetheless, unless a symmetry forces active-sterile mixing to vanish exactly, the DW mechanism will unavoidably affect the sterile neutrino DM population created by any other production mechanism. We present a semi-analytic approach to the DW mechanism acting on an arbitrary initial abundance of sterile neutrinos, allowing to combine DW with any other preceding production mechanism in a physical and precise way. While previous analyses usually assumed that the spectra produced by DW and another mechanism can simply be added, we use our semi-analytic results to discuss the validity of this assumption and to quantify its accurateness, thereby also scrutinising the DW spectrum and the derived mass bounds. We then map our results to the case of sterile neutrino DM from the decay of a real SM singlet coupled to the Higgs. Finally, we will investigate aspects of structure formation beyond the usual simple free-streaming estimates in order to judge on the effects of the DW modification on the sterile neutrino DM spectra generated by scalar decay.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Production of Sterile Neutrino Dark Matter and the 3.5 keV line

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Structure formation with suppressed small-scale perturbations

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 nonlinear 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The same with less: The cosmic web of warm versus cold dark matter dwarf galaxies

We explore fundamental properties of the distribution of low mass dark matter halos within the cosmic web using warm dark matter (WDM) and cold dark matter (CDM) cosmological simulations. Using self abundance-matched mock galaxy catalogs, we show that the distribution of dwarf galaxies in a WDM universe, wherein low mass halo formation is heavily suppressed, is nearly indistinguishable to that of a CDM universe whose low mass halos are not seen because galaxy formation is suppressed below some threshold halo mass. However, if the scatter between dwarf galaxy luminosity and halo properties is large enough, low mass CDM halos would sometimes host relatively bright galaxies thereby populating CDM voids with the occasional isolated galaxy and reducing the numbers of completely empty voids. Otherwise, without high mass to light scatter, all mock galaxy clustering statistics that we consider--the auto-correlation function, the numbers and radial profiles of satellites, the numbers of isolated galaxies, and the PDF of small voids--are nearly identical in CDM and WDM. WDM voids are neither larger nor emptier than CDM voids, when constructed from abundance-matched halo catalogs. It is thus a challenge to determine whether the CDM problem of the over-abundance of small halos with respect to the number density of observed dwarf galaxies has a cosmological solution or an astrophysical solution. However, some clues about the dark matter particle and the scatter between the properties of dwarf galaxies and their dark matter halo hosts might be found in the cosmic web of galaxies in future surveys of the local volume.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sussing Merger Trees : The Impact of Halo Merger Trees on Galaxy Properties in a Semi-Analytic Model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sussing Merger Trees: the influence of the halo finder

Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of \subhalos) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The life and death of cosmic voids

We investigate the formation, growth, merger history, movement, and destruction of cosmic voids detected via the watershed transform code VIDE in a cosmological N-body dark matter ΛCDM simulation. By adapting a method used to construct halo merger trees, we are able to trace individual voids back to their initial appearance and record the merging and evolution of their progenitors at high redshift. For the scales of void sizes captured in our simulation, we find that the void formation rate peaks at scale factor 0.3, which coincides with a growth in the void hierarchy and the emergence of dark energy. Voids of all sizes appear at all scale factors, though the median initial void size decreases with time. When voids become detectable they have nearly their present-day volumes. Almost all voids have relatively stable growth rates and suffer only infrequent minor mergers. Dissolution of a void via merging is very rare. Instead, most voids maintain their distinct identity as annexed subvoids of a larger parent. The smallest voids are collapsing at the present epoch, but void destruction ceases after scale factor 0.3. In addition, voids centers tend to move very little, less than 0.01 of their effective radii per ln a, over their lifetimes. Overall, most voids exhibit little radical dynamical evolution; their quiet lives make them pristine probes of cosmological initial conditions and the imprint of dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Warm dark matter does not do better than cold dark matter in solving small-scale inconsistencies

Over the last decade, warm dark matter (WDM) has been repeatedly proposed as an alternative scenario to the standard cold dark matter (CDM) one, potentially resolving several disagreements between the CDM model and observations on small scales. Here, we reconsider the most important CDM small-scale discrepancies in the light of recent observational constraints on WDM. As a result, we find that a conventional thermal (or thermal-like) WDM cosmology with a particle mass in agreement with Lyman-$α$ is nearly indistinguishable from CDM on the relevant scales and therefore fails to alleviate any of the small-scale problems. The reason for this failure is that the power spectrum of conventional WDM falls off too rapidly. To maintain WDM as a significantly different alternative to CDM, more evolved production mechanisms leading to multiple dark matter components or a gradually decreasing small-scale power spectrum have to be considered.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Halo Mass Function and the Free Streaming Scale

The nature of structure formation around the particle free streaming scale is still far from understood. Many attempts to simulate hot, warm, and cold dark matter cosmologies with a free streaming cutoff have been performed with cosmological particle-based simulations, but they all suffer from spurious structure formation at scales below their respective free streaming scales -- i.e. where the physics of halo formation is most affected by free streaming. We perform a series of high resolution numerical simulations of different WDM models, and develop an approximate method to subtract artificial structures in the measured halo mass function. The corrected measurements are then used to construct and calibrate an extended Press-Schechter (EPS) model with sharp-$k$ window function and adequate mass assignment. The EPS model gives accurate predictions for the low redshift halo mass function of CDM and WDM models, but it significantly under-predicts the halo abundance at high redshifts. By taking into account the ellipticity of the initial patches and connecting the characteristic filter scale to the smallest ellipsoidal axis, we are able to eliminate this inconsistency and obtain an accurate mass function over all redshifts and all dark matter particle masses covered by the simulations. As an additional application we use our model to predict the microhalo abundance of the standard neutralino-CDM scenario and we give the first quantitative prediction of the mass function over the full range of scales of CDM structure formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sussing Merger Trees: The Merger Trees Comparison Project

Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the SUSSING MERGER TREES Workshop in which ten different tree-building algorithms were applied to the same set of halo catalogues and their results compared. Although many of these codes were similar in nature, all algorithms produced distinct results. Our main conclusions are that a useful merger-tree code should possess the following features: (i) the use of particle IDs to match haloes between snapshots; (ii) the ability to skip at least one, and preferably more, snapshots in order to recover subhaloes that are temporarily lost during merging; (iii) the ability to cope with (and ideally smooth out) large, temporary flucuations in halo mass. Finally, to enable different groups to communicate effectively, we defined a common terminology that we used when discussing merger trees and we encourage others to adopt the same language. We also specified a minimal output format to record the results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Galactic Halo in Mixed Dark Matter Cosmologies

A possible solution to the small scale problems of the cold dark matter (CDM) scenario is that the dark matter consists of two components, a cold and a warm one. We perform a set of high resolution simulations of the Milky Way halo varying the mass of the WDM particle ($m_{\rm WDM}$) and the cosmic dark matter mass fraction in the WDM component ($\bar{f}_{\rm W}$). The scaling ansatz introduced in combined analysis of LHC and astroparticle searches postulates that the relative contribution of each dark matter component is the same locally as on average in the Universe (e.g. $f_{\rm W,\odot} = \bar{f}_{\rm W}$). Here we find however, that the normalised local WDM fraction ($f_{\rm W,\odot}$ / $\bar{f}_{\rm W}$) depends strongly on $m_{\rm WDM}$ for $m_{\rm WDM} <$ 1 keV. Using the scaling ansatz can therefore introduce significant errors into the interpretation of dark matter searches. To correct this issue a simple formula that fits the local dark matter densities of each component is provided.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Toward an accurate mass function for precision cosmology

Cosmological surveys aim to use the evolution of the abundance of galaxy clusters to accurately constrain the cosmological model. In the context of LCDM, we show that it is possible to achieve the required percent level accuracy in the halo mass function with gravity-only cosmological simulations, and we provide simulation start and run parameter guidelines for doing so. Some previous works have had sufficient statistical precision, but lacked robust verification of absolute accuracy. Convergence tests of the mass function with, for example, simulation start redshift can exhibit false convergence of the mass function due to counteracting errors, potentially misleading one to infer overly optimistic estimations of simulation accuracy. Percent level accuracy is possible if initial condition particle mapping uses second order Lagrangian Perturbation Theory, and if the start epoch is between 10 and 50 expansion factors before the epoch of halo formation of interest. The mass function for halos with fewer than ~1000 particles is highly sensitive to simulation parameters and start redshift, implying a practical minimum mass resolution limit due to mass discreteness. The narrow range in converged start redshift suggests that it is not presently possible for a single simulation to capture accurately the cluster mass function while also starting early enough to model accurately the numbers of reionisation era galaxies, whose baryon feedback processes may affect later cluster properties. Ultimately, to fully exploit current and future cosmological surveys will require accurate modeling of baryon physics and observable properties, a formidable challenge for which accurate gravity-only simulations are just an initial step.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cores in warm dark matter haloes: a Catch 22 problem

The free streaming of warm dark matter particles dampens the fluctuation spectrum, flattens the mass function of haloes and imprints a fine grained phase density limit for dark matter structures. The phase space density limit is expected to imprint a constant density core at the halo center on the contrary to what happens for cold dark matter. We explore these effects using high resolution simulations of structure formation in different warm dark matter scenarios. We find that the size of the core we obtain in simulated haloes is in good agreement with theoretical expectations based on Liouville's theorem. However, our simulations show that in order to create a significant core, (r_c~1 kpc), in a dwarf galaxy (M~1e10 Msun), a thermal candidate with a mass as low as 0.1 keV is required. This would fully prevent the formation of the dwarf galaxy in the first place. For candidates satisfying large scale structure constrains (m_wdm larger than 1-2 keV) the expected size of the core is of the order of 10 (20) pc for a dark matter halo with a mass of 1e10 (1e8) Msun. We conclude that "standard" warm dark matter is not viable solution for explaining the presence of cored density profiles in low mass galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear Evolution of Cosmological Structures in Warm Dark Matter Models

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 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 ~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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Microhalos and Dark Matter Detection

Cosmological structure formation predicts that our galactic halo contains an enormous hierarchy of substructures and streams, the remnants of the merging hierarchy that began with tiny Earth mass microhalos. If these structures persist until the present time, they could influence dramatically the detection signatures of weakly interacting elementary particle dark matter (WIMP). Using numerical simulations that follow the tidal disruption within the Galactic potential and heating from stellar encounters, we find that neither microhalos nor streams have significant impact on direct detection, implying that dark matter constraints derived using simple smooth halo models are relatively robust. We also find that many dense central cusps survive, yielding a small enhancement in the signal for indirect detection experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the Stability of Tidal Streams

We explore the stability of tidal streams to perturbations, motivated by recent claims that the clumpy structure of the stellar streams surrounding the globular cluster Palomar 5 are the result of gravitational instability. We calculate the Jeans length of tidal streams by treating them as a thin expanding cylinder of collisionless matter. We also find a general relation between the density and the velocity dispersion inside a stream, which is used to determine the longitudinal Jeans criterion. Our analytic results are checked by following the time evolution of the phase space density within streams using numerical simulations. We conclude that tidal streams within our galactic halo are stable on all length scales and over all timescales.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Impact of Dark Matter Microhalos on Signatures for Direct and Indirect Detection

Detecting dark matter as it streams through detectors on Earth relies on knowledge of its phase space density on a scale comparable to the size of our solar system. Numerical simulations predict that our Galactic halo contains an enormous hierarchy of substructures, streams and caustics, the remnants of the merging hierarchy that began with tiny Earth mass microhalos. If these bound or coherent structures persist until the present time, they could dramatically alter signatures for the detection of weakly interacting elementary particle dark matter (WIMP). Using numerical simulations that follow the coarse grained tidal disruption within the Galactic potential and fine grained heating from stellar encounters, we find that microhalos, streams and caustics have a negligible likelihood of impacting direct detection signatures implying that dark matter constraints derived using simple smooth halo models are relatively robust. We also find that many dense central cusps survive, yielding a small enhancement in the signal for indirect detection experiments.

0 citations0 reviewsNo code linkedOpen access