The cold dark matter model has become the leading theoretical picture for the formation of structure in the Universe. This model, together with the theory of cosmic inflation, makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a simulation of the growth of dark matter structure using 2,160(3) particles, following them from redshift z = 127 to the present in a cube-shaped region 2.230 billion lightyears on a side. In postprocessing, we also follow the formation and evolution of the galaxies and quasars. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with future generations of observational surveys of galaxies.
We present the results of a large library of cosmological N‐body simulations, using power‐law initial spectra. The non‐linear evolution of the matter power spectra is compared with the predictions of existing analytic scaling formulae based on the work of Hamilton et al. The scaling approach has assumed that highly non‐linear structures obey ‘stable clustering’ and are frozen in proper coordinates. Our results show that, when transformed under the self‐similarity scaling, the scale‐free spectra define a non‐linear locus that is clearly shallower than would be required under stable clustering. Furthermore, the small‐scale non‐linear power increases as both the power spectrum index n and the density parameter Ω decrease, and this evolution is not well accounted for by the previous scaling formulae. This breakdown of stable clustering can be understood as resulting from the modification of dark matter haloes by continuing mergers. These effects are naturally included in the analytic ‘halo model’ for non‐linear structure; we use this approach to fit both our scale‐free results and also our previous cold dark matter data. This method is more accurate than the commonly used Peacock–Dodds formula and should be applicable to more general power spectra. Code to evaluate non‐linear power spectra using this method is available from http://as1.chem.nottingham.ac.uk/~res/software.html. Following publication, we will make the power‐law simulation data publically available through the Virgo website http://www.mpa-garching.mpg.de/Virgo/.
We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.Comment: 15 pages, 13 figures, to be submitted to MNRAS. For high-resolution figures, please see http://www-theorie.physik.unizh.ch/~agertz
We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.
A B S T R A C TWe investigate the importance of several numerical artefacts such as lack of resolution on spectral properties of the Lya forest as computed from cosmological hydrodynamic simulations in a standard cold dark matter universe. We use a new simulation code which is based on a combination of a hierarchical particle-particle-particle-mesh (P3M) scheme for gravity and smoothed particle hydrodynamics (SPH) for gas dynamics. We have performed extensive comparisons between this new code and a modified version of the HYDRA code of Couchman et al. and find excellent agreement. We have also rerun the TREESPH simulations of Hernquist et al. using our new codes and find very good agreement with their published results. This shows that results from hydrodynamical simulations that include cooling are reproducible with different numerical algorithms. We then use our new code to investigate several numerical effects, such as resolution, on spectral statistics deduced from Voigt profile fitting of lines by running simulations with gas particle masses of 1:4 × 10 8 , 1:8 × 10 7 , 2:2 × 10 6 and 2:1 × 10 5 M ᭪ . When we increase the numerical resolution the mean effective hydrogen optical depth converges and so does the column density distribution. However, higher resolution simulations produce narrower lines and consequently the b parameter (velocity width) distribution has only marginally converged in our highest resolution run. Obtaining numerical convergence for the mean He ii transmission is demanding. When progressively smaller haloes are resolved at better resolution, a larger fraction of low-density gas contracts to moderate overdensities in which He ii is already optically thick, and this increases the net transmission, making it difficult to simulate He ii reliably. Our highest resolution simulation gives a mean effective optical depth in He ii 5 per cent lower than the simulation with eight times lower mass resolution, illustrating the degree to which the He ii optical depth has converged. In contrast, the hydrogen mean optical depth for these runs is identical. As many properties of the simulated Lya forest depend on resolution, one should be careful when deducing physical parameters from a comparison of the simulated forest with the observed one. We compare predictions from our highest resolution simulation in a cold dark matter universe, with a photoionizing background inferred from quasars as computed by Haardt & Madau, with observations. The simulation reproduces both the H I column density and b parameter distribution when we assume a high baryon density, Q B h 2 տ 0:028. In addition we need to impose a higher intergalactic medium (IGM) temperature than predicted within our basic set of assumptions. We argue that such a higher temperature could be caused by differences between the assumed and true reionization history. The simulated H I optical depth is in good agreement with observations, but the He ii optical depth is lower than observed. Fitting the He ii optical depth requires a larger jump, ϳ14, b...
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.