No Arabic abstract
We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with mass, from expansion to contraction, with up to factor ~10 changes in the enclosed dark matter mass at one per cent of the virial radius. Confirming previous studies, the halo response is correlated with the integrated efficiency of star formation: e_SF=(M_{star}/M_{200})/(Omega_b/Omega_m). In addition we report a new correlation with the compactness of the stellar system: e_R=r_{1/2}/R_{200}. We provide an analytic formula depending on e_SF and e_R for the response of cold dark matter haloes to baryonic processes. An observationally testable prediction is that, at fixed mass, larger galaxies experience more halo expansion, while the smaller galaxies more halo contraction. This diversity of dark halo response is captured by a toy model consisting of cycles of adiabatic inflow (causing contraction) and impulsive gas outflow (causing expansion). For net outflow, or equal inflow and outflow fractions, f, the overall effect is expansion, with more expansion with larger f. For net inflow, contraction occurs for small f (large radii), while expansion occurs for large f (small radii), recovering the phenomenology seen in our simulations. These regularities in the galaxy formation process provide a step towards a fully predictive model for the structure of cold dark matter haloes.
We perform a stacking analysis of the neutral ad,$lambdalambda$5889,5895,AA ISM doublet using the SDSS DR7 spectroscopic data set to probe the prevalence and characteristics of cold (T,$lesssim$,10$^{4}$,K) galactic-scale gas flows in local (0.025$leqslant zleqslant$0.1) inactive and AGN-host galaxies across the SFR-M$_{*}$ plane. We find low-velocity outflows to be prevalent in regions of high SFRs and stellar masses (10 $lesssim$log M$_{*}$/M$_{odot}$ $lesssim$ 11.5), however we do not find any detections in the low mass (log M$_{*}$/M$_{odot}$ $lesssim$ 10) regime. We also find tentative detections of inflowing gas in high mass galaxies across the star-forming population. We derive mass outflow rates in the range of 0.14-1.74,M$_{odot}$yr$^{-1}$ and upper limits on inflow rates <1,M$_{odot}$yr$^{-1}$, allowing us to place constraints on the mass loading factor ($eta$=$dot{M}_{text{out}}$/SFR) for use in simulations of the local Universe. We discuss the fate of the outflows by comparing the force provided by the starburst to the critical force needed to push the outflow outward, and find the vast majority of the outflows unlikely to escape the host system. Finally, as outflow detection rates and central velocities do not vary strongly with the presence of a (weak) active supermassive black hole, we determine that star formation appears to be the primary driver of outflows at $zsim$0.
We use cosmological hydrodynamical galaxy formation simulations from the NIHAO project to investigate the impact of the threshold for star formation on the response of the dark matter (DM) halo to baryonic processes. The fiducial NIHAO threshold, $n=10, {rm cm}^{-3}$, results in strong expansion of the DM halo in galaxies with stellar masses in the range $10^{7.5} < M_{star} < 10^{9.5} M_{odot}$. We find that lower thresholds such as $n=0.1$ (as employed by the EAGLE/APOSTLE and Illustris/AURIGA projects) do not result in significant halo expansion at any mass scale. Halo expansion driven by supernova feedback requires significant fluctuations in the local gas fraction on sub-dynamical times (i.e., < 50 Myr at galaxy half-light radii), which are themselves caused by variability in the star formation rate. At one per cent of the virial radius, simulations with $n=10$ have gas fractions of $simeq 0.2$ and variations of $simeq 0.1$, while $n=0.1$ simulations have order of magnitude lower gas fractions and hence do not expand the halo. The observed DM circular velocities of nearby dwarf galaxies are inconsistent with CDM simulations with $n=0.1$ and $n=1$, but in reasonable agreement with $n=10$. Star formation rates are more variable for higher $n$, lower galaxy masses, and when star formation is measured on shorter time scales. For example, simulations with $n=10$ have up to 0.4 dex higher scatter in specific star formation rates than simulations with $n=0.1$. Thus observationally constraining the sub-grid model for star formation, and hence the nature of DM, should be possible in the near future.
We use cosmological hydrodynamical galaxy formation simulations from the NIHAO project to investigate the response of cold dark matter (CDM) haloes to baryonic processes. Previous work has shown that the halo response is primarily a function of the ratio between galaxy stellar mass and total virial mass, and the density threshold above which gas is eligible to form stars, $n [{rm cm}^{-3}]$. At low $n$ all simulations in the literature agree that dwarf galaxy haloes are cuspy, but at high $nge 100$ there is no consensus. We trace halo contraction in dwarf galaxies with $nge 100$ reported in some previous simulations to insufficient spatial resolution. Provided the adopted star formation threshold is appropriate for the resolution of the simulation, we show that the halo response is remarkably stable for $nge 5$, up to the highest star formation threshold that we test, $n=500$. This free parameter can be calibrated using the observed clustering of young stars. Simulations with low thresholds $nle 1$ predict clustering that is too weak, while simulations with high star formation thresholds $nge 5$, are consistent with the observed clustering. Finally, we test the CDM predictions against the circular velocities of nearby dwarf galaxies. Low thresholds predict velocities that are too high, while simulations with $nsim 10$ provide a good match to the observations. We thus conclude that the CDM model provides a good description of the structure of galaxies on kpc scales provided the effects of baryons are properly captured.
We use N-body simulations of dark matter haloes in cold dark matter (CDM) and a large set of different warm dark matter (WDM) cosmologies to demonstrate that the spherically averaged density profile of dark matter haloes has a shape that depends on the power spectrum of matter perturbations. Density profiles are steeper in WDM but become shallower at scales less than one percent of the virial radius. Virialization isotropizes the velocity dispersion in the inner regions of the halo but does not erase the memory of the initial conditions in phase space. The location of the observed deviations from CDM in the density profile and in phase space can be directly related to the ratio between the halo mass and the filtering mass and are most evident in small mass haloes, even for a 34 keV thermal relic WDM. The rearrangement of mass within the haloes supports analytic models of halo structure that include angular momentum. We also find evidence of a dependence of the slope of the inner density profile in CDM cosmologies on the halo mass with more massive haloes exhibiting steeper profiles, in agreement with the model predictions and with previous simulation results. Our work complements recent studies of microhaloes near the filtering scale in CDM and strongly argue against a universal shape for the density profile.
We study the hot and cold circum-galactic medium (CGM) of 86 galaxies of the cosmological, hydrodynamical simulation suite NIHAO. NIHAO allows a study of how the $z=0$ CGM varies across 5 orders of magnitude of stellar mass using OVI and HI as proxies for hot and cold gas. The cool HI covering fraction and column density profiles match observations well, particularly in the inner CGM. OVI shows increasing column densities with mass, a trend seemingly echoed in the observations. As in multiple previous simulations, the OVI column densities in simulations are lower than observed and optically thick HI does not extend as far out as in observations. We take a look at the collisional ionisation fraction of OVI as a function of halo mass. We make observable predictions of the bipolarity of outflows and their effect on the general shape of the CGM. Bipolar outflows can be seen out to around 40 kpc in intermediate and low mass halos ($M_{mathrm{halo}}<10^{11}M_{mathrm{sun}}$), but outside that radius, the CGM is too well mixed to detect an elongated shape. Larger halos have extended gas discs beyond the stellar disc that dominate the shape of the inner CGM. The simulated CGM is remarkably spherical even in low mass simulations. The chemical enrichment of both halo and disc gas follow expected increasing trends as a function of halo mass that are well fit with power laws. These relations can be used in non-hydrodynamic models, such as semi-analytic models.