No Arabic abstract
We present the Millennium-II Simulation (MS-II), a very large N-body simulation of dark matter evolution in the concordance LCDM cosmology. The MS-II assumes the same cosmological parameters and uses the same particle number and output data structure as the original Millennium Simulation (MS), but was carried out in a periodic cube one-fifth the size (100 Mpc/h) with 5 times better spatial resolution (a Plummer equivalent softening of 1.0 kpc/h) and with 125 times better mass resolution (a particle mass of 6.9 times 10^6 Msun/h). By comparing results at MS and MS-II resolution, we demonstrate excellent convergence in dark matter statistics such as the halo mass function, the subhalo abundance distribution, the mass dependence of halo formation times, the linear and nonlinear autocorrelations and power spectra, and halo assembly bias. Together, the two simulations provide precise results for such statistics over an unprecedented range of scales, from halos similar to those hosting Local Group dwarf spheroidal galaxies to halos corresponding to the richest galaxy clusters. The Milky Way halos of the Aquarius Project were selected from a lower resolution version of the MS-II and were then resimulated at much higher resolution. As a result, they are present in the MS-II along with thousands of other similar mass halos. A comparison of their assembly histories in the MS-II and in resimulations of 1000 times better resolution shows detailed agreement over a factor of 100 in mass growth. We publicly release halo catalogs and assembly trees for the MS-II in the same format within the same archive as those already released for the MS.
Certain configurations of massive structures projected along the line of sight maximize the number of detections of gravitationally lensed $zsim10$ galaxies. We characterize such lines of sight with the etendue $sigma_mu$, the area in the source plane magnified over some threshold $mu$. We use the Millennium I and Millennium XXL cosmological simulations to determine the frequency of high $sigma_mu$ beams on the sky, their properties, and efficient selection criteria. We define the best beams as having $sigma_{mu>3} >2000$ arcsec$^2$, for a $zsim10$ source plane, and predict $477 pm 21$ such beams on the sky. The total mass in the beam and $sigma_{mu>3}$ are strongly correlated. After controlling for total mass, we find a significant residual correlation between $sigma_{mu>3}$ and the number of cluster-scale halos ($>10^{14} M_odot h^{-1}$) in the beam. Beams with $sigma_{mu>3} >2000$ arcsec$^2$, which should be best at lensing $zsim10$ galaxies, are ten times more likely to contain multiple cluster-scale halos than a single cluster-scale halo. Beams containing an Abell 1689-like massive cluster halo often have additional structures along the line of sight, including at least one additional cluster-scale ($M_{200}>10^{14}M_odot h^{-1}$) halo 28% of the time. Selecting beams with multiple, massive structures will lead to enhanced detection of the most distant and intrinsically faint galaxies.
We detected 10 compact galaxy groups (CGs) at $z=0$ in the semi-analytic galaxy catalog of Guo et al. (2011) for the milli-Millennium Cosmological Simulation (sCGs in mGuo2010a). We aimed to identify potential canonical pathways for compact group evolution and thus illuminate the history of observed nearby compact groups. By constructing merger trees for $z=0$ sCG galaxies, we studied the cosmological evolution of key properties, and compared them with $z=0$ Hickson CGs (HCGs). We found that, once sCG galaxies come within 1 (0.5) Mpc of their most massive galaxy, they remain within that distance until $z=0$, suggesting sCG birth redshifts. At $z=0$ stellar masses of sCG most-massive galaxies are within $10^{10} lesssim M_{ast}/M_{odot} lesssim 10^{11}$. In several cases, especially in the two 4- and 5-member systems, the amount of cold gas mass anti-correlates with stellar mass, which in turn correlates with hot gas mass. We define the angular difference between group members 3D velocity vectors, $Deltatheta_{rm vel}$, and note that many of the groups are long-lived because their small values of $Deltatheta_{rm vel}$ indicate a significant parallel component. For triplets in particular, $Deltatheta_{rm vel}$ values range between $20^{circ}$ and $40^{circ}$ so that galaxies are coming together along roughly parallel paths, and pairwise separations do not show large pronounced changes after close encounters. The best agreement between sCG and HCG physical properties is for $M_{ast}$ galaxy values, but HCG values are higher overall, including for SFRs. Unlike HCGs, due to a tail at low SFR and $M_{ast}$, and a lack of $M_{ast}gtrsim 10^{11}M_{odot}$ galaxies, only a few sCG galaxies are on the star-forming main sequence.
Kinetic Field Theory (KFT) is a statistical field theory for an ensemble of point-like classical particles in or out of equilibrium. We review its application to cosmological structure formation. Beginning with the construction of the generating functional of the theory, we describe in detail how the theory needs to be adapted to reflect the expanding spatial background and the homogeneous and isotropic, correlated initial conditions for cosmic structures. Based on the generating functional, we develop three main approaches to non-linear, late-time cosmic structures, which rest either on the Taylor expansion of an interaction operator, suitable averaging procedures for the interaction term, or a resummation of perturbation terms. We show how an analytic, parameter-free equation for the non-linear cosmic power spectrum can be derived. We explain how the theory can be used to derive the density profile of gravitationally bound structures and use it to derive power spectra of cosmic velocity densities. We further clarify how KFT relates to the BBGKY hierarchy. We then proceed to apply kinetic field theory to fluids, introduce a reformulation of KFT in terms of macroscopic quantities which leads to a resummation scheme, and use this to describe mixtures of gas and dark matter. We discuss how KFT can be applied to study cosmic structure formation with modified theories of gravity. As an example for an application to a non-cosmological particle ensemble, we show results on the spatial correlation function of cold Rydberg atoms derived from KFT.
Future galaxy surveys require realistic mock catalogues to understand and quantify systematics in order to make precise cosmological measurements. We present a halo lightcone catalogue and halo occupation distribution (HOD) galaxy catalogue built using the Millennium-XXL (MXXL) simulation. The halo catalogue covers the full sky, extending to z = 2 with a mass resolution of ~1e11 Msun/h . We use this to build a galaxy catalogue, which has an r-band magnitude limit of r < 20.0, with a median redshift of z~0.2. A Monte Carlo HOD method is used to assign galaxies to the halo lightcone catalogue, and we evolve the HODs to reproduce a target luminosity function; by construction, the luminosity function of galaxies in the mock is in agreement with the Sloan Digital Sky Survey (SDSS) at low redshifts and the Galaxy and Mass Assembly (GAMA) survey at high redshifts. A Monte Carlo method is used to assign a 0.1(g-r) colour to each galaxy, and the colour distribution of galaxies at different redshifts agrees with measurements from GAMA. The clustering of galaxies in the mock for galaxies in different magnitude and redshift bins is in good agreement with measurements from SDSS and GAMA, and the colour-dependent clustering is in reasonable agreement. We show that the baryon acoustic oscillation (BAO) can be measured in the mock catalogue, and the redshift space distortions (RSDs) are in agreement with measurements from SDSS, illustrating that this catalogue will be useful for upcoming surveys.
The Millennium N-body simulation and the Sloan Digital Sky Survey seventh data release (SDSS DR7) galaxy and galaxy group catalogues are compared to study the properties of galaxy groups and the distribution of galaxies in groups. We construct mock galaxy group catalogues for a Millennium semi-analytical galaxy catalogue by using the same friends-of-friends method, which was used by Tago et al to analyse the SDSS data. We analyse in detail the group luminosities, group richnesses, virial radii, sizes of groups and their rms velocities for four volume-limited samples from observations and simulations. Our results show that the spatial densities of groups agree within one order of magnitude in all samples with a rather good agreement between the mock catalogues and observations. All group property distributions have similar shapes and amplitudes for richer groups. For galaxy pairs and small groups, the group properties for observations and simulations are clearly different. In addition, the spatial distribution of galaxies in small groups is different: at the outskirts of the groups the galaxy number distributions do not agree, although the agreement is relatively good in the inner regions. Differences in the distributions are mainly due to the observational limitations in the SDSS sample and to the problems in the semi-analytical methods that produce too compact and luminous groups.