Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe effect of dark matter resolution on the collapse of baryons in high redshift numerical simulations
114
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We examine the impact of dark matter particle resolution on the formation of a baryonic core in high resolution adaptive mesh refinement simulations. We test the effect that both particle smoothing and particle splitting have on the hydrodynamic properties of a collapsing halo at high redshift (z > 20). Furthermore, we vary the background field intensity, with energy below the Lyman limit (< 13.6 eV), as may be relevant for the case of metal-free star formation and super-massive black hole seed formation. We find that using particle splitting methods greatly increases our particle resolution without introducing any numerical noise and allows us to achieve converged results over a wide range of external background fields. Additionally, we find that for lower values of the background field a lower dark matter particle mass is required. We define the radius of the core as the point at which the enclosed baryonic mass dominates over the enclosed dark matter mass. For our simulations this results in $rm{R_{core} sim 5 pc}$. We find that in order to produce converged results which are not affected by dark matter particles requires that the relationship ${M_{rm{core}} / M_{rm{DM}}} > 100.0$ be satisfied, where ${M_{rm{core}}}$ is the enclosed baryon mass within the core and $M_{rm{DM}}$ is the minimum dark matter particle mass. This ratio should provide a very useful starting point for conducting convergence tests before any production run simulations. We find that dark matter particle smoothing is a useful adjunct to already highly resolved simulations.
rate research
Read More
The spatial and velocity distributions of dark matter particles in the Milky Way Halo affect the signals expected to be observed in searches for dark matter. Results from direct detection experiments are often analyzed assuming a simple isothermal distribution of dark matter, the Standard Halo Model (SHM). Yet there has been skepticism regarding the validity of this simple model due to the complicated gravitational collapse and merger history of actual galaxies. In this paper we compare the SHM to the results of cosmological hydrodynamical simulations of galaxy formation to investigate whether or not the SHM is a good representation of the true WIMP distribution in the analysis of direct detection data. We examine two Milky Way-like galaxies from the MaGICC cosmological simulations (a) with dark matter only and (b) with baryonic physics included. The inclusion of baryons drives the shape of the DM halo to become more spherical and makes the velocity distribution of dark matter particles less anisotropic especially at large heliocentric velocities, thereby making the SHM a better fit. We also note that we do not find a significant disk-like rotating dark matter component in either of the two galaxy halos with baryons that we examine, suggesting that dark disks are not a generic prediction of cosmological hydrodynamical simulations. We conclude that in the Solar neighborhood, the SHM is in fact a good approximation to the true dark matter distribution in these cosmological simulations (with baryons) which are reasonable representations of the Milky Way, and hence can also be used for the purpose of dark matter direct detection calculations.
The James Webb Space Telescope (JWST) will revolutionise 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. We find that while at fixed star formation efficiency, $varepsilon$, 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 $varepsilon$ in the same dark matter model. We find that a high $varepsilon$ WDM and a low $varepsilon$ 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 $varepsilon$ 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 $varepsilon$. 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.
Self-gravitating astronomical objects often show a central plateau in the density profile (core) whose physical origin is hotly debated. Cores are theoretically expected in N-body systems of maximum entropy, however, they are not present in the canonical N-body numerical simulations of cold dark matter (CDM). Our work shows that despite this apparent contradiction between theory and numerical simulations, they are fully consistent. Simply put, cores are characteristic of systems in thermodynamic equilibrium, but thermalizing collisions are purposely suppressed in CDM simulations. When collisions are allowed, N-body numerical simulations develop cored density profiles, in perfect agreement with the theoretical expectation. We compare theory and two types of numerical simulations: (1) when DM particles are self-interacting (SIDM) with enough cross-section, then the effective two-body relaxation timescale becomes shorter than the Hubble time resulting in cored DM haloes. The haloes thus obtained, with masses from dwarf galaxies to galaxy clusters, collapse to a single shape after normalization, and this shape agrees with the polytropic density profile theoretically expected. (2) The inner radii in canonical N-body numerical simulations are always discarded because the use of finite-mass DM particles artificially increases the two-body collision rate. We show that the discarded radii develop cores that are larger than the employed numerical softening and have polytropic shapes independently of halo mass. Our work suggests that the presence of cores in simulated (or observed) density profiles can used as evidence for systems in thermodynamic equilibrium.
We study the density structures of dark matter subhalos for both cold dark matter and self-interacting dark matter models using high-resolution cosmological $N$-body simulations. We quantify subhalos central density at 150 pc from the center of each subhalo at the classical dwarf spheroidal and ultrafaint dwarf scales. By comparing them with observations, we find that the self-interacting scattering cross-section of $sigma/m<3 rm{cm^{2}g^{-1}}$ is favored. Due to the combination of hosts tide and self-interactions, the central density of subhalos with small pericenter shows a noticeable difference between the cold and the self-interacting models, indicating that dwarf satellites with small pericenter are ideal sites to further constrain the nature of dark matter by future large spectroscopic surveys.
Resolving faint galaxies in large volumes is critical for accurate cosmic reionisation simulations. While less demanding than hydrodynamical simulations, semi-analytic reionisation models still require very large N-body simulations in order to resolve the atomic cooling limit across the whole reionisation history within box sizes $gtrsim 100 , h^{-1} {rm Mpc}$. To facilitate this, we extend the mass resolution of N-body simulations using a Monte Carlo algorithm. We also propose a method to evolve positions of Monte Carlo halos, which can be an input for semi-analytic reionisation models. To illustrate, we present an extended halo catalogue that reaches a mass resolution of $M_text{halo} = 3.2 times 10^7 , h^{-1} text{M}_odot$ in a $105 , h^{-1} {rm Mpc}$ box, equivalent to an N-body simulation with $sim 6800^3$ particles. The resulting halo mass function agrees with smaller volume N-body simulations with higher resolution. Our results also produce consistent two-point correlation functions with analytic halo bias predictions. The extended halo catalogues are applied to the textsc{meraxes} semi-analytic reionisation model, which improves the predictions on stellar mass functions, star formation rate densities and volume-weighted neutral fractions. Comparison of high resolution large volume simulations with both small volume or low resolution simulations confirms that both low resolution and small volume simulations lead to reionisation ending too rapidly. Lingering discrepancies between the star formation rate functions predicted with and without our extensions can be traced to the uncertain contribution of satellite galaxies.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
