We present the measurements and modelling of the small-to-intermediate scale (0.1--25 Mpc/h) projected and three-dimensional (3D) redshift-space two-point correlation functions (2PCFs) of local galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7. We find a clear dependence of galaxy clustering on luminosity in both projected and redshift spaces, generally being stronger for more luminous samples. The measurements are successfully interpreted within the halo occupation distribution (HOD) framework with central and satellite velocity bias parameters to describe galaxy kinematics inside haloes and to model redshift-space distortion (RSD) effects. In agreement with previous studies, we find that more luminous galaxies reside in more massive haloes. Including the redshift-space 2PCFs helps tighten the HOD constraints. Moreover, we find that luminous central galaxies are not at rest at the halo centres, with the velocity dispersion about 30% that of the dark matter. Such a relative motion may reflect the consequence of galaxy and halo mergers, and we find that central galaxies in lower mass haloes tend to be more relaxed with respect to their host haloes. The motion of satellite galaxies in luminous samples is consistent with their following that of the dark matter. For faint samples, satellites tends to have slower motion, with velocity dispersion inside haloes about 85% that of the dark matter. We discuss possible applications of the velocity bias constraints on studying galaxy evolution and cosmology. In the appendix, we characterize the distribution of galaxy redshift measurement errors, which is well described by a Gaussian-convolved double exponential distribution.

It is well known that satellite galaxies are not isotropically distributed among their host galaxies as suggested by most interpretations of the $\Lambda$CDM model. One type of anisotropy recently detected in the SDSS (and seen when examining the distribution of satellites in the Local Group and in the Centaurus group) is a tendency to be so-called "lopsided". Namely, in pairs of galaxies (like Andromeda and the Milky Way) the satellites are more likely to inhabit the region in between the pair, rather than on opposing sides. Although recent studies found a similar set up when comparing pairs of galaxies in $\Lambda$CDM simulations indicating that such a set up is not inconsistent with $\Lambda$CDM, the origin has yet to be explained. Here we examine the origin of such lopsided setups by first identifying such distributions in pairs of galaxies in numerical cosmological simulations, and then tracking back the orbital trajectories of satellites (which at $z=0$ display the effect). We report two main results: first, the lopsided distribution was stronger in the past and weakens towards $z=0$. Second, the weakening of the signal is due to the interaction of satellite galaxies with the pair. Finally, we show that the $z=0$ signal is driven primarily by satellites that are on first approach, who have yet to experience a "flyby". This suggests that the signal seen in the observations is also dominated by dynamically young accretion events.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

Today's galaxies experienced cosmic reionization at different times in different locations. For the first time, reionization ($50\%$ ionized) redshifts, $z_R$, at the location of their progenitors are derived from new, fully-coupled radiation-hydrodynamics simulation of galaxy formation and reionization at $z > 6$, matched to N-body simulation to z = 0. Constrained initial conditions were chosen to form the well-known structures of the local universe, including the Local Group and Virgo, in a (91 Mpc)$^3$ volume large enough to model both global and local reionization. Reionization simulation CoDa I-AMR, by CPU-GPU code EMMA, used (2048)$^3$ particles and (2048)$^3$ initial cells, adaptively-refined, while N-body simulation CoDa I-DM2048, by Gadget2, used (2048)$^3$ particles, to find reionization times for all galaxies at z = 0 with masses $M(z=0)\ge 10^8 M_\odot$. Galaxies with $M(z=0) \gtrsim 10^{11} M_\odot$ reionized earlier than the universe as a whole, by up to $\sim$ 500 Myrs, with significant scatter. For Milky-Way-like galaxies, $z_R$ ranged from 8 to 15. Galaxies with $M(z=0) \lesssim 10^{11} M_\odot$ typically reionized as late or later than globally-averaged $50\%$ reionization at $\langle z_R\rangle =7.8$, in neighborhoods where reionization was completed by external radiation. The spread of reionization times within galaxies was sometimes as large as the galaxy-to-galaxy scatter. The Milky Way and M31 reionized earlier than global reionization but later than typical for their mass, neither dominated by external radiation. Their most massive progenitors at $z>6$ had $z_R$ = 9.8 (MW) and 11 (M31), while their total masses had $z_R$ = 8.2 (both).

The Zone of Avoidance (ZOA), whose emptiness is an artifact of our Galaxy dust, has been challenging observers as well as theorists for many years. Multiple attempts have been made on the observational side to map this region in order to better understand the local flows. On the theoretical side, however, this region is often simply statistically populated with structures but no real attempt has been made to confront theoretical and observed matter distributions. This paper takes a step forward using constrained realizations of the local Universe shown to be perfect substitutes of local Universe-like simulations for smoothed high density peak studies. Far from generating completely `random' structures in the ZOA, the reconstruction technique arranges matter according to the surrounding environment of this region. More precisely, the mean distributions of structures in a series of constrained and random realizations differ: while densities annihilate each other when averaging over 200 random realizations, structures persist when summing 200 constrained realizations. The probability distribution function of ZOA grid cells to be highly overdense is a Gaussian with a 15% mean in the random case, while that of the constrained case exhibits large tails. This implies that areas with the largest probabilities host most likely a structure. Comparisons between these predictions and observations, like those of the Puppis 3 cluster, show a remarkable agreement and allow us to assert the presence of the, recently highlighted by observations, Vela supercluster at about 180 Mpc/h, right behind the thickest dust layers of our Galaxy.

Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while resolving millions of galaxy sources above ~10^8 Msun , combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the local universe today. Cosmic Dawn ("CoDa") is the first such fully-coupled, radiation-hydrodynamics simulation of reionization of the local universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc)^3 with 4096^3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10^9 Msun , decreasing the abun- dance of faint galaxies around MAB_1600 = [-10,-12]. For most of reionization, star formation was dominated by haloes between 10^10 - 10^11 Msun , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. (Abridged)

We present the measurements of the luminosity-dependent redshift-space three-point correlation functions (3PCFs) for the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample. We compare the 3PCF measurements to the predictions from three different halo and subhalo models. One is the halo occupation distribution (HOD) model and the other two are extensions of the subhalo abundance matching (SHAM) model by allowing the central and satellite galaxies to have different occupation distributions in the host halos and subhalos. Parameters in all the models are chosen to best describe the projected and redshift-space two-point correlation functions (2PCFs) of the same set of galaxies. All three model predictions agree well with the 3PCF measurements for the most luminous galaxy sample, while the HOD model better performs in matching the 3PCFs of fainter samples (with luminosity threshold below $L^*$), which is similar in trend to the case of fitting the 2PCFs. The decomposition of the model 3PCFs into contributions from different types of galaxy triplets shows that on small scales the dependence of the 3PCFs on triangle shape is driven by nonlinear redshift-space distortion (and not by the intrinsic halo shape) while on large scales it reflects the filamentary structure. The decomposition also reveals more detailed differences in the three models, which are related to the radial distribution, the mean occupation function, and the velocity distribution of satellite galaxies inside halos. The results suggest that galaxy 3PCFs can further help constrain the above galaxy-halo relation and test theoretical models.

While the Virgo cluster is the nearest galaxy cluster and therefore the best observed one, little is known about its formation history. In this paper, a set of cosmological simulations that resemble the Local Universe is used to shed the first light on this mystery. The initial conditions for these simulations are constrained with galaxy peculiar velocities of the second catalog of the Cosmicflows project using algorithms developed within the Constrained Local UniversE Simulation project. Boxes of 500 Mpc/h on a side are set to run a series of dark matter only constrained simulations. In each simulation, a unique dark matter halo can be reliably identified as Virgo's counterpart. The properties of these Virgo halos are in agreement at a 10-20% level with the global properties of the observed Virgo cluster. Their zero-velocity masses agree at one-sigma with the observational mass estimate. In all the simulations, the matter falls onto the Virgo objects along a preferential direction that corresponds to the observational filament and the slowest direction of collapse. A study of the mass accretion history of the Virgo candidates reveals the most likely formation history of the Virgo cluster, namely a quiet accretion over the last 7 Gigayears.

We model the luminosity-dependent projected and redshift-space two-point correlation functions (2PCFs) of the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample, using the halo occupation distribution (HOD) model and the subhalo abundance matching (SHAM) model and its extension. All the models are built on the same high-resolution $N$-body simulations. We find that the HOD model generally provides the best performance in reproducing the clustering measurements in both projected and redshift spaces. The SHAM model with the same halo-galaxy relation for central and satellite galaxies (or distinct haloes and subhaloes), when including scatters, has a best-fitting $\chi^2/\rm{dof}$ around $2$--$3$. We therefore extend the SHAM model to the subhalo clustering and abundance matching (SCAM) by allowing the central and satellite galaxies to have different galaxy--halo relations. We infer the corresponding halo/subhalo parameters by jointly fitting the galaxy 2PCFs and abundances and consider subhaloes selected based on three properties, the mass $M_{\rm acc}$ at the time of accretion, the maximum circular velocity $V_{\rm acc}$ at the time of accretion, and the peak maximum circular velocity $V_{\rm peak}$ over the history of the subhaloes. The three subhalo models work well for luminous galaxy samples (with luminosity above $L_*$). For low-luminosity samples, the $V_{\rm acc}$ model stands out in reproducing the data, with the $V_{\rm peak}$ model slightly worse, while the $M_{\rm acc}$ model fails to fit the data. We discuss the implications of the modeling results.

Determining the precise value of the tangential component of the velocity of M31 is a non trivial astrophysical issue, that relies on complicated modeling. This has recently lead to con- flicting estimates, obtained by several groups that used different methodologies and assump- tions. This letter addresses the issue by computing a Bayesian posterior distribution function of this quantity, in order to measure the compatibility of those estimates with LambdaCDM. This is achieved using an ensemble of local group (LG) look-alikes collected from a set of Con- strained Simulations (CSs) of the local Universe, and a standard unconstrained LambdaCDM. The latter allows us to build a control sample of LG-like pairs and to single out the influence of the environment in our results. We find that neither estimate is at odds with LambdaCDM; how- ever, whereas CSs favour higher values of vtan , the reverse is true for estimates based on LG samples gathered from unconstrained simulations, overlooking the environmental element

We use a cosmological simulation of the formation of the Local Group to explore the origin of age and metallicity gradients in dwarf spheroidal galaxies. We find that a number of simulated dwarfs form "outside-in", with an old, metal-poor population that surrounds a younger, more concentrated metal-rich component, reminiscent of dwarf spheroidals like Sculptor or Sextans. We focus on a few examples where stars form in two populations distinct in age in order to elucidate the origin of these gradients. The spatial distributions of the two components reflect their diverse origin; the old stellar component is assembled through mergers, but the young population forms largely in situ. The older component results from a first episode of star formation that begins early but is quickly shut off by the combined effects of stellar feedback and reionization. The younger component forms when a late accretion event adds gas and reignites star formation. The effect of mergers is to disperse the old stellar population, increasing their radius and decreasing their central density relative to the young population. We argue that dwarf-dwarf mergers offer a plausible scenario for the formation of systems with multiple distinct populations and, more generally, for the origin of age and metallicity gradients in dwarf spheroidals.

This paper combines observational datasets and cosmological simulations to generate realistic numerical replicas of the nearby Universe. These latter are excellent laboratories for studies of the non-linear process of structure formation in our neighborhood. With measurements of radial peculiar velocities in the Local Universe (cosmicflows-2) and a newly developed technique, we produce Constrained Local UniversE Simulations (CLUES). To assess the quality of these constrained simulations, we compare them with random simulations as well as with local observations. The cosmic variance, defined as the mean one-sigma scatter of cell-to-cell comparison between two fields, is significantly smaller for the constrained simulations than for the random simulations. Within the inner part of the box where most of the constraints are, the scatter is smaller by a factor 2 to 3 on a 5 Mpc/h scale with respect to that found for random simulations. This one-sigma scatter obtained when comparing the simulated and the observation-reconstructed velocity fields is only 104 +/- 4 km/s i.e. the linear theory threshold. These two results demonstrate that these simulations are in agreement with each other and with the observations of our neighborhood. For the first time, simulations constrained with observational radial peculiar velocities resemble the Local Universe up to a distance of 150 Mpc/h on a scale of a few tens of megaparsecs. When focusing on the inner part of the box, the resemblance with our cosmic neighborhood extends to a few megaparsecs (< 5 Mpc/h). The simulations provide a proper Large Scale environment for studies of the formation of nearby objects.

Recent observational studies have demonstrated that the majority of satellite galaxies tend to orbit their hosts on highly flattened, vast, possibly co-rotating planes. Two nearly parallel planes of satellites have been confirmed around the M31 galaxy and around the Centaurus A galaxy, while the Milky Way also sports a plane of satellites. It has been argued that such an alignment of satellites on vast planes is unexpected in the standard ({\Lambda}CDM) model of cosmology if not even in contradiction to its generic predictions. Guided by {\Lambda}CDM numerical simulations, which suggest that satellites are channeled towards hosts along the axis of the slowest collapse as dictated by the ambient velocity shear tensor, we re-examine the planes of local satellites systems within the framework of the local shear tensor derived from the Cosmicflows-2 dataset. The analysis reveals that the Local Group and Centaurus A reside in a filament stretched by the Virgo cluster and compressed by the expansion of the Local Void. Four out of five thin planes of satellite galaxies are indeed closely aligned with the axis of compression induced by the Local Void. Being the less massive system, the moderate misalignment of the Milky Way's satellite plane can likely be ascribed to its greater susceptibility to tidal torques, as suggested by numerical simulations. The alignment of satellite systems in the local universe with the ambient shear field is thus in general agreement with predictions of the {\Lambda}CDM model.

We study the two main constituent galaxies of a constrained simulation of the Local Group as candidates for the Milky Way (MW) and Andromeda (M31). We focus on the formation of the stellar discs and its relation to the formation of the group as a rich system with two massive galaxies, and investigate the effects of mergers and accretion as drivers of morphological transformations. We use a state-of-the-art hydrodynamical code which includes star formation, feedback and chemical enrichment to carry out our study. We run two simulations, where we include or neglect the effects of radiation pressure from stars, to investigate the impact of this process on the morphologies and star formation rates of the simulated galaxies. We find that the simulated M31 and MW have different formation histories, even though both inhabit, at z=0, the same environment. These differences directly translate into and explain variations in their star formation rates, in-situ fractions and final morphologies. The M31 candidate has an active merger history, as a result of which its stellar disc is unable to survive unaffected until the present time. In contrast, the MW candidate has a smoother history with no major mergers at late times, and forms a disc that grows steadily; at z=0 the simulated MW has an extended, rotationally-supported disc which is dominant over the bulge. Our two feedback implementations predict similar evolution of the galaxies and their discs, although some variations are detected, the most important of which is the formation time of the discs: in the model with weaker/stronger feedback the discs form earlier/later. In summary, by comparing the formation histories of the two galaxies, we conclude that the particular merger/accretion history of a galaxy rather than its environment at the LG-scales is the main driver of the formation and subsequent growth or destruction of galaxy discs.

In this letter we present, for the first time, a study of star formation rate, gas fraction and galaxy morphology of a constrained simulation of the Milky Way (MW) and Andromeda (M31) galaxies, compared to other MW-mass galaxies. By combining with unconstrained simulations we cover a sufficient volume to compare these galaxies environmental densities ranging from the field to that of the Local Group (LG). This is particularly relevant as it has been shown that, quite generally, galaxy properties depend intimately upon their environment, most prominently when galaxies in clusters are compared to those in the field. For galaxies in loose groups such as the LG, however, environmental effects have been less clear. We consider the galaxy's environmental density in spheres of 1200 kpc (comoving) and find that whilst environment does not appear to directly affect morphology, there is a positive trend with star formation rates. This enhancement in star formation occurs systematically for galaxies in higher density environments, regardless whether they are part of the LG or in filaments. Our simulations suggest that the richer environment at Mpc-scales may help replenish the star-forming gas, allowing higher specific star formation rates in galaxies such as the MW.

We examine the angular infall pattern of subhaloes onto host haloes in the context of the large-scale structure. We find that this infall pattern is essentially driven by the shear tensor of the ambient velocity field. Dark matter subhaloes are found to be preferentially accreted along the principal axis of the shear tensor which corresponds to the direction of weakest collapse. We examine the dependence of this preferential infall on subhalo mass, host halo mass and redshift. Although strongest for the most massive hosts and the most massive subhaloes at high redshift, the preferential infall of subhaloes is effectively universal in the sense that its always aligned with the axis of weakest collapse of the velocity shear tensor. It is the same shear tensor that dictates the structure of the cosmic web and hence the shear field emerges as the key factor that governs the local anisotropic pattern of structure formation. Since the small (sub-Mpc) scale is strongly correlated with the mid-range (~10 Mpc) scale - a scale accessible by current surveys of peculiar velocities - it follows that findings presented here open a new window into the relation between the observed large scale structure unveiled by current surveys of peculiar velocities and the preferential infall direction of the Local Group. This may shed light on the unexpected alignments of dwarf galaxies seen in the Local Group

We study the gas distribution in the Milky Way and Andromeda using a constrained cosmological simulation of the Local Group (LG) within the context of the CLUES (Constrained Local UniversE Simulations) project. We analyse the properties of gas in the simulated galaxies at $z=0$ for three different phases: `cold', `hot' and HI, and compare our results with observations. The amount of material in the hot halo ($M_{hot}\approx 4-5\times10^{10}\,$M$_{\odot}$), and the cold ($M_{cold}(r\lesssim10\,$kpc$)\approx10^{8}\,$M$_{\odot}$) and HI ($M_{HI}(r\lesssim50\,$kpc$)\approx 3-4\times10^8\,$M$_{\odot}$) components display a reasonable agreement with observations. We also compute the accretion/ejection rates together with the HI (radial and all-sky) covering fractions. The integrated HI accretion rate within $r=50\,$kpc gives $\sim$$0.2-0.3\,$M$_{\odot}\,$yr$^{-1}$, i.e. close to that obtained from high-velocity clouds in the Milky Way. We find that the global accretion rate is dominated by hot material, although ionized gas with $T\lesssim10^5\,$K can contribute significantly too. The $net$ accretion rates of $all$ material at the virial radii are $6-8\,$M$_{\odot}\,$yr$^{-1}$. At $z=0$, we find a significant gas excess between the two galaxies, as compared to any other direction, resulting from the overlap of their gaseous haloes. In our simulation, the gas excess first occurs at $z\sim1$, as a consequence of the kinematical evolution of the LG.

We explore the impact of cosmic reionization on nearby isolated dwarf galaxies using a compilation of SFHs estimated from deep HST data and a cosmological hydrodynamical simulation of the Local Group. The nearby dwarfs show a wide diversity of star formation histories; from ancient systems that have largely completed their star formation $\sim 10$ Gyr ago to young dwarfs that have formed the majority of their stars in the past $\sim 5$ Gyr to two-component systems characterized by the overlap of comparable numbers of old and young stars. Taken as an ensemble, star formation in nearby dwarfs dips to lower-than-average rates at intermediate times ($4<t$/Gyr $<8$), a feature that we trace in the simulation to the effects of cosmic reionization. Reionization heats the gas and drives it out of the shallow potential wells of low mass halos, affecting especially those below a sharp mass threshold that corresponds to a virial temperature of $\sim 2 \times 10^4 $ $\mathrm{K}$ at $z_{\rm reion}$. The loss of baryons leads to a sharp decline in the star forming activity of early-collapsing systems, which, compounded by feedback from early star formation, empties halos of gas and leaves behind systems where a single old stellar component prevails. In halos below the threshold at $z_{\rm reion}$, reionization heating leads to a delay in the onset of star formation that lasts until the halo grows massive enough to allow some of the remaining gas to cool and form stars. Young stellar components therefore dominate in dwarfs whose halos assemble late and thus form few stars before reionization. Two-component systems may be traced to late mergers of individual examples of the two aforementioned cases. The relative dearth of intermediate-age stars in nearby dwarfs might thus be the clearest signature yet identified of the imprint of cosmic reionization on the star formation history of dwarf galaxies.

We present the results of a numerical study based on the analysis of the MUSIC-2 simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many lines-of-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of un-relaxed systems. For such a sample we measure an average 2D concentration which is ~11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3-6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.

The increasing precision in the determination of the Hubble parameter has reached a per cent level at which large-scale cosmic flows induced by inhomogeneities of the matter distribution become non-negligible. Here we use large-scale cosmological N-body simulations to study statistical properties of the local Hubble parameter as measured by local observers. We show that the distribution of the local Hubble parameter depends not only on the scale of inhomogeneities, but also on how one defines the positions of observers in the cosmic web and what reference frame is used. Observers located in random dark matter haloes measure on average lower expansion rates than those at random positions in space or in the centres of cosmic voids, and this effect is stronger from the halo rest frames compared to the CMB rest frame. We compare the predictions for the local Hubble parameter with observational constraints based on type Ia supernovae (SNIa) and CMB observations. Due to cosmic variance, for observers located in random haloes we show that the Hubble constant determined from nearby SNIa may differ from that measured from the CMB by 0.8 per cent at 1sigma statistical significance. This scatter is too small to significantly alleviate a recently claimed discrepancy between current measurements assuming a flat LCDM model. However, for observers located in the centres of the largest voids permitted by the standard LCDM model, we find that Hubble constant measurements from SNIa would be biased high by 5 per cent, rendering this tension inexistent in this extreme case.

We review how dark matter is distributed in our local neighbourhood from an observational and theoretical perspective. We will start by describing first the dark matter halo of our own galaxy and in the Local Group. Then we proceed to describe the dark matter distribution in the more extended area known as the Local Universe. Depending on the nature of dark matter, numerical simulations predict different abundances of substructures in Local Group galaxies, in the number of void regions and in the abundance of low rotational velocity galaxies in the Local Universe. By comparing these predictions with the most recent observations, strong constrains on the physical properties of the dark matter particles can be derived. We devote particular attention to the results from the Constrained Local UniversE Simulations (CLUES) project, a special set of simulations whose initial conditions are constrained by observational data from the Local Universe. The resulting simulations are designed to reproduce the observed structures in the nearby universe. The CLUES provides a numerical laboratory for simulating the Local Group of galaxies and exploring the physics of galaxy formation in an environment designed to follow the observed Local Universe. It has come of age as the numerical analogue of Near-Field Cosmology.

Peculiar velocities, obtained from direct distance measurements, are data of choice to achieve constrained simulations of the Local Universe reliable down to a scale of a few Megaparsecs. Unlike redshift surveys, peculiar velocities are direct tracers of the underlying gravitational field as they trace both baryonic and dark matter. This paper presents the first attempt to use solely observational peculiar velocities to constrain cosmological simulations of the nearby universe. In order to set up Initial Conditions, a Reverse Zel'dovich Approximation (RZA) is used to displace constraints from their positions at z=0 to their precursors' locations at higher redshifts. An additional new feature replaces original observed radial peculiar velocity vectors by their full 3D reconstructions provided by the Wiener-Filter (WF) estimator. Subsequently, the Constrained Realization of Gaussian fields technique (CR) is applied to build various realizations of the Initial Conditions. The WF/RZA/CR method is first tested on realistic mock catalogs built from a reference simulation similar to the Local Universe. These mocks include errors on peculiar velocities, on data-point positions and a large continuous zone devoid of data in order to mimic galactic extinction. Large scale structures are recovered with a typical accuracy of 5 Megaparsecs/h in position, the best realizations reaching a 2-3 Mpc/h precision, the limit imposed by the RZA linear theory. Then, the method is applied to the first observational radial peculiar velocity catalog of the project Cosmicflows. This paper is a proof of concept that the WF/RZA/CR method can be applied to observational peculiar velocities to successfully build constrained Initial Conditions.

The detection of extremely massive clusters at $z>1$ such as SPT-CL J0546-5345, SPT-CL J2106-5844, and XMMU J2235.3-2557 has been considered by some authors as a challenge to the standard \LCDM$\;$cosmology. In fact,assuming Gaussian initial conditions, the theoretical expectation of detecting such objects is as low as $\leq 1%$. In this \textit{Letter} we discuss the probability of the existence of such objects in the light of the Vector Dark Energy (VDE) paradigm, showing by means of a series of $N$-body simulations that chances of detection are substantially enhanced in this non-standard framework.

We present the results of a series of cosmological $N$-body simulations of a Vector Dark Energy (VDE) model, performed using a suitably modified version of the publicly available \texttt{GADGET}-2 code. The setups of our simulations were calibrated pursuing a twofold aim: 1) to analyze the large scale distribution of massive objects and 2) to determine the properties of halo structure in this different ramework.We observe that structure formation is enhanced in VDE, since the mass function at high redshift is boosted up to a factor of ten with respect to \LCDM, possibly alleviating tensions with the observations of massive clusters at high redshifts and early reionization epoch. Significant differences can also be found for the value of the growth factor, that in VDE shows a completely different behaviour, and in the distribution of voids, which in this cosmology are on average smaller and less abundant. We further studied the structure of dark matter haloes more massive than $5\times10^{13}$\hMsun, finding that no substantial difference emerges when comparing spin parameter, shape, triaxiality and profiles of structures evolved under different cosmological pictures. Nevertheless, minor differences can be found in the concentration-mass relation and the two point correlation function; both showing different amplitudes and steeper slopes.Using an additional series of simulations of a \LCDM\ scenario with the same $\Omega_M$ and $\sigma_8$ used in the VDE cosmology, we have been able to establish whether the modifications induced in the new cosmological picture were due to the particular nature of the dynamical dark energy or a straightforward consequence of the cosmological parameters.

In this work we present cosmological N-body simulations of the Local Universe with initial conditions constrained by the Two-Micron Redshift Survey (2MRS) within a cubic volume of 180 Mpc/h side-length centred at the Local Group. We use a self-consistent Bayesian based approach to explore the joint parameter space of primordial density fluctuations and peculiar velocity fields, which are compatible with the 2MRS galaxy distribution after cosmic evolution. This method (the KIGEN-code) includes the novel ALPT (Augmented Lagrangian Perturbation Theory) structure formation model which combines second order LPT (2LPT) on large scales with the spherical collapse model on small scales. Furthermore we describe coherent flows with 2LPT and include a dispersion term to model fingers-of-god (fogs). These implementations are crucial to avoid artificial filamentary structures, which appear when using a structure formation model with 2LPT and data with compressed fogs. We assume LCDM cosmology throughout our method. The recovered initial Gaussian fields are used to perform a set of 25 constrained simulations. Statistically this ensemble of simulations is in agreement with a reference set of 25 simulations based on randomly seeded Gaussian fluctuations in terms of matter statistics, power-spectra and mass functions. Considering the entire volume, we obtain correlation coefficients of about 98.3% for the comparison between the simulated density fields and the galaxy density field in log-space with Gaussian smoothing scales of r_S=3.5 Mpc/h (74% for r_S=1.4 Mpc/h). The cross power-spectra show correlations with the galaxy distribution which weaken towards smaller length scales until they vanish at scales of 2.2-3.0 Mpc/h. The simulations we present provide a fully nonlinear density and velocity field with a high level of correlation with the observed galaxy distribution at scales of a few Mpc.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

The velocity anisotropy of particles inside dark matter (DM) haloes is an important physical quantity, which is required for the accurate modelling of mass profiles of galaxies and clusters of galaxies. It is typically measured using the ratio of the radial-to-tangential velocity dispersions at a given distance from the halo centre. However, this measure is insufficient to describe the dynamics of realistic haloes, which are typically quite elongated. Studying the velocity distribution in massive DM haloes in cosmological simulations, we find that in the inner parts of the haloes the local velocity ellipsoids are strongly aligned with the major axis of the halo, the alignment being stronger for more relaxed haloes. In the outer regions of the haloes, the alignment becomes gradually weaker and the orientation is more random. These two distinct regions of different degree of the alignment coincide with two characteristic regimes of the DM density profile: shallower and steeper than \rho r^{-2}. This alignment of the local velocity ellipsoids requires reinterpretation of features found in measurements based on the spherically averaged ratio of the radial-to-tangential velocity dispersions. In particular, we show that the velocity distribution in the central halo regions is highly anisotropic. For cluster-size haloes with mass 10^{14}-10^{15} h^-1 Msun, the velocity anisotropy along the major axis is nearly independent of radius and is equal to \beta=1-\sigma^2_{perp}/\sigma^2_{radial}=0.4, which is significantly larger than the previously estimated spherically averaged velocity anisotropy. The alignment of density and velocity anisotropies, and the radial trends may also have some implications for the mass modelling based on kinematical data of such objects as galaxy clusters or dwarf spheroidals, where the orbital anisotropy is a key element in an unbiased mass inference.

The differences between cold (CDM) and warm (WDM) dark matter in the formation of a group of galaxies is examined by running two identical simulations where in the WDM case the initial power spectrum has been altered to mimic a 1keV dark matter particle. The CDM initial conditions were constrained to reproduce at z = 0 the correct local environment within which a "Local Group" (LG) of galaxies may form. Two significant differences between the two simulations are found. While in the CDM case a group of galaxies that resembles the real LG forms, the WDM run fails to reproduce a viable LG, instead forming a diffuse group which is still expanding at z = 0. This is surprising since, due to the suppression of small scale power in its power spectrum, WDM is naively expected to only affect the collapse of small haloes and not necessarily the dynamics on a scale of a group of galaxies. Furthermore the concentration of baryons in halo center's is greater in CDM than in WDM and the properties of the disks differ.

We present results on the clustering of 282,068 galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) sample of massive galaxies with redshifts 0.4<z<0.7 which is part of the Sloan Digital Sky Survey III project. Our results cover a large range of scales from ~0.5 to ~90 Mpc/h. We compare these estimates with the expectations of the flat LCDM cosmological model with parameters compatible with WMAP7 data. We use the MultiDark cosmological simulation together with a simple halo abundance matching technique, to estimate galaxy correlation functions, power spectra, abundance of subhaloes and galaxy biases. We find that the LCDM model gives a reasonable description to the observed correlation functions at z~0.5, which is a remarkably good agreement considering that the model, once matched to the observed abundance of BOSS galaxies, does not have any free parameters. However, we find a deviation (>~10%) in the correlation functions for scales less than ~1 Mpc/h and ~10-40 Mpc/h. A more realistic abundance matching model and better statistics from upcoming observations are needed to clarify the situation. We also estimate that about 12% of the "galaxies" in the abundance-matched sample are satellites inhabiting central haloes with mass M>~1e14 M_sun/h. Using the MultiDark simulation we also study the real space halo bias b(r) of the matched catalogue finding that b=2.00+/-0.07 at large scales, consistent with the one obtained using the measured BOSS projected correlation function. Furthermore, the linear large-scale bias depends on the number density n of the abundance-matched sample as b=-0.048-(0.594+/-0.02)*log(n/(h/Mpc)^3). Extrapolating these results to BAO scales we measure a scale-dependent damping of the acoustic signal produced by non-linear evolution that leads to ~2-4% dips at ~3 sigma level for wavenumbers k>~0.1 h/Mpc in the linear large-scale bias.

Recent observations constrained the tangential velocity of M31 with respect to the Milky Way (MW) to be v_tan<34.4 km/s and the radial velocity to be in the range v_rad=-109+/- 4.4 km/s (van der Marel et al. 2012). In this study we use a large volume high resolution N-body cosmological simulation (Bolshoi) together with three constrained simulations to statistically study this kinematics in the context of the LCDM. The comparison of the ensembles of simulated pairs with the observed LG at the 1-sigma level in the uncertainties has been done with respect to the radial and tangential velocities, the reduced orbital energy (e_tot), angular momentum (l_orb) and the dimensionless spin parameter, lambda. Our main results are: (i) the preferred radial and tangential velocities for pairs in LCDM are v_rad=-80+/-20 km/s, v_tan=50+/-10 km/s, (ii) pairs around that region are 3 to 13 times more common than pairs within the observational values, (iii) 15%to 24% of LG-like pairs in LCDM have energy and angular momentum consistent with observations while (iv) 9% to 13% of pairs in the same sample show similar values in the inferred dimensionless spin parameter. It follows that within current observational uncertainties the quasi-conserved quantities that characterize the orbit of the LG, i.e. e_tot, r_orb and lambda, do not challenge the standard LCDM model, but the model is in tension with regard to the actual values of the radial and tangential velocities. This might hint to a problem of the LCDM model to reproduce the observed LG.

We use dark matter only and full hydrodynamical Constrained Local UniversE Simulations (CLUES) of the formation of the Local Group to study the density profile of subhaloes of the simulated Milky Way and Andromeda galaxies. We show that the Einasto model provides the best description of the subhaloes' density profile, as opposed to the more commonly used NFW profile or any generalisation of it. We further find that the Einasto shape parameter \nEin\ is strongly correlated with the total subhalo mass, pointing towards the notion of a non-universality of the subhaloes' density profile. We observe that the effect of mass loss due to tidal stripping, in both the dark matter only and the hydrodynamical run, is the reduction of the shape parameter \nEin\ between the infall and the present time. Assuming now that the dSphs of our Galaxy follow the Einasto profile and using the maximum and minimum values of \nEin\ from our hydrodynamical simulation as a gauge, we can improve the observational constraints on the \Rmax-\Vmax\ pairs obtained for the brightest satellite galaxies of the Milky Way. When considering only the subhaloes with $-13.2\lesssim M_V\lesssim-8.8$, i.e. the range of luminosity of the classical dwarfs, we find that all our simulated objects are consistent with the observed dSphs if their haloes follow the Einasto model with $1.6\lesssim n_{\rm E} \lesssim5.3$. The numerically motivated Einasto profile for the observed dSphs will alleviate the recently presented "massive failures" problem.

The Reverse Zeldovich Approximation (RZA) is a reconstruction method which allows to estimate the cosmic displacement field from galaxy peculiar velocity data and to constrain initial conditions for cosmological simulations of the Local Universe. In this paper, we investigate the effect of different observational errors on the reconstruction quality of this method. For this, we build a set of mock catalogues from a cosmological simulation, varying different error sources like the galaxy distance measurement error (0 - 20%), the sparseness of the data points, and the maximum catalogue radius (3000 - 6000 km/s). We perform the RZA reconstruction of the initial conditions on these mock catalogues and compare with the actual initial conditions of the simulation. We also investigate the impact of the fact that only the radial part of the peculiar velocity is observationally accessible. We find that the sparseness of a dataset has the highest detrimental effect on RZA reconstruction quality. Observational distance errors also have a significant influence, but it is possible to compensate this relatively well with Wiener Filter reconstruction. We also investigate the effect of different object selection criteria and find that distance catalogues distributed randomly and homogeneously across the sky (such as spiral galaxies selected for the Tully-Fisher method) allow for a higher reconstruction quality than if when data is preferentially drawn from massive objects or dense environments (such as elliptical galaxies). We find that the error of estimating the initial conditions with RZA is always dominated by the inherent non-linearity of data observed at z=0 rather than by the combined effect of the observational errors. Even an extremely sparse dataset with high observational errors still leads to a good reconstruction of the initial conditions on a scale of about 5 Mpc/h.

We propose a new method to recover the cosmological initial conditions of the presently observed galaxy distribution, which can serve to run constrained simulations of the Local Universe. Our method, the Reverse Zeldovich Approximation (RZA), can be applied to radial galaxy peculiar velocity data and extends the previously used Constrained Realizations (CR) method by adding a Lagrangian reconstruction step. The RZA method consists of applying the Zeldovich approximation in reverse to galaxy peculiar velocities to estimate the cosmic displacement field and the initial linear matter distribution from which the present-day Local Universe evolved.We test our method with a mock survey taken from a cosmological simulation. We show that the halo peculiar velocities at z = 0 are close to the linear prediction of the Zeldovich approximation, if a grouping is applied to the data to remove virial motions. We find that the addition of RZA to the CR method significantly improves the reconstruction of the initial conditions. The RZA is able to recover the correct initial positions of the velocity tracers with a median error of only 1.36 Mpc/h in our test simulation. For realistic sparse and noisy data, this median increases to 5 Mpc/h. This is a significant improvement over the previous approach of neglecting the displacement field, which introduces errors on a scale of 10 Mpc/h or even higher. Applying the RZA method to the upcoming high-quality observational peculiar velocity catalogues will generate much more precise constrained simulations of the Local Universe.

In previous works we proposed the Reverse Zeldovich Approximation (RZA) method, which can be used to estimate the cosmological initial conditions underlying the galaxy distribution in the Local Universe using peculiar velocity data. In this paper, we apply the technique to run constrained cosmological simulations from the RZA-reconstructed initial conditions, designed to reproduce the large-scale structure of the Local Universe. We test the method with mock peculiar velocity catalogues extracted from a reference simulation. We first reconstruct the initial conditions of this reference simulation using the mock data, and then run the reconstructed initial conditions forward in time until z=0. We compare the resulting constrained simulations with the original simulation at z=0 to test the accuracy of this method. We also compare them with constrained simulations run from the mock data without the addition of RZA, i.e. using only the currently established constrained realizations (CR) method. Our re-simulations are able to correctly recover the evolution of the large-scale structure underlying the data. The results show that the addition of RZA to the CR method significantly improves both the reconstruction of the initial conditions and the accuracy of the obtained constrained resimulations. Haloes from the original simulation are recovered in the re-simulations with an average accuracy of about 2 Mpc/h on their position and a factor of 2 in mass, down to haloes with a mass of approx 10^14 M_odot/h. In comparison, without RZA the re-simulations recover only the most massive haloes with masses of about 5x10^14 M_odot/h and higher, and with a systematic shift on their position of about approx 10 Mpc/h due to the cosmic displacement field. We show that with the additional Lagrangian reconstruction step introduced by the RZA, this shift can be removed.

In the standard model of cosmology, structure emerges out of non-rotational flow and the angular momentum of collapsing halos is induced by tidal torques. The growth of halo angular momentum in the linear and quasi-linear phases is associated with a shear, curl-free, flow and it is well described within the linear framework of tidal torque theory (TTT). However, TTT is rendered irrelevant as haloes approach turn around and virialization. At that stage the flow field around halos has non-zero vorticity. Using a cosmological simulation, we have examined the importance of the curl of the velocity field (vorticity) in determining halo spin, finding a strong alignment between the two. We have also examined the alignment of vorticity with the principle axes of the shear tensor, finding that it tends to be perpendicular to the axis along which material is collapsing fastest (e1). This behavior is independent of halo masses and cosmic web environment. Our results agree with previous findings on the tendency of halo spin to be perpendicular to e1, and of the spin of (simulated) halos and (observed) galaxies to be aligned with the large-scale structure. Our results imply that angular momentum growth proceeds in two distinct phases. In the first phase angular momentum emerges out of a shear, curl-free, potential flow, as described by TTT. In the second phase, in which haloes approach virialization, the angular momentum emerges out of a vortical flow and halo spin becomes strongly aligned with the vorticity of the ambient flow field.

We introduce the Marenostrum-MultiDark SImulations of galaxy Clusters (MUSIC) Dataset, one of the largest sample of hydrodynamically simulated galaxy clusters with more than 500 clusters and 2000 groups. The objects have been selected from two large N-body simulations and have been resimulated at high resolution using SPH together with relevant physical processes (cooling, UV photoionization, star formation and different feedback processes). We focus on the analysis of the baryon content (gas and star) of clusters in the MUSIC dataset both as a function of aperture radius and redshift. The results from our simulations are compared with the most recent observational estimates of the gas fraction in galaxy clusters at different overdensity radii. When the effects of cooling and stellar feedbacks are included, the MUSIC clusters show a good agreement with the most recent observed gas fractions quoted in the literature. A clear dependence of the gas fractions with the total cluster mass is also evident. The impact of the aperture radius choice, when comparing integrated quantities at different redshifts, is tested: the standard definition of radius at a fixed overdensity with respect to critical density is compared with a definition based on the redshift dependent overdensity with respect to background density. We also present a detailed analysis of the scaling relations of the thermal SZ (Sunyaev Zel'dovich) Effect derived from MUSIC clusters. The integrated SZ brightness, Y, is related to the cluster total mass, M, as well as, the M-Y counterpart, more suitable for observational applications. Both laws are consistent with predictions from the self-similar model, showing a very low scatter. The effects of the gas fraction on the Y-M scaling and the presence of a possible redshift dependence on the Y-M scaling relation are also explored.

We present a cosmography analysis of the Local Universe based on the recently released Two-Micron All-Sky Redshift Survey (2MRS). Our method is based on a Bayesian Networks Machine Learning algorithm (the Kigen-code) which self-consistently samples the initial density fluctuations compatible with the observed galaxy distribution and a structure formation model given by second order Lagrangian perturbation theory (2LPT). From the initial conditions we obtain an ensemble of reconstructed density and peculiar velocity fields which characterize the local cosmic structure with high accuracy unveiling nonlinear structures like filaments and voids in detail. Coherent redshift space distortions are consistently corrected within 2LPT. From the ensemble of cross-correlations between the reconstructions and the galaxy field and the variance of the recovered density fields we find that our method is extremely accurate up to k ~ 1 h Mpc^-1 and still yields reliable results down to scales of about 3-4 h^-1 Mpc. The motion of the local group we obtain within ~ 80 h^-1 Mpc (v_LG=522+-86 km s^-1, l_LG=291^o +- 16^o, b_LG=34^o+-8^o) is in good agreement with measurements derived from the CMB and from direct observations of peculiar motions and is consistent with the predictions of LambdaCDM.

We use a cosmological simulation of the formation of the Local Group of Galaxies to identify a mechanism that enables the removal of baryons from low-mass halos without appealing to feedback or reionization. As the Local Group forms, matter bound to it develops a network of filaments and pancakes. This moving web of gas and dark matter drifts and sweeps a large volume, overtaking many halos in the process. The dark matter content of these halos is unaffected but their gas can be efficiently removed by ram-pressure. The loss of gas is especially pronounced in low-mass halos due to their lower binding energy and has a dramatic effect on the star formation history of affected systems. This "cosmic web stripping" may help to explain the scarcity of dwarf galaxies compared with the numerous low-mass halos expected in \Lambda CDM and the large diversity of star formation histories and morphologies characteristic of faint galaxies. Although our results are based on a single high-resolution simulation, it is likely that the hydrodynamical interaction of dwarf galaxies with the cosmic web is a crucial ingredient so far missing from galaxy formation models.

There are thousands of confirmed detections of star forming galaxies at high redshift (z > 4). These observations rely primarily on the detection of the spectral Lyman Break and the Lyman-alpha emission line. Theoretical modelling of these sources helps to interpret the observations in the framework of the standard cosmological paradigm. We present results from the High-z MareNostrum Project, aimed at constructing a panchromatic picture of the high redshift galaxy evolution that will improve our understanding of young star forming galaxies. Our simulation successfully reproduces the observational constraints from Lyman Break Galaxies and Lyman-alpha emitters at 5 < z < 7 . Based on this model we make predictions on the expected Far Infrared (FIR) emission that should be observed for LAEs. These predictions will help to settle down the question on the dust content of massive high-z galaxies, an issue that will be feasible to probe observationally with the Atacama Large Millimetre Array (ALMA).

We use the recently proposed scale-free mass estimators to determine the masses of the Milky Way (MW) and Andromeda (M31) galaxy in a dark matter only Constrained Local UniversE Simulation (CLUES). While these mass estimators work rather well for isolated spherical host systems, we examine here their applicability to a simulated binary system with a unique satellite population similar to the observed satellites of MW and M31. We confirm that the scale-free estimators work also very well in our simulated Local Group galaxies with the right number of satellites which follow the observed radial distribution. In the isotropic case and under the assumption that the satellites are tracking the total gravitating mass, the power-law index of the radial satellite distribution $N(<r)\propto r^{3-\gamma}$ is directly related to the host's mass profile $M(<r)\propto r^{1-\alpha}$ as $\alpha=\gamma-2$. The use of this relation for any given $\gamma$ leads to highly accurate mass estimations which is a crucial point for observer, since they do not know a priori the mass profile of the MW and M31 haloes. We discuss possible bias in the mass estimators and conclude that the scale-free mass estimators can be satisfactorily applied to the real MW and M31 system.

Radio relics have been discovered in many galaxy clusters. They are believed to trace shock fronts induced by cluster mergers. Cosmological simulations allow us to study merger shocks in detail since the intra-cluster medium is heated by shock dissipation. Using high resolution cosmological simulations, identifying shock fronts and applying a parametric model for the radio emission allows us to simulate the formation of radio relics. We analyze a simulated shock front in detail. We find a rather broad Mach number distribution. The Mach number affects strongly the number density of relativistic electrons in the downstream area, hence, the radio luminosity varies significantly across the shock surface. The abundance of radio relics can be modeled with the help of the radio power probability distribution which aims at predicting radio relic number counts. Since the actual electron acceleration efficiency is not known, predictions for the number counts need to be normalized by the observed number of radio relics. For the characteristics of upcoming low frequency surveys we find that about thousand relics are awaiting discovery.

We present all-sky simulated Fermi maps of gamma-rays from dark matter decay and annihilation in the Local Universe. The dark matter distribution is obtained from a constrained cosmological simulation of the neighboring large-scale structure provided by the CLUES project. The dark matter fields of density and density squared are then taken as an input for the Fermi observation simulation tool to predict the gamma-ray photon counts that Fermi would detect in 5 years of all-sky survey for given dark matter models. Signal-to-noise sky maps have also been obtained by adopting the current Galactic and isotropic diffuse background models released by the Fermi collaboration. We point out the possibility for Fermi to detect a dark matter gamma-ray signal in local extragalactic structures. In particular, we conclude here that Fermi observations of nearby clusters (e.g. Virgo and Coma) and filaments are expected to give stronger constraints on decaying dark matter compared to previous studies. As an example, we find a significant signal-to-noise ratio in dark matter models with a decay rate fitting the positron excess as measured by PAMELA. This is the first time that dark matter filaments are shown to be promising targets for indirect detection of dark matter. On the other hand, the prospects for detectability of annihilating dark matter in local extragalactic structures are less optimistic even with extreme cross-sections. We make the dark matter density and density squared maps available online at http://www.clues-project.org/articles/darkmattermaps.html

The size of the horizon at the matter-radiation equality is a key scale of the Big Bang cosmology that is directly related to the energy-matter content of the Universe. In this letter, we argue that this scale can be accurately measured from the observed clustering of galaxies in new large scale surveys. We demonstrate that the zero-crossing, r_c, of the 2-point galaxy correlation function is closely related to the horizon size at matter-radiation equality for a large variety of flat LCDM models. Using large-volume cosmological simulations, we also show that the pristine zero-crossing is unaltered by non-linear evolution of density fluctuations, redshift distortions and galaxy biases. This makes r_c a very powerful standard ruler that can be accurately measured, at a percent level, in upcoming experiments that will collect redshifts of millions of galaxies and quasars.

Upcoming radio telescopes will allow to study the radio sky at low frequencies with unprecedent sensitivity and resolution. New surveys are expected to discover a large number of new radio sources. Here we investigate the abundance of radio relics, i.e. steep-spectrum diffuse radio emission coming from the periphery of galaxy clusters, which are believed to trace shock waves induced by cluster mergers. With the advent of comprehensive relic samples a framework is needed to analyze statistically the relic abundance. To this end, we introduce the probability to find a relic located in a galaxy cluster with given mass and redshift allowing us to relate the halo mass function of the Universe with the radio relic number counts. Up to date about 45 relics have been reported and we compile the resulting counts, N(>S_1.4). In principle, the parameters of the distribution could be determined using a sufficiently large relic sample. However, since the number of known relics is still small for that purpose we use the MareNostrum Universe simulation to determine the relic radio power scaling with cluster mass and redshift. Our model is able to reproduce the recently found tentative evidence for an increase in the fraction of clusters hosting relics, both with X-ray luminosity and redshift, using an X-ray flux limited cluster sample. Moreover, we find that a considerable fraction of faint relics (S_1.4 < ~10 mJy) reside in clusters with an X-ray flux below ~3e-12 erg/s/cm^2. Finally, we estimate the number of radio relics which await discovery by future low frequency surveys proposed for LOFAR and WSRT. We estimate that the WODAN survey proposed for WSRT may discover 900 relics and that the LOFAR-Tier 1-120 MHz survey may discover about 2500 relics. However, the actual number of newly discovered relics will crucially depend on the existence of sufficiently complete galaxy cluster catalogues.

We present the online MultiDark Database -- a Virtual Observatory-oriented, relational database for hosting various cosmological simulations. The data is accessible via an SQL (Structured Query Language) query interface, which also allows users to directly pose scientific questions, as shown in a number of examples in this paper. Further examples for the usage of the database are given in its extensive online documentation (www.multidark.org). The database is based on the same technology as the Millennium Database, a fact that will greatly facilitate the usage of both suites of cosmological simulations. The first release of the MultiDark Database hosts two 8.6 billion particle cosmological N-body simulations: the Bolshoi (250/h Mpc simulation box, 1/h kpc resolution) and MultiDark Run1 simulation (MDR1, or BigBolshoi, 1000/h Mpc simulation box, 7/h kpc resolution). The extraction methods for halos/subhalos from the raw simulation data, and how this data is structured in the database are explained in this paper. With the first data release, users get full access to halo/subhalo catalogs, various profiles of the halos at redshifts z=0-15, and raw dark matter data for one time-step of the Bolshoi and four time-steps of the MultiDark simulation. Later releases will also include galaxy mock catalogs and additional merging trees for both simulations as well as new large volume simulations with high resolution. This project is further proof of the viability to store and present complex data using relational database technology. We encourage other simulators to publish their results in a similar manner.

We present theoretical results for the expected fraction of Lyman Break Galaxies (LBGs) to be detected as strong Lyman alpha emitters (LAEs) in the redshift range 5 < z < 7. We base our analysis on the 2-billion particle SPH simulation MareNostrum High-z Universe. We approximate galaxies as static dusty slabs with an additional clumpy dust distribution affecting stellar populations younger than 25 Myr. The model for the Lyman alpha escape fraction is based on the results of our Monte-Carlo radiative transfer code (CLARA) for a slab configuration. We also fix the transmission of Lyman alpha photons through the intergalactic medium to a constant value of 50% at all redshifts. From the results of this model we calculate xLya, the fraction of Lyman Break Galaxies with Lyman alpha equivalent width (EW) larger than 50 Angstrom. We find a remarkable agreement with observational data at 4.5 < z < 6. For bright (-22 < MUV < -20.5) and faint (-20.5 < MUV < -18.5) galaxies our model predicts xLya = 0.02 \pm 0.01 and xLya = 0.47 \pm 0.01 while observers report xLya = 0.08 \pm 0.02 and xLya = 0.47 \pm 0.16, respectively. Additional evolution of the extinction model at redshift z \sim 7, that decreases the intensity of transmitted Lyman alpha radiation by a factor of f_T = 0.4 as to match the LAE luminosity function at z \sim 6.5, naturally provides a good match for the recently reported xLya fractions at z > 6.3. Exploring different toy models for the Lyman alpha escape fraction, we show that a decreasing Lyman alpha escape fraction with increasing UV galaxy luminosity is a key element in our model to explain the of larger xLya fractions for fainter LBGs.

Using Constrained Local UniversE Simulations (CLUES) of the formation of the Local Group in a cosmological context we investigate the recently highlighted problem that the majority of the most massive dark subhaloes of the Milky Way are too dense to host any of its bright satellites. In particular, we examine the influence of baryonic processes and find that they leave a twofold effect on the relation between the peak of the rotation curve and its position (Vmax and Rmax). Satellites with a large baryon fraction experience adiabatic contraction thus decreasing Rmax while leaving Vmax more or less unchanged. Subhaloes with smaller baryon fractions undergo a decrease in Vmax possibly due to outflows of material. Furthermore, the situation of finding subhaloes in simulations that lie outside the confidence interval for possible hosts of the bright MW dwarf spheroidals, appears to be far more prominent in cosmologies with a high sigma8 normalisation and depends on the mass of the host. We conclude that the problem cannot be simply solved by including baryonic processes and hence demands further investigations.

Using a dark matter only Constrained Local UniversE Simulation (CLUES) we examine the existence of subhaloes that change their affiliation from one of the two prominent hosts in the Local Group (i.e. the Milky Way and the Andromeda galaxy) to the other, and call these objects "renegade subhaloes". In light of recent claims that the two Magellanic Clouds (MCs) may have originated from another region (or even the outskirts) of the Local Group or that they have been spawned by a major merger in the past of the Andromeda galaxy, we investigate the nature of such events. However, we cannot confirm that renegade subhaloes enter as deep into the potential well of their present host nor that they share the most simplest properties with the MCs, namely mass and relative velocity. Our simulation rather suggests that these renegade subhaloes appear to be flying past one host before being pulled into the other. A merger is not required to trigger such an event, it is rather the distinct environment of our simulated Local Group facilitating such behavior. Since just a small fraction of the full z=0 subhalo population are renegades, our study indicates that it will be intrinsically difficult to distinguish them despite clear differences in their velocity, radial distribution, shape and spin parameter distributions.

Using CLARA (Code for Lyman Alpha Radiation Analysis) we constrain the escape fraction of Lyman-Alpha radiation in galaxies in the redshift range 5<z<7, based on the MareNostrum High-z Universe, a SPH cosmological simulation with more than 2 billion particles. We approximate Lyman-Alpha Emitters (LAEs) as dusty gaseous slabs with Lyman-Alpha radiation sources homogeneously mixed in the gas. Escape fractions for such a configuration and for different gas and dust contents are calculated using our newly developed radiative transfer code CLARA. The results are applied to the MareNostrum High-z Universe numerical galaxies. The model shows a weak redshift evolution and good agreement with estimations of the escape fraction as a function of reddening from observations at z \sim 2.2 and z \sim 3. We extend the slab model by including additional dust in a clumpy component in order to reproduce the UV con- tinuum luminosity function and UV colours at redshifts z>~5. The LAE Luminosity Function (LF) based on the extended clumpy model reproduces broadly the bright end of the LF derived from observations at z \sim 5 and z \sim 6. At z \sim 7 our model over-predicts the LF by roughly a factor of four, presumably because the effects of the neutral intergalactic medium are not taken into account. The remaining tension between the observed and simulated faint end of the LF, both in the UV-continuum and Lyman-Alpha at redshifts z \sim 5 and z \sim 6 points towards an overabundance of simulated LAEs hosted in haloes of masses 1.0x10^10h-1Msol < Mh < 4.0x10^10h-1Msol. Given the difficulties in explaining the observed overabundance by dust absorption, a probable origin of the mismatch are the high star formation rates in the simulated haloes around the quoted mass range. A more efficient supernova feedback should be able to regulate the star formation process in the shallow potential wells of these haloes.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

A long-standing problem of strong lensing by galaxy clusters regards the observed high rate of giant gravitational arcs as compared to the predictions in the framework of the "standard" cosmological model. Recently, few other inconsistencies between theoretical expectations and observations have been claimed which regard the large size of the Einstein rings and the high concentrations of few clusters with strong lensing features. All of these problems consistently indicate that observed galaxy clusters may be gravitational lenses stronger than expected. We use clusters extracted from the MareNostrum Universe to build up mock catalogs of galaxy clusters selected through their X-ray flux. We use these objects to estimate the probability distributions of lensing cross sections, Einstein rings, and concentrations for the sample of 12 MACS clusters at $z>0.5$ presented in Ebeling et al. (2007) and discussed in Zitrin et al. (2010). We find that simulated clusters produce $\sim 50%$ less arcs than observed clusters do. The medians of the distributions of the Einstein ring sizes differ by $\sim 25%$ between simulations and observations. We estimate that, due to cluster triaxiality and orientation biases affecting the lenses with the largest cross sections, the concentrations of the individual MACS clusters inferred from the lensing analysis should be up to a factor of $\sim 2$ larger than expected from the $\Lambda$CDM model. The arc statistics, the Einstein ring, and the concentration problems in strong lensing clusters are mitigated but not solved on the basis of our analysis. Nevertheless, due to the lack of redshifts for most of the multiple image systems used for modeling the MACS clusters, the results of this work will need to be verified with additional data. The upcoming CLASH program will provide an ideal sample for extending our comparison (abridged).

We identify shocked gas in simulated galaxy clusters extracted from the MareNostrum Universe simulation (Gottloeber et al. 2006) assuming that shock waves are regions of electron acceleration. We perform flux number counts within the framework of the non-thermal emission model developed by Hoeft et al. (2008). Results are presented at two different observing frequencies, i.e. 1.4 GHz and 120 MHz, posing interesting constraints for LOFAR and upcoming radio telescopes.

We study the differences and similarities in the luminosities of bound, infalling and the so-called backsplash (Gill et al. 2005) galaxies of the Milky Way and M31 using a hydrodynamical simulation performed within the Constrained Local UniversE Simulation (CLUES) project. The simulation models the formation of the Local Group within a self-consistent cosmological framework. We find that even though backsplash galaxies passed through the virial radius of their host halo and hence may have lost a (significant) fraction of their mass, their stellar populations are hardly affected. This leaves us with comparable luminosity functions for infalling and backsplash galaxies and hence little hope to decipher their past (and different) formation and evolutionary histories by luminosity measurements alone. Nevertheless, due to the tidal stripping of dark matter we find that the mass-to-light ratios have changed when comparing the various populations against each other: they are highest for the infalling galaxies and lowest for the bound satellites with the backsplash galaxies in-between.

Using a high resolution DM simulation of the Local Group (LG), conducted within the framework of the Constrained Local UniversE Simulation (CLUES) project, we investigate the nature of how satellites of the MW and M31 are accreted. Satellites of these galaxies are accreted anisotropically, entering the virial radius of their hosts from specific "spots" with respect to the large scale structure. Furthermore, material which is tidally stripped from these satellites is also, at z=0, distributed anisotropically and is characterized by an ellipsoidal sub-volume embedded in the halo. The angular pattern created by the locus of satellite infall points and the projected stripped DM is investigated within a coordinate system determined by the location of the LG companion and the simulated Virgo cluster across concentric shells ranging from 0.1 to 5 r_vir. Remarkably, the principal axis of the ellipsoidal sub-volume shows a coherent alignment extending from well within the halo to a few r_vir. A spherical harmonics transform applied to the angular distributions confirms the visual impression: namely, the angular distributions of both the satellites entry points and stripped DM is dominated by the l=2 quadrupole term, whose major principal axis is aligned across the shells considered. It follows that the structure of the outer (r>0.5 r_vir) regions of the main halos is closely related to the cosmic web within which it is embedded. Given the hypothesis that a large fraction of the stellar halo of the Milky Way has been accreted from satellites, the present results can be applied to the stellar halo of the MW and M31. We predict that the remnants of tidally stripped satellites should be embedded in streams of material composed of dark matter and stars. The present results can therefore shed light on the existence of satellites embedded within larger streams of matter, such as the Segue 2 satellite.

We estimate the expected distribution of displacements between the dark matter and gas cores in simulated clusters. We use the MareNostrum Universe, one of the largest non radiative, SPH Lambda CDM cosmological simulations. We find that projected 2-D displacements between dark matter and gas, equal or larger than the observed in the Bullet Cluster, are expected in 1% to 2% of the clusters with masses larger than 10^{14} Msun. The 2-D displacement distribution is roughly the same between redshifts 0<z<0.5 when multiplied by a factor of (1+z)^{-1/2}. We conclude that the separations between dark matter and gas as observed in the bullet cluster can be easily found in a Lambda CDM universe. Furthermore we find that the displacement distribution is not very sensitive to the normalization of the power spectrum. Upcoming surveys could extend the measurements of these displacements between dark matter and gas into large samples of hundreds of clusters, providing a potential test for Lambda CDM.

The strong lensing events that are observed in compact clusters of galaxies can, both statistically and individually, return important clues about the structural properties of the most massive structures in the Universe. Substantial work is ongoing in order to understand the degree of similarity between the lensing cluster population and the population of clusters as a whole, with members of the former being likely more massive, compact, and substructured than members of the latter. In this work we exploit synthetic clusters extracted from the {\sc MareNostrum Universe} cosmological simulation in order to estimate the correlation between the strong lensing efficiency and other bulk properties of lensing clusters, such as the virial mass and the bolometric X-ray luminosity. We found that a positive correlation exist between all these quantities, with the substantial scatter being smaller for the luminosity-cross section relation. We additionally used the relation between the lensing efficiency and the virial mass in order to construct a synthetic optical depth that agrees well with the true one, while being extremely faster to be evaluated. We finally estimated what fraction of the total giant arc abundance is recovered when galaxy clusters are selected according to their dynamical activity or their X-ray luminosity. Our results show that there is a high probability for high-redshift strong lensing clusters to be substantially far away from dynamical equilibrium, and that $30-40\%$ of the total amount of giant arcs are lost if looking only at very X-ray luminous objects.

We present the first detailed structure formation and radiative transfer simulations of the reionization history of our cosmic neighbourhood. To this end, we follow the formation of the Local Group of galaxies and nearby clusters by means of constrained simulations, which use the available observational constraints to construct a representation of those structures which reproduces their actual positions and properties at the present time. We find that the reionization history of the Local Group is strongly dependent on the assumed photon production efficiencies of the ionizing sources, which are still poorly constrained. If sources are relatively efficient, i.e. the process is 'photon-rich', the Local Group is primarily ionized externally by the nearby clusters. Alternatively, if the sources are inefficient, i.e. reionization is 'photon-poor' the Local Group evolves largely isolated and reionizes itself. The mode of reionization, external vs. internal, has important implications for the evolution of our neighbourhood, in terms of e.g. its satellite galaxy populations and primordial stellar populations. This therefore provides an important avenue for understanding the young universe by detailed studies of our nearby structures.

The local universe is the best known part of our universe. Within the CLUES project (http://clues-project.org - Constrained Local UniversE Simulations) we perform numerical simulations of the evolution of the local universe. For these simulations we construct initial conditions based on observational data of the galaxy distribution in the local universe. Here we review the technique of these constrained simulations. In the second part we summarize our predictions of a possible Warm Dark Matter cosmology for the observed local distribution of galaxies and the local spectrum of mini-voids as well as a study of the satellite dynamics in a simulated Local Group.

Strong lensing is one of the most direct probes of the mass distribution in the inner regions of galaxy clusters. It can be used to constrain the density profiles and to measure the mass of the lenses. Moreover, the abundance of strong lensing events can be used to constrain the structure formation and the cosmological parameters through the so-called "arc-statistics" approach. However, several issues related to the usage of strong lensing clusters in cosmological applications are still controversial, leading to the suspect that several biases may affect this very peculiar class of objects. With this study we aim at better understanding the properties of galaxy clusters which can potentially act as strong lenses. We do so by investigating the properties of a large sample of galaxy clusters extracted from the N-body/hydrodynamical simulation MareNostrum Universe. We explore the correlation between the cross section for lensing and many properties of clusters, like the mass, the three-dimensional and projected shapes, their concentrations, the X-ray luminosity and the dynamical activity. We find that the probability of strong alignments between the major axes of the lenses and the line of sight is a growing function of the lensing cross section. In projection, the strong lenses appear rounder within R200, but we find that their cores tend to be more elliptical as the lensing cross section increases. We also find that the cluster concentrations estimated from the projected density profiles tend to be biased high. The X-ray luminosity of strong lensing clusters is higher than that of normal lenses of similar mass and redshift. This is particular significant for the least massive lenses. Finally, we find that the strongest lenses generally exhibit an excess of kinetic energy within the virial radius, indicating that they are more dynamically active than usual clusters.

We use two simulations performed within the Constrained Local UniversE Simulation (CLUES) project to study both the shape and radial alignment of (the dark matter component of) subhaloes; one of the simulations is a dark matter only model while the other run includes all the relevant gas physics and star formation recipes. We find that the involvement of gas physics does not have a statistically significant effect on either property -- at least not for the most massive subhaloes considered in this study. However, we observe in both simulations including and excluding gasdynamics a (pronounced) evolution of the dark matter shapes of subhaloes as well as of the radial alignment signal since infall time. Further, this evolution is different when positioned in the central and outer regions of the host halo today; while subhaloes tend to become more aspherical in the central 50% of their host's virial radius, the radial alignment weakens in the central regime while strengthening in the outer parts. We confirm that this is due to tidal torquing and the fact that subhaloes at pericentre move too fast for the alignment signal to respond.

We use a simulation performed within the Constrained Local UniversE Simulation (CLUES) project to study a realistic Local Group-like object. We employ this group as a numerical laboratory for studying the evolution of the population of its subhaloes from the point of view of the effects it may have on the origin of different types of dwarf galaxies. We focus on the processes of tidal stripping of the satellites, their interaction, merging and grouping before infall. The tidal stripping manifests itself in the transition between the phase of mass accretion and mass loss seen in most subhaloes, which occurs at the moment of infall on to the host halo, and the change of the shape of their mass function with redshift. Although the satellites often form groups, they are loosely bound within them and do not interact with each other. The infall of a large group could however explain the observed peculiar distribution of the Local Group satellites, but only if it occurred recently. Mergers between prospective subhaloes are significant only during an early stage of evolution, i.e. more than 7 Gyr ago, when they are still outside the host haloes. Such events could thus contribute to the formation of more distant early type Milky Way companions. Once the subhaloes enter the host halo the mergers become very rare.

Galaxy surveys have shown that luminous galaxies are mainly distributed in large filaments and galaxy clusters. The remaining large volumes are virtually devoid of luminous galaxies. This is in concordance with the formation of the large-scale structure in Universe as derived from cosmological simulations. However, the numerical results indicate that cosmological voids are abundantly populated with dark matter haloes which may in principle host dwarf galaxies. Observational efforts have in contrast revealed, that voids are apparently devoid of dwarf galaxies. We investigate the formation of dwarf galaxies in voids by hydrodynamical cosmological simulations. Due to the cosmic ultra-violet background radiation low-mass haloes show generally are reduced baryon fraction. We determine the characteristic mass below which dwarf galaxies are baryon deficient. We show that the circular velocity below which the accretion of baryons is suppressed is approximately 40 km/s. The suppressed baryon accretion is caused by the photo-heating due to the UV-background. We set up a spherical halo model and show that the effective equation of state of the gas in the periphery of dwarf galaxies determines the characteristic mass. This implies that any process which heats the gas around dwarf galaxies increases the characteristic mass and thus reduces the number of observable dwarf galaxies.

We discuss the rest-frame UV emission between 5< z < 7 from the MareNostrum High-z Universe, a SPH simulation done with more than 2 billion particles. Cosmological simulations of galaxy formation generally overpredict the UV restframe luminosity function at high redshift, both at the bright and faint ends. In this Letter we explore a dust attenuation model where a larger extinction is applied to star populations younger than a given age, mimicking the effect of a clumpy interstellar medium. We show that this scenario fits reasonably well both the UV luminosity functions and the UV-continuum slopes derived from observations. The model assumes a large obscuration for stars younger than 25 Myr from the gas clouds where they should be embedded at their formation time. We find that the optical depth in these clouds should be between 30 and 100 times larger than the mean optical depth for the homogeneous part of the interstellar medium. These values are one order of magnitude larger than those estimated in local galaxies. Therefore, we conclude that LambdaCDM predictions for the high-z UV emission can accommodate the current observations if we consider a dust extinction model based on the assumption of a clumpy environment at high redshift.

We examine the properties of satellites found in high resolution simulations of the local group. We use constrained simulations designed to reproduce the main dynamical features that characterize the local neighborhood, i.e. within tens of Mpc around the Local Group (LG). Specifically, a LG-like object is found located within the 'correct' dynamical environment and consisting of three main objects which are associated with the Milky Way, M31 and M33. By running two simulations of this LG from identical initial conditions - one with and one without baryons modeled hydrodynamically - we can quantify the effect of gas physics on the $z=0$ population of subhaloes in an environment similar to our own. We find that above a certain mass cut, $M_{\rm sub} > 2\times10^{8}h^{-1} M_{\odot}$ subhaloes in hydrodynamic simulations are more radially concentrated than those in simulations with out gas. This is caused by the collapse of baryons into stars that typically sit in the central regions of subhaloes, making them denser. The increased central density of such a subhalo, results in less mass loss due to tidal stripping than the same subhalo simulated with only dark matter. The increased mass in hydrodynamic subhaloes with respect to dark matter ones, causes dynamical friction to be more effective, dragging the subhalo towards the centre of the host. This results in these subhaloes being effectively more radially concentrated then their dark matter counterparts.

Using high-resolution simulations within the Cold and Warm Dark Matter models we study the evolution of small scale structure in the Local Volume, a sphere of 8 Mpc radius around the Local Group. We compare the observed spectrum of mini-voids in the Local Volume with the spectrum of mini-voids determined from the simulations. We show that the \LWDM model can easily explain both the observed spectrum of mini-voids and the presence of low-mass galaxies observed in the Local Volume, provided that all haloes with circular velocities greater than 20 km/s host galaxies. On the contrary within the LCDM model the distribution of the simulated mini-voids reflects the observed one if haloes with maximal circular velocities larger than 35 km/s host galaxies. This assumption is in contradiction with observations of galaxies with circular velocities as low as 20 km/s in our Local Universe. A potential problem of the LWDM model could be the late formation of the haloes in which the gas can be efficiently photo-evaporated. Thus star formation is suppressed and low-mass haloes might not host any galaxy at all.

We present a new method of constraining the mass and velocity anisotropy profiles of galaxy clusters from kinematic data. The method is based on a model of the phase space density which allows the anisotropy to vary with radius between two asymptotic values. The characteristic scale of transition between these asymptotes is fixed and tuned to a typical anisotropy profile resulting from cosmological simulations. The model is parametrized by two values of anisotropy, at the centre of the cluster and at infinity, and two parameters of the NFW density profile, the scale radius and the scale mass. In order to test the performance of the method in reconstructing the true cluster parameters we analyze mock kinematic data for 20 relaxed galaxy clusters generated from a cosmological simulation of the standard LCDM model. We use Bayesian methods of inference and the analysis is carried out following the Markov Chain Monte Carlo approach. The parameters of the mass profile are reproduced quite well, but we note that the mass is typically underestimated by 15 percent, probably due to the presence of small velocity substructures. The constraints on the anisotropy profile for a single cluster are in general barely conclusive. Although the central asymptotic value is determined accurately, the outer one is subject to significant systematic errors caused by substructures at large clustercentric distance. The anisotropy profile is much better constrained if one performs joint analysis of at least a few clusters. In this case it is possible to reproduce the radial variation of the anisotropy over two decades in radius inside the virial sphere.

The paper has been withdrawn because double checking and comparison with other data sets after the original submission showed that a broken R-band filter at the Maidanak telescope had affected our quasar monitoring observations in the years 2004 and 2005. They had led to partially spurious measurements, hence our original analysis and conclusions are not reliable.

A dynamical classification of the cosmic web is proposed. The large scale environment is classified into four web types: voids, sheets, filaments and knots. The classification is based on the evaluation of the deformation tensor, i.e. the Hessian of the gravitational potential, on a grid. The classification is based on counting the number of eigenvalues above a certain threshold, lambda_th at each grid point, where the case of zero, one, two or three such eigenvalues corresponds to void, sheet, filament or a knot grid point. The collection of neighboring grid points, friends-of-friends, of the same web attribute constitutes voids, sheets, filaments and knots as web objects. A simple dynamical consideration suggests that lambda_th should be approximately unity, upon an appropriate scaling of the deformation tensor. The algorithm has been applied and tested against a suite of (dark matter only) cosmological N-body simulations. In particular, the dependence of the volume and mass filling fractions on lambda_th and on the resolution has been calculated for the four web types. Also, the percolation properties of voids and filaments have been studied. Our main findings are: (a) Already at lambda_th = 0.1 the resulting web classification reproduces the visual impression of the cosmic web. (b) Between 0.2 < lambda_th < 0.4, a system of percolated voids coexists with a net of interconected filaments. This suggests a reasonable choice for lambda_th as the parameter that defines the cosmic web. (c) The dynamical nature of the suggested classification provides a robust framework for incorporating environmental information into galaxy formation models, and in particular the semi-analytical ones.

Large-scale diffuse radio emission is observed in some clusters of galaxies. There is ample of evidence that the emission has its origin in synchrotron losses of relativistic electrons, accelerated in the course of clusters mergers. In a cosmological simulation we locate the structure formation shocks and estimate their radio emission. We proceed as follows: Introducing a novel approach to identify strong shock fronts in an SPH simulation, we determine the Mach number as well as the downstream density and temperature in the MareNostrum Universe simulation which has 2x1024^3 particles in a 500 Mpc/h box and was carried out with non-radiative physics. Then, we estimate the radio emission using the formalism derived in Hoeft & Brueggen (2007) and produce artificial radio maps of massive clusters. Several of our clusters show radio objects with similar morphology to large-scale radio relics found in the sky, whereas about half of the clusters show only very little radio emission. In agreement with observational findings, the maximum diffuse radio emission of our clusters depends strongly on their X-ray temperature. We find that the so-called accretion shocks cause only very little radio emission. We conclude that a moderate efficiency of shock acceleration, namely xi_e <= 0.005, and moderate magnetic fields in the region of the relics, namely 0.07 to 0.8 muGauss are sufficient to reproduce the number density and luminosity of radio relics.

We study the distribution function (DF) of dark matter particles in haloes of mass range 10^{14}--10^{15}\Msun. In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard \LambdaCDM model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L, a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form f_{E}(E)[1+L^{2}/(2L_{0}^{2})]^{-\beta_{\infty}+\beta_{0}}L^{-2\beta_{0}}. This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (\beta_{0} and \beta_{\infty}) and the scale of transition between them (L_{0}). The energy part f_{E}(E) is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase space density Q=\rho(r)/\sigma^{3}(r).

This technical paper describes a software package that was designed to produce initial conditions for large cosmological simulations in the context of the Horizon collaboration. These tools generalize E. Bertschinger's Grafic1 software to distributed parallel architectures and offer a flexible alternative to the Grafic2 software for ``zoom'' initial conditions, at the price of large cumulated cpu and memory usage. The codes have been validated up to resolutions of 4096^3 and were used to generate the initial conditions of large hydrodynamical and dark matter simulations. They also provide means to generate constrained realisations for the purpose of generating initial conditions compatible with, e.g. the local group, or the SDSS catalog.

Despite a history that dates back at least a quarter of a century studies of voids in the large--scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void--finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recent increased interest in voids, both in very large galaxy surveys and in detailed simulations of cosmic structure formation, this situation is very unfortunate. We here present the first systematic comparison study of thirteen different void finders constructed using particles, haloes, and semi--analytical model galaxies extracted from a subvolume of the Millennium simulation. The study includes many groups that have studied voids over the past decade. We show their results and discuss their differences and agreements. As it turns out, the basic results of the various methods agree very well with each other in that they all locate a major void near the centre of our volume. Voids have very underdense centres, reaching below 10 percent of the mean cosmic density. In addition, those void finders that allow for void galaxies show that those galaxies follow similar trends. For example, the overdensity of void galaxies brighter than $m_B = -20 $ is found to be smaller than about -0.8 by all our void finding algorithms.

We describe the basic ideas of MPI parallelization of the N-body Adaptive Refinement Tree (ART) code. The code uses self-adaptive domain decomposition where boundaries of the domains (parallelepipeds) constantly move -- with many degrees of freedom -- in the search of the minimum of CPU time. The actual CPU time spent by each MPI task on previous time-step is used to adjust boundaries for the next time-step. For a typical decomposition of 5^3 domains, the number of possible changes in boundaries is 3^{84}. We describe two algorithms of finding minimum of CPU time for configurations with a large number of domains. Each MPI task in our code solves the N-body problem where the large-scale distribution of matter outside of the boundaries of a domain is represented by relatively few temporary large particles created by other domains. At the beginning of a zero-level time-step, domains create and exchange large particles. Then each domain advances all its particles for many small time-steps. At the end of the large step, the domains decide where to place new boundaries and re-distribute particles. The scheme requires little communications between processors and is very efficient for large cosmological simulations.

Thanks to the ever increasing computational power and the development of more sophisticated algorithms, numerical N-body simulations are now uncovering several phenomenological relations between the physical properties of dark matter haloes in position and velocity space. It is the aim of the present work to investigate in detail the dynamical structure of dark matter haloes, as well as its possible dependence on mass and its evolution with redshift up to z=5. We use high-resolution cosmological simulations of individual objects to compute the radially-averaged profiles of several quantities, scaled by the radius Rmax at which the circular velocity attains its maximum value, Vmax. No systematic dependence on mass or cosmic epoch are found within Rmax, and all the different radial profiles are well fit by simple analytical models. However, our results suggest that several properties are not `universal' outside this radius. [Abridged]

We have used GADGET2 to simulate the formation of an elliptical galaxy in a cosmological dark matter halo with mass 3x10^12M_Sun/h. Using a stellar population synthesis model has allowed us to compute magnitudes, colours and surface brightness profiles. We have included a model to follow the growth of a central black hole and we have compared the results of simulations with and without feedback from AGNs. We have studied the interplay between cold gas accretion and merging in the development of galactic morphologies, the link between colour and morphology evolution, the effect of AGN feedback on the photometry of early type galaxies, the redshift evolution in the properties of quasar hosts, and the impact of AGN winds on the chemical enrichment of the intergalactic medium (IGM). We have found that the early phases of galaxy formation are driven by the accretion of cold filamentary flows, which form a disc at the centre of the dark matter halo. When the dark matter halo is sufficiently massive to support the propagation of a stable shock, cold accretion is shut down, and the star formation rate begins to decline. Mergers transform the disc into an elliptical galaxy, but also bring gas into the galaxy. Without a mechanism that removes gas from the merger remnants, the galaxy ends up with blue colours, atypical for its elliptical morphology. AGN feedback can solve this problem even with a fairly low heating efficiency. We have also demonstrated that AGN winds are potentially important for the metal enrichment of the IGM a high redshift.(abridged)

We explore the dependence of the radial alignment of subhalos on the mass of the host halo they orbit in. As the effect is seen on a broad range of scales including massive clusters as well as galactic systems it only appears natural to explore this phenomenon by means of cosmological simulations covering the same range in masses. We have 25 well resolved host dark matter halos at our disposal ranging from 10^15 Msun/h down to 10^12 Msun/h each consisting of order of a couple of million particles within the virial radius. We observe that subhalos tend to be more spherical than isolated objects. Both the distributions of sphericity and triaxiality of subhalos are Gaussian distributed with peak values of s approx. 0.80 and T approx. 0.56, irrespective of host mass. Interestingly we note that the radial alignment is independent of host halo mass and the distribution of \cos\theta (i.e. the angle between the major axis E_a of each subhalo and the radius vector of the subhalo in the reference frame of the host) is well fitted by a simple power law P(\cos\theta) proportional to \cos^4\theta with the same fitting parameters for all host halos.

To study the local Hubble flow, we have run constrained dark matter (DM) simulations of the Local Group (LG) in the concordance LCDM and OCDM cosmologies, with identical cosmological parameters apart from the Lambda term. The simulations were performed within a computational box of 64 h^{-1}Mpc centred on the LG. The initial conditions were constrained by the observed peculiar velocities of galaxies and positions of X-ray nearby clusters of galaxies. The simulations faithfully reproduce the nearby large scale structure, and in particular the Local Supercluster and the Virgo cluster. LG-like objects have been selected from the DM halos so as to closely resemble the dynamical properties of the LG. Both the LCDM and OCDM simulations show very similar local Hubble flow around the LG-like objects. It follows that, contrary to recent statements, the dark energy (DE) does not manifest itself in the local dynamics.

Using the high resolution cosmological N-body simulation MareNostrum Universe we study the orientation of shape and angular momentum of galaxy-size dark matter haloes around large voids. We find that haloes located on the shells of the largest cosmic voids have angular momenta that tend to be preferentially perpendicular to the direction that joins the centre of the halo and the centre of the void. This alignment has been found in spiral galaxies around voids using galaxy redshift surveys. We measure for the first time the strength of this alignment, showing how it falls off with increasing distance to the centre of the void. We also confirm the correlation between the intensity of this alignment and the halo mass. The analysis of the orientation of the halo main axes confirms the results of previous works. Moreover, we find a similar alignment for the baryonic matter inside dark matter haloes, which is much stronger in their inner parts.

We summarize the method of mass modelling of galaxy clusters based on reproducing the dispersion and kurtosis of the projected velocity distribution of galaxies. The models are parametrized within the framework of the NFW density profile, characterized by the virial mass and concentration, together with the constant anisotropy of galaxy orbits. The use of velocity dispersion alone does not allow to constrain all the three parameters from kinematic data due to the mass-anisotropy degeneracy. The degeneracy is broken by introducing the fourth velocity moment, the kurtosis. We tested the method based on fitting both moments on mock data sets drawn from simulated dark matter haloes and showed it to reproduce reliably the properties of the haloes. The method has been applied to estimate the mass, concentration and anisotropy of more than 20 clusters which allowed us to confirm, for the first time using kinematic data, the mass-concentration relation found in N-body simulations.

We investigate the origin and evolution of fossil groups in a concordance LCDM cosmological simulation. We consider haloes with masses between $(1-5)\times10^{13} \hMsun$ and study the physical mechanisms that lead to the formation of the large gap in magnitude between the brightest and the second most bright group member, which is typical for these fossil systems. Fossil groups are found to have high dark matter concentrations, which we can relate to their early formation time. The large magnitude-gaps arise after the groups have build up half of their final mass, due to merging of massive group members. We show that the existence of fossil systems is primarily driven by the relatively early infall of massive satellites, and that we do not find a strong environmental dependence for these systems. In addition, we find tentative evidence for fossil group satellites falling in on orbits with typically lower angular momentum, which might lead to a more efficient merger onto the host. We find a population of groups at higher redshifts that go through a ``fossil phase'': a stage where they show a large magnitude-gap, which is terminated by renewed infall from their environment.

The cosmic ultraviolet background (UVB) heats the intergalactic medium (IGM), as a result the gas in dark matter halos below a certain mass is too hot to cool within a Hubble time. The UVB effectively suppresses the formation of dwarf galaxies. Using high resolution cosmological hydrodynamical simulations we show that photo heating leads to small baryon fractions in halos below ~ 6x10^9 h^{-1}M_sun, independent of the cosmic environment. The simulations are carried out assuming a homogeneous UVB with flux densities as given by Haardt & Madau (1996). A halo may stop to condense gas significantly after the universe is reionised, namely when its mass falls below the characteristic mass scale set by the photo heating. Assuming a spherical halo model we derive this characteristic mass analytically and identify the main mechanisms that prevent the gas from cooling in small halos. The theoretically derived characteristic mass is smaller than the one obtained from observations. Increasing the energy per ionising photon by a factor between four and eight would be sufficient to reconcile both. This is equivalent to an average temperature of the IGM of ~ 10^4 K. In this sense the faint end of the luminosity function may serve as a calorimeter for the IGM.

We present the mass and X-ray temperature functions derived from a sample of more than 15,000 galaxy clusters of the MareNostrum Universe cosmological SPH simulations. In these simulations, we follow structure formation in a cubic volume of 500/h Mpc on a side assuming cosmological parameters consistent with either the first or third year WMAP data and gaussian initial conditions. We compare our numerical predictions with the most recent observational estimates of the cluster X-ray temperature functions and find that the low normalization cosmological model inferred from the 3 year WMAP data results is barely compatible with the present epoch X-ray cluster abundances. We can only reconcile the simulations with the observational data if we assume a normalization of the Mass-Temperature relation which is a factor of 2.5--3 smaller than our non-radiative simulations predict. This deviation seems to be too large to be accounted by the effects of star formation or cooling in the ICM, not taken into account in these simulations.

On the basis of a large scale 'adiabatic', namely non-radiative and non-dissipative, cosmological smooth particle hydrodynamic simulation we compare the entropy profiles of the gas and the dark matter (DM) in galaxy clusters. The quantity K_g = T_g \rho_g^{-2/3} provides a measure for the entropy of the intra-cluster gas. By analogy with the thermodynamic variables of the gas the velocity dispersion of the DM is associated with a formal temperature and thereby K_DM = \sigma_DM^2 \rho_DM^{-2/3} is defined. This DM entropy is related to the DM phase space density by K_DM \propto Q_DM^{-2/3}. In accord with other studies the DM phase space density follows a power law behaviour, Q_DM \propto r^{-1.82}, which corresponds to K_DM \propto r^{1.21}. The simulated intra-cluster gas has a flat entropy core within (0.8 \pm 0.4) R_s, where R_s is the NFW scale radius. The outer profile follows the DM behaviour, K_g \propto r^{1.21}, in close agreement with X-ray observations. Upon scaling the DM and gas densities by their mean cosmological values we find that outside the entropy core a constant ratio of K_g / K_{DM} = 0.71 \pm 0.18 prevails. By extending the definition of the gas temperature to include also the bulk kinetic energy the ratio of the DM and gas extended entropy is found to be unity for r > 0.8 R_s. The constant ratio of the gas thermal entropy to that of the DM implies that observations of the intra-cluster gas can provide an almost direct probe of the DM.

The MareNostrum Universe is one of the largest cosmological SPH simulation done so far. It consists of $1024^3$ dark and $1024^3$ gas particles in a box of 500 $h^{-1}$ Mpc on a side. Here we study the shapes and spins of the dark matter and gas components of the 10,000 most massive objects extracted from the simulation as well as the gas fraction in those objects. We find that the shapes of objects tend to be prolate both in the dark matter and gas. There is a clear dependence of shape on halo mass, the more massive ones being less spherical than the less massive objects. The gas distribution is nevertheless much more spherical than the dark matter, although the triaxiality parameters of gas and dark matter differ only by a few percent and it increases with cluster mass. The spin parameters of gas and dark matter can be well fitted by a lognormal distribution function. On average, the spin of gas is 1.4 larger than the spin of dark matter. We find a similar behavior for the spins at higher redshifts, with a slightly decrease of the spin ratios to 1.16 at $z=1.$ The cosmic normalized baryon fraction in the entire cluster sample ranges from $Y_b = 0.94$, at $z=1$ to $Y_b = 0.92$ at $z=0$. At both redshifts we find a slightly, but statistically significant decrease of $Y_b$ with cluster mass.

The aim of this paper is to study the efficiency of different approaches to interloper treatment in dynamical modelling of galaxy clusters. Using cosmological N-body simulation of standard LCDM model, we select 10 massive dark matter haloes and use their particles to emulate mock kinematic data in terms of projected galaxy positions and velocities as they would be measured by a distant observer. Taking advantage of the full 3D information available from the simulation, we select samples of interlopers defined with different criteria. The interlopers thus selected provide means to assess the efficiency of different interloper removal schemes found in the literature. We study direct methods of interloper removal based on dynamical or statistical restrictions imposed on ranges of positions and velocities available to cluster members. In determining these ranges, we use either the velocity dispersion criterion or a maximum velocity profile. We also generalize the common approaches taking into account both the position and velocity information. Another criterion is based on the dependence of the commonly used virial mass and projected mass estimators on the presence of interlopers. We find that the direct methods exclude on average 60-70 percent of unbound particles producing a sample with contamination as low as 2-4 percent. Next, we consider indirect methods of interloper treatment which are applied to the data stacked from many objects. In these approaches, interlopers are treated in a statistical way as a uniform background which modifies the distribution of cluster members. Using a Bayesian approach, we reproduce the properties of composite clusters and estimate the probability of finding an interloper as a function of distance from the object centre.

Within the cosmic framework clusters of galaxies are relatively young objects. Many of them have recently experienced major mergers. Here we investigate an equal mass merging event at z = ~0.6 resulting in a dark matter haloe of ~2.2 times 10^{14} Msol/h at z=0. The merging process is covered by 270 outputs of a high resolution cosmological N-body simulation performed with the ART (adaptive refinement tree) code. Some 2 Gyrs elapse between the first peri-centre passage of the progenitor cores and their final coalescence. During that phase the cores experience six peri-centre passages with minimal distances declining from ~30 to ~2 kpc/h. The time intervals between the peri-centre passages continuously decrease from 9 to 1 times 10^8 yrs. We follow the mean density, the velocity dispersion and the entropy of the two progenitors within a set of fixed proper radii (25, 50, 100, 250, 500, 1000 kpc/h). During the peri-centre passages we find sharp peaks of the mean densities within these radii, which exceed the sum of the corresponding progenitor densities. In addition to the intermixing of the merging haloes, the densities increase due to contraction caused by the momentary deepening of the potential well. At the end of the oscillatory relaxation phase the material originating from the less concentrated of the two equal mass progenitors is deposited at larger radii and shows a slightly more radially anisotropic velocity dispersion compared to the material coming from the more concentrated progenitor. Every peri-centre passage is accompanied by a substantial drop of the central potential well. We briefly discuss the possibility that AGN outbursts are triggered by the periodically changing potential.

We report some results from one of the largest hydrodynamical cosmological simulations of large scale structures that has been done up to date. The MareNostrum Universe SPH simulation consists of 2 billion particles (2 times 1024^3) in a cubic box of 500 h^-1 Mpc on a side. This simulation has been done in the MareNostrum parallel supercomputer at the Barcelona SuperComputer Center. Due to the large simulated volume and good mass resolution, our simulated catalog of dark matter halos comprises more than half a million objects with masses larger than a typical Milky Way galaxy halo. From this dataset we have studied several statistical properties such as the evolution of the halo mass function, the void distribution, the shapes of dark and gas halos and the large scale distribution of baryons.

We test the Secondary Infall Model (SIM) by direct comparison with the results of N-body simulations. Eight cluster-size and six galactic-size dark matter haloes have been selected at $z=0$ and re-simulated with high resolution. Based on their density profiles at the initial redshift, we compute their evolution by the SIM, assuming a simple prescription for the angular momentum. A comparison of the density profiles obtained by the SIM and the numerical experiments at $z=5$, 1 and 0 shows that, for most of the haloes at most epochs, the SIM reproduces the simulated mater distribution with a typical fractional deviation of less than 40 per cent over more than six order of magnitudes in the density. It is also found that, within the SIM framework, most of the diversity in the shape of the density profiles at $z=0$ arises from the scatter in the primordial initial conditions rather than the scatter in the angular momentum distribution. A crude optimization shows that a similar degree of agreement is obtained for galactic and cluster haloes, but the former seem to require slightly higher amounts of angular momentum. Our main conclusion is that the SIM provides a viable dynamical model for predicting the structure and evolution of the density profile of dark matter haloes.

We show by means of a high-resolution N-body simulation how the mass assembly histories of galaxy-size cold dark matter (CDM) halos depend on environment. Halos in high density environments form earlier and a higher fraction of their mass is assembled in major mergers, compared to low density environments. The distribution of the present-day specific mass aggregation rate is strongly dependent on environment. While in low density environments only ~20% of the halos are not accreting mass at the present epoch, this fraction rises to ~80% at high densities. At z=1 the median of the specific aggregation rate is ~4 times larger than at z=0 and almost independent on environment. All the dependences on environment found here are critically enhanced by local processes associated to subhalos because the fraction of subhalos increases as the environment gets denser. The distribution of the halo specific mass aggregation rate as well as its dependence on environment resemble the relations for the specific star formation rate distribution of galaxies. An analogue of the morphology-density relation is also present at the level of CDM halos, being driven by the halo major merging history. Nevertheless, baryonic processes are necessary in order to explain further details and the evolution of the star formation rate-, color- and morphology-environment relations.

The MareNostrum Universe is one of the biggest SPH cosmological simulations done so far. It contains more than 2 billion particles (2 times 1024^3) in a 500 Mpc/h cubic volume. This simulation has been performed on the MareNostrum supercomputer at the Barcelona Supercomputer Center. We have obtained more than 0.5 million halos with masses greater than a typical Milky Way galaxy halo. We report results about the halo mass function, the shapes of dark matter and gas distributions in halos, the baryonic fraction in galaxy clusters and groups, baryon oscillations in the dark matter and the halo power spectra as well as the distribution and evolution of the gas fraction at large scales.

We study structure formation in cosmological void regions using high-resolution hydrodynamical simulations. Despite being significantly underdense, voids are populated abundantly with small dark matter halos which should appear as dwarf galaxies if their star formation is not suppressed significantly. We here investigate to which extent the cosmological UV-background photo-evaporates baryons out of halos of dwarf galaxies, and thereby limits their cooling and star formation rates. Assuming a Haardt & Madau UV-background with reionisation at redshift z=6, our samples of simulated galaxies show that halos with masses below a characteristic mass of M_c(z=0) = 6.5 x 10^9 h^{-1} M_sun are baryon-poor, but in general not completely empty, because baryons that are in the condensed cold phase or are already locked up in stars resist evaporation. In halos with mass M < M_c, we find that photo-heating suppresses further cooling of gas. The redshift and UV-background dependent characteristic mass M_c(z) can be understood from the equilibrium temperature between heating and cooling at a characteristic overdensity of \delta ~ 1000. If a halo is massive enough to compress gas to this density despite the presence of the UV background, gas is free to `enter' the condensed phase and cooling continues in the halo, otherwise it stalls. By analysing the mass accretion histories of dwarf galaxies in voids, we show that they can build up a significant amount of condensed mass at early times before the epoch of reionisation. Later on, the amount of mass in this phase remains roughly constant, but the masses of the dark matter halos continue to increase. (abridged)

The aim of the present work is to show that, contrary to popular belief, galaxy clusters are **not** expected to be self-similar, even when the only energy sources available are gravity and shock-wave heating. In particular, we investigate the scaling relations between mass, luminosity and temperature of galaxy groups and clusters in the absence of radiative processes. Theoretical expectations are derived from a polytropic model of the intracluster medium and compared with the results of high-resolution adiabatic gasdynamical simulations. It is shown that, in addition to the well-known relation between the mass and concentration of the dark matter halo, the effective polytropic index of the gas also varies systematically with cluster mass, and therefore neither the dark matter nor the gas profiles are exactly self-similar. It is remarkable, though, that the effects of concentration and polytropic index tend to cancel each other, leading to scaling relations whose logarithmic slopes roughly match the predictions of the most basic self-similar models. We provide a phenomenological fit to the relation between polytropic index and concentration, as well as a self-consistent scheme to derive the non-linear scaling relations expected for any cosmology and the best-fit normalizations of the M-T, L-T and F-T relations appropriate for a Lambda-CDM universe. The predicted scaling relations reproduce observational data reasonably well for massive clusters, where the effects of cooling and star formation are expected to play a minor role.

We study the mass distribution in six nearby (z<0.06) relaxed Abell clusters of galaxies A0262, A0496, A1060, A2199, A3158 and A3558. Given the dominance of dark matter in galaxy clusters we approximate their total density distribution by the NFW formula characterized by virial mass and concentration. We also assume that the anisotropy of galactic orbits is reasonably well described by a constant and that galaxy distribution traces that of the total density. Using the velocity and position data for 120-420 galaxies per cluster we calculate, after removal of interlopers, the profiles of the lowest-order even velocity moments, dispersion and kurtosis. We then reproduce the velocity moments by jointly fitting the moments to the solutions of the Jeans equations. Including the kurtosis in the analysis allows us to break the degeneracy between the mass distribution and anisotropy and constrain the anisotropy as well as the virial mass and concentration. The method is tested in detail on mock data extracted from N-body simulations of dark matter haloes. We find that the best-fitting galactic orbits are remarkably close to isotropic in most clusters. Using the fitted pairs of mass and concentration parameters for the six clusters we conclude that the trend of decreasing concentration for higher masses found in cosmological N-body simulations is consistent with the data. By scaling the individual cluster data by mass we combine them to create a composite cluster with 1465 galaxies and perform a similar analysis on such sample. The estimated concentration parameter then lies in the range 1.5 < c < 14 and the anisotropy parameter in the range -1.1 < \beta < 0.5 at the 95 percent confidence level.

Dynamical Dark Energy (DE) is a viable alternative to the cosmological constant. Yet, constructing tests to discriminate between Lambda and dynamical DE models is difficult because the differences are not large. In this paper we explore tests based on the galaxy mass function, the void probability function (VPF), and the number of galaxy clusters. At high z the number density of clusters shows large differences between DE models, but geometrical factors reduce the differences substantially. We find that detecting a model dependence in the cluster redshift distribution is a hard challenge. We show that the galaxy redshift distribution is potentially a more sensitive characteristics. We do so by populating dark matter halos in Nbody simulations with galaxies using well-tested Halo Occupation Distribution (HOD). We also estimate the Void Probability Function and find that, in samples with the same angular surface density of galaxies in different models, the VPF is almost model independent and cannot be used as a test for DE. Once again, geometry and cosmic evolution compensate each other. By comparing VPF's for samples with fixed galaxy mass limits, we find measurable differences.

A large scale SPH+N-body simulation (GADGET) of the concordance LCDM universe is used to investigate orientation and angular momentum of galaxy clusters at z=0 in connection with their recent accretion histories. The basic cluster sample comprises the 3000 most massive friends-of-friends halos found in the 500 Mpc/h simulation box. Two disjoint sub-samples are constructed, using the mass ratio of the two most massive progenitors at z=0.5 m_2 / m_1 (m_1 < m_2), namely a recent major merger sample and a steady accretion mode sample. The mass of clusters in the merger sample is on average ~43% larger than the mass of the two progenitors (m_1 + m_2), whereas in the steady accretion mode sample a smaller increase of ~25% is found. The separation vector connecting the two most massive progenitor halos at z=0.5 is strongly correlated with the orientation of the cluster at z=0. The angular momentum of the clusters in the recent major merger sample tends to be parallel to orbital angular momentum of the two progenitors, whereas the angular momentum of the steady accretion mode sample is mainly determined by the angular momentum of the most massive progenitor. The long range correlations for the major and the minor principal axes of cluster pairs extend to distances of ~100 Mpc/h. Weak angular momentum correlations are found for distances < 20 Mpc/h. Within these ranges the major axes tend to be aligned with the connecting line of the cluster pairs whereas minor axes and angular momenta tend to be perpendicular to this line. A separate analysis of the two sub-samples reveals that the long range correlations are independent of the mass accretion mode. Thus orientation and angular momentum of galaxy clusters is mainly determined by the accretion along the filaments independently of the particular accretion mode.

Using high resolution DM simulations we study the shape of dark matter halos. Halos become more spherical with decreasing mass. This trend is even more pronounced for the inner part of the halo. Angular momentum and shape are correlated. The angular momenta of neighboring halos are correlated.

The Abell 1689 galaxy cluster has recently become a subject of intensive study. Thanks to its intermediate redshift (z=0.183) its mass distribution can be reconstructed using numerous methods including gravitational lensing, galaxy kinematics and X-ray imaging. The methods used to yield conflicting mass estimates in the past and recently the cluster mass distribution has been claimed to be in conflict with standard CDM scenarios due to rather large concentration and steep mass profile obtained from detailed studies of Broadhurst et al. using lensing. By studying in detail the kinematics of about 200 galaxies with measured redshifts in the vicinity of the cluster we show that the cluster is probably surrounded by a few structures, quite distant from each other, but aligned along the line of sight. We support our arguments by referring to cosmological N-body simulations and showing explicitly that distant, non-interacting haloes can produce entangled multi-peak line-of-sight velocity distributions similar to that in A1689. We conclude that it is difficult to estimate the cluster mass reliably from galaxy kinematics, but the value we obtain after applying a simple cut-off in velocity agrees roughly with the mass estimated from lensing. The complicated mass distribution around the cluster may however increase the uncertainty in the determination of the density profile shape obtained with weak lensing.

Using cosmological N-body simulations we study the line-of-sight velocity distribution of dark matter haloes focusing on the lowest-order even moments, dispersion and kurtosis, and their application to estimate the mass profiles of cosmological structures. For each of the ten massive haloes selected from the simulation box we determine the virial mass, concentration and the anisotropy parameter. In order to emulate observations from each halo we choose randomly 300 particles and project their velocities and positions along the line of sight and on the surface of the sky, respectively. After removing interlopers we calculate the profiles of the line-of-sight velocity moments and fit them with the solutions of the Jeans equations. The estimates of virial mass, concentration parameter and velocity anisotropy obtained in this way are in good agreement with the values found from the full 3D analysis.

High-resolution LCDM cosmological N-body simulations are used to study the properties of galaxy-size dark halos in different environments (cluster, void, and "field"). Halos in clusters and their surroundings have a median spin parameter ~1.3 times lower, and tend to be more spherical and to have less aligned internal angular momentum than halos in voids and the field. For halos in clusters the concentration parameters decrease on average with mass with a slope of ~0.1; for halos in voids these concentrations do not change with mass. For masses <5 10^11 M_sh^-1, halos in clusters are on average ~30-40% more concentrated and have ~2 times higher central densities than halos in voids. When comparing only parent halos, the differences are less pronounced but they are still significant. The Vmax-and Vrms-mass relations are shallower and more scattered for halos in clusters than in voids, and for a given Vmax or Vrms, the mass is smaller at z=1 than at z=0 in all the environments. At z=1, the differences in the halo properties with environment almost dissapear, suggesting this that the differences were stablished mainly after z~1. The halos in clusters undergo more dramatic changes than those in the field or the voids. The differences with environment are owing to (i) the dependence of halo formation time on environment, and (ii) local effects as tidal stripping and the tumultuos histories that halos suffer in high-density regions. We calculate seminumerical models of disk galaxy evolution in halos with the properties found for the different environments. For a given disk mass, the galaxy disks have higher surface density, larger Vd,max and secular bulge-to-disk ratio, lower gas fraction, and are redder as one goes from cluster to void environments, in rough agreement with observations. (abridged)

We study the density profiles of collapsed galaxy-size dark matter halos with masses 1e11-5e12 Msun focusing mostly on the halo outer regions from the formal virial radius Rvir up to 5-7Rvir. We find that isolated halos in this mass range extend well beyond Rvir exhibiting all properties of virialized objects up to 2-3Rvir: relatively smooth density profiles and no systematic infall velocities. The dark matter halos in this mass range do not grow as one naively may expect through a steady accretion of satellites, i.e., on average there is no mass infall. This is strikingly different from more massive halos, which have large infall velocities outside of the virial radius. We provide accurate fit for the density profile of these galaxy-size halos. For a wide range (0.01-2)Rvir of radii the halo density profiles are fit with the approximation rho=rho_s exp(-2n[x^{1/n}-1])+rho_m, where x=r/r_s, rho_m is the mean matter density of the Universe, and the index n is in the range n=6-7.5. These profiles do not show a sudden change of behavior beyond the virial radius. For larger radii we combine the statistics of the initial fluctuations with the spherical collapse model to obtain predictions for the mean and most probable density profiles for halos of several masses. The model give excellent results beyond 2-3 formal virial radii.

We present the results of a monitoring campaign of the double quasar SBS1520+530 at Maidanak observatory from April 2003 to August 2004. We obtained light curves in V and R filters that show small-amplitude \Delta m~0.1 mag intrinsic variations of the quasar on time scales of about 100 days. The data set is consistent with the previously determined time delay of \Delta t=(130+-3) days by Burud et al. (2002). We find that the time delay corrected magnitude difference between the quasar images is now larger by (0.14+-0.03) mag than during the observations by Burud et al. (2002). This confirms the presence of gravitational microlensing variations in this system.

Using cosmological N-body simulations we study the radial velocity distribution in dark matter haloes focusing on the lowest-order even moments, dispersion and kurtosis. We determine the properties of ten massive haloes in the simulation box approximating their density distribution by the NFW formula characterized by the virial mass and concentration. We also calculate the velocity anisotropy parameter of the haloes and find it mildly radial and increasing with distance from the halo centre. The radial velocity dispersion of the haloes shows a characteristic profile with a maximum, while the radial kurtosis profile decreases with distance starting from a value close to Gaussian near the centre. We therefore confirm that dark matter haloes possess intrinsically non-Gaussian, flat-topped velocity distributions. We find that the radial velocity moments of the simulated haloes are quite well reproduced by the solutions of the Jeans equations obtained for the halo parameters with the anisotropy measured in the simulations. We also study the radial velocity moments for a composite cluster made of ten haloes out to ten virial radii. In this region the velocity dispersion decreases systematically to reach the value of the background, while kurtosis increases from below to above the Gaussian value of 3 signifying a transition from a flat-topped to a strongly peaked velocity distribution with respect to the Gaussian, which can be interpreted as the dominance of ordered flow with a small dispersion. We illustrate the transition by showing explicitly the velocity distribution of the composite cluster in a few radial bins.

We study motions of galaxies in galaxy clusters formed in the concordance LCDM cosmology. We use high-resolution cosmological simulations that follow dynamics of dark matter and gas and include various physical processes critical for galaxy formation: gas cooling, heating and star formation. Analysing motions of galaxies and the properties of intracluster gas in the sample of eight simulated clusters at z=0, we study velocity dispersion profiles of the dark matter, gas, and galaxies. We measure the mean velocity of galaxy motions and gas sound speed as a function of radius and calculate the average Mach number of galaxy motions. The simulations show that galaxies, on average, move supersonically with the average Mach number of ~1.4, approximately independent of the cluster-centric radius. The supersonic motions of galaxies may potentially provide an important source of heating for the intracluster gas by driving weak shocks and via dynamical friction, although these heating processes appear to be inefficient in our simulations. We also find that galaxies move faster than the dark matter particles in clusters. The magnitude of the velocity bias, b_v~1.1, is, however, smaller than the bias estimated for subhalos in dissipationless simulations. Interestingly, we find velocity bias in the tangential component of the velocity dispersion, but not in the radial component. Finally, we find significant random bulk motions of gas. The typical gas velocities are of order ~20-30% of the gas sound speed. These random motions provide about 10% of the total pressure support in our simulated clusters. The non-thermal pressure support, if neglected, will bias measurements of the total mass in the hydrostatic analyses of the X-ray cluster observations.

We analyze the mass accretion histories (MAHs) and density profiles of cluster- size halos with virial masses of 0.6-2.5x10^14/h Msun in a flat LCDM cosmology. In agreement with previous studies,we find that the concentration of the density distribution is tightly correlated with the halo's MAH and its formation redshift.During the period of fast mass growth the concentration remains approximately constant and low c_v ~ 3-4,while during the slow accretion stages it increases with decreasing redshift as c_v ~ (1+z)^-1.We consider fits of three widely discussed analytic density profiles to the simulated clusters focusing on the most relaxed inner regions.We find that there is no unique best fit analytic profile for all the systems.If,however,a cluster is best fit by a particular analytic profile at z=0,the same is usually true at earlier epochs out to z ~ 1-2.The local logarithmic slope of the density profiles at 3% of the virial radius ranges from -1.2 to -2.0,a remarkable diversity for the relatively narrow mass range of our cluster sample.For all the studied clusters the logarithmic slope becomes shallower with decreasing radius without reaching an asymptotic value down to the smallest resolved scale < ~1% of the virial radius.We do not find a clear correlation of the inner slope with the formation redshift or the shape of the halo's MAH.We do find, however,that during the period of rapid mass growth the density profiles can be well described by a single power law rho(r) ~ r^-gamma with gamma ~ 1.5-2.The relatively shallow power law slopes result in low concentrations at these stages of evolution.This indicates that the inner power law like density distribution of halos is built up during the periods of rapid mass accretion and active merging,while outer steeper profile is formed when the mass accretion is slow.

We analyze the halo occupation distribution (HOD), the probability for a halo of mass M to host a number of subhalos N, and two-point correlation function of galaxy-size dark matter halos using high-resolution dissipationless simulations of the concordance flat LCDM model. The halo samples include both the host halos and the subhalos, distinct gravitationally-bound halos within the virialized regions of larger host systems. We find that the first moment of the HOD, <N>(M), has a complicated shape consisting of a step, a shoulder, and a power law high-mass tail. The HOD can be described by a Poisson statistics at high halo masses but becomes sub-Poisson for <N><4. We show that the HOD can be understood as a combination of the probability for a halo of mass M to host a central galaxy and the probability to host a given number Ns of satellite galaxies. The former can be approximated by a step-like function, while the latter can be well approximated by a Poisson distribution, fully specified by its first moment <Ns>(M). We find that <Ns>~M^b with b~1 for a wide range of number densities, redshifts, and different power spectrum normalizations. This formulation provides a simple but accurate model for the halo occupation distribution found in simulations. At z=0, the two-point correlation function (CF) of galactic halos can be well fit by a power law down to ~100/h kpc with an amplitude and slope similar to those of observed galaxies. At redshifts z>~1, we find significant departures from the power-law shape of the CF at small scales. If the deviations are as strong as indicated by our results, the assumption of the single power law often used in observational analyses of high-redshift clustering is likely to bias the estimates of the correlation length and slope of the correlation function.

Simple self-consistent models of galaxy groups and clusters are tested against the results of high-resolution adiabatic gasdynamical simulations. We investigate two models based on the existence of a 'universal' dark matter density profile and two versions of the beta-model. The mass distribution of relaxed clusters can be fitted by phenomenological formulae proposed in the literature. Haloes that have experienced a recent merging event are systematically less concentrated and show steeper profiles than relaxed objects near the centre. The hot X-ray emitting gas is found to be in approximate hydrostatic equilibrium with the dark matter potential, and it is well described by a polytropic equation of state. Analytic formulae for the gas density and temperature can be derived from these premises. Though able to reproduce the X-ray surface brightness, the beta-model is shown to provide a poor description of our numerical clusters. We find strong evidence of a 'universal' temperature profile that decreases by a factor of 2-3 from the centre to the virial radius, whereas baryon fraction and entropy are monotonically increasing functions. Numerical resolution and entropy conservation play a key role in the shapes of the profiles at small radii.

The radial mass distribution of dark matter haloes is investigated within the framework of the spherical infall model. We present a new formulation of spherical collapse including non-radial motions, and compare the analytical profiles with a set of high-resolution N-body simulations ranging from galactic to cluster scales. We argue that the dark matter density profile is entirely determined by the initial conditions, which are described by only two parameters: the height of the primordial peak and the smoothing scale. These are physically meaningful quantities in our model, related to the mass and formation time of the halo. Angular momentum is dominated by velocity dispersion, and it is responsible for the shape of the density profile near the centre. The phase-space density of our simulated haloes is well described by a power-law profile, rho/sigma^3 = 10^{1.46\pm0.04} (rho_c/Vvir^3) (r/Rvir)^{-1.90\pm0.05}. Setting the eccentricity of particle orbits according to the numerical results, our model is able to reproduce the mass distribution of individual haloes.

Numerous numerical studies indicate that dark matter halos show an almost universal radial density profile. The origin of the profile is still under debate. We investigate this topic and pay particular attention to the velocity dispersion profile. To this end we have performed high-resolution simulations with two independent codes, ART and {\sc Gadget}. The radial velocity dispersion can be approximated as function of the potential by $\sigma_r^2 = a (\Phi / \Phi_{\rm{out}})^\kappa (\Phi_{\rm{out}} - \Phi)$, where $\Phi_{\rm{out}}$ is the outer potential of the halo. For the parameters $a$ and $\kappa$ we find $a=0.29\pm0.04$ and $\kappa=0.41\pm0.03$. We find that the power-law asymptote $\sigma^2 \propto \Phi^\kappa$ is valid out to much larger distances from the halo center than any power asymptote for the density profile $\rho \propto r^{-n}$. The asymptotic slope $n(r \to 0)$ of the density profile is related to the exponent $\kappa$ via $n=2\kappa/(1+\kappa)$. Thus the value obtained for $\kappa$ from the available simulation data can be used to obtain an estimate of the density profile below presently resolved scales. We predict a continuously decreasing $n$ towards the halo center with the asymptotic value $n \lesssim 0.58$ at $r=0$.

Using high resolution N-body simulations we address the problem of emptiness of giant 20 Mpc/h diameter voids found in the distribution of bright galaxies. Are the voids filled by dwarf galaxies? Do cosmological models predict too many small dark matter haloes inside the voids? Can the problems of cosmological models on small scales be addressed by studying the abundance of dwarf galaxies inside voids? We find that voids in the distribution of 10^12 Msun/h haloes (expected galactic magnitudes ~ M_*) are almost the same as the voids in 10^11 Msun/h haloes. Yet, much smaller haloes with masses 10^9 Msun/h and circular velocities v_circ about 20 km/s readily fill the voids: there should be almost 1000 of these haloes in a 20 Mpc/h void. A typical void of diameter 20 Mpc/h contains about 50 haloes with v_circ > 50 km/s. The haloes are arranged in a pattern, which looks like a miniature Universe: it has the same structural elements as the large-scale structure of the galactic distribution of the Universe. There are filaments and voids; larger haloes are at the intersections of filaments. The only difference is that all masses are four orders of magnitude smaller. There is severe (anti)bias in the distribution of haloes, which depends on halo mass and on the distance from the centre of the void. Large haloes are more antibiased and have a tendency to form close to void boundaries. The mass function of haloes in voids is different from the ``normal'' mass function. It is much steeper for high masses resulting in very few M33-type galaxies (v_circ about 100 km/s). We present an analytical approximation for the mass function of haloes in voids.

The relative orientation of clusters' major elongation axes and clusters' angular momentum is studied using a large N-body simulation in a box of 500 Mpc/h base length for a standard LambdaCDM model. Employing the technique of mark correlation functions, we successfully separated the correlations in the orientation from the well known clustering signal traced by the two-point correlation function. The correlations in the orientation are highly significant for our sample of 3000 clusters. We found an alignment of neighboring clusters, i.e. an enhanced probability of the major elongation axes of neighboring cluster pairs to be in parallel with each other. At 10 Mpc/h separation the amplitude of this signal is ~10% above the value expected from random orientations, and it vanishes on scales larger than 15 Mpc/h. The ``filamentary'' alignment between clusters' major elongation axes and the lines pointing towards neighboring clusters shows even stronger deviations from random orientation, which can be detected out to scales of 100 Mpc/h, both in 2D and 3D analyses. Similarly, strong correlations of the angular momentum were seen. Also a clear signal in the scalar correlation of the absolute value of the angular momentum, the spin parameter and the mass was found. They extend up to 50 Mpc/h and have an amplitude of 40%, 15%, and 10% above a random distribution at 10 Mpc/h separation, respectively.

The theory of structure formation predicts that galaxies form within extended massive dark matter halos built from smaller pieces that collided and merged, resulting in the hierarchy of galaxies, groups, and clusters observed today. Here we present constrained numerical simulations designed to match the observed local universe as well as possible and to study the formation, evolution and present day properties of such dark matter halos in different environments. Simulations have been done with different versions of the Adaptive Refinement Tree (ART) N-body code, some of them including eulerian gasdynamics.

Using numerical simulations and the Sheth-Tormen approximation we study the mass function of dark matter halos in voids. We find that the void mass function is significantly lower and its shape is different than that of the field halos. We predict that in the standard LCDM model a void with radius 10 h^-1 Mpc should have 50 halos with circular velocity v_c > 50 km/s and 600 halos with v_c > 20 km/s.

We use a novel statistical tool, the mark correlation functions (MCFs), to study clustering of galaxy-size halos as a function of their properties and environment in a high-resolution numerical simulation of the LambdaCDM cosmology. We applied MCFs using several types of continuous and discrete marks: maximum circular velocity of halos, merger mark indicating whether halos experienced or not a major merger in their evolution history, and a stripping mark indicating whether the halo underwent a tidal stripping. We find that halos which experienced a relatively early (z>1) major merger or mass loss (due to tidal stripping) in their evolution histories are over-abundant in halo pairs with separations < 3 Mpc/h. This result can be interpreted as spatial segregation of halos with different merger histories, qualitatively similar to the morphological segregation in the observed galaxy distribution. The analysis presented in this paper demonstrate that MCFs provide powerful, yet algorithmically simple, quantitative measures of segregation in the spatial distribution of objects with respect to their various properties (marks).

The cosmic star formation history in Cold Dark Matter dominated cosmological scenarios is studied by means of hydrodynamical numerical simulations. In particular, we explore a low density model with a Lambda-term and two high density models with different power spectra, all of them being spatially flat. Our simulations employ a fully nonlinear N-body and Eulerian hydrodynamics algorithm with a model for star formation and supernovae feedback that depends on two phenomenological parameters determined in agreement with previous papers. We find a nearly constant star formation rate beyond z=1, and we discuss which facts may determine the decrease in the SFR from z=1 to the present epoch. The Lambda-term cosmology with realistic parameters for star formation and feedback best reproduces the observed star formation history.

We present cosmological simulations which closely mimic the real Universe within 100Mpc of the Local Group. The simulations, called Constrained Simulations, reproduce the large-scale density field with major nearby structures, including the Local Group, the Coma and Virgo clusters, the Great Attractor, the Perseus-Pices, and the Local Supercluster, in approximately correct locations. The MARK III survey of peculiar velocities of the observed structures inside 80Mpc/h sphere is used to constrain the initial conditions. Fourier modes on scales larger then 5Mpc/h are dominated by the constraints, while small scale waves are random. The main aim of this paper is the structure of the Local Supercluster (LSC; 30Mpc/h around the Virgo cluster) and the Local Group environment. We find that at the current epoch most of the mass (7.5e14Msun/h) of the LSC is located in a filament roughly centered on the Virgo cluster and extending over 40Mpc/h. The simulated Local Group (LG) is located in an adjacent smaller filament, which is not a part of the main body of the LSC, and has a peculiar velocity of 250kms toward the Virgo cluster. The peculiar velocity field in the LSC region is complicated and is drastically different from the field assumed in the Virgocentric infall models. The peculiar velocity flow in the vicinity of the LG in the simulation is ``cold'': the peculiar line-of-sight velocity dispersion within 7Mpc/h of the LG is less than 60km/s, comparable to the observed velocity dispersion of nearby galaxies.

According to the hierarchical scenario, galaxies form via merging and accretion of small objects. Using N-body simulations, we study the frequency of merging events in the history of the halos. We find that at z<~2 the merging rate of the overall halo population can be described by a simple power law (1+z)^3. The main emphasis of the paper is on the effects of environment of halos at the present epoch (z=0). We find that the halos located inside clusters have formed earlier (dz \approx 1) than isolated halos of the same mass. At low redshifts (z<1), the merger rate of cluster halos is 3 times lower than that of isolated halos and 2 times lower than merger rate of halos that end up in groups by z=0. At higher redshifts (z~1-4), progenitors of cluster and group halos have 3--5 times higher merger rates than isolated halos. We briefly discuss implications of our results for galaxy evolution in different environments.

The possibility of determining cosmological parameters on the basis of a wide set of observational data including the Abell-ACO cluster power spectrum and mass function, peculiar velocities of galaxies, the distribution of Ly-$\alpha$ clouds and CMB temperature fluctuations is analyzed. Using a $\chi^2$ minimization method, assuming $\Omega_{\Lambda}+\Omega_m =1$ and no contribution from gravity waves, we found that a tilted $\Lambda$MDM model with one sort of massive neutrinos and the parameters $n\approx 1.12$, $\Omega_m\approx 0.4$ ($\Omega_{\Lambda}\approx 0.6$), $\Omega_{cdm}\approx 0.3$, $\Omega_{\nu}\approx 0.06$, $\Omega_b\approx 0.04$ and $h\approx 0.7$ matches observational data best. $\Omega_{\nu}$ is higher for more species of massive neutrinos, $\sim 0.1$ for two and $\sim 0.13$ for three species. $\Omega_m$ raises by $\sim 0.08$ and $\sim 0.15$ respectively. Varying only a subset of parameters and fixing the others shows also that the observational data set used here rules out pure CDM models with $h\ge 0.5$, scale invariant primordial power spectrum, zero cosmological constant and spatial curvature at a very high confidence level, $>99.99%$. The corresponding class of MDM models are ruled out at $\sim 95%$ C.L. It is notable also that this data set determines the amplitude of scalar fluctuations approximately at the same level as COBE four-year data. It indicates that a possible tensor component in the COBE data cannot be very substantial.

The possibility of determining cosmological parameters on the basis of a wide set of observational data including the Abell-ACO cluster power spectrum and mass function, peculiar velocities of galaxies, the distribution of Ly-$\alpha$ clouds and CMB temperature fluctuations is analyzed. Using a $\chi^2$ minimization method, assuming $\Omega_{\Lambda}+\Omega_{\rm{matter}} =1 $ and no contribution from gravity waves, we found that a tilted $\Lambda$MDM model with one sort of massive neutrinos and the parameters $n=1.12\pm 0.10$, $\Omega_m=0.41\pm 0.11$ ($\Omega_{\Lambda}=0.59\pm0.11$), $\Omega_{cdm}=0.31\pm 0.15$, $\Omega_{\nu}=0.059\pm 0.028$, $\Omega_b=0.039\pm 0.014$ and $h=0.70\pm 0.12$ matches observational data best. The 1$\sigma$ (68.3%) confidence limits on each cosmological parameter, which are obtained by marginalizing over the other parameters, are $0.82\le n\le1.39$, $0.19\le\Omega_m\le 1$ ($0\le\Omega_{\Lambda}\le 0.81$), $0\le\Omega_{\nu}\le 0.17$, $0.021\le \Omega_b\le 0.13$ and $0.38\le h\le 0.85$ $1.5\le b_{cl}\le 3.5$. Here $b_{cl}$ is the cluster bias parameter. The best-fit parameters for 31 models which are inside of $1\sigma$ range of the best model are presented. It is shown also that observational data set used here rules out the class of CDM models with $h\ge 0.5$, scale invariant primordial power spectrum, zero cosmological constant and spatial curvature at very high confidence level, $>99.99%$. The corresponding class of MDM models are ruled out at $\sim 95%$ C.L.

We present a study of numerical effects in dissipationless cosmological simulations. The numerical effects are evaluated and studied by comparing results of a series of 64^3-particle simulations of varying force resolution and number of time steps, performed using three different N-body techniques: the Particle Mesh (PM), the adaptive P3M (AP3M) code, and the Adaptive Refinement Tree (ART) code. This study can therefore be interesting both as an analysis of numerical effects and as a systematic comparison of different codes. We find that the AP3M and the ART codes produce similar results, given that convergence is reached within the code type. We also find that numerical effects may affect the high-resolution simulations in ways that have not been discussed before. In particular, our study revealed the presence of two-body scattering, effects of which can be greatly amplified by inaccuracies of time integration. This process appears to affect the correlation function of matter, mass function and inner density of dark matter halos and other statistics at scales much larger than the force resolution, although different statistics may be affected in different fashion. We discuss the conditions at which strong two-body scattering is possible and discuss the choice of the force resolution and integration time step. Furthermore, we discuss recent claims that simulations with force softening smaller than the mean interparticle separation are not trustworthy and argue that this claim is incorrect in general and applies only to the phase-sensitive statistics. Our conclusion is that, depending on the choice of mass and force resolution and integration time step, a force resolution as small as 0.01 of the mean interparticle separation can be justified.

Using high resolution N-body simulations we have studied the formation and evolution of clusters and groups in a LCDM cosmological model. Groups of galaxies have been formed already before z=4. The total number of small bound systems and the total number of galaxies in these small systems rapidly decreases after z=1.5. The fraction of isolated galaxies remains approximately constant after z =1, whereas the fraction of galaxies in groups decreases. Merging of groups and accretion leads to cluster formation at z < 2. Some of the groups merge into large isolated halos.

We discuss the impact of the cosmological environment on the evolution of dark matter halos using a high-resolution simulation within a spatially flat LCDM cosmology.

We use a high-resolution dissipationless simulation to study the evolution of the dark matter and halo distributions in a spatially flat cosmological model dominated by a cosmological constant $\Lambda$ and cold dark matter ($\Lambda$CDM). In order to quantify the evolution of structure, we calculate the Minkowski functionals of the halos and the dark matter component at various redshifts. A comparison of Minkowski functionals and the more standard correlation function analysis shows that the Minkowski functionals contain information about correlation functions of arbitrary order. While little evolution of the Minkowski functionals of halos between $z=4$ and $z=0$ is observed, we find that the Minkowski functionals of the dark matter evolve rapidly with time. The difference of the Minkowski functionals of halos and dark matter can be interpreted as a scale dependent bias. This implies that scale-dependent bias is a property of not only the two-point halo correlation function, but also of correlation functions of higher order.

We present analysis of the evolution of dark matter halos in dense environments of groups and clusters in dissipationless cosmological simulations. The premature destruction of halos in such environments, known as the overmerging, reduces the predictive power of N-body simulations and makes difficult any comparison between models and observations. We analyze the possible processes that cause the overmerging and assess the extent to which this problem can be cured with current computer resources and codes. Using both analytic estimates and high resolution numerical simulations, we argue that the overmerging is mainly due to the lack of numerical resolution. We find that the force and mass resolution required for a simulated halo to survive in galaxy groups and clusters is extremely high and was almost never reached before: ~1-3 kpc and 10^8-10^9 Msun, respectively. We use the high-resolution Adaptive Refinement Tree (ART) N-body code to run cosmological simulations with the particle mass of \approx 2x10^8/h Msun} and the spatial resolution of \approx 1-2/h kpc, and show that in these simulations the halos do survive in regions that would appear overmerged with lower force resolution. Nevertheless, the halo identification in very dense environments remains a challenge even with the resolution this high. We present two new halo finding algorithms developed to identify both isolated and satellite halos that are stable (existed at previous moments) and gravitationally bound. To illustrate the use of the satellite halos that survive the overmerging, we present a series of halo statistics, that can be compared with those of observed galaxies. (Abridged)

We study the formation and evolution of the dark matter halos in a LCDM model. The dynamics of 16.8 million DM particles is followed numerically in a box of 60 Mpc/h with the dynamic range of 32,000 in spatial resolution. The high resolution of the simulation allows us to examine evolution of both isolated and satellite halos in a representative volume. We present results on the evolution of velocity function of DM halos, compare it with the Press-Schechter function and discuss the evolution of power spectrum of matter and halo distributions. The mass evolution of isolated virialized objects determined from the simulation is in good agreement with predictions of the extended Press-Schechter models. However, satellite halos evolve very different: for some of them the mass decreases with time, which happens if the halos fall into clusters or groups. We discuss the dependence of the correlation function of halo populations on their environment and merging history.

A significant fraction of mass in the universe is believed to be in the form of dark matter (DM). Due to gravitational instability, the DM collapses hierarchically into DM halos. In this contribution we present a study of the formation and evolution of such DM halos in a COBE-normalized spatially flat LCDM model (Omega_0=1-Omega_L=0.3; h=0.7) using high-resolution N-body simulations. The novelty of this study is use of the newly developed halo-finding algorithms to study the evolution of both isolated and satellite (located inside virial radii of larger group- and cluster-size systems) halos. The force and mass resolution required for a simulated halo to survive in the high-density environments typical of groups and clusters is high: ~1-3 kpc and \~10^9 M_sun, respectively. We use the high-resolution Adaptive Refinement Tree (ART) N-body code to follow the evolution of 256^3 dark matter particles with dynamic range in spatial resolution of 32,000 in a box of 60/h Mpc. We show that the correlation function of these halos is anti-biased with respect to the dark matter correlation function and is high and steeper than the correlation function of the isolated virialized objects. The correlation function evolves only mildly between z=3 and z=1. The mass evolution of isolated virialized objects determined from the simulation is in good agreement with prediction of semi-analytical models. The differences exist, however, if we include satellite halos in the halo catalogs.

Using extensive N-body simulations we estimate redshift space power spectra of clusters of galaxies for different cosmological models (SCDM, TCDM, CHDM, Lambda-CDM, OCDM, BSI, tau-CDM) and compare the results with observational data for Abell-ACO clusters. Our mock samples of galaxy clusters have the same geometry and selection functions as the observational sample which contains 417 clusters of galaxies in a double cone of galactic latitude |b| > 30 degrees up to a depth of 240 Mpc/h. The power spectrum has been estimated for wave numbers k in the range 0.03 < k < 0.2 h/Mpc. For k > k_max ~ 0.05 h/Mpc the power spectrum of the Abell-ACO clusters has a power-law shape, P(k)\propto k^n, with n ~ -1.9, while it changes sharply to a positive slope at k < k_max. By comparison with the mock catalogues SCDM, TCDM (n=0.9), and also OCDM with Omega_0 = 0.35 are rejected. Better agreement with observation can be found for the Lambda-CDM model with Omega_0 = 0.35 and h = 0.7 and the CHDM model with two degenerate neutrinos and Omega_HDM = 0.2 as well as for a CDM model with broken scale invariance (BSI) and the tau-CDM model. As for the peak in the Abell-ACO cluster power spectrum, we find that it does not represent a very unusual finding within the set of mock samples extracted from our simulations.

Three halo finder algorithms - hierarchical friends-of-friends, bound density maxima, and 6D minimum spanning tree - are discussed.

The large scale gravitational potential distribution and its influence on the large-scale matter clustering is considered on the basis of six simulations. It is found that the mean separation between zero levels of the potential along random straight lines coincides with the theoretical expectations, but it scatters largely. A strong link of the initial potential and the structure evolution is shown. It is found that the under-dense and over-dense regions correlate with regions of positive and negative gravitational potential at large redshifts. The over-dense regions arise due to a slow matter flow into the negative potential regions, where more pronounced non-linear structures appear. Such regions are related to the formation of huge super-large scale structures seen in the galaxy distribution.

We compare the distribution of internal velocity dispersions of galaxy clusters for an observational sample to those obtained from a set of N-body simulations of seven COBE-normalised cosmological scenarios: the standard CDM (SCDM) and a tilted (n=0.85) CDM (TCDM) model, a CHDM model with 25% of massive neutrinos, two low-density LCDM models with Omega_0=0.3 and 0.5, two open OCDM models with Omega_0=0.4 and 0.6. Simulated clusters are observed in projection so as to reproduce the main observational biases and are analysed by applying the same algorithm for interlopers removal and velocity dispersion estimate as for the reference observational sample. Velocity dispersions for individual clusters can be largely affected by observational biases in a model-dependent way: models in which clusters had less time to virialize show larger discrepancies between 3D and projected velocity dispersions. From the comparison with real clusters we find that both SCDM and TCDM largely overproduce clusters. The CHDM model marginally overproduces clusters and requires a somewhat larger sigma_8 than a purely CDM model in order to produce the same cluster abundance. The LCDM model with Omega_0=0.3 agrees with data, while the open model with Omega_0=0.4 and 0.6 underproduces and marginally overproduces clusters, respectively.

Presented are measurements of the observed redshift-space galaxy-galaxy autocorrelation function, xi(s), for the Las Campanas Redshift Survey (LCRS). For separations 2.0/h Mpc < s < 16.4/h Mpc, xi(s) can be approximated by a power law with slope of -1.52 +/- 0.03 and a correlation length of s_0 = (6.28 +\- 0.27)/h Mpc. A zero-crossing occurs on scales of roughly 30 - 40/h Mpc. On larger scales, xi(s) fluctuates closely about zero, indicating a high level of uniformity in the galaxy distribution on these scales. In addition, two aspects of the LCRS selection criteria - a variable field-to-field galaxy sampling rate and a 55 arcsec galaxy pair separation limit - are tested and found to have little impact on the measurement of xi(s). Finally, the LCRS xi(s) is compared with those from numerical simulations; it is concluded that, although the LCRS xi(s) does not discriminate sharply among modern cosmological models, redshift-space distortions in the LCRS xi(s) will likely provide a strong test of theory.

The large-scale matter distribution represents a complex network of structure elements such as voids, clusters, filaments, and sheets. This network is spanned by a point distribution. The global properties of the point process can be measured by different statistical methods, which, however, do not describe directly the structure elements. The morphology of structure elements is an important property of the point distribution. Here we apply the core-sampling method to various Voronoi tessellations. Using the core-sampling method we identify one- and two-dimensional structure elements (filaments and sheets) in these Voronoi tessellations and reconstruct their mean separation along random straight lines. We compare the results of the core-sampling method with the a priori known structure elements of the Voronoi tessellations under consideration and find good agreement between the expected and found structure parameters, even in the presence of substantial noise. We conclude that the core-sampling method is a potentially powerful tool to investigate the distribution of such structure elements like filaments and walls of galaxies.

We use the void probability function to compare the redshift--space galaxy distribution in the Perseus--Pisces survey with artificial samples from N-body simulations of standard cold dark matter (CDM) and broken scale invariance (BSI) models. Galaxies are identified as residing in peaks of the evolved density field in such a way as to reproduce both the observed luminosity and two--point correlation functions. Using a similar approach, it was recently shown that the VPF can discriminate between CDM and a cold+hot dark matter (CHDM) model with cold/hot/baryon fractions = 0.6/0.3/0.1. Our main result is that both CDM (as expected from a previous analysis) and BSI fit observational data. The robustness of the result is checked against changing the observer's position in the simulations and the galaxy identification in the evolved density field. Therefore, while the void statistics is sensitive to the passage from CDM to CHDM (different spectrum and different nature of dark matter), it is not to the passage from CDM to BSI (different spectrum but same dark matter). On such a basis, we conjecture that the distribution of voids could be directly sensitive to the nature of dark matter, but scarcely sensitive to the shape of the transfer function.

We report on the performance of Lagrangian perturbation theory up to the second order for the standard cold dark matter (SCDM) and broken scale invariance (BSI) scenarios. We normalize both models to the COBE data, the BSI model serves as an example of models which fit the small-scale power of galaxy surveys. We optimize Lagrangian perturbation solutions by removing small-scale power from the initial data and compare the results with those of numerical simulations. We find an excellent performance of the optimized Lagrangian schemes down to scales around the correlation length or smaller, depending on the statistics used for the comparison. The optimization scheme can be expressed in a way which is independent of the type of fluctuation spectrum and of the size of the simulations.

We investigate the performance of Lagrangian perturbation theory up to the second order for two scenarios of cosmological large-scale structure formation, SCDM (standard cold dark matter) and BSI (broken scale invariance). The latter model we study as a representative of COBE-normalized CDM models which fit the small-scale power of galaxy surveys. In this context we optimize the performance of the Lagrangian perturbation schemes by smoothing the small-scale fluctuations in the initial data. The results of the so obtained Lagrangian mappings are computed for a set of COBE-normalized SCDM and BSI initial data of different sizes and at different times. We compare these results against those obtained with a numerical PM-code. We find an excellent performance of the optimized Lagrangian schemes down to scales close to the correlation length. This is explained by the counterintuitive fact that nonlinearities in the model can produce more small-scale power, if initially such power is removed. The optimization scheme can be expressed in a way which is independent of the type of fluctuation spectrum and of the size of the simulations.

A quantitative theory, based on the Zeldovich approximation, to provide an approximate description of the evolution of structure is presented. We give an expression for the characteristic scale of superlarge-scale structure which can also be applied to power spectra with a Harrison-Zeldovich asymptotic. The evolution of the characteristic scale of large-scale structure is described with an analytical expression. It corresponds to the mean separation of filaments or chains of galaxies. We apply the theory to some simulations and show that the theoretically defined characteristic scale for LSS does also correspond to the 'physical measurement' of the mean separation of dark matter LSS elements. This supports our claim for a new theory for the quantitative description of the formation and evolution of large-scale structure in the universe.

Observational constraints on standard CDM spectra and perturbation spectra with broken scale invariance are discussed.

To test the primordial power spectra predicted by a double inflationary model with a break of amplitude $\Delta=3$ at a scale of $2\pi/k\approx 10 \hm$ and CDM as dominant matter content, we perform PM simulations with 128$^3$ particles on a 256$^3$ grid. The broken scale invariance of the power spectra explains the extra power observed in the large-scale matter distribution. COBE-normalized spectra and a linear biasing with $b\approx 2$ are shown to reproduce the reconstructed power spectra from the CfA catalog. Identifying galactic halos with overdensity of approximately two times the cell variance, we can fit the angular correlation function using both the Limber equation and creating a APM-like angular projection with the observed luminosity function. Finally, the higher order moments of the galaxy distribution are shown to fit reasonably well the observed values.

CDM models with non-scale-free step-like spectra of adiabatic perturbations produced in a realistic double inflationary model are compared with recent observational data. The model contains two additional free parameters relatively to the standard CDM model with the flat ($n=1$) initial spectrum. Results of the COBE experiment are used for the determination of a free overall spectrum normalization. Then predictions for the galaxy biasing parameter, the variance for "counts in cells", the galaxy angular correlation function, bulk flow peculiar velocities and the Mach number test are obtained. Also considered are conditions for galaxy and quasar formation. Observational data strongly restricts allowed values for the two remaining model parameters. However, a non-empty region for them satisfying all considered tests is found.