We present the evolution of the structure of relaxed cold dark matter (CDM) haloes in the cosmology from the Planck satellite. Our simulations cover five decades in halo mass, from dwarf galaxies to ...galaxy clusters. Because of the increased matter density and power spectrum normalization the concentration–mass relation in the Planck
cosmology has a ∼20 per cent higher normalization at redshift z = 0 compared to Wilkinson Microwave Anisotropy Probe cosmology. We confirm that CDM haloes are better described by the Einasto profile; for example, at scales near galaxy half-light radii CDM haloes have significantly steeper density profiles than implied by Navarro–Frenk–White (NFW) fits. There is a scatter of ∼0.2 dex in the Einasto shape parameter at fixed halo mass, adding further to the diversity of CDM halo profiles. The evolution of the concentration–mass relation in our simulations is not reproduced by any of the analytic models in the literature. We thus provide a simple fitting formula that accurately describes the evolution between redshifts z = 5 and 0 for both NFW and Einasto fits. Finally, the observed concentrations and halo masses of spiral galaxies, groups and clusters of galaxies at low redshifts are in good agreement with our simulations, suggesting only mild halo response to galaxy formation on these scales.
ABSTRACT
We study ultra-diffuse galaxies (UDGs) in zoom in cosmological simulations, seeking the origin of UDGs in the field versus galaxy groups. We find that while field UDGs arise from dwarfs in a ...characteristic mass range by multiple episodes of supernova feedback (Di Cintio et al.), group UDGs may also form by tidal puffing up and they become quiescent by ram-pressure stripping. The field and group UDGs share similar properties, independent of distance from the group centre. Their dark-matter haloes have ordinary spin parameters and centrally dominant dark-matter cores. Their stellar components tend to have a prolate shape with a Sérsic index n ∼ 1 but no significant rotation. Ram pressure removes the gas from the group UDGs when they are at pericentre, quenching star formation in them and making them redder. This generates a colour/star-formation-rate gradient with distance from the centre of the dense environment, as observed in clusters. We find that ∼20 per cent of the field UDGs that fall into a massive halo survive as satellite UDGs. In addition, normal field dwarfs on highly eccentric orbits can become UDGs near pericentre due to tidal puffing up, contributing about half of the group-UDG population. We interpret our findings using simple toy models, showing that gas stripping is mostly due to ram pressure rather than tides. We estimate that the energy deposited by tides in the bound component of a satellite over one orbit can cause significant puffing up provided that the orbit is sufficiently eccentric. We caution that while the simulations produce UDGs that match the observations, they under-produce the more compact dwarfs in the same mass range, possibly because of the high threshold for star formation or the strong feedback.
Abstract
Simulating thin and extended galactic disks has long been a challenge in computational astrophysics. We introduce the NIHAO-UHD suite of cosmological hydrodynamical simulations of Milky Way ...mass galaxies and study stellar disk properties such as stellar mass, size and rotation velocity which agree well with observations of the Milky Way and local galaxies. In particular, the simulations reproduce the age-velocity dispersion relation and a multi-component stellar disk as observed for the Milky Way. Half of our galaxies show a double exponential vertical profile, while the others are well described by a single exponential model which we link to the disk merger history. In all cases, mono-age populations follow a single exponential whose scale height varies monotonically with stellar age and radius. The scale length decreases with stellar age while the scale height increases. The general structure of the stellar disks is already set at time of birth as a result of the inside-out and upside-down formation. Subsequent evolution modifies this structure by increasing both the scale length and height of all mono-age populations. Thus, our results put tight constraints on how much dynamical memory stellar disks can retain over cosmological timescales. Our simulations demonstrate that it is possible to form thin galactic disks in cosmological simulations provided there are no significant stellar mergers at low redshifts. Most of the stellar mass is formed in-situ with only a few percent ($\lesssim 5\%$) brought in by merging satellites at early times. Redshift zero snapshots and halo catalogues are publicly available.
Recent studies have shown that massive elliptical galaxies have total mass density profiles within an effective radius that can be approximated as
, with mean slope 〈γ′〉 = 2.08 ± 0.03 and scatter
. ...The small scatter of the slope (known as the bulge-halo conspiracy) is not generic in Λ cold dark matter (ΛCDM) based models and therefore contains information about the galaxy formation process. We compute the distribution of γ′ for ΛCDM-based models that reproduce the observed correlations between stellar mass, velocity dispersion, and effective radius of early-type galaxies in the Sloan Digital Sky Survey. The models have a range of stellar initial mass functions (IMFs) and dark halo responses to galaxy formation. The observed distribution of γ′ is well reproduced by a model with cosmologically motivated but uncontracted dark matter haloes, and a Salpeter-type IMF. Other models are on average ruled out by the data, even though they may happen in individual cases. Models with adiabatic halo contraction (and lighter IMFs) predict too small values of γ′. Models with halo expansion, or mass-follows-light predict too high values of γ′. Our study shows that the non-homologous structure of massive early-type galaxies can be precisely reproduced by ΛCDM models if the IMF is not universal and if mechanisms, such as feedback from active galactic nuclei, or dynamical friction, effectively on average counterbalance the contraction of the halo expected as a result of baryonic cooling.
Abstract
We address the origin of ultra-diffuse galaxies (UDGs), which have stellar masses typical of dwarf galaxies but effective radii of Milky Way-sized objects. Their formation mechanism, and ...whether they are failed L
⋆ galaxies or diffuse dwarfs, are challenging issues. Using zoom-in cosmological simulations from the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project, we show that UDG analogues form naturally in dwarf-sized haloes due to episodes of gas outflows associated with star formation. The simulated UDGs live in isolated haloes of masses 1010–11 M⊙, have stellar masses of 107–8.5 M⊙, effective radii larger than 1 kpc and dark matter cores. They show a broad range of colours, an average Sérsic index of 0.83, a typical distribution of halo spin and concentration, and a non-negligible H i gas mass of 107 − 9 M⊙, which correlates with the extent of the galaxy. Gas availability is crucial to the internal processes which form UDGs: feedback-driven gas outflows, and subsequent dark matter and stellar expansion, are the key to reproduce faint, yet unusually extended, galaxies. This scenario implies that UDGs represent a dwarf population of low surface brightness galaxies and should exist in the field. The largest isolated UDGs should contain more H i gas than less extended dwarfs of similar M
⋆.
ABSTRACT We present a summary of recent assessments of the mass distribution in disk galaxies. In order to characterize galaxy formation models the relative fraction of baryons and dark matter at all ...radii in galaxies must be determined. For disk galaxies, various measurements of the mass distribution in galaxies based on vertical kinematics, strong lensing, residuals of scaling relations, fluid dynamical modeling, bar strength and pattern speed, warps, and others have called for either a maximal or sub-maximal contribution of the baryons in the inner parts of the disk. We propose a global picture whereby all galaxies are typically baryon-dominated (maximal) at the center and dark-matter dominated (sub-maximal) in their outskirts. Using as a fiducial radius the peak of the rotation curve of a pure baryonic exponential disk at , where is the disk scale length, the transition from maximal to sub-maximal baryons occurs within or near for low-mass disk galaxies (with 200 ) and beyond for more massive systems. The mass fraction's inverse dependence on circular velocity in disk galaxies is largely insensitive to the presence of a bar at . The mean mass fractions of late- and early-type galaxies are shown to be different at the same fiducial radius and circular velocity, pointing to different galaxy formation mechanisms. Feedback, dynamical friction, size, and IMF variations are key ingredients for understanding these differences.
Abstract
Most of the baryons in the Universe are not in the form of stars and cold gas in galaxies. Galactic outflows driven by supernovae/stellar winds are the leading mechanisms for explaining this ...fact. The scaling relation between galaxy mass and outer rotation velocity (also known as the baryonic Tully-Fisher relation, BTF) has recently been used as evidence against this viewpoint. We use a Λ cold dark matter (ΛCDM)-based semi-analytic disc galaxy formation model to investigate these claims. In our model, galaxies with less efficient star formation and higher gas fractions are more efficient at ejecting gas from galaxies. This somewhat counter intuitive result is due to the (observational) fact that galaxies with less efficient star formation and higher gas fractions tend to live in dark matter haloes with lower circular velocities, from which less energy is required to escape the potential well. In our model the intrinsic scatter in the BTF is ≃0.15 dex, and mostly reflects scatter in dark halo concentration. The scatter is largely independent of galaxy structure because of the large radius within which galaxy rotation velocities are measured. The observed scatter, equal to ≃0.24 dex, is dominated by measurement errors. The best estimate for the intrinsic scatter is that it is less than 0.15 dex, and thus our ΛCDM-based model (which does not include all possible sources of scatter) is only just consistent with this. Future observations of the BTF scatter could be made with a more stringent measurement of the intrinsic scatter, and thus provide a strong constraint to galaxy formation models. In our model, gas-rich galaxies, at fixed virial velocity (V
vir), with lower stellar masses have lower baryonic masses. This is consistent with the expectation that galaxies with lower stellar masses have had less energy available to drive an outflow. However, when the outer rotation velocity (V
flat) is used the correlation has the opposite sign, with a slope in agreement with observations. This is due to the fact there is scatter in the relation between V
flat and V
vir. In summary, contrary to some previous claims, we show that basic features of the BTF are consistent with a ΛCDM-based model in which the low efficiency of galaxy formation is determined by galactic outflows.
Abstract
State-of-the-art cosmological hydrodynamical simulations of galaxy formation have reached the point at which their outcomes result in galaxies with ever more realism. Still, the employed ...sub-grid models include several free parameters such as the density threshold, n, to localize the star-forming gas. In this work, we investigate the possibilities to utilize the observed clustered nature of star formation (SF) in order to refine SF prescriptions and constrain the density threshold parameter. To this end, we measure the clustering strength, correlation length, and power-law index of the two-point correlation function of young (τ < 50 Myr) stellar particles and compare our results to observations from the HST Legacy Extragalactic UV Survey (LEGUS). Our simulations reveal a clear trend of larger clustering signal and power-law index and lower correlation length as the SF threshold increases with only mild dependence on galaxy properties such as stellar mass or specific star formation rate. In conclusion, we find that the observed clustering of SF is inconsistent with a low threshold for SF (n < 1 cm−3) and strongly favours a high value for the density threshold of SF (n > 10 cm−3), as, for example employed in the NIHAO project.
We use a suite of 31 simulated galaxies drawn from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter haloes. The sample covers a wide mass ...range: 9.4 × 109 < M
halo/ M < 7.8 × 1011, hosting galaxies with stellar masses in the range 5.0 × 105 < M
*/ M < 8.3 × 1010, i.e. from dwarf to L*. The galaxies are simulated with blastwave supernova feedback and, for some of them, an additional source of energy from massive stars is included. Within this feedback scheme we vary several parameters, such as the initial mass function, the density threshold for star formation, and energy from supernovae and massive stars. The main result is a clear dependence of the inner slope of the dark matter density profile, α in ρ ∝ r
α, on the stellar-to-halo mass ratio, M
*/M
halo. This relation is independent of the particular choice of parameters within our stellar feedback scheme, allowing a prediction for cusp versus core formation. When M
*/M
halo is low, 0.01 per cent, energy from stellar feedback is insufficient to significantly alter the inner dark matter density, and the galaxy retains a cuspy profile. At higher stellar-to-halo mass ratios, feedback drives the expansion of the dark matter and generates cored profiles. The flattest profiles form where M
*/M
halo ∼ 0.5 per cent. Above this ratio, stars formed in the central regions deepen the gravitational potential enough to oppose the supernova-driven expansion process, resulting in cuspier profiles. Combining the dependence of α on M
*/M
halo with the empirical abundance matching relation between M
* and M
halo provides a prediction for how α varies as a function of stellar mass. Further, using the Tully-Fisher relation allows a prediction for the dependence of the dark matter inner slope on the observed rotation velocity of galaxies. The most cored galaxies are expected to have V
rot ∼ 50 km s−1, with α decreasing for more massive disc galaxies: spirals with V
rot ∼ 150 km s−1 have central slopes α ≤ −0.8, approaching again the Navarro-Frenk-White profile. This novel prediction for the dependence of α on disc galaxy mass can be tested using observational data sets and can be applied to theoretical modelling of mass profiles and populations of disc galaxies.