We examine the thermodynamic state and cooling of the low-$z$ Circum-Galactic Medium (CGM) in five FIRE-2 galaxy formation simulations of Milky Way-mass galaxies. We find that the CGM in these simulations is generally multiphase and dynamic, with a wide spectrum of largely nonlinear density perturbations sourced by the accretion of gas from the IGM and winds from both the central and satellite galaxies. We investigate the origin of the multiphase structure of the CGM with a particle tracking analysis and find that most of the low entropy gas has cooled from the hot halo as a result of thermal instability triggered by these perturbations. The ratio of cooling to free-fall timescales $t_{\rm cool}/t_{\rm ff}$ in the hot component of the CGM spans a wide range $\sim 1-100$ at a given radius, but exhibits approximately constant median values $\sim 5-20$ at all radii $0.1 R_{\rm vir} < r < R_{\rm vir}$. These are similar to the $\approx 10-20$ value typically adopted as the thermal instability threshold in "precipitation" models of the ICM. Consequently, a one-dimensional model based on the assumption of a constant $t_{\rm cool}/t_{\rm ff}$ and hydrostatic equilibrium approximately reproduces the simulation number density and entropy profiles, but only if it assumes the metallicity profiles taken directly from the simulations. We explicitly show that the $t_{\rm cool}/t_{\rm ff}$ value of a gas parcel in the hot component of the CGM does not predict its probability of cooling and subsequently accreting onto the central galaxy. This suggests that the value of $t_{\rm cool}/t_{\rm ff}$ is a poor predictor of thermal stability in gaseous halos in which large-amplitude density perturbations are prevalent.

The population of Milky Way (MW) satellites contains the faintest known galaxies and thus provides essential insight into galaxy formation and dark matter microphysics. Here we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high Galactic latitude sky. In particular, we use cosmological zoom-in simulations of MW-like halos that include realistic Large Magellanic Cloud (LMC) analogs to fit the position-dependent MW satellite luminosity function. We report decisive evidence for the statistical impact of the LMC on the MW satellite population due to an estimated $6\pm 2$ observed LMC-associated satellites, consistent with the number of LMC satellites inferred from Gaia proper-motion measurements, confirming the predictions of cold dark matter models for the existence of satellites within satellite halos. Moreover, we infer that the LMC fell into the MW within the last $2\ \rm{Gyr}$ at high confidence. Based on our detailed full-sky modeling, we find that the faintest observed satellites inhabit halos with peak virial masses below $3.2\times 10^{8}\ M_{\rm{\odot}}$ at $95\%$ confidence, and we place the first robust constraints on the fraction of halos that host galaxies in this regime. We predict that the faintest detectable satellites occupy halos with peak virial masses above $10^{6}\ M_{\rm{\odot}}$, highlighting the potential for powerful galaxy formation and dark matter constraints from future dwarf galaxy searches.

We quantify the shape and overdensity of the galaxy distribution in the `Local Sheet' within a sphere of $R=8$ Mpc, and compare these properties with the expectations of the $\Lambda$CDM model. We measure ellipsoidal axis ratios of $c/a\approx0.16$ and $b/a\approx0.79$, indicating that the distribution of galaxies in the Local Volume can be approximated by a flattened oblate ellipsoid, consistent with the `sheet'-like configuration noted in previous studies. In contrast with previous estimates that the Local Sheet has a density close to average, we find that the number density of faint and bright galaxies in the Local Volume is $\approx1.7$ and $\approx5.2$ times denser, respectively, than the mean number density of galaxies of the same luminosity. Comparison with simulations shows that the number density contrasts of bright and faint galaxies within $8$ Mpc alone make the Local Volume a $\approx 2.5\sigma$ outlier in the $\Lambda$CDM cosmology. Our results indicate that the cosmic neighbourhood of the Milky Way may be unusual for galaxies of similar luminosity. The impact of the peculiar properties of our neighbourhood on the properties of the Milky Way and other nearby galaxies is not yet understood and warrants further study.

We present the analysis of deep optical imaging of the galaxy cluster Abell 133 with the IMACS instrument on Magellan. Our multi-band photometry enables stellar mass measurements in the cluster member galaxies down to a mass limit of $M_\star=3\times10^8\,M_\odot$ ($\approx 0.1$ of the Large Magellanic Cloud stellar mass). We observe a clear difference in the spatial distribution of large and dwarf galaxies within the cluster. Modeling these galaxies populations separately, we can confidently track the distribution of stellar mass locked in the galaxies to the cluster's virial radius. The extended envelope of the cluster's brightest galaxy can be tracked to $\sim 200$ kpc. The central galaxy contributes $\sim 1/3$ of the the total cluster stellar mass within the radius $r_{500c}$.

Non-gravitational interactions between dark matter particles with strong scattering, but relatively small annihilation and dissipation, has been proposed to match various observables on cluster and group scales. In this paper, we present the results from large cosmological simulations which include the effects of different self-interaction scenarios. In particular we explore a model with the differential cross section that can depend on both the relative velocity of the interacting particles and the angle of scattering. We focus on how quantities, such as the stacked density profiles, subhalo counts and the splashback radius change as a function of different forms of self-interaction. We find that self-interactions not only affect the central region of the cluster, the effect well known from previous studies, but also significantly alter the distribution of subhalos and the density of particles out to the splashback radius. Our results suggest that current weak lensing data can already put constraints on the self-interaction cross-section that are only slightly weaker than the Bullet Cluster constraints ($\sigma/m \lesssim 2$ cm$^2/$g), and future lensing surveys should be able to tighten them even further making halo profiles on cluster scales a competitive probe for DM physics.

We present a detection of the splashback feature around galaxy clusters selected using their Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, $r_{\rm sp}$, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that $r_{\rm sp}$ inferred from the observed radial distribution of galaxies is affected by selection effects related to the optical cluster-finding algorithms. We test this possibility by measuring the splashback feature in clusters selected via the SZ effect in data from the South Pole Telescope SZ survey and the Atacama Cosmology Telescope Polarimeter survey. The measurement is accomplished by correlating these clusters with galaxies detected in the Dark Energy Survey Year 3 data. The SZ observable used to select clusters in this analysis is expected to have a tighter correlation with halo mass and to be more immune to projection effects and aperture-induced biases than optically selected clusters. We find that the measured $r_{\rm sp}$ for SZ-selected clusters is consistent with the expectations from simulations, although the small number of SZ-selected clusters makes a precise comparison difficult. In agreement with previous work, when using optically selected redMaPPer clusters, $r_{\rm sp}$ is $\sim$ $2\sigma$ smaller than in the simulations. These results motivate detailed investigations of selection biases in optically selected cluster catalogs and exploration of the splashback feature around larger samples of SZ-selected clusters. Additionally, we investigate trends in the galaxy profile and splashback feature as a function of galaxy color, finding that blue galaxies have profiles close to a power law with no discernible splashback feature, which is consistent with them being on their first infall into the cluster.

The majority of baryons reside beyond the optical extent of a galaxy in the circumgalactic and intergalactic media (CGM/IGM). Gaseous halos are inextricably linked to the appearance of their host galaxies through a complex story of accretion, feedback, and continual recycling. The energetic processes, which define the state of gas in the CGM, are the same ones that 1) regulate stellar growth so that it is not over-efficient, and 2) create the diversity of today's galaxy colors, SFRs, and morphologies spanning Hubble's Tuning Fork Diagram. They work in concert to set the speed of growth on the star-forming Main Sequence, transform a galaxy across the Green Valley, and maintain a galaxy's quenched appearance on the Red Sequence. Most baryons in halos more massive than 10^12 Msolar along with their high-energy physics and dynamics remain invisible because that gas is heated above the UV ionization states. We argue that information on many of the essential drivers of galaxy evolution is primarily contained in this "missing" hot gas phase. Completing the picture of galaxy formation requires uncovering the physical mechanisms behind stellar and SMBH feedback driving mass, metals, and energy into the CGM. By opening galactic hot halos to new wavebands, we not only obtain fossil imprints of >13 Gyrs of evolution, but observe on-going hot-mode accretion, the deposition of superwind outflows into the CGM, and the re-arrangement of baryons by SMBH feedback. A description of the flows of mass, metals, and energy will only be complete by observing the thermodynamic states, chemical compositions, structure, and dynamics of T>=10^6 K halos. These measurements are uniquely possible with a next-generation X-ray observatory if it provides the sensitivity to detect faint CGM emission, spectroscopic power to measure absorption lines and gas motions, and high spatial resolution to resolve structures.

We use a statistical sample of galaxy clusters from a large cosmological $N$-body$+$hydrodynamics simulation to examine the relation between morphology, or shape, of the X-ray emitting intracluster medium (ICM) and the mass accretion history of the galaxy clusters. We find that the mass accretion rate (MAR) of a cluster is correlated with the ellipticity of the ICM. The correlation is largely driven by material accreted in the last $\sim 4.5$~Gyr, indicating a characteristic time-scale for relaxation of cluster gas. Furthermore, we find that the ellipticity of the outer regions ($R\sim R_{\rm 500c}$) of the ICM is correlated with the overall MAR of clusters, while ellipticity of the inner regions ($\lesssim 0.5 R_{\rm 500c}$) is sensitive to recent major mergers with mass ratios of $\geq 1:3$. Finally, we examine the impact of variations in cluster mass accretion history on the X-ray observable-mass scaling relations. We show that there is a {\it continuous\/} anti-correlation between the residuals in the $T_x-M$ relation and cluster MARs, within which merging and relaxed clusters occupy extremes of the distribution rather than form two peaks in a bi-modal distribution, as was often assumed previously. Our results indicate the systematic uncertainties in the X-ray observable-mass relations can be mitigated by using the information encoded in the apparent ICM ellipticity.

We present a detailed analysis of the physical processes that cause halo assembly bias -- the dependence of halo clustering on proxies of halo formation time. We focus on the origin of assembly bias in the mass range corresponding to the hosts of typical galaxies and use halo concentration as our chief proxy of halo formation time. We also repeat our key analyses across a broad range of halo masses and for alternative formation time definitions. We show that splashback subhaloes are responsible for two thirds of the assembly bias signal, but do not account for the entire effect. After splashback subhaloes have been removed, we find that the remaining assembly bias signal is due to a relatively small fraction ($\lesssim 10\%$) of haloes in dense regions. We test a number of additional physical processes thought to contribute to assembly bias and demonstrate that the two key processes are the slowing of mass growth by large-scale tidal fields and by the high velocities of ambient matter in sheets and filaments. We also rule out several other proposed physical causes of halo assembly bias. Based on our results, we argue that there are three processes that contribute to assembly bias of low-mass halos: large-scale tidal fields, gravitational heating due to the collapse of large-scale structures, and splashback subhaloes located outside the virial radius.

We present a suite of 18 synthetic sky catalogs designed to support science analysis of galaxies in the Dark Energy Survey Year 1 (DES Y1) data. For each catalog, we use a computationally efficient empirical approach, ADDGALS, to embed galaxies within light-cone outputs of three dark matter simulations that resolve halos with masses above ~5x10^12 h^-1 m_sun at z <= 0.32 and 10^13 h^-1 m_sun at z~2. The embedding method is tuned to match the observed evolution of galaxy counts at different luminosities as well as the spatial clustering of the galaxy population. Galaxies are lensed by matter along the line of sight --- including magnification, shear, and multiple images --- using CALCLENS, an algorithm that calculates shear with 0.42 arcmin resolution at galaxy positions in the full catalog. The catalogs presented here, each with the same LCDM cosmology (denoted Buzzard), contain on average 820 million galaxies over an area of 1120 square degrees with positions, magnitudes, shapes, photometric errors, and photometric redshift estimates. We show that the weak-lensing shear catalog, redMaGiC galaxy catalogs and redMaPPer cluster catalogs provide plausible realizations of the same catalogs in the DES Y1 data by comparing their magnitude, color and redshift distributions, angular clustering, and mass-observable relations, making them useful for testing analyses that use these samples. We make public the galaxy samples appropriate for the DES Y1 data, as well as the data vectors used for cosmology analyses on these simulations.

The surface densities of molecular gas, $\Sigma_{\rm H_2}$, and the star formation rate (SFR), $\dot\Sigma_\star$, correlate almost linearly on kiloparsec scales in observed star-forming (non-starburst) galaxies. We explore the origin of the linear slope of this correlation using a suite of isolated $L_\star$ galaxy simulations. We show that in simulations with efficient feedback, the slope of the $\dot\Sigma_\star$-$\Sigma_{\rm H_2}$ relation on kiloparsec scales is insensitive to the slope of the $\dot\rho_\star$-$\rho$ relation assumed at the resolution scale. We also find that the slope on kiloparsec scales depends on the criteria used to identify star-forming gas, with a linear slope arising in simulations that identify star-forming gas using a virial parameter threshold. This behavior can be understood using a simple theoretical model based on conservation of interstellar gas mass as the gas cycles between atomic, molecular, and star-forming states under the influence of feedback and dynamical processes. In particular, we show that the linear slope emerges when feedback efficiently regulates and stirs the evolution of dense, molecular gas. We show that the model also provides insights into the likely origin of the relation between the SFR and molecular gas in real galaxies on different scales.

We use mock galaxy survey simulations designed to resemble the Dark Energy Survey Year 1 (DES Y1) data to validate and inform cosmological parameter estimation. When similar analysis tools are applied to both simulations and real survey data, they provide powerful validation tests of the DES Y1 cosmological analyses presented in companion papers. We use two suites of galaxy simulations produced using different methods, which therefore provide independent tests of our cosmological parameter inference. The cosmological analysis we aim to validate is presented in DES Collaboration et al. (2017) and uses angular two-point correlation functions of galaxy number counts and weak lensing shear, as well as their cross-correlation, in multiple redshift bins. While our constraints depend on the specific set of simulated realisations available, for both suites of simulations we find that the input cosmology is consistent with the combined constraints from multiple simulated DES Y1 realizations in the $\Omega_m-\sigma_8$ plane. For one of the suites, we are able to show with high confidence that any biases in the inferred $S_8=\sigma_8(\Omega_m/0.3)^{0.5}$ and $\Omega_m$ are smaller than the DES Y1 $1-\sigma$ uncertainties. For the other suite, for which we have fewer realizations, we are unable to be this conclusive; we infer a roughly 70% probability that systematic biases in the recovered $\Omega_m$ and $S_8$ are sub-dominant to the DES Y1 uncertainty. As cosmological analyses of this kind become increasingly more precise, validation of parameter inference using survey simulations will be essential to demonstrate robustness.

Splashback refers to the process of matter that is accreting onto a dark matter halo reaching its first orbital apocenter and turning around in its orbit. The cluster-centric radius at which this process occurs, r_sp, defines a halo boundary that is connected to the dynamics of the cluster. A rapid decline in the halo profile is expected near r_sp. We measure the galaxy number density and weak lensing mass profiles around redMaPPer galaxy clusters in the first year Dark Energy Survey (DES) data. For a cluster sample with mean M_200m mass ~2.5 x 10^14 M_sun, we find strong evidence of a splashback-like steepening of the galaxy density profile and measure r_sp=1.13 +/- 0.07 Mpc/h, consistent with earlier SDSS measurements of More et al. (2016) and Baxter et al. (2017). Moreover, our weak lensing measurement demonstrates for the first time the existence of a splashback-like steepening of the matter profile of galaxy clusters. We measure r_sp=1.34 +/- 0.21 Mpc/h from the weak lensing data, in good agreement with our galaxy density measurements. For different cluster and galaxy samples, we find that consistent with LCDM simulations, r_sp scales with R_200m and does not evolve with redshift over the redshift range of 0.3--0.6. We also find that potential systematic effects associated with the redMaPPer algorithm may impact the location of r_sp. We discuss progress needed to understand the systematic uncertainties and fully exploit forthcoming data from DES and future surveys, emphasizing the importance of more realistic mock catalogs and independent cluster samples.

We explore the multiphase structure of the circumgalactic medium (CGM) probed by synthetic spectra through a cosmological zoom-in galaxy formation simulation. We employ a Bayesian method for modelling a combination of absorption lines to derive physical properties of absorbers with a formal treatment of detections, including saturated systems, and non-detections in a uniform manner. We find that in the lines of sight passing through localized density structures, absorption lines of low, intermediate and high ions are present in the spectrum and overlap in velocity space. Low, intermediate and high ions can be combined to derive the mass-weighted properties of a density-varying peak, although the ions are not co-spatial within the structure. By contrast, lines of sight that go through the hot halo only exhibit detectable HI and high ions. In such lines of sight, the absorption lines are typically broad due to the complex velocity fields across the entire halo. We show that the derived gas density, temperature, and metallicity match closely the corresponding HI mass-weighted averages along the LOS. We also show that when the data quality allows, our Bayesian technique allows one to recover the underlying physical properties of LOS by incorporating both detections and non-detections. It is especially useful to include non-detections, of species such as NV or NeVIII, when the number of detections of strong absorbers, such as HI and OVI, is smaller than the number of model parameters (density, temperature, and metallicity).

Using a suite of isolated $L_\star$ galaxy simulations, we show that global depletion times and star-forming gas mass fractions in simulated galaxies exhibit systematic and well-defined trends as a function of the local star formation efficiency per freefall time, $\epsilon_{\rm ff}$, strength of stellar feedback, and star formation threshold. We demonstrate that these trends can be reproduced and explained by a simple physical model of global star formation in galaxies. Our model is based on mass conservation and the idea of gas cycling between star-forming and non-star-forming states on certain characteristic time scales under the influence of dynamical and feedback processes. Both the simulation results and our model predictions exhibit two limiting regimes with rather different dependencies of global galactic properties on the local parameters. When $\epsilon_{\rm ff}$ is small and feedback is inefficient, the total star-forming mass fraction, $f_{\rm sf}$, is independent of $\epsilon_{\rm ff}$ and the global depletion time, $\tau_{\rm dep}$, scales inversely with $\epsilon_{\rm ff}$. When $\epsilon_{\rm ff}$ is large or feedback is very efficient, these trends are reversed: $f_{\rm sf} \propto \epsilon_{\rm ff}^{-1}$ and $\tau_{\rm dep}$ is independent of $\epsilon_{\rm ff}$ but scales linearly with the feedback strength. We also compare our results with the observed depletion times and mass fractions of star-forming and molecular gas and show that they provide complementary constraints on $\epsilon_{\rm ff}$ and the feedback strength. We show that useful constraints on $\epsilon_{\rm ff}$ can also be obtained using measurements of the depletion time and its scatter on different spatial scales.

We find, using cosmological simulations of galaxy clusters, that the hot X-ray emitting intra-cluster medium (ICM) enclosed within the outer accretion shock extends out to $R_{\rm shock}\sim(2 - 3) R_{\rm vir}$, where $R_{\rm vir}$ is the standard virial radius of the halo. Using a simple analytic model for satellite galaxies in the cluster, we evaluate the effect of ram pressure stripping on the gas in the inner discs and in the haloes at different distances from the cluster centre. We find that significant removal of star-forming disc gas occurs only at $r \lesssim 0.5 R_{\rm vir}$, while gas removal from the satellite halo is more effective and can occur when the satellite is found between $ R_{\rm vir}$ and $R_{\rm shock}$. Removal of halo gas sets the stage for quenching of the star formation by starvation over $2\textrm{--}3\,\mathrm{Gyr}$, prior to the satellite entry to the inner cluster halo. This scenario explains the presence of quenched galaxies, preferentially discs, at the outskirts of galaxy clusters, and the delayed quenching of satellites compared to central galaxies.

Cold Fronts and shocks are hallmarks of the complex intra-cluster medium (ICM) in galaxy clusters. They are thought to occur due to gas motions within the ICM and are often attributed to galaxy mergers within the cluster. Using hydro-cosmological simulations of clusters of galaxies, we show that collisions of inflowing gas streams, seen to penetrate to the very centre of about half the clusters, offer an additional mechanism for the formation of shocks and cold fronts in cluster cores. Unlike episodic merger events, a gas stream inflow persists over a period of several Gyrs and it could generate a particular pattern of multiple cold fronts and shocks.

An optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the violations can be considered to be "weak" in a certain sense and the corresponding numerical solution may be stable, such calculation does not guarantee the correct propagation speed for arbitrary waves. We use an experimental fluid dynamics code that allows cubic "patches" of grid cells to step with independent, locally constrained time steps to demonstrate how the CFL condition can be enforced by imposing a condition on the time steps of neighboring patches. We perform several numerical tests that illustrate errors introduced in the numerical solutions by weak CFL condition violations and show how strict enforcement of the CFL condition eliminates these errors. In all our tests the strict enforcement of the CFL condition does not impose a significant performance penalty.

We present the umbrella sampling (US) technique and show that it can be used to sample extremely low probability areas of the posterior distribution that may be required in statistical analyses of data. In this approach sampling of the target likelihood is split into sampling of multiple biased likelihoods confined within individual umbrella windows. We show that the US algorithm is efficient and highly parallel and that it can be easily used with other existing MCMC samplers. The method allows the user to capitalize on their intuition and define umbrella windows and increase sampling accuracy along specific directions in the parameter space. Alternatively, one can define umbrella windows using an approach similar to parallel tempering. We provide a public code that implements umbrella sampling as a standalone python package. We present a number of tests illustrating the power of the US method in sampling low probability areas of the posterior and show that this ability allows a considerably more robust sampling of multi-modal distributions compared to the standard sampling methods. We also present an application of the method in a real world example of deriving cosmological constraints using the supernova type Ia data. We show that umbrella sampling can sample the posterior accurately down to the $\approx 15\sigma$ credible region in the $\Omega_{\rm m}-\Omega_\Lambda$ plane, while for the same computational work the affine-invariant MCMC sampling implemented in the {\tt emcee} code samples the posterior reliably only to $\approx 3\sigma$.

We empirically constrain how galaxy size relates to halo virial radius using new measurements of the size- and stellar mass-dependent clustering of galaxies in the Sloan Digital Sky Survey. We find that small galaxies cluster much more strongly than large galaxies of the same stellar mass. The magnitude of this clustering difference increases on small scales, and decreases with increasing stellar mass. Using Halotools to forward model the observations, we test an empirical model in which present-day galaxy size is proportional to the size of the virial radius at the time the halo reached its maximum mass. This simple model reproduces the observed size-dependence of galaxy clustering in striking detail. The success of this model provides strong support for the conclusion that satellite galaxies have smaller sizes relative to central galaxies of the same halo mass. Our findings indicate that satellite size is set prior to the time of infall, and that a remarkably simple, linear size--virial radius relation emerges from the complex physics regulating galaxy size. We make quantitative predictions for future measurements of galaxy-galaxy lensing, including dependence upon size, scale, and stellar mass, and provide a scaling relation of the ratio of mean sizes of satellites and central galaxies as a function of their halo mass that can be used to calibrate hydrodynamical simulations and semi-analytic models.

We introduce a Bayesian approach for modeling Voigt profiles in absorption spectroscopy and its implementation in the python package, BayesVP, publicly available at https://github.com/cameronliang/BayesVP. The code fits the absorption line profiles within specified wavelength ranges and generates posterior distributions for the column density, Doppler parameter, and redshifts of the corresponding absorbers. The code uses publicly available efficient parallel sampling packages to sample posterior and thus can be run on parallel platforms. BayesVP supports simultaneous fitting for multiple absorption components in high-dimensional parameter space. We provide other useful utilities in the package, such as explicit specification of priors of model parameters, continuum model, Bayesian model comparison criteria, and posterior sampling convergence check.

We present a model that explains why galaxies form stars on a time scale significantly longer than the time scales of processes governing the evolution of interstellar gas. We show that gas evolves from a non-star-forming to a star-forming state on a relatively short time scale and thus the rate of this evolution does not limit the star formation rate. Instead, the star formation rate is limited because only a small fraction of star-forming gas is converted into stars before star-forming regions are dispersed by feedback and dynamical processes. Thus, gas cycles into and out of star-forming state multiple times, which results in a long time scale on which galaxies convert gas into stars. Our model does not rely on the assumption of equilibrium and can be used to interpret trends of depletion times with the properties of observed galaxies and the parameters of star formation and feedback recipes in simulations. In particular, the model explains how feedback self-regulates the star formation rate in simulations and makes it insensitive to the local star formation efficiency. We illustrate our model using the results of an isolated $L_*$-sized galaxy simulation that reproduces the observed Kennicutt-Schmidt relation for both molecular and atomic gas. Interestingly, the relation for molecular gas is almost linear on kiloparsec scales, although a nonlinear relation is adopted in simulation cells. We discuss how a linear relation emerges from non-self-similar scaling of the gas density PDF with the average gas surface density.

The splashback radius $R_{\rm sp}$, the apocentric radius of particles on their first orbit after falling into a dark matter halo, has recently been suggested as a physically motivated halo boundary that separates accreting from orbiting material. Using the SPARTA code presented in Paper I, we analyze the orbits of billions of particles in cosmological simulations of structure formation and measure $R_{\rm sp}$ for a large sample of halos that span a mass range from dwarf galaxy to massive cluster halos, reach redshift 8, and include WMAP, Planck, and self-similar cosmologies. We analyze the dependence of $R_{\rm sp}/R_{\rm 200m}$ and $M_{\rm sp}/M_{\rm 200m}$ on the mass accretion rate $\Gamma$, halo mass, redshift, and cosmology. The scatter in these relations varies between 0.02 and 0.1 dex. While we confirm the known trend that $R_{\rm sp}/R_{\rm 200m}$ decreases with $\Gamma$, the relationships turn out to be more complex than previously thought, demonstrating that $R_{\rm sp}$ is an independent definition of the halo boundary that cannot trivially be reconstructed from spherical overdensity definitions. We present fitting functions for $R_{\rm sp}/R_{\rm 200m}$ and $M_{\rm sp}/M_{\rm 200m}$ as a function of accretion rate, peak height, and redshift, achieving an accuracy of 5% or better everywhere in the parameter space explored. We discuss the physical meaning of the distribution of particle apocenters and show that the previously proposed definition of $R_{\rm sp}$ as the radius of the steepest logarithmic density slope encloses roughly three-quarters of the apocenters. Finally, we conclude that no analytical model presented thus far can fully explain our results.

The density field in the outskirts of dark matter halos is discontinuous due to a caustic formed by matter at its first apocenter after infall. In this paper, we present an algorithm to identify the "splashback shell" formed by these apocenters in individual simulated halos using only a single snapshot of the density field. We implement this algorithm in the code SHELLFISH (SHELL Finding In Spheroidal Halos) and demonstrate that the code identifies splashback shells correctly and measures their properties with an accuracy of $<5\%$ for halos with more than 50,000 particles and mass accretion rates of $\Gamma_\textrm{DK14}>0.5$. Using SHELLFISH, we present the first estimates for several basic properties of individual splashback shells, such as radius, $R_\textrm{sp}$, mass, and overdensity, and provide fits to the distribution of these quantities as functions of $\Gamma_\textrm{DK14}$, $\nu_\textrm{200m}$, and $z.$ We confirm previous findings that $R_\textrm{sp}$ decreases with increasing $\Gamma_\textrm{DK14}$, but we show that independent of accretion rate, it also decreases with increasing $\nu_\textrm{200m}$. We also study the 3D structures of these shells and find that they generally have non-ellipsoidal oval shapes. We find that splashback radii estimated by SHELLFISH are $20\%-30\%$ larger than those estimated in previous studies from stacked density profiles at high accretion rates. We demonstrate that the latter are biased low due to the contribution of high-mass subhalos to these profiles and show that using the median instead of mean density in each radial bin mitigates the effect of substructure on density profiles and removes the bias.

We present the first constraints on cosmology from the Dark Energy Survey (DES), using weak lensing measurements from the preliminary Science Verification (SV) data. We use 139 square degrees of SV data, which is less than 3\% of the full DES survey area. Using cosmic shear 2-point measurements over three redshift bins we find $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.81 \pm 0.06$ (68\% confidence), after marginalising over 7 systematics parameters and 3 other cosmological parameters. We examine the robustness of our results to the choice of data vector and systematics assumed, and find them to be stable. About $20$\% of our error bar comes from marginalising over shear and photometric redshift calibration uncertainties. The current state-of-the-art cosmic shear measurements from CFHTLenS are mildly discrepant with the cosmological constraints from Planck CMB data; our results are consistent with both datasets. Our uncertainties are $\sim$30\% larger than those from CFHTLenS when we carry out a comparable analysis of the two datasets, which we attribute largely to the lower number density of our shear catalogue. We investigate constraints on dark energy and find that, with this small fraction of the full survey, the DES SV constraints make negligible impact on the Planck constraints. The moderate disagreement between the CFHTLenS and Planck values of $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ is present regardless of the value of $w$.

This paper describes digital radio astronomical receivers used for decameter and meter wavelength observations. This paper describes digital radio astronomical receivers used for decameter and meter wavelength observations. Since 1998, digital receivers performing on-the-fly dynamic spectrum calculations or waveform data recording without data loss have been used at the UTR-2 radio telescope, the URAN VLBI system, and the GURT new generation radio telescope. Here we detail these receivers developed for operation in the strong interference environment that prevails in the decameter wavelength range. Data collected with these receivers allowed us to discover numerous radio astronomical objects and phenomena at low frequencies, a summary of which is also presented.

Mass around dark matter halos can be divided into "infalling" material and "collapsed" material that has passed through at least one pericenter. Analytical models and simulations predict a rapid drop in the halo density profile associated with the transition between these two regimes. Using data from SDSS, we explore the evidence for such a feature in the density profiles of galaxy clusters and investigate the connection between this feature and a possible phase space boundary. We first estimate the steepening of the outer galaxy density profile around clusters: the profiles show an abrupt steepening, providing evidence for truncation of the halo profile. Next, we measure the galaxy density profile around clusters using two sets of galaxies selected based on color. We find evidence of an abrupt change in the galaxy colors that coincides with the location of the steepening of the density profile. Since galaxies are likely to be quenched of star formation and turn red inside of clusters, this change in the galaxy color distribution can be interpreted as the transition from an infalling regime to a collapsed regime. We also measure this transition using a model comparison approach which has been used recently in studies of the "splashback" phenomenon, but find that this approach is not a robust way to quantify the significance of detecting a splashback-like feature. Finally, we perform measurements using an independent cluster catalog to test for potential systematic errors associated with cluster selection. We identify several avenues for future work: improved understanding of the small-scale galaxy profile, lensing measurements, identification of proxies for the halo accretion rate, and other tests. With upcoming data from the DES, KiDS and HSC surveys, we can expect significant improvements in the study of halo boundaries.

We present a new implementation of star formation in cosmological simulations, by considering star clusters as a unit of star formation. Cluster particles grow in mass over several million years at the rate determined by local gas properties, with high time resolution. The particle growth is terminated by its own energy and momentum feedback on the interstellar medium. We test this implementation for Milky Way-sized galaxies at high redshift, by comparing the properties of model clusters with observations of young star clusters. We find that the cluster initial mass function is best described by a Schechter function rather than a single power law. In agreement with observations, at low masses the logarithmic slope is $\alpha\approx 1.8-2$, while the cutoff at high mass scales with the star formation rate. A related trend is a positive correlation between the surface density of star formation rate and fraction of stars contained in massive clusters. Both trends indicate that the formation of massive star clusters is preferred during bursts of star formation. These bursts are often associated with major merger events. We also find that the median timescale for cluster formation ranges from 0.5 to 4 Myr and decreases systematically with increasing star formation efficiency. Local variations in the gas density and cluster accretion rate naturally lead to the scatter of the overall formation efficiency by an order of magnitude, even when the instantaneous efficiency is kept constant. Comparison of the formation timescale with the observed age spread of young star clusters provides an additional important constraint on the modeling of star formation and feedback schemes.

We utilize cosmological simulations of 16 galaxy clusters at redshifts \zeq{0} and \zeq{0.6} to study the effect of inflowing streams on the properties of the X-ray emitting intracluster medium. We find that the mass accretion occurs predominantly along streams that originate from the cosmic web and consist of heated gas. Clusters that are unrelaxed in terms of their X-ray morphology are characterized by higher mass inflow rates and deeper penetration of the streams, typically into the inner third of the virial radius. The penetrating streams generate elevated random motions, bulk flows and cold fronts. The degree of penetration of the streams may change over time such that clusters can switch from being unrelaxed to relaxed over a time-scale of several giga years.

We present measurements of weak gravitational lensing cosmic shear two-point statistics using Dark Energy Survey Science Verification data. We demonstrate that our results are robust to the choice of shear measurement pipeline, either ngmix or im3shape, and robust to the choice of two-point statistic, including both real and Fourier-space statistics. Our results pass a suite of null tests including tests for B-mode contamination and direct tests for any dependence of the two-point functions on a set of 16 observing conditions and galaxy properties, such as seeing, airmass, galaxy color, galaxy magnitude, etc. We furthermore use a large suite of simulations to compute the covariance matrix of the cosmic shear measurements and assign statistical significance to our null tests. We find that our covariance matrix is consistent with the halo model prediction, indicating that it has the appropriate level of halo sample variance. We compare the same jackknife procedure applied to the data and the simulations in order to search for additional sources of noise not captured by the simulations. We find no statistically significant extra sources of noise in the data. The overall detection significance with tomography for our highest source density catalog is 9.7sigma. Cosmological constraints from the measurements in this work are presented in a companion paper (DES et al. 2015).

We present a study of a star formation prescription in which star formation efficiency depends on local gas density and turbulent velocity dispersion, as suggested by direct simulations of SF in turbulent giant molecular clouds (GMCs). We test the model using a simulation of an isolated Milky Way-sized galaxy with a self-consistent treatment of turbulence on unresolved scales. We show that this prescription predicts a wide variation of local star formation efficiency per free-fall time, $\epsilon_{\rm ff} \sim 0.1 - 10\%$, and gas depletion time, $t_{\rm dep} \sim 0.1 - 10$ Gyr. In addition, it predicts an effective density threshold for star formation due to suppression of $\epsilon_{\rm ff}$ in warm diffuse gas stabilized by thermal pressure. We show that the model predicts star formation rates in agreement with observations from the scales of individual star-forming regions to the kiloparsec scales. This agreement is non-trivial, as the model was not tuned in any way and the predicted star formation rates on all scales are determined by the distribution of the GMC-scale densities and turbulent velocities $\sigma$ in the cold gas within the galaxy, which is shaped by galactic dynamics. The broad agreement of the star formation prescription calibrated in the GMC-scale simulations with observations, both gives credence to such simulations and promises to put star formation modeling in galaxy formation simulations on a much firmer theoretical footing.

We analyze circumgalactic medium (CGM) in a suite of high-resolution cosmological re-simulations of a Milky-Way size galaxy and show that CGM properties are quite sensitive to details of star formation--feedback loop modelling. The simulation that produces a realistic late-type galaxy, fails to reproduce existing observations of the CGM. In contrast, simulation that does not produce a realistic galaxy has the predicted CGM in better agreement with observations. This illustrates that properties of galaxies and properties of their CGM provide strong ${\it complementary}$ constraints on the processes governing galaxy formation. Our simulations predict that column density profiles of ions are well described by an exponential function of projected distance $d$: $N \propto e^{-d/h_s}$. Simulations thus indicate that the sharp drop in absorber detections at larger distances in observations does not correspond to a "boundary" of an ion, but reflects the underlying steep exponential column density profile. Furthermore, we find that ionization energy of ions is tightly correlated with the scale height $h_s$: $h_s \propto E_{\rm ion}^{0.74}$. At $z \approx 0$, warm gas traced by low-ionization species (e.g., Mg II and C IV) has $ h_s \approx 0.03-0.07 R_{\rm vir}$, while higher ionization species (O VI and Ne VIII) have $h_s \approx 0.32-0.45R_{\rm vir}$. Finally, the scale heights of ions in our simulations evolve slower than the virial radius for $z\leq 2$, but similarly to the halo scale radius, $r_s$. Thus, we suggest that the column density profiles of galaxies at different redshifts should be scaled by $r_s$ rather than the halo virial radius.

The fact that the clustering of dark matter halos depends not only on their mass, but also the formation epoch, is a prominent, albeit subtle, feature of the cold dark matter structure formation theory, and is known as assembly bias. At low mass scales ($\sim 10^{12}\,h^{-1}M_\odot$), early-forming halos are predicted to be more strongly clustered than the late-forming ones. In this study we aim to robustly detect the signature of assembly bias observationally, making use of formation time indicators of central galaxies in low mass halos as a proxy for the halo formation history. Weak gravitational lensing is employed to ensure our early- and late-forming halo samples have similar masses, and are free of contamination of satellites from more massive halos. For the two formation time indicators used (resolved star formation history and current specific star formation rate), we do not find convincing evidence of assembly bias. For a pair of early- and late-forming galaxy samples with mean mass $M_{200c} \approx 9\times10^{11}\,h^{-1}M_\odot$, the relative bias is $1.00\pm 0.12$. We attribute the lack of detection to the possibilities that either the current measurements of these indicators are too noisy, or they do not correlate well with the halo formation history. Alternative proxies for the halo formation history that should perform better are suggested for future studies.

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$\sigma$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $\Lambda$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies.

We use cosmological zoom-in simulations of galaxy formation in a Milky Way (MW)-sized halo started from identical initial conditions to investigate the evolution of galaxy sizes, baryon fractions, morphologies and angular momenta in runs with different parameters of the star formation--feedback cycle. Our fiducial model with a high local star formation efficiency, which results in efficient feedback, produces a realistic late-type galaxy that matches the evolution of basic properties of late-type galaxies: stellar mass, disk size, morphology dominated by a kinematically cold disk, stellar and gas surface density profiles, and specific angular momentum. We argue that feedback's role in this success is twofold: (1) removal of low-angular momentum gas and (2) maintaining a low disk-to-halo mass fraction which suppresses disk instabilities that lead to angular momentum redistribution and a central concentration of baryons. However, our model with a low local star formation efficiency, but large energy input per supernova, chosen to produce a galaxy with a similar star formation history as our fiducial model, leads to a highly irregular galaxy with no kinematically cold component, overly extended stellar distribution and low angular momentum. This indicates that only when feedback is allowed to become vigorous via locally efficient star formation in dense cold gas, resulting galaxy sizes, gas/stellar surface density profiles and stellar disk angular momenta agree with observed $z=0$ galaxies.

The boundaries of cold dark matter halos are commonly defined to enclose a density contrast $\Delta$ relative to a reference (mean or critical) density. We argue that a more physical boundary of halos is the radius at which accreted matter reaches its first orbital apocenter after turnaround. This splashback radius, $R_{sp}$, manifests itself as a sharp density drop in the halo outskirts, at a location that depends upon the mass accretion rate. We present calibrations of $R_{sp}$ and the enclosed mass, $M_{sp}$, as a function of the accretion rate and alternatively peak height. We find that $R_{sp}$ varies between $\approx0.8-1R_{200m}$ for rapidly accreting halos and $\approx1.5R_{200m}$ for slowly accreting halos. The extent of a halo and its associated environmental effects can thus extend well beyond the conventionally defined "virial" radius. We show that $M_{sp}$ and $R_{sp}$ evolve relatively strongly compared to other commonly used definitions. In particular, $M_{sp}$ evolves significantly even for the smallest dwarf-sized halos at $z=0$. We also contrast $M_{sp}$ with the mass enclosed within four scale radii of the halo density profile, $M_{<4rs}$, which characterizes the inner halo. During the early stages of halo assembly, $M_{sp}$ and $M_{<4rs}$ evolve similarly, but in the late stages $M_{<4rs}$ stops increasing while $M_{sp}$ continues to grow significantly. This illustrates that halos at low $z$ can have "quiet" interiors while continuing to accrete mass in their outskirts. We discuss potential observational estimates of the splashback radius and show that it may already have been detected in galaxy clusters.

We investigate the star formation-feedback cycle in cosmological galaxy formation simulations, focusing on progenitors of Milky Way (MW)-sized galaxies. We find that in order to reproduce key properties of the MW progenitors, such as semi-empirically derived star formation histories and the shape of rotation curves, our implementation of star formation and stellar feedback requires 1) a combination of local early momentum feedback via radiation pressure and stellar winds and subsequent efficient supernovae feedback, and 2) efficacy of feedback that results in self-regulation of the global star formation rate on kiloparsec scales. We show that such feedback-driven self-regulation is achieved globally for a local star formation efficiency per free fall time of $\epsilon_{\rm ff}\approx 10\%$. Although this value is larger that the $\epsilon_{\rm ff}\sim 1\%$ value usually inferred from the Kennicutt-Schmidt (KS) relation, we show that it is consistent with direct observational estimates of $\epsilon_{\rm ff}$ in molecular clouds. Moreover, we show that simulations with local efficiency of $\epsilon_{\rm ff}\approx 10\%$ reproduce the global observed KS relation. Such simulations also reproduce the cosmic star formation history of the Milky Way sized galaxies and satisfy a number of other observational constraints. Conversely, we find that simulations that a priori assume an inefficient mode of star formation, instead of achieving it via stellar feedback regulation, fail to produce sufficiently vigorous outflows and do not reproduce observations. This illustrates the importance of understanding the complex interplay between star formation and feedback and the detailed processes that contribute to the feedback-regulated formation of galaxies.

We revisit previous suggestions that the heating provided by the winds of dying low-mass stars plays an important role in preventing star formation in quiescent galaxies. At the end of their asymptotic giant branch phase, intermediate and low-mass stars eject their envelopes rapidly in a super-wind phase, usually giving rise to planetary nebulae. In spheroidal galaxies with high stellar velocity dispersions, the interaction of these ejected envelopes with the ambient diffuse gas can lead to significant, isotropic and steady-state heating that scales as $\dot{M}_\ast\sigma_\ast^2$. We show that cooling of the central regions of the hot diffuse halo gas can be delayed for a Hubble time for halos more massive than $\sim10^{12.5}M_{\odot}$ at $0<z<2$, although stellar heating alone is unlikely to forestall cooling in the most massive clusters at $z=0$. This mechanism provides a natural explanation for the strong trend of galaxy quiescence with stellar surface density and velocity dispersion. In addition, since the ejected material will thermalize to $kT\sim\sigma_\ast^2$, this mechanism provides an explanation for the observed similarity between the central temperature of the hot diffuse gas and $\sigma_\ast^2$, a result which is not trivial in light of the short inferred cooling times of the hot gas. The main uncertainty in this analysis is the ultimate fate of the stellar ejecta. Preventing accumulation of the ejecta in the central regions may require energy input from another source, such as Type Ia supernovae. Detailed simulations of the interaction of the stellar wind with the ambient gas are required to better quantify the net effect of AGB heating.

We present a numerical study of dark matter halo concentrations in $\Lambda$CDM and self-similar cosmologies. We show that the relation between concentration, $c$, and peak height, $\nu$, exhibits the smallest deviations from universality if halo masses are defined with respect to the critical density of the universe. These deviations can be explained by the residual dependence of concentration on the local slope of the matter power spectrum, $n$, which affects both the normalization and shape of the $c$-$\nu$ relation. In particular, there is no well-defined floor in the concentration values. Instead, the minimum concentration depends on redshift: at fixed $\nu$, halos at higher $z$ experience steeper slopes $n$, and thus have lower minimum concentrations. We show that the concentrations in our simulations can be accurately described by a universal seven-parameter function of only $\nu$ and $n$. This model matches our $\Lambda$CDM results to $\lesssim 5\%$ accuracy up to $z = 6$, and matches scale-free $\Omega_{\rm m} = 1$ models to $\lesssim 15\%$. The model also reproduces the low concentration values of Earth--mass halos at $z \approx 30$, and thus correctly extrapolates over $16$ orders of magnitude in halo mass. The predictions of our model differ significantly from all models previously proposed in the literature at high masses and redshifts. Our model is in excellent agreement with recent lensing measurements of cluster concentrations.

The stellar mass-halo mass relation is a key constraint in all semi-analytic, numerical, and semi-empirical models of galaxy formation and evolution. However, its exact shape and redshift dependence remain debated. Several recent works support a relation in the local Universe steeper than previously thought. Based on the comparisons with a variety of data on massive central galaxies, we show that this steepening holds up to z~1, for stellar masses Mstar>2e11 Msun. Specifically, we find significant evidence for a high-mass end slope of \beta>0.35-0.70, instead of the usual \beta~0.20-0.30 reported by a number of previous results. When including the independent constraints from the recent BOSS clustering measurements, the data, independent of any systematic errors in stellar masses, tend to favor a model with a very small scatter (< 0.15 dex) in stellar mass at fixed halo mass, in the redshift range z < 0.8 and for Mstar>3e11 Msun, suggesting a close connection between massive galaxies and host halos even at relatively recent epochs. We discuss the implications of our results with respect to the evolution of the most massive galaxies since z~1.

We use recent proper motion measurements of the tangential velocity of M31, along with its radial velocity and distance, to derive the likelihood of the sum of halo masses of the Milky Way and M31. This is done using a sample halo pairs in the Bolshoi cosmological simulation of $\Lambda$CDM cosmology selected to match properties and environment of the Local Group. The resulting likelihood gives estimate of the sum of masses of $M_{\rm MW,200}+M_{\rm M31,200}=$ $2.40_{-1.05}^{+1.95}\times10^{12}\,M_{\odot}$ ($90\%$ confidence interval). This estimate is consistent with individual mass estimates for the Milky Way and M31 and is consistent, albeit somewhat on the low side, with the mass estimated using the timing argument. We show that although the timing argument is unbiased on average for all pairs, for pairs constrained to have radial and tangential velocities similar to that of the Local Group the argument overestimates the sum of masses by a factor of $1.6$. Using similar technique we estimate the total dark matter mass enclosed within $1$ Mpc from the Local Group barycenter to be $M_{\rm LG}(r<1\, {\rm Mpc})=4.2_{-2.0}^{+3.4}\times10^{12}\,M_{\odot}$ ($90\%$ confidence interval).

We present a systematic study of the density profiles of LCDM halos, focusing on the outer regions, 0.1 < r/Rvir < 9. We show that the median and mean profiles of halo samples of a given peak height exhibit significant deviations from the universal analytic profiles discussed previously in the literature, such as the Navarro-Frenk-White and Einasto profiles, at radii r > 0.5 R200m. In particular, at these radii the logarithmic slope of the median density profiles of massive or rapidly accreting halos steepens more sharply than predicted. The steepest slope of the profiles occurs at r ~ R200m, and its absolute value increases with increasing peak height or mass accretion rate, reaching slopes of -4 and steeper. Importantly, we find that the outermost density profiles at r > R200m are remarkably self-similar when radii are rescaled by R200m. This self-similarity indicates that radii defined with respect to the mean density are preferred for describing the structure and evolution of the outer profiles. However, the inner density profiles are most self-similar when radii are rescaled by R200c. We propose a new fitting formula that describes the median and mean profiles of halo samples selected by their peak height or mass accretion rate with accuracy < 10% at all radii, redshifts and masses we studied, r < 9 Rvir, 0 < z < 6 and Mvir > 1.7E10 Msun/h. We discuss observational signatures of the profile features described above, and show that the steepening of the outer profile should be detectable in future weak-lensing analyses of massive clusters. Such observations could be used to estimate the mass accretion rate of cluster halos.

We study relation between stellar mass and halo mass for high-mass halos using a sample of galaxy clusters with accurate measurements of stellar masses from optical and IR data and total masses from X-ray observations. We find that stellar mass of the brightest cluster galaxies (BCGs) scales as M*BCG\propto M500^a_BCG with the best fit slope of a_BCG~0.35+-0.1 and scatter of M*BCG at a fixed M500 of ~0.2 dex. We show that M*-M relations from abundance matching or halo modelling reported in recent studies underestimate stellar masses of BCGs by a factor of ~2-4, because these studies used stellar mass functions (SMF) based on photometry that severely underestimates the outer surface brightness profiles of massive galaxies. We show that M*-M relation derived using abundance matching with the recent SMF calibration by Bernardi et al. (2013) based on improved photometry is in a much better agreement with the relation we derive. The total stellar mass of galaxies correlates with total mass M500 with the slope of \approx 0.6+-0.1 and scatter of 0.1 dex. This indicates that efficiency with which baryons are converted into stars decreases with increasing cluster mass. We show that for a fixed choice of the initial mass function (IMF) total stellar fraction in clusters is only a factor of ~3-5 lower than the peak stellar fraction reached in M\approx 10^12 Msun halos, and only a factor of ~1.5-3 if the IMF becomes progressively more bottom heavy with increasing mass in early type galaxies, as indicated by recent observations. The larger normalization and slope of the M*-M relation derived in this study shows that the overall efficiency of star formation in massive halos is suppressed much less than was thought before and that feedback and associated suppression of star formation in massive halos should be weaker than assumed in most of the current semi-analytic models and simulations.

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.

Understanding the evolution of scaling relations between the observable properties of clusters and their total mass is key to realizing their potential as cosmological probes. In this study, we investigate whether the evolution of cluster scaling relations is affected by the spurious evolution of mass caused by the evolving reference density with respect to which halo masses are defined (pseudo-evolution). We use the relation between mass, M, and velocity dispersion, sigma, as a test case, and show that the deviation from the M-sigma relation of cluster-sized halos caused by pseudo-evolution is smaller than 10% for a wide range of mass definitions. The reason for this small impact is a tight relation between the velocity dispersion and mass profiles, sigma(<r) = const * (GM(<r) / r)^1/2, which holds across a wide range of radii. We show that such a relation is generically expected for a variety of density profiles, as long as halos are in approximate Jeans equilibrium. Thus, as the outer "virial" radius used to define the halo mass, R, increases due to pseudo-evolution, halos approximately preserve their M-sigma relation. This result highlights the fact that tight scaling relations are the result of tight equilibrium relations between radial profiles of physical quantities. We find exceptions at very small and very large radii, where the profiles deviate from the relations they exhibit at intermediate radii. We discuss the implications of these results for other cluster scaling relations, and argue that pseudo-evolution should have a small effect on most scaling relations, except for those that involve the stellar masses of galaxies. In particular, we show that the relation between stellar-mass fraction and total mass is affected by pseudo-evolution and is largely shaped by it for halo masses of less than 10^14 solar masses.

We assemble a statistical set of global mass models for ~2,000 nearly spherical SDSS galaxies at a mean redshift of 0.12 based on their aperture velocity dispersions and newly derived luminosity profiles in conjunction with published velocity dispersion profiles and empirical properties and relations of galaxy and halo parameters. When two-component (i.e. stellar plus dark) mass models are fitted to the SDSS aperture velocity dispersions, the predicted velocity dispersion profile (VP) slopes within the effective radius R_eff match well the distribution in observed elliptical galaxies. In contrast, the single-component models cannot reproduce the VP slope distribution. From a number of input variations the models exhibit for the radial range 0.1 R_eff < r < R_eff a tight correlation <gamma_e>=(1.865+/-0.008)+(-4.93+/-0.15)<eta> where <gamma_e> is the mean slope absolute value of the total mass density and <eta> is the mean slope of the velocity dispersion profile, which leads to a super-isothermal <gamma_e> = 2.15+/-0.04 for <eta>=-0.058+/-0.008 in observed elliptical galaxies. Furthermore, the successful two-component models appear to imply a typical slope curvature pattern in the total mass profile because for the observed steep luminosity (stellar mass) profile and the weak lensing inferred halo profile at large radii a total mass profile with monotonically varying slope would require too high DM density in the optical region giving rise to too large aperture velocity dispersion and too shallow VP.

We present results from high-resolution hydrodynamic simulations of isolated SMC- and Milky Way-sized galaxies that include a model for feedback from galactic cosmic rays (CRs). We find that CRs are naturally able to drive winds with mass loading factors of up to ~10 in dwarf systems. The scaling of the mass loading factor with circular velocity between the two simulated systems is consistent with \propto v_c^{1-2} required to reproduce the faint end of the galaxy luminosity function. In addition, simulations with CR feedback reproduce both the normalization and the slope of the observed trend of wind velocity with galaxy circular velocity. We find that winds in simulations with CR feedback exhibit qualitatively different properties compared to SN driven winds, where most of the acceleration happens violently in situ near star forming sites. In contrast, the CR-driven winds are accelerated gently by the large-scale pressure gradient established by CRs diffusing from the star-forming galaxy disk out into the halo. The CR-driven winds also exhibit much cooler temperatures and, in the SMC-sized system, warm (T~10^4 K) gas dominates the outflow. The prevalence of warm gas in such outflows may provide a clue as to the origin of ubiquitous warm gas in the gaseous halos of galaxies detected via absorption lines in quasar spectra.

A dark matter halo is commonly defined as a spherical overdensity of matter with respect to a reference density, such as the critical density or the mean matter density of the Universe. Such definitions can lead to a spurious pseudo-evolution of halo mass simply due to redshift evolution of the reference density, even if its physical density profile remains constant over time. We estimate the amount of such pseudo-evolution of mass between z=1 to 0 for halos identified in a large N-body simulation, and show that it accounts for almost the entire mass evolution of the majority of halos with M200 of about 1E12 solar masses and can be a significant fraction of the apparent mass growth even for cluster-sized halos. We estimate the magnitude of the pseudo-evolution assuming that halo density profiles remain static in physical coordinates, and show that this simple model predicts the pseudo-evolution of halos identified in numerical simulations to good accuracy, albeit with significant scatter. We discuss the impact of pseudo-evolution on the evolution of the halo mass function and show that the non-evolution of the low-mass end of the halo mass function is the result of a fortuitous cancellation between pseudo-evolution and the absorption of small halos into larger hosts. We also show that the evolution of the low mass end of the concentration-mass relation observed in simulations is almost entirely due to the pseudo-evolution of mass. Finally, we discuss the implications of our results for the interpretation of the evolution of various scaling relations between the observable properties of galaxies and galaxy clusters and their halo masses.

X-ray observations of galaxy clusters provide emission measure weighted spectra, arising from a range of density and temperature fluctuations in the intra-cluster medium (ICM). This is fitted to a single temperature plasma emission model to provide an estimate of the gas density and temperature, which are sensitive to the gas inhomogeneities. Therefore, X-ray observations yield a potentially biased estimate of the thermal gas pressure, P_X. At the same time Sunyaev-Zeldovich (SZ) observations directly measure the integrated gas pressure, P_SZ. If the X-ray pressure profiles are strongly biased with respect to the SZ, then one has the possibility to probe the gas inhomogeneities, even at scales unresolved by the current generation of telescopes. At the same time, a weak bias has implications for the interchangeable use of mass proxies like Y_SZ and Y_X as cosmological probes. In this paper we investigate the dependence of the bias, defined as b_P(r)=P_X(r)/P_SZ(r)-1, on the characteristics of fluctuations in the ICM taking into account the correlation between temperature and density fluctuations. We made a simple prediction of the irreducible bias in idealised X-ray vs SZ observations using multi-temperature plasma emission model. We also provide a simple fitting form to estimate the bias given the distribution of fluctuations. Analysing a sample of 16 clusters extracted from hydrodynamical simulations, we find that the median value of bias is within +/-3% within R_500, it decreases to -5% at R_500<r<1.5R_500 and then rises back to ~0% at r >~ 2R_500. The scatter of b_P(r) between individual relaxed clusters is small -- at the level of <0.03 within R_500, but turns significantly larger (0.25) and highly skewed at r>~ 1.5R_500. For any relaxed cluster we find |b_P(r)| < 15% within R_500, across different implementations of input physics in the simulations. [abridged]

Previous studies showed that an estimate of the likelihood distribution of the Milky Way halo mass can be derived using the properties of the satellites similar to the Large and Small Magellanic Clouds (LMC and SMC). However, it would be straightforward to interpret such an estimate only if the properties of the Magellanic Clouds (MCs) are fairly typical and are not biased by the environment. In this study we explore whether the environment of the Milky Way affects the properties of the SMC and LMC such as their velocities. To test for the effect of the environment, we compare velocity distributions for MC-sized subhalos around Milky Way hosts in a sample selected simply by mass and in the second sample of such halos selected with additional restrictions on the distance to the nearest cluster and the local galaxy density, designed to mimic the environment of the Local Group (LG). We find that satellites in halos in the LG-like environments do have somewhat larger velocities, as compared to the halos of similar mass in the sample without environmental constraints. We derive the host halo likelihood distribution for the samples in the LG-like envirionment and in the control sample and find that the environment does not significantly affect the derived likelihood. We use the updated properties of the SMC and LMC to derive the constraint on the MW halo mass $\log{({\rm M}_{200} /\msol)}=12.06^{+0.31}_{-0.19}$ (90% confidence interval). We also explore the incidence of close pairs with relative velocities and separations similar to those of the LMC and SMC and find that such pairs are quite rare among $\Lambda$CDM halos. Taking into account the close separation of the MCs in the Busha et al.\ 2011 method results in the shift of the MW halo mass estimate to smaller masses, with the peak shifting approximately by a factor of two.[Abridged]

We study the high column density regime of the HI column density distribution function and argue that there are two distinct features: a turnover at NHI ~ 10^21 cm^-2 which is present at both z=0 and z ~ 3, and a lack of systems above NHI ~ 10^22 cm^-2 at z=0. Using observations of the column density distribution, we argue that the HI-H2 transition does not cause the turnover at NHI ~ 10^21 cm^-2, but can plausibly explain the turnover at NHI > 10^22 cm^-2. We compute the HI column density distribution of individual galaxies in the THINGS sample and show that the turnover column density depends only weakly on metallicity. Furthermore, we show that the column density distribution of galaxies, corrected for inclination, is insensitive to the resolution of the HI map or to averaging in radial shells. Our results indicate that the similarity of HI column density distributions at z=3 and z=0 is due to the similarity of the maximum HI surface densities of high-z and low-z disks, set presumably by universal processes that shape properties of the gaseous disks of galaxies. Using fully cosmological simulations, we explore other candidate physical mechanisms that could produce a turnover in the column density distribution. We show that while turbulence within GMCs cannot affect the DLA column density distribution, stellar feedback can affect it significantly if the feedback is sufficiently effective in removing gas from the central 2-3 kpc of high-redshift galaxies. Finally, we argue that it is meaningful to compare column densities averaged over ~ kpc scales with those estimated from quasar spectra which probe sub-pc scales due to the steep power spectrum of HI column density fluctuations observed in nearby galaxies.

Sizes of galaxies are an important diagnostic for galaxy formation models. In this study I use the abundance matching ansatz, which has proven to be successful in reproducing galaxy clustering and other statistics, to derive estimates of the virial radius, R200, for galaxies of different morphological types and wide range of stellar mass. I show that over eight of orders of magnitude in stellar mass galaxies of all morphological types follow an approximately linear relation between 3D half-mass radius of their stellar distribution, rhalf and virial radius, rhalf~0.015R200 with a scatter of ~0.2 dex. Such scaling is in remarkable agreement with expectation of models which assume that galaxy sizes are controlled by halo angular momentum, which implies rhalf\propto lambda R200, where lambda is the spin of galaxy parent halo. The scatter about the relation is comparable with the scatter expected from the distribution of $\lambda$ and normalization of the relation agrees with that predicted by the model of Mo, Mao & White (1998), if galaxy sizes were set on average at z~1-2. Moreover, I show that when stellar and gas surface density profiles of galaxies of different morphological types are rescaled using radius r_n= 0.015 R200, the rescaled surface density profiles follow approximately universal exponential (for late types) and de Vaucouleurs (for early types) profiles with scatter of only 30-50% at R~1-3r_n. Remarkably, both late and early type galaxies have similar mean stellar surface density profiles at R>r_n. The main difference between their stellar distributions is thus at R<r_n. The results of this study imply that galaxy sizes and radial distribution of baryons are shaped primarily by properties of their parent halo and that sizes of both late type disks and early type spheroids are controlled by halo angular momentum.

The next generation of wide-area sky surveys offer the power to place extremely precise constraints on cosmological parameters and to test the source of cosmic acceleration. These observational programs will employ multiple techniques based on a variety of statistical signatures of galaxies and large-scale structure. These techniques have sources of systematic error that need to be understood at the percent-level in order to fully leverage the power of next-generation catalogs. Simulations of large-scale structure provide the means to characterize these uncertainties. We are using XSEDE resources to produce multiple synthetic sky surveys of galaxies and large-scale structure in support of science analysis for the Dark Energy Survey. In order to scale up our production to the level of fifty 10^10-particle simulations, we are working to embed production control within the Apache Airavata workflow environment. We explain our methods and report how the workflow has reduced production time by 40% compared to manual management.

We present a new method to identify and characterize the structure of the intracluster medium (ICM) in simulated galaxy clusters. The method uses the median of gas properties, such as density and pressure, which we show to be very robust to the presence of gas inhomogeneities. In particular, we show that the radial profiles of median gas properties are smooth and do not exhibit fluctuations at locations of massive clumps in contrast to mean and mode properties. It is shown that distribution of gas properties in a given radial shell can be well described by a log-normal PDF and a tail. The former corresponds to a nearly hydrostatic bulk component, accounting for ~99% of the volume, while the tail corresponds to high density inhomogeneities. We show that this results in a simple and robust separation of the diffuse and clumpy components of the ICM. The FWHM of the density distribution grows with radius and varies from ~0.15 dex in cluster centre to ~0.5 dex at 2r_500 in relaxed clusters. The small scatter in the width between relaxed clusters suggests that the degree of inhomogeneity is a robust characteristic of the ICM. It broadly agrees with the amplitude of density perturbations in the Coma cluster. We discuss the origin of ICM density variations in spherical shells and show that less than 20% of the width can be attributed to the triaxiality of the cluster gravitational potential. As a link to X-ray observations of real clusters we evaluated the ICM clumping factor with and without high density inhomogeneities. We argue that these two cases represent upper and lower limits on the departure of the observed X-ray emissivity from the median value. We find that the typical value of the clumping factor in the bulk component of relaxed clusters varies from ~1.1-1.2 at r_500 up to ~1.3-1.4 at r_200, in broad agreement with recent observations.

Stellar feedback plays a key role in galaxy formation by regulating star formation, driving interstellar turbulence and generating galactic scale outflows. Although modern simulations of galaxy formation can resolve scales of 10-100 pc, star formation and feedback operate on smaller, "subgrid" scales. Great care should therefore be taken in order to properly account for the effect of feedback on global galaxy evolution. We investigate the momentum and energy budget of feedback during different stages of stellar evolution, and study its impact on the interstellar medium using simulations of local star forming regions and galactic disks at the resolution affordable in modern cosmological zoom-in simulations. In particular, we present a novel subgrid model for the momentum injection due to radiation pressure and stellar winds from massive stars during early, pre-supernova evolutionary stages of young star clusters. Early injection of momentum acts to clear out dense gas in star forming regions, hence limiting star formation. The reduced gas density mitigates radiative losses of thermal feedback energy from subsequent supernova explosions, leading to an increased overall efficiency of stellar feedback. The detailed impact of stellar feedback depends sensitively on the implementation and choice of parameters. Somewhat encouragingly, we find that implementations in which feedback is efficient lead to approximate self-regulation of global star formation efficiency. We compare simulation results using our feedback implementation to other phenomenological feedback methods, where thermal feedback energy is allowed to dissipate over time scales longer than the formal gas cooling time. We find that simulations with maximal momentum injection suppress star formation to a similar degree as is found in simulations adopting adiabatic thermal feedback.

With SDSS galaxy data and halo data from up-to-date N-body simulations we construct a semi-empirical catalog (SEC) of early-type systems by making a self-consistent bivariate statistical match of stellar mass (M_star) and velocity dispersion (sigma) with halo virial mass (M_vir). We then assign stellar mass profile and velocity dispersion profile parameters to each system in the SEC using their observed correlations with M_star and sigma. Simultaneously, we solve for dark matter density profile of each halo using the spherical Jeans equation. The resulting dark matter density profiles deviate in general from the dissipationless profile of NFW or Einasto and their mean inner density slope and concentration vary systematically with M_vir. Statistical tests of the distribution of profiles at fixed M_vir rule out the null hypothesis that it follows the distribution predicted by N-body simulations for M_vir ~< 10^{13.5-14.5} M_solar. These dark matter profiles imply that dark matter density is, on average, enhanced significantly in the inner region of halos with M_vir ~< 10^{13.5-14.5} M_solar supporting halo contraction. The main characteristics of halo contraction are: (1) the mean dark matter density within the effective radius has increased by a factor varying systematically up to ~ 3-4 at M_vir = 10^{12} M_solar, and (2) the inner density slope has a mean of <alpha> ~ 1.3 with rho(r) ~ r^{-alpha} and a halo-to-halo rms scatter of rms(alpha) ~ 0.4-0.5 for 10^{12} M_solar ~< M_vir ~< 10^{13-14} M_solar steeper than the NFW profile (alpha=1). Based on our results we predict that halos of nearby elliptical and lenticular galaxies can, in principle, be promising targets for gamma-ray emission from dark matter annihilation.

Simulations of cluster formation have demonstrated that condensation of baryons into central galaxies during cluster formation can drive the shape of the gas distribution in galaxy clusters significantly rounder, even at radii as large as half of the virial radius. However, such simulations generally predict stellar fractions within cluster virial radii that are ~2 to 3 times larger than the stellar masses deduced from observations. In this work we compare ellipticity profiles of clusters simulated with and without baryonic cooling to the cluster ellipticity profiles derived from Chandra and ROSAT observations in an effort to constrain the fraction of gas that cools and condenses into the central galaxies within clusters. We find that the observed ellipticity profiles are fairly constant with radius, with an average ellipticity of 0.18 +/- 0.05. The observed ellipticity profiles are in good agreement with the predictions of non-radiative simulations. On the other hand, the ellipticity profiles of the clusters in simulations that include radiative cooling, star formation, and supernova feedback (but no AGN feedback) deviate significantly from the observed ellipticity profiles at all radii. The simulations with cooling overpredict (underpredict) ellipticity in the inner (outer) regions of galaxy clusters. By comparing the simulations with and without cooling, we show that the cooling of gas via cooling flows in the central regions of simulated clusters causes the gas distribution to be more oblate in the central regions, but makes the outer gas distribution more spherical. We find that late-time gas cooling and star formation are responsible for the significantly oblate gas distributions in cluster cores, but the gas shapes outside of cluster cores are set primarily by baryon dissipation at high redshift z > 2.

In this review, we describe our current understanding of cluster formation: from the general picture of collapse from initial density fluctuations in an expanding Universe to detailed simulations of cluster formation including the effects of galaxy formation. We outline both the areas in which highly accurate predictions of theoretical models can be obtained and areas where predictions are uncertain due to uncertain physics of galaxy formation and feedback. The former includes the description of the structural properties of the dark matter halos hosting cluster, their mass function and clustering properties. Their study provides a foundation for cosmological applications of clusters and for testing the fundamental assumptions of the standard model of structure formation. The latter includes the description of the total gas and stellar fractions, the thermodynamical and non-thermal processes in the intracluster plasma. Their study serves as a testing ground for galaxy formation models and plasma physics. In this context, we identify a suitable radial range where the observed thermal properties of the intra-cluster plasma exhibit the most regular behavior and thus can be used to define robust observational proxies for the total cluster mass. We put particular emphasis on examining assumptions and limitations of the widely used self-similar model of clusters. Finally, we discuss the formation of clusters in non-standard cosmological models, such as non-Gaussian models for the initial density field and models with modified gravity, along with prospects for testing these alternative scenarios with large cluster surveys in the near future.

There is ample observational evidence that the star formation rate (SFR) surface density, Sigma_SFR, is closely correlated with the surface density of molecular hydrogen, Sigma_H2. This empirical relation holds both for galaxy-wide averages and for individual >=kpc sized patches of the interstellar medium (ISM), but appears to degrade substantially at a sub-kpc scale. Identifying the physical mechanisms that determine the scale-dependent properties of the observed Sigma_H2-Sigma_SFR relation remains a challenge from a theoretical perspective. To address this question, we analyze the slope and scatter of the Sigma_H2-Sigma_SFR relation using a set of cosmological, galaxy formation simulations with a peak resolution of ~100 pc. These simulations include a chemical network for molecular hydrogen, a model for the CO emission, and a simple, stochastic prescription for star formation that operates on ~100 pc scales. Specifically, star formation is modeled as a Poisson process in which the average SFR is directly proportional to the present mass of H2. The predictions of our numerical model are in good agreement with the observed Kennicutt-Schmidt and Sigma_H2-Sigma_SFR relations. We show that observations based on CO emission are ill suited to reliably measure the slope of the latter relation at low (<20 M_sun pc^-2) H2 surface densities on sub-kpc scales. Our models also predict that the inferred Sigma_H2-Sigma_SFR relation steepens at high H2 surface densities as a result of the surface density dependence of the CO/H2 conversion factor. Finally, we show that on sub-kpc scales most of the scatter in the relation is a consequence of discreteness effects in the star formation process. In contrast, variations of the CO/H2 conversion factor are responsible for most of the scatter measured on super-kpc scales.

Future X-ray observations of galaxy clusters by high spectral resolution missions will provide spatially resolved measurements of the energy and width for the brightest emission lines in the intracluster medium (ICM) spectrum. In this paper we discuss various ways of using these high resolution data to constrain velocity power spectrum in galaxy clusters. We argue that variations of these quantities with the projected distance R in cool core clusters contain important information on the velocity field length scales in the ICM. The effective length $l_{\rm eff}$ along the line of sight, which provides dominant contribution to the line flux, increases with R, allowing one to probe the amplitude of the velocity variations at different spatial scales. In particular, we show that the width of the line as a function of R is closely linked to the structure function of the 3D velocity field. Yet another easily obtainable proxy of the velocity field length scales is the ratio of the amplitude of the projected velocity field (line energy) variations to the dispersion of the velocity along the line of sight (line width). Finally the projected velocity field can be easily converted into 3D velocity field, especially for clusters like Coma with an extended flat core in the surface brightness. Under assumption of a homogeneous isotropic Gaussian 3D velocity field we derived simple expressions relating the power spectrum of the 3D velocity field (or structure function) and the observables. The uncertainties in the observables, caused by stochastic nature of the velocity field, are estimated by making multiple realizations of the random Gaussian velocity field and evaluating the scatter in observables. If large scale motions are present in the ICM these uncertainties may dominate the statistical errors of line width and shift measurements.

Characterizing the conversion factor between CO emission and column density of molecular hydrogen, X_CO, is crucial in studying the gaseous content of galaxies, its evolution, and relation to star formation. In most cases the conversion factor is assumed to be close to that of giant molecular clouds (GMCs) in the Milky Way, except possibly for mergers and star-bursting galaxies. However, there are physical grounds to expect that it should also depend on the gas metallicity, surface density, and strength of the interstellar radiation field. The XCO factor may also depend on the scale on which CO emission is averaged due to effects of limited resolution. We study the dependence of X_CO on gas properties and averaging scale using a model that is based on a combination of results of sub-pc scale magneto-hydrodynamic simulations and on the gas distribution from self-consistent cosmological simulations of galaxy formation. Our model predicts a value of X_CO that is consistent with the Galactic value for interstellar medium conditions typical for the Milky Way. For such conditions the predicted X_CO varies by only a factor of two for gas surfaced densities in the range \sim 50 - 500 M_sun / pc^2. However, the model also predicts that more generally on the scale of GMCs, X_CO is a strong function of metallicity, and depends on the column density and the interstellar UV flux. We show explicitly that neglecting these dependencies in observational estimates can strongly bias the inferred distribution of H2 column densities of molecular clouds to have a narrower and offset range compared to the true distribution. We find that when averaged on \sim kpc scales the X-factor depends only weakly on radiation field and column density, but is still a strong function of metallicity. The predicted metallicity dependence can be approximated as X_CO \sim Z^{-{\gamma}} with {\gamma} ~ 0.5 - 0.8.

We study the relationship between two major baryonic components in galaxy clusters, namely the stars in galaxies, and the ionized gas in the intracluster medium (ICM), using 94 clusters that span the redshift range 0-0.6. Accurately measured total and ICM masses from Chandra observations, and stellar masses derived from the Wide-field Infrared Survey Explorer and the Two-Micron All-Sky Survey allow us to trace the evolution of cluster baryon content in a self-consistent fashion. We find that, within r_{500}, the evolution of the ICM mass--total mass relation is consistent with the expectation of self-similar model, while there is no evidence for redshift evolution in the stellar mass--total mass relation. This suggests that the stellar mass and ICM mass in the inner parts of clusters evolve differently.

The tight relation of star formation with molecular gas indicated by observations and assumed in recent models implies that the efficiency with which galaxies convert their gas into stars depends on gas metallicity. This is because the abundance of molecular hydrogen is sensitive to the abundance of dust, which catalyzes the formation of H_2 and helps to shield it from dissociating radiation. In this study we point out that in the absence of significant pre-enrichment by Population III stars forming out of zero metallicity gas, such H_2-based star formation is expected to leave an imprint in the form of bi-modality in the metallicity distribution among dwarf galaxies and in the metallicity distribution of stars within individual galaxies. The bi-modality arises because when gas metallicity (and dust abundance) is low, formation of molecular gas is inefficient, the gas consumption time scale is long, and star formation and metal enrichment proceed slowly. When metallicity reaches a critical threshold value star formation and enrichment accelerate, which leads to rapid increase in both stellar mass and metallicity of galaxies. We demonstrate this process both using a simple analytical model and full cosmological simulations. In contrast, observed metallicity distributions of dwarf galaxies or stars within them are not bi-modal. We argue that this discrepancy points to substantial early stochastic pre-enrichment by population III stars to levels Z ~ 0.01 Z_sun in dense, star forming regions of early galaxies.

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 study the detailed structure of galaxies at redshifts z > 2 using cosmological simulations with improved modeling of the interstellar medium and star formation. The simulations follow the formation and dissociation of molecular hydrogen, and include star formation only in cold molecular gas. The molecular gas is more concentrated towards the center of galaxies than the atomic gas, and as a consequence, the resulting stellar distribution is very compact. For halos with total mass above 10^{11} Mo, the median half-mass radius of the stellar disks is 0.8 kpc at z = 3. The vertical structure of the molecular disk is much thinner than that of the atomic neutral gas. Relative to the non-radiative run, the inner regions of the dark matter halo change shape from prolate to mildly oblate and align with the stellar disk. However, we do not find evidence for a significant dark disk of dark matter around the stellar disk. The outer halo regions retain the orientation acquired during accretion and mergers, and are significantly misaligned with the inner regions. The radial profile of the dark matter halo contracts in response to baryon dissipation, establishing an approximately isothermal profile throughout most of the halo. This effect can be accurately described by a modified model of halo contraction. The angular momentum of a fixed amount of inner dark matter is approximately conserved over time, while in the dissipationless case most of it is transferred outward during mergers. The conservation of the dark matter angular momentum provides supporting evidence for the validity of the halo contraction model in a hierarchical galaxy formation process.

We present one of the most precise measurement to date of the spatial clustering of X-ray selected AGNs using a sample derived from the Chandra X-ray Observatory survey in the Bootes field. The real-space two-point correlation function over a redshift interval from z=0.17 to z~3 is well described by the power law, xi(r)=(r/r0)^-gamma, for comoving separations r<~20h^-1 Mpc. We find gamma=1.84+-0.12 and r0 consistent with no redshift trend within the sample (varying between r0=5.5+-0.6 h^-1 Mpc for <z>=0.37 and r0=6.9+-1.0 h^-1 Mpc for <z>=1.28). Further, we are able to measure the projections of the two-point correlation function both on the sky plane and in the line of sight. We use these measurements to show that the Chandra/Bootes AGNs are predominantly located at the centers of dark matter halos with the circular velocity Vmax>320 km/s or M_200 > 4.1e12 h^-1 Msun, and tend to avoid satellite galaxies in halos of this or higher mass. The halo occupation properties inferred from the clustering properties of Chandra/Bootes AGNs --- the mass scale of the parent dark matter halos, the lack of significant redshift evolution of the clustering length, and the low satellite fraction --- are broadly consistent with the Hopkins et al. scenario of quasar activity triggered by mergers of similarly-sized galaxies.

We study the bias and scatter in mass measurements of galaxy clusters resulting from fitting a spherically-symmetric Navarro, Frenk & White model to the reduced tangential shear profile measured in weak lensing observations. The reduced shear profiles are generated for ~10^4 cluster-sized halos formed in a LCDM cosmological N-body simulation of a 1 Gpc/h box. In agreement with previous studies, we find that the scatter in the weak lensing masses derived using this fitting method has irreducible contributions from the triaxial shapes of cluster-sized halos and uncorrelated large-scale matter projections along the line-of-sight. Additionally, we find that correlated large-scale structure within several virial radii of clusters contributes a smaller, but nevertheless significant, amount to the scatter. The intrinsic scatter due to these physical sources is ~20% for massive clusters, and can be as high as ~30% for group-sized systems. For current, ground-based observations, however, the total scatter should be dominated by shape noise from the background galaxies used to measure the shear. Importantly, we find that weak lensing mass measurements can have a small, ~5%-10%, but non-negligible amount of bias. Given that weak lensing measurements of cluster masses are a powerful way to calibrate cluster mass-observable relations for precision cosmological constraints, we strongly emphasize that a robust calibration of the bias requires detailed simulations which include more observational effects than we consider here. Such a calibration exercise needs to be carried out for each specific weak lensing mass estimation method, as the details of the method determine in part the expected scatter and bias. We present an iterative method for estimating mass M500c that can eliminate the bias for analyses of ground-based data.

A basic property of objects, like galaxies and halos that form in cosmological structure formation simulations, is their shape. Here, we critically investigate shape determination methods that are commonly used in the literature. It is found that using an enclosed integration volume and weight factors r^{-2} and r_{ell}^{-2} (elliptical radius) for the contribution of each particle or volume element in the shape tensor leads to biased axis ratios and smoothing of details when calculating the local shape as a function of distance from the center. To determine the local shape of matter distributions as a function of distance for well resolved objects (typically more than O(10^4) particles), we advocate a method that (1) uses an ellipsoidal shell (homoeoid) as an integration volume without any weight factors in the shape tensor and (2) removes subhalos.

We present analysis of the three-dimensional shape of intracluster gas in clusters formed in cosmological simulations of the Lambda-CDM cosmology and compare it to the shape of dark matter distribution and the shape of the overall isopotential surfaces. We find that in simulations with radiative cooling, star formation and stellar feedback (CSF), intracluster gas outside the cluster core is more spherical compared to non-radiative (NR) simulations, while in the core the gas in the CSF runs is more triaxial and has a distinctly oblate shape. The latter reflects the ongoing cooling of gas, which settles into a thick oblate ellipsoid as it loses thermal energy. The shape of the gas in the inner regions of clusters can therefore be a useful diagnostic of gas cooling. We find that gas traces the shape of the underlying potential rather well outside the core, as expected in hydrostatic equilibrium. At smaller radii, however, the gas and potential shapes differ significantly. In the CSF runs, the difference reflects the fact that gas is partly rotationally supported. Interestingly, we find that in NR simulations the difference between gas and potential shape at small radii is due to random gas motions, which make the gas distribution more spherical than the equipotential surfaces. Finally, we use mock Chandra X-ray maps to show that the differences in shapes observed in three-dimensional distribution of gas are discernible in the ellipticity of X-ray isophotes. Contrasting the ellipticities measured in simulated clusters against observations can therefore constrain the amount of cooling of the intracluster medium and the presence of random gas motions in cluster cores.

As emphasized by previous studies, proper treatment of the density fluctuation on the fundamental scale of a cosmological simulation volume - the "DC mode" - is critical for accurate modeling of spatial correlations on scales > 10% of simulation box size. We provide further illustration of the effects of the DC mode on the abundance of halos in small boxes and show that it is straightforward to incorporate this mode in cosmological codes that use the "supercomoving" variables. The equations governing evolution of dark matter and baryons recast with these variables are particularly simple and include the expansion factor, and hence the effect of the DC mode, explicitly only in the Poisson equation.

We examine the importance of secular stellar mass loss for fueling ongoing star formation in disk galaxies during the late stages of their evolution. For a galaxy of a given stellar mass, we calculate the total mass loss rate of its entire stellar population using star formation histories derived from the observed evolution of the M*-star formation rate relation, along with the predictions of standard stellar evolution models for stellar mass loss for a variety of initial stellar mass functions. Our model shows that recycled gas from stellar mass loss can provide most or all of the fuel required to sustain the current level of star formation in late type galaxies. Stellar mass loss can therefore remove the tension between the low gas infall rates that are derived from observations and the relatively rapid star formation occurring in disk galaxies. For galaxies where cold gas infall rates have been estimated, we demonstrate explicitly that stellar mass loss can account for most of the deficit between their star formation and infall rates.

It is a well established empirical fact that the surface density of the star formation rate, Sigma_SFR, strongly correlates with the surface density of molecular hydrogen, Sigma_H2, at least when averaged over large (~kpc) scales. Much less is known, however, if (and how) the Sigma_SFR-Sigma_H2 relation depends on environmental parameters, such as the metallicity or the UV radiation field in the interstellar medium (ISM). Furthermore, observations indicate that the scatter in the Sigma_SFR-Sigma_H2 relation increases rapidly with decreasing averaging scale. How the scale-dependent scatter is generated and how one recovers a tight ~ kpc scale Sigma_SFR-Sigma_H2 relation in the first place is still largely debated. Here, these questions are explored with hydrodynamical simulations that follow the formation and destruction of H2, include radiative transfer of UV radiation, and resolve the ISM on ~60 pc scales. We find that within the considered range of H2 surface densities (10-100 Msun/pc^2) the Sigma_SFR-Sigma_H2 relation is steeper in environments of low metallicity and/or high radiation fields (compared to the Galaxy), that the star formation rate at a given H2 surface density is larger, and the scatter is increased. Deviations from a "universal" Sigma_SFR-Sigma_H2 relation should be particularly relevant for high redshift galaxies or for low-metallicity dwarfs at z~0. We also find that the use of time-averaged SFRs produces a large, scale dependent scatter in the Sigma_SFR-Sigma_H2 relation. Given the plethora of observational data expected from upcoming surveys such as ALMA the scale-scatter relation may indeed become a valuable tool for determining the physical mechanisms connecting star formation and H2 formation.

The friends-of-friends algorithm (hereafter, FOF) is a percolation algorithm which is routinely used to identify dark matter halos from N-body simulations. We use results from percolation theory to show that the boundary of FOF halos does not correspond to a single density threshold but to a range of densities close to a critical value that depends upon the linking length parameter, b. We show that for the commonly used choice of b = 0.2, this critical density is equal to 81.62 times the mean matter density. Consequently, halos identified by the FOF algorithm enclose an average overdensity which depends on their density profile (concentration) and therefore changes with halo mass contrary to the popular belief that the average overdensity is ~180. We derive an analytical expression for the overdensity as a function of the linking length parameter b and the concentration of the halo. Results of tests carried out using simulated and actual FOF halos identified in cosmological simulations show excellent agreement with our analytical prediction. We also find that the mass of the halo that the FOF algorithm selects crucially depends upon mass resolution. We find a percolation theory motivated formula that is able to accurately correct for the dependence on number of particles for the mock realizations of spherical and triaxial Navarro-Frenk-White halos. However, we show that this correction breaks down when applied to the real cosmological FOF halos due to presence of substructures. Given that abundance of substructure depends on redshift and cosmology, we expect that the resolution effects due to substructure on the FOF mass and halo mass function will also depend on redshift and cosmology and will be difficult to correct for in general. Finally, we discuss the implications of our results for the universality of the mass function.

We calculate the probability distribution function (PDF) of the expected annihilation luminosities of dark matter subhalos as a function of subhalo mass and distance from the Galactic center using a semi-analytical model of halo evolution. We find that the PDF of luminosities is relatively broad, exhibiting a spread of as much as an order of magnitude at fixed subhalo mass and halo-centric distance. The luminosity PDF allows for simple construction of mock samples of gamma-ray luminous subhalos and assessment of the variance in among predicted gamma-ray signals from dark matter annihilation. Other applications include quantifying the variance among the expected luminosities of dwarf spheroidal galaxies, assessing the level at which dark matter annihilation can be a contaminant in the expected gamma-ray signal from other astrophysical sources, as well as estimating the level at which nearby subhalos can contribute to the antimatter flux.

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 present optical emission-line spectra for outlying HII regions in the extended neutral gas disk surrounding the blue compact dwarf galaxy NGC 2915. Using a combination of strong-line R23 and direct oxygen abundance measurements, we report a flat, possibly increasing, metallicity gradient out to 1.2 times the Holmberg radius. We find the outer-disk of NGC 2915 to be enriched to a metallicity of 0.4 Z_solar. An analysis of the metal yields shows that the outer disk of NGC 2915 is overabundant for its gas fraction, while the central star-foming core is similarly under-abundant for its gas fraction. Star formation rates derived from very deep ~14 ks GALEX FUV exposures indicate that the low-level of star formation observed at large radii is not sufficient to have produced the measured oxygen abundances at these galactocentric distances. We consider 3 plausible mechanisms that may explain the metal-enriched outer gaseous disk of NGC 2915: radial redistribution of centrally generated metals, strong galactic winds with subsequent fallback, and galaxy accretion. Our results have implications for the physical origin of the mass-metallicity relation for gas-rich dwarf galaxies.

Dwarf galaxies provide opportunities for drawing inferences about the processes in the early universe by observing our "cosmological backyard"-the Local Group and its vicinity. This special issue of the open-access journal Advances in Astronomy is a snapshot of the current state of the art of dwarf-galaxy cosmology.

We present a detailed description of a phenomenological H2 formation model and local star formation prescription based on the density of molecular (rather than total) gas. Such approach allows us to avoid the arbitrary density and temperature thresholds typically used in star formation recipes. We present results of the model based on realistic cosmological simulations of high-z galaxy formation for a grid of numerical models with varied dust-to-gas ratios and interstellar far UV (FUV) fluxes. Our results show that both the atomic-to-molecular transition on small, ~10 pc scales and the Kennicutt-Schmidt (KS) relation on ~kpc scales are sensititive to the dust-to-gas ratio and the FUV flux. The atomic-to-molecular transition as a function of gas density or column density has a large scatter but is rather sharp and shifts to higher densities with decreasing dust-to-gas ratio and/or increasing FUV flux. Consequently, star formation is concentrated to higher gas surface density regions, resulting in steeper slope and lower amplitude of the KS relation at a given gas surface density, in less dusty and/or higher FUV flux environments. These trends should have a particularly strong effect on the evolution of low-mass, low surface brightness galaxies which typically have low dust content and anemic star formation, but are also likely to be important for evolution of the Milky Way-sized systems. We parameterize the dependencies observed in our simulations in convenient fitting formulae, which can be used to model the dependence of the KS relation on the dust-to-gas ratio and FUV flux in semi-analytic models and in cosmological simulations that do not include radiative transfer and H2 formation.

We present results of self-consistent, high-resolution cosmological simulations of galaxy formation at z~3. The simulations employ recently developed recipe for star formation based on the local abundance of molecular hydrogen, which is tracked self-consistently during the course of simulation. The phenomenological H2 formation model accounts for the effects of dissociating UV radiation of stars in each galaxy, as well as self-shielding and shielding of H2 by dust, and therefore allows us to explore effects of lower metallicities and higher UV fluxes prevalent in high redshift galaxies on their star formation. We compare stellar masses, metallicities, and star formation rates of the simulated galaxies to available observations of the Lyman Break Galaxies (LBGs) and find a reasonable agreement. We find that the Kennicutt-Schmidt (KS) relation exhibited by our simulated galaxies at z~3 is substantially steeper and has a lower amplitude than the z=0 relation at Sigma_gas < 100 Msun/pc^2. The predicted relation, however, is consistent with existing observational constraints for the z~3 Damped Lyman $\alpha$ (DLA) and LBGs. Our tests show that the main reason for the difference from the local KS relation is lower metallicity of the ISM in high redshift galaxies. We discuss several implications of the metallicity-dependence of the KS relation for galaxy evolution and interpretation of observations.

We investigate effects of baryon dissipation on the dark matter virial scaling relation between total mass and velocity dispersion and the velocity bias of galaxies in groups and clusters using self-consistent cosmological simulations. We show that the baryon dissipation increases the velocity dispersion of dark matter within the virial radius by 5% - 10%. The effect is mainly driven by the change in density and gravitational potential in inner regions of cluster, and it is larger in lower mass systems where gas cooling and star formation are more efficient. We also show that the galaxy velocity bias depends on how galaxies are selected. Galaxies selected based on their stellar mass exhibit no velocity bias, while galaxies selected based on their total mass show positive bias of ~10%, consistent with previous results based on collisionless dark matter- only simulations. We further find that observational estimates of galaxy velocity dispersion are unbiased with respect to the velocity dispersion of dark matter, provided galaxies are selected using their stellar masses and and their velocity dispersions are computed with more than twenty most massive galaxies. Velocity dispersions estimated with fewer galaxies, on the other hand, can lead to significant underestimate of dynamical masses. Results presented in this paper should be useful in interpretating high-redshift groups and clusters as well as cosmological constraints derived from upcoming optical cluster surveys.

Eulerian hydrodynamical simulations are a powerful and popular tool for modeling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.

We present analysis of bulk and random gas motions in the intracluster medium using high-resolution Eulerian cosmological simulations of sixteen simulated clusters, including both very relaxed and unrelaxed systems and spanning a virial mass range of 5*10^13 - 2*10^15 Msun/h. We investigate effects of the residual subsonic gas motions on the hydrostatic estimates of mass profiles and concentrations of galaxy clusters. In agreement with previous studies we find that the gas motions contribute up to ~ 5%-15% of the total pressure support in relaxed clusters with contribution increasing with cluster-centric radius. The fractional pressure support is higher in unrelaxed systems. This contribution would not be accounted for in hydrostatic estimates of the total mass profile and would lead to systematic underestimate of mass. We demonstrate that total mass can be recovered accurately if pressure due to gas motions measured in simulations is explicitly taken into account in the equation of hydrostatic equilibrium. Given that the underestimate of mass is increasing at larger radii, where gas is less relaxed and contribution of gas motions to pressure is larger, the total density profile derived from hydrostatic analysis is more concentrated than the true profile. This may at least partially explain some high values of concentrations of clusters estimated from hydrostatic analysis of X-ray data.

We study the evolution of phase-space density during the hierarchical structure formation of LCDM halos. We compute both a spherically-averaged surrogate for phase-space density (Q) and the coarse-grained distribution function f(x,v) for dark matter particles that lie within~2 virial radii of four Milky-Way-sized dark matter halos. The estimated f(x,v) spans over four decades at any radius. Dark matter particles that end up within two virial radii of a Milky-Way-sized DM halo at $z=0$ have an approximately Gaussian distribution in log(f) at early redshifts, but the distribution becomes increasingly skewed at lower redshifts. The value corresponding to the peak of the Gaussian decreases as the evolution progresses and is well described by a power-law in (1+z). The highest values of f are found at the centers of dark matter halos and subhalos, where f can be an order of magnitude higher than in the center of the main halo. The power-law Q(r) profile likely reflects the distribution of entropy (K = sigma^2/rho^{2/3} \propto r^{1.2}), which dark matter acquires as it is accreted onto a growing halo. The estimated f(x, v), on the other hand, exhibits a more complicated behavior. Although the median coarse-grained phase-space density profile F(r) can be approximated by a power-law in the inner regions of halos and at larger radii the profile flattens significantly. This is because phase-space density averaged on small scales is sensitive to the high-f material associated with surviving subhalos, as well as relatively unmixed material (probably in streams) resulting from disrupted subhalos, which contribute a sizable fraction of matter at large radii. (ABRIDGED)

We review recent progress in the description of the formation and evolution of galaxy clusters in a cosmological context by using numerical simulations. We focus our presentation on the comparison between simulated and observed X-ray properties, while we will also discuss numerical predictions on properties of the galaxy population in clusters. Many of the salient observed properties of clusters, such as X-ray scaling relations, radial profiles of entropy and density of the intracluster gas, and radial distribution of galaxies are reproduced quite well. In particular, the outer regions of cluster at radii beyond about 10 per cent of the virial radius are quite regular and exhibit scaling with mass remarkably close to that expected in the simplest case in which only the action of gravity determines the evolution of the intra-cluster gas. However, simulations generally fail at reproducing the observed cool-core structure of clusters: simulated clusters generally exhibit a significant excess of gas cooling in their central regions, which causes an overestimate of the star formation and incorrect temperature and entropy profiles. The total baryon fraction in clusters is below the mean universal value, by an amount which depends on the cluster-centric distance and the physics included in the simulations, with interesting tensions between observed stellar and gas fractions in clusters and predictions of simulations. Besides their important implications for the cosmological application of clusters, these puzzles also point towards the important role played by additional physical processes, beyond those already included in the simulations. We review the role played by these processes, along with the difficulty for their implementation, and discuss the outlook for the future progress in numerical modeling of clusters.

A decade ago cosmological simulations of increasingly higher resolution were used to demonstrate that virialized regions of Cold Dark Matter (CDM) halos are filled with a multitude of dense, gravitationally-bound clumps. These dark matter subhalos are central regions of halos that survived strong gravitational tidal forces and dynamical friction during the hierarchical sequence of merging and accretion via which the CDM halos form. Comparisons with observations revealed that there is a glaring discrepancy between abundance of subhalos and luminous satellites of the Milky Way and Andromeda as a function of their circular velocity or bound mass within a fixed aperture. This large discrepancy, which became known as the ``substructure'' or the ``missing satellites'' problem, begs for an explanation. In this paper I review the progress made during the last several years both in quantifying the problem and in exploring possible scenarios in which it could be accommodated and explained in the context of galaxy formation in the framework of the CDM paradigm of structure formation. In particular, I show that the observed luminosity function, radial distribution, and the remarkable similarity of the inner density profiles of luminous satellites can be understood within hierarchical CDM framework using a simple model in which efficiency of star formation monotonically decreases with decreasing virial mass satellites had before their accretion without any actual sharp galaxy formation threshold.

We use dissipationless N-body simulations to investigate the evolution of the true coarse-grained phase-space density distribution f(x,v) in equal-mass mergers between dark matter (DM) halos. The halo models are constructed with various asymptotic power-law indices ranging from steep cusps to core-like profiles and we employ the phase-space density estimator ``Enbid'' developed by Sharma & Steinmetz to compute f(x,v). The adopted force resolution allows robust phase-space density profile estimates in the inner ~1% of the virial radii of the simulated systems. We confirm that mergers result in a decrease of the coarse-grained phase-space density in accordance with expectations from Mixing Theorems for collisionless systems. We demonstrate that binary mergers between identical DM halos produce remnants that retain excellent memories of the inner slopes and overall shapes of the phase-space density distribution of their progenitors. The robustness of the phase-space density profiles holds for a range of orbital energies, and a variety of encounter configurations including sequences of several consecutive merger events, designed to mimic hierarchical merging, and collisions occurring at different cosmological epochs. If the progenitor halos are constructed with appreciably different asymptotic power-law indices, we find that the inner slope and overall shape of the phase-space density distribution of the remnant are substantially closer to that of the initial system with the steepest central density cusp. These results explicitly demonstrate that mixing is incomplete in equal-mass mergers between DM halos, as it does not erase memory of the progenitor properties. Our results also confirm the recent analytical predictions of Dehnen (2005) regarding the preservation of merging self-gravitating central density cusps.

(Abridged) We perform dissipationless N-body simulations to elucidate the dynamical response of thin disks to bombardment by cold dark matter (CDM) substructure. Our method combines (1) cosmological simulations of the formation of Milky Way (MW)-sized CDM halos to derive the properties of substructure and (2) controlled numerical experiments of consecutive subhalo impacts onto an initially-thin, fully-formed MW type disk galaxy. The present study is the first to account for the evolution of satellite populations over cosmic time in such an investigation of disk structure. We find that accretions of massive subhalos onto the central regions of host halos, where the galactic disks reside, since z~1 should be common. One host halo accretion history is used to initialize the controlled simulations of satellite-disk encounters. We show that these accretion events severely perturb the thin galactic disk and produce a wealth of distinctive dynamical signatures on its structure and kinematics. These include (1) considerable thickening and heating at all radii, with the disk thickness and velocity ellipsoid nearly doubling at the solar radius; (2) prominent flaring associated with an increase in disk thickness greater than a factor of 4 in the disk outskirts; (3) surface density excesses at large radii, beyond ~5 disk scale lengths, resembling those of observed antitruncated disks; (4) lopsidedness at levels similar to those measured in observational samples of disk galaxies; and (5) substantial tilting. The interaction with the most massive subhalo drives the disk response while subsequent bombardment is much less efficient at disturbing the disk. We conclude that substructure-disk encounters of the kind expected in the LCDM paradigm play a significant role in setting the structure of disk galaxies and driving galaxy evolution.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

For over four decades, synchrotron-radiating sources have played a series of pathfinding roles in the study of galaxy clusters and large scale structure. Such sources are uniquely sensitive to the turbulence and shock structures of large-scale environments, and their cosmic rays and magnetic fields often play important dynamic and thermodynamic roles. They provide essential complements to studies at other wavebands. Over the next decade, they will fill essential gaps in both cluster astrophysics and the cosmological growth of structure in the universe, especially where the signatures of shocks and turbulence, or even the underlying thermal plasma itself, are otherwise undetectable. Simultaneously, synchrotron studies offer a unique tool for exploring the fundamental question of the origins of cosmic magnetic fields. This work will be based on the new generation of m/cm-wave radio telescopes now in construction, as well as major advances in the sophistication of 3-D MHD simulations.

We are learning much about how structure forms, in particular how clusters as nodes in the cosmic web evolve and accrete matter, and about the physical processes within these objects. In the next decade, the study of clusters will enable us to tackle important questions regarding the nature of Dark Matter and Dark Energy, how clusters co-evolve with super-massive black holes at their centers, and to advance our knowledge about fundamental plasma astrophysics. This science white paper outlines the key questions and research opportunities in cluster astrophysics that are emerging in the coming decade and beyond, and serves as an overview to other cluster related white papers.

We describe a phenomenological model for molecular hydrogen formation suited for applications in galaxy formation simulations, which includes on-equilibrium formation of molecular hydrogen on dust and approximate treatment of both its self-shielding and shielding by dust from the dissociating UV radiation. The model is applicable in simulations in which individual star forming regions - the giant molecular complexes - can be identified (resolution of tens of pc) and their mean internal density estimated reliably, even if internal structure is not resolved. In agreement with previous studies, calculations based on our model show that the transition from atomic to fully molecular phase depends primarily on the metallicity, which we assume is directly related to the dust abundance, and clumpiness of the interstellar medium. The clumpiness simply boosts the formation rate of molecular hydrogen, while dust serves both as a catalyst of molecular hydrogen formation and as an additional shielding from dissociating UV radiation. The upshot is that it is difficult to form fully-shielded giant molecular clouds while gas metallicity is low. However, once the gas is enriched to Z ~ 0.01-0.1 solar, the subsequent star formation and enrichment can proceed at a much faster rate. This may keep star formation efficiency in the low-mass, low-metallicity progenitors of galaxies very low for a certain period of time with the effect similar to a strong "feedback" mechanism. [abridged]

Groups and clusters of galaxies occupy a special position in the hierarchy of large-scale cosmic structures because they are the largest and the most massive (from ~10^13 Msun to over 10^15 Msun) objects in the universe that have had time to undergo gravitational collapse. The large masses of clusters imply that their contents have been accreted from regions of ~8-40 comoving Mpc in size and should thus be representative of the mean matter content of the universe. During the next decade sensitive multi-wavelength observations should be able to map the radial distributions of all main mass components (stars, cold, warm, and hot gas and total mass) at z<~ 1 out to the virial radius. At the same time, comparative studies of real and simulated cluster samples sould allow us to use clusters as veritable astrophysical laboratories for studying galaxy formation, as well as testing our theoretical models of structure formation and underlying assumptions about fundamental physics governing the universe.

The hot, diffuse gas that fills the largest overdense structures in the Universe -- clusters of galaxies and a web of giant filaments connecting them -- provides us with tools to address a wide array of fundamental astrophysical and cosmological questions via observations in the X-ray band. Clusters are sensitive cosmological probes. To utilize their full potential for precision cosmology in the following decades, we must precisely understand their physics -- from their cool cores stirred by jets produced by the central supermassive black hole (itself fed by inflow of intracluster gas), to their outskirts, where the infall of intergalactic medium (IGM) drives shocks and accelerates cosmic rays. Beyond the cluster confines lies the virtually unexplored warm IGM, believed to contain most of the baryonic matter in the present-day Universe. As a depository of all the matter ever ejected from galaxies, it carries unique information on the history of energy and metal production in the Universe. Currently planned major observatories, such as Astro-H and IXO, will make deep inroads into these areas, but to see the most interesting parts of the picture will require an almost science-fiction-grade facility with tens of m^2 of effective area, subarcsecond angular resolution, a matching imaging calorimeter and a super high-dispersion spectrograph, such as Generation-X.

Our heuristic understanding of the abundance of dark matter halos centers around the concept of a density threshold, or "barrier", for gravitational collapse. If one adopts the ansatz that regions of the linearly evolved density field smoothed on mass scale M with an overdensity that exceeds the barrier will undergo gravitational collapse into halos of mass M, the corresponding abundance of such halos can be estimated simply as a fraction of the mass density satisfying the collapse criterion divided by the mass M. The key ingredient of this ansatz is therefore the functional form of the collapse barrier as a function of mass M or, equivalently, of the variance sigma^2(M). Several such barriers based on the spherical, Zel'dovich, and ellipsoidal collapse models have been extensively discussed. Using large scale cosmological simulations, we show that the relation between the linear overdensity and the mass variance for regions that collapse to form halos by the present epoch resembles expectations from dynamical models of ellipsoidal collapse. However, we also show that using such a collapse barrier with the excursion set ansatz predicts a halo mass function inconsistent with that measured directly in cosmological simulations. This inconsistency demonstrates a failure of the excursion set ansatz as a physical model for halo collapse. We discuss implications of our results for understanding the collapse epoch for halos as a function of mass, and avenues for improving consistency between analytical models for the collapse epoch and the results of cosmological simulations.

We discuss the measurements of the galaxy cluster mass functions at z=~0.05 and z=~0.5 using high-quality Chandra observations of samples derived from the ROSAT PSPC All-Sky and 400deg^2 surveys. We provide a full reference for the data analysis procedures, present updated calibration of relations between the total cluster mass and its X-ray indicators (T_X, Mgas, and Y_X) based on a subsample of low-z relaxed clusters, and present a first measurement of the evolving L_X-Mtot relation (with Mtot estimated from Y_X) obtained from a well-defined statistically complete cluster sample and with appropriate corrections for the Malmquist bias applied. Finally, we present the derived cluster mass functions, estimate the systematic uncertainties in this measurement, and discuss the calculation of the likelihood function. We confidently measure the evolution in the cluster comoving number density at a fixed mass threshold, e.g., by a factor of 5.0 +- 1.2 at M_500=2.5e14 h^-1 Msun between z=0 and 0.5. This evolution reflects the growth of density perturbations and can be used for the cosmological constraints complementing those from the distance-redshift relation.

Chandra observations of large samples of galaxy clusters detected in X-rays by ROSAT provide a new, robust determination of the cluster mass functions at low and high redshifts. Statistical and systematic errors are now sufficiently small, and the redshift leverage sufficiently large for the mass function evolution to be used as a useful growth of structure based dark energy probe. In this paper, we present cosmological parameter constraints obtained from Chandra observations of 36 clusters with <z>=0.55 derived from 400deg^2 ROSAT serendipitous survey and 49 brightest z=~0.05 clusters detected in the All-Sky Survey. Evolution of the mass function between these redshifts requires Omega_Lambda>0 with a ~5sigma significance, and constrains the dark energy equation of state parameter to w0=-1.14+-0.21, assuming constant w and flat universe. Cluster information also significantly improves constraints when combined with other methods. Fitting our cluster data jointly with the latest supernovae, WMAP, and baryonic acoustic oscillations measurements, we obtain w0=-0.991+-0.045 (stat) +-0.039 (sys), a factor of 1.5 reduction in statistical uncertainties, and nearly a factor of 2 improvement in systematics compared to constraints that can be obtained without clusters. The joint analysis of these four datasets puts a conservative upper limit on the masses of light neutrinos, Sum m_nu<0.33 eV at 95% CL. We also present updated measurements of Omega_M*h and sigma_8 from the low-redshift cluster mass function.

(Abridged) We conduct a series of high-resolution, dissipationless N-body simulations to investigate the cumulative effect of substructure mergers onto thin disk galaxies in the context of the LCDM paradigm of structure formation. Our simulation campaign is based on a hybrid approach. Substructure properties are culled directly from cosmological simulations of galaxy-sized cold dark matter (CDM) halos. In contrast to what can be inferred from statistics of the present-day substructure populations, accretions of massive subhalos onto the central regions of host halos, where the galactic disk resides, since z~1 should be common occurrences. One host halo merger history is subsequently used to seed controlled numerical experiments of repeated satellite impacts on an initially-thin Milky Way-type disk galaxy. We show that these accretion events produce several distinctive observational signatures in the stellar disk including: a ring-like feature in the outskirts; a significant flare; a central bar; and faint filamentary structures that (spuriously) resemble tidal streams. The final distribution of disk stars exhibits a complex vertical structure that is well-described by a standard ``thin-thick'' disk decomposition. We conclude that satellite-disk encounters of the kind expected in LCDM models can induce morphological features in galactic disks that are similar to those being discovered in the Milky Way, M31, and in other disk galaxies. These results highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution. Upcoming galactic structure surveys and astrometric satellites may be able to distinguish between competing cosmological models by testing whether the detailed structure of galactic disks is as excited as predicted by the CDM paradigm.

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.

(ABRIDGED) We use hydrodynamical simulations of disk galaxies to study relations between star formation and properties of the molecular interstellar medium (ISM). We implement a model for the ISM that includes low-temperature (T<10^4K) cooling, directly ties the star formation rate to the molecular gas density, and accounts for the destruction of H2 by an interstellar radiation field from young stars. We demonstrate that the ISM and star formation model simultaneously produces a spatially-resolved molecular-gas surface density Schmidt-Kennicutt relation of the form Sigma_SFR \propto Sigma_Hmol^n_mol with n_mol~1.4 independent of galaxy mass, and a total gas surface density -- star formation rate relation Sigma_SFR \propto Sigma_gas^n_tot with a power-law index that steepens from n_tot~2 for large galaxies to n_tot>~4 for small dwarf galaxies. We show that deviations from the disk-averaged Sigma_SFR \propto Sigma_gas^1.4 correlation determined by Kennicutt (1998) owe primarily to spatial trends in the molecular fraction f_H2 and may explain observed deviations from the global Schmidt-Kennicutt relation.

Using a hydrodynamic adaptive mesh refinement code, we simulate the growth and evolution of a galaxy, which could potentially host a supermassive black hole, within a cosmological volume. Reaching a dynamical range in excess of 10 million, the simulation follows the evolution of the gas structure from super-galactic scales all the way down to the outer edge of the accretion disk. Here, we focus on global instabilities in the self-gravitating, cold, turbulence-supported, molecular gas disk at the center of the model galaxy, which provide a natural mechanism for angular momentum transport down to sub-pc scales. The gas density profile follows a power-law scaling as r^-8/3, consistent with an analytic description of turbulence in a quasi-stationary circumnuclear disk. We analyze the properties of the disk which contribute to the instabilities, and investigate the significance of instability for the galaxy's evolution and the growth of a supermassive black hole at the center.

We present detailed comparisons of the intracluster medium (ICM) in cosmological Eulerian cluster simulations with deep Chandra observations of nearby relaxed clusters. To assess the impact of galaxy formation, we compare two sets of simulations, one performed in the non-radiative regime and another with radiative cooling and several physical processes critical to various aspects of galaxy formation: star formation, metal enrichment and stellar feedback. We show that the observed ICM properties outside cluster cores are well-reproduced in the simulations that include cooling and star formation, while the non-radiative simulations predict an overall shape of the ICM profiles inconsistent with observations. In particular, we find that the ICM entropy in our runs with cooling is enhanced to the observed levels at radii as large as half of the virial radius. We also find that outside cluster cores entropy scaling with the mean ICM temperature in both simulations and Chandra observations is consistent with being self-similar within current error bars. We find that the pressure profiles of simulated clusters are also close to self-similar and exhibit little cluster-to-cluster scatter. The X-ray observable-total mass relations for our simulated sample agree with the Chandra measurements to \~10%-20% in normalization. We show that this systematic difference could be caused by the subsonic gas motions, unaccounted for in X-ray hydrostatic mass estimates. The much improved agreement of simulations and observations in the ICM profiles and scaling relations is encouraging and the existence of tight relations of X-ray observables, such as Yx, and total cluster mass and the simple redshift evolution of these relations hold promise for the use of clusters as cosmological probes.

Nearby dwarf galaxies exhibit tight correlations between their global stellar and dynamical properties, such as circular velocity, mass-to-light ratio, stellar mass, surface brightness, and metallicity. Such correlations have often been attributed to gas or metal-rich outflows driven by supernova energy feedback to the interstellar medium. We use high-resolution cosmological simulations of high-redshift galaxies with and without energy feedback, as well as analytic modeling, to investigate whether the observed correlations can arise without supernova-driven outflows. We find that the simulated dwarf galaxies exhibit correlations similar to those observed as early as z~10, regardless of whether supernova feedback is included. We also show that the correlations can be well reproduced by our analytic model that accounts for realistic gas inflow but assumes no outflows, and star formation rate obeying the Kennicutt-Schmidt law with a critical density threshold. We argue that correlations in simulated galaxies arise due to the increasingly inefficient conversion of gas into stars in low-mass dwarf galaxies rather than supernova-driven outflows. We also show that the decrease of the observed effective yield in low-mass objects, often used as an indicator of gas and metal outflows, can be reasonably reproduced in our simulations without outflows. We show that this trend can arise if a significant fraction of metals in small galaxies is spread to the outer regions of the halo outside the stellar extent via mixing. In this case the effective yield can be significantly underestimated if only metals within the stellar radius are taken into account. Measurements of gas metallicity in the outskirts of gaseous disks of dwarfs would thus provide a key test of such explanation.

We model the escape of ionizing radiation from high-redshift galaxies using high-resolution Adaptive Mesh Refinement N-body + hydrodynamics simulations. Our simulations include time-dependent and spatially-resolved transfer of ionizing radiation in three dimensions, including effects of dust absorption. For galaxies of total mass M > 10^11 Msun and star formation rates SFR ~ 1-5 Msun/yr, we find angular averaged escape fractions of 0.01-0.03 over the entire redshift interval studied (3<z<9). In addition, we find that the escape fraction varies by more than an order of magnitude along different lines-of-sight within individual galaxies, from the largest values near galactic poles to the smallest along the galactic disk. The escape fraction declines steeply at lower masses and SFR. We show that the low values of escape fractions are due to a small fraction of young stars located just outside the edge of HI disk. We compare our predicted escape fraction of ionizing photons with previous results, and find a general agreement with both other simulation results and available direct detection measurements at z ~ 3. We also compare our simulations with a novel method to estimate the escape fraction in galaxies from the observed distribution of neutral hydrogen column densities along the lines of sights to long duration gamma-ray bursts. Using this method we find escape fractions of the GRB host galaxies of 2-3%, consistent with our theoretical predictions. [abridged]

(Abridged) We present a semi-analytic model for Cold Dark Matter halo substructure that can be used as a framework for studying galaxy formation and an ingredient in halo models of galaxy clustering. We perform a comprehensive comparison of the model to the results of a suite of high-resolution cosmological simulations. The comparisons reveal that subhalo statistics, such as velocity and mass functions, radial distributions, and the halo occupation distributions, agree well over three orders of magnitude in host halo mass at various redshifts. Both in the simulations and in our model, the radial distributions of subhalos are significantly shallower than that of the dark matter density. The abundance of subhalos in a host is set by competition between destruction and new accretion. High mass halos and halos at high redshift tend to host more subhalos because the subhalos have, on average, been accreted more recently. Similarly, at fixed mass, halos that formed more recently host more subhalos. Observed "fossil groups" may represent an extreme tail of this correlation. We find a related correlation between host halo concentration and subhalo abundance at fixed host mass, Nsat ~ c^-a, where a changes with redshift and host-to-subhalo mass ratio. Lastly, we use our model to populate host halos in one of our high-resolution cosmological simulations, replacing the subhalos resolved in the simulation with subhalos computed according to our model. We show that the resulting 2pt correlation function of such a "hybrid" halo ensemble is indistinguishable from that measured directly in the simulation. This supports one of the key tenets of the standard halo model -- the assumption that the halo occupation distribution is statistically independent of host halo environment.

We study the importance of baryonic physics on predictions of the matter power spectrum as it is relevant for forthcoming weak lensing surveys. We quantify the impact of baryonic physics using a set of three cosmological numerical simulations. Each simulation has the same initial density field, but models a different set of physical processes. The first simulation evolves the density field using gravity alone, the second includes non-radiative gasdynamics, and the third includes radiative heating and cooling of baryons, star formation, and supernova feedback. We find that baryonic processes alter predictions for the matter power spectrum significantly relative to models that include only gravitational interactions. Our results imply that future weak lensing experiments such as LSST and SNAP will be very sensitive to the poorly-understood physics governing the nonlinear evolution of the baryonic component of the universe. The net effect is significantly larger in the case of the model with cooling and star formation, in which case our results imply that contemporary surveys such as the CFHT Wide survey may also be sensitive to baryonic processes. In particular, this effect could be important for forecasts of the constraining power of future surveys if information from scales larger than l ~ 1000 is included in the analysis. We find that deviations are caused primarily by the rearrangement of matter within individual dark matter halos relative to the gravity-only case, rather than a large-scale rearrangement of matter. Consequently, we propose a simple model, based on the phenomenological halo model of dark matter clustering, for baryonic effects that can be used to aid in the interpretation of forthcoming weak lensing data.

We use a set of simulation-based models for the dissipationless evolution of galaxies since z=1 to constrain the fate of accreted satellites embedded in dark matter subhalos. These models assign stellar mass to dark matter halos at z=1 by relating the observed galaxy stellar mass function (GSMF) to the halo+subhalo mass function monotonically. The evolution of the stellar mass content is then followed using halo merger trees extracted from N-body simulations. Our models are differentiated only in the fate assigned to satellite galaxies once subhalos, within which satellites are embedded, disrupt. These models are confronted with the observed evolution in the massive end of the GSMF, the z~0 brightest cluster galaxy (BCG)-cluster mass relation, and the combined BCG and intracluster light (ICL) luminosity distribution -- all observables expected to evolve approximately dissipationlessly since z=1. The combined observational constraints favor a model in which the vast majority (>80%) of satellite stars from disrupted subhalos go into the ICL. Conversely, models that leave behind a significant population of satellite galaxies once the subhalo has disrupted are strongly disfavored, as are models that put a significant fraction of satellite stars into the BCG. Our results show that observations of the ICL provide useful and unique constraints on models of galaxy merging and the dissipationless evolution of galaxies in groups and clusters.

We analyze N-body simulations of halo mergers to investigate the mechanisms responsible for driving mixing in phase-space and the evolution to dynamical equilibrium. We focus on mixing in energy and angular momentum and show that mixing occurs in step-like fashion following pericenter passages of the halos. This makes mixing during a merger unlike other well known mixing processes such as phase mixing and chaotic mixing whose rates scale with local dynamical time. We conclude that the mixing process that drives the system to equilibrium is primarily a response to energy and angular momentum redistribution that occurs due to impulsive tidal shocking and dynamical friction rather than a result of chaotic mixing in a continuously changing potential. We also analyze the merger remnants to determine the degree of mixing at various radii by monitoring changes in radius, energy and angular momentum of particles. We confirm previous findings that show that the majority of particles retain strong memory of their original kinetic energies and angular momenta but do experience changes in their potential energies owing to the tidal shocks they experience during pericenter passages. Finally, we show that a significant fraction of mass (~ 40%) in the merger remnant lies outside its formal virial radius and that this matter is ejected roughly uniformly from all radii outside the inner regions. This highlights the fact that mass, in its standard virial definition, is not additive in mergers. We discuss the implications of these results for our understanding of relaxation in collisionless dynamical systems.

X-ray observations of galaxy clusters potentially provide powerful cosmological probes if systematics due to our incomplete knowledge of the intracluster medium (ICM) physics are understood and controlled. In this paper, we study the effects of galaxy formation on the properties of the ICM and X-ray observable-mass relations using high-resolution self-consistent cosmological simulations of galaxy clusters and comparing their results with recent Chandra X-ray observations. We show that despite complexities of their formation and uncertainties in their modeling, clusters of galaxies both in observations and numerical simulations are remarkably regular outside of their cores, which holds great promise for their use as cosmological probes.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.

Using ray tracing for simple analytic profiles, we demonstrate that the lensing cross section for producing giant arcs has distinct contributions due to arcs formed through image distortion only, and arcs form from the merging of two or three images. We investigate the dependence of each of these contributions on halo ellipticity and on the slope of the density profile, and demonstrate that at fixed Einstein radius, the lensing cross section increases as the halo profile becomes steeper. We then compare simulations with and without baryonic cooling of the same cluster for a sample of six clusters, and demonstrate that cooling can increase the overall abundance of giant arcs by factors of a few. The net boost to the lensing probability for individual clusters is mass dependent, and can lower the effective low mass limit of lensing clusters. This last effect can potentially increase the number of lensing clusters by an extra 50%. While the magnitude of these effects may be overestimated due to the well known overcooling problem in simulations, it is evident that baryonic cooling has a non-negligible impact on the expected abundance of giant arcs, and hence cosmological constraints from giant arc abundances may be subject to large systematic errors.

X-ray observations of galaxy clusters potentially provide powerful cosmological probes if systematics due to our incomplete knowledge of the intracluster medium (ICM) physics are understood and controlled. In this paper, we present mock Chandra analyses of cosmological cluster simulations and assess X-ray measurements of galaxy cluster properties using a model and procedure essentially identical to that used in real data analysis. We show that reconstruction of three-dimensional ICM density and temperature profiles is excellent for relaxed clusters, but still reasonably accurate for unrelaxed systems. The total ICM mass is measured quite accurately (<6%) in all clusters, while the hydrostatic estimate of the gravitationally bound mass is biased low by about 5%-20% through the virial region, primarily due to additional pressure support provided by subsonic bulk motions in the ICM, ubiquitous in our simulations even in relaxed systems. Gas fraction determinations are therefore biased high; the bias increases toward cluster outskirts and depends sensitively on its dynamical state, but we do not observe significant trends of the bias with cluster mass or redshift. We also find that different average ICM temperatures, such as the X-ray spectroscopic Tspec and gas-mass-weighted Tmg, are related to each other by a constant factor with a relatively small object-to-object scatter and no systematic trend with mass, redshift or the dynamical state of clusters. We briefly discuss direct applications of our results for different cluster-based cosmological tests.

This presentation is a Moriond version of our recent paper (Kravtsov, Vikhlinin & Nagai astro-ph/0603205) where we discussed X-ray proxies for the total cluster mass, including the spectral temperature (Tx), gas mass measured within r500 (Mg), and the new proxy, Yx, which is a simple product of Tx and Mg. We use mock Chandra images constructed for a sample of clusters simulated with high resolution in the concordance LambdaCDM cosmology. The simulated clusters exhibit tight correlations between the considered observables and total mass. The normalizations of the M500-Tx, Mg-Tx, and M500-Yx relations agree to better than =~ 10-15% with the current observational measurements of these relations. Our results show that Yx is the best mass proxy with a remarkably low scatter of only =~ 5-7% in M500 for a fixed Yx, at both low and high redshifts and regardless of whether clusters are relaxed or not. In addition, we show that redshift evolution of the Yx-M500 relation is close to the self-similar prediction, which makes Yx a very attractive mass indicator for measurements of the cluster mass function from X-ray selected samples.

I review recent progress in understanding and modeling galaxy clustering in cosmological simulations, with emphasis on models based on high-resolution dissipationless simulations. During the last decade, significant advances in our understanding of abundance and clustering of dark matter halos allowed construction of accurate, quantitative models of galaxy clustering both in linear and non-linear regimes. Results of several recent studies show that dissipationless simulations with a simple, non-parametric model for the relation between halo circular velocity and luminosity of the galaxy they host predict the shape, amplitude, and luminosity dependence of the two-point correlation function in excellent agreement with the observed galaxy clustering in the SDSS data at z~0 and in the DEEP2 samples at z~1 over the entire probed range of projected separations. In particular, the small-scale upturn of the correlation function from the power-law form in the SDSS and DEEP2 luminosity-selected samples is reproduced very well. At z~3-5, predictions also match the observed shape and amplitude of the angular two-point correlation function of Lyman-break galaxies (LBGs) on both large and small scales, including the theoretically predicted strong upturn at small scales. This suggests that, like galaxies in lower redshift samples, the LBGs are fair tracers of the overall halo population and that their luminosity is tightly correlated with the circular velocity (and hence mass) of their dark matter halos.

We investigate the dependence of dark matter halo clustering on halo formation time, density profile concentration, and subhalo occupation number, using high-resolution numerical simulations of a LCDM cosmology. We confirm results that halo clustering is a function of halo formation time, and that this trend depends on halo mass. For the first time, we show unequivocally that halo clustering is a function of halo concentration and show that the dependence of halo bias on concentration, mass, and redshift can be accurately parameterized in a simple way: b(c,M|z) = b(M|z) b(c|M/M*). The scaling between bias and concentration changes sign with the value of M/M*: high concentration (early forming) objects cluster more strongly for M <~ M* while low concentration (late forming) objects cluster more strongly for rare high-mass halos, M >~ M*. We show the first explicit demonstration that host dark halo clustering depends on the halo occupation number (of dark matter subhalos) and discuss implications for halo model calculations of dark matter power spectra and galaxy clustering statistics. The effect of these halo properties on clustering is strongest for early-forming dwarf-mass halos, which are significantly more clustered than typical halos of their mass. Our results suggest that isolated low-mass galaxies (e.g. low surface-brightness dwarfs) should have more slowly-rising rotation curves than their clustered counterparts, and may have consequences for the dearth of dwarf galaxies in voids. They also imply that self calibrating richness-selected cluster samples with their clustering properties might overestimate cluster masses and bias cosmological parameter estimation.

In this paper we examine how well galaxies and intra-cluster gas trace the gravitational potential of clusters. Utilizing mass profiles derived from gravitational lensing and X-ray observations, coupled with measured galaxy velocities, we solve for the velocity anisotropy parameter using the anisotropic Jeans equation. This is done for five clusters, three at low redshift: A2199, A496 and A576 and two at high redshifts: A2390 and MS1358. With X-ray temperature profiles obtained from Chandra and ASCA/ROSAT data, we estimate the ratio of energy in the galaxies compared to the X-ray gas. We find that none of these clusters is strictly in hydro-static equilibrium. We compare the properties of our sample with clusters that form in high-resolution cosmological N-body simulations that include baryonic physics. Simulations and data show considerable scatter both these profiles. We demonstrate the future feasibility and potential for directly comparing the orbital structure of clusters inferred from multi-wavelength observations with high resolution simulated clusters.

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.

We present comparison of X-ray proxies for the total cluster mass, including the spectral temperature (Tx), gas mass measured within r500 (Mg), and the new proxy, Yx, which is a simple product of Tx and Mg and is related to the total thermal energy of the ICM. We use mock Chandra images constructed for a sample of clusters simulated with the eulerian N-body+gasdynamics adaptive mesh refinement ART code in the concordance LCDM cosmology. The simulations achieve high spatial and mass resolution and include radiative cooling, star formation, and other processes accompanying galaxy formation. Our analysis shows that simulated clusters exhibit a high degree of regularity and tight correlations between the considered observables and total mass. The normalizations of the M-Tx, Mg-Tx, and M-Yx relations agree to better than 10-15% with the current observational measurements of these relations. Our results show that Yx is the best mass proxy with a remarkably low scatter of only ~5-7% in M500 for a fixed Yx, at both low and high redshifts and regardless of whether clusters are relaxed or not. In addition, we show that redshift evolution of the Yx-M500 relation is close to the self-similar prediction, which makes Yx a very attractive mass indicator for measurements of the cluster mass function from X-ray selected samples.

We employ high-resolution dissipationless simulations of the concordance LCDM cosmology to model the observed luminosity dependence and evolution of galaxy clustering through most of the age of the universe, from z~5 to z~0. We use a simple, non-parametric model which monotonically relates galaxy luminosities to the maximum circular velocity of dark matter halos (V_max) by preserving the observed galaxy luminosity function in order to match the halos in simulations with observed galaxies. The novel feature of the model is the use of the maximum circular velocity at the time of accretion, V_max,acc, for subhalos, the halos located within virial regions of larger halos. We argue that for subhalos in dissipationless simulations, V_max,acc reflects the luminosity and stellar mass of the associated galaxies better than the circular velocity at the epoch of observation, V_max,now. The simulations and our model L-V_max relation predict the shape, amplitude, and luminosity dependence of the two-point correlation function in excellent agreement with the observed galaxy clustering in the SDSS data at z~0 and in the DEEP2 samples at z~1 over the entire probed range of projected separations, 0.1<r_p/(Mpc/h)<10.0. In particular, the small-scale upturn of the correlation function from the power-law form in the SDSS and DEEP2 luminosity-selected samples is reproduced very well. At z~3-5, our predictions also match the observed shape and amplitude of the angular two-point correlation function of Lyman-break galaxies (LBGs) on both large and small scales, including the small-scale upturn.

We use a combination of high-resolution gasdynamics simulations of high-redshift dwarf galaxies and dissipationless simulations of a Milky Way sized halo to estimate the expected abundance and spatial distribution of the dwarf satellite galaxies that formed most of their stars around z~8 and evolved only little since then. Such galaxies can be considered as fossils of the reionization era, and studying their properties could provide a direct window into the early, pre-reionization stages of galaxy formation. We show that 5-15% of the objects existing at z~8 do indeed survive until the present in the MW like environment without significant evolution. This implies that it is plausible that the fossil dwarf galaxies do exist in the Local Group. Because such galaxies form their stellar systems early during the period of active merging and accretion, they should have spheroidal morphology regardless of their current distance from the host galaxy. We show that both the expected luminosity function and spatial distribution of dark matter halos which are likely to host fossil galaxies agree reasonably well with the observed distributions of the luminous (L_V>10^6 Lsun) Local Group fossil candidates near the host galaxy (d<200 kpc). However, the predicted abundance is substantially larger (by a factor of 2-3) for fainter galaxies (L_V<10^6 Lsun) at larger distances (d>300 kpc). We discuss several possible explanations for this discrepancy.

We present a comprehensive series of dissipationless N-body simulations to investigate the evolution of density distribution in equal-mass mergers between dark matter (DM) halos and multicomponent galaxies. The DM halo models are constructed with various asymptotic power-law indices ranging from steep cusps to core-like profiles and the structural properties of the galaxy models are motivated by the LCDM paradigm of structure formation. The adopted force resolution allows robust density profile estimates in the inner ~1% of the virial radii of the simulated systems. We demonstrate that the central slopes and overall shapes of the remnant density profiles are virtually identical to those of the initial systems suggesting that the remnants retain a remarkable memory of the density structure of their progenitors, despite the relaxation that accompanies merger activity. We also find that halo concentrations remain approximately constant through hierarchical merging involving identical systems and show that remnants contain significant fractions of their bound mass well beyond their formal virial radii. These conclusions hold for a wide variety of initial asymptotic density slopes, orbital energies, and encounter configurations, including sequences of consecutive merger events, simultaneous mergers of severals ystems, and mergers of halos with embedded cold baryonic components in the form of disks, spheroids, or both. As an immediate consequence, the net effect of gas cooling, which contracts and steepens the inner density profiles of DM halos, should be preserved through a period of dissipationless major merging. Our results imply that the characteristic universal shape of DM density profiles may be set early in the evolution of halos.

The Dark Energy Survey (DES) will use a new imaging camera on the Blanco 4-m telescope at CTIO to image 5000 square degrees of sky in the South Galactic Cap in four optical bands, and to carry out repeat imaging over a smaller area to identify and measure lightcurves of Type Ia supernovae. The main imaging area overlaps the planned Sunyaev-Zel'dovich survey of the South Pole Telescope. The idea behind DES is to use four distinct and largely independent methods to probe the properties of dark energy: baryon oscillations of the power spectrum, abundance and spatial distribution of clusters, weak gravitational lensing, and Type Ia supernovae. This white paper outlines, in broad terms, some of the theoretical issues associated with the first three of these probes (the issues for supernovae are mostly different in character), and with the general task of characterizing dark energy and distinguishing it from alternative explanations for cosmic acceleration. A companion white paper discusses the kind of numerical simulations and other theoretical tools that will be needed to address the these issues and to create mock catalogs that allow end-to-end tests of analysis procedures. Although we have been thinking about these problems in the specific context of DES, many of them are also relevant to other planned dark energy studies.

The ability to test the nature of dark mass-energy components in the universe through large-scale structure studies hinges on accurate predictions of sky survey expectations within a given world model. Numerical simulations predict key survey signatures with varying degrees of confidence, limited mainly by the complex astrophysics of galaxy formation. As surveys grow in size and scale, systematic uncertainties in theoretical modeling can become dominant. Dark energy studies will challenge the computational cosmology community to critically assess current techniques, develop new approaches to maximize accuracy, and establish new tools and practices to efficiently employ globally networked computing resources.

Galaxy-galaxy lensing is a powerful probe of the relation between galaxies and dark matter halos, but its theoretical interpretation requires a careful modeling of various contributions, such as the contribution from central and satellite galaxies. For this purpose, a phenomenological approach based on the halo model has been developed, allowing for fast exploration of the parameter space of models. In this paper, we investigate the ability of the halo model to extract information from the g-g weak lensing signal by comparing it to high-resolution dissipationless simulations that resolve subhalos. We find that the halo model reliably determines parameters such as the host halo mass of central galaxies, the fraction of galaxies that are satellites, and their radial distribution inside larger halos. If there is a significant scatter present in the central galaxy host halo mass distribution, then the mean and median mass of that distribution can differ significantly from one another, and the halo model mass determination lies between the two. This result suggests that when analyzing the data, galaxy subsamples with a narrow central galaxy halo mass distribution, such as those based on stellar mass, should be chosen for a simpler interpretation of the results.

(Abridged) We study predictions for galaxy cluster observables that can test the statistics of dark matter halo shapes expected in a flat LCDM universe. We present a simple analytical model for the prediction of cluster-scale X-ray observations, approximating clusters as isothermal systems in hydrostatic equilibrium, and dark matter haloes as ellipsoids with uniform axial ratios. We test the model against high-resolution, hydrodynamic cluster simulations to gauge its reliability. We find that this simple prescription does a good job of predicting the distribution of cluster X-ray ellipticities compared to the simulations as long as one focuses on cluster regions that are less sensitive to recent mergers. Based on this simple model, the distribution of cluster-size halo shapes expected in the concordance LCDM cosmology implies an X-ray ellipticity distribution with a mean of 0.32 +- 0.01 and a scatter of 0.14 +- 0.01 for the mass range (1-4)x10^{14} Msun/h. We find it important to include the mass dependence of halo shape to make comparisons to observational samples that contain many, very massive clusters. We analyse the systematics of four observational samples of cluster ellipticities and find that our results are statistically compatible with observations. In particular, we find remarkably good agreement between two recent ROSAT samples and LCDM predictions that DO NOT include gas cooling. We also test how well our analytical model can predict Sunyaev-Zel'dovich decrement maps and find that it is less successful although still useful; the model does not perform as well as a function of flux level in this case because of the changing triaxiality of dark matter haloes as a function of radial distance. Both this effect and the changing alignment of isodensity shells of dark matter haloes leave an imprint on cluster gas...

Using six high resolution dissipationless simulations with a varying box size in a flat LCDM universe, we study the mass and redshift dependence of dark matter halo shapes for M_vir = 9.0e11 - 2.0e14, over the redshift range z=0-3, and for two values of sigma_8=0.75 and 0.9. Remarkably, we find that the redshift, mass, and sigma_8 dependence of the mean smallest-to-largest axis ratio of halos is well described by the simple power-law relation <s> = (0.54 +- 0.02)(M_vir/M_*)^(-0.050 +- 0.003), where s is measured at 0.3 R_vir and the z and sigma_8 dependences are governed by the characteristic nonlinear mass, M_*=M_*(z,sigma_8). We find that the scatter about the mean s is well described by a Gaussian with sigma ~ 0.1, for all masses and redshifts. We compare our results to a variety of previous works on halo shapes and find that reported differences between studies are primarily explained by differences in their methodologies. We address the evolutionary aspects of individual halo shapes by following the shapes of the halos through ~100 snapshots in time. We determine the formation scalefactor a_c as defined by Wechsler et al. (2002) and find that it can be related to the halo shape at z = 0 and its evolution over time.

We explore the effect of substructure in dark matter halos on the power spectrum and bispectrum of matter fluctuations and weak lensing shear. By experimenting with substructure in a cosmological N = 512^3 simulation, we find that when a larger fraction of the host halo mass is in subhalos, the resulting power spectrum has less power at 1 < k < 100 h Mpc^{-1} and more power at k > 100 h Mpc^{-1}. We explain this effect using an analytic halo model including subhalos, which shows that the 1 < k < 100 h Mpc^{-1} regime depends sensitively on the radial distribution of subhalo centers while the interior structure of subhalos is important at k > 100 h Mpc^{-1}. The corresponding effect due to substructures on the weak lensing power spectrum is up to about 11% at angular scale l < 10^4. Predicting the nonlinear power spectrum to a few percent accuracy for future surveys would therefore require large cosmological simulations that also have exquisite numerical resolution to model accurately the survivals of dark matter subhalos in the tidal fields of their hosts.

We present gas and total mass profiles for 13 low-redshift, relaxed clusters spanning a temperature range 0.7-9 keV, derived from all available Chandra data of sufficient quality. In all clusters, gas temperature profiles are measured to large radii (Vikhlinin et al.) so that direct hydrostatic mass estimates are possible to nearly r_500 or beyond. The gas density was accurately traced to larger radii; its profile is not described well by a beta-model, showing continuous steepening with radius. The derived rho_tot profiles and their scaling with mass generally follow the Navarro-Frenk-White model with concentration expected for dark matter halos in LambdaCDM cosmology. In the inner region (r<0.1r_500), the gas density and temperature profiles exhibit significant scatter and trends with mass, but they become nearly self-similar at larger radii. Correspondingly, we find that the slope of the mass-temperature relation for these relaxed clusters is in good agreement with the simple self-similar behavior, M_500 ~ T^alpha, where alpha=(1.5-1.6)+-0.1, if the gas temperatures are measured excluding the central cool cores. The normalization of this M-T relation is significantly, by =~ 30%, higher than most previous X-ray determinations. We derive accurate gas mass fraction profiles, which show increase both with radius and cluster mass. The enclosed f_gas profiles within r_2500 =~ 0.4 r_500 have not yet reached any asymptotic value and are still far (by a factor of 1.5-2) from the Universal baryon fraction according to the CMB observations. The f_gas trends become weaker and its values closer to Universal at larger radii, in particular, in spherical shells r_2500<r<r_500.

We study the effects of dissipation on the baryon fractions in clusters using high-resolution cosmological simulations of nine clusters that resolve formation of cluster galaxies. The simulations of each cluster are performed with the shock-capturing eulerian adaptive mesh refinement N-body+gasdynamics ART code with and without radiative cooling. We show that dissipation and associated galaxy formation increase the total baryon fractions within radii as large as the virial radius. The effect is the strongest within cluster cores, where the simulations with cooling have baryon fractions larger than the universal value, while the fraction of baryons in adiabatic simulations are smaller than universal. At larger radii (r >~ r_500) the cumulative baryon fractions in simulations with cooling are close to, while those in the adiabatic runs remain below than, the universal value. The gas fractions in simulations with dissipation are reduced by ~20-40% at r<0.3r_vir and ~10% at larger radii compared to the adiabatic runs, because a fraction of gas is converted into stars. There is an indication that gas fractions within different radii increase with increasing cluster mass as f_gas ~ M_vir^0.2. We find that the total baryon fraction within the virial radius does not evolve with time in both adiabatic simulations and in simulations with cooling. Finally, to evaluate systematic uncertainties in the baryon fraction in cosmological simulations we present a comparison of gas fractions in our adiabatic simulations to re-simulations of the same objects with the entropy-conserving SPH code Gadget. The cumulative gas fraction profiles in the two sets of simulations on average agree to <~3% at r/r_vir>0.2, but differ systematically by up to 10% at small radii.

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.

We study the formation of globular clusters (GCs) in a Milky Way-size galaxy using a high-resolution cosmological simulation. The clusters in our model form in the dense cores of supergiant molecular clouds in the gaseous disks of high-redshift galaxies. The properties of clusters are estimated using a physically-motivated subgrid model of the isothermal cloud collapse. The first clusters in the simulation form at z ~ 12, while we conjecture that the best conditions for GC formation appear to be at z ~ 3-5. Most clusters form in the progenitor galaxies of the virial mass >10^9 Msun and the total mass of the cluster population is strongly correlated with the mass of its host galaxy with a fraction ~2x10^-4 of the galactic baryons being in the form of GCs. In addition, the mass of the GC population and the maximum cluster mass in a given region strongly correlate with the local average star formation rate. We find that the mass, size, and metallicity distributions of the cluster population identified in the simulation are remarkably similar to the corresponding distributions of the Milky Way globulars. We find no clear mass-metallicity or age-metallicity correlations for the old clusters. The zero-age cluster mass function can be approximated by a power-law, dN/dM ~ M^-alpha, with alpha ~ 2, in agreement with the mass function of young stellar clusters in starbursting galaxies. However, the shape of the zero-age mass function may be better described by the high-mass tail of a lognormal distribution which peaks at ~10^3 Msun. We discuss in detail the origin, universality, and dynamical evolution of the globular cluster mass function. Our results indicate that globular clusters with properties similar to those of observed clusters can form naturally within young dense gaseous disks at z >~ 3 in the LCDM cosmology.

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

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 study the radial distribution of subhalos and galaxies using high-resolution cosmological simulations of galaxy clusters formed in the concordance LCDM cosmology. In agreement with previous studies, we find that the radial distribution of subhalos is significantly less concentrated than that of the dark matter, when subhalos are selected using their present-day gravitationally bound mass. We show that the difference in the radial distribution is not a numerical artifact and is due to tidal stripping. The subhalos in the cluster core lose more than 70% of their initial mass since accretion, while the average tidal mass loss for halos near the virial radius is ~30%. This introduces a radial bias in the spatial distribution of subhalos when they are selected using their tidally truncated mass. We demonstrate that the radial bias disappears almost entirely if subhalos are selected using their mass or circular velocity at the accretion epoch. The comparisons of the results of dissipationless simulations to the observed distribution of galaxies in clusters are therefore sensitive to the selection criteria used to select subhalo samples. Using the simulations that include cooling and starformation, we show that the radial distribution of subhalos is in reasonable agreement with the observed radial distribution of galaxies in clusters for 0.1<R/R200<2.0, if subhalos are selected using the stellar mass of galaxies. The radial bias is minimized in this case because the stars are located in the centers of dark matter subhalos and are tightly bound. The stellar mass of an object is therefore approximately conserved as the dark matter is stripped from the outer regions. Nevertheless, the concentration of the radial distribution of galaxies is systematically lower than that of the dark matter.

We analyze the effect of dissipation on the shapes of dark matter (DM) halos using high-resolution cosmological gasdynamics simulations of clusters and galaxies in the LCDM cosmology. We find that halos formed in simulations with gas cooling are significantly more spherical than corresponding halos formed in adiabatic simulations. Gas cooling results in an average increase of the principle axis ratios of halos by ~ 0.2-0.4 in the inner regions. The systematic difference decreases slowly with radius but persists almost to the virial radius. We argue that the differences in simulations with and without cooling arise both during periods of quiescent evolution, when gas cools and condenses toward the center, and during major mergers. We perform a series of high-resolution N-body simulations to study the shapes of remnants in major mergers of DM halos and halos with embedded stellar disks. In the DM halo-only mergers, the shape of the remnants depends only on the orbital angular momentum of the encounter and not on the internal structure of the halos. However, significant shape changes in the DM distribution may result if stellar disks are included. In this case the shape of the DM halos is correlated with the morphology of the stellar remnants.

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 use high-resolution dissipationless simulations of the concordance flat LCDM model to make predictions for the galaxy--mass (GM) correlations and compare them to the recent SDSS weak lensing measurements.We use a simple observationally motivated scheme to assign luminosities and colors to the halos.This allows us to closely match the selection criteria used to define observational samples.The simulations reproduce the observed GM correlation function and its observed dependencies on luminosity and color.The luminosity dependence of the correlation function is primarily determined by the changing relative contribution of central and satellite galaxies at different luminosities. The color dependence of the GM correlations reflects the difference in the typical environments of blue and red galaxies. We also find agreement between the predicted and observed cross-bias, b_x=b/r,at all probed scales.The GM correlation coefficient, r, is close to unity on scales >1/h Mpc.The cross bias is thus expected to measure the actual bias of galaxy clustering on these scales.The aperture mass-to-light ratio is independent of galaxy color.The aperture mass scales approximately linearly with luminosity at L_r>10^{10}h^{-2} Lsun, while at lower luminosities the scaling is shallower: L_r^{0.5}. We show that most of the luminous galaxies (M_r<-21) are near the centers of their halos and their GM correlation function at r<100/h kpc can therefore be interpreted as the average dark matter density profile of these galaxies. We find that for galaxies in a given narrow luminosity range, there is a broad and possibly non-gaussian distribution of halo virial masses. Therefore, the average relation between mass and luminosity derived from the weak lensing analyses should be interpreted with caution.

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 the effects of radiative cooling, star formation and stellar feedback on the properties and evolution of galaxy clusters using high-resolution Adaptive Mesh Refinement N-body+gasdynamics simulations of clusters forming in the LCDM universe. Cooling leads to the condensation of gas in the inner regions of clusters, which in turn leads to steepening of the dark matter profile. The cooling gas is replaced by the higher-entropy gas from the outer regions, which raises the entropy and temperature of gas in the cluster core. The magnitude of these effects is likely overestimated in the current simulations because they suffer from the overcooling problem: a much larger fraction of baryons is in the form of cold gas and stars than is observed. We find that the thermal stellar feedback alone does not remedy this problem. Additional ad-hoc preheating can lower the amount of cold gas but a simple uniform preheating results in incorrect star formation history, as it delays the bulk of star formation until z<1. Our analysis shows that the overcooling in a cluster as a whole is really the overcooling in the central galaxy and its progenitors at high redshifts. This indicates that an additional heating mechanism that can continuously heat the gas in the cluster core is required to reproduce the observed cluster properties. Energy injection by the Active Galactic Nuclei, which may provide such heating, may thus be an important missing ingredient in the current theoretical models of cluster formation.

Multiply-imaged quasar lenses can be used to constrain the substructure mass fraction in galaxy-sized dark matter halos via anomalous flux ratios in lensed images. The flux ratios, however, can be affected by both the substructure in the lens halo and by isolated small-mass halos along the entire line-of-sight to the lensed source. While lensing by dark matter clumps near the lens galaxy is more efficient than elsewhere, the cumulative effect of all objects along the line-of-sight could be significant. Here we estimate the potential contribution of isolated clumps to the substructure lensing signal using a simple model motivated by cosmological simulations. We find that the contribution of isolated clumps to the total lensing optical depth ranges from a few to tens percent, depending on assumptions and the particular configuration of the lens. Therefore, although the contribution of isolated clumps to the lensing signal is not dominant, it should not be neglected in detailed analyses of substructure lensing. For the currently favored Lambda-CDM model, the total calculated optical depth for lensing is high, tau ~ 0.2-20 and could, therefore, naturally explain the high frequency of anomalous flux ratios in observed lenses. The prediction, however, is highly sensitive to the spatial distribution of substructure halos in the innermost regions of the lens halo, which is still very uncertain. Therefore, constraints on the properties of the substructure population or accurate cosmological constraints, such as the mass of the warm dark matter particle, are difficult if not impossible to derive at this point.

One of the most puzzling properties of observed galaxies is the universality of the empirical correlation between the star formation rate and average gas surface density on kiloparsec scales (the Schmidt law). In this study I present results of self-consistent cosmological simulations of high-redshift galaxy formation that reproduce the Schmidt law naturally, without assuming it, and provide some clues to this puzzle. The simulations incorporate the main physical processes critical to various aspects of galaxy formation and have a dynamic range high enough to identify individual star forming regions. The results indicate that the global Schmidt law is a manifestation of the overall density distribution of the interstellar medium (ISM). In particular, the density probability distribution function (PDF) in the simulated disks is similar to that observed in recent state-of-the-art modeling of the turbulent ISM and has a well-defined generic shape. The shape of the PDF in a given region of the disk depends on the local average surface density Sigma_g. The dependence is such that the fraction of gas mass in the high-density tail of the distribution scales as Sigma_g^{n-1} with n~1.4, which gives rise to the Schmidt-like correlation. The high-density tail of the PDF is remarkably insensitive to the inclusion of feedback and details of the cooling and heating processes. This indicates that the global star formation rate is determined by the supersonic turbulence driven by gravitational instabilities on large scales, rather than stellar feedback or thermal instabilities on small scales.

Recently, high-resolution Chandra observations revealed the existence of very sharp features in the X-ray surface brightness and temperature maps of several clusters (Vikhlinin et. al., 2001). These features, called ``cold fronts'', are characterized by an increase in surface brightness by a factor >2 over 10-50 kpc, accompanied by a drop in temperature of a similar magnitude. The existence of such sharp gradients can be used to put interesting constraints on the physics of the intracluster medium (ICM), if their mechanism and longevity are well understood. Here, we present results of a search for cold fronts in high-resolution simulations of galaxy clusters in cold dark matter (CDM) models. We show that sharp gradients with properties similar to those of observed cold fronts naturally arise in cluster mergers when the shocks heat gas surrounding the merging sub-cluster, while its dense core remains relatively cold. The compression induced by supersonic motions and shock heating during the merger enhance the amplitude of gas density and temperature gradients across the front. Our results indicate that cold fronts are non-equilibrium transient phenomena and can be observed for a period of less than a billion years. We show that the velocity and density fields of gas surrounding the cold front can be very irregular which would complicate analyses aiming to put constraints on the physical conditions of the intracluster medium in the vicinity of the front.

Galaxy clusters are potentially powerful probes of the large-scale velocity field in the Universe because their peculiar velocity can be estimated directly via the kinematic Sunyaev-Zeldovich effect (kSZ). Using high-resolution cosmological simulations of an evolving cluster of galaxies, we evaluate how well the average velocity obtained via a kSZ measurement reflects the actual cluster peculiar velocity. We find that the internal velocities in the intracluster gas are comparable to the overall cluster peculiar velocity, 20 to 30% of the sound speed even when a cluster is relatively relaxed. Nevertheless, the velocity averaged over the kSZ map inside a circular aperture matched to the cluster virial region provides an unbiased estimate of a cluster's radial peculiar velocity with a dispersion of 50 to 100 km/s, depending on the line of sight and dynamical state of the cluster. This dispersion puts a lower limit on the accuracy with which cluster peculiar velocity can be measured. Although the dispersion of the average is relatively small, the velocity distribution is broad; regions of low signal must be treated with care to avoid bias. We discuss the extent to which systematic errors might be modelled, and the resulting limitations on using galaxy clusters as cosmological velocity tracers.

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.

We present a general statistical framework for describing the effect of sample variance in the number counts of virialized objects and examine its effect on cosmological parameter estimation. Specifically, we consider effects of sample variance on the power spectrum normalization and properties of dark energy extracted from current and future local and high-redshift samples of clusters. We show that for future surveys that probe ever lower cluster masses and temperatures, sample variance is generally comparable to or greater than shot noise and thus cannot be neglected in deriving precision cosmological constraints. For example, sample variance is usually more important than shot variance in constraints on the equation of state of the dark energy from z < 1 clusters. Although we found that effects of sample variance on the sigma_8-Omega_m constraints from the current flux and temperature limited X-ray surveys are not significant, they may be important for future studies utilizing the shape of the temperature function to break the sigma_8-Omega_m degeneracy. We also present numerical tests clarifying the definition of cluster mass employed in cosmological modelling and accurate fitting formula for the conversion between different definitions of halo mass (e.g., virial vs. fixed overdensity).

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

(abridged) We study the relation between the density profiles of dark matter halos and their mass assembly histories, using a statistical sample of halos in a high-resolution N-body simulation of the LCDM cosmology. For each halo at z=0, we identify its merger-history tree, and determine concentration parameters c_vir for all progenitors, thus providing a structural merger tree for each halo. We fit the mass accretion histories by a universal function with one parameter, the formation epoch a_c, defined when the log mass accretion rate dlogM/dloga falls below a critical value S. We find that late forming galaxies tend to be less concentrated, such that c_vir ``observed'' at any epoch a_o is strongly correlated with a_c via c_vir=c_1*a_o/a_c. Scatter about this relation is mostly due to measurement errors in c_v and a_c, implying that the actual spread in c_vir for halos of a given mass can be mostly attributed to scatter in a_c. We demonstrate that this relation can also be used to predict the mass and redshift dependence of c_v, and the scatter about the median c_vir(M,z), using accretion histories derived from the Extended Press-Schechter (EPS) formalism, after adjusting for a constant offset between the formation times as predicted by EPS and as measured in the simulations;this new ingredient can thus be easily incorporated into semi-analytic models of galaxy formation. The correlation found between halo concentration and mass accretion rate suggests a physical interpretation: for high mass infall rates the central density is related to the background density; when the mass infall rate slows, the central density stays approximately constant and the halo concentration just grows as R_vir. The tight correlation demonstrated here provides an essential new ingredient for galaxy formation modeling.

We study the relation between the density profiles of dark matter halos and their mass assembly histories, using a statistical sample of halos in a high-resolution N-body simulation of the LCDM cosmology. For each halo at z=0, we identify its merger-history tree, and determine concentration parameters c_vir for all progenitors, thus providing a structural merger tree for each halo. We fit the mass accretion histories by a universal function with one parameter, the formation epoch a_c, defined when the log mass accretion rate dlog M/doga falls below a critical value. We find that late forming galaxies tend to be less concentrated, such that c_vir ``observed'' at any epoch a_obs is strongly correlated with a_c via c_vir=c_1 a_obs/a_c. Scatter about this relation is mostly due to measurement errors in c_vir and a_c, implying that the actual spread in c_vir for halos of a given mass can be mostly attributed to scatter in a_c. Because of the direct connection between halo concentration and velocity rotation curves, and because of probable connections between halo mass assembly history and star formation history, the tight correlation between these properties provides an essential new ingredient for galaxy formation modeling.

We compare the specific angular momentum profiles of virialized dark halos in cold dark matter (CDM) and warm dark matter (WDM) models using high-resolution dissipationless simulations. The simulations were initialized using the same set of modes, except on small scales, where the power was suppressed in WDM below the filtering length. Remarkably, WDM as well as CDM halos are well-described by the two-parameter angular momentum profile of Bullock et al. (2001), even though the halo masses are below the filtering scale of the WDM. Although the best-fit shape parameters change quantitatively for individual halos in the two simulations, we find no systematic variation in profile shapes as a function of the dark matter type. The scatter in shape parameters is significantly smaller for the WDM halos, suggesting that substructure and/or merging history plays a role producing scatter about the mean angular momentum distribution, but that the average angular momentum profiles of halos originate from larger-scale phenomena or a mechanism associated with the virialization process. The known mismatch between the angular momentum distributions of dark halos and disk galaxies is therefore present in WDM as well as CDM models. Our WDM halos tend to have a less coherent (more misaligned) angular momentum structure and smaller spin parameters than do their CDM counterparts, although we caution that this result is based on a small number of halos.

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

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

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 develop an explicit model for the formation of the stellar halo from tidally disrupted, accreted dwarf satellites in the cold dark matter (CDM) framework, focusing on predictions testable with the Sloan Digital Sky Survey (SDSS) and other wide-field surveys. Subhalo accretion and orbital evolution are calculated using a semi-analytic approach within the Press-Schechter formalism. Motivated by our previous work, we assume that low-mass subhalos (v < 30 km/s) can form significant populations of stars only if they accreted a substantial fraction of their mass before the epoch of reionization. With this assumption, the model reproduces the observed velocity function of galactic satellites in the Local Group, solving the ``dwarf satellite problem'' without modifying the popular LCDM cosmology. The disrupted satellites yield a stellar distribution with a total mass and radial density profile consistent with those observed for the Milky Way stellar halo. Most significantly, the model predicts the presence of many large-scale, coherent substructures in the outer halo. These substructures are remnants of individual, tidally disrupted dwarf satellite galaxies. Substructure is more pronounced at large galactocentric radii because of the smaller number density of tidal streams and the longer orbital times. This model provides a natural explanation for the coherent structures in the outer stellar halo found in the SDSS commissioning data, and it predicts that many more such structures should be found as the survey covers more of the sky. The detection (or non-detection) and characterization of such structures could eventually test variants of the CDM scenario, especially those that aim to solve the dwarf satellite problem by enhancing satellite disruption.

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 present estimates of the energy input from supernovae (SNe) into the intergalactic medium using (i) recent measurements of Si and Fe abundances in the intracluster medium (ICM) and (ii) self-consistent gasdynamical galaxy formation simulations that include processes of cooling, star formation, SNe feedback, and a multi-phase model of the interstellar medium. We estimate the energy input from observed abundances using two different assumptions: (i) spatial uniformity of metal abundances in the ICM and (ii) radial abundance gradients. We show that these two cases lead to energy input estimates which are different by an order of magnitude, highlighting a need for observational data on large-scale abundance gradients in clusters. Analysis of galaxy formation results and estimates from observed Fe and Si abundances indicates that the SNe energy input can be important for heating of the entire ICM (providing energy of ~1 keV per particle) only if the ICM abundances are uniform and the efficiency of gas heating by SN explosions is close to 100% (implying that all of the initial kinetic energy of the explosion goes into heating of the ICM). We conclude that unless these most favorable conditions are met, SNe alone are unlikely to provide sufficient energy input to heat all of the cluster ICM and may need to be supplemented or even substituted by some other heating process(es). (Abridged)

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.

One of the main challenges facing standard hierarchical structure formation models is that the predicted abundance of galactic subhalos with circular velocities of 10-30 km/s is an order of magnitude higher than the number of satellites actually observed within the Local Group. Using a simple model for the formation and evolution of dark halos, based on the extended Press-Schechter formalism and tested against N-body results, we show that the theoretical predictions can be reconciled with observations if gas accretion in low-mass halos is suppressed after the epoch of reionization. In this picture, the observed dwarf satellites correspond to the small fraction of halos that accreted substantial amounts of gas before reionization. The photoionization mechanism naturally explains why the discrepancy between predicted halos and observed satellites sets in at about 30 km/s, and for reasonable choices of the reionization redshift (z_re = 5-12) the model can reproduce both the amplitude and shape of the observed velocity function of galactic satellites. If this explanation is correct, then typical bright galaxy halos contain many low-mass dark matter subhalos. These might be detectable through their gravitational lensing effects, through their influence on stellar disks, or as dwarf satellites with very high mass-to-light ratios. This model also predicts a diffuse stellar component produced by large numbers of tidally disrupted dwarfs, perhaps sufficient to account for most of the Milky Way's stellar halo.

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

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.

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)

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

We have used integral two-point spatial correlation function and its second derivative to analyze the distribution of quasars in three very deep surveys published in the literature. Statistically significant ($\sim 2-3\sigma $) correlations were found at scales of $\sim 50-100h^{-1}$ Mpc in all of the analyzed surveys. We have used the friends-of-friends cluster analysis to show that these correlations can be explained by the presence of relatively small quasar clusters (3-6 objects) which may possibly belong to larger structures such as Large Quasar Groups found in the bigger surveys. The sizes of these clusters along the redshift direction and distances between them are similar to those for structures found recently in studies of CIV absorption systems. These results present further evidence for the existence of large-scale structures at redshifts $z\sim 1-2$.

Recently, it was suggested that large concentrations or groups of quasars may trace sites of enhanced matter density at medium and high redshifts analogous to how galaxy clusters trace them in nearby space. We have checked existing quasar data for the presence of such groups. Large Quasar Groups (LQGs) were identified using a well-known cluster analysis technique. Our final list consists of 12 groups, including one reported previously. Most of the quasars in these groups come from deep homogeneous surveys. We have found that: (i) the probability that the detected groups are random is rather small (generally a few per cent); (ii) their sizes range from $\sim 70$ to $\sim 160h^{-1} Mpc$, which is comparable to sizes of nearby rich superclusters; (iii) the detected groups all have redshifts $0.5<z\leq 2$; (iv) the abundance of the LQGs is comparable with the abundance of large superclusters at $z\sim 0$, which is consistent with the idea that quasar groups and superclusters can be evolutionarily related. We argue that quasar groups could be a common feature of the spatial distribution of medium redshift quasars, and that the quasars in groups may belong to concentrations of young galaxy clusters and groups (distant superclusters) and hence be biased tracers of large-scale structure of matter distribution in the early Universe. Theoretical implications as well as other observations needed to test this point are discussed.

We report the investigation of spatial distribution of quasars using two different methods: the statistical study by means of 2-point correlation function and direct search for structures such as previously reported Large Quasar Groups. Having analyzed the combined sample of eight homogeneous quasar surveys (totalling near 2200 objects), we detected strong clustering at scales $r\leq 20h^{-1}$ Mpc ($H_0=100h$ and $q_0=0.5$) which seems, however, to be redshift dependent. Using reasonable selection criteria we have also sought for Large Quasar Groups (LQG) similar to those which were earlier found observationally. We discuss the properties of the detected groups and their implications for cosmology. We also discuss the suitability of the present data for getting the reliable statistical inference on quasar clustering.