We test four commonly used astrophysical simulation codes, enzo, flash, gadget and hydra, using a suite of numerical problems with analytic initial and final states. Situations similar to the ...conditions of these tests, a Sod shock, a Sedov blast, and both a static and translating King sphere, occur commonly in astrophysics, where the accurate treatment of shocks, sound waves, supernovae explosions and collapsed haloes is a key condition for obtaining reliable validated simulations. We demonstrate that comparable results can be obtained for Lagrangian and Eulerian codes by requiring that approximately one particle exists per grid cell in the region of interest. We conclude that adaptive Eulerian codes, with their ability to place refinements in regions of rapidly changing density, are well suited to problems where physical processes are related to such changes. Lagrangian methods, on the other hand, are well suited to problems where large density contrasts occur and the physics are related to the local density itself rather than the local density gradient.
The baryon fraction of ΛCDM haloes Crain, Robert A.; Eke, Vincent R.; Frenk, Carlos S. ...
Monthly notices of the Royal Astronomical Society,
20/May , Letnik:
377, Številka:
1
Journal Article
Recenzirano
We investigate the baryon fraction in dark matter haloes formed in non-radiative gas-dynamical simulations of the Λ cold dark matter (CDM) cosmogony. By combining a realization of the Millennium ...Simulation with a simulation of a smaller volume focusing on dwarf haloes, our study spans five decades in halo mass, from 1010 to 1015 h−1 M⊙. We find that the baryon fraction within the halo virial radius is typically 90 per cent of the cosmic mean, with a rms scatter of 6 per cent, independently of redshift and of halo mass down to the smallest resolved haloes. Our results show that, contrary to the proposal of Mo et al., pre-virialization gravitational heating is unable to prevent the collapse of gas within galactic and protogalactic haloes, and confirm the need for non-gravitational feedback in order to reduce the efficiency of gas cooling and star formation in dwarf galaxy haloes. Simulations including a simple photoheating model (where a gas temperature floor of Tfloor= 2 × 104 K is imposed from z= 11) confirm earlier suggestions that photoheating can only prevent the collapse of baryons in systems with virial temperatures T200≲ 2.2 Tfloor≈ 4.4 × 104 K (corresponding to a virial mass of M200∼ 1010 h−1 M⊙ and a circular velocity of V200∼ 35 km s−1). Photoheating may thus help regulate the formation of dwarf spheroidals and other galaxies at the extreme faint end of the luminosity function, but it cannot, on its own, reconcile the abundance of sub-L★ galaxies with the vast number of dwarf haloes expected in the ΛCDM cosmogony. The lack of evolution or mass dependence seen in the baryon fraction augurs well for X-ray cluster studies that assume a universal and non-evolving baryon fraction to place constraints on cosmological parameters.
Abstract
We present a clustering comparison of 12 galaxy formation models including semi-analytic models (SAMs) and halo occupation distribution (HOD) models all run on halo catalogues and merger ...trees extracted from a single Λ cold dark matter N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the two-point correlation functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 h−1 Mpc. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs.
Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using ...a wide variety of such codes, here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of subhaloes) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced.
We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated ... cold dark matter galaxy cluster studied by Elahi et al. at z = 0. ...The 1.1 x 10... h... M... galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and sub-grid schemes. All models completed a dark-matter-only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution mDM = 9.01 x 10 super( 8) h... M... and mgas = 1.9 x 10... h... M..., respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60 per cent of haloes in the infall region with mass ~1012.5-10... h... M..., including 2-3 group-sized haloes (>10... h... M...). However, we find that only ~10 per cent of objects in the infall region are sub-haloes residing in haloes, which may suggest that there is not much ongoing pre-processing occurring in the infall region at z = 0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in sub-grid schemes and calibration procedures that each model uses. Models that do not include active galactic nucleus feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude. (ProQuest: ... denotes formulae/symbols omitted.)
Sussing merger trees: stability and convergence Wang, Yang; Pearce, Frazer R; Knebe, Alexander ...
Monthly notices of the Royal Astronomical Society,
06/2016, Letnik:
459, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Merger trees are routinely used to follow the growth and merging history of dark matter haloes and subhaloes in simulations of cosmic structure formation. Srisawat et al. compared a wide range of ...merger-tree-building codes. Here we test the influence of output strategies and mass resolution on tree-building. We find that, somewhat surprisingly, building the tree from more snapshots does not generally produce more complete trees; instead, it tends to shorten them. Significant improvements are seen for patching schemes that attempt to bridge over occasional dropouts in the underlying halo catalogues or schemes that combine the halo-finding and tree-building steps seamlessly. The adopted output strategy does not affect the average number of branches (bushiness) of the resultant merger trees. However, mass resolution has an influence on both main branch length and the bushiness. As the resolution increases, a halo with the same mass can be traced back further in time and will encounter more small progenitors during its evolutionary history. Given these results, we recommend that, for simulations intended as precursors for galaxy formation models where of the order of 100 or more snapshots are analysed, the tree-building routine should be integrated with the halo finder, or at the very least be able to patch over multiple adjacent snapshots.
Using the Millennium N-body Simulation we explore how the shape and angular momentum of galaxy dark matter haloes surrounding the largest cosmological voids are oriented. We find that the major and ...intermediate axes of the haloes tend to lie parallel to the surface of the voids, whereas the minor axis points preferentially in the radial direction. We have quantified the strength of these alignments at different radial distances from the void centres. The effect of these orientations is still detected at distances as large as 2.2 Rvoid from the void centre. Taking a subsample of haloes expected to contain disc-dominated galaxies at their centres we detect, at the 99.9 per cent confidence level, a signal that the angular momentum of those haloes tends to lie parallel to the surface of the voids. Contrary to the alignments of the inertia axes, this signal is only detected in shells at the void surface (1 < R < 1.07 Rvoid) and disappears at larger distances. This signal, together with the similar alignment observed using real spiral galaxies, strongly supports the prediction of the Tidal Torque theory that both dark matter haloes and baryonic matter have acquired, conjointly, their angular momentum before the moment of turnaround.
Cosmological simulations of the intracluster medium Kay, Scott T.; Thomas, Peter A.; Jenkins, Adrian ...
Monthly notices of the Royal Astronomical Society,
December 2004, Letnik:
355, Številka:
4
Journal Article
Recenzirano
Odprti dostop
Abstract
We investigate the properties of the intracluster medium (ICM) that forms within N-body/hydrodynamical simulations of galaxy clusters in a ΛCDM cosmology. When radiative cooling and a simple ...model for galactic feedback are included, our clusters have X-ray luminosities and temperatures in good agreement with observed systems, demonstrating the required excess entropy in their cores. More generally, cooling and feedback increases the entropy of the ICM everywhere, albeit without significantly affecting the slope of the profile (S∝r) at large radii. The temperature of the ICM is only modestly increased by these processes, with projected temperature profiles being in reasonable agreement with the observations. Star/galaxy formation is still too efficient in our simulations, however, and so our gas mass fractions are around 60 per cent of the observed value at r
2500. Finally, we examine the reliability of using the hydrostatic equilibrium equation to estimate cluster masses and find that it underpredicts the true mass of our clusters by up to 20 per cent, due to incomplete thermalization of the gas. Feedback reduces this discrepancy, however, with estimates being accurate to within 10 per cent out to r
500.