Search

The Voids-within-Voids-within-Voids project used dark-matter-only (DMO) simulations to study the abundance and structure of dark matter (DM) haloes over the full mass range populated in the standard Lambda cold dark matter cosmology. Here, we explore how baryonic effects modify these results for |$z=0$| halo masses in the range |$10^4$| – |$10^7~\mathrm{M_\odot }$|⁠ , below the threshold for galaxy formation. Our main study focuses on three simulations from identical initial conditions at |$z=127$|⁠ , one following DMO, one including non-radiative gas, and one additionally including the baryonic physics relevant in this halo mass range (cooling and photoheating). In the non-radiative simulation, above |$10^{5.5}~\mathrm{M_\odot }$|⁠ , halo abundance and internal structure are very similar to the DMO simulation, and the baryon to DM ratio is everywhere close to the cosmic value. At lower mass, this ratio drops and haloes are less concentrated and less massive in the non-radiative case. Test simulations at higher resolution show this to be mainly a resolution effect; the expected drop in baryon content due to residual pressure effects only becomes substantial for |$z=0$| haloes below |${\sim}10^{2.7}~\mathrm{M_\odot }$|⁠. However, gas is heated by reionization at |$z=6$| in our 'full physics' run, and this results in almost complete expulsion of gas from all haloes in our simulated mass range. This suppresses the halo mass function by |${\sim}30{{\ \rm per\ cent}}$|⁠ , lowers halo concentration, and consequently weakens the DM annihilation signal by |${\sim}40{-}60{{\ \rm per\ cent}}$|⁠. [ABSTRACT FROM AUTHOR]

Dwarf satellites in galaxy groups are distributed in an anisotropic and asymmetric manner, which is called the 'lopsided satellite distribution'. This lopsided signal has been observed not only in galaxy pairs but also in isolated systems. However, the physical origin of the lopsided signal in isolated systems is still unknown. In this work, we investigate this in the state-of-the-art hydrodynamical simulation of the MillenniumTNG Project by tracing each system back to high redshift. We find that the lopsided signal is dominated by satellites located in the outer regions of the halo and is also dominated by recently accreted satellites. The lopsided signal originates from the anisotropic accretion of galaxies from the surrounding large-scale structure and that, after accretion, the non-linear evolution of satellites inside the dark matter halo weakens the lopsidedness. The signal decreases as cosmic time passes because of a competition between anisotropic accretion and internal evolution within dark matter haloes. Our findings provide a useful perspective for the study of galaxy evolution, especially for the origin of the spatial satellite galaxy distributions. [ABSTRACT FROM AUTHOR]

We use the voids-within-voids-within-voids (VVV) simulations, a suite of successive nested N -body simulations with extremely high resolution (denoted, from low to high resolution, by L0 to L7), to test the Press–Schechter (PS), Sheth–Tormen (ST), and extended Press–Schechter (EPS) formulae for the halo abundance over the entire mass range, from minihaloes of 10−6 M⊙, to cluster haloes of 1015 M⊙, at different redshifts, from z = 30 to the present. We find that at z = 0 and z = 2, ST best reproduces the results of L0, which has the mean cosmic density (overdensity δ = 0), at 1011−15 M⊙. The higher resolution levels (L1–L7) are biased underdense regions (δ < −0.6). The EPS formalism takes this into account since it gives the mass function of a region conditioned, in this case, on having a given underdensity. EPS provides good matches to these higher levels, with deviations ≲20 per cent, at 10−6−12.5 M⊙. At z ∼ 7−15, the ST predictions for L0 and the EPS for L1–L7 show somewhat larger deviations from the simulation results. However, at even higher redshifts, z ∼ 30, EPS fits the simulations well again. We confirm our results by picking more subvolumes from the L0 simulation, finding that our conclusions depend only weakly on the size and overdensity of the region. The good agreement of EPS with the higher level simulations implies that PS (or ST) gives an accurate description of the total halo mass function in representative regions of the universe. [ABSTRACT FROM AUTHOR]

We take advantage of the unprecedented dynamical range provided by the 'Cosmic-Zoom' project to study the mass accretion history (MAH) of present-day dark matter haloes over the entire mass range present in the Lambda cold dark matter paradigm when the dark matter is made of weakly interacting massive particles of mass 100 GeV. In particular, we complement previous studies by exploring the MAHs of haloes with mass from |$10^8\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$| down to Earth mass, |$10^{-6}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$|⁠. The formation redshift of low-mass haloes anticorrelates weakly with mass, peaking at z = 3 for haloes of mass |$10^{-4}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$|⁠. Even lower masses are affected by the free-streaming cut-off in the primordial spectrum of density fluctuations and form at lower redshift. We compare MAHs in our simulations with predictions from two analytical models based on the extended Press–Schechter theory (EPS), and three empirical models derived by fitting and extrapolating either results from cosmological N -body simulations or Monte Carlo realizations of halo growth. All models fit our simulations reasonably well over the mass range for which they were calibrated. While the empirical models match better for more massive haloes, |$M\gt 10^{10}\ h^{-1}\,\mathrm{{\rm M}_{\odot }}$|⁠ , the analytical models do better when extrapolated down to Earth mass. At the higher masses, we explore the correlation between local environment density and MAH, finding that biases are relatively weak, with typical MAHs for haloes in extremely low-density and in typical regions differing by less than 20 per cent at high redshift. If this result can be extrapolated to lower halo masses, we conclude that EPS theory predicts the hierarchical build up of dark matter haloes quite well over the entire halo mass range. [ABSTRACT FROM AUTHOR]

Upcoming large galaxy surveys will subject the standard cosmological model, Lambda Cold Dark Matter, to new precision tests. These can be tightened considerably if theoretical models of galaxy formation are available that can predict galaxy clustering and galaxy–galaxy lensing on the full range of measurable scales, throughout volumes as large as those of the surveys, and with sufficient flexibility that uncertain aspects of the underlying astrophysics can be marginalized over. This, in particular, requires mock galaxy catalogues in large cosmological volumes that can be directly compared to observation, and can be optimized empirically by Monte Carlo Markov Chains or other similar schemes, thus eliminating or estimating parameters related to galaxy formation when constraining cosmology. Semi-analytic galaxy formation methods implemented on top of cosmological dark matter simulations offer a computationally efficient approach to construct physically based and flexibly parametrized galaxy formation models, and as such they are more potent than still faster, but purely empirical models. Here, we introduce an updated methodology for the semi-analytic L-Galaxies code, allowing it to be applied to simulations of the new MillenniumTNG project, producing galaxies directly on fully continuous past lightcones, potentially over the full sky, out to high redshift, and for all galaxies more massive than |$\sim 10^8\, {\rm M}_\odot$|⁠. We investigate the numerical convergence of the resulting predictions, and study the projected galaxy clustering signals of different samples. The new methodology can be viewed as an important step towards more faithful forward-modelling of observational data, helping to reduce systematic distortions in the comparison of theory to observations. [ABSTRACT FROM AUTHOR]

Cosmological inference with large galaxy surveys requires theoretical models that combine precise predictions for large-scale structure with robust and flexible galaxy formation modelling throughout a sufficiently large cosmic volume. Here, we introduce the millenniumTNG (MTNG) project which combines the hydrodynamical galaxy formation model of illustrisTNG with the large volume of the millennium simulation. Our largest hydrodynamic simulation, covering |$(500 \, h^{-1}{\rm Mpc})^3 \simeq (740\, {\rm Mpc})^3$|⁠ , is complemented by a suite of dark-matter-only simulations with up to 43203 dark matter particles (a mass resolution of |$1.32\times 10^8 \, h^{-1}{\rm M}_\odot$|⁠) using the fixed-and-paired technique to reduce large-scale cosmic variance. The hydro simulation adds 43203 gas cells, achieving a baryonic mass resolution of |$2\times 10^7 \, h^{-1}{\rm M}_\odot$|⁠. High time-resolution merger trees and direct light-cone outputs facilitate the construction of a new generation of semi-analytic galaxy formation models that can be calibrated against both the hydro simulation and observation, and then applied to even larger volumes – MTNG includes a flagship simulation with 1.1 trillion dark matter particles and massive neutrinos in a volume of |$(3000\, {\rm Mpc})^3$|⁠. In this introductory analysis we carry out convergence tests on basic measures of non-linear clustering such as the matter power spectrum, the halo mass function and halo clustering, and we compare simulation predictions to those from current cosmological emulators. We also use our simulations to study matter and halo statistics, such as halo bias and clustering at the baryonic acoustic oscillation scale. Finally we measure the impact of baryonic physics on the matter and halo distributions. [ABSTRACT FROM AUTHOR]

We investigate cosmological structure formation in fuzzy dark matter (FDM) with the attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum 'pressure' and alter structure formation. We find in our simulations that the SI can enhance small-scale structure formation, and soliton cores above a critical mass undergo a phase transition, transforming from dilute to dense solitons. [ABSTRACT FROM AUTHOR]

Observations of the high redshift universe provide a promising avenue for constraining the nature of the dark matter (DM). This will be even more true with the advent of the James Webb Space Telescope (JWST). We run cosmological simulations of galaxy formation as part of the Effective Theory of Structure Formation (ETHOS) project to compare high redshift galaxies in cold dark matter (CDM) and alternative DM models which have varying relativistic coupling and self-interaction strengths. The interacting DM scenarios produce a cutoff in the linear power spectrum on small-scales, followed by a series of 'dark acoustic oscillations'. We find that DM interactions suppress the abundance of galaxies below |$M_\star \sim 10^8\, {\rm M}_\odot$| for the models considered. The cutoff in the power spectrum delays structure formation relative to CDM. Objects in ETHOS that end up at the same final masses as their CDM counterparts are characterized by a more vigorous phase of early star formation. While galaxies with |$M_\star \lesssim 10^6\, {\rm M_\odot }$| make up more than 60 per cent of star formation in CDM at z ≈ 10, they contribute only about half the star formation density in ETHOS. These differences diminish with decreasing redshift. We find that the effects of DM self-interactions are negligible compared to effects of relativistic coupling (i.e. the effective initial conditions for galaxy formation) in all properties of the galaxy population we examine. Finally, we show that the clustering strength of galaxies at high redshifts depends sensitively on DM physics, although these differences are manifest on scales that may be too small to be measurable by JWST. [ABSTRACT FROM AUTHOR]

We employ the hydrodynamical simulation illustrisTNG to inform the galaxy–halo connection of the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples of the Dark Energy Spectroscopic Instrument (DESI) survey at redshift z ∼ 0.8. Specifically, we model the galaxy colours of illustrisTNG and apply sliding DESI colour–magnitude cuts, matching the DESI target densities. We study the halo occupation distribution (HOD) model of the selected samples by matching them to their corresponding dark matter haloes in the illustrisTNG dark matter run. We find the HOD of both the LRG and ELG samples to be consistent with their respective baseline models, but also we find important deviations from common assumptions about the satellite distribution, velocity bias, and galaxy secondary biases. We identify strong evidence for concentration-based and environment-based occupational variance in both samples, an effect known as 'galaxy assembly bias'. The central and satellite galaxies have distinct dependencies on secondary halo properties, showing that centrals and satellites have distinct evolutionary trajectories and should be modelled separately. These results serve to inform the necessary complexities in modelling galaxy–halo connection for DESI analyses and also prepare for building high-fidelity mock galaxies. Finally, we present a shuffling-based clustering analysis that reveals a 10–15  |${{\ \rm per\ cent}}$| excess in the LRG clustering of modest statistical significance due to secondary galaxy biases. We also find a similar excess signature for the ELGs, but with much lower statistical significance. When a larger hydrodynamical simulation volume becomes available, we expect our analysis pipeline to pinpoint the exact sources of such excess clustering signatures. [ABSTRACT FROM AUTHOR]

Tracking the formation and evolution of dark matter haloes is a critical aspect of any analysis of cosmological N -body simulations. In particular, the mass assembly of a halo and its progenitors, encapsulated in the form of its merger tree, serves as a fundamental input for constructing semi-analytic models of galaxy formation and, more generally, for building mock catalogues that emulate galaxy surveys. We present an algorithm for constructing halo merger trees from abacussummit , the largest suite of cosmological N -body simulations performed to date consisting of nearly 60 trillion particles, and which has been designed to meet the Cosmological Simulation Requirements of the Dark Energy Spectroscopic Instrument (DESI) survey. Our method tracks the cores of haloes to determine associations between objects across multiple time slices, yielding lists of halo progenitors and descendants for the several tens of billions of haloes identified across the entire suite. We present an application of these merger trees as a means to enhance the fidelity of abacussummit halo catalogues by flagging and 'merging' haloes deemed to exhibit non-monotonic past merger histories. We show that this cleaning technique identifies portions of the halo population that have been deblended due to choices made by the halo finder, but which could have feasibly been part of larger aggregate systems. We demonstrate that by cleaning halo catalogues in this post-processing step, we remove potentially unphysical features in the default halo catalogues, leaving behind a more robust halo population that can be used to create highly accurate mock galaxy realizations from abacussummit. [ABSTRACT FROM AUTHOR]

We use N -body cosmological simulations and empirical galaxy models to study the merger history of dwarf-mass galaxies (with |$M_{\rm halo}\sim 10^{10}\, \mathrm{M}_\odot$|⁠). Our input galaxy models describe the stellar mass–halo mass relation, and the galaxy occupation fraction. The number of major and minor mergers depends on the type of dark matter; in particular, minor mergers are greatly suppressed in warm dark matter models. In addition, the number of mergers that bring in stars is strongly dependent on the galaxy occupation model. For example, minor mergers are negligible for stellar halo growth in models with a high mass threshold for galaxy formation (i.e. |$10^{9.3}\, \mathrm{M}_\odot$| at z = 0). Moreover, this threshold for galaxy formation can also determine the relative difference (if any) between the stellar haloes of satellite and field dwarfs. Using isolated simulations of dwarf–dwarf mergers, we show that the relative frequency of major and minor mergers predict very different stellar haloes: Typically, 'intermediate' dark matter merger ratios (∼1:5) maximize the growth of distant stellar haloes. We discuss the observability of dwarf stellar haloes and find that the surface brightness of these features are incredibly faint. However, when several dwarfs are stacked together, models that form particularly rich stellar haloes could be detectable. Finally, we show that stellar streams in the Galactic halo overlapping in phase space with known dwarf satellites are likely remnants of their stripped stellar haloes. The mere existence of dwarf stellar haloes can already put constraints on some small-scale models, and thus observational probes should be a high priority. [ABSTRACT FROM AUTHOR]

We introduce the AbacusHOD model and present two applications of AbacusHOD and the AbacusSummit simulations to observations. AbacusHOD is a Halo Occupation Distribution (HOD) framework written in Python that is particle-based, multitracer, highly generalized, and highly efficient. It is designed specifically with multitracer/cosmology analyses for next-generation large-scale structure surveys in mind, and takes advantage of the volume and precision offered by the new state-of-the-art AbacusSummit cosmological simulations. The model is also highly customizable and should be broadly applicable to any upcoming surveys and a diverse range of cosmological analyses. In this paper, we demonstrate the capabilities of the AbacusHOD framework through two example applications. The first example demonstrates the high efficiency and the large HOD extension feature set through an analysis of full-shape redshift-space clustering of BOSS galaxies at intermediate to small scales (⁠|$\lt 30\, h^{-1}$|  Mpc), assessing the necessity of introducing secondary galaxy biases (assembly bias). We find strong evidence for using halo environment instead of concentration to trace secondary galaxy bias, a result which also leads to a moderate reduction in the 'lensing is low' tension. The second example demonstrates the multitracer capabilities of the AbacusHOD package through an analysis of the extended Baryon Oscillation Spectroscopic Survey cross-correlation measurements between three different galaxy tracers: luminous red galaxies, emission-line galaxies, and quasi-stellar objects. We expect the AbacusHOD framework, in combination with the AbacusSummit simulation suite, to play an important role in a simulation-based analysis of the upcoming Dark Energy Spectroscopic Instrument data sets. [ABSTRACT FROM AUTHOR]

The James Webb Space Telescope (JWST) will revolutionize our understanding of early galaxy formation, and could potentially set stringent constraints on the nature of dark matter. We use a semi-empirical model of galaxy formation to investigate the extent to which uncertainties in the implementation of baryonic physics may be degenerate with the predictions of two different models of dark matter – cold dark matter (CDM) and a 7 keV sterile neutrino, which behaves as warm dark matter (WDM). Our models are calibrated to the observed UV luminosity function at z  = 4 using two separate dust attenuation prescriptions, which manifest as high and low star formation efficiency in low-mass haloes. These efficiencies capture the net effect of processes that regulate star formation. We find that while at fixed star formation efficiency, ε, there are marked differences in the abundance of faint galaxies in the two dark matter models at high- z ; these differences are mimicked easily by varying ε in the same dark matter model. We find that a high ε WDM and a low ε CDM model – which provide equally good fits to the z  = 4 UV luminosity function – exhibit nearly identical evolution in the cosmic stellar mass and star formation rate densities. We show that differences in the star formation rate at fixed stellar mass are larger for variations in ε in a given dark matter model than they are between dark matter models; however, the scatter in star formation rates is larger between the two models than they are when varying ε. Our results suggest that JWST will likely be more informative in constraining baryonic processes operating in high- z galaxies than it will be in constraining the nature of dark matter. [ABSTRACT FROM AUTHOR]

We investigate the morphology of the stellar distribution (SD) in a sample of Milky Way–like galaxies in the TNG50 simulation. Using a local in shell iterative method as the main approach, we explicitly show evidence of twisting (in about 52% of halos) and stretching (in 48% of them) in real space. This is matched with the reorientation observed in the eigenvectors of the inertia tensor and gives us a clear picture of having a reoriented SD. We make a comparison between the shape profile of the dark matter (DM) halo and SD and quite remarkably see that their radial profiles are fairly close, especially at small galactocentric radii, where the stellar disk is located. This implies that the DM halo is somewhat aligned with stars in response to the baryonic potential. The level of alignment mostly decreases away from the center. We study the impact of substructures in the orbital circularity parameter. It is demonstrated that in some cases, faraway substructures are counterrotating compared with the central stars and may flip the sign of total angular momentum and thus the orbital circularity parameter. Truncating them above 150 kpc, however, retains the disky structure of the galaxy as per initial selection. Including the impact of substructures in the shape of stars, we explicitly show that their contribution is subdominant. Overlaying our theoretical results on the observational constraints from previous literature, we establish fair agreement. [ABSTRACT FROM AUTHOR]

We present a comprehensive analysis of the shape of dark matter (DM) halos in a sample of 25 Milky Way-like galaxies in TNG50 simulation. Using an enclosed volume iterative method, we infer an oblate-to-triaxial shape for the DM halo with median T ≃ 0.24. We group DM halos into three different categories. Simple halos (32% of the population) establish principal axes whose ordering in magnitude does not change with radius and whose orientations are almost fixed throughout the halo. Twisted halos (32%) experience levels of gradual rotations throughout their radial profiles. Finally, stretched halos (36%) demonstrate a stretching in the lengths of their principal axes where the ordering of different eigenvalues changes with radius. Subsequently, the halo experiences a "rotation" of ∼90° where the stretching occurs. Visualizing the 3D ellipsoid of each halo, for the first time, we report signs of a reorienting ellipsoid in twisted and stretched halos. We examine the impact of baryonic physics on DM halo shape through a comparison to dark matter only (DMO) simulations. This suggests a triaxial (prolate) halo. We analyze the impacts of substructure on DM halo shape in both hydrodynamical and DMO simulations and confirm that they are subdominant. We study the distribution of satellites in our sample. In simple and twisted halos, the angle between satellites' angular momentum and the galaxy's angular momentum grows with radius. However, stretched halos show a flat distribution of angles. Overlaying our theoretical outcome on the observational results presented in the literature establishes a fair agreement. [ABSTRACT FROM AUTHOR]

We use the IllustrisTNG (TNG) simulations of galaxy formation to measure the velocity dispersion profiles of dark matter and stars in Milky Way–mass, galaxy group, and cluster-scale dark matter halos. The mean profiles calculated from both tracers are similar in shape, exhibiting a large halo-to-halo scatter around the average profile. The so-called "splashback" radius demarcates the outer boundary of the halo, and manifests as a kink in the velocity dispersion profile, located on average between ∼1.0–1.5r200m, where r200m is the radius within which the density of the halo equals 200 times the background density of the universe. We find that this location may also be identified as the radius at which the (stacked) dispersion profile drops to 60% of its peak value (for line-of-sight motions in TNG halos). We further show that the scatter in the dispersion profiles may be attributed to the variations in the assembly history of the host halos. In particular, this segregates the profile into two regimes: one within ∼0.1r200m, where the scatter is set by the early assembly history of the halo; and the other beyond this radius, where the scatter is influenced more strongly by its late-time assembly. Finally, we show that a two-parameter model can be used to fit the measured velocity dispersion profiles and the fit parameters can be related directly to two fundamental halo properties: mass and concentration. We describe a simple model that allows us to express the stellar velocity dispersion profile in terms of these halo properties only. [ABSTRACT FROM AUTHOR]

Building accurate and flexible galaxy–halo connection models is crucial in modelling galaxy clustering on non-linear scales. Recent studies have found that halo concentration by itself cannot capture the full galaxy assembly bias effect and that the local environment of the halo can be an excellent indicator of galaxy assembly bias. In this paper, we propose an extended halo occupation distribution (HOD) model that includes both a concentration-based assembly bias term and an environment-based assembly bias term. We use this model to achieve a good fit (χ2/degrees of freedom = 1.35) on the 2D redshift-space two-point correlation function (2PCF) of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy sample. We find that the inclusion of both assembly bias terms is strongly favoured by the data and the standard five-parameter HOD model is strongly rejected. More interestingly, the redshift-space 2PCF drives the assembly bias parameters in a way that preferentially assigns galaxies to lower mass haloes. This results in galaxy–galaxy lensing predictions that are within 1σ agreement with the observation, alleviating the perceived tension between galaxy clustering and lensing. We also showcase a consistent 3σ–5σ preference for a positive environment-based assembly bias that persists over variations in the fit. We speculate that the environmental dependence might be driven by underlying processes such as mergers and feedback, but might also be indicative of a larger halo boundaries such as the splashback radius. Regardless, this work highlights the importance of building flexible galaxy–halo connection models and demonstrates the extra constraining power of the redshift-space 2PCF. [ABSTRACT FROM AUTHOR]

Bose–Einstein condensate dark matter (BECDM, also known as fuzzy dark matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established cold dark matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts (z ≳ 5) for a boson mass m  = 2.5 × 10−22 eV, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered – a warm dark matter-like ('WDM') model often used as a proxy for BECDM. Our simulations confirm that 'WDM' is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/'WDM' compared to the CDM case: in BECDM/'WDM' first stars form at z ∼ 13/13.5, while in CDM star formation starts at z ∼ 35. The signature of BECDM interference, not present in 'WDM', is seen in the evolved dark matter power spectrum: although the small-scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM. [ABSTRACT FROM AUTHOR]

We use the IllustrisTNG (TNG) simulations to explore the galaxy–halo connection as inferred from state-of-the-art cosmological, magnetohydrodynamical simulations. With the high-mass resolution and large volume achieved by combining the 100 Mpc (TNG100) and 300 Mpc (TNG300) volumes, we establish the mean occupancy of central and satellite galaxies and their dependence on the properties of the dark matter haloes hosting them. We derive best-fitting HOD parameters from TNG100 and TNG300 for target galaxy number densities of |$\bar{n}_g = 0.032\,$|   and |$\bar{n}_g = 0.016\, h^3$|  Mpc−3, respectively, corresponding to a minimum galaxy stellar mass of |$M_\star \sim 1.9\times 10^9\, $| and |$M_\star \sim 3.5\times 10^9\, {\rm M}_\odot$|⁠ , respectively, in hosts more massive than |$10^{11}\, {\rm M}_\odot$|⁠. Consistent with previous work, we find that haloes located in dense environments, with low concentrations, later formation times, and high angular momenta are richest in their satellite population. At low mass, highly concentrated haloes and those located in overdense regions are more likely to contain a central galaxy. The degree of environmental dependence is sensitive to the definition adopted for the physical boundary of the host halo. We examine the extent to which correlations between galaxy occupancy and halo properties are independent and demonstrate that HODs predicted by halo mass and present-day concentration capture the qualitative dependence on the remaining halo properties. At fixed halo mass, concentration is a strong predictor of the stellar mass of the central galaxy, which may play a defining role in the fate of the satellite population. The radial distribution of satellite galaxies, which exhibits a universal form across a wide range of host halo mass, is described accurately by the best-fitting NFW density profile of their host haloes. [ABSTRACT FROM AUTHOR]

Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamical simulations of galaxy formation – including the APOSTLE and C-EAGLE simulations – we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55 keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the dark matter decay line for Perseus-analogue galaxy clusters lies in the range of 1300–1700  |${\rm km\, s^{-1}}$|  , and exceptionally up to 3000  |${\rm km\, s^{-1}}$|  . Therefore, the non-detection of the 3.55 keV line in the centre of the Perseus cluster by the Hitomi collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with XMM–Newton and with the planned X-ray telescopes XRISM and ATHENA. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter. [ABSTRACT FROM AUTHOR]

We study cosmic voids in the normal-branch Dvali–Gabadadze–Porrati (nDGP) braneworld models that are representative of a class of modified gravity theories where deviations from General Relativity are usually hidden by the Vainshtein screening in high-density environments. This screening is less efficient away from these environments, which makes voids ideally suited for testing this class of models. We use N -body simulations of Λ-cold dark matter (ΛCDM) and nGDP universes, where dark matter haloes are populated with mock galaxies that mimic the clustering and number densities of the boss cmass galaxy sample. We measure the force, density, and weak lensing profiles around voids identified with six different algorithms. Compared to ΛCDM, voids in nDGP are more underdense due to the action of the fifth force that arises in these models, which leads to a faster evacuation of matter from voids. This leaves an imprint on the weak lensing tangential shear profile around nDGP voids, an effect that is particularly strong for 2D underdensities that are identified in the plane of the sky. We make predictions for the feasibility of distinguishing between nDGP and ΛCDM using void lensing in upcoming large-scale surveys such as lsst and euclid. We compare with the analysis of voids in chameleon gravity theories and find that the weak lensing signal for 3D voids is similar to nDGP, whereas for 2D voids the differences with ΛCDM are much stronger for the chameleon gravity case, a direct consequence of the different screening mechanisms operating in these theories. [ABSTRACT FROM AUTHOR]

We study galaxy formation in sterile neutrino dark matter models that differ significantly from both cold and from 'warm thermal relic' models. We use the EAGLE code to carry out hydrodynamic simulations of the evolution of pairs of galaxies chosen to resemble the Local Group, as part of the APOSTLE simulations project. We compare cold dark matter (CDM) with two sterile neutrino models with 7 keV mass: one, the warmest among all models of this mass (LA120) and the other, a relatively cold case (LA10). We show that the lower concentration of sterile neutrino subhaloes compared to their CDM counterparts makes the inferred inner dark matter content of galaxies like Fornax (or Magellanic Clouds) less of an outlier in the sterile neutrino cosmologies. In terms of the galaxy number counts, the LA10 simulations are indistinguishable from CDM when one takes into account halo-to-halo (or 'simulation-to-simulation') scatter. In order for the LA120 model to match the number of Local Group dwarf galaxies, a higher fraction of low-mass haloes is required to form galaxies than is predicted by the EAGLE simulations. As the census of the Local Group galaxies nears completion, this population may provide a strong discriminant between cold and warm dark matter models. [ABSTRACT FROM AUTHOR]

N-body dark matter simulations of structure formation in the Λ cold dark matter (ΛCDM) model predict a population of subhaloes within Galactic haloes that have higher central densities than inferred for the Milky Way satellites, a tension known as the 'too big to fail' problem. Proposed solutions include baryonic effects, a smaller mass for the Milky Way halo and warm dark matter (WDM). We test these possibilities using a semi-analytic model of galaxy formation to generate luminosity functions for Milky Way halo-analogue satellite populations, the results of which are then coupled to the Jiang & van den Bosch model of subhalo stripping to predict the subhalo Vmax functions for the 10 brightest satellites. We find that selecting the brightest satellites (as opposed to the most massive) and modelling the expulsion of gas by supernovae at early times increases the likelihood of generating the observed Milky Way satellite Vmax function. The preferred halo mass is 6 × 1011 M⊙, which has a 14 per cent probability to host a Vmax function like that of the Milky Way satellites. We conclude that the Milky Way satellite Vmax function is compatible with a CDM cosmology, as previously found by Sawala et al. using hydrodynamic simulations. Sterile neutrino-WDM models achieve a higher degree of agreement with the observations, with a maximum 50 per cent chance of generating the observed Milky Way satellite Vmax function. However, more work is required to check that the semi-analytic stripping model is calibrated correctly for each sterile neutrino cosmology. [ABSTRACT FROM AUTHOR]

We use the Copernicus Complexio (COCO) high-resolution N-body simulations to investigate differences in the properties of small-scale structures in the standard cold dark matter (CDM) model and in a model with a cutoff in the initial power spectrum of density fluctuations consistent with both a thermally produced warm dark matter (WDM) particle with a rest mass of 3.3 keV and a sterile neutrino with mass 7 keV and leptogenesis parameter L6 = 8.7. The latter corresponds to the 'coldest' model with this sterile neutrino mass compatible with the identification of the recently detected 3.5 keV X-ray line as resulting from particle decay. CDM and WDM predict very different number densities of subhaloes with mass ≲ 109 h-1Mʘ although they predict similar, nearly universal, normalized subhalo radial density distributions. Haloes and subhaloes in both models have cuspy Navarro-Frenk-White profiles, but WDM subhaloes below the cut-off scale in the power spectrum (corresponding to maximum circular velocities Vz=0max ≤ 50 kms-1) are less concentrated than their CDM counterparts. We make predictions for observable properties using the GALFORM semi-analytic model of Galaxy formation. Both models predict Milky Way satellite luminosity functions consistent with observations, although the WDM model predicts fewer very faint satellites. This model, however, predicts slightly more UV bright galaxies at redshift z > 7 than CDM, but both are consistent with observations. Gravitational lensing offers the best prospect of distinguishing between the models. [ABSTRACT FROM AUTHOR]

We investigate the process of reionization in a model in which the dark matter is a warm elementary particle such as a sterile neutrino. We focus on models that are consistent with the dark matter decay interpretation of the recently detected line at 3.5 keV in the X-ray spectra of galaxies and clusters. In warm dark matter models, the primordial spectrum of density perturbations has a cut-off on the scale of dwarf galaxies. Structure formation therefore begins later than in the standard cold dark matter (CDM) model and very few objects form below the cut-off mass scale. To calculate the number of ionizing photons, we use the Durham semi-analytic model of galaxy formation, GALFORM.We find that even the most extreme 7 keV sterile neutrino we consider is able to reionize the Universe early enough to be compatible with the bounds on the epoch of reionization from Planck. This, perhaps surprising, result arises from the rapid build-up of high redshift galaxies in the sterile neutrino models which is also reflected in a faster evolution of their far-UV luminosity function between 10 > z > 7 than in CDM. The dominant sources of ionizing photons are systematically more massive in the sterile neutrino models than in CDM. As a consistency check on the models, we calculate the present-day luminosity function of satellites of Milky Way-like galaxies. When the satellites recently discovered in the Dark Energy Survey are taken into account, strong constraints are placed on viable sterile neutrino models. [ABSTRACT FROM AUTHOR]

Theoretical models of galaxy formation based on the cold dark matter cosmogony typically require strong feedback from supernova (SN) explosions in order to reproduce the Milky Way satellite galaxy luminosity function and the faint end of the field galaxy luminosity function. However, too strong a SN feedback also leads to the universe reionizing too late, and the metallicities of Milky Way satellites being too low. The combination of these four observations therefore places tight constraints on SN feedback. We investigate these constraints using the semi-analytical galaxy formation model GALFORM. We find that these observations favour a SN feedback model in which the feedback strength evolves with redshift. We find that, for our best-fitting model, half of the ionizing photons are emitted by galaxies with rest-frame far-UV absolute magnitudes MAB(1500Å) < -17.5, which implies that already observed galaxy populations contribute about half of the photons responsible for reionization. The z = 0 descendants of these galaxies are mainly galaxies with stellar mass M* > 1010M⊚ and preferentially inhabit haloes with mass Mhalo > 1013M⊚. [ABSTRACT FROM AUTHOR]

The sterile neutrino is a viable dark matter candidate that can be produced in the early Universe via non-equilibrium processes, and would therefore possess a highly non-thermal spectrum of primordial velocities. In this paper we analyse the process of structure formation with this class of dark matter particles. To this end we construct primordial dark matter power spectra as a function of the lepton asymmetry, L6, that is present in the primordial plasma and leads to resonant sterile neutrino production. We compare these power spectra with those of thermally produced dark matter particles and show that resonantly produced sterile neutrinos are much colder than their thermal relic counterparts. We also demonstrate that the shape of these power spectra is not determined by the free-streaming scale alone. We then use the power spectra as an input for semi-analytic models of galaxy formation in order to predict the number of luminous satellite galaxies in a Milky Way-like halo. By assuming that the mass of the Milky Way halo must be no more than 2 × 1012 M⊙ (the adopted upper bound based on current astronomical observations) we are able to constrain the value of L6 for Ms ≤ 8 keV. We also show that the range of L6 that is in best agreement with the 3.5 keV line (if produced by decays of 7 keV sterile neutrino) requires that the Milky Way halo has a mass no smaller than 1.5 × 1012 M⊙. Finally, we compare the power spectra obtained by direct integration of the Boltzmann equations for a non-resonantly produced sterile neutrino with the fitting formula of Viel et al. and find that the latter significantly underestimates the power amplitude on scales relevant to satellite galaxies. [ABSTRACT FROM AUTHOR]

We use a suite of cosmological simulations to study the mass-concentration-redshift relation, c(M, z), of dark matter haloes. Our simulations include standard Λ-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM haloes are selfsimilar and well approximated by the Einasto profile. The c(M, z) relation of CDM haloes is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM haloes are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM haloes: their MAHs are not scale-free because of the characteristic scale imposed by the power spectrum suppression. Further, the WDM c(M, z) relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the 'collapsed mass history'; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass-concentration-redshift relation of both CDM and WDM haloes over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations. [ABSTRACT FROM AUTHOR]

The cold dark matter (CDM) cosmological model unambiguously predicts that a large number of haloes should survive as subhaloes when they are accreted into a larger halo. The CDM model would be ruled out if such substructures were shown not to exist. By contrast, if the dark matter consists of Warm Dark Matter (WDM) particles, then below a threshold mass that depends on the particle mass far fewer substructures would be present. Finding subhaloes below a certain mass would then rule out warm particle masses below some value. Strong gravitational lensing provides a clean method to measure the subhalo mass function through distortions in the structure of Einstein rings and giant arcs. Using mock lensing observations constructed from high-resolution N-body simulations, we show that measurements of approximately 100 strong lens systems with a detection limit of Mlow = 107 h-1M⊙ would clearly distinguish CDM from WDM in the case where this consists of 7 keV sterile neutrinos such as those that might be responsible for the 3.5 keV X-ray emission line recently detected in galaxies and clusters. [ABSTRACT FROM AUTHOR]

We introduce Copernicus Complexio (coco), a high-resolution cosmological N-body simulation of structure formation in the ΛCDM model. coco follows an approximately spherical region of radius ~17.4 h−1 Mpc embedded in a much larger periodic cube that is followed at lower resolution. The high-resolution volume has a particle mass of 1.135 × 105 h−1 M⊙ (60 times higher than the Millennium-II simulation). coco gives the dark matter halo mass function over eight orders of magnitude in halo mass; it forms ~60 haloes of galactic size, each resolved with about 10 million particles. We confirm the power-law character of the subhalo mass function, .... (>μ)∝μ−s down to a reduced subhalo mass Msub/M200 ≡ μ = 10−6, with a best-fitting power-law index, s = 0.94, for hosts of mass 〈M200〉 = 1012 h−1 M⊙. The concentration-mass relation of coco haloes deviates from a single power law for masses M200 < afew × 108 h−1 M⊙, where it flattens, in agreement with results by Sanchez-Conde et al. The host mass invariance of the reduced maximum circular velocity function of subhaloes, ν ≡ Vmax/V200, hinted at in previous simulations, is clearly demonstrated over five orders of magnitude in host mass. Similarly, we find that the average, normalized radial distribution of subhaloes is approximately universal (i.e. independent of subhalo mass), as previously suggested by the Aquarius simulations of individual haloes. Finally, we find that at fixed physical subhalo size, subhaloes in lower mass hosts typically have lower central densities than those in higher mass hosts. [ABSTRACT FROM AUTHOR]

The recent detection of a 3.5 keV X-ray line from the centres of galaxies and clusters by Bulbul et al. and Boyarsky et al. has been interpreted as emission from the decay of 7 keV sterile neutrinos which could make up the (warm) dark matter (WDM). As part of the Copernicus Complexio (COCO) programme, we investigate the properties of dark matter haloes formed in a high-resolution cosmological N-body simulation from initial conditions similar to those expected in a universe in which the dark matter consists of 7 keV sterile neutrinos. This simulation and its cold dark matter (CDM) counterpart have ~13.4 bn particles, each of mass ~105 h-1M☉, providing detailed information about halo structure and evolution down to dwarf galaxy mass scales. Non-linear structure formation on small scales (M200 ≲2 × 109 h-1M☉) begins slightly later in COCO-WARM than in COCO-COLD. The halo mass function at the present day in the WDM model begins to drop below its CDM counterpart at a mass ~2 × 109 h-1M☉ and declines very rapidly towards lower masses so that there are five times fewer haloes of mass M200 = 108 h-1M☉ in COCO-WARM than in COCO-COLD. Halo concentrations on dwarf galaxy scales are correspondingly smaller in COCO-WARM, and we provide a simple functional form that describes its evolution with redshift. The shapes of haloes are similar in the two cases, but the smallest haloes in COCO-WARM rotate slightly more slowly than their CDM counterparts. [ABSTRACT FROM AUTHOR]

The detection of planar structures within the satellite systems of both the Milky Way (MW) and Andromeda (M31) has been reported as being in stark contradiction to the predictions of the standard cosmological model (Λ cold dark matter - ΛCDM). Given the ambiguity in defining a planar configuration, it is unclear how to interpret the low incidence of the MW and M31 planes in ΛCDM. We investigate the prevalence of satellite planes around galactic mass haloes identified in high-resolution cosmological simulations. We find that planar structures are very common, and that ~10 percent of ΛCDM haloes have even more prominent planes than those present in the Local Group. While ubiquitous, the planes of satellite galaxies show a large diversity in their properties. This precludes using one or two systems as small-scale probes of cosmology, since a large sample of satellite systems is needed to obtain a good measure of the object-to-object variation. This very diversity has been misinterpreted as a discrepancy between the satellite planes observed in the Local Group and ΛCDM predictions. In fact, ~10 percent of ΛCDM galactic haloes have planes of satellites that are as infrequent as the MW and M31 planes. The look-elsewhere effect plays an important role in assessing the detection significance of satellite planes and accounting for it leads to overestimating the significance level by a factor of 30 and 100 for the MW and M31 systems, respectively. [ABSTRACT FROM AUTHOR]

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter, we investigate how the first galaxies are assembled in a "fuzzy" dark matter (FDM) cosmology where dark matter is an ultralight ~10-22 eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical solitonlike cores, which are unstable under gravity and collapse into kiloparsec-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter--a smoking gun signature of FDM. [ABSTRACT FROM AUTHOR]