The relationship between galaxy properties and environment is a widely discussed topic within astrophysics. Here we use galaxy samples from hydrodynamical re-simulations to examine this relationship. We use the over-density ($\delta_1$) within a $1 h^{-1}{\rm Mpc}$ sphere around a galaxy to evaluate its environment. Then the relations between galaxy properties, such as specific star formation rate(sSFR), fraction of star forming galaxies, $g-r$ colour and $\delta_1$ are examined within three galactic samples formed from galaxies within large clusters, those in the vicinity of large clusters and those in the field. We find tight environmental correlations for these galaxy properties. In brief, galaxies in denser environments tend to be redder and are more likely to be quenched. This is consistent with observations. We find that although the sSFR decreases with $\delta_1$, this is mainly because that galaxies with higher stellar mass reside in environment with higher overdensity. At fixed over-density a galaxy's color is also independent of whether it lives within a cluster or within the field, but the relative fractions of the two samples varies dramatically with over-density and this drives an apparent evolution.

Accurately predicting structural properties of dark matter halos is one of the fundamental goals of modern cosmology. We use the new suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies. The MultiDark simulations cover a large range of masses 1e10-1e15Msun and volumes upto 50Gpc**3. The total number of dark matter halos in all the simulations exceeds 60 billion. We find that in order to understand the structure of dark matter halos and to make ~1% accurate predictions for density profiles, one needs to realize that halo concentration is more complex than the traditional ratio of the virial radius to the core radius in the NFW profile. For massive halos the averge density profile is far from the NFW shape and the concentration is defined by both the core radius and the shape parameter alpha in the Einasto approximation. Combining results from different redshifts, masses and cosmologies, we show that halos progress through three stages of evolution. (1) They start as rare density peaks that experience very fast and nearly radial infall. This radial infall brings mass closer to the center producing a high concentrated halo. Here, the halo concentration increases with the increasing halo mass and the concentration is defined by the alpha parameter with nearly constant core radius. Later halos slide into (2) the plateau regime where the accretion becomes less radial, but frequent mergers still affect even the central region. Now the concentration does not depend on halo mass. (3) Once the rate of accretion slows down, halos move into the domain of declining concentration-mass relation because new accretion piles up mass close to the virial radius while the core radius is staying constant. We provide accurate analytical fits to the numerical results for halo density profiles and concentrations.

We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.

Future dark energy experiments will require better and more accurate theoretical predictions for the baryonic acoustic oscillations (BAO) signature in the spectrum of cosmological perturbations. Here, we use large N-body simulations of the \LambdaCDM Planck cosmology to study any possible systematic shifts and damping in BAO due to the impact of nonlinear gravitational growth of structure, scale dependent and non-local bias, and redshift-space distortions. The effect of cosmic variance is largely reduced by dividing the tracer power spectrum by that from a BAO-free simulation starting with the same phases. This permits us to study with unprecedented accuracy (better than 0.02% for dark matter and 0.07% for low-bias halos) small shifts of the pristine BAO wavenumbers towards larger k, and non-linear damping of BAO wiggles in the power spectrum of dark matter and halo populations in the redshift range z=0-1. For dark matter, we provide an accurate parametrization of the evolution of \alpha as a function of the linear growth factor D(z). For halo samples, with bias ranging from 1.2 to 2.8, we measure a typical BAO shift of ~0.25%, observed in real-space, which does not show an appreciable evolution with redshift within the uncertainties. Moreover, we report a constant shift as a function of halo bias. We find a different evolution of the damping of the acoustic feature in all halo samples as compared to dark matter with haloes suffering less damping, and also find some weak dependence on bias. A larger BAO shift and damping is measured in redshift-space which can be well explained by linear theory due to redshift-space distortions. A clear modulation in phase with the acoustic scale is observed in the scale-dependent halo bias due to the presence of the baryonic acoustic oscillations.

Accurate predictions of the abundance and clustering of dark matter haloes play a key role in testing the standard cosmological model. Here, we investigate the accuracy of one of the leading methods of connecting the simulated dark matter haloes with observed galaxies -- the Halo Abundance Matching (HAM) technique. We show how to choose the optimal values of the mass and force resolution in large-volume N-body simulations so that they provide accurate estimates for correlation functions and circular velocities for haloes and their subhaloes -- crucial ingredients of the HAM method. At the 10% accuracy, results converge for 50 particles for haloes and 150 particles for progenitors subhaloes. In order to achieve this level of accuracy a number of conditions should be satisfied. The force resolution for the smallest resolved (sub)haloes should be in the range (0.1-0.3)rs, where rs is the scale radius of (sub)haloes. The number of particles for progenitors of subhaloes should be 150. We also demonstrate that the two-body scattering plays a minor role for the accuracy of N-body simulations thanks to the relatively small number of crossing-times of dark matter in haloes, and the limited force resolution of cosmological simulations.

This paper presents an analysis of the local peculiar velocity field based on the Wiener Filter reconstruction method. We used our currently available catalog of distance measurements containing 1,797 galaxies within 3000 km/s: Cosmicflows-1. The Wiener Filter method is used to recover the full 3D peculiar velocity field from the observed map of radial velocities and to recover the underlying linear density field. The velocity field within a data zone of 3000 km/s is decomposed into a local component that is generated within the data zone and a tidal one that is generated by the mass distribution outside that zone. The tidal component is characterized by a coherent flow toward the Norma-Hydra-Centaurus (Great Attractor) region while the local component is dominated by a flow toward the Virgo Cluster and away from the Local Void. A detailed analysis shows that the local flow is predominantly governed by the Local Void and the Virgo Cluster plays a lesser role. The analysis procedure was tested against a mock catalog. It is demonstrated that the Wiener Filter accurately recovers the input velocity field of the mock catalog on the scale of the extraction of distances and reasonably recovers the velocity field on significantly larger scales. The Bayesian Wiener Filter reconstruction is carried out within the ?CDM WMAP5 framework. The Wiener Filter reconstruction draws particular attention to the importance of voids in proximity to our neighborhood. The prominent structure of the Local Supercluster is wrapped in a horseshoe collar of under density with the Local Void as a major component.

We use dark matter haloes identified in the MareNostrum Universe and galaxy groups identified in the Sloan Data Release 7 galaxy catalogue, to study the relation between halo shape and halo dynamics, parametrizing out the mass of the systems. A strong shape-dynamics, independent of mass, correlation is present in the simulation data, which we find it to be due to different halo formation times. Early formation time haloes are, at the present epoch, more spherical and have higher velocity dispersions than late forming-time haloes. The halo shape-dynamics correlation, albeit weaker, survives the projection in 2D (ie., among projected shape and 1-D velocity dispersion). A similar shape-dynamics correlation, independent of mass, is also found in the SDSS DR7 groups of galaxies and in order to investigate its cause we have tested and used, as a proxy of the group formation time, a concentration parameter. We have found, as in the case of the simulated haloes, that less concentrated groups, corresponding to late formation times, have lower velocity dispersions and higher elongations than groups with higher values of concentration, corresponding to early formation times.

We measure the clustering of dark matter halos in a large set of collisionless cosmological simulations of the flat LCDM cosmology. Halos are identified using the spherical overdensity algorithm, which finds the mass around isolated peaks in the density field such that the mean density is Delta times the background. We calibrate fitting functions for the large scale bias that are adaptable to any value of Delta we examine. We find a ~6% scatter about our best fit bias relation. Our fitting functions couple to the halo mass functions of Tinker et. al. (2008) such that bias of all dark matter is normalized to unity. We demonstrate that the bias of massive, rare halos is higher than that predicted in the modified ellipsoidal collapse model of Sheth, Mo, & Tormen (2001), and approaches the predictions of the spherical collapse model for the rarest halos. Halo bias results based on friends-of-friends halos identified with linking length 0.2 are systematically lower than for halos with the canonical Delta=200 overdensity by ~10%. In contrast to our previous results on the mass function, we find that the universal bias function evolves very weakly with redshift, if at all. We use our numerical results, both for the mass function and the bias relation, to test the peak-background split model for halo bias. We find that the peak-background split achieves a reasonable agreement with the numerical results, but ~20% residuals remain, both at high and low masses.

Using constrained simulations of the local Universe for generic cold dark matter and for 1keV warm dark matter, we investigate the difference in the abundance of dark matter halos in the local environment. We find that the mass function within 20 Mpc/h of the Local Group is ~2 times larger than the universal mass function in the 10^9-10^13 M_odot/h mass range. Imposing the field of view of the on-going HI blind survey ALFALFA in our simulations, we predict that the velocity function in the Virgo-direction region exceeds the universal velocity function by a factor of 3. Furthermore, employing a scheme to translate the halo velocity function into a galaxy velocity function, we compare the simulation results with a sample of galaxies from the early catalog release of ALFALFA. We find that our simulations are able to reproduce the velocity function in the 80-300 km/s velocity range, having a value ~10 times larger than the universal velocity function in the Virgo-direction region. In the low velocity regime, 35-80 km/s, the warm dark matter simulation reproduces the observed flattening of the velocity function. On the contrary, the simulation with cold dark matter predicts a steep rise in the velocity function towards lower velocities; for V_max=35 km/s, it forecasts ~10 times more sources than the ones observed. If confirmed by the complete ALFALFA survey, our results indicate a potential problem for the cold dark matter paradigm or for the conventional assumptions about energetic feedback in dwarf galaxies.

We investigate the large-scale distribution of galaxy clusters taken from several X-ray catalogs. Different statistics of clustering like the conditional correlation function (CCF) and the minimal spanning tree (MST) as well as void statistics were used. Clusters show two distinct regimes of clustering: 1) on scales of superclusters (~40/h Mpc) the CCF is represented by a power law; 2) on larger scales a gradual transition to homogeneity (~100/h Mpc) is observed. We also present the correlation analysis of the galaxy distribution taken from DR6 SDSS main galaxy database. In case of galaxies the limiting scales of the different clustering regimes are 1)10-15/h Mpc; 2) 40-50/h Mpc. The differences in the characteristic scales and scaling exponents of the cluster and galaxy distribution can be naturally explained within the theory of biased structure formation. We compared the density contrasts of inhomogeneities in the cluster and galaxy distributions in the SDSS region. The estimation of the relative cluster-galaxy bias gives the value b = 5 +/- 2. The distribution of real clusters is compared to that of simulated (model) clusters (the MareNostrum Universe simulations). We selected a cluster sample from 500/h Mpc simulation box with WMAP3 cosmological parameters and sigma_8 = 0.8. We found a general agreement between the distribution of observed and simulated clusters. The differences are mainly due to the presences of the Shapley supercluster in the observed sample. On the basis of SDSS galaxy sample we study properties of the power law behavior showed by the CCF on small scales. We show that this phenomenon is quite complex, with significant scatter in scaling properties, and characterized by a non-trivial dependence on galaxy properties and environment.

We measure the mass function of dark matter halos in a large set of collisionless cosmological simulations of flat LCDM cosmology and investigate its evolution at z<~2. Halos are identified as isolated density peaks, and their masses are measured within a series of radii enclosing specific overdensities. We argue that these spherical overdensity masses are more directly linked to cluster observables than masses measured using the friends-of-friends algorithm (FOF), and are therefore preferable for accurate forecasts of halo abundances. Our simulation set allows us to calibrate the mass function at z=0 for virial masses in the range 10^{11} Msol/h < M < 10^{15} Msol/h, to <~ 5%. We derive fitting functions for the halo mass function in this mass range for a wide range of overdensities, both at z=0 and earlier epochs. In addition to these formulae, which improve on previous approximations by 10-20%, our main finding is that the mass function cannot be represented by a universal fitting function at this level of accuracy. The amplitude of the "universal" function decreases monotonically by ~20-50%, depending on the mass definition, from z=0 to 2.5. We also find evidence for redshift evolution in the overall shape of the mass function.

We study properties of dark matter halos at high redshifts z=2-10 for a vast range of masses with the emphasis on dwarf halos with masses 10^7-10^9 Msun/h. We find that the density profiles of relaxed dwarf halos are well fitted by the NFW profile and do not have cores. We compute the halo mass function and the halo spin parameter distribution and find that the former is very well reproduced by the Sheth & Tormen model while the latter is well fitted by a lognormal distribution with lambda_0 = 0.042 and sigma_lambda = 0.63. We estimate the distribution of concentrations for halos in mass range that covers six orders of magnitude from 10^7 Msun/h to 10^13} Msun/h, and find that the data are well reproduced by the model of Bullock et al. The extrapolation of our results to z = 0 predicts that present-day isolated dwarf halos should have a very large median concentration of ~ 35. We measure the subhalo circular velocity functions for halos with masses that range from 4.6 x 10^9 Msun/h to 10^13 Msun/h and find that they are similar when normalized to the circular velocity of the parent halo. Dwarf halos studied in this paper are many orders of magnitude smaller than well-studied cluster- and Milky Way-sized halos. Yet, in all respects the dwarfs are just down-scaled versions of the large halos. They are cuspy and, as expected, more concentrated. They have the same spin parameter distribution and follow the same mass function that was measured for large halos.

We present results from a set of high (512^3 effective resolution), and ultra-high (1024^3) SPH adiabatic cosmological simulations of cluster formation aimed at studying the internal structure of the intracluster medium (ICM). We derive a self-consistent analytical model of the structure of the intracluster medium (ICM). We discuss the radial structure and scaling relations expected from purely gravitational collapse, and show that the choice of a particular halo model can have important consequences on the interpretation of observational data. The validity of the approximations of hydrostatic equilibrium and a polytropic equation of state are checked against results of our simulations. The properties of the ICM are fully specified when a 'universal' profile is assumed for either the dark or the baryonic component. We also show the first results from an unprecedented large-scale simulation of 500 Mpc/h and 2 times 512^3 gas and dark matter particles. This experiment will make possible a detailed study of the large-scale distribution of clusters as a function of their X-ray properties.

We present some of the results of an ongoing collaboration to sudy the dynamical properties of galaxy clusters by means of high resolution adiabatic SPH cosmological simulations. Results from our numerical clusters have been tested against analytical models often used in X-ray observations: $\beta$ model (isothermal and polytropic) and those based on universal dark matter profiles. We find a universal temperature profile, in agreement with AMR gasdynamical simulations of galaxy clusters. Temperature decreases by a factor 2-3 from the center to virial radius. Therefore, isothermal models (e.g. $\beta$ model) give a very poor fit to simulated data. Moreover, gas entropy profiles deviate from a power law near the center, which is also in very good agreement with independent AMR simulations. Thus, if future X-ray observations confirm that gas in clusters has an extended isothermal core, then non-adiabatic physics would be required in order to explain it.

Formation of luminous matter in the Universe is a complicated process, which includes many processes and components. It is the vastly different scales involved in the process (from star formation on few parsec scales to galaxy clusters and superclusters on megaparsecs scales) and numerous ill-understood processes, which make the whole field a maze of unsolved, but exciting problems. We present new approximations for numerical treatment of multiphase ISM forming stars. The approximations were tested and calibrated using N-body+fluid numerical simulations. We specifically target issues related with effects of unresolved lumpinesses of the gas. We show that the degree of freedom is much smaller than naively expected because of self-regulating nature of the process of global star formation. One of the problems of numerical simulations is related with the uncertainties of approximation of the supernovae (SN) feedback. It is often assumed that the feedback is mainly due to momentum transferred by SN in to the ISM. We argue that this may not be true. We present a realistic example of gas actively forming stars with short cooling time for which the SF feedback is important, but the kinetic energy of the gas motion due to SN is negligibly small as compared with the thermal energy of gas.