Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNIHAO XX: The impact of the star formation threshold on the cusp-core transformation of cold dark matter haloes
65
0
0.0
(
0
)
Authors
Aaron A. Dutton
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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 quantify the impact of galaxy formation on dark matter halo shapes using cosmological simulations at redshift $z=0$. The haloes are drawn from the IllustrisTNG project, a suite of magneto-hydrodynamic simulations of galaxies. We focus on haloes of mass $10^{10-14} M_odot$ from the 50-Mpc (TNG50) and 100-Mpc (TNG100) boxes, and compare them to dark matter-only (DMO) analogues and other simulations e.g. NIHAO and Eagle. We further quantify the prediction uncertainty by varying the baryonic feedback models in a series of smaller 25 Mpc $h^{-1}$ boxes. We find that: (i) galaxy formation results in rounder haloes compared to the DMO simulations, in qualitative agreement with past hydrodynamic models. Haloes of mass $approx 2times 10^{12} M_odot$ are most spherical, with an average minor-to-major axis ratio of $left< s right> approx 0.75$ in the inner halo, an increase of 40 per cent compared to their DMO counterparts. No significant change in halo shape is found for low-mass $10^{10} M_odot$ haloes; (ii) stronger feedback, e.g. increasing galactic wind speed, reduces the impact of baryons; (iii) the inner halo shape correlates with the stellar mass fraction, which can explain the dependence of halo shapes on different feedback models; (iv) the fiducial and weaker feedback models are most consistent with observational estimates of the Milky Way halo shape. Yet, at fixed halo mass, very diverse and possibly unrealistic feedback models all predict inner halo shapes that are closer to one another than to the DMO results. This implies that a larger observational sample would be required to statistically distinguish between different baryonic prescriptions due to large halo-to-halo variation in halo shapes.
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 address the issue of numerical convergence in cosmological smoothed particle hydrodynamics simulations using a suite of runs drawn from the EAGLE project. Our simulations adopt subgrid models that produce realistic galaxy populations at a fiducial mass and force resolution, but systematically vary the latter in order to study their impact on galaxy properties. We provide several analytic criteria that help guide the selection of gravitational softening for hydrodynamical simulations, and present results from runs that both adhere to and deviate from them. Unlike dark matter-only simulations, hydrodynamical simulations exhibit a strong sensitivity to gravitational softening, and care must be taken when selecting numerical parameters. Our results--which focus mainly on star formation histories, galaxy stellar mass functions and sizes--illuminate three main considerations. First, softening imposes a minimum resolved escape speed, $v_epsilon$, due to the binding energy between gas particles. Runs that adopt such small softening lengths that $v_epsilon gt 10,{rm km s^{-1}}$ (the sound speed in ionised $sim 10^4,{rm K}$ gas) suffer from reduced effects of photo-heating. Second, feedback from stars or active galactic nuclei may suffer from numerical over-cooling if the gravitational softening length is chosen below a critical value, $epsilon_{rm eFB}$. Third, we note that small softening lengths exacerbate the segregation of stars and dark matter particles in halo centres, often leading to the counter-intuitive result that galaxy sizes {em increase} as softening is reduced. The structure of dark matter haloes in hydrodynamical runs respond to softening in a way that reflects the sensitivity of their galaxy populations to numerical parameters.
The Local Group is a unique environment in which to study the astrophysics of galaxy formation. The proximity of the Milky Way and M31 causes a large fraction of the low-mass halo population to interact with more massive dark matter haloes, which increases their concentrations and strips them of gas and other material. Some low-mass haloes pass through the haloes of the Milky Way or M31 and are either ejected into the field or exchanged between the two primary hosts. We use high resolution gas-dynamical simulations to describe a new class of field halo that passed through the haloes of both the Milky Way and M31 at early times and is almost twice as concentrated as isolated field haloes. These Hermeian haloes are distributed anisotropically at greater distances from the Local Group barycentre than the primary haloes and appear to cluster close to the Milky Way and M31 in projection. We show that some Hermeian haloes can host galaxies that are promising targets for indirect dark matter searches and are competitive with signals from other dwarf galaxies. Hermeian galaxies in the Local Group should be detectable by forthcoming wide-field imaging surveys.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
