Search

We discuss the mechanism of cluster formation in a numerical simulation of a molecular cloud (MC) undergoing global hierarchical collapse (GHC). The global nature of the collapse implies that the SFR increases over time. The hierarchical nature of the collapse consists of small-scale collapses within larger-scale ones. The large-scale collapses culminate a few Myr later than the small-scale ones and consist of filamentary flows that accrete onto massive central clumps. The small-scale collapses form clumps that are embedded in the filaments and falling onto the large-scale collapse centers. The stars formed in the early, small-scale collapses share the infall motion of their parent clumps. Thus, the filaments feed both gaseous and stellar material to the massive central clump. This leads to the presence of a few older stars in a region where new protostars are forming, and also to a self-similar structure, in which each unit is composed of smaller-scale sub-units that approach each other and may merge. Because the older stars formed in the filaments share the infall motion of the gas onto the central clump, they tend to have larger velocities and to be distributed over larger areas than the younger stars formed in the central clump. Finally, interpreting the IMF at face-value as a probability distribution implies that massive stars only form once the {\it local} SFR is large enough to sample the IMF up to high masses. In combination with the increase of the SFR, this implies that massive stars tend to appear late in the evolution of the MC, and only in the central massive clumps. We discuss the correspondence of these features with observed properties of young stellar clusters, finding very good qualitative agreement, thus providing support to the scenario of global, hierarchical collapse of MCs, while explaining the origin of the observed cluster structure., Comment: 17 pages, 12 figures. Submitted to MNRAS. Comments welcome

We introduce a new set of eight Milky Way-sized cosmological simulations performed using the AMR code ART + Hydrodynamics in a LCDM cosmology. The set of zoom-in simulations covers present-day virial masses in the 0.83-1.56 x 10^12 msun range and is carried out with our simple but effective deterministic star formation (SF) and ``explosive' stellar feedback prescriptions. The work is focused on showing the goodness of the simulated set of ``field' Milky Way-sized galaxies. Our results are as follows. (a) The circular velocity curves of our simulated galaxies are nearly flat. (b) Runs ending with a significant disk component, for their stellar masses, have V_max, radius, SF rate, gas fraction, and specific angular momentum values consistent with observations of late-type galaxies. (C) The two most spheroid-dominated galaxies formed in halos with late active merger histories, but other run that ends also as spheroid-dominated, never had major mergers. (d) Our simulations are consistent with the empirical stellar-to-halo mass correlation, and those that end as disk-dominated, evolve mostly along the low-mass branch of this correlation. (e) Moreover, since the last 6.5-10 Gyr, the baryonic/stellar and halo mass growth histories are proportional. (f) Within Rvir ~ 25-50% of the baryons are missed. (g) The z ~ 0 gas velocity dispersion profiles, sigma_z(r), are nearly flat and can be mostly explained by the kinetic energy injected by stars. (h) The average values of sigma_z increase at higher redshifts, following roughly the shape of the SF history., Comment: 17 pages, 14 figures, ApJ accepted

We discuss the mechanism of cluster formation in hierarchically collapsing molecular clouds. Recent evidence, both observational and numerical, suggests that molecular clouds (MCs) may be undergoing global, hierarchical gravitational collapse. The "hierarchical" regime consists of small-scale collapses within larger-scale ones. The latter implies that the star formation rate increases systematically during the early stages of evolution, and occurs via filamentary flows onto "hubs" of higher density, mass, and velocity dispersion, and culminates a few Myr after than the small-scale collapses have started to form stars. In turn, the small-scale collapses occur in clumps embedded in the filaments, and are themselves falling into the larger potential well of the still-ongoing large-scale collapse. The stars formed in the early, small-scale collapses share the infall motion of their parent clumps towards the larger potential trough, so that the filaments feed both gaseous and stellar material to the hubs. This leads to the presence of older stars in a region where new protostars are still forming, to a scale-free or fractal structure of the clusters, in which each unit is composed of smaller-scale ones, and to the eventual merging of the subunits, explaining the observed structural features of open clusters., Comment: 6 pages, 7 figures, to appear in the proceedings of IAU Symposium 316: Formation, Evolution, and Survival of Massive Star Clusters

By means of N-body+Hydrodynamic zoom-in simulations we study the evolution of the inner dark matter and stellar mass distributions of central dwarf galaxies formed in halos of virial masses Mv=2-3x10^10 Msun at z=0, both in a WDM and CDM cosmology. The half-mode mass in the WDM power spectrum of our simulations is Mf= 2x 10^10 Msun. In the dark matter only simulations halo density profiles are well described by the NFW parametric fit in both cosmologies, though the WDM halos have concentrations lower by factors 1.5--2.0 than their CDM counterparts. In the hydrodynamic simulations, the effects of baryons significantly flatten the inner density, velocity dispersion, and pseudo phase-space density profiles of the WDM halos but not of the CDM ones. The density slope measured at ~ 0.02xRv, alpha, becomes shallow in periods of 2 to 5 Gyr in the WDM runs. We explore whether this flattening process correlates with the global SF, Ms/Mv ratio, gas outflow, and internal specific angular momentum histories. We do not find any clear trends but when alpha is shallower than -0.5, Ms/Mv is always between 0.25 and 1%. We conclude that the main reason of the formation of the shallow core is the presence of strong gas mass fluctuations inside the inner halo, which are a consequence of the feedback driven by a very bursty and sustained SF history in shallow gravitational potentials. Our WDM halos, which assemble late and are less concentrated than the CDM ones, obey these conditions. There are also (rare) CDM systems with extended mass assembly histories that obey these conditions and form indeed shallow cores. The dynamical heating and expansion processes, behind the DM core flattening, apply also to the stars in a such a way that the stellar age and metallicity gradients of the dwarfs are softened, their stellar half-mass radii strongly grow with time, and their central surface densities decrease., Comment: 12 pages, 8 figures. Published in ApJ

We introduce a new set of simulations of Milky Way-sized galaxies using the AMR code ART + hydrodynamics in a $\Lambda$CDM cosmogony. The simulation series is named GARROTXA and follow the formation of a halo/galaxy from z~$=$~60 to z~$=$~0. The final virial mass of the system is $\sim$7.4$\times$10$^{11}$M$_{\odot}$. Our results are as follows: (a) contrary to many previous studies, the circular velocity curve shows no central peak and overall agrees with recent MW observations. (b) Other quantities, such as M$\_{*}$(6$\times$10$^{10}$M$_\odot$) and R$_d$ (2.56 kpc), fall well inside the observational MW range. (c) We measure the disk-to-total ratio kinematically and find that D/T=0.42. (d) The cold gas fraction and star formation rate (SFR) at z=0, on the other hand, fall short from the values estimated for the Milky Way. As a first scientific exploitation of the simulation series, we study the spatial distribution of the hot X-ray luminous gas. We have found that most of this X-ray emitting gas is in a halo-like distribution accounting for an important fraction but not all of the missing baryons. An important amount of hot gas is also present in filaments. In all our models there is not a massive disk-like hot gas distribution dominating the column density. Our analysis of hot gas mock observations reveals that the homogeneity assumption leads to an overestimation of the total mass by factors 3 to 5 or to an underestimation by factors $0.7-0.1$, depending on the used observational method. Finally, we confirm a clear correlation between the total hot gas mass and the dark matter halo mass of galactic systems., Comment: Accepted for publication, ApJ, April 2016

We present zoom-in N-body + Hydrodynamic simulations of dwarf central galaxies formed in Warm Dark Matter (WDM) halos with masses at present-day of $2-4\times 10^{10}$ \msun. Two different cases are considered, the first one when halo masses are close to the corresponding half-mode filtering scale \Mhm\ (\mwdm =1.2 keV), and the second when they are 20 to 30 times the corresponding \Mhm\ (\mwdm = 3.0 keV). The WDM simulations are compared with the respective Cold Dark Matter (CDM) simulations. The dwarfs formed in halos of masses (20-30)\Mhm have roughly similar properties and evolution than their CDM counterparts; on the contrary, those formed in halos of masses around \Mhm, are systematically different from their CDM counterparts. As compared to the CDM dwarfs, they assemble the dark and stellar masses later, having mass-weighted stellar ages 1.4--4.8 Gyr younger; their circular velocity profiles are shallower, with maximal velocities 20--60% lower; their stellar distributions are much less centrally concentrated and with larger effective radii, by factors 1.3--3. The WDM dwarfs at the filtering scale (\mwdm =1.2 keV) have disk-like structures, and end in most cases with higher gas fractions and lower stellar-to-total mass ratios than their CDM counterparts. The late halo assembly, low halo concentrations, and the absence of satellites of the former with respect to the latter, are at the basis of the differences., Comment: Accepted for publication in ApJ. 15 pages, 9 figures, and 1 table

Motivated by recent numerical simulations of molecular cloud (MC)evolution, in which the clouds engage in global gravitational contraction, and local collapse events culminate significantly earlier than the global collapse, we investigate the growth of density perturbations embedded in a collapsing background, to which we refer as an Inverse Hubble Flow (IHF). We use the standard procedure for the growth of perturbations in a universe that first expands (the usual Hubble Flow) and then recollapses (the IHF). We find that linear density perturbations immersed in an IHF grow faster than perturbations evolving in a static background (the standard Jeans analysis). A fundamental distinction between the two regimes is that, in the Jeans case, the time $\tau_\mathrm{nl}$ for a density fluctuation to become nonlinear increases without limit as its initial value approaches zero, while in the IHF case $\tau_\mathrm{nl} \le \tau_\mathrm{ff}$ always, where $\tau_\mathrm{ff}$ is the free-fall time of the background density. We suggest that this effect, although moderate, implies that small-scale density fluctuations embedded in globally-collapsing clouds must collapse earlier than their parent cloud, regardless of whether the initial amplitude of the fluctuations is moderate or strongly nonlinear, thus allowing the classical mechanism of Hoyle fragmentation to operate in multi-Jeans-mass MCs. More fundamentally, our results show that, contrary to the standard paradigm that fluctuations of all scales grow at the same rate in the linear regime, the hierarchical nesting of the fluctuations of different scales does affect their growth even in the linear stage., Comment: 7 pages, 3 figures; published in MNRAS

Supermassive black hole (BH) binaries would comprise the strongest sources of gravitational waves (GW) once they reach <<1 pc separations, for both pulsar timing arrays (PTAs) and space based (SB) detectors. While BH binaries coalescences constitute a natural outcome of the cosmological standard model and galaxy mergers, their dynamical evolution is still poorly understood and therefore their abundances at different stages. We use a dynamical model for the decay of BH binaries coupled with a cosmological simulation and semi-empirical approaches to the occupation of haloes by galaxies and BHs, in order to follow the evolution of the properties distribution of galaxies hosting BH binaries candidates to decay due to GWs emission. Our models allow us to relax simplifying hypothesis about the binaries occupation in galaxies and their mass, as well as redshift evolution. Following previously proposed electromagnetic (EM) signatures of binaries in the subpc regime, that include spectral features and variability, we model possible distributions of such signatures and also set upper limits to their lifespan. We found a bimodal distribution of hosts properties, corresponding to BH binaries suitable to be detected by PTA and the ones detectable only from space missions, as eLISA. Although it has been discussed that the peak of eLISA sources may happen at high z, we show that there must be a population of such sources in the nearby Universe that might show detectable EM signatures, representing an important laboratory for multimessenger astrophysics. We found a weak dependence of galaxy host properties on the binaries occupation, that can be traced back to the BH origin. The combination of the host correlations reported here with the expected EM signal, may be helpful to verify the presence of nearby GW candidates, and to distinguish them from 'regular' intrinsic AGN variability., Comment: 16 pages, 17 figures, submitted, comments are welcome

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

We present a numerical study of the evolution of molecular clouds, from their formation by converging flows in the warm ISM, to their destruction by the ionizing feedback of the massive stars they form. We improve with respect to our previous simulations by including a different stellar-particle formation algorithm, which allows them to have masses corresponding to single stars rather than to small clusters, and with a mass distribution following a near-Salpeter stellar IMF. We also employ a simplified radiative-transfer algorithm that allows the stellar particles to feed back on the medium at a rate that depends on their mass and the local density. Our results are as follows: a) Contrary to the results from our previous study, where all stellar particles injected energy at a rate corresponding to a star of ~ 10 Msun, the dense gas is now completely evacuated from 10-pc regions around the stars within 10-20 Myr, suggesting that this feat is accomplished essentially by the most massive stars. b) At the scale of the whole numerical simulations, the dense gas mass is reduced by up to an order of magnitude, although star formation (SF) never shuts off completely, indicating that the feedback terminates SF locally, but new SF events continue to occur elesewhere in the clouds. c) The SF efficiency (SFE) is maintained globally at the ~ 10% level, although locally, the cloud with largest degree of focusing of its accretion flow reaches SFE ~ 30%. d) The virial parameter of the clouds approaches unity before the stellar feedback begins to dominate the dynamics, becoming much larger once feedback dominates, suggesting that clouds become unbound as a consequence of the stellar feedback. e) The erosion of the filaments that feed the star-forming clumps produces chains of isolated dense blobs reminiscent of those observed in the vicinity of the dark globule B68., Comment: Accepted for publication in MNRAS. This version include two more figures and four more pages with a more extended discusion

We study the effect that the dark matter background (DMB) has on the gravitational energy content and, in general, on the star formation efficiency of a molecular cloud (MC). We first analyze the effect that a dark matter halo, described by the Navarro et al. (1996) density profile, has on the energy budget of a spherical, homogeneous, cloud located at different distances from the halo center. We found that MCs located in the innermost regions of a massive galaxy can feel a contraction force greater than their self-gravity due to the incorporation of the potential of the galaxy's dark matter halo. We also calculated analytically the gravitational perturbation that a MC produces over a uniform DMB (uniform at the scales of a MC) and how this perturbation will affect the evolution of the MC itself. The study shows that the star formation in a MC will be considerably enhanced if the cloud is located in a dense and low velocity dark matter environment. We confirm our results by measuring the star formation efficiency in numerical simulations of the formation and evolution of MCs within different DMBs. Our study indicates that there are situations where the dark matter's gravitational contribution to the evolution of the molecular clouds should not be neglected., Comment: Same as previous version, but accepted by the Astrophysical Journal; 12 pages, 12 figures

It has been recently shown that molecular clouds do not exhibit a unique shape for the column density probability distribution function (Npdf). Instead, clouds without star formation seem to possess a lognormal distribution, while clouds with active star formation develope a power-law tail at high column densities. The lognormal behavior of the Npdf has been interpreted in terms of turbulent motions dominating the dynamics of the clouds, while the power-law behavior occurs when the cloud is dominated by gravity. In the present contribution we use thermally bi-stable numerical simulations of cloud formation and evolution to show that, indeed, these two regimes can be understood in terms of the formation and evolution of molecular clouds: a very narrow lognormal regime appears when the cloud is being assembled. However, as the global gravitational contraction occurs, the initial density fluctuations are enhanced, resulting, first, in a wider lognormal Npdf, and later, in a power-law Npdf. We thus suggest that the observed Npdf of molecular clouds are a manifestation of their global gravitationally contracting state. We also show that, contrary to recent suggestions, the exact value of the power-law slope is not unique, as it depends on the projection in which the cloud is being observed., Comment: 9 pages, 4 figures. MNRAS, accepted

We present an idealized, semi-empirical model for the evolution of gravitationally contracting molecular clouds (MCs) and their star formation rate (SFR) and efficiency (SFE). The model assumes that the instantaneous SFR is given by the mass above a certain density threshold divided by its free-fall time. The instantaneous number of massive stars is computed assuming a Kroupa IMF. These stars feed back on the cloud through ionizing radiation, eroding it. The main controlling parameter of the evolution turns out to be the maximum cloud mass, $\Mmax$. This allows us to compare various properties of the model clouds against their observational counterparts. A giant molecular cloud (GMC) model ($\Mmax \sim 10^5 \Msun$) adheres very well to the evolutionary scenario recently inferred by Kawamura et al. (2009) for GMCs in the Large Magellanic Cloud. A model cloud with $\Mmax \approx 2000 \Msun$ evolves in the Kennicutt-Schmidt diagram first passing through the locus of typical low- to-intermediate mass star-forming clouds, and then moving towards the locus of high-mass star-forming ones over the course of $\sim 10$ Myr. Also, the stellar age histograms for this cloud a few Myr before its destruction agree very well with those observed in the $\rho$-Oph stellar association, whose parent cloud has a similar mass, and imply that the SFR of the clouds increases with time. Our model thus agrees well with various observed properties of star-forming MCs, suggesting that the scenario of gravitationally collapsing MCs, with their SFR regulated by stellar feedback, is entirely feasible and in agreement with key observed properties of molecular clouds., Comment: Version accepted for publication in ApJ. At referee's suggestion, includes comparison with numerical models in addition to comparison with observational data

We numerically investigate the effect of feedback from the ionizing radiation heating from massive stars on the evolution of giant molecular clouds (GMCs) and their star formation efficiency (SFE). We find that the star-forming regions within the GMCs are invariably formed by gravitational contraction. After an initial period of contraction, the collapsing clouds begin forming stars, whose feedback evaporates part of the clouds' mass, opposing the continuing accretion from the infalling gas. The competition of accretion against dense gas consumption by star formation (SF) and evaporation by the feedback, regulates the clouds' mass and energy balance, as well as their SFE. We find that, in the presence of feedback, the clouds attain levels of the SFE that are consistent at all times with observational determinations for regions of comparable SF rates (SFRs). However, we observe that the dense gas mass is larger in general in the presence of feedback, while the total (dense gas + stars) is nearly insensitive to the presence of feedback, suggesting that the total mass is determined by the accretion, while the feedback inhibits mainly the conversion of dense gas to stars. The factor by which the SFE is reduced upon the inclusion of feedback is a decreasing function of the cloud's mass, for clouds of size ~ 10 pc. This naturally explains the larger observed SFEs of massive-star forming regions. We also find that the clouds may attain a pseudo-virialized state, with a value of the virial mass very similar to the actual cloud mass. However, this state differs from true virialization in that the clouds are the center of a large-scale collapse, continuously accreting mass, rather than being equilibrium entities., Comment: Submitted to ApJ (abstract abridged)

We present numerical simulations aimed at exploring the effects of varying the sub-grid physics parameters on the evolution and the properties of the galaxy formed in a low-mass dark matter halo (~7 times 10^10 Msun/h at redshift z=0). The simulations are run within a cosmological setting with a nominal resolution of 218 pc comoving and are stopped at z = 0.43. In all of our simulations, an extended old/intermediate-age stellar halo and a more compact younger stellar disk are formed. We found that a non negligible fraction of the halo stars are formed in situ in a spheroidal distribution. Changes in the sub-grid physics parameters affect significantly and in a complex way the evolution and properties of the galaxy: (i) Lower threshold densities nsf produce larger stellar effective radii Re, less peaked circular velocity curves V_c(R), and greater amounts of low-density and hot gas in the disk mid-plane; (ii) When stellar feedback is modeled by temporarily switching off radiative cooling in the star forming regions, Re increases (by a factor of ~ 2 in our particular model), the circular velocity curve becomes flatter, and a complex multi-phase gaseous disk structure develops; (iii) A more efficient local conversion of gas mass to stars, measured by a stellar particle mass distribution biased toward larger values, increases the strength of the feedback energy injection -driving outflows and inducing burstier SF histories; iv) If feedback is too strong, gas loss by galactic outflows -which are easier to produce in low-mass galaxies- interrupts SF, whose history becomes episodic. The simulations exhibit two important shortcomings: the baryon fractions are higher, and the specific SF rates are much smaller, than observationally inferred values for redshifts ~ 0.4-1., Comment: Accepted by ApJ after referee's revision (abstract abridged)

Based on the photometric properties of NGC 6822 we derive a robust star formation history. Adopting this history we compute 15 models of galactic chemical evolution. All of them match present-day photometric properties. The dark halo mass in all models evolves according to the mass assembly history predicted by a $\Lambda$CDM cosmology. We model the evolution of the baryonic mass aggregation history in this cosmological context assuming that part of the gas available for accretion never falls into the system due to two different physical processes. For seven models we assume that during accretion the universal baryon fraction is reduced by reionization only. The best model of this first group, a complex model with an early outflow, fits the observed gaseous mass, and the O/H, C/O, and Fe/O present-day values. This model requires a lower upper mass limit for the IMF than that of the solar vicinity, in agreement with recent results by other authors. We have also computed eight models where, in addition to reionization, the accreted baryon fraction is reduced by large-scale shock heating.The best model of this series, that also requires and early outflow and a lower upper mass limit for the IMF,can marginally fit the gaseous mass and the O/H and C/O observed values., Comment: 25 pages, ApJ accepted

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., Comment: 5 pages, 1 figure, Submitted to ApJL

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., Comment: 8 pages, 6 figures, Proceedings of "Cosmic Flows", Victoria, UVic, July 1999, paspconf and epsf

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., Comment: 13 pages, Latex, AASTeXv5 and natbib

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., Comment: Contributions to the Proceedings of the MPA/ESO Conference "Evolution of LSS: from Recombination to Garching", Garching, August 2-7, 1998. LaTeX (uses tp.sty, psfig.sty), 5 pages, 4 figs. Full-resolution version of Fig.2 is available at ftp://charon.nmsu.edu/pub/kravtsov/PAPERS/Garching98/

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), Comment: submitted to the Astrophys.Journal; 21 pages, LaTeX (uses emulateapj.sty); full resolution versions of figs.1 and 2 are available at http://astro.nmsu.edu/~akravtso/GROUP/group_publications.html or at ftp://charon.nmsu.edu/pub/kravtsov/PAPERS/Bias/

A newly developed isochrone synthesis algorithm for the photometric evolution of galaxies is described. Two initial mass functions, IMFs, in particular, the recent IMF determined by Kroupa, Tout, and Gilmore, three photometric transformations, and a 1-Gyr-burst star formation rate, SFR, are used to compute the $B-V$ and $V-K$ color index evolution. Non-negligible differences are observed among model results. In the framework of the galaxy count model by Col\'in, Schramm, and Peimbert a simple merging scenario is considered to account for the excess of galaxies observed in the blue band counts. The excess is explained by the number and luminosity evolution of a group of galaxies called interacting, I. It is assumed that the number of I galaxies increases as $(1+z)^{\eta}$ due to mergers. Moreover, it is proposed that their characteristic luminosity increases as $(1+z)^3$ due to starbursts driven by galaxy-galaxy collision and decreases as $(1+z)^{-\eta}$ due to the change in the size of the galaxies. Not much number evolution is needed to account for the excess; for example, a model with $\eta = 4.0$ predicts that about 17 \% of the galaxies at $z = 0.4$ are interacting. Number evolution models with a rather high value of $\eta$ fit better the data; in particular, the model with $\eta = 4.0$ predicts that about 13 \% of the galaxies have $z > 0.7$ in the $21.0 < m_{b_J} < 22.5$ interval, this contrasts with the upper bound of 5 \% obtained with the sample of 78 galaxies by Colless et al. The excess of high redshift galaxies can not be simply explained by changing reasonably the parameters of the luminosity function of I galaxies. This result could indicate that mergers are not the whole story. Our best-fit model produces the following values for the, Comment: 16 pages uuencoded gzip-compressed Postscript including 6 figures, to appear in RevMexAA