Lyman break galaxies at z= 4–6 in cosmological smoothed particle hydrodynamics simulations

We perform a spectrophotometric analysis of galaxies at redshifts z = 4–6 in cosmological smoothed particle hydrodynamics simulations of a Λ cold dark matter universe. Our models include radiative cooling and heating by a uniform ultraviolet (UV) background, star formation, supernova feedback, and a phenomenological model for galactic winds. Analysing a series of simulations of varying box size and particle number allows us to isolate the impact of numerical resolution on our results. Specifically, we determine the luminosity functions in B , V , R , i ' and z ' filters, and compare the results with observational surveys of Lyman break galaxies (LBGs) performed with the Subaru telescope and the Hubble Space Telescope . We find that the simulated galaxies have UV colours consistent with observations and fall in the expected region of the colour–colour diagrams used by the Subaru group. The stellar masses of the most massive galaxies in our largest simulation increase their stellar mass from   M _★∼ 10¹¹ M_⊙  at z = 6 to   M _★∼ 10^11.7 M_⊙  at z = 3. Assuming a uniform extinction of E ( B − V ) = 0.15, we also find reasonable agreement between simulations and observations in the space density of UV bright galaxies at z = 3–6, down to the magnitude limit of each survey. For the same moderate extinction level of E ( B − V ) ∼ 0.15, the simulated luminosity functions match observational data, but have a steep faint-end slope with α∼−2.0. We discuss the implications of the steep faint-end slope found in the simulations. Our results confirm the generic conclusion from earlier numerical studies that UV bright LBGs at z ≥ 3 are the most massive galaxies with E ( B − V ) ∼ 0.15 at each epoch.     

Two-body relaxation in cosmological simulations

It is logically possible that early two-body relaxation in simulations of cosmological clustering influences the final structure of massive clusters. Convergence studies in which mass and spatial resolution are simultaneously increased cannot eliminate this possibility. We test the importance of two-body relaxation in cosmological simulations with simulations in which there are two species of particles. The cases of two mass ratios, √2:1 and 4:1, are investigated. Simulations are run with both a spatially fixed softening length and adaptive softening using the publicly available codes gadget and mlapm , respectively.
The effects of two-body relaxation are detected in both the density profiles of haloes and the mass function of haloes. The effects are more pronounced with a fixed softening length, but even in this case they are not so large as to suggest that results obtained with one mass species are significantly affected by two-body relaxation.
The simulations that use adaptive softening are less affected by two-body relaxation and produce slightly higher central densities in the largest haloes. They run about three times faster than the simulations that use a fixed softening length.     

Gas in simulations of high-redshift galaxies and minihaloes

We study the gas content of haloes in the early universe using high-resolution hydrodynamical simulations. We extract from the simulations and also predict, based on linear theory, the halo mass for which the enclosed baryon fraction equals half of the mean cosmic fraction. We find a rough agreement between the simulations and the predictions, which suggests that during the high-redshift era before stellar heating, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was  ∼3 × 10⁴ M_⊙  . We also carry out a detailed resolution analysis and show that in order to determine a halo's gas fraction even to 20 per cent accuracy, the halo must be resolved into at least 500 dark matter particles.     

ΛCDM predictions for galaxy protoclusters – I. The relation between galaxies, protoclusters and quasars at z∼ 6

Motivated by recent observational studies of the environment of   z ∼ 6  QSOs, we have used the Millennium Run (MR) simulations to construct a very large  (∼4°× 4°)  mock redshift survey of star-forming galaxies at   z ∼ 6  . We use this simulated survey to study the relation between density enhancements in the distribution of i ₇₇₅-dropouts and Lyα emitters, and their relation to the most massive haloes and protocluster regions at   z ∼ 6  . Our simulation predicts significant variations in surface density across the sky with some voids and filaments extending over scales of 1°, much larger than probed by current surveys. Approximately one-third of all   z ∼ 6  haloes hosting i -dropouts brighter than   z = 26.5  mag  (≈ M *_{UV, z =6})  become part of   z = 0  galaxy clusters. i -dropouts associated with protocluster regions are found in regions where the surface density is enhanced on scales ranging from a few to several tens of arcminutes on the sky. We analyse two structures of i -dropouts and Lyα emitters observed with the Subaru Telescope and show that these structures must be the seeds of massive clusters in formation. In striking contrast, six   z ∼ 6  QSO fields observed with Hubble Space Telescope show no significant enhancements in their i ₇₇₅-dropout number counts. With the present data, we cannot rule out the QSOs being hosted by the most massive haloes. However, neither can we confirm this widely used assumption. We conclude by giving detailed recommendations for the interpretation and planning of observations by current and future ground- and space-based instruments that will shed new light on questions related to the large-scale structure at   z ∼ 6  .     

Alignment between luminous red galaxies and surrounding structures at z∼ 0.5

We analyse a high-redshift sample (0.4 < z < 0.5) of luminous red galaxies (LRGs) extracted from the Sloan Digital Sky Survey data release 4 and their surrounding structures to explore the presence of alignment effects of these bright galaxies with neighbour objects. In order to avoid projection effects, we compute photometric redshifts for galaxies within 3  h ⁻¹ Mpc in projection of LRGs and calculate the relative angle between the LRG major axis and the direction to neighbours within 1000 km s⁻¹. We find a clear signal of alignment between LRG orientations and the distribution of galaxies within 1.5  h ⁻¹ Mpc. The alignment effects are present only for the red population of tracers; LRG orientation is uncorrelated to the blue population of neighbour galaxies. These results add evidence to the alignment effects between primaries and satellites detected at low redshifts. We conclude that such alignments were already present at z ∼ 0.5.     

Parameter tests within cosmological simulations of galaxy formation

Numerical simulations of galaxy formation require a number of parameters. Some of these are intrinsic to the numerical integration scheme (e.g., the time-step), while others describe the physical model (e.g., the gas metallicity). In this paper we present results of a systematic exploration of the effects of varying a subset of these parameters on simulations of galaxy formation. We use N -body and 'Smoothed Particle Hydrodynamics' techniques to follow the evolution of cold dark matter and gas in a small volume. We compare a fiducial model with 24 different simulations, in which one parameter at a time is varied, focusing on properties such as the relative fraction of hot and cold gas, and the abundance and masses of galaxies. We find that for reasonable choices of numerical values, many parameters have relatively little effect on the galaxies, with the notable exception of the parameters that control the resolution of the simulation and the efficiency with which gas cools.     

Satellite systems around galaxies in hydrodynamic simulations

We investigate the properties of satellite galaxies formed in N -body/SPH simulations of galaxy formation in the ΛCDM cosmology. The simulations include the main physical effects thought to be important in galaxy formation and, in several cases, produce realistic spiral discs. In total, a sample of nine galaxies of luminosity comparable to the Milky Way was obtained. At magnitudes brighter than the resolution limit,   M_V =−12  , the luminosity function of the satellite galaxies in the simulations is in excellent agreement with data for the Local Group. The radial number density profile of the model satellites, as well as their gas fractions also match observations very well. In agreement with previous N -body studies, we find that the satellites tend to be distributed in highly flattened configurations whose major axis is aligned with the major axis of the (generally triaxial) dark halo. In two out of three systems with sufficiently large satellite populations, the satellite system is nearly perpendicular to the plane of the galactic disc, a configuration analogous to that observed in the Milk Way. The discs themselves are perpendicular to the minor axis of their host haloes in the inner parts, and the correlation between the orientation of the galaxy and the shape of the halo persists even out to the virial radius. However, in one case the disc's minor axis ends up, at the virial radius, perpendicular to the minor axis of the halo. The angular momenta of the galaxies and their host halo tend to be well aligned.     

Multi-level adaptive particle mesh (MLAPM): a c code for cosmological simulations

We present a computer code written in c that is designed to simulate structure formation from collisionless matter. The code is purely grid-based and uses a recursively refined Cartesian grid to solve Poisson's equation for the potential, rather than obtaining the potential from a Green's function. Refinements can have arbitrary shapes and in practice closely follow the complex morphology of the density field that evolves. The time-step shortens by a factor of 2 with each successive refinement.
Competing approaches to N -body simulation are discussed from the point of view of the basic theory of N -body simulation. It is argued that an appropriate choice of softening length ε is of great importance and that ε should be at all points an appropriate multiple of the local interparticle separation. Unlike tree and P³M codes, multigrid codes automatically satisfy this requirement. We show that at early times and low densities in cosmological simulations, ε needs to be significantly smaller relative to the interparticle separation than in virialized regions. Tests of the ability of the code's Poisson solver to recover the gravitational fields of both virialized haloes and Zel'dovich waves are presented, as are tests of the code's ability to reproduce analytic solutions for plane-wave evolution. The times required to conduct a ΛCDM cosmological simulation for various configurations are compared with the times required to complete the same simulation with the ART, AP³M and GADGET codes. The power spectra, halo mass functions and halo–halo correlation functions of simulations conducted with different codes are compared.
The code is available from http://www-thphys.physics.ox.ac.uk/users/MLAPM .     

The properties of Lyα emitting galaxies in hierarchical galaxy formation models

We present detailed predictions for the properties of Lyα-emitting galaxies in the framework of the Λ cold dark matter cosmology, calculated using the semi-analytical galaxy formation model galform . We explore a model that assumes a top-heavy initial mass function in starbursts and that has previously been shown to explain the sub-millimetre number counts and the luminosity function of Lyman-break galaxies at high redshift. We show that this model, with the simple assumption that a fixed fraction of Lyα photons escape from each galaxy, is remarkably successful at explaining the observed luminosity function of Lyα emitters (LAEs) over the redshift range  3 < z < 6.6  . We also examine the distribution of Lyα equivalent widths and the broad-band continuum magnitudes of emitters, which are in good agreement with the available observations. We look more deeply into the nature of LAEs, presenting predictions for fundamental properties such as the stellar mass and radius of the emitting galaxy and the mass of the host dark matter halo. The model predicts that the clustering of LAEs at high redshifts should be strongly biased relative to the dark matter, in agreement with observational estimates. We also present predictions for the luminosity function of LAEs at   z > 7  , a redshift range that is starting to be be probed by near-infrared surveys and using new instruments such as the Dark Ages Z Lyman Explorer (DAzLE).