Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn efficient hybrid method to produce high resolution large volume dark matter simulations for semi-analytic models of reionisation
58
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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 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.
We study the correlation between the specific star formation rate of central galaxies and neighbour galaxies, also known as galactic conformity, out to 20 Mpc/h using three semi-analytic models (SAMs, one from L-GALAXIES and other two from GALFORM). The aim is to establish whether SAMs are able to show galactic conformity using different models and selection criteria. In all the models, when the selection of primary galaxies is based on an isolation criterion in real space, the mean fraction of quenched galaxies around quenched primary galaxies is higher than that around star-forming primary galaxies of the same stellar mass. The overall signal of conformity decreases when we remove satellites selected as primary galaxies, but the effect is much stronger in GALFORM models compared with the L-GALAXIES model. We find this difference is partially explained by the fact that in GALFORM once a galaxy becomes a satellite remains as such, whereas satellites can become centrals at a later time in L-GALAXIES. The signal of conformity decreases down to 60% in the L-GALAXIES model after removing central galaxies that were ejected from their host halo in the past. Galactic conformity is also influenced by primary galaxies at fixed stellar mass that reside in dark matter haloes of different masses. Finally, we explore a proxy of conformity between distinct haloes. In this case the conformity is weak beyond ~ 3 Mpc/h (<3% in L-GALAXIES, <1-2% in GALFORM models). Therefore, it seems difficult that conformity is directly related with a long-range effect.
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.
We argue that dark radiation is naturally generated from the decay of the overall volume modulus in the LARGE volume scenario. We consider both sequestered and non-sequestered cases, and find that the axionic superpartner of the modulus is produced by the modulus decay and it can account for the dark radiation suggested by observations, while the modulus decay through the Giudice-Masiero term gives the dominant contribution to the total decay rate. In the sequestered case, the lightest supersymmetric particles produced by the modulus decay can naturally account for the observed dark matter density. In the non-sequestered case, on the other hand, the supersymmetric particles are not produced by the modulus decay, since the soft masses are of order the heavy gravitino mass. The QCD axion will then be a plausible dark matter candidate.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
