LCDM cosmological models with Early Dark Energy (EDE) have been proposed to resolve tensions between the Hubble constant H0 = 100h km/s/Mpc measured locally, giving h ~ 0.73, and H0 deduced from Planck cosmic microwave background (CMB) and other early universe measurements plus LCDM, giving h ~ 0.67. EDE models do this by adding a scalar field that temporarily adds dark energy equal to about 10% of the cosmological energy density at the end of the radiation-dominated era at redshift z ~ 3500. Here we compare linear and nonlinear predictions of a Planck-normalized LCDM model including EDE giving h = 0.728 with those of standard Planck-normalized LCDM with h = 0.678. We find that nonlinear evolution reduces the differences between power spectra of fluctuations at low redshifts. As a result, at z = 0 the halo mass functions on galactic scales are nearly the same, with differences only 1-2%. However, the differences dramatically increase at high redshifts. The EDE model predicts 50% more massive clusters at z = 1 and twice more galaxy-mass halos at z = 4. Even greater increases in abundances of galaxy-mass halos at higher redshifts may make it easier to reionize the universe with EDE. Predicted galaxy abundances and clustering will soon be tested by JWST observations. Positions of baryonic acoustic oscillations (BAOs) and correlation functions differ by about 2% between the models -- an effect that is not washed out by nonlinearities. Both standard LCDM and the EDE model studied here agree well with presently available acoustic-scale observations, but DESI and Euclid measurements will provide stringent new tests.

Luminous red galaxies (LRGs) are one of the key tracers of the large-scale structure of the Universe used by galaxy surveys. Hence, it is important to make accurate predictions for their properties and clustering, including the errors on these statistics. Here, we describe a novel technique which uses the semi-analytical model of galaxy formation {\sc Galform}, embedded in the high-resolution $N$-body Planck-Millennium simulation, to populate a thousand halo catalogues generated using the Parallel-PM $N$-body {\sc glam} code. Our hybrid scheme allows us to make clustering predictions on scales that cannot be modelled in the original $N$-body simulation. LRGs are selected in the redshift range $z=0.6-1$ from the {\sc Galform} output using similar colour-magnitude cuts in the $r$, $z$ and $W1$ bands to those that will be applied in the Dark Energy Spectroscopic Instrument (DESI) survey. We find that the LRG-halo connection is non-trivial, leading to the prediction of a non-standard halo occupation distribution; in particular, the occupation of central galaxies does not reach unity for the most massive haloes, and drops with increasing mass. The {\sc glam} catalogues reproduce the abundance and clustering of the LRGs predicted by {\sc Galform}, and show good agreement with recent measurements of the clustering of DESI-like LRGs using photometric redshifts. We use the \glam{} mocks to compute the covariance matrices for the two-point correlation function and power spectrum of the LRGs and their background dark matter density field, revealing important differences. We also make predictions for the linear-growth rate and the baryon acoustic oscillations distances at $z=0.6$, $0.74$ and $0.93$. All DESI-like LRG catalogues are made publicly available.

We derive cosmological constraints on the matter density, $\Omega_m$, and the amplitude of fluctuations, $\sigma_8$, using $\mathtt{GalWCat19}$, a catalog of 1800 galaxy clusters we identified in the Sloan Digital Sky Survey-DR13 spectroscopic data set using our GalWeight technique to determine cluster membership (Abdullah et al. 2018; Abdullah et al. 2020). By analyzing a subsample of 843 clusters in the redshift range $0.01 \leq z \leq 0.125$ with virial masses of $M\geq 0.8\times10^{14}$ $h^{-1} \ M_{\odot}$, we obtain $\Omega_m=0.305^{+0.037}_{-0.042}$ and $\sigma_8=0.810^{+0.053}_{-0.056}$, with a cluster normalization relation of $\sigma_8= 0.44 \Omega_m^{-0.52}$. There are several unique aspects to our approach: we use the largest spectroscopic data set currently available, and we assign membership using the GalWeight technique which we have shown to be very effective at simultaneously maximizing the number of {\it{bona fide}} cluster members while minimizing the number of contaminating interlopers. Moreover, rather than employing scaling relations, we calculate cluster masses individually using the virial mass estimator. Since $\mathtt{GalWCat19}$ is a low-redshift cluster catalog we do not need to make any assumptions about evolution either in cosmological parameters or in the properties of the clusters themselves. Our constraints on $\Omega_m$ and $\sigma_8$ are consistent and very competitive with those obtained from non-cluster abundance cosmological probes such as Cosmic Microwave Background (CMB), Baryonic Acoustic Oscillation (BAO), and supernovae (SNe). The joint analysis of our cluster data with Planck18+BAO+Pantheon gives $\Omega_m=0.310^{+0.013}_{-0.010}$ and $\sigma_8=0.810^{+0.010}_{-0.013}$.

Utilizing the SDSS-DR13 spectroscopic dataset, we create a new publicly-available catalog of 1,800 galaxy clusters (GalWeight cluster catalog, $\mathtt{GalWCat19}$) and a corresponding catalog of 34,471 identified member galaxies. The clusters are identified from overdensities in redshift-phase space. The GalWeight technique introduced in Abdullah, Wilson and Klypin (AWK18) is then applied to identify cluster members. The completeness of the cluster catalog ($\mathtt{GalWCat19}$) and the procedure followed to determine cluster mass are tested on the Bolshoi N-body simulations. The 1,800 $\mathtt{GalWCat19}$ clusters range in redshift between $0.01 - 0.2$ and in mass between $(0.4 - 14) \times 10^{14}h^{-1}M_{\odot}$. The cluster catalog provides a large number of cluster parameters including sky position, redshift, membership, velocity dispersion, and mass at overdensities $\Delta = 500, 200, 100, 5.5$. The 34,471 member galaxies are identified within the radius at which the density is 200 times the critical density of the Universe. The galaxy catalog provides the coordinates of each galaxy and the ID of the cluster that the galaxy belongs to. The cluster velocity dispersion scales with mass as $\log(\sigma_{200})=\log(946\pm52~ \mbox{km} ~ \mbox{s}^{-1}) +(0.349\pm0.142)\log\left[h(z) ~ M_{200}/10^{15}M_\odot\right]$ with scatter of $\delta_{\log\sigma} = 0.06$. The catalogs are publicly available at the following website\footnote{\url{https://mohamed-elhashash-94.webself.net/galwcat/}}.

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.

In order to capture as much information as possible large galaxy surveys have been increasing their volume and redshift depth. To face this challenge theory has responded by making cosmological simulations of huge computational volumes with equally increasing the number of dark matter particles and supercomputing resources. Thus, it is taken for granted that the ideal situation is when a single computational box encompasses the whole effective volume of the observational survey, e.g., ~50 Gpch^3 for the DESI and Euclid surveys. Here we study the effects of missing long-waves in a finite volume using several relevant statistics: the abundance of dark matter halos, the PDF, the correlation function and power spectrum, and covariance matrices. Finite volume effects can substantially modify the results if the computational volumes are less than ~(500Mpch)^3. However, the effects become extremely small and practically can be ignored when the box-size exceeds ~1Gpch^3. We find that the average power spectra of dark matter fluctuations show remarkable lack of dependence on the computational box-size with less than 0.1% differences between 1Gpch and 4Gpch boxes. No measurable differences are expected for the halo mass functions for these volumes. The covariance matrices are scaled trivially with volume, and small corrections due to super-sample modes can be added. We conclude that there is no need to make those extremely large simulations when a box-size of 1-1.5Gpch is sufficient to fulfil most of the survey science requirements.

Making cosmological inferences from the observed galaxy clustering requires accurate predictions for the mean clustering statistics and their covariances. Those are affected by cosmic variance -- the statistical noise due to the finite number of harmonics. The cosmic variance can be suppressed by fixing the amplitudes of the harmonics instead of drawing them from a Gaussian distribution predicted by the inflation models. Initial realizations also can be generated in pairs with 180 degrees flipped phases to further reduce the variance. Here, we compare the consequences of using paired-and-fixed vs Gaussian initial conditions on the average dark matter clustering and covariance matrices predicted from N-body simulations. As in previous studies, we find no measurable differences between paired-and-fixed and Gaussian simulations for the average density distribution function, power spectrum and bispectrum. Yet, the covariances from paired-and-fixed simulations are suppressed in a complicated scale- and redshift-dependent way. The situation is particularly problematic on the scales of Baryon Acoustic Oscillations where the covariance matrix of the power spectrum is lower by only 20% compared to the Gaussian realizations, implying that there is not much of a reduction of the cosmic variance. The non-trivial suppression, combined with the fact that paired-and-fixed covariances are noisier than from Gaussian simulations, suggests that there is no path towards obtaining accurate covariance matrices from paired-and-fixed simulations. Because the covariances are crucial for the observational estimates of galaxy clustering statistics and cosmological parameters, paired-and-fixed simulations, though useful for some applications, cannot be used for the production of mock galaxy catalogs.

The most massive and luminous galaxies in the Universe serve as powerful probes to study the formation of structure, the assembly of mass, and cosmology. However, their detailed formation and evolution is still barely understood. Here we extract a sample of massive mock galaxies from the semi-analytical model of galaxy formation (SAM) GALACTICUS from the MultiDark-Galaxies, by replicating the CMASS photometric selection from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). The comparison of the GALACTICUS CMASS-mock with BOSS-CMASS data allows us to explore different aspects of the massive galaxy population at $0.5<z<0.6$, including the galaxy-halo connection and the galaxy clustering. We find good agreement between our modelled galaxies and observations regarding the galaxy-halo connection, but our CMASS-mock over-estimates the clustering amplitude of the 2-point correlation function, due to a smaller number density compared to BOSS, a lack of blue objects, and a small intrinsic scatter in stellar mass at fixed halo mass of $<0.1$ dex. To alleviate this problem, we construct an alternative mock catalogue mimicking the CMASS colour-magnitude distribution by randomly down-sampling the SAM catalogue. This CMASS-mock reproduces the clustering of CMASS galaxies within 1$\sigma$ and shows some environmental dependency of star formation properties that could be connected to the quenching of star formation and the assembly bias.

We present the first galaxy-OVI absorption kinematic study for 20 absorption systems (EW>0.1~{\AA}) associated with isolated galaxies (0.15$<z<$0.55) that have accurate redshifts and rotation curves obtained using Keck/ESI. Our sample is split into two azimuthal angle bins: major axis ($\Phi<25^{\circ}$) and minor axis ($\Phi>33^{\circ}$). OVI absorption along the galaxy major axis is not correlated with galaxy rotation kinematics, with only 1/10 systems that could be explained with rotation/accretion models. This is in contrast to co-rotation commonly observed for MgII absorption. OVI along the minor axis could be modeled by accelerating outflows but only for small opening angles, while the majority of the OVI is decelerating. Along both axes, stacked OVI profiles reside at the galaxy systemic velocity with the absorption kinematics spanning the entire dynamical range of their galaxies. The OVI found in AMR cosmological simulations exists within filaments and in halos of ~50 kpc surrounding galaxies. Simulations show that major axis OVI gas inflows along filaments and decelerates as it approaches the galaxy while increasing in its level of co-rotation. Minor axis outflows in the simulations are effective within 50-75 kpc beyond that they decelerate and fall back onto the galaxy. Although the simulations show clear OVI kinematic signatures they are not directly comparable to observations. When we compare kinematic signatures integrated through the entire simulated galaxy halo we find that these signatures are washed out due to full velocity distribution of OVI throughout the halo. We conclude that OVI alone does not serve as a useful kinematic indicator of gas accretion, outflows or star-formation and likely best probes the halo virial temperature.

Halo bias is the main link between the matter distribution and dark matter halos. In its simplest form, halo bias is determined by halo mass, but there are known additional dependencies on other halo properties which are of consequence for accurate modeling of galaxy clustering. Here we present the most precise measurement of these secondary-bias dependencies on halo age, concentration, and spin, for a wide range of halo masses spanning from 10$^{10.7}$ to 10$^{14.7}$ $h^{-1}$ M$_{\odot}$. At the high-mass end, we find no strong evidence of assembly bias for masses above M$_{vir}$ $\sim10^{14}$ $h^{-1}$ M$_{\odot}$. Secondary bias exists, however, for halo concentration and spin, up to cluster-size halos, in agreement with previous findings. For halo spin, we report, for the first time, two different regimes: above M$_{vir}\sim$10$^{11.5}$ $h^{-1}$ M$_{\odot}$, halos with larger values of spin have larger bias, at fixed mass, with the effect reaching almost a factor 2. This trend reverses below this characteristic mass. In addition to these results, we test, for the first time, the performance of a multi-tracer method for the determination of the relative bias between different subsets of halos. We show that this method increases significantly the signal-to-noise of the secondary-bias measurement as compared to a traditional approach. This analysis serves as the basis for follow-up applications of our multi-tracer method to real data.

Cosmological simulations are fundamental tools to study structure formation and the astrophysics of evolving structures, in particular clusters of galaxies. While hydrodynamical simulations cannot sample efficiently large volumes and explore different cosmologies at the same time, N-body simulations lack the baryonic physics that is crucial to determine the observed properties of clusters. One solution is to use (semi-)analytical models to implement the needed baryonic physics. In this way, we can generate the many mock universes that will be required to fully exploit future large sky surveys, such as that from the upcoming eROSITA X-ray telescope. We developed a phenomenological model based on observations of clusters to implement gas density and temperature information on the dark-matter-only halos of the MultiDark simulations. We generate several full-sky mock light-cones of clusters for the WMAP and Planck cosmologies, adopting different parameters in our phenomenological model of the intra-cluster medium. For one of these simulations and models, we also generate 100 light-cones corresponding to 100 random observers and explore the variance among them in several quantities. In this first paper on MultiDark mock galaxy cluster light-cones, we focus on presenting our methodology and discuss predictions for eROSITA, in particular, exploring the potential of angular power spectrum analyses of its detected (and undetected) cluster population to study X-ray scaling relations, the intra-cluster medium, and the composition of the cosmic X-ray background. We make publicly available on-line more than 400 GB of light-cones, which include the expected eROSITA count rate, on Skies & Universes (http://www.skiesanduniverses.org).

We introduce GalWeight, a new technique for assigning galaxy cluster membership. This technique is specifically designed to simultaneously maximize the number of bona fide cluster members while minimizing the number of contaminating interlopers. The GalWeight technique can be applied to both massive galaxy clusters and poor galaxy groups. Moreover, it is effective in identifying members in both the virial and infall regions with high efficiency. We apply the GalWeight technique to MDPL2 \& Bolshoi N-body simulations, and find that it is $> 98\%$ accurate in correctly assigning cluster membership. We show that GalWeight compares very favorably against four well-known existing cluster membership techniques (shifting gapper, den Hartog, caustic, SIM). We also apply the GalWeight technique to a sample of twelve Abell clusters (including the Coma cluster) using observations from the Sloan Digital Sky Survey. We end by discussing GalWeight's potential for other astrophysical applications.

In Montero-Dorta et al. 2017, we show that luminous red galaxies (LRGs) from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) at $z\sim0.55$ can be divided into two groups based on their star formation histories. So-called fast-growing LRGs assemble $80\%$ of their stellar mass at $z\sim5$, whereas slow-growing LRGs reach the same evolutionary state at $z\sim1.5$. We further demonstrate that these two subpopulations present significantly different clustering properties on scales of $\sim1 - 30 \mathrm{Mpc}$. Here, we measure the mean halo mass of each subsample using the galaxy-galaxy lensing technique, in the $\sim190\deg^2$ overlap of the LRG catalogue and the CS82 and CFHTLenS shear catalogues. We show that fast- and slow-growing LRGs have similar lensing profiles, which implies that they live in haloes of similar mass: $\log\left(M_{\rm halo}^{\rm fast}/h^{-1}\mathrm{M}_{\odot}\right) = 12.85^{+0.16}_{-0.26}$ and $\log\left(M_{\rm halo}^{\rm slow}/h^{-1}\mathrm{M}_{\odot}\right) =12.92^{+0.16}_{-0.22}$. This result, combined with the clustering difference, suggests the existence of galaxy assembly bias, although the effect is too subtle to be definitively proven given the errors on our current weak-lensing measurement. We show that this can soon be achieved with upcoming surveys like DES.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

Numerical simulations play a key important role in modern cosmology. Examples are plenty including the cosmic web - large scale structure of the distribution of galaxies in space - which was first observed in N-body simulations and later discovered in observations. The cuspy dark matter halo profiles, the overabundance of satellites, the Too-Big-Too-Fail problem are other examples of theoretical predictions that have a dramatic impact on recent developments in cosmology and galaxy formation. Large observational surveys such as e.g. SDSS, Euclid, and LSST are intimately connected with extensive cosmological simulations that provide statistical errors and tests for systematics. Accurate predictions for baryonic acoustic oscillations and redshift space distortions from high-resolution and large-volume cosmological simulations are required for interpretation of these large-scale galaxy/qso surveys. However, most of the results from extensive computer simulations, that would be greatly beneficial if publicly available, are still in hands of few research groups. Even when the simulation data is available, sharing vast amounts of data can be overwhelming. We argue that there is an effective and simple path to expand the data access and dissemination of numerous results from different cosmological models. Here we demonstrate that public access can be effectively provided with relatively modest resources. Among different results, we release for the astronomical community terabytes of raw data of th popular Bolshoi and MultiDark simulations. We also provide numerous results that are focused on mimicking observational data and galaxy surveys for major projects. Skies and Universes is a community effort: data are produced and shared by many research groups. We offer to other cosmologists and astronomers to host their data products in the skiesanduniverses.org space.

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available.

We study the galaxy clustering dependence on the $\left[\mathrm{O\scriptsize{II}}\right]$ emission line luminosity in the SDSS DR7 Main galaxy sample at mean redshift $z\sim0.1$. We select volume-limited samples of galaxies with different $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity thresholds and measure their projected, monopole and quadrupole two-point correlation functions. We model these observations using the 1$h^{-1}\rm{Gpc}$ MultiDark Planck cosmological simulation and generate light-cones with the SUrvey GenerAtoR algorithm. To interpret our results, we adopt a modified (Sub)Halo Abundance Matching scheme, accounting for the stellar mass incompleteness of the emission line galaxies. The satellite fraction constitutes an extra parameter in this model and allows to optimize the clustering fit on both small and intermediate scales (i.e. $r_p\lesssim 30h^{-1}\rm{Mpc})$, with no need of any velocity bias correction. We find that, in the local Universe, the $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity correlates with all the clustering statistics explored and with the galaxy bias. This latter quantity correlates more strongly with the SDSS $r$-band magnitude than $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity. In conclusion, we propose a straightforward method to produce reliable clustering models, entirely built on the simulation products, which provides robust predictions of the typical ELG host halo masses and satellite fraction values. The SDSS galaxy data, MultiDark mock catalogues and clustering results are made publicly available.

The one-point probability distribution function (PDF) of the matter density field in the universe is a fundamental property that plays an essential role in cosmology for estimates such as gravitational weak lensing, non-linear clustering, massive production of mock galaxy catalogs, and testing predictions of cosmological models. Here we make a comprehensive analysis of the dark matter PDF using a suite of 7000 N-body simulations that covers a wide range of numerical and cosmological parameters. We find that the PDF has a simple shape: it declines with density as a power-law P~rho**(-2), which is exponentially suppressed on both small and large densities. The proposed double-exponential approximation provides an accurate fit to all our N-body results for small filtering scales R< 5Mpc/h with rms density fluctuations sigma>1. In combination with the spherical infall model that works well for small fluctuations sigma<1, the PDF is now approximated with just few percent errors over the range of twelve orders of magnitude -- a remarkable example of precision cosmology. We find that at 5-10% level the PDF explicitly depends on redshift (at fixed sigma) and on cosmological density parameter Omega_m. We test different existing analytical approximations and find that the often used log-normal approximation is always 3-5 times less accurate than either the double-exponential approximation or the spherical infall model.

$N$-body cosmological simulations are an essential tool to understand the observed distribution of galaxies. We use the MultiDark simulation suite, run with the Planck cosmological parameters, to revisit the mass and velocity functions. At redshift $z=0$, the simulations cover four orders of magnitude in halo mass from $\sim10^{11}M_\odot$ with 8,783,874 distinct halos and 532,533 subhalos. The total volume used is $\sim$515 Gpc$^3$, more than 8 times larger than in previous studies. We measure and model the halo mass function, its covariance matrix w.r.t halo mass and the large scale halo bias. With the formalism of the excursion-set mass function, we explicit the tight interconnection between the covariance matrix, bias and halo mass function. We obtain a very accurate ($<2\%$ level) model of the distinct halo mass function. We also model the subhalo mass function and its relation to the distinct halo mass function. The set of models obtained provides a complete and precise framework for the description of halos in the concordance Planck cosmology. Finally, we provide precise analytical fits of the $V_{max}$ maximum velocity function up to redshift $z<2.3$ to push for the development of halo occupation distribution using $V_{max}$. The data and the analysis code are made publicly available in the \textit{Skies and Universes} database.

We analyze the spectra of 300,000 luminous red galaxies (LRGs) with stellar masses $M_* \gtrsim 10^{11} M_{\odot}$ from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). By studying their star-formation histories, we find two main evolutionary paths converging into the same quiescent galaxy population at $z\sim0.55$. Fast-growing LRGs assemble $80\%$ of their stellar mass very early on ($z\sim5$), whereas slow-growing LRGs reach the same evolutionary state at $z\sim1.5$. Further investigation reveals that their clustering properties on scales of $\sim$1-30 Mpc are, at a high level of significance, also different. Fast-growing LRGs are found to be more strongly clustered and reside in overall denser large-scale structure environments than slow-growing systems, for a given stellar-mass threshold. Our results imply a dependence of clustering on stellar-mass assembly history (naturally connected to the mass-formation history of the corresponding halos) for a homogeneous population of similar mass and color, which constitutes a strong observational evidence of galaxy assembly bias.

We report and provide fitting functions for the abundance of dark matter halos and subhalos as a function of mass, circular velocity, and redshift from the new Bolshoi-Planck and MultiDark-Planck $\Lambda$CDM cosmological simulations, based on the Planck cosmological parameters. We also report the halo mass accretion rates, which may be connected with galaxy star formation rates. We show that the higher cosmological matter density of the Planck parameters compared with the WMAP parameters leads to higher abundance of massive halos at high redshifts. We find that the median halo spin parameter $\lambda_{\rm B} = J(2M_{\rm vir}R_{\rm vir}V_{\rm vir})^{-1}$ is nearly independent of redshift, leading to predicted evolution of galaxy sizes that is consistent with observations, while the significant decrease with redshift in median $\lambda_{\rm P} = J|E|^{-1/2}G^{-1}M^{-5/2}$ predicts more decrease in galaxy sizes than is observed. Using the Tully-Fisher and Faber-Jackson relations between galaxy velocity and mass, we show that a simple model of how galaxy velocity is related to halo maximum circular velocity leads to increasing overprediction of cosmic stellar mass density as redshift increases beyond redshifts $z\sim1$, implying that such velocity-mass relations must change at redshifts $z>1$. By making a realistic model of how observed galaxy velocities are related to halo circular velocity, we show that recent optical and radio observations of the abundance of galaxies are in good agreement with our $\Lambda$CDM simulations. Our halo demographics are based on updated versions of the \rockstar\ and \ctrees\ codes, and this paper includes appendices explaining all of their outputs. This paper is an introduction to a series of related papers presenting other analyses of the Bolshoi-Planck and MultiDark-Planck simulations.

In current and future surveys, quasars play a key role. The new data will extend our knowledge of the Universe as it will be used to better constrain the cosmological model at redshift $z>1$ via baryon acoustic oscillation and redshift space distortion measurements. Here, we present the first clustering study of quasars observed by the extended Baryon Oscillation Spectroscopic Survey. We measure the clustering of $\sim 70,000$ quasars located in the redshift range $0.9<z<2.2$ that cover 1,168 deg$^2$. We model the clustering and produce high-fidelity quasar mock catalogues based on the BigMultiDark Planck simulation. Thus, we use a modified (Sub)Halo Abundance Matching model to account for the specificities of the halo population hosting quasars. We find that quasars are hosted by halos with masses $\sim10^{12.7}M_\odot$ and their bias evolves from 1.54 ($z=1.06$) to 3.15 ($z=1.98$). Using the current eBOSS data, we cannot distinguish between models with different fractions of satellites. The high-fidelity mock light-cones, including properties of halos hosting quasars, are made publicly available.

Upcoming and existing large-scale surveys of galaxies require accurate theoretical predictions of the dark matter clustering statistics for thousands of mock galaxy catalogs. We demonstrate that this goal can be achieve with our new Parallel Particle-Mesh (PM) Nbody code (PPM-GLAM) at a very low computational cost. We run about 15,000 simulations with ~2 billion particles that provide ~1% accuracy of the dark matter power spectra P(k) for wave-numbers up to k~ 1h/Mpc. Using this large data-set we study the power spectrum covariance matrix, the stepping stone for producing mock catalogs. In contrast to many previous analytical and numerical results, we find that the covariance matrix normalised to the power spectrum C(k,k')/P(k)P(k') has a complex structure of non-diagonal components. It has an upturn at small k, followed by a minimum at k=0.1-0.2h/Mpc. It also has a maximum at k=0.5-0.6h/Mpc. The normalised covariance matrix strongly evolves with redshift: C(k,k')~delta(t)^alpha P(k)P(k'), where delta is the linear growth factor and alpha ~ 1-1.25, which indicates that the covariance matrix depends on cosmological parameters. We also show that waves longer than 1Gpc have very little impact on the power spectrum and covariance matrix. This significantly reduces the computational costs and complexity of theoretical predictions: relatively small volume ~ (1Gpc)^3 simulations capture the necessary properties of dark matter clustering statistics. All the power spectra obtained from many thousands of our simulations are publicly available.

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.

Current and future large redshift surveys, as the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), will use emission-line galaxies (ELG) to probe cosmological models by mapping the large-scale structure of the Universe in the redshift range $0.6 < z < 1.7$. With current data, we explore the halo-galaxy connection by measuring three clustering properties of $g$-selected ELGs as matter tracers in the redshift range $0.6 < z < 1$: (i) the redshift-space two-point correlation function using spectroscopic redshifts from the BOSS ELG sample and VIPERS; (ii) the angular two-point correlation function on the footprint of the CFHT-LS; (iii) the galaxy-galaxy lensing signal around the ELGs using the CFHTLenS. We interpret these observations by mapping them onto the latest high-resolution MultiDark Planck N-body simulation, using a novel (Sub)Halo-Abundance Matching technique that accounts for the ELG incompleteness. ELGs at $z\sim0.8$ live in halos of $(1\pm 0.5)\times10^{12}\,h^{-1}$M$_{\odot}$ and 22.5$\pm2.5$% of them are satellites belonging to a larger halo. The halo occupation distribution of ELGs indicates that we are sampling the galaxies in which stars form in the most efficient way, according to their stellar-to-halo mass ratio.

We present a study of the clustering and halo occupation distribution of BOSS CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a Halo Abundance Matching (HAM) clustering model that assigns galaxies to dark matter halos selected from the large BigMultiDark $N$-body simulation of a flat $\Lambda$CDM Planck cosmology. We compare the observational data with the simulated ones on a light-cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fiber collisions are included in the model, which reproduces within 1-$\sigma$ errors the observed monopole of the 2-point correlation function at all relevant scales: from the smallest scales, 0.5 $h^{-1}$ Mpc, up to scales beyond the Baryonic Acoustic Oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to $k\sim1$ $h$ Mpc$^{-1}$), and the Three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar to halo mass relation for the CMASS sample measured using weak lensing in the CFHT Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1-$\sigma$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available.

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.

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 new measurements of the abundance of galaxies with a given circular velocity in the Local Volume: a region centered on the Milky Way Galaxy and extending to distance 10Mpc. The sample of 750 mostly dwarf galaxies provides a unique opportunity to study the abundance and properties of galaxies down to absolute magnitudes MB= -10, and virial masses Mvir= 1e9Msun. We find that the standard LCDM model gives remarkably accurate estimates for the velocity function of galaxies with circular velocities V>60kms and corresponding virial masses Mvir> 3e10Msun, but it badly fails by over-predicting 5 times the abundance of large dwarfs with velocities V= 30-50kms. The Warm Dark Matter models cannot explain the data either, regardless of mass of the WDM particle. Just as in previous observational studies, we find a shallow asymptotic slope dN/dlog V = V**alpha, alpha =-1 of the velocity function, which is inconsistent with the standard LCDM model that predicts the slope alpha =-3. Though reminiscent to the known overabundance of satellites problem, the overabundance of field galaxies is a much more difficult problem. For the LCDM model to survive, in the 10Mpc radius of the Milky Way there should be 1000 dark galaxies with virial mass Mvir= 1e10Msun, extremely low surface brightness and no detectable HI gas. So far none of this type of galaxies have been discovered.

We reproduce the galaxy clustering catalogue from the SDSS-III Baryon Oscillation Spectroscopic Survey Final Data Release (BOSS DR11&DR12) with high fidelity on all relevant scales in order to allow a robust analysis of baryon acoustic oscillations and redshift space distortions. We have generated (6,000) 12,288 MultiDark PATCHY BOSS (DR11) DR12 light-cones corresponding to an effective volume of $\sim192,000\,[h^{-1}\,{\rm Gpc}]^3$ (the largest ever simulated volume), including cosmic evolution in the redshift range from 0.15 to 0.75. The mocks have been calibrated using a reference galaxy catalogue based on the halo abundance matching modelling of the BOSS DR11&DR12 galaxy clustering data and on the data themselves. The production follows three steps. First, we apply the PATCHY code to generate a dark matter field and an object distribution including nonlinear stochastic galaxy bias. Secondly, we run the halo/stellar distribution reconstruction HADRON code to assign masses to the various objects. This step uses the mass distribution as a function of local density and non-local indicators (i.e., tidal field tensor eigenvalues and relative halo exclusion separation for massive objects) from the reference simulation applied to the corresponding PATCHY dark matter and galaxy distribution. Finally, we apply the SUGAR code to build the light cones. The resulting MultiDark PATCHY mock light cones reproduce the number density, selection function, survey geometry, and in general within 1 $\sigma$, for arbitrary stellar mass bins, the power spectrum up to $k=0.3\,h\,{\rm Mpc}^{-1}$, the two-point correlation functions down to a few Mpc scales, and the three-point statistics of the BOSS DR11&DR12 galaxy samples.

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.

Using 3D spectroscopy with a scanning Fabry-Perot interferometer, we study the ionized gas kinematics in 59 nearby dwarf galaxies. Combining our results with data from literature, we provide a global relation between the gas velocity dispersion (sigma) and the star formation rate (SFR) and H\alpha luminosity for galaxies in a very broad range of star formation rates SFR=0.001-300 Msun/yr. We find that the SFR-sigma relation for the combined sample of dwarf galaxies, star forming, local luminous, and ultra-luminous infrared galaxies can be fitted as sigma~ SFR^(5.3+-0.2). This implies that the slope of the L-sigma relation inferred from the sample of rotation supported disc galaxies (including mergers) is similar to the L-sigma relation of individual giant HII regions. We present arguments that the velocity dispersion of the ionized gas does not reflect the virial motions in the gravitational potential of dwarf galaxies, and instead is mainly determined by the energy injected into the interstellar medium by the ongoing star formation.

We study the circumgalactic medium (CGM) of a z=0.54 simulated dwarf galaxy using hydroART simulations. We present our analysis methods, which emulate observations, including objective absorption line detection, apparent optical depth (AOD) measurements, Voigt profile (VP) decomposition, and ionization modeling. By comparing the inferred CGM gas properties from the absorption lines directly to the gas selected by low ionization HI and MgII, and by higher ionization CIV and OVI absorption, we examine how well observational analysis methods recover the "true" properties of CGM gas. In this dwarf galaxy, low ionization gas arises in sub-kiloparsec "cloud" structures, but high ionization gas arises in multiple extended structures spread over 100 kpc; due to complex velocity fields, highly separated structures give rise to absorption at similar velocities. We show that AOD and VP analysis fails to accurately characterize the spatial, kinematic, and thermal conditions of high ionization gas. We find that HI absorption selected gas and OVI absorption gas arise in totally distinct physical gas structures, calling into question current observational techniques employed to infer metallicities and the total mass of "warm-hot" CGM gas. We present a method to determine whether CIV and OVI absorbing gas is photo or collisionally ionized and whether the assumption of ionization equilibrium is sound. As we discuss, these and additional findings have strong implications for how accurately currently employed observational absorption line methods recover the true gas properties, and ultimately, our ability to understand the CGM and its role in galaxy evolution.

Despite recent success in forming realistic present-day galaxies, simulations still form the bulk of their stars earlier than observations indicate. We investigate the process of stellar mass assembly in low-mass field galaxies, a dwarf and a typical spiral, focusing on the effects of radiation from young stellar clusters on the star formation (SF) histories. We implement a novel model of SF with a deterministic low efficiency per free-fall time, as observed in molecular clouds. Stellar feedback is based on observations of star-forming regions, and includes radiation pressure from massive stars, photoheating in H II regions, supernovae and stellar winds. We find that stellar radiation has a strong effect on the formation of low-mass galaxies, especially at z > 1, where it efficiently suppresses SF by dispersing cold and dense gas, preventing runaway growth of the stellar component. This behaviour is evident in a variety of observations but had so far eluded analytical and numerical models without radiation feedback. Compared to supernovae alone, radiation feedback reduces the SF rate by a factor of ~100 at z < 2, yielding rising SF histories which reproduce recent observations of Local Group dwarfs. Stellar radiation also produces bulgeless spiral galaxies and may be responsible for excess thickening of the stellar disc. The galaxies also feature rotation curves and baryon fractions in excellent agreement with current data. Lastly, the dwarf galaxy shows a very slow reduction of the central dark matter density caused by radiation feedback over the last ~7 Gyr of cosmic evolution.

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.

Any successful model of galaxy formation needs to explain the low rate of star formation in the small progenitors of today's galaxies. This inefficiency is necessary for reproducing the low stellar-to-virial mass fractions, suggested by current abundance matching models. A possible driver of this low efficiency is the radiation pressure exerted by ionizing photons from massive stars. The effect of radiation pressure in cosmological, zoom-in galaxy formation simulations is modeled as a non-thermal pressure that acts only in dense and optically thick star-forming regions. We also include photoionization and photoheating by massive stars. The full photoionization of hydrogen reduces the radiative cooling in the $10^{4-4.5}$ K regime. The main effect of radiation pressure is to regulate and limit the high values of gas density and the amount of gas available for star formation. This maintains a low star formation rate of $\sim 1 \ {\rm M_\odot} \ {\rm yr}^{-1}$ in halos with masses about $10^{11} \ {M_\odot}$ at $z\simeq3$. Infrared trapping and photoionization/photoheating processes are secondary effects in this mass range. The galaxies residing in these low-mass halos contain only $\sim0.6\%$ of the total virial mass in stars, roughly consistent with abundance matching. Radiative feedback maintains an extended galaxy with a rising circular velocity profile.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

Dwarf spheroidal galaxies (dSphs) are extremely gas-poor, dark matter-dominated galaxies, which make them ideal to test the predictions of the cold dark matter (CDM) model. We argue that the removal of the baryonic component from gas-rich dwarf irregular galaxies, the progenitors of dSphs, can substantially reduce their central density. Thus, it may play an important role in alleviating one of the problems of the CDM model related to the structure of relatively massive satellite galaxies of the Milky Way (MW). Traditionally, collisionless cosmological N-body simulations are used when confronting theoretical predictions with observations. However, these simulations assume that the baryon fraction everywhere in the Universe is equal to the cosmic mean, which can be incorrect for dSphs. We find that the combination of (i) the lower baryon fraction in dSphs compared to the cosmic mean and (ii) the concentration of baryons in the inner part of the MW halo can go a long way towards explaining the observed circular velocity profiles of dSphs. We find that the blowing away of baryons by ram pressure, when the dwarfs fall into larger galaxies, decreases the circular velocity profile of the satellite and reduces the density in the central ~200-500 pc by a factor of (1 - f_b)^4 ~ 0.5, where f_b is the cosmological fraction of baryons. Additionally, the enhanced baryonic mass in the central regions of the parent galaxy generates tidal forces, which are larger than those experienced by subhaloes in traditional N-body simulations. Increased tidal forces substantially alter circular velocity profiles for satellites with pericentres less than 50 kpc. We show that these two effects are strong enough to bring the predictions of subhaloes from CDM simulations into agreement with the observed structure of MW dSphs, regardless of the details of the baryonic processes.

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.

We find that infalling dark matter halos (i.e., the progenitors of satellite halos) begin losing mass well outside the virial radius of their eventual host halos. The peak mass occurs at a range of clustercentric distances, with median and 68th percentile range of 1.8 +2.3/-1.0 R_(vir,host) for progenitors of z=0 satellites. The peak circular velocity for infalling halos occurs at significantly larger distances (3.7 +3.3/-2.2 R_(vir,host) at z=0). This difference arises because different physical processes set peak circular velocity (typically, ~1:5 and larger mergers which cause transient circular velocity spikes) and peak mass (typically, smooth accretion) for infalling halos. We find that infalling halos also stop having significant mergers well before they enter the virial radius of their eventual hosts. Mergers larger than a 1:40 ratio in halo mass end for infalling halos at similar clustercentric distances (~ 1.9 R_(vir,host)) as the end of overall mass accretion. However, mergers larger than 1:3 typically end for infalling halos at more than 4 virial radial away from their eventual hosts. This limits the ability of mergers to affect quenching and morphology changes in clusters. We also note that the transient spikes which set peak circular velocity may lead to issues with abundance matching on that parameter, including unphysical galaxy stellar mass growth profiles near clusters; we propose a simple observational test to check if a better halo proxy for galaxy stellar mass exists.

We derive constraints on parameters of generic dark matter candidates by comparing theoretical predictions with the gamma-ray emission observed by the Fermi-LAT from the region around the Galactic Center. Our analysis is conservative since it simply requires that the expected dark matter signal does not exceed the observed emission. The constraints obtained in the likely case that the collapse of baryons to the Galactic Center is accompanied by the contraction of the dark matter are strong. In particular, we find that for bb and \tau+\tau- or W+W- dark matter annihilation channels, the upper limits on the annihilation cross section imply that the thermal cross section is excluded for a Weakly Interacting Massive Particle (WIMP) mass smaller than about 700 and 500 GeV, respectively. For the \mu+ \mu- channel, where the effect of the inverse Compton scattering is important, depending on models of the Galactic magnetic field the exclusion of the thermal cross-section is for a WIMP mass smaller than about 150 to 400 GeV. The upper limits on the annihilation cross section of dark matter particles obtained are two orders of magnitude stronger than without contraction. In the latter case our results are compatible with the upper limits from the Galactic halo analysis reported by the Fermi-LAT collaboration for the case in which the same conservative approach without modeling of the astrophysical background is employed.

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.

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.

We present a new algorithm for generating merger trees and halo catalogs which explicitly ensures consistency of halo properties (mass, position, and velocity) across timesteps. Our algorithm has demonstrated the ability to improve both the completeness (through detecting and inserting otherwise missing halos) and purity (through detecting and removing spurious objects) of both merger trees and halo catalogs. In addition, our method is able to robustly measure the self-consistency of halo finders; it is the first to directly measure the uncertainties in halo positions, halo velocities, and the halo mass function for a given halo finder based on consistency between snapshots in cosmological simulations. We use this algorithm to generate merger trees for two large simulations (Bolshoi and Consuelo) and evaluate two halo finders (ROCKSTAR and BDM). We find that both the ROCKSTAR and BDM halo finders track halos extremely well; in both, the number of halos which do not have physically consistent progenitors is at the 1-2% level across all halo masses. Our code is publicly available at http://code.google.com/p/consistent-trees . Our trees and catalogs are publicly available at http://hipacc.ucsc.edu/Bolshoi/ .

The alignment of DM halos and the surrounding large scale structure (LSS) is examined in the context of the cosmic web. Halo spin, shape and the orbital angular momentum of subhaloes is investigated relative to the LSS using the eigenvectors of the velocity shear tensor evaluated on a grid with a scale of 1 Mpc/h, deep within the non-linear regime. Knots, filaments, sheets and voids are associated with regions that are collapsing along 3, 2, 1 or 0 principal directions simultaneously. Each halo is tagged with a web classification (i.e. knot halo, filament halo, etc) according to the nature of the collapse at the halo's position. The full distribution of shear eigenvalues is found to be substantially different from that tagged to haloes, indicating that the observed velocity shear is significantly biased. We find that larger mass haloes live in regions where the shear is more isotropic, namely the expansion or collapse is more spherical. A correlation is found between the halo's shape and the eigenvectors of the shear tensor, with the longest (shortest) axis of the halo's shape being aligned with the slowest (fastest) collapsing eigenvector. This correlation is web independent, suggesting that the velocity shear is a fundamental tracer of the halo alignment. A similar result is found for the alignment of halo spin with the cosmic web. It has been shown that high mass haloes exhibit a spin flip with respect to the LSS: we find the mass at which this spin flip occurs is web dependent and not universal as suggested previously. Although weaker than haloes, subhalo orbits too exhibit an alignment with the LSS, providing a possible insight into the highly correlated co-rotation of the Milky Way's satellite system. The present study suggests that the velocity shear tensor constitutes the natural framework for studying the directional properties of the non-linear LSS and of halos and galaxies.

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

We calculate the probability that a Milky-Way-like halo in the standard cosmological model has the observed number of Magellanic Clouds (MCs). The statistics of the number of MCs in the LCDM model are in good agreement with observations of a large sample of SDSS galaxies. Under the sub-halo abundance matching assumption of a relationship with small scatter between galaxy r-band luminosities and halo internal velocities v_max, we make detailed comparisons to similar measurements using SDSS DR7 data by Liu et al. (2010). Models and observational data give very similar probabilities for having zero, one, and two MC-like satellites. In both cases, Milky Way-luminosity hosts have just a \sim 10% chance of hosting two satellites similar to the Magellanic Clouds. In addition, we present a prediction for the probability for a host galaxy to have Nsats satellite galaxies as a function of the magnitudes of both the host and satellite. This probability and its scaling with host properties is significantly different from that of mass-selected objects because of scatter in the mass- luminosity relation and because of variations in the star formation efficiency with halo mass.

We present a new measurement of the mass of the Milky Way (MW) based on observed properties of its largest satellite galaxies, the Magellanic Clouds (MCs), and an assumed prior of a {\Lambda}CDM universe. The large, high-resolution Bolshoi cosmological simulation of this universe provides a means to statistically sample the dynamical properties of bright satellite galaxies in a large population of dark matter halos. The observed properties of the MCs, including their circular velocity, distance from the center of the MW, and velocity within the MW halo, are used to evaluate the likelihood that a given halo would have each or all of these properties; the posterior PDF for any property of the MW system can thus be constructed. This method provides a constraint on the MW virial mass, 1.2 +0.7 -0.3(stat.) +0.3 -0.4 (sys.) x 10^12 M\odot (68% confidence), which is consistent with recent determinations that involve very different assumptions. In addition, we calculate the posterior PDF for the density profile of the MW and its satellite accretion history. Although typical satellites of 10^12 M\odot halos are accreted over a wide range of epochs over the last 10 Gyr, we find a \sim72% probability that the Magellanic Clouds were accreted within the last Gyr, and a 50% probability that they were accreted together.

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.

It has long been regarded as difficult for a cosmological model to account simultaneously for the galaxy luminosity, mass, and velocity distributions. We revisit this issue using a modern compilation of observational data along with the best available large-scale cosmological simulation of dark matter. We find that the standard cosmological model, used in conjunction with halo abundance matching (HAM) and simple dynamical corrections, fits all basic statistics of galaxies with circular velocities Vcirc > 80 km/s. Our observational constraint is the luminosity-velocity relation which allows all types of galaxies to be included. We have compiled data for a variety of galaxies ranging from dwarf irregulars to giant ellipticals. The data present a clear monotonic luminosity-velocity relation from 50 km/s to 500 km/s, with a bend below 80 km/s and a systematic offset between late- and early-type galaxies. For comparison to theory, we employ our LCDM "Bolshoi" simulation of dark matter, which has unprecedented mass and force resolution. We use halo abundance matching to assign rank-ordered galaxy luminosities to the dark matter halos. The resulting predictions for the luminosity-velocity relation are in excellent agreement with the available data on both early-type and late-type galaxies for the luminosity range from Mr = -14-22. We also compare our predictions for the "cold" baryon mass (i.e., stars and cold gas) of galaxies as a function of circular velocity with the available observations, again finding a very good agreement. The predicted circular velocity function is in agreement with the galaxy velocity function for 80-400 km/s. However, we find that the dark matter halos with Vcirc < 80 km/s are much more abundant than observed galaxies with the same Vcirc . We find that the two-point correlation function of galaxies in our model matches very well the results from the SDSS.

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.

The condensation of gas and stars in the inner regions of dark matter halos leads to a more concentrated dark matter distribution. While this effect is based on simple gravitational physics, the question of its validity in hierarchical galaxy formation has led to an active debate in the literature. We use a collection of several state-of-the-art cosmological hydrodynamic simulations to study the halo contraction effect in systems ranging from dwarf galaxies to clusters of galaxies, at high and low redshift. The simulations are run by different groups with different codes and include hierarchical merging, gas cooling, star formation, and stellar feedback. We show that in all our cases the inner dark matter density increases relative to the matching simulation without baryon dissipation, at least by a factor of several. The strength of the contraction effect varies from system to system and cannot be reduced to a simple prescription. We present a revised analytical model that describes the contracted mass profile to an rms accuracy of about 10%. The model can be used to effectively bracket the response of the dark matter halo to baryon dissipation. The halo contraction effect is real and must be included in modeling of the mass distribution of galaxies and galaxy clusters.

We present the first results from the new Bolshoi N-body cosmological LCDM simulation that uses cosmological parameters favored by current observations. The Bolshoi simulation was done in a volume 250Mpc on a side using 8billion particles with mass and force resolution adequate to follow subhalos down to a completeness limit of Vcirc=50km/ s circular velocity. Using excellent statistics of halos and subhalos (10M at every moment and 50M over the whole history) we present accurate approximations for statistics such as the halo mass function, the concentrations for distinct halos and subhalos, abundance of halos as function of their circular velocity, the abundance and the spatial distribution of subhalos. We find that at high redshifts the concentration falls to a minimum of about 3.8 and then rises slightly for higher values of halo mass. We find that while the Sheth-Tormen approximation for the mass function of halos found by spherical overdensity is accurate at low redshifts, it over-predicts the abundance of halos by nearly an order of magnitude by z=10. We find that the number of subhalos scales with the circular velocity of the host halo as Vhost**0.5, and that subhalos have nearly the same radial distribution as dark matter particles at radii 0.3-2 times the host halo virial radius. The subhalo velocity function n(>V) behaves as V**(-3). We give normalization of this relation for different masses and redshifts. Finally, we use an abundance-matching procedure to assign r-band luminosities to dark matter halos as a function of halo Vcirc, and find that the luminosity-velocity relation is in remarkably good agreement with the observed Tully-Fisher relation for galaxies in the range 50-200km/s.

We make detailed theoretical predictions for the assembly properties of the Local Group (LG) in the standard LambdaCDM cosmological model. We use three cosmological N-body dark matter simulations from the CLUES project, which are designed to reproduce the main dynamical features of the matter distribution down to the scale of a few Mpc around the LG. Additionally, we use the results of an unconstrained simulation with a sixty times larger volume to calibrate the influence of cosmic variance. We characterize the Mass Aggregation History (MAH) for each halo by three characteristic times, the formation, assembly and last major merger times. A major merger is defined by a minimal mass ratio of 10:1. We find that the three LGs share a similar MAH with formation and last major merger epochs placed on average \approx 10 - 12 Gyr ago. Between 12% and 17% of the halos in the mass range 5 x 10^11 Msol/h < M_h < 5 x 10^12 Msol/h have a similar MAH. In a set of pairs of halos within the same mass range, a fraction of 1% to 3% share similar formation properties as both halos in the simulated LG. An unsolved question posed by our results is the dynamical origin of the MAH of the LGs. The isolation criteria commonly used to define LG-like halos in unconstrained simulations do not narrow down the halo population into a set with quiet MAHs, nor does a further constraint to reside in a low density environment. The quiet MAH of the LGs provides a favorable environment for the formation of disk galaxies like the Milky Way and M31. The timing for the beginning of the last major merger in the Milky Way dark matter halo matches with the gas rich merger origin for the thick component in the galactic disk. Our results support the view that the specific large and mid scale environment around the Local Group play a critical role in shaping its MAH and hence its baryonic structure at present.

We study the concentration of dark matter halos and its evolution in N-body simulations of the standard LCDM cosmology. The results presented in this paper are based on 4 large N-body simulations with about 10 billion particles each: the Millennium-I and II, Bolshoi, and MultiDark simulations. The MultiDark (or BigBolshoi) simulation is introduced in this paper. This suite of simulations with high mass resolution over a large volume allows us to compute with unprecedented accuracy the concentration over a large range of scales (about six orders of magnitude in mass), which constitutes the state-of-the-art of our current knowledge on this basic property of dark matter halos in the LCDM cosmology. We find that there is consistency among the different simulation data sets. We confirm a novel feature for halo concentrations at high redshifts: a flattening and upturn with increasing mass. The concentration c(M,z) as a function of mass and the redshift and for different cosmological parameters shows a remarkably complex pattern. However, when expressed in terms of the linear rms fluctuation of the density field sigma(M,z), the halo concentration c(sigma) shows a nearly-universal simple U-shaped behaviour with a minimum at a well defined scale at sigma=0.71. Yet, some small dependences with redshift and cosmology still remain. At the high-mass end (sigma < 1) the median halo kinematic profiles show large signatures of infall and highly radial orbits. This c-sigma(M,z) relation can be accurately parametrized and provides an analytical model for the dependence of concentration on halo mass. When applied to galaxy clusters, our estimates of concentrations are substantially larger -- by a factor up to 1.5 -- than previous results from smaller simulations, and are in much better agreement with results of observations. (abridged)

[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).

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.

We obtained ESI/Keck rotation curves of 10 MgII absorption selected galaxies (0.3 < z < 1.0) for which we have WFPC-2/HST images and high resolution HIRES/Keck and UVES/VLT quasar spectra of the MgII absorption profiles. We perform a kinematic comparison of these galaxies and their associated halo MgII absorption. For all 10 galaxies, the majority of the absorption velocities lie in the range of the observed galaxy rotation velocities. In 7/10 cases, the absorption velocities reside fully to one side of the galaxy systemic velocity and usually align with one arm of the rotation curve. In all cases, a constant rotating thick-disk model poorly reproduces the full spread of observed MgII absorption velocities when reasonably realistic parameters are employed. In 2/10 cases, the galaxy kinematics, star formation surface densities, and absorption kinematics have a resemblance to those of high redshift galaxies showing strong outflows. We find that MgII absorption velocity spread and optical depth distribution may be dependent on galaxy inclination. To further aid in the spatial-kinematic relationships of the data, we apply quasar absorption line techniques to a galaxy (v_c=180 km/s) embedded in LCDM simulations. In the simulations, MgII absorption selects metal enriched "halo" gas out to roughly 100 kpc from the galaxy, tidal streams, filaments, and small satellite galaxies. Within the limitations inherent in the simulations, the majority of the simulated MgII absorption arises in the filaments and tidal streams and is infalling towards the galaxy with velocities between -200 < v_r < -180 km/s. The MgII absorption velocity offset distribution (relative to the simulated galaxy) spans ~200 km/s with the lowest frequency of detecting MgII at the galaxy systematic velocity.

We study the luminosity function, the peculiar velocities, and the sizes of voids in the Local Volume (LV) in observational samples of galaxies which contain galaxies down to MB=-10-12 within the distance 4-8Mpc. Comparing the results with the predictions of the standard cosmological LCDM model, we find that the theory faces a sever problem: it predicts a factor of ten more dwarf haloes as compared with the observed number of dwarf galaxies. In the LV we identify voids with sizes ranging from 1 to 4.5Mpc and compare the observational distribution of void sizes with the voids in very high resolution simulations of the LCDM model with WMAP1 and WMAP3 parameters. The theoretical void function matches the observations remarkably well only if we use haloes with circular velocities Vc larger than 40-45 km/s (Mvir=(1-2)e10msun) for models with with sigma8 = 0.9 and Vc>35km/s (Mvir=(6-8)e9msun) for sigma8 = 0.75. We exclude the possibility that in the LCDM model haloes with circular velocities < 35km/s can host galaxies as bright as MB=-12: there are too many small haloes in the LCDM model resulting in voids being too small as compared with the observations. The problem is that many of the observed dwarf galaxies have HI rotational velocities below 25km/s that strictly contradicts the LCDM predictions. Thus, the LCDM model faces the same overabundance problem, which it had with the number of satellites in the LG. We also estimate the rms deviations from the Hubble flow sigmaH and find that in most of our model LV-candidates the rms peculiar velocities are consistent with observational values sigmaH=50km/s for distances less than 3Mpc and sigmaH=80km/s for distances less than 8Mpc. At distances 4 (8) Mpc, the observed overdensities of galaxies are 3.5-5.5 (1.3-1.6) -- significantly larger than typically assumed.

We develop a new realistic prescription for modeling the stellar feedback, which minimizes any ad hoc assumptions about sub-grid physics. We start with developing high resolution models of the ISM and formulate the conditions required for its realistic functionality: formation of multi-phase medium with hot chimneys, super-bubbles, cold molecular phase, and very slow consumption of gas. Another important ingredient is the runaway stars. They greatly facilitate the feedback. Once these effects are implemented into cosmological simulations, we do not have the overcooling problem and the angular momentum problem (resulting in a too massive bulge) is also reduced substantially: the rotation curves are nearly flat. Just as it is often observed in QSO absorption lines, our models produce substantial outflows from forming and active galaxies. At high redshifts we routinely find gas with few hundred km/s and occasionally 1000-2000 km/s. The density profile of dark matter is still consistent with a cuspy profile. The simulations reproduce this picture only if the resolution is very high: better than 50 pc, which is much better than the typical resolution in previous cosmological simulations. Our simulations of galaxy formation reach the resolution of 35 pc. (Abridged)

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.

The motion of satellite galaxies around normal galaxies at distances 50-500 kpc provides a sensitive test for the theories. We study the surface density and the velocities of satellites around isolated galaxies in the Sloan Digital Sky Survey. We find that the surface number-density of satellites declines with the projected distance as a power law with the slope -1.5-2. The rms velocities gradually decline: observations exclude constant velocities at a 10 sigma level. We show that observational data strongly favor the standard model: all three major statistics of satellites - the number-density profile, the line-of-sight velocity dispersion, and the distribution function of the velocities -- agree remarkably well with the predictions of the standard cosmological model. Thus, that the success of the standard model extends to scales (50-500) kpc, much lower than what was previously considered. MOND fails on these scales for models which assume any single power-law number-density profile of satellites and any constant velocity anisotropy by predicting nearly constant rms velocities of satellites. Satellite data can be fit by fine-tuned models, which require (1) specific non-power-law density profile, (2) very radial orbits at large distances (velocity anisotropy beta =0.6-0.7 at 200-300 kpc), and (3) 2-2.5 times more stellar mass than what is found in the galaxies. The external gravity force - a necessary component for MOND -- makes the situation even worse. We argue that a combination of satellite data and observational constraints on stellar masses make these models very problematic.

We study dynamics of bars in models of disk galaxies embeded in realistic dark matter halos. We find that disk thickness plays an important, if not dominant, role in the evolution and structure of the bars. We also make extensive numerical tests of different N-body codes used to study bar dynamics. Models with thick disks typically used in this type of modeling (height-to-length ratio hz/Rd=0.2) produce slowly rotating, and very long, bars. In contrast, more realistic thin disks with the same parameters as in our Galaxy (hz/Rd= 0.1) produce bars with normal length Rbar approx R_d, which rotate quickly with the ratio of the corotation radius to the bar radius 1.2-1.4 compatible with observations. Bars in these models do not show a tendency to slow down, and may lose as little as 2-3 percent of their angular momentum due to dynamical friction with the dark matter over cosmological time. We attribute the differences between the models to a combined effect of high phase-space density and smaller Jeans mass in the thin disk models, which result in the formation of a dense central bulge. Special attention is paid to numerical effects such as the accuracy of orbital integration, force and mass resolution. Using three N-body codes -- Gadget, ART, and Pkdgrav -- we find that numerical effects are very important and, if not carefully treated, may produce incorrect and misleading results. Once the simulations are performed with sufficiently small time-steps and with adequate force and mass resolution, all the codes produce nearly the same results: we do not find any systematic deviations between the results obtained with TREE codes (Gadget and Pkdgrav) and with the Adaptive-Mesh-Refinement (ART) code.

We present an approximation for the average density profile of dark matter haloes in the LambdaCDM cosmological model, which is accurate to within 10--15% even for large radial distances from 0.05R_{vir} up to 10R_{vir} for halo masses ranging from 10^{11.5} to 10^{15.0} h^{-1}M_{Sun}. We propose a modified form of the Navarro, Frenk & White (NFW) approximation: \rho(r)=\rho_{NFW}(r)+ A (r/R_{vir})^{-1}+B (1+r/R_{vir})^{-1}. This generalized expression, which is applicable to the external regions of dark matter haloes, only very slightly affects the density in the inner regions of haloes. The strong correlation among the different parameters in the model allows us to describe the profile in terms of just one parameter: the virial mass. We integrate our density profile to derive the enclosed mass in a sphere of a given radius and compare it with the NFW results. We find that the NFW underestimates the enclosed mass by more than 50% at 10R_{vir}, whereas our model reproduces the results from numerical simulations to within 2% accuracy even at this distance. We also use this new approximation to study the weak gravitational lensing and to obtain an analytic expression for the tangential shear. This allows us to quantify the contribution to the shear from the outer regions of the density profile. For the first time we calculate the difference between the tangential shear calculated via the NFW profile and the corresponding result when the external regions of haloes in cosmological simulations are taken into account. We find a 4% difference for all the mass ranges under study.

(Abridged) Virial mass is used as an estimator for the mass of a dark matter halo. However, the commonly used constant overdensity criterion does not reflect the dynamical structure of haloes. Here we analyze dark matter cosmological simulations in order to obtain properties of haloes of different masses focusing on the size of the region with zero mean radial velocity. Dark matter inside this region is stationary, and thus the mass of this region is a much better approximation for the virial mass. We call this mass the static mass to distinguish from the commonly used constant overdensity mass. We also study the relation of this static mass with the traditional virial mass, and we find that the matter inside galaxy-size haloes is underestimated by the virial mass by nearly a factor of two. At redshift zero the virial mass is close to the static mass for cluster-size haloes. The same pattern - large haloes having M_vir > M_static - exists at all redshifts, but the transition mass M_0 = M_vir = M_static decreases dramatically with increasing redshift. When rescaled to the same M_0 haloes clearly demonstrate a self-similar behaviour, which in a statistical sense gives a relation between the static and virial mass. To our surprise we find that the abundance of haloes with a given static mass, i.e. the static mass function, is very accurately fitted by the Press & Schechter approximation at z=0, but this approximation breaks at higher redshifts. Instead, the virial mass function is well fitted as usual by the Sheth & Tormen approximation. We find an explanation why the static radius can be 2-3 times larger as compared with the constant overdensity estimate. Applying the non-stationary Jeans equation we find that the role of the pressure gradients is significantly larger for small haloes.

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).

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.

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 have explored the galaxy disk/extended halo gas kinematic relationship using rotation curves (Keck/ESI) of ten intermediate redshift galaxies which were selected by MgII halo gas absorption observed in quasar spectra. Previous results of six edge-on galaxies, probed along their major axis, suggest that observed halo gas velocities are consistent with extended disk-like halo rotation at galactocentric distances of 25-72 kpc. Using our new sample, we demonstrate that the gas velocities are by and large not consistent with being directly coupled to the galaxy kinematics. Thus, mechanisms other than co-rotation dynamics (i.e., gas inflow, feedback, galaxy-galaxy interactions, etc.) must be invoked to account for the overall observed kinematics of the halo gas. In order to better understand the dynamic interaction of the galaxy/halo/cosmic web environment, we performed similar mock observations of galaxies and gaseous halos in Lambda-CDM cosmological simulations. We discuss an example case of a z=0.92 galaxy with various orientations probing halo gas at a range of positions. The gas dynamics inferred using simulated quasar absorption lines are consistent with observational data.

Although supernova explosions and stellar winds happens at scales bellow 100 pc, they affect the interstellar medium(ISM) and galaxy formation. We use cosmological N-body+Hydrodynamics simulations of galaxy formation, as well as simulations of the ISM to study the effect of stellar feedback on galactic scales. Stellar feedback maintains gas with temperatures above a million degrees. This gas fills bubbles, super-bubbles and chimneys. Our model of feedback, in which 10%-30% of the feedback energy is coming from runaway stars, reproduces this hot gas only if the resolution is better than 50 pc. This is 10 times better than the typical resolution in cosmological simulations of galaxy formation. Only with this resolution, the effect of stellar feedback in galaxy formation is resolved without any assumption about sub-resolution physics. Stellar feedback can regulate the formation of bulges and can shape the inner parts of the rotation curve.

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.

Current explanation of the overabundance of dark matter subhalos in the Local Group (LG) indicates that there maybe a limit on mass of a halo, which can host a galaxy. This idea can be tested using voids in the distribution of galaxies: at some level small voids should not contain any (even dwarf) galaxies. We use observational samples complete to M_B=-12 with distances less than 8 Mpc to construct the void function (VF): the distribution of sizes of voids empty of any galaxies. There are ~ 30 voids with sizes ranging from 1 to 5 Mpc. We also study the distribution of dark matter halos in very high resolution simulations of the LCDM model. The theoretical VF matches the observations remarkably well only if we use halos with circular velocities larger than 45 +/- 10 km/s. This agrees with the Local Group predictions. Small voids look quite similar to heir giant cousins: the density has a minimum at the center of a void and it increases as we get closer to the border. Thus, both the Local Group data and the nearby voids indicate that isolated halos below 45 +/- 10 km/s must not host galaxies and that small (few Mpc) voids are truly dark.

Current explanation of the overabundance of dark matter subhalos in the Local Group (LG) indicates that there maybe a limit on mass of a halo, which can host a galaxy. This idea can be tested using voids in the distribution of galaxies: at some level small voids should not contain any (even dwarf) galaxies. We use observational samples complete to M_B = -12 with distances less than 8 Mpc to construct the void function (VF): the distribution of sizes of voids empty of any galaxies. There are ~30 voids with sizes ranging from 1 to 5 Mpc. We then study the distribution of dark matter halos in very high resolution simulations of the LCDM model. The theoretical VF matches the observations remarkably well only if we use halos with circular velocities larger than 45 +/- 10 km/s. This agrees with the Local Group predictions. There are smaller halos in the voids, but they should not produce any luminous matter. Small voids look quite similar to their giant cousins: the density has a minimum at the center of a void and it increases as we get closer to the border. Small nonluminous halos inside the void form a web of tiny filaments. Thus, both the Local Group data and the nearby voids indicate that isolated halos below 45 +/- 10 km/s must not host galaxies and that small (few Mpc) voids are truly dark.

The inner parts of many spiral galaxies are dominated by bars. These are strong non-axisymmetric features which significantly affect orbits of stars and dark matter particles. One of the main effects is the dynamical resonances between galactic material and the bar. We detect and characterize these resonances in N-body models of barred galaxies by measuring angular and radial frequencies of individual orbits. We found narrow peaks in the distribution of orbital frequencies with each peak corresponding to a specific resonance. We found five different resonances in the stellar disk and two in the dark matter. The corotation resonance and the inner and outer Lindblad resonances are the most populated. The spatial distributions of particles near resonances are wide. For example, the inner Lindblad resonance is not localized at a given radius. Particles near this resonance are mainly distributed along the bar and span a wide range of radii. On the other hand, particles near the corotation resonance are distributed in two broad areas around the two stable Lagrange points. The distribution resembles a wide ring at the corotation radius. Resonances capture disk and halo material in near-resonant orbits. Our analysis of orbits in both N-body simulations and in simple analytical models indicates that resonances tend to prevent the dynamical evolution of this trapped material. Only if the bar evolves as a whole, resonances drift through the phase space. In this case particles anchored near resonant orbits track the resonance shift and evolve. The criteria to ensure a correct resonant behavior discussed in Weinberg & Katz (2007a) can be achieved with few millions particles because the regions of trapped orbits near resonances are large and evolving.

We present a statistical analysis of voids in the 2dF galaxy redshift survey (2dFGRS). In order to detect the voids, we have developed two robust algorithms. We define voids as non-overlapping maximal spheres empty of halos or galaxies with mass or luminosity above a given one. We search for voids in cosmological $N$-Body simulations to test the performance of our void finders. We obtain and analyze the void statistics for several volume-limited samples for the North Galactic Strip (NGP) and the South Galactic Strip (SGP) constructed from the 2dFGRS full data release. We find that the results obtained from the NGP and the SGP are statistically compatible. From the results of several statistical tests we conclude that voids are essentially uncorrelated, with at most a mild anticorrelation and that there is a dependence of the void number density on redshift at least at the 99.5% confidence level. We develop a technique to correct the distortion caused by the fact that we use the redshift as the radial coordinate. We calibrate this technique with mock catalogues and find that the correction might be of some relevance to carry out accurate inferences from void statistics. We study the statistics of the galaxies inside nine nearby voids. We find that galaxies in voids are not randomly distributed: they form structures like filaments. We also obtain the galaxy number density profile in voids. This profile follow a similar but steeper trend to that follow by halos in voids.

We investigate the orientation of the axes and angular momentum of dark matter halos with respect to their neighboring voids using high resolution N-body cosmological simulations. We find that the minor axis of halos tends to be aligned along the line joining the halo with the center of the void, and that the major axis is perpendicular to this line. However, we find that the angular momentum of halos does not have any particular orientation. These results may provide information about the mechanisms whereby the large-scale structure of the Universe affects galaxy formation and cast light upon the issue of the orientation of galaxy disks with respect to their host halos.

Two well studied dwarf galaxies -- NGC 3109 and NGC 6822 -- present some of the strongest observational support for a flat core at the center of galactic dark matter (DM) halos. We use detailed cosmologically motivated numerical models to investigate the systematics and the accuracy of recovering parameters of the galaxies. Some of our models match the observed structure of the two galaxies remarkably well. Our analysis shows that the rotation curves of these two galaxies are instead quite compatible with their DM halos having steep cuspy density profiles. The rotation curves in our models are measured using standard observational techniques. The models reproduce the rotation curves of both galaxies, the disk surface brightness profiles as well as the profile of isophotal ellipticity and position angle. The models are centrally dominated by baryons; however, the dark matter component is globally dominant. The simulated disk mass is marginally consistent with a stellar mass-to-light ratio in agreement with the observed colors. We show that non-circular motions combined with gas pressure support and projection effects results in a large underestimation of the circular velocity in the central $\sim 1$ kpc region, creating the illusion of a constant density core. Although the systematic effects mentioned above are stronger in barred systems, they are also present in axisymetric disks. Our results strongly suggest that there is no contradiction between the observed rotation curves in dwarf galaxies and the cuspy central dark matter density profiles predicted by Cold Dark Matter models.

We measure the evolution in the virial mass-to-light ratio (M_{200}/L_B) and virial-to-stellar mass ratio (M_{200}/M_\ast) for isolated ~ L* galaxies between z~1 and z~0 by combining data from the DEEP2 Galaxy Redshift Survey and the Sloan Digital Sky Survey. Utilizing the motions of satellite galaxies around isolated galaxies, we measure line-of-sight velocity dispersions and derive dark matter halo virial masses for these host galaxies. At both epochs the velocity dispersion of satellites correlates with host galaxy stellar mass, \sigma\propto M_\ast^{0.4+/-0.1}, while the relation between satellite velocity dispersion and host galaxy B-band luminosity may grow somewhat shallower from \sigma\propto L_B^{0.6+/-0.1} at z~1 to \sigma\propto L_B^{0.4+/-0.1} at z~0. The evolution in M_200/M_\ast from z~1 to z~0 displays a bimodality insofar as host galaxies with stellar mass below M_\ast ~10^{11} M_Sun/h maintain a constant ratio (the intrinsic increase is constrained to a factor of 1.1+/-0.7) while host galaxies above M_\ast ~10^{11} M_Sun/h experience a factor of 4+/-3 increase in their virial-to-stellar mass ratio. This result can be easily understood if galaxies below this stellar mass scale continue to form stars while star formation in galaxies above this scale is quenched and the dark matter halos of galaxies both above and below this scale grow in accordance with LCDM cosmological simulations. Host galaxies that are red in U-B color have larger satellite dispersions and hence reside on average in more massive halos than blue galaxies at both z~1 and z~0. The redshift and host galaxy stellar mass dependence of M_200/M_\ast agrees qualitatively with the Millennium Run semi-analytic model of galaxy formation. (ABRIDGED)

We present a comparison of the properties of galaxies in the most underdense regions of the Universe, where the galaxy number density is less than 10% of the mean density, with galaxies from more typical regions. We have compiled a sample of galaxies in 46 large nearby voids that were identified using the Sloan Digital Sky Survey DR4, which provides the largest coverage of the sky. We study the u-r color distribution, morphology, specific star formation rate, and radial number density profiles for a total of 495 galaxies fainter than M_r=-20.4 +5logh located inside the voids and compare these properties with a control sample of field galaxies. We show that there is an excess of blue galaxies inside the voids. However, inspecting the properties of blue and red galaxies separately, we find that galaxy properties such as color distribution, bulge-to-total ratios, and concentrations are remarkably similar between the void and overall sample. The void galaxies also show the same specific star formation rate at fixed color as the control galaxies. We compare our results with the predictions of cosmological simulations of galaxy formation using the Millennium Run semi-analytic galaxy catalog. We show that the properties of the simulated galaxies in large voids are in reasonably good agreement with those found in similar environments in the real Universe. To summarize, in spite of the fact that galaxies in voids live in the least dense large-scale environment, this environment makes very little impact on properties of galaxies.

We use measurements of the stellar velocity dispersion profile of the Fornax dwarf spheroidal galaxy to derive constraints on its dark matter distribution. Though the data are unable to distinguish between models with small cores and those with cusps, we show that a large > 1 kpc dark matter core in Fornax is highly implausible. Irrespective of the origin of the core, reasonable dynamical limits on the mass of the Fornax halo constrain its core radius to be no larger than 700 pc. We derive an upper limit core radius of 300 pc by demanding that the central phase space density of Fornax not exceed that directly inferred from the rotation curves of low-mass spiral galaxies. Further, if the halo is composed of warm dark matter then phase-space constraints force the core to be quite small in order to avoid conservative limits from the Ly alpha forest power spectrum, implying a core radius < 85 pc. We discuss our results in the context of the idea that the extended globular cluster distribution in Fornax can be explained by the presence of a large 1.5 kpc core. A self-consistent core of this size would be drastically inconsistent with the expectations of standard warm or cold dark matter models, and would also require an unreasonably massive dark matter halo, with a maximum circular velocity of 200 km/s.

We use the Sloan Digital Sky Survey (SDSS) spectroscopic sample to constrain the projected radial distribution of satellites around isolated ~ L* galaxies. We employ mock galaxy catalogs derived from high-resolution cosmological simulations to investigate the effects of interloper contamination and show that interlopers significantly bias the estimated slope of the projected radial distribution of satellites. We also show that the distribution of interlopers around galaxies is expected to be non-uniform in velocity space because galaxies are clustered and reside in crowded environments. Successful methods of interloper contamination correction should therefore take into account environments of the host galaxies. Two such new methods are presented and the most reliable of them is used to correct for interloper contamination in analyses of the SDSS galaxy sample. The best fit power-law slope of the interloper-corrected surface density distribution of satellites, Sigma(R) ~ R^alpha, in the volume-limited SDSS sample is alpha ~ -1.7 +/- 0.1, independent of the galaxy and satellite luminosities. Comparison with Lambda-CDM simulations shows that the radial distribution of the SDSS satellites is more concentrated than that of subhalos around galaxy-sized halos, especially at R < 100 kpc/h. The predicted dark matter radial distribution is somewhat more concentrated than the profile of the SDSS satellites, but the difference is not statistically significant for our sample.

The central part of a dark matter halo reacts to the presence and evolution of a bar. Not only does the halo absorb angular momentum from the disk, it can also be compressed and have its shape modified. We study these issues in a series of cosmologically motivated, highly resolved N-body simulations of barred galaxies run under different initial conditions. In all models we find that the inner halo's central density increases. We model this density increase using the standard adiabatic approximation and the modified formula by Gnedin et al. and find that halo mass profiles are better reproduced by this latter. In models with a strong bar, the dark matter in the central region forms a bar-like structure (``dark matter bar''), which rotates together with the normal bar formed by the stellar component (``stellar bar''). The minor-to-major axial ratio of a halo bar changes with radius with a typical value 0.7 in the central disk region. DM bar amplitude is mostly a function of the stellar bar strength. Models in which the bar amplitude increases or stays roughly constant with time, initially large (40%-60%) misalignment between the halo and disk bars quickly decreases with time as the bar grows. The halo bar is nearly aligned with the stellar bar (~10 degrees lag for the halo) after ~2 Gyr. The torque, which the halo bar exerts on the stellar bar, can serve as a mechanism to regulate the angular momentum transfer from the disk to the halo.

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.

We analyze the angular distribution and the orbital rotation directions of a sample of carefully-selected satellite galaxies extracted from the Sloan Digital Sky Survey (SDSS). We also study these statistics in an N-body simulation of cosmological structure formation set within the LCDM paradigm under various assumptions for the orientations of disk angular momenta. Assuming that the angular momenta of the disks are aligned with the angular momenta of the inner regions of their host dark matter halos, we find that the fraction of simulated satellite halos that exhibit prograde motion is 0.55-0.60, with larger satellites more likely to be prograde. In our observational sample, approximately 60% of the satellites exhibit prograde motion, a result that is broadly consistent with the simulated sample. Contrary to several recent studies, our observational sample of satellite galaxies show no evidence for being anisotropically distributed about their primary disks. Again, this result is broadly consistent with our simulated sample of satellites under the assumption that disk and halo angular momenta are aligned. However, the small size of our observational sample does not yet allow us to distinguish between various assumptions regarding the orientations of disks in their halos. Finally, we assessed the importance of contamination by interlopers on the measured prograde and retrograde statistics.

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.

It was argued in the past that bulges of galaxies cannot be formed through collisionless secular evolution because that would violate constraints on the phase-space density: the phase-space density in bulges is several times larger than in the inner parts of discs. We show that these arguments against secular evolution are not correct. Observations give estimates of the coarsely grained phase-space densities of galaxies, f'=rho_s/(sigma_R sigma_phi sigma_z), where rho_s is stellar density and sigma_R, sigma_phi, sigma_z are the radial, tangential, and vertical rms velocities of stars. Using high-resolution N-body simulations, we study the evolution of f' in discs of Galaxy-size models. During the secular evolution, the discs, which are embedded in live CDM haloes, form a bar and then a thick, dynamically hot, central mass concentration. In the course of evolution f' declines at all radii, not just in the central region. However, the decline is different in different parts of the disc. In the inner disc, f'(R) develops a valley with a minimum around the end of the central mass concentration. The final result is that the values of f' in the central regions are significantly larger than those in the inner disc. The minimum, which gets deeper with time, seems to be due to a large phase mixing produced by the outer bar. We find that the shape and the amplitude of f'(R) for different simulations agree qualitatively with the observed f'(R) in our Galaxy. Curiously enough, the fact that the coarsely grained phase-space density of the bulge is significantly larger than the one of the inner disc turns out to be an argument in favor of secular formation of bulges, not against it.

We present a study of the spatial distribution of subhalos in galactic dark matter halos using dissipationless cosmological simulations of the concordance LCDM model. We find that subhalos are distributed anisotropically and are preferentially located along the major axes of the triaxial mass distributions of their hosts. The Kolmogorov-Smirnov probability for drawing our simulated subhalo sample from an isotropic distribution is P_KS \simeq 1.5 \times 10^{-4}. An isotropic distribution of subhalos is thus not the correct null hypothesis for testing the CDM paradigm. The nearly planar distribution of observed Milky Way (MW) satellites is marginally consistent (probability \simeq 0.02) with being drawn randomly from the subhalo distribution in our simulations. Furthermore, if we select the subhalos likely to be luminous, we find a distribution that is consistent with the observed MW satellites. In fact, we show that subsamples of the subhalo population with a centrally-concentrated radial distribution, similar to that of the MW dwarfs, typically exhibit a comparable degree of planarity. We explore the origin of the observed subhalo anisotropy and conclude that it is likely due to (1) preferential accretion of subhalos along filaments, often closely aligned with the major axis of the host halo, and (2) evolution of satellite orbits within the prolate, triaxial potentials typical of CDM halos. Agreement between predictions and observations requires the major axis of the outer dark matter halo of the Milky Way to be nearly perpendicular to the disk. We discuss possible observational tests of such disk-halo alignment with current large galaxy surveys.

The dynamics of the inner disk of many galaxies is dominated by bars. As a result, resonances between the bar and the disk become an important factor. In order to detect resonances, we measure angular and radial frequencies of individual orbits in N-body models of barred galaxies. In our simulations with several million particles, stellar disks are immersed in live halos of dark matter. A bar forms spontaneously and interacts with the disk through resonances. We detect 9 of them. The corotation (CR) and the inner Lindblad resonance (ILR) are the most significant. The spatial distributions of particles trapped at these resonances show that ILR is not localized at a given radius. Particles trapped at ILR lie along the bar spanning a wide range of radii. On the other hand, particles at the corotation resonance form a ring around the corotation radius. We study how orbits evolve near corotation and how they are trapped. When an orbit is getting trapped, it evolves strongly, but after that there is little evolution. Those trapped orbits circulate along the corotation ring or librate around a stable lagragian point. The same technique is applied to study resonances in the halo. The bar interacts strongly with the halo particles that stay close to the disk. We find the corotation resonance and the inner Lindblad resonance for halo particles.

The cooling of gas in the centers of dark matter halos is expected to lead to a more concentrated dark matter distribution. The response of dark matter to the condensation of baryons is usually calculated using the model of adiabatic contraction, which assumes spherical symmetry and circular orbits. In contrast, halos in the hierarchical structure formation scenarios grow via multiple violent mergers and accretion along filaments, and particle orbits in the halos are highly eccentric. We study the effects of the cooling of gas in the inner regions of halos using high-resolution cosmological simulations which include gas dynamics, radiative cooling, and star formation. We find that the dissipation of gas indeed increases the density of dark matter and steepens its radial profile in the inner regions of halos compared to the case without cooling. For the first time, we test the adiabatic contraction model in cosmological simulations and find that the standard model systematically overpredicts the increase of dark matter density in the inner 5% of the virial radius. We show that the model can be improved by a simple modification of the assumed invariant from M(r)r to M(r_av)r, where r and r_av are the current and orbit-averaged particle positions. This modification approximately accounts for orbital eccentricities of particles and reproduces simulation profiles to within 10-20%. We present analytical fitting functions that accurately describe the transformation of the dark matter profile in the modified model and can be used for interpretation of observations.

We address the issue of accuracy in recovering density profiles from observations of rotation curves of galaxies. We ``observe'' and analyze our models in much the same way as observers do the real galaxies. We find that the tilted ring model analysis produces an underestimate of the central rotational velocity. In some cases the galaxy halo density profile seems to have a flat core, while in reality it does not. We identify three effects, which explain the systematic biases: (1) inclination (2), small bulge, and (3) bar. The presence of even a small non-rotating bulge component reduces the rotation velocity. In the case of a disk with a bar, the underestimate of the circular velocity is larger due to a combination of non-circular motions and random velocities. Signatures of bars can be difficult to detect in the surface brightness profiles of the model galaxies. The variation of inclination angle and isophote position angle with radius are more reliable indicators of bar presence than the surface brightness profiles. The systematic biases in the central ~ 1 kpc of galaxies are not large. Each effect separately gives typically a few kms error, but the effects add up. In some cases the error in circular velocity was a factor of two, but typically we get about 20 percent. The result is the false inference that the density profile of the halo flattens in the central parts. Our observations of real galaxies show that for a large fraction of galaxies the velocity of gas rotation (as measured by emission lines) is very close to the rotation of stellar component (as measured by absorption lines). This implies that the systematic effects discussed in this paper are also applicable both for the stars and emission-line gas.

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.

Hierarchical Cold Dark Matter (CDM) models predict that Milky Way (MW) sized halos contain hundreds of dense low-mass dark satellites, an order of magnitude more than the number of observed satellites in the Local Group (LG). If the CDM paradigm is correct, we need to understand why most of these halos failed to form stars and become galaxies. We analyze the dynamical evolution of the satellite halos in a high-resolution cosmological simulation of MW sized halos in the LCDM cosmology. We find that about 10% of the substructure halos with the present masses <10^8-10^9 Msun (circular velocities Vmax<30 km/s) had considerably larger masses and circular velocities when they formed at redshifts z>2. After the initial period of mass accretion in isolation, these objects experience dramatic mass loss due to tidal stripping, in some cases even before they are accreted by their host halo. This can explain how the smallest dwarf spheroidal galaxies of the LG were able to build up a sizable stellar mass in their seemingly shallow potential wells. We propose a new model in which all of the luminous dwarf spheroidals in the Local Group are descendants of the relatively massive (>10^9 Msun) high-redshift systems, in which the gas could cool efficiently by atomic line emission and which were not significantly affected by the extragalactic ultraviolet radiation. We present a simple galaxy formation model based on the trajectories extracted from the simulation, which accounts for the bursts of star formation after strong tidal shocks and the inefficiency of gas cooling in halos with virial temperatures Tvir<~10^4 K. Our model reproduces the abundance, spatial distribution, and morphological segregation of the observed Galactic satellites. The results are insensitive to the redshift of reionization.

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 a method for computing the 6-dimensional coarse-grained phase-space density $f(x,v)$ in an N-body system, and derive its distribution function $v(f)$. The method is based on Delaunay tessellation, where $v(f)$ is obtained with an effective fixed smoothing window over a wide $f$ range. The errors are estimated, and $v(f)$ is found to be insensitive to the sampling resolution or the simulation technique. We find that in gravitationally relaxed haloes built by hierarchical clustering, $v(f)$ is well approximated by a robust power law, $v(f) \propto f^{-2.5 \pm 0.05}$, over more than 4 decades in $f$, from its virial level to the numerical resolution limit. This is tested to be valid in the $\Lambda$CDM cosmology for haloes with masses $10^9-10^{15}\msun$, indicating insensitivity to the slope of the initial fluctuation power spectrum. By mapping the phase-space density in position space, we find that the high-$f$ end of $v(f)$ is dominated by the ``cold'' subhaloes rather than the parent-halo central region and its global spherical profile. The value of $f$ in subhaloes near the virial radius is typically $>100$ times higher than its value at the halo centre, and it decreases gradually from outside in toward its value at the halo centre. This seems to reflect phase mixing due to mergers and tidal effects involving puffing up and heating. The phase-space density can thus provide a sensitive tool for studying the evolution of subhaloes during the hierarchical buildup of haloes. It remains to be understood why the evolved substructure adds up to the actual universal power law of $v(f) \propto f^{-2.5}$. It seems that this behaviour results from the hierarchical clustering process and is not a general result of violent relaxation.

If dark matter (DM), which is considered to constitute most of the mass of galaxies, is made of supersymmetric (SUSY) particles, the centers of galaxies should emit gamma-rays produced by their self-annihilation. We present accurate estimates of continuum gamma-ray fluxes due to neutralino annihilation in the central regions of the Milky Way. We use detailed models of our Galaxy, which satisfy available observational data, and include some important physical processes, which were previously neglected. Our models predict that spatially extended annihilation signal should be detected at high confidence levels by incoming experiments assuming that neutralinos make up most of the DM in the Universe and that they annihilate according to current SUSY models.

We perform N-body simulations for models with a DE component. Besides of DE with constant negative state parameter w, we consider DE due to scalar fields, self-interacting through RP or SUGRA potentials. According to our post-linear analysis, at z=0, DM power spectra and halo mass functions do not depend on DE nature. This is welcome, as LCDM fits observations. Halo profiles, instead, are denser than LCDM. For example, the density at 15 kpc of a DE halo with M=10^13 exceeds LCDM by ~45%. Differences, therefore, are small but, however, DE does not ease the problem with cuspy DM profiles. We study also subhalos and find that, at $z=0$, the number of satellites coincides in all DE models. At higher z, DE models show increasing differences from LCDM and among themselves (i.e. in the mass function evolution); this is the obvious pattern to distinguish between different DE state equations.

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.

Using high resolution N-body simulations of stellar disks embedded in cosmologically motivated dark matter halos, we study the evolution of bars and the transfer of angular momentum between halos and bars. We find that dynamical friction results in some transfer of angular momentum to the halo, but the effect is much smaller than previously found in low resolution simulations and is incompatible with early analytical estimates. In simulations with millions of particles reaching a resolution of 20-40 pc, the pattern speed may not change over billions of years. Our high resolution models produce bars which are fast rotators, where the ratio of the corotation radius to the bar major semi-axis lies in the range R =1.2-1.7, marginally compatible with observational results. In contrast to many previous simulations, we find that bars are relatively short. As in many observed cases, the bar major semi-axis is close to the exponential length of the disk. The transfer of angular momentum between inner and outer parts of the disk plays a very important role in the secular evolution of the disk and the bar. The bar formation increases the exponential length of the disk by a factor of 1.2 -1.5. The transfer substantially increases the stellar mass in the centre of the galaxy and decreases the dark matter-to-baryons ratio. As the result, the central 2 kpc region is always strongly dominated by the baryonic component. At intermediate 3-10kpc scales the disk is sub-dominant. In summary, realistic models produce bar structure in striking agreement with observational results.

Cosmological models with different types of Dark Energy are becoming viable alternatives for standard models with the cosmological constant. Yet, such models are more difficult to analyze and to simulate. We present analytical approximations and discuss ways of making simulations for two families of models, which cover a wide range of possibilities and include models with both slow and fast changing ratio w=p\rho. More specifically, we give analytical expressions for the evolution of the matter density parameter Omega_m(z) and the virial density contrast Delta_c at any redshift z. The latter is used to identify halos and to find their virial masses. We also provide an approximation for the linear growth factor of linear fluctuations between redshift z=40 and z=0. This is needed to set the normalization of the spectrum of fluctuations. Finally, we discuss the expected behavior of the halo mass function and its time evolution.

We study properties of dark matter halos in a variety of models which include Dark Energy (DE). We consider both DE due to a scalar field self--interacting through Ratra-Peebles or SUGRA potentials, and DE with constant negative w=p\rho >-1. We find that at redshift zero the nonlinear power spectrum of the dark matter, and the mass function of halos, practically do not depend on DE state equation and are almost indistinguishable from predictions of the LCDM model. This is consistent with the nonlinear analysis presented in the accompanying paper. It is also a welcome feature because LCDM models fit a large variety of data. On the other hand, at high redshifts DE models show substantial differences from LCDM and substantial differences among themselves. Halo profiles differ even at z=0. DE halos are denser than LCDM in their central parts because the DE halos collapse earlier. Nevertheless, differences between the models are not so large. For example, the density at 10 kpc of a DE \~10^{13}Msun halo deviates from LCDM by not more than 50%. This, however, means that DE is not a way to ease the problem with cuspy dark matter profiles. Addressing another cosmological problem - abundance of subhalos -- we find that the number of satellites of halos in various DE models does not change relative to the LCDM, when normalized to the same circular velocity of the parent halo. To summarize, the best way to find which DE model fits the observed Universe is to look for evolution of halo properties. For example, the abundance of galaxy groups with mass larger than 10^{13}\Msun at z> 2 can be used to discriminate between the models, and, thus, to constrain the nature of DE.

Using the Sloan Digital Sky Survey (SDSS), we probe the halo mass distribution by studying the velocities of satellites orbiting isolated galaxies. In a subsample that covers 2500 sq. degrees on the sky, we detect about 3000 satellites with absolute blue magnitudes going down to M_B = -14; most of the satellites have M_B=-16 to -18, comparable to the magnitudes of M32 and the Magellanic Clouds. After a careful, model-independent removal of interlopers, we find that the line-of-sight velocity dispersion of satellites declines with distance to the primary. For an L* galaxy the r.m.s. line-of-sight velocity changes from ~120 km/s at 20 kpc to ~60 km/s at 350 kpc. This decline agrees remarkably well with theoretical expectations, as all modern cosmological models predict that the density of dark matter in the peripheral parts of galaxies declines as rho_DM propto r^{-3}. Thus, for the first time we find direct observational evidence of the density decline predicted by cosmological models; we also note that this result contradicts alternative theories of gravity such as MOND. We also find that the velocity dispersion of satellites within 100 kpc scales with the absolute magnitude of the central galaxy as sigma propto L^{0.3}; this is very close to the Tully-Fisher relation for normal spiral galaxies.

We propose a new explanation for the origin of angular momentum in galaxies and their dark halos, in which the halos obtain their spin through the cumulative acquisition of angular momentum from satellite accretion. In our model, the build-up of angular momentum is a random walk process associated with the mass assembly history of the halo's major progenitor. We assume no correlation between the angular momenta of accreted objects. Using the extended Press-Schechter approximation, we calculate the growth of mass, angular momentum, and spin parameter $\lambda$ for many halos. Our random walk model reproduces the key features of the angular momentum of halos found in N-body simulations: a lognormal distribution in $\lambda$ with an average of $<\lambda> \approx 0.04$, independent of mass and redshift. The evolution of the spin parameter in individual halos in this model is quite different from the steady increase with time of angular momentum in the tidal torque picture. We find both in N-body simulations and in our random walk model that the value of $\lambda$ changes significantly with time for a halo's major progenitor. It typically has a sharp increase due to major mergers, and a steady decline during periods of gradual accretion of small satellites. The model predicts that on average the $\lambda$ of halos which had major mergers after redshift $z=2$ should be substantially larger than the $\lambda$ of those which did not. Perhaps surprisingly, this suggests that halos that host late-forming elliptical galaxies should rotate faster than halos of spiral galaxies.

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).

We apply standard disk formation theory with adiabatic contraction within cuspy halo models predicted by the standard LCDM cosmology. The resulting models score remarkably well when confronted with the broad range of observational data available for the Milky Way and M31 galaxies, giving a Milky Way virial mass of 1-2x10^12Msun and concentration C=12-17. We consider two types of models, in which: (A) baryons conserve angular momentum and (B) some of the angular momentum of the baryons is transferred to the dark matter. Type-A models produce good agreement with observed rotation curves and obey constraints in the solar neighborhood, but may have too much dark matter in the center to allow a fast rotating bar. The type-B models with angular momentum transport have a slightly more massive disk and less dark matter in the central part, allowing a fast rotating bar to persist. Both classes of models probably have sufficient baryonic mass in the central 3.5kpc to reproduce recent observational values of the optical depth to microlensing events towards the Galactic center. All models require that about 1/2 of all baryons expected inside the virial radius must not be in the disk or bulge. We investigate whether the range of virial masses allowed by our dynamical models is compatible with constraints from the galaxy luminosity function, and find a range of parameter space that is allowed by this constraint. We conclude that rotation curves and dynamical properties of ``normal'' high surface brightness spiral galaxies appear to be consistent with standard LCDM.

We present results of gasdynamics+N-body constrained cosmological simulations of the Local Supercluster region (LSC; about 30/h Mpc around the Virgo cluster), which closely mimic the real Universe within 100 Mpc by imposing constraints from Mark III catalog of galaxy peculiar velocities. The simulations are used to study the properties and possible observational signatures of intergalactic medium in the LSC region. We find that, in agreement with previous unconstrained simulations, about 30% of the gas in this region is in the warm/hot phase at T~10^5-10^7 K, and about 40% in the diffuse phase at T<10^5 K in low-density regions. The best prospects for detection of the warm/hot intergalactic medium of the LSC located in filaments and in the vicinity of virialized regions of groups and clusters are through absorption in resonant lines of OVII and OVIII in soft X-rays and in the OVI doublet in UV. If intergalactic gas in filaments (rho/<rho>~ 1-10) is enriched to typical metallicities of >0.05, the column densities of OVI, OVII, and OVIII along a random line of sight near the North Galactic Pole, especially near the supergalactic plane, have a significant probability to be in the range detectable by the current (FUSE,XMM) and future (Constellation-X) instruments. (Abridged)

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.

We study the effects of mass and force resolution on the density profiles of galaxy-size Cold Dark Matter (CDM) halos in a flat, low-density cosmological model with vacuum energy. We show thatalthough increasing the mass and force resolution allows us to probe deeper into the inner halo regions, it does not lead to steeper inner density profiles. Instead, the halo profiles converge at scales larger than four times the formal resolution or the radius containing more than 200 particles, whichever is larger. In the simulations presented in this paper, we are able to probe density profile of a relaxed isolated galaxy-size halo at scales r=(0.005-1)r_vir. We find that the density distribution can be well approximated by the profile suggested by Moore etal (1998): rho= x^{-1.5}(1+x^{1.5})^-1, where x=r/r_s and r_s is the characteristic radius. The analytical profile proposed by Navarro et al. (1996) rho= x^{-1}(1+x)^-2, also provides a good fit, with the same relative errors of about 10% for radii larger than 1% of the virial radius. Both analytical profiles fit equally well because for high-concentration galaxy-size halos the differences between these profiles become significant only at scales well below 0.01r_vir. We also find that halos of similar mass may have somewhat different parameters (characteristic radius, maximum rotation velocity, etc.) and shapes of their density profiles. We associate this scatter in properties with differences in halo merger histories and the amount of substructure present in the analyzed halos.

[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv \mu j/(j_0+j). This profile is defined by one shape parameter (\mu or j_0) in addition to the global spin parameter \lambda. It follows a power-law over most of the mass, and flattens at large j, with the flattening more pronounced for small values of \mu. Compared to a uniform sphere in solid-body rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. The spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each halo for ~80% of the halos. We investigate two ideas for the origin of this profile. The first is based on a revised version of linear tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on j transport in minor mergers. Finally, we briefly explore implications of the M(<j) profile on the formation of galactic disks assuming that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent of the "angular-momentum problem" seen in hydrodynamic simulations, even though we have assumed perfect j conservation. A possible solution is to associate the central excesses with bulge components and the outer regions with extended gaseous disks.

We summarize recent developments in the study of the origin of halo spin profiles and preliminary implications on disk formation. The specific angular-momentum distributions within halos in N-body simulations match a universal shape, M(<j) \propto j/(j_0+j). It is characterized by a power law over most of the mass, and one shape parameter in addition to the spin parameter \lambda. The angular momentum tends to be aligned throughout the halo and of cylindrical symmetry. Even if angular momentum is conserved during baryonic infall, the resultant disk density profile is predicted to deviate from exponential, with a denser core and an extended tail. A slightly corrected version of the scaling relation due to linear tidal-torque theory is used to explain the origin of a typical power-law profile in shells, j(M) \propto M^s with s \gsim 1. While linear theory crudely predicts the amplitudes of halo spins, it is not a good predictor of their directions. Independently, mergers of halos are found to produce a similar profile due to j transfer from the orbit to the product halo via dynamical friction and tidal stripping. The halo spin is correlated with having a recent major merger, though this correlation is weakened by mass loss. These two effects are due to a correlation between the spins of neighboring halos and their orbit, leading to prograde mergers.

[Abridged] We analyze the exchange of dark matter between halos, subhalos, and their environments in a high-resolution cosmological N-body simulation of a Lambda CDM cosmology. At each analyzed redshift z we divide the dark matter particles into 4 components: (i) isolated galactic halos, (ii) subhalos, (iii) the diffuse medium of group and cluster halos, and (iv) the background outside of virialized halos. We follow the time evolution of the mass distribution and flows between these components and provide fitting functions for the exchange rates. We use our derived exchange rates to gauge the importance of metal redistribution in the universe due solely to gravity-induced interactions. The diffuse metallicity in clusters is predicted to be ~40% that in isolated galaxies (~55% of groups) at z=0, and should be lower only slightly by z=1, consistent with observations. The metallicity of the diffuse media in poor groups is expected to be lower by a factor of 5 by z~2, in agreement with the observed metallicity of damped Ly$\alpha$ systems. The metallicity of the background IGM is predicted to be (1-3)x10^{-4} that of z=0 clusters, also consistent with observations. The agreement of predicted and observed trends indicates that gravitational interaction alone may play an important role in metal enrichment of the intra-cluster and intergalactic media.

[Abridged] The interaction rates of dark-matter halos and subhalos, are computed using high-resolution cosmological N-body simulations of the Lambda-CDM model. Although the number fraction of subhalos of mass $>2x10^{11}\hsolmass$ is only ~10%, we find that the interaction rate of such subhalos is relatively high because they reside in high density environments. At low redshift, the subhalo collisions dominate the total collision rate, and even at z=3 they are involved in more than 30% of all collisions. About 40% of the "major" mergers are between subhalos. Therefore subhalo interactions must be incorporated in models of structure formation. We find that the collision rate between halos in physical density units, is $\propto (1+z)^\delta $, with $\delta = 3-4$, in agreement with earlier simulations and most observational data. We test previous analytic estimates of the interaction rates of subhalos based on statistical models, which could be very inaccurate because of the small number of subhalos and the variation of conditions within small host halos. We find that, while such statistical estimates may severely overestimate the rate within hosts $< 10^{13} \hsolmass$, typical of high redshifts, they are valid for larger hosts regardless of the number of subhalos in them. We find the Makino & Hut (1997) estimate of the subhalo merger rate to be valid for hosts $ > 6x10^{11}\hsolmass$ at all redshifts. The collision rate between subhalos and the central object of their host halo is approximated relatively well using the timescale for dynamical friction in circular orbits. This approximation fails in ~40% of the cases, partly because of deviations from the assumption of circular orbits and partly because of the invalidity of the assumption that the subhalo mass is negligible.

We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the redshift dependence of the median concentration is c_v ~ 1/(1+z), corresponding to a roughly constant r_s with redshift. We present an improved analytic treatment of halo formation that fits the measured relations between halo parameters and their redshift dependence. The implications are that high-redshift galaxies are predicted to be more extended and dimmer than expected before. Second, we find that the scatter in log(c_v) is ~0.18, corresponding to a scatter in maximum rotation velocities of dV/V ~ 0.12. We discuss implications for modelling the Tully-Fisher relation, which has a smaller reported intrinsic scatter. Third, haloes in dense environments tend to be more concentrated than isolated haloes. These results suggest that c_v is an essential parameter for the theory of galaxy modelling, and we briefly discuss implications for the universality of the Tully-Fisher relation, the formation of low surface brightness galaxies, and the origin of the Hubble sequence.

We give a summary of recent results on spatial and velocity biases in cosmological models. Progress in numerical techniques made it possible to simulate halos in large volumes with a such accuracy that halos survive in dense environments of groups and clusters of galaxies. Halos in simulations look like real galaxies, and, thus, can be used to study the biases -- differences between galaxies and the dark matter. The biases depend on scale, redshift, and circular velocities of selected halos. Two processes seem to define the evolution of the spatial bias: (1) statistical bias and (2) merger bias (merging of galaxies, which happens preferentially in groups, reduces the number of galaxies, but does not affect the clustering of the dark matter). There are two kinds of velocity bias. The pair-wise velocity bias is b_12=0.6-0.8 at r< 5Mpc/h, z=0. This bias mostly reflects the spatial bias and provides almost no information on the relative velocities of the galaxies and the dark matter. One-point velocity bias is a better measure of the velocities. Inside clusters the galaxies should move slightly faster (b_v =1.1-1.3) than the dark matter. Qualitatively this result can be understood using the Jeans equations of the stellar dynamics. For the standard LCDM model we find that the correlation function and the power spectrum of galaxy-size halos at z=0 are antibiased on scales r<5Mpc/h and k=(0.15-30)h/Mpc.

Properties of dark matter halos are reviewed. Taken from different publications, we present results on (1) the mass and velocity functions, (2) density and velocity profiles, and (3) concentration of halos. In the range of radii r=(0.005-1)rvir the density profile for a quiet isolated halo is very accurately approximated by a fit suggested by Moore etal (1997): rho=1/x^1.5(1+x^1.5), where x=r/rs and rs is a characteristic radius. The fit suggested by Navarro et al (1995) rho= 1/x(1+x)^2, also gives a very satisfactory approximation with relative errors of about 10% for radii not smaller than 1% of the virial radius. The mass function of z=0 halos with mass below 10^{13}Msun/h is approximated by a power law with slope alpha =-1.85. The slope increases with the redshift. The velocity function of halos with Vmax< 500km/s is also a power law with the slope beta= -3.8-4. The power-law extends to halos at least down to 10km/s. It is also valid for halos inside larger virialized halos. The concentration of halos depends on mass (more massive halos are less concentrated) and environment, with isolated halos being less concentrated than halos of the same mass inside clusters. Halos have intrinsic scatter of concentration: at 1sigma level halos with the same mass have Delta(log{c_vir})=0.18 or, equivalently, DeltaVmax/Vmax =0.12. Velocity anisotropy for both subhalos and the dark matter is approximated by beta(r) =0.15+2x/[x^2+4], where x is radius in units of the virial radius.

This lecture gives a short introduction to different methods used in cosmology. The focus is on major features of N-body simulations: equations, main numerical techniques, effects of resolution, and methods of halo identification.

We study the distribution functions of mass and circular velocity for dark matter halos in N-body simulations of the $\Lambda$CDM cosmology, addressing redshift and environmental dependence. The dynamical range enables us to resolve subhalos and distinguish them from "distinct" halos. The mass function is compared to analytic models, and is used to derive the more observationally relevant circular velocity function. The distribution functions in the velocity range 100--500 km/s are well fit by a power-law with two parameters, slope and amplitude. We present the parameter dependence on redshift and provide useful fitting formulae. The amplitudes of the mass functions decrease with z, but, contrary to naive expectation, the comoving density of halos of a fixed velocity ~200 km/s actually increases out to z=5. This is because high-z halos are denser, so a fixed velocity corresponds to a smaller mass. The slope of the velocity function at z=0 is as steep as ~ -4, and the mass and velocity functions of distinct halos steepen with increasing z, while the functions of subhalos do not steepen with $z$, and become even flatter at z>2. A simple observable prediction is that the slope of the velocity function of isolated galaxies is steeper than that of galaxies in groups by as much as unity, reflecting the density biasing of high-velocity halos. We confirm that the Press-Schechter approximation typically overestimates the halo mass function by a factor of ~2, while modified approximations provide improved predictions.

We have simulated the evolution of a large sample of X-ray clusters in a Lambda-CDM universe at high spatial resolution using adaptive mesh refinement and placed the results in an online archive for public access. The Simulated Cluster Archive website http://sca.ncsa.uiuc.edu provides tools for interactive 2D and 3D analysis of gas and dark matter fields, X-ray and SZ imaging, and data export. We encourage community use and solicit their feedback.

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.

We present results of extensive model calculations of disk galaxy evolution within an hierarchical inside-out formation scenario. We first compare properties of the dark halos identified in a cosmological N-body simulation with predictions of a seminumerical method based on an extended collapse model and find a good agreement. We also study the properties of the halos in dependence on their environment. We then describe detailed modelling of the formation and evolution of disks within the growing isolated cold dark matter halos and predictions for the main properties, correlations and evolutionary features of normal disk galaxies. The possible shortcomings of the scenario are discussed.

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.

We make available to the public 4 dark matter simulations that cover 4 CDM cosmological models and were run with the AP^3M code (Couchamn 1991). These simulations are described in Colin et al. (ApJ, 523, 32, 1999). There are about 20 epochs available per model. They were run in a 30 Mpc/h box and have 128^3 particles in a 128^3 mesh. To download the binary data please contact me at: colin@astroscu.unam.mx

We give a summary of our recent studies of spatial and velocity biases of galaxy-size halos in cosmological models. Recent progress in numerical techniques made it possible to simulate halos in large volumes with a such accuracy that halos survive in dense environments of groups and clusters of galaxies. Dark matter halos in simulations look like real galaxies, and, thus, can be used to study the biases - differences between galaxies and the dark matter. For the standard LCDM model we find that the correlation function and the power spectrum of galaxy-size halos at z=0 are antibiased on scales r<5Mpc/h and k=(0.15-30)h/Mpc. The biases depend on scale, redshift, and circular velocities of selected halos. Two processes seem to define the evolution of the spatial bias: (1) statistical bias (or peak bias) and merger bias (merging of galaxies, which happens preferentially in groups, reduces the number of galaxies, but does not affect the clustering of the dark matter). There are two kinds of velocity bias. The pair-wise velocity bias is b_12=0.6-0.8 at r<5Mpc/h, z=0. This bias mostly reflects the spatial bias and provides almost no information on the relative velocities of the galaxies and the dark matter. One-point velocity bias is a better measure of the velocities. Inside clusters the galaxies should move slightly faster (b_v=1.1-1.3) than the dark matter. Qualitatively this result can be understood using the Jeans equations of the stellar dynamics.

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 N-body simulations to study the velocity bias of dark matter halos, the difference in the velocity fields of dark matter and halos, in a flat low- density LCDM model. The high force, 2kpc/h, and mass, 10^9Msun/h, resolution allows dark matter halos to survive in very dense environments of groups and clusters making it possible to use halos as galaxy tracers. We find that the velocity bias pvb measured as a ratio of pairwise velocities of the halos to that of the dark matter evolves with time and depends on scale. At high redshifts (z ~5) halos move generally faster than the dark matter almost on all scales: pvb(r)~1.2, r>0.5Mpc/h. At later moments the bias decreases and gets below unity on scales less than r=5Mpc/h: pvb(r)~(0.6-0.8) at z=0. We find that the evolution of the pairwise velocity bias follows and probably is defined by the spatial antibias of the dark matter halos at small scales. One-point velocity bias b_v, defined as the ratio of the rms velocities of halos and dark matter, provides a more direct measure of the difference in velocities because it is less sensitive to the spatial bias. We analyze b_v in clusters of galaxies and find that halos are ``hotter'' than the dark matter: b_v=(1.2-1.3) for r=(0.2-0.8)r_vir, where r_vir is the virial radius. At larger radii, b_v decreases and approaches unity at r=(1-2)r_vir. We argue that dynamical friction may be responsible for this small positive velocity bias b_v>1 found in the central parts of clusters. We do not find significant difference in the velocity anisotropy of halos and the dark matter. The dark matter the velocity anisotropy can be approximated as beta(x)=0.15 +2x/(x^2+4), where x is measured in units of the virial radius.

(Abridged) We study the outer density profiles of dark matter haloes predicted by a generalized secondary infall model and observed in a N-body cosmological simulation of a \Lambda CDM model. We find substantial systematic variations in shapes and concentrations of the halo profiles as well as a strong correlation of the profiles with the environment. In the N-body simulation, the average outer slope of the density profiles, \beta (\rho\propto r^{-\beta}), of isolated haloes is \approx 2.9; 68% of these haloes have values of \beta between 2.5 and 3.8. Haloes in dense environments of clusters are more concentrated and exhibit a broad distribution of \beta with values larger than for isolated haloes . Contrary to what one may expect, the haloes contained within groups and galaxy systems are less concentrated and have flatter outer density profiles than the isolated haloes. The concentration decreases with M_h, but its scatter for a given mass is substantial. The mass and circular velocity of the haloes are strongly correlated: M_h \propto V_m^{\alpha} with \alpha ~ 3.3 (isolated) and ~3.5 (haloes in clusters). For M_h=10^12M_sun the rms deviations from these relations are \Delta logM_h=0.12 and 0.18, respectively. Approximately 30% of the haloes are contained within larger haloes or have massive companions (larger than ~0.3 the mass of the current halo) within 3 virial radii. The remaining 70% of the haloes are isolated objects. The distribution of \beta as well as the concentration-mass and M_h-V_m relations for the isolated haloes agree very well with the predictions of our seminumerical approach which is based on a generalization of the secondary infall model and on the extended Press-Schechter formalism.

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 have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.

Lyman break galaxies (LBGs) at $z \sim 3$ exhibit number densities and clustering similar to local $L_*$ galaxies; however, their star formation rates (SFRs) are much higher. We explore the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an explanation for these key properties. The relative paucity of starburst galaxies at low redshift can be attributed to a much higher collision rate in the past. We use high-resolution cosmological $N$-body simulations and a hierarchical halo finder to estimate the galaxy collision rate as a function of time in the popular $\Lambda$CDM cosmological model. We find that bright collisional starbursts are frequent enough to account for most of the high-$z$ ($2.5-4.5$) LBGs. Moreover, many of the objects are of relatively small mass, but they cluster about large-mass halos. They therefore exhibit strong clustering, similar to that observed and stronger than that of the relevant massive halos.

Using published data, we have compiled the circular velocity (Vc) distribution function (VDF) of galaxy satellites in the Local Group. We find that within the volumes of radius of 570 kpc (400/h kpc for h=0.7) centered on the Milky Way and Andromeda, the average VDF is roughly approximated as n(>Vc)~ 45(Vc/10 km/s)^{-1} h^3 Mpc^{-3} for Vc in the range ~10-70 km/s. The observed VDF is compared with results of high-resolution cosmological simulations. We find that the VDF in models is very different from the observed one: n(>Vc)~1200(Vc/10 km/s)^{-2.75}h^3 Mpc^{-3}. Cosmological models thus predict that a halo of the size of our Galaxy should have about 50 dark matter satellites with circular velocity >20 km/s and mass >3x10^8/h Msun within a 570 kpc radius. This number is significantly higher than the approximate dozen satellites actually observed around our Galaxy. The observed and predicted VDFs cross at ~50 km/s, indicating that the predicted abundance of satellites with Vc> 50 km/s is in reasonably good agreement with observations. We conclude, therefore, that unless a large fraction of the Local Group satellites has been missed in observations, there is a dramatic discrepancy between observations and hierarchical models, regardless of the model parameters. We discuss several possible explanations for this discrepancy including identification of some satellites with the High Velocity Clouds observed in the Local Group, and the existence of dark satellites that failed to accrete gas and form stars due either to the expulsion of gas in the supernovae-driven winds or to gas heating by the intergalactic ionizing background. (Abridged)

We present results from a study of bias and its evolution for galaxy-size halos in a large, high-resolution simulation of a LCDM model. We consider the evolution of bias estimated using two-point correlation function (b_xi), power spectrum (b_P), and a direct correlation of smoothed halo and matter overdensity fields (b_d). We present accurate estimates of the evolution of the matter power spectrum probed deep into the stable clustering regime (k~[0.1-200]h/Mpc at z=0). The halo power spectrum evolves much slower than the power spectrum of matter and has a different shape which indicates that the bias is time- and scale-dependent. At z=0, the halo power spectrum is anti-biased with respect to the matter power spectrum at wavenumbers k~[0.15-30]h/Mpc, and provides an excellent match to the power spectrum of the APM galaxies at all probed k. In particular, it nicely matches the inflection observed in the APM power spectrum at k~0.15h/Mpc. We complement the power spectrum analysis with a direct estimate of bias using smoothed halo and matter overdensity fields and show that the evolution observed in the simulation in linear and mildly non-linear regimes can be well described by the analytical model of Mo & White (1996), if the distinction between formation redshift of halos and observation epoch is introduced into the model. We present arguments and evidence that at higher overdensities, the evolution of bias is significantly affected by dynamical friction and tidal stripping operating on the satellite halos in high-density regions of clusters and groups; we attribute the strong anti-bias observed in the halo correlation function and power spectrum to these effects. (Abridged)

This article concerns the formation and structure of dark matter halos, including (1) their radial density profiles, (2) their abundance, and (3) their merger rates. The last topic may be relevant to the nature of the small, bright, high-redshift galaxies discovered by the Lyman break technique. (1) Study of a statistical sample of galaxy-mass dark halos in high-resolution Adaptive Refinement Tree simulations shows that they have a central density profile p(r) \propto 1/r^g with g \approx 0.2, in agreement with data on dark-matter-dominated disk galaxies. We present recent, higher resolution results on this. (2) Another important new result is that the Press-Schechter approximation predicts about twice as many galaxy-mass halos at z=0 as are present in large dissipationless N-body simulations; more generally, PS overpredicts the abundance of M ~< 0.1M_* halos at all redshifts. (3) Finally, we discuss the assembly of these halos, in particular the merger rate of (sub-)halos at high redshift and the distribution of the starbursts that these mergers are likely to trigger. If most of the Lyman-break galaxies are such starbursts, this perhaps resolves the apparent paradox that these galaxies appear to cluster like massive halos (~10^{12} M_\odot), while their relatively low linewidths and their spectral energy distributions suggest that they have relatively low mass (few X 10^{10} M_\odot) and young ages (few X 10^8 yr). It also predicts much more star formation at high redshift in CDM-type hierarchical models for structure formation than if only quiescent star formation is included.

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.

We study the evolution of the halo-halo correlation function and bias in a LCDM model using very high-resolution N-body simulations. The high force and mass resolution allows dark matter (DM) halos to survive in the tidal fields in high-density regions and thus prevents the ambiguities related with the ``overmerging problem.'' Numbers of galaxy-size halos in cluster-like objects in our simulation are similar to the numbers of galaxies observed in real clusters. This allows us to estimate for the first time the evolution of the correlation function and bias at small (down to ~100/h kpc) scales. We compare particle distribution, dark matter correlation function, density profiles, and halo mass function produced with our N-body code and corresponding results of the AP3M simulations. We find that at all epochs the 2-point correlation function of galaxy-size halos xihh is well approximated by a power-law with slope ~1.6-1.8. The difference between the shape of xihh and the shape of the correlation function of matter results in the scale-dependent bias at scales <7/h Mpc, which we find to be a generic prediction of the hierarchical models, independent of the epoch and of the model details. We find that our results agree well with existing clustering data at different redshifts, indicating the general success of the picture of structure formation in which galaxies form inside the host DM halos. Particularly, we find an excellent agreement in both slope and the amplitude between xihh(z=0) in our simulation and the galaxy correlation function measured using the APM galaxy survey. At high redshifts, the observed clustering of the Lyman-break galaxies is also reasonably well reproduced by the models.

We study the evolution of the halo-halo correlation function and bias in four cosmological models (LCDM, OCDM, tauCDM, and SCDM) using very high-resolution N-body simulations. The high force and mass resolution allows dark matter (DM) halos to survive in the tidal fields of high-density regions and thus prevents the ambiguities related with the ``overmerging problem.'' This allows us to estimate for the first time the evolution of the correlation function and bias at small (down to ~100/h kpc) scales. We find that at all epochs the 2-point correlation function of galaxy-size halos xi_hh is well approximated by a power-law with slope ~1.6-1.8. The difference between the shape of xi_hh and the shape of the correlation function of matter results in the scale-dependent bias at scales <7/h Mpc, which we find to be a generic prediction of the hierarchical models. The bias evolves rapidly from a high value of ~2-5 at z~3-7 to the anti-bias of b~0.5-1 at small <5/h Mpc scales at z=0. We find that our results agree well with existing clustering data at different redshifts. Particularly, we find an excellent agreement in both slope and the amplitude between xi_hh(z=0) in our LCDM simulation and the galaxy correlation function measured using the APM galaxy survey. At high redshifts, the observed clustering of the Lyman-break galaxies is also well reproduced by the models. The agreement with the data at high and low z indicates the general success of the hierarchical models of structure formation in which galaxies form inside the host DM halos. (Abridged)

Using N-body+hydro simulations we study relations between the local environments of galaxies on 0.5 Mpc scale and properties of the luminous components of galaxies. Our numerical simulations include effects of star formation and supernova feedback in different cosmological scenarios: the standard Cold Dark Matter model, the Broken Scale Invariance model (BSI), and a model with cosmological constant (LCDM). In this paper, we concentrate on the effects of environment on colors and morphologies of galaxies, on the star formation rate and on the relation between the total luminosity of a galaxy and its circular velocity. We demonstrate a statistically significant theoretical relationship between morphology and environment. In particular, there is a strong tendency for high-mass galaxies and for elliptical galaxies to form in denser environments, in agreement with observations. We find that in models with denser environments (CDM scenario) ~ 13 % of the galactic halos can be identified as field ellipticals, according to their colors. In simulations with less clustering (BSI and LCDM), the fraction of ellipticals is considerably lower (~ 2-3 %). The strong sensitivity of morphological type to environment is rather remarkable because our results are applicable to ``field'' galaxies and small groups. If all galaxies in our simulations are included, we find a statistically significant dependence of the galaxy luminosity - circular velocity relation on dark matter overdensity within spheres of radius 0.5 Mpc, for the CDM simulations. But if we remove ``elliptical'' galaxies from our analysis to mimic the Tully-Fisher relation for spirals, then no dependence is found in any model.

Using 3D hydrodynamical simulations of galaxy formation with supernova feedback and a multiphase medium, we derive theoretical relations analogous to the observed Tully-Fisher (TF) relations in various photometric bands. This paper examines the influence of self-regulation mechanisms including supernova feedback on galaxy luminosities and the TF relation in three cosmological scenarios (CDM, \Lambda CDM and BSI (broken scale invariance)). The galaxy catalogs derived from our hydrodynamical simulations lead to an acceptably small scatter in the theoretical TF relation amounting to \Delta M =0.2-0.4 in the I band, and increasing by 0.1 magnitude from the I-band to the B-band. Our results give strong evidence that the tightness of the TF relation cannot be attributed to supernova feedback alone. However, although eliminating supernova feedback hardly affects the scatter, it does influence the slope of the TF relation quite sensitively. With supernova feedback, L \propto V_c^{3-3.5} (depending on the strength of supernova feedback). Without it, L \propto V_c^{2} as predicted by the virial theorem with constant M/L. The luminosity functions in the B and K bands are quite sensitive to supernova feedback at the faint end studied here. We find that the faint end of the B-band luminosity function (-18 \leq M_B \leq -15) has a slope that is steeper than the Stromlo-APM estimate, but in rough agreement with the recent ESO Slice Project estimates.

We use the rotation curves of a sample of dark matter dominated dwarf and low-surface brightness (LSB) late-type galaxies to study their radial mass distributions. We find that the shape of the rotation curves is remarkably similar for all (both dwarf and LSB) galaxies in the sample, suggesting a self-similar density distribution of their dark matter (DM) halos. This shape can be reproduced well by a density profile with a shallow central cusp (rho(r)\propto 1/r^gamma, gamma\approx 0.2-0.4), corresponding to a steeply rising velocity curve (v(r)\propto r^g, g\approx 0.9-0.8). We further show that the observed shape of the rotation curves is well matched by the density profiles of DM halos formed in very high resolution simulations of the CDM, LCDM, and CHDM models of structure formation. This is surprising in light of several previous studies which suggested that the structure of simulated DM halos is inconsistent with the dynamics of dwarf galaxies. We discuss possible explanations for this discrepancy. We show that although the mass distribution in the hierarchically formed halos is on average consistent with the shape of rotation curves of dark matter dominated galaxies, the scatter of the individual profiles around the average is substantial and should not be neglected in comparisons with the data. Finally, we show that the DM halos in our hierarchical simulations and the real galaxies in our sample exhibit a similar decrease in their characteristic densities with increasing characteristic radial scales, and increase in their maximum rotation velocities with increasing radii at which their maximum velocity occurs. (Abridged)

We present two matched sets of five simulations each, covering five presently favored simple modifications to the standard cold dark matter (CDM) scenario. One simulation suite, with a linear box size of 75 Mpc/h, is designed for high resolution and good statistics on the group/poor cluster scale, and the other, with a box size of 300 Mpc/h, is designed for good rich cluster statistics. All runs had 57 million cold particles, and models with massive neutrinos had an additional 113 million hot particles. We consider separately models with massive neutrinos, tilt, curvature, and a nonzero cosmological constant in addition to the standard CDM model. We find that our tilted Omega+Omega_Lambda=1 (TLCDM) model produces too much small-scale power by a factor of ~3, and our open Lambda=0 (OCDM) model also exceeds observed small-scale power by a factor of 2. In addition, we take advantage of the large dynamic range in detectable halo masses our simulations allow to check the shape of the Press-Schechter approximation. We find good fits at cluster masses for delta_c=1.27--1.35 for a Gaussian filter and delta_c=1.57--1.73 for a tophat filter. However, Press-Schechter overpredicts the number density of halos compared to the simulations in the high resolution suite by a weakly cosmology-dependent factor of 1.5--2 at galaxy and group masses, which cannot be fixed by adjusting delta_c within reasonable bounds. An appendix generalizes the spherical collapse model to any isotropic cosmology.

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.

Particle-Mesh (PM) codes are still very useful tools for testing predictions of cosmological models in cases when extra high resolution is not very important. We release for public use a cosmological PM N-body code. We provide a complete package of routines needed to set initial conditions, to run the code, and to analyze the results. The package allows you to simulate models with numerous combinations of parameters: open/flat/closed background, with or without the cosmological constant, different values of the Hubble constant, with or without hot neutrinos, tilted or non-tilted initial spectra, different amount of baryons. Routines are included to measure the power spectrum and the density distribution function in your simulations, and a bound-density-maxima code for halo finding. We also provide results of test runs. A simulation with 256^3 mesh and 128^3 particles can be done in a couple of days on a typical workstation (70Mb of RAM are needed). To run simulations with 800^3 mesh and 256^3 particles one needs a computer with 1Gb memory and 1Gb disk space. The code has been successfully tested on an HP workstation and on a Sun workstation running Solaris. Most of the files (not tests) can be obtained from ftp://astro.nmsu.edu/pub/aklypin/PMCODE The package can be downloaded from http://astro.nmsu.edu/~aklypin/PM/pmcode/index.html We provide this tool as a service to the astronomical community, but we cannot guarantee results or publications.

We present a new high-resolution N-body algorithm for cosmological simulations. The algorithm employs a traditional particle-mesh technique on a cubic grid and successive multilevel relaxations on the finer meshes, introduced recursively in a fully adaptive manner in the regions where the density exceeds a predefined threshold. The mesh is generated to effectively match an arbitrary geometry of the underlying density field -- a property particularly important for cosmological simulations. In a simulation the mesh structure is not created at every time step but is properly adjusted to the evolving particle distribution. The algorithm is fast and effectively parallel. We present a detailed description of the methodology, implementation, and tests of the code. We further use the code to study the structure of dark matter halos in high-resolution (2/h kpc) simulations of standard CDM (Omega=1, h=0.5, sigma_8=0.63) and LCDM (Omega_Lambda=1-Omega_0=0.7, h=0.7, sigma_8=1.0) models. We find that halo density profiles in both CDM and LCDM models are well fitted by the analytical model presented recently by Navarro et al., which predicts a singular [$\rho(r)\propto r^{-1}$] behavior of the halo density profiles at small radii. We therefore conclude that halos formed in the $\Lambda$CDM model have structure similar to CDM halos and thus cannot explain the dynamics of the central parts of dwarf spiral galaxies, as inferred from the galaxies' rotation curves.

Reconstruction of the linear power spectrum from observational data provides a way to compare cosmological models to a large amount of data, as Peacock & Dodds (1994, 1996) have shown. By applying the appropriate corrections to the observational power spectrum it is possible to recover the underlying linear power spectrum for any cosmological model. Using extensive N-body simulations we demonstrate that the method is applicable to a wide range of cosmological models. However, we find that the recovery of the linear power spectrum from observations following PD94 is misleading because the corrections are model- dependent. When we apply the proper corrections for a given model to the observational power spectrum, we find that no model in our test group recovers the linear power spectrum well for the bias suggested by PD94 between Abell, Radio, Optical, and IRAS catalogs 4.5:1.9:1.3:1, with b_IRAS=1. When we allow b_IRAS to vary we find that: (i)CHDM models give very good fits to observations if optically-selected galaxies are slightly biased b_Opt=1.1 (ii) Most LCDM models give worse but acceptable fits if blue galaxies are considerably antibiased: 0.6<b_Opt<0.9 and fail if optical galaxies are biased. (iii)There is a universal shape of the recovered linear power spectrum of all LCDM models over their entire range of explored wavenumbers,0.01<k<0.6h\Mpc. Recovered spectra of CDM and CHDM models are nearly the same as that of LCDM in the region 0.01<k<0.2h/Mpc but diverge from this spectrum at higher k.

(Shortened Abstract). We analyze a redshift sample of Abell/ACO clusters and compare them with numerical simulations based on the truncated Zel'dovich approximation (TZA), for a list of eleven dark matter (DM) models. For each model we run several realizations, on which we estimate cosmic variance effects. We analyse correlation statistics, the probability density function, and supercluster properties from percolation analysis. As a general result, we find that the distribution of galaxy clusters provides a constraint only on the shape of the power spectrum, but not on its amplitude: a shape parameter 0.18 < \Gamma < 0.25 and an effective spectral index at 20Mpc/h in the range [-1.1,-0.9] are required by the Abell/ACO data. In order to obtain complementary constraints on the spectrum amplitude, we consider the cluster abundance as estimated using the Press--Schechter approach, whose reliability is explicitly tested against N--body simulations. We conclude that, of the cosmological models considered here, the only viable models are either Cold+Hot DM ones with \Omega_\nu = [0.2-0.3], better if shared between two massive neutrinos, and flat low-density CDM models with \Omega_0 = [0.3-0.5].

We apply several statistical estimators to high-resolution N-body simulations of two currently viable cosmological models: a mixed dark matter model, having $\Omega_\nu=0.2$ contributed by two massive neutrinos (C+2\nuDM), and a Cold Dark Matter model with Cosmological Constant (\LambdaCDM) with $\Omega_0=0.3$ and h=0.7. Our aim is to compare simulated galaxy samples with the Perseus-Pisces redshift survey (PPS). We consider the n-point correlation functions (n=2-4), the N-count probability functions P_N, including the void probability function P_0, and the underdensity probability function U_\epsilon (where \epsilon fixes the underdensity threshold in percentage of the average). We find that P_0 (for which PPS and CfA2 data agree) and P_1 distinguish efficiently between the models, while U_\epsilon is only marginally discriminatory. On the contrary, the reduced skewness and kurtosis are, respectively, S_3\simeq 2.2 and S_4\simeq 6-7 in all cases, quite independent of the scale, in agreement with hierarchical scaling predictions and estimates based on redshift surveys. Among our results, we emphasize the remarkable agreement between PPS data and C+2\nuDM in all the tests performed. In contrast, the above \LambdaCDM model has serious difficulties in reproducing observational data if galaxies and matter overdensities are related in a simple way.

A new class of geometric statistics for analyzing galaxy catalogs is presented. Filament statistics quantify filamentarity and planarity in large scale structure in a manner consistent with catalog visualizations. These statistics are based on sequences of spatial links which follow local high-density structures. From these link sequences we compute the discrete curvature, planarity, and torsion. Filament statistics are applied to CDM and CHDM ($\Omega_\nu = 0.3$) simulations of Klypin \etal (1996), the CfA1-like mock redshift catalogs of Nolthenius, Klypin and Primack (1994, 1996), and the CfA1 catalog. We also apply the moment-based shape statistics developed by Babul \& Starkman (1992), Luo \& Vishniac (1995), and Robinson \& Albrecht (1996) to these same catalogs, and compare their robustness and discriminatory power versus filament statistics. For 100 Mpc periodic simulation boxes ($H_0 = 50$ km s$^{-1}$ Mpc$^{-1}$), we find discrimination of $\sim 4\sigma$ (where $\sigma$ represents resampling errors) between CHDM and CDM for selected filament statistics and shape statistics, including variations in the galaxy identification scheme. Comparing the CfA1 data versus the models does not yield a conclusively favored model; no model is excluded at more than a $\sim 2\sigma$ level for any statistic, not including cosmic variance which could further degrade the discriminatory power. We find that CfA1 discriminates between models poorly mainly due to its sparseness and small number of galaxies, not due to redshift distortion, magnitude limiting, or geometrical effects. We anticipate that the proliferation of large redshift surveys and simulations will enable the statistics presented here to provide robust discrimination between large-scale structure in various cosmological models.

Here we discuss what are perhaps the two most popular variants of CDM that might agree with the data: \lcdm\ and CHDM. While the predictions of COBE-normalized \lcdm\ and CHDM both agree well with the available data on scales of $\sim 10$ to $100 \hMpc$, each has potential virtues and defects. \lcdm\ with $\Omega_0 \sim 0.3$ has the possible virtue of allowing a higher expansion rate $H_0$ for a given cosmic age $t_0$, but the defect of predicting too much fluctuation power on small scales. CHDM has less power on small scales, so its predictions appear to be in good agreement with data on the galaxy distribution, but it remains to be seen whether it predicts early enough galaxy formation to be compatible with the latest high-redshift data. Also, two very recent observational results favor high cosmic density, and thus favor $\Omega=1$ models such as CHDM over \lcdm\ --- (1) the positive deceleration parameter $q_0>0$ measured using high-redshift Type Ia supernovae, and (2) the low primordial deuterium/hydrogen ratio measured in two different quasar absorption spectra. We try to identify ``best'' variants of both \lcdm\ and CHDM, and discuss critical observational tests for both models.

We present the first X-ray Luminosity Function (XLF) for an optically-selected sample of 49 nearby poor clusters of galaxies and a sample of 67 Abell clusters with z < 0.15. We have extended the measured cluster XLF by more than a factor of 10 in X-ray luminosity. Our poor cluster sample was drawn from an optical catalog of groups with 0.01<z<0.03 composed of Zwicky galaxies. The X-ray emission was measured from the ROSAT all-sky survey. About 45% of the poor clusters were detected with 0.5-2.0 keV luminosities from (1.7-65)x 10^{41} h^{-2} ergs/sec. These are among the X-ray brightest, optically-selected poor clusters in the northern hemisphere. For this sample, the poor cluster XLF was found to be smooth extrapolation of the rich cluster XLF. A new Hydro/N-body simulation of a Hot +Cold dark matter model with \Omega_{total}=1, \Omega_{\nu}=0.2, and a baryon fraction of 7.5% was used to model and understand our observational selection effects. We found that the observed cluster Gas Mass Function was consistent with our model.

The purpose of these lectures is to give a short introduction into a very vast field of numerical simulations for cosmological applications. I focus on major features of the simulations: the equations, main numerical techniques, effects of resolution, and methods of halo identification.

We numerically simulate some of the most critical physical processes in galaxy formation: The supernova feedback, in conjunction with gasdynamics and gravity, plays a crucial role in determining how galaxies arise within the context of a model for large-scale structure. Our treatment incorporates a multi-phase model of the interstellar medium and includes the effects of cooling, heating and metal enrichment by supernovae, and evaporation of cold clouds. The star formation happens inside the clouds of cold gas, which are produced via thermal instability. We simulate the galaxy formation in standard biased CDM model for a variety of parameters and for several resolutions in the range 2--20$h^{-1}$kpc. In our picture, supernova feedback regulates the evolution of the gas components and star formation. The efficiency of cloud evaporation by supernova strongly influences star formation rates. This feedback results in a steady rate of star formation in large galaxies (mass larger than $2-3x10^{11}\Msun$) at a level of $(1-10)\Msun\yr$ for $z<3$. Supernova feedback has an even stronger effect on the evolution of dwarf galaxies, most of which have a small fraction of stars and extremely low luminosities: $M_R>-15$. In the case of both large and small galaxies, the distribution of luminous matter (stars) is strongly BIASED with respect to the dark matter. We find an approximate biasing measure of the form $\rho_{lum}= (\rho_{dm}/133)^{1.7}$ for z=0 and overdensities exceeding 1000. Deviations from this relation (a factor 2-3) depend on the environment. For halo masses exceeding $2x10^{10}\Msun$, the dependence of the absolute magnitude on the total mass can be approximated as $M_V=-18.5-4\log(M_{tot}/10^{11}\Msun)$, with a scatter of less than 0.5mag.

Cosmological models with a positive cosmological constant and $\Omega_0<1$ have a number of attractive features. A larger Hubble constant, which can be compatible with the recent HST estimate, and a large fraction of baryon density in galaxy clusters make them current favorites. Early galaxy formation also is considered as a welcome feature of these models But early galaxy formation implies that fluctuations on few megaparsec scales spent more time in the nonlinearregime, as compared with standard Cold Dark Matter (CDM) or Cold+Hot Dark Matter models. This results in excessive clustering on small scales. We show that a typical LCDM model with $H_0=70$ km/s/Mpc, $\Omega_0=0.3$ normalized to COBE on large scales and compatible with the number-density of galaxy clusters, predicts a power spectrum of galaxy clustering in real space which is too high: {\it at least} twice larger than CfA estimates (Park \etal 1994) and 3 times larger than APM estimates (Baugh \& Efstathiou 1994) for wavenumbers $k=(0.4-1)h/{\rm Mpc}$. This conclusion holds if we assume either that galaxies trace the dark matter ($\sigma_8\approx 1.1$ for this model) or just that a region with higher density produces more galaxies than a region with lower density. The only way to reconcile the model with the observed power spectrum is to assume that regions with high dark matter density produce fewer galaxies than regions with low density. Theoretically this is possible, but it seems very unlikely: X-ray emission from groups and clusters indicates that places with a large density of dark matter produce a large number of galaxies. Since it follows that the low- $\Omega$ LCDM models are in serious trouble, we discuss which LCDM models have the best hope of surviving the confrontation with available observational data.

We discuss the structure of clusters in a class of flat cosmological models with the fraction of mass \Omega_{CDM} ~0.8 in cold dark matter, and the rest in hot dark matter in the form of massive neutrinos. We consider such Cold+Hot Dark Matter (CHDM) models with one, two or three massive neutrinos, with total mass ~4.6eV. Neutrinos of such low mass cannot constitute halos of galaxies and groups, but only of clusters of galaxies. The limit on the density of neutrinos in the central parts of galaxy clusters is estimated from the phase space den- sity constraints. The ratio of the neutrino density to that of CDM through the cluster is found analytically. It appears that the density of neutrinos is suppressed within the Abell radius. However, neutrinos contribute ~20% of the mass density to the cluster halo. Our numerical simulations match analytical results. The simulations indicate that the cluster halo dark matter density profile has the power-law slope ~-2.5 which is close to that in the model with cosmological constant. We also found that in the CHDM models the velocity dispersion is almost constant across the cluster. This is quite different from the model with cosmological constant or the open model where the velocity dispersion falls in the cluster outskirts. We discuss X-ray emission and weak gravitational lensing by clusters in the model. We input the found spherically symmetrical fit to the CHDM mass density profile and the X-ray surface brightness for the cluster A2256 into simple equation of hydrostatic equilibrium of the hot gas. X-ray temperature derived this way departures from both the data and actual prediction of the model, which give almost constant temperature. We found also that the problem of high baryonic fraction in clusters is not resolved in the CHDM models.

We use N-body simulations to study evolution of galaxy clusters over the redshift interval 0 <= z <= 0.5 in cosmological models with a mixture of cold and hot dark matter (CHDM). Four different techniques are utilized: the cluster-cluster correlation function, axial ratios and quadrupoles of the dark matter distribution in individual clusters, and virial properties. We find that the correlation function for clusters of the same mass limit was larger and steeper at high redshifts. The slope increases from 1.8 at z=0 to 2.1 at z=0.5. Comoving correlation length r_c scales with the mass limit M within comoving radius 1.5 h^-1 Mpc and the redshift z as r_c ~= 20(1+z)(M/M_*)^1/3, where M_* = 3*10^14 h^-1 M_sun. When the correlation length is normalized to the mean cluster separation d_c, it remains almost constant: r_c~=(0.45-0.5) d_c. For small masses (M_clust <2*10^14 h^-1 M_sun) there is an indication that r_c goes slightly above the relation with the constant of proportionality being ~= 0.55-0.6. Anisotropy of density distribution in a cluster shows no change over redshift with axial ratios remaining constant around 1.2. In other words, clusters at present are as elongated as they were at the epoch of their first appearance. While the anisotropy of clusters does not change with time, the density profile shows visible evolution: the slope of density profile changes from gamma ~= -3.5 at z=0.5 to gamma ~= -2.5 at the present. We find that the core of a cluster remains essentially the same over time, but the density of the outlying regions increases noticeably. The virial relation M ~ v^2 is a good approximation, but there is a large fraction of clusters with peculiar velocities greater than given by this relation, and clusters with the same rms velocities have smaller masses in the past, a factor of 2 at z=0.5.

A new class of geometric statistics for analyzing galaxy catalogs is presented. {\it Filament statistics} quantify filamentarity and planarity in large scale structure in a manner consistent with catalog visualizations. These statistics are based on sequences of spatial links which follow local high-density structures. From these link sequences we compute the discrete curvature, planarity, and torsion. Filament statistics are applied to CDM and CHDM ($\Omega_\nu = 0.3$) simulations of Klypin \etal (1995), the CfA1-like sky catalogs of Nolthenius, Klypin and Primack (1994, 1995), and the CfA1 catalog. For 100 Mpc periodic simulation boxes ($H_0 = 50$ km s$^{-1}$ Mpc$^{-1}$), we find robust discrimination of over 4$\sigma$ (where $\sigma$ represents resampling errors) between CHDM and CDM. The {\it reduced filament statistics} show that CfA1 data is intermediate between CHDM and CDM, but more consistent with the CHDM models. Filament statistics provide robust and discriminatory shape statistics with which to test cosmological simulations of various models against present and future redshift survey data.

This paper presents an angular correlation analysis of the Green Bank 4.85 GHz radio catalog (Gregory \& Condon 1991) of 54,579 sources (S \gsim 25 mJy). The Green Bank catalog is found to be complete to S $\geq$ 35 mJy over 20$^{\circ} \leq \delta < 74^{\circ}$, 0$^h \leq \alpha < 24^h$, and Galactic latitude $|b| \geq 10^{\circ}$. The 2-point angular correlation function shows evidence for the clustering of radio sources, with a power-law distribution consistent with a slope $\gamma = -0.8$. This may well provide the $first$ detection of an angular correlation in a large area, complete deep radio survey.

We present a series of four simulations of Cold Dark Matter (CDM) and Cold + Hot Dark Matter (CHDM) cosmologies. We discuss the power spectrum and correlation functions in real and redshift space, with comparisons to the CfA2 and IRAS redshift data, the pairwise velocity of galaxies, and the distribution of hot and cold particles in CHDM simulations. We confirm that CHDM with cold/hot/baryon density ratios $\Omega_c/\Omega_\nu/\Omega_b= 0.6/0.3/0.1$ is a good fit to a wide variety of present-epoch data, much better than CDM. In particular, with reasonable assumptions about identification of galaxies and biasing, we find that the power spectrum agrees rather well with both the CfA2 and IRAS power specta in both the nonlinear and linear regimes. New CHDM models (e.g., with $\Omega_\nu=0.2$ or with two massive neutrinos) predict a significantly larger rate of formation of galaxies at high redshift, which may be needed to explain some observational data. At the same time, the difference between the variants is rather small at $z=0$. The results presented in this paper are interesting for two purposes: (i) For a comparison with other classes of models (like CDM and $\Lambda$CDM) at $z=0$.(ii) As a reference point for comparison between different variants of the CHDM model.

We use the moments of counts of neighbors as given by the Generalized Correlation Integrals, to study the clustering properties of Dark Matter Halos (DH) in Cold Dark Matter (CDM) and Cold+Hot Dark Matter (CHDM) models. We compare the results with those found in the CfA and SSRS galaxy catalogs. We show that if we apply the analysis in redshift space, both models reproduce equally well the observed clustering of galaxies. Mass segregation is also found in the models: more massive DHs are more clustered compared with less massive ones. In redshift space, this mass segregation is reduced by a factor 2-3 due to the peculiar velocities. Observational catalogs give an indication of luminosity and size segregation, which is consistent with the predictions of the models. Because the mass segregation is smaller in redshift space, it is suggestive that the real luminosity or size segregation of galaxies could be significantly larger than what it is found in redshift catalogs.

Cold + Hot Dark Matter (CHDM) $\Omega=1$ cosmological models require a total neutrino mass $\sim 5$ eV. Because recent data support the $\nu_\mu \to \nu_\tau$ oscillation explanation of the cosmic ray $\nu_\mu$ deficit, which requires that $m(\nu_\mu) \approx m(\nu_\tau)$, this suggests that $m(\nu_\mu) \approx m(\nu_\tau) \approx 2.4$ eV. The linear calculations and N-body simulation reported here indicate that an $\Omega=1$ CHDM model with two 2.4 eV neutrinos (designated C$\nu^2$DM) agrees remarkably well with all available observations, but only if the Hubble parameter $h \approx 0.5$. We also show that even one 2.4 eV neutrino raises serious difficulties for low-$\Omega$ flat CDM models.

We present techniques for identifying and analyzing galaxy groups and show that they provide a powerful discriminator between cosmological models. We apply these to high-resolution PM N-body simulations of structure formation in CHDM (cold/hot/baryon fractions=.6/.3/.1, b=1.5 (COBE norm) and two CDM models; b=1.5 and b=1.0 (COBE norm). Groups are identified with the Nolthenius (1993) algorithm. We find: (1) Properties of groups are a powerful and robust discriminator between these Gaussian models whose spatial and velocity properties differ on Mpc scales. We test robustness against several methods for assigning luminosity to dark matter halos, for merging CfA1 galaxies, and for breaking up massive over-mergers (2) When allowance is made for the high large scale power P(k) present in the CfA1 data, CHDM at our $\Omnu=0.30$ produces slightly too many groups and groups too high a fraction of the galaxies, while the fraction grouped in CDM is far too low. Lowering $\Omnu$ would appear to produce close agreement with all measures, save one: for all simulations, median group sizes are x1.4-1.7 larger than equivalently selected CfA1 groups. (3) The standard group M/L method gives $\Omega\simeq 0.1$ for CHDM and CfA1, and $\Omega\simeq 0.35$ for CDM. A conspiracy of three effects produce these very low values, and we conclude that low observed $\Omega$'s need not argue for a low $\Omega$ universe. When overmergers are broken up, the median virial- to-DM mass of 3D selected groups is $\sim1$ for all simulations. $M_{DM}>10^{14}\Msun$ groups appear virialized in all simulations. Global velocity biases $b_v$ are similar to previous studies. Within 3D-selected groups, CHDM and CDM b=1.5 show a stronger bias of $b_v=0.7-0.8$, while CDM b=1.0 shows group $b_v$'s $\simeq 1$.

Damped Ly$\alpha$ systems provide possibly the most significant evidence for early structure formation, and thus a stringent constraint on the Cold + Hot Dark Matter (CHDM) cosmology. Using the numbers of halos in N-body simulations to normalize Press-Schechter (PS) estimates of the number densities of protogalaxies as a function of redshift, we find that CHDM with $\Omega_c/\Omega_\nu/\Omega_b = 0.6/0.3/0.1$ is compatible with the damped Ly$\alpha$ data only at $z<2.5$, but that it is probably incompatible with data at $z>3$. The predictions of CHDM are quite sensitive to the neutrino fraction. We find that $\Omega_c/\Omega_\nu/\Omega_b = 0.725/0.20/0.075$ is compatible with the $z>3$ data. With one massive neutrino species, this corresponds to lowering the neutrino mass from 7.0 to 4.7 eV, for $H_0=50\kmsMpc$ and $T=2.726$ K. By analysing our numerical simulations with different resolutions and box sizes as well as those of Ma \& Bertchinger (1994), we show that for the CHDM models with $\Omega_\nu$=0.2--0.3 the PS approximation should be used with Gaussian filter with $\delta_c=1.3-1.4$ if one tries to recover the total mass of a collapsed halo and to include nonlinear effects, due to waves both longer and shorter than those within the simulation box.

We present the first self-consistently computed results on the evolution of X-ray properties of galaxy clusters in a Cold + Hot Dark Matter (CHDM) model. We have performed a hydrodynamic plus N-body simulation for the COBE-compatible CHDM model with standard mass components: Omega(hot) = 0.3, Omega(cold) = 0.6 and Omega(baryon) = 0.1 (h = 0.5). In contrast with the CDM model, which fails to reproduce the observed temperature distribution function dN/dT (Bryan et al. 1994b), the CHDM model fits the observational dN/dT quite well. Our results on X-ray luminosity are less firm but even more intriguing. We find that the resulting X-ray luminosity functions at redshifts z = 0.0, 0.2, 0.4, 0.7 are well fit by observations, where they overlap. The fact that both temperatures and luminosities provide a reasonable fit to the available observational data indicates that, unless we are missing some essential physics, there is neither room nor need for a large fraction of gas in rich clusters: 10% (or less) in baryons is sufficient to explain their X-ray properties. We also see a tight correlation between X-ray luminosity and gas temperature.

We use the void probability statistics to study the redshift-space galaxy distribution as described by a volume-limited subsample of the Perseus-Pisces survey. We compare the results with the same analysis realized on artificial samples, extracted from high-resolution N-body simulations by reproducing the observational biases of the real data set. Simulations are run for the Cold+HotDM model (CHDM) and for unbiased and biased (b=1.5) CDM models in a 50 Mpc/h box. We identify galaxies as residing in peaks of the evolved density field. We fragment overmerged structures into individual galaxies so as to reproduce both the correct luminosity function (after assuming M/ L values for the resulting galaxy groups) and the two-point correlation function. Our main result is that a void-probability function (VPF) from the standard CHDM model with fractions 60% cold, 30% hot, 10% barions, exceeds the observational VPF with a high confidence level. CDM models produce smaller VPF independent of the biasing parameter. We verify the robustness of this result against changing the observer position in the simulations and the galaxy identification in the evolved density field.

We estimate the angular correlation function for the standard CDM, tilted $n=0.7$ CDM and hybrid (CHDM) models, and compare with observations. When compared with the APM observational results scaled to the Lick depth, there appears to be fair agreement with the estimate from the CHDM model. But a more detailed comparison using the unscaled APM data for the five magnitude slices with $\Delta b_j=0.5$ shows that, in fact, none of the models can actually fit w(theta) for all the slices simultaneously. As $n=0.7$ tilted CDM and SCDM bracket all possible tilted models, we conclude that none of the tilted models is consistent with the APM results. With CHDM, a \wtheta amplitude that is 30\%--40\% too high on scales $\theta < 0.\deg 5$ is predicted for the deepest slices ($b_j\approx 20$). We find that no reasonable simple variation of the parameters for the luminosity function or for the evolution of the correlation function with redshift could change the situation. Thus, the discrepancy between the APM data and the CHDM model, though small, seems to be real.

We discuss the effects of redshift--space distortions on the estimate of high order galaxy correlation functions. We use both the Perseus--Pisces redshift survey and the results of high--resolution N--body simulations to explore the consequences of working in redshift space on the detection of deviations from hierarchical clustering. Both for real and simulated data, significant deviations from hierarchical clustering seem to be present in real space. Their behaviour is coherent with that expected from an initially biased galaxy field, once the displacement from the sites where galaxies formed, due to the nonlinear gravitational evolution, is taken into account. The passage to redshift space has the net effect of filtering out the higher powers required to fit the distribution in real space. This magnifies the distortions operated by nonlinear evolution on the initial distribution, erasing the residual narrow scale range (2--5~$h^{-1}$Mpc) where deviations from hierarchical clustering are still detectable in real space. We conclude that such deviations can hardly be estimated using data in redshift space.

Although the Cold + Hot Dark Matter (CHDM) cosmology provides perhaps the best fit of any model to all the available data at the current epoch, CHDM produces structure at relatively low redshifts and thus could be ruled out if there were evidence for formation of massive objects at high redshifts. Damped Ly$\alpha$ systems are abundant in quasar absorption spectra and thus provide possibly the most significant evidence for early structure formation, and thus perhaps the most stringent constraint on CHDM. Using the numbers of halos in N-body simulations to normalize Press-Schechter estimates of the number densities of protogalaxies as a function of redshift, we find that CHDM with $\Omega_c/\Omega_\nu/\Omega_b = 0.6/0.3/0.1$ is compatible with the damped Ly$\alpha$ data at $\le 2.5$, but that it is probably incompatible with the limited $z>3$ damped Ly$\alpha$ data. The situation is uncertain because there is very little data for $z>3$, and also it is unclear whether all damped Ly$\alpha$ systems are associated with collapsed protogalaxies. The predictions of CHDM are quite sensitive to the hot (neutrino) fraction, and we find that $\Omega_c/\Omega_\nu/\Omega_b = 0.675/0.25/0.075$ is compatible even with the $z>3$ data. This corresponds to lowering the neutrino mass from 6.8 to 5.7 eV, for $H_0=50\kmsMpc$. In CHDM, the higher redshift damped Ly$\alpha$ systems are predicted to have lower masses, which can be checked by measuring the velocity widths of the associated metal line systems.

This letter presents results of new high resolution $\Omega=1$ Cold + Hot Dark Matter (CHDM) and Cold Dark Matter (CDM) simulations. Properties of groups in these simulations reflect the lower small-scale velocities and the greater tendency to form distinct filaments on both small and large scales in CHDM as compared to CDM. The fraction of galaxies in groups and the median group rms velocity are found to be powerful discriminators between models. We combine these two features into a very robust statistic, median group rms velocity $v_{\rm gr}(f_{\rm gr})$ as a function of the fraction $f_{\rm gr}$ of galaxies in groups. Using this statistic, we compare ``observed'' simulations to CfA data in redshift space in a careful and consistent way. We find that CHDM remains a promising model, with for example $v_{\rmgr}(0.45) \approx 125 \pm 25 \kms$ in agreement with the CfA data, while CDM with bias b=1.0 (COBE-compatible) or b=1.5, both giving $v_{\rm gr}(0.45) \approx 400 \pm 25 \kms$, can be virtually ruled out. Using median $M/L$, the observed value of $\Omega$ is $0.10$ (CHDM) to $0.38$ (CDM).

We use high resolution PM N-body simulations to follow the development of non-linear clustering in a flat Universe, dominated by Cold + Hot Dark Matter (CHDM) with 60% of CDM, 30% of HDM and 10% of baryons; a simulation box of 100 Mpc a side ($h=0.5$) is used. We analyze two CHDM simulations with $b =1.5$ (COBE normalization). We also compare them with CDM simulations with $b=1.5$ and $b=1$. We evaluate high-order correlation functions and the void-probability-function (VPF). Correlation functions are obtained both from counts in cells and counts of neighbors. The analysis is made for DM particles and for galaxies, identified as massive halos in the evolved density field. We also check the effects of dynamical evolution and redshift space distortions. We find that clustering of DM particles exhibits deviations from the hierarchical scaling, which decrease somewhat in redshift space. Galaxies follow hierarchical scaling far more closely, with coefficients S_3=2.5 and S_4=7.5, in general agreement with observational results, just marginally affected by redshift distortions and on the choices for the initial spectra. The hierarchical scaling of galaxy clustering is confirmed by the VPF analysis and in all the cases considered a good agreement with observational results is obtained.

We propose the zero-point of the cluster-cluster correlation function as a sensitive test for the shape of the power spectrum of initial fluctuations. It is now possible to go beyond the power law description to measure the point at which the correlation function becomes zero. Four independent measurements indicate that the zero point, $r_0$, should be in the range (40-60)$h^{-1}$Mpc. The large value of $r_0$ at which the zero-point occurs rules out conventional CDM models independently of the assumed amplitude. Models with $\Omega < 0.25$ should be rejected because they predict too large $r_0$. If the age of the Universe is assumed to be larger than 15 Gyr, models with either $\Omega<0.5$ or $h> 0.55$ are rejected. We present the results of numerical simulations of clusters in an $\Omega=1$ cosmological model with a mixture of cold plus hot dark matter (CHDM). The correlation function we determined for the simulated clusters has a zero-point, $r_0=55h^{-1}$Mpc that accurately matches the zero point of the observed function.

We report results from high-resolution particle-mesh (PM) N-body simulations of structure formation in an $\Omega=1$ cosmological model with a mixture of Cold plus Hot Dark Matter (C+HDM) having $\Omega_{\rm cold}=0.6$, $\Omega_\nu=0.3$, and $\Omega_{\rm baryon}=0.1$. We present analytic fits to the C+HDM power spectra for both cold and hot ($\nu$) components, which provide initial conditions for our nonlinear simulations. In order to sample the neutrino velocities adequately, these simulations included six times as many neutrino particles as cold particles. Our simulation boxes were 14, 50, and 200~Mpc cubes (with $H_0=50$ km s$^{-1}$ Mpc$^{-1}$); we also did comparison simulations for Cold Dark Matter (CDM) in a 50~Mpc box. C+HDM with linear bias factor $b=1.5$ is consistent both with the COBE data and with the galaxy correlations we calculate. We find the number of halos as a function of mass and redshift in our simulations; our results for both CDM and C+HDM are well fit by a Press-Schechter model. The number density of galaxy-mass halos is smaller than for CDM, especially at redshift $z>2$, but the numbers of cluster-mass halos are comparable. We also find that on galaxy scales the neutrino velocities and flatter power spectrum in C+HDM result in galaxy pairwise velocities that are in good agreement with the data, and about 30\% smaller than in CDM with the same biasing factor. On scales of several tens of Mpc, the C+HDM streaming velocities are considerably larger than CDM. Thus C+HDM looks promising as a model of structure formation.